Cis- and Trans-regulatory Effects on Gene Expression in a Natural Population of Drosophila melanogaster

Intraspecific gene regulation in cis- and trans

Changes in gene expression underlie much of evolution and occur via either cis-acting mutations, which lie near the affected gene and act in a context-specific manner, or trans-acting mutations, which may be far from the affected gene and act through diffusible molecules such as transcription factors. A commonly held view is that most expression variation within species is controlled in trans- while expression differences between species are largely controlled in cis-. Here, we summarize recent intraspecific gene regulation studies and find, contrary to this widely held view, that many studies in diverse taxa have revealed a large role for cis-acting mutations underlying expression variation within species. A review of the existing literature also shows that preparations using whole organisms rather than individual tissues may be biased toward identifying trans-regulation. Moreover, we note several examples of predominantly cis-acting regulation in recently diverged populations adapted to different environments. We highlight the challenges of drawing general conclusions from comparisons among studies that use different methodologies and we offer suggestions for studies that will address outstanding questions concerning the evolution of gene regulation.

Genotype-by-Environment Effects of Cis-Variations in the Atgl Promoter Mediate the Divergent Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species

Phenotypic plasticity plays an essential role in adaptive evolution. However, the molecular mechanisms of how genotype-by-environment interaction (G × E) effects shape phenotypic plasticity in marine organisms remain poorly understood. The crucial temperature-responsive trait triacylglycerol (TAG) content and its major gene adipose triglyceride lipase (Atgl) expression have divergent plastic patterns in two congeneric oyster species (Crassostrea gigas and Crassostrea angulata) to adapt to relative-cold/northern and relative-warm/southern habitats, respectively. In this study, eight putative loci were identified in the Atgl promoter region (cis-variations) between wild C. gigas and C. angulata that exhibited differential environmental responsiveness (G × E). The G and G × E effects of each locus were further dissected by measuring the Atgl gene expression of different genotypes in response to temperature changes at the cellular and organismal levels. Two transcription factors, non-environmentally responsive non-POU domain-containing octamer-binding protein (Nono) and environmentally responsive heterogeneous nuclear ribonucleoprotein K (Hnrnpk), were screened for binding to g.-1804 (G locus) and g.-1919 (G + G × E locus), respectively. The specificity of Nono binding to the C. angulata allele mediated the G effects of g.-1804, and the lower environmental sensitivity of Hnrnpk in C. angulata mediated the G × E effects of g.-1919, jointly regulating the trade-offs between higher constitutive and lower plastic expression of Atgl gene expression in C. angulata. This study served as an experimental case to reveal how the genetic variations with G and (or) G × E effects propagate into the divergent pattern of plasticity in environmental adaptive traits, which provides new insights into predicting the adaptability of marine organisms to future climate changes.

Cis-regulatory Variation in Relation to Sex and Sexual Dimorphism in Drosophila melanogaster

Much of sexual dimorphism is likely due to sex-biased gene expression, which results from differential regulation of a genome that is largely shared between males and females. Here, we use allele-specific expression to explore cis-regulatory variation in Drosophila melanogaster in relation to sex. We develop a Bayesian framework to infer the transcriptome-wide joint distribution of cis-regulatory effects across the sexes. We also examine patterns of cis-regulatory variation with respect to two other levels of variation in sexual dimorphism: (i) across genes that vary in their degree of sex-biased expression and (ii) among tissues that vary in their degree of dimorphism (e.g. relatively low dimorphism in heads vs. high dimorphism in gonads). We uncover evidence of widespread cis-regulatory variation in all tissues examined, with female-biased genes being especially enriched for this variation. A sizeable proportion of cis-regulatory variation is inferred to have sex-specific effects, with sex-dependent cis effects being much more frequent in gonads than in heads. Finally, we find some genes where 1 allele contributes to more than 50% of a gene's expression in heterozygous males but <50% of its expression in heterozygous females. Such variants could provide a mechanism for sex-specific dominance reversals, a phenomenon important for sexually antagonistic balancing selection. However, tissue differences in allelic imbalance are approximately as frequent as sex differences, perhaps suggesting that sexual conflict may not be particularly unique in shaping patterns of expression variation.

Unifying framework explaining how parental regulatory divergence can drive gene expression in hybrids and allopolyploids

Hybridization and polyploidy are powerful evolutionary forces, inducing a range of phenotypic outcomes, including non-additive expression, subgenome dominance, deviations in genomic dosage, and transcriptome downsizing. The reasons for these patterns and whether they are universal adaptive responses to genome merging and doubling remain debated. To address this, we develop a thermodynamic model of gene expression based on transcription factor (TF)-promoter binding. Applied to hybridization between species with divergent gene expression levels, cell volumes, or euchromatic ratios, this model distinguishes the effects of hybridization from those of polyploidy. Our results align with empirical observations, suggesting that gene regulation patterns in hybrids and polyploids often stem from the constrained interplay between inherited diverged regulatory networks rather than from subsequent adaptive evolution. In addition, occurrence of certain phenotypic traits depend on specific assumptions about promoter-TF coevolution and their distribution within the hybrid's nucleoplasm, offering new research avenues to understand the underlying mechanisms. In summary, our model explains how the legacy of divergent species directly influences the phenotypic traits of hybrids and allopolyploids.

Pervasive tissue-, genetic background-, and allele-specific gene expression effects in Drosophila melanogaster

The pervasiveness of gene expression variation and its contribution to phenotypic variation and evolution is well known. This gene expression variation is context dependent, with differences in regulatory architecture often associated with intrinsic and environmental factors, and is modulated by regulatory elements that can act in cis (linked) or in trans (unlinked) relative to the genes they affect. So far, little is known about how this genetic variation affects the evolution of regulatory architecture among closely related tissues during population divergence. To address this question, we analyzed gene expression in the midgut, hindgut, and Malpighian tubule as well as microbiome composition in the two gut tissues in four Drosophila melanogaster strains and their F1 hybrids from two divergent populations: one from the derived, European range and one from the ancestral, African range. In both the transcriptome and microbiome data, we detected extensive tissue- and genetic background-specific effects, including effects of genetic background on overall tissue specificity. Tissue-specific effects were typically stronger than genetic background-specific effects, although the two gut tissues were not more similar to each other than to the Malpighian tubules. An examination of allele specific expression revealed that, while both cis and trans effects were more tissue-specific in genes expressed differentially between populations than genes with conserved expression, trans effects were more tissue-specific than cis effects. Despite there being highly variable regulatory architecture, this observation was robust across tissues and genetic backgrounds, suggesting that the expression of trans variation can be spatially fine-tuned as well as or better than cis variation during population divergence and yielding new insights into cis and trans regulatory evolution.

Genetic, Environmental, and Stochastic Components of Lifespan Variability: The Drosophila Paradigm

Lifespan is a complex quantitative trait involving genetic and non-genetic factors as well as the peculiarities of ontogenesis. As with all quantitative traits, lifespan shows considerable variation within populations and between individuals. Drosophila, a favourite object of geneticists, has greatly advanced our understanding of how different forms of variability affect lifespan. This review considers the role of heritable genetic variability, phenotypic plasticity and stochastic variability in controlling lifespan in Drosophila melanogaster. We discuss the major historical milestones in the development of the genetic approach to study lifespan, the breeding of long-lived lines, advances in lifespan QTL mapping, the environmental factors that have the greatest influence on lifespan in laboratory maintained flies, and the mechanisms, by which individual development affects longevity. The interplay between approaches to study ageing and lifespan limitation will also be discussed. Particular attention will be paid to the interaction of different types of variability in the control of lifespan.

Transforming our understanding of species-specific gene regulation

Mechanisms underlying phenotypic divergence across species remain unresolved. In this issue of Cell Genomics, Hansen, Fong, et al.1 systematically dissect human and rhesus macaque gene expression divergence by screening tens of thousands of orthologous elements for enhancer activity in lymphoblastoid cell lines, revealing a much greater role for trans divergence at levels equal to those of cis effects, counter to the prevailing consensus in the field.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Tissue-specific and cis-regulatory changes underlie parallel, adaptive gene expression evolution in house mice

Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence.

Genome-wide expression QTL mapping reveals the highly dynamic regulatory landscape of a major wheat pathogen

BACKGROUND

In agricultural ecosystems, outbreaks of diseases are frequent and pose a significant threat to food security. A successful pathogen undergoes a complex and well-timed sequence of regulatory changes to avoid detection by the host immune system; hence, well-tuned gene regulation is essential for survival. However, the extent to which the regulatory polymorphisms in a pathogen population provide an adaptive advantage is poorly understood.

RESULTS

We used Zymoseptoria tritici, one of the most important pathogens of wheat, to generate a genome-wide map of regulatory polymorphism governing gene expression. We investigated genome-wide transcription levels of 146 strains grown under nutrient starvation and performed expression quantitative trait loci (eQTL) mapping. We identified cis-eQTLs for 65.3% of all genes and the majority of all eQTL loci are within 2kb upstream and downstream of the transcription start site (TSS). We also show that polymorphism in different gene elements contributes disproportionally to gene expression variation. Investigating regulatory polymorphism in gene categories, we found an enrichment of regulatory variants for genes predicted to be important for fungal pathogenesis but with comparatively low effect size, suggesting a separate layer of gene regulation involving epigenetics. We also show that previously reported trait-associated SNPs in pathogen populations are frequently cis-regulatory variants of neighboring genes with implications for the trait architecture.

CONCLUSIONS

Overall, our study provides extensive evidence that single populations segregate large-scale regulatory variation and are likely to fuel rapid adaptation to resistant hosts and environmental change.

A Quantitative, Genome-Wide Analysis in Drosophila Reveals Transposable Elements' Influence on Gene Expression Is Species-Specific

Transposable elements (TEs) are parasite DNA sequences that are able to move and multiply along the chromosomes of all genomes. They can be controlled by the host through the targeting of silencing epigenetic marks, which may affect the chromatin structure of neighboring sequences, including genes. In this study, we used transcriptomic and epigenomic high-throughput data produced from ovarian samples of several Drosophila melanogaster and Drosophila simulans wild-type strains, in order to finely quantify the influence of TE insertions on gene RNA levels and histone marks (H3K9me3 and H3K4me3). Our results reveal a stronger epigenetic effect of TEs on ortholog genes in D. simulans compared with D. melanogaster. At the same time, we uncover a larger contribution of TEs to gene H3K9me3 variance within genomes in D. melanogaster, which is evidenced by a stronger correlation of TE numbers around genes with the levels of this chromatin mark in D. melanogaster. Overall, this work contributes to the understanding of species-specific influence of TEs within genomes. It provides a new light on the considerable natural variability provided by TEs, which may be associated with contrasted adaptive and evolutionary potentials.

A computational investigation of cis-gene regulation in evolution

In biological processes involving gene networks, genes regulate other genes that determine the phenotypic traits. Gene regulation plays an important role in evolutionary dynamics. In a genetic algorithm, a trans-gene regulatory mechanism was shown to speed up adaptation and evolution. Here, we examine the effect of cis-gene regulation on an adaptive system. The model is haploid. A chromosome is partitioned into regulatory loci and structural loci. The regulatory genes regulate the expression and functioning of structural genes via the cis-elements in a probabilistic manner. In the simulation, the change in the allele frequency, the mean population fitness and the efficiency of phenotypic selection are monitored. Cis-gene regulation increases adaption and accelerates the evolutionary process in comparison with the case involving absence of gene regulation. Some special features of the simulation results are as follows. A low ratio of regulatory loci and structural loci gives higher adaptation for fixed total number of loci. Plasticity is advantageous beyond a threshold value. Adaptation is better for large number of total loci when the ratio of regulatory loci to structural loci is one. However, it reaches a saturation beyond which the increase in the total loci is not advantageous. Efficiency of the phenotypic selection is higher for larger value of the initial plasticity.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using ATAC-STARR-seq. In addition to thousands of cis changes, we discover an unexpected number (~10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie >50% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-to regulatory differences between species.

Cis- and Trans-variations of Stearoyl-CoA Desaturase Provide New Insights into the Mechanisms of Diverged Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species

The evolution of phenotypic plasticity plays an essential role in adaptive responses to climate change; however, its regulatory mechanisms in marine organisms which exhibit high phenotypic plasticity still remain poorly understood. The temperature-responsive trait oleic acid content and its major gene stearoyl-CoA desaturase (Scd) expression have diverged in two allopatric congeneric oyster species, cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. In this study, genetic and molecular methods were used to characterize fatty acid desaturation and membrane fluidity regulated by oyster Scd. Sixteen causative single-nucleotide polymorphisms (SNPs) were identified in the promoter/cis-region of the Scd between wild C. gigas and C. angulata. Further functional experiments showed that an SNP (g.-333C [C. gigas allele] >T [C. angulata allele]) may influence Scd transcription by creating/disrupting the binding motif of the positive trans-factor Y-box factor in C. gigas/C. angulata, which mediates the higher/lower constitutive expression of Scd in C. gigas/C. angulata. Additionally, the positive trans-factor sterol-regulatory element-binding proteins (Srebp) were identified to specifically bind to the promoter of Scd in both species, and were downregulated during cold stress in C. gigas compared to upregulated in C. angulata. This partly explains the relatively lower environmental sensitivity (plasticity) of Scd in C. gigas. This study serves as an experimental case to reveal that both cis- and trans-variations shape the diverged pattern of phenotypic plasticity, which provides new insights into the formation of adaptive traits and the prediction of the adaptive potential of marine organisms to future climate change.

Environmental factors have a major effect in shaping the gene expression of Siberian larch in the Altai Mountains of China

The differentiation of gene expression is an important link between genotype and phenotype and has important contributions to species adaptation and ecosystem evolution. As a major component of the world's forests, boreal forests play an important role in regulating the global climate, and the phenology of tree species has been and is undergoing changes during global warming. Here, to understand the impact of global warming on gene expression in boreal forest species, we used PacBio and Illumina sequencing methods to study the transcriptome of natural populations of Siberian larch (Larix sibirica Ledeb.) from the Altai Mountains in Xinjiang, China. We found that populations in this area had low genetic differentiation, but individuals were genetically clustered together when they had close geographic distance. Environmental factors, especially temperature, dominated differential gene expression of Siberian larch, while the contribution of genetic variation is relatively small. We speculate that Siberian larch adapts to changes in temperature and precipitation by altering its own gene expression. These results not only predict the tolerance of boreal forests to higher temperatures in the future, but also inform forest management strategies under global climate change.

Assessing Loss of Regulatory Divergence, Genome-Transcriptome Incongruence, and Preferential Expression Switching in Abaca × Banana Backcrosses

The Musa textilis var. Abuab has high fiber quality (FQ) but is susceptible to abaca bunchy top virus (AbBTV); the Musa balbisiana var. Pacol has low FQ but is resistant against AbBTV. Their backcrosses (BC2 and BC3) possess both desirable traits. Analysis using RNA-seq showed that the regulatory divergence of Abuab and Pacol is largely explained by cis differences with 27.4% and 22.3% if we are to assess it using BC2 and BC3, respectively. Cis differences between the two genotypes are significantly reduced from BC2 to BC3 due to changes in genomic constitution. Trans, on the other hand, is robust to changes in allelic composition. All these are attributed to the loss of heterozygosity in BC3 relative to BC2. Further analysis showed that both backcrosses exhibited genome-wide preferential expression of Pacol- over Abuab-specific alleles, despite the wider genetic presence of the latter in the hybrids. The ratio of the two genotype-specific expressed transcripts and the ratio of their corresponding genetic make-up are significantly disproportionate, a phenomenon that we refer to here as "genome-transcriptome incongruence". We also observed preferential expression switching in which several genes prefer the Abuab- (or Pacol-) specific allele in BC2 but switched to the Pacol- (or Abuab-) specific allele in the BC3 genome.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Alternative splicing in seasonal plasticity and the potential for adaptation to environmental change

Seasonal plasticity is accomplished via tightly regulated developmental cascades that translate environmental cues into trait changes. Little is known about how alternative splicing and other posttranscriptional molecular mechanisms contribute to plasticity or how these mechanisms impact how plasticity evolves. Here, we use transcriptomic and genomic data from the butterfly Bicyclus anynana, a model system for seasonal plasticity, to compare the extent of differential expression and splicing and test how these axes of transcriptional plasticity differ in their potential for evolutionary change. Between seasonal morphs, we find that differential splicing affects a smaller but functionally unique set of genes compared to differential expression. Further, we find strong support for the novel hypothesis that spliced genes are more susceptible than differentially expressed genes to erosion of genetic variation due to selection on seasonal plasticity. Our results suggest that splicing plasticity is especially likely to experience genetic constraints that could affect the potential of wild populations to respond to rapidly changing environments.

Evolution of a New Testis-Specific Functional Promoter Within the Highly Conserved Map2k7 Gene of the Mouse

Map2k7 (synonym Mkk7) is a conserved regulatory kinase gene and a central component of the JNK signaling cascade with key functions during cellular differentiation. It shows complex transcription patterns, and different transcript isoforms are known in the mouse (Mus musculus). We have previously identified a newly evolved testis-specific transcript for the Map2k7 gene in the subspecies M. m. domesticus. Here, we identify the new promoter that drives this transcript and find that it codes for an open reading frame (ORF) of 50 amino acids. The new promoter was gained in the stem lineage of closely related mouse species but was secondarily lost in the subspecies M. m. musculus and M. m. castaneus. A single mutation can be correlated with its transcriptional activity in M. m. domesticus, and cell culture assays demonstrate the capability of this mutation to drive expression. A mouse knockout line in which the promoter region of the new transcript is deleted reveals a functional contribution of the newly evolved promoter to sperm motility and the spermatid transcriptome. Our data show that a new functional transcript (and possibly protein) can evolve within an otherwise highly conserved gene, supporting the notion of regulatory changes contributing to the emergence of evolutionary novelties.

Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations

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 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 superabundance 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 coexpressed 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.

Regulatory regions in natural transposable element insertions drive interindividual differences in response to immune challenges in Drosophila

Background

Variation in gene expression underlies interindividual variability in relevant traits including immune response. However, the genetic variation responsible for these gene expression changes remains largely unknown. Among the non-coding variants that could be relevant, transposable element insertions are promising candidates as they have been shown to be a rich and diverse source of cis-regulatory elements.

Results

In this work, we use a population genetics approach to identify transposable element insertions likely to increase the tolerance of Drosophila melanogaster to bacterial infection by affecting the expression of immune-related genes. We identify 12 insertions associated with allele-specific expression changes in immune-related genes. We experimentally validate three of these insertions including one likely to be acting as a silencer, one as an enhancer, and one with a dual role as enhancer and promoter. The direction in the change of gene expression associated with the presence of several of these insertions is consistent with an increased survival to infection. Indeed, for one of the insertions, we show that this is the case by analyzing both natural populations and CRISPR/Cas9 mutants in which the insertion is deleted from its native genomic context.

Conclusions

We show that transposable elements contribute to gene expression variation in response to infection in D. melanogaster and that this variation is likely to affect their survival capacity. Because the role of transposable elements as regulatory elements is not restricted to Drosophila, transposable elements are likely to play a role in immune response in other organisms as well.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02471-3.

Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain

Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F₁ hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.

Unraveling regulatory divergence, heterotic malleability, and allelic imbalance switching in rice due to drought stress

The indica ecotypes, IR64, an elite drought-susceptible variety adapted to irrigated ecosystem, and Apo (IR55423-01 or NSIC RC9), a moderate drought-tolerant upland genotype together with their hybrid (IR64 × Apo) were exposed to non- and water-stress conditions. By sequencing (RNA-seq) these genotypes, we were able to map genes diverging in cis and/or trans factors. Under non-stress condition, cis dominantly explains (11.2%) regulatory differences, followed by trans (8.9%). Further analysis showed that water-limiting condition largely affects trans and cis + trans factors. On the molecular level, cis and/or trans regulatory divergence explains their genotypic differences and differential drought response. Between the two parental genotypes, Apo appears to exhibit more photosynthetic efficiency even under water-limiting condition and is ascribed to trans. Statistical analyses showed that regulatory divergence is significantly influenced by environmental conditions. Likewise, the mode of parental expression inheritance which drives heterosis (HET) is significantly affected by environmental conditions indicating the malleability of heterosis to external factors. Further analysis revealed that the HET class, dominance, was significantly enriched under water-stress condition. We also identified allelic imbalance switching in which several genes prefer IR64- (or Apo-) specific allele under non-stress condition but switched to Apo- (or IR64-) specific allele when exposed to water-stress condition.

The mode of expression divergence in Drosophila fat body is infection-specific

Transcription is controlled by interactions of cis-acting DNA elements with diffusible trans-acting factors. Changes in cis or trans factors can drive expression divergence within and between species, and their relative prevalence can reveal the evolutionary history and pressures that drive expression variation. Previous work delineating the mode of expression divergence in animals has largely used whole-body expression measurements in one condition. Because cis-acting elements often drive expression in a subset of cell types or conditions, these measurements may not capture the complete contribution of cis-acting changes. Here, we quantify the mode of expression divergence in the Drosophila fat body, the primary immune organ, in several conditions, using two geographically distinct lines of D. melanogaster and their F1 hybrids. We measured expression in the absence of infection and in infections with Gram-negative S. marcescens or Gram-positive E. faecalis bacteria, which trigger the two primary signaling pathways in the Drosophila innate immune response. The mode of expression divergence strongly depends on the condition, with trans-acting effects dominating in response to Gram-negative infection and cis-acting effects dominating in Gram-positive and preinfection conditions. Expression divergence in several receptor proteins may underlie the infection-specific trans effects. Before infection, when the fat body has a metabolic role, there are many compensatory effects, changes in cis and trans that counteract each other to maintain expression levels. This work shows that within a single tissue, the mode of expression divergence varies between conditions and suggests that these differences reflect the diverse evolutionary histories of host-pathogen interactions.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.

Genetic association study of a novel indel polymorphism in HSPA1B with the risk of sudden cardiac death in the Chinese populations

Sudden cardiac death (SCD) has become a global problem due to its high mortality in the general population. Identification of genetic factors predisposed to SCD is significant since it enables genetic testing that would contribute to molecular diagnosis and risk stratification of SCD. It has been reported that HSPA1B gene mutations might be related with SCD. In this study, based on candidate-gene-based approach and systematic screening strategy, a 5-base pair insertion/deletion (Indel) polymorphism (rs3036297) in the 3'UTR of HSPA1B gene was selected to perform a case-control study aiming to investigate its association with SCD susceptibility in Chinese populations. Logistic regression analysis showed that the insertion allele of rs3036297 was correlated with a comparatively lower risk for SCD [OR=0.58, 95%CI=0.43-0.77, P=1.28×10^-4] compared with the deletion allele. Luciferase activity assay indicated that HSPA1B expression could be regulated by rs3036297 through interfering binding with miR-134-5p. Furthermore, analysis of database from Haploreg and GTEx revealed that the rs3036297 variant was involved in potential cis-regulatory element with the promoter of HLA-DRB5 through a long-range interaction and the deletion allele of rs3036297 increased HLA-DRB5 expression. In conclusion, the rs3036297 variant may regulate HSPA1B expression via a mechanism of miRNA binding and HLA-DRB5 expression via a long-range promoter interaction through which contributed to SCD susceptibility. Therefore, rs3036297 would be a potential marker for molecular diagnosis and genetic counseling of SCD.

Parallel evolution of apetalous lineages within the buttercup family (Ranunculaceae): outward expansion of AGAMOUS1, rather than disruption of APETALA3-3

Complete loss of petals, or becoming apetalous, has occurred independently in many flowering plant lineages. However, the mechanisms underlying the parallel evolution of naturally occurring apetalous lineages remain largely unclear. Here, by sampling representatives of all nine apetalous genera/tribes of the family Ranunculaceae and conducting detailed morphological, expression, molecular evolutionary and functional studies, we investigate the mechanisms underlying parallel petal losses. We found that while non-expression/downregulation of the petal identity gene APETALA3-3 (AP3-3) is tightly associated with complete petal losses, disruptions of the AP3-3 orthologs were unlikely to be the real causes for the parallel evolution of apetalous lineages. We also found that, compared with their close petalous relatives, naturally occurring apetalous taxa usually bear slightly larger numbers of stamens, whereas the number of sepals remains largely unchanged, suggestive of petal-to-stamen rather than petal-to-sepal transformations. In addition, in the recently originated apetalous genus Enemion, the petal-to-stamen transformations have likely been caused by the mutations that led to the elevation and outward expansion of the expression of the C-function gene, AGAMOUS1 (AG1). Our results not only provide a general picture of parallel petal losses within the Ranunculaceae but also help understand the mechanisms underlying the independent originations of other apetalous lineages.

microRNA expression variation as a potential molecular mechanism contributing to adaptation to hydrogen sulphide

microRNAs (miRNAs) are post-transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression in Poecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H₂ S). We found that P. mexicana has a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute to P. mexicana adaptation to sulphidic environments. In silico identification of differentially expressed miRNA-mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H₂ S-rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H₂ S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H₂ S adaptation by promoting functions needed for survival and reducing functions affected by H₂ S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H₂ S. Overall, this study provides a critical stepping-stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

Polymorphism and Divergence of Novel Gene Expression Patterns in Drosophila melanogaster

Transcriptomes may evolve by multiple mechanisms, including the evolution of novel genes, the evolution of transcript abundance, and the evolution of cell, tissue, or organ expression patterns. Here, we focus on the last of these mechanisms in an investigation of tissue and organ shifts in gene expression in Drosophila melanogaster. In contrast to most investigations of expression evolution, we seek to provide a framework for understanding the mechanisms of novel expression patterns on a short population genetic timescale. To do so, we generated population samples of D. melanogaster transcriptomes from five tissues: accessory gland, testis, larval salivary gland, female head, and first-instar larva. We combined these data with comparable data from two outgroups to characterize gains and losses of expression, both polymorphic and fixed, in D. melanogaster We observed a large number of gain- or loss-of-expression phenotypes, most of which were polymorphic within D. melanogaster Several polymorphic, novel expression phenotypes were strongly influenced by segregating cis-acting variants. In support of previous literature on the evolution of novelties functioning in male reproduction, we observed many more novel expression phenotypes in the testis and accessory gland than in other tissues. Additionally, genes showing novel expression phenotypes tend to exhibit greater tissue-specific expression. Finally, in addition to qualitatively novel expression phenotypes, we identified genes exhibiting major quantitative expression divergence in the D. melanogaster lineage.

Cis and trans effects differentially contribute to the evolution of promoters and enhancers

BACKGROUND

Gene expression differences between species are driven by both cis and trans effects. Whereas cis effects are caused by genetic variants located on the same DNA molecule as the target gene, trans effects are due to genetic variants that affect diffusible elements. Previous studies have mostly assessed the impact of cis and trans effects at the gene level. However, how cis and trans effects differentially impact regulatory elements such as enhancers and promoters remains poorly understood. Here, we use massively parallel reporter assays to directly measure the transcriptional outputs of thousands of individual regulatory elements in embryonic stem cells and measure cis and trans effects between human and mouse.

RESULTS

Our approach reveals that cis effects are widespread across transcribed regulatory elements, and the strongest cis effects are associated with the disruption of motifs recognized by strong transcriptional activators. Conversely, we find that trans effects are rare but stronger in enhancers than promoters and are associated with a subset of transcription factors that are differentially expressed between human and mouse. While we find that cis-trans compensation is common within promoters, we do not see evidence of widespread cis-trans compensation at enhancers. Cis-trans compensation is inversely correlated with enhancer redundancy, suggesting that such compensation may often occur across multiple enhancers.

CONCLUSIONS

Our results highlight differences in the mode of evolution between promoters and enhancers in complex mammalian genomes and indicate that studying the evolution of individual regulatory elements is pivotal to understand the tempo and mode of gene expression evolution.

Speciation and changes in male gene expression in Drosophila

It has long been acknowledged that changes in the regulation of gene expression may account for major organismal differences. However, we still do not fully understand how changes in gene expression evolve and how do such changes influence organisms' differences. We are even less aware of the impact such changes might have in restricting gene flow between species. Here, we focus on studies of gene expression and speciation in the Drosophila model. We review studies that have identified gene interactions in post-mating reproductive isolation and speciation, particularly those that modulate male gene expression. We also address studies that have experimentally manipulated changes in gene expression to test their effect in post-mating reproductive isolation. We highlight the need for a more in-depth analysis of the role of selection causing disrupted gene expression of such candidate genes in sterile/inviable hybrids. Moreover, we discuss the relevance to incorporate more routinely assays that simultaneously evaluate the potential effects of environmental factors and genetic background in modulating plastic responses in male genes and their potential role in speciation.

The contribution of cis- and trans-acting variants to gene regulation in wild and domesticated barley under cold stress and control conditions

Barley, like other crops, has experienced a series of genetic changes that have impacted its architecture and growth habit to suit the needs of humans, termed the domestication syndrome. Domestication also resulted in a concomitant bottleneck that reduced sequence diversity in genes and regulatory regions. Little is known about regulatory changes resulting from domestication in barley. We used RNA sequencing to examine allele-specific expression in hybrids between wild and domesticated barley. Our results show that most genes have conserved regulation. In contrast to studies of allele-specific expression in interspecific hybrids, we find almost a complete absence of trans effects. We also find that cis regulation is largely stable in response to short-term cold stress. Our study has practical implications for crop improvement using wild relatives. Genes regulated in cis are more likely to be expressed in a new genetic background at the same level as in their native background.

Gene duplication and transposition of mobile elements drive evolution of the Rpv3 resistance locus in grapevine

A wild grape haplotype (Rpv3-1) confers resistance to Plasmopara viticola. We mapped the causal factor for resistance to an interval containing a TIR-NB-LRR (TNL) gene pair that originated 1.6-2.6 million years ago by a tandem segmental duplication. Transient coexpression of the TNL pair in Vitis vinifera leaves activated pathogen-induced necrosis and reduced sporulation compared with control leaves. Even though transcripts of the TNL pair from the wild haplotype appear to be partially subject to nonsense-mediated mRNA decay, mature mRNA levels in a homozygous resistant genotype were individually higher than the mRNA trace levels observed for the orthologous single-copy TNL in sensitive genotypes. Allelic expression imbalance in a resistant heterozygote confirmed that cis-acting regulatory variation promotes expression in the wild haplotype. The movement of transposable elements had a major impact on the generation of haplotype diversity, altering the DNA context around similar TNL coding sequences and the GC-content in their proximal 5'-intergenic regions. The wild and domesticated haplotypes also diverged in conserved single-copy intergenic DNA, but the highest divergence was observed in intraspecific and not in interspecific comparisons. In this case, introgression breeding did not transgress the genetic boundaries of the domesticated species, because haplotypes present in modern varieties sometimes predate speciation events between wild and cultivated species.

Cloudy with a Chance of Insights: Context Dependent Gene Regulation and Implications for Evolutionary Studies

Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic evolution. An exceptional contribution of cis-regulatory divergence has been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and stage-specific expression. We argue that the regulatory context must be considered in comparative expression studies.

Comparative transcriptomics of 3 high-altitude passerine birds and their low-altitude relatives

High-altitude environments present strong stresses for living organisms, which have driven striking phenotypic and genetic adaptations. While previous studies have revealed multiple genetic adaptations in high-altitude species, how evolutionary history (i.e., phylogenetic background) contributes to similarity in genetic adaptations to high-altitude environments is largely unknown, in particular in a group of birds. We explored this in 3 high-altitude passerine birds from the Qinghai-Tibet Plateau and their low-altitude relatives in lowland eastern China. We generated transcriptomic data for 5 tissues across these species and compared sequence changes and expression shifts between high- and low-altitude pairs. Sequence comparison revealed that similarity in all 3 high-altitude species was high for genes under positive selection (218 genes) but low in amino acid substitutions (only 4 genes sharing identical amino acid substitutions). Expression profiles for all genes identified a tissue-specific expression pattern (i.e., all species clustered by tissue). By contrast, an altitude-related pattern was observed in genes differentially expressed between all 3 species pairs and genes associated with altitude, suggesting that the high-altitude environment may drive similar expression shifts in the 3 high-altitude species. Gene expression level, gene connectivity, and the interactions of these 2 factors with altitude were correlated with evolutionary rates. Our results provide evidence for how gene sequence changes and expression shifts work in a concerted way in a group of high-altitude birds, leading to similar evolution routes in response to high-altitude environmental stresses.

Chromatin Evolution-Key Innovations Underpinning Morphological Complexity

The history of life consists of a series of major evolutionary transitions, including emergence and radiation of complex multicellular eukaryotes from unicellular ancestors. The cells of multicellular organisms, with few exceptions, contain the same genome, however, their organs are composed of a variety of cell types that differ in both structure and function. This variation is largely due to the transcriptional activity of different sets of genes in different cell types. This indicates that complex transcriptional regulation played a key role in the evolution of complexity in eukaryotes. In this review, we summarize how gene duplication and subsequent evolutionary innovations, including the structural evolution of nucleosomes and chromatin-related factors, contributed to the complexity of the transcriptional system and provided a basis for morphological diversity.

Determining Allele-Specific Protein Expression (ASPE) Using a Novel Quantitative Concatamer Based Proteomics Method

Measuring allele-specific expression (ASE) is a powerful approach for identifying cis-regulatory genetic variants. Here, we developed a novel targeted proteomics method for the quantification of allele-specific protein expression (ASPE) based on scheduled parallel reaction monitoring (PRM) with a heavy stable isotope-labeled quantitative concatamer (QconCAT) internal protein standard. This strategy was applied to the determination of the ASPE of UGT2B15 in human livers using the common UGT2B15 nonsynonymous variant rs1902023 (i.e., Y85D) as the marker to differentiate expressions from the two alleles. The QconCAT standard contains both the wild-type tryptic peptide and the Y85D mutant peptide at a ratio of 1:1 to ensure accurate measurement of the ASPE of UGT2B15. The results from 18 UGT2B15 Y85D heterozygotes revealed that the ratios between the wild-type Y allele and the mutant D allele varied from 0.60 to 1.46, indicating the presence of cis-regulatory variants. In addition, we observed no significant correlations between the ASPE and mRNA ASE of UGT2B15, suggesting the involvement of different cis-acting variants in regulating the transcription and translation processes of the gene. This novel ASPE approach provides a powerful tool for capturing cis-genetic variants involved in post-transcription processes, an important yet understudied area of research.

Gene Regulatory Variation in Drosophila melanogaster Renal Tissue

Genetic variation influencing levels of gene expression is abundant in natural populations, and may exert its effects through complex mechanisms that depend on an organism's genetic background and the tissue in which expression is measured. We investigated natural variation in gene expression in the Malpighian tubules of three inbred Drosophila melanogaster strains and their F₁ hybrids. One of the strains was from a population in the species' ancestral range (Zambia), while the other two were from a more recently derived population (Sweden). Although closely related, the two Swedish strains differed greatly in terms of their expression inheritance when hybridized with the Zambian strain, with one Swedish strain showing a large excess of genes with recessive expression inheritance, as well as a large number of genes with overdominant inheritance. Although most expression variation could be attributed to trans-regulation, there were ∼200 genes that showed allele-specific expression differences in each of the between-population hybrids, indicating that cis-regulation contributes as well. The cis-regulated genes were enriched with cytochrome P450 genes, and the upstream regions of six of these genes were incorporated into transgenic reporter gene constructs to test their effects on expression. Differential expression was observed for five of the six reporter genes in the Malpighian tubule, suggesting that a large proportion of cis-regulatory variation lies directly upstream of the affected gene. In most cases, the differential expression was specific to the Malpighian tubule or greater in this tissue than in the rest of the body, highlighting the importance of single-tissue studies of gene expression variation.

The Evolution of Gene Expression in cis and trans

There is abundant variation in gene expression between individuals, populations, and species. The evolution of gene regulation and expression within and between species is thought to frequently contribute to adaptation. Yet considerable evidence suggests that the primary evolutionary force acting on variation in gene expression is stabilizing selection. We review here the results of recent studies characterizing the evolution of gene expression occurring in cis (via linked polymorphisms) or in trans (through diffusible products of other genes) and their contribution to adaptation and response to the environment. We review the evidence for buffering of variation in gene expression at the level of both transcription and translation, and the possible mechanisms for this buffering. Lastly, we summarize unresolved questions about the evolution of gene regulation.

Convergent evolution of gene expression in two high-toothed stickleback populations

Changes in developmental gene regulatory networks enable evolved changes in morphology. These changes can be in cis regulatory elements that act in an allele-specific manner, or changes to the overall trans regulatory environment that interacts with cis regulatory sequences. Here we address several questions about the evolution of gene expression accompanying a convergently evolved constructive morphological trait, increases in tooth number in two independently derived freshwater populations of threespine stickleback fish (Gasterosteus aculeatus). Are convergently evolved cis and/or trans changes in gene expression associated with convergently evolved morphological evolution? Do cis or trans regulatory changes contribute more to gene expression changes accompanying an evolved morphological gain trait? Transcriptome data from dental tissue of ancestral low-toothed and two independently derived high-toothed stickleback populations revealed significantly shared gene expression changes that have convergently evolved in the two high-toothed populations. Comparing cis and trans regulatory changes using phased gene expression data from F1 hybrids, we found that trans regulatory changes were predominant and more likely to be shared among both high-toothed populations. In contrast, while cis regulatory changes have evolved in both high-toothed populations, overall these changes were distinct and not shared among high-toothed populations. Together these data suggest that a convergently evolved trait can occur through genetically distinct regulatory changes that converge on similar trans regulatory environments.