This content will become publicly available on September 26, 2026

Title: Craniofacial-specific transcriptomics uncovers novel genes underlying jaw divergence in dietary specialist pupfishes
Abstract Changes in gene expression underlie most phenotypic differences among closely related species. While previous studies in model systems have identified conserved genes and pathways involved in craniofacial evolution, less is known about gene expression differences associated with craniofacial divergence in rapidly evolving species. Here, we investigate craniofacial-specific gene expression in a nascent adaptive radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas, which includes 3 trophic specialists with highly divergent craniofacial morphologies (two scale-eaters and a molluscivore) derived from an ancestral Caribbean-wide generalist. We compared gene expression in the most morphologically divergent craniofacial region with the relatively conserved caudal region across 5 Cyprinodon species and 9 populations. We focused on the hatchling stage, the earliest developmental stage at which craniofacial differences among species are evident. Our approach revealed a large proportion of differentially expressed genes (DEGs) found exclusively in the craniofacial region of the specialists only. By intersecting these specialist-specific craniofacial-exclusive genes with genomic regions harboring fixed single-nucleotide variants under selection in the specialists, we identified 14 candidate genes. We confirmed novel craniofacial expression for 2 of these candidates, pycr3 and atp8a1, genes not previously associated with craniofacial development or function, in hatchlings using in situ mRNA hybridization and observed species-specific differences in the pharyngeal arches and craniofacial muscles, respectively. Our findings demonstrate how an “evolutionary mutant” model can reveal novel gene expression patterns, highlighting the power of integrating tissue-species transcriptomics with speciation genomics to identify novel regulators of craniofacial evolution.  more » « less
Award ID(s):
1938571
PAR ID:
10668706
Author(s) / Creator(s):
; ;
Editor(s):
Peichel, C
Publisher / Repository:
Genetics Society of America
Date Published:
Journal Name:
GENETICS
Volume:
231
Issue:
4
ISSN:
1943-2631
Format(s):
Medium: X
Associated Dataset(s):
View Associated Dataset(s) >>
Sponsoring Org:
National Science Foundation
More Like this
  1. Wittkopp, Patricia (Ed.)
    Abstract Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists—a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species—only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development. 
    more » « less
  2. Understanding the genetic basis of novel adaptations in new species is a fundamental question in biology. Here we demonstrate a new role for galr2 in vertebrate craniofacial development using an adaptive radiation of trophic specialist pupfishes endemic to San Salvador Island, Bahamas. We confirmed the loss of a putative Sry transcription factor binding site upstream of galr2 in scale-eating pupfish and found significant spatial differences in galr2 expression among pupfish species in Meckel's cartilage usingin situhybridization chain reaction (HCR). We then experimentally demonstrated a novel role for Galr2 in craniofacial development by exposing embryos to Garl2-inhibiting drugs. Galr2-inhibition reduced Meckel's cartilage length and increased chondrocyte density in both trophic specialists but not in the generalist genetic background. We propose a mechanism for jaw elongation in scale-eaters based on the reduced expression of galr2 due to the loss of a putative Sry binding site. Fewer Galr2 receptors in the scale-eater Meckel's cartilage may result in their enlarged jaw lengths as adults by limiting opportunities for a circulating Galr2 agonist to bind to these receptors during development. Our findings illustrate the growing utility of linking candidate adaptive SNPs in non-model systems with highly divergent phenotypes to novel vertebrate gene functions. 
    more » « less
  3. Wittkopp, Patricia (Ed.)
    Abstract Genes involved in spermatogenesis tend to evolve rapidly, but we lack a clear understanding of how protein sequences and patterns of gene expression evolve across this complex developmental process. We used fluorescence-activated cell sorting (FACS) to generate expression data for early (meiotic) and late (postmeiotic) cell types across 13 inbred strains of mice (Mus) spanning ∼7 My of evolution. We used these comparative developmental data to investigate the evolution of lineage-specific expression, protein-coding sequences, and expression levels. We found increased lineage specificity and more rapid protein-coding and expression divergence during late spermatogenesis, suggesting that signatures of rapid testis molecular evolution are punctuated across sperm development. Despite strong overall developmental parallels in these components of molecular evolution, protein and expression divergences were only weakly correlated across genes. We detected more rapid protein evolution on the X chromosome relative to the autosomes, whereas X-linked gene expression tended to be relatively more conserved likely reflecting chromosome-specific regulatory constraints. Using allele-specific FACS expression data from crosses between four strains, we found that the relative contributions of different regulatory mechanisms also differed between cell types. Genes showing cis-regulatory changes were more common late in spermatogenesis, and tended to be associated with larger differences in expression levels and greater expression divergence between species. In contrast, genes with trans-acting changes were more common early and tended to be more conserved across species. Our findings advance understanding of gene evolution across spermatogenesis and underscore the fundamental importance of developmental context in molecular evolutionary studies. 
    more » « less
  4. Evolutionary transitions to a social lifestyle in insects are associated with lineage-specific changes in gene expression, but the key nodes that drive these regulatory changes are unknown. We examined the relationship between social organization and lineage-specific microRNAs (miRNAs). Genome scans across 12 bee species showed that miRNA copy-number is mostly conserved and not associated with sociality. However, deep sequencing of small RNAs in six bee species revealed a substantial proportion (20–35%) of detected miRNAs had lineage-specific expression in the brain, 24–72% of which did not have homologues in other species. Lineage-specific miRNAs disproportionately target lineage-specific genes, and have lower expression levels than shared miRNAs. The predicted targets of lineage-specific miRNAs are not enriched for genes with caste-biased expression or genes under positive selection in social species. Together, these results suggest that novel miRNAs may coevolve with novel genes, and thus contribute to lineage-specific patterns of evolution in bees, but do not appear to have significant influence on social evolution. Our analyses also support the hypothesis that many new miRNAs are purged by selection due to deleterious effects on mRNA targets, and suggest genome structure is not as influential in regulating bee miRNA evolution as has been shown for mammalian miRNAs. 
    more » « less
  5. Skeletal muscles are essential for movement, supporting a wide range of locomotor behaviors. Muscle tissue is composed of multiple cell types including “fast” and “slow” myofibers, whose contractile properties are largely influenced by selective expression of myosin heavy chain (MyHC) isoforms. While ‘super-enhancers’ regulating MyHC gene clusters have been identified, the cis-regulatory elements (CREs) controlling non-MyHC genes important to myofiber physiology remain less defined. Here, we profile the regulatory landscape of two pairs of mouse hind limb muscles differing in MyHC expression at a late embryonic (E18.5) and adult time point to identify candidate CREs that may regulate genes important to myofiber type. Gene expression and chromatin accessibility analyses revealed that epigenetic differences at E18.5 largely reflect limb patterning, whereas adult differences reflect myofiber differentiation. We identified thousands of differentially accessible regions that may regulate genes important for muscle development, muscle biology, and myofiber identity. Among these, twelve conserved, muscle-specific CREs associated with myofiber type were tested for regulatory activity. Nine enhanced and three reduced gene activity in vitro, although their phenotypic effects remain unknown. By profiling multiple muscles across two time points, our study extends current understanding of conserved, muscle-specific CREs that regulate gene expression during myogenesis. 
    more » « less