The evolutionary genomics of adaptation to stress in wild rhizobium bacteria.

Microbiota comprise the bulk of life's diversity, yet we know little about how populations of microbes accumulate adaptive diversity across natural landscapes. Adaptation to stressful soil conditions in plants provides seminal examples of adaptation in response to natural selection via allelic substitution. For microbes symbiotic with plants however, horizontal gene transfer allows for adaptation via gene gain and loss, which could generate fundamentally different evolutionary dynamics. We use comparative genomics and genetics to elucidate the evolutionary mechanisms of adaptation to physiologically stressful serpentine soils in rhizobial bacteria in western North American grasslands. In vitro experiments demonstrate that the presence of a locus of major effect, the nre operon, is necessary and sufficient to confer adaptation to nickel, a heavy metal enriched to toxic levels in serpentine soil, and a major axis of environmental soil chemistry variation. We find discordance between inferred evolutionary histories of the core genome and nreAXY genes, which often reside in putative genomic islands. This suggests that the evolutionary history of this adaptive variant is marked by frequent losses, and/or gains via horizontal acquisition across divergent rhizobium clades. However, different nre alleles confer distinct levels of nickel resistance, suggesting allelic substitution could also play a role in rhizobium adaptation to serpentine soil. These results illustrate that the interplay between evolution via gene gain and loss and evolution via allelic substitution may underlie adaptation in wild soil microbiota. Both processes are important to consider for understanding adaptive diversity in microbes and improving stress-adapted microbial inocula for human use.

Topological data analysis reveals a core gene expression backbone that defines form and function across flowering plants.

Since they emerged approximately 125 million years ago, flowering plants have evolved to dominate the terrestrial landscape and survive in the most inhospitable environments on earth. At their core, these adaptations have been shaped by changes in numerous, interconnected pathways and genes that collectively give rise to emergent biological phenomena. Linking gene expression to morphological outcomes remains a grand challenge in biology, and new approaches are needed to begin to address this gap. Here, we implemented topological data analysis (TDA) to summarize the high dimensionality and noisiness of gene expression data using lens functions that delineate plant tissue and stress responses. Using this framework, we created a topological representation of the shape of gene expression across plant evolution, development, and environment for the phylogenetically diverse flowering plants. The TDA-based Mapper graphs form a well-defined gradient of tissues from leaves to seeds, or from healthy to stressed samples, depending on the lens function. This suggests that there are distinct and conserved expression patterns across angiosperms that delineate different tissue types or responses to biotic and abiotic stresses. Genes that correlate with the tissue lens function are enriched in central processes such as photosynthetic, growth and development, housekeeping, or stress responses. Together, our results highlight the power of TDA for analyzing complex biological data and reveal a core expression backbone that defines plant form and function.

Hosts winnow symbionts with multiple layers of absolute and conditional discrimination mechanisms.

In mutualism, hosts select symbionts via partner choice and preferentially direct more resources to symbionts that provide greater benefits via sanctions. At the initiation of symbiosis, prior to resource exchange, it is not known how the presence of multiple symbiont options (i.e. the symbiont social environment) impacts partner choice outcomes. Furthermore, little research addresses whether hosts primarily discriminate among symbionts via sanctions, partner choice or a combination. We inoculated the legume, Acmispon wrangelianus, with 28 pairs of fluorescently labelled Mesorhizobium strains that vary continuously in quality as nitrogen-fixing symbionts. We find that hosts exert robust partner choice, which enhances their fitness. This partner choice is conditional such that a strain's success in initiating nodules is impacted by other strains in the social environment. This social genetic effect is as important as a strain's own genotype in determining nodulation and has both transitive (consistent) and intransitive (idiosyncratic) effects on the probability that a symbiont will form a nodule. Furthermore, both absolute and conditional partner choice act in concert with sanctions, among and within nodules. Thus, multiple forms of host discrimination act as a series of sieves that optimize host benefits and select for costly symbiont cooperation in mixed symbiont populations.

Negotiating mutualism: A locus for exploitation by rhizobia has a broad effect size distribution and context-dependent effects on legume hosts.

In mutualisms, variation at genes determining partner fitness provides the raw material upon which coevolutionary selection acts, setting the dynamics and pace of coevolution. However, we know little about variation in the effects of genes that underlie symbiotic fitness in natural mutualist populations. In some species of legumes that form root nodule symbioses with nitrogen-fixing rhizobial bacteria, hosts secrete nodule-specific cysteine-rich (NCR) peptides that cause rhizobia to differentiate in the nodule environment. However, rhizobia can cleave NCR peptides through the expression of genes like the plasmid-borne Host range restriction peptidase (hrrP), whose product degrades specific NCR peptides. Although hrrP activity can confer host exploitation by depressing host fitness and enhancing symbiont fitness, the effects of hrrP on symbiosis phenotypes depend strongly on the genotypes of the interacting partners. However, the effects of hrrP have yet to be characterised in a natural population context, so its contribution to variation in wild mutualist populations is unknown. To understand the distribution of effects of hrrP in wild rhizobia, we measured mutualism phenotypes conferred by hrrP in 12 wild Ensifer medicae strains. To evaluate context dependency of hrrP effects, we compared hrrP effects across two Medicago polymorpha host genotypes and across two experimental years for five E. medicae strains. We show for the first time in a natural population context that hrrP has a wide distribution of effect sizes for many mutualism traits, ranging from strongly positive to strongly negative. Furthermore, we show that hrrP effect size varies across host genotypes and experiment years, suggesting that researchers should be cautious about extrapolating the role of genes in natural populations from controlled laboratory studies of single genetic variants.

Previously unmeasured genetic diversity explains part of Lewontin's paradox in a k-mer-based meta-analysis of 112 plant species.

At the molecular level, most evolution is expected to be neutral. A key prediction of this expectation is that the level of genetic diversity in a population should scale with population size. However, as was noted by Richard Lewontin in 1974 and reaffirmed by later studies, the slope of the population size-diversity relationship in nature is much weaker than expected under neutral theory. We hypothesize that one contributor to this paradox is that current methods relying on single nucleotide polymorphisms (SNPs) called from aligning short reads to a reference genome underestimate levels of genetic diversity in many species. To test this idea, we calculated nucleotide diversity ( π ) and k-mer-based metrics of genetic diversity across 112 plant species, amounting to over 205 terabases of DNA sequencing data from 27,488 individual plants. We then compared how these different metrics correlated with proxies of population size that account for both range size and population density variation across species. We found that our population size proxies scaled anywhere from about 3 to over 20 times faster with k-mer diversity than nucleotide diversity after adjusting for evolutionary history, mating system, life cycle habit, cultivation status, and invasiveness. The relationship between k-mer diversity and population size proxies also remains significant after correcting for genome size, whereas the analogous relationship for nucleotide diversity does not. These results suggest that variation not captured by common SNP-based analyses explains part of Lewontin's paradox in plants.

Weaker selection on genes with treatment-specific expression consistent with a limit on plasticity evolution in Arabidopsis thaliana.

Differential gene expression between environments often underlies phenotypic plasticity. However, environment-specific expression patterns are hypothesized to relax selection on genes, and thus limit plasticity evolution. We collated over 27 terabases of RNA-sequencing data on Arabidopsis thaliana from over 300 peer-reviewed studies and 200 treatment conditions to investigate this hypothesis. Consistent with relaxed selection, genes with more treatment-specific expression have higher levels of nucleotide diversity and divergence at nonsynonymous sites but lack stronger signals of positive selection. This result persisted even after controlling for expression level, gene length, GC content, the tissue specificity of expression, and technical variation between studies. Overall, our investigation supports the existence of a hypothesized trade-off between the environment specificity of a gene's expression and the strength of selection on said gene in A. thaliana. Future studies should leverage multiple genome-scale datasets to tease apart the contributions of many variables in limiting plasticity evolution.

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.