Genomic mechanisms and consequences of diverse postzygotic barriers between monkeyflower species.

The evolution of genomic incompatibilities causing postzygotic barriers to hybridization is a key step in species divergence. Incompatibilities take 2 general forms-structural divergence between chromosomes leading to severe hybrid sterility in F1 hybrids and epistatic interactions between genes causing reduced fitness of hybrid gametes or zygotes (Dobzhansky-Muller incompatibilities). Despite substantial recent progress in understanding the molecular mechanisms and evolutionary origins of both types of incompatibility, how each behaves across multiple generations of hybridization remains relatively unexplored. Here, we use genetic mapping in F2 and recombinant inbred line (RIL) hybrid populations between the phenotypically divergent but naturally hybridizing monkeyflowers Mimulus cardinalis and M. parishii to characterize the genetic basis of hybrid incompatibility and examine its changing effects over multiple generations of experimental hybridization. In F2s, we found severe hybrid pollen inviability (<50% reduction vs parental genotypes) and pseudolinkage caused by a reciprocal translocation between Chromosomes 6 and 7 in the parental species. RILs retained excess heterozygosity around the translocation breakpoints, which caused substantial pollen inviability when interstitial crossovers had not created compatible heterokaryotypic configurations. Strong transmission ratio distortion and interchromosomal linkage disequilibrium in both F2s and RILs identified a novel 2-locus genic incompatibility causing sex-independent gametophytic (haploid) lethality. The latter interaction eliminated 3 of the expected 9 F2 genotypic classes via F1 gamete loss without detectable effects on the pollen number or viability of F2 double heterozygotes. Along with the mapping of numerous milder incompatibilities, these key findings illuminate the complex genetics of plant hybrid breakdown and are an important step toward understanding the genomic consequences of natural hybridization in this model system.

Supergene potential of a selfish centromere.

Selfishly evolving centromeres bias their transmission by exploiting the asymmetry of female meiosis and preferentially segregating to the egg. Such female meiotic drive systems have the potential to be supergenes, with multiple linked loci contributing to drive costs or enhancement. Here, we explore the supergene potential of a selfish centromere (D) in Mimulus guttatus, which was discovered in the Iron Mountain (IM) Oregon population. In the nearby Cone Peak population, D is still a large, non-recombining and costly haplotype that recently swept, but shorter haplotypes and mutational variation suggest a distinct population history. We detected D in five additional populations spanning more than 200 km; together, these findings suggest that selfish centromere dynamics are widespread in M. guttatus. Transcriptome comparisons reveal elevated differences in expression between driving and non-driving haplotypes within, but not outside, the drive region, suggesting large-scale cis effects of D's spread on gene expression. We use the expression data to refine linked candidates that may interact with drive, including Nuclear Autoantigenic Sperm Protein (NASPSIM3), which chaperones the centromere-defining histone CenH3 known to modify Mimulus drive. Together, our results show that selfishly evolving centromeres may exhibit supergene behaviour and lay the foundation for future genetic dissection of drive and its costs. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Population genomic consequences of life-history and mating system adaptation to a geothermal soil mosaic in yellow monkeyflowers.

Local selection can promote phenotypic divergence despite gene flow across habitat mosaics, but adaptation itself may generate substantial barriers to genetic exchange. In plants, life-history, phenology, and mating system divergence have been proposed to promote genetic differentiation in sympatry. In this study, we investigate phenotypic and genetic variation in Mimulus guttatus (yellow monkeyflowers) across a geothermal soil mosaic in Yellowstone National Park (YNP). Plants from thermal annual and nonthermal perennial habitats were heritably differentiated for life-history and mating system traits, consistent with local adaptation to the ephemeral thermal-soil growing season. However, genome-wide genetic variation primarily clustered plants by geographic region, with little variation sorting by habitat. The one exception was an extreme thermal population also isolated by a 200 m geographical gap of no intermediate habitat. Individual inbreeding coefficients (FIS ) were higher (and predicted by trait variation) in annual plants and annual pairs showed greater isolation by distance at local (<1 km) scales. Finally, YNP adaptation does not reuse a widespread inversion that underlies M. guttatus life-history ecotypes range-wide, suggesting a novel genetic mechanism. Overall, this work suggests that life-history and mating system adaptation strong enough to shape individual mating patterns does not necessarily generate incipient speciation without geographical barriers.

Quantitative trait locus mapping reveals an independent genetic basis for joint divergence in leaf function, life-history, and floral traits between scarlet monkeyflower (Mimulus cardinalis) populations.

PREMISE: Across taxa, vegetative and floral traits that vary along a fast-slow life-history axis are often correlated with leaf functional traits arrayed along the leaf economics spectrum, suggesting a constrained set of adaptive trait combinations. Such broad-scale convergence may arise from genetic constraints imposed by pleiotropy (or tight linkage) within species, or from natural selection alone. Understanding the genetic basis of trait syndromes and their components is key to distinguishing these alternatives and predicting evolution in novel environments. METHODS: We used a line-cross approach and quantitative trait locus (QTL) mapping to characterize the genetic basis of twenty leaf functional/physiological, life history, and floral traits in hybrids between annualized and perennial populations of scarlet monkeyflower (Mimulus cardinalis). RESULTS: We mapped both single and multi-trait QTLs for life history, leaf function and reproductive traits, but found no evidence of genetic co-ordination across categories. A major QTL for three leaf functional traits (thickness, photosynthetic rate, and stomatal resistance) suggests that a simple shift in leaf anatomy may be key to adaptation to seasonally dry habitats. CONCLUSIONS: Our results suggest that the co-ordination of resource-acquisitive leaf physiological traits with a fast life-history and more selfing mating system results from environmental selection rather than functional or genetic constraint. Independent assortment of distinct trait modules, as well as a simple genetic basis to leaf physiological traits associated with drought escape, may facilitate adaptation to changing climates.

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.

Ancient and recent introgression shape the evolutionary history of pollinator adaptation and speciation in a model monkeyflower radiation (Mimulus section Erythranthe).

Inferences about past processes of adaptation and speciation require a gene-scale and genome-wide understanding of the evolutionary history of diverging taxa. In this study, we use genome-wide capture of nuclear gene sequences, plus skimming of organellar sequences, to investigate the phylogenomics of monkeyflowers in Mimulus section Erythranthe (27 accessions from seven species). Taxa within Erythranthe, particularly the parapatric and putatively sister species M. lewisii (bee-pollinated) and M. cardinalis (hummingbird-pollinated), have been a model system for investigating the ecological genetics of speciation and adaptation for over five decades. Across >8000 nuclear loci, multiple methods resolve a predominant species tree in which M. cardinalis groups with other hummingbird-pollinated taxa (37% of gene trees), rather than being sister to M. lewisii (32% of gene trees). We independently corroborate a single evolution of hummingbird pollination syndrome in Erythranthe by demonstrating functional redundancy in genetic complementation tests of floral traits in hybrids; together, these analyses overturn a textbook case of pollination-syndrome convergence. Strong asymmetries in allele sharing (Patterson's D-statistic and related tests) indicate that gene tree discordance reflects ancient and recent introgression rather than incomplete lineage sorting. Consistent with abundant introgression blurring the history of divergence, low-recombination and adaptation-associated regions support the new species tree, while high-recombination regions generate phylogenetic evidence for sister status for M. lewisii and M. cardinalis. Population-level sampling of core taxa also revealed two instances of chloroplast capture, with Sierran M. lewisii and Southern Californian M. parishii each carrying organelle genomes nested within respective sympatric M. cardinalis clades. A recent organellar transfer from M. cardinalis, an outcrosser where selfish cytonuclear dynamics are more likely, may account for the unexpected cytoplasmic male sterility effects of selfer M. parishii organelles in hybrids with M. lewisii. Overall, our phylogenomic results reveal extensive reticulation throughout the evolutionary history of a classic monkeyflower radiation, suggesting that natural selection (re-)assembles and maintains species-diagnostic traits and barriers in the face of gene flow. Our findings further underline the challenges, even in reproductively isolated species, in distinguishing re-use of adaptive alleles from true convergence and emphasize the value of a phylogenomic framework for reconstructing the evolutionary genetics of adaptation and speciation.

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

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

Standard Deviations: The Biological Bases of Transmission Ratio Distortion.

The rule of Mendelian inheritance is remarkably robust, but deviations from the equal transmission of alternative alleles at a locus [a.k.a. transmission ratio distortion (TRD)] are also commonly observed in genetic mapping populations. Such TRD reveals locus-specific selection acting at some point between the diploid heterozygous parents and progeny genotyping and therefore can provide novel insight into otherwise-hidden genetic and evolutionary processes. Most of the classic selfish genetic elements were discovered through their biasing of transmission, but many unselfish evolutionary and developmental processes can also generate TRD. In this review, we describe methodologies for detecting TRD in mapping populations, detail the arenas and genetic interactions that shape TRD during plant and animal reproduction, and summarize patterns of TRD from across the genetic mapping literature. Finally, we point to new experimental approaches that can accelerate both detection of TRD and characterization of the underlying genetic mechanisms.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity.

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

Making a Murderer: The Evolutionary Framing of Hybrid Gamete-Killers.

Recent molecular investigations of hybrid incompatibilities have revealed fascinating patterns of genetic interactions that have been interpreted as the remnants of a history of selfish evolution. Instead of framing hybrid incompatibilities in light of genetic conflict, we advocate assuming their innocence. Researchers must build a strong theory for each case, supported by population genetic evidence, such that the role of conflict in the evolution of a hybrid incompatibility can be proven beyond reasonable doubt. This will require careful investigation of the evolutionary history of these incompatibilities, a reckoning of how the reproductive biology of study organisms impacts on the likelihood of genetic conflict, and molecular evidence of the rapid selfish spread of these alleles.

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers.

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness-relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate-frequency three-copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over-represented in late-flowering and large-flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250-300 copies of RLG1a totalling ~5.7 Mb (20-40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR-based copy counts indicate that high+ is a single tandem array unlinked to the single-copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single-copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress-responsive life-history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.

When Two Rights Make a Wrong: The Evolutionary Genetics of Plant Hybrid Incompatibilities.

Hybrids between flowering plant species often exhibit reduced fitness, including sterility and inviability. Such hybrid incompatibilities create barriers to genetic exchange that can promote reproductive isolation between diverging populations and, ultimately, speciation. Additionally, hybrid breakdown opens a window into hidden molecular and evolutionary processes occurring within species. Here, we review recent work on the mechanisms and origins of hybrid incompatibility in flowering plants, including both diverse genic interactions and chromosomal incompatibilities. Conflict and coevolution among and within plant genomes contributes to the evolution of some well-characterized genic incompatibilities, but duplication and drift also play important roles. Inversions, while contributing to speciation by suppressing recombination, rarely cause underdominant sterility. Translocations cause severe F1 sterility by disrupting meiosis in heterozygotes, making their fixation in outcrossing sister species a paradox. Evolutionary genomic analyses of both genic and chromosomal incompatibilities, in the context of population genetic theory, can explicitly test alternative scenarios for their origins.

Selfish evolution of cytonuclear hybrid incompatibility in Mimulus.

Intraspecific coevolution between selfish elements and suppressors may promote interspecific hybrid incompatibility, but evidence of this process is rare. Here, we use genomic data to test alternative models for the evolution of cytonuclear hybrid male sterility in Mimulus In hybrids between Iron Mountain (IM) Mimulus guttatus × Mimulus nasutus, two tightly linked M. guttatus alleles (Rf1/Rf2) each restore male fertility by suppressing a local mitochondrial male-sterility gene (IM-CMS). Unlike neutral models for the evolution of hybrid incompatibility loci, selfish evolution predicts that the Rf alleles experienced strong selection in the presence of IM-CMS. Using whole-genome sequences, we compared patterns of population-genetic variation in Rf at IM to a neighbouring population that lacks IM-CMS. Consistent with local selection in the presence of IM-CMS, the Rf region shows elevated FST, high local linkage disequilibrium and a distinct haplotype structure at IM, but not at Cone Peak (CP), suggesting a recent sweep in the presence of IM-CMS. In both populations, Rf2 exhibited lower polymorphism than other regions, but the low-diversity outliers were different between CP and IM. Our results confirm theoretical predictions of ubiquitous cytonuclear conflict in plants and provide a population-genetic mechanism for the evolution of a common form of hybrid incompatibility.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus.

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 with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription 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 of Mimulus guttatus adaptation 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.

Genetic loci with parent-of-origin effects cause hybrid seed lethality in crosses between Mimulus species.

In flowering plants, F1 hybrid seed lethality is a common outcome of crosses between closely related diploid species, but the genetic basis of this early-acting and potentially widespread form of postzygotic reproductive isolation is largely unknown. We intercrossed two closely related species of monkeyflower, Mimulus guttatus and Mimulus tilingii, to characterize the mechanisms and strength of postzygotic reproductive isolation. Then, using a reciprocal backcross design, we performed high-resolution genetic mapping to determine the genetic architecture of hybrid seed lethality and directly test for loci with parent-of-origin effects. We found that F1 hybrid seed lethality is an exceptionally strong isolating barrier between Mimulus species, with reciprocal crosses producing < 1% viable seeds. This form of postzygotic reproductive isolation appears to be highly polygenic, indicating that multiple incompatibility loci have accumulated rapidly between these closely related Mimulus species. It is also primarily caused by genetic loci with parent-of-origin effects, suggesting a possible role for imprinted genes in the evolution of Mimulus hybrid seed lethality. Our findings suggest that divergence in loci with parent-of-origin effects, which is probably driven by genomic coevolution within lineages, might be an important source of hybrid incompatibilities between flowering plant species.

The Ecology and Evolutionary Dynamics of Meiotic Drive.

Meiotic drivers are genetic variants that selfishly manipulate the production of gametes to increase their own rate of transmission, often to the detriment of the rest of the genome and the individual that carries them. This genomic conflict potentially occurs whenever a diploid organism produces a haploid stage, and can have profound evolutionary impacts on gametogenesis, fertility, individual behaviour, mating system, population survival, and reproductive isolation. Multiple research teams are developing artificial drive systems for pest control, utilising the transmission advantage of drive to alter or exterminate target species. Here, we review current knowledge of how natural drive systems function, how drivers spread through natural populations, and the factors that limit their invasion.

A Segregating Inversion Generates Fitness Variation in Yellow Monkeyflower (Mimulus guttatus).

Polymorphic chromosomal rearrangements can bind hundreds of genes into single genetic loci with diverse effects. Rearrangements are often associated with local adaptation and speciation and may also be an important component of genetic variation within populations. We genetically and phenotypically characterize a segregating inversion (inv6) in the Iron Mountain (IM) population of Mimulus guttatus (yellow monkeyflower). We initially mapped inv6 as a region of recombination suppression in three F2 populations resulting from crosses among IM plants. In each case, the F1 parent was heterozygous for a derived haplotype, homogenous across markers spanning over 5 Mb of chromsome 6. In the three F2 populations, inv6 reduced male and female fitness components. In addition,i nv6 carriers suffered an ∼30% loss of pollen viability in the field. Despite these costs, inv6 exists at moderate frequency (∼8%) in the natural population, suggesting counterbalancing fitness benefits that maintain the polymorphism. Across 4 years of monitoring in the field, inv6 had an overall significant positive effect on seed production (lifetime female fitness) of carriers. This benefit was particularly strong in harsh years and may be mediated (in part) by strong positive effects on flower production. These data suggest that opposing fitness effects maintain an intermediate frequency, and as a consequence, inv6 generates inbreeding depression and high genetic variance. We discuss these findings in relation to the theory of inbreeding depression and the maintenance of fitness variation.

Duplication and Adaptive Evolution of a Key Centromeric Protein in Mimulus, a Genus with Female Meiotic Drive.

The fundamental asymmetry of female meiosis creates an arena for genetic elements to compete for inclusion in the egg, promoting the selfish evolution of centromere variants that maximize their transmission to the future egg. Such "female meiotic drive" has been hypothesized to explain the paradoxically complex and rapidly evolving nature of centromeric DNA and proteins. Although theoretically widespread, few cases of active drive have been observed, thereby limiting the opportunities to directly assess the impact of centromeric drive on molecular variation at centromeres and binding proteins. Here, we characterize the molecular evolutionary patterns of CENH3, the centromere-defining histone variant, in Mimulus monkeyflowers, a genus with one of the few known cases of active centromere-associated female meiotic drive. First, we identify a novel duplication of CENH3 in diploid Mimulus, including in lineages with actively driving centromeres. Second, we demonstrate long-term adaptive evolution at several sites in the N-terminus of CENH3, a region with some meiosis-specific functions that putatively interacts with centromeric DNA. Finally, we infer that the paralogs evolve under different selective regimes; some sites in the N-terminus evolve under positive selection in the pro-orthologs or only one paralog (CENH3_B) and the paralogs exhibit significantly different patterns of polymorphism within populations. Our finding of long-term, adaptive evolution at CENH3 in the context of centromere-associated meiotic drive supports an antagonistic, coevolutionary battle for evolutionary dominance between centromeric DNA and binding proteins.

Centromere-associated meiotic drive and female fitness variation in Mimulus.

Female meiotic drive, in which chromosomal variants preferentially segregate to the egg pole during asymmetric female meiosis, is a theoretically pervasive but still mysterious form of selfish evolution. Like other selfish genetic elements, driving chromosomes may be maintained as balanced polymorphisms by pleiotropic or linked fitness costs. A centromere-associated driver (D) with a ∼58:42 female-specific transmission advantage occurs at intermediate frequency (32-40%) in the Iron Mountain population of the yellow monkeyflower, Mimulus guttatus. Previously determined male fertility costs are sufficient to prevent the fixation of D, but predict a higher equilibrium frequency. To better understand the dynamics and effects of D, we developed a new population genetic model and measured genotype-specific lifetime female fitness in the wild. In three of four years, and across all years, D imposed significant recessive seedset costs, most likely due to hitchhiking by deleterious mutations. With both male and female costs as measured, and 58:42 drive, our model predicts an equilibrium frequency of D (38%) very close to the observed value. Thus, D represents a rare selfish genetic element whose local population genetic dynamics have been fully parameterized, and the observation of equilibrium sets the stage for investigations of coevolution with suppressors.

The genetic architecture of traits associated with the evolution of self-pollination in Mimulus.

Quantitative trait locus (QTL) mapping is a first step toward understanding the genetic basis of adaptive evolution and may also reveal reproductive incompatibilities unique to hybrids. In plants, the shift from outcrossing to self-pollination is common, providing the opportunity for comparisons of QTL architecture among parallel evolutionary transitions. We used QTL mapping in hybrids between the bee-pollinated monkeyflower Mimulus lewisii and the closely related selfer Mimulus parishii to determine the genetic basis of divergence in floral traits and flowering time associated with mating-system evolution, and to characterize hybrid anther sterility. We found a moderately polygenic and highly directional basis for floral size evolution, suggesting adaptation from standing variation or in pursuit of a moving optimum, whereas only a few major loci accounted for substantial flowering-time divergence. Cytonuclear incompatibilities caused hybrid anther sterility, confounding estimation of reproductive organ QTLs. The genetic architecture of floral traits associated with selfing in M. parishii was primarily polygenic, as in other QTL studies of this transition, but in contrast to the previously characterized oligogenic basis of a pollinator shift in close relatives. Hybrid anther sterility appeared parallel at the molecular level to previously characterized incompatibilities, but also raised new questions about cytonuclear co-evolution in plants.