Search

Heritable variation in gene expression is common within and among species and contributes to phenotypic diversity. Mutations affecting either cis - or trans -regulatory sequences controlling gene expression give rise to variation in gene expression, and natural selection acting on this variation causes some regulatory variants to persist in a population for longer than others. To understand how mutation and selection interact to produce the patterns of regulatory variation we see within and among species, my colleagues and I have been systematically determining the effects of new mutations on expression of the TDH3 gene in Saccharomyces cerevisiae and comparing them to the effects of polymorphisms segregating within this species. We have also investigated the molecular mechanisms by which regulatory variants act. Over the past decade, this work has revealed properties of cis - and trans -regulatory mutations including their relative frequency, effects, dominance, pleiotropy and fitness consequences. Comparing these mutational effects to the effects of polymorphisms in natural populations, we have inferred selection acting on expression level, expression noise and phenotypic plasticity. Here, I summarize this body of work and synthesize its findings to make inferences not readily discernible from the individual studies alone. This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

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 theSaccharomyces cerevisiaegeneTDH3by its paralogsTDH1andTDH2. TDH1andTDH2are upregulated in a dose-dependent manner in response to reductions inTDH3by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar toTDH2, suggesting that the active compensation byTDH3paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

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.

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.

Variation in gene expression arises from cis- and trans-regulatory mutations, which contribute differentially to expression divergence. We compare the impacts on gene expression and fitness resulting from cis- and trans-regulatory mutations in Saccharomyces cerevisiae , with a focus on the TDH3 gene. We use the effects of cis-regulatory mutations to infer effects of trans-regulatory mutations attributable to impacts beyond the focal gene, revealing a distribution of pleiotropic effects. Cis- and trans-regulatory mutations had different effects on gene expression with pleiotropic effects of trans-regulatory mutants affecting expression of genes both in parallel to and downstream of the focal gene. The more widespread and deleterious effects of trans-regulatory mutations we observed are consistent with their decreasing relative contribution to expression differences over evolutionary time.

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Heritable variation in gene expression is common within species. Much of this variation is due to genetic differences outside of the gene with altered expression and is trans-acting. This trans-regulatory variation is often polygenic, with individual variants typically having small effects, making the genetic architecture and evolution of trans-regulatory variation challenging to study. Consequently, key questions about trans-regulatory variation remain, including the variability of trans-regulatory variation within a species, how selection affects trans-regulatory variation, and how trans-regulatory variants are distributed throughout the genome and within a species. To address these questions, we isolated and measured trans-regulatory differences affecting TDH3 promoter activity among 56 strains of Saccharomyces cerevisiae, finding that trans-regulatory backgrounds varied approximately twofold in their effects on TDH3 promoter activity. Comparing this variation to neutral models of trans-regulatory evolution based on empirical measures of mutational effects revealed that despite this variability in the effects of trans-regulatory backgrounds, stabilizing selection has constrained trans-regulatory differences within this species. Using a powerful quantitative trait locus mapping method, we identified ∼100 trans-acting expression quantitative trait locus in each of three crosses to a common reference strain, indicating that regulatory variation is more polygenic than previous studies have suggested. Loci altering expression were located throughout the genome, and many loci were strain specific. This distribution and prevalence of alleles is consistent with recent theories about the genetic architecture of complex traits. In all mapping experiments, the nonreference strain alleles increased and decreased TDH3 promoter activity with similar frequencies, suggesting that stabilizing selection maintained many trans-acting variants with opposing effects. This variation may provide the raw material for compensatory evolution and larger scale regulatory rewiring observed in developmental systems drift among species.

Evolutionary changes in enhancers, which control gene expression, often contribute to phenotypic evolution. Here, Kalay et al. examine how enhancer activities are encoded within non-coding sequences surrounding the pigmentation gene yellow from three Drosophila species... Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5′ intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni. We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Heritable variation in a gene’s expression arises from mutations impactingcis- andtrans-acting components of its regulatory network, with expression variation often derived fromtrans-regulatory mutations within species. Here, we investigate howtrans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations withtrans-regulatory effects on expression of the same focal gene inSaccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that thesetrans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of thetrans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism and signal transduction. Finally, we find that the genomic distribution of these new regulatory mutations significantly overlaps with the genomic distribution oftrans-regulatory variants segregating withinS. cerevisiae.

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. In this Review, Hill et al. discuss how high-throughput methods for creating and characterizing mutations are providing insight into how regulatory variation is generated and evolves.

Understanding the genetic basis of species differences is a major goal in evolutionary biology. Pigmentation divergence betweenDrosophilaspecies often involves genetic changes in pigmentation candidate genes that pattern the body and wings, but it remains unclear how these changes affect pigmentation evolution in multiple body parts between the same diverging species.Drosophila elegansandD. gunungcolashow pigmentation differences in the thorax, legs, and wings, withD. elegansexhibiting male-specific wing spots andD. gunungcolalacking wing spots with intensely dark thoraces and legs. Here, we performed QTL mapping to identify the genetic architecture of these differences. We find a large effect QTL on the X chromosome for all three body parts. QTL on Muller Element E were found for thorax pigmentation in both backcrosses but were only marginally significant in one backcross for the legs and wings. Consistent with this observation, we isolated the effects of the Muller Element E QTL by introgressingD. gunungcolaalleles into aD. elegansgenetic background and found thatD. gunungcolaalleles linked near the pigmentation candidate geneebonycaused intense darkening of the thorax, minimal darkening of legs, and minimal shrinking of wing spots.D. elegans ebonymutants showed changes in pigmentation consistent with Ebony having different effects on pigmentation in different tissues. Our results suggest that multiple genes have evolved differential effects on pigmentation levels in different body regions.

Phenotypic variation within a species is often structured geographically in clines. InDrosophila americana, a longitudinal cline for body color exists within North America that appears to be due to local adaptation. Thetanandebonygenes have been hypothesized to contribute to this cline, with alleles of both genes that lighten body color found inD. americana. These alleles are similar in sequence and function to the allele fixed inD. americana’smore lightly pigmented sister species,Drosophila novamexicana. To test this hypothesis, we examined the frequency and geographic distribution ofD. novamexicana-like alleles oftanandebonyinD. americana. Among alleles from over 100 strains ofD. americanaisolated from 21 geographic locations, we failed to identify additional alleles oftanorebonywith as much sequence similarity toD. novamexicanaas the alleles previously described. However, using genetic analysis of 51D. americanastrains derived from 20 geographic locations, we identified one new allele ofebonyand one new allele oftansegregating inD. americanathat are functionally equivalent to theD. novamexicanaallele. An additional 5 alleles oftanalso 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 analyzed indicate a more complex genetic architecture, with diversity in the number and effects of loci altering pigmentation observed both within and among populations ofD. americana. This genetic heterogeneity poses a challenge to association studies and genomic scans for clinal variation, but might be common in natural populations.

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.

Variation in gene expression is widespread within and between species, but fitness consequences of this variation are generally unknown. Here we use mutations in theSaccharomyces cerevisiae TDH3promoter to assess how changes inTDH3expression affect cell growth. From these data, we predict the fitness consequences ofde novomutations and natural polymorphisms in theTDH3promoter. Nearly all mutations and polymorphisms in theTDH3promoter 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 newcis-regulatory mutations.

Heritable changes in gene expression are important contributors to phenotypic differences within and between species and are caused by mutations in cis-regulatory elements and trans-regulatory factors. Although previous work has suggested that cis-regulatory differences preferentially accumulate with time, technical restrictions to closely related species and limited comparisons have made this observation difficult to test. To address this problem, we used allele-specific RNA-seq data from Saccharomyces species and hybrids to expand both the evolutionary timescale and number of species in which the evolution of regulatory divergence has been investigated. We find that as sequence divergence increases, cis-regulatory differences do indeed become the dominant type of regulatory difference between species, ultimately becoming a better predictor of expression divergence than trans-regulatory divergence. When both cis- and trans-regulatory differences accumulate for the same gene, they more often have effects in opposite directions than in the same direction, indicating widespread compensatory changes underlying the evolution of gene expression. The frequency of compensatory changes within and between species and the magnitude of effect for the underlying cis- and trans-regulatory differences suggests that compensatory changes accumulate primarily due to selection against divergence in gene expression as a result of weak stabilizing selection on gene expression levels. These results show that cis-regulatory differences and compensatory changes in regulation play increasingly important roles in the evolution of gene expression as time increases.

Drosophila sechellia is a species of fruit fly endemic to the Seychelles islands. Unlike its generalist sister species, D. sechellia has evolved to be a specialist on the host plant Morinda citrifolia. This specialization is interesting because the plant’s fruit contains secondary defense compounds, primarily octanoic acid (OA), that are lethal to most other Drosophilids. Although ecological and behavioral adaptations to this toxic fruit are known, the genetic basis for evolutionary changes in OA resistance are not. Prior work showed that a genomic region on chromosome 3R containing 18 genes has the greatest contribution to differences in OA resistance between D. sechellia and D. simulans. To determine which gene(s) in this region might be involved in the evolutionary change in OA resistance, we knocked-down expression of each gene in this region in D. melanogaster with RNA interference (RNAi) (i) ubiquitously throughout development, (ii) during only the adult stage, and (iii) within specific tissues. We identified three neighboring genes in the Osiris family, Osiris 6 (Osi6), Osi7, and Osi8, that lead to decreased OA resistance when ubiquitously knocked-down. Tissue specific RNAi, however, showed that decreasing expression of Osi6 and Osi7 specifically in the fat body and/or salivary glands increased OA resistance. Gene expression analyses of Osi6 and Osi7 revealed that while standing levels of expression are higher in D. sechellia, Osi6 expression is significantly downregulated in salivary glands in response to OA exposure, suggesting that evolved tissue-specific environmental plasticity of Osi6 expression may be responsible for OA resistance in D. sechellia.

The evolution of sexual traits often involves correlated changes in morphology and behavior. For example, in Drosophila, divergent mating displays are often accompanied by divergent pigment patterns. To better understand how such traits co-evolve, we investigated the genetic basis of correlated divergence in wing pigmentation and mating display between the sibling species Drosophila elegans and Drosophila gunungcola. Drosophila elegans males have an area of black pigment on their wings known as a wing spot and appear to display this spot to females by extending their wings laterally during courtship. By contrast, D. gunungcola lost both of these traits. Using Multiplexed Shotgun Genotyping (MSG), we identified a ∼440 kb region on the X chromosome that behaves like a genetic switch controlling the presence or absence of male-specific wing spots. This region includes the candidate gene optomotor-blind (omb), which plays a critical role in patterning the Drosophila wing. The genetic basis of divergent wing display is more complex, with at least two loci on the X chromosome and two loci on autosomes contributing to its evolution. Introgressing the X-linked region affecting wing spot development from D. gunungcola into D. elegans reduced pigmentation in the wing spots but did not affect the wing display, indicating that these are genetically separable traits. Consistent with this observation, broader sampling of wild D. gunungcola populations confirmed that the wing spot and wing display are evolving independently: some D. gunungcola males performed wing displays similar to D. elegans despite lacking wing spots. These data suggest that correlated selection pressures rather than physical linkage or pleiotropy are responsible for the coevolution of these morphological and behavioral traits. They also suggest that the change in morphology evolved prior to the change in behavior.

The evolution of sexual traits often involves correlated changes in morphology and behavior. For example, in Drosophila, divergent mating displays are often accompanied by divergent pigment patterns. To better understand how such traits co-evolve, we investigated the genetic basis of correlated divergence in wing pigmentation and mating display between the sibling speciesDrosophila elegansandD. gunungcola.Drosophila elegansmales have an area of black pigment on their wings known as a wing spot and appear to display this spot to females by extending their wings laterally during courtship. By contrast,D. gunungcolalacks both of these traits. Using Multiplexed Shotgun Genotyping (MSG), we identified a ∼440 kb region on the X chromosome that behaves like a genetic switch controlling the presence or absence of male-specific wing spots. This region includes the candidate geneoptomotor-blind(omb), which plays a critical role in patterning theDrosophilawing. The genetic basis of divergent wing display is more complex, with at least two loci on the X chromosome and two loci on autosomes contributing to its evolution. Introgressing the X-linked region affecting wing spot development fromD. gunungcolaintoD. elegansreduced pigmentation in the wing spots but did not affect the wing display, indicating that these are genetically separable traits. Consistent with this observation, broader sampling of wildD. gunungcolapopulations confirmed the wing spot and wing display are evolving independently: someD. gunungcolamales preformed wing displays similar toD. elegansdespite lacking wing spots. These data suggest that correlated selection pressures rather than physical linkage or pleiotropy are responsible for the coevolution of these morphological and behavioral traits. They also suggest that the change in morphology evolved prior to the change in behavior.

Drosophila melanogastermales perform a series of courtship behaviors that, when successful, result in copulation with a female. For over a century, mutations in theyellowgene, named for its effects on pigmentation, have been known to reduce male mating success. Prior work has suggested thatyellowinfluences mating behavior through effects on wing extension, song, and/or courtship vigor. Here, we rule out these explanations, as well as effects on the nervous system more generally, and find instead that the effects ofyellowon male mating success are mediated by its effects on pigmentation of male-specific leg structures called sex combs. Loss ofyellowexpression in these modified bristles reduces their melanization, which changes their structure and causes difficulty grasping females prior to copulation. These data illustrate why the mechanical properties of anatomy, and not just neural circuitry, must be considered to fully understand the development and evolution of behavior.

Heritable variation in gene expression is common within species. Much of this variation is due to genetic changes at loci other than the affected gene and is thus trans-acting. This trans-regulatory variation is often polygenic, with individual variants typically having small effects, making the genetic architecture of trans-regulatory variation challenging to study. Consequently, key questions about trans-regulatory variation remain, including how selection affects this variation and how trans-regulatory variants are distributed throughout the genome and within species. Here, we show that trans-regulatory variation affecting TDH3 promoter activity is common among strains of Saccharomyces cerevisiae. Comparing this variation to neutral models of trans-regulatory evolution based on empirical measures of mutational effects revealed that stabilizing selection has constrained this variation. Using a powerful quantitative trait locus (QTL) mapping method, we identified ∼100 loci altering expression between a reference strain and each of three genetically distinct strains. In all three cases, the non-reference strain alleles increased and decreased TDH3 promoter activity with similar frequencies, suggesting that stabilizing selection maintained many trans-acting variants with opposing effects. Loci altering expression were located throughout the genome, with many loci being strain specific and others being shared among multiple strains. These findings are consistent with theory showing stabilizing selection for quantitative traits can maintain many alleles with opposing effects, and the wide-spread distribution of QTL throughout the genome is consistent with the omnigenic model of complex trait variation. Furthermore, the prevalence of alleles with opposing effects might provide raw material for compensatory evolution and developmental systems drift.Significance statementGene expression varies among individuals in a population due to genetic differences in regulatory components. To determine how this variation is distributed within genomes and species, we used a powerful genetic mapping approach to examine multiple strains of Saccharomyces cerevisiae. Despite evidence of stabilizing selection maintaining gene expression levels among strains, we find hundreds of loci that affect expression of a single gene. These loci vary among strains and include similar frequencies of alleles that increase and decrease expression. As a result, each strain contains a unique set of compensatory alleles that lead to similar levels of gene expression among strains. This regulatory variation might form the basis for large scale regulatory rewiring observed between distantly related species.

Pleiotropic genes are genes that affect more than one trait. For example, many genes required for pigmentation in the fruit flyDrosophila melanogasteralso affect traits such as circadian rhythms, vision, and mating behavior. Here, we present evidence that two pigmentation genes,ebonyandtan, which encode enzymes catalyzing reciprocal reactions in the melanin biosynthesis pathway, also affect cuticular hydrocarbon (CHC) composition inD. melanogasterfemales. More specifically, we report thatebonyloss-of-function mutants have a CHC profile that is biased toward long (>25C) chain CHCs, whereastanloss-of-function mutants have a CHC profile that is biased toward short (ebonymutants, making the CHC profiles similar to those seen intanmutants. These observations suggest that genetic variation affectingebonyand/ortanactivity might cause correlated changes in pigmentation and CHC composition in natural populations. We tested this possibility using theDrosophilaGenetic Reference Panel (DGRP) and found that CHC composition covaried with pigmentation as well as levels ofebonyandtanexpression in newly eclosed adults in a manner consistent with theebonyandtanmutant phenotypes. These data suggest that the pleiotropic effects ofebonyandtanmight contribute to covariation of pigmentation and CHC profiles inDrosophila.

The regulation of gene expression controls development, and changes in this regulation often contribute to phenotypic evolution. Drosophila pigmentation is a model system for studying evolutionary changes in gene regulation, with differences in expression of pigmentation genes such as yellow that correlate with divergent pigment patterns among species shown to be caused by changes in cis- and trans-regulation. Currently, much more is known about the cis-regulatory component of divergent yellow expression than the trans-regulatory component, in part because very few trans-acting regulators of yellow expression have been identified. This study aims to improve our understanding of the trans-acting control of yellow expression by combining yeast-one-hybrid and RNAi screens for transcription factors binding to yellow cis-regulatory sequences and affecting abdominal pigmentation in adults, respectively. Of the 670 transcription factors included in the yeast-one-hybrid screen, 45 showed evidence of binding to one or more sequence fragments tested from the 5′ intergenic and intronic yellow sequences from D. melanogaster, D. pseudoobscura, and D. willistoni, suggesting that they might be direct regulators of yellow expression. Of the 670 transcription factors included in the yeast-one-hybrid screen, plus another TF previously shown to be genetically upstream of yellow, 125 were also tested using RNAi, and 32 showed altered abdominal pigmentation. Nine transcription factors were identified in both screens, including four nuclear receptors related to ecdysone signaling (Hr78, Hr38, Hr46, and Eip78C). This finding suggests that yellow expression might be directly controlled by nuclear receptors influenced by ecdysone during early pupal development when adult pigmentation is forming.

Gene expression noise is an evolvable property of biological systems that describes differences in expression among genetically identical cells in the same environment. Prior work has shown that expression noise is heritable and can be shaped by selection, but the impact of variation in expression noise on organismal fitness has proven difficult to measure. Here, we quantify the fitness effects of altering expression noise for the TDH3 gene in Saccharomyces cerevisiae. We show that increases in expression noise can be deleterious or beneficial depending on the difference between the average expression level of a genotype and the expression level maximizing fitness. We also show that a simple model relating single-cell expression levels to population growth produces patterns consistent with our empirical data. We use this model to explore a broad range of average expression levels and expression noise, providing additional insight into the fitness effects of variation in expression noise., eLife digest Single-celled organisms that reproduce by dividing, like yeast, can create whole populations of genetically identical cells. However, some differences will exist among such cells, even when they have all experienced the same environment. These differences are known as “noise”. By definition, noise is not caused by differences in DNA sequence, but some DNA sequences are noisier than others (i.e. they cause more differences among cells). Because the amount of noise can be under genetic control, noise could evolve due to natural selection. Scientists often study noise at the level of gene expression – in other words, how many RNA or protein molecules are produced from each gene within each cell. Prior work has suggested that this type of noise can affect how often individual cells divide in a population, which is a component of that population’s fitness. Yet directly measuring these effects has proven challenging. Different studies have in the past reached opposite conclusions about whether a change in gene expression noise would increase or decrease fitness. One major reason for the lack of clear results is that most mutations that alter gene expression noise also alter the average level of expression of that gene. To find DNA sequences that produced the same average amount of protein but different levels of expression noise, Duveau et al. compared the effects of hundreds of mutations in the DNA sequence regulating the expression of a gene in baker’s yeast. Experiments focused on 43 DNA sequences then showed that increased expression noise could either speed up or slow down the growth of the population by affecting how long it took each cell to divide. More specifically, the effects of increasing expression noise depended on the average amount of protein produced among the cells in the population. If the average expression level was close to the optimum amount at which cells divided as fast as possible, increasing expression noise reduced the growth of the whole population. If, however, the average protein level caused cells to divide slower than their maximum rate, increasing expression noise resulted in faster growth of the population as a whole. Duveau et al. explain their results as follows: more expression noise in a population that is already making the optimal amount of protein can reduce fitness because it increases the fraction of that population making a suboptimal amount of the protein. However, when the average expression level is not optimal, more expression noise would mean more cells producing an amount of protein that is closer to the optimum and thus having higher fitness. These findings provide conceptual tools needed to understand how genetic variation affecting expression noise evolves. They could also help understand how expression noise might contribute to biological processes that depend upon cell division, such as diseases like cancer.

Gene expression noise is an evolvable property of biological systems that describes differences in gene expression among genetically identical cells in the same environment. Prior work has shown that expression noise is heritable and can be shaped by natural selection, but the impact of variation in expression noise on organismal fitness has proven difficult to measure. Here, we quantify the fitness effects of altering expression noise for theTDH3gene inSaccharomyces cerevisiae. We show that increases in expression noise can be deleterious or beneficial depending on the difference between the average expression level of a genotype and the expression level maximizing fitness. We also show that a simple model relating single-cell expression levels to population growth produces patterns that are consistent with our empirical data. We use this model to explore a broad range of average expression levels and expression noise, providing additional insight into the fitness effects of variation in expression noise.

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.

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.

From seasonal variation in the color of butterfly wings to trees bending toward the light, organisms often change in response to their environment. These changes, known as phenotypic plasticity, can result from differences in how genes are expressed among environments. Mutations causing environment-specific changes in gene expression provide raw material for phenotypic plasticity, but their frequency, effect size, and direction of effects among environments are not well understood. This study shows that mutations in the promoter of a yeast metabolic gene often display environment-dependent effects on gene expression and that these environment-dependent effects have been shaped by selection in natural populations.

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.

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 theSaccharomyces cerevisiae TDH3promoter (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 to the effects of 30PTDH3polymorphisms on expression plasticity in the same environments provided evidence of natural selection acting to prevent the plastic response inPTDH3activity 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, however, 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.Significance StatementFrom seasonal variation in the color of butterfly wings to trees bending toward the light, organisms often change in response to their environment. These changes, known as phenotypic plasticity, can result from differences in how genes are expressed among environments. Mutations causing environment-specific changes in gene expression provide raw material for phenotypic plasticity, but their frequency, effect size and direction of effects among environments are not well understood. This study shows that mutations in the promoter of a yeast metabolic gene often display environment-dependent effects on gene expression and that these environment-dependent effects have been shaped by selection in natural populations.

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.

Genetic variants identified by mapping are biased toward large phenotypic effects because of methodologic challenges for detecting genetic variants with small phenotypic effects. Recently, bulk segregant analysis combined with next-generation sequencing (BSA-seq) was shown to be a powerful and cost-effective way to map small effect variants in natural populations. Here, we examine the power of BSA-seq for efficiently mapping small effect mutations isolated from a mutagenesis screen. Specifically, we determined the impact of segregant population size, intensity of phenotypic selection to collect segregants, number of mitotic generations between meiosis and sequencing, and average sequencing depth on power for mapping mutations with a range of effects on the phenotypic mean and standard deviation as well as relative fitness. We then used BSA-seq to map the mutations responsible for three ethyl methanesulfonate−induced mutant phenotypes in Saccharomyces cerevisiae. These mutants display small quantitative variation in the mean expression of a fluorescent reporter gene (−3%, +7%, and +10%). Using a genetic background with increased meiosis rate, a reliable mating type marker, and fluorescence-activated cell sorting to efficiently score large segregating populations and isolate cells with extreme phenotypes, we successfully mapped and functionally confirmed a single point mutation responsible for the mutant phenotype in all three cases. Our simulations and experimental data show that the effects of a causative site not only on the mean phenotype, but also on its standard deviation and relative fitness should be considered when mapping genetic variants in microorganisms such as yeast that require population growth steps for BSA-seq.

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.

One of the central questions in evolutionary genetics is how much of the genome is involved in the early stages of divergence between populations, causing them to be reproductively isolated. In this article, we investigate genomic differentiation in a pair of closely related field crickets (Gryllus firmus and G. pennsylvanicus). These two species are the result of allopatric divergence and now interact along an extensive hybrid zone in eastern North America. Genes encoding seminal fluid proteins (SFPs) are often divergent between species, and it has been hypothesized that these proteins may play a key role in the origin and maintenance of reproductive isolation between diverging lineages. Hence, we chose to scan the accessory gland transcriptome to enable direct comparisons of differentiation for genes known to encode SFPs with differentiation in a much larger set of genes expressed in the same tissue. We have characterized differences in allele frequency between two populations for >6000 SNPs and >26,000 contigs. About 10% of all SNPs showed nearly fixed differences between the two species. Genes encoding SFPs did not have significantly elevated numbers of fixed SNPs per contig, nor did they seem to show larger differences than expected in their average allele frequencies. The distribution of allele frequency differences across the transcriptome is distinctly bimodal, but the relatively high proportion of fixed SNPs does not necessarily imply “ancient” divergence between these two lineages. Further studies of linkage disequilibrium and introgression across the hybrid zone are needed to direct our attention to those genome regions that are important for reproductive isolation.

Pigmentation is a model trait for evolutionary and developmental analysis that is particularly amenable to molecular investigation in the genus Drosophila. To better understand how this phenotype evolves, we examined divergent pigmentation and gene expression over developmental time in the dark-bodied D. americana and its light-bodied sister species D. novamexicana. Prior genetic analysis implicated two enzyme-encoding genes, tan and ebony, in pigmentation divergence between these species, but questions remain about the underlying molecular mechanisms. Here, we describe stages of pupal development in both species and use this staging to determine when pigmentation develops and diverges between D. americana and D. novamexicana. For the developmental stages encompassing pigment divergence, we compare mRNA expression of tan and ebony over time and between species. Finally, we use allele-specific expression assays to determine whether interspecific differences in mRNA abundance have a cis-regulatory basis and find evidence of cis-regulatory divergence for both tan and ebony. cis-regulatory divergence affecting tan had a small effect on mRNA abundance and was limited to a few developmental stages, yet previous data suggests that this divergence is likely to be biologically meaningful. Our study suggests that small and developmentally transient expression changes may contribute to phenotypic diversification more often than commonly appreciated. Recognizing the potential phenotypic impact of such changes is important for a scientific community increasingly focused on dissecting quantitative variation, but detecting these types of changes will be a major challenge to elucidating the molecular basis of complex traits.