Is there a latitudinal diversity gradient for symbiotic microbes? A case study with sensitive partridge peas.

Mutualism is thought to be more prevalent in the tropics than temperate zones and may therefore play an important role in generating and maintaining high species richness found at lower latitudes. However, results on the impact of mutualism on latitudinal diversity gradients are mixed, and few empirical studies sample both temperate and tropical regions. We investigated whether a latitudinal diversity gradient exists in the symbiotic microbial community associated with the legume Chamaecrista nictitans. We sampled bacteria DNA from nodules and the surrounding soil of plant roots across a latitudinal gradient (38.64-8.68 °N). Using 16S rRNA sequence data, we identified many non-rhizobial species within C. nictitans nodules that cannot form nodules or fix nitrogen. Species richness increased towards lower latitudes in the non-rhizobial portion of the nodule community but not in the rhizobial community. The microbe community in the soil did not effectively predict the non-rhizobia community inside nodules, indicating that host selection is important for structuring non-rhizobia communities in nodules. We next factorially manipulated the presence of three non-rhizobia strains in greenhouse experiments and found that co-inoculations of non-rhizobia strains with rhizobia had a marginal effect on nodule number and no effect on plant growth. Our results suggest that these non-rhizobia bacteria are likely commensals-species that benefit from associating with a host but are neutral for host fitness. Overall, our study suggests that temperate C. nictitans plants are more selective in their associations with the non-rhizobia community, potentially due to differences in soil nitrogen across latitude.

Quantifying the role of genome size and repeat content in adaptive variation and the architecture of flowering time in Amaranthus tuberculatus.

Genome size variation, largely driven by repeat content, is poorly understood within and among populations, limiting our understanding of its significance for adaptation. Here we characterize intraspecific variation in genome size and repeat content across 186 individuals of Amaranthus tuberculatus, a ubiquitous native weed that shows flowering time adaptation to climate across its range and in response to agriculture. Sequence-based genome size estimates vary by up to 20% across individuals, consistent with the considerable variability in the abundance of transposable elements, unknown repeats, and rDNAs across individuals. The additive effect of this variation has important phenotypic consequences-individuals with more repeats, and thus larger genomes, show slower flowering times and growth rates. However, compared to newly-characterized gene copy number and polygenic nucleotide changes underlying variation in flowering time, we show that genome size is a marginal contributor. Differences in flowering time are reflected by genome size variation across sexes and marginally, habitats, while polygenic variation and a gene copy number variant within the ATP synthesis pathway show consistently stronger environmental clines than genome size. Repeat content nonetheless shows non-neutral distributions across the genome, and across latitudinal and environmental gradients, demonstrating the numerous governing processes that in turn influence quantitative genetic variation for phenotypes key to plant adaptation.

Strong selection is poorly aligned with genetic variation in Ipomoea hederacea: implications for divergence and constraint.

The multivariate evolution of populations is the result of the interactions between natural selection, drift, and the underlying genetic structure of the traits involved. Covariances among traits bias responses to selection, and the multivariate axis which describes the greatest genetic variation is expected to be aligned with patterns of divergence across populations. An exception to this expectation is when selection acts on trait combinations lacking genetic variance, which limits evolutionary change. Here we used a common garden field experiment of individuals from 57 populations of Ipomoea hederacea to characterize linear and nonlinear selection on 5 quantitative traits in the field. We then formally compare patterns of selection to previous estimates of within population genetic covariance structure (the G-matrix) and population divergence in these traits. We found that selection is poorly aligned with previous estimates of genetic covariance structure and population divergence. In addition, the trait combinations favored by selection were generally lacking genetic variation, possessing approximately 15%-30% as much genetic variation as the most variable combination of traits. Our results suggest that patterns of population divergence are likely the result of the interplay between adaptive responses, correlated responses, and selection favoring traits lacking genetic variation.

G-matrix stability in clinally diverging populations of an annual weed.

How phenotypic and genetic divergence among populations is influenced by the genetic architecture of those traits, and how microevolutionary changes in turn affect the within-population patterns of genetic variation, are of major interest to evolutionary biology. Work on Ipomoea hederacea, an annual vine, has found genetic clines in the means of a suite of ecologically important traits, including flowering time, growth rate, seed mass, and corolla width. Here we investigate the genetic (co)variances of these clinally varying traits in two northern range-edge and two central populations of I. hederacea to evaluate the influence of the genetic architecture on divergence across the range. We find (1) limited evidence for clear differentiation between Northern and Southern populations in the structure of G, suggesting overall stability of G across the range despite mean trait divergence and (2) that the axes of greatest variation (gmax) were unaligned with the axis of greatest multivariate divergence. Together these results indicate the role of the quantitative genetic architecture in constraining evolutionary response and divergence among populations across the geographic range.

Nonsymbiotic legumes are more invasive, but only if polyploid.

Both mutualism and polyploidy are thought to influence invasion success in plants, but few studies have tested their joint effects. Mutualism can limit range expansion when plants cannot find a compatible partner in a novel habitat, or facilitate range expansion when mutualism increases a plant's niche breadth. Polyploids are also expected to have greater niche breadth because of greater self-compatibility and phenotypic plasticity, increasing invasion success. For 847 legume species, we compiled data from published sources to estimate ploidy, symbiotic status with rhizobia, specificity on rhizobia, and the number of introduced ranges. We found that diploid species have had limited spread around the globe regardless of whether they are symbiotic or how many rhizobia partners they can host. Polyploids, by contrast, have been successfully introduced to many new ranges, but interactions with rhizobia constrain their range expansion. In a hidden state model of trait evolution, we also found evidence of a high rate of re-diploidization in symbiotic legume lineages, suggesting that symbiosis and ploidy may interact at macroevolutionary scales. Overall, our results suggest that symbiosis with rhizobia limits range expansion when legumes are polyploid but not diploid.

Rapid weed adaptation and range expansion in response to agriculture over the past two centuries.

North America has experienced a massive increase in cropland use since 1800, accompanied more recently by the intensification of agricultural practices. Through genome analysis of present-day and historical samples spanning environments over the past two centuries, we studied the effect of these changes in farming on the extent and tempo of evolution across the native range of the common waterhemp (Amaranthus tuberculatus), a now pervasive agricultural weed. Modern agriculture has imposed strengths of selection rarely observed in the wild, with notable shifts in allele frequency trajectories since agricultural intensification in the 1960s. An evolutionary response to this extreme selection was facilitated by a concurrent human-mediated range shift. By reshaping genome-wide diversity across the landscape, agriculture has driven the success of this weed in the 21st century.

The relative role of plasticity and demographic history in Capsella bursa-pastoris: a common garden experiment in Asia and Europe.

The colonization success of a species depends on the interplay between its phenotypic plasticity, adaptive potential and demographic history. Assessing their relative contributions during the different phases of a species range expansion is challenging, and requires large-scale experiments. Here, we investigated the relative contributions of plasticity, performance and demographic history to the worldwide expansion of the shepherd's purse, Capsella bursa-pastoris. We installed two large common gardens of the shepherd's purse, a young, self-fertilizing, allopolyploid weed with a worldwide distribution. One common garden was located in Europe, the other in Asia. We used accessions from three distinct genetic clusters (Middle East, Europe and Asia) that reflect the demographic history of the species. Several life-history traits were measured. To explain the phenotypic variation between and within genetic clusters, we analysed the effects of (i) the genetic clusters, (ii) the phenotypic plasticity and its association to fitness and (iii) the distance in terms of bioclimatic variables between the sampling site of an accession and the common garden, i.e. the environmental distance. Our experiment showed that (i) the performance of C. bursa-pastoris is closely related to its high phenotypic plasticity; (ii) within a common garden, genetic cluster was a main determinant of phenotypic differences; and (iii) at the scale of the experiment, the effect of environmental distance to the common garden could not be distinguished from that of genetic clusters. Phenotypic plasticity and demographic history both play important role at different stages of range expansion. The success of the worldwide expansion of C. bursa-pastoris was undoubtedly influenced by its strong phenotypic plasticity.

Repeated origins, widespread gene flow, and allelic interactions of target-site herbicide resistance mutations.

Causal mutations and their frequency in agricultural fields are well-characterized for herbicide resistance. However, we still lack understanding of their evolutionary history: the extent of parallelism in the origins of target-site resistance (TSR), how long these mutations persist, how quickly they spread, and allelic interactions that mediate their selective advantage. We addressed these questions with genomic data from 19 agricultural populations of common waterhemp (Amaranthus tuberculatus), which we show to have undergone a massive expansion over the past century, with a contemporary effective population size estimate of 8 x 107. We found variation at seven characterized TSR loci, two of which had multiple amino acid substitutions, and three of which were common. These three common resistance variants show extreme parallelism in their mutational origins, with gene flow having shaped their distribution across the landscape. Allele age estimates supported a strong role of adaptation from de novo mutations, with a median age of 30 suggesting that most resistance alleles arose soon after the onset of herbicide use. However, resistant lineages varied in both their age and evidence for selection over two different timescales, implying considerable heterogeneity in the forces that govern their persistence. Two such forces are intra- and inter-locus allelic interactions; we report a signal of extended haplotype competition between two common TSR alleles, and extreme linkage with genome-wide alleles with known functions in resistance adaptation. Together, this work reveals a remarkable example of spatial parallel evolution in a metapopulation, with important implications for the management of herbicide resistance.

Selective ancestral sorting and de novo evolution in the agricultural invasion of Amaranthus tuberculatus.

The relative role of hybridization, de novo evolution, and standing variation in weed adaptation to agricultural environments is largely unknown. In Amaranthus tuberculatus, a widespread North American agricultural weed, adaptation is likely influenced by recent secondary contact and admixture of two previously isolated lineages. We characterized the extent of adaptation and phenotypic differentiation accompanying the spread of A. tuberculatus into agricultural environments and the contribution of ancestral divergence. We generated phenotypic and whole-genome sequence data from a manipulative common garden experiment, using paired samples from natural and agricultural populations. We found strong latitudinal, longitudinal, and sex differentiation in phenotypes, and subtle differences among agricultural and natural environments that were further resolved with ancestry inference. The transition into agricultural environments has favored southwestern var. rudis ancestry that leads to higher biomass and treatment-specific phenotypes: increased biomass and earlier flowering under reduced water availability, and reduced plasticity in fitness-related traits. We also detected de novo adaptation in individuals from agricultural habitats independent of ancestry effects, including marginally higher biomass, later flowering, and treatment-dependent divergence in time to germination. Therefore, the invasion of A. tuberculatus into agricultural environments has drawn on adaptive variation across multiple timescales-through both preadaptation via the preferential sorting of var. rudis ancestry and de novo local adaptation.

Selection on Accessible Chromatin Regions in Capsella grandiflora.

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.

Multiple Mutualism Effects generate synergistic selection and strengthen fitness alignment in the interaction between legumes, rhizobia and mycorrhizal fungi.

Nearly all organisms participate in multiple mutualisms, and complementarity within these complex interactions can result in synergistic fitness effects. However, it remains largely untested how multiple mutualisms impact eco-evolutionary dynamics in interacting species. We tested how multiple microbial mutualists-N-fixing bacteria and mycorrrhizal fungi-affected selection and heritability of traits in their shared host plant (Medicago truncatula), as well as fitness alignment between partners. Our results demonstrate for the first time that multiple mutualisms synergistically affect the selection and heritability of host traits and enhance fitness alignment between mutualists. Specifically, we found interaction with multiple microbial symbionts doubled the strength of natural selection on a plant architectural trait, resulted in 2- to 3-fold higher heritability of plant reproductive success, and more than doubled fitness alignment between N-fixing bacteria and plants. These findings show synergism generated by multiple mutualisms extends to key components of microevolutionary change, emphasising the importance of multiple mutualism effects on evolutionary trajectories.

Introduced populations of ragweed show as much evolutionary potential as native populations.

Invasive species are a global economic and ecological problem. They also offer an opportunity to understand evolutionary processes in a colonizing context. The impacts of evolutionary factors, such as genetic variation, on the invasion process are increasingly appreciated, but there remain gaps in the empirical literature. The adaptive potential of populations can be quantified using genetic variance-covariance matrices (G), which encapsulate the heritable genetic variance in a population. Here, we use a multivariate Bayesian approach to assess the adaptive potential of invasive populations of ragweed (Ambrosia artemisiifolia), a serious allergen and agricultural weed. We compared several aspects of genetic architecture and the structure of G matrices between three native and three introduced populations, based on phenotypic data collected in a field common garden experiment. We found moderate differences in the quantitative genetic architecture among populations, but we did not find that introduced populations suffer from a limited adaptive potential or increased genetic constraint compared with native populations. Ragweed has an annual life history, is an obligate outcrosser, and produces very large numbers of seeds and pollen grains. These characteristics, combined with the significant additive genetic variance documented here, suggest ragweed will be able to respond quickly to selection pressures in both its native and introduced ranges.

The genetic architecture and population genomic signatures of glyphosate resistance in Amaranthus tuberculatus.

Much of what we know about the genetic basis of herbicide resistance has come from detailed investigations of monogenic adaptation at known target-sites, despite the increasingly recognized importance of polygenic resistance. Little work has been done to characterize the broader genomic basis of herbicide resistance, including the number and distribution of genes involved, their effect sizes, allele frequencies and signatures of selection. In this work, we implemented genome-wide association (GWA) and population genomic approaches to examine the genetic architecture of glyphosate (Round-up) resistance in the problematic agricultural weed Amaranthus tuberculatus. A GWA was able to correctly identify the known target-gene but statistically controlling for two causal target-site mechanisms revealed an additional 250 genes across all 16 chromosomes associated with non-target-site resistance (NTSR). The encoded proteins had functions that have been linked to NTSR, the most significant of which is response to chemicals, but also showed pleiotropic roles in reproduction and growth. Compared to an empirical null that accounts for complex population structure, the architecture of NTSR was enriched for large effect sizes and low allele frequencies, suggesting the role of pleiotropic constraints on its evolution. The enrichment of rare alleles also suggested that the genetic architecture of NTSR may be population-specific and heterogeneous across the range. Despite their rarity, we found signals of recent positive selection on NTSR-alleles by both window- and haplotype-based statistics, and an enrichment of amino acid changing variants. In our samples, genome-wide single nucleotide polymorphisms explain a comparable amount of the total variation in glyphosate resistance to monogenic mechanisms, even in a collection of individuals where 80% of resistant individuals have large-effect TSR mutations, indicating an underappreciated polygenic contribution to the evolution of herbicide resistance in weed populations.

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.

The Potential for Genotype-by-Environment Interactions to Maintain Genetic Variation in a Model Legume-Rhizobia Mutualism.

The maintenance of genetic variation in mutualism-related traits is key for understanding mutualism evolution, yet the mechanisms maintaining variation remain unclear. We asked whether genotype-by-environment (G×E) interaction is a potential mechanism maintaining variation in the model legume-rhizobia system, Medicago truncatula-Ensifer meliloti. We planted 50 legume genotypes in a greenhouse under ambient light and shade to reflect reduced carbon availability for plants. We found an expected reduction under shaded conditions for plant performance traits, such as leaf number, aboveground and belowground biomass, and a mutualism-related trait, nodule number. We also found G×E for nodule number, with ∼83% of this interaction due to shifts in genotype fitness rank order across light environments, coupled with strong positive directional selection on nodule number regardless of light environment. Our results suggest that G×E can maintain genetic variation in a mutualism-related trait that is under consistent positive directional selection across light environments.

Population genomics of parallel adaptation.

Parallel evolution is one of the striking patterns in nature. The presence of repeated evolution of the same phenotypes, suites of traits, and adaptations suggests a strong role for natural selection in shaping biological diversity. The reasoning is straightforward: each instance of repeated evolution makes it less likely that these features evolved neutrally or due to stochastic forces in each population or species. With the growing sequencing capability, we are now poised to examine the genetic basis of parallel evolution in model and nonmodel systems. On pages 4102-4117 of this issue of Molecular Ecology, van Boheemen and Hodgins (2020) provide an exemplar study of this kind, using common ragweed (Ambrosia artemisiifolia; Figure 1a). Their study is noteworthy and ambitious in many respects, and we think will serve as a model for studying parallel adaptation, even in nonmodel species.

Parallel flowering time clines in native and introduced ragweed populations are likely due to adaptation.

As introduced species expand their ranges, they often encounter differences in climate which are often correlated with geography. For introduced species, encountering a geographically variable climate sometimes leads to the re-establishment of clines seen in the native range. However, clines can also be caused by neutral processes, and so it is important to gather additional evidence that population differentiation is the result of selection as opposed to nonadaptive processes. Here, we examine phenotypic and genetic differences in ragweed from the native (North America) and introduced (European) ranges. We used a common garden to assess phenotypic differentiation in size and flowering time in ragweed populations. We found significant parallel clines in flowering time in both North America and Europe. Height and branch number had significant clines in North America, and, while not statistically significant, the patterns in Europe were the same. We used SNP data to assess population structure in both ranges and to compare phenotypic differentiation to neutral genetic variation. We failed to detect significant patterns of isolation by distance, geographic patterns in population structure, or correlations between the major axes of SNP variation and phenotypes or latitude of origin. We conclude that the North American clines in size and the parallel clines seen for flowering time are most likely the result of adaptation.

The Evolutionary Forces Shaping Cis- and Trans-Regulation of Gene Expression within a Population of Outcrossing Plants.

Understanding the persistence of genetic variation within populations has long been a goal of evolutionary biology. One promising route toward achieving this goal is using population genetic approaches to describe how selection acts on the loci associated with trait variation. Gene expression provides a model trait for addressing the challenge of the maintenance of variation because it can be measured genome-wide without information about how gene expression affects traits. Previous work has shown that loci affecting the expression of nearby genes (local or cis-eQTLs) are under negative selection, but we lack a clear understanding of the selective forces acting on variants that affect the expression of genes in trans. Here, we identify loci that affect gene expression in trans using genomic and transcriptomic data from one population of the obligately outcrossing plant, Capsella grandiflora. The allele frequencies of trans-eQTLs are consistent with stronger negative selection acting on trans-eQTLs than cis-eQTLs, and stronger negative selection acting on trans-eQTLs associated with the expression of multiple genes. However, despite this general pattern, we still observe the presence of a trans-eQTL at intermediate frequency that affects the expression of a large number of genes in the same coexpression module. Overall, our work highlights the different selective pressures shaping variation in cis- and trans-regulation.

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.