Mimulus is an emerging model system for the integration of ecological and genomic studies

Within-species floral evolution reveals convergence in adaptive walks during incipient pollinator shift

Understanding how evolution proceeds from molecules to organisms to interactions requires integrative studies spanning biological levels. Linking phenotypes with associated genes and fitness illuminates how adaptive walks move organisms between fitness peaks. Floral evolution can confer rapid reproductive isolation, often converging in association with pollinator guilds. Within the monkeyflowers (Mimulus sect. Erythranthe), yellow flowers within red hummingbird-pollinated species have arisen at least twice, suggesting possible pollinator shifts. We compare two yellow-flowered forms of M. cardinalis and M. verbenaceus to their red counterparts in floral phenotypes, biochemistry, transcriptomic and genomic variation, and pollinator interactions. We find convergence in ongoing adaptive walks of both yellow morphs, with consistent changes in traits of large effect (floral pigments, associated gene expression), resulting in strong preference for yellow flowers by bumblebees. Shifts in scent emission and floral opening size also favor bee adaptation, suggesting smaller-effect steps from hummingbird to bee pollination. By examining intraspecific, incipient pollinator shifts in two related species, we elucidate adaptive walks at early stages, revealing how convergent large effect mutations (floral color) may drive pollinator attraction, followed by smaller effect changes for mechanical fit and reward access. Thus, ongoing adaptive walks may impact reproductive isolation and incipient speciation via convergent evolution.

Progeny array analysis to estimate outcrossing rates, inbreeding coefficients, and inbreeding depression among native, naturalized, and invasive populations of Mimulus guttatus (Phrymaceae)

The mating system of non-native plant populations plays a role in determining the colonizing success following introduction into locations outside of the native distribution. For plant species capable of mixed-mating, both selfing and outcrossing can be advantageous and promote the establishment, persistence, and spread of newly arrived populations. To investigate how mating systems may contribute to the invasion process we estimated mating system parameters in perennial populations of the model plant species, Mimulus guttatus from its native range (West coast USA), non-native populations that are established but have not become invasive (East coast USA, >50 years), and populations in invasive regions (UK >200 years). Studies that include mating system data across the continuum of the invasion process are rare, thus here we utilize molecular markers to estimate outcrossing rates (t), inbreeding coefficients (F), and inbreeding depression in native, naturalized, and invasive populations. Overall, we found support for the persistence of mixed-mating across populations, variability in the relationship between outcrossing rates and inbreeding depression across populations, and evidence for the bridgehead process, where non-native populations may be the sources for the further establishment or reinforcement of additional non-native populations. The methodology we deployed had its own assumptions and sampling design constraints, that contributed to the variation in the parameter estimates. All M. guttatus populations likely rely on selfing through both within clone, and within flower and plant mating in addition to vegetative propagation. The results underscore the importance of introduction history in determining the role of both sexual and asexual reproduction in the successful establishment of a plant species outside its native range.

Unveiling the Behavior of an Endangered Facultative Cuprophyte Coincya Species in an Abandoned Copper Mine (Southeast Portugal)

Plant-soil interactions of endangered species with a high-priority conservation status are important to define in situ and ex situ conservation and restoration projects. The threatened endemic Coincya transtagana, thriving in the southwest of the Iberian Peninsula, can grow in metalliferous soils. The main goal of this study was to investigate the behavior of this species in soils rich in potentially toxic elements in the abandoned Aparis Cu mine. Soil samples were characterized for physicochemical properties and multielemental composition, as well as biological activity, through an analysis of enzymatic activities. Plant biomass was assessed, and multielemental analysis of the plants was also performed. The mine soils had slightly basic pH values and were non-saline and poor in mineral N-NH4, with medium-to-high organic matter concentration and medium cation-exchange capacity. In these soils, dehydrogenase had the highest activity, whereas protease had the lowest activity. The total concentrations of Cu (1.3-5.9 g/kg) and As (37.9-118 mg/kg) in soils were very high, and the available fraction of Cu in the soil also had high concentration values (49-491 mg/kg). Moreover, this study shows for the first time that C. transtagana had high uptake and translocation capacities from roots to shoots for Cu, Ni, and Cr. Although Cu in the plants' aerial parts (40-286 mg/kg) was considered excessive/toxic, no signs of plant toxicity disorders or P uptake reduction were detected. This preliminary study revealed that C. transtagana is Cu-tolerant, and it could be used for phytoremediation of soils contaminated with potentially toxic elements, while also contributing to its conservation.

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

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

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.

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

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

Suillus: an emerging model for the study of ectomycorrhizal ecology and evolution

Research on mycorrhizal symbiosis has been slowed by a lack of established study systems. To address this challenge, we have been developing Suillus, a widespread ecologically and economically relevant fungal genus primarily associated with the plant family Pinaceae, into a model system for studying ectomycorrhizal (ECM) associations. Over the last decade, we have compiled extensive genomic resources, culture libraries, a phenotype database, and protocols for manipulating Suillus fungi with and without their tree partners. Our efforts have already resulted in a large number of publicly available genomes, transcriptomes, and respective annotations, as well as advances in our understanding of mycorrhizal partner specificity and host communication, fungal and plant nutrition, environmental adaptation, soil nutrient cycling, interspecific competition, and biological invasions. Here, we highlight the most significant recent findings enabled by Suillus, present a suite of protocols for working with the genus, and discuss how Suillus is emerging as an important model to elucidate the ecology and evolution of ECM interactions.

Hybrid Mimulus flowers attract a new pollinator

Hybridization is common in flowering plants and is believed to be an important force driving adaptation and speciation. The flowers of hybrids often exhibit new trait combinations, which, theoretically, could attract new species of pollinators. In this study, we found that the hybrids between a hummingbird-pollinated species Mimulus cardinalis and a self-pollinated species Mimulus parishii attract bumblebees (Bombus impatiens), a pollinator not attracted to either of the progenitor species. This novel attraction is explained by new combinations of floral traits in hybrids, including, most importantly, petal color, in addition to nectar concentration and corolla size. To understand how petal color variation is perceived by bumblebees, we performed reflectance spectroscopy and multispectral imaging to model the flower appearance in bee vision. This analysis showed that color variation would impact the ease of detection. We also found that YUP, the genetic locus responsible for a large portion of floral color variation and previously shown to be important in bee interactions with other Mimulus species, also played an important role in this novel attraction. These results together suggest that the attraction of new pollinators to hybrid plants could be an underexplored avenue for pollinator shift and speciation.

Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

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

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

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

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

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

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

PREMISE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Genetic diversity and population structure among native, naturalized, and invasive populations of the common yellow monkeyflower, Mimulus guttatus (Phrymaceae)

An ongoing controversy in invasion biology is the prevalence of colonizing plant populations that are able to establish and spread, while maintaining limited amounts of genetic variation. Invasive populations can be established through several routes including from a single source or from multiple introductions. The aim of this study was to examine genetic diversity in populations of Mimulus guttatus in the United Kingdom, where the species is considered invasive, and compare this diversity to that in native populations on the west coast of North America. Additionally, we looked at diversity in non-native populations that have not yet become invasive (naturalized populations) in eastern North America. We investigated population structure among populations in these three regions and attempted to uncover the sources for populations that have established in the naturalized and invasive regions. We found that genetic diversity was, on average, relatively high in populations from the invasive UK region and comparable to native populations. Contrastingly, two naturalized M. guttatus populations were low in both genetic and genotypic diversity, indicating a history of asexual reproduction and self-fertilization. A third naturalized population was found to be a polyploid Mimulus hybrid of unknown origin. Our results demonstrate that M. guttatus has likely achieved colonization success outside of its native western North America distribution by a variety of establishment pathways, including those with genetic and demographic benefits resulting from multiple introductions in the UK, reproductive assurance through selfing, and asexual reproduction in eastern North America, and possible polyploidization in one Canadian population.

Spatially and temporally varying selection influence species boundaries in two sympatric Mimulus

Spatially and temporally varying selection can maintain genetic variation within and between populations, but it is less well known how these forces influence divergence between closely related species. We identify the interaction of temporal and spatial variation in selection and their role in either reinforcing or eroding divergence between two closely related Mimulus species. Using repeated reciprocal transplant experiments with advanced generation hybrids, we compare the strength of selection on quantitative traits involved in adaptation and reproductive isolation in Mimulus guttatus and Mimulus laciniatus between two years with dramatically different water availability. We found strong divergent habitat-mediated selection on traits in the direction of species differences during a drought in 2013, suggesting that spatially varying selection maintains species divergence. However, a relaxation in divergent selection on most traits in an unusually wet year (2019), including flowering time, which is involved in pre-zygotic isolation, suggests that temporal variation in selection may weaken species differences. Therefore, we find evidence that temporally and spatially varying selection may have opposing roles in mediating species boundaries. Given our changing climate, future growing seasons are expected to be more similar to the dry year, suggesting that in this system climate change may actually increase species divergence.

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

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

Populations of western North American monkeyflowers accrue niche breadth primarily via genotypic divergence in environmental optima

Niche breadth, the range of environments that individuals, populations, and species can tolerate, is a fundamental ecological and evolutionary property, yet few studies have examined how niche breadth is partitioned across biological scales. We use a published dataset of thermal performance for a single population from each of 10 closely related species of western North American monkeyflowers (genus Mimulus) to investigate whether populations achieve broad thermal niches through general purpose genotypes, specialized genotypes with divergent environmental optima, and/or variation among genotypes in the degree of generalization. We found the strongest relative support for the hypothesis that populations with greater genetic variation for thermal optimum had broader thermal niches, and for every unit increase in among-family variance in thermal optimum, population-level thermal breadth increased by 0.508°C. While the niche breadth of a single genotype represented up to 86% of population-level niche breadth, genotype-level niche breadth had a weaker positive effect on population-level breadth, with every 1°C increase in genotypic thermal breadth resulting in a 0.062°C increase in population breadth. Genetic variation for thermal breadth was not predictive of population-level thermal breadth. These findings suggest that populations of Mimulus species have achieved broad thermal niches primarily through genotypes with divergent thermal optima and to a lesser extent via general-purpose genotypes. Future work examining additional biological hierarchies would provide a more comprehensive understanding of how niche breadth partitioning impacts the vulnerabilities of individuals, populations, and species to environmental change.

Fitness effects of somatic mutations accumulating during vegetative growth

The unique life form of plants promotes the accumulation of somatic mutations that can be passed to offspring in the next generation, because the same meristem cells responsible for vegetative growth also generate gametes for sexual reproduction. However, little is known about the consequences of somatic mutation accumulation for offspring fitness. We evaluate the fitness effects of somatic mutations in Mimulus guttatus by comparing progeny from self-pollinations made within the same flower (autogamy) to progeny from self-pollinations made between stems on the same plant (geitonogamy). The effects of somatic mutations are evident from this comparison, as autogamy leads to homozygosity of a proportion of somatic mutations, but progeny from geitonogamy remain heterozygous for mutations unique to each stem. In two different experiments, we find consistent fitness effects of somatic mutations from individual stems. Surprisingly, several progeny groups from autogamous crosses displayed increases in fitness compared to progeny from geitonogamy crosses, likely indicating that beneficial somatic mutations occurred in some stems. These results support the hypothesis that somatic mutations accumulate during vegetative growth, but they are filtered by different forms of selection that occur throughout development, resulting in the culling of expressed deleterious mutations and the retention of beneficial 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.

Chromosome-Scale Genome Assembly of Gilia yorkii Enables Genetic Mapping of Floral Traits in an Interspecies Cross

Substantial morphological variation in land plants remains inaccessible to genetic analysis because current models lack variation in important ecological and agronomic traits. The genus Gilia was historically a model for biosystematics studies and includes variation in morphological traits that are poorly understood at the genetic level. We assembled a chromosome-scale reference genome of G. yorkii and used it to investigate genome evolution in the Polemoniaceae. We performed QTL (quantitative trait loci) mapping in a G. yorkii×G. capitata interspecific population for traits related to inflorescence architecture and flower color. The genome assembly spans 2.75 Gb of the estimated 2.80-Gb genome, with 96.7% of the sequence contained in the nine largest chromosome-scale scaffolds matching the haploid chromosome number. Gilia yorkii experienced at least one round of whole-genome duplication shared with other Polemoniaceae after the eudicot paleohexaploidization event. We identified QTL linked to variation in inflorescence architecture and petal color, including a candidate for the major flower color QTL-a tandem duplication of flavanol 3',5'-hydroxylase. Our results demonstrate the utility of Gilia as a forward genetic model for dissecting the evolution of development in plants including the causal loci underlying inflorescence architecture transitions.

Variation in seasonal timing traits and life history along a latitudinal transect in Mimulus ringens

Seasonal timing traits are commonly under recurrent, spatially variable selection, and are therefore predicted to exhibit clinal variation. Temperate perennial plants often require vernalization to prompt growth and reproduction; however, little is known about whether vernalization requirements change across the range of a broadly distributed species. We performed a critical vernalization duration study in Mimulus ringens, coupled with population genomic analysis. Plants from eight populations spanning the latitudinal range were exposed to varying durations of 4°C vernalization between 0 and 56 days, and flowering response was assessed. RADSeq was also performed to generate 1179 polymorphic SNPs, which were used to examine population structure. We found unexpected life history variation, with some populations lacking vernalization requirement. Population genomic analyses show that these life history variants are highly divergent from perennials, potentially revealing a cryptic species. For perennial populations, minimum vernalization time was surprisingly consistent. However, once vernalized, northern populations flowered almost 3 weeks faster than southern. Furthermore, southern populations exhibited sensitivity to vernalization times beyond flowering competency, suggesting an ability to respond adaptively to different lengths of winter. Mimulus ringens, therefore, reveals evidence of clinal variation, and provides opportunities for future studies addressing mechanistic and ecological hypotheses both within and between incipient species.

Rescaling Biology: Increasing Integration Across Biological Scales and Subdisciplines to Enhance Understanding and Prediction

The organization of the living world covers a vast range of spatiotemporal scales, from molecules to the biosphere, seconds to centuries. Biologists working within specialized subdisciplines tend to focus on different ranges of scales. Therefore, developing frameworks that enable testing questions and predictions of scaling require sufficient understanding of complex processes across biological subdisciplines and spatiotemporal scales. Frameworks that enable scaling across subdisciplines would ideally allow us to test hypotheses about the degree to which explicit integration across spatiotemporal scales is needed for predicting the outcome of biological processes. For instance, how does genomic variation within populations allow us to explain community structure? How do the dynamics of cellular metabolism translate to our understanding of whole-ecosystem metabolism? Do patterns and processes operate seamlessly across biological scales, or are there fundamental laws of biological scaling that limit our ability to make predictions from one scale to another? Similarly, can sub-organismal structures and processes be sufficiently understood in isolation of potential feedbacks from the population, community, or ecosystem levels? And can we infer the sub-organismal processes from data on the population, community, or ecosystem scale? Concerted efforts to develop more cross-disciplinary frameworks will open doors to a more fully integrated field of biology. In this paper we discuss how we might integrate across scales, specifically by 1. Identifying scales and boundaries, 2. Determining analogous units and processes across scales, 3. Developing frameworks to unite multiple scales, and 4. Extending frameworks to new empirical systems.

New genomic resources and comparative analyses reveal differences in floral gene expression in selfing and outcrossing Collinsia sister species

The evolutionary transition from outcross- to self-fertilization is one of the most common in angiosperms and is often associated with a parallel shift in floral morphological and developmental traits, such as reduced flower size and pollen to ovule ratios, known as the "selfing syndrome." How these convergent phenotypes arise, the extent to which they are shaped by selection, and the nature of their underlying genetic basis are unsettled questions in evolutionary biology. The genus Collinsia (Plantaginaceae) includes seven independent transitions from outcrossing or mixed mating to high selfing rates accompanied by selfing syndrome traits. Accordingly, Collinsia represents an ideal system for investigating this parallelism, but requires genomic resource development. We present a high quality de novo genome assembly for the highly selfing species Collinsia rattanii. To begin addressing the basis of selfing syndrome developmental shifts, we evaluate and contrast patterns of gene expression from floral transcriptomes across three stages of bud development for C. rattanii and its outcrossing sister species Collinsia linearis. Relative to C. linearis, total gene expression is less variable among individuals and bud stages in C. rattanii. In addition, there is a common pattern among differentially expressed genes: lower expression levels that are more constant across bud development in C. rattanii relative to C. linearis. Transcriptional regulation of enzymes involved in pollen formation specifically in early bud development may influence floral traits that distinguish selfing and outcrossing Collinsia species through pleiotropic functions. Future work will include additional Collinsia outcrossing-selfing species pairs to identify genomic signatures of parallel evolution.

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

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

Model Systems in Ecology, Evolution, and Behavior: A Call for Diversity in Our Model Systems and Discipline

AbstractEcologists and evolutionary biologists are fascinated by life's variation but also seek to understand phenomena and mechanisms that apply broadly across taxa. Model systems can help us extract generalities from amid all the wondrous diversity, but only if we choose and develop them carefully, use them wisely, and have a range of model systems from which to choose. In this introduction to the Special Feature on Model Systems in Ecology, Evolution, and Behavior (EEB), we begin by grappling with the question, What is a model system? We then explore where our model systems come from, in terms of the skills and other attributes required to develop them and the historical biases that influence traditional model systems in EEB. We emphasize the importance of communities of scientists in the success of model systems-narrow scientific communities can restrict the model organisms themselves. We also consider how our discipline was built around one type of "model scientist"-a history still reflected in the field. This lack of diversity in EEB is unjust and also narrows the field's perspective, including by restricting the questions asked and talents used to answer them. Increasing diversity, equity, and inclusion will require acting at many levels, including structural changes. Diversity in EEB, in both model systems and the scientists who use them, strengthens our discipline.

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

Rapid local adaptation in both sexual and asexual invasive populations of monkeyflowers (Mimulus spp.)

BACKGROUND AND AIMS

Traditionally, local adaptation has been seen as the outcome of a long evolutionary history, particularly with regard to sexual lineages. By contrast, phenotypic plasticity has been thought to be most important during the initial stages of population establishment and in asexual species. We evaluated the roles of adaptive evolution and phenotypic plasticity in the invasive success of two closely related species of invasive monkeyflowers (Mimulus) in the UK that have contrasting reproductive strategies: M. guttatus combines sexual (seeds) and asexual (clonal growth) reproduction while M. × robertsii is entirely asexual.

METHODS

We compared the clonality (number of stolons), floral and vegetative phenotype, and phenotypic plasticity of native (M. guttatus) and invasive (M. guttatus and M. × robertsii) populations grown in controlled environment chambers under the environmental conditions at each latitudinal extreme of the UK. The goal was to discern the roles of temperature and photoperiod on the expression of phenotypic traits. Next, we tested the existence of local adaptation in the two species within the invasive range with a reciprocal transplant experiment at two field sites in the latitudinal extremes of the UK, and analysed which phenotypic traits underlie potential local fitness advantages in each species.

KEY RESULTS

Populations of M. guttatus in the UK showed local adaptation through sexual function (fruit production), while M. × robertsii showed local adaptation via asexual function (stolon production). Phenotypic selection analyses revealed that different traits are associated with fitness in each species. Invasive and native populations of M. guttatus had similar phenotypic plasticity and clonality. M. × robertsii presents greater plasticity and clonality than native M. guttatus, but most populations have restricted clonality under the warm conditions of the south of the UK.

CONCLUSIONS

This study provides experimental evidence of local adaptation in a strictly asexual invasive species with high clonality and phenotypic plasticity. This indicates that even asexual taxa can rapidly (<200 years) adapt to novel environmental conditions in which alternative strategies may not ensure the persistence of populations.

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.

Population genomic and historical analysis suggests a global invasion by bridgehead processes in Mimulus guttatus

Imperfect historical records and complex demographic histories present challenges for reconstructing the history of biological invasions. Here, we combine historical records, extensive worldwide and genome-wide sampling, and demographic analyses to investigate the global invasion of Mimulus guttatus from North America to Europe and the Southwest Pacific. By sampling 521 plants from 158 native and introduced populations genotyped at >44,000 loci, we determined that invasive M. guttatus was first likely introduced to the British Isles from the Aleutian Islands (Alaska), followed by admixture from multiple parts of the native range. We hypothesise that populations in the British Isles then served as a bridgehead for vanguard invasions worldwide. Our results emphasise the highly admixed nature of introduced M. guttatus and demonstrate the potential of introduced populations to serve as sources of secondary admixture, producing novel hybrids. Unravelling the history of biological invasions provides a starting point to understand how invasive populations adapt to novel environments.

Vallejo-Marín et al. combine historical records, extensive worldwide and genome-wide sampling, and demographic analyses to investigate the global invasion of Mimulus guttatus from North America to Europe and the Southwest Pacific. They found that M. guttatus was first likely introduced to the British Isles from the Aleutian Islands (Alaska), followed by admixture from multiple parts of the native range, and hypothesise that populations in the British Isles then served as a bridgehead for vanguard invasions worldwide.

Evolutionary history and ecology shape the diversity and abundance of phytochemical arsenals across monkeyflowers

We examine the extent to which phylogenetic effects and ecology are associated with macroevolutionary patterns of phytochemical defence production across the Mimulus phylogeny. We grew plants from 21 species representing the five major sections of the Mimulus phylogeny in a common garden to assess how the arsenals (NMDS groupings) and abundances (concentrations) of a phytochemical defence, phenylpropanoid glycosides (PPGs), vary across the phylogeny. Very few PPGs are widespread across the genus, but many are common to multiple sections of the genus. Phytochemical arsenals cluster among sections in an NMDS and are not associated with total concentration of PPGs. There is a strong phylogenetic signal for phytochemical arsenal composition across the Mimulus genus, whereas ecological variables such as growing season length, latitude, and elevation do not significantly influence arsenal. In contrast, there is little phylogenetic signal for total PPG concentration, and this trait is significantly influenced by several ecological factors. Phytochemical arsenals and abundances are influenced by plant life history form. Both phylogenetic effects and ecology are related to phytochemical patterns across species, albeit in different ways. The independence of phytochemical defence concentrations from arsenal compositions indicates that these aspects of defence may continue to evolve independently of one another.

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.

Predicting evolutionary change at the DNA level in a natural Mimulus population

Evolution by natural selection occurs when the frequencies of genetic variants change because individuals differ in Darwinian fitness components such as survival or reproductive success. Differential fitness has been demonstrated in field studies of many organisms, but it remains unclear how well we can quantitatively predict allele frequency changes from fitness measurements. Here, we characterize natural selection on millions of Single Nucleotide Polymorphisms (SNPs) across the genome of the annual plant Mimulus guttatus. We use fitness estimates to calibrate population genetic models that effectively predict allele frequency changes into the next generation. Hundreds of SNPs experienced "male selection" in 2013 with one allele at each SNP elevated in frequency among successful male gametes relative to the entire population of adults. In the following generation, allele frequencies at these SNPs consistently shifted in the predicted direction. A second year of study revealed that SNPs had effects on both viability and reproductive success with pervasive trade-offs between fitness components. SNPs favored by male selection were, on average, detrimental to survival. These trade-offs (antagonistic pleiotropy and temporal fluctuations in fitness) may be essential to the long-term maintenance of alleles. Despite the challenges of measuring selection in the wild, the strong correlation between predicted and observed allele frequency changes suggests that population genetic models have a much greater role to play in forward-time prediction of evolutionary change.

Multi-level patterns of genetic structure and isolation by distance in the widespread plant Mimulus guttatus

An understanding of genetic structure is essential for answering many questions in population genetics. However, complex population dynamics and scale-dependent processes can make it difficult to detect if there are distinct genetic clusters present in natural populations. Inferring discrete population structure is particularly challenging in the presence of continuous genetic variation such as isolation by distance. Here, we use the plant species Mimulus guttatus as a case study for understanding genetic structure at three spatial scales. We use reduced-representation sequencing and marker-based genotyping to understand dispersal dynamics and to characterise genetic structure. Our results provide insight into the spatial scale of genetic structure in a widespread plant species, and demonstrate how dispersal affects spatial genetic variation at the local, regional, and range-wide scale. At a fine-spatial scale, we show dispersal is rampant with little evidence of spatial genetic structure within populations. At a regional-scale, we show continuous differentiation driven by isolation by distance over hundreds of kilometres, with broad geographic genetic clusters that span major barriers to dispersal. Across Western North America, we observe geographic genetic structure and the genetic signature of multiple postglacial recolonisation events, with historical gene flow linking isolated populations. Our genetic analyses show M. guttatus is highly dispersive and maintains large metapopulations with high intrapopulation variation. This high diversity and dispersal confounds the inference of genetic structure, with multi-level sampling and spatially-explicit analyses required to understand population history.

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.

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.

Senecio as a model system for integrating studies of genotype, phenotype and fitness

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

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.

Patterns of Hybrid Seed Inviability in the Mimulus guttatus sp. Complex Reveal a Potential Role of Parental Conflict in Reproductive Isolation

Genomic conflicts may play a central role in the evolution of reproductive barriers. Theory predicts that early-onset hybrid inviability may stem from conflict between parents for resource allocation to offspring. Here, we describe M. decorus: a group of cryptic species within the M. guttatus species complex that are largely reproductively isolated by hybrid seed inviability (HSI). HSI between M. guttatus and M. decorus is common and strong, but populations of M. decorus vary in the magnitude and directionality of HSI with M. guttatus. Patterns of HSI between M. guttatus and M. decorus, as well as within M. decorus, conform to the predictions of parental conflict: first, reciprocal F1s exhibit size differences and parent-of-origin-specific endosperm defects; second, the extent of asymmetry between reciprocal F1 seed size is correlated with asymmetry in HSI; and third, inferred differences in the extent of conflict predict the extent of HSI between populations. We also find that HSI is rapidly evolving, as populations that exhibit the most HSI are each others' closest relative. Lastly, although all populations appear largely outcrossing, we find that the differences in the inferred strength of conflict scale positively with π, suggesting that demographic or life history factors other than transitions to self-fertilization may influence the rate of parental-conflict-driven evolution. Overall, these patterns suggest the rapid evolution of parent-of-origin-specific resource allocation alleles coincident with HSI within and between M. guttatus and M. decorus. Parental conflict may therefore be an important evolutionary driver of reproductive isolation.

Divergence in drought-response traits between sympatric species of Mimulus

Differential adaptation to local environmental conditions is thought to be an important driver of speciation. Plants, whose sedentary lifestyle necessitates fine-tuned adaptation to edaphic conditions such as water availability, are often distributed based on these conditions. Populations occupying water-limited habitats may employ a variety of strategies, involving numerous phenotypes, to prevent and withstand desiccation. In sympatry, two closely related Mimulus species-M. guttatus and M. nasutus-occupy distinct microhabitats that differ in seasonal water availability. In a common garden experiment, we characterized natural variation within and between sympatric M. guttatus and M. nasutus in the ability to successfully set seed under well-watered and drought conditions. We also measured key phenotypes for drought adaptation, including developmental timing, plant size, flower size, and stomatal density. Consistent with their microhabitat associations in nature, M. nasutus set seed much more successfully than M. guttatus under water-limited conditions. This divergence in reproductive output under drought was due to differences in mortality after the onset of flowering, with M. nasutus surviving at a much higher rate than M. guttatus. Higher seed set in M. nasutus was mediated, at least in part, by a plastic increase in the rate of late-stage development (i.e., fruit maturation), consistent with the ability of this species to inhabit more ephemeral habitats in the field. Our results suggest adaptation to water availability may be an important factor in species maintenance of these Mimulus taxa in sympatry.

Rapid experimental evolution of reproductive isolation from a single natural population

Ecological speciation occurs when local adaptation generates reproductive isolation as a by-product of natural selection. Although ecological speciation is a fundamental source of diversification, the mechanistic link between natural selection and reproductive isolation remains poorly understood, especially in natural populations. Here, we show that experimental evolution of parasite body size over 4 y (approximately 60 generations) leads to reproductive isolation in natural populations of feather lice on birds. When lice are transferred to pigeons of different sizes, they rapidly evolve differences in body size that are correlated with host size. These differences in size trigger mechanical mating isolation between lice that are locally adapted to the different sized hosts. Size differences among lice also influence the outcome of competition between males for access to females. Thus, body size directly mediates reproductive isolation through its influence on both intersexual compatibility and intrasexual competition. Our results confirm that divergent natural selection acting on a single phenotypic trait can cause reproductive isolation to emerge from a single natural population in real time.

A hot topic: the genetics of adaptation to geothermal vents in Mimulus guttatus

Identifying the individual loci and mutations that underlie adaptation to extreme environments has long been a goal of evolutionary biology. However, finding the genes that underlie adaptive traits is difficult for several reasons. First, because many traits and genes evolve simultaneously as populations diverge, it is difficult to disentangle adaptation from neutral demographic processes. Second, finding the individual loci involved in any trait is challenging given the respective limitations of quantitative and population genetic methods. In this issue of Molecular Ecology, Hendrick et al. (2016) overcome these difficulties and determine the genetic basis of microgeographic adaptation between geothermal vent and nonthermal populations of Mimulus guttatus in Yellowstone National Park. The authors accomplish this by combining population and quantitative genetic techniques, a powerful, but labour-intensive, strategy for identifying individual causative adaptive loci that few studies have used (Stinchcombe & Hoekstra ). In a previous common garden experiment (Lekberg et al. 2012), thermal M. guttatus populations were found to differ from their closely related nonthermal neighbours in various adaptive phenotypes including trichome density. Hendrick et al. (2016) combine quantitative trait loci (QTL) mapping, population genomic scans for selection and admixture mapping to identify a single genetic locus underlying differences in trichome density between thermal and nonthermal M. guttatus. The candidate gene, R2R3 MYB, is homologous to genes involved in trichome development across flowering plants. The major trichome QTL, Tr14, is also involved in trichome density differences in an independent M. guttatus population comparison (Holeski et al. 2010) making this an example of parallel genetic evolution.

Monkeyflowers (Mimulus): new model for plant developmental genetics and evo-devo

Contents Summary 694 I. Introduction 694 II. The system 695 III. Regulation of carotenoid pigmentation 695 IV. Formation of periodic pigmentation patterns 696 V. Developmental genetics of corolla tube formation and elaboration 697 VI. Molecular basis of floral trait variation underlying pollinator shift 698 VII. Outlook 699 Acknowledgements 699 References 699 SUMMARY: Monkeyflowers (Mimulus) have long been recognized as a classic ecological and evolutionary model system. However, only recently has it been realized that this system also holds great promise for studying the developmental genetics and evo-devo of important plant traits that are not found in well-established model systems such as Arabidopsis. Here, I review recent progress in four different areas of plant research enabled by this new model, including transcriptional regulation of carotenoid biosynthesis, formation of periodic pigmentation patterns, developmental genetics of corolla tube formation and elaboration, and the molecular basis of floral trait divergence underlying pollinator shift. These examples suggest that Mimulus offers ample opportunities to make exciting discoveries in plant development and evolution.

Patterns of Transmission Ratio Distortion in Interspecific Lettuce Hybrids Reveal a Sex-Independent Gametophytic Barrier

Interspecific crosses can result in progeny with reduced vitality or fertility due to genetic incompatibilities between species, a phenomenon known as hybrid incompatibility (HI). HI is often caused by a bias against deleterious allele combinations, which results in transmission ratio distortion (TRD). Here, we determined the genome-wide distribution of HI between wild lettuce, Lactuca saligna, and cultivated lettuce, L. sativa, in a set of backcross inbred lines (BILs) with single introgression segments from L. saligna introgressed into a L. sativa genetic background. Almost all BILs contained an introgression segment in a homozygous state except a few BILs, for which we were able to obtain only a single heterozygous introgression. Their inbred progenies displayed severe TRD with a bias toward the L. sativa allele and complete nontransmission of the homozygous L. saligna introgression, i.e., absolute HI. These HI might be caused by deleterious heterospecific allele combinations at two loci. We used an multilocus segregating interspecific F2 population to identify candidate conspecific loci that can nullify the HI in BILs. Segregation analysis of developed double-introgression progenies showed nullification of three HI and proved that these HI are explained by nuclear pairwise incompatibilities. One of these digenic HI showed 29% reduced seed set and its pattern of TRD pointed to a sex-independent gametophytic barrier. Namely, this HI was caused by complete nontransmission of one heterospecific allele combination at the haploid stage, surprisingly in both male and female gametophytes. Our study shows that two-locus incompatibility systems contribute to reproductive barriers among Lactuca species.

Cloning and characterization of a Mimulus lewisii NPR1 gene involved in regulating plant resistance to Rhizoctonia solani

The monkey flower Mimulus lewisii is a new emerging model plant for the study in corolla tube formation, pigmentation patterns and pollinator selection, etc. However, the cultivation and management of this plant are difficult due to its susceptibility to a wide range of pathogens and the lack of rigid varieties with high levels of resistance to pathogens. In this regard, genetic engineering is a promising tool that may possibly allow us to enhance the M. lewisii disease resistance against pathogens. Here, we reported the isolation and characterization of non-expressor of pathogenesis related gene 1 (NPR1) gene from M. lewisii. The phylogenetic tree constructed based on the deduced sequence of MlNPR1 with homologs from other species revealed that MlNPR1 grouped together with other known NPR1 proteins of Scrophulariaceae family, and was nearest to Mimulus guttatus. Furthermore, expression analysis showed that MlNPR1 was upregulated after SA treatment and fungal infection. To understand the defensive role of this gene, we overexpressed MlNPR1 in M. lewisii. The transgenic lines showed slight phenotypic abnormalities, but constitutive expression of MlNPR1 activates defense signaling pathways by priming the expression of antifungal PR genes. Moreover, MlNPR1 transgenic lines showed enhanced resistance to Rhizoctonia solani there was delay in symptoms and reduced disease severity than non-transgenic plants. Altogether, the present study suggests that increasing the expression level of MlNPR1 may be a promising approach for development of monkey flower cultivars with enhanced resistance to diseases.