- Short genome report
- Open Access
High-quality permanent draft genome sequence of Rhizobium sullae strain WSM1592; a Hedysarum coronarium microsymbiont from Sassari, Italy
Standards in Genomic Sciences volume 10, Article number: 44 (2015)
Rhizobium sullae strain WSM1592 is an aerobic, Gram-negative, non-spore-forming rod that was isolated from an effective nitrogen (N2) fixing root nodule formed on the short-lived perennial legume Hedysarum coronarium (also known as Sulla coronaria or Sulla). WSM1592 was isolated from a nodule recovered from H. coronarium roots located in Ottava, bordering Sassari, Sardinia in 1995. WSM1592 is highly effective at fixing nitrogen with H. coronarium, and is currently the commercial Sulla inoculant strain in Australia. Here we describe the features of R. sullae strain WSM1592, together with genome sequence information and its annotation. The 7,530,820 bp high-quality permanent draft genome is arranged into 118 scaffolds of 118 contigs containing 7.453 protein-coding genes and 73 RNA-only encoding genes. This rhizobial genome is sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.
The accessibility and supply of nitrogen fertilizer is an ever-increasing challenge that world agriculture faces . Despite the fact the Earth’s atmosphere consists of approximately 78 % dinitrogen, it is in a form that must be converted before it can be utilised by plants . Conversion of N2 can be achieved by the chemical synthesis of natural gas but these methods can be considered unsustainable because of the use of exhaustible and costly fossil fuel resources . In addition, the manufacturing process not only increases the greenhouse gas emissions but also field N fertiliser application have been directly linked to contaminating and leading to detrimental effects in ecosystems and waterways. Alternatively, a more sustainable and environmentally friendly process of acquiring N is through the biological process of N fixation by diazotrophs . Most biological fixation in world agriculture is provided from the process of symbiotic nitrogen fixation, which occurs following the successful formation of an effective symbiosis by leguminous plants and bacterial microsymbionts .
The productive efficiencies of SNF in agricultural areas rely on considerable efforts by researchers and producers in matching suitable legume hosts with their compatible microsymbionts . Some agricultural areas farm with indigenous legumes, while others embark on introducing exotic legumes and their compatible microsymbionts from different geographical locations that are edaphically and climatically suited to their own . In Australia for instance, selection programs have enabled the domestication of new Mediterranean legume species and their microsymbionts . One such grazing legume species commercially introduced into Australian and New Zealand agriculture includes the Papilionoid legume Hedysarum coronarium (also known as Sulla coronaria or Sulla). Sulla is a deep-rooted, short-lived perennial pasture legume that is grown throughout Mediterranean countries where it is fed green, used for silage or as hay . It is noted that the microsymbionts of Sulla display a high level of specificity for nodulation and nitrogen fixation . However, when effectively nodulated Sulla plants have the ability to biologically fix large amounts of nitrogen for increased paddock fertility .
Rhizobium sullae strain WSM1592 is the current Australian commercial inoculant for Sulla after replacing strain CC1335 in 2006. This strain has also been deposited in the Western Australian Soil Microbiology collection and is available for research. WSM1592 was isolated in 1995 from a nodule collected from a Sulla plant sampled on a roadside in calcareous loamy sand near the Ottava agriculture research farm, east of Sassari in Sardinia, Italy. The location has a Mediterranean climate with a long-term mean seasonal rainfall of 547 mm. Here we present a preliminary description of the general features for Rhizobium sullae strain WSM1592 together with its genome sequence and annotation.
Classification and features
R. sullae strain WSM1592 is a motile, Gram-negative rod (Fig. 1 Left and Center) in the order Rhizobiales of the class Alphaproteobacteria . It is fast growing, forming colonies within 3–4 days when grown on half strength Lupin Agar (½LA) , tryptone-yeast extract agar (TY)  or a modified yeast-mannitol agar (YMA)  at 28 °C. Colonies on ½LA are white-opaque, slightly domed and moderately mucoid with smooth margins (Fig. 1 Right).
Figure 2 shows the phylogenetic relationship of R. sullae strain WSM1592 in a 16S rRNA gene sequence based tree. This strain is phylogenetically the most related to Rhizobium sullae IS 123T, Rhizobium leguminosarum USDA 2370T and Rhizobium phaseoli ATCC 14482T with sequence identities to the WSM1592 16S rRNA gene sequence of 100 %, 99.84 % and 99.84 %, respectively, as determined using the EzTaxon-e server . Rhizobium sullae IS 123T was isolated from a Hedysarum coronarium root nodule discovered in Southern Spain . In contrast, R. leguminosarum USDA 2370T was isolated from an effective nodule of Pisum sativum and is also able to nodulate Trifolium repens and Phaseolus vulgaris . R. phaseoli ATCC 14482T was originally isolated from nodules of Phaseolus vulgaris and has been shown to nodulate Trifolium repens , but not Pisum sativum .
Hedysarum coronarium is a short-lived perennial pasture legume native to the Mediterranean basin and throughout the Hedysarum genus there is a large degree of specificity in symbiotic compatibility within this region . Rhizobium sullae WSM1592 nodulates (Nod+) and fixes nitrogen effectively (Fix+) with Hedysarum coronarium . However, inoculation of H. spinosissimum, H. flexuosum and H. carnosum with WSM1592 results in mostly Nod- but always Fix-.
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 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
R. sullae WSM1592 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 CTAB (Cetyl trimethyl ammonium bromide) bacterial genomic DNA isolation method .
Genome sequencing and assembly
The draft genome of R. sullae strain WSM1592 was generated at the DOE Joint Genome Institute using state of the art technology . An Illumina Std shotgun library was constructed and sequenced using the Illumina HiSeq 2000 platform which generated 29,255,624 reads. All general aspects of library construction and sequencing performed at the JGI can be found at the JGI’s 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). Following steps were then performed for assembly: (1) filtered 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) Illumina reads were assembled with simulated read pairs using Allpaths–LG (version r39750) . Parameters for assembly steps were: 1) Velvet (velveth: 63 –shortPaired and velvetg: −very_clean yes –exportFiltered yes –min_contig_lgth 500 –scaffolding no–cov_cutoff 10) 2) wgsim (−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 118 contigs in 118 scaffolds. The total size of the genome is 7.5 Mbp and the final assembly is based on 2,498,075,850 bp of Illumina data, which provides an average of 877× coverage of the genome.
Genes were identified using Prodigal , as part of the DOE-JGI genome annotation pipeline [27, 28]. The predicted CDSs were translated and used to search the National Centre 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 7,530,820 nucleotides 59.87 % GC content (Table 3 and comprised of 118 scaffolds of 118 contigs. From a total of 7,526 genes, 7,453 were protein encoding and 73 RNA only encoding genes. The majority of genes (78.42 %) 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.
Rhizobium sullae WSM1592 was isolated from a root nodule of Hedysarum coronarium (also known as Sulla coronaria ). Phylogenetic analysis revealed that WSM1592 is the most closely related to Hedysarum coronarium IS 123T, which was also isolated from Hedysarum coronarium growing in Southern Spain. The genome of WSM1592 is the first to be described for a strain of Rhizobium sullae and is 7.5 Mbp, with a GC content of 59.87 %. As expected this genome contains the nitrogenase-RXN MetaCyc pathway characterized by the multiprotein nitrogenase complex and has been shown to fix effectively with Hedysarum coronarium . The genome attributes of WSM1592 will be important for the characterisation of the genetic determinants required for the establishment of an effective symbiosis with Hedysarum .
Genomic Encyclopedia of Bacteria and Archaea – Root Nodule Bacteria
Joint Genome Institute
Cetyl trimethyl ammonium bromide
Western Australian Soil Microbiology
Symbiotic Nitrogen Fixation
Herridge DF, Peoples MB, Boddey RM. Global inputs of biological nitrogen fixation in agricultural systems. Plant Soil. 2008;311:1–18.
Sprent JI. Legume nodulation: A global perspective. Oxford, U.K: Wiley-Blackwell; 2009.
Vance CP. Symbiotic nitrogen fixation and phosphorus acquisition. Plant nutrition in a world of declining renewable resources. Plant Physiol. 2001;127:390–7.
Howieson JG, Yates RJ, Foster K, Real D, Besier B. Prospects for the future use of legumes. In: Dilworth MJ, James EK, Sprent JI, Newton WE, editors. Leguminous nitrogen-fixing symbioses. London, UK: Elsevier; 2008. p. 363–94.
Sessitsch A, Howieson JG, Perret X, Antoun H, Martinez-Romero E. Advances in Rhizobium research. Crit Rev Plant Sci. 2002;21:323–78.
Howieson JG, Loi A, Carr SJ. Biserrula pelecinus L. - a legume pasture species with potential for acid, duplex soils which is nodulated by unique root-nodule bacteria. Aust J Agric Res. 1995;46:997–1009.
Foote A. Local cultivar adaptation of Mediterranean Sulla. N Z J Exp Agric. 1988;153:25–7.
Kishinevsky BD, Nandasena KG, Yates RJ, Nemas C, Howieson JG. Phenotypic and genetic diversity among rhizobia isolated from three Hedysarum species: H. spinosissimum, H. coronarium and H. flexuosum. Plant Soil. 2003;251:143–53.
Sulas L, Seddaiu G, Muresu R, Roggero PP. Nitrogen fixation of Sulla under mediterranean conditions. Agron J. 2009;101:1470–8.
Howieson JG, Ewing MA, D‘antuono MF. Selection for acid tolerance in Rhizobium meliloti. Plant Soil. 1988;105:179–88.
Beringer JE. R factor transfer in Rhizobium leguminosarum. J Gen Microbiol. 1974;84:188–98.
Terpolilli JJ. Why are the symbioses between some genotypes of Sinorhizobium and Medicago suboptimal for N2 fixation? PhD thesis. Murdoch University: Microbiology; 2009.
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.
Squartini A, Struffi P, Döring H, Selenska-Pobell S, Tola E, Giacomini A, et al. Rhizobium sullae sp. nov. (formerly ’Rhizobium hedysari’), the root-nodule microsymbiont of Hedysarum coronarium. Int J Syst Evol Microbiol. 2002;52:1267–76.
Ramirez-Bahena MH, Garcia-Fraile P, Peix A, Valverde A, Rivas R, Igual JM, et al. Revision of the taxonomic status of the species Rhizobium leguminosarum (Frank 1879) Frank 1889(AL), Rhizobium phaseoli Dangeard 1926(AL) and Rhizobium trifolii Dangeard 1926AL. R-trifolii is a later synonym of R-leguminosarum. Reclassification of the strain R-leguminosarum DSM 30132 (= NCIMB 11478) as Rhizobium pisi sp nov. Int J Syst Evol Microbiol. 2008;58:2484–90.
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.
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]
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.
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.
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.
Chen WX, Wang ET, Kuykendall LD. The Proteobacteria. 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.
Garrity GM, Bell JA, Lilburn T, Class I. Alphaproteobacteria. In: Garrity GM, Brenner DJ, Kreig NR, Staley JT, editors. Bergey‘s manual of systematic bacteriology. Secondth ed. New York: Springer - Verlag; 2005.
Kuykendall LD: Order VI. Rhizobiales. 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.
Kuykendall LD: Family I. Rhizobiaceae. In: Brenner DJ, Krieg NR, Staley JT, editors. Bergey’s manual of systematic bacteriology, vol. 2. 2nd ed. New York: Springer; 2005. p. 324–58. part C.
Kuykendall LD, Young JM, Martínez Romero E, Kerr A, Sawada H. Genus I. Rhizobium. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey’s Manual of Systematic Bacteriology, vol. 2. Second ed. New York: Springer - Verlag; 2005.
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.
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.
Reeve W, Ardley J, Tian R, Eshraghi 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. 2015;10:14.
Reddy TB, Thomas AD, Stamatis D, Bertsch J, Isbandi M, Jansson J, Mallajosyula J, Pagani I, Lobos EA, Kyrpides NC: The Genomes OnLine Database (GOLD) v.5: a metadata management system based on a four level (meta)genome project classification. Nucleic Acids Res 2014.
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.
CTAB DNA extraction protocol [http://jgi.doe.gov/collaborate-with-jgi/pmo-overview/protocols-sample-preparation-information/]
JGI Website [http://www.jgi.doe.gov]
Zerbino D, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.
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 U S A. 2011;108:1513–8.
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, SDM and WR drafted the paper, RY provided financial support and all other authors were involved in sequencing the genome and editing the final paper. All authors read and approved the final manuscript.