Cell Cycle Control by the Master Regulator CtrA in Sinorhizobium meliloti.

In all domains of life, proper regulation of the cell cycle is critical to coordinate genome replication, segregation and cell division. In some groups of bacteria, e.g. Alphaproteobacteria, tight regulation of the cell cycle is also necessary for the morphological and functional differentiation of cells. Sinorhizobium meliloti is an alphaproteobacterium that forms an economically and ecologically important nitrogen-fixing symbiosis with specific legume hosts. During this symbiosis S. meliloti undergoes an elaborate cellular differentiation within host root cells. The differentiation of S. meliloti results in massive amplification of the genome, cell branching and/or elongation, and loss of reproductive capacity. In Caulobacter crescentus, cellular differentiation is tightly linked to the cell cycle via the activity of the master regulator CtrA, and recent research in S. meliloti suggests that CtrA might also be key to cellular differentiation during symbiosis. However, the regulatory circuit driving cell cycle progression in S. meliloti is not well characterized in both the free-living and symbiotic state. Here, we investigated the regulation and function of CtrA in S. meliloti. We demonstrated that depletion of CtrA cause cell elongation, branching and genome amplification, similar to that observed in nitrogen-fixing bacteroids. We also showed that the cell cycle regulated proteolytic degradation of CtrA is essential in S. meliloti, suggesting a possible mechanism of CtrA depletion in differentiated bacteroids. Using a combination of ChIP-Seq and gene expression microarray analysis we found that although S. meliloti CtrA regulates similar processes as C. crescentus CtrA, it does so through different target genes. For example, our data suggest that CtrA does not control the expression of the Fts complex to control the timing of cell division during the cell cycle, but instead it negatively regulates the septum-inhibiting Min system. Our findings provide valuable insight into how highly conserved genetic networks can evolve, possibly to fit the diverse lifestyles of different bacteria.

Evolution of Intra-specific Regulatory Networks in a Multipartite Bacterial Genome.

Reconstruction of the regulatory network is an important step in understanding how organisms control the expression of gene products and therefore phenotypes. Recent studies have pointed out the importance of regulatory network plasticity in bacterial adaptation and evolution. The evolution of such networks within and outside the species boundary is however still obscure. Sinorhizobium meliloti is an ideal species for such study, having three large replicons, many genomes available and a significant knowledge of its transcription factors (TF). Each replicon has a specific functional and evolutionary mark; which might also emerge from the analysis of their regulatory signatures. Here we have studied the plasticity of the regulatory network within and outside the S. meliloti species, looking for the presence of 41 TFs binding motifs in 51 strains and 5 related rhizobial species. We have detected a preference of several TFs for one of the three replicons, and the function of regulated genes was found to be in accordance with the overall replicon functional signature: house-keeping functions for the chromosome, metabolism for the chromid, symbiosis for the megaplasmid. This therefore suggests a replicon-specific wiring of the regulatory network in the S. meliloti species. At the same time a significant part of the predicted regulatory network is shared between the chromosome and the chromid, thus adding an additional layer by which the chromid integrates itself in the core genome. Furthermore, the regulatory network distance was found to be correlated with both promoter regions and accessory genome evolution inside the species, indicating that both pangenome compartments are involved in the regulatory network evolution. We also observed that genes which are not included in the species regulatory network are more likely to belong to the accessory genome, indicating that regulatory interactions should also be considered to predict gene conservation in bacterial pangenomes.

Exploring the dynamics of bacterial community composition in soil: the pan-bacteriome approach.

We performed a longitudinal study (repeated observations of the same sample over time) to investigate both the composition and structure of temporal changes of bacterial community composition in soil mesocosms, subjected to three different treatments (water and 5 or 25 mg kg(-1) of dried soil Cd(2+)). By analogy with the pan genome concept, we identified a core bacteriome and an accessory bacteriome. Resident taxa were assigned to the core bacteriome, while occasional taxa were assigned to the accessory bacteriome. Core and accessory bacteriome represented roughly 35 and 50 % of the taxa detected, respectively, and were characterized by different taxonomic signatures from phylum to genus level while 15 % of the taxa were found to be unique to a particular sample. In particular, the core bacteriome was characterized by higher abundance of members of Planctomycetes, Actinobacteria, Verrucomicrobia and Acidobacteria, while the accessory bacteriome included more members of Firmicutes, Clamydiae and Proteobacteria, suggesting potentially different responses to environmental changes of members from these phyla. We conclude that the pan-bacteriome model may be a useful approach to gain insight for modeling bacterial community structure and inferring different abilities of bacteria taxa.

DuctApe: a suite for the analysis and correlation of genomic and OmniLog™ Phenotype Microarray data.

Addressing the functionality of genomes is one of the most important and challenging tasks of today's biology. In particular the ability to link genotypes to corresponding phenotypes is of interest in the reconstruction and biotechnological manipulation of metabolic pathways. Over the last years, the OmniLog™ Phenotype Microarray (PM) technology has been used to address many specific issues related to the metabolic functionality of microorganisms. However, computational tools that could directly link PM data with the gene(s) of interest followed by the extraction of information on gene-phenotype correlation are still missing. Here we present DuctApe, a suite that allows the analysis of both genomic sequences and PM data, to find metabolic differences among PM experiments and to correlate them with KEGG pathways and gene presence/absence patterns. As example, an application of the program to four bacterial datasets is presented. The source code and tutorials are available at http://combogenomics.github.io/DuctApe/.

The DivJ, CbrA and PleC system controls DivK phosphorylation and symbiosis in Sinorhizobium meliloti.

Sinorhizobium meliloti is a soil bacterium that invades the root nodules it induces on Medicago sativa, whereupon it undergoes an alteration of its cell cycle and differentiates into nitrogen-fixing, elongated and polyploid bacteroid with higher membrane permeability. In Caulobacter crescentus, a related alphaproteobacterium, the principal cell cycle regulator, CtrA, is inhibited by the phosphorylated response regulator DivK. The phosphorylation of DivK depends on the histidine kinase DivJ, while PleC is the principal phosphatase for DivK. Despite the importance of the DivJ in C. crescentus, the mechanistic role of this kinase has never been elucidated in other Alphaproteobacteria. We show here that the histidine kinases DivJ together with CbrA and PleC participate in a complex phosphorylation system of the essential response regulator DivK in S. meliloti. In particular, DivJ and CbrA are involved in DivK phosphorylation and in turn CtrA inactivation, thereby controlling correct cell cycle progression and the integrity of the cell envelope. In contrast, the essential PleC presumably acts as a phosphatase of DivK. Interestingly, we found that a DivJ mutant is able to elicit nodules and enter plant cells, but fails to establish an effective symbiosis suggesting that proper envelope and/or low CtrA levels are required for symbiosis.

Exploring the plant-associated bacterial communities in Medicago sativa L.

BACKGROUND: Plant-associated bacterial communities caught the attention of several investigators which study the relationships between plants and soil and the potential application of selected bacterial species in crop improvement and protection. Medicago sativa L. is a legume crop of high economic importance as forage in temperate areas and one of the most popular model plants for investigations on the symbiosis with nitrogen fixing rhizobia (mainly belonging to the alphaproteobacterial species Sinorhizobium meliloti). However, despite its importance, no studies have been carried out looking at the total bacterial community associated with the plant. In this work we explored for the first time the total bacterial community associated with M. sativa plants grown in mesocosms conditions, looking at a wide taxonomic spectrum, from the class to the single species (S. meliloti) level. RESULTS: Results, obtained by using Terminal-Restriction Fragment Length Polymorphism (T-RFLP) analysis, quantitative PCR and sequencing of 16 S rRNA gene libraries, showed a high taxonomic diversity as well as a dominance by members of the class Alphaproteobacteria in plant tissues. Within Alphaproteobacteria the families Sphingomonadaceae and Methylobacteriaceae were abundant inside plant tissues, while soil Alphaproteobacteria were represented by the families of Hyphomicrobiaceae, Methylocystaceae, Bradyirhizobiaceae and Caulobacteraceae. At the single species level, we were able to detect the presence of S. meliloti populations in aerial tissues, nodules and soil. An analysis of population diversity on nodules and soil showed a relatively low sharing of haplotypes (30-40%) between the two environments and between replicate mesocosms, suggesting drift as main force shaping S. meliloti population at least in this system. CONCLUSIONS: In this work we shed some light on the bacterial communities associated with M. sativa plants, showing that Alphaproteobacteria may constitute an important part of biodiversity in this system, which includes also the well known symbiont S. meliloti. Interestingly, this last species was also found in plant aerial part, by applying cultivation-independent protocols, and a genetic diversity analysis suggested that population structure could be strongly influenced by random drift.

Effect of Non-Lethal Selection on Spontaneous Revertants of Frameshift Mutations: The Escherichia coli hisF Case

Microorganisms possess the potential to adapt to fluctuations in environmental parameters, and their evolution is driven by the continuous generation of mutations. The reversion of auxotrophic mutations has been widely studied; however, little is known about the reversion of frameshift mutations resulting in amino acid auxotrophy and on the structure and functioning of the protein encoded by the revertant mutated gene. The aims of this work were to analyze the appearance of reverse mutations over time and under different selective pressures and to investigate revertant enzymes' three-dimensional structures and their correlation with a different growth ability. Escherichia coli FB182 strain, carrying the hisF892 single nucleotide deletion resulting in histidine auxotrophy, was subjected to different selective pressures, and revertant mutants were isolated and characterized. The obtained results allowed us to identify different indels of different lengths located in different positions in the hisF gene, and relations with the incubation time and the selective pressure applied were observed. Moreover, the structure of the different mutant proteins was consistent with the respective revertant ability to grow in absence of histidine, highlighting a correlation between the mutations and the catalytic activity of the mutated HisF enzyme.

Exploring the symbiotic pangenome of the nitrogen-fixing bacterium Sinorhizobium meliloti.

BACKGROUND: Sinorhizobium meliloti is a model system for the studies of symbiotic nitrogen fixation. An extensive polymorphism at the genetic and phenotypic level is present in natural populations of this species, especially in relation with symbiotic promotion of plant growth. AK83 and BL225C are two nodule-isolated strains with diverse symbiotic phenotypes; BL225C is more efficient in promoting growth of the Medicago sativa plants than strain AK83. In order to investigate the genetic determinants of the phenotypic diversification of S. meliloti strains AK83 and BL225C, we sequenced the complete genomes for these two strains. RESULTS: With sizes of 7.14 Mbp and 6.97 Mbp, respectively, the genomes of AK83 and BL225C are larger than the laboratory strain Rm1021. The core genome of Rm1021, AK83, BL225C strains included 5124 orthologous groups, while the accessory genome was composed by 2700 orthologous groups. While Rm1021 and BL225C have only three replicons (Chromosome, pSymA and pSymB), AK83 has also two plasmids, 260 and 70 Kbp long. We found 65 interesting orthologous groups of genes that were present only in the accessory genome, consequently responsible for phenotypic diversity and putatively involved in plant-bacterium interaction. Notably, the symbiosis inefficient AK83 lacked several genes required for microaerophilic growth inside nodules, while several genes for accessory functions related to competition, plant invasion and bacteroid tropism were identified only in AK83 and BL225C strains. Presence and extent of polymorphism in regulons of transcription factors involved in symbiotic interaction were also analyzed. Our results indicate that regulons are flexible, with a large number of accessory genes, suggesting that regulons polymorphism could also be a key determinant in the variability of symbiotic performances among the analyzed strains. CONCLUSIONS: In conclusions, the extended comparative genomics approach revealed a variable subset of genes and regulons that may contribute to the symbiotic diversity.

Plasmid electroporation of Sinorhizobium strains: The role of the restriction gene hsdR in type strain Rm1021.

Although horizontal gene transfer mediated by plasmids is important to the generation of the genetic variability of Sinorhizobium strains, the barriers which can reduce horizontal gene transfer between bacteria have not yet been studied in Sinorhizobium. We studied the plasmid transfer by electroporation and its restriction in strains of Sinorhizobium meliloti and S. medicae. After conditions for electroporation were established, three S. meliloti strains (including the sequenced type strain Rm1021) and two S. medicae strains were electroporated with plasmid DNA extracted from strains of both species. The efficiency of transformation was found to be variable among different strains. The acquisition of plasmid DNA was found to be donor-dependent in S. meliloti strain Rm1021 that prefers self-DNA more than the DNA from other Sinorhizobium strains. All other strains tested did not show a preference for self-DNA. In strain Rm1021, the inactivation of the hsdR gene, coding for a putative type-I restriction enzyme, increased the efficiency of transformation and conjugation with non-self DNA; the transformation capability was again reduced in hsdR mutant when the cloned hsdR gene was expressed from a lac promoter. Phylogenetic analysis of the hsdR gene clearly indicated that this gene was horizontally transferred to strain Rm1021, explaining its absence in the other strains tested.

Metabolic capacity of Sinorhizobium (Ensifer) meliloti strains as determined by phenotype MicroArray analysis.

Sinorhizobium meliloti is a soil bacterium that fixes atmospheric nitrogen in plant roots. The high genetic diversity of its natural populations has been the subject of extensive analysis. Recent genomic studies of several isolates revealed a high content of variable genes, suggesting a correspondingly large phenotypic differentiation among strains of S. meliloti. Here, using the Phenotype MicroArray (PM) system, hundreds of different growth conditions were tested in order to compare the metabolic capabilities of the laboratory reference strain Rm1021 with those of four natural S. meliloti isolates previously analyzed by comparative genomic hybridization (CGH). The results of PM analysis showed that most phenotypic differences involved carbon source utilization and tolerance to osmolytes and pH, while fewer differences were scored for nitrogen, phosphorus, and sulfur source utilization. Only the variability of the tested strain in tolerance to sodium nitrite and ammonium sulfate of pH 8 was hypothesized to be associated with the genetic polymorphisms detected by CGH analysis. Colony and cell morphologies and the ability to nodulate Medicago truncatula plants were also compared, revealing further phenotypic diversity. Overall, our results suggest that the study of functional (phenotypic) variability of S. meliloti populations is an important and complementary step in the investigation of genetic polymorphism of rhizobia and may help to elucidate rhizobial evolutionary dynamics, including adaptation to diverse environments.

Plant-by-plant variations of bacterial communities associated with leaves of the nickel hyperaccumulator Alyssum bertolonii Desv.

Bacteria associated with tissues of metal-hyperaccumulating plants are of great interest due to the multiple roles they may play with respect to plant growth and resistance to heavy metals. The variability of bacterial communities associated with plant tissues of three populations of Alyssum bertolonii, a Ni hyperaccumulator endemic of serpentine outcrops of Central Italy, was investigated. Terminal-restriction fragment length polymorphism (T-RFLP) analysis of bacterial 16S rRNA genes was applied to DNA extracted from leaf tissues of 30 individual plants from three geographically separated serpentine outcrops. Moreover, T-RFLP fingerprinting was also performed on DNA extracted from the same soils from which the plants were collected. Fifty-nine unique terminal-restriction fragments (TRFs) were identified, with more than half of the taxonomically interpreted TRFs assigned to Alpha- and Gamma-Proteobacteria and Clostridia. Data were then used to define the extent of variation of bacterial communities due to single plants or to plant populations. Results indicated a very high plant-by-plant variation of leaf-associated community (more than 93% of total variance observed). However, a core (numerically small) of plant-specific TRFs was found. This work demonstrates that plant-associated bacterial communities represent a large reservoir of biodiversity and that the high variability existing between plants, even from the same population, should be taken into account in future studies on association between bacteria and metal-hyperaccumulating plants.

Analysis of plasmid genes by phylogenetic profiling and visualization of homology relationships using Blast2Network.

BACKGROUND: Phylogenetic methods are well-established bioinformatic tools for sequence analysis, allowing to describe the non-independencies of sequences because of their common ancestor. However, the evolutionary profiles of bacterial genes are often complicated by hidden paralogy and extensive and/or (multiple) horizontal gene transfer (HGT) events which make bifurcating trees often inappropriate. In this context, plasmid sequences are paradigms of network-like relationships characterizing the evolution of prokaryotes. Actually, they can be transferred among different organisms allowing the dissemination of novel functions, thus playing a pivotal role in prokaryotic evolution. However, the study of their evolutionary dynamics is complicated by the absence of universally shared genes, a prerequisite for phylogenetic analyses. RESULTS: To overcome such limitations we developed a bioinformatic package, named Blast2Network (B2N), allowing the automatic phylogenetic profiling and the visualization of homology relationships in a large number of plasmid sequences. The software was applied to the study of 47 completely sequenced plasmids coming from Escherichia, Salmonella and Shigella spps. CONCLUSION: The tools implemented by B2N allow to describe and visualize in a new way some of the evolutionary features of plasmid molecules of Enterobacteriaceae; in particular it helped to shed some light on the complex history of Escherichia, Salmonella and Shigella plasmids and to focus on possible roles of unannotated proteins.The proposed methodology is general enough to be used for comparative genomic analyses of bacteria.

Strategy for in situ detection of natural transformation-based horizontal gene transfer events.

A strategy is described that enables the in situ detection of natural transformation in Acinetobacter baylyi BD413 by the expression of a green fluorescent protein. Microscale detection of bacterial transformants growing on plant tissues was shown by fluorescence microscopy and indicated that cultivation-based selection of transformants on antibiotic-containing agar plates underestimates transformation frequencies.

Diversity of Sinorhizobium strains nodulating Medicago sativa from different Iranian regions.

Alfalfa is believed to have originated in north-western Iran and has a long history of coexistence with its bacterial symbiont Sinorhizobium in soils of Iran. However, little is known about the diversity of Sinorhizobium strains nodulating Iranian alfalfa genotypes. In this study, Sinorhizobium populations were sampled from eight different Iranian sites using three cultivars of Medicago sativa as trap host plants. A total of 982 rhizobial strains were isolated and species were identified showing a large prevalence of Sinorhizobium meliloti over Sinorhizobium medicae. Analysis of salt tolerance demonstrated a great phenotypic diversity. The genetic diversity of the Sinorhizobium isolates was analysed using BOX-PCR and enterobacterial repetitive intergenic consensus (ERIC)-PCR. Patterns ofBOX-PCR fingerprinting were statistically analysed with AMOVA to evaluate the role of plant variety and site of origin in the genetic variance observed. Results indicated that most of the total molecular variance was attributable to divergence among strains isolated from different sites and cultivars (intrapopulation, strain-by-strain variance). Moreover, the analysis showed the presence of two geographic populations (west and northwest), indicating that the effect of the site of origin could be more relevant in shaping population genetic diversity than the effect of cultivar or individual plant.

Comparison of Highly and Weakly Virulent Dickeya solani Strains, With a View on the Pangenome and Panregulon of This Species.

Bacteria belonging to the genera Dickeya and Pectobacterium are responsible for significant economic losses in a wide variety of crops and ornamentals. During last years, increasing losses in potato production have been attributed to the appearance of Dickeya solani. The D. solani strains investigated so far share genetic homogeneity, although different virulence levels were observed among strains of various origins. The purpose of this study was to investigate the genetic traits possibly related to the diverse virulence levels by means of comparative genomics. First, we developed a new genome assembly pipeline which allowed us to complete the D. solani genomes. Four de novo sequenced and ten publicly available genomes were used to identify the structure of the D. solani pangenome, in which 74.8 and 25.2% of genes were grouped into the core and dispensable genome, respectively. For D. solani panregulon analysis, we performed a binding site prediction for four transcription factors, namely CRP, KdgR, PecS and Fur, to detect the regulons of these virulence regulators. Most of the D. solani potential virulence factors were predicted to belong to the accessory regulons of CRP, KdgR, and PecS. Thus, some differences in gene expression could exist between D. solani strains. The comparison between a highly and a low virulent strain, IFB0099 and IFB0223, respectively, disclosed only small differences between their genomes but significant differences in the production of virulence factors like pectinases, cellulases and proteases, and in their mobility. The D. solani strains also diverge in the number and size of prophages present in their genomes. Another relevant difference is the disruption of the adhesin gene fhaB2 in the highly virulent strain. Strain IFB0223, which has a complete adhesin gene, is less mobile and less aggressive than IFB0099. This suggests that in this case, mobility rather than adherence is needed in order to trigger disease symptoms. This study highlights the utility of comparative genomics in predicting D. solani traits involved in the aggressiveness of this emerging plant pathogen.

Creation and Characterization of a Genomically Hybrid Strain in the Nitrogen-Fixing Symbiotic Bacterium Sinorhizobium meliloti.

Many bacteria, often associated with eukaryotic hosts and of relevance for biotechnological applications, harbor a multipartite genome composed of more than one replicon. Biotechnologically relevant phenotypes are often encoded by genes residing on the secondary replicons. A synthetic biology approach to developing enhanced strains for biotechnological purposes could therefore involve merging pieces or entire replicons from multiple strains into a single genome. Here we report the creation of a genomic hybrid strain in a model multipartite genome species, the plant-symbiotic bacterium Sinorhizobium meliloti. We term this strain as cis-hybrid, since it is produced by genomic material coming from the same species' pangenome. In particular, we moved the secondary replicon pSymA (accounting for nearly 20% of total genome content) from a donor S. meliloti strain to an acceptor strain. The cis-hybrid strain was screened for a panel of complex phenotypes (carbon/nitrogen utilization phenotypes, intra- and extracellular metabolomes, symbiosis, and various microbiological tests). Additionally, metabolic network reconstruction and constraint-based modeling were employed for in silico prediction of metabolic flux reorganization. Phenotypes of the cis-hybrid strain were in good agreement with those of both parental strains. Interestingly, the symbiotic phenotype showed a marked cultivar-specific improvement with the cis-hybrid strains compared to both parental strains. These results provide a proof-of-principle for the feasibility of genome-wide replicon-based remodelling of bacterial strains for improved biotechnological applications in precision agriculture.

Trade, Diplomacy, and Warfare: The Quest for Elite Rhizobia Inoculant Strains.

Rhizobia form symbiotic nitrogen-fixing nodules on leguminous plants, which provides an important source of fixed nitrogen input into the soil ecosystem. The improvement of symbiotic nitrogen fixation is one of the main challenges facing agriculture research. Doing so will reduce the usage of chemical nitrogen fertilizer, contributing to the development of sustainable agriculture practices to deal with the increasing global human population. Sociomicrobiological studies of rhizobia have become a model for the study of the evolution of mutualistic interactions. The exploitation of the wide range of social interactions rhizobia establish among themselves, with the soil and root microbiota, and with the host plant, could constitute a great advantage in the development of a new generation of highly effective rhizobia inoculants. Here, we provide a brief overview of the current knowledge on three main aspects of rhizobia interaction: trade of fixed nitrogen with the plant; diplomacy in terms of communication and possible synergistic effects; and warfare, as antagonism and plant control over symbiosis. Then, we propose new areas of investigation and the selection of strains based on the combination of the genetic determinants for the relevant rhizobia symbiotic behavioral phenotypes.

Role and Regulation of ACC Deaminase Gene in Sinorhizobium meliloti: Is It a Symbiotic, Rhizospheric or Endophytic Gene?

Plant-associated bacteria exhibit a number of different strategies and specific genes allow bacteria to communicate and metabolically interact with plant tissues. Among the genes found in the genomes of plant-associated bacteria, the gene encoding the enzyme 1-aminocyclopropane-1-carboxylate (ACC) deaminase (acdS) is one of the most diffused. This gene is supposed to be involved in the cleaving of plant-produced ACC, the precursor of the plant stress-hormone ethylene toning down the plant response to infection. However, few reports are present on the actual role in rhizobia, one of the most investigated groups of plant-associated bacteria. In particular, still unclear is the origin and the role of acdS in symbiotic competitiveness and on the selective benefit it may confer to plant symbiotic rhizobia. Here we present a phylogenetic and functional analysis of acdS orthologs in the rhizobium model-species Sinorhizobium meliloti. Results showed that acdS orthologs present in S. meliloti pangenome have polyphyletic origin and likely spread through horizontal gene transfer, mediated by mobile genetic elements. When acdS ortholog from AK83 strain was cloned and assayed in S. meliloti 1021 (lacking acdS), no modulation of plant ethylene levels was detected, as well as no increase in fitness for nodule occupancy was found in the acdS-derivative strain compared to the parental one. Surprisingly, AcdS was shown to confer the ability to utilize formamide and some dipeptides as sole nitrogen source. Finally, acdS was shown to be negatively regulated by a putative leucine-responsive regulator (LrpL) located upstream to acdS sequence (acdR). acdS expression was induced by root exudates of both legumes and non-leguminous plants. We conclude that acdS in S. meliloti is not directly related to symbiotic interaction, but it could likely be involved in the rhizospheric colonization or in the endophytic behavior.

Mixed Nodule Infection in Sinorhizobium meliloti-Medicago sativa Symbiosis Suggest the Presence of Cheating Behavior.

In the symbiosis between rhizobia and legumes, host plants can form symbiotic root nodules with multiple rhizobial strains, potentially showing different symbiotic performances in nitrogen fixation. Here, we investigated the presence of mixed nodules, containing rhizobia with different degrees of mutualisms, and evaluate their relative fitness in the Sinorhizobium meliloti-Medicago sativa model symbiosis. We used three S. meliloti strains, the mutualist strains Rm1021 and BL225C and the non-mutualist AK83. We performed competition experiments involving both in vitro and in vivo symbiotic assays with M. sativa host plants. We show the occurrence of a high number (from 27 to 100%) of mixed nodules with no negative effect on both nitrogen fixation and plant growth. The estimation of the relative fitness as non-mutualist/mutualist ratios in single nodules shows that in some nodules the non-mutualist strain efficiently colonized root nodules along with the mutualist ones. In conclusion, we can support the hypothesis that in S. meliloti-M. sativa symbiosis mixed nodules are formed and allow non-mutualist or less-mutualist bacterial partners to be less or not sanctioned by the host plant, hence allowing a potential form of cheating behavior to be present in the nitrogen fixing symbiosis.

Metabolic modelling reveals the specialization of secondary replicons for niche adaptation in Sinorhizobium meliloti.

The genome of about 10% of bacterial species is divided among two or more large chromosome-sized replicons. The contribution of each replicon to the microbial life cycle (for example, environmental adaptations and/or niche switching) remains unclear. Here we report a genome-scale metabolic model of the legume symbiont Sinorhizobium meliloti that is integrated with carbon utilization data for 1,500 genes with 192 carbon substrates. Growth of S. meliloti is modelled in three ecological niches (bulk soil, rhizosphere and nodule) with a focus on the role of each of its three replicons. We observe clear metabolic differences during growth in the tested ecological niches and an overall reprogramming following niche switching. In silico examination of the inferred fitness of gene deletion mutants suggests that secondary replicons evolved to fulfil a specialized function, particularly host-associated niche adaptation. Thus, genes on secondary replicons might potentially be manipulated to promote or suppress host interactions for biotechnological purposes.