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.

Herbivory modifies plant symbiont number and impact on host plant performance in the field

Species interactions are a unifying theme in ecology and evolution. Both fields are currently moving beyond their historical focus on isolated pairwise relationships to understand how ecological communities affect focal interactions. Additional species can modify both the number of interactions and the fitness consequences of each interaction (i.e., selection). Although only selection affects the evolution of the focal interaction, the two are often conflated, limiting our understanding of the evolution of multispecies interactions. We manipulated aboveground herbivory on the legume Medicago lupulina in the field and quantified its effect on number of symbionts and the per-symbiont impact on plant performance in two belowground symbioses: mutualistic rhizobia bacteria (Ensifer meliloti) and parasitic root-knot nematodes (Meloidogyne hapla). We found that herbivores modified the number of rhizobia nodules, as well as the benefit per nodule. However, each effect was specific to a distinct herbivory regime: natural herbivory affected nodule number, whereas leafhoppers (Cicadellidae) weakened the per nodule benefit. We did not detect any effect of herbivory on nematode gall number or the cost of infection. Our data demonstrate that distinguishing between symbiont number from the fitness consequences of symbiosis is crucial to accurately infer how pairwise interactions will evolve in a community.

Plant‐associate interactions and diversification across trophic levels

Interactions between species are widely understood to have promoted the diversification of life on Earth, but how interactions spur the formation of new species remains unclear. Interacting species often become locally adapted to each other, but they may also be subject to shared dispersal limitations and environmental conditions. Moreover, theory predicts that different kinds of interactions have different effects on diversification. To better understand how species interactions promote diversification, we compiled population genetic studies of host plants and intimately associated herbivores, parasites, and mutualists. We used Bayesian multiple regressions and the BEDASSLE modeling framework to test whether host and associate population structures were correlated over and above the potentially confounding effects of geography and shared environmental variation. We found that associates' population structure often paralleled their hosts' population structure, and that this effect is robust to accounting for geographic distance and climate. Associate genetic structure was significantly explained by plant genetic structure somewhat more often in antagonistic interactions than in mutualistic ones. This aligns with a key prediction of coevolutionary theory that antagonistic interactions promote diversity through local adaptation of antagonists to hosts, while mutualistic interactions more often promote diversity via the effect of hosts' geographic distribution on mutualists' dispersal.

Unraveling coevolutionary dynamics using ecological genomics

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

The soil microbiome increases plant survival and modifies interactions with root endosymbionts in the field

Evidence is accumulating that the soil microbiome-the community of microorganisms living in soils-has a major effect on plant traits and fitness. However, most work to date has taken place under controlled laboratory conditions and has not experimentally disentangled the effect of the soil microbiome on plant performance from the effects of key endosymbiotic constituents. As a result, it is difficult to extrapolate from existing data to understand the role of the soil microbiome in natural plant populations. To address this gap, we performed a field experiment using the black medick Medicago lupulina to test how the soil microbiome influences plant performance and colonization by two root endosymbionts (the mutualistic nitrogen-fixing bacteria Ensifer spp. and the parasitic root-knot nematode Meloidogyne hapla) under natural conditions. We inoculated all plants with nitrogen-fixing bacteria and factorially manipulated the soil microbiome and nematode infection. We found that plants grown in microbe-depleted soil exhibit greater mortality, but that among the survivors, there was no effect of the soil microbiome on plant performance (shoot biomass, root biomass, or shoot-to-root ratio). The soil microbiome also impacted parasitic nematode infection and affected colonization by mutualistic nitrogen-fixing bacteria in a plant genotype-dependent manner, increasing colonization in some plant genotypes and decreasing it in others. Our results demonstrate the soil microbiome has complex effects on plant-endosymbiont interactions and may be critical for survival under natural conditions.

Priority effects alter interaction outcomes in a legume-rhizobium mutualism

Priority effects occur when the order of species arrival affects the final community structure. Mutualists often interact with multiple partners in different orders, but if or how priority effects alter interaction outcomes is an open question. In the field, we paired the legume Medicago lupulina with two nodulating strains of Ensifer bacteria that vary in nitrogen-fixing ability. We inoculated plants with strains in different orders and measured interaction outcomes. The first strain to arrive primarily determined plant performance and final relative abundances of rhizobia on roots. Plants that received effective microbes first and ineffective microbes second grew larger than plants inoculated with the same microbes in the opposite order. Our results show that mutualism outcomes can be influenced not just by partner identity, but by the interaction order. Furthermore, hosts receiving high-quality mutualists early can better tolerate low-quality symbionts later, indicating that priority effects may help explain the persistence of ineffective symbionts.

Decreased coevolutionary potential and increased symbiont fecundity during the biological invasion of a legume-rhizobium mutualism

Although most invasive species engage in mutualism, we know little about how mutualism evolves as partners colonize novel environments. Selection on cooperation and standing genetic variation for mutualism traits may differ between a mutualism's invaded and native ranges, which could alter cooperation and coevolutionary dynamics. To test for such differences, we compare mutualism traits between invaded- and native-range host-symbiont genotype combinations of the weedy legume, Medicago polymorpha, and its nitrogen-fixing rhizobium symbiont, Ensifer medicae, which have coinvaded North America. We find that mutualism benefits for plants are indistinguishable between invaded- and native-range symbioses. However, rhizobia gain greater fitness from invaded-range mutualisms than from native-range mutualisms, and this enhancement of symbiont fecundity could increase the mutualism's spread by increasing symbiont availability during plant colonization. Furthermore, mutualism traits in invaded-range symbioses show lower genetic variance and a simpler partitioning of genetic variance between host and symbiont sources, compared to native-range symbioses. This suggests that biological invasion has reduced mutualists' potential to respond to coevolutionary selection. Additionally, rhizobia bearing a locus (hrrP) that can enhance symbiotic fitness have more exploitative phenotypes in invaded-range than in native-range symbioses. These findings highlight the impacts of biological invasion on the evolution of mutualistic interactions.

A need to consider the evolutionary genetics of host-symbiont mutualisms

Despite the ubiquity and importance of mutualistic interactions, we know little about the evolutionary genetics underlying their long-term persistence. As in antagonistic interactions, mutualistic symbioses are characterized by substantial levels of phenotypic and genetic diversity. In contrast to antagonistic interactions, however, we, by and large, do not understand how this variation arises, how it is maintained, nor its implications for future evolutionary change. Currently, we rely on phenotypic models to address the persistence of mutualistic symbioses, but the success of an interaction almost certainly depends heavily on genetic interactions. In this review, we argue that evolutionary genetic models could provide a framework for understanding the causes and consequences of diversity and why selection may favour processes that maintain variation in mutualistic interactions.

Soil origin and plant genotype structure distinct microbiome compartments in the model legume Medicago truncatula

BACKGROUND

Understanding the genetic and environmental factors that structure plant microbiomes is necessary for leveraging these interactions to address critical needs in agriculture, conservation, and sustainability. Legumes, which form root nodule symbioses with nitrogen-fixing rhizobia, have served as model plants for understanding the genetics and evolution of beneficial plant-microbe interactions for decades, and thus have added value as models of plant-microbiome interactions. Here we use a common garden experiment with 16S rRNA gene amplicon and shotgun metagenomic sequencing to study the drivers of microbiome diversity and composition in three genotypes of the model legume Medicago truncatula grown in two native soil communities.

RESULTS

Bacterial diversity decreased between external (rhizosphere) and internal plant compartments (root endosphere, nodule endosphere, and leaf endosphere). Community composition was shaped by strong compartment × soil origin and compartment × plant genotype interactions, driven by significant soil origin effects in the rhizosphere and significant plant genotype effects in the root endosphere. Nevertheless, all compartments were dominated by Ensifer, the genus of rhizobia that forms root nodule symbiosis with M. truncatula, and additional shotgun metagenomic sequencing suggests that the nodulating Ensifer were not genetically distinguishable from those elsewhere in the plant. We also identify a handful of OTUs that are common in nodule tissues, which are likely colonized from the root endosphere.

CONCLUSIONS

Our results demonstrate strong host filtering effects, with rhizospheres driven by soil origin and internal plant compartments driven by host genetics, and identify several key nodule-inhabiting taxa that coexist with rhizobia in the native range. Our results set the stage for future functional genetic experiments aimed at expanding our pairwise understanding of legume-rhizobium symbiosis toward a more mechanistic understanding of plant microbiomes. Video Abstract.

Environmental variation impacts trait expression and selection in the legume-rhizobium symbiosis

PREMISE

The ecological outcomes of mutualism are well known to shift across abiotic or biotic environments, but few studies have addressed how different environments impact evolutionary responses, including the intensity of selection on and the expression of genetic variance in key mutualism-related traits.

METHODS

We planted 30 maternal lines of the legume Medicago lupulina in four field common gardens and compared our measures of selection on and genetic variance in nodulation, a key trait reflecting legume investment in the symbiosis, with those from a previous greenhouse experiment using the same 30 M. lupulina lines.

RESULTS

We found that both the mean and genetic variance for nodulation were much greater in the greenhouse than in the field and that the form of selection on nodulation significantly differed across environments. We also found significant genotype-by-environment (G × E) effects for fitness-related traits that were generated by differences in the rank order of plant lines among environments.

CONCLUSIONS

Overall, our results suggest that the expression of genotypic variation and selection on nodulation differ across environments. In the field, significant rank-order changes for plant fitness potentially help maintain genetic variation in natural populations, even in the face of directional or stabilizing selection.

Exploring Genetic Diversity and Signatures of Horizontal Gene Transfer in Nodule Bacteria Associated with Lotus japonicus in Natural Environments

To investigate the genetic diversity and understand the process of horizontal gene transfer (HGT) in nodule bacteria associated with Lotus japonicus, we analyzed sequences of three housekeeping and five symbiotic genes using samples from a geographically wide range in Japan. A phylogenetic analysis of the housekeeping genes indicated that L. japonicus in natural environments was associated with diverse lineages of Mesorhizobium spp., whereas the sequences of symbiotic genes were highly similar between strains, resulting in remarkably low nucleotide diversity at both synonymous and nonsynonymous sites. Guanine-cytosine content values were lower in symbiotic genes, and relative frequencies of recombination between symbiotic genes were also lower than those between housekeeping genes. An analysis of molecular variance showed significant genetic differentiation among populations in both symbiotic and housekeeping genes. These results confirm that the Mesorhizobium genes required for symbiosis with L. japonicus behave as a genomic island (i.e., a symbiosis island) and suggest that this island has spread into diverse genomic backgrounds of Mesorhizobium via HGT events in natural environments. Furthermore, our data compilation revealed that the genetic diversity of symbiotic genes in L. japonicus-associated symbionts was among the lowest compared with reports of other species, which may be related to the recent population expansion proposed in Japanese populations of L. japonicus.

Genetic divergence among Bradyrhizobium strains nodulating wild and cultivated Kummerowia spp. in China

Distribution of rhizobial species is affected by geographical isolation and selected by leguminous hosts, however, little is known about the molecular evolution of rhizobia nodulating the same legume in different eco-environments. In present study, the microevolution of Bradyrhizobium associated with the leguminous grass Kummerowia grown in exurban areas and cultivated in urban areas in China was investigated. Total 14 genospecies, including seven new groups, were identified based on a concatenated sequence analysis of taxonomic markers (SMc00019, truA and thrA) for 94 representative strains. Results demonstrated that lower levels of nucleotide diversity were found in the strains isolated from urban areas compared with those isolated from exurban areas, based on the evolutional analyses of three housekeeping genes (atpD, glnII and recA), two symbiosis-related genes (nodC and nifH), and the taxonomic markers. Moreover, compared with urban areas, gene exchange and recombination occurred more frequently among the genospecies isolated from exurban areas, regardless of the geographical distribution. Finally, the evolutionary lineage of Bradyrhizobium strains isolated from urban areas was independent of that of the strains isolated from exurban areas. In summary, the evolutionary history of Kummerowia bradyrhizobia may have been gradually segregated to different evolutionary lineages, irrespective of distinct biogeography.

Dynamic genomic architecture of mutualistic cooperation in a wild population of Mesorhizobium

Research on mutualism seeks to explain how cooperation can be maintained when uncooperative mutants co-occur with cooperative kin. Gains and losses of the gene modules required for cooperation punctuate symbiont phylogenies and drive lifestyle transitions between cooperative symbionts and uncooperative free-living lineages over evolutionary time. Yet whether uncooperative symbionts commonly evolve from within cooperative symbiont populations or from within distantly related lineages with antagonistic or free-living lifestyles (i.e., third-party mutualism exploiters or parasites), remains controversial. We use genomic data to show that genotypes that differ in the presence or absence of large islands of symbiosis genes are common within a single wild recombining population of Mesorhizobium symbionts isolated from host tissues and are an important source of standing heritable variation in cooperation in this population. In a focal population of Mesorhizobium, uncooperative variants that lack a symbiosis island segregate at 16% frequency in nodules, and genome size and symbiosis gene number are positively correlated with cooperation. This finding contrasts with the genomic architecture of variation in cooperation in other symbiont populations isolated from host tissues in which the islands of genes underlying cooperation are ubiquitous and variation in cooperation is primarily driven by allelic substitution and individual gene gain and loss events. Our study demonstrates that uncooperative mutants within mutualist populations can comprise a significant component of genetic variation in nature, providing biological rationale for models and experiments that seek to explain the maintenance of mutualism in the face of non-cooperators.

Genetic conflict with a parasitic nematode disrupts the legume-rhizobia mutualism

Genetic variation for partner quality in mutualisms is an evolutionary paradox. One possible resolution to this puzzle is that there is a tradeoff between partner quality and other fitness‐related traits. Here, we tested whether susceptibility to parasitism is one such tradeoff in the mutualism between legumes and nitrogen‐fixing bacteria (rhizobia). We performed two greenhouse experiments with the legume Medicago truncatula. In the first, we inoculated each plant with the rhizobia Ensifer meliloti and with one of 40 genotypes of the parasitic root‐knot nematode Meloidogyne hapla. In the second experiment, we inoculated all plants with rhizobia and half of the plants with a genetically variable population of nematodes. Using the number of nematode galls as a proxy for infection severity, we found that plant genotypes differed in susceptibility to nematode infection, and nematode genotypes differed in infectivity. Second, we showed that there was a genetic correlation between the number of mutualistic structures formed by rhizobia (nodules) and the number of parasitic structures formed by nematodes (galls). Finally, we found that nematodes disrupt the rhizobia mutualism: nematode‐infected plants formed fewer nodules and had less nodule biomass than uninfected plants. Our results demonstrate that there is genetic conflict between attracting rhizobia and repelling nematodes in Medicago. If genetic conflict with parasitism is a general feature of mutualism, it could account for the maintenance of genetic variation in partner quality and influence the evolutionary dynamics of positive species interactions.

No evidence for adaptation to local rhizobial mutualists in the legume Medicago lupulina

Local adaptation is a common but not ubiquitous feature of species interactions, and understanding the circumstances under which it evolves illuminates the factors that influence adaptive population divergence. Antagonistic species interactions dominate the local adaptation literature relative to mutualistic ones, preventing an overall assessment of adaptation within interspecific interactions. Here, we tested whether the legume Medicago lupulina is adapted to the locally abundant species of mutualistic nitrogen-fixing rhizobial bacteria that vary in frequency across its eastern North American range. We reciprocally inoculated northern and southern M. lupulina genotypes with the northern (Ensifer medicae) or southern bacterium (E. meliloti) in a greenhouse experiment. Despite producing different numbers of root nodules (the structures in which the plants house the bacteria), neither northern nor southern plants produced more seeds, flowered earlier, or were more likely to flower when inoculated with their local rhizobia. We then used a pre-existing dataset to perform a genome scan for loci that showed elevated differentiation between field-collected plants that hosted different bacteria. None of the loci we identified belonged to the well-characterized suite of legume-rhizobia symbiosis genes, suggesting that the rhizobia do not drive genetic divergence between M. lupulina populations. Our results demonstrate that symbiont local adaptation has not evolved in this mutualism despite large-scale geographic variation in the identity of the interacting species.