Skip to main content

Advertisement

High-quality permanent draft genome sequence of Bradyrhizobium sp. Tv2a.2, a microsymbiont of Tachigali versicolor discovered in Barro Colorado Island of Panama

Abstract

Bradyrhizobiumsp. Tv2a.2 is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated from an effective nitrogen-fixing root nodule of Tachigali versicolor collected in Barro Colorado Island of Panama. Here we describe the features of Bradyrhizobiumsp. Tv2a.2, together with high-quality permanent draft genome sequence information and annotation. The 8,496,279 bp high-quality draft genome is arranged in 87 scaffolds of 87 contigs, contains 8,109 protein-coding genes and 72 RNA-only encoding genes. This rhizobial genome was sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.

Introduction

Legumes engage in nitrogen-fixation symbioses with bacterial partners from at least 13 genera of Proteobacteria [1-4]. Despite the high extent of phylogenetic diversity of root nodule bacteria, the very broad distribution of one particular genus ( Bradyrhizobium ) across host legume clades suggests that bacteria in this genus may have been the first legume symbionts [5]. Bradyrhizobium interacts with the widest diversity of legume clades (at least 24 of ca. 33 nodule-forming legume tribes; [6]) and is associated with nodulating groups that represent early branching lineages [7] in all three legume subfamilies [8,9]. Analysis of basal Bradyrhizobium lineages that are associated with early-diverging legume groups may thus shed light on the origins of this symbiosis.

Here we report the genome sequence of one such organism, Bradyrhizobium strain Tv2a.2. Strain Tv2a.2 was sampled in 1997 from the tree Tachigali versicolor on Barro Colorado Island, Panama, a biological preserve with an old-growth moist tropical forest [10]. Tachigali is one of just a handful of nodule-forming genera in the legume Subfamily Caesalpinioideae [11], which is comprised of the earliest branching lineages in the legume family [7]. Tachigali versicolor is a large canopy tree with an unusual monocarpic life history, in which trees grow for decades without flowering. They produce just a single crop of seeds, and then die [12].

Strain Tv2a.2 is a typical representative of the nodule symbionts that are associated with Tachigali in this tropical forest habitat [13], and appears to represent a unique early-diverging lineage of Bradyrhizobium . Phylogenetic analyses have placed Tv2a.2 somewhere near the early split in the genus between two large superclades represented by B. diazoefficiens USDA 110 and B. elkanii USDA 76. However, its exact position near the base of the Bradyrhizobium tree varies to some extent in different analyses, depending on the loci, the strains included, and the method of tree analysis [5,13]. For example, a Bayesian analysis of 16S rRNA sequences from the type strains of 21 Bradyrhizobium species and strain ORS278 placed Tv2a.2 as the earliest diverging Bradyrhizobium lineage [14].

Here we provide an analysis of the complete genome sequence of Tv2a.2, one of the rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project proposal [15], whose properties should help to clarify early events in the diversification of the genus Bradyrhizobium as a whole.

Organism information

Classification and features

Bradyrhizobium sp. Tv2a.2 is a motile, non-sporulating, non-encapsulated, Gram-negative strain in the order Rhizobiales of the class Alphaproteobacteria . The rod shaped form (Figure 1 Left, Center) has dimensions of approximately 0.5 μm in width and 1.5-2.0 μm in length. It is relatively slow growing, forming colonies after 6–7 days when grown on half strength Lupin Agar (½LA) [16], tryptone-yeast extract agar (TY) [17] or a modified yeast-mannitol agar (YMA) [18] at 28°C. Colonies on ½LA are opaque, slightly domed and moderately mucoid with smooth margins (Figure 1 Right).

Figure 1
figure1

Images of Bradyrhizobium sp. Tv2a.2 using scanning (Left) and transmission (Center) electron microscopy as well as light microscopy to visualize colony morphology on solid media (Right).

Figure 2 shows the phylogenetic relationship of Bradyrhizobium sp. Tv2a.2 in a 16S rRNA gene sequence based tree. This strain is phylogenetically the most related to Bradyrhizobiumsp. EC3.3 based on a 16S rRNA gene sequence identity of 99.31% as determined using BLAST analysis [19]. Tv2a.2 is also related to the type strains Bradyrhizobium ingae BR 10250T and Bradyrhizobium iriomotense EK05T with 16S rRNA gene sequence identities of 99.16 % and 99.08%, respectively, based on results from the EzTaxon-e server [20,21].

Figure 2
figure2

Phylogenetic tree highlighting the position of Bradyrhizobium sp. Tv2a.2 (shown in blue print) relative to other type and non-type strains in the Bradyrhizobium genus using a 1,310 bp intragenic sequence of the 16S rRNA gene. Azorhizobium caulinodans ORS 571T sequence was used as an outgroup. All sites were informative and there were no gap-containing sites. Phylogenetic analyses were performed using MEGA, version 5.05 [41]. The tree was built using the maximum likelihood method with the General Time Reversible model. Bootstrap analysis with 500 replicates was performed to assess the support of the clusters. Type strains are indicated with a superscript T. Strains with a genome sequencing project registered in GOLD [22] have the GOLD ID mentioned after the strain number and are represented in bold, otherwise the NCBI accession number is provided.

Minimum Information about the Genome Sequence (MIGS) of Tv2a.2 is provided in Table 1 and Additional file 1: Table S1.

Table 1 Classification and general features of Bradyrhizobium sp. Tv2a.2 in accordance with the MIGS recommendations [42] published by the Genome Standards Consortium [43]

Symbiotaxonomy

Bradyrhizobium strain Tv2a.2 was isolated from nodules of Tachigali versicolor found in a tropical forest on Barro Colorado Island, Panama [10]. Due to the highly erratic pattern of seed production from this host, no seeds of this legume were available to authenticate the symbiotic proficiency of strain Tv2a.2. Nodulation and nitrogen fixation was therefore tested on two promiscuous legumes ( Vigna unguiculata , Macroptilium atropurpureum ) and revealed that nodules could only develop on M. atropurpureum. Acetylene reduction assays also showed that these nodules lacked nitrogenase activity [13]. A further indication that Tv2a.2 may be relatively host-specific is the fact that extensive sampling of other legume hosts in Panama (and elsewhere in the Neotropics) have never recovered strains belonging to the Tv2a.2 lineage from any legume taxa other than T. versicolor [9].

Genome sequencing and annotation information

Genome project history

This organism was selected for sequencing on the basis of its environmental and agricultural relevance to issues in global carbon cycling, alternative energy production, and biogeochemical importance, and is part of the Genomic Encyclopedia of Bacteria and Archaea, Root Nodulating Bacteria (GEBA-RNB) project at the U.S. Department of Energy, Joint Genome Institute (JGI). The genome project is deposited in the Genomes OnLine Database [22] and a high-quality permanent draft genome sequence in IMG [23]. Sequencing, finishing and annotation were performed by the JGI using state of the art sequencing technology [24]. A summary of the project information is shown in Table 2.

Table 2 Project information

Growth conditions and genomic DNA preparation

Bradyrhizobium sp. Tv2a.2 was cultured to mid logarithmic phase in 60 ml of TY rich media on a gyratory shaker at 28°C [25]. DNA was isolated from the cells using a CTAB (Cetyl trimethyl ammonium bromide) bacterial genomic DNA isolation method [26].

Genome sequencing and assembly

The draft genome of Bradyrhizobium sp. Tv2a.2 was generated at the DOE Joint Genome Institute (JGI) using the Illumina technology [27]. An Illumina standard shotgun library was constructed and sequenced using the Illumina HiSeq 2000 platform which generated 8,336,316 reads totaling 1250.45 Mbp. All general aspects of library construction and sequencing were performed at the JGI and details can be found on the JGI website [28]. All raw Illumina sequence data was passed through DUK, a filtering program developed at JGI, which removes known Illumina sequencing and library preparation artifacts (Mingkun L, Copeland A, Han J, Unpublished). Following steps were then performed for assembly: (1) filtered Illumina reads were assembled using Velvet (version 1.1.04) [29], (2) 1–3 Kbp simulated paired end reads were created from Velvet contigs using wgsim [30], (3) Illumina reads were assembled with simulated read pairs using Allpaths–LG (version r39750) [31]. Parameters for the assembly steps were 1) velveth: −-v --s 51 --e 71 --i 2 --t 1 --f “-shortPaired -fastq $FASTQ” --o “-ins_length 250 -min_contig_lgth 500” for Velvet and 2) wgsim: −e 0–1 76–2 76 -r 0 -R 0 -X 0. The final draft assembly contained 87 contigs in 87 scaffolds. The total size of the genome is 8.5 Mb with an average of 109.04x coverage of the genome.

Genome annotation

Genes were identified using Prodigal [32], as part of the DOE-JGI genome annotation pipeline [33,34]. The predicted CDSs were translated and used to search the National Center for Biotechnology Information non-redundant database, UniProt, TIGRFam, Pfam, KEGG, COG, and InterPro databases. The tRNAScanSE tool [35] was used to find tRNA genes, whereas ribosomal RNA genes were found by searches against models of the ribosomal RNA genes built from SILVA [36]. Other non–coding RNAs such as the RNA components of the protein secretion complex and the RNase P were identified by searching the genome for the corresponding Rfam profiles using INFERNAL [37]. Additional gene prediction analysis and manual functional annotation was performed within the Integrated Microbial Genomes-Expert Review (IMG-ER) system [38] developed by the Joint Genome Institute, Walnut Creek, CA, USA.

Genome properties

The genome is 8,496,279 nucleotides with 62.20% GC content (Table 3) and comprised of 87 scaffolds. From a total of 8,181 genes, 8,109 were protein encoding and 72 RNA only encoding genes. The majority of genes (72.94%) were assigned a putative function whilst the remaining genes were annotated as hypothetical. The distribution of genes into COGs functional categories is presented in Table 4.

Table 3 Genome statistics for Bradyrhizobium sp. Tv2a.2
Table 4 Number of genes associated with the general COG functional categories

Conclusions

Bradyrhizobium sp. Tv2a.2 was collected in 1997 from a nodule of the tree Tachigali versicolor on Barro Colorado Island, Panama. Based on 16S rRNA gene analyses, Tv2a.2 is phylogenetically the most closely related to Bradyrhizobium sp. EC3.3 (a strain isolated from a nodule of Erythrina costaricensis collected from Barro Colorado Island, Panama) and to the type strains Bradyrhizobium ingae BR 10250T and Bradyrhizobium iriomotense EK05T isolated from Inga laurina (Sw.) Willd. growing in the Cerrado Amazon region, State of Roraima, Brazil [39] and from Entada koshunensis, a legume available in Okinawa, Japan [40], respectively. Strain Tv2a.2 is one of 25 Bradyrhizobium genomes that were sequenced within the GEBA-RNB project [15]; of these, the Tv2a.2 genome has the fifth lowest genome size (8.5 Mbp), gene count (8,181) and Pfam percentage (74.32%) amongst these strains. The specific genome attributes of Bradyrhizobium sp. Tv2a.2 compared to the other Bradyrhizobium genomes will be important to understand the interactions required for the successful establishment of an effective symbiosis with the host Tachigali versicolor .

Abbreviations

GEBA-RNB:

Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria

JGI:

Joint Genome Institute

½LA:

Half strength Lupin Agar

TY:

Tryptone yeast

YMA:

Yeast mannitol agar

CTAB:

Cetyl trimethyl ammonium bromide

References

  1. 1.

    Willems A. The taxonomy of rhizobia; an overview. Plant Soil. 2006;287:3–14.

  2. 2.

    Gyaneshwar P, Hirsch AM, Moulin L, Chen WM, Elliott GN, Bontemps C, et al. Legume-nodulating betaproteobacteria: diversity, host range, and future prospects. Mol Plant Microbe In. 2011;24:1276–88.

  3. 3.

    Ardley JK, Parker MA, De Meyer SE, Trengove RD, O’Hara GW, Reeve WG, et al. Microvirga lupini sp. nov., Microvirga lotononidis sp. nov. and Microvirga zambiensis sp. nov. are alphaproteobacterial root-nodule bacteria that specifically nodulate and fix nitrogen with geographically and taxonomically separate legume hosts. Int J Syst Evol Microbiol. 2012;62:2579–88.

  4. 4.

    Maynaud G, Willems A, Soussou S, Vidal C, Maure L, Moulin L, et al. Molecular and phenotypic characterization of strains nodulating Anthyllis vulneraria in mine tailings, and proposal of Aminobacter anthyllidis sp. nov., the first definition of Aminobacter as legume-nodulating bacteria. Syst Appl Microbiol. 2012;35:65–72.

  5. 5.

    Parker MA. The spread of Bradyrhizobium lineages across host legume clades: from Abarema to Zygia. Microb Ecol. 2014;69:630–640.

  6. 6.

    Lewis GP. Legumes of the world. Royal Botanic Gardens, Kew: Richmond, UK; 2005.

  7. 7.

    Doyle JJ. Phylogenetic perspectives on the origins of nodulation. Mol Plant Microbe Int. 2011;24:1289–95.

  8. 8.

    Fonseca MB, Peix A, de Faria SM, Mateos PF, Rivera LP, Simoes-Araujo JL, et al. Nodulation in Dimorphandra wilsonii Rizz. (Caesalpinioideae), a threatened species native to the Brazilian Cerrado. PLoS One. 2012;7:e49520.

  9. 9.

    Parker MA, Rousteau A. Mosaic origins of Bradyrhizobium legume symbionts on the Caribbean island of Guadeloupe. Mol Phylogenet Evol. 2014;77:110–5.

  10. 10.

    Foster RB, Brokaw NVL. Structure and history of the vegetation on Barro Colorado Island. The ecology of a tropical forest: seasonal rhythms and long-term changes. Smithsonian Institution: Washington, DC; 1996.

  11. 11.

    Sprent JI, Ardley JK, James EK. From North to South: a latitudinal look at legume nodulation processes. South Afr J Bot. 2013;89:31–41.

  12. 12.

    Foster RB. Tachigali versicolor is a suicidal neotropical tree. Nature. 1977;268:624–26.

  13. 13.

    Parker MA. Divergent Bradyrhizobium symbionts on Tachigali versicolor from Barro Colorado Island, Panama. Syst Appl Microbiol. 2000;23:585–90.

  14. 14.

    Giraud E, Moulin L, Vallenet D, Barbe V, Cytryn E, Avarre JC, et al. Legumes symbioses: absence of Nod genes in photosynthetic bradyrhizobia. Science. 2007;316:1307–12.

  15. 15.

    Reeve W, Ardley J, Tian R, Eshragi L, Yoon JW, Ngamwisetkun P, et al. A genomic encyclopedia of the root nodule bacteria: assessing genetic diversity through a systematic biogeographic survey. Stand Genomic Sci. 2015;10:14.

  16. 16.

    Howieson JG, Ewing MA, D’antuono MF. Selection for acid tolerance in Rhizobium meliloti. Plant Soil. 1988;105:179–88.

  17. 17.

    Beringer JE. R factor transfer in Rhizobium leguminosarum. J Gen Microbiol. 1974;84:188–98.

  18. 18.

    Vincent JM. A manual for the practical study of the root-nodule bacteria. International biological programme. UK: Blackwell Scientific Publications, Oxford; 1970.

  19. 19.

    NCBI BLAST [http://blast.ncbi.nlm.nih.gov/Blast.cgi]

  20. 20.

    Kim O-S, Cho Y-J, Lee K, Yoon S-H, Kim M, Na H, et al. Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species. Int J Syst Evol Microbiol. 2012;62:716–21.

  21. 21.

    EZTaxon [http://eztaxon-e.ezbiocloud.net/]

  22. 22.

    Pagani I, Liolios K, Jansson J, Chen IM, Smirnova T, Nosrat B, et al. The Genomes OnLine Database (GOLD) v. 4: status of genomic and metagenomic projects and their associated metadata. Nucleic Acids Res. 2012;40:D571–9.

  23. 23.

    Markowitz VM, Chen I-MA, Palaniappan K, Chu K, Szeto E, Pillay M, et al. IMG 4 version of the integrated microbial genomes comparative analysis system. Nucleic Acids Res. 2014;42:D560–7.

  24. 24.

    Mavromatis K, Land ML, Brettin TS, Quest DJ, Copeland A, Clum A, et al. The fast changing landscape of sequencing technologies and their impact on microbial genome assemblies and annotation. PLoS One. 2012;7, e48837.

  25. 25.

    Reeve WG, Tiwari RP, Worsley PS, Dilworth MJ, Glenn AR, Howieson JG. Constructs for insertional mutagenesis, transcriptional signal localization and gene regulation studies in root nodule and other bacteria. Microbiology. 1999;145:1307–16.

  26. 26.

    Protocols and sample preparation information [http://jgi.doe.gov/collaborate-with-jgi/pmo-overview/protocols-sample-preparation-information/]

  27. 27.

    Bennett S. Solexa Ltd. Pharmacogenomics. 2004;5:433–8.

  28. 28.

    JGI:Joint Genome Institute [http://www.jgi.doe.gov]

  29. 29.

    Zerbino D, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.

  30. 30.

    Reads simulator wgsim [https://github.com/lh3/wgsim]

  31. 31.

    Gnerre S, MacCallum I, Przybylski D, Ribeiro FJ, Burton JN, Walker BJ, et al. High-quality draft assemblies of mammalian genomes from massively parallel sequence data. P Natl A Sci. 2011;108:1513–8.

  32. 32.

    Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.

  33. 33.

    Mavromatis K, Ivanova NN, Chen IM, Szeto E, Markowitz VM, Kyrpides NC. The DOE-JGI Standard operating procedure for the annotations of microbial genomes. Stand Genomic Sci. 2009;1:63–7.

  34. 34.

    Chen IM, Markowitz VM, Chu K, Anderson I, Mavromatis K, Kyrpides NC, et al. Improving microbial genome annotations in an integrated database context. PLoS One. 2013;8, e54859.

  35. 35.

    Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.

  36. 36.

    Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35:7188–96.

  37. 37.

    Infernal: inference of RNA alignments [http://infernal.janelia.org/]

  38. 38.

    Markowitz VM, Mavromatis K, Ivanova NN, Chen IM, Chu K, Kyrpides NC. IMG ER: a system for microbial genome annotation expert review and curation. Bioinformatics. 2009;25:2271–8.

  39. 39.

    da Silva K, De Meyer SE, Rouws LF, Farias EN, dos Santos MA, O’Hara G, et al. Bradyrhizobium ingae sp. nov., isolated from effective nodules of Inga laurina grown in Cerrado soil. Int J Syst Evol Microbiol. 2014;64:3395–401.

  40. 40.

    Islam MS, Kawasaki H, Muramatsu Y, Nakagawa Y, Seki T. Bradyrhizobium iriomotense sp. nov., isolated from a tumor-like root of the legume Entada koshunensis from Iriomote Island in Japan. Biosci Biotechnol Biochem. 2008;72:1416–29.

  41. 41.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.

  42. 42.

    Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, et al. Towards a richer description of our complete collection of genomes and metagenomes “Minimum Information about a Genome Sequence” (MIGS) specification. Nature Biotechnol. 2008;26:541–7.

  43. 43.

    Field D, Amaral-Zettler L, Cochrane G, Cole JR, Dawyndt P, Garrity GM, et al. The Genomic Standards Consortium. PLoS Biol. 2011;9, e1001088.

  44. 44.

    Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. P Natl A Sci USA. 1990;87:4576–9.

  45. 45.

    Garrity GM, Bell JA, Lilburn T. Phylum XIV. Proteobacteria phyl. nov. In: Garrity GM, Brenner DJ, Kreig NR, Staley JT, editors. Bergey's Manual of Systematic Bacteriology. Volume 2. Second ed. New York: Springer - Verlag; 2005. p. 1.

  46. 46.

    Validation of publication of new names and new combinations previously effectively published outside the IJSEM. List no. 106. Int J Syst Evol Microbiol. 2005;55:223538.

  47. 47.

    Garrity GM, Bell JA, Lilburn T. Class I. Alphaproteobacteria class. In: Garrity GM, Brenner DJ, Kreig NR, Staley JT, editors. Bergey's Manual of Systematic Bacteriology. Volume 2. Second ed. New York: Springer - Verlag; 2005. p. 1.

  48. 48.

    Garrity GM, Bell JA, Lilburn T. Class I. Alphaproteobacteria class. In: Garrity GM, Brenner DJ, Kreig NR, Staley JT, editors. Bergey's Manual of Systematic Bacteriology. Volume 2. Second ed. New York: Springer - Verlag; 2005. p. 1.

  49. 49.

    Kuykendall LD. Order VI. Rhizobiales ord. nov. In: Garrity GM, Brenner DJ, Kreig NR, Staley JT, editors. Bergey's Manual of Systematic Bacteriology. Secondth ed. New York: Springer - Verlag; 2005. p. 324.

  50. 50.

    Garrity GM, Bell JA, Lilburn T. Family VII. Bradyrhizobiaceae fam. nov. In Bergey’s Manual of Systematic Bacteriology. Volume 2. Second edition. Edited by Brenn DJ; New York: Springer - Verlag; 2005: 438

  51. 51.

    Jordan DC. Transfer of Rhizobium japonicum Buchanan 1980 to Bradyrhizobium gen. nov., a genus of slow-growing, root nodule bacteria from leguminous plants. Int J Syst Bacteriol. 1982;32:136–9.

  52. 52.

    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.

  53. 53.

    Guide to GO evidence codes [http://www.geneontology.org/GO.evidence.shtml]

  54. 54.

    GOLD ID for Bradyrhizobium sp.Tv2a.2 [https://gold.jgi-psf.org/projects?id=9814]

Download references

Acknowledgements

This work was performed under the auspices of the US Department of Energy’s Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Berkeley National Laboratory under contract No. DE-AC02-05CH11231. We thank Gordon Thompson (Murdoch University) for the preparation of SEM and TEM photos. We would also like to thank the Center of Nanotechnology at King Abdulaziz University for their support.

Author information

Correspondence to Wayne Reeve.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MP supplied the strain and background information for this project and the DNA to the JGI, TR performed all imaging, TR and WR drafted the paper, MNB and NAB provided financial support and all other authors were involved in sequencing the genome and/or editing the final paper. All authors read and approved the final manuscript.

Additional file

Additional file 1:

Associated MIGS record.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Root-nodule bacteria
  • Nitrogen fixation
  • Symbiosis
  • Alphaproteobacteria
  • GEBA-RNB