Antagonistic pleiotropy (AP) is a genetic trade-off between different fitness components. In annual plants, a trade-off between days to flower (DTF) and reproductive capacity often determines how many individuals survive to flower in a short growing season, and also influences the seed set of survivors. We develop a model of viability and fecundity selection informed by many experiments on the yellow monkeyflower, Mimulus guttatus, but applicable to many annual species. A viability/fecundity trade-off maintains stable polymorphism under surprisingly general conditions. We also introduce both spatial heterogeneity and temporal stochasticity in environmental parameters. Neither is necessary for polymorphism, but spatial heterogeneity allows polymorphism while also generating the often observed non-negative correlations in fitness components.

The estimation of outcrossing rates in hermaphroditic species has been a major focus in the evolutionary study of reproductive strategies, and is also essential for plant breeding and conservation. Surprisingly, genomics has thus far minimally influenced outcrossing rate studies. In this article, we generalize a Bayesian inference method (BORICE) to accommodate genomic data from multiple subpopulations of a species. As an empirical demonstration, BORICE is applied to 115 maternal families of Mimulus guttatus. The analysis shows that low‐level whole genome sequencing of parents and offspring is sufficient for individualized mating system estimation: 208 offspring (88.5%) were definitively called as outcrossed, 23 (9.8%) as selfed. After mating system parameters are established (each offspring as outcrossed or selfed and the inbreeding level of maternal plants), BORICE outputs posterior genotype probabilities for each SNP genomewide. Individual SNP calls are often burdened with considerable uncertainty and distilling information from closely linked sites (within genomic windows) can be a useful strategy. For the Mimulus data, principal components based on window statistics were sufficient to diagnose inversion polymorphisms and estimate their effects on spatial structure, phenotypic and fitness measures. More generally, mating system estimation with BORICE can set the stage for population and quantitative genomic analyses, particularly researchers collect phenotypic or fitness data from maternal individuals.

Most flowering plants are hermaphroditic and experience strong pressures to evolve self‐pollination (automatic selection and reproductive assurance). Inbreeding depression (ID) can oppose selection for selfing, but it remains unclear if ID is typically strong enough to maintain outcrossing. To measure the full cost of sustained inbreeding on fitness, and its genomic basis, we planted highly homozygous, fully genome‐sequenced inbred lines of yellow monkeyflower (Mimulus guttatus) in the field next to outbred plants from crosses between the same lines. The cost of full homozygosity is severe: 65% for survival and 86% for lifetime seed production. Accounting for the unmeasured effect of lethal and sterile mutations, we estimate that the average fitness of fully inbred genotypes is only 3–4% that of outbred competitors. The genome sequence data provide no indication of simple overdominance, but the number of rare alleles carried by a line, especially within rare allele clusters nonrandomly distributed across the genome, is a significant negative predictor of fitness measurements. These findings are consistent with a deleterious allele model for ID. High variance in rare allele load among lines and the genomic distribution of rare alleles both suggest that migration might be an important source of deleterious alleles to local populations.

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the Allele Frequency Spectrum (AFS) of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with “rarity disequilibrium” produce a super-abundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of co-expressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high “connectedness” (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

Leopard complex spotting (LP), the result of an incompletely dominant mutation in TRPM1, produces a collection of unique depigmentation patterns in the horse. Although the LP mutation allows for expression of the various patterns, other loci are responsible for modification of the extent of white. Pedigree analysis of families segregating for high levels of patterning indicated a single dominant gene, named Pattern-1 (PATN1), as a major modifier of LP. Linkage analysis in two half-sibling families segregating for PATN1 identified a 15-Mb region on ECA3p that warranted further investigation. Whole transcriptome sequencing of skin samples from horses with and without the PATN1 allele was performed to identify genic SNPs for fine mapping. Two Sequenom assays were utilized to genotype 192 individuals from five LP-carrying breeds. The initial panel highlighted a 1.6-Mb region without a clear candidate gene. In the second round of fine mapping, SNP ECA3:23 658 447T>G in the 3'-UTR of RING finger and WD repeat domain 3 (RFWD3) reached a significance level of P = 1.063 × 10(-39). Sequencing of RFWD3 did not identify any coding polymorphisms specific to PATN1 horses. Genotyping of the RFWD3 3'-UTR SNP in 54 additional LP animals and 327 horses from nine breeds not segregating for LP further supported the association (P = 4.17 × 10(-115)). This variant is a strong candidate for PATN1 and may be particularly useful for LP breeders to select for high levels of white patterning.

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the Allele Frequency Spectrum (AFS) of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with rarity disequilibrium produce a super-abundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of co-expressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high connectedness (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

Summary Angiosperm diversity has been shaped by mating system evolution, with the most common transition from outcrossing to self-fertilizing.To investigate the genetic basis of this transition, we performed crosses between two species endemic to the Canary Islands, the self-compatible (SC) species Tolpis coronopifolia and its self-incompatible (SI) relative T. santosii. We scored self-compatibility as self-seedset of recombinant plants within two F2 populations.To map and genetically characterize the breakdown of SI, we built a draft genome sequence of T. coronopifolia, genotyped F2 plants using multiplexed shotgun genotyping (MSG), and located MSG markers to the genome sequence. We identified a single quantitative trait locus (QTL) that explains nearly all variation in self-seedset in both F2 populations.To identify putative causal genetic variants within the QTL, we performed transcriptome sequencing on mature floral tissue from both SI and SC species, constructed a transcriptome for each species, and then located each predicted transcript to the T. coronopifolia genome sequence. We annotated each predicted gene within the QTL and found two strong candidates for SI breakdown. Each gene has a coding sequence insertion/deletion mutation within the SC species that produces a truncated protein. Homologues of each gene have been implicated in pollen development, pollen germination, and pollen tube growth in other species.

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 (NASP SIM3 ), 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’.