Search

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.

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 (

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.

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.

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.

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. Author summary: Adaptive radiations, which involve both divergent evolution of new traits and recurrent trait evolution, provide insight into the processes that generate and maintain organismal diversity. However, rapid radiations also generate particular challenges for inferring the evolutionary history and mechanistic basis of adaptation and speciation, as multiple processes can cause different parts of the genome to have distinct phylogenetic trees. Thus, inferences about the mode and timing of divergence and the causes of parallel trait evolution require a fine-grained understanding of the flow of genomic variation through time. In this study, we used genome-wide sampling of thousands of genes to re-construct the evolutionary histories of a model plant radiation, the monkeyflowers of Mimulus section Erythranthe. Work over the past half-century has established the parapatric and putatively sister species M. lewisii (bee-pollinated) and M. cardinalis (hummingbird-pollinated, as are three other species in the section) as textbook examples of both rapid speciation via shifts in pollination syndrome and convergent evolution of floral syndromes. Our phylogenomic analyses re-write both of these stories, placing M. cardinalis in a clade with other hummingbird-pollinated taxa and demonstrating that abundant introgression between ancestral lineages as well as in areas of current sympatry contributes to the real (but misleading) affinities between M. cardinalis and M. lewisii. This work illustrates the pervasive influence of gene flow and introgression during adaptive radiation and speciation, and underlines the necessity of a gene-scale and genome-wide phylogenomics framework for understanding trait divergence, even among well-established species. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]