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.

Evolution of the Highly Repetitive PEVK Region of Titin Across Mammals.

The protein titin plays a key role in vertebrate muscle where it acts like a giant molecular spring. Despite its importance and conservation over vertebrate evolution, a lack of high quality annotations in non-model species makes comparative evolutionary studies of titin challenging. The PEVK region of titin-named for its high proportion of Pro-Glu-Val-Lys amino acids-is particularly difficult to annotate due to its abundance of alternatively spliced isoforms and short, highly repetitive exons. To understand PEVK evolution across mammals, we developed a bioinformatics tool, PEVK_Finder, to annotate PEVK exons from genomic sequences of titin and applied it to a diverse set of mammals. PEVK_Finder consistently outperforms standard annotation tools across a broad range of conditions and improves annotations of the PEVK region in non-model mammalian species. We find that the PEVK region can be divided into two subregions (PEVK-N, PEVK-C) with distinct patterns of evolutionary constraint and divergence. The bipartite nature of the PEVK region has implications for titin diversification. In the PEVK-N region, certain exons are conserved and may be essential, but natural selection also acts on particular codons. In the PEVK-C, exons are more homogenous and length variation of the PEVK region may provide the raw material for evolutionary adaptation in titin function. The PEVK-C region can be further divided into a highly repetitive region (PEVK-CA) and one that is more variable (PEVK-CB). Taken together, we find that the very complexity that makes titin a challenge for annotation tools may also promote evolutionary adaptation.

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.

Genes Integral to the Reproductive Function of Male Reproductive Tissues Drive Heterogeneity in Evolutionary Rates in Japanese Quail.

Early comparative genomics studies originally uncovered a nonintuitive pattern; genes involved in reproduction appeared to evolve more rapidly than other classes of genes. Currently, the emerging consensus is that genes encoding reproductive proteins evolve under variable selective pressures, producing more heterogeneous divergence patterns than previously appreciated. Here, we investigate a facet of that heterogeneity and explore the factors that drive male reproductive tissue-based heterogeneity in evolutionary rates. In Japanese quail (Coturnix japonica), genes with enriched expression in the testes evolve much more rapidly than those enriched in the foam gland (FG), a novel gland that secretes an airy foam that males transfer to females during mating. We compared molecular evolutionary patterns among (1) genes with induced expression in breeding vs. wintering conditions for both tissues and (2) genes that encode foam proteins (FPs) vs. those with varying degrees of expression specificity in the FG. We report two major findings. First, genes upregulated in breeding condition testes evolve exceptionally rapidly, while those induced in breeding condition FGs evolve slowly. These differences hold even after correcting for hormonally-dependent gene expression and chromosomal location. Second, genes encoding FPs are extremely conserved in terms of gene identity and sequence. Together, these finding suggest that genes involved in the reproductive function of each tissue drive the marked rate of heterogeneity.

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.

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.

A comparison of next generation sequencing technologies for transcriptome assembly and utility for RNA-Seq in a non-model bird.

De novo assembled transcriptomes, in combination with RNA-Seq, are powerful tools to explore gene sequence and expression level in organisms without reference genomes. Investigators must first choose which high throughput sequencing platforms will provide data most suitable for their experimental goals. In this study, we explore the utility of 454 and Illumina sequences for de novo transcriptome assembly and downstream RNA-Seq applications in a reproductive gland from a non-model bird species, the Japanese quail (Coturnix japonica). Four transcriptomes composed of either pure 454 or Illumina reads or mixtures of read types were assembled and evaluated for the same cost. Illumina assemblies performed best for de novo transcriptome characterization in terms of contig length, transcriptome coverage, and complete assembly of gene transcripts. Improvements over the Hybrid assembly were marginal, with the exception that the addition of 454 data significantly increased the number of genes annotated. The Illumina assembly provided the best reference to align an independent set of RNA-Seq data as ∼84% of reads mapped to single genes in the transcriptome. Contigs constructed solely from 454 data may impose problems for RNA-Seq as our 454 transcriptome revealed a high number of indels and many ambiguously mapped reads. Correcting the 454 transcriptome with Illumina reads was an effective strategy to deal with indel and frameshift errors inherent to the 454 transcriptome, but at the cost of transcriptome coverage. In the absence of a reference genome, we find that Illumina reads alone produced a high quality transcriptome appropriate for RNA-Seq gene expression analyses.

Selective constraint dominates the evolution of genes expressed in a novel reproductive gland.

One striking pattern in molecular evolution is that genes encoding proteins involved in reproduction tend to evolve rapidly. Seminal fluid proteins frequently exhibit this pattern and directly affect multiple reproductive processes including enhancing sperm performance and mediating postmating sexual selection. Here, we investigate molecular evolutionary patterns of genes expressed in the foam gland of Japanese quail (Coturnix japonica), a novel reproductive phenotype. Foam provides an interesting contrast to seminal fluid because it plays a similar functional role, yet is produced, stored, and transferred to females independent of semen. We combined RNA-Seq and comparative genomics to examine evolutionary rates of genes with enriched expression in the foam gland of Japanese quail and those that exhibit enriched expression in two other tissues (testis and liver) and with broadly expressed genes. Overall, we found pronounced heterogeneity in evolutionary rates. Foam gland genes evolved under strong evolutionary constraint, whereas testis genes evolved rapidly and sometimes adaptively. These striking differences were robust to variation in gene expression. Genes with enriched expression in the foam gland did not show major shifts in selective pressure after the quail and chicken lineages split; in contrast, testis-expressed genes experienced a burst of accelerated evolution specifically along the Coturnix lineage. Our work demonstrates that, as a class, genes expressed in the novel foam gland experience different selection regimes than genes expressed in many other tissues producing seminal fluid proteins. Our results also highlight the importance of selective constraint in shaping the evolution of male reproductive genes.