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.

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.

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.

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.

Population responses to a historic drought across the range of the common monkeyflower (Mimulus guttatus).

PREMISE: Due to climate change, more frequent and intense periodic droughts are predicted to increasingly pose major challenges to the persistence of plant populations. When a severe drought occurs over a broad geographical region, independent responses by individual populations provide replicated natural experiments for examining the evolution of drought resistance and the potential for evolutionary rescue. METHODS: We used a resurrection approach to examine trait evolution in populations of the common monkeyflower, Mimulus guttatus, exposed to a record drought in California from 2011 to 2017. Specifically, we compared variation in traits related to drought escape and avoidance from seeds collected from 37 populations pre- and post-drought in a common garden. In a parallel experiment, we evaluated fitness in two populations, one which thrived and one which was nearly extirpated during the drought, under well-watered and dry-down conditions. RESULTS: We observed substantial variation among populations in trait evolution. In the subset of populations where phenotypes changed significantly, divergence proceeded along trait correlations with some populations flowering rapidly with less vegetative tissue accumulation and others delaying flowering with greater vegetative tissue accumulation. The degree of trait evolution was only weakly correlated with drought intensity but strongly correlated with initial levels of standing variation. Fitness was higher in the post-drought than pre-drought accessions in both treatments for the thriving population, but lower in both treatments for the nearly extirpated population. CONCLUSIONS: Together, our results indicate that evolutionary responses to drought are context dependent and reflect the standing genetic variation and genetic correlations present within populations.

The Genetic Architecture of Plant Defense Trade-offs in a Common Monkeyflower.

Determining how adaptive combinations of traits arose requires understanding the prevalence and scope of genetic constraints. Frequently observed phenotypic correlations between plant growth, defenses, and/or reproductive timing have led researchers to suggest that pleiotropy or strong genetic linkage between variants affecting independent traits is pervasive. Alternatively, these correlations could arise via independent mutations in different genes for each trait and extensive correlational selection. Here we evaluate these alternatives by conducting a quantitative trait loci (QTL) mapping experiment involving a cross between 2 populations of common monkeyflower (Mimulus guttatus) that differ in growth rate as well as total concentration and arsenal composition of plant defense compounds, phenylpropanoid glycosides (PPGs). We find no evidence that pleiotropy underlies correlations between defense and growth rate. However, there is a strong genetic correlation between levels of total PPGs and flowering time that is largely attributable to a single shared QTL. While this result suggests a role for pleiotropy/close linkage, several other QTLs also contribute to variation in total PPGs. Additionally, divergent PPG arsenals are influenced by a number of smaller-effect QTLs that each underlie variation in 1 or 2 PPGs. This result indicates that chemical defense arsenals can be finely adapted to biotic environments despite sharing a common biochemical precursor. Together, our results show correlations between defense and life-history traits are influenced by pleiotropy or genetic linkage, but genetic constraints may have limited impact on future evolutionary responses, as a substantial proportion of variation in each trait is controlled by independent loci.

The role of cold cues at different life stages on germination and flowering phenology.

PREMISE OF THE STUDY: The timing of major phenological transitions is critical to lifetime fitness, and life history theory predicts differences for annual and perennial plants. To correctly time these transitions, many plants rely on environmental cues such as exposure to extended periods of cold, which may occur at different stages throughout their lifetime. METHODS: We studied the role of cold at different life stages, by jointly exposing seed (stratification) and rosettes (vernalization) to cold. We used 23 populations of Mimulus guttatus, which vary from annuals to perennials, and investigated how cold at one or both stages affected germination, flowering, growth, and biomass. KEY RESULTS: We found that stratification and vernalization interact to affect life cycle transitions, and that cold at either stage could synchronize flowering phenology. For perennials, either stratification or vernalization is necessary for maximum flowering. We also found that germination timing covaried with later traits. Moreover, plants from environments with dissimilar climates displayed different phenological responses to stratification or vernalization. CONCLUSIONS: In general, cold is more important for seed germination in annuals and plants from environments with warm temperatures and variable precipitation. In contrast, cold is more important for flowering in perennials: it accelerates flowering in plants from lower precipitation environments, and it increases flowering proportion in plants from cooler, more stable precipitation environments. We discuss our findings in the context of the variable environments plants experience within a population and the variation encountered across the biogeographic native range of the species.

The Climatic Variability Hypothesis and trade-offs in thermal performance in coastal and inland populations of Mimulus guttatus.

Ecologists and evolutionary biologists have long predicted that organisms in more climatically variable environments should be adapted to handle a wider range of conditions. This intuitive idea, known as the Climatic Variability Hypothesis (CVH), has gained mixed support from empirical studies. We tested the CVH in a novel system by comparing the thermal breadth of coastal and inland populations of Mimulus guttatus. To quantify thermal breadth, we performed a thermal performance experiment and built performance curves. Using these performance curves, we also evaluated evidence for a breadth-performance trade-off and the Hotter-is-Better hypothesis. We did not find support for the CVH; coastal and inland populations did not differ in thermal breadth. However, we found evidence for a breadth-performance trade-off and the Hotter-is-Better hypothesis. Surprisingly, the two most inland populations differed the most in the thermal performance traits we evaluated. Our results highlight the importance of explicitly measuring thermal performance to test explanations of species distribution patterns and the need to examine alternative mechanisms by which organisms occupy different climatic regimes.

Life history and chemical defense interact to drive patterns of adaptation in an annual monkeyflower.

Although chemical defenses and herbivore pressure are widely established as key targets and agents of selection, their roles in local adaptation and determining potential evolutionary responses to changing climates are often neglected. Here, we explore fitness differences between 11 rangewide M. guttatus populations in a field common garden experiment and assess the agents and targets of selection driving relative fitness patterns. We use piecewise structural equation models to disentangle associations between chemical defenses, (phenylpropanoid glycosides; PPGs), and life history traits with herbivory and fitness. While the historical environment of populations is not predictive of fitness differences between populations, >90% of variation in fitness can be predicted by the flowering time and foliar PPG defense arsenal of a population. Piecewise structural equation models indicate that life history traits, particularly earlier flowering time, are strongly and directly linked to fitness. However, herbivory, particularly fruit predation, is also an important agent of selection that creates indirect links between fitness and both chemical defenses and life history traits. Our results emphasize the multivariate nature of the agents and targets of selections in producing adaptation and suggest that future responses to selection must navigate a complex fitness landscape.

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.

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 (FIS ) 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.

Abiotic Environment Predicts Micro- but Not Macroevolutionary Patterns of Flower Color in Monkeyflowers (Phrymaceae).

Anthocyanin pigments are responsible for many of the vivid pink, purple, red, and blue flower colors across angiosperms and frequently vary within and between closely related species. While anthocyanins are well known to influence pollinator attraction, they are also associated with tolerance to abiotic stressors such as extreme temperatures, reduced precipitation, and ultraviolet radiation. Using a comparative approach, we tested whether abiotic variables predict floral anthocyanin in monkeyflowers (Phrymaceae) across western North America. Within two polymorphic species, we found that abiotic variables predicted flower color across their geographic ranges. In Erythranthe discolor, the frequency of pink flowered (anthocyanin producing) individuals was greater in populations with reduced precipitation. In Diplacus mephiticus, the frequency of pink flowered individuals was greater at higher elevations that had reduced precipitation and lower temperatures but less ultraviolet radiation. At the macroevolutionary scale, across two parallel radiations of North American monkeyflowers, species with floral anthocyanins (pink, purple, or red corollas) occupied areas with reduced precipitation in Erythranthe but not Diplacus. However, after accounting for phylogenetic relatedness, we found no evidence for the joint evolution of flower color and environmental affinity in either clade. We conclude that although abiotic stressors may play a role in the evolution of flower color within polymorphic species, we found no evidence that these processes lead to macroevolutionary patterns across monkeyflowers.

Regional differences in rapid evolution during severe drought.

Climate change is increasing drought intensity, threatening biodiversity. Rapid evolution of drought adaptations might be required for population persistence, particularly in rear-edge populations that may already be closer to physiological limits. Resurrection studies are a useful tool to assess adaptation to climate change, yet these studies rarely encompass the geographic range of a species. Here, we sampled 11 populations of scarlet monkeyflower (Mimulus cardinalis), collecting seeds across the plants' northern, central, and southern range to track trait evolution from the lowest to the greatest moisture anomaly over a 7-year period. We grew families generated from these populations across well-watered and terminal drought treatments in a greenhouse and quantified five traits associated with dehydration escape and avoidance. When considering pre-drought to peak-drought phenotypes, we find that later date of flowering evolved across the range of M. cardinalis, suggesting a shift away from dehydration escape. Instead, traits consistent with dehydration avoidance evolved, with smaller and/or thicker leaves evolving in central and southern regions. The southern region also saw a loss of plasticity in these leaf traits by the peak of the drought, whereas flowering time remained plastic across all regions. This observed shift in traits from escape to avoidance occurred only in certain regions, revealing the importance of geographic context when examining adaptations to climate change.

Effect of Whole-Genome Duplication on the Evolutionary Rescue of Sterile Hybrid Monkeyflowers.

Hybridization is a creative evolutionary force, increasing genomic diversity and facilitating adaptation and even speciation. Hybrids often face significant challenges to establishment, including reduced fertility that arises from genomic incompatibilities between their parents. Whole-genome duplication in hybrids (allopolyploidy) can restore fertility, cause immediate phenotypic changes, and generate reproductive isolation. Yet the survival of polyploid lineages is uncertain, and few studies have compared the performance of recently formed allopolyploids and their parents under field conditions. Here, we use natural and synthetically produced hybrid and polyploid monkeyflowers (Mimulus spp.) to study how polyploidy contributes to the fertility, reproductive isolation, phenotype, and performance of hybrids in the field. We find that polyploidization restores fertility and that allopolyploids are reproductively isolated from their parents. The phenotype of allopolyploids displays the classic gigas effect of whole-genome duplication, in which plants have larger organs and are slower to flower. Field experiments indicate that survival of synthetic hybrids before and after polyploidization is intermediate between that of the parents, whereas natural hybrids have higher survival than all other taxa. We conclude that hybridization and polyploidy can act as sources of genomic novelty, but adaptive evolution is key in mediating the establishment of young allopolyploid lineages.

A resurrection study reveals limited evolution of phenology in response to recent climate change across the geographic range of the scarlet monkeyflower.

PREMISE OF THE STUDY: As global climate change alters drought regimes, rapid evolution of traits that facilitate adaptation to drought can rescue populations in decline. The evolution of phenological advancement can allow plant populations to escape drought, but evolutionary responses in phenology can vary across a species' range due to differences in drought intensity and standing genetic variation. METHODS: Mimulus cardinalis, a perennial herb spanning a broad climatic gradient, recently experienced a period of record drought. Here, we used a resurrection study comparing flowering time and stem height at first flower of pre-drought ancestors and post-drought descendants from northern-edge, central, and southern-edge populations in a common environment to examine the evolution of drought escape across the latitudinal range. KEY RESULTS: Contrary to the hypothesis of the evolution of advanced phenology in response to recent drought, flowering time did not advance between ancestors and descendants in any population, though storage condition and maternal effects could have impacted these results. Stem height was positively correlated with flowering time, such that plants that flowered earlier were shorter at first flower. This correlation could constrain the evolution of earlier flowering time if selection favors flowering early at a large size. CONCLUSIONS: These findings suggest that rapid evolution of phenology will not rescue these populations from recent climate change. Future work is needed to examine the potential for the evolution of alternative drought strategies and phenotypic plasticity to buffer M. cardinalis populations from changing climate.

Populations Are Differentiated in Biological Rhythms without Explicit Elevational Clines in the Plant Mimulus laciniatus.

Environmental variation along an elevational gradient can yield phenotypic differentiation resulting from varying selection pressures on plant traits related to seasonal responses. Thus, genetic clines can evolve in a suite of traits, including the circadian clock, that drives daily cycling in varied traits and that shares its genetic background with adaptation to seasonality. We used populations of annual Mimulus laciniatus from different elevations in the Sierra Nevada in California to explore among-population differentiation in the circadian clock, flowering responses to photoperiod, and phenological traits (days to cotyledon emergence, days to flowering, and days to seed ripening) in controlled common-garden conditions. Further, we examined correlations of these traits with environmental variables related to temperature and precipitation. We observed that the circadian period in leaf movement was differentiated among populations sampled within about 100 km, with population means varying by 1.6 h. Significant local genetic variation occurred within 2 populations in which circadian period among families varied by up to 1.8 h. Replicated treatments with variable ecologically relevant photoperiods revealed marked population differentiation in critical day length for flowering that ranged from 11.0 to 14.1 h, corresponding to the time period between late February and mid-May in the wild. Flowering time varied among populations in a 14-h photoperiod. Regardless of this substantial population-level diversity, obvious linear clinality in trait variability across elevations could not be determined based on our genotypic sample; it is possible that more complex spatial patterns of variation arise in complex terrains such as those in the Sierra Nevada. Moreover, we did not find statistically significant bivariate correlations between population means of different traits. Our research contributes to the understanding of genetic variation in the circadian clock and in seasonal responses in natural populations, highlighting the need for more comprehensive investigations on the association between the clock and other adaptive traits in plants.

Two MYB Proteins in a Self-Organizing Activator-Inhibitor System Produce Spotted Pigmentation Patterns.

Many organisms exhibit visually striking spotted or striped pigmentation patterns. Developmental models predict that such spatial patterns can form when a local autocatalytic feedback loop and a long-range inhibitory feedback loop interact. At its simplest, this self-organizing network only requires one self-activating activator that also activates a repressor, which inhibits the activator and diffuses to neighboring cells. However, the molecular activators and inhibitors fully fitting this versatile model remain elusive in pigmentation systems. Here, we characterize an R2R3-MYB activator and an R3-MYB repressor in monkeyflowers (Mimulus). Through experimental perturbation and mathematical modeling, we demonstrate that the properties of these two proteins correspond to an activator-inhibitor pair in a two-component, reaction-diffusion system, explaining the formation of dispersed anthocyanin spots in monkeyflower petals. Notably, disrupting this pattern impacts pollinator visitation. Thus, subtle changes in simple activator-inhibitor systems are likely essential contributors to the evolution of the remarkable diversity of pigmentation patterns in flowers.

Assortative mating by flowering time and its effect on correlated traits in variable environments.

Reproductive timing is a key life-history trait that impacts the pool of available mates, the environment experienced during flowering, and the expression of other traits through genetic covariation. Selection on phenology, and its consequences on other life-history traits, has considerable implications in the context of ongoing climate change and shifting growing seasons. To test this, we grew field-collected seed from the wildflower Mimulus guttatus in a greenhouse to assess the standing genetic variation for flowering time and covariation with other traits. We then created full-sib families through phenological assortative mating and grew offspring in three photoperiod treatments representing seasonal variation in daylength. We find substantial quantitative genetic variation for the onset of flowering time, which covaried with vegetative traits. The assortatively-mated offspring varied in their critical photoperiod by over two hours, so that families differed in their probability of flowering across treatments Allocation to flowering and vegetative growth changed across the daylength treatments, with consistent direction and magnitude of covariation among flowering time and other traits. Our results suggest that future studies of flowering time evolution should consider the joint evolution of correlated traits and shifting seasonal selection to understand how environmental variation influences life histories.

Mimulus peregrinus (Phrymaceae): A new British allopolyploid species.

Polyploidization plays an important role in species formation as chromosome doubling results in strong reproductive isolation between derivative and parental taxa. In this note I describe a new species, Mimulus peregrinus (Phrymaceae), which represents the first recorded instance of a new British polyploid species of Mimulus (2n = 6x = 92) that has arisen since the introduction of this genus into the United Kingdom in the 1800's. Mimulus peregrinus presents floral and vegetative characteristics intermediate between Mimulus guttatus and Mimulus luteus, but can be distinguished from all naturalized British Mimulus species and hybrids based on a combination of reproductive and vegetative traits. Mimulus peregrinus displays high pollen and seed fertility as well as traits usually associated with genome doubling such as increased pollen and stomata size. The intermediate characteristics of Mimulus peregrinus between Mimulus guttatus (2n = 2x = 28)and Mimulus luteus (2n = 4x = 60-62), and its close affinity with the highly sterile, triploid (2n = 3x = 44-45) hybrid taxon Mimulus × robertsii (Mimulus guttatus × Mimulus luteus), suggests that Mimulus peregrinus mayconstitute an example of recent allopolyploid speciation.