Is local adaptation in Mimulus guttatus caused by trade-offs at individual loci?

Ocean exposure and latitude drive multiple clines within the coastal perennial ecotype of the yellow monkeyflower, Mimulus guttatus

PREMISE

A key goal of evolutionary biologists is to understand how and why genetic variation is partitioned within species. In the yellow monkeyflower, Mimulus guttatus (syn. Erythranthe guttata), coastal perennial populations constitute a single genetically and morphologically differentiated ecotype compared to inland M. guttatus populations. While the coastal ecotype's distinctiveness has now been well documented, there is also environmental variation across the ecotype's range that could drive more continuous differentiation among its component populations.

METHODS

Based on previous observations of a potential cline within this ecotype, we quantified plant height, among other traits, across coastal perennial accessions from 74 populations in a greenhouse common garden experiment. To evaluate potential drivers of the relationship between trait variation and latitude, we regressed height against multiple climatic factors, including temperature, precipitation, and coastal wind speeds. We also accounted for exposure to the open ocean in all analyses.

RESULTS

Multiple traits were correlated with latitude of origin, but none more than plant height. Height was negatively correlated with latitude, and plants directly exposed to the open ocean were shorter than those protected from coastal winds. Further analyses revealed that height was correlated with climatic factors (precipitation, temperature, and wind speeds) that were autocorrelated with latitude. We hypothesize that one or more of these climatic factors drove the evolution of latitudinal clinal variation within the coastal ecotype.

CONCLUSIONS

Overall, our study illustrates the complexity of how the distribution of environmental variation can simultaneously drive the evolution of both distinct ecotypes and continuous clines within those ecotypes.

Conditionally Deleterious Mutation Load Accumulates in Genomic Islands of Local Adaptation but Can Be Purged with Sufficient Genotypic Redundancy

AbstractLocal adaptation frequently evolves in patches or environments that are connected via migration. In these cases, genomic regions that are linked to a locally adapted locus experience reduced effective migration rates. Via individual-based simulations of a two-patch system, we show that this reduced effective migration results in the accumulation of conditionally deleterious mutations, but not universally deleterious mutations, adjacent to adaptive loci. When there is redundancy in the genetic basis of local adaptation (i.e., genotypic redundancy), turnover of locally adapted polymorphisms allows conditionally deleterious mutation load to be purged. The amount of mutational load that accumulates adjacent to locally adapted loci is dependent on redundancy, recombination rate, migration rate, population size, strength of selection, and the phenotypic effect size of adaptive alleles. Our results highlight the need to be cautious when interpreting patterns of local adaptation at the level of phenotype or fitness, as the genetic basis of local adaptation can be transient, and evolution may confer a degree of maladaptation to nonlocal environments.

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation

Identifying the genetic basis of local adaptation and fitness trade-offs across environments is a central goal of evolutionary biology. Cold acclimation is an adaptive plastic response for surviving seasonal freezing, and costs of acclimation may be a general mechanism for fitness trade-offs across environments in temperate zone species. Starting with locally adapted ecotypes of Arabidopsis thaliana from Italy and Sweden, we examined the fitness consequences of a naturally occurring functional polymorphism in CBF2. This gene encodes a transcription factor that is a major regulator of cold-acclimated freezing tolerance and resides within a locus responsible for a genetic trade-off for long-term mean fitness. We estimated the consequences of alternate genotypes of CBF2 on 5-y mean fitness and fitness components at the native field sites by comparing near-isogenic lines with alternate genotypes of CBF2 to their genetic background ecotypes. The effects of CBF2 were validated at the nucleotide level using gene-edited lines in the native genetic backgrounds grown in simulated parental environments. The foreign CBF2 genotype in the local genetic background reduced long-term mean fitness in Sweden by more than 10%, primarily via effects on survival. In Italy, fitness was reduced by more than 20%, primarily via effects on fecundity. At both sites, the effects were temporally variable and much stronger in some years. The gene-edited lines confirmed that CBF2 encodes the causal variant underlying this genetic trade-off. Additionally, we demonstrated a substantial fitness cost of cold acclimation, which has broad implications for potential maladaptive responses to climate change.

Ecological genetics of local adaptation in Arabidopsis: An 8-year field experiment

There is considerable evidence for local adaptation in nature, yet important questions remain regarding its genetic basis. How many loci are involved? What are their effect sizes? What is the relative importance of conditional neutrality versus genetic trade-offs? Here we address these questions in the self-pollinating, annual plant Arabidopsis thaliana. We used 400 recombinant inbred lines (RILs) derived from two locally adapted populations in Italy and Sweden, grew the RILs and parents at the parental locations, and mapped quantitative trait loci (QTL) for mean fitness (fruits/seedling planted). We previously published results from the first 3 years of the study, and here add five additional years, providing a unique opportunity to assess how temporal variation in selection might affect QTL detection and classification. We found 10 adaptive and one maladaptive QTL in Italy, and six adaptive and four maladaptive QTL in Sweden. The discovery of maladaptive QTL at both sites suggests that even locally adapted populations are not always at their genotypic optimum. Mean effect sizes for adaptive QTL, 0.97 and 0.55 fruits in Italy and Sweden, respectively, were large relative to the mean fitness of the RILs (approximately 8 fruits/seedling planted at both sites). Both genetic trade-offs (four cases) and conditional neutrality (seven cases) contribute to local adaptation in this system. The 8-year dataset provided greater power to detect QTL and to estimate their locations compared to our previous 3-year study, identifying one new genetic trade-off and resolving one genetic trade-off into two conditionally adaptive QTL.

Local adaptation: Causal agents of selection and adaptive trait divergence

Divergent selection across the landscape can favor the evolution of local adaptation in populations experiencing contrasting conditions. Local adaptation is widely observed in a diversity of taxa, yet we have a surprisingly limited understanding of the mechanisms that give rise to it. For instance, few have experimentally confirmed the biotic and abiotic variables that promote local adaptation, and fewer yet have identified the phenotypic targets of selection that mediate local adaptation. Here, we highlight critical gaps in our understanding of the process of local adaptation and discuss insights emerging from in-depth investigations of the agents of selection that drive local adaptation, the phenotypes they target, and the genetic basis of these phenotypes. We review historical and contemporary methods for assessing local adaptation, explore whether local adaptation manifests differently across life history, and evaluate constraints on local adaptation.

The Boechera model system for evolutionary ecology

Model systems in biology expand the research capacity of individuals and the community. Closely related to Arabidopsis, the genus Boechera has emerged as an important ecological model owing to the ability to integrate across molecular, functional, and eco-evolutionary approaches. Boechera species are broadly distributed in relatively undisturbed habitats predominantly in western North America and provide one of the few experimental systems for identification of ecologically important genes through genome-wide association studies and investigations of selection with plants in their native habitats. The ecologically, evolutionarily, and agriculturally important trait of apomixis (asexual reproduction via seeds) is common in the genus, and field experiments suggest that abiotic and biotic environments shape the evolution of sex. To date, population genetic studies have focused on the widespread species B. stricta, detailing population divergence and demographic history. Molecular and ecological studies show that balancing selection maintains genetic variation in ~10% of the genome, and ecological trade-offs contribute to complex trait variation for herbivore resistance, flowering phenology, and drought tolerance. Microbiome analyses have shown that host genotypes influence leaf and root microbiome composition, and the soil microbiome influences flowering phenology and natural selection. Furthermore, Boechera offers numerous opportunities for investigating biological responses to global change. In B. stricta, climate change has induced a shift of >2 weeks in the timing of first flowering since the 1970s, altered patterns of natural selection, generated maladaptation in previously locally-adapted populations, and disrupted life history trade-offs. Here we review resources and results for this eco-evolutionary model system and discuss future research directions.

A Pleiotropic Flowering Time QTL Exhibits Gene-by-Environment Interaction for Fitness in a Perennial Grass

Appropriate flowering time is a crucial adaptation impacting fitness in natural plant populations. Although the genetic basis of flowering variation has been extensively studied, its mechanisms in nonmodel organisms and its adaptive value in the field are still poorly understood. Here, we report new insights into the genetic basis of flowering time and its effect on fitness in Panicum hallii, a native perennial grass. Genetic mapping in populations derived from inland and coastal ecotypes identified flowering time quantitative trait loci (QTL) and many exhibited extensive QTL-by-environment interactions. Patterns of segregation within recombinant hybrids provide strong support for directional selection driving ecotypic divergence in flowering time. A major QTL on chromosome 5 (q-FT5) was detected in all experiments. Fine-mapping and expression studies identified a gene with orthology to a rice FLOWERING LOCUS T-like 9 (PhFTL9) as the candidate underlying q-FT5. We used a reciprocal transplant experiment to test for local adaptation and the specific impact of q-FT5 on performance. We did not observe local adaptation in terms of fitness tradeoffs when contrasting ecotypes in home versus away habitats. However, we observed that the coastal allele of q-FT5 conferred a fitness advantage only in its local habitat but not at the inland site. Sequence analyses identified an excess of low-frequency polymorphisms at the PhFTL9 promoter in the inland lineage, suggesting a role for either selection or population expansion on promoter evolution. Together, our findings demonstrate the genetic basis of flowering variation in a perennial grass and provide evidence for conditional neutrality underlying flowering time divergence.

A polygenic architecture with habitat-dependent effects underlies ecological differentiation in Silene

Ecological differentiation can drive speciation but it is unclear how the genetic architecture of habitat-dependent fitness contributes to lineage divergence. We investigated the genetic architecture of cumulative flowering, a fitness component, in second-generation hybrids between Silene dioica and Silene latifolia transplanted into the natural habitat of each species. We used reduced-representation sequencing and Bayesian sparse linear mixed models (BSLMMs) to analyze the genetic control of cumulative flowering in each habitat. Our results point to a polygenic architecture of cumulative flowering. Allelic effects were mostly beneficial or deleterious in one habitat and neutral in the other. Positive-effect alleles often were derived from the native species, whereas negative-effect alleles, at other loci, tended to originate from the non-native species. We conclude that ecological differentiation is governed and maintained by many loci with small, habitat-dependent effects consistent with conditional neutrality. This pattern may result from differences in selection targets in the two habitats and from environmentally dependent deleterious load. Our results further suggest that selection for native alleles and against non-native alleles acts as a barrier to gene flow between species.

Three problems in the genetics of speciation by selection

Speciation is the process by which barriers to gene flow evolve between populations. Although we now know that speciation is largely driven by natural selection, knowledge of the agents of selection and the genetic and genomic mechanisms that facilitate divergence is required for a satisfactory theory of speciation. In this essay, we highlight three advances/problems in our understanding of speciation that have arisen from studies of the genes and genomic regions that underlie the evolution of reproductive isolation. First, we describe how the identification of “speciation” genes makes it possible to identify the agents of selection causing the evolution of reproductive isolation, while also noting that the link between the genetics of phenotypic divergence and intrinsic postzygotic reproductive barriers remains tenuous. Second, we discuss the important role of recombination suppressors in facilitating speciation with gene flow, but point out that the means and timing by which reproductive barriers become associated with recombination cold spots remains uncertain. Third, we establish the importance of ancient genetic variation in speciation, although we argue that the focus of speciation studies on evolutionarily young groups may bias conclusions in favor of ancient variation relative to new mutations.

Phylogenomic analyses in Phrymaceae reveal extensive gene tree discordance in relationships among major clades

PREMISE

Phylogenomic datasets using genomes and transcriptomes provide rich opportunities beyond resolving bifurcating phylogenetic relationships. Monkeyflower (Phrymaceae) is a model system for evolutionary ecology. However, it lacks a well-supported phylogeny as a basis for a stable taxonomy and for macroevolutionary comparisons.

METHODS

We sampled 24 genomes and transcriptomes in Phrymaceae and closely related families, including eight newly sequenced transcriptomes. We reconstructed the phylogeny using IQ-TREE and ASTRAL, evaluated gene tree discordance using PhyParts, Quartet Sampling, and a cloudogram, and carried out reticulation analyses using PhyloNet and HyDe. We searched for whole genome duplication (WGD) events using chromosome numbers, synonymous distances, and gene duplication events as evidence.

RESULTS

Most gene trees support the monophyly of Phrymaceae and each of its tribes. Most gene trees also support tribe Mimuleae being sister to Phrymeae + Diplaceae + Leucocarpeae, with extensive gene tree discordance among the latter three. Despite the discordance, the monophyly of Mimulus s.l. is rejected, and no individual reticulation event among the Phrymaceae tribes is well-supported. Reticulation likely occurred among Erythranthe bicolor and closely related species. No ancient WGD was detected in Phrymaceae. Instead, small-scale duplications are among potential drivers of macroevolutionary diversification of Phrymaceae.

CONCLUSIONS

We show that analysis of reticulate evolution is sensitive to taxon sampling and methods used. We also demonstrate that phylogenomic datasets using genomes and transcriptomes present rich opportunities to investigate gene family evolution and genome duplication events involved in lineage diversification and adaptation.

Multivariate selection mediated by aridity predicts divergence of drought-resistant traits along natural aridity gradients of an invasive weed

Geographical variation in the environment underpins selection for local adaptation and evolutionary divergence among populations. Because many environmental conditions vary across species' ranges, identifying the specific environmental variables underlying local adaptation is profoundly challenging. We tested whether natural selection mediated by aridity predicts clinal divergence among invasive populations of capeweed (Arctotheca calendula) that established and spread across southern Australia during the last two centuries. Using common garden experiments with two environmental treatments (wet and dry) that mimic aridity conditions across capeweed's invasive range, we estimated clinal divergence and effects of aridity on fitness and multivariate phenotypic selection in populations sampled along aridity gradients in Australia. We show that: (1) capeweed populations have relatively high fitness in aridity environments similar to their sampling locations; (2) the magnitude and direction of selection strongly differs between wet and dry treatments, with drought stress increasing the strength of selection; and (3) differences in directional selection between wet and dry treatments predict patterns of clinal divergence across the aridity gradient, particularly for traits affecting biomass, flowering phenology and putative antioxidant expression. Our results suggest that aridity-mediated selection contributes to trait diversification among invasive capeweed populations, possibly facilitating the expansion of capeweed across southern Australia.

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.

A B73×Palomero Toluqueño mapping population reveals local adaptation in Mexican highland maize

Generations of farmer selection in the central Mexican highlands have produced unique maize varieties adapted to the challenges of the local environment. In addition to possessing great agronomic and cultural value, Mexican highland maize represents a good system for the study of local adaptation and acquisition of adaptive phenotypes under cultivation. In this study, we characterize a recombinant inbred line population derived from the B73 reference line and the Mexican highland maize variety Palomero Toluqueño. B73 and Palomero Toluqueño showed classic rank-changing differences in performance between lowland and highland field sites, indicative of local adaptation. Quantitative trait mapping identified genomic regions linked to effects on yield components that were conditionally expressed depending on the environment. For the principal genomic regions associated with ear weight and total kernel number, the Palomero Toluqueño allele conferred an advantage specifically in the highland site, consistent with local adaptation. We identified Palomero Toluqueño alleles associated with expression of characteristic highland traits, including reduced tassel branching, increased sheath pigmentation and the presence of sheath macrohairs. The oligogenic architecture of these three morphological traits supports their role in adaptation, suggesting they have arisen from consistent directional selection acting at distinct points across the genome. We discuss these results in the context of the origin of phenotypic novelty during selection, commenting on the role of de novo mutation and the acquisition of adaptive variation by gene flow from endemic wild relatives.

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.

Genetic Divergence between Two Sympatric Ecotypes of Phalaenopsis pulcherrima on Hainan Island

Ecotypes are the result of ecological differentiation at the early stages of speciation. Adaptation to soil conditions offers arguably the best examples of local adaptation in plants. Two sympatric ecotypes, with either a red or green abaxial leaf surface, were found without clear geographical isolation in Phalaenopsis pulcherrima, a Southeast Asia endemic and endangered orchid. The soil of the red leaf ecotype has a higher water content and nutrient content than the green ecotype. What is the genetic structure of the two ecotypes? Is there complete or partial reproductive isolation between the two ecotypes? In this work, leaf reflection of the two ecotypes in P. pulcherrima were compared, to illustrate their difference in leaf color. The genetic differentiation between two ecotypes was examined, using ISSR and SRAP markers to determine the genetic structure of the populations. Our results showed that the green ecotype had reflectance spectrum peaks at 530 nm and 620 nm, while in the red ecotype, the peak at 530 nm was absent. A total of 165 ISSR and SRAP loci showed a high level of genetic diversity within the green ecotype, and analyses of the population structure revealed two genetic clusters that corresponded to the red and green ecotypes. The percentage of variation between the two ecotypes (24.55%) was greater than the percentage of variation among the populations (16.54%)—indicating partial reproductive isolation, high genetic differentiation, and that ecological differentiation has been more important than geographical barriers among populations within ecotypes. Most pairwise FST values between the populations within either ecotype on Hainan Island were less than 0.15; however, the FST between both the Thai and Malaysian populations and the Hainan Island population was greater than 0.25, due to South China sea isolation. Ecotypic differentiation is an important part of speciation; therefore, we must take into account the axes along which lineages sort, when formulating protection strategies.

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.

Rapid Parallel Adaptation to Anthropogenic Heavy Metal Pollution

The impact of human-mediated environmental change on the evolutionary trajectories of wild organisms is poorly understood. In particular, capacity of species to adapt rapidly (in hundreds of generations or less), reproducibly and predictably to extreme environmental change is unclear. Silene uniflora is predominantly a coastal species, but it has also colonized isolated, disused mines with phytotoxic, zinc-contaminated soils. To test whether rapid, parallel adaptation to anthropogenic pollution has taken place, we used reduced representation sequencing (ddRAD) to reconstruct the evolutionary history of geographically proximate mine and coastal population pairs and found largely independent colonization of mines from different coastal sites. Furthermore, our results show that parallel evolution of zinc tolerance has occurred without gene flow spreading adaptive alleles between mine populations. In genomic regions where signatures of selection were detected across multiple mine-coast pairs, we identified genes with functions linked to physiological differences between the putative ecotypes, although genetic differentiation at specific loci is only partially shared between mine populations. Our results are consistent with a complex, polygenic genetic architecture underpinning rapid adaptation. This shows that even under a scenario of strong selection and rapid adaptation, evolutionary responses to human activities (and other environmental challenges) may be idiosyncratic at the genetic level and, therefore, difficult to predict from genomic data.

Life-history trade-offs and the genetic basis of fitness in Arabidopsis thaliana

Resources allocated to survival cannot be used to increase fecundity, but the extent to which this trade-off constrains adaptation depends on overall resource status. Adaptation to local environmental conditions may therefore entail the evolution of traits that increase the amount of resources available to individuals (their resource status or 'condition'). We examined the relative contribution of trade-offs and increased condition to adaptive evolution in a recombinant inbred line population of Arabidopsis thaliana planted at the native sites of the parental ecotypes in Italy and Sweden in 2 years. We estimated genetic correlations among fitness components based on genotypic means and explored their causes with QTL mapping. The local ecotype produced more seeds per fruit than did the non-local ecotype, reflected in stronger adaptive differentiation than was previously shown based on survival and fruit number only. Genetic correlations between survival and overall fecundity, and between number of fruits and number of seeds per fruit, were positive, and there was little evidence of a trade-off between seed size and number. Quantitative trait loci for these traits tended to map to the same regions of the genome and showed positive pleiotropic effects. The results indicate that adaptive differentiation between the two focal populations largely reflects the evolution of increased ability to acquire resources in the local environment, rather than shifts in the relative allocation to different life-history traits. Differentiation both in phenology and in tolerance to cold is likely to contribute to the advantage of the local genotype at the two sites.

The genetic architecture and evolution of life-history divergence among perennials in the Mimulus guttatus species complex

Ecological divergence is a fundamental source of phenotypic diversity between closely related species, yet the genetic architecture of most ecologically relevant traits is poorly understood. Differences in elevation can impose substantial divergent selection on both complex, correlated suites of traits (such as life-history), as well as novel adaptations. We use the Mimulus guttatus species complex to assess if the divergence in elevation is accompanied by trait divergence in a group of closely related perennials and determine the genetic architecture of this divergence. We find that divergence in elevation is associated with differences in life-history, as well as a unique trait, the production of rhizomes. The divergence between two perennials is largely explained by few mid-to-large effect quantitative trait loci (QTLs). However, the presence of QTLs with correlated, but opposing effects on multiple traits leads to some hybrids with transgressive trait combinations. Lastly, we find that the genetic architecture of the ability to produce rhizomes changes through development, wherein most hybrids produce rhizomes, but only later in development. Our results suggest that elevational differences may shape life-history divergence between perennials, but aspects of the genetic architecture of divergence may have implications for hybrid fitness in nature.

Asymmetric, but opposing reductions in immigrant viability and fecundity promote reproductive isolation among host-associated populations of an insect herbivore

Immigrant inviability can contribute to reproductive isolation (RI) during ecological speciation by reducing the survival of immigrants in non-native environments. However, studies that assess the fitness consequence of immigrants moving from native to non-native environments typically fail to explore the potential role of concomitant reductions in immigrant fecundity despite recent evidence suggesting its prominent role during local adaptation. Here, we evaluate the directionality and magnitude of both immigrant viability and fecundity to RI in a host-specific gall-forming wasp, Belonocnema treatae. Using reciprocal transplant experiments replicated across sites, we measure immigrant viability and fecundity by comparing differences in the incidence of gall formation (viability) and predicted the number of eggs per female (fecundity) between residents and immigrants in each of two host-plant environments. Reduced immigrant viability was found in one host environment while reduced immigrant fecundity was found in the other. Such habitat-dependent barriers resulted in asymmetric RI between populations. By surveying recent literature on local adaptation, we find that asymmetry in immigrant viability and fecundity are widespread across disparate taxa, which highlights the need to combine estimates of both common and overlooked barriers in cases of potential bi-directional gene flow to create a more comprehensive view of the evolution of RI.

Gene regulatory effects of a large chromosomal inversion in highland maize

Chromosomal inversions play an important role in local adaptation. Inversions can capture multiple locally adaptive functional variants in a linked block by repressing recombination. However, this recombination suppression makes it difficult to identify the genetic mechanisms underlying an inversion's role in adaptation. In this study, we used large-scale transcriptomic data to dissect the functional importance of a 13 Mb inversion locus (Inv4m) found almost exclusively in highland populations of maize (Zea mays ssp. mays). Inv4m was introgressed into highland maize from the wild relative Zea mays ssp. mexicana, also present in the highlands of Mexico, and is thought to be important for the adaptation of these populations to cultivation in highland environments. However, the specific genetic variants and traits that underlie this adaptation are not known. We created two families segregating for the standard and inverted haplotypes of Inv4m in a common genetic background and measured gene expression effects associated with the inversion across 9 tissues in two experimental conditions. With these data, we quantified both the global transcriptomic effects of the highland Inv4m haplotype, and the local cis-regulatory variation present within the locus. We found diverse physiological effects of Inv4m across the 9 tissues, including a strong effect on the expression of genes involved in photosynthesis and chloroplast physiology. Although we could not confidently identify the causal alleles within Inv4m, this research accelerates progress towards understanding this inversion and will guide future research on these important genomic features.

Life-History Plasticity and Water-Use Trade-Offs Associated with Drought Resistance in a Clade of California Jewelflowers

Water limitation is a primary driver of plant geographic distributions and individual plant fitness. Drought resistance is the ability to survive and reproduce despite limited water, and numerous studies have explored its physiological basis in plants. However, it is unclear how drought resistance and trade-offs associated with drought resistance evolve within plant clades. We quantified the relationship between water availability and fitness for 13 short-lived plant taxa in the Streptanthus clade that vary in their phenology and the availability of water in the environments where they occur. We derived two parameters from these relationships: plant fitness when water is not limiting and the water inflection point (WIF), the watering level at which additional water is most efficiently turned into fitness. We used phylogenetic comparative methods to explore trade-offs related to drought resistance and trait plasticity and the degree to which water relationship parameters are conserved. Taxa from drier climates produced fruits at the lowest water levels, had a lower WIF, flowered earlier, had shorter life spans, had greater plastic water-use efficiency (WUE), and had lower fitness at nonlimiting water. In contrast, later-flowering Streptanthus taxa from less xeric climates experienced high fitness at nonlimiting water but had no fitness at the lowest water levels. Across the clade, we found a trade-off between drought resistance and fitness at high water, though a single ruderal species was an outlier in this relationship. Our results suggest that drought escape trades off with maximal fitness under nonlimiting water, and both are tied to phenology. We also found that variation in trait plasticity determines how different plant species produce fitness over a water gradient.

Life history evolution, species differences, and phenotypic plasticity in hemiparasitic eyebrights (Euphrasia)

PREMISE

Species delimitation in parasitic organisms is challenging because traits used to identify species are often plastic and vary depending on the host. Here, we use species from a recent radiation of generalist hemiparasitic Euphrasia to investigate trait variation and trait plasticity. We tested whether Euphrasia species show reliable trait differences, investigated whether these differences correspond to life history trade-offs between growth and reproduction, and quantified plasticity in response to host species.

METHODS

Common garden experiments were used to evaluate trait differences between 11 Euphrasia taxa grown on a common host, document phenotypic plasticity when a single Euphrasia species is grown on eight different hosts, and relate observations to trait differences recorded in the wild.

RESULTS

Euphrasia exhibited variation in life history strategies; some individuals transitioned rapidly to flowering at the expense of early season growth, while others invested in vegetative growth and delayed flowering. Life history differences were present between some species, though many related taxa lacked clear trait differences. Species differences were further blurred by phenotypic plasticity-many traits were plastic and changed with host type or between environments.

CONCLUSIONS

Phenotypic plasticity in response to host and environment confounds species delimitation in Euphrasia. When grown in a common garden environment, some morphologically distinct taxa can be identified, though others represent morphologically similar shallow segregates. Trait differences present between some species and populations demonstrate the rapid evolution of distinct life history strategies in response to local ecological conditions.

Contrasting environmental factors drive local adaptation at opposite ends of an environmental gradient in the yellow monkeyflower (Mimulus guttatus)

PREMISE

Identifying the environmental factors responsible for natural selection across different habitats is crucial for understanding the process of local adaptation in plants. Despite its importance, few studies have successfully isolated the environmental factors driving local adaptation in nature. In this study, we evaluated the agents of selection responsible for local adaptation of the monkeyflower Mimulus guttatus to California's coastal and inland habitats.

METHODS

We implemented a manipulative reciprocal transplant experiment at coastal and inland sites, where we excluded aboveground stressors in an effort to elucidate their role in the evolution of local adaptation.

RESULTS

Excluding aboveground stressors, most likely a combination of salt spray and herbivory, completely rescued inland annual plant fitness when transplanted to coastal habitat. The exclosures in inland habitat provided a benefit to the performance of coastal perennial plants. However, the exclosures are unlikely to provide much fitness benefit to the coastal plants at the inland site because of their general inability to flower in time to escape from the summer drought.

CONCLUSIONS

Our study demonstrates that a distinct set of selective agents (aboveground vs. belowground) are responsible for local adaptation at opposite ends of an environmental gradient.

Genetic and physiological mechanisms of freezing tolerance in locally adapted populations of a winter annual

PREMISE

Despite myriad examples of local adaptation, the phenotypes and genetic variants underlying such adaptive differentiation are seldom known. Recent work on freezing tolerance and local adaptation in ecotypes of Arabidopsis thaliana from Italy and Sweden provides an essential foundation for uncovering the genotype-phenotype-fitness map for an adaptive response to a key environmental stress.

METHODS

We examined the consequences of a naturally occurring loss-of-function (LOF) mutation in an Italian allele of the gene that encodes the transcription factor CBF2, which underlies a major freezing-tolerance locus. We used four lines with a Swedish genetic background, each containing a LOF CBF2 allele. Two lines had introgression segments containing the Italian CBF2 allele, and two contained deletions created using CRISPR-Cas9. We used a growth chamber experiment to quantify freezing tolerance and gene expression before and after cold acclimation.

RESULTS

Freezing tolerance was lower in the Italian (11%) compared to the Swedish (72%) ecotype, and all four experimental CBF2 LOF lines had reduced freezing tolerance compared to the Swedish ecotype. Differential expression analyses identified 10 genes for which all CBF2 LOF lines, and the IT ecotype had similar patterns of reduced cold responsive expression compared to the SW ecotype.

CONCLUSIONS

We identified 10 genes that are at least partially regulated by CBF2 that may contribute to the differences in cold-acclimated freezing tolerance between the Italian and Swedish ecotypes. These results provide novel insight into the molecular and physiological mechanisms connecting a naturally occurring sequence polymorphism to an adaptive response to freezing conditions.

QTL × environment interactions underlie adaptive divergence in switchgrass across a large latitudinal gradient

Understanding how individual genetic loci contribute to trait variation across geographic space is of fundamental importance for understanding evolutionary adaptations. Our study demonstrates that most loci underlying locally adaptive trait variation have beneficial effects in some geographic regions while conferring little or no detectable cost in other parts of the geographic range of switchgrass over two field seasons of study. Thus, loci that contribute to local adaptation vary in the degree to which they are costly in alternative environments but typically confer greater benefits than costs. Further, our study suggests that breeding locally adapted varieties of switchgrass will be a boon to the biofuel industry, as locally adaptive loci could be combined to increase local yields in switchgrass.

Local adaptation is the process by which natural selection drives adaptive phenotypic divergence across environmental gradients. Theory suggests that local adaptation results from genetic trade-offs at individual genetic loci, where adaptation to one set of environmental conditions results in a cost to fitness in alternative environments. However, the degree to which there are costs associated with local adaptation is poorly understood because most of these experiments rely on two-site reciprocal transplant experiments. Here, we quantify the benefits and costs of locally adaptive loci across 17° of latitude in a four-grandparent outbred mapping population in outcrossing switchgrass (Panicum virgatum L.), an emerging biofuel crop and dominant tallgrass species. We conducted quantitative trait locus (QTL) mapping across 10 sites, ranging from Texas to South Dakota. This analysis revealed that beneficial biomass (fitness) QTL generally incur minimal costs when transplanted to other field sites distributed over a large climatic gradient over the 2 y of our study. Therefore, locally advantageous alleles could potentially be combined across multiple loci through breeding to create high-yielding regionally adapted cultivars.

Unpacking Conditional Neutrality: Genomic Signatures of Selection on Conditionally Beneficial and Conditionally Deleterious Mutations

It is common to look for signatures of local adaptation in genomes by identifying loci with extreme levels of allele frequency divergence among populations. This approach to finding genes associated with local adaptation often assumes antagonistic pleiotropy, wherein alternative alleles are strongly favored in alternative environments. Conditional neutrality has been proposed as an alternative to antagonistic pleiotropy, but conditionally neutral polymorphisms are transient, and it is unclear how much outlier signal would be maintained under different forms of conditional neutrality. Here, we use individual-based simulations and a simple analytical heuristic to show that a pattern that mimics local adaptation at the phenotypic level, where each genotype has the highest fitness in its home environment, can be produced by the accumulation of mutations that are neutral in their home environment and deleterious in nonlocal environments. Because conditionally deleterious mutations likely arise at a rate many times higher than conditionally beneficial mutations, they can have a significant cumulative effect on fitness even when individual effect sizes are small. We show that conditionally deleterious mutations driving nonlocal maladaptation may be undetectable by even the most powerful genome scans, as differences in allele frequency between populations are typically small. We also explore the evolutionary effects of conditionally beneficial mutations and find that they can maintain significant signals of local adaptation, and they would be more readily detectable than conditionally deleterious mutations using conventional genome scan approaches. We discuss implications for interpreting outcomes of transplant experiments and genome scans that are used to study the genetic basis of local adaptation.

Gene regulatory divergence between locally adapted ecotypes in their native habitats

Local adaptation is a key driver of ecological specialization and the formation of new species. Despite its importance, the evolution of gene regulatory divergence among locally adapted populations is poorly understood, especially how that divergence manifests in nature. Here, we evaluate gene expression divergence and allele-specific gene expression responses for locally adapted coastal perennial and inland annual accessions of the yellow monkeyflower, Mimulus guttatus, in a field reciprocal transplant experiment. Overall, 6765 (73%) of surveyed genes were differentially expressed between coastal and inland habitats, while 7213 (77%) were differentially expressed between the coastal perennial and inland annual accessions. Cis-regulatory variation was pervasive, affecting 79% (5532) of differentially expressed genes. We detected trans effects for 52% (3611) of differentially expressed genes. Expression plasticity of alleles across habitats (G × E interactions) appears to be relatively common (affecting 18% of transcripts) and could minimize fitness trade-offs at loci that contribute to local adaptation. We also found evidence that at least one chromosomal inversion may act as supergene by holding together haplotypes of differentially expressed genes, but this pattern depends on habitat context. Our results highlight multiple key patterns regarding the relationship between gene expression and the evolution of locally adapted populations.

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less-fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co-occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre-date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.

Large-effect mutations generate trade-off between predatory and locomotor ability during arms race coevolution with deadly prey

Adaptive evolution in response to one selective challenge may disrupt other important aspects of performance. Such evolutionary trade‐offs are predicted to arise in the process of local adaptation, but it is unclear if these phenotypic compromises result from the antagonistic effects of simple amino acid substitutions. We tested for trade‐offs associated with beneficial mutations that confer tetrodotoxin (TTX) resistance in the voltage‐gated sodium channel (Na_V1.4) in skeletal muscle of the common garter snake (Thamnophis sirtalis). Separate lineages in California and the Pacific Northwest independently evolved TTX‐resistant changes to the pore of Na_V1.4 as a result of arms race coevolution with toxic prey, newts of the genus Taricha. Snakes from the California lineage that were homozygous for an allele known to confer large increases in toxin resistance (Na_V1.4^LVNV) had significantly reduced crawl speed compared to individuals with the ancestral TTX‐sensitive channel. Heterologous expression of native snake Na_V1.4 proteins demonstrated that the same Na_V1.4^LVNV allele confers a dramatic increase in TTX resistance and a correlated decrease in overall channel excitability. Our results suggest the same mutations that accumulate during arms race coevolution and beneficially interfere with toxin‐binding also cause changes in electrophysiological function of the channel that may affect organismal performance. This trade‐off was only evident in the predator lineage where coevolution has led to the most extreme resistance phenotype, determined by four critical amino acid substitutions. If these biophysical changes also translate to a fitness cost—for example, through the inability of T. sirtalis to quickly escape predators—then pleiotropy at this single locus could contribute to observed variation in levels of TTX resistance across the mosaic landscape of coevolution.

Environmentally induced development costs underlie fitness tradeoffs

Local adaptation can lead to genotype-by-environment interactions, which can create fitness tradeoffs in alternative environments, and govern the distribution of biodiversity across geographic landscapes. Exploring the ecological circumstances that promote the evolution of fitness tradeoffs requires identifying how natural selection operates and during which ontogenetic stages natural selection is strongest. When organisms disperse to areas outside their natural range, tradeoffs might emerge when organisms struggle to reach key life history stages, or alternatively, die shortly after reaching life history stages if there are greater risks of mortality associated with costs to developing in novel environments. We used multiple populations from four ecotypes of an Australian native wildflower (Senecio pinnatifolius) in reciprocal transplants to explore how fitness tradeoffs arise across ontogeny. We then assessed whether the survival probability for plants from native and foreign populations was contingent on reaching key developmental stages. We found that fitness tradeoffs emerged as ontogeny progressed when native plants were more successful than foreign plants at reaching seedling establishment and maturity. Native and foreign plants that failed to reach seedling establishment died at the same rate, but plants from foreign populations died quicker than native plants after reaching seedling establishment, and died quicker regardless of whether they reached sexual maturity or not. Development rates were similar for native and foreign populations, but changed depending on the environment. Together, our results suggest that natural selection for environment-specific traits early in life history created tradeoffs between contrasting environments. Plants from foreign populations were either unable to develop to seedling establishment, or they suffered increased mortality as a consequence of reaching seedling establishment. The observation of tradeoffs together with environmentally dependent changes in development rate suggest that foreign environments induce organisms to develop at a rate different from their native habitat, incurring consequences for lifetime fitness and population divergence.

Trait differentiation and adaptation of plants along elevation gradients

Studies of genetic adaptation in plant populations along elevation gradients in mountains have a long history, but there has until now been neither a synthesis of how frequently plant populations exhibit adaptation to elevation nor an evaluation of how consistent underlying trait differences across species are. We reviewed studies of adaptation along elevation gradients (i) from a meta-analysis of phenotypic differentiation of three traits (height, biomass and phenology) from plants growing in 70 common garden experiments; (ii) by testing elevation adaptation using three fitness proxies (survival, reproductive output and biomass) from 14 reciprocal transplant experiments; (iii) by qualitatively assessing information at the molecular level, from 10 genomewide surveys and candidate gene approaches. We found that plants originating from high elevations were generally shorter and produced less biomass, but phenology did not vary consistently. We found significant evidence for elevation adaptation in terms of survival and biomass, but not for reproductive output. Variation in phenotypic and fitness responses to elevation across species was not related to life history traits or to environmental conditions. Molecular studies, which have focussed mainly on loci related to plant physiology and phenology, also provide evidence for adaptation along elevation gradients. Together, these studies indicate that genetically based trait differentiation and adaptation to elevation are widespread in plants. We conclude that a better understanding of the mechanisms underlying adaptation, not only to elevation but also to environmental change, will require more studies combining the ecological and molecular approaches.

Integrative Regulation of Drought Escape through ABA-Dependent and -Independent Pathways in Rice

Many plants have evolved a drought escape (DE) mechanism to shorten their life cycle when facing water-deficit conditions. While drought tolerance has been intensively investigated, the genetic and molecular mechanisms of DE remain elusive. In this study, we found that low water-deficit treatment (LWT) at the early stage of rice development can trigger early flowering and reduced tiller numbers. LWT induced the accumulation of abscisic acid (ABA), which in turn has feed-back effects on light perception and circadian clock by synchronously regulating many flowering-related genes to promote early flowering. Moreover, some of light receptors, circadian components, and flowering-related genes including OsTOC1, Ghd7, and PhyB were found to be involved in LWT in an ABA-dependent manner, whereas some of the other flowering-related genes including OsGI, OsELF3, OsPRR37, and OsMADS50 were involved in the regulation of DE independent of ABA. In addition, we found that strigolactones and OsTB1 are involved in the tillering inhibition under LWT, which is independent of the flowering pathway in rice. Taken together, our findings provide compelling evidence that DE in rice is coordinately regulated by multiple pathways during the reproduction (flowering) switch.

The effect of artificial selection on phenotypic plasticity in maize

Remarkable productivity has been achieved in crop species through artificial selection and adaptation to modern agronomic practices. Whether intensive selection has changed the ability of improved cultivars to maintain high productivity across variable environments is unknown. Understanding the genetic control of phenotypic plasticity and genotype by environment (G × E) interaction will enhance crop performance predictions across diverse environments. Here we use data generated from the Genomes to Fields (G2F) Maize G × E project to assess the effect of selection on G × E variation and characterize polymorphisms associated with plasticity. Genomic regions putatively selected during modern temperate maize breeding explain less variability for yield G × E than unselected regions, indicating that improvement by breeding may have reduced G × E of modern temperate cultivars. Trends in genomic position of variants associated with stability reveal fewer genic associations and enrichment of variants 0–5000 base pairs upstream of genes, hypothetically due to control of plasticity by short-range regulatory elements.

Breeding has increased crop productivity, but whether it has also changed phenotypic plasticity is unclear. Here, the authors find maize genomic regions selected for high productivity show reduced contribution to genotype by environment variation and provide evidence for regulatory control of phenotypic stability.

Transcriptome predictors of coral survival and growth in a highly variable environment

Concern over rapid environmental shifts associated with climate change has led to a search for molecular markers of environmental tolerance. Climate-associated gene expression profiles exist for a number of systems, but have rarely been tied to fitness outcomes, especially in nonmodel organisms. We reciprocally transplanted corals between two backreef locations with more and less variable temperature regimes to disentangle effects of recent and native environment on survival and growth. Coral growth over 12 months was largely determined by local environment. Survival, however, was impacted by native environment; corals from the more variable environment had 22% higher survivorship. By contrast, corals native to the less variable environment had more variable survival. This might represent a "selective sieve" where poor survivors are filtered from the more stressful environment. We also find a potential fitness trade-off-corals with high survival under stressful conditions grew less in the more benign environment. Transcriptome samples taken a year before transplantation were used to examine gene expression patterns that predicted transplant survival and growth. Two separate clusters of coexpressed genes were predictive of survival in the two locations. Genes from these clusters are candidate biomarkers for predicting persistence of corals under future climate change scenarios.

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site-specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole-genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree-based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.

Adaptation to an invasive host is driving the loss of a native ecotype

Locally adapted populations are often used as model systems for the early stages of ecological speciation, but most of these young divergent populations will never become complete species. The maintenance of local adaptation relies on the strength of natural selection overwhelming the homogenizing effects of gene flow; however, this balance may be readily upset in changing environments. Here I show that soapberry bugs (Jadera haematoloma) have lost adaptations to their native host plant (Cardiospermum corindum) and are regionally specializing on an invasive host (Koelreuteria elegans), collapsing a classic and well-documented example of local adaptation. All populations that were adapted to the native host-including those still found on that host today-are now better adapted to the invasive host in multiple phenotypes. Weak differentiation remains in two traits, suggesting that homogenization across the region is incomplete. This study highlights the potential for adaptation to invasive species to disrupt native communities by swamping adaptation to native conditions through maladaptive gene flow.