Intrapopulation genomics in a model mutualist: Population structure and candidate symbiosis genes under selection in Medicago truncatula

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.

Phenotypic and genomic signatures of interspecies cooperation and conflict in naturally occurring isolates of a model plant symbiont

Given the need to predict the outcomes of (co)evolution in host-associated microbiomes, whether microbial and host fitnesses tend to trade-off, generating conflict, remains a pressing question. Examining the relationships between host and microbe fitness proxies at both the phenotypic and genomic levels can illuminate the mechanisms underlying interspecies cooperation and conflict. We examined naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti, paired with each of two host Medicago truncatula genotypes in single- or multi-strain experiments to determine how multiple proxies of microbial and host fitness were related to one another and test key predictions about mutualism evolution at the genomic scale, while also addressing the challenge of measuring microbial fitness. We found little evidence for interspecies fitness conflict; loci tended to have concordant effects on both microbe and host fitnesses, even in environments with multiple co-occurring strains. Our results emphasize the importance of quantifying microbial relative fitness for understanding microbiome evolution and thus harnessing microbiomes to improve host fitness. Additionally, we find that mutualistic coevolution between hosts and microbes acts to maintain, rather than erode, genetic diversity, potentially explaining why variation in mutualism traits persists in nature.

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.

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.

Genome-Wide SNPs Provide Insights on the Cryptic Genetic Structure and Signatures of Climate Adaption in Amorphophallus albus Germplasms

Domesticated species represent unique systems in which the evolutionary genomic consequences of intensive selective breeding and adaptation can be thoroughly investigated. Amorphophallus albus occurs naturally and is in cultivation throughout the downstream region of the Jinshajiang River in Southwest China. This species is characterised by high konjac glucomannan content, and has been cultivated in China for nearly 2,000 years. To study genetic differentiation and local adaption of A. albus, we sampled 13 distinct local cultivated populations of this species. Restriction site-associated DNA sequencing was conducted with 87 samples, resulting in 24,225 SNPs. The population structure analyses suggest two main genetic groups: one in the relatively upstream region, and one downstream. We found evidence of additional sub-structure within the upstream group, demonstrating the statistical power of genomic SNPs in discovering subtle genetic structure. The environmental and geographic factors were all identified as significant in shaping the genetic differentiation of this species. Notably, the proportion of environmental factors was larger than geographic factors in influencing the population genetic patterns of A. albus. We also discovered loci that were associated with local adaptation. These findings will help us understand the genetic differentiation of this newly domesticated species, thereby informing future breeding programs of A. albus.

Distinct genetic basis for root responses to lipo-chitooligosaccharide signal molecules from different microbial origins

Lipo-chitooligosaccharides (LCOs) were originally found as symbiotic signals called Nod Factors (Nod-LCOs) controlling the nodulation of legumes by rhizobia. More recently, LCOs were also found in symbiotic fungi and, more surprisingly, very widely in the kingdom Fungi, including in saprophytic and pathogenic fungi. The LCO-V(C18:1, fucosylated/methyl fucosylated), hereafter called Fung-LCOs, are the LCO structures most commonly found in fungi. This raises the question of how legume plants such as Medicago truncatula can discriminate between Nod-LCOs and Fung-LCOs. To address this question, we performed a genome-wide association study on 173 natural accessions of M. truncatula, using a root branching phenotype and a newly developed local score approach. Both Nod-LCOs and Fung-LCOs stimulated root branching in most accessions, but the root responses to these two types of LCO molecules were not correlated. In addition, the heritability of the root response was higher for Nod-LCOs than for Fung-LCOs. We identified 123 loci for Nod-LCO and 71 for Fung-LCO responses, of which only one was common. This suggests that Nod-LCOs and Fung-LCOs both control root branching but use different molecular mechanisms. The tighter genetic constraint of the root response to Fung-LCOs possibly reflects the ancestral origin of the biological activity of these molecules.

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.

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.

Reference-based RADseq resolves robust relationships among closely related species of lichen-forming fungi using metagenomic DNA

Despite increasing availability of phylogenomic datasets, strategies to generate genome-scale data from organisms involved in symbiotic relationships remains challenging. Restriction site-associated DNA sequencing (RADseq) can effectively generated reduced representation genomic loci. However, when using metagenomic DNA from inseparable symbiotic organisms, RADseq loci may belong to any number of the organisms involved in these intimate associations. In this study, we explored the potential for a reference-based RADseq approach to generate data for lichen-forming fungi from metagenomic DNA extracted from intact lichens. We simulated RAD data from draft genomes of closely related lichenized fungi to test if RADseq can reconstruct robust evolutionary relationships. Subsequently, we generated empirical RADseq data from metagenomic lichen DNA, with RADseq loci mapped back to a reference genome to exclude loci from other lichen symbionts that are represented in metagenomic libraries. In all cases, phylogenetic reconstructions using RADseq loci recovered diversification histories consistent with a previous study based on more comprehensive genome sampling. Furthermore, RADseq loci were found to resolve relationships among closely related species, which were otherwise indistinguishable using a phylogenetic species recognition criterion. Our studies revealed that a modified, reference-based RADseq approach can successfully be implemented to generate symbiont-specific phylogenomic data from metagenomic reads.

Policing the legume-Rhizobium symbiosis: a critical test of partner choice

In legume-Rhizobium symbioses, specialised soil bacteria fix atmospheric nitrogen in return for carbon. However, ineffective strains can arise, making discrimination essential. Discrimination can occur via partner choice, where legumes prevent ineffective strains from entering, or via sanctioning, where plants provide fewer resources. Several studies have inferred that legumes exercise partner choice, but the rhizobia compared were not otherwise isogenic. To test when and how plants discriminate ineffective strains we developed sets of fixing and non-fixing strains that differed only in the expression of nifH - essential for nitrogen fixation - and could be visualised using marker genes. We show that the plant is unable to select against the non-fixing strain at the point of entry, but that non-fixing nodules are sanctioned. We also used the technique to characterise mixed nodules (containing both a fixing and a non-fixing strain), whose frequency could be predicted using a simple diffusion model. We discuss that sanctioning is likely to evolve in preference to partner choice in any symbiosis where partner quality cannot be adequately assessed until goods or services are actively exchanged.