Molecular and evolutionary processes generating variation in gene expression.

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Active compensation for changes in TDH3 expression mediated by direct regulators of TDH3 in Saccharomyces cerevisiae.

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralog TDH2. TDH2 is upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. TDH1, a second and more distantly related paralog of TDH3, has diverged in its regulation and is upregulated by another mechanism. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

Network Topology Can Explain Differences in Pleiotropy Between Cis- and Trans-regulatory Mutations.

A mutation's degree of pleiotropy (i.e., the number of traits it alters) is predicted to impact the probability of the mutation being detrimental to fitness. For mutations that impact gene expression, mutations acting in cis have been hypothesized to generally be less pleiotropic than mutations affecting the same gene's expression in trans, suggesting that cis-regulatory mutations should be less deleterious and more likely to fix over evolutionary time. Here, we use expression and fitness data from Saccharomyces cerevisiae gene deletion strains to test these hypotheses. By treating deletion of each gene as a cis-regulatory mutation affecting its own expression and deletions of other genes affecting expression of this focal gene as trans-regulatory mutations, we find that cis-acting mutations do indeed tend to be less pleiotropic than trans-acting mutations affecting expression of the same gene. This pattern was observed for the vast majority of genes in the data set and could be explained by the topology of the regulatory network controlling gene expression. Comparing the fitness of cis- and trans-acting mutations affecting expression of the same gene also confirmed that trans-acting deletions tend to be more deleterious. These findings provide strong support for pleiotropy playing a role in the preferential fixation of cis-regulatory alleles over evolutionary time.

Differential Grainy head binding correlates with variation in chromatin structure and gene expression in Drosophila melanogaster.

Phenotypic evolution is often caused by variation in gene expression resulting from altered gene regulatory mechanisms. Genetic variation affecting chromatin remodeling has been identified as a potential source of variable gene expression; however, the roles of specific chromatin remodeling factors remain unclear. Here, we address this knowledge gap by examining the relationship between variation in gene expression, variation in chromatin structure, and variation in binding of the pioneer factor Grainy head between imaginal wing discs of two divergent strains of Drosophila melanogaster and their F1 hybrid. We find that (1) variation in Grainy head binding is mostly due to sequence changes that act in cis but are located outside of the canonical Grainy head binding motif, (2) variation in Grainy head binding correlates with changes in chromatin accessibility, and (3) this variation in chromatin accessibility, coupled with variation in Grainy head binding, correlates with variation in gene expression in some cases but not others. Interactions among these three molecular layers is complex, but these results suggest that genetic variation affecting the binding of pioneer factors contributes to variation in chromatin remodeling and the evolution of gene expression.

Empirical measures of mutational effects define neutral models of regulatory evolution in Saccharomyces cerevisiae.

Understanding how phenotypes evolve requires disentangling the effects of mutation generating new variation from the effects of selection filtering it. Tests for selection frequently assume that mutation introduces phenotypic variation symmetrically around the population mean, yet few studies have tested this assumption by deeply sampling the distributions of mutational effects for particular traits. Here, we examine distributions of mutational effects for gene expression in the budding yeast Saccharomyces cerevisiae by measuring the effects of thousands of point mutations introduced randomly throughout the genome. We find that the distributions of mutational effects differ for the 10 genes surveyed and are inconsistent with normality. For example, all 10 distributions of mutational effects included more mutations with large effects than expected for normally distributed phenotypes. In addition, some genes also showed asymmetries in their distribution of mutational effects, with new mutations more likely to increase than decrease the gene's expression or vice versa. Neutral models of regulatory evolution that take these empirically determined distributions into account suggest that neutral processes may explain more expression variation within natural populations than currently appreciated.

Distinct genetic architectures underlie divergent thorax, leg, and wing pigmentation between Drosophila elegans and D. gunungcola.

Pigmentation divergence between Drosophila species has emerged as a model trait for studying the genetic basis of phenotypic evolution, with genetic changes contributing to pigmentation differences often mapping to genes in the pigment synthesis pathway and their regulators. These studies of Drosophila pigmentation have tended to focus on pigmentation changes in one body part for a particular pair of species, but changes in pigmentation are often observed in multiple body parts between the same pair of species. The similarities and differences of genetic changes responsible for divergent pigmentation in different body parts of the same species thus remain largely unknown. Here we compare the genetic basis of pigmentation divergence between Drosophila elegans and D. gunungcola in the wing, legs, and thorax. Prior work has shown that regions of the genome containing the pigmentation genes yellow and ebony influence the size of divergent male-specific wing spots between these two species. We find that these same two regions of the genome underlie differences in leg and thorax pigmentation; however, divergent alleles in these regions show differences in allelic dominance and epistasis among the three body parts. These complex patterns of inheritance can be explained by a model of evolution involving tissue-specific changes in the expression of Yellow and Ebony between D. elegans and D. gunungcola.

Selection on noise constrains variation in a eukaryotic promoter.

Genetic variation segregating within a species reflects the combined activities of mutation, selection, and genetic drift. In the absence of selection, polymorphisms are expected to be a random subset of new mutations; thus, comparing the effects of polymorphisms and new mutations provides a test for selection. When evidence of selection exists, such comparisons can identify properties of mutations that are most likely to persist in natural populations. Here we investigate how mutation and selection have shaped variation in a cis-regulatory sequence controlling gene expression by empirically determining the effects of polymorphisms segregating in the TDH3 promoter among 85 strains of Saccharomyces cerevisiae and comparing their effects to a distribution of mutational effects defined by 236 point mutations in the same promoter. Surprisingly, we find that selection on expression noise (that is, variability in expression among genetically identical cells) appears to have had a greater impact on sequence variation in the TDH3 promoter than selection on mean expression level. This is not necessarily because variation in expression noise impacts fitness more than variation in mean expression level, but rather because of differences in the distributions of mutational effects for these two phenotypes. This study shows how systematically examining the effects of new mutations can enrich our understanding of evolutionary mechanisms. It also provides rare empirical evidence of selection acting on expression noise.

Genetic architecture of a body colour cline in Drosophila americana.

Phenotypic variation within a species is often structured geographically in clines. In Drosophila americana, a longitudinal cline for body colour exists within North America that appears to be due to local adaptation. The tan and ebony genes have been hypothesized to contribute to this cline, with alleles of both genes that lighten body colour found in D. americana. These alleles are similar in sequence and function to the allele fixed in D. americana's more lightly pigmented sister species, Drosophila novamexicana. Here, we examine the frequency and geographic distribution of these D. novamexicana-like alleles in D. americana. Among alleles from over 100 strains of D. americana isolated from 21 geographic locations, we failed to identify additional alleles of tan or ebony with as much sequence similarity to D. novamexicana as the D. novamexicana-like alleles previously described. However, using genetic analysis of 51 D. americana strains derived from 20 geographic locations, we identified one new allele of ebony and one new allele of tan segregating in D. americana that are functionally equivalent to the D. novamexicana allele. An additional 5 alleles of tan also showed marginal evidence of functional similarity. Given the rarity of these alleles, however, we conclude that they are unlikely to be driving the pigmentation cline. Indeed, phenotypic distributions of the 51 backcross populations analysed indicate a more complex genetic architecture, with diversity in the number and effects of loci altering pigmentation observed both within and among populations of D. americana. This genetic heterogeneity poses a challenge to association studies and genomic scans for clinal variation, but might be common in natural populations.

Effects of mutation and selection on plasticity of a promoter activity in Saccharomyces cerevisiae.

Phenotypic plasticity is an evolvable property of biological systems that can arise from environment-specific regulation of gene expression. To better understand the evolutionary and molecular mechanisms that give rise to plasticity in gene expression, we quantified the effects of 235 single-nucleotide mutations in the Saccharomyces cerevisiae TDH3 promoter (PTDH3 ) on the activity of this promoter in media containing glucose, galactose, or glycerol as a carbon source. We found that the distributions of mutational effects differed among environments because many mutations altered the plastic response exhibited by the wild-type allele. Comparing the effects of these mutations with the effects of 30 PTDH3 polymorphisms on expression plasticity in the same environments provided evidence of natural selection acting to prevent the plastic response in PTDH3 activity between glucose and galactose from becoming larger. The largest changes in expression plasticity were observed between fermentable (glucose or galactose) and nonfermentable (glycerol) carbon sources and were caused by mutations located in the RAP1 and GCR1 transcription factor binding sites. Mutations altered expression plasticity most frequently between the two fermentable environments, with mutations causing significant changes in plasticity between glucose and galactose distributed throughout the promoter, suggesting they might affect chromatin structure. Taken together, these results provide insight into the molecular mechanisms underlying gene-by-environment interactions affecting gene expression as well as the evolutionary dynamics affecting natural variation in plasticity of gene expression.

Structure of the Transcriptional Regulatory Network Correlates with Regulatory Divergence in Drosophila.

Transcriptional control of gene expression is regulated by biochemical interactions between cis-regulatory DNA sequences and trans-acting factors that form complex regulatory networks. Genetic changes affecting both cis- and trans-acting sequences in these networks have been shown to alter patterns of gene expression as well as higher-order organismal phenotypes. Here, we investigate how the structure of these regulatory networks relates to patterns of polymorphism and divergence in gene expression. To do this, we compared a transcriptional regulatory network inferred for Drosophila melanogaster to differences in gene regulation observed between two strains of D. melanogaster as well as between two pairs of closely related species: Drosophila sechellia and Drosophila simulans, and D. simulans and D. melanogaster. We found that the number of transcription factors predicted to directly regulate a gene ("in-degree") was negatively correlated with divergence in both gene expression (mRNA abundance) and cis-regulation. This observation suggests that the number of transcription factors directly regulating a gene's expression affects the conservation of cis-regulation and gene expression over evolutionary time. We also tested the hypothesis that transcription factors regulating more target genes (higher "out-degree") are less likely to evolve changes in their cis-regulation and expression (presumably due to increased pleiotropy), but found little support for this predicted relationship. Taken together, these data show how the architecture of regulatory networks can influence regulatory evolution.

Fitness Effects of Cis-Regulatory Variants in the Saccharomyces cerevisiae TDH3 Promoter.

Variation in gene expression is widespread within and between species, but fitness consequences of this variation are generally unknown. Here, we use mutations in the Saccharomyces cerevisiae TDH3 promoter to assess how changes in TDH3 expression affect cell growth. From these data, we predict the fitness consequences of de novo mutations and natural polymorphisms in the TDH3 promoter. Nearly all mutations and polymorphisms in the TDH3 promoter were found to have no significant effect on fitness in the environment assayed, suggesting that the wild-type allele of this promoter is robust to the effects of most new cis-regulatory mutations.

Contrasting Frequencies and Effects of cis- and trans-Regulatory Mutations Affecting Gene Expression.

Heritable differences in gene expression are caused by mutations in DNA sequences encoding cis-regulatory elements and trans-regulatory factors. These two classes of regulatory change differ in their relative contributions to expression differences in natural populations because of the combined effects of mutation and natural selection. Here, we investigate how new mutations create the regulatory variation upon which natural selection acts by quantifying the frequencies and effects of hundreds of new cis- and trans-acting mutations altering activity of the TDH3 promoter in the yeast Saccharomyces cerevisiae in the absence of natural selection. We find that cis-regulatory mutations have larger effects on expression than trans-regulatory mutations and that while trans-regulatory mutations are more common overall, cis- and trans-regulatory changes in expression are equally abundant when only the largest changes in expression are considered. In addition, we find that cis-regulatory mutations are skewed toward decreased expression while trans-regulatory mutations are skewed toward increased expression. We also measure the effects of cis- and trans-regulatory mutations on the variability in gene expression among genetically identical cells, a property of gene expression known as expression noise, finding that trans-regulatory mutations are much more likely to decrease expression noise than cis-regulatory mutations. Because new mutations are the raw material upon which natural selection acts, these differences in the frequencies and effects of cis- and trans-regulatory mutations should be considered in models of regulatory evolution.

The roles of cis- and trans-regulation in the evolution of regulatory incompatibilities and sexually dimorphic gene expression.

Evolutionary changes in gene expression underlie many aspects of phenotypic diversity within and among species. Understanding the genetic basis for evolved changes in gene expression is therefore an important component of a comprehensive understanding of the genetic basis of phenotypic evolution. Using interspecific introgression hybrids, we examined the genetic basis for divergence in genome-wide patterns of gene expression between Drosophila simulans and Drosophila mauritiana. We find that cis-regulatory and trans-regulatory divergences differ significantly in patterns of genetic architecture and evolution. The effects of cis-regulatory divergence are approximately additive in heterozygotes, quantitatively different between males and females, and well predicted by expression differences between the two parental species. In contrast, the effects of trans-regulatory divergence are associated with largely dominant introgressed alleles, have similar effects in the two sexes, and generate expression levels in hybrids outside the range of expression in both parental species. Although the effects of introgressed trans-regulatory alleles are similar in males and females, expression levels of the genes they regulate are sexually dimorphic between the parental D. simulans and D. mauritiana strains, suggesting that pure-species genotypes carry unlinked modifier alleles that increase sexual dimorphism in expression. Our results suggest that independent effects of cis-regulatory substitutions in males and females may favor their role in the evolution of sexually dimorphic phenotypes, and that trans-regulatory divergence is an important source of regulatory incompatibilities.

Evolution of splicing regulatory networks in Drosophila.

The proteome expanding effects of alternative pre-mRNA splicing have had a profound impact on eukaryotic evolution. The events that create this diversity can be placed into four major classes: exon skipping, intron retention, alternative 5' splice sites, and alternative 3' splice sites. Although the regulatory mechanisms and evolutionary pressures among alternative splicing classes clearly differ, how these differences affect the evolution of splicing regulation remains poorly characterized. We used RNA-seq to investigate splicing differences in D. simulans, D. sechellia, and three strains of D. melanogaster. Regulation of exon skipping and tandem alternative 3' splice sites (NAGNAGs) were more divergent than other splicing classes. Splicing regulation was most divergent in frame-preserving events and events in noncoding regions. We further determined the contributions of cis- and trans-acting changes in splicing regulatory networks by comparing allele-specific splicing in F1 interspecific hybrids, because differences in allele-specific splicing reflect changes in cis-regulatory element activity. We find that species-specific differences in intron retention and alternative splice site usage are primarily attributable to changes in cis-regulatory elements (median ∼80% cis), whereas species-specific exon skipping differences are driven by both cis- and trans-regulatory divergence (median ∼50% cis). These results help define the mechanisms and constraints that influence splicing regulatory evolution and show that networks regulating the four major classes of alternative splicing diverge through different genetic mechanisms. We propose a model in which differences in regulatory network architecture among classes of alternative splicing affect the evolution of splicing regulation.

Tempo and mode of regulatory evolution in Drosophila.

Genetic changes affecting gene expression contribute to phenotypic divergence; thus, understanding how regulatory networks controlling gene expression change over time is critical for understanding evolution. Prior studies of expression differences within and between species have identified properties of regulatory divergence, but technical and biological differences among these studies make it difficult to assess the generality of these properties or to understand how regulatory changes accumulate with divergence time. Here, we address these issues by comparing gene expression among strains and species of Drosophila with a range of divergence times and use F1 hybrids to examine inheritance patterns and disentangle cis- and trans-regulatory changes. We find that the fixation of compensatory changes has caused the regulation of gene expression to diverge more rapidly than gene expression itself. Specifically, we observed that the proportion of genes with evidence of cis-regulatory divergence has increased more rapidly with divergence time than the proportion of genes with evidence of expression differences. Surprisingly, the amount of expression divergence explained by cis-regulatory changes did not increase steadily with divergence time, as was previously proposed. Rather, one species (Drosophila sechellia) showed an excess of cis-regulatory divergence that we argue most likely resulted from positive selection in this lineage. Taken together, this work reveals not only the rate at which gene expression evolves, but also the molecular and evolutionary mechanisms responsible for this evolution.

Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence.

Cis-regulatory sequences, such as enhancers and promoters, control development and physiology by regulating gene expression. Mutations that affect the function of these sequences contribute to phenotypic diversity within and between species. With many case studies implicating divergent cis-regulatory activity in phenotypic evolution, researchers have recently begun to elucidate the genetic and molecular mechanisms that are responsible for cis-regulatory divergence. Approaches include detailed functional analysis of individual cis-regulatory elements and comparing mechanisms of gene regulation among species using the latest genomic tools. Despite the limited number of mechanistic studies published to date, this work shows how cis-regulatory activity can diverge and how studies of cis-regulatory divergence can address long-standing questions about the genetic mechanisms of phenotypic evolution.

Molecular Mechanisms and Evolutionary Processes Contributing to Accelerated Divergence of Gene Expression on the Drosophila X Chromosome.

In species with a heterogametic sex, population genetics theory predicts that DNA sequences on the X chromosome can evolve faster than comparable sequences on autosomes. Both neutral and nonneutral evolutionary processes can generate this pattern. Complex traits like gene expression are not predicted to have accelerated evolution by these theories, yet a "faster-X" pattern of gene expression divergence has recently been reported for both Drosophila and mammals. Here, we test the hypothesis that accelerated adaptive evolution of cis-regulatory sequences on the X chromosome is responsible for this pattern by comparing the relative contributions of cis- and trans-regulatory changes to patterns of faster-X expression divergence observed between strains and species of Drosophila with a range of divergence times. We find support for this hypothesis, especially among male-biased genes, when comparing different species. However, we also find evidence that trans-regulatory differences contribute to a faster-X pattern of expression divergence both within and between species. This contribution is surprising because trans-acting regulators of X-linked genes are generally assumed to be randomly distributed throughout the genome. We found, however, that X-linked transcription factors appear to preferentially regulate expression of X-linked genes, providing a potential mechanistic explanation for this result. The contribution of trans-regulatory variation to faster-X expression divergence was larger within than between species, suggesting that it is more likely to result from neutral processes than positive selection. These data show how accelerated evolution of both coding and noncoding sequences on the X chromosome can lead to accelerated expression divergence on the X chromosome relative to autosomes.

The significance and scope of evolutionary developmental biology: a vision for the 21st century.

Evolutionary developmental biology (evo-devo) has undergone dramatic transformations since its emergence as a distinct discipline. This paper aims to highlight the scope, power, and future promise of evo-devo to transform and unify diverse aspects of biology. We articulate key questions at the core of eleven biological disciplines-from Evolution, Development, Paleontology, and Neurobiology to Cellular and Molecular Biology, Quantitative Genetics, Human Diseases, Ecology, Agriculture and Science Education, and lastly, Evolutionary Developmental Biology itself-and discuss why evo-devo is uniquely situated to substantially improve our ability to find meaningful answers to these fundamental questions. We posit that the tools, concepts, and ways of thinking developed by evo-devo have profound potential to advance, integrate, and unify biological sciences as well as inform policy decisions and illuminate science education. We look to the next generation of evolutionary developmental biologists to help shape this process as we confront the scientific challenges of the 21st century.

Contrasting properties of gene-specific regulatory, coding, and copy number mutations in Saccharomyces cerevisiae: frequency, effects, and dominance.

Genetic variation within and between species can be shaped by population-level processes and mutation; however, the relative impact of "survival of the fittest" and "arrival of the fittest" on phenotypic evolution remains unclear. Assessing the influence of mutation on evolution requires understanding the relative rates of different types of mutations and their genetic properties, yet little is known about the functional consequences of new mutations. Here, we examine the spectrum of mutations affecting a focal gene in Saccharomyces cerevisiae by characterizing 231 novel haploid genotypes with altered activity of a fluorescent reporter gene. 7% of these genotypes had a nonsynonymous mutation in the coding sequence for the fluorescent protein and were classified as "coding" mutants; 2% had a change in the S. cerevisiae TDH3 promoter sequence controlling expression of the fluorescent protein and were classified as "cis-regulatory" mutants; 10% contained two copies of the reporter gene and were classified as "copy number" mutants; and the remaining 81% showed altered fluorescence without a change in the reporter gene itself and were classified as "trans-acting" mutants. As a group, coding mutants had the strongest effect on reporter gene activity and always decreased it. By contrast, 50%-95% of the mutants in each of the other three classes increased gene activity, with mutants affecting copy number and cis-regulatory sequences having larger median effects on gene activity than trans-acting mutants. When made heterozygous in diploid cells, coding, cis-regulatory, and copy number mutant genotypes all had significant effects on gene activity, whereas 88% of the trans-acting mutants appeared to be recessive. These differences in the frequency, effects, and dominance among functional classes of mutations might help explain why some types of mutations are found to be segregating within or fixed between species more often than others.