skip to main content

Title: 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 wildflowerMimulus guttatusand selfingM. 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 pureM. guttatus, demonstrating ongoing hybridization and asymmetric introgression fromM. nasutusintoM. 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 photoperiodQTL, 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 againstM. nasutusancestry inM. guttatusgenetic backgrounds. Together, our results suggest thatMimulusspeciation 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.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Date Published:
Journal Name:
Molecular Ecology
Page Range / eLocation ID:
p. 2499-2517
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The genetic architecture of local adaptation has been of central interest to evolutionary biologists since the modern synthesis. In addition to classic theory on the effect size of adaptive mutations by Fisher, Kimura and Orr, recent theory addresses the genetic architecture of local adaptation in the face of ongoing gene flow. This theory predicts that with substantial gene flow between populations local adaptation should proceed primarily through mutations of large effect or tightly linked clusters of smaller effect loci. In this study, we investigate the genetic architecture of divergence in flowering time, mating system‐related traits, and leaf shape betweenMimulus laciniatusand a sympatric population of its close relativeM. guttatus. These three traits are probably involved inM. laciniatus’adaptation to a dry, exposed granite outcrop environment. Flowering time and mating system differences are also reproductive isolating barriers making them ‘magic traits’. Phenotypic hybrids in this population provide evidence of recent gene flow. Using next‐generation sequencing, we generate denseSNPmarkers across the genome and map quantitative trait loci (QTLs) involved in flowering time, flower size and leaf shape. We find that interspecific divergence in all three traits is due to fewQTLof large effect including a highly pleiotropicQTLon chromosome 8. ThisQTLregion contains the pleiotropic candidate gene TCP4 and is involved in ecologically important phenotypes in otherMimulusspecies. Our results are consistent with theory, indicating that local adaptation and reproductive isolation with gene flow should be due to few loci with large and pleiotropic effects.

    more » « less
  2. Abstract

    The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers,Mimulus guttatusandMimulus nasutus,from the sympatric Catherine Creek population. We discover that the threeM. guttatusfounders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed fromM. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among theM. guttatusfounders, two due to admixture. We find strong, genome‐wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter‐chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky‐Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of threeM. guttatuslines, we discover abundant segregating variation for hybrid incompatibilities withM. nasutus,suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities betweenMimulusspecies persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

    more » « less
  3. Abstract

    The reuse of old genetic variation can promote rapid diversification in evolutionary radiations, but in most cases, the historical events underlying this divergence are not known. For example, ancient hybridization can generate new combinations of alleles that sort into descendant lineages, potentially providing the raw material to initiate divergence. In the Mimulus aurantiacus species complex, there is evidence for widespread gene flow among members of this radiation. In addition, allelic variation in the MaMyb2 gene is responsible for differences in flower color between the closely related ecotypes of subspecies puniceus, contributing to reproductive isolation by pollinators. Previous work suggested that MaMyb2 was introgressed into the red-flowered ecotype of puniceus. However, additional taxa within the radiation have independently evolved red flowers from their yellow-flowered ancestors, raising the possibility that this introgression had a more ancient origin. In this study, we used repeated tests of admixture from whole-genome sequence data across this diverse radiation to demonstrate that there has been both ancient and recurrent hybridization in this group. However, most of the signal of this ancient introgression has been removed due to selection, suggesting that widespread barriers to gene flow are in place between taxa. Yet, a roughly 30 kb region that contains the MaMyb2 gene is currently shared only among the red-flowered taxa. Patterns of admixture, sequence divergence, and extended haplotype homozygosity across this region confirm a history of ancient hybridization, where functional variants have been preserved due to positive selection in red-flowered taxa but lost in their yellow-flowered counterparts. The results of this study reveal that selection against gene flow can reduce genomic signatures of ancient hybridization, but that historical introgression can provide essential genetic variation that facilitates the repeated evolution of phenotypic traits between lineages.

    more » « less
  4. Abstract

    Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches – quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals – to investigate one trait associated withMimulus guttatus(yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a majorQTLcausing dense leaf trichomes in thermally adapted plants to a <50‐kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independentM. guttatuspopulations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferredQTLpeak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3MYBtranscription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r2 = 0.27) with trichome phenotype in analyses of wild‐collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis ofMimulus guttatusadaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

    more » « less
  5. Burbrink, Frank (Ed.)
    Abstract In cryptic amphibian complexes, there is a growing trend to equate high levels of genetic structure with hidden cryptic species diversity. Typically, phylogenetic structure and distance-based approaches are used to demonstrate the distinctness of clades and justify the recognition of new cryptic species. However, this approach does not account for gene flow, spatial, and environmental processes that can obfuscate phylogenetic inference and bias species delimitation. As a case study, we sequenced genome-wide exons and introns to evince the processes that underlie the diversification of Philippine Puddle Frogs—a group that is widespread, phenotypically conserved, and exhibits high levels of geographically based genetic structure. We showed that widely adopted tree- and distance-based approaches inferred up to 20 species, compared to genomic analyses that inferred an optimal number of five distinct genetic groups. Using a suite of clustering, admixture, and phylogenetic network analyses, we demonstrate extensive admixture among the five groups and elucidate two specific ways in which gene flow can cause overestimations of species diversity: 1) admixed populations can be inferred as distinct lineages characterized by long branches in phylograms; and 2) admixed lineages can appear to be genetically divergent, even from their parental populations when simple measures of genetic distance are used. We demonstrate that the relationship between mitochondrial and genome-wide nuclear $p$-distances is decoupled in admixed clades, leading to erroneous estimates of genetic distances and, consequently, species diversity. Additionally, genetic distance was also biased by spatial and environmental processes. Overall, we showed that high levels of genetic diversity in Philippine Puddle Frogs predominantly comprise metapopulation lineages that arose through complex patterns of admixture, isolation-by-distance, and isolation-by-environment as opposed to species divergence. Our findings suggest that speciation may not be the major process underlying the high levels of hidden diversity observed in many taxonomic groups and that widely adopted tree- and distance-based methods overestimate species diversity in the presence of gene flow. [Cryptic species; gene flow; introgression; isolation-by-distance; isolation-by-environment; phylogenetic network; species delimitation.] 
    more » « less