Skip to main content

The complete genome sequence of Ensifer meliloti strain CCMM B554 (FSM-MA), a highly effective nitrogen-fixing microsymbiont of Medicago truncatula Gaertn


Strain CCMM B554, also known as FSM-MA, is a soil dwelling and nodule forming, nitrogen-fixing bacterium isolated from the nodules of the legume Medicago arborea L. in the Maamora Forest, Morocco. The strain forms effective nitrogen fixing nodules on species of the Medicago, Melilotus and Trigonella genera and is exceptional because it is a highly effective symbiotic partner of the two most widely used accessions, A17 and R108, of the model legume Medicago truncatula Gaertn. Based on 16S rRNA gene sequence, multilocus sequence and average nucleotide identity analyses, FSM-MA is identified as a new Ensifer meliloti strain. The genome is 6,70 Mbp and is comprised of the chromosome (3,64 Mbp) harboring 3574 predicted genes and two megaplasmids, pSymA (1,42 Mbp) and pSymB (1,64 Mbp) with respectively 1481 and 1595 predicted genes. The average GC content of the genome is 61.93%. The FSM-MA genome structure is highly similar and co-linear to other E. meliloti strains in the chromosome and the pSymB megaplasmid while, in contrast, it shows high variability in the pSymA plasmid. The large number of strain-specific sequences in pSymA as well as strain-specific genes on pSymB involved in the biosynthesis of the lipopolysaccharide and capsular polysaccharide surface polysaccharides may encode novel symbiotic functions explaining the high symbiotic performance of FSM-MA.


To secure their nitrogen supply, legumes such as alfalfa, pea, (soy−/faba-)bean establish an endosymbiotic interaction with soil bacteria collectively called rhizobia that can reduce atmospheric nitrogen gas and produce reduced nitrogen molecules metabolizable by the plants. This symbiosis between legumes and rhizobia is of ecological and economic importance because of its contribution to the global nitrogen cycle, its impact on sustainable agriculture and its biotechnological potential to ensure nitrogen supply in agriculture [1].

The reduction of atmospheric nitrogen by rhizobia takes place in a specific niche, within the cells of de novo formed organs called nodules found usually on the roots and in some cases on the stem of the plants. Nodule development is initiated when flavonoids released by the plants induce the expression of the bacterial nodulation (nod) genes resulting in the production of the lipo-chitooligosaccharide signal molecules, the Nod factors. Nod factors cause a change in the direction of polar growth in developing root hairs and simultaneously induce cell division in the root cortex cells. As a result, a nodule primordium is formed that turns into meristematic tissue to produce the cells of the nodule and bacteria become entrapped in the curled root hair where they form an infection pocket. From the site of the infection pocket, a tubular structure, called infection thread, is formed in the root hair that grows toward the cells of the developing nodule. In the infection thread, bacteria multiply and finally they are released into the cytoplasm of the nodule cells via a mechanism resembling endocytosis resulting in organelle-like structures called symbiosomes. Symbiosomes have a membrane of plant origin which surrounds one or more bacteria. After bacterial release, the cells of both partners differentiate into mature symbiotic cells. The nodule cells become enlarged polyploid cells which host several tens of thousands of bacteria that are themselves differentiated into a nitrogen-fixing form called bacteroid [2,3,4]. Interestingly, in Medicago and closely related species like Pisum and Vicia , the host imposes a terminal differentiation on the bacterial partner that is accompanied by the increase in the DNA content and size of the bacteroids and results in the loss of their cell division capacity [5]. This terminal differentiation is orchestrated by nodule-specific cysteine-rich peptides that are expressed exclusively in the infected cells of the nodule [6, 7].

To effectively investigate these interactions, two genetic model legume species, Lotus japonicus (Regel) K. Larsen (bird’s-foot trefoil) and Medicago truncatula Gaertn. (barrel clover/barrel medic) have been chosen for which structural and functional genomics tools and databases have been developed [8, 9]. M. truncatula is a diploid, self-pollinating annual plant belonging to the Medicago genus, which contains species that are among the most extensively cultivated forage and pasture plants. Medicago plants establish symbiosis only with a limited number of bacterial species, mainly with Ensifer (synonym Sinorhizobium ) meliloti and Ensifer medicae , and with certain Ensifer fredii strains and Rhizobium mongolense [10,11,12]. However, some combinations of wild-type plants (species, sub-species and ecotypes) and bacterial strains of the most-studied bacterial species, E. meliloti and E. medicae , often lead to incompatible interactions [13,14,15,16,17], i.e. nodule formation is initiated but bacteria cannot invade nodules or cannot persist and fix nitrogen in the symbiotic organ. The incompatibility can be caused by functions/proteins encoded by genes in the accessory genome of the bacteria [14] such as the strain-specific HrrP peptidase [18], strain specific exopolysaccharide production [19] and/or allelic variants of the host genes like the NFS1 and NFS2 genes encoding NCR peptides in M. truncatula [20, 21]. Strikingly, the model bacterium E. meliloti strain 1021 (with the reference genome and most of the available mutants) is poorly matched for nitrogen fixation with the most widely used M. truncatula accessions Jemalong A17 and M. truncatula ssp. tricycla R108 [22, 23].

E. meliloti strain FSM-MA (first catalogued as E. arboris strain CCMM B554, also known as LMG-R33403 and MR372) was isolated from the nodules of Medicago arborea L. (moontrefoil/tree medic) in Maamora Forest between Rabat and Meknes, Morocco, and is stored in The Moroccan Coordinated Collections of Microorganisms as CCMM B554. Recently, Kazmierczak et al. [22] identified E. meliloti strain FSM-MA as a highly effective symbiotic partner of the two most widely used M. truncatula ecotypes, A17 and R108, as well as all tested Medicago sativa L. (alfalfa) cultivars. To gain the potential to identify novel bacterial symbiotic genes and genes associated with FSM-MA’s exceptional symbiotic performance, we sequenced the genome of the strain FSM-MA. Here we present a summary classification and a set of general features for E. meliloti strain FSM-MA, together with a description of its genome sequence and annotation.

Organism information

Classification and features

E. meliloti FSM-MA is a motile, non-sporulating, Gram-negative strain (Fig. 1) in the order Rhizobiales of the class Alphaproteobacteria . This fast growing strain forms colonies within 3 days on YEB agar plates [22] at 30 °C. The colonies (Fig. 1a, b) are light beige colored on YEB plates, slightly doomed, mucoid and have a smooth margin. The rod shaped free-living form (Fig. 1c, d) has dimensions of 1.0–2.0 μm in length and approximately 0.5 μm in width, while bacteroids in M. truncatula Jemalong A17 nodules (Fig. 1e, f) have the same width and are elongated to 5–8 μm. A summary of the classification is provided in Table 1.

Fig. 1

Colony morphology of E. meliloti strain FSM-MA on solid medium (a,b) at 5× (A) and 50× (b) magnifications as well as scanning electron microscopic images at 4000× (c,e) and 20,000× magnifications (d,f) of free-living cells (c,d) and bacteroids isolated from M. truncatula Jemalong A17 nodules (e,f)

Table 1 Classification and general features of E. meliloti strain FSM-MA

Extended feature descriptions

Phylogenetic analysis of E. meliloti strain FSM-MA was performed by aligning the 16S rRNA sequence to the 16S rRNA sequences (consensus sequence length of 1346 basepairs (bp)) of other Ensifer strains (Fig. 2). The FSM-MA 16S rRNA sequence has 100% sequence identity with those of the widely used E. meliloti strains such as strain 1021 or Rm41, while four mismatches can be observed with the E. medicae strain WSM419 sequence. Moreover, there are five mismatches between the 16S rRNA sequence fragments of strain FSM-MA and E. arboris strain LMG14919 T. A Multilocus Sequence Analysis (Additional file 1: Figure S1) using 14 chromosomal genes further confirmed FSM-MA as an E. meliloti strain and clearly separated it from E. arboris strain LMG14919 T, E. medicae strain WSM419 and the E. fredii strains NGR234, USDA257 and HH103. Among the E. meliloti strains, strain FSM-MA is most closely related to strains BO21CC and BL225C which were isolated from M. sativa nodules in Lodi, Italy [24]. Finally, the two-way average nucleotide identity [25] was calculated between genomes using the default settings of the ANI calculator. The genome of strain FSM-MA showed 99,42% identity with the genome of E. meliloti strain 1021 and only 90,09% identity with the genome of E. arboris strain LMG14919 T , 87,09% identity with the genome of E. medicae strain WSM419 and 83,16% and 83,31% identity with the genomes of E. fredii strains NGR234 and HH103, respectively. Once more this analysis showed that FSM-MA is an E. meliloti strain and not an E. arboris strain, considering a cut-off for species delineation at 95% identity [25].

Fig. 2

Phylogenetic tree showing the relationship of E. meliloti strain FSM-MA. The evolutionary history was inferred by using the Maximum Likelihood method based on the General Time Reversible model [30]. The tree with the highest log likelihood (−2208.71) is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach, and then selecting the topology with superior log likelihood value. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. Evolutionary analyses were conducted in MEGA7 [31]


Strain FSM-MA forms effective nitrogen fixing nodules on Medicago species M. sativa L., M. truncatula, M. arborea L., M. sativa subsp. x varia, M. ruthenica (L.) Trautv. as well as on Trigonella calliceras Fisch., Melilotus albus (L.) Lam. (white sweetclover) and Melilotus officinalis (L.) Lam. (yellow sweetclover). Moreover, in agreement with its classification as E. meliloti , it nodulates Medicago polymorpha L. (burclover) –that forms nitrogen-fixing symbiosis with E. medicae strains – but there is no nitrogen fixation in the formed nodules.

Genome sequencing information

Genome project history

This organism was selected for sequencing on the basis of its superior symbiotic performance [22] with the most widely used accessions (A17 and R108) of the model legume M. truncatula. The genome project and the sequence of the three replicons are deposited in the National Center for Biotechnology Information (NCBI; accession numbers: CP019584, CP019585, CP019586). Genome sequencing and sequence assembling were performed at the University of Malaya (Kuala Lumpur, Malaysia) and at the Seqomics Biotechnology Ltd. (Mórahalom, Hungary). Annotation was carried out at Seqomics Biotechnology Ltd. A summary of the project information can be found in Table 2.

Table 2 Genome sequencing project information for E. meliloti strain FSM-MA

Growth condition and genomic DNA preparation

E. meliloti strain FSM-MA was grown on solid YEB medium (0.5% beef extract; 0.1% yeast extract; 0.5%peptone; 0.5% sucrose; 0.04% MgSO4.7H2O; pH 7.5) for 3 days and a single colony was used to inoculate 3 ml YEB broth medium. The culture was grown for 24 h on a gyratory shaker at 225 rpm at 30 °C, then 0.5 ml of the starter culture was used to inoculate 50 ml YEB broth medium. The culture was grown at 30 °C at 225 rpm until OD600 = 0.6 was reached. DNA was isolated from the cells with the MasterPure Complete DNA and RNA Purification Kit (Epicentre). The integrity of the extracted genomic DNA was analyzed by 0.7% agarose gel electrophoresis. The final concentration of the DNA, estimated with the help of a Qubit Fluorometer (ThermoFisher Scientific), was 0.45 mg ml−1.

Genome sequencing and assembly

The genome sequence of E. meliloti strain FSM-MA was generated using Pacific BioScience (PacBio) and Illumina technologies. An Illumina Mate Paired library (average insert length 7 kbp) was constructed and sequenced using the Illumina MiSeq platform, which generated 3,387,162 reads. Similarly, a PacBio SMRTbell library was constructed and sequenced on the PacBio RS II platform to generate 254,443 filtered reads (N50 value at 8643 bp and total bases at 1,726,776,880 bp). Assembly was then carried out using HGAP version 3 [26] yielding three contigs with an average coverage of 186.71×. Subsequently, Illumina reads were aligned to the PacBio assembly with the help of the CLC Genomics Workbench version 9.5 and the observed 17 InDels were corrected. The final assembly contains three circular contigs corresponding to the three replicons (the chromosome and the pSymA and pSymB megaplasmids) totaling 6,703,999 bp and total input read coverage was at 249.2×.

Genome annotation

Genes were identified and annotated using the NCBI Prokaryotic Genome Annotation Pipeline. The NCBI non-redundant database, UniProt, TIGR/Fam, Pfam, PRIAM, KEGG, COG, and InterPro databases were used to analyse the predicted coding sequences after translation. HMMER [27] and tRNAscan-SE [28] were used to identify the rRNA and tRNA genes, respectively.

Genome properties

The genome is 6,703,999 bp and comprised of three replicons (Table 3) with the size of 3,641,423 bp (chromosome), 1,422,736 bp (pSymA) and 1,639,840 bp (pSymB). The average GC content is 61.93%. Three rRNA operons, 67 RNA only genes were identified and 6583 protein coding genes were predicted in the genome. Five thousand thirty-two protein-coding genes were assigned a putative function and 1551 genes were predicted to code for hypothetical proteins (Table 4). The distribution of genes in COG functional categories is presented in Table 5.

Table 3 Summary of genome: one chromosome and 2 plasmids
Table 4 Genome statistics for E. meliloti strain FSM-MA
Table 5 Number of genes of Ensifer meliloti FSM-MA associated with general COG functional categories

Insights from the genome sequence

The genome size of FSM-MA falls within the expected size range of 6.65–8.94 Mbp observed in the 33 sequenced E. meliloti genomes that have been deposited in the Integrated Microbial Genomes (IMG) database. The genome of all E. meliloti strains is composed of a circular chromosome and two megaplasmids/chromids, however, certain strains harbour additional replicons too. In strain FSM-MA, however, no additional plasmid was detected. The strain contains three rRNA gene clusters as other E. meliloti strains do. Similarly to other Ensifer strains, the highest number of genes is assigned to the COG functional categories amino acid transport and metabolism (9.46%), carbohydrate transport and metabolism (8.65%) and transcription (7.76%). An enrichment of the COG functional categories amino acid transport and metabolism, transcription and signal transduction mechanisms is observed in pSymA, while carbohydrate transport and metabolism and cell wall/membrane/envelope biogenesis are overrepresented on pSymB (Table 5).

Extended insights

Comparing the FSM-MA genome structure with that of other E. meliloti strains using the Mauve software [29] revealed high co-linearity of the chromosomes and the pSymB megaplasmids in contrast to the pSymA plasmids that are highly variable. For example, the average sequence identity between FSM-MA and strain 1021 is 99.4% and their chromosomes and pSymB plasmids are essentially co-linear. The major differences between the chromosomes originated from the insertion of three putative prophages/insertion elements into the FSM-MA genome at genes coding for tRNAs (SMB554__06910: tRNA-Thr, SMB554_09150: tRNA-Lys, SMB554_16265: tRNA-Met). These inserted elements are of approximately 48, 43 and 44 kbp and contains 70, 54 and 34 predicted ORFs, respectively. In the putative prophages at tRNA-Thr and tRNA-Lys, among hypothetical proteins, a number of phage related functions such as terminase, phage portal and capsid proteins (both prophages) as well as ORFs encoding endonucleases, transcriptional regulators, site-specific integrase, DNA ligase, peptidase or peptidoglycan-binding protein are encoded (prophage at tRNA-Lys). The inserted sequence at tRNA-Met seems to contain genes coding for type I restriction-modification system elements, an N6-DNA-methylase, chromosome segregation and AAA family ATPases as well as transcriptional regulators among hypothetical proteins. On the other hand, one putative prophage in the 1021 genome at a tRNA-Ser_CGA gene and the SMc01989-SMc02032 gene cluster coding for transcriptional regulators, membrane transporter and oxido-reductase elements are missing from the FSM-MA genome. The differences between the pSymB plasmids are mainly attributed to mobile genetic elements (IS elements, transposons) that are associated with strain-specific genes, essentially coding for proteins involved in the biosynthesis and transport of strain-specific LPS (lipopolysaccharide) and K-antigen (capsular polysaccharide) surface polysaccharides (discussed later). The pSymA plasmids – that are the carriers of major symbiotic functions such as genes encoding Nod factor biosynthesis and the nitrogenase enzyme and co-factor biosynthesis – have a number of co-linear blocks but have about 80 kbp size difference (FSM-MA > 1021), and more than 200 kbp (>1/7) of the sequences are absent in the other strain.

As the FSM-MA strain is interesting from the symbiotic point of view, we analysed those genes that are important for the development and functioning of the nitrogen-fixing symbioses. The initiation of the symbiotic interaction requires the production of Nod factors with proper chemical structure via the activity of the so-called Nod, Noe and Nol proteins. The FSM-MA genome contains all the known nod, noe and nol genes described in E. meliloti . The nif and fix genes code for the structural elements of the nitrogenase complex (nitrogenase, nitrogenase reductase, electron transport proteins) performing the reduction of atmospheric nitrogen as well as for proteins required for the biosynthesis of co-factors and the assembly of the complex. All these genes – including the ones that are present in multiple copies such as the three fixNOQP operons – can be found in the FSM-MA genome. Notably, despite the high diversity of the E. meliloti pSymA plasmids harbouring these symbiotic genes, the arrangement and the genomic environment of the nodulation and nitrogen fixation genes in FSM-MA and strain 1021 are the same.

Surface polyasaccharides play an essential role during the infection process [4] when bacteria enter the cells of the developing nodules via the infection threads. In the Medicago - E. meliloti symbiosis, the production of the succinoglycan exopolysaccharide is required for the continuous growth of the infection threads and its lack can be suppressed by the production of galactoglycan or certain capsular polysaccharides. Lipopolysaccharides might also affect both the infection and bacteroid differentiation processes. The organization and genomic environment of genes for the production and transport of the species-specific polysaccahrides EPS I (exo and exs genes) and EPS II (exp genes) as well as of the conserved part (lipidA and O-antigene core) of LPS (chromosomal and pSymB-born genes) and the KPS transporters are the same in the two E. meliloti strains. In contrast, the genes responsible for the production of the strain-specific polysaccharide moieties of LPS and KPS (Additional file 2: Figure S2), located on pSymB, are unique for the given strains.


The genome sequence of FSM-MA is of particular interest because the strain is highly effective with the most widely used ecotypes, Jemalong and R108 of the model legume M. truncatula. Comparative genomics with less and similarly effective strains as well as the creation and use of genomic libraries from FSM-MA has the potential to identify novel symbiotic genes and genes/operons that contribute to the exceptional symbiotic performance of the strain.







Capsular polysaccharide




Nodule-specific cysteine-rich


  1. 1.

    Oldroyd GE, Dixon R. Biotechnological solutions to the nitrogen problem. Curr Opin Biotechnol. 2014;26:19–24.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Oldroyd GE, Murray JD, Poole PS, Downie JA. The rules of engagement in the legume-rhizobial symbiosis. Annu Rev Genet. 2011;45:119–44.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Kereszt A, Mergaert P, Kondorosi E. Bacteroid development in legume nodules: evolution of mutual benefit or of sacrificial victims? Mol Plant-Microbe Interact. 2011;24:1300–9.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Gibson KE, Kobayashi H, Walker GC. Molecular determinants of a symbiotic chronic infection. Annu Rev Genet. 2008;42:413–41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Mergaert P, Uchiumi T, Alunni B, Evanno G, Cheron A, Catrice O, et al. Eukaryotic control on bacterial cell cycle and differentiation in the Rhizobium-legume symbiosis. Proc Natl Acad Sci U S A. 2006;103:5230–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Van de Velde W, Zehirov G, Szatmari A, Debreczeny M, Ishihara H, Kevei Z, et al. Plant peptides govern terminal differentiation of bacteria in symbiosis. Science. 2010;327:1122–6.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Guefrachi I, Nagymihaly M, Pislariu CI, Van de Velde W, Ratet P, Mars M, et al. Extreme specificity of NCR gene expression in Medicago truncatula. BMC Genomics. 2014;15:712.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Stacey G, Libault M, Brechenmacher L, Wan J, May GD. Genetics and functional genomics of legume nodulation. Curr Opin Plant Biol. 2006;9:110–21.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Suzaki T, Yoro E, Kawaguchi M. Leguminous plants: inventors of root nodules to accommodate symbiotic bacteria. Int Rev Cell Mol Biol. 2015;316:111–58.

    Article  PubMed  Google Scholar 

  10. 10.

    Kuykendall LD, Hashem FM, Bauchan GR, Devine TE, Dadson RB. Symbiotic competence of Sinorhizobium fredii on twenty alfalfa cultivars of diverse dormancy. Symbiosis. 1999;27:1–16.

    Google Scholar 

  11. 11.

    Rome S, Fernandez MP, Brunel MP, Normand P, Cleyet-Marel JC. Evidence that two genomic species of Rhizobium are associated with Medicago truncatula. Int J Syst Bacteriol. 1996;46:972–80.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Van Berkum P, Beyene D, Bao G, Campbell TA, Eardly BD. Rhizobium mongolense sp. nov. is one of three rhizobial genotypes identified which nodulate and form nitrogen-fixing symbioses with Medicago ruthenica (L.) Ledebour. Int J Syst Bacteriol. 1998;48:13–22.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Bena G, Lyet A, Huguet T, Olivieri I. Medicago-Sinorhizobium symbiotic specificity evolution and the geographic expansion of Medicago. J Evol Biol. 2005;18:1547–58.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Crook MB, Lindsay DP, Biggs MB, Bentley JS, Price JC, Clement SC, et al. Rhizobial plasmids that cause impaired symbiotic nitrogen fixation and enhanced host invasion. Mol Plant-Microbe Interact. 2012;25:1026–33.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Liu J, Yang S, Zheng Q, Zhu H. Identification of a dominant gene in Medicago truncatula that restricts nodulation by Sinorhizobium meliloti strain Rm41. BMC Plant Biol. 2014;14:167.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Tirichine L, de Billy F, Huguet T. Mtsym6, a gene conditioning Sinorhizobium strain-specific nitrogen fixation in Medicago truncatula. Plant Physiol. 2000;123:845–51.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Sugawara M, Epstein B, Badgley BD, Unno T, Xu L, Reese J, et al. Comparative genomics of the core and accessory genomes of 48 Sinorhizobium strains comprising five genospecies. Genome Biol. 2013;14:R17.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Price PA, Tanner HR, Dillon BA, Shabab M, Walker GC, Griffitts JS. Rhizobial peptidase HrrP cleaves host-encoded signaling peptides and mediates symbiotic compatibility. Proc Natl Acad Sci U S A. 2015;112:15244–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Simsek S, Ojanen-Reuhs T, Stephens SB, Reuhs BL. Strain-ecotype specificity in Sinorhizobium meliloti-Medicago truncatula symbiosis is correlated to succinoglycan oligosaccharide structure. J Bacteriol. 2007;189:7733–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Wang Q, Yang S, Liu J, Terecskei K, Abraham E, Gombar A, et al. Host-secreted antimicrobial peptide enforces symbiotic selectivity in Medicago truncatula. Proc Natl Acad Sci U S A. 2017;114:6854–9.

    CAS  PubMed  Google Scholar 

  21. 21.

    Yang S, Wang Q, Fedorova E, Liu J, Qin Q, Zheng Q, et al. Microsymbiont discrimination mediated by a host-secreted peptide in Medicago truncatula. Proc Natl Acad Sci U S A. 2017;114:6848–53.

    CAS  PubMed  Google Scholar 

  22. 22.

    Kazmierczak T, Nagymihaly M, Lamouche F, Barriere Q, Guefrachi I, Alunni B, et al. Specific host-responsive associations between Medicago truncatula accessions and Sinorhizobium strains. Mol Plant-Microbe Interact. 2017;30:399–409.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Terpolilli JJ, O'Hara GW, Tiwari RP, Dilworth MJ, Howieson JG. The model legume Medicago truncatula A17 is poorly matched for N2 fixation with the sequenced microsymbiont Sinorhizobium meliloti 1021. New Phytol. 2008;179:62–6.

    Article  PubMed  Google Scholar 

  24. 24.

    Carelli M, Gnocchi S, Fancelli S, Mengoni A, Paffetti D, Scotti C, et al. Genetic diversity and dynamics of Sinorhizobium meliloti populations nodulating different alfalfa cultivars in Italian soils. Appl Environ Microbiol. 2000;66:4785–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007;57:81–91.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10:563–9.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39:W29–37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Darling AC, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14:1394–403.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Nei M, Kumar S. Molecular evolution and phylogenetics. New York: Oxford University Press; 2000.

    Google Scholar 

  31. 31.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Garrity GM, Bell JA, Lilburn T. Phylum XIV. Proteobacteria phyl. Nov. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey's manual of systematic bacteriology. 2nd ed. New York: Springer; 2005.

    Google Scholar 

  34. 34.

    Garrity GM, Bell JA, Lilburn T. Class I. Alphaproteobacteria class. Nov. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey's manual of systematic bacteriology. 2nd ed. New York: Springer; 2005.

    Google Scholar 

  35. 35.

    Validation EJ, No L. 107. List of new names and new combinations previously effectively, but not validly, published. Int J Syst Evol Microbiol. 2006;56:1–6.

    Article  Google Scholar 

  36. 36.

    Kuykendall LD. Order VI. Rhizobiales ord. Nov. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey's manual of systematic bacteriology. 2nd ed. New York: Springer; 2005.

    Google Scholar 

  37. 37.

    Skerman VBD, McGowan V, Sneath PHA. Approved lists of bacterial names. Int J Syst Bacteriol. 1980;30:225–420.

    Article  Google Scholar 

  38. 38.

    Kuykendall LD. Family I. Rhizobiaceae. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey's manual of systematic bacteriology. 2nd ed. New York: Springer; 2005.

    Google Scholar 

  39. 39.

    Casida LE. Ensifer adhaerens gen. Nov., sp. nov.: a bacterial predator of bacteria in soil. Int J Syst Bacteriol. 1982;32:339–45.

    Article  Google Scholar 

  40. 40.

    de Lajudie P, Willems A, Pot B, Dewettinck D, Maestrojuan G, Neyra M, et al. Polyphasic taxonomy of rhizobia: emendation of the genus Sinorhizobium and description of Sinorhizobium meliloti comb. nov., Sinorhizobium saheli sp. nov., and Sinorhizobium teranga sp. nov. Int J Syst Bacteriol. 1994;44:715–33.

    Article  Google Scholar 

  41. 41.

    Kuykendall LD, Hashem FM, Wang ET. Genus VII. Sinorhizobium. In: Garrity GM, Brenner DJ, Krieg NR, Staley JT, editors. Bergey's manual of systematic bacteriology. New York: Springer; 2005.

    Google Scholar 

  42. 42.

    Judicial Commission of the International Committee on Systematics of Prokaryotes. The genus name Sinorhizobium Chen et al. 1988 is a later synonym of Ensifer Casida 1982 and is not conserved over the latter genus name, and the species name ‘Sinorhizobium adhaerens’ is not validly published. Opinion 84. Int J Syst Evol Microbiol. 2008;58:1973.

    Article  Google Scholar 

  43. 43.

    Chen WX, Yan GH, Li JL. Numerical taxonomic study of fast-growing soybean rhizobia and a proposal that Rhizobium fredii be assigned to Sinorhizobium gen. Nov. Int J Syst Bacteriol. 1988;28:392–7.

    Article  Google Scholar 

  44. 44.

    Geddes BA, Oresnik IJ. Physiology, genetics, and biochemistry of carbon metabolism in the alphaproteobacterium Sinorhizobium meliloti. Can J Microbiol. 2014;60:491–507.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Biological Agents: Technical rules for biological agents. [].

  46. 46.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors thank Anne Willems (Laboratory of Microbiology, University Gent, Belgium) for initial characterization of the strain CCMM B554/FSM-MA.


This work was supported by a grant from the National Research, Development and Innovation Office of Hungary to AK (grant number: 120,122), by a grant from the Agence National de la Recherche to PM (grant number ANR-13-BSV7–0013) and by a CNRS MEDICRABAT PRAD grant (PRAD N° 06–10) to PR and MA.

Author information




JI isolated the strain; MA performed the taxonomic characterization of the strain; PR performed the first characterization of the strain’s symbiotic performance and with JI and MA assembled and provided the strain and the background information. Symbiotaxonomic evaluation of the strain was performed by QB (Melilotus and Trigonella species) with the supervision of VG and PM and by MN (Medicago species) with the supervision of EK and AK. AF performed all the imaging. MN also isolated the genomic DNA for sequencing. PacBio and Illumina sequencing were performed by TMC supervised by KGC and by PB, respectively. PacBio reads were analysed and assembled by KWH under the supervision of KGC. BMV created the final assembly supervised by BB while BH carried out the annotation of the genome supervised by BB and AK. PM and AK devised the experiments, analysed the data and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Attila Kereszt.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Figure S1.

Multilocus Sequence Analysis of 14 genes, recA, gltA, glnA, ctrA, ftsA, ftsZ1, ftsZ2, gyrB, dnaK, pnp, rpoB, thrC, atpD and gap in E. meliloti strains FSM-MA, Sm1021, Su47, Rm41, AK58, AK83, SM11, GR4, BO21CC and BL225C, E. arboris strain LMG14919, E. medicae strain WSM419 and E. fredii strains USDA257, NGR234 and HH103. The concatenated gene sequences (total 23,220 bp) were aligned by ClustalW and a maximum likelihood tree was inferred from the aligned sequences using MEGA ver. 6.0.6 software (Tamura et al., 2007). The tree was estimated using the Tamura-Nei substitution model (Tamura and Nei, 1993). Bootstrap tests were performed with 1000 replications. The inset shows the topology of the maximum likelihood tree. Tamura K, Nei M. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10: 512–526. Tamura K, Dudley J, Nei M, Kumar S. (2007). MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599. (TIFF 14278 kb)

Additional file 2: Figure S2.

Comparison of the organization of genes responsible for the production of the strain-specifc KPS in E. meliloti strains FSM-MA, 1021 and Rm41. The gene clusters are located between conserved genes (red arrow) coding for a 3-methyl-2-oxobutanoate-hydroxymethyl transferase (MOBHMT) and a nucleotidyl transferase. Genes determining conserved functions in KPS production such as transport (RkpR, RkpS, RkpS) or chain-length determination (RkpZ) are drawn as solid blue boxes. Open arrows with blue line indicate strain-specific rkp genes. Mustard arrows indicate genes conserved between two strains in the region. Open arrows with black line show genes with unknown function or function that could not be related to KPS synthesis. The genes are not drawn to scale. HypProt: hypothetical protein; pAcetylT: putative acetyl transferase; pMethylT: putative methyl transferase; GlycosylT: glycosyl transferase; pLysozime: putative lysozyme; SecCaBProt: putative secreted calcium-binding protein; pMembProt: putative membrane protein. (TIFF 75 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nagymihály, M., Vásarhelyi, B.M., Barrière, Q. et al. The complete genome sequence of Ensifer meliloti strain CCMM B554 (FSM-MA), a highly effective nitrogen-fixing microsymbiont of Medicago truncatula Gaertn. Stand in Genomic Sci 12, 75 (2017).

Download citation


  • Ensifer meliloti
  • Root nodule bacteria
  • Nitrogen-fixation
  • Symbiosis