Natural variation for a hybrid incompatibility between two species of Mimulus

Understanding the process by which hybrid incompatibility alleles become established in natural populations remains a major challenge to evolutionary biology. Previously, we discovered a two-locus Dobzhansky-Muller incompatibility that causes severe hybrid male sterility between two inbred lines of the incompletely isolated wildflower species, Mimulus guttatus and M. nasutus. An interspecific cross between these two inbred lines revealed that the M. guttatus (IM62) allele at hybrid male sterility 1 (hms1) acts dominantly in combination with recessive M. nasutus (SF5) alleles at hybrid male sterility 2 (hms2) to cause nearly complete hybrid male sterility. In this report, we extend these genetic analyses to investigate intraspecific variation for the hms1-hms2 incompatibility in natural populations of M. nasutus and M. guttatus, performing a series of interspecific crosses between individuals collected from a variety of geographic locales. Our results suggest that hms2 incompatibility alleles are common and geographically widespread within M. nasutus, but absent or rare in M. guttatus. In contrast, the hms1 locus is polymorphic within M. guttatus and the incompatibility allele appears to be extremely geographically restricted. We found evidence for the presence of the hms1 incompatibility allele in only two M. guttatus populations that exist within a few kilometers of each other. The restricted distribution of the hms1 incompatibility allele might currently limit the potential for the hms1-hms2 incompatibility to act as a species barrier between sympatric populations of M. guttatus and M. nasutus. Extensive sampling within a single M. guttatus population revealed that the hms1 locus is polymorphic and that the incompatibility allele appears to segregate at intermediate frequency, a pattern that is consistent with either genetic drift or natural selection.     

Patterns of nucleotide variation and reproductive isolation between a Mimulus allotetraploid and its progenitor species

Here we report our characterization of a widespread, highly selfing Mimulus allotetraploid formed by interspecific hybridization between M. nasutus and M. guttatus. Nucleotide variation at two nuclear loci (mCYCA and mAP3) within and among tetraploid populations resolves two haplotype clusters for each locus: one shares near identity with sequences from M. nasutus and the other group shares substantial variation with M. guttatus. With respect to the two loci studied, each allotetraploid individual is a 'fixed heterozygote' carrying sequences from both clusters. Moreover, mCYCA variation is consistent with at least two evolutionary origins for the Mimulus allotetraploid. We show that the allotetraploid is strongly reproductively isolated from M. nasutus and M. guttatus; interploidy crosses produce almost no viable seeds. By extension, we infer strong triploid block and argue that Mimulus allotetraploid formation might proceed in one step via the union of unreduced gametes in an M. nasutus-M. guttatus F(1) hybrid. We also discuss the potential roles of mating system and flowering asynchrony in allotetraploid establishment.     

Cytoplasmic male sterility in Mimulus hybrids has pleiotropic effects on corolla and pistil traits

The mechanisms underlying genetic associations have important consequences for evolutionary outcomes, but distinguishing linkage from pleiotropy is often difficult. Here, we use a fine mapping approach to determine the genetic basis of association between cytonuclear male sterility and other floral traits in Mimulus hybrids. Previous work has shown that male sterility in hybrids between Mimulus guttatus and Mimulus nasutus is due to interactions between a mitochondrial gene from M. guttatus and two tightly linked nuclear restorer alleles on Linkage Group 7, and that male sterility is associated with reduced corolla size. In the present study, we generated a set of nearly isogenic lines segregating for the restorer region and male sterility, but with unique flanking introgressions. Male-sterile flowers had significantly smaller corollas, longer styles and greater stigmatic exsertion than fertile flowers. Because these effects were significant regardless of the genotypic composition of introgressions flanking the restorer region, they suggest that these floral differences are a direct byproduct of the genetic incompatibility causing anther abortion. In addition, we found a non-significant but intriguing trend for male-sterile plants to produce more seeds per flower than fertile siblings after supplemental pollination. Such pleiotropic effects may underlie the corolla dimorphism frequently observed in gynodioecious taxa and may affect selection on cytoplasmic male sterility genes when they initially arise.     

Minor quantitative trait loci underlie floral traits associated with mating system divergence in Mimulus

The genetic basis of species differences provides insight into the mode and tempo of phenotypic divergence. We investigate the genetic basis of floral differences between two closely related plant taxa with highly divergent mating systems, Mimulus guttatus (large-flowered outcrosser) and M. nasutus (small-flowered selfer). We had previously constructed a framework genetic linkage map of the hybrid genome containing 174 markers spanning approximately 1800 cM on 14 linkage groups. In this study, we analyze the genetics of 16 floral, reproductive, and vegetative characters measured in a large segregating M. nasutus x M. guttatus F2 population (N = 526) and in replicates of the parental lines and F1 hybrids. Phenotypic analyses reveal strong genetic correlations among floral traits and epistatic breakdown of male and female fertility traits in the F2 hybrids. We use multitrait composite interval mapping to jointly locate and characterize quantitative trait loci (QTLs) underlying interspecific differences in seven floral traits. We identified 24 floral QTLs, most of which affected multiple traits. The large number of QTLs affecting each trait (mean = 13, range = 11-15) indicates a strikingly polygenic basis for floral divergence in this system. In general, QTL effects are small relative to both interspecific differences and environmental variation within genotypes, ruling out QTLs of major effect as contributors to floral divergence between M. guttatus and M. nasutus. QTLs show no pattern of directional dominance. Floral characters associated with pollinator attraction (corolla width) and self-pollen deposition (stigma-anther distance) share several pleiotropic or linked QTLs, but unshared QTLs may have allowed selfing to evolve independently from flower size. We discuss the polygenic nature of divergence between M. nasutus and M. guttatus in light of theoretical work on the evolution of selfing, genetics of adaptation, and maintenance of variation within populations.     

Pollen tube competition as a mechanism of prezygotic reproductive isolation between Mimulus nasutus and its presumed progenitor M. guttatus

This paper considers the extent to which differences in pollen tube growth rates can provide prezygotic reproductive isolation between Mimulus nasutus and its presumed progenitor, Mimulus guttatus . Mimulus nasutus is partially cleistogamous, but its larger chasmogamous flowers offer appreciable opportunity for outcrossing. Mimulus nasutus was found to have smaller pollen grains and shorter styles than M. guttatus . No differences were observed in pollen grain germination on conspecific and heterospecific stigmas. However, pollen tube growth rates of M. nasutus were found to be much slower than those of M. guttatus in the styles of that species. Consequently, any M. nasutus pollen transferred to an M. guttatus stigma was found to be competitively disadvantaged in an M. gutattus style. By contrast, no difference in pollen tube growth rate was detected between the species when growing in M. nasutus styles, possibly because M. nasutus styles are unable to support fast pollen tube growth. We tested the prediction from the pollen tube studies that a 50:50 mix of M. guttatus and M. nasutus pollen would produce 50% hybrid seeds when M. nasutus was the maternal parent, and near to 0% hybrid seed when M. guttatus was the maternal parent. The results were found to support this prediction. We conclude that pollen–pistil interactions can effect strong reproductive isolation between these species, as M. guttatus pollen tubes have a competitive advantage over those of M. nasutus in an M. guttatus style, but not in an M. nasutus style.     

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.     

A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions.

As part of a study of the genetics of floral adaptation and speciation in the Mimulus guttatus species complex, we constructed a genetic linkage map of an interspecific cross between M. guttatus and M. nasutus. We genotyped an F(2) mapping population (N = 526) at 255 AFLP, microsatellite, and gene-based markers and derived a framework map through repeated rounds of ordering and marker elimination. The final framework map consists of 174 marker loci on 14 linkage groups with a total map length of 1780 cM Kosambi. Genome length estimates (2011-2096 cM) indicate that this map provides thorough coverage of the hybrid genome, an important consideration for QTL mapping. Nearly half of the markers in the full data set (49%) and on the framework map (48%) exhibited significant transmission ratio distortion (alpha = 0.05). We localized a minimum of 11 transmission ratio distorting loci (TRDLs) throughout the genome, 9 of which generate an excess of M. guttatus alleles and a deficit of M. nasutus alleles. This pattern indicates that the transmission ratio distortion results from particular interactions between the heterospecific genomes and suggests that substantial genetic divergence has occurred between these Mimulus species. We discuss possible causes of the unequal representation of parental genomes in the F(2) generation.     

A cytonuclear incompatibility causes anther sterility in Mimulus hybrids

Multilocus interactions (also known as Dobzhansky-Muller incompatibilities) are thought to be the major source of hybrid inviability and sterility. Because cytoplasmic and nuclear genomes have conflicting evolutionary interests and are often highly coevolved, cytonuclear incompatibilities may be among the first to develop in incipient species. Here, we report the discovery of cytoplasm-dependent anther sterility in hybrids between closely related Mimulus species, outcrossing M. guttatus and selfing M. nasutus. A novel pollenless anther phenotype was observed in F2 hybrids with the M. guttatus cytoplasm (F2G) but not in the reciprocal F2N hybrids, F1 hybrids or parental genotypes. The pattern of phenotypic segregation in the F2G hybrids and two backcross populations fit a Mendelian single-locus recessive model, allowing us to map the underlying nuclear locus to a small region on LG7 of the Mimulus linkage map. Anther sterility was associated with a 20% reduction in flower size in backcross hybrids and we mapped a major cytoplasm-dependent corolla width QTL with its peak at the anther sterility locus. We argue that the cytonuclear anther sterility seen in hybrids reflects the presence of a cryptic cytoplasmic male sterility (CMS) and restorer system within the hermaphroditic M. guttatus population and therefore name the anther sterility locus restorer-of-male-fertility (RMF). The genetic mapping of RMF is a first step toward testing hypotheses about the molecular basis, individual fitness consequences, and ecological context of CMS and restoration in a system without stable CMS-restorer polymorphism (i.e., gynodioecy). The discovery of cryptic CMS in a hermaphroditic wildflower further suggests that selfish cytoplasmic evolution may play an important, but often undetected, role in shaping patterns of hybrid incompatibility and interspecific introgression in plants.     

Toward the evolutionary genomics of gametophytic divergence: patterns of transmission ratio distortion in monkeyflower (Mimulus) hybrids reveal a complex genetic basis for conspecific pollen precedence

Conspecific pollen precedence (CPP) is a major component of reproductive isolation between many flowering plant taxa and may reveal mechanisms of gametophytic evolution within species, but little is known about the genetic basis and evolutionary history of CPP. We systematically investigated the genetic architecture of CPP using patterns of transmission ratio distortion (TRD) in F2 and backcross hybrids between closely related species of Mimulus (Phrymaceae) with divergent mating systems. We found that CPP in Mimulus hybrids was polygenic and was the majority source of interspecific TRD genome-wide, with at least eight genomic regions contributing to the transmission advantage of M. guttatus pollen grains on M. guttatus styles. In aggregate, these male-specific transmission ratio distorting loci (TRDLs) were more than sufficient to account for the 100% precedence of pure M. guttatus pollen over M. nasutus pollen in mixed pollinations of M. guttatus. All but one of these pollen TRDLs were style-dependent; that is, we observed pollen TRD in F(1) and/or M. guttatus styles, but not in M. nasutus styles. These findings suggest that species-specific differences in pollen tube performance accumulate gradually and may have been driven by coevolution between pollen and style in the predominantly outcrossing M. guttatus.     

Patterns of nucleotide diversity in two species of Mimulus are affected by mating system and asymmetric introgression

The evolutionary transition from outcrossing to self-fertilization has far-reaching implications for patterns of intraspecific genetic diversity and the potential for speciation. Using DNA sequence variation at two nuclear loci, we examined the divergence history of two closely related species of Mimulus. To investigate the effects of mating system and introgressive hybridization on the outcrossing M. guttatus and the selfing M. nasutus, we inspected nucleotide diversity within and between natural populations spanning the species' geographic ranges. High sequence similarity among populations of the selfing M. nasutus points to a single evolutionary origin for the species. Consistent with their distinct mating systems, all genetic variation in M. nasutus is distributed among populations, whereas M. guttatus exhibits appreciable levels of nucleotide diversity within populations. Silent genetic diversity is extensive in M. guttatus (mean theta(sil)/site = 0.077) and greatly exceeds the predicted twofold elevation in neutral variation for outcrossers relative to selfers. The finding of several M. guttatus sequences that share complete identity with sequences from M. nasutus suggests that recent asymmetric introgression may have occurred. We argue that exceptionally high nucleotide diversity in M. guttatus is consistent with a long-term history of directional introgression from M. nasutus to M. guttatus throughout the divergence of these two species.     

Ecological divergence associated with mating system causes nearly complete reproductive isolation between sympatric Mimulus species

Speciation often involves the evolution of numerous prezygotic and postzygotic isolating barriers between divergent populations. Detailed knowledge of the strength and nature of those barriers provides insight into ecological and genetic factors that directly or indirectly influenced their origin, and may help predict whether they will be maintained in the face of sympatric hybridization and introgression. We estimated the magnitude of pre- and postzygotic barriers between naturally occurring sympatric populations of Mimulus guttatus and M. nasutus. Prezygotic barriers, including divergent flowering phenologies, differential pollen production, mating system isolation, and conspecific pollen precedence, act asymmetrically to completely prevent the formation of F(1) hybrids among seeds produced by M. guttatus (F(1)g), and reduce F(1) hybrid production among seeds produced by M. nasutus (F(1)n) to only about 1%. Postzygotic isolation is also asymmetric: in field experiments, F(1)g but not F(1)n hybrids had significantly reduced germination rates and survivorship compared to parental species. Both hybrid classes had flower, pollen, and seed production values within the range of parental values. Despite the moderate degree of F(1)g hybrid inviability, postzygotic isolation contributes very little to the total isolation between these species in the wild. We also found that F(1) hybrid flowering phenology overlapped more with M. guttatus than M. nasutus. These results, taken together, suggest greater potential for introgression from M. nasutus to M. guttatus than for the reverse direction. We also address problems with commonly used indices of isolation, discuss difficulties in calculating meaningful measures of reproductive isolation when barriers are asymmetric, and propose novel measures of prezygotic isolation that are consistent with postzygotic measures.     

A novel meiotic drive locus almost completely distorts segregation in mimulus (monkeyflower) hybrids

We report the discovery, mapping, and characterization of a meiotic drive locus (D) exhibiting nearly 100% nonrandom transmission in hybrids between two species of yellow monkeyflowers, outcrossing Mimulus guttatus and selfing M. nasutus. Only 1% of F(2) hybrids were M. nasutus homozygotes at the marker most tightly linked to D. We used a set of reciprocal backcrosses to distinguish among male-specific, female-specific, and zygote-specific sources of transmission ratio distortion. Transmission was severely distorted only when the heterozygous F(1) acted as the female parent in crosses to either parental species, ruling out pollen competition and zygote mortality as potential sources of drive. After four generations of backcrossing to M. nasutus, nearly isogenic lines were still >90% heterozygous at markers linked to D, suggesting that heterozygosity at the drive locus alone is sufficient for nonrandom transmission. A lack of dramatic female fitness costs in these lines rules out alternatives involving ovule or seed mortality and points to a truly meiotic mechanism of drive. The strength and direction of drive in this system is consistent with population genetic theory of selfish element evolution under different mating systems. These results are the first empirical demonstration of the strong female-specific drive predicted by new models of selfish centromere turnover.     

Transmission ratio distortion in intraspecific hybrids of Mimulus guttatus: implications for genomic divergence

We constructed a genetic linkage map between two divergent populations of Mimulus guttatus. We genotyped an F(2) mapping population (N = 539) at 154 AFLP, microsatellite, and gene-based markers. A framework map was constructed consisting of 112 marker loci on 14 linkage groups with a total map length of 1518 cM Kosambi. Nearly half of all markers (48%) exhibited significant transmission ratio distortion (alpha = 0.05). By using a Bayesian multipoint mapping method and visual inspection of significantly distorted markers, we detected 12 transmission ratio distorting loci (TRDL) throughout the genome. The high degree of segregation distortion detected in this intraspecific map indicates substantial genomic divergence that perhaps suggests genomic incompatibilities between these two populations. We compare the pattern of transmission ratio distortion in this map to an interspecific map constructed between M. guttatus and M. nasutus. A similar level of segregation distortion is detected in both maps. Collinear regions between maps are compared to determine if there are shared genetic patterns of non-Mendelian segregation distortion within and among Mimulus species.     

Polymorphic microsatellite loci in Mimulus guttatus and related species

The evolution of self-fertilization from outcrossing is one of the most common evolutionary transitions in flowering plants. One goal of our research is to develop a linkage map based upon microsatellite loci in Mimulus guttatus and the selfing derivative M. nasutus in order to examine the genetics of one such transition. Here we present 19 primer pairs that detect polymorphisms between these two species. Within-population allelic surveys demonstrate the value of these markers for detailed analysis of population structure. We have found also that many of the microsatellites are conserved in more distantly related Mimulus species, and may thus serve as useful genetic markers for population genetic and comparative mapping studies within this diverse genus.     

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 F₁ 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 F₁ 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 F₁s showed highly restored fertility (83.4% pollen fertility) relative to diploids F₁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.     

Pleiotropic quantitative trait loci contribute to population divergence in traits associated with life-history variation in Mimulus guttatus

Evolutionary biologists seek to understand the genetic basis for multivariate phenotypic divergence. We constructed an F2 mapping population (N = 539) between two distinct populations of Mimulus guttatus. We measured 20 floral, vegetative, and life-history characters on parents and F1 and F2 hybrids in a common garden experiment. We employed multitrait composite interval mapping to determine the number, effect, and degree of pleiotropy in quantitative trait loci (QTL) affecting divergence in floral, vegetative, and life-history characters. We detected 16 QTL affecting floral traits; 7 affecting vegetative traits; and 5 affecting selected floral, vegetative, and life-history traits. Floral and vegetative traits are clearly polygenic. We detected a few major QTL, with all remaining QTL of small effect. Most detected QTL are pleiotropic, implying that the evolutionary shift between these annual and perennial populations is constrained. We also compared the genetic architecture controlling floral trait divergence both within (our intraspecific study) and between species, on the basis of a previously published analysis of M. guttatus and M. nasutus. Eleven of our 16 floral QTL map to approximately the same location in the interspecific map based on shared, collinear markers, implying that there may be a shared genetic basis for floral divergence within and among species of Mimulus.     

Mapping and characterization of a 'speciation gene' in Drosophila.

Almost nothing is known about the identity of the genes causing reproductive isolation between species. As a first step towards molecular isolation of a 'speciation gene', I mapped and partly characterized a gene causing hybrid male sterility in Drosophila. This analysis shows that sterility of D. melanogaster males who carry the 'dot' fourth chromosome from D. simulans is due entirely to a very small region of the D. simulans chromosome (including only about 5 salivary gland bands or approximately 250 kb of DNA). Thus the hybrid sterility effect of the D. simulans fourth chromosome is almost surely due to a single gene of very large effect (here named hms, hybrid male sterile). Hms is zygotically acting, and the D. simulans allele of hms is completely recessive. Furthermore, complementation tests suggest that hms is not an allele of any known locus in D. melanogaster.     

Epistasis among Drosophila persimilis factors conferring hybrid male sterility with D. pseudoobscura bogotana

The Bateson-Dobzhansky-Muller model posits that hybrid incompatibilities result from genetic changes that accumulate during population divergence. Indeed, much effort in recent years has been devoted to identifying genes associated with hybrid incompatibilities, often with limited success, suggesting that hybrid sterility and inviability are frequently caused by complex interactions between multiple loci and not by single or a small number of gene pairs. Our previous study showed that the nature of epistasis between sterility-conferring QTL in the Drosophila persimilis-D. pseudoobscura bogotana species pair is highly specific. Here, we further dissect one of the three QTL underlying hybrid male sterility between these species and provide evidence for multiple factors within this QTL. This result indicates that the number of loci thought to contribute to hybrid dysfunction may have been underestimated, and we discuss how linkage and complex epistasis may be characteristic of the genetics of hybrid incompatibilities. We further pinpoint the location of one locus that confers hybrid male sterility when homozygous, dubbed "mule-like", to roughly 250 kilobases.     

Five anthocyanin polymorphisms are associated with an R2R3-MYB cluster in Mimulus guttatus (Phrymaceae)

? Premise of study: Botanists have long been interested in the reasons for genetic variation among individuals, populations, and species of plants. The anthocyanin pathway is ideal for studying the evolution of such phenotypic variation. ? Methods: We used a combination of quantitative trait loci mapping and association studies to understand the genetic basis of variation in five anthocyanin phenotypes including calyx, corolla, and leaf coloration patterns that vary within and among populations of Mimulus guttatus. We then examined what genes might be responsible for this phenotypic variation and whether one of the traits, calyx spotting, is randomly distributed across the geographic range of the species. ? Key results: All five phenotypes in M. guttatus were primarily controlled by the same major locus (PLA1), which contains a tandem array of three R2R3-MYB genes known to be involved in the evolution of flower color in a related species of Mimulus. Calyx spotting was nonrandomly distributed across the range of M. guttatus and correlated with multiple climate variables. ? Conclusions: The results of this study suggest that variation in R2R3-MYB genes is the primary cause of potentially important anthocyanin phenotypic variation within and among populations of M. guttatus, a finding consistent with recent theoretical and empirical research on flower color evolution.