Regulatory Evolution: Yeast rises to the occasion

  1. Mark A Ragan  Is a corresponding author
  1. University of Queensland, Australia

One of the less-anticipated outcomes from the past decade of genomics is how poorly the complexity of organisms correlates with their gene number. Even allowing for a large measure of anthropocentrism, it remains puzzling that humans have only about 20,800 protein-coding genes, whereas water fleas have about 30,900 and the rice plant has 40,000 or so. Working with this pedestrian number of genes, the human body nonetheless gives rise to more than 400 specialised cell types (Vickaryous and Hall, 2006), and the remarkable diversity of these cell types in terms of form and function arises from different sets of genes (known as modules) being co-expressed at different times.

We understand quite a lot about how genes themselves evolve and diversify, but we know very much less about the evolution of the processes that regulate the expression of genes. Are modules stable over evolutionary timescales, or are they assembled opportunistically as required? Are duplicate copies of genes retained in the ancestral module or are they reassigned to another module? And if they are reassigned, do they tend to be reassigned to the same module or to different modules, and does this happen shortly after duplication or does it continue over a much longer time? Does it matter whether the duplicates were generated sporadically or via whole-genome duplication? Is regulatory evolution driven by natural selection, and does it correlate with changes in lifestyle or the copy number of chromosomes? Indeed, are there general principles of modular gene regulation, or is the story of gene regulation one of contingency and anecdote?

Writing in eLife, Dawn Thompson, Aviv Regev and co-workers—including Thompson and Sushmita Roy as joint first authors—report the results of a large-scale assault on these questions. Thompson et al. focused on 15 species of yeast for which complete genome sequences are known (Figure 1). Yeasts offer many advantages for studying the evolution of gene regulation, over and above their well known genetics. In particular, their evolutionary history over the past 300 million years is well known and is comparable to those of teleost fishes (Near et al., 2012) and seed plants (Clarke et al., 2012) in terms of its depth. Their physiology is also well understood, but it is also surprisingly diverse: for example, different species of yeast colonise different ecological niches, utilise a range of different carbon sources, and differ in their preference for oxidative phosphorylation vs a more fermentative lifestyle.

Given gene expression profiles for a number of species, and a gene tree and a species tree for these species, an algorithm called Arboretum (Roy et al., 2013) can be used to determine how sets of genes called modules have evolved over the period covered by the these trees.

(A) Schematic diagram showing gene expression profiles (bottom) for two gene modules in three different species (X, Y and Z); the left module contains 5 or 6 genes in these species, while the right module contains 2, 3 or 4 genes. A and B are the inferred ancestral states of these modules. A local gene duplication event along the lineage AB results in duplicates (paralogs) of the red gene in the left module: one or both of these duplicates can be retained or assigned to a different module; both are retained in the module along BX, and one is reassigned along B → Y. Genes can also be lost (AZ) or gained by modules. (B) Thompson et al. studied 15 species of yeast; 13 of these are descended from a single Saccharomycotina species, and a whole genome duplication event (yellow star) resulted in seven of the species. See Figure 1A of Thompson et al. for full details.

Moreover, as a bonus, the 15-yeast dataset compiled by Thompson et al. spans a whole-genome duplication event that has affected seven of the 15 species (Figure 1; Wolfe and Shields, 1997; Kellis et al., 2004). By simultaneously duplicating all genes and their regulatory elements, even the most ancestral elements, whole-genome duplication events make it possible for a lineage to explore modes of gene regulation that would not become accessible as a result of sporadic, localised duplication events (Lynch and Katju, 2004).

Patterns of gene co-expression are sometimes conserved over substantial timescales, despite significant turnover in the associated transcription factors and chromatin organisation. Sometimes this turnover is coupled to adaptive changes in lifestyle, whereas other changes in regulation may be neutral, analogous to the genetic drift that happens at the sequence level (Tsankov et al., 2011; Baker et al., 2012). However, the intrinsic technical difficulty of these experiments, coupled with physiological diversity displayed by different species of yeast, has until now made it difficult to test the generality of these findings.

Focusing on growth in glucose and its depletion in batch culture, Thompson, Regev and co-workers—who are based at the Broad Institute of MIT and Harvard, and also at MIT—began by devising a medium that supports the growth of all 15 yeasts at comparable rates (Thompson et al., 2013). They next identified six physiologically comparable time-points along the growth curve of each yeast. Gene expression profiles confirmed that these six time-points were indeed physiologically comparable. Thompson et al. then used oligomer arrays to profile the transcriptome of each yeast species at each time point.

To compare expression profiles and track the assignment of each gene to one or another module along the phylogenetic tree, the Broad-MIT team developed a probabilistic algorithm (called Arboretum) that delineates the modules (which can change size and composition over time) and computes the trajectory of every gene through a module in each extant and inferred ancestral species (see Figure 1A; Roy et al., 2013). Thompson et al. used Arboretum to map the evolution of functional annotation, cis-regulatory motifs and nucleosome-free regions within and across modules for single-copy genes, and also for duplicates arising from sporadic or whole-genome duplication.

Many key questions yield to this systematic approach. Two-thirds of the variation in transcriptional response is captured by five expression modules. Genes are conserved within these modules in a way that is inversely proportional to evolutionary time, with two modules (those related to growth and stress-response functions) being more conservative than the other three. Gene reassignment between modules is often consistent with changes in lifestyle or the copy number of chromosomes. Duplicates are more likely to be reassigned than single-copy genes, with this reassignment often occurring in a brief ‘window of opportunity’ after duplication; however, duplicates that arise from the whole-genome duplication continue to be reassigned over a much longer time. Neo-functionalization (where one copy is retained, the other reassigned) and symmetric divergence (both copies reassigned to the same module) are more frequent than asymmetric divergence (reassignment to different modules).

Do any of these observations look like a principle? They all seem to apply regardless of gene function, lifestyle or evolutionary distance on the yeast tree. In an experiment with eight of these species, many of the same responses were seen in response to heat shock. Finally, Thompson, Regev and co-workers point to fascinating similarities between the up-regulation of genes for nucleotide salvage and glycine synthesis at one of the six time-points they study (when the available glucose has been depleted) in some yeasts, and regulatory changes involving the same pathways in cancer cells that show an effect called the Warburg effect.

Yeast is one of the less-complex eukaryotes, but the emergence of evolutionary principles for gene regulation in these experiments represents another addition to the list of unanticipated outcomes of genomic biology.


Article and author information

Author details

  1. Mark A Ragan

    Institute for Molecular Bioscience and the School of Information Technology and Electrical Engineering, University of Queensland, Brisbane, Australia
    For correspondence
    Competing interests
    The author declares that no competing interests exist.

Publication history

  1. Version of Record published: June 18, 2013 (version 1)


© 2013, Ragan

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.


  • 344
  • 23
  • 0

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Mark A Ragan
Regulatory Evolution: Yeast rises to the occasion
eLife 2:e00933.
  1. Further reading

Further reading

    1. Chromosomes and Gene Expression
    Rupam Choudhury, Anuroop Venkateswaran Venkatasubramani ... Axel Imhof
    Research Article Updated

    Eukaryotic chromatin is organized into functional domains, that are characterized by distinct proteomic compositions and specific nuclear positions. In contrast to cellular organelles surrounded by lipid membranes, the composition of distinct chromatin domains is rather ill described and highly dynamic. To gain molecular insight into these domains and explore their composition, we developed an antibody-based proximity biotinylation method targeting the RNA and proteins constituents. The method that we termed antibody-mediated proximity labelling coupled to mass spectrometry (AMPL-MS) does not require the expression of fusion proteins and therefore constitutes a versatile and very sensitive method to characterize the composition of chromatin domains based on specific signature proteins or histone modifications. To demonstrate the utility of our approach we used AMPL-MS to characterize the molecular features of the chromocenter as well as the chromosome territory containing the hyperactive X chromosome in Drosophila. This analysis identified a number of known RNA-binding proteins in proximity of the hyperactive X and the centromere, supporting the accuracy of our method. In addition, it enabled us to characterize the role of RNA in the formation of these nuclear bodies. Furthermore, our method identified a new set of RNA molecules associated with the Drosophila centromere. Characterization of these novel molecules suggested the formation of R-loops in centromeres, which we validated using a novel probe for R-loops in Drosophila. Taken together, AMPL-MS improves the selectivity and specificity of proximity ligation allowing for novel discoveries of weak protein–RNA interactions in biologically diverse domains.

    1. Chromosomes and Gene Expression
    2. Developmental Biology
    F Javier DeHaro-Arbona, Charalambos Roussos ... Sarah Bray
    Research Article

    Developmental programming involves the accurate conversion of signalling levels and dynamics to transcriptional outputs. The transcriptional relay in the Notch pathway relies on nuclear complexes containing the co-activator Mastermind (Mam). By tracking these complexes in real time, we reveal that they promote the formation of a dynamic transcription hub in Notch ON nuclei which concentrates key factors including the Mediator CDK module. The composition of the hub is labile and persists after Notch withdrawal conferring a memory that enables rapid reformation. Surprisingly, only a third of Notch ON hubs progress to a state with nascent transcription, which correlates with polymerase II and core Mediator recruitment. This probability is increased by a second signal. The discovery that target-gene transcription is probabilistic has far-reaching implications because it implies that stochastic differences in Notch pathway output can arise downstream of receptor activation.