Short genome report | Open | Published:
Draft genome sequences for the obligate bacterial predators Bacteriovorax spp. of four phylogenetic clusters
Standards in Genomic Sciencesvolume 10, Article number: 11 (2015)
Bacteriovorax is the halophilic genus of the obligate bacterial predators, Bdellovibrio and like organisms. The predators are known for their unique biphasic life style in which they search for and attack their prey in the free living phase; penetrate, grow, multiply and lyse the prey in the intraperiplasmic phase. Bacteriovorax isolates representing four phylogenetic clusters were selected for genomic sequencing. Only one type strain genome has been published so far from the genus Bacteriovorax. We report the genomes from non-type strains isolated from aquatic environments. Here we describe and compare the genomic features of the four strains, together with the classification and annotation.
As a member of the highly diverse Deltaproteobacteria class, the obligate bacterial predators Bdellovibrio and like organisms possess unique ecological features that are worth exploring. They are the only known predatory bacteria that exhibit a life cycle alternating between an extracellular free-living phase and an intraperiplasmic phase and are capable of invading the periplasmic space of prey cells, resulting in the lysis of the prey and release of new progeny . Based on their small size, about 1/5th that of a typical bacterium cell, BALOs have been called “the world’s smallest hunters”. Nevertheless, their genomes are larger than expected, more than 3.98 Mb in Bdellovibrio. bacteriovorus Tiberius , 3.78 Mb in B. bacteriovorus HD100  and 3.44 Mb in Bacteriovorax marinus SJ . Despite the uniqueness , and increasing understanding, of the potential of these organisms in various applications [5–7], their phylogeny and unique predatory features are only beginning to be understood.
Systematics has played a most important role in advancing the study of the BALOs. Based on systematic genomic molecular techniques, the original BALO genus, Bdellovibrio, has been subdivided into four genera: Bdellovibrio,Bacteriolyticum,Peredibacter, andBacteriovorax[8–10]. Being an exclusive saltwater genus, Bacteriovorax is distinct from the freshwater/terrestrial members of BALOs in many ways. It is ubiquitous in salt-water environments , requires at least 0.5% NaCl for growth, prefers saltwater prey , thrive at a lower temperature range  and has a lower % GC ratio of ca. 37%  compared to the 50.65% of the freshwater Bdellovibrio bacteriovorus HD100. Currently, Bacteriovorax marinus SJ is the only strain from the genus Bacteriovorax of which the complete genome has been sequenced and reported.
To date, variations in the 16S rRNA sequences have yielded approximately eight Bacteriovorax clusters or OTUs. The previously sequenced Bacteriovorax marinus SJT is one of the representatives that belong to phylogenetic Cluster III. This classification scheme has enabled for the first time the detection of specific Bacteriovorax strains in environmental/ecological studies. The validity of using the 16S rRNA gene was tested by comparison with the rpoB gene . The results of recent studies monitoring the activities and distribution of specific phylogenetic clusters have yielded new discoveries on the distribution, predation patterns, prey preferences, and ecology of this bacterial predator [14–16].
Here we present a description of the draft genomes of Bacteriovorax isolates of four phylogenetic clusters isolated from estuarine systems, together with the description of the genomic sequencing and annotation.
A 16S rRNA phylogenetic tree was constructed showing the phylogenetic neighborhood of the four newly sequenced Bacteriovorax strains within the family of Bdellovibrionaceae (Figure 1). As expected, Bacteriovoraxsp. strain BSW11_IV was grouped together with cluster IV, strain SEQ25 _V with cluster V, Strain DB6_IX with Cluster IX and lastly strain BAL6_X with cluster X.
General features of Bacteriovoraxspp. are summarized in Table 1. Individual features of Bacteriovorax isolates have not been sufficiently explored and are largely unknown. Micrographs generated by both transmission electron microscopy and scanning microscopy (Figure 2) suggest that Bacteriovoraxspp. employ similar predation strategies as other BALO members to attack and reside in the periplasic space of its prey.
Genome sequencing information
Genome project history
The four genomes were selected for sequencing on the basis of their phylogenetic position and isolation source. Low salt Bacteriovoraxsp. BSW11_IV was isolated from Breton Sound, Louisiana (salinity 0.6 ppt; Temperature 26.4°C) and SEQ25_V was obtained from water samples of Barataria Bay, Louisiana (salinity 5.2 ppt; Temperature 19.2°C). High salt DB6_IX (Salinity 32.4 ppt; Temperature 24.1°C) and BAL6_ X (Salinity 30.9 ppt; Temperature 25.2°C) were obtained from Apalachicola Bay, Florida. The genome sequences were deposited in GenBank. Sequencing and annotation were performed at the J. Craig Venter Institute. Table 2 presents the project information and its association with MIGS version 2.0 compliance .
Growth conditions and DNA isolation
Bacteriovorax cultures were grown separately in 70% artificial sea water (ASW) (Instant Ocean, Aquarium Systems, Inc., Mentor, Ohio) (pH 8, salinity 22 ppt.) amended with prey, Vibrio. vulnificus CMCP6 (for Bx sp.BSW11_IV and SEQ25_V), or V. parahaemolyticusRIMD 2210633 (for Bacteriovoraxsp. DB6_IX and BAL6_X). The genomes of both prey bacteria have been sequenced previously [28, 29]. When cultures became clear (2–3 days after inoculation of the prey), which indicated the majority of the prey cells were lysed by the predators, 300 ml suspensions were filtered consecutively through 0.45 and 0.22 μm sterile syringe filters (Corning, NY, USA) to remove any remaining prey. Filtrates containing high concentrations of Bacteriovorax cells (ca. 4 × 108 PFU ml−1) were centrifuged at 27,485 × g for 20 min at 4°C. The pellets were then re-suspended in 1 ml of ASW respectively. To test that the concentrated Bacteriovorax suspensions were free of prey cell contamination, aliquots of 0.1 ml of the filtrate were spread-plated onto LB agar and incubated at 37°C for two days.
Subsequently, total DNA from the cell pellets were extracted using the QIAGEN Kit (QIAamp DNA Mini Kit), according to the manufacturer’s protocol. The concentration and purity of DNA was measured by a NanoDrop Spectrophotometer (ND 1000, Thermo Fisher Scientific, DE). To reconfirm the phylotype of the isolations, the DNA was PCR amplified using Bacteriovorax specific primers, Bac-676 F (5′-ATT TCG CAT GTA GGG GTA-3′) and Bac-1442R (5′-GCC ACG GCT TCA GGT AAG-3′)  by puReTaq Ready-To-Go PCR Beads (GE Healthcare Bio-Sciences). PCR products were purified with the QIAquick PCR-Purification Kit (QIAGEN) and sequenced with Bac-676 F primer at the DNA Sequencing Laboratory at Florida State University. 16S DNA sequences were analyzed with the Basic Local Alignment Search Tool (BLAST) server from the National Center of Biotechnology Information .
Genome sequencing and assembly
Genome sequencing of the four Bacteriovorax isolates was conducted at the J. Craig Venter Institute employing a combination of Illumina and 454 sequencing platforms. The 454 data consisted of a half plate of 454 FLX per genome from 3 KB mate paired libraries. The Illumina data consisted of one-quarter lane of 2 × 100 bp Illumina HiSeq data per genome. On average, 300,000 454 reads (average length trimmed 300 bp) and 10 million Illumina sequences (average length trimmed 100 bp) were generated per genome. To incorporate a hybrid assembly using both 454 and Illumina sequence libraries, one million reads were randomly sampled (with their mates) from the Illumina library using Celera , which was sufficient to provide high coverage in the initial assemblies.
Genes were identified using GLIMMER3  as part of the JCVI prokaryotic annotation pipeline followed by manual curation using the Manatee annotation-editing platform. The JCVI automated pipeline incorporates HMM3 searches against Pfam  and TIGRFAMs  and BLASTP against UniProt , JCVI’s database of experimentally characterized proteins CharProt DB , and PIR .
The Bacteriovoraxsp. BSW11_ IV draft genome contains 3,650,096 bp with a GC content of 37%. The hybrid assembly was scanned for contamination using BlastP and the appropriate contigs were filtered out. The final assembly comprised of 3 scafolds, 30 RNAs and 3457 CDS. For the CDSs, 2591 (75%) proteins had a BLASTP hit with an e-value of 1e-9 or better to Bacteriovorax marinus SJ, and an additional 151 (4%) CDSs had a hit within the genus Bdellovibrio.
The Bacteriovoraxsp. SEQ25_V draft genome contains 3,450,786 bp with a GC content of 37%. The sequences were assembled into 29 contigs comprised of 35 RNAs and 3,292 CDSs. Among the CDSs, 2,456 (75%) of proteins had a BLASTP hit with an e-value of 1e-9 or better to Bacteriovorax marinusSJ, and an additional 131 (4%) CDSs had a hit within the genus Bdellovibrio.
The Bacteriovoraxsp. DB6_IX draft genome contains 2,969,235 bp with a GC content of 38%. Sequences were assembled into 10 scaffolds with 30 RNAs and 3192 CDSs. Among the CDSs 2,253 (71%) proteins had a BLASTP hit with an e-value of 1e-9 or better to Bx marinus SJ, and an additional 97 (3%) CDSs had a hit within the genus Bdellovibrio.
The Bacteriovoraxsp. BAL6_ X draft genome contains 3,233,679 bp with a GC content of 36%. The reads were assembled into 9 contigs with 37 RNAs and 3,065 CDSs. Among the CDSs, 2,298 (72%) proteins had a BLASTP hit with an e-value of 1e-9 or better to Bacteriovorax marinusSJ, and an additional 92 (3%) CDSs had a hit within the genus Bdellovibrio.
It is noteworthy to point out that three phage tail fiber proteins were identified within the Bx sp. BSW11_IV genome but were absent from all the other BALO genome including the completed Bacteriovorax marinus SJ and Bdellovibrio bacteriovorus HD100 genomes. A staphylococcal phi-Mu50B-like prophage element was present in both SJ and HD 100 genomes but was not found in the genomes of the four newly sequenced Bacteriovorax isolates. The properties and the statistics of the genome are summarized in Tables 3, 4 and 5 and (Additional file 1: Table S1).
Insights from the genome sequences
Genome Comparisons between BALO Members
Crossman et al.,  reported that the genomic sequences of Bacteriovorax marinus SJ were unique with about one third of predicted genes over 500 bp in length having no significant hit in the databases. No genomic synteny was found between SJ and its closest whole genome sequenced relative at that time, Bdellovibrio bacteriovorus HD100.
We found that even within the genus Bacteriovorax, the genomic sequences were highly divergent with an average identity of 70%. A Venn diagram summarizing the comparison of the four Bacteriovorax isolates is presented in Figure 3. As shown in the diagram, a core of 1,513 proteins is shared by all four Bacteriovorax genomes and each encodes many proteins without orthologs in the other three (Figure 3A). When compared to the freshwater/terrestrial Bdellovibrio bacteriovorus HD 100, only a total of 843 genes were shared between all BALO members (Figure 3B). The calculated ANI  for BALO members (Additional file 2: Table S2) is below 75%, which is the threshold for the scores to be reliable. The AAI among the five Bacteriovorax genomes ranged between 50% to 60% (Additional file 3: Table S3), also significantly lower than the typical values found for species within a genus (73%-99.5%) [42, 43]. Currently, several proposals to clarify and revise the systematics of BALOs are under consideration.
Comparisons of BALOs and non-predatory bacteria
Phylogenetically, most genera of BALOs (including Bacteriovorax) are classified as Deltaproteobacteria. Members of this class are found in diverse environments with various lifestyles such as Myxococcus xanthus which is characterized by its gliding motility and wolf pack predatory strategy to prey on other bacteria , Pelobacter carbinolicus which grows by using iron and sulfur as electron acceptors , and the focus of this study, the obligate predators Bacteriovorax spp. which replicate within the periplasmic space of prey bacteria. Although their ecological features are distinct, the genomes of Deltaproteobacteria were found to exhibit some common characteristics. For example, most Deltaproteobacteria, including the Bdellovibrio bacteriovorus HD 100, typically possess two giant S1 ribosomal protein genes and high numbers of TonB receptors and ferric siderophore receptors which facilitate metal uptake and removal . In contrast, only one giant S1 protein was found in the Bacteriovorax marinus SJ genome , and our study confirmed that this is the case for the other four Bacteriovorax genomes. Bacteriovorax genomes also encodes multiple TonB receptor proteins (6–11 copies) and ferric siderophore receptors (2–4 copies) that they may use for predation.
Using a reciprocal best match analysis with e-value cutoff of 10–9, 843 core genes were found to contain orthologs in all six BALO genomes including previously sequenced SJ and HD 100 genomes (see center of Figure 3B). Fifty nine of these genes (Additional file 4: Table S4) have no homologs with an E-value of 10–9 or lower to proteins from any non-predatory bacterium in the NCBI “nr” database. These genes, including periplasmic proteins, a radical activating enzyme and an outer membrane channel protein, may represent a core set of unique genes involved in the predatory process and prey interactions such as locating the prey, degradation and consumption of prey cellular content, formation of bdelloplast, synchronous nonbinary septation or release of progeny from the ghost cell.
The genomes of four Bacteriovorax phylogenetic clusters isolated from the environment were sequenced. The genome sizes of the four strains were comparable with Bacteriovorax SJ but were slightly smaller than the two freshwater BALOs, Bdellovibrio bacteriovorus Tiberius and B. bacteriovorus HD100. Fifty-nine genes were identified that are conserved among BALOs, but not present in other organisms, that may be responsible for their predatory life style. The unique genomic features of Bacteriovorax that are essential for their ecological function were also reported.
Bdellovibrio and Like Organisms
Average nucleotide identity
Average amino acid identity.
Rendulic S, Jagtap P, Rosinus A, Eppinger M, Baar C, Lanz C: A predator unmasked: life cycle of Bdellovibrio bacteriovorus from a genomic perspective. Science 2004, 303:689–92. 10.1126/science.1093027
Hobley L, Lerner TR, Williams LE, Lambert C, Till R, Milner DS, et al.: Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria. BMC Genomics 2012, 13:670. 10.1186/1471-2164-13-670
Baer ML, Ravel J, Pineiro SA, Guether-Borg D, Williams HN: Reclassification of salt-water Bdellovibrio sp. as Bacteriovorax marinus sp. nov. and Bacteriovorax litoralis sp. nov. Int J Syst Evol Microbiol 2004, 54:1011–6. 10.1099/ijs.0.02458-0
Chen H, Williams HN: Sharing of prey: coinfection of a bacterium by a virus and a prokaryotic predator. MBio 2012,3(2):e00051–12.
Atterbury RJ, Hobley L, Till R, Lambert C, Capeness MJ, Lerner TR, et al.: Effects of orally administered Bdellovibrio bacteriovorus on the well-being and Salmonella colonization of young chicks. Appl Environ Microbiol 2011, 77:5794–803. 10.1128/AEM.00426-11
Dashiff A, Junka R, Libera M, Kadouri D: Predation of human pathogens by the predatory bacteria Micavibrio aeruginosavorus and Bdellovibrio bacteriovorus . J Appl Microbiol 2011, 110:431–44. 10.1111/j.1365-2672.2010.04900.x
Cao H, Hou S, He S, Lu L, Yang X: Identification of a Bacteriovorax sp. isolate as a potential biocontrol bacterium against snakehead fish‒pathogenic Aeromonas veronii . J Fish Dis 2013,37(3):283–289.
Baer ML, Ravel J, Chun J, Hill RT, Williams HN: A proposal for the reclassification of Bdellovibrio stolpii and Bdellovibrio starrii into a new genus, Bacteriovorax gen. nov. as Bacteriovorax stolpii comb. nov. and Bacteriovorax starrii comb. nov., respectively. Int J Syst Evol Microbiol 2000, 50 Pt 1:219–24.
Davidov Y, Jurkevitch E: Diversity and evolution of Bdellovibrio -and-like organisms (BALOs), reclassification of Bacteriovorax starrii as Peredibacter starrii gen. nov., comb. nov., and description of the Bacteriovorax-Peredibacter clade as Bacteriovoracaceae fam. nov. Int J Syst Evol Microbiol 2004, 54:1439–52. 10.1099/ijs.0.02978-0
Pineiro SA, Williams HN, Stine OC: Phylogenetic relationships amongst the saltwater members of the genus Bacteriovorax using rpoB sequences and reclassification of Bacteriovorax stolpii as Bacteriolyticum stolpii gen. nov., comb. nov. Int J Syst Evol Microbiol 2008, 58:1203–9. 10.1099/ijs.0.65710-0
Schoeffield AJ, Williams HN: Efficiencies of Recovery of Bdellovibrios from Brackish- Water Environments by Using Various Bacterial Species as Prey. Appl Environ Microbiol 1990, 56:230–6.
Williams HN: Cultural, immunologic and ecologic studies of marine bdellovibrios isolated from the Atlantic Ocean and the Chesapeake Bay. Baltimore, MD, USA: University of Maryland; 1979.
Marbach A, Varon M, Shilo M: Properties of marine bdellovibrios. Microb Ecol 1976, 2:284–95.
Chen H, Young S, Berhane T-K, Williams HN: Predatory Bacteriovorax Communities Ordered by Various Prey Species. PLoS ONE 2012, 7:e34174. 10.1371/journal.pone.0034174
Chauhan A, Cherrier J, Williams HN: Impact of sideways and bottom-up control factors on bacterial community succession over a tidal cycle. Proc Natl Acad Sci U S A 2009, 106:4301–6. 10.1073/pnas.0809671106
Chen H, Athar R, Zheng G, Williams HN: Prey bacteria shape the community structure of their predators. ISME J 2011, 5:1314–22. 10.1038/ismej.2011.4
Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML web servers. Syst Biol 2008, 57:758–71. 10.1080/10635150802429642
Pineiro S, Stine O, Chauhan A, Steyert S, Smith R, Williams H: Global survey of diversity among environmental saltwater Bacteriovoracaceae . Environ Microbiol 2007, 9:2441–50. 10.1111/j.1462-2920.2007.01362.x
Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, et al.: The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol 2008, 26:541–7. 10.1038/nbt1360
Woese CR, Kandler O, Wheelis ML: Towards a natural system of organisms: proposal for the domains Archaea , Bacteria, and Eucarya. Proc Natl Acad Sci U S A 1990, 87:4576–9. 10.1073/pnas.87.12.4576
Garrity GM, Holt JG: The Road Map to the Manual. In Bergey’s Manual of Systematic Bacteriology. Volume 1. 2nd edition. Edited by: Garrity GM, Boone DR, Castenholz RW. New York: Springer; 2001:119–69.
Validation List No. 107: List of new names and new combinations previously effectively, but not validly, published Int J Syst Evol Microbiol 2006, 56:1–6.
Kuever J, Rainey FA, Widdel F: Class VI. Deltaproteobacteria class nov. In Bergey’s Manual of Systematic Bacteriology. Volume 2, part C. 2nd edition. Edited by: Brenner DJ, Krieg NR, Staley JT, Garrity GM. New York: Springer; 2005.
Garrity GM, Bell JA, Lilburn T: Order VII. Bdellovibrionales ord. nov. In Bergey’s Manual of Systematic Bacteriology. Volume 2, Part C. 2nd edition. Edited by: Brenner DJ, Krieg NR, Staley JT, Garrity GM. New York: Springer; 2005:1040.
Garrity GM, Bell JA, Lilburn T: Family I. Bdellovibrionaceae fam. nov. In Bergey’s Manual of Systematic Bacteriology. Volume 2, Part C. 2nd edition. Edited by: Brenner DJ, Krieg NR, Staley JT, Garrity GM. New York: Springer; 2005:1040–1.
Williams H, Baer M, Tudor J: Bdellovibrio Stolp and Starr 1963, 243 AL . Bergey’s Manual® of Systematic Bacteriology 2005, 1041–53.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25:25–9. 10.1038/75556
Makino K, Oshima K, Kurokawa K, Yokoyama K, Uda T, Tagomori K, et al.: Genome sequence of Vibrio parahaemolyticus : a pathogenic mechanism distinct from that of V. cholerae . Lancet 2003, 361:743–9. 10.1016/S0140-6736(03)12659-1
Kim YR, Lee SE, Kim CM, Kim SY, Shin EK, Shin DH, et al.: Characterization and pathogenic significance of Vibrio vulnificus antigens preferentially expressed in septicemic patients. Infect Immun 2003, 71:5461. 10.1128/IAI.71.10.5461-5471.2003
Davidov Y, Friedjung A, Jurkevitch E: Structure analysis of a soil community of predatory bacteria using culture-dependent and culture-independent methods reveals a hitherto undetected diversity of Bdellovibrio -and-like organisms. Environ Microbiol 2006, 8:1667–73. 10.1111/j.1462-2920.2006.01052.x
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990,215(3):403–410. 10.1016/S0022-2836(05)80360-2
Miller JR, Delcher AL, Koren S, Venter E, Walenz BP, Brownley A, et al.: Aggressive assembly of pyrosequencing reads with mates. Bioinformatics 2008, 24:2818–24. 10.1093/bioinformatics/btn548
Delcher AL, Harmon D, Kasif S, White O, Salzberg SL: Improved microbial gene identification with GLIMMER. Nucleic Acids Res 1999, 27:4636–41. 10.1093/nar/27.23.4636
Eddy SR: Accelerated profile HMM searches. PLoS Comput Biol 2011, 7:e1002195. 10.1371/journal.pcbi.1002195
Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res 2011, 39:W29–37. 10.1093/nar/gkr367
Haft DH, Selengut JD, Richter RA, Harkins D, Basu MK, Beck E: TIGRFAMs and genome properties in 2013. Nucleic Acids Res 2013, 41:D387–95. 10.1093/nar/gks1234
Magrane M: UniProt Knowledgebase: a hub of integrated protein data. Database: J Biol Dat Curation 2011, 2011:..
Madupu R, Richter A, Dodson RJ, Brinkac L, Harkins D, Durkin S, et al.: CharProtDB: a database of experimentally characterized protein annotations. Nucleic Acids Res 2011,40(D1):D237.
Barker WC, Garavelli JS, Huang H, McGarvey PB, Orcutt BC, Srinivasarao GY, et al.: The Protein Information Resource (PIR). Nucleic Acids Res 2000, 28:41–4. 10.1093/nar/28.1.41
Crossman LC, Chen H, Cerdeño-Tárraga A-M, Brooks K, Quail MA, Pineiro SA, et al.: A small predatory core genome in the divergent marine Bacteriovorax marinus SJ and the terrestrial Bdellovibrio bacteriovorus . ISME J 2012, 7:148–60.
Konstantinidis KT, Tiedje JM: Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci U S A 2005, 102:2567–72. 10.1073/pnas.0409727102
Konstantinidis K, Tiedje J: Towards a genome-based taxonomy for prokaryotes. J Bacteriol 2005, 187:6258. 10.1128/JB.187.18.6258-6264.2005
Thompson C, Vicente A, Souza R, Vasconcelos A, Vesth T, Alves N, et al.: Genomic taxonomy of vibrios. BMC Evol Biol 2009, 9:258. 10.1186/1471-2148-9-258
Reichenbach H: The ecology of the myxobacteria. Environ Microbiol 1999, 1:15–21. 10.1046/j.1462-2920.1999.00016.x
Lovley DR, Phillips E, Lonergan DJ, Widman PK: Fe (III) and S 0 reduction by Pelobacter carbinolicus . Appl Environ Microbiol 1995, 61:2132.
Karlin S, Brocchieri L, Mrázek J, Kaiser D: Distinguishing features of δ-proteobacterial genomes. Proc Natl Acad Sci U S A 2006, 103:11352. 10.1073/pnas.0604311103
This work was financially supported by grants from the National Science Foundation HBCU-RISE (#0531523) and DDIG (DEB-1110620). We thank the assistance of Drs Jill W Verlander and Sharon W Matthews at the University of Florida, College of Medicine Electron Microscopy Facility; Megan Lamb at Apalachicola Bay, National Estuarine Research Reserve for research vessel operation and Dr. Aixin Hou at the Louisiana State University for assisting in sample collection and providing laboratory space for initial processing of samples.
The authors declare that they have no competing interests.
JHB and HNW initiated and supervised the study. HC draft the manuscript, conducted wetlab work and performed electron microscopy. HC, LMB, DSL, TLD and NG annotated the genome. HC, PM, LMB and JHB worked on genome sequencing and assembly. HC, NL, JHB, PM and HNW discussed, analyzed the data and revised the manuscript. All authors read and approved the final manuscript.