Rewiring of the Yeast Transcriptional Network Through the Evolution of Motif Usage

See allHide authors and affiliations

Science  05 Aug 2005:
Vol. 309, Issue 5736, pp. 938-940
DOI: 10.1126/science.1113833


Recent experiments revealed large-scale differences in the transcription programs of related species, yet little is known about the genetic basis underlying the evolution of gene expression and its contribution to phenotypic diversity. Here we describe a large-scale modulation of the yeast transcription program that is connected to the emergence of the capacity for rapid anaerobic growth. Genes coding for mitochondrial and cytoplasmic ribosomal proteins display a strongly correlated expression pattern in Candida albicans, but this correlation is lost in the fermentative yeast Saccharomyces cerevisiae. We provide evidence that this change in gene expression is connected to the loss of a specific cis-regulatory element from dozens of genes following the apparent whole-genome duplication event. Our results shed new light on the genetic mechanisms underlying the large-scale evolution of transcriptional networks.

Evolution of gene expression plays a prominent role in generating phenotypic diversity (13), but little is known about the genetic basis underlying broad modulations of the genome-wide transcription program. Here we describe a rewiring of the yeast transcriptional network that is linked to the massive loss of a conserved cis-regulatory element.

With the recent sequencing of many yeast genomes, the hemiascomycete lineage emerged as a central paradigm for studying the genetic basis of phenotypic diversity (4). Despite a strong conservation of gene content, yeast species exhibit major phenotypic differences. Most striking, perhaps, is their differential requirement for oxygen. Whereas growth of most yeast species requires oxygen, Saccharomyces cerevisiae grows rapidly in its absence and prefers to ferment glucose anaerobically even when oxygen is present. The emergence of this capacity is linked to the apparent whole-genome-duplication event (57), which allowed the specialization of enzymes to aerobic versus hypoxic conditions (4, 8) (Fig. 1A).

Fig. 1.

(A) Phylogenetic tree of the yeast species used in this study. (B) Expression similarity between pairs of genes in the two yeast species. Each matrix represents pairwise correlation coefficients between genes associated with ribosomal proteins (RP), rRNA processing (rRNA), mitochondrial ribosomal proteins (MRP), and the stress-related genes (STR). Gene sets were defined as in (10, 14, 15). The number of genes in each set for S. cerevisiae are 175 (RP), 50 (rRNA), 72 (MRP), and 139 (STR), and for C. albicans 90 (RP), 50 (rRNA), 39 (MRP), and 33 (STR).

The preference for fermentation is reflected in the organization of the S. cerevisiae transcriptional network. For example, genes involved in oxidative phosphorylation are repressed in the presence of glucose (9). To examine the transcription program of S. cerevisiae in more detail, we compiled a database of over 1000 published genome-wide expression profiles (10). In accordance with its anaerobic growth phenotype, we found that the expression of S. cerevisiae genes encoding mitochondrial functions, most notably the mitochondrial ribosomal proteins (MRP), is not correlated (and perhaps is even weakly anticorrelated) with that of genes coding for cytoplasmic ribosomal proteins (RP) and with ribosomal RNA (rRNA) processing genes (10) (Fig. 1B). Instead, expression of the MRP genes exhibits a distinct correlation with that of genes induced in response to environmental stress conditions. At least in part, this correlation reflects the induction of the stress-related genes during the relatively slow growth in nonfermentable carbon sources, such as glycerol or ethanol, which requires mitochondrial function.

To examine the manifestation of this relationship between the expression of MRP genes and genes required for rapid growth in aerobic yeast species, we assembled a data set of 198 genome-wide expression profiles of Candida albicans, a human pathogen that primarily grows aerobically. In accordance with the requirement for respiration to support rapid growth, we found that in C. albicans the expression of MRP genes is strongly correlated with the expression of genes involved in ribsome assembly, including RP genes and rRNA processing genes (Fig. 1B). Thus, the organization of the S. cerevisiae and C. albicans transcriptional networks reflects their respective physiologies.

Gene expression can evolve through mutations in trans-acting regulatory proteins or in cis-acting DNA elements (compare fig. S4). For example, the regulation of a common transcription factor could diverge between the species. Alternatively, modifications of the individual cis-regulatory elements of each downstream gene could change their regulation, potentially coopting or losing regulation by a transcription factor. To distinguish between these scenarios, we searched the regions 600 base pairs (bp) upstream of open reading frames of both organisms for overrepresented sequence elements. Specific regulatory sequences that appear to be conserved between S. cerevisiae and C. albicans were identified (11, 12). Most notably, the PAC (GATGAG) sequence (13) was associated with rRNA processing genes in both species. In agreement with previous reports (10), no overrepresented sequence was associated with the MRP genes in S. cerevisiae.

In contrast, in C. albicans, the MRP genes were clearly associated with an overrepresented upstream sequence motif (AATTTT). This sequence was previously implicated in the regulation of rRNA processing genes in S. cerevisiae (13), although its functional role was not demonstrated experimentally. To examine whether this motif contributes to the regulation of MRP genes in C. albicans, we fused the promoter of the MRP7 gene to a yellow fluorescent protein (YFP) reporter. As expected, strong induction of the YFP reporter was observed during exponential growth, whereas expression in stationary phase was much weaker (Fig. 2A). Mutating the AATTTT sequence in the MRP7 promoter to either a GC-rich or to an AT-rich sequence drastically reduced expression of the YPF (Fig. 2A), indicating that it is the AATTTT sequence itself, and not its AT-rich nature, that is important for its function. Similar results were obtained with two additional mitochondrial ribosomal protein gene promoters (fig. S5). We conclude that the AATTTT sequence functions as a cis-regulatory element of MRP genes in C. albicans.

Fig. 2.

(A) The AATTTT sequence is required to induce high-level transcription of a reporter protein fused to the C. albicans MRP7 promoter in exponentially growing cells. (B) Frequency of occurrence of RGE in a 50-bp window of the gene groups indicated. Background genomic frequency is indicated in red (15). Multiple occurrences of the RGE were included. The same quantitative result is obtained when only a single occurrence per gene is considered (fig. S2).

Examining the appearance of this motif in more detail, we found that its position relative to that of the ORF start codon is highly confined in both S. cerevisiae and C. albicans, although the precise position is somewhat different (160 versus 110 bp in S. cerevisiae and C. albicans, respectively). In both organisms, this sequence is significantly overrepresented in genes involved in rRNA processing and ribosomal proteins, but not in genes associated with the environmental stress response. In addition, only in C. albicans is it overrepresented also in MRP genes (Fig. 2B). Because this conserved element is associated with genes required for rapid growth in both organisms, we denote it as a rapid growth element (RGE).

This RGE could have been coopted to C. albicans MRP genes, or alternatively, could have been lost in S. cerevisiae. To distinguish between these alternatives, we examined the appearance of the RGE in nine sequenced yeast species that are considered to be intermediate in the evolution between S. cerevisiae and C. albicans (Fig. 3). In all species examined, the RGE is significantly overrepresented in genes involved in rRNA processing. Notably, overrepresentation of the RGE in MRP promoters was found in all genomes that diverged from the S. cerevisiae lineage before the whole-genome duplication event (Fig. 3). This indicates a widespread loss of the RGE in MRP genes following the genome duplication.

Fig. 3.

Frequency of occurrence of RGE in intermediate yeast species (15). The gene sets for S. cerevisiae were defined as in Fig. 1B. Gene sets for the remaining species consist of orthologs of the S. cervevisiae genes.

Taken together, it appears that the emergence of anaerobic growth capacity in yeast is associated with a global rewiring of the yeast transcriptional network. This rewiring involved changes in the promoter regions of dozens of genes, manifested by the loss of a specific regulatory motif from MRP gene promoters. It would be interesting to examine which transcription factor binds to the RGE site and whether its coding sequence was also modified in association with the emergence of anaerobic growth. Intriguingly, the most parsimonious scenario for the described promoter adaptation appears to be intimately linked with the whole-genome duplication event. Our results suggest that gene duplication can facilitate the evolution of new function not only by specialization of coding sequences but also by facilitating the evolution of gene expression.

Supporting Online Material


Figs. S1 to S5

Tables S1 to S4

References and Notes

View Abstract

Navigate This Article