Selective trade-offs maintain alleles underpinning complex trait variation in plants

Fluctuating reproductive isolation and stable ancestry structure in a fine-scaled mosaic of hybridizing Mimulus monkeyflowers

Hybridization among taxa impacts a variety of evolutionary processes from adaptation to extinction. We seek to understand both patterns of hybridization across taxa and the evolutionary and ecological forces driving those patterns. To this end, we use whole-genome low-coverage sequencing of 458 wild-grown and 1565 offspring individuals to characterize the structure, stability, and mating dynamics of admixed populations of Mimulus guttatus and Mimulus nasutus across a decade of sampling. In three streams, admixed genomes are common and a M. nasutus organellar haplotype is fixed in M. guttatus, but new hybridization events are rare. Admixture is strongly unidirectional, but each stream has a unique distribution of ancestry proportions. In one stream, three distinct cohorts of admixed ancestry are spatially structured at ~20-50m resolution and stable across years. Mating system provides almost complete isolation of M. nasutus from both M. guttatus and admixed cohorts, and is a partial barrier between admixed and M. guttatus cohorts. Isolation due to phenology is near-complete between M. guttatus and M. nasutus. Phenological isolation is a strong barrier in some years between admixed and M. guttatus cohorts, but a much weaker barrier in other years, providing a potential bridge for gene flow. These fluctuations are associated with differences in water availability across years, supporting a role for climate in mediating the strength of reproductive isolation. Together, mating system and phenology accurately predict fluctuations in assortative mating across years, which we estimate directly using paired maternal and offspring genotypes. Climate-driven fluctuations in reproductive isolation may promote the longer-term stability of a complex mosaic of hybrid ancestry, preventing either complete isolation or complete collapse of species barriers.

The genetic architecture of floral trait divergence between hummingbird- and self-pollinated monkeyflower (Mimulus) species

Pollination syndromes are a key component of flowering plant diversification, prompting questions about the architecture of single traits and genetic coordination among traits. Here, we investigate the genetics of extreme floral divergence between naturally hybridizing monkeyflowers, Mimulus parishii (self-pollinated) and M. cardinalis (hummingbird-pollinated). We mapped quantitative trait loci (QTLs) for 18 pigment, pollinator reward/handling, and dimensional traits in parallel sets of F2 hybrids plus recombinant inbred lines and generated nearly isogenic lines (NILs) for two dimensional traits, pistil length and corolla size. Our multi-population approach revealed a highly polygenic basis (n = 190 QTLs total) for pollination syndrome divergence, capturing minor QTLs even for pigment traits with leading major loci. There was significant QTL overlap within pigment and dimensional categories. Nectar volume QTLs clustered with those for floral dimensions, suggesting a partially shared module. The NILs refined two pistil length QTLs, only one of which has tightly correlated effects on other dimensional traits. An overall polygenic architecture of floral divergence is partially coordinated by genetic modules formed by linkage (pigments) and likely pleiotropy (dimensions plus nectar). This work illuminates pollinator syndrome diversification in a model radiation and generates a robust framework for molecular and ecological genomics.

Among- and within-population variation in germination response shapes ecological resilience in the Mediterranean cliff species Brassica incana

BACKGROUND AND AIMS

Understanding how plant species respond to extreme conditions is crucial for predicting their ecological resilience under climate change. Here, we aimed to forecast the ecological resilience of the Mediterranean cliff species Brassica incana (Brassicaceae) by estimating population variation in germination response under novel extreme environmental conditions.

METHODS

We investigated the thermal germination responses in 14 populations of B. incana by exposing seeds to temperatures within and outside conditions experienced in their local environment. Then, we quantified among- and within-population variation in germination response to extreme temperatures, estimated genotype-by-environment interactions (G × E) and tested if population performance at extreme temperatures is explained by local climate.

KEY RESULTS

We found significant among-population differences in germination response, a different level of within-population variability and different mechanisms underlying G × E patterns. Also, populations experiencing higher temperatures in their local environment showed a better performance at both cold and hot extremes while populations experiencing lower temperatures showed a limited ability to germinate under extreme conditions.

CONCLUSIONS

Our results suggest that populations experiencing higher temperatures in their local environment have a greater potential to face future thermal extreme conditions and their role is thus crucial to promote species ecological resilience.

Museum genomics reveals temporal genetic stasis and global genetic diversity in Arabidopsis thaliana

Global patterns of population genetic variation through time offer a window into evolutionary processes that maintain diversity. Over time, lineages may expand or contract their distribution, causing turnover in population genetic composition. At individual loci, migration, drift, environmental change (among other processes) may affect allele frequencies. Museum specimens of widely distributed species offer a unique window into the genetics of understudied populations and changes over time. Here, we sequenced genomes of 130 herbarium specimens and 91 new field collections of Arabidopsis thaliana and combined these with published genomes. We sought a broader view of genomic diversity across the species, and to test if population genomic composition is changing through time. We documented extensive and previously uncharacterized diversity in a range of populations in Africa, populations that are under threat from anthropogenic climate change. Through time, we did not find dramatic changes in genomic composition of populations. Instead, we found a pattern of genetic change every 100 years of the same magnitude seen when comparing Eurasian populations that are 185 km apart, potentially due to a combination of drift and changing selection. We found only mixed signals of polygenic adaptation at phenology and physiology QTL. We did find that genes conserved across eudicots show altered levels of directional allele frequency change, potentially due to variable purifying and background selection. Our study highlights how museum specimens can reveal new dimensions of population diversity and show how wild populations are evolving in recent history.

Testing the evolutionary theory of inversion polymorphisms in the yellow monkeyflower (Mimulus guttatus)

Chromosomal inversions have been implicated in a remarkable range of natural phenomena, but it remains unclear how much they contribute to standing genetic variation. Here, we evaluate 64 inversions that segregate within a single natural population of the yellow monkeyflower (Mimulus guttatus). Nucleotide diversity patterns confirm low internal variation for the derived orientation (predicted by recent origin), elevated diversity between orientations (predicted by natural selection), and localized fluctuations (predicted by gene flux). Sequence divergence between orientations varies idiosyncratically by position, not following the suspension bridge pattern predicted if the breakpoints are the targets of selection. Genetic variation in gene expression is not inflated close to inversion breakpoints but is clearly partitioned between orientations. Like sequence variation, the pattern of expression variation suggests that the capture of coadapted alleles is more important than the breakpoints for the fitness effects of inversions. This work confirms several evolutionary predictions for inversion polymorphisms, but clarity emerges only by synthesizing estimates across many loci.

Pervasive fitness trade-offs revealed by rapid adaptation in large experimental populations of Drosophila melanogaster

Life-history trade-offs are an inherent feature of organismal biology that evolutionary theory posits play a key role in patterns of divergence within and between species. Efforts to quantify trade-offs are largely confined to phenotypic measurements and the identification of negative genetic-correlations among fitness-relevant traits. Here, we use time-series genomic data collected during experimental evolution in large, genetically diverse populations of Drosophila melanogaster to directly measure the manifestation of trade-offs in response to temporally fluctuating selection pressures on ecological timescales. Specifically, we quantify the genome-wide signal of antagonistic pleiotropy suggestive of trade-offs between reproduction and stress tolerance. We further identify a putative role of two cosmopolitan inversions in these trade-offs, and show that loci responding to selection during lab-based, reproduction selection exhibit signals of fluctuating selection in an outdoor mesocosm exposed to natural environmental conditions. Our results demonstrate the utility of time-series genomic data in revealing the presence and genomic architecture underlying fitness trade-offs, and add credence to models positing a role of generic life history trade-offs in the maintenance of variation in natural populations.

36-year study reveals stability of a wild wheat population across microhabitats

Long-term genetic studies of wild populations are very scarce, but are essential for connecting ecological and population genetics models, and for understanding the dynamics of biodiversity. We present a study of a wild wheat population sampled over a 36-year period at high spatial resolution. We genotyped 832 individuals from regular sampling along transects during the course of the experiment. Genotypes were clustered into ecological microhabitats over scales of tens of metres, and this clustering was remarkably stable over the 36 generations of the study. Simulations show that it is difficult to determine whether this spatial and temporal stability reflects extremely limited dispersal or fine-scale local adaptation to ecological parameters. Using a common-garden experiment, we showed that the genotypes found in distinct microhabitats differ phenotypically. Our results provide a rare insight into the population genetics of a natural population over a long monitoring period.

Early Life History Divergence Mediates Elevational Adaptation in a Perennial Alpine Plant

Spatially divergent natural selection can drive adaptation to contrasting environments and thus the evolution of ecotypes. In perennial plants, selection shapes life history traits by acting on subsequent life stages, each contributing to fitness. While evidence of adaptation in perennial plants is common, the expression of life history traits is rarely characterized, limiting our understanding of their role in adaptive evolution. We conducted a multi-year reciprocal transplant experiment with seedlings from low and high elevation populations of the alpine carnation Dianthus carthusianorum to test for adaptation linked to contrasting climates and inferred specific contributions of early life stages to fitness. We assessed genotype by environment interactions in single fitness components, applied matrix population models to achieve an integrated estimate of fitness through population growth rates, and performed trade-off analyses to investigate the advantage of alternate life history traits across environments. We found evidence of genotype by environment interactions consistent with elevational adaptation at multiple stages of the early life cycle. Estimates of population growth rates corroborated a strong advantage of the local genotype. Early reproduction and survival are alternate major contributors to adaptation at low and high elevation, respectively, and are linked by trade-offs that underlie the evolution of divergent life history traits across environments. While these traits have a strong genetic basis, foreign populations express co-gradient plasticity, reflecting the adaptive strategy of the local populations. Our study reveals that selection associated to climate has driven the evolution of divergent life histories and the formation of elevational ecotypes. While the high energy environment and strong competition favor investment in early reproduction at low elevation, limiting resources favor a more conservative strategy relying on self-maintenance at high elevation. The co-gradient plasticity expressed by high-elevation populations may facilitate their persistence under warming climatic conditions.

Short-term fluctuating and long-term divergent selection on sympatric Monkeyflowers: insights from decade-spanning reciprocal transplants

Sympatric species are often locally adapted to distinct microhabitats. However, temporal variation may cause local maladaptation and species boundary breakdown, especially during extreme climatic events leading to episodic selection. Repeated reciprocal transplants can reveal the interplay between short and long-term patterns of natural selection. To examine evolutionary trajectories of sympatric Monkeyflowers adapted to different niches, Mimulus guttatus and M. laciniatus, we performed three replicated transplants and combined them with previous experiments to leverage a dataset of five transplants spanning 10 years. We performed phenotypic selection analyses on parents and hybrids in parental habitats in Yosemite NP, CA during years of drastically differing snowpack. If there is ecological isolation, then we predicted divergent phenotypic selection between habitats in line with species' differences and local adaptation. We found interannual fluctuations in phenotypic selection, often in unpredicted directions. However, a combined-year analysis detected longer-term divergent selection on flowering time, a key temporally isolating and adaptative trait, suggesting that selection may reinforce species boundaries despite short-term fluctuations. Finally, we found temporal variation in local adaptation with M. laciniatus locally adapted in low snowpack years, while an extremely high snowpack year contributed to average local maladaptation of M. guttatus.

Experimental evolution suggests rapid assembly of the 'selfing syndrome' from standing variation in Mimulus guttatus

Ecological and evolutionary changes are likely to occur rapidly when outcrossing populations experience pollinator loss. However, the number and identify of plant traits that will respond to this form of selection, as well as the overall predictability of evolutionary responses, remain unclear. We experimentally evolved 20 large replicate populations of Mimulus guttatus for 10 generations under three treatments: pure outcrossing, mixed mating (10% outcrossing) and pure selfing. These populations were founded from the same genetically diverse and outcrossing natural population. After 10 generations, all measured traits evolved with flower size, phenology, and reproductive traits diverging consistently among mating system treatments. Autogamy increased dramatically in the selfing treatment, but the magnitude of adaptation only becomes clear once inbreeding depression is factored out. Selfing treatment plants evolved reduced stigma-anther separation, and also exhibited declines in flower size and per-flower reproductive capacity. Flower size also declined in selfing populations but this was driven mainly by inbreeding depression and cannot be attributed to adaptation towards the selfing syndrome. Generally, the mixed mating populations evolved trait values intermediate to the fully selfing and outcrossing populations. Overall, our experimental treatments reiterated differences that have been documented in interspecific comparisons between selfing and outcrossing species pairs. Given that such contrasts involve species separated by thousands or even millions of generations, it is noteworthy that large evolutionary responses were obtained from genetic variation segregating within a single natural population.

A tradeoff between bacteriophage resistance and bacterial motility is mediated by the Rcs phosphorelay in Escherichia coli

Across the tree of life, pleiotropy is thought to constrain adaptation through evolutionary tradeoffs. However, few examples of pleiotropy exist that are well explained at the genetic level, especially for pleiotropy that is mediated by multiple genes. Here, we describe a set of pleiotropic mutations that mediate two key fitness components in bacteria: parasite resistance and motility. We subjected Escherichia coli to strong selection by phage U136B to obtain 27 independent mucoid mutants. Mucoidy is a phenotype that results from excess exopolysaccharide and can act as a barrier against viral infection but can also interfere with other cellular functions. We quantified the mutants' phage resistance using efficiency of plaquing assays and swimming motility using swim agar plates, and we sequenced the complete genomes of all mutants to identify mucoid-causing mutations. Increased phage resistance co-occurred with decreased motility. This relationship was mediated by highly parallel (27/27) mutations to the Rcs phosphorelay pathway, which senses membrane stress to regulate exopolysaccharide production. Together, these results provide an empirical example of a pleiotropic relationship between two traits with intermediate genetic complexity.

Polygenic selection to a changing optimum under self-fertilisation

Many traits are polygenic, affected by multiple genetic variants throughout the genome. Selection acting on these traits involves co-ordinated allele-frequency changes at these underlying variants, and this process has been extensively studied in random-mating populations. Yet many species self-fertilise to some degree, which incurs changes to genetic diversity, recombination and genome segregation. These factors cumulatively influence how polygenic selection is realised in nature. Here, we use analytical modelling and stochastic simulations to investigate to what extent self-fertilisation affects polygenic adaptation to a new environment. Our analytical solutions show that while selfing can increase adaptation to an optimum, it incurs linkage disequilibrium that can slow down the initial spread of favoured mutations due to selection interference, and favours the fixation of alleles with opposing trait effects. Simulations show that while selection interference is present, high levels of selfing (at least 90%) aids adaptation to a new optimum, showing a higher long-term fitness. If mutations are pleiotropic then only a few major-effect variants fix along with many neutral hitchhikers, with a transient increase in linkage disequilibrium. These results show potential advantages to self-fertilisation when adapting to a new environment, and how the mating system affects the genetic composition of polygenic selection.

Biotic interactions promote local adaptation to soil in plants

Although different ecological factors shape adaptative evolution in natural habitats, we know little about how their interactions impact local adaptation. Here we used eight generations of experimental evolution with outcrossing Brassica rapa plants as a model system, in eight treatment groups that varied in soil type, herbivory (with/without aphids), and pollination mode (hand- or bumblebee-pollination), to study how biotic interactions affect local adaptation to soil. First, we show that several plant traits evolved in response to biotic interactions in a soil-specific way. Second, using a reciprocal transplant experiment, we demonstrate that significant local adaptation to soil-type evolved in the "number of open flowers", a trait used as a fitness proxy, but only in plants that evolved with herbivory and bee pollination. Whole genome re-sequencing of experimental lines revealed that biotic interactions caused a 10-fold increase in the number of SNPs across the genome with significant allele frequency change, and that alleles with opposite allele frequency change in different soil types (antagonistic pleiotropy) were most common in plants with an evolutionary history of herbivory and bee pollination. Our results demonstrate that the interaction with mutualists and antagonists can facilitate local adaptation to soil type through antagonistic pleiotropy.

Weak response to selection on stigma-anther distance in a primarily selfing population of yellow monkeyflower

Stebbins hypothesized that selfing lineages are evolutionary dead ends because they lack adaptive potential. While selfing populations often possess limited nucleotide variability compared with closely related outcrossers, reductions in the genetic variability of quantitative characters remain unclear, especially for key traits determining selfing rates. Yellow monkeyflower (Mimulus guttatus) populations generally outcross and maintain extensive quantitative genetic variation in floral traits. Here, we study the Joy Road population (Bodega Bay, CA, USA) of M. guttatus, where individuals exhibit stigma-anther distances (SAD) typical of primarily selfing monkeyflowers. We show that this population is closely related to nearby conspecifics on the Pacific Coast with a modest 33% reduction in genome-wide variation compared with a more highly outcrossing population. A five-generation artificial selection experiment challenged the hypothesis that the Joy Road population harbours comparatively low evolutionary potential in stigma-anther distance, a critical determinant of selfing rate in Mimulus. Artificial selection generated a weak phenotypic response, with low realized heritabilities (0.020-0.028) falling 84% below those measured for floral characters in more highly outcrossing M. guttatus. These results demonstrate substantial declines in evolutionary potential with a transition toward selfing. Whether these findings explain infrequent reversals to outcrossing or general limits on adaptation in selfers requires further investigation.

Strategies to identify and dissect trade-offs in plants

Trade-offs between traits arise and reflect constraints imposed by the environment and physicochemical laws. Trade-off situations are expected to be highly relevant for sessile plants, which have to respond to changes in the environment to ensure survival. Despite increasing interest in determining the genetic and molecular basis of plant trade-offs, there are still gaps and differences with respect to how trade-offs are defined, how they are measured, and how their genetic architecture is dissected. The first step to fill these gaps is to establish what is meant by trade-offs. In this review we provide a classification of the existing definitions of trade-offs according to: (1) the measures used for their quantification, (2) the dependence of trade-offs on environment, and (3) experimental designed used (i.e. a single individual across different environments or a population of individuals in single or multiple environments). We then compare the approaches for quantification of trade-offs based on phenotypic, between-individual, and genetic correlations, and stress the need for developing further quantification indices particularly for trade-offs between multiple traits. Lastly, we highlight the genetic mechanisms underpinning trade-offs and experimental designs that facilitate their discovery in plants, with focus on usage of natural variability. This review also offers a perspective for future research aimed at identification of plant trade-offs, dissection of their genetic architecture, and development of strategies to overcome trade-offs, with applications in crop breeding.

The quantitative genetics of gene expression in Mimulus guttatus

Gene expression can be influenced by genetic variants that are closely linked to the expressed gene (cis eQTLs) and variants in other parts of the genome (trans eQTLs). We created a multiparental mapping population by sampling genotypes from a single natural population of Mimulus guttatus and scored gene expression in the leaves of 1,588 plants. We find that nearly every measured gene exhibits cis regulatory variation (91% have FDR < 0.05). cis eQTLs are usually allelic series with three or more functionally distinct alleles. The cis locus explains about two thirds of the standing genetic variance (on average) but varies among genes and tends to be greatest when there is high indel variation in the upstream regulatory region and high nucleotide diversity in the coding sequence. Despite mapping over 10,000 trans eQTL / affected gene pairs, most of the genetic variance generated by trans acting loci remains unexplained. This implies a large reservoir of trans acting genes with subtle or diffuse effects. Mapped trans eQTLs show lower allelic diversity but much higher genetic dominance than cis eQTLs. Several analyses also indicate that trans eQTLs make a substantial contribution to the genetic correlations in expression among different genes. They may thus be essential determinants of "gene expression modules," which has important implications for the evolution of gene expression and how it is studied by geneticists.

Dominance reversals: the resolution of genetic conflict and maintenance of genetic variation

Beneficial reversals of dominance reduce the costs of genetic trade-offs and can enable selection to maintain genetic variation for fitness. Beneficial dominance reversals are characterized by the beneficial allele for a given context (e.g. habitat, developmental stage, trait or sex) being dominant in that context but recessive where deleterious. This context dependence at least partially mitigates the fitness consequence of heterozygotes carrying one non-beneficial allele for their context and can result in balancing selection that maintains alternative alleles. Dominance reversals are theoretically plausible and are supported by mounting empirical evidence. Here, we highlight the importance of beneficial dominance reversals as a mechanism for the mitigation of genetic conflict and review the theory and empirical evidence for them. We identify some areas in need of further research and development and outline three methods that could facilitate the identification of antagonistic genetic variation (dominance ordination, allele-specific expression and allele-specific ATAC-Seq (assay for transposase-accessible chromatin with sequencing)). There is ample scope for the development of new empirical methods as well as reanalysis of existing data through the lens of dominance reversals. A greater focus on this topic will expand our understanding of the mechanisms that resolve genetic conflict and whether they maintain genetic variation.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

Transposition, duplication, and divergence of the telomerase RNA underlies the evolution of Mimulus telomeres

Telomeres are nucleoprotein complexes with a crucial role of protecting chromosome ends. It consists of simple repeat sequences and dedicated telomere-binding proteins. Because of its vital functions, components of the telomere, for example its sequence, should be under strong evolutionary constraint. But across all plants, telomere sequences display a range of variation and the evolutionary mechanism driving this diversification is largely unknown. Here, we discovered in Monkeyflower (Mimulus) the telomere sequence is even variable between species. We investigated the basis of Mimulus telomere sequence evolution by studying the long noncoding telomerase RNA (TR), which is a core component of the telomere maintenance complex and determines the telomere sequence. We conducted total RNA-based de novo transcriptomics from 16 Mimulus species and analyzed reference genomes from 6 species, and discovered Mimulus species have evolved at least three different telomere sequences: (AAACCCT)n, (AAACCCG)n, and (AAACCG)n. Unexpectedly, we discovered several species with TR duplications and the paralogs had functional consequences that could influence telomere evolution. For instance, M. lewisii had two sequence-divergent TR paralogs and synthesized a telomere with sequence heterogeneity, consisting of AAACCG and AAACCCG repeats. Evolutionary analysis of the M. lewisii TR paralogs indicated it had arisen from a transposition-mediate duplication process. Further analysis of the TR from multiple Mimulus species showed the gene had frequently transposed and inserted into new chromosomal positions during Mimulus evolution. From our results, we propose the TR transposition, duplication, and divergence model to explain the evolutionary sequence turnovers in Mimulus and potentially all plant telomeres.

Adaptation of perennial flowering phenology across the European range of Arabis alpina

Flowering phenology is important in the adaptation of many plants to their local environment, but its adaptive value has not been extensively studied in herbaceous perennials. We used Arabis alpina as a model system to determine the importance of flowering phenology to fitness of a herbaceous perennial with a wide geographical range. Individual plants representative of local genetic diversity (accessions) were collected across Europe, including in Spain, the Alps and Scandinavia. The flowering behaviour of these accessions was documented in controlled conditions, in common-garden experiments at native sites and in situ in natural populations. Accessions from the Alps and Scandinavia varied in whether they required exposure to cold (vernalization) to induce flowering, and in the timing and duration of flowering. By contrast, all Spanish accessions obligately required vernalization and had a short duration of flowering. Using experimental gardens at native sites, we show that an obligate requirement for vernalization increases survival in Spain. Based on our analyses of genetic diversity and flowering behaviour across Europe, we propose that in the model herbaceous perennial A. alpina, an obligate requirement for vernalization, which is correlated with short duration of flowering, is favoured by selection in Spain where the plants experience a long growing season.

The quantitative genetics of gene expression in Mimulus guttatus

Gene expression can be influenced by genetic variants that are closely linked to the expressed gene (cis eQTLs) and variants in other parts of the genome (trans eQTLs). We created a multiparental mapping population by sampling genotypes from a single natural population of Mimulus guttatus and scored gene expression in the leaves of 1,588 plants. We find that nearly every measured gene exhibits cis regulatory variation (91% have FDR < 0.05) and that cis eQTLs are usually allelic series with three or more functionally distinct alleles. The cis locus explains about two thirds of the standing genetic variance (on average) but varies among genes and tends to be greatest when there is high indel variation in the upstream regulatory region and high nucleotide diversity in the coding sequence. Despite mapping over 10,000 trans eQTL / affected gene pairs, most of the genetic variance generated by trans acting loci remains unexplained. This implies a large reservoir of trans acting genes with subtle or diffuse effects. Mapped trans eQTLs show lower allelic diversity but much higher genetic dominance than cis eQTLs. Several analyses also indicate that trans eQTL make a substantial contribution to the genetic correlations in expression among different genes. They may thus be essential determinants of "gene expression modules", which has important implications for the evolution of gene expression and also how it is studied by geneticists.

Demographic consequences of an extreme heat wave are mitigated by spatial heterogeneity in an annual monkeyflower

Heat waves are becoming more frequent and intense with climate change, but the demographic and evolutionary consequences of heat waves are rarely investigated in herbaceous plant species. We examine the consequences of a short but extreme heat wave in Oregon populations of the common yellow monkeyflower (Mimulus guttatus) by leveraging a common garden experiment planted with range-wide populations and observational studies of 11 local populations. In the common garden, 89% of seedlings died during the heat wave including >96% of seedlings from geographically local populations. Some populations from hotter and drier environments had higher fitness, however, others from comparable environments performed poorly. Observational studies of local natural populations drastically differed in the consequences of the heat wave-one population was completely extirpated and nearly half had a >50% decrease in fitness. However, a few populations had greater fitness during the heat wave year. Differences in mortality corresponded to the impact of the heat wave on soil moisture-retention of soil moisture throughout the heat wave led to greater survivorship. Our results suggest that not all populations experience the same intensity or degree of mortality during extreme events and such heterogeneity could be important for genetic rescue or to facilitate the distribution of adaptive variants throughout the region.

Unique drought resistance strategies occur among monkeyflower populations spanning an aridity gradient

PREMISE

Annual plants often exhibit drought-escape and avoidance strategies to cope with limited water availability. Determining the extent of variation and factors underlying the evolution of divergent strategies is necessary for determining population responses to more frequent and severe droughts.

METHODS

We leveraged five Mimulus guttatus populations collected across an aridity gradient within manipulative drought and quantitative genetics experiments to examine constitutive and terminal-drought induced responses in drought resistance traits.

RESULTS

Populations varied considerably in drought-escape- and drought-avoidance-associated traits. The most mesic population demonstrated a unique resource conservative strategy. Xeric populations exhibited extreme plasticity when exposed to terminal drought that included flowering earlier at shorter heights, increasing water-use efficiency, and shifting C:N ratios. However, plasticity responses also differed between populations, with two populations slowing growth rates and flowering at earlier nodes and another population increasing growth rate. While nearly all traits were heritable, phenotypic correlations differed substantially between treatments and often, populations.

CONCLUSIONS

Our results suggest drought resistance strategies of populations may be finely adapted to local patterns of water availability. Substantial plastic responses suggest that xeric populations can already acclimate to drought through plasticity, but populations not frequently exposed to drought may be more vulnerable.

The design of mapping populations: Impacts of geographic scale on genetic architecture and mapping efficacy for defense and immunity

Genome-wide association studies (GWAS) have yielded tremendous insight into the genetic architecture of trait variation. However, the collections of loci they uncover are far from exhaustive. As many of the complicating factors that confound or limit the efficacy of GWAS are exaggerated over broad geographic scales, a shift toward more analyses using mapping panels sampled from narrow geographic localities ("local" populations) could provide novel, complementary insights. Here, we present an overview of the major complicating factors, review mounting evidence from genomic analyses that these factors are pervasive, and synthesize theoretical and empirical evidence for the power of GWAS in local populations.

No evidence for schistosome parasite fitness trade-offs in the intermediate and definitive host

BACKGROUND

The trematode parasite Schistosoma mansoni uses an aquatic snail intermediate and a vertebrate definitive host to complete its life cycle. We previously showed that a key transmission trait-the number of cercariae larvae shed from infected Biomphalaria spp. snails-varies significantly within and between different parasite populations and is genetically controlled by five loci. We investigated the hypothesis that the success of parasite genotypes showing high propagative fitness in the intermediate snail host may be offset by lower reproductive fitness in the definitive vertebrate host.

METHODS

We investigated this trade-off hypothesis by selecting parasite progeny producing high or low number of larvae in the snail and then comparing fitness parameters and virulence in the rodent host. We infected inbred BALB/c mice using two Schistosoma mansoni parasite lines [high shedder (HS) and low shedder (LS) lines] isolated from F2 progeny generated by genetic crosses between SmLE (HS parent) and SmBRE (LS parent) parasites. We used the F3 progeny to infect two populations of inbred Biomphalaria glabrata snails. We then compared life history traits and virulence of these two selected parasite lines in the rodent host to understand pleiotropic effects of genes determining cercarial shedding in parasites infecting the definitive host.

RESULTS

HS parasites shed high numbers of cercariae, which had a detrimental impact on snail physiology (measured by laccase-like activity and hemoglobin rate), regardless of the snail genetic background. In contrast, selected LS parasites shed fewer cercariae and had a lower impact on snail physiology. Similarly, HS worms have a higher reproductive fitness and produced more viable F3 miracidia larvae than LS parasites. This increase in transmission is correlated with an increase in virulence toward the rodent host, characterized by stronger hepato-splenomegaly and hepatic fibrosis.

CONCLUSIONS

These experiments revealed that schistosome parasite propagative and reproductive fitness was positively correlated in intermediate and definitive host (positive pleiotropy). Therefore, we rejected our trade-off hypothesis. We also showed that our selected schistosome lines exhibited low and high shedding phenotype regardless of the intermediate snail host genetic background. ​.

Re-thinking the environment in landscape genomics

Detecting the extrinsic selective pressures shaping genomic variation is critical for a better understanding of adaptation and for forecasting evolutionary responses of natural populations to changing environmental conditions. With increasing availability of geo-referenced environmental data, landscape genomics provides unprecedented insights into how genomic variation and underlying gene functions affect traits potentially under selection. Yet, the robustness of genotype-environment associations used in landscape genomics remains tempered due to various limitations, including the characteristics of environmental data used, sampling designs employed, and statistical frameworks applied. Here, we argue that using complementary or new environmental data sources and well-informed sampling designs may help improve the detection of selective pressures underlying patterns of local adaptation in various organisms and environments.

Genotype accounts for intraspecific variation in the timing and duration of multiple, sequential life-cycle events in a willow species

PREMISE

Phenological variation among individuals within populations is common and has a variety of ecological and evolutionary consequences, including forming the basis for population-level responses to environmental change. Although the timing of life-cycle events has genetic underpinnings, whether intraspecific variation in the duration of life-cycle events reflects genetic differences among individuals is poorly understood.

METHODS

We used a common garden experiment with 10 genotypes of Salix hookeriana (coastal willow) from northern California, United States to investigate the extent to which genetic variation explains intraspecific variation in the timing and duration of multiple, sequential life-cycle events: flowering, leaf budbreak, leaf expansion, fruiting, and fall leaf coloration. We used seven clones of each genotype, for a total of 70 individual trees.

RESULTS

Genotype affected each sequential life-cycle event independently and explained on average 62% of the variation in the timing and duration of vegetative and reproductive life-cycle events. All events were significantly heritable. A single genotype tended to be "early" or "late" across life-cycle events, but for event durations, there was no consistent response within genotypes.

CONCLUSIONS

This research demonstrates that genetic variation can be a major component underlying intraspecific variation in the timing and duration of life-cycle events. It is often assumed that the environment affects durations, but we show that genetic factors also play a role. Because the timing and duration of events are independent of one another, our results suggest that the effects of environmental change on one event will not necessarily cascade to subsequent events.

Taking quantitative genomics into the wild

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.

The potential for floral evolution in response to competing selection pressures following the loss of hawkmoth pollination in Ruellia humilis

PREMISE

In the absence of hawkmoth pollinators, chasmogamous (CH) flowers of Ruellia humilis self-pollinate by two secondary mechanisms. Other floral visitors might exert selection on CH floral traits to restore outcrossing, but at the same time preferential predation of CH seeds generates selection to increase the allocation of resources to cleistogamous (CL) flowers.

METHODS

To assess the potential for an evolutionary response to these competing selection pressures, we estimated additive genetic variances ( σ A 2 ${\sigma }_{{\rm{A}}}^{2}$ ) and covariances for 14 reproductive traits and three fitness components in a Missouri population lacking hawkmoth pollinators.

RESULTS

We found significant σ A 2 ${\sigma }_{{\rm{A}}}^{2}$ for all 11 floral traits and two measures of resource allocation to CL flowers, indicating the potential for a short-term response to selection on most reproductive traits. Selection generated by seed predators is predicted to increase the percentage of CL flowers by 0.24% per generation, and mean stigma-anther separation is predicted to decrease as a correlated response, increasing the fraction of plants that engage in prior selfing. However, the initial response to this selection is opposed by strong directional dominance.

CONCLUSIONS

The predicted evolutionary decrease in the number of CH flowers available for potential outcrossing, combined with the apparent preclusion of potential diurnal pollinators by the pollen-harvesting activities of sweat bees, suggest that 100% cleistogamy is the likely outcome of evolution in the absence of hawkmoths. However, rare mutations with large effects, such as delaying budbreak until after sunrise, could provide pathways for the restoration of outcrossing that are not reachable by gradual quantitative-genetic evolution.

Climate and the biotic community structure plant resistance across biogeographic groups of yellow monkeyflower

Characterizing correlates of phytochemical resistance trait variation across a landscape can provide insight into the ecological factors that have shaped the evolution of resistance arsenals. Using field-collected data and a greenhouse common garden experiment, we assessed the relative influences of abiotic and biotic drivers of genetic-based defense trait variation across 41 yellow monkeyflower populations from western and eastern North America and the United Kingdom. Populations experience different climates, herbivore communities, and neighboring vegetative communities, and have distinct phytochemical resistance arsenals. Similarities in climate as well as herbivore and vegetative communities decline with increasing physical distance separating populations, and phytochemical resistance arsenal composition shows a similarly decreasing trend. Of the abiotic and biotic factors examined, temperature and the neighboring vegetation community had the strongest relative effects on resistance arsenal differentiation, whereas herbivore community composition and precipitation have relatively small effects. Rather than simply controlling for geographic proximity, we jointly assessed the relative strengths of both geographic and ecological variables on phytochemical arsenal compositional dissimilarity. Overall, our results illustrate how abiotic conditions and biotic interactions shape plant defense traits in natural populations.

The influence of genetic architecture on responses to selection under drought in rice

Accurately predicting responses to selection is a major goal in biology and important for successful crop breeding in changing environments. However, evolutionary responses to selection can be constrained by such factors as genetic and cross-environment correlations, linkage, and pleiotropy, and our understanding of the extent and impact of such constraints is still developing. Here, we conducted a field experiment to investigate potential constraints to selection for drought resistance in rice (Oryza sativa) using phenotypic selection analysis and quantitative genetics. We found that traits related to drought response were heritable, and some were under selection, including selection for earlier flowering, which could allow drought escape. However, patterns of selection generally were not opposite under wet and dry conditions, and we did not find individual or closely linked genes that influenced multiple traits, indicating a lack of evidence that antagonistic pleiotropy, linkage, or cross-environment correlations would constrain selection for drought resistance. In most cases, genetic correlations had little influence on responses to selection, with direct and indirect selection largely congruent. The exception to this was seed mass under drought, which was predicted to evolve in the opposite direction of direct selection due to correlations. Because of this indirect effect on selection on seed mass, selection for drought resistance was not accompanied by a decrease in seed mass, and yield increased with fecundity. Furthermore, breeding lines with high fitness and yield under drought also had high fitness and yield under wet conditions, indicating that there was no evidence for a yield penalty on drought resistance. We found multiple genes in which expression influenced both water use efficiency (WUE) and days to first flowering, supporting a genetic basis for the trade-off between drought escape and avoidance strategies. Together, these results can provide helpful guidance for understanding and managing evolutionary constraints and breeding stress-resistant crops.

Mutualist- and antagonist-mediated selection contribute to trait diversification of flowers

Flowers are generally short-lived, and they all face a multidimensional challenge because they have to attract mutualists, compel them to vector pollen with minimal investment in rewards, and repel floral enemies during this short time window. Their displays are under complex selection, either consistent or conflicting, to maximize reproductive fitness under heterogeneous environments. The phenological or morphological mismatches between flowers and visitors will influence interspecific competition, resource access, mating success and, ultimately, population and community dynamics. To better understand the effects of the plant visitors on floral traits, it is necessary to determine the functional significance of specific floral traits for the visitors; how plants respond to both mutualists and antagonists through adaptive changes; and to evaluate the net fitness effects of biological mutualisms and antagonism on plants. In this review, we bring together insights from fields as diverse as floral biology, insect behavioral responses, and evolutionary biology to explain the processes and patterns of floral diversity evolution. Then, we discuss the ecological significance of plant responses to mutualists and antagonists from a community perspective, and propose a set of research questions that can guide the research field to integrate studies of plant defense and reproduction.

Pollinator loss causes rapid adaptive evolution of selfing and dramatically reduces genome-wide genetic variability

Although selfing populations harbor little genetic variation limiting evolutionary potential, the causes are unclear. We experimentally evolved large, replicate populations of Mimulus guttatus for nine generations in greenhouses with or without pollinating bees and studied DNA polymorphism in descendants. Populations without bees adapted to produce more selfed seed yet exhibited striking reductions in DNA polymorphism despite large population sizes. Importantly, the genome-wide pattern of variation cannot be explained by a simple reduction in effective population size, but instead reflects the complicated interaction between selection, linkage, and inbreeding. Simulations demonstrate that the spread of favored alleles at few loci depresses neutral variation genome wide in large populations containing fully selfing lineages. It also generates greater heterogeneity among chromosomes than expected with neutral evolution in small populations. Genome-wide deviations from neutrality were documented in populations with bees, suggesting widespread influences of background selection. After applying outlier tests to detect loci under selection, two genome regions were found in populations with bees, yet no adaptive loci were otherwise mapped. Large amounts of stochastic change in selfing populations compromise evolutionary potential and undermine outlier tests for selection. This occurs because genetic draft in highly selfing populations makes even the largest changes in allele frequency unremarkable.

Supergene potential of a selfish centromere

Selfishly evolving centromeres bias their transmission by exploiting the asymmetry of female meiosis and preferentially segregating to the egg. Such female meiotic drive systems have the potential to be supergenes, with multiple linked loci contributing to drive costs or enhancement. Here, we explore the supergene potential of a selfish centromere (D) in Mimulus guttatus, which was discovered in the Iron Mountain (IM) Oregon population. In the nearby Cone Peak population, D is still a large, non-recombining and costly haplotype that recently swept, but shorter haplotypes and mutational variation suggest a distinct population history. We detected D in five additional populations spanning more than 200 km; together, these findings suggest that selfish centromere dynamics are widespread in M. guttatus. Transcriptome comparisons reveal elevated differences in expression between driving and non-driving haplotypes within, but not outside, the drive region, suggesting large-scale cis effects of D's spread on gene expression. We use the expression data to refine linked candidates that may interact with drive, including Nuclear Autoantigenic Sperm Protein (NASP^SIM3), which chaperones the centromere-defining histone CenH3 known to modify Mimulus drive. Together, our results show that selfishly evolving centromeres may exhibit supergene behaviour and lay the foundation for future genetic dissection of drive and its costs. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Population genetics of transposable element load: A mechanistic account of observed overdispersion

In an empirical analysis of transposable element (TE) abundance within natural populations of Mimulus guttatus and Drosophila melanogaster, we found a surprisingly high variance of TE count (e.g., variance-to-mean ratio on the order of 10 to 300). To obtain insight regarding the evolutionary genetic mechanisms that underlie the overdispersed population distributions of TE abundance, we developed a mathematical model of TE population genetics that includes the dynamics of element proliferation and purifying selection on TE load. The modeling approach begins with a master equation for a birth-death process and extends the predictions of the classical theory of TE dynamics in several ways. In particular, moment-based analyses of population distributions of TE load reveal that overdispersion is likely to arise via copy-and-paste proliferation dynamics, especially when the elementary processes of proliferation and excision are approximately balanced. Parameter studies and analytic work confirm this result and further suggest that overdispersed population distributions of TE abundance are probably not a consequence of purifying selection on total element load.

Population genomic consequences of life-history and mating system adaptation to a geothermal soil mosaic in yellow monkeyflowers

Local selection can promote phenotypic divergence despite gene flow across habitat mosaics, but adaptation itself may generate substantial barriers to genetic exchange. In plants, life-history, phenology, and mating system divergence have been proposed to promote genetic differentiation in sympatry. In this study, we investigate phenotypic and genetic variation in Mimulus guttatus (yellow monkeyflowers) across a geothermal soil mosaic in Yellowstone National Park (YNP). Plants from thermal annual and nonthermal perennial habitats were heritably differentiated for life-history and mating system traits, consistent with local adaptation to the ephemeral thermal-soil growing season. However, genome-wide genetic variation primarily clustered plants by geographic region, with little variation sorting by habitat. The one exception was an extreme thermal population also isolated by a 200 m geographical gap of no intermediate habitat. Individual inbreeding coefficients (F_IS ) were higher (and predicted by trait variation) in annual plants and annual pairs showed greater isolation by distance at local (<1 km) scales. Finally, YNP adaptation does not reuse a widespread inversion that underlies M. guttatus life-history ecotypes range-wide, suggesting a novel genetic mechanism. Overall, this work suggests that life-history and mating system adaptation strong enough to shape individual mating patterns does not necessarily generate incipient speciation without geographical barriers.

Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the allele frequency spectrum of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with "rarity disequilibrium" produce a superabundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of coexpressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high "connectedness" (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

Standard metabolic rate does not associate with age-at-maturity genotype in juvenile Atlantic salmon

Atlantic salmon (Salmo salar) is a species with diverse life-history strategies, to which the timing of maturation contributes considerably. Recently, the genome region including the gene vgll3 has gained attention as a locus with a large effect on Atlantic salmon maturation timing, and recent studies on the vgll3 locus in salmon have indicated that its effect might be mediated through body condition and accumulation of adipose tissue. However, the cellular and physiological pathways leading from vgll3 genotype to phenotype are still unknown. Standard metabolic rate is a potentially important trait for resource acquisition and assimilation and we hypothesized that this trait, being a proxy for the maintenance energy expenditure of an individual, could be an important link in the pathway from vgll3 genotype to maturation timing phenotype. As a first step to studying links between vgll3 and the metabolic phenotype of Atlantic salmon, we measured the standard metabolic rate of 150 first-year Atlantic salmon juveniles of both sexes, originating from 14 different families with either late-maturing or early-maturing vgll3 genotypes. No significant difference in mass-adjusted standard metabolic rate was detected between individuals with different vgll3 genotypes, indicating that juvenile salmon of different vgll3 genotypes have similar maintenance energy requirements in the experimental conditions used and that the effects of vgll3 on body condition and maturation are not strongly related to maintenance energy expenditure in either sex at this life stage.

Domestication reshaped the genetic basis of inbreeding depression in a maize landrace compared to its wild relative, teosinte

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. The extent to which the genetic load of mutations contributing to inbreeding depression is due to large-effect mutations versus variants with very small individual effects is unknown and may be affected by population history. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits, congruent with the greater segregating genetic load in the maize population that we predicted from sequence data. Parental breeding values were highly consistent between outcross and selfed offspring, indicating that additive effects determine most of the genetic value even in the presence of strong inbreeding depression. We developed a novel linkage scan to identify quantitative trait loci (QTL) representing large-effect rare variants carried by only a single parent, which were more important in teosinte than maize. Teosinte also carried more putative juvenile-acting lethal variants identified by segregation distortion. These results suggest a mixture of mostly polygenic, small-effect partially recessive effects in linkage disequilibrium underlying inbreeding depression, with an additional contribution from rare larger-effect variants that was more important in teosinte but depleted in maize following the domestication bottleneck. Purging associated with the maize domestication bottleneck may have selected against some large effect variants, but polygenic load is harder to purge and overall segregating mutational burden increased in maize compared to teosinte.

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. Mating of close relatives increases the probability that an individual inherits two non-functioning mutations at the same gene, resulting in lower fitness of such matings. We do not know the extent to which inbreeding depression is due to mutations with large-effects versus small-effect polygenic variants. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits and we found that this was consistent with higher predicted ‘genetic load’ in maize based solely on the evolutionary conservation of the sequence variants observed in the population. We also mapped genome positions associated with inbreeding depression, identifying more and larger-effect genetic variants in teosinte than maize. These results suggest that during domestication, some of the rare large-effect variants in teosinte were bred out, but many genetic variants of small effects on inbreeding depression increased in frequency maize.

Identifying loci under selection via explicit demographic models

Adaptive genetic variation is a function of both selective and neutral forces. To accurately identify adaptive loci, it is thus critical to account for demographic history. Theory suggests that signatures of selection can be inferred using the coalescent, following the premise that genealogies of selected loci deviate from neutral expectations. Here, we build on this theory to develop an analytical framework to identify loci under selection via explicit demographic models (LSD). Under this framework, signatures of selection are inferred through deviations in demographic parameters, rather than through summary statistics directly, and demographic history is accounted for explicitly. Leveraging the property of demographic models to incorporate directionality, we show that LSD can provide information on the environment in which selection acts on a population. This can prove useful in elucidating the selective processes underlying local adaptation, by characterizing genetic trade-offs and extending the concepts of antagonistic pleiotropy and conditional neutrality from ecological theory to practical application in genomic data. We implement LSD via approximate Bayesian computation and demonstrate, via simulations, that LSD (a) has high power to identify selected loci across a large range of demographic-selection regimes, (b) outperforms commonly applied genome-scan methods under complex demographies and (c) accurately infers the directionality of selection for identified candidates. Using the same simulations, we further characterize the behaviour of isolation-with-migration models conducive to the study of local adaptation under regimes of selection. Finally, we demonstrate an application of LSD by detecting loci and characterizing genetic trade-offs underlying flower colour in Antirrhinum majus.

The promise and deceit of genomic selection component analyses

Selection component analyses (SCA) relate individual genotype to fitness components such as viability, fecundity and mating success. SCA are based on population genetic models and yield selection estimates directly in terms of predicted allele frequency change. This paper explores the statistical properties of gSCA: experiments that apply SCA to genome-wide scoring of SNPs in field sampled individuals. Computer simulations indicate that gSCA involving a few thousand genotyped samples can detect allele frequency changes of the magnitude that has been documented in field experiments on diverse taxa. To detect selection, imprecise genotyping from low-level sequencing of large samples of individuals provides much greater power than precise genotyping of smaller samples. The simulations also demonstrate the efficacy of 'haplotype matching', a method to combine information from a limited collection of whole genome sequence (the reference panel) with the much larger sample of field individuals that are measured for fitness. Pooled sequencing is demonstrated as another way to increase statistical power. Finally, I discuss the interpretation of selection estimates in relation to the Beavis effect, the overestimation of selection intensities at significant loci.

Genomically correlated trait combinations and antagonistic selection contributing to counterintuitive genetic patterns of adaptive diapause divergence in Rhagoletis flies

Adaptation to novel environments can result in unanticipated genomic responses to selection. Here, we illustrate how multifarious, correlational selection helps explain a counterintuitive pattern of genetic divergence between the recently derived apple- and ancestral hawthorn-infesting host races of Rhagoletis pomonella (Diptera: Tephritidae). The apple host race terminates diapause and emerges as adults earlier in the season than the hawthorn host race, to coincide with the earlier fruiting phenology of their apple hosts. However, alleles at many loci associated with later emergence paradoxically occur at higher frequencies in sympatric populations of the apple compared to the hawthorn race. We present genomic evidence that historical selection over geographically varying environmental gradients across North America generated genetic correlations between two life history traits, diapause intensity and diapause termination, in the hawthorn host race. Moreover, the loci associated with these life history traits are concentrated in genomic regions in high linkage disequilibrium (LD). These genetic correlations are antagonistic to contemporary selection on local apple host race populations that favours increased initial diapause depth and earlier, not later, diapause termination. Thus, the paradox of apple flies appears due, in part, to pleiotropy or linkage of alleles associated with later adult emergence and increased initial diapause intensity, the latter trait strongly selected for by the earlier phenology of apples. Our results demonstrate how understanding of multivariate trait combinations and the correlative nature of selective forces acting on them can improve predictions concerning adaptive evolution and help explain seemingly counterintuitive patterns of genetic diversity in nature.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

Growth-defense trade-offs shape population genetic composition in an iconic forest tree species

All organisms experience fundamental conflicts between divergent metabolic processes. In plants, a pivotal conflict occurs between allocation to growth, which accelerates resource acquisition, and to defense, which protects existing tissue against herbivory. Trade-offs between growth and defense traits are not universally observed, and a central prediction of plant evolutionary ecology is that context-dependence of these trade-offs contributes to the maintenance of intraspecific variation in defense [Züst and Agrawal, Annu. Rev. Plant Biol., 68, 513-534 (2017)]. This prediction has rarely been tested, however, and the evolutionary consequences of growth-defense trade-offs in different environments are poorly understood, especially in long-lived species [Cipollini et al., Annual Plant Reviews (Wiley, 2014), pp. 263-307]. Here we show that intraspecific trait trade-offs, even when fixed across divergent environments, interact with competition to drive natural selection of tree genotypes corresponding to their growth-defense phenotypes. Our results show that a functional trait trade-off, when coupled with environmental variation, causes real-time divergence in the genetic architecture of tree populations in an experimental setting. Specifically, competitive selection for faster growth resulted in dominance by fast-growing tree genotypes that were poorly defended against natural enemies. This outcome is a signature example of eco-evolutionary dynamics: Competitive interactions affected microevolutionary trajectories on a timescale relevant to subsequent ecological interactions [Brunner et al., Funct. Ecol. 33, 7-12 (2019)]. Eco-evolutionary drivers of tree growth and defense are thus critical to stand-level trait variation, which structures communities and ecosystems over expansive spatiotemporal scales.

Ecological factors influence balancing selection on leaf chemical profiles of a wildflower

Balancing selection is frequently invoked as a mechanism that maintains variation within and across populations. However, there are few examples of balancing selection operating on loci underpinning complex traits, which frequently display high levels of variation. We investigated mechanisms that may maintain variation in a focal polymorphism-leaf chemical profiles of a perennial wildflower (Boechera stricta, Brassicaceae)-explicitly interrogating multiple ecological and genetic processes including spatial variation in selection, antagonistic pleiotropy and frequency-dependent selection. A suite of common garden and greenhouse experiments showed that the alleles underlying variation in chemical profile have contrasting fitness effects across environments, implicating two ecological drivers of selection on chemical profile: herbivory and drought. Phenotype-environment associations and molecular genetic analyses revealed additional evidence of past selection by these drivers. Together, these data are consistent with balancing selection on chemical profile, probably caused by pleiotropic effects of secondary chemical biosynthesis genes on herbivore defence and drought response.

Ecological Interactions, Environmental Gradients, and Gene Flow in Local Adaptation

Despite long-standing interest in local adaptation of plants to their biotic and abiotic environment, existing theory, and many case studies, little work to date has addressed within-species evolution of concerted strategies and how these might contrast with patterns across species. Here we consider the interactions between pollinators, herbivores, and resource availability in shaping plant local adaptation, how these interactions impact plant phenotypes and gene flow, and the conditions where multiple traits align along major environmental gradients such as latitude and elevation. Continued work in emerging model systems will benefit from the melding of classic experimental approaches with novel population genetic analyses to reveal patterns and processes in plant local adaptation.

Adaptive and nonadaptive causes of heterogeneity in genetic differentiation across the Mimulus guttatus genome

Genetic diversity becomes structured among populations over time due to genetic drift and divergent selection. Although population structure is often treated as a uniform underlying factor, recent resequencing studies of wild populations have demonstrated that diversity in many regions of the genome may be structured quite dissimilar to the genome-wide pattern. Here, we explored the adaptive and nonadaptive causes of such genomic heterogeneity using population-level, whole genome resequencing data obtained from annual Mimulus guttatus individuals collected across a rugged environment landscape. We found substantial variation in how genetic differentiation is structured both within and between chromosomes, although, in contrast to other studies, known inversion polymorphisms appear to serve only minor roles in this heterogeneity. In addition, much of the genome can be clustered into eight among-population genetic differentiation patterns, but only two of these clusters are particularly consistent with patterns of isolation by distance. By performing genotype-environment association analysis, we also identified genomic intervals where local adaptation to specific climate factors has accentuated genetic differentiation among populations, and candidate genes in these windows indicate climate adaptation may proceed through changes affecting specialized metabolism, drought resistance, and development. Finally, by integrating our findings with previous studies, we show that multiple aspects of plant reproductive biology may be common targets of balancing selection and that variants historically involved in climate adaptation among populations have probably also fuelled rapid adaptation to microgeographic environmental variation within sites.

Larval habitats impose trait-dependent limits on the direction and rate of adult evolution in dragonflies

Natural selection on juveniles is often invoked as a constraint on adult evolution, but it remains unclear when such restrictions will have their greatest impact. Selection on juveniles could, for example, mainly limit the evolution of adult traits that mostly develop prior to maturity. Alternatively, selection on juveniles might primarily constrain the evolution of adult traits that experience weak or context-dependent selection in the adult stage. Using a comparative study of dragonflies, I tested these hypotheses by examining how a species' larval habitat was related to the evolution of two adult traits that differ in development and exposure to selection: adult size and male ornamentation. Whereas adult size is fixed at metamorphosis and experiences consistent positive selection in the adult stage, ornaments develop throughout adulthood and provide context-dependent fitness benefits. My results show that species that develop in less stable larval habitats have smaller adult sizes and slower rates of adult size evolution. However, these risky larval habitats do not limit ornament expression or rates of ornament evolution. Selection on juveniles may therefore primarily affect the evolution of adult traits that mostly develop prior to maturity.

Quantitative trait locus mapping reveals an independent genetic basis for joint divergence in leaf function, life-history, and floral traits between scarlet monkeyflower (Mimulus cardinalis) populations

PREMISE

Across taxa, vegetative and floral traits that vary along a fast-slow life-history axis are often correlated with leaf functional traits arrayed along the leaf economics spectrum, suggesting a constrained set of adaptive trait combinations. Such broad-scale convergence may arise from genetic constraints imposed by pleiotropy (or tight linkage) within species, or from natural selection alone. Understanding the genetic basis of trait syndromes and their components is key to distinguishing these alternatives and predicting evolution in novel environments.

METHODS

We used a line-cross approach and quantitative trait locus (QTL) mapping to characterize the genetic basis of twenty leaf functional/physiological, life history, and floral traits in hybrids between annualized and perennial populations of scarlet monkeyflower (Mimulus cardinalis).

RESULTS

We mapped both single and multi-trait QTLs for life history, leaf function and reproductive traits, but found no evidence of genetic co-ordination across categories. A major QTL for three leaf functional traits (thickness, photosynthetic rate, and stomatal resistance) suggests that a simple shift in leaf anatomy may be key to adaptation to seasonally dry habitats.

CONCLUSIONS

Our results suggest that the co-ordination of resource-acquisitive leaf physiological traits with a fast life-history and more selfing mating system results from environmental selection rather than functional or genetic constraint. Independent assortment of distinct trait modules, as well as a simple genetic basis to leaf physiological traits associated with drought escape, may facilitate adaptation to changing climates.