Divergent population structure and climate associations of a chromosomal inversion polymorphism across the Mimulus guttatus species complex

Chromosomal rearrangement polymorphisms are common and increasingly found to be associated with adaptive ecological divergence and speciation. Rearrangements, such as inversions, reduce recombination in heterozygous individuals and thus can protect favourable allelic combinations at linked loci, facilitating their spread in the presence of gene flow. Recently, we identified a chromosomal inversion polymorphism that contributes to ecological adaptation and reproductive isolation between annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus. Here we evaluate the population genetic structure of this inverted region in comparison with the collinear regions of the genome across the M. guttatus species complex. We tested whether annual and perennial M. guttatus exhibit different patterns of divergence for loci in the inverted and noninverted regions of the genome. We then evaluated whether there are contrasting climate associations with these genomic regions through redundancy analysis. We found that the inversion exhibits broadly different patterns of divergence among annual and perennial M. guttatus and is associated with environmental variation across population accessions. This study is the first widespread population genetic survey of the diversity of the M. guttatus species complex. Our findings contribute to a greater understanding of morphological, ecological, and genetic evolutionary divergence across this highly diverse group of closely related ecotypes and species. Finally, understanding species relationships among M. guttatus sp. has hitherto been stymied by accumulated evidence of substantial gene flow among populations as well as designated species. Nevertheless, our results shed light on these relationships and provide insight into adaptation in life history traits within the complex.     

Niche partitioning between close relatives suggests trade‐offs between adaptation to local environments and competition

Niche partitioning among close relatives may reflect trade‐offs underlying species divergence and coexistence (e.g., between stress tolerance and competitive ability). We quantified the effects of habitat and congeneric species interactions on fitness for two closely related herbaceous plant species, Mimulus guttatus and Mimulus laciniatus, in three common habitat types within their sympatric range. Drought stress strongly reduced survival of M. guttatus in fast‐drying seeps occupied by M. laciniatus, suggesting that divergent habitat adaptation maintains this niche boundary. However, neither seedling performance nor congeneric competition explained the absence of M. laciniatus from shady streams where M. guttatus thrives. M. laciniatus may be excluded from this habitat by competition with other species in the community or mature M. guttatus. Species performance and competitive ability were similar in sympatric meadows where plant community stature and the growing season length are intermediate between seeps and streams. Stochastic effects (e.g., dispersal among habitats or temporal variation) may contribute to coexistence in this habitat. Habitat adaptation, species interactions, and stochastic mechanisms influence sympatric distributions for these recently diverged species.     

REPRODUCTIVE ASSURANCE MECHANISMS IN THREE TAXA OF THE MIMULUS GUTTATUS COMPLEX (SCROPHULARIACEAE)

The evolution of inbreeding in plants has often been attributed to selection for the ability to set seed in the absence of mates or pollinators. Mechanisms of reproductive assurance in five populations of mixed mating Mimulus guttatus, three populations of inbreeding M. platycalyx, and two populations of inbreeding M. nasutus were examined in a pollinator-free greenhouse. Reproductive assurance was manifested in all populations by autofertility, vegetative reproduction, or both. The inbreeding taxa had significantly greater levels of autofertility and less vegetative reproduction. Three modes of autofertility were identified: 1) due to corolla abscission only, occurring in three M. guttatus populations; 2) due to both corolla abscission and direct anther-stigma contact by curling of the lower stigmatic lobe into the anthers, occurring in two M. guttatus populations; and 3) direct stigma-anther contact by stigma curling alone prior to corolla abscission, found in each M. platycalyx and M. nasutus population. Stigma-anther distance and its interaction with stigma curling contributed to differences in autofertility among populations. Significant levels of intrapopulation quantitative genetic variation were found for seven of ten traits examined; average levels were similar between inbreeding and mixed mating populations. Genetic variation within populations for autofertility per se was not detected, but significant levels controlling stigma-anther distance were found in two M. guttatus populations. These results show that evolution of inbreeding by natural selection for reproductive assurance is possible in Mimulus, and illustrate the complex changes in floral dynamics and morphology it may involve.     

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 2008 ). 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.     

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

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

CROSS ABILITY OF MIMULUS GUTTATUS IN RELATION TO COMPONENTS OF GENE FIXATION

The joint effects of parental gene fixation and consanguinity of mates upon the fitness of matings was examined in Mimulus guttatus. Plants from four populations were crossed at five levels of genetic relatedness, and five viability characters were scored in progeny. Parental gene fixation at 12 polymorphic allozyme loci was partitioned into local, subpopulation, and population components. A model is proposed wherein parental gene fixation influences distance-dependent crossing success. At a fixed locus, inbreeding is favored if natural selection caused allele fixation, or is disfavored if gene fixation was random. The distance between mates required to eliminate gene fixation depends upon patch size of fixation. When selective fixation and patch size differ among loci, an optimal crossing distance is possible. In M. guttatus, progeny viability generally decreased with greater relatedness between mates, but this decrease was often heterogeneous among populations. The highest pollen viability and the lowest seed set were found at an intermediate relatedness between mates. To determine whether parental gene fixation influences these crossing patterns, a type of mutational-load analysis was performed. Progeny fitness was regressed on parent F and fitness estimated at F = 1. This was done for each component of F, for a) crosses that maintain gene fixation and b) crosses that eliminate gene fixation. A multiplicative, composite measure of fitness indicates that, in M. guttatus, genes fixed during local or population differentiation favor outbreeding, while genes fixed during subpopulation differentiation favor inbreeding. This predicts that random mating within subpopulations confers highest progeny fitness, exclusive of between-population matings. However, predictions did not fit the observed patterns of crossing success very well, perhaps because gene fixation was relatively low or was not adequately measured at loci influencing fitness.     

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.     

Speciation and Introgression between Mimulus nasutus and Mimulus guttatus

Mimulus guttatus and M. nasutus are an evolutionary and ecological model sister species pair differentiated by ecology, mating system, and partial reproductive isolation. Despite extensive research on this system, the history of divergence and differentiation in this sister pair is unclear. We present and analyze a population genomic data set which shows that M. nasutus budded from a central Californian M. guttatus population within the last 200 to 500 thousand years. In this time, the M. nasutus genome has accrued genomic signatures of the transition to predominant selfing, including an elevated proportion of nonsynonymous variants, an accumulation of premature stop codons, and extended levels of linkage disequilibrium. Despite clear biological differentiation, we document genomic signatures of ongoing, bidirectional introgression. We observe a negative relationship between the recombination rate and divergence between M. nasutus and sympatric M. guttatus samples, suggesting that selection acts against M. nasutus ancestry in M. guttatus.     

Natural variation for drought-response traits in the Mimulus guttatus species complex

Soil moisture is a key factor affecting plant abundance and distribution, both across and within species. In response to water limitation, plants have evolved numerous morphological, physiological, and phenological adaptations. In both well-watered and water-limited conditions, we identified considerable natural variation in drought-related whole-plant and leaf-level traits among closely related members of the Mimulus guttatus species complex that occupy a diversity of habitats in the field. The self-fertilizing Mimulus nasutus and serpentine-endemic Mimulus nudatus demonstrated the overall greatest tolerance to soil water limitation, exhibiting the smallest reduction in seed set relative to well-watered conditions. This may be due in part to early flowering, faster fruit development, and low stomatal density. In contrast, flowering of coastal M. guttatus was so delayed that it precluded any seed production in water-limited conditions. This range of phenotypic responses to soil water deficit in Mimulus, coupled with developing genomic resources, holds considerable promise for identifying genomic variation responsible for adaptive responses to soil water availability.     

The extent and genetic basis of phenotypic divergence in life history traits in Mimulus guttatus

Differential natural selection acting on populations in contrasting environments often results in adaptive divergence in multivariate phenotypes. Multivariate trait divergence across populations could be caused by selection on pleiotropic alleles or through many independent loci with trait‐specific effects. Here, we assess patterns of association between a suite of traits contributing to life history divergence in the common monkey flower, Mimulus guttatus, and examine the genetic architecture underlying these correlations. A common garden survey of 74 populations representing annual and perennial strategies from across the native range revealed strong correlations between vegetative and reproductive traits. To determine whether these multitrait patterns arise from pleiotropic or independent loci, we mapped QTLs using an approach combining high‐throughput sequencing with bulk segregant analysis on a cross between populations with divergent life histories. We find extensive pleiotropy for QTLs related to flowering time and stolon production, a key feature of the perennial strategy. Candidate genes related to axillary meristem development colocalize with the QTLs in a manner consistent with either pleiotropic or independent QTL effects. Further, these results are analogous to previous work showing pleiotropy‐mediated genetic correlations within a single population of M. guttatus experiencing heterogeneous selection. Our findings of strong multivariate trait associations and pleiotropic QTLs suggest that patterns of genetic variation may determine the trajectory of adaptive divergence.     

Differential adaptation to a harsh granite outcrop habitat between sympatric Mimulus species

Understanding which environmental variables and traits underlie adaptation to harsh environments is difficult because many traits evolve simultaneously as populations or species diverge. Here, we investigate the ecological variables and traits that underlie Mimulus laciniatus’ adaptation to granite outcrops compared to its sympatric, mesic‐adapted progenitor, Mimulus guttatus. We use fine‐scale measurements of soil moisture and herbivory to examine differences in selective forces between the species’ habitats, and measure selection on flowering time, flower size, plant height, and leaf shape in a reciprocal transplant using M. laciniatus × M. guttatus F₄ hybrids. We find that differences in drought and herbivory drive survival differences between habitats, that M. laciniatus and M. guttatus are each better adapted to their native habitat, and differential habitat selection on flowering time, plant stature, and leaf shape. Although early flowering time, small stature, and lobed leaf shape underlie plant fitness in M. laciniatus’ seasonally dry environment, increased plant size is advantageous in a competitive mesic environment replete with herbivores like M. guttatus’. Given that we observed divergent selection between habitats in the direction of species differences, we conclude that adaptation to different microhabitats is an important component of reproductive isolation in this sympatric species pair.     

Leaf shape evolution has a similar genetic architecture in three edaphic specialists within the Mimulus guttatus species complex

Background and Aims The genetic basis of leaf shape has long interested botanists because leaf shape varies extensively across the plant kingdom and this variation is probably adaptive. However, knowledge of the genetic architecture of leaf shape variation in natural populations remains limited. This study examined the genetic architecture of leaf shape diversification among three edaphic specialists in the Mimulus guttatus species complex. Lobed and narrow leaves have evolved from the entire, round leaves of M. guttatus in M. laciniatus, M. nudatus and a polymorphic serpentine M. guttatus population (M2L).Methods Bulk segregant analysis and next-generation sequencing were used to map quantitative trait loci (QTLs) that underlie leaf shape in an M. laciniatus × M. guttatus F₂ population. To determine whether the same QTLs contribute to leaf shape variation in M. nudatus and M2L, F₂s from M. guttatus × M. nudatus and lobed M2L × unlobed M. guttatus crosses were genotyped at QTLs from the bulk segregant analysis.Key Results Narrow and lobed leaf shapes in M. laciniatus, M. nudatus and M. guttatus are controlled by overlapping genetic regions. Several promising leaf shape candidate genes were found under each QTL.Conclusions The evolution of divergent leaf shape has taken place multiple times in the M. guttatus species complex and is associated with the occupation of dry, rocky environments. The genetic architecture of elongated and lobed leaves is similar across three species in this group. This may indicate that parallel genetic evolution from standing variation or new mutations is responsible for the putatively adaptive leaf shape variation in Mimulus.     

Floral volatile alleles can contribute to pollinator‐mediated reproductive isolation in monkeyflowers (Mimulus)

Pollinator‐mediated reproductive isolation is a major factor in driving the diversification of flowering plants. Studies of floral traits involved in reproductive isolation have focused nearly exclusively on visual signals, such as flower color. The role of less obvious signals, such as floral scent, has been studied only recently. In particular, the genetics of floral volatiles involved in mediating differential pollinator visitation remains unknown. The bumblebee‐pollinated Mimulus lewisii and hummingbird‐pollinated Mimulus cardinalis are a model system for studying reproductive isolation via pollinator preference. We have shown that these two species differ in three floral terpenoid volatiles – d ‐limonene, β‐myrcene, and E‐β‐ocimene – that are attractive to bumblebee pollinators. By genetic mapping and in vitro analysis of enzyme activity we demonstrate that these interspecific differences are consistent with allelic variation at two loci, LIMONENE‐MYRCENE SYNTHASE (LMS) and OCIMENE SYNTHASE (OS). Mimulus lewisii LMS (MlLMS) and OS (MlOS) are expressed most strongly in floral tissue in the last stages of floral development. Mimulus cardinalis LMS (McLMS) is weakly expressed and has a nonsense mutation in exon 3. Mimulus cardinalis OS (McOS) is expressed similarly to MlOS, but the encoded McOS enzyme produces no E‐β‐ocimene. Recapitulating the M. cardinalis phenotype by reducing the expression of MlLMS by RNA interference in transgenic M. lewisii produces no behavioral difference in pollinating bumblebees; however, reducing MlOS expression produces a 6% decrease in visitation. Allelic variation at the OCIMENE SYNTHASE locus is likely to contribute to differential pollinator visitation, and thus promote reproductive isolation between M. lewisii and M. cardinalis. OCIMENE SYNTHASE joins a growing list of ‘speciation genes’ (‘barrier genes’) in flowering plants.     

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

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

Digest: Linking coordinated shifts in plant resource allocation to a chromosomal inversion*

Local adaptation in plants often requires coordinated shifts among resources. Lowry et al. provide evidence for physiological and genomic mechanisms underpinning adaptive shifts in yellow monkeyflower (Mimulus guttatus), such as the transition between annual and perennial life histories. In M. guttatus, differential activity of gibberellins, governed partially by a chromosomal inversion, is responsible for shifts between growth, reproduction, and herbivore defense (secondary compound production).     

Inbreeding depression in an asexual population of Mimulus guttatus

The reproductive mechanism, that is whether an organism outcrosses, selfs or asexually reproduces, has a substantial impact on the amount and pattern of genetic variation. In this study, we estimate genetic variation and genetic load for a predominately asexual population of Mimulus guttatus, and then compare our results to other studies of predominately sexually reproducing (outcrossing and selfing) populations of M. guttatus. The asexual population had low levels of heterozygosity (H_e = 0.03) and low (but significantly non‐zero) inbreeding load, especially when compared with other M. guttatus populations. This differs greatly from the sexual populations of Mimulus that display substantial inbreeding depression. We discuss a variety of reasons why we see such low load in this study and suggest future research projects to further explore the questions.     

Speciation on a local geographic scale: the evolution of a rare rock outcrop specialist in Mimulus

Speciation can occur on both large and small geographical scales. In plants, local speciation, where small populations split off from a large-ranged progenitor species, is thought to be the dominant mode, yet there are still few examples to verify speciation has occurred in this manner. A recently described morphological species in the yellow monkey flowers, Mimulus filicifolius, is an excellent candidate for local speciation because of its highly restricted geographical range. Mimulus filicifolius was formerly identified as a population of M. laciniatus due to similar lobed leaf morphology and rocky outcrop habitat. To investigate whether M. filicifolius is genetically divergent and reproductively isolated from M. laciniatus, we examined patterns of genetic diversity in ten nuclear and eight microsatellite loci, and hybrid fertility in M. filicifolius and its purported close relatives: M. laciniatus, M. guttatus and M. nasutus. We found that M. filicifolius is genetically divergent from the other species and strongly reproductively isolated from M. laciniatus. We conclude that M. filicifolius is an independent rock outcrop specialist despite being morphologically and ecologically similar to M. laciniatus, and that its small geographical range nested within other wide-ranging members of the M. guttatus species complex is consistent with local speciation.     

Spatial population genetic structure of a bacterial parasite in close coevolution with its host

Knowledge of a species’ population genetic structure can provide insight into fundamental ecological and evolutionary processes including gene flow, genetic drift and adaptive evolution. Such inference is of particular importance for parasites, as an understanding of their population structure can illuminate epidemiological and coevolutionary dynamics. Here, we describe the population genetic structure of the bacterium Pasteuria ramosa, a parasite that infects planktonic crustaceans of the genus Daphnia. This system has become a model for investigations of host–parasite interactions and represents an example of coevolution via negative frequency‐dependent selection (aka “Red Queen” dynamics). To sample P. ramosa, we experimentally infected a panel of Daphnia hosts with natural spore banks from the sediments of 25 ponds throughout much of the species range in Europe and western Asia. Using 12 polymorphic variable number tandem repeat loci (VNTR loci), we identified substantial genetic diversity, both within and among localities, that was structured geographically among ponds. Genetic diversity was also structured among host genotypes within ponds, although this pattern varied by locality, with P. ramosa at some localities partitioned into distinct host‐specific lineages, and other localities where recombination had shuffled genetic variation among different infection phenotypes. Across the sample range, there was a pattern of isolation by distance, and principal components analysis coupled with Procrustes rotation identified congruence between patterns of genetic variation and geography. Our findings support the hypothesis that Pasteuria is an endemic parasite coevolving closely with its host. These results provide important context for previous studies of this model system and inform hypotheses for future research.     

Evidence of Natural Selection Acting on a Polymorphic Hybrid Incompatibility Locus in Mimulus

As a common cause of reproductive isolation in diverse taxa, hybrid incompatibilities are fundamentally important to speciation. A key question is which evolutionary forces drive the initial substitutions within species that lead to hybrid dysfunction. Previously, we discovered a simple genetic incompatibility that causes nearly complete male sterility and partial female sterility in hybrids between the two closely related yellow monkeyflower species Mimulus guttatus and M. nasutus. In this report, we fine map the two major incompatibility loci—hybrid male sterility 1 (hms1) and hybrid male sterility 2 (hms2)—to small nuclear genomic regions (each <70 kb) that include strong candidate genes. With this improved genetic resolution, we also investigate the evolutionary dynamics of hms1 in a natural population of M. guttatus known to be polymorphic at this locus. Using classical genetic crosses and population genomics, we show that a 320-kb region containing the hms1 incompatibility allele has risen to intermediate frequency in this population by strong natural selection. This finding provides direct evidence that natural selection within plant species can lead to hybrid dysfunction between species.