The complete replicons of 16 Ensifer meliloti strains offer insights into intra- and inter-replicon gene transfer, transposon-associated loci, and repeat elements

Guidelines for the description of rhizobial symbiovars

Rhizobia are bacteria that form nitrogen-fixing nodules in legume plants. The sets of genes responsible for both nodulation and nitrogen fixation are carried in plasmids or genomic islands that are often mobile. Different strains within a species sometimes have different host specificities, while very similar symbiosis genes may be found in strains of different species. These specificity variants are known as symbiovars, and many of them have been given names, but there are no established guidelines for defining or naming them. Here, we discuss the requirements for guidelines to describe symbiovars, propose a set of guidelines, provide a list of all symbiovars for which descriptions have been published so far, and offer a mechanism to maintain a list in the future.

Impacts of Natural Selection on Evolution of Core and Symbiotically Specialized (sym) Genes in the Polytypic Species Neorhizobium galegae

Nodule bacteria (rhizobia) represent a suitable model to address a range of fundamental genetic problems, including the impacts of natural selection on the evolution of symbiotic microorganisms. Rhizobia possess multipartite genomes in which symbiotically specialized (sym) genes differ from core genes in their natural histories. Diversification of sym genes is responsible for rhizobia microevolution, which depends on host-induced natural selection. By contrast, diversification of core genes is responsible for rhizobia speciation, which occurs under the impacts of still unknown selective factors. In this paper, we demonstrate that in goat's rue rhizobia (Neorhizobium galegae) populations collected at North Caucasus, representing two host-specific biovars orientalis and officianalis (N2-fixing symbionts of Galega orientalis and G. officinalis), the evolutionary mechanisms are different for core and sym genes. In both N. galegae biovars, core genes are more polymorphic than sym genes. In bv. orientalis, the evolution of core genes occurs under the impacts of driving selection (dN/dS > 1), while the evolution of sym genes is close to neutral (dN/dS ≈ 1). In bv. officinalis, the evolution of core genes is neutral, while for sym genes, it is dependent on purifying selection (dN/dS < 1). A marked phylogenetic congruence of core and sym genes revealed using ANI analysis may be due to a low intensity of gene transfer within and between N. galegae biovars. Polymorphism in both gene groups and the impacts of driving selection on core gene evolution are more pronounced in bv. orientalis than in bv. officianalis, reflecting the diversities of their respective host plant species. In bv. orientalis, a highly significant (P0 < 0.001) positive correlation is revealed between the p-distance and dN/dS values for core genes, while in bv. officinalis, this correlation is of low significance (0.05 < P0 < 0.10). For sym genes, the correlation between p-distance and dN/dS values is negative in bv. officinalis but is not revealed in bv. orientalis. These data, along with the functional annotation of core genes implemented using Gene Ontology tools, suggest that the evolution of bv. officinalis is based mostly on adaptation for in planta niches while in bv. orientalis, evolution presumably depends on adaptation for soil niches. New insights into the tradeoff between natural selection and genetic diversity are presented, suggesting that gene nucleotide polymorphism may be extended by driving selection only in ecologically versatile organisms capable of supporting a broad spectrum of gene alleles in their gene pools.

Discordant population structure among rhizobium divided genomes and their legume hosts

Symbiosis often occurs between partners with distinct life history characteristics and dispersal mechanisms. Many bacterial symbionts have genomes comprising multiple replicons with distinct rates of evolution and horizontal transmission. Such differences might drive differences in population structure between hosts and symbionts and among the elements of the divided genomes of bacterial symbionts. These differences might, in turn, shape the evolution of symbiotic interactions and bacterial evolution. Here we use whole genome resequencing of a hierarchically structured sample of 191 strains of Sinorhizobium meliloti collected from 21 locations in southern Europe to characterize population structures of this bacterial symbiont, which forms a root nodule symbiosis with the host plant Medicago truncatula. S. meliloti genomes showed high local (within-site) variation and little isolation by distance. This was particularly true for the two symbiosis elements, pSymA and pSymB, which have population structures that are similar to each other, but distinct from both the bacterial chromosome and the host plant. Given limited recombination on the chromosome, compared to the symbiosis elements, distinct population structures may result from differences in effective gene flow. Alternatively, positive or purifying selection, with little recombination, may explain distinct geographical patterns at the chromosome. Discordant population structure between hosts and symbionts indicates that geographically and genetically distinct host populations in different parts of the range might interact with genetically similar symbionts, potentially minimizing local specialization.

Host-Associated Rhizobial Fitness: Dependence on Nitrogen, Density, Community Complexity, and Legume Genotype

The environmental context of the nitrogen-fixing mutualism between leguminous plants and rhizobial bacteria varies over space and time. Variation in resource availability, population density, and composition likely affect the ecology and evolution of rhizobia and their symbiotic interactions with hosts. We examined how host genotype, nitrogen addition, rhizobial density, and community complexity affected selection on 68 rhizobial strains in the Sinorhizobium meliloti–Medicago truncatula mutualism. As expected, host genotype had a substantial effect on the size, number, and strain composition of root nodules (the symbiotic organ). The understudied environmental variable of rhizobial density had a stronger effect on nodule strain frequency than the addition of low nitrogen levels. Higher inoculum density resulted in a nodule community that was less diverse and more beneficial but only in the context of the more selective host genotype. Higher density resulted in more diverse and less beneficial nodule communities with the less selective host. Density effects on strain composition deserve additional scrutiny as they can create feedback between ecological and evolutionary processes. Finally, we found that relative strain rankings were stable across increasing community complexity (2, 3, 8, or 68 strains). This unexpected result suggests that higher-order interactions between strains are rare in the context of nodule formation and development. Our work highlights the importance of examining mechanisms of density-dependent strain fitness and developing theoretical predictions that incorporate density dependence. Furthermore, our results have translational relevance for overcoming establishment barriers in bioinoculants and motivating breeding programs that maintain beneficial plant-microbe interactions across diverse agroecological contexts.

IMPORTANCE Legume crops establish beneficial associations with rhizobial bacteria that perform biological nitrogen fixation, providing nitrogen to plants without the economic and greenhouse gas emission costs of chemical nitrogen inputs. Here, we examine the influence of three environmental factors that vary in agricultural fields on strain relative fitness in nodules. In addition to manipulating nitrogen, we also use two biotic variables that have rarely been examined: the rhizobial community's density and complexity. Taken together, our results suggest that (i) breeding legume varieties that select beneficial strains despite environmental variation is possible, (ii) changes in rhizobial population densities that occur routinely in agricultural fields could drive evolutionary changes in rhizobial populations, and (iii) the lack of higher-order interactions between strains will allow the high-throughput assessments of rhizobia winners and losers during plant interactions.

Combining GWAS and population genomic analyses to characterize coevolution in a legume-rhizobia symbiosis

The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two-species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single-strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype-by-genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation-selection balance. Consistent with the lack of positive selection acting on host genes, we found that among-host variation in growth was similar whether plants were grown with rhizobia or N-fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

Diversification of plasmids in a genus of pathogenic and nitrogen-fixing bacteria

Members of the agrobacteria-rhizobia complex (ARC) have multiple and diverse plasmids. The extent to which these plasmids are shared and the consequences of their interactions are not well understood. We extracted over 4000 plasmid sequences from 1251 genome sequences and constructed a network to reveal interactions that have shaped the evolutionary histories of oncogenic virulence plasmids. One newly discovered type of oncogenic plasmid is a mosaic with three incomplete, but complementary and partially redundant virulence loci. Some types of oncogenic plasmids recombined with accessory plasmids or acquired large regions not known to be associated with pathogenicity. We also identified two classes of partial virulence plasmids. One class is potentially capable of transforming plants, but not inciting disease symptoms. Another class is inferred to be incomplete and non-functional but can be found as coresidents of the same strain and together are predicted to confer pathogenicity. The modularity and capacity for some plasmids to be transmitted broadly allow them to diversify, convergently evolve adaptive plasmids and shape the evolution of genomes across much of the ARC. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Pervasive RNA Regulation of Metabolism Enhances the Root Colonization Ability of Nitrogen-Fixing Symbiotic α-Rhizobia

The rhizosphere and rhizoplane are nutrient-rich but selective environments for the root microbiome. Here, we deciphered a posttranscriptional network regulated by the homologous trans-small RNAs (sRNAs) AbcR1 and AbcR2, which rewire the metabolism of the nitrogen-fixing α-rhizobium Sinorhizobium meliloti during preinfection stages of symbiosis with its legume host alfalfa. The LysR-type regulator LsrB, which transduces the cell redox state, is indispensable for AbcR1 expression in actively dividing bacteria, whereas the stress-induced transcription of AbcR2 depends on the alternative σ factor RpoH1. MS2 affinity purification coupled with RNA sequencing unveiled exceptionally large and overlapping AbcR1/2 mRNA interactomes, jointly representing ⁓6% of the S. meliloti protein-coding genes. Most mRNAs encode transport/metabolic proteins whose translation is silenced by base pairing to two distinct anti-Shine Dalgarno motifs that function independently in both sRNAs. A metabolic model-aided analysis of the targetomes predicted changes in AbcR1/2 expression driven by shifts in carbon/nitrogen sources, which were confirmed experimentally. Low AbcR1/2 levels in some defined media anticipated overexpression growth phenotypes linked to the silencing of specific mRNAs. As a proof of principle, we confirmed AbcR1/2-mediated downregulation of the l-amino acid AapQ permease. AbcR1/2 interactomes are well represented in rhizosphere-related S. meliloti transcriptomic signatures. Remarkably, a lack of AbcR1 specifically compromised the ability of S. meliloti to colonize the root rhizoplane. The AbcR1 regulon likely ranks the utilization of available substrates to optimize metabolism, thus conferring on S. meliloti an advantage for efficient rhizosphere/rhizoplane colonization. AbcR1 regulation is predicted to be conserved in related α-rhizobia, which opens unprecedented possibilities for engineering highly competitive biofertilizers. IMPORTANCE Nitrogen-fixing root nodule symbioses between rhizobia and legume plants provide more than half of the combined nitrogen incorporated annually into terrestrial ecosystems, rendering plant growth independent of environmentally unfriendly chemical fertilizers. The success of symbiosis depends primarily on the capacity of rhizobia to establish competitive populations in soil and rhizosphere environments. Here, we provide insights into the regulation and architecture of an extensive RNA posttranscriptional network that fine-tunes the metabolism of the alfalfa symbiont S. meliloti, thereby enhancing the ability of this beneficial bacterium to colonize nutrient-rich but extremely selective niches, such as the rhizosphere of its host plant. This pervasive RNA regulation of metabolism is a major adaptive mechanism, predicted to operate in diverse rhizobial species. Because RNA regulation relies on modifiable base-pairing interactions, our findings open unexplored avenues for engineering the legumes rhizobiome within sustainable agricultural practices.

Complete Genome Sequence of Sinorhizobium meliloti Strain AK21, a Salt-Tolerant Isolate from the Aral Sea Region

We report here the complete genome sequence of the salt-tolerant Sinorhizobium meliloti strain AK21, isolated from nodules of Medicago sativa L. subsp. ambigua inhabiting the northern Aral Sea Region. This genome (7.36 Mb) consists of a chromosome and four accessory plasmids, two of which are the symbiotic megaplasmids pSymA and pSymB.

We report here the complete genome sequence of the salt-tolerant Sinorhizobium meliloti strain AK21, isolated from nodules of Medicago sativa L. subsp. ambigua inhabiting the northern Aral Sea Region. This genome (7.36 Mb) consists of a chromosome and four accessory plasmids, two of which are the symbiotic megaplasmids pSymA and pSymB.

A Select and Resequence Approach Reveals Strain-Specific Effects of Medicago Nodule-Specific PLAT-Domain Genes

Genetic studies of legume symbiosis with nitrogen-fixing rhizobial bacteria have traditionally focused on nodule and nitrogen-fixation phenotypes when hosts are inoculated with a single rhizobial strain. These approaches overlook the potential effect of host genes on rhizobial fitness (i.e. how many rhizobia are released from host nodules) and strain-specific effects of host genes (i.e. genome × genome interactions). Using Medicago truncatula mutants in the recently described nodule-specific PLAT domain (NPD) gene family, we show how inoculating plants with a mixed inoculum of 68 rhizobial strains (Ensifer meliloti) via a select-and-resequence approach can be used to efficiently assay host mutants for strain-specific effects of late-acting host genes on interacting bacteria. The deletion of a single NPD gene (npd2) or all five members of the NPD gene family (npd1-5) differentially altered the frequency of rhizobial strains in nodules even though npd2 mutants had no visible nodule morphology or N-fixation phenotype. Also, npd1-5 nodules were less diverse and had larger populations of colony-forming rhizobia despite their smaller size. Lastly, NPD mutations disrupt a positive correlation between strain fitness and wild-type host biomass. These changes indicate that the effects of NPD proteins are strain dependent and that NPD family members are not redundant with regard to their effects on rhizobial strains. Association analyses of the rhizobial strains in the mixed inoculation indicate that rhizobial genes involved in chromosome segregation, cell division, GABA metabolism, efflux systems, and stress tolerance play an important role in the strain-specific effects of NPD genes.

Legacy of prior host and soil selection on rhizobial fitness in planta

Measuring selection acting on microbial populations in natural or even seminatural environments is challenging because many microbial populations experience variable selection. The majority of rhizobial bacteria are found in the soil. However, they also live symbiotically inside nodules of legume hosts and each nodule can release thousands of daughter cells back into the soil. We tested how past selection (i.e., legacies) by two plant genotypes and by the soil alone affected selection and genetic diversity within a population of 101 strains of Ensifer meliloti. We also identified allelic variants most strongly associated with soil- and host-dependent fitness. In addition to imposing direct selection on rhizobia populations, soil and host environments had lasting effects across host generations. Host presence and genotype during the legacy period explained 22% and 12% of the variance in the strain composition of nodule communities in the second cohort, respectively. Although strains with high host fitness in the legacy cohort tended to be enriched in the second cohort, the diversity of the strain community was greater when the second cohort was preceded by host rather than soil legacies. Our results indicate the potential importance of soil selection driving the evolution of these plant-associated microbes.

Codon Usage Heterogeneity in the Multipartite Prokaryote Genome: Selection-Based Coding Bias Associated with Gene Location, Expression Level, and Ancestry

Prokaryotes represent an ancestral lineage in the tree of life and constitute optimal resources for investigating the evolution of genomes in unicellular organisms. Many bacterial species possess multipartite genomes offering opportunities to study functional variations among replicons, how and where new genes integrate into a genome, and how genetic information within a lineage becomes encoded and evolves. To analyze these issues, we focused on the model soil bacterium Sinorhizobium meliloti, which harbors a chromosome, a chromid (pSymB), a megaplasmid (pSymA), and, in many strains, one or more accessory plasmids. The analysis of several genomes, together with 1.4 Mb of accessory plasmid DNA that we purified and sequenced, revealed clearly different functional profiles associated with each genomic entity. pSymA, in particular, exhibited remarkable interstrain variation and a high density of singletons (unique, exclusive genes) featuring functionalities and modal codon usages that were very similar to those of the plasmidome. All this evidence reinforces the idea of a close relationship between pSymA and the plasmidome. Correspondence analyses revealed that adaptation of codon usages to the translational machinery increased from plasmidome to pSymA to pSymB to chromosome, corresponding as such to the ancestry of each replicon in the lineage. We demonstrated that chromosomal core genes gradually adapted to the translational machinery, reminiscent of observations in several bacterial taxa for genes with high expression levels. Such findings indicate a previously undiscovered codon usage adaptation associated with the chromosomal core information that likely operates to improve bacterial fitness. We present a comprehensive model illustrating the central findings described here, discussed in the context of the changes occurring during the evolution of a multipartite prokaryote genome.IMPORTANCE Bacterial genomes usually include many thousands of genes which are expressed with diverse spatial-temporal patterns and intensities. A well-known evidence is that highly expressed genes, such as the ribosomal and other translation-related proteins (RTRPs), have accommodated their codon usage to optimize translation efficiency and accuracy. Using a bioinformatic approach, we identify core-genes sets with different ancestries, and demonstrate that selection processes that optimize codon usage are not restricted to RTRPs but extended at a genome-wide scale. Such findings highlight, for the first time, a previously undiscovered adaptation strategy associated with the chromosomal-core information. Contrasted with the translationally more adapted genes, singletons (i.e., exclusive genes, including those of the plasmidome) appear as the gene pool with the less-ameliorated codon usage in the lineage. A comprehensive summary describing the inter- and intra-replicon heterogeneity of codon usages in a complex prokaryote genome is presented.

Genome-Wide Association Analyses in the Model Rhizobium Ensifer meliloti

Genome-wide association studies (GWAS) can identify genetic variants responsible for naturally occurring and quantitative phenotypic variation. Association studies therefore provide a powerful complement to approaches that rely on de novo mutations for characterizing gene function. Although bacteria should be amenable to GWAS, few GWAS have been conducted on bacteria, and the extent to which nonindependence among genomic variants (e.g., linkage disequilibrium [LD]) and the genetic architecture of phenotypic traits will affect GWAS performance is unclear. We apply association analyses to identify candidate genes underlying variation in 20 biochemical, growth, and symbiotic phenotypes among 153 strains of Ensifer meliloti For 11 traits, we find genotype-phenotype associations that are stronger than expected by chance, with the candidates in relatively small linkage groups, indicating that LD does not preclude resolving association candidates to relatively small genomic regions. The significant candidates show an enrichment for nucleotide polymorphisms (SNPs) over gene presence-absence variation (PAV), and for five traits, candidates are enriched in large linkage groups, a possible signature of epistasis. Many of the variants most strongly associated with symbiosis phenotypes were in genes previously identified as being involved in nitrogen fixation or nodulation. For other traits, apparently strong associations were not stronger than the range of associations detected in permuted data. In sum, our data show that GWAS in bacteria may be a powerful tool for characterizing genetic architecture and identifying genes responsible for phenotypic variation. However, careful evaluation of candidates is necessary to avoid false signals of association.IMPORTANCE Genome-wide association analyses are a powerful approach for identifying gene function. These analyses are becoming commonplace in studies of humans, domesticated animals, and crop plants but have rarely been conducted in bacteria. We applied association analyses to 20 traits measured in Ensifer meliloti, an agriculturally and ecologically important bacterium because it fixes nitrogen when in symbiosis with leguminous plants. We identified candidate alleles and gene presence-absence variants underlying variation in symbiosis traits, antibiotic resistance, and use of various carbon sources; some of these candidates are in genes previously known to affect these traits whereas others were in genes that have not been well characterized. Our results point to the potential power of association analyses in bacteria, but also to the need to carefully evaluate the potential for false associations.