How yeast re-programmes its transcriptional profile in response to different nutrient impulses
© Dikicioglu et al; licensee BioMed Central Ltd. 2011
Received: 5 May 2011
Accepted: 25 September 2011
Published: 25 September 2011
A microorganism is able to adapt to changes in its physicochemical or nutritional environment and this is crucial for its survival. The yeast, Saccharomyces cerevisiae, has developed mechanisms to respond to such environmental changes in a rapid and effective manner; such responses may demand a widespread re-programming of gene activity. The dynamics of the re-organization of the cellular activities of S. cerevisiae in response to the sudden and transient removal of either carbon or nitrogen limitation has been studied by following both the short- and long-term changes in yeast's transcriptomic profiles.
The study, which spans timescales from seconds to hours, has revealed the hierarchy of metabolic and genetic regulatory switches that allow yeast to adapt to, and recover from, a pulse of a previously limiting nutrient. At the transcriptome level, a glucose impulse evoked significant changes in the expression of genes concerned with glycolysis, carboxylic acid metabolism, oxidative phosphorylation, and nucleic acid and sulphur metabolism. In ammonium-limited cultures, an ammonium impulse resulted in the significant changes in the expression of genes involved in nitrogen metabolism and ion transport. Although both perturbations evoked significant changes in the expression of genes involved in the machinery and process of protein synthesis, the transcriptomic response was delayed and less complex in the case of an ammonium impulse. Analysis of the regulatory events by two different system-level, network-based approaches provided further information about dynamic organization of yeast cells as a response to a nutritional change.
The study provided important information on the temporal organization of transcriptomic organization and underlying regulatory events as a response to both carbon and nitrogen impulse. It has also revealed the importance of a long-term dynamic analysis of the response to the relaxation of a nutritional limitation to understand the molecular basis of the cells' dynamic behaviour.
The ability of a microorganism to adapt changes in its physicochemical (e.g. temperature , pH ) or nutritional [3, 4] environment is crucial for its survival. The yeast, Saccharomyces cerevisiae, has developed mechanisms to respond to such environmental changes in a rapid and effective manner; such responses may demand a widespread re-programming of gene activity [2, 5]. Among other transcription factors, Msn2p and Msn4p regulate the expression of ~200 genes in response to several stresses, including heat shock, osmotic shock, oxidative stress, low pH, glucose starvation, sorbic acid and high ethanol concentrations . This is especially true of changes in the nutrient environment and the ability to sense and respond to changes in nutrient availability is essential for cells from both unicellular and multicellular organisms. Glucose is the most abundant monosaccharide on earth and is the preferred carbon source for most organisms and, accordingly, changes in glucose availability often have profound consequences in many types of cell . A negative regulator of the glucose-sensing signal transduction pathway, Mth1p, is required for repression of transcription by Rgt1p. Mth1p interacts with Rgt1p and the glucose sensors Snf3p and Rgt2p to play one of the key regulatory roles in response to the amount of available glucose in the environment . Gcr1p and Gcr2p are transcription factors that activate genes involved in glycolysis, and are also among the major regulators mediating carbon catabolite repression . The introduction of glucose to a culture of S. cerevisiae cells growing by respiration evokes changes at both the level of gene expression and of metabolism, with several proteins being activated or deactivated and gene expression being completely re-programmed to accommodate the switch from respiration to fermentation. Many transcriptional regulators are involved in the process including the HAP complex, which is a transcriptional activator and a global regulator of respiratory gene expression . Carbon catabolite repression down-regulates the expression of genes that encode enzymes involved in gluconeogenesis, the Krebs cycle, respiration, mitochondrial development, and the utilization of carbon sources other than glucose, fructose or mannose . While the main effect of glucose is exerted at the transcriptional level , changes in mRNA and protein stability are also involved in the process [13, 14].
Ammonium assimilation in yeast occurs through its incorporation into glutamate, the source of nearly 80% of all cellular nitrogen . Growth on ammonium causes a decrease in the activities of the enzymes used to assimilate less favourable nitrogen sources. This phenomenon is termed nitrogen catabolite repression, although the effect is not as well characterised as its carbon counterpart, particularly with respect to sudden changes in ammonium availability. Much less is known of the cellular response to sudden changes in the concentration of ammonium or other nitrogen sources available to the cell. A complex regulatory scheme is invoked in diploid yeast cells when they are deprived of nitrogen, this can result in pseudohyphal growth, a process regulated by a set of transcriptional activators and repressors including Phd1p, Hms1p, Mga1p and Msn1p . It should be noted that, while ammonium is not one of the most preferred nitrogen sources for S. cerevisiae, the yeast grows well on ammonium and its presence evokes nitrogen catabolite repression . Ammonium is taken via two high-affinity permeases (Mep1p and Mep2p) as well as by a low-affinity permease (Mep3p). The expression of the GDH1, GLN1, and GAP1 genes is regulated by the concentration of ammonium present in the growth medium [17, 18]. The expression of nitrogen-regulated genes is controlled by both positively (Gln3p and Nil1p) and negatively acting proteins (Nil2p and Dal80p). In addition, it has been shown that the TOR kinases play an essential role in preventing the expression of nitrogen-regulated genes , and they probably have an important integrative role.
Several investigations of the transient responses of yeast metabolism to a sudden change in nutritional availability have been carried out. Kresnowati et al.  have investigated the transient short-term transcriptome and metabolome response of yeast cells to glucose perturbation in chemostats and have indicated that both the transcriptomic and metabolomic changes mediate two kinds of response - one concerned with the transition from fully respiratory to respiro-fermentative metabolism and the other with the increase in growth rate that is the consequence of an increase in nutrient supply. Ronen and Botstein  have investigated the transient transcriptional response to switching carbon sources between galactose and glucose and concluded that experimental designs that involve short-term transient perturbations may be useful in understanding dynamic metabolic regulatory networks. The transient response to nitrogen catabolite repression was investigated by introducing an ammonium pulse into a glutamine-limited culture  and showed that the ammonium-induced repression was not due to a generalised stress response but, instead, represented a specific signal for nitrogen catabolite regulation. The effect of sulphate or phosphate limitation in the growth medium, together with uracil and leucine deficiency, was also investigated and it has been deduced that the cells adjust their growth rate to nutrient availability and maintain homeostasis in the same way in both batch and steady-state conditions .
In this study, the dynamic re-organization of yeast's cellular activity was analyzed by following the short- and long-term transcriptomic response to a sudden relaxation of either carbon and nitrogen limitation by an impulse of glucose or ammonium, respectively. The experimental design was such that the specific perturbation was uniquely introduced into an otherwise carefully controlled environment. Thus a glucose impulse was given to a steady-state glucose-limited culture and an ammonium impulse to a corresponding ammonium-limited steady-state culture. The response of the yeast cells was monitored at the transcriptomic level until the steady state was re-established. Thus the time-scale of this investigation ranged from seconds to hours, allowing the elucidation of both the metabolic and regulatory switches that enable yeast cells to adapt to, and recover from, a transient change in nutrient availability. We believe that this study makes a significant contribution to our understanding of nutritional control in yeast since the response is studied over both short and long time-scales for two different nutrients under well-controlled physiological conditions.
Results and discussion
The immediate, as well as the adaptive (long-term), response to the release from nutritional limitation, followed by the system's slow return to the nutrient-limited steady state was investigated using a systems biology approach. Glucose (as a carbon source) and ammonium (as a nitrogen source) were injected into their respective nutrient-limited cultures in two matched fermenters operated in fully controlled chemostat mode. Samples for transcriptome analysis were taken at different time intervals, ranging from seconds to hours, until the culture had reached a second steady state.
Correlation analyses of genome-wide expression profiles
The gravimetrically determined biomass values and the optical densities at two steady states that were approximately 100 hr apart indicated that a significant change (loss or gain) of fitness was not detected (p-value > 0.01). The correlation between the expression profiles of two steady states was determined to be the highest among the individual profiles at each time point and the first steady state. There are ca. 50 generations between collecting the initial sample in the first steady state and the final sample in the second steady state. Thus, these observations indicate that spontaneous mutations during the course of fermentation are not a significant confounding factor population and so justify a detailed investigation of the transient response of the transcriptome.
Temporal organization of the global transcriptional response
The release from ammonium limitation, by providing an ammonium impulse, revealed a very different transcriptional response to that observed upon release from glucose limitation, in that the re-programming of gene expression started later and took longer both to complete and to return to the steady-state profile. Thus, the transcriptome profile recorded 20 sec. after the ammonium impulse was not significantly different to that of the preceding steady state. The profiles of cells collected 40 and 60 sec post-impulse were clustered with those from the 8 and 16 min samples and were still closely related to the steady state. The main impact of the ammonium impulse on gene transcription is seen in the period between 24 min. and 3 h. post-impulse, while the period 4-7 h. post-impulse represents a slow return to the steady-state profile (Figure 2B).
Genes showing a significant change in expression level in response to the nutrient impulses
Gene Ontology (GO) Annotations to Differentially Expressed Genes
Significantly Associated Process GO Term
Fraction of Differentially Expressed Subset Associated with the GO Term
Fraction of Transcriptome Associated with the GO Term
Carbon Catabolite Repression
Carboxylic metabolic processes
4.07 × 10-21
Aspartate family amino acid metabolic processes
1.09 × 10-6
Glutamine family amino acid metabolic processes
2.35 × 10-4
Methionine family amino acid metabolic processes
1.69 × 10-3
Serine family amino acid metabolic processes
9.85 × 10-3
Purine metabolic processes
2.37 × 10-9
4.20 × 10-8
4.79 × 10-8
Alcohol catabolic processes
1.19 × 10-7
Energy coupled proton transport
3.81 × 10-7
Nitrogen Catabolite Repression
3.95 × 10-16
5.85 × 10-6
8.02 × 10-6
2.98 × 10-4
Aspartate family amino acid metabolic process
5.38 × 10-4
Amino acid and derivative metabolic process
4.19 × 10-3
Response to the glucose impulse
A glucose impulse was found to elicit significant changes in the transcript levels of 372 genes which are associated with the following biological process terms: 'carboxylic acid metabolic processes'; 'aspartate, glutamine, methionine and serine family amino acid metabolic processes'; 'purine metabolic processes'; 'glycolysis'; 'oxidative phosphorylation'; 'alcohol catabolic processes'; 'energy-coupled proton transport' (Table 1) (Additional file 1). Kresnowati et al.  studied the changes in transcript levels during the first 6 min following a glucose impulse and have also reported significant changes in transcript levels belonging to energy, purine ribonucleotide, amino-acid metabolism, and signal transduction functional categories in the MIPS classification as a short-term response to shifting from glucose limitation to conditions where glucose was in excess.
Introduction of glucose also immediately stimulated the expression of genes associated with the sulphate assimilation pathway (c5) but, as the glucose levels started to decline again, the genes associated with this pathway were down-regulated gradually after the first hour and later recovered to levels similar to that of the preceding steady state. Expression of genes for transcription factors related to sulphur metabolism were also up-regulated within 5 minutes following the introduction of glucose . This is most likely to be associated with methylation, reflecting the huge demand for the post-transcriptional processing of rRNA to sustain the transiently boosted growth rate.
The transcripts in clusters c1 and c6 (Figure 4A) displayed a sustained response of either up-(c1) or down-regulation (c6) throughout the experiment after the introduction of the glucose pulse. However, analysis of the genes in these clusters failed to reveal their significant (p-value < 10-4) enrichment for any GO biological process category. Among the genes in c1 (the cluster displaying sustained up-regulation following the glucose impulse) were a sub-group of transcripts that were related to methylation: SAM1 and SAM2, whose products are S-adenosylmethionine synthetases. It has previously been reported that an increase in growth rate requires Sam1p, and further increases results in yet more demand for methyl donors to sustain rRNA modification, also requiring higher levels of Sam2p, a close homolog of Sam1p . This supports the idea that the initial stimulation of the expression of genes concerned with sulphur metabolism is associated with the increased demand for methyl donors. Five members of the 'alcohol catabolic' process, namely PFK1, PFK2, ENO2, TKL1 and CTS1 were also members of the up-regulated cluster c1. Cluster 6 contains genes that displayed sustained down-regulation following the glucose impulse and included several amino-acid metabolism genes: CIT2, CPA2, IDP2, ARG1 and CPA1 in the glutamine family amino-acid metabolic process; LYS20, LYS21, LYS9 and HOM3 in the aspartate family amino acid metabolic process; HOM3, CYS4 and FPR1 in homoserine metabolic process, as well as four members of the nicotinamide nucleotide metabolic process, PYC1, PYC2, ADH2 and ALD4.
Interestingly, the transcript levels of genes for glucose transporters did not go through any major change in response to a sudden shift from glucose-limited to glucose-abundant conditions. Expression of the high-affinity glucose transporters would be expected to be fully derepressed during the preceding glucose-limited steady state. However, of the genes encoding high-affinity glucose-repressible hexose transporters, only HXT7 displayed a significant down-regulation in the level of its transcript immediately following the pulse. Transcript levels for the other three genes encoding high-affinity glucose transporters (HXT2, HXT4, and HXT7) are up-regulated from 1 h post-impulse as the glucose concentration in the growth medium starts to fall.
Published values  for the poly (A) tail lengths of all mRNA molecules were checked in order to identify any possible differences in mRNA degradation since no direct measurement was available. The down-regulated transcripts were not found to be significantly enriched with short poly (A) tails neither for carbon or nitrogen catabolite repression with the distribution of poly (A) tail length among up-and down-regulated transcripts appearing to be random. Since the shortest mRNA half-lives in yeast were in the range of 3 to 6 minutes , even the transcripts of the samples taken within the first minute are likely to be the result of an increase in transcription activity rather than an effect of mRNA degradation.
Response to the ammonium impulse
Relieving nitrogen limitation in the fermentation with an ammonium impulse resulted in significant changes in the transcription levels of 369 genes (Additional file 3). The members of this gene set are significantly enriched for GO bioprocess annotations associated with: 'central carbon metabolism', including 'glycolysis'; 'gluconeogenesis'; 'proton transport' and 'oxidative phosphorylation'; as well as amino acid production pathways, such as 'aspartate family amino acid metabolic process' and 'amino acid and derivative metabolic process' (Table 1).
A similar clustering of genes with a significant change in their transcript levels following the ammonium impulse, using self-organizing maps, produced 9 groups with bioprocess GO terms that can be significantly associated with each subset (Figure 4B) (Additional file 4). It was observed that the cells respond to the ammonium impulse more slowly than they do to a carbon impulse.
Down-regulation of the transcripts clustered in c0 started after the first minute displaying a sharp decrease in the expression levels after the first hour. This cluster was significantly enriched for glycolytic genes whose expression levels recovered towards the second steady state. Another cluster (c1), which was also significantly enriched with glycolytic genes, exhibited a delayed up-regulated transcriptional profile. This indicated that recovery from nitrogen limitation allowed the yeast cells also to utilize glucose better, thus resulting in down-regulation of glycolytic genes in mid-length response periods and then an up-regulation towards the cessation of the effect of the pulse as the re-establishment of high glucose concentrations resulted in the cells switching back into fermentative metabolism. Clusters significantly enriched with oxidative phosphorylation and trans-membrane ion transport processes were observed to display a down-regulation trend having the most distinct down-regulation between the 3rd and the 5th hours, recovering towards the second steady state (c2 and c4, respectively). This might have been due to the presence of excess glucose repressing respiration-related events during this latter period.
The transient abundance of ammonium led to an up-regulation of genes concerned with the process, and regulation, of translation. This up-regulated expression profile was displayed at the minute and hour timescales in c2, c5, c7, and c8. Clusters that were significantly enriched with 'translation process' terms (c5, c7 and c8) were also significantly enriched for 'cellular biosynthetic process' (p-value < 10-25 (c2, c5), and p-value < 10-37 (c8)), which indicated an up-regulation of growth-related events following a release from ammonium limitation. The induction of growth also required a higher demand for the methylation of tRNAs and rRNAs. The expression level of SAM1 in c2 was also observed to be up-regulated as is the case for the glucose impulse.
The expression levels of transcripts that were enriched with 'cation transport process' were sharply turned off around the first hour following the pulse (c4). Among the members of this cluster, an ammonium permease, Mep2p, works in conjunction with Pmp1p, Pmp3p and Pma1p to facilitate the trans-membrane transport of the slightly acidic ammonium during the uptake of the nitrogen source. This might have been due to the fast consumption of ammonium at that time, altering the intracellular pH, which resulted in the down-regulation of the relevant genes, only to be up-regulated again at later time points.
Dynamic transcriptional reprogramming of the cell during the transition created by a nutrient impulse
Temporal organization of dynamic regulatory events within the transcriptional response of yeast cells to a nutrient impulse was investigated using two different systems-based approaches, namely Dynamic Regulatory Event Miner (DREM)  and Negative Positive Network Analysis (NP) . A hidden-input/hidden-output Markov model integration of protein-DNA interactions and dynamic transcriptome data has been used in the Dynamic Regulatory Events Miner for identifying bifurcation events in the time series, where sets of genes which previously had roughly similar expression level diverge under the regulation of transcription factors that are selectively responsible for the controlling the expression of a certain subset of genes. A network-based modular approach was adapted to incorporate dynamic transcriptome data with protein-protein and protein-DNA interactions in the Negative-Positive Network Analysis in order to identify modular activity in response to catabolite repression.
The dynamic programming of the cells in response to nitrogen and carbon catabolite repression was identified using DREM . The dynamic reprogramming of the cells in response to a perturbation causing a change in nutrient availability exhibits a more complex pattern when carbon limitation is relieved than when ammonium was added to nitrogen-limited culture.
Gene Ontology Enrichment of Transcripts Following Bifurcation
Branch of Down-regulated Transcripts
Branch of Up-regulated Transcripts
Bifurcation Controlling TFs
Process Gene Ontology Terms
Bifurcation Controlling TFs
Process Gene Ontology Terms
Carbon Catabolite Repression
Mth1p, Msn2p, Msn4p, Hap2p, Swi4p, Ume6p
ion trans-membrane transporter activity
5.3 × 10-5
Mth1p, Msn2p, Msn4p, Hap2p, Swi4p, Ume6p
Endoplasmic reticulum membrane
4.7 × 10-5
G1-specific transcription in mitotic cell cycle
5.0 × 10-3
Microtubule associated complex
2.8 × 10-3
8.4 × 10-5
1.0 × 10-3
8.1 × 10-3
1.3 × 10-3
Nitrogen Catabolite Repression
Hms1p, Mga1p, Msn1p, Phd1p, Hap2p
ion trans-membrane transport
5.9 × 10-3
ribosome biogenesis and assembly
2.9 × 10-3
Pulse injection of ammonium sulphate into fermentation medium resulted in a bifurcation of transcripts into up-regulated and down-regulated branches (Figure 5B). However, the observed response was delayed for 20 seconds, similar to what has been observed by the hierarchical clustering of the 20th second time point with the steady state samples. A single transcription factor, Gcr2p; glycolysis regulatory protein, was significantly responsible for the split of the up-regulated branch with transcripts significantly enriched for ribosome biogenesis and assembly. Five transcription factors four of which were related to nitrogen starvation and pseudohyphal growth directly (Hms1p, Mga1p, Msn1p and Phd1p), as well as another transcription factor activating respiratory gene expression and a member of the complex facilitating the cross-pathway regulation of carbon and nitrogen metabolisms ; Hap2p were significantly associated with the down-regulated branch of transcripts which were enriched for ion trans-membrane transport (Table 2).
Network reduction via NP-analysis
GO biological process terms of modules with top 10 highest detection parameters
regulation of transcription, DNA-dependent
rRNA metabolic process
regulation of translation
cellular response to heat
regulation of transcription, DNA-dependent
ubiquitin-dependent protein catabolic process
regulation of translation
regulation of transcription, DNA-dependent
regulation of transcription, DNA-dependent
ubiquitin-dependent protein catabolic process
protein amino acid phosphorylation
regulation of translation
The transcription factors among the module interface genes, which had inter-modular interactions with genes from other modules, were investigated in terms of their significance in inter-modular communication. 28 out of 32 transcription factors in module P were interface genes that were interacting with module D, indicating the importance of inter-modular regulation of cellular growth and cellular trafficking. Growth and energy metabolism, on the other hand were observed not to possess enriched regulatory communication. Rather, the regulation of cellular trafficking enabled the indirect communication between these two processes. Only two TFs were identified in Module S, which was associated with cellular proliferation and none were identified as interface genes.
We have examined, over both short and long time-scales, the dynamic re-organization of gene expression in S. cerevisiae cells in response to a sudden relaxation of either carbon or nitrogen limitation using a system-based integrative approach. The observation of the genome-wide response at both levels, in a wide-ranging time span from seconds to hours, revealed metabolic and regulatory switches of yeast cells to adapt to and recover from an impulse-like perturbation.
The transcriptional response to impulse like addition of glucose was immediate whereas the response to ammonium was delayed for approximately 20 seconds. A larger change in the magnitude of expression of differentially transcribed genes was observed in response to carbon catabolite repression than to nitrogen catabolite repression. The transcriptional response was time-scale dependent. The response to glucose perturbation was different in the first minute, the first hour, the early hours and the late hours post-impulse both from the glucose-limited steady-state expression levels and from each other. On the other hand, the 20 seconds delay in response to ammonium perturbation caused a shift in the timescale with the remaining seconds-scale samples responding in a similar fashion to those from the early minutes, and the remaining minutes-scale samples responding in a similar fashion to those from the early hours.
The temporal distribution of priority of the biological processes was observed from the expression levels of the genes that were significantly associated with these processes. The most immediate response to the impulse-like addition of glucose was observed as the immediate down-regulation of the respiration-related transcripts within the first minute of induction. The up-regulation of growth-related transcripts that were significantly associated with translation and sulphate assimilation processes followed this response, reaching a maximum level of expression within the first half hour post-impulse, indicating that the increased demand for methyl donors, which are required for translational machinery, was met almost simultaneously with translation. The up-regulation of the glycolytic genes was the most delayed response. The expression of only one high-affinity hexose transporter, HXT7, was decreased in response to a surplus of glucose, possibly indicating a preferential use of this particular gene as a response to sudden changes in the availability of glucose. HXT7 was previously reported as having the highest affinity for glucose  and being the high affinity transporter, whose repression was not mediated through Mig1p . The selective down-regulation of this gene might be an indicator of a Mig1p-independent transient response to changes in the extracellular availability of glucose. Interestingly, the expression levels of the genes involved with central carbon metabolism were significantly affected by the impulse-like addition of ammonium. This might have resulted from the cross-pathway regulation of carbon and nitrogen metabolisms in yeast through common transcription factors including the HAP complex. In fact, dynamic regulatory analysis identified Hap2p as the sole common transcriptional regulator of the response to either one of the two different nutritional perturbations. The down-regulation of the expression of glycolytic genes could indicate that the yeast might actually perceive the relief from ammonium limitation as a release from an environment of excess glucose, possibly explaining this combinatorial change as a relief actually from carbon catabolite repression in addition to being subjected to nitrogen catabolite repression.
Dynamic regulatory analysis was able to reveal divergence patterns in the temporal profiles of gene expression on the basis of cellular component and localization criteria, specifically highlighting different patterns for microtubule-associated complex, endoplasmic reticulum and the proteasome complex as a response to carbon catabolite repression. The presence of a relationship between the expression levels of the genes associated with retrograde transport and aerobic respiration as well as the G1-specific transcription in the mitotic cell cycle could be identified although the detailed mechanics and the identification of common transcription factors regulating these processes require further investigation.
Dynamic correlated network analysis revealed that yeast's response to nitrogen catabolite repression could be represented as a single global module, but as four complementary and interacting modules in charge of growth, energy, transport/trafficking, and cellular proliferation in response to carbon catabolite repression. The interactions among the modules were regulated by inter-modular transcription factors. Specifically, the regulation of cellular trafficking and transport acted as a bridge in the regulation of energy and growth metabolisms. The down-regulation in the expression of a set of unknown genes as a response to growth-inducing perturbations requires further investigation from an evolutionary perspective as these genes might be candidates in phenotypic analysis studies for the determination of haploinsufficiency or complex haploinsufficiency.
The present study revealed the importance of long-term analysis of the response to the relaxation from nutritional deprivation to understand the molecular basis of the dynamic behaviour of the cells. A further detailed systems-based study that integrates additional levels of functional genomics analyses may provide further information on the dynamic re-organization of yeast cells to changing environmental conditions.
Strain and Growth Conditions
Wild type BY4743 (MAT a /MATα his3Δ/his3Δ leu2Δ/leu2Δ LYS2/lys2Δ MET15/met15Δ ura3Δ/ura3Δ; ) was cultivated in 2L fermenters (Applikon®) with 1L working volume under aerobic conditions in glucose-or ammonium-limited F1 media  in chemostat mode at a dilution rate of 0.1hr-1. Temperature and pH were controlled to 30°C and pH 4.5, respectively. Fermenters were stirred at 800 rpm which, together with constant air flow at a rate of 0.1 vvm, provided dissolved oxygen at ≥ 80% dO2 saturation at all times during cultivation.
Pulse injections and sampling
Samples for transcriptome analyses were collected at all time points. Biomass concentration was determined gravimetrically.
Sampling for transcriptome analysis, RNA isolation and transcriptome analysis
A culture sample (20 ml) was centrifuged at 4000 rpm for 3 min. Most of the supernatant was discarded, allowing re-suspension of cells in a small volume of growth medium. The cell suspension was released into liquid nitrogen and stored at -80°C until RNA isolation. Total RNA was isolated as described by . Total RNA was qualitatively assessed on an Agilent 2100 Bioanalyser (Agilent Technologies) and quantified using Nanodrop ultra-low-volume spectrophotometer (Nanodrop Technologies). cDNA was synthesised, and double-stranded cDNA was retrieved from ca.15 μg of total RNA as described in the Affymetrix GeneChip®Expression Analysis Technical Manual, using appropriate kits. cDNA was checked for quality using the Agilent 2100 Bioanalyser and was quantified using Nanodrop. Biotin-labelled cRNA was synthesized and was purified using clean up kits and it was quantified using the Nanodrop spectrophotometer before hybridization. Hybridization and loading onto Affymetrix Yeast2 arrays were carried out as described in the GeneChip®Expression Analysis Technical Manual. The chips were then loaded into a fluidics station for washing and staining using Microarray Suite 5 with EukGe W S2v4 programme. Lastly, the chips were loaded onto the Agilent GeneArray scanner 2500 and another quality check was performed using Microarray Suite 5 .
Microarray data acquisition and analysis
The raw data files were assessed with dChip software for outliers at the array level as well as at the probe-set level . Different nutritional conditions (glucose and ammonium pulse experiments) were treated as different sets and were assessed for their quality control separately. RMA Express software was then used to normalize the data, again as two separate data sets . The data was log2 transformed prior to analysis. In compliance with MIAME guidelines , the microarray data from this study has been submitted to ArrayExpress at the European Bioinformatics Institute under accession number [E-MTAB-643].
In order to identify transcripts whose expression significantly differed from steady-state levels following the nutrient pulse [40, 41], the software package EDGE  was used. The 'within-class' analysis carried out by EDGE was used to look for any differential expression over time within a single group. A significance measure was then assigned to the transcript level of each gene via the Q-value methodology, using the Benjamini-Hochberg correction for the calculation of false-discovery rates. Microsoft Excel Built-In commands were used to calculate the Pearson correlation coefficients of the transcriptome at the first steady state and the rest of the sampling times. All p-values were corrected for the false-positives introduced by the multiple testing presented by many time points using Bonferroni correction and 10-3 was selected as the cut-off threshold for p-values. GeneCluster 2.0  was used for clustering of significant transcripts via self-organizing maps and Hierarchical Clustering Explorer (HCE) 3.0  was used for hierarchical clustering purposes. The significantly enriched functional categories and the process ontology terms of the genes falling into the same cluster were determined by Saccharomyces Genome Database GO Term Finder tool  or AmiGO Term Finder tool . The threshold p-value was selected as 10-3. Transcription factors (TF) were taken from two sources TRANSFAC Professional Gene Transcription Factor Database  and YEASTRACT  and TFs that are common in both databases were considered in further analyses.
Dynamic Regulatory Events Miner analysis
Bifurcation points were determined using a hidden-input/hidden-output Markov model based software, DREM - The Dynamic Regulatory Events Miner as described by the authors . Different nutritional perturbations, where glucose and ammonium sulphate were introduced into their corresponding limited cultures, were analyzed independently. The chromatin immunoprecipitation experiments (chip-CHIP experiments) from which the TF-gene interactions were acquired were garnered from [49, 50].
Negative-Positive Network Analysis
The Negative-Positive (NP) Network was constructed as described by . The protein-protein and the protein-DNA interactome for Saccharomyces cerevisiae was obtained from the compilation available in the MATISSE software . A total of 5539 genes for which both expression and interaction data was available were considered for analysis. The threshold for the Pearson correlation coefficient (PCC) was selected as |PCC| ≥ 0.7. Following network reduction, Hierarchical Clustering Explorer (HCE) 3.0  was used for module identification. The largest possible modules with less than 1% intra-cluster anti-correlated links were manually dissected. MATLAB7.0 (R2006b) (The MathWorks Inc.) and Python 2.6 (Python Software Foundation) were utilized for computational purposes.
where n was the number of genes associated with the GO term and t was the total number of genes in the module. A GO term was used iff dp > (0.05/g) where g was the total number of GO terms associated with all the genes in a module . GO Slim terms were also used in conjunction with dp. The definition of the interface and the core members were also adapted from .
The authors greatly acknowledge the financial support for the research from the BBSRC (Grant BB/C505140/1 to SGO), and the travel grants for DD kindly provided by the Research Council of Turkey (TUBITAK) through the BDP programme and the Turkish State Planning Organization DPT09K120520. The research was also financially supported by Bogazici University Research Fund through Project No 631 and TUBITAK through Project No 106M444. Further support came from European Commission though the Coordination Action Project YSBN (Contract No.018942 to both BK and SGO) and UNICELLSYS Collaborative Project (No. 201142 to SGO).
- Shin DY, Matsumoto K, Iida H, Uno I, Ishikawa T: Heat shock response of Saccharomyces cerevisiae mutants altered in cyclic AMP-dependent protein phosphorylation. Mol Cell Biol. 1987, 7 (1): 244-250.PubMed CentralView ArticlePubMedGoogle Scholar
- Viladevall L, Serrano R, Ruiz A, Domenech G, Giraldo J, Barcelo A, Arino J: Characterization of the Calcium-mediated response to alkaline stress in Saccharomyces cerevisiae. J Biol Chem. 2004, 279 (42): 43614-43624. 10.1074/jbc.M403606200View ArticlePubMedGoogle Scholar
- Kresnowati MTAP, van Winden WA, Almering MJH, 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-PubMed CentralView ArticlePubMedGoogle Scholar
- Ronen M, Botstein D: Transcriptional response of steady-state yeast cultures to transient perturbations in carbon source. P Natl Acad Sci USA. 2006, 103 (2): 389-394. 10.1073/pnas.0509978103.View ArticleGoogle Scholar
- Causton HC, Ren B, Koh SS, Harbison CT, Kanin E, Jennings EG, Lee TI, True HL, Lander ES, Young RA: Remodelling of yeast genome expression in response to environmental changes. Mol Biol Cell. 2001, 12: 323-337.PubMed CentralView ArticlePubMedGoogle Scholar
- Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11: 4241-4257.PubMed CentralView ArticlePubMedGoogle Scholar
- Moriya H, Johnston M: Glucose sensing and signaling in Saccharomyces cerevisiae through the Rgt2 glucose sensor and casein kinase I. P Natl Acad Sci. 2004, 101 (6): 1572-1577. 10.1073/pnas.0305901101.View ArticleGoogle Scholar
- Lafuente MJ, Gancedo C, Jaunlaux J-C, Gancedo JM: Mth1 receives the signal given by the glucose sensors Snf3 and Rgt2 in Saccharomyces cerevisiae. Mol Microbiol. 2000, 35 (1): 161-172. 10.1046/j.1365-2958.2000.01688.xView ArticlePubMedGoogle Scholar
- Uemra H, Fraenkel DG: Glucose metabolism in gcr mutants of Saccharomyces cerevisiae. J Bacteriol. 1999, 181 (15): 4719-4723.Google Scholar
- Forsburg SL, Guarente L: Identification and characterization of HAP4: a third component of the CCAAT-bound HAP2/HAP3 heteromer. Genes Dev. 1989, 3: 1166-1178. 10.1101/gad.3.8.1166View ArticlePubMedGoogle Scholar
- Meijer MCM, Boonstr J, Verkleij AJ, Verrips CT: Glucose repression in Saccharomyces cerevisiae is related to the glucose concentration rather than the glucose flux. J Biol Chem. 1998, 273 (37): 24102-24107. 10.1074/jbc.273.37.24102View ArticlePubMedGoogle Scholar
- Gancedo JM: Yeast carbon catabolite repression. Microbiol Mol Biol R. 1998, 62 (2): 334-361.Google Scholar
- Foat BC, Houshmandi SS, Olivas WM, Bussemaker HJ: Profiling condition specific, genome-wide regulation of mRNA stability in yeast. P Natl Acad Sci. 2005, 102 (49): 17675-17680. 10.1073/pnas.0503803102.View ArticleGoogle Scholar
- Kim SY, Herbst A, Tworkowski KA, Salghetti SE, Tansey WP: Skp2 regulates Myc protein stability and activity. Mol Cell. 2003, 11: 1177-1188. 10.1016/S1097-2765(03)00173-4View ArticlePubMedGoogle Scholar
- Magasnik B: Ammonia Assimilation by Saccharomyces cerevisiae. Eucaryot Cell. 2003, 2 (5): 827-829. 10.1128/EC.2.5.827-829.2003.View ArticleGoogle Scholar
- Lorenz MC, Heitman J: Regulators of pseudohyphal differentiation in Saccharomyces cerevisiae identified through multicopy suppressor analysis in ammonium permease mutant strains. Genetics. 1998, 150: 1443-1457.PubMed CentralPubMedGoogle Scholar
- Magasnik B, Keiser CA: Nitrogen regulation in Saccharomyces cerevisiae. Gene. 2002, 290: 1-18. 10.1016/S0378-1119(02)00558-9View ArticleGoogle Scholar
- ter Schure EG, Sillje HHW, Vermeulen EE, Kalhorn JW, Verkleij AJ, Boonstra J, Verrips CT: Repression of nitrogen catabolic genes by ammonia and glutamine in nitrogen-limited continuous cultures of Saccharomyces cerevisiae. Microbiology. 1998, 144: 1451-1462. 10.1099/00221287-144-5-1451View ArticlePubMedGoogle Scholar
- Saldanha AJ, Brauer MJ, Botstein D: Nutritional homeostasis in batch and steady-state culture of yeast. Mol Biol Cell. 2004, 15 (9): 4089-4104. 10.1091/mbc.E04-04-0306PubMed CentralView ArticlePubMedGoogle Scholar
- Castrillo JI, Zeef LA, Hoyle DC, Zhang N, Hayes A, Gardner DJC, Cornell MJ, Petty J, Hakes L, Wardleworth L, Rhash B, Brown M, Dunn WB, Broadhurst D, O'Donoghue K, Hester SS, Dunkley TPJ, Hart SR, Swainston N, Li P, Gaskell SJ, Paton NW, Lilley KS, Kell DB, Oliver SG: Growth control of the eukaryote cell: a systems biology study in yeast. Journal of Biology. 2007, 6: 4- 10.1186/jbiol54PubMed CentralView ArticlePubMedGoogle Scholar
- Gutteridge A, Pir P, Castrillo JI, Charles PD, Lilley KS, Oliver SG: Nutrient control of eukaryote cell growth: a systems biology study in yeast. BMC Biol. 2010, 8: 68- 10.1186/1741-7007-8-68PubMed CentralView ArticlePubMedGoogle Scholar
- Tamayo P, Slonim D, Mesirov J, Zhu Q, Dmitrovsky E, Lander ES, Golub TR: Interpreting gene expression with self-organizing maps: Methods and application to hematopoeitic differentiation. P Natl Acad Sci. 1999, 96: 2907-2912. 10.1073/pnas.96.6.2907.View ArticleGoogle Scholar
- Beilharz TH, Preiss T: Widespread use of poly-A tail length control to accentuate expression of the yeast transcriptome. RNA. 2007, 13: 982-997. 10.1261/rna.569407PubMed CentralView ArticlePubMedGoogle Scholar
- Carlini DB: Context-dependent codon bias and mRNA longevity in the yeast transcriptome. Math Biosci Eng. 2005, 22 (6): 1403-1411.Google Scholar
- Ernst J, Vainas O, Harbison CT, Simon I, Bar-Josef Z: Reconstructing Dynamic Regulatory Maps. Mol Syst Biol. 2007, 3: 74-PubMed CentralView ArticlePubMedGoogle Scholar
- Xia K, Dong D, Xue H, Zhu S, Wang J, Zhang Q, Hou L, Chen H, Tao R, Huang Z, Fu Z, Chen YG, Han JD: Identification of the proliferation/differentiation switch in the cellular network of multicellular organisms. PLoS Comput Biol. 2006, 2: e145- 10.1371/journal.pcbi.0020145PubMed CentralView ArticlePubMedGoogle Scholar
- Lang T, Schaeffeler E, Bernreuther D, Bredschneider M, Wolf DH, Thumm M: Aut2p and Aut7p, two novel microtubule-associated proteins are essential for delivery for autophagic vesicles to the vacuole. EMBO J. 1998, 17: 3597-3607. 10.1093/emboj/17.13.3597PubMed CentralView ArticlePubMedGoogle Scholar
- Dang VD, Bohn C, Bolotin-Fukuhara M, Daignan-Fornier B: The CCAAT box-binding factor stimulates ammonium assimilation in Saccharomyces cerevisiae, defining a new cross-pathway regulation between nitrogen and carbon metabolisms. J Biol. 1996, 178 (7): 1842-1849.Google Scholar
- Ohashi Y, Munro S: Membrane delivery to the yeast autophagosome from the Golgi-endosomal system. Mol Biol Cell. 2011, 15; 21 (22): 3998-4008.Google Scholar
- Fazio A, Jewett MC, Daran-Lapujade P, Mustacchi R, Usaite R, Pronk JT, Workman CT, Nielsen J: Transcription factor control of growth rate dependent genes in Saccharomyces cerevisiae: A three factor design. BMC Genomics. 2008, 9: 341- 10.1186/1471-2164-9-341PubMed CentralView ArticlePubMedGoogle Scholar
- Kasahara T, Kasahara M: Identification of a key residue determining substrate affinity in the yeast glucose transporter Hxt7 - A two dimensional comprehensive study. J Biol Chem. 2010, 285: 26263-26268. 10.1074/jbc.M110.149716PubMed CentralView ArticlePubMedGoogle Scholar
- Geladé R, de Velde SV, Van Dijck P, Thevelein JM: Multi-level response of the yeast genome to glucose. Genome Biol. 2003, 4 (11): 233.33-View ArticleGoogle Scholar
- Brachmann CB, Davies A, Cost GJ, Caputo E, Li J, Hieter P, Boeke JD: Designer deletion strains derived from Saccharomyces cerevisiae S288C: a useful set of strains and plasmids for PCR-mediated gene disruption and other applications. Yeast. 1998, 14: 115-32. 10.1002/(SICI)1097-0061(19980130)14:2<115::AID-YEA204>3.0.CO;2-2View ArticlePubMedGoogle Scholar
- Baganz F, Hayes A, Marren D, Gardner DCJ, Oliver SG: Suitability of replacement markers for functional analysis studies in Saccharomyces cerevisiae. Yeast. 1997, 13: 1563-1573. 10.1002/(SICI)1097-0061(199712)13:16<1563::AID-YEA240>3.0.CO;2-6View ArticlePubMedGoogle Scholar
- Hayes A, Zhang N, Wu J, Butler PR, Hauser NC, Hoheisel JD, Lim FL, Sharrocks AD, Oliver SG: Hybridization array technology coupled with chemostat culture: Tools to interrogate gene expression in Saccharomyces cerevisiae. Methods. 2002, 26: 281-290. 10.1016/S1046-2023(02)00032-4View ArticlePubMedGoogle Scholar
- Wishart JA, Hayes A, Wardleworth L, Zhang N, Oliver SG: Doxycycline, the drug used to control the tet-regulatable promoter system, has no effect on global gene expression in Saccharomyces cerevisiae. Yeast. 2005, 22: 565-569. 10.1002/yea.1225View ArticlePubMedGoogle Scholar
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: Expression index computation and outlier detection. Proc Natl Acad Sci. 2001, 98: 31-36. 10.1073/pnas.011404098PubMed CentralView ArticlePubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on bias and variance. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185View ArticlePubMedGoogle Scholar
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FPC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME) - toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365View ArticlePubMedGoogle Scholar
- Storey JD, Xiao W, Leek JT, Tompkins RG, Davis RW: Significance analysis of time course microarray experiments. P Natl Acad Sci. 2005, 102: 12837-12842. 10.1073/pnas.0504609102.View ArticleGoogle Scholar
- Storey JD, Dai JY, Leek JT: The optimal discovery procedure for large-scale significance testing, with applications to comparative microarray experiments. Biostatistics. 2007, 8: 414-432.View ArticlePubMedGoogle Scholar
- Leek JT, Monsen EC, Dabney AR, Storey JD: EDGE: Extraction and analysis of differential gene expression. Bioinformatics. 2006, 22: 507-508. 10.1093/bioinformatics/btk005View ArticlePubMedGoogle Scholar
- Reich M, Ohm K, Tamayo P, Angelo M, Mesirov JP: GeneCluster 2.0: An advanced toolset for bioarray analysis. Bioinformatics. 2004, 20 (11): 1797-1798. 10.1093/bioinformatics/bth138View ArticlePubMedGoogle Scholar
- Seo J, Bakay M, Chen YW, Hilmer S, Shneiderman B, Hoffman EP: Interactively optimizing signal-to-noise ratios in expression profiling: project-specific algorithm selection and detection p-value weighting in Affymetrix microarrays. Bioinformatics. 2004, 20: 2534-2544. 10.1093/bioinformatics/bth280View ArticlePubMedGoogle Scholar
- Saccharomyces Genome Database GO Term Finder. http://db.yeastgenome.org/cgi-bin/GO/goTermFinder.pl
- Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S: The AmiGO Hub; the Web Presence Working Group. AmiGO: online access to ontology and annotation data. Bioinformatics. 2008,Google Scholar
- Matys V, Fricke1 E, Geffers R, Goßling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel1 AE, Kel-Margoulis OV, Kloos DU, Land S, Lewicki-Potapov B, Michael H, Munch R, Reuter I, Rotert S, Saxel H, Scheer M, Thiele S, Wingender E: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 1:31 (1): 374-378.View ArticlePubMedGoogle Scholar
- Teixeira MC, Monteiro P, Jain P, Tenreiro S, Fernandes AR, Mira NP, Alenquer M, Freitas AT, Oliveira AL, Sá-Correia I: The YEASTRACT database: a tool for the analysis of transcription regulatory associations in Saccharomyces cerevisiae. Nucleic Acids Res. 2006, 34: 446-451. 10.1093/nar/gkj013.View ArticleGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hanett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takgusakawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99-104. 10.1038/nature02800PubMed CentralView ArticlePubMedGoogle Scholar
- MacIsaak KD, Wang T, Gordon DB, Gifford DK, Stormo GD, Fraenkel E: An improved map of conserved regulatory sites for Saccharomyces cerevisiae. BMC Bioinformatics. 2006, 7: 113- 10.1186/1471-2105-7-113View ArticleGoogle Scholar
- Ulitsky I, Shamir R: Identification of functional modules using network topology and high-throughput data. BMC Syst Biol. 2007, 1: 8- 10.1186/1752-0509-1-8PubMed CentralView ArticlePubMedGoogle Scholar