New insights into the Saccharomyces cerevisiae fermentation switch: Dynamic transcriptional response to anaerobicity and glucose-excess
BMC Genomics volume 9, Article number: 100 (2008)
The capacity of respiring cultures of Saccharomyces cerevisiae to immediately switch to fast alcoholic fermentation upon a transfer to anaerobic sugar-excess conditions is a key characteristic of Saccharomyces cerevisiae in many of its industrial applications. This transition was studied by exposing aerobic glucose-limited chemostat cultures grown at a low specific growth rate to two simultaneous perturbations: oxygen depletion and relief of glucose limitation.
The shift towards fully fermentative conditions caused a massive transcriptional reprogramming, where one third of all genes within the genome were transcribed differentially. The changes in transcript levels were mostly driven by relief from glucose-limitation. After an initial strong response to the addition of glucose, the expression profile of most transcriptionally regulated genes displayed a clear switch at 30 minutes. In this respect, a striking difference was observed between the transcript profiles of genes encoding ribosomal proteins and those encoding ribosomal biogenesis components. Not all regulated genes responded with this binary profile. A group of 87 genes showed a delayed and steady increase in expression that specifically responded to anaerobiosis.
Our study demonstrated that, despite the complexity of this multiple-input perturbation, the transcriptional responses could be categorized and biologically interpreted. By comparing this study with public datasets representing dynamic and steady conditions, 14 up-regulated and 11 down-regulated genes were determined to be anaerobic specific. Therefore, these can be seen as true "signature" transcripts for anaerobicity under dynamic as well as under steady state conditions.
In the majority of industrial fermentation applications of bakers' yeast (Saccharomyces cerevisiae), a high initial and sustained capacity to ferment the available sugar is a highly important characteristic, especially when the biomass is introduced in an application environment with high sugar concentrations and/or absence of oxygen. Despite several attempts [1–3], quantitative data concerning the dynamics of the adaptation to such industrially relevant fermentative conditions have not been obtained. The majority of studies published to date on fermentative capacity under defined conditions rely on the use of batch or chemostat cultures [4–6]. The high specific growth rate in batch cultures does not reflect typical industrial conditions for aerobic cultivation of yeast biomass and has a drastic impact on fermentative capacity . In chemostat cultures, which can be used in physiological studies to specifically investigate the effect of individual culture parameters, several physiological and transcriptional responses to the availability of oxygen and/or glucose have been identified [8–12]. However, in steady state chemostats dynamic responses to change in culture parameters can not be observed. A perturbation of one parameter in a chemostat cultivation results in a reproducible dynamic response from a defined constant culture . By means of such experimental set-up, short and long term dynamics have been studied to pulses of low glucose concentrations [14, 15].
The goal of the present study was to investigate the dynamic adaptation of S. cerevisiae to the industrially relevant transition from aerobic, sugar-limited and respiratory growth to fully fermentative (i.e., anaerobic glucose-excess) conditions and to dissect responses to the glucose up-shift and onset of anaerobicity. To this end, aerobic glucose-limited chemostat cultures grown at a moderate specific growth rate (0.10 h-1) were exposed to two simultaneous perturbations: a rapid depletion of oxygen and an increase of glucose concentration to a high value (40 g·l-1). Physiological analysis confirmed that the chemostat culture was fully respiratory before, and fully fermentative after the shift. Global dynamic responses to this combined perturbation were analyzed through genome-wide transcription analysis.
Results and discussion
To invoke rapid and full induction of fermentative capacity, respiratory, aerobic glucose-limited chemostat cultures (D = 0.1·h-1) were shifted to fully fermentative conditions by sudden depletion of oxygen and addition of glucose. The glucose was added two min after sparging the continuous culture with pure nitrogen, when the dissolved oxygen concentration had decreased from 75–80% to 10–15% of air saturation (Fig. 1). This raised the glucose concentration to 200 mM and ensured that the residual glucose concentration after 2 h of cultivation would still be above 100 mM, thus maintaining strong glucose catabolite repression throughout the experiment. (Fig. 2A) . Indeed, the sudden shift to fermentative conditions resulted in fully fermentative metabolism within the first 5 min with CO2, ethanol and glycerol as the major metabolic products (Fig. 2). This metabolic shift coincided with an increasing specific glucose consumption rate, up to 12-fold, over 2 h following the perturbation (Fig. 2B). The specific ethanol production rate, which under these anaerobic glucose-excess conditions reflects the culture's fermentative capacity, steadily increased to 19.6 mmol ethanol·g-1·h-1 (Fig. 2B). While the overall metabolic response was rapid and strong, the biomass concentration, the cell count and the cellular protein content did not change significantly throughout the experiment (Fig. 2C). During the two hours of the experiment, the growth rate did not exceed the starting growth rate of 0.1 h-1. The biomass therefore only contributed to 5% of the total carbon flux, while the main metabolic products (i.e. carbon dioxide, ethanol and glycerol) accounted for ca. 90% of the total carbon produced.
Microarray data processing and general transcriptional response
To identify genome-wide transcriptional changes connected to the induced metabolic adaptation, micro-array analysis was performed on samples from two independent replicate steady-state chemostat cultures and on samples taken 5, 10, 30, 60 and 120 min after glucose addition. The coefficient of variation between replicates was below 20%, which is comparable with previous chemostat-based transcriptome analyses [10, 11].
A first main concern was with normalization of these microarray data from non-steady-state culture samples. In previous transcriptome studies on steady-state chemostat cultures using Affymetrix microarrays, setting the average signal intensity of all probe-sets to a fixed value (also called global scaling) provided a good normalization method [10, 11]. As this normalization method might not be appropriate for dynamic cultivation conditions, we considered transcript levels of a few so-called 'house-keeping' genes commonly used as loading standards for Northern analysis and quantitative RT-PCR. After global scaling, the expression of ACT1, HHT2 and SHR3 (encoding respectively, for actin, histone and endoplasmic reticulum packaging chaperone protein) remained constant throughout the experiment with a variation coefficient around or below 20%. The stable transcript levels of house-keeping gene expression obtained with a global scaling approach indicated that no major changes in the total mRNA pools occurred during the experiment, which would require another type of normalization.
After global scaling, the significance of the changes in transcript levels during the dynamic experiment was estimated using the EDGE software (p-value threshold 0.005, , see Methods section for details). A set of 1923 genes was thus identified as being transcriptionally regulated in response to combined oxygen depletion and glucose addition (Additional file 1). This large group of genes was divided in several subgroups according to their expression profiles. 607 genes whose transcript levels increased after the perturbation were separated into four clusters according to their initial and later response (Clusters A-D; Fig. 3). 1316 genes with reduced transcripts responded rapidly to the perturbations (within 10 min) and were clustered according to their secondary response (Clusters 1–6; Fig. 3). All clusters were subsequently searched for overrepresentation of specific functional categories (as defined by MIPS ), and of promoter elements corresponding to specific transcriptional regulation networks (see Methods section). Upon a first inspection, some of clusters, despite subtle differences in their time-dependent transcript profiles, showed an overrepresentation of genes from the same functional categories. These were pooled to further improve the enrichment analysis. Thus cluster A and B, as well as 2 and 3, and also 4, 5 and 6 were pooled (Table 1 and 2), resulting in a final set of six different clusters.
Sudden relief from glucose limitation enables yeast cells to accelerate to a higher specific growth rate. Although faster growth was not observed in the 2 h after the relief of glucose limitation, over one third of the initially up-regulated genes were related to protein synthesis (Fig. 3; cluster A, B and C). This included a massive and fast up-regulation of genes within clusters A and B that encode components of the translational machinery, including 126 genes involved in rRNA synthesis, processing and modification and 49 genes involved in ribosomal biogenesis (Table 1). Taking into account that the total RNA pool mainly consists of rRNA , an up-regulation of rRNA synthesis was confirmed by an increase of the RNA content of the biomass after the relief from glucose limitation (Fig. 4).
Genes in cluster C displayed a sustained, slower increase of their transcript levels than those in clusters A and B. 37% of the genes in cluster C encoded ribosomal proteins. The delay between the expression of ribosomal biogenesis/rRNA genes and ribosomal protein genes is in line with previous observations indicating the existence of different regulatory mechanisms for these two groups of genes [20–22]. Accordingly, PAC and RRPE regulatory elements were enriched in the promoter regions of genes in clusters A & B, whereas Rap1p/Sfp1p and Fhl1p motifs were overrepresented in the promoter regions of cluster C genes (Table 2). In addition to the translational machinery, 57 genes involved in amino acid metabolism and 46 genes involved in nucleotide metabolism were up-regulated. This was consistent with the overrepresentation of Met32p , Gcn4p  and Bas1p  binding sites in the promoter regions of these genes, and indicated the need for synthesis of building blocks for transcription and translation.
Among the 1316 genes with reduced expression, one cluster comprising 122 genes showed rapid and strong repression (Cluster 1, Fig. 3). Although this cluster appeared relatively heterogeneous, one functional category was clearly enriched. It consists of seven transcription factor genes (ACE2, PRP45, OAF1, GTS1, SWI5, MSN1 and STB1) involved in various cellular functions, like fatty acid oxidation, stress response and cell cycle progression [26–31]. A large number of known targets of these transcription factors were also present in the down-regulated clusters (Additional file 2). Most of the remaining 1194 down-regulated genes were associated to metabolism and energy generation. In addition, a large number of genes involved in protein degradation (97 genes in total) were down-regulated, indicating a decreased requirement for proteolytic activity. Interestingly, 73 genes involved in stress response were down-regulated, including 18 related to oxidative stress response. This observation suggests that anaerobicity per se does not evoke an immediate stress for yeast.
As expected, the initial response to fully fermentative conditions showed quite some overlap with published datasets for glucose pulses to aerobic cultures [14, 32], including induction of the translational machinery and repression of the respiratory chain [33–35]. With this study, we aimed to go beyond the primary response to see how yeast adjusted to its altered growth environment.
We did not identify genes whose transcript levels continuously increased or decreased in the 2 h following the perturbation. At 30 min after the shift, a pivotal point appeared to be reached at which the transcript profiles either indicated a reverse regulation mode (clusters A, B, 4, 5 and 6) or a stable mRNA level (clusters C, 2 and 3). In this respect, a striking difference was observed between the transcript profiles of genes encoding ribosomal proteins and those encoding ribosomal biogenesis components (Fig. 5). A steady transcript level after 30 min of ribosomal proteins was indicative for a constitutive requirement for translational building blocks to support faster growth. In contrast, transcriptional up-regulated genes involved in the synthesis, processing and modification of the translational machinery appeared only to be temporarily required for a rapid adaptation to the new environmental conditions. In addition to the ribosomal protein genes, genes involved in de novo purine biosynthesis, methionine metabolism, and tetrahydrofolate-dependent C1 metabolism were continuously transcribed at an elevated level after 30 min. All three functional categories have previously been correlated with each other, and with a response to the decrease in the adenine nucleotide pool .
Also the initially down-regulated genes with a turning point after 30 min could be divided in two groups: steady pattern after 30 min (clusters 2 & 3) or again up-regulated after 30 min (clusters 4, 5 & 6). During the course of the experiment, the glucose concentration remained high and hence, functional categories known to be repressed by glucose were enriched among the clusters in which the transcript level remained low after 30 min. Regulatory factors involved in regulation of the respiratory chain (HAP2, HAP4 and HAP5) were down-regulated together with their targets . Stress-response genes also maintained low transcript level during the experiment. In contrast, transcripts of genes involved in lipid biosynthesis, reserve carbohydrate metabolism and protein degradation tended to increase again after 30 min. The large and coordinated transcriptional up-regulation of the translational machinery, specifically ribosomal proteins, complemented an opposite transcriptional regulation pattern of genes related to proteolytic activity. The down-regulation of target genes of the Mbp1/Swi6 complex, involved in G1 to S transition , correlated with a delay in cell cycle progression and correspondingly, a constant cell number over the two h monitored.
Many genes involved in the metabolism of storage carbohydrates (trehalose and glycogen) showed a decreased transcript level after the perturbation. To further investigate the observed changes in trehalose and glycogen metabolism, intracellular levels of trehalose and glycogen were measured. Both reserve carbohydrates were completely degraded within 30 min (Fig. 2D), consistent with a post-transcriptional activation of trehalose and glycogen phosphorylases [38, 39]. Physiological interpretation of the trehalose and glycogen degradation however, is less straightforward, since trehalose and glycogen are known to be involved in flux regulation, stress response and cell cycle .
Delayed responses related to anaerobiosis
Eighty-seven of the 1923 genes that showed a significantly altered transcript level after the combined glucose pulse and oxygen depletion only showed an increased transcript level after 30 min (cluster D). One of the few functional categories enriched within this group involved modification by glycosylation (ALG7, GNT1, MNT4, OST5, PMT2, PMT4, PMT5, SEC53 and SWP1). PMT2, PMT4 and PMT5 are specifically involved in O-linked mannosyl glycosylation, which is indispensable for cell wall integrity . In addition, this 'delayed response' cluster contained five of the nine genes encoding anaerobically induced mannoproteins (DAN1, DAN4, TIR1, TIR2 and TIR4) . Two other anaerobically induced mannoproteins (DAN2 and DAN3) were initially down-regulated, whereas transcript levels of the gene encoding the major cell wall mannoprotein (TIP1) did not significantly change at all.
A strongly anaerobiosis-related character of the genes in cluster D was not only suggested by the presence of the abovementioned genes involved in cell wall maintenance, but additionally by the presence of several genes involved in lipid transport (AUS1, FAA4 and DNF2), heme biosynthesis (HEM13; Rox1p repressed), sterol metabolism and regulation (ARE1, HES1 and NCP1), and cell wall biosynthesis (EXG2). Accordingly, a high number of genes contained AR1 elements in their promoter (Table 2), indicating a role of Upc2p [9, 42]. The delayed up-regulation of these 'anaerobic genes' indicated that the response to anaerobiosis is slow compared to the fast response to the relief from glucose limitation (clusters A, B and C).
Dissecting the response to anaerobiosis
The response to the anaerobic shift described in this study was compared with a dataset from a previously published study [20, 21], in which the transcriptional response of batch cultures was monitored for several generations after a shift from aerobic to anaerobic conditions. Surprisingly, only 51 genes were overlapping with the significant up-regulated genes of our study. Half of these resided in our delayed response cluster D, which contains many anaerobiosis-related genes. The absence of a glucose pulse in the study of Lai et al. [20, 21] explains the absence of genes encoding components of the translational machinery among the up-regulated genes in their dataset. Similarly, the large group of genes related to protein degradation found in the present study was not observed among the down-regulated genes identified by Lai et al. [20, 21]. A strong overlap (464 genes) was found between the down-regulated genes identified in the two studies. Most of this overlap resided in the constitutively low expressed clusters 2 & 3 of our study (45% of the genes overlapped), which include many genes related to oxidative stress response. The majority of genes within the functional category Stress Response responded slower in the anaerobic shift study of Lai et al. [20, 21] than in our study which included a step-up of the glucose concentration. Hence, we conclude that the observed regulation of stress response correlated with the relief from growth limitation rather than with a mere depletion of oxygen.
Anaerobic "signature" transcripts
In an attempt to identify robust 'signature transcripts' that show a consistent response to anaerobiosis, the set of significantly responding genes in this dynamic study was compared with several datasets from glucose pulses and aerobic-to-anaerobic shift experiments (Fig. 6)[14, 20, 32]. 457 genes of the 607 genes up-regulated in this study were previously identified in two other glucose-induced studies. Sixty-seven of the 150 non-overlapping genes resided within the delayed response of cluster D, indicating that more than 70% of the genes within cluster D were not responding to glucose. Twenty of these 150 genes are also up-regulated in the aerobic-to-anaerobic shift study of Lai et al. and can therefore been seen as specifically anaerobiosis-responsive (Fig. 6). These 20 anaerobic genes are involved in cell wall maintenance (DAN/TIR genes and related glycolysations), in membrane composition (WSC4, DAL5 and FET4) and metabolism (HEM13, MET13, ARE1, AUS1 and NCP1). Interestingly, 17 of those 20 genes resided in cluster D, and 12 genes contained the Upc2-binding promoter element (TCGTTTA), which earlier was associated with about 1/3 of anaerobic genes . Transcription factor Upc2 has been reported to be strictly regulated by heme and sterol levels . The delayed response of the anaerobic genes was likely due to the almost complete absence of growth during the experiment, thus sterol levels may not have been depleted rapidly through dilution over newly formed cells. A similar comparison with the down-regulated genes in our study resulted in 46 commonly responding genes (Fig. 6). The majority of the anaerobic down-regulated genes had functions related to mitochondrial function or oxidative stress response, of which 13 genes contained a Hap1 or Hap2/3/4/5 binding element [36, 44]. Heme levels are likely to respond rapidly to the depletion of oxygen from the culture, consistent with the fast response of the anaerobically down-regulated genes.
The use of transcripts as a diagnostic tool for biotechnological applications has been proposed previously [8, 11, 45]. Based on steady-state cultivation experiments, a consistent response to anaerobiosis had been determined by analyzing aerobic and anaerobic chemostat cultures grown under different nutrient limitations (carbon-, nitrogen-, phosphorus-, and sulfur limitation), resulting in 65 anaerobically up-regulated and 90 down-regulated genes . Surprisingly, only 14 up-regulated genes and 11 down-regulated genes were also found in both our dynamic study and the dynamic study of Lai et al.. Therefore, these can be seen as true "signature" transcripts for anaerobicity both within dynamic and steady state conditions (Table 3).
We have studied the induction of yeast fermentative capacity by switching a fully respiratory culture to fully fermentative conditions. The aerobic glucose-limited chemostat culture with a low specific growth rate became, as seen in the physiology measures, fully fermentative for the entire experiment due to a rapid depletion of oxygen and addition of a high glucose concentration (40 g·l-1). The shift caused a massive transcriptional reprogramming, where one third of all genes within the genome were transcribed differentially. Our study demonstrates that, despite the complexity of this multiple-input perturbation, the transcriptional responses could be categorized and biologically interpreted. This required clustering of genes that shared discernable time-dependent responses to the perturbation, followed by a systematic analysis of overrepresented gene categories and upstream regulatory elements. This approach revealed that this reprogramming of the transcriptome was mostly driven by relief from the glucose-limitation, exemplified by preparation for faster growth (induction of ribosomes, nucleotide biosynthesis and amino acids biosynthesis) and glucose repression of various metabolic pathways. Contrary to previous observations [9, 12], but as argued by Lai et al [20, 21], the apparent relief from stress clearly indicates that anaerobicity per se does not evoke a stress in yeast.
A recent study by our group  studied transcriptional responses in the first five min after a glucose pulse to aerobic, glucose limited chemostat cultures. While that study revealed important and virtually instantaneous transcriptional events after imposition of a relief from glucose limitation, the present study shows that transcriptional reprogramming continues well beyond this 5 min period. Interestingly, most responses changed character after the 30 minutes point. This is clearly illustrated by the difference between the expression pattern of genes encoding ribosomal proteins versus genes encoding components for ribosomal biogenesis. Therefore, we have used this experimental set-up for further studying molecular details governing the observed differences in the regulatory mechanisms of the various groups of genes (manuscript in preparation).
One exception to the binary response mechanism observed around 30 minutes is presented by the anaerobic induction response, which appears only after the initial response to the glucose pulse. Most of the genes specifically induced by anaerobiosis are related to cell wall and plasma membrane remodeling. This is in contrast with Lai et al. where this response was only apparent after one generation . The time span of anaerobic remodeling is therefore significantly shorter during a shift to complete fermentative metabolism on high glucose. By comparing this study with public datasets representing dynamic and steady conditions, the determined group of anaerobic "signature transcripts" will be better suited for use as a diagnostic tool in biotechnological applications.
Most of the transcriptional changes were due to sensitivity to the carbon supply. Still, the observed minor changes in transcripts for glycolytic enzymes cannot explain the 12-fold increase in flux through glycolysis under these conditions. Therefore we are presently studying the central carbon metabolism under such dynamic conditions by a multilevel approach, where transcripts, enzyme activities, metabolites and fluxes will be integrated. Hence, we will try to understand in more detail the regulatory mechanisms controlling fermentative capacity in yeast.
Strain and media
The S. cerevisiae strain used in this study was a prototrophic haploid reference strain CEN.PK113-7D (MATa) . Stock cultures were grown at 30°C in shake flasks containing 100 ml of synthetic medium with 20 g of glucose per liter.
The synthetic medium contained per liter of demineralized water 5 g of (NH4)2SO4, 3 g of KH2PO4, 0.5 g of MgSO4·7H2O, 0.15 ml of silicon antifoam (BDH), and trace element concentrations according to Verduyn et al. . After heat sterilization of the medium for 20 min at 120°C, a filter-sterilized vitamin solution  was added. The concentration of glucose in the reservoir medium was 7.5 g·l-1. This glucose was added to the synthetic medium after separate heat sterilization at 110°C.
CEN.PK113-7D (MATa) was grown at 30°C in 2-l bioreactors (Applikon) with a working volume of 1.5 l via an electrical level sensor. Removal of effluent from the center of the culture ensured that biomass concentrations in the effluent line differed by less than 1% from those in the culture . The dilution rate was set at 0.10 h-1. The pH was measured on-line and kept constant at 5.0 by the automatic addition of 2 M KOH using an Applikon ADI 1030 Biocontroller. A stirrer speed of 800 rpm and air flow of 0.75 liter·min-1 were applied to keep the dissolved-oxygen concentration, as measured with an oxygen electrode, above 60% of air saturation in all chemostat cultivations performed. Steady-state samples were taken after ~10 volume changes to avoid strain adaptation due to long-term cultivation [49, 50]. Biomass dry weight, metabolite, dissolved oxygen, and gas profiles were constant over at least three volume changes.
Anaerobic glucose-pulse experiments were started by sparging the medium reservoir of the fermentor of a steady-state glucose-limited aerobic chemostat culture (airflow of 0.5 liter·min-1) with pure nitrogen gas (Hoek-Loos, Schiedam, <5 ppm O2). Norprene™ tubing and butyl septa were used to minimize oxygen diffusion into the anaerobic cultures . Two min after nitrogen sparging and just before adding the glucose, the medium-supply and effluent-removal pump was switched off. The 200 mM (60 g of glucose monohydrate in 60 ml water) glucose pulse was injected aseptically through a rubber septum. Samples were taken 5, 10, 30, 60 and 120 min following glucose addition.
The exhaust gas was cooled by a condenser connected to a cryostat set at 2°C and dried with a Permapure™ dryer (Inacom Instruments) before analysis of the O2 and CO2 concentrations with a Rosemount NGA 2000 analyzer. The gas flow rate was determined with an Ion Science Saga digital flow meter.
Acetate, ethanol, glycerol, and glucose concentrations in supernatants were determined by HPLC analysis with a Bio-Rad Aminex HPX-87H column at 60°C. The column was eluted with 5 mM sulfuric acid at a flow rate of 0.6 ml min-1. Acetate was detected by a Waters 2487 dual-wavelength absorbance detector at 214 nm. Glucose, ethanol and glycerol were detected by a Waters 2410 refractive index detector.
Culture dry weights were determined as described in  while whole cell protein determination was carried out as described in . Cell numbers were counted by a Coulter counter (Multisizer II; Beckman Coulter) by using a 50 μm aperture.
Trehalose and glycogen
Trehalose and glycogen concentration measurements were performed as described previously  in duplicate measurements on two independent replicate cultures. Glucose was determined using the UV-method based on Roche kit no. 0716251.
Samples were collected during the pulse, washed three times with cold 5% trichloroacetic acid and the pellet is stored at -20°C. The samples were resuspended in 3% perchloric acid and heated at 90°C for 30 min. After centrifugation, the supernatant was mixed with 37% hydrochloric acid, containing 10 g l-1 orcinol monohydrate (crystalline, Sigma-Aldrich, Germany) and 5 g l-1 iron(III) chloride hexahydrate. The mixture was heated at 90°C for 20 min before measuring absorbance at 660 nm . Absorbance values were related to a concentration (expressed as μg·ml-1) using a calibration curve of a standard yeast RNA solution (Sigma-Aldrich, Germany).
Microarrays processing and analysis
Sampling of cells from chemostats, probe preparation, and hybridization to Affymetrix Genechip® microarrays were performed as described previously . The results for each time point after the perturbation (5, 10, 30, 60 and 120 min) were derived from two independently cultured replicates, while steady state data were derived from three independent chemostats. The complete dataset therefore comprised 13 arrays.
Acquisition and quantification of array images and data filtering were performed using Affymetrix GeneChip® Operating Software version 1.2. Before comparison, all arrays were globally scaled to a target value of 150 using the average signal from all gene features. To eliminate insignificant variations, genes with expression values below 12 were set to 12 and genes for which maximum expression was 20 over the 13 arrays were discarded. From the 9335 transcript features on the YG-S98 arrays, a filter was applied to extract 6383 yeast open reading frames, as previously described . To represent the variation in the measurements, the coefficient of variation was calculated as the mean deviation divided by the mean . The array data used in this study can be retrieved at Genome Expression Omnibus  with series number GSE8187.
For additional statistical analyses, Microsoft Excel running the EDGE (version 1.1.208) add-in was used  for a time course differential expression analysis. To determine the genes called significantly changed according to EDGE a p-value of 0.005 was used. K-means clustering of the genes with significantly changed expression levels was subsequently performed using Genedata Expressionist® Pro (version 3.1). The k-means algorithm used positive correlation as distance metric. The maximum number of iterations was set to 1000. Initially, the algorithm was run with k equal to 2, dividing the genes into an up- and a down-regulated cluster. Each cluster was then clustered again using k-means with k ranging from 2 to 10. The optimal k-value, i.e. 4 for the initially up-regulated and 6 for initially down-regulated genes, were based on the explained variance between clusters and the overrepresentation of functional categories (for detailed explanation please refer to Additional file 3).
Each cluster was consulted for enrichment in functional annotation and significant transcription factor (TF) binding (experimentally identified by Harbison et al. ) as described previously . In addition, specific TF binding sites not present in the Harbison dataset were analyzed by using web-based Regulatory Sequence Analysis Tools [11, 59].
Tanaka F, Ando A, Nakamura T, Takagi H, Shima J: Functional genomic analysis of commercial baker's yeast during initial stages of model dough-fermentation. Food Microbiol. 2006, 23: 717-728. 10.1016/j.fm.2006.02.003.
Higgins VJ, Beckhouse AG, Oliver AD, Rogers PJ, Dawes IW: Yeast genome-wide expression analysis identifies a strong ergosterol and oxidative stress response during the initial stages of an industrial lager fermentation. Appl Environ Microbiol. 2003, 69: 4777-4787. 10.1128/AEM.69.8.4777-4787.2003.
Novo M, Beltran G, Rozes N, Guillamon JM, Sokol S, Leberre V, Francois J, Mas A: Early transcriptional response of wine yeast after rehydration: osmotic shock and metabolic activation. FEMS Yeast Res. 2006, 7: 304-316. 10.1111/j.1567-1364.2006.00175.x.
van Hoek P, van Dijken JP, Pronk JT: Regulation of fermentative capacity and levels of glycolytic enzymes in chemostat cultures of Saccharomyces cerevisiae. Enzyme Microb Technol. 2000, 26: 724-736. 10.1016/S0141-0229(00)00164-2.
Thomsson E, Larsson C, Albers E, Nilsson A, Franzen CJ, Gustafsson L: Carbon starvation can induce energy deprivation and loss of fermentative capacity in Saccharomyces cerevisiae. Appl Environ Microbiol. 2003, 69: 3251-3257. 10.1128/AEM.69.6.3251-3257.2003.
van Hoek P, de Hulster E, van Dijken JP, Pronk JT: Fermentative capacity in high-cell-density fed-batch cultures of baker's yeast. Biotechnol Bioeng. 2000, 68: 517-523. 10.1002/(SICI)1097-0290(20000605)68:5<517::AID-BIT5>3.0.CO;2-O.
Weusthuis RA, Pronk JT, van den Broek PJ, van Dijken JP: Chemostat cultivation as a tool for studies on sugar transport in yeasts. Microbiol Rev. 1994, 58: 616-630.
Boer VM, de Winde JH, Pronk JT, Piper MD: The genome-wide transcriptional responses of Saccharomyces cerevisiae grown on glucose in aerobic chemostat cultures limited for carbon, nitrogen, phosphorus, or sulfur. J Biol Chem. 2003, 278: 3265-3274. 10.1074/jbc.M209759200.
Kwast KE, Lai LC, Menda N, James DT, Aref S, Burke PV: Genomic analyses of anaerobically induced genes in Saccharomyces cerevisiae: functional roles of Rox1 and other factors in mediating the anoxic response. J Bacteriol. 2002, 184: 250-265. 10.1128/JB.184.1.250-265.2002.
Piper MD, Daran-Lapujade P, Bro C, Regenberg B, Knudsen S, Nielsen J, Pronk JT: Reproducibility of oligonucleotide microarray transcriptome analyses. An interlaboratory comparison using chemostat cultures of Saccharomyces cerevisiae. J Biol Chem. 2002, 277: 37001-37008. 10.1074/jbc.M204490200.
Tai SL, Boer VM, Daran-Lapujade P, Walsh MC, de Winde JH, Daran JM, Pronk JT: Two-dimensional transcriptome analysis in chemostat cultures. Combinatorial effects of oxygen availability and macronutrient limitation in Saccharomyces cerevisiae. J Biol Chem. 2005, 280: 437-447. 10.1074/jbc.M501243200.
ter Linde JJ, Liang H, Davis RW, Steensma HY, van Dijken JP, Pronk JT: Genome-wide transcriptional analysis of aerobic and anaerobic chemostat cultures of Saccharomyces cerevisiae. J Bacteriol. 1999, 181: 7409-7413.
Flikweert MT, Kuyper M, van Maris AJ, Kotter P, van Dijken JP, Pronk JT: Steady-state and transient-state analysis of growth and metabolite production in a Saccharomyces cerevisiae strain with reduced pyruvate-decarboxylase activity. Biotechnol Bioeng. 1999, 66: 42-50. 10.1002/(SICI)1097-0290(1999)66:1<42::AID-BIT4>3.0.CO;2-L.
Kresnowati MT, van Winden WA, Almering MJ, ten Pierick A, Ras C, Knijnenburg TA, Daran-Lapujade P, Pronk JT, Heijnen JJ, Daran JM: When transcriptome meets metabolome: fast cellular responses of yeast to sudden relief of glucose limitation. Mol Syst Biol. 2006, 2 (): 49-10.1038/msb4100083. Epub 2006 Sep 12.
Ronen M, Botstein D: Transcriptional response of steady-state yeast cultures to transient perturbations in carbon source. Proc Natl Acad Sci U S A. 2006, 103: 389-394. 10.1073/pnas.0509978103.
Bisson LF, Fraenkel DG: Expression of kinase-dependent glucose uptake in Saccharomyces cerevisiae. J Bacteriol. 1984, 159: 1013-1017.
Storey JD, Xiao W, Leek JT, Tompkins RG, Davis RW: Significance analysis of time course microarray experiments. Proc Natl Acad Sci U S A. 2005, 102: 12837-12842. 10.1073/pnas.0504609102.
MIPS Functional catalogue. 2008, [http://mips.gsf.de/proj/funcatDB/]
Warner JR: The economics of ribosome biosynthesis in yeast. Trends Biochem Sci. 1999, 24: 437-440. 10.1016/S0968-0004(99)01460-7.
Lai LC, Kosorukoff AL, Burke PV, Kwast KE: Dynamical remodeling of the transcriptome during short-term anaerobiosis in Saccharomyces cerevisiae: differential response and role of Msn2 and/or Msn4 and other factors in galactose and glucose media. Mol Cell Biol. 2005, 25: 4075-4091. 10.1128/MCB.25.10.4075-4091.2005.
Lai LC, Kosorukoff AL, Burke PV, Kwast KE: Metabolic-state-dependent remodeling of the transcriptome in response to anoxia and subsequent reoxygenation in Saccharomyces cerevisiae. Eukaryot Cell. 2006, 5: 1468-1489. 10.1128/EC.00107-06.
Wade CH, Umbarger MA, McAlear MA: The budding yeast rRNA and ribosome biosynthesis (RRB) regulon contains over 200 genes. Yeast. 2006, 23: 293-306. 10.1002/yea.1353.
Blaiseau PL, Isnard AD, Surdin-Kerjan Y, Thomas D: Met31p and Met32p, two related zinc finger proteins, are involved in transcriptional regulation of yeast sulfur amino acid metabolism. Mol Cell Biol. 1997, 17: 3640-3648.
Hinnebusch AG: Translational regulation of GCN4 and the general amino acid control of yeast. Annu Rev Microbiol. 2005, 59: 407-450. 10.1146/annurev.micro.59.031805.133833.
Zhang F, Kirouac M, Zhu N, Hinnebusch AG, Rolfes RJ: Evidence that complex formation by Bas1p and Bas2p (Pho2p) unmasks the activation function of Bas1p in an adenine-repressible step of ADE gene transcription. Mol Cell Biol. 1997, 17: 3272-3283.
Albers M, Diment A, Muraru M, Russell CS, Beggs JD: Identification and characterization of Prp45p and Prp46p, essential pre-mRNA splicing factors. RNA. 2003, 9: 138-150. 10.1261/rna.2119903.
Ho Y, Costanzo M, Moore L, Kobayashi R, Andrews BJ: Regulation of transcription at the Saccharomyces cerevisiae start transition by Stb1, a Swi6-binding protein. Mol Cell Biol. 1999, 19: 5267-5278.
Karpichev IV, Small GM: Global regulatory functions of Oaf1p and Pip2p (Oaf2p), transcription factors that regulate genes encoding peroxisomal proteins in Saccharomyces cerevisiae. Mol Cell Biol. 1998, 18: 6560-6570.
McBride HJ, Yu Y, Stillman DJ: Distinct regions of the Swi5 and Ace2 transcription factors are required for specific gene activation. J Biol Chem. 1999, 274: 21029-21036. 10.1074/jbc.274.30.21029.
Rep M, Reiser V, Gartner U, Thevelein JM, Hohmann S, Ammerer G, Ruis H: Osmotic stress-induced gene expression in Saccharomyces cerevisiae requires Msn1p and the novel nuclear factor Hot1p. Mol Cell Biol. 1999, 19: 5474-5485.
Yaguchi S, Mitsui K, Kawabata K, Xu Z, Tsurugi K: The pleiotropic effect of the GTS1 gene product on heat tolerance, sporulation and the life span of Saccharomyces cerevisiae. Biochem Biophys Res Commun. 1996, 218: 234-237. 10.1006/bbrc.1996.0041.
Wang Y, Pierce M, Schneper L, Guldal CG, Zhang X, Tavazoie S, Broach JR: Ras and Gpa2 mediate one branch of a redundant glucose signaling pathway in yeast. PLoS Biol. 2004, 2 (5): E128-10.1371/journal.pbio.0020128. Epub 2004 May 11.
Gancedo JM: Yeast carbon catabolite repression. Microbiol Mol Biol Rev. 1998, 62: 334-361.
Johnston M: Feasting, fasting and fermenting. Glucose sensing in yeast and other cells. Trends Genet. 1999, 15: 29-33. 10.1016/S0168-9525(98)01637-0.
Warner JR: Synthesis of ribosomes in Saccharomyces cerevisiae. Microbiol Rev. 1989, 53: 256-271.
McNabb DS, Pinto I: Assembly of the Hap2p/Hap3p/Hap4p/Hap5p-DNA complex in Saccharomyces cerevisiae. Eukaryot Cell. 2005, 4: 1829-1839. 10.1128/EC.4.11.1829-1839.2005.
Alberghina L, Smeraldi C, Ranzi BM, Porro D: Control by nutrients of growth and cell cycle progression in budding yeast, analyzed by double-tag flow cytometry. J Bacteriol. 1998, 180: 3864-3872.
Francois J, Neves MJ, Hers HG: The control of trehalose biosynthesis in Saccharomyces cerevisiae: evidence for a catabolite inactivation and repression of trehalose-6-phosphate synthase and trehalose-6-phosphate phosphatase. Yeast. 1991, 7: 575-587. 10.1002/yea.320070605.
Francois J, Parrou JL: Reserve carbohydrates metabolism in the yeast Saccharomyces cerevisiae. FEMS Microbiol Rev. 2001, 25: 125-145.
Girrbach V, Strahl S: Members of the evolutionarily conserved PMT family of protein O-mannosyltransferases form distinct protein complexes among themselves. J Biol Chem. 2003, 278: 12554-12562. 10.1074/jbc.M212582200.
Abramova N, Sertil O, Mehta S, Lowry CV: Reciprocal regulation of anaerobic and aerobic cell wall mannoprotein gene expression in Saccharomyces cerevisiae. J Bacteriol. 2001, 183: 2881-2887. 10.1128/JB.183.9.2881-2887.2001.
Cohen BD, Sertil O, Abramova NE, Davies KJ, Lowry CV: Induction and repression of DAN1 and the family of anaerobic mannoprotein genes in Saccharomyces cerevisiae occurs through a complex array of regulatory sites. Nucleic Acids Res. 2001, 29: 799-808. 10.1093/nar/29.3.799.
Davies BS, Rine J: A role for sterol levels in oxygen sensing in Saccharomyces cerevisiae. Genetics. 2006, 174: 191-201. 10.1534/genetics.106.059964.
Hickman MJ, Winston F: Heme levels switch the function of Hap1 of Saccharomyces cerevisiae between transcriptional activator and transcriptional repressor. Mol Cell Biol. 2007, 27: 7414-7424. 10.1128/MCB.00887-07.
Higgins VJ, Rogers PJ, Dawes IW: Application of genome-wide expression analysis to identify molecular markers useful in monitoring industrial fermentations. Appl Environ Microbiol. 2003, 69: 7535-7540. 10.1128/AEM.69.12.7535-7540.2003.
van Dijken JP, Bauer J, Brambilla L, Duboc P, Francois JM, Gancedo C, Giuseppin ML, Heijnen JJ, Hoare M, Lange HC, Madden EA, Niederberger P, Nielsen J, Parrou JL, Petit T, Porro D, Reuss M, van Riel N, Rizzi M, Steensma HY, Verrips CT, Vindelov J, Pronk JT: An interlaboratory comparison of physiological and genetic properties of four Saccharomyces cerevisiae strains. Enzyme Microb Technol. 2000, 26: 706-714. 10.1016/S0141-0229(00)00162-9.
Verduyn C, Postma E, Scheffers WA, van Dijken JP: Effect of benzoic acid on metabolic fluxes in yeasts: a continuous-culture study on the regulation of respiration and alcoholic fermentation. Yeast. 1992, 8: 501-517. 10.1002/yea.320080703.
van den Berg MA, Jong-Gubbels P, Kortland CJ, van Dijken JP, Pronk JT, Steensma HY: The two acetyl-coenzyme A synthetases of Saccharomyces cerevisiae differ with respect to kinetic properties and transcriptional regulation. J Biol Chem. 1996, 271: 28953-28959. 10.1074/jbc.271.49.31243.
Ferea TL, Botstein D, Brown PO, Rosenzweig RF: Systematic changes in gene expression patterns following adaptive evolution in yeast. Proc Natl Acad Sci U S A. 1999, 96: 9721-9726. 10.1073/pnas.96.17.9721.
Jansen ML, Daran-Lapujade P, de Winde JH, Piper MD, Pronk JT: Prolonged maltose-limited cultivation of Saccharomyces cerevisiae selects for cells with improved maltose affinity and hypersensitivity. Appl Environ Microbiol. 2004, 70: 1956-1963. 10.1128/AEM.70.4.1956-1963.2004.
Visser W, Scheffers WA, WH BV, van Dijken JP: Oxygen requirements of yeasts. Appl Environ Microbiol. 1990, 56: 3785-3792.
Postma E, Kuiper A, Tomasouw WF, Scheffers WA, van Dijken JP: Competition for glucose between the yeasts Saccharomyces cerevisiae and Candida utilis. Appl Environ Microbiol. 1989, 55: 3214-3220.
Verduyn C, Postma E, Scheffers WA, van Dijken JP: Physiology of Saccharomyces cerevisiae in anaerobic glucose-limited chemostat cultures. J Gen Microbiol. 1990, 136: 395-403.
Parrou JL, Francois J: A simplified procedure for a rapid and reliable assay of both glycogen and trehalose in whole yeast cells. Anal Biochem. 1997, 248: 186-188. 10.1006/abio.1997.2138.
Dische Z: Qualitative and quantitative colorimetric determination of heptoses. J Biol Chem. 1953, 204: 983-997.
Genome Expression Omnibus. 2008, [http://www.ncbi.nlm.nih.gov/geo/]
Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takusagawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99-104. 10.1038/nature02800.
Knijnenburg TA, de Winde JH, Daran JM, Daran-Lapujade P, Pronk JT, Reinders MJ, Wessels LF: Exploiting combinatorial cultivation conditions to infer transcriptional regulation. BMC Genomics. 2007, 8: 25-10.1186/1471-2164-8-25.
van Helden J, Andre B, Collado-Vides J: A web site for the computational analysis of yeast regulatory sequences. Yeast. 2000, 16: 177-187. 10.1002/(SICI)1097-0061(20000130)16:2<177::AID-YEA516>3.0.CO;2-9.
Yeast Protein Database. 2008, [http://www.proteome.com]
We thank Erwin Suir for technical assistance, and specifically thank Theo Knijnenburg for valuable statistical advice and technical assistance. This project was financially supported by the IOP Genomics program of Senter Novem, The Netherlands.
JvdB carried out the experimental work, interpreted the results and drafted the manuscript. PDL, JTP and JHdW contributed to conception and design. PDL interpreted the results and assisted in structuring the manuscript. JTP and JHdW assisted in structuring the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table with all significant changed genes listed percluster. The Affymetrix ID, systematic name and the standard name were given for each gene. (XLS 154 KB)
Additional file 2: Figure with expression patterns of transcription factors Oaf1, Ace2, Swi5 and Gts1 (▲) and average expression patterns of their significant changed targets (□). Targets were all experimentally determined, as described at the proteome database of BIOBASE . Expression values were given in fold change between each time point and time t = 0 (steady state). (XLS 36 KB)
Additional file 3: Quality measurements of the K-means clustering. The k-values in a range from 2 to 10 were analyzed for its explained variance and the overrepresentation of functional categories. The quality of the individual clusters was measured by the same pairs proportion. (DOC 258 KB)
Authors’ original submitted files for images
About this article
Cite this article
van den Brink, J., Daran-Lapujade, P., Pronk, J.T. et al. New insights into the Saccharomyces cerevisiae fermentation switch: Dynamic transcriptional response to anaerobicity and glucose-excess. BMC Genomics 9, 100 (2008). https://doi.org/10.1186/1471-2164-9-100
- Chemostat Culture
- Fermentative Condition
- Translational Machinery
- Glucose Limitation