- Short genome report
- Open Access
High-quality permanent draft genome sequence of Bradyrhizobium sp. strain WSM1743 - an effective microsymbiont of an Indigofera sp. growing in Australia
Standards in Genomic Sciences volume 10, Article number: 87 (2015)
Bradyrhizobium sp. strain WSM1743 is an aerobic, motile, Gram-negative, non-spore-forming rod that can exist as a soil saprophyte or as a legume microsymbiont of an Indigofera sp. WSM1743 was isolated from a nodule recovered from the roots of an Indigofera sp. growing 20 km north of Carnarvon in Australia. It is slow growing, tolerates up to 1 % NaCl and is capable of growth at 37 °C. Here we describe the features of Bradyrhizobium sp. strain WSM1743, together with genome sequence information and its annotation. The 8,341,956 bp high-quality permanent draft genome is arranged into 163 scaffolds and 167 contigs, contains 7908 protein-coding genes and 75 RNA-only encoding genes and was sequenced as part of the Root Nodule Bacteria chapter of the Genomic Encyclopedia of Bacteria and Archaea project.
Rhizobia are soil-dwelling bacteria that have acquired the ability to establish associations with leguminous plants to symbiotically fix nitrogen. After infection of the plant, the rhizobia become established within root nodules and can fix atmospheric dinitrogen gas into ammonia using a reaction that is catalyzed by the nitrogenase enzyme . The export of fixed nitrogen to the plant improves growth and productivity under N-limiting environmental conditions. The effective use of the symbiosis leads to sustainable cropping systems with a net positive impact on the environment . In Australia, the majority of productive legumes and their rhizobia in agricultural systems have been deliberately, or accidentally, introduced since European settlement . However, recently, there has been an interest in the diversity of Australian native legumes and their microsymbionts .
The northwest of Western Australia is an ideal landscape to discover rhizobia nodulating indigenous legume flora  and is an area low in introduced legumes and inoculants. In 1996, an extensive survey was conducted of the area revealing a range of indigenous legume genera including a number of Indigofera spp. . In Australia, this species has been found at dispersed locations in the Northern Territory, Queensland and Western Australia on dark brown clay loams and frequently on lands under cultivation . The Australian Indigofera spp., based on their habitat, can be placed into three categories; i) shrubs, including I. brevidens, I. australis, I. adesmiifolia and some members of the I. pratensis group, which occur mainly on better soil types in the east coast, ii) perennial herbs, such as I. baileyi, I. efoliata, I. triflora, I. georgei, I. rugosa and members of the I. triflora and I. pratensis groups, which occur in the more arid, or seasonally dry, parts of Australia, iii) annual herbs, uncommon amongst the endemic species, including the two annual species, I. haplophylla and I. ammobia, which occur in the monsoon tropics and the Tanami and Great Sandy Deserts, respectively . The native species with wide extra-Australian distributions (particularly I. colutea, I. hirsuta, I. linnaei and I. linifolia) occur in a variety of habitats, mostly towards the northern parts of Australia. It is likely that these taxa now inhabit a greater range than they did before European settlement, and the Australian populations of these species may have been augmented by the introduction of seed from non-Australian sources .
Since there is a paucity of information regarding microsymbionts of Indigofera , a collection of root nodules was therefore obtained from the most prevalent Indigofera spp. present in northwest Australia and the microsymbionts from these nodules were then isolated. One microsymbiont, Bradyrhizobium sp. strain WSM1743, was isolated from a nodule recovered from the roots from an indigenous Indigofera sp. growing in red-brown sandy loam 40 m above sea level. The plant was located in natural bush land, approximately 20 km Northeast of the town Carnarvon in Western Australia . The collection area has a warm semi-arid climate with a long-term mean seasonal rainfall of 226 mm.
Strain WSM1743 was identified as a Bradyrhizobium sp. based on 16S rRNA typing . Most Bradyrhizobium spp., including WSM1743, cannot grow on sucrose or lactose, which may indicate the lack of a disaccharide uptake system . However, WSM1743, unlike other Bradyrhizobium spp., is able to grow at 37 °C and this ability could be a specific adaptation to the high soil temperatures experienced in the northwest of Western Australia . Here we present a summary classification and a set of general features for this microsymbiont together with a description of its genome sequence and annotation done as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria project  (Additional file 2).
Classification and features
Bradyrhizobium sp. strain WSM1743 is a motile, Gram-negative non-spore-forming rod (Fig. 1 Left and Center) in the order Rhizobiales of the class Alphaproteobacteria . It is slow growing, forming colonies within 7–10 days when grown on half strength Lupin Agar  at 28 °C. Colonies on ½LA are white-opaque, slightly domed and moderately mucoid with smooth margins (Fig. 1 Right). This strain was isolated together with 7 other bacteria from native Indigofera plants and physiologically characterised. Strain WSM1743 was identified as slow growing with poor growth in 1 % NaCl, no growth at 2 to 3 % NaCl and average growth in pH5 to 9 . Bradyrhizobium type strains have a slow generation time (9 to 18 h) and fail to grow in media containing 2 % NaCl , which indicates that WSM1743 belongs to this genus. The maximal growth temperature for most Bradyrhizobium strains is 33 to 35 °C, with many strains failing to grow above 34 °C . However, WSM1743 was able to grow at 37 °C on ½LA medium , and therefore extending the temperature range for Bradyrhizobium .
Figure 2 shows the phylogenetic relationship of Bradyrhizobium sp. strain WSM1743 in a 16S rRNA gene sequence based tree. This strain is phylogenetically the most related to the RNB type strains B. japonicum USDA 6T , B. lupini DSM30140T and B. yuanmingense LMG 21827T with sequence identities to the WSM1743 16S rRNA gene sequence of 99.78 %, 99.71 % and 99.63 %, respectively, as determined using the EzTaxon-e server . B. japonicum USDA6T was originally isolated in Japan from Glycine max root nodules and is able to nodulate and fix nitrogen effectively with several other Glycine species and Macroptillium atropurpureum . B. lupini DSM30140T is a microsymbiont of Lupinus luteus and Lupinus angustifolius . B. yuanmingense B071T was isolated from Lespedeza cuneata root nodules from China but is also able to nodulate and fix nitrogen effectively with Vigna unguiculata and Glycyrrhiza uralensis . Additionally, a recent report showed that B. yuanmingense and B. japonicum are the preferred microsymbionts of Vigna unguiculata and Vigna radiate in the subtropical region of China .
Bradyrhizobium sp. strain WSM1743 was isolated from Indigofera sp. nodules collected at site 32, north of Carnarvon, Western Australia . The site of collection contained several Australian native legumes, with a soil pH of 7.5. Symbiotic interactions of Bradyrhizobium sp. strain WSM1743 were assessed on three annual, one biennial and nine perennial exotic legume species that have agricultural use, or potential use, in southern Australia. WSM1743 consistently nodulated with the exotic legume species, Macroptilium atropurpureum and Phaseolus vulgaris , inconsistently with Ononis natrix and did not form nodules with Argyrolobium uniflorum , Chamaecytisus proliferus, Sutherlandia microphylla , Hedysarum coronarium , Medicago sativa , Ornithopus sativus , O. compressus, Trifolium burchellianum , T. polymorphum and T. uniflorum. Strain WSM1743 was able to consistently nodulate the Australian native legumes, Acacia saligna , Kennedia prorepens and K. coccinea, but could not nodulate Swainsona pterostylis , S. formosa and S. macculochiana . Additionally it was noted that the isolate could not nodulate Indigofera brevidens , an indigenous Indigofera found in the same location as the host of WSM1743 .
Genome sequencing 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, The Root Nodulating Bacteria chapter (GEBA-RNB) project at the U.S. Department of Energy, Joint Genome Institute . The genome project is deposited in the Genomes OnLine Database  and the high-quality permanent draft genome sequence in IMG . Sequencing, finishing and annotation were performed by the JGI using state of the art sequencing technology . A summary of the project information is shown in Table 2.
Growth conditions and genomic DNA preparation
Bradyrhizobium sp. strain WSM1743 was cultured to mid logarithmic phase in 60 ml of TY rich media on a gyratory shaker at 28 °C . DNA was isolated from the cells using a Cetyl trimethyl ammonium bromide bacterial genomic DNA isolation method .
Genome sequencing and assembly
The draft genome of Bradyrhizobium sp. strain WSM1743 was generated at the DOE Joint Genome Institute . An Illumina standard shotgun library was constructed and sequenced using the Illumina HiSeq 2000 platform, which generated 14,683,452 reads totaling 2.2 Gbp. All general aspects of library construction and sequencing performed at the JGI can be found at the JGI web site . 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). Artifact filtered sequence data was then screened and trimmed according to the k–mers present in the dataset (Mingkun L, Copeland A, Han J. unpublished). High–depth k–mers, presumably derived from MDA amplification bias, cause problems in the assembly, especially if the k–mer depth varies in orders of magnitude for different regions of the genome. Reads with high k–mer coverage (>30x average k–mer depth) were normalized to an average depth of 30x. Reads with an average kmer depth of less than 2x were removed. Following steps were then performed for assembly: (1) normalized Illumina reads were assembled using Velvet version 1.1.04  (2) 1–3 Kbp simulated paired end reads were created from Velvet contigs using wgsim  (3) normalized Illumina reads were assembled with simulated read pairs using Allpaths–LG (version r39750) . Parameters for assembly steps were: 1) Velvet optimizing parameters (−−v --s 51 --e 71 --i 2 --t 1 --f "-shortPaired -fastq $FASTQ" --o "-ins_length 250 -min_contig_lgth 500) 2) wgsim version 0.3.0 (−e 0–1 76–2 76 -r 0 -R 0 -X 0) 3) Allpaths–LG (PrepareAllpathsInputs: PHRED 64 = 1 PLOIDY = 1 FRAG COVERAGE = 125 JUMP COVERAGE = 25 LONG JUMP COV = 50, RunAllpathsLG: THREADS = 8 RUN = std_shredpairs TARGETS = standard VAPI_WARN_ONLY = True OVERWRITE = True). The final draft assembly contained 167 contigs in 163 scaffolds. The total size of the genome is 8.3 Mbp and the final assembly is based on 1238 Mbp of Illumina data, which provides an average of 440x coverage.
Genes were identified using Prodigal , as part of the DOE-JGI genome annotation pipeline [36, 37]. 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  was used to find tRNA genes, whereas ribosomal RNA genes were found by searches against models of the ribosomal RNA genes built from SILVA . 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 . Additional gene prediction analysis and manual functional annotation was performed within the Integrated Microbial Genomes-Expert Review system  developed by the Joint Genome Institute, Walnut Creek, CA, USA.
The genome is 8,341,956 nucleotides with 63.37 % GC content (Table 3) and comprised of 163 scaffolds of 167 contigs. From a total of 7983 genes, 7908 were protein encoding and 75 RNA only encoding genes. The majority of genes (71.51 %) were assigned a putative function whilst the remaining genes were annotated as hypothetical. The distribution of genes into COG functional categories is presented in Table 4.
Bradyrhizobium sp. WSM1743 belongs to a group of Alpha-rhizobia microsymbionts from native Australian legumes and was isolated from a nodule of an Indigofera species growing 20 km north of Carnarvon in northwestern Australia. Phylogenetic analysis revealed that WSM1743 is most closely related to B. japonicum USDA 6T , which was obtained from Glycine max root nodules from Japan and is able to nodulate and fix nitrogen effectively with several other Glycine species and Macroptillium atropurpureum . Strain WSM1743 has been shown to nodulate with Macroptilium atropurpureum and endemic Australian legumes including Acacia saligna , Kennedia prorepens and K. coccinea .
A comparison of the genome of WSM1743 to that of USDA 6T  reveals that WSM1743 has a lower GC content, gene count, coding base count %, rRNAcount, COG % and transmembrane %. In contrast, the paralogs % is much higher for WSM1743 than for USDA 6T (81.75 % versus 42.54 %). These two genomes are included within a group of 54 Bradyrhizobium genomes that have been deposited into the IMG database . Within this group, strains known to symbiotically fix nitrogen all contain the nitrogenase-RXN MetaCyc pathway that is characterized by the multiprotein nitrogenase complex. The genome of Bradyrhizobium sp. WSM1743, in conjunction with the other Bradyrhizobium genomes, will be important for on-going comparative and functional analyses of the plant microbe interactions required for the successful establishment of native Australian legume symbioses.
Genomic encyclopedia of bacteria and archaea – root nodule bacteria
Joint genome institute
Cetyl trimethyl ammonium bromide
Western Australian soil microbiology
Integrated microbial genomes-expert review
National centre for biotechnology information
- N2 :
Terpolilli JJ, Hood GA, Poole PS. What determines the efficiency of N2-fixing Rhizobium-Legume symbioses. Adv Microb Physiol. 2012;60:325–89.
Bohlool BB, Ladha JK, Garrity DP, George T. Biological nitrogen fixation for sustainable agriculture - A perspective. Plant Soil. 1992;141:1–11.
Donald CM. Temperate pasture species. Canberra: Australian National University Press; 1970.
Yates RJ, Howieson JG, Nandasena KG, O'Hara GW. Root-nodule bacteria from indigenous legumes in the north-west of Western Australia and their interaction with exotic legumes. Soil Biol Biochem. 2004;36:1319–29.
Wilson PG, Rowe R. A revision of the Indigofereae (Fabaceae) in Australia. 2. Indigofera species with trifoliolate and alternately pinnate leaves. Telopea. 2008;12:293–307.
Wilson PG, Rowe R. A revision of the Indigofereae (Fabaceae) in Australia. 1. Indigastrum and the simple or unifoliolate species of Indigofera. Telopea. 2004;10:651–82.
Glenn AR, Reeve WG, Tiwari RP, Dilworth MJ. Acid tolerance in root nodule bacteria. Novartis Found Symp. 1999;221:112–26. discussion 126–130.
Reeve W, Ardley J, Tian R, Eshragi L, Yoon J, Ngamwisetkun P, et al. A genomic encyclopedia of the root nodule bacteria: Assessing genetic diversity through a systematic biogeographic survey. Stand Genomic Sci In press.
Howieson JG, Ewing MA, D'antuono MF. Selection for acid tolerance in Rhizobium meliloti. Plant Soil. 1988;105:179–88.
Kuykendall LD. Genus I. Bradyrhizobium. In: Brenner DJ, Krieg NR, Staley JT, editors. Bergey's Manual of Systematic Bacteriology, vol. 2 part C. 2nd ed. New York: Springer; 2005. p. 438–43.
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.
Jordan DC. NOTES: 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.
Eckhardt MM, Baldwin IL, Fred EB. Studies of the Root-Nodule Organism of Lupinus. J Bacteriol. 1931;21:273–85.
Yao ZY, Kan FL, Wang ET, Wei GH, Chen WX. Characterization of rhizobia that nodulate legume species of the genus Lespedeza and description of Bradyrhizobium yuanmingense sp nov. Int J Syst Evol Microbiol. 2002;52:2219–30.
Zhang YF, Wang ET, Tian CF, Wang FQ, Han LL, Chen WF, et al. Bradyrhizobium elkanii, Bradyrhizobium yuanmingense and Bradyrhizobium japonicum are the main rhizobia associated with Vigna unguiculata and Vigna radiata in the subtropical region of China. FEMS Microbiol Lett. 2008;285:146–54.
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. Nat Biotechnol. 2008;26:541–7.
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.
Reddy TB, Thomas AD, Stamatis D, Bertsch J, Isbandi M, Jansson J, et al. The Genomes OnLine Database (GOLD) v. 5: a metadata management system based on a four level (meta)genome project classification. Nucleic Acids Res. 2015;43(Database issue):D1099-106. doi: 10.1093/nar/gku950.
Field D, Amaral-Zettler L, Cochrane G, Cole JR, Dawyndt P, Garrity GM, et al. The Genomic Standards Consortium. PLoS Biol. 2011;9:e1001088.
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.
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. Secondth ed. New York: Springer - Verlag; 2005.
Validation of publication of new names and new combinations previously effectively published outside the IJSEM. Int J Syst Evol Microbiol 2005, 55:2235–2238
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.
Garrity GM, Bell JA, Liburn T. Family VII. Bradyrhizobiaceae. In: Brenner DJ, Krieg NR, Staley JT, editors. Bergey's Manual of Systematic Bacteriology. Volume 2 part C. New York: Springer; 2005. p. 438–75.
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.
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.
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.
GOLD ID Bradyrhizobium sp. WSM1743 [https://gold.jgi-psf.org/project?id=9884]. Access date 23/10/2015.
Beringer JE. R factor transfer in Rhizobium leguminosarum. J Gen Microbiol. 1974;84:188–98.
CTAB DNA extraction protocol [http://jgi.doe.gov/collaborate-with-jgi/pmo-overview/protocols-sample-preparation-information/]. Access date 23/10/2015.
JGI Website [http://www.jgi.doe.gov]. Access date 23/10/2015.
Zerbino D, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.
wgsim [https://github.com/lh3/wgsim]. Access date 23/10/2015.
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. Proc Natl Acad Sci. 2011;108:1513–8.
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.
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.
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.
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.
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.
INFERNAL. Inference of RNA alignments [http://infernal.janelia.org]. Access date 23/10/2015.
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.
Kaneko T, Maita H, Hirakawa H, Uchiike N, Minamisawa K, Watanabe A, et al. Complete genome sequence of the soybean symbiont Bradyrhizobium japonicum strain USDA6T. Genes. 2011;2:763–87.
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, Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344, and Los Alamos National Laboratory under contract No. DE-AC02-06NA25396.
The authors declare that they have no competing interests.
RY and JH supplied the strain and background information for this project, TR supplied DNA to JGI and performed all imaging, LE and SDM drafted the paper, RT provided financial support and all other authors were involved in sequencing the genome and editing the final manuscript. All authors read and approved the final manuscript.