Selfish chromosomal drive shapes recent centromeric histone evolution in monkeyflowers.

Centromeres are essential mediators of chromosomal segregation, but both centromeric DNA sequences and associated kinetochore proteins are paradoxically diverse across species. The selfish centromere model explains rapid evolution by both components via an arms-race scenario: centromeric DNA variants drive by distorting chromosomal transmission in female meiosis and attendant fitness costs select on interacting proteins to restore Mendelian inheritance. Although it is clear than centromeres can drive and that drive often carries costs, female meiotic drive has not been directly linked to selection on kinetochore proteins in any natural system. Here, we test the selfish model of centromere evolution in a yellow monkeyflower (Mimulus guttatus) population polymorphic for a costly driving centromere (D). We show that the D haplotype is structurally and genetically distinct and swept to a high stable frequency within the past 1500 years. We use quantitative genetic mapping to demonstrate that context-dependence in the strength of drive (from near-100% D transmission in interspecific hybrids to near-Mendelian in within-population crosses) primarily reflects variable vulnerability of the non-driving competitor chromosomes, but also map an unlinked modifier of drive coincident with kinetochore protein Centromere-specific Histone 3 A (CenH3A). Finally, CenH3A exhibits a recent (<1000 years) selective sweep in our focal population, implicating local interactions with D in ongoing adaptive evolution of this kinetochore protein. Together, our results demonstrate an active co-evolutionary arms race between DNA and protein components of the meiotic machinery in Mimulus, with important consequences for individual fitness and molecular divergence.

Inbreeding depression contributes to the maintenance of habitat segregation between closely related monkeyflower species.

Incompletely reproductively isolated species often segregate into different microhabitats, even when they are able to survive and reproduce in both habitats. Longer term evolutionary factors may contribute to this lack of cross-habitat persistence. When reproductive interference reduces immigrant fitness, assortative mating, including self-fertilization, increases immigrants' fitness in a single generation, but longer term, inbreeding depression may reduce the chance of population persistence. Two California monkeyflower species repeatedly segregate into drier and wetter areas in their zone of sympatry. To test whether inbreeding depression may contribute to the maintenance of this segregation pattern, we transplanted outbred and successively inbred Mimulus guttatus and Mimulus nudatus into their native habitats and heterospecific habitats. We measured germination, survival, and seed set and found that recurrent selfing reduced all aspects of fitness in both species, most strongly in foreign habitats. A simulation model, parameterized from the transplant experiment, found that inbreeding reduced fitness to such an extent that sequentially inbred populations of either species would be unable to persist in heterospecific-occupied habitats in the absence of continued gene flow. These results demonstrate that individual immigrants are unlikely to form persistent populations and thus, inbreeding depression contributes to the absence of fine-scale coexistence in this species pair.

Evolution of multiple postzygotic barriers between species of the Mimulus tilingii complex.

Species are often defined by their ability to interbreed (i.e., Biological Species Concept), but determining how and why reproductive isolation arises between new species can be challenging. In the Mimulus tilingii species complex, three species (M. caespitosa, M. minor, and M. tilingii) are largely allopatric and grow exclusively at high elevations (>2000 m). The extent to which geographic separation has shaped patterns of divergence among the species is not well understood. In this study, we determined that the three species are morphologically and genetically distinct, yet recently diverged. Additionally, we performed reciprocal crosses within and between the species and identified several strong postzygotic reproductive barriers, including hybrid seed inviability, F1 hybrid necrosis, and F1 hybrid male and female sterility. In this study, such postzygotic barriers are so strong that a cross between any species pair in the M. tilingii complex would cause nearly complete reproductive isolation. We consider how geographical and topographical patterns may have facilitated the evolution of several postzygotic barriers and contributed to speciation of closely related members within the M. tilingii species complex.

Widespread selection and gene flow shape the genomic landscape during a radiation of monkeyflowers.

Speciation genomic studies aim to interpret patterns of genome-wide variation in light of the processes that give rise to new species. However, interpreting the genomic "landscape" of speciation is difficult, because many evolutionary processes can impact levels of variation. Facilitated by the first chromosome-level assembly for the group, we use whole-genome sequencing and simulations to shed light on the processes that have shaped the genomic landscape during a radiation of monkeyflowers. After inferring the phylogenetic relationships among the 9 taxa in this radiation, we show that highly similar diversity (π) and differentiation (FST) landscapes have emerged across the group. Variation in these landscapes was strongly predicted by the local density of functional elements and the recombination rate, suggesting that the landscapes have been shaped by widespread natural selection. Using the varying divergence times between pairs of taxa, we show that the correlations between FST and genome features arose almost immediately after a population split and have become stronger over time. Simulations of genomic landscape evolution suggest that background selection (BGS; i.e., selection against deleterious mutations) alone is too subtle to generate the observed patterns, but scenarios that involve positive selection and genetic incompatibilities are plausible alternative explanations. Finally, tests for introgression among these taxa reveal widespread evidence of heterogeneous selection against gene flow during this radiation. Combined with previous evidence for adaptation in this system, we conclude that the correlation in FST among these taxa informs us about the processes contributing to adaptation and speciation during a rapid radiation.

A two-locus hybrid incompatibility is widespread, polymorphic, and active in natural populations of Mimulus.

Reproductive isolation, which is essential for the maintenance of species in sympatry, is often incomplete between closely related species. In these taxa, reproductive barriers must evolve within species, without being degraded by ongoing gene flow. To better understand this dynamic, we investigated the frequency and geographic distribution of alleles underlying a two-locus, hybrid lethality system between naturally hybridizing species of monkeyflower (Mimulus guttatus and M. nasutus). We found that M. guttatus typically carries hybrid lethality alleles at one locus (hl13) and M. nasutus typically carries hybrid lethality alleles at the other locus (hl14). As a result, natural hybrids carry incompatible alleles at both loci, and express hybrid lethality in later generations. We also discovered considerable polymorphism at both hl13 and hl14 within both species. For M. guttatus, polymorphism at both loci occurs within populations, meaning that incompatible allele pairings likely arise through intraspecific gene flow. Genetic variation at markers linked to hl13 and hl14 suggest that introgression from M. nasutus is the primary driver of this polymorphism within M. guttatus. Additionally, patterns of introgression at the two hybrid lethality loci suggest that natural selection eliminates incompatible allele pairings, suggesting that even weak reproductive barriers might promote genomic divergence between species.

Losing one's touch: Evolution of the touch-sensitive stigma in the Mimulus guttatus species complex.

PREMISE OF THE STUDY: The stigmas of several species are touch sensitive and respond to pressure by closing. Previous research suggests that stigma closure could prevent self pollination within a flower during a pollinator's visit or enhance male function by increasing pollen export. Both factors could be favored in outcrossers, and neither would be beneficial in selfers. METHODS: We investigated variation in stigma-closing and the duration of closure in annual and perennial populations of the variable species Mimulus guttatus and whether four closely related selfing species (M. cupriphilus, M. laciniatus, M. nasutus, and M. pardalis) have lost their touch sensitivity. We grew plants in a controlled environment and performed experiments with and without the addition of pollen to the stigma. KEY RESULTS: In M. guttatus, the speed of stigma-closing was rapid and unaffected by the deposition of pollen. Populations varied significantly in closing speed, which may reflect their geographic location. For annual populations only, anther-stigma separation significantly affected closing speed. Also, stigmas that closed quickly stayed closed longer, and stigmas that received pollen remained closed longer. Finally, in the selfing species, stigma-closing was more variable; some populations have entirely lost the ability to respond to touch. CONCLUSIONS: We discuss our results in the context of traits that promote outcrossing and traits that are under selection during the evolution of selfing. This is the first characterization of variation in touch responses across multiple populations within a species and the first to demonstrate the loss of touch sensitivity in selfing lineages.

Reproductive isolation and introgression between sympatric Mimulus species.

Incompletely isolated species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain distinct species in the face of ongoing gene flow. Here, we use field surveys and reduced representation sequencing to characterize the patterns of reproductive isolation, admixture and genomic divergence between populations of the outcrossing wildflower Mimulus guttatus and selfing M. nasutus. Focusing on a single site where these two species have come into secondary contact, we find that phenological isolation is strong, although incomplete, and is likely driven by divergence in response to photoperiod. In contrast to previous field studies, which have suggested that F1 -hybrid formation might be rare, we discover patterns of genomic variation consistent with ongoing introgression. Strikingly, admixed individuals vary continuously from highly admixed to nearly pure M. guttatus, demonstrating ongoing hybridization and asymmetric introgression from M. nasutus into M. guttatus. Patterns of admixture and divergence across the genome show that levels of introgression are more variable than expected by chance. Some genomic regions show a reduced introgression, including one region that overlaps a critical photoperiod QTL, whereas other regions show elevated levels of interspecific gene flow. In addition, we observe a genome-wide negative relationship between absolute divergence and the local recombination rate, potentially indicating natural selection against M. nasutus ancestry in M. guttatus genetic backgrounds. Together, our results suggest that Mimulus speciation is both ongoing and dynamic and that a combination of divergence in phenology and mating system, as well as selection against interspecific alleles, likely maintains these sympatric species.

Comparative Transcriptomics Indicates a Role for SHORT VEGETATIVE PHASE (SVP) Genes in Mimulus guttatus Vernalization Response.

The timing of reproduction in response to variable environmental conditions is critical to plant fitness, and is a major driver of taxon differentiation. In the yellow monkey flower, Mimulus guttatus, geographically distinct North American populations vary in their photoperiod and chilling (vernalization) requirements for flowering, suggesting strong local adaptation to their surroundings. Previous analyses revealed quantitative trait loci (QTL) underlying short-day mediated vernalization responsiveness using two annual M. guttatus populations that differed in their vernalization response. To narrow down candidate genes responsible for this variation, and to reveal potential downstream genes, we conducted comparative transcriptomics and quantitative PCR (qPCR) in shoot apices of parental vernalization responsive IM62, and unresponsive LMC24 inbred lines grown under different photoperiods and temperatures. Our study identified several metabolic, hormone signaling, photosynthetic, stress response, and flowering time genes that are differentially expressed between treatments, suggesting a role for their protein products in short-day-mediated vernalization responsiveness. Only a small subset of these genes intersected with candidate genes from the previous QTL study, and, of the main candidates tested with qPCR under nonpermissive conditions, only SHORT VEGETATIVE PHASE (SVP) gene expression met predictions for a population-specific short-day-repressor of flowering that is repressed by cold.

Genomic studies on the nature of species: adaptation and speciation in Mimulus.

Evolutionary biology is in an exciting era, in which powerful genomic tools make the answers accessible to long-standing questions about variation, adaptation and speciation. The availability of a suite of genomic resources, a shared knowledge base and a long history of study have made the phenotypically diverse plant genus Mimulus an important system for understanding ecological and evolutionary processes. An international Mimulus Research Meeting was held at Duke University in June 2014 to discuss developments in ecological and evolutionary genetic studies in Mimulus. Here, we report major recent discoveries presented at the meeting that use genomic approaches to advance our understanding of three major themes: the parallel genetic basis of adaptation; the ecological genomics of speciation; and the evolutionary significance of structural genetic variation. We also suggest future research directions for studies of Mimulus and highlight challenges faced when developing new ecological and evolutionary model systems.

Chromosomal rearrangements directly cause underdominant F1 pollen sterility in Mimulus lewisii-Mimulus cardinalis hybrids.

Chromosomal rearrangements can contribute to the evolution of postzygotic reproductive isolation directly, by disrupting meiosis in F1 hybrids, or indirectly, by suppressing recombination among genic incompatibilities. Because direct effects of rearrangements on fertility imply fitness costs during their spread, understanding the mechanism of F1 hybrid sterility is integral to reconstructing the role(s) of rearrangements in speciation. In hybrids between monkeyflowers Mimulus cardinalis and Mimulus lewisii, rearrangements contain all quantitative trait loci (QTLs) for both premating barriers and pollen sterility, suggesting that they may have facilitated speciation in this model system. We used artificial chromosome doubling and comparative mapping to test whether heterozygous rearrangements directly cause underdominant male sterility in M. lewisii-M. cardinalis hybrids. Consistent with a direct chromosomal basis for hybrid sterility, synthetic tetraploid F1 s showed highly restored fertility (83.4% pollen fertility) relative to diploids F1 s (36.0%). Additional mapping with Mimulus parishii-M. cardinalis and M. parishii-M. lewisii hybrids demonstrated that underdominant male sterility is caused by one M. lewisii specific and one M. cardinalis specific reciprocal translocation, but that inversions had no direct effects on fertility. We discuss the importance of translocations as causes of reproductive isolation, and consider models for how underdominant rearrangements spread and fix despite intrinsic fitness costs.

The evolution of environmental tolerance and range size: a comparison of geographically restricted and widespread Mimulus.

The geographic ranges of closely related species can vary dramatically, yet we do not fully grasp the mechanisms underlying such variation. The niche breadth hypothesis posits that species that have evolved broad environmental tolerances can achieve larger geographic ranges than species with narrow environmental tolerances. In turn, plasticity and genetic variation in ecologically important traits and adaptation to environmentally variable areas can facilitate the evolution of broad environmental tolerance. We used five pairs of western North American monkeyflowers to experimentally test these ideas by quantifying performance across eight temperature regimes. In four species pairs, species with broader thermal tolerances had larger geographic ranges, supporting the niche breadth hypothesis. As predicted, species with broader thermal tolerances also had more within-population genetic variation in thermal reaction norms and experienced greater thermal variation across their geographic ranges than species with narrow thermal tolerances. Species with narrow thermal tolerance may be particularly vulnerable to changing climatic conditions due to lack of plasticity and insufficient genetic variation to respond to novel selection pressures. Conversely, species experiencing high variation in temperature across their ranges may be buffered against extinction due to climatic changes because they have evolved tolerance to a broad range of temperatures.

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.

Comparative linkage maps suggest that fission, not polyploidy, underlies near-doubling of chromosome number within monkeyflowers (Mimulus; Phrymaceae).

Changes in chromosome number and structure are important contributors to adaptation, speciation and macroevolution. In flowering plants, polyploidy and subsequent reductions in chromosome number by fusion are major sources of chromosomal evolution, but chromosome number increase by fission has been relatively unexplored. Here, we use comparative linkage mapping with gene-based markers to reconstruct chromosomal synteny within the model flowering plant genus Mimulus (monkeyflowers). Two sections of the genus with haploid numbers ≥ 14 have been inferred to be relatively recent polyploids because they are phylogenetically nested within numerous taxa with low base numbers (n=8-10). We combined multiple data sets to build integrated genetic maps of the M. guttatus species complex (section Simiolus, n=14) and the M. lewisii group (section Erythranthe; n=8), and then aligned the two integrated maps using >100 shared markers. We observed strong segmental synteny between M. lewisii and M. guttatus maps, with essentially 1-to-1 correspondence across each of 16 chromosomal blocks. Assuming that the M. lewisii (and widespread) base number of 8 is ancestral, reconstruction of 14 M. guttatus chromosomes requires at least eight fission events (likely shared by Simiolus and sister section Paradanthus (n=16)), plus two fusion events. This apparent burst of fission in the yellow monkeyflower lineages raises new questions about mechanisms and consequences of chromosomal fission in plants. Our comparative maps also provide insight into the origins of a chromosome exhibiting centromere-associated female meiotic drive and create a framework for transferring M. guttatus genome resources across the entire genus.