Population mixing of Rhizobium leguminosarum bv. viciae nodulating Vicia faba: the role of recombination and lateral gene transfer

Recombination facilitates adaptive evolution in rhizobial soil bacteria

Homologous recombination is expected to increase natural selection efficacy by decoupling the fate of beneficial and deleterious mutations and by readily creating new combinations of beneficial alleles. Here, we investigate how the proportion of amino acid substitutions fixed by adaptive evolution (α) depends on the recombination rate in bacteria. We analyze 3,086 core protein-coding sequences from 196 genomes belonging to five closely related species of the genus Rhizobium. These genes are found in all species and do not display any signs of introgression between species. We estimate α using the site frequency spectrum (SFS) and divergence data for all pairs of species. We evaluate the impact of recombination within each species by dividing genes into three equally sized recombination classes based on their average level of intragenic linkage disequilibrium. We find that α varies from 0.07 to 0.39 across species and is positively correlated with the level of recombination. This is both due to a higher estimated rate of adaptive evolution and a lower estimated rate of nonadaptive evolution, suggesting that recombination both increases the fixation probability of advantageous variants and decreases the probability of fixation of deleterious variants. Our results demonstrate that homologous recombination facilitates adaptive evolution measured by α in the core genome of prokaryote species in agreement with studies in eukaryotes.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Phylogeography and Symbiotic Effectiveness of Rhizobia Nodulating Chickpea (Cicer arietinum L.) in Ethiopia

Chickpea (Cicer arietinum L.) used to be considered a restrictive host that nodulated and fixed nitrogen only with Mesorhizobium ciceri and M. mediterraneum. Recent analysis revealed that chickpea can also establish effective symbioses with strains of several other Mesorhizobium species such as M. loti, M. haukuii, M. amorphae, M. muleiense, etc. These strains vary in their nitrogen fixation potential inviting further exploration. We characterized newly collected mesorhizobial strains isolated from various locations in Ethiopia to evaluate genetic diversity, biogeographic structure and symbiotic effectiveness. Symbiotic effectiveness was evaluated in Leonard Jars using a locally released chickpea cultivar "Nattoli". Most of the new isolates belonged to a clade related to M. plurifarium, with very few sequence differences, while the total collection of strains contained three additional mesorhizobial genospecies associated with M. ciceri, M. abyssinicae and an unidentified Mesorhizobium species isolated from a wild host in Eritrea. The four genospecies identified represented a subset of the eight major Mesorhizobium clades recently reported for Ethiopia based on metagenomic data. All Ethiopian strains had nearly identical symbiotic genes that grouped them in a single cluster with M. ciceri, M. mediterraneum and M. muleiense, but not with M. plurifarium. Some phylogeographic structure was observed, with elevation and geography explaining some of the genetic differences among strains, but the relation between genetic identity and symbiotic effectiveness was observed to be weak.

Phylogeographic distribution of rhizobia nodulating common bean (Phaseolus vulgaris L.) in Ethiopia

Rhizobia are soilborne bacteria that form symbiotic relations with legumes and fix atmospheric nitrogen. The nitrogen fixation potential depends on several factors such as the type of host and symbionts and on environmental factors that affect the distribution of rhizobia. We isolated bacteria nodulating common bean in Southern Ethiopia to evaluate their genetic diversity and phylogeography at nucleotide, locus (gene/haplotype) and species levels of genetic hierarchy. Phylogenetically, eight rhizobial genospecies (including previous collections) were determined that had less genetic diversity than found among reference strains. The limited genetic diversity of the Ethiopian collections was due to absence of many of the Rhizobium lineages known to nodulate beans. Rhizobium etli and Rhizobiumphaseoli were predominant strains of bean-nodulating rhizobia in Ethiopia. We found no evidence for a phylogeographic pattern in strain distribution. However, joint analysis of the current and previous collections revealed differences between the two collections at nucleotide level of genetic hierarchy. The differences were due to genospecies Rhizobium aethiopicum that was only isolated in the earlier collection.

Host-specific competitiveness to form nodules in Rhizobium leguminosarum symbiovar viciae

Fabeae legumes such as pea and faba bean form symbiotic nodules with a large diversity of soil Rhizobium leguminosarum symbiovar viciae (Rlv) bacteria. However, bacteria competitive to form root nodules (CFN) are generally not the most efficient to fix dinitrogen, resulting in a decrease in legume crop yields. Here, we investigate differential selection by host plants on the diversity of Rlv. A large collection of Rlv was collected by nodule trapping with pea and faba bean from soils at five European sites. Representative genomes were sequenced. In parallel, diversity and abundance of Rlv were estimated directly in these soils using metabarcoding. The CFN of isolates was measured with both legume hosts. Pea/faba bean CFN were associated to Rlv genomic regions. Variations of bacterial pea and/or faba bean CFN explained the differential abundance of Rlv genotypes in pea and faba bean nodules. No evidence was found for genetic association between CFN and variations in the core genome, but variations in specific regions of the nod locus, as well as in other plasmid loci, were associated with differences in CFN. These findings shed light on the genetic control of CFN in Rlv and emphasise the importance of host plants in controlling Rhizobium diversity.

New Insight into the Evolution of Symbiotic Genes in Black Locust-Associated Rhizobia

Nitrogen fixation in legumes occurs via symbiosis with rhizobia. This process involves packages of symbiotic genes on mobile genetic elements that are readily transferred within or between rhizobial species, furnishing the recipient with the ability to interact with plant hosts. However, it remains elusive whether plant host migration has played a role in shaping the current distribution of genetic variation in symbiotic genes. Herein, we examined the genetic structure and phylogeographic pattern of symbiotic genes in 286 symbiotic strains of Mesorhizobium nodulating black locust (Robinia pseudoacacia), a cross-continental invasive legume species that is native to North America. We conducted detailed phylogeographic analysis and approximate Bayesian computation to unravel the complex demographic history of five key symbiotic genes. The sequencing results indicate an origin of symbiotic genes in Germany rather than North America. Our findings provide strong evidence of prehistoric lineage splitting and spatial expansion events resulting in multiple radiations of descendent clones from founding sequence types worldwide. Estimates of the timescale of divergence in North American and Chinese subclades suggest that black locust-specific symbiotic genes have been present in these continent many thousands of years before recent migration of plant host. Although numerous crop plants, including legumes, have found their centers of origin as centers of evolution and diversity, the number of legume-specific symbiotic genes with a known geographic origin is limited. This work sheds light on the coevolution of legumes and rhizobia.

Electrophoretic profiles of lipopolysaccharides from Rhizobium strains nodulating Pisum sativum do not reflect phylogenetic relationships between these strains

Rhizobia that nodulate peas comprise a heterogeneous group of bacteria. The aim of this study was to investigate the relationship between phylogeny and electrophoretic and hydroxy fatty acid lipopolysaccharide (LPS) profiles of pea microsymbionts. Based on amplified fragment length polymorphism (AFLP) fingerprinting data, the pea microsymbionts were grouped into two clusters distinguished at 58% similarity level. Based on the concatenated 16S rRNA, recA, and atpD housekeeping gene data, the microsymbionts appeared to be most closely related to Rhizobium leguminosarum biovars viciae and trifolii. Applying cluster analysis to their LPS electrophoretic profiles, the strains were assigned to two major groups with different banding patterns. All hydroxy fatty acids common to R. leguminosarum and R. etli were detected in each examined strain. Differences in the proportions of 3- to ω-1 hydroxy fatty acids allowed us to distinguish two groups of strains. This classification did not overlap with one based on LPS electrophoretic profiles. No clear correlation was apparent between the genetic traits and LPS profiles of the pea nodule isolates.

Genetic Diversity and Symbiotic Efficiency of Nodulating Rhizobia Isolated from Root Nodules of Faba Bean in One Field

Thirty-one nodulating rhizobium strains were collected from root nodules of spring and winter type faba bean cultivars grown in micro ecoarea, i.e. the same field in Chengdu plain, China. The symbiotic efficiency and phylogeny of these strains were studied. Effectively nitrogen fixing strains were isolated from both winter type and spring type cultivars. Based on phylogenetic analysis of 16S rRNA gene and concatenated sequence of atpD, glnII and recA genes, the isolates were assigned as Rhizobium anhuiense and a potential new Rhizobium species. The isolates were diverse on symbiosis related gene level, carrying five, four and three variants of nifH, nodC and nodD, respectively. Strains carrying similar gene combinations were trapped by both winter and spring cultivars, disagreeing with the specificity of symbiotic genotypes to reported earlier faba bean ecotypes.

The Plasmid Mobilome of the Model Plant-Symbiont Sinorhizobium meliloti: Coming up with New Questions and Answers

Rhizobia are Gram-negative Alpha- and Betaproteobacteria living in the underground which have the ability to associate with legumes for the establishment of nitrogen-fixing symbioses. Sinorhizobium meliloti in particular-the symbiont of Medicago, Melilotus, and Trigonella spp.-has for the past decades served as a model organism for investigating, at the molecular level, the biology, biochemistry, and genetics of a free-living and symbiotic soil bacterium of agricultural relevance. To date, the genomes of seven different S. meliloti strains have been fully sequenced and annotated, and several other draft genomic sequences are also available. The vast amount of plasmid DNA that S. meliloti frequently bears (up to 45% of its total genome), the conjugative ability of some of those plasmids, and the extent of the plasmid diversity has provided researchers with an extraordinary system to investigate functional and structural plasmid molecular biology within the evolutionary context surrounding a plant-associated model bacterium. Current evidence indicates that the plasmid mobilome in S. meliloti is composed of replicons varying greatly in size and having diverse conjugative systems and properties along with different evolutionary stabilities and biological roles. While plasmids carrying symbiotic functions (pSyms) are known to have high structural stability (approaching that of chromosomes), the remaining plasmid mobilome (referred to as the non-pSym, functionally cryptic, or accessory compartment) has been shown to possess remarkable diversity and to be highly active in conjugation. In light of the modern genomic and current biochemical data on the plasmids of S. meliloti, the current article revises their main structural components, their transfer and regulatory mechanisms, and their potential as vehicles in shaping the evolution of the rhizobial genome.

Recombination and horizontal transfer of nodulation and ACC deaminase (acdS) genes within Alpha- and Betaproteobacteria nodulating legumes of the Cape Fynbos biome

The goal of this work is to study the evolution and the degree of horizontal gene transfer (HGT) within rhizobial genera of both Alphaproteobacteria (Mesorhizobium, Rhizobium) and Betaproteobacteria (Burkholderia), originating from South African Fynbos legumes. By using a phylogenetic approach and comparing multiple chromosomal and symbiosis genes, we revealed conclusive evidence of high degrees of horizontal transfer of nodulation genes among closely related species of both groups of rhizobia, but also among species with distant genetic backgrounds (Rhizobium and Mesorhizobium), underscoring the importance of lateral transfer of symbiosis traits as an important evolutionary force among rhizobia of the Cape Fynbos biome. The extensive exchange of symbiosis genes in the Fynbos is in contrast with a lack of significant events of HGT among Burkholderia symbionts from the South American Cerrado and Caatinga biome. Furthermore, homologous recombination among selected housekeeping genes had a substantial impact on sequence evolution within Burkholderia and Mesorhizobium. Finally, phylogenetic analyses of the non-symbiosis acdS gene in Mesorhizobium, a gene often located on symbiosis islands, revealed distinct relationships compared to the chromosomal and symbiosis genes, suggesting a different evolutionary history and independent events of gene transfer. The observed events of HGT and incongruence between different genes necessitate caution in interpreting topologies from individual data types.

The Diversity of Pea Microsymbionts in Various Types of Soils and Their Effects on Plant Host Productivity

The growth and yield of peas cultivated on eight different soils, as well as the diversity of pea microsymbionts derived from these soils were investigated in the present study. The experimental plot was composed of soils that were transferred from different parts of Poland more than a century ago. The soils were located in direct vicinity of each other in the experimental plot. All soils examined contained pea microsymbionts, which were suggested to belong to Rhizobium leguminosarum sv. viciae based on the nucleotide sequence of the partial 16S rRNA gene. PCR-RFLP analyses of the 16S-23S rRNA gene ITS region and nodD alleles revealed the presence of numerous and diversified groups of pea microsymbionts and some similarities between the tested populations, which may have been the result of the spread or displacement of strains. However, most populations retained their own genetic distinction, which may have been related to the type of soil. Most of the tested populations comprised low-effective strains for the promotion of pea growth. No relationships were found between the characteristics of soil and symbiotic effectiveness of rhizobial populations; however, better seed yield was obtained for soil with medium biological productivity inhabited by high-effective rhizobial populations than for soil with high agricultural quality containing medium-quality pea microsymbionts, and these results showed the importance of symbiosis for plant hosts.

Bacterial genospecies that are not ecologically coherent: population genomics of Rhizobium leguminosarum

Biological species may remain distinct because of genetic isolation or ecological adaptation, but these two aspects do not always coincide. To establish the nature of the species boundary within a local bacterial population, we characterized a sympatric population of the bacterium Rhizobium leguminosarum by genomic sequencing of 72 isolates. Although all strains have 16S rRNA typical of R. leguminosarum, they fall into five genospecies by the criterion of average nucleotide identity (ANI). Many genes, on plasmids as well as the chromosome, support this division: recombination of core genes has been largely within genospecies. Nevertheless, variation in ecological properties, including symbiotic host range and carbon-source utilization, cuts across these genospecies, so that none of these phenotypes is diagnostic of genospecies. This phenotypic variation is conferred by mobile genes. The genospecies meet the Mayr criteria for biological species in respect of their core genes, but do not correspond to coherent ecological groups, so periodic selection may not be effective in purging variation within them. The population structure is incompatible with traditional ‘polyphasic taxonomy′ that requires bacterial species to have both phylogenetic coherence and distinctive phenotypes. More generally, genomics has revealed that many bacterial species share adaptive modules by horizontal gene transfer, and we envisage a more consistent taxonomic framework that explicitly recognizes this. Significant phenotypes should be recognized as ‘biovars' within species that are defined by core gene phylogeny.

Genetic divergence of bradyrhizobium strains nodulating soybeans as revealed by multilocus sequence analysis of genes inside and outside the symbiosis island

The genus Bradyrhizobium has been considered to be a taxonomically difficult group. In this study, phylogenetics and evolutionary genetics analyses were used to investigate divergence levels among Bradyrhizobium strains nodulating soybeans in China. Eleven genospecies were identified by sequence analysis of three phylogenetic and taxonomic markers (SMc00019, thrA, and truA). This was also supported by analyses of eight genes outside the symbiosis island ("off-island" genes; SMc00019, thrA, truA, fabB, glyA, phyR, exoN, and hsfA). However, seven genes inside the symbiosis island ("island" genes; nifA, nifH, nodC, nodV, fixA, trpD, and rhcC2) showed contrasting lower levels of nucleotide diversity and recombination rates than did off-island genes. Island genes had significantly incongruent gene phylogenies compared to the species tree. Four phylogenetic clusters were observed in island genes, and the epidemic cluster IV (harbored by Bradyrhizobium japonicum, Bradyrhizobium diazoefficiens, Bradyrhizobium huanghuaihaiense, Bradyrhizobium liaoningense, Bradyrhizobium daqingense, Bradyrhizobium sp. I, Bradyrhizobium sp. III, and Bradyrhizobium sp. IV) was not found in Bradyrhizobium yuanmingense, Bradyrhizobium sp. II, or Bradyrhizobium elkanii. The gene flow level of island genes among genospecies is discussed in the context of the divergence level of off-island genes.

Replicon-dependent differentiation of symbiosis-related genes in Sinorhizobium strains nodulating Glycine max

In order to investigate the genetic differentiation of Sinorhizobium strains nodulating Glycine max and related microevolutionary mechanisms, three housekeeping genes (SMc00019, truA, and thrA) and 16 symbiosis-related genes on the chromosome (7 genes), pSymA (6 genes), and pSymB (3 genes) were analyzed. Five distinct species were identified among the test strains by calculating the average nucleotide identity (ANI) of SMc00019-truA-thrA: Sinorhizobium fredii, Sinorhizobium sojae, Sinorhizobium sp. I, Sinorhizobium sp. II, and Sinorhizobium sp. III. These species assignments were also supported by population genetics and phylogenetic analyses of housekeeping genes and symbiosis-related genes on the chromosome and pSymB. Different levels of genetic differentiation were observed among these species or different replicons. S. sojae was the most divergent from the other test species and was characterized by its low intraspecies diversity and limited geographic distribution. Intergenic recombination dominated the evolution of 19 genes from different replicons. Intraspecies recombination happened frequently in housekeeping genes and symbiosis-related genes on the chromosome and pSymB, whereas pSymA genes showed a clear pattern of lateral-transfer events between different species. Moreover, pSymA genes were characterized by a lower level of polymorphism and recombination than those on the chromosome and pSymB. Taken together, genes from different replicons of rhizobia might be involved in the establishment of symbiosis with legumes, but these symbiosis-related genes might have evolved differently according to their corresponding replicons.

Microevolution of symbiotic Bradyrhizobium populations associated with soybeans in east North America

Microevolution and origins of Bradyrhizobium populations associated with soybeans at two field sites (A and B, 280 km apart in Canada) with contrasting histories of inoculation was investigated using probabilistic analyses of six core (housekeeping) gene sequences. These analyses supported division of 220 isolates in five lineages corresponding either to B. japonicum groups 1 and 1a or to one of three novel lineages within the genus Bradyrhizobium. None of the isolates from site A and about 20% from site B (the only site with a recent inoculation history) were attributed to inoculation sources. The data suggest that most isolates were of indigenous origin based on sequence analysis of 148 isolates of soybean-nodulating bacteria from native legumes (Amphicarpaea bracteata and Desmodium canadense). Isolates from D. canadense clustered with B. japonicum group 1, whereas those from A. bracteata were placed in two novel lineages encountered at soybean field sites. One of these novel lineages predominated at soybean sites and exhibited a significant clonal expansion likely reflecting selection by the plant host. Homologous recombination events detected in the 35 sequence types from soybean sites had an effect on genetic diversification that was approximately equal to mutation. Interlineage transfer of core genes was infrequent and mostly attributable to gyrB that had a history of frequent recombination. Symbiotic gene sequences (nodC and nifH) of isolates from soybean sites and native legumes clustered in two lineages corresponding to B. japonicum and B. elkani with the inheritance of these genes appearing predominantly by vertical transmission. The data suggest that soybean-nodulating bacteria associated with native legumes represent a novel source of ecologically adapted bacteria for soybean inoculation.

Robust markers reflecting phylogeny and taxonomy of rhizobia

Genomic ANI (Average Nucleotide Identity) has been found to be able to replace DNA-DNA hybridization in prokaryote taxonomy. The ANI of each of the core genes that has a phylogeny congruent with the reference species tree of rhizobia was compared to the genomic ANI. This allowed us to identify three housekeeping genes (SMc00019-truA-thrA) whose ANI reflected the intraspecies and interspecies genomic ANI among rhizobial strains, revealing an ANI gap (≥2%) between the inter- and intra-species comparisons. The intraspecies (96%) and interspecies (94%) ANI boundaries calculated from three genes (SMc00019-truA-thrA) provided a criterion for bacterial species definition and confirmed 621/629 of known interspecies relationships within Bradyrhizobium, Mesorhizobium, Sinorhizobium and Rhizobium. Some widely studied strains should be renamed. The SMc00019-truA-thrA ANI also correlates well with the genomic ANI of strains in Agrobacterium, Methylobacterium, Ralstonia, Rhodopseudomonas, Cupriavidus and Burkholderia, suggesting their wide applicability in other bacteria.

Comparative genomics of rhizobia nodulating soybean suggests extensive recruitment of lineage-specific genes in adaptations

The rhizobium-legume symbiosis has been widely studied as the model of mutualistic evolution and the essential component of sustainable agriculture. Extensive genetic and recent genomic studies have led to the hypothesis that many distinct strategies, regardless of rhizobial phylogeny, contributed to the varied rhizobium-legume symbiosis. We sequenced 26 genomes of Sinorhizobium and Bradyrhizobium nodulating soybean to test this hypothesis. The Bradyrhizobium core genome is disproportionally enriched in lipid and secondary metabolism, whereas several gene clusters known to be involved in osmoprotection and adaptation to alkaline pH are specific to the Sinorhizobium core genome. These features are consistent with biogeographic patterns of these bacteria. Surprisingly, no genes are specifically shared by these soybean microsymbionts compared with other legume microsymbionts. On the other hand, phyletic patterns of 561 known symbiosis genes of rhizobia reflected the species phylogeny of these soybean microsymbionts and other rhizobia. Similar analyses with 887 known functional genes or the whole pan genome of rhizobia revealed that only the phyletic distribution of functional genes was consistent with the species tree of rhizobia. Further evolutionary genetics revealed that recombination dominated the evolution of core genome. Taken together, our results suggested that faithfully vertical genes were rare compared with those with history of recombination including lateral gene transfer, although rhizobial adaptations to symbiotic interactions and other environmental conditions extensively recruited lineage-specific shell genes under direct or indirect control through the speciation process.

The mosaicism of plasmids revealed by atypical genes detection and analysis

BACKGROUND

From an evolutionary viewpoint, prokaryotic genomes are extremely plastic and dynamic, since large amounts of genetic material are continuously added and/or lost through promiscuous gene exchange. In this picture, plasmids play a key role, since they can be transferred between different cells and, through genetic rearrangement(s), undergo gene(s) load, leading, in turn, to the appearance of important metabolic innovations that might be relevant for cell life. Despite their central position in bacterial evolution, a massive analysis of newly acquired functional blocks [likely the result of horizontal gene transfer (HGT) events] residing on plasmids is still missing.

RESULTS

We have developed a computational, composition-based, pipeline to scan almost 2000 plasmids for genes that differ significantly from their hosting molecule. Plasmids atypical genes (PAGs) were about 6% of the total plasmids ORFs and, on average, each plasmid possessed 4.4 atypical genes. Nevertheless, conjugative plasmids were shown to possess an amount of atypical genes than that found in not mobilizable plasmids, providing strong support for the central role suggested for conjugative plasmids in the context of HGT. Part of the retrieved PAGs are organized into (mainly short) clusters and are involved in important biological processes (detoxification, antibiotic resistance, virulence), revealing the importance of HGT in the spreading of metabolic pathways within the whole microbial community. Lastly, our analysis revealed that PAGs mainly derive from other plasmid (rather than coming from phages and/or chromosomes), suggesting that plasmid-plasmid DNA exchange might be the primary source of metabolic innovations in this class of mobile genetic elements.

CONCLUSIONS

In this work we have performed the first large scale analysis of atypical genes that reside on plasmid molecules to date. Our findings on PAGs function, organization, distribution and spreading reveal the importance of plasmids-mediated HGT within the complex bacterial evolutionary network and in the dissemination of important biological traits.