Abstract
Saccharomyces cerevisiae SFP1 is required for nutrient-dependent regulation of ribosome biogenesis and cell size. A mutant deleted for SFP1 shows specific traits, including a slow growth phenotype, especially when growing on glucose. We recently analysed the physiology of an sfp1Δ mutant and its isogenic reference strain in chemostat cultures. This approach was successful in revealing the effects of nutrients on the activity of Sfp1 independent of growth rate-related feedback. In the present work we exposed carbon-limited cultures of an sfp1Δ mutant and its reference strain to sudden glucose excess. This allowed us to study the effect of SFP1 deletion on cell physiology when the cells are forced to exploit their maximum growth potential; this is similar to what happens in shake-flask cultures but with no bias due to growth rate differences. We show that nutrients differentiallly affect the role of Sfp1 in cell-size modulation and in transcriptional control. Furthermore, we report that while Sfp1 is necessary for the efficient glucose-dependent regulation of ribosome biogenesis genes, it is not required for the proper induction of ribosomal protein genes in response to glucose excess. Finally, our data suggest a role for Sfp1 in the regulation of glycolysis, further underlining its involvement in the network that links ribosome biogenesis and cell metabolism.
- CmM, mM carbon (content)
- FDR, false discovery rate
- MIPS, Munich Institute for Protein Sequences
- RP, ribosomal protein
- SAM, significance analysis of microarrays
-
↵†Present address: Department of Biotechnology, Delft University of Technology, Julianalaan 67, 2628 BC, Delft, The Netherlands.
-
The array data discussed in this publication have been deposited in the NCBI Gene Expression Omnibus (GEO; ) and are accessible through GEO series accession number GSE9644.
-
Three supplementary tables, listing genes showing significantly different expression in the sfp1 null mutant compared with the reference strain during ethanol/glucose-limited growth, significant differences in the transcriptional response between the reference and the mutant strain following the glucose pulse, and genes in particular clusters belonging to specified functional categories, and five supplementary figures, showing the cell-size distributions, extracellular glucose concentration and specific activity of glycolytic enzymes of the two strains, the Z score profiles of the 12 clusters showing differential expression, and modulation of the budding index and cell concentration after the glucose pulse, are available with the online version of this paper.
Edited by: H. A. B. Wösten
INTRODUCTION
The budding yeast Saccharomyces cerevisiae has evolved to survive sudden changes in the environment resulting from natural fluctuations in nutrient availability, temperature, osmotic pressure and other factors (Alberghina et al., 1998; Gasch & Werner-Washburne, 2002; Kresnowati et al., 2006). As a consequence, S. cerevisiae can rapidly adjust its physiology in response to most different external perturbations. A number of pathways exist in yeast which successfully sense external conditions and generate a signal that is transmitted from the cell surface to the target biosynthetic machineries (Crespo & Hall, 2002; Estruch, 2000; Rolland et al., 2002). Besides external stimuli, metabolic activity itself is an important player in the generation of such signals (Newcomb et al., 2003; Rolland et al., 2002).
In budding yeast, the TOR and PKA kinases are the major upstream regulators for the control of growth, metabolism and cell cycle progression in response to nutrients (Crespo & Hall, 2002; Rolland et al., 2002). One of the most important targets controlled by TOR and PKA is the ribosome biosynthesis machinery (Klein & Struhl, 1994; Martin et al., 2004; Powers & Walter, 1999; Zurita-Martinez & Cardenas, 2005). During optimal growth, a yeast cell spends most of its energy in making ribosomes. All three RNA polymerases are involved in this process and the activity of the RNA polymerase I, responsible for the transcription of the rRNA, accounts for up to 60 % of the transcriptional activity of a growing cell (Warner, 1999). Ribosomal protein (RP) mRNAs monopolize the cell translational capacity during fast proliferation, and a network of more than 100 factors (encoded by the so called RiBi gene cluster) is required for the correct assembly of the mature ribosomal particles (Fatica & Tollervey, 2002; Jorgensen et al., 2002). Such a massive utilization of cell resources explains why the control of ribosome biogenesis is a key issue for the economy of the cell. This is indeed finely tuned depending on the nutritional/environmental conditions (Martin et al., 2004; Zurita-Martinez & Cardenas, 2005).
The control of rRNA, RP and RiBi gene expression by TOR and PKA occurs via parallel and partially redundant pathways within a coordinated regulatory network (Crespo & Hall, 2002; Zurita-Martinez & Cardenas, 2005). Although the topology of this network has not yet been fully elucidated, in recent years some important downstream targets, such as Fhl1, Sch9, Hmo1 and Sfp1, have been discovered (Berger et al., 2007; Jorgensen et al., 2002; Martin et al., 2004; Urban et al., 2007).
Sfp1 is a zinc finger protein for which a role in the control of RiBi and RP gene expression in response to nutrient availability has been proposed (Fingerman et al., 2003; Jorgensen et al., 2004; Marion et al., 2004). The localization and activity of Sfp1 appear to be modulated by TOR via a PKA-dependent and a PKA-independent pathway (Marion et al., 2004). It has been suggested that a condition of respiro-fermentative growth on glucose stimulates Sfp1 activity (Cipollina et al., 2005). Despite extensive research, the mechanisms of action and regulation of Sfp1 remain elusive, and the identification of its specific transcriptional targets is still controversial (Cipollina et al., 2008; Fingerman et al., 2003; Jorgensen et al., 2004; Marion et al., 2004).
Besides its role in ribosome biogenesis, Sfp1 is involved in the modulation of cell-size setting; a mutant deleted for SFP1 is characterized by a carbon source-dependent small size (whi) phenotype (Cipollina et al., 2005; Jorgensen et al., 2002). This dual role has led several researchers to suspect Sfp1 to be at the interface between ribosome biogenesis and condition-dependent cell-size control. Whether the setting of cell size is a readout of the growth events or is dictated upstream by ribosome biogenesis still remains an open question (Jorgensen et al., 2002; Schneider et al., 2004).
Most of the current evidence on the possible function of SFP1 has so far been derived from the comparison of sfp1Δ mutants with their isogenic reference strains growing in shake flasks. However, the mutant strain has a slow growth phenotype which hinders the separation of the direct effects due to SFP1 deletion from secondary growth rate-related effects (Cipollina et al., 2005; Fingerman et al., 2003; Jorgensen et al., 2004; Marion et al., 2004).
We recently grew an sfp1Δ strain and its isogenic reference strain at the same specific growth rate by using chemostat cultures (Hoskisson & Hobbs, 2005; Novick & Szilard, 1950). Chemostat cultivation has proved to be a powerful tool to study the involvement of Sfp1 in cell-size modulation and transcriptional control in response to nutrients, independent of the growth rate-related phenotype (Cipollina et al., 2008).
Nevertheless, this chemostat-based approach suffers from several shortcomings. For instance, to maintain steady growth in chemostats, the growth rates are slower than the maximum specific growth rates (μmax) of the strains. In addition, the cultures were run under carbon limitation with slow feeding rates, leading to relatively slow in vivo carbon fluxes. It could therefore be surmised that the effect of the deletion of SFP1 on S. cerevisiae physiology would be weaker under these conditions than when the cells have to exploit their full potential for fast growth and metabolism (as in exponential growth in shake-flask cultivation). A new set of experiments was therefore designed to circumvent these drawbacks while preventing growth-rate bias.
In the present work, sfp1Δ and its reference strain cultivated at the same growth rate in aerobic carbon-limited chemostats were exposed to a sudden increase in glucose concentration coupled with a switch to anaerobiosis. Using this experimental set-up we explored the role of Sfp1 when cells have to sense the change in nutrient abundance and rapidly adjust their metabolism and growth to this new condition, exploiting their maximum growth potential. The study of the early response following glucose addition allowed us to quantify the effect of SFP1 deletion on cell physiology, cell-size control and transcriptional regulation, with no bias due to the slow-growth phenotype. Our results shed new light on the role of Sfp1 at the interface between nutrient sensing, ribosome biogenesis, size setting and carbon metabolism.
METHODS
Strains and growth conditions.
The prototrophic reference S. cerevisiae strain CEN.PK113.7D (MATa MAL2-8c SUC2) and the isogenic strain CEN.PK111.32D sfp1Δ (MATa MAL2-8c SUC2 leu2-3,112 sfp1Δ : : KlLEU2) were used in this study. Cells were grown at 30 °C in 2 l chemostats (Applikon), with a working volume of 1.5 l. The medium composition was based on that described by Verduyn et al. (1992). Cultures were fed with a defined mineral medium that limited growth on a mixture of glucose and ethanol. The substrate concentrations were 200 and 50 mM carbon (content) (CmM) for glucose and ethanol, respectively. The addition of ethanol successfully prevented the sfp1Δ cultures from oscillating (Mashego et al., 2007). The dilution rate (D) was set at 0.10 h−1 and the airflow at 0.5 l min−1. The stirrer speed was 800 r.p.m. The volume and pH control (pH 5) and off-gas analysis were performed as previously described (van Maris et al., 2003). Chemostat cultures were assumed to be in a steady-state when, after at least five volume changes, the culture dry-weight, specific CO2 production rate and oxygen consumption rate changed by less than 2 % during 24 h. Steady-state samples were taken from cultures after 10 volume changes. We will refer to the moment right before the onset of the pulse as time 0 (t=0).
Glucose anaerobic pulse.
After 10 volume changes, a few seconds after the samples for the steady-state analysis were collected, the anaerobic glucose pulse experiments were started by sparging the medium reservoir and the fermenter with pure nitrogen gas (airflow 0.5 l min−1, Hoek-Loos, <5 p.p.m. O2). Norprene tubing and butyl septa were used to minimize oxygen diffusion into the anaerobic culture. Two minutes after nitrogen sparging and just before the addition of glucose, the medium and the effluent pumps were switched off. At this time point, 60 ml 5.0 M glucose was injected aseptically through a rubber septum giving a glucose concentration of ∼190 mM in the fermenter. CO2 formation was monitored online by continuously analysing the off-gas composition. Samples were collected at different time points depending on the specific analysis and the strain.
Analytical methods.
Extracellular metabolite concentration and culture dry weight were determined as described previously (Cipollina et al., 2008; Postma et al., 1989). The specific consumption/production rates for glucose and ethanol were calculated by first fitting the concentration profile of the given metabolite to a second-order polynomial. The derivative of the polynomial was then calculated for each time point and divided by the biomass concentration measured at that time point.
In vitro measurements of enzymatic activities.
The enzyme assays were performed on fresh cell-free extracts. A 60 ml volume of culture broth was sampled from the fermenter and washed twice in ice-cold freeze buffer (10 mM potassium phosphate, 2 mM EDTA, pH 7.5). The pellet was resuspended in freeze buffer (10 ml) and stored at −20 °C in 2 ml aliquots. Frozen samples were thawed at room temperature, washed with sonication buffer (100 mM potassium phosphate, 2 mM MgCl2, pH 7.5) and centrifuged (5 min at 4600 g, 4 °C). The pellet was resuspended in 2 ml ice-cold sonication buffer to which DTT was added to a final concentration of 1 mM. The suspension was added to 1.5 ml Safe-Lock tubes (Eppendorf) filled beforehand with 0.75 g cold glass beads (Sigma). The tubes were placed in a FastPrep120 machine (Qbiogene) and shaken in four bursts of 20 s, at speed 6, for the efficient breakdown of the cell membranes; between the shaking cycles the tubes were cooled on ice for at least 60 s. Unbroken cells and debris were removed by centrifugation (20 min at 47 000 g, 4 °C). The cell-free extracts obtained were used to perform the enzyme assays as described previously (Jansen et al., 2005).
Absorbance measurements were performed with a Tecan GENios Pro microtitre plate reader. All determinations were performed at 30 °C and 340 nm (ϵNAD(P)H=6.3 mM−1 cm−1). The assays were performed in duplicate and with two different dilutions of cell extracts. Protein concentration in the cell-free extracts was measured by the Lowry method.
Microarray analysis.
Sampling from chemostats, total RNA extraction, probe preparation and hybridization to Affymetrix Genechip microarrays were performed as described previously (Cipollina et al., 2008). Samples were collected at steady state and then at 5, 10, 30, 60 and 120 min after the pulse. The results related to steady-state samples were derived from three independent cultures; those related to the time-course analysis were derived from two independent cultures.
Acquisition and quantification of array images, global scaling and data filtering were performed using Affymetrix GeneChip Operating Software (GCOS) version 1.2, as described previously (Cipollina et al., 2008; Piper et al., 2002). To represent the variation in triplicate and duplicate measurements, the coefficient of variation (CV) was calculated as the standard deviation divided by the mean (Boer et al., 2003). For all the replicates the CV varied by less than 13 %, indicating a good reproducibility of the results. After global scaling (Piper et al., 2002), the expression of the housekeeping genes ACT1, HHT2 and SHR3 (encoding, respectively, actin, histone and an endoplasmic reticulum packaging chaperone) remained constant throughout the experiment, showing a CV below or around 30 %.
For statistical analyses of steady-state data, Microsoft Excel running the Significance Analysis of Microarrays (SAM, version 2.23) add-in was used for pair-wise comparisons between the reference and the sfp1Δ mutant datasets (Boer et al., 2003; Tusher et al., 2001). Enrichment for specific Munich Institute for Protein Sequences (MIPS) functional categories among the differentially expressed genes compared with their genome-wide occurrence was computed by considering a hypergeometric distribution. The obtained P-value indicated the statistical significance of the results.
The statistical analysis of dynamic data was performed with the Extraction of Differential Gene Expression (EDGE) software (Leek et al., 2006). We compared the time-course transcriptional profile between the sfp1Δ mutant and the reference strain (EDGE command: ‘Identify Differentially Expressed Genes – between class’), and 1524 genes were classified as significant (q value cutoff: 0.005). This list of genes was then used as input for the k-means clustering. Each gene was represented by a vector of 12 elements: elements 1–6 represented the expression levels of the gene in the reference strain; elements 7–12 were the expression levels of the same gene in the mutant. The expression levels were normalized so that they had zero mean and unit variance. The coefficient of correlation was used as measure of the distance between vectors. The k-means determination was performed so as to allow the k value (i.e. the number of clusters) to range from 2 up to 30 (in steps of two). A k value of 12, i.e. 12 clusters, was chosen as the one giving the best compromise between good partitioning of the dataset (the sd of each cluster was used as measure for this) and enrichment for MIPS functional categories.
To have a quantitative measure of the significance of the expression changes for a given cluster, the Z score was computed for each time point using the following equation (Knijnenburg et al., 2008):
Quantification of RNA contents and flow cytometric analysis.
Samples were collected at steady state and at the selected time points during the pulse. RNA content was chemically determined by the orcinol method (Popolo et al., 1982). The flow cytometric analysis for protein determination was performed on FITC-stained cells, as described by Vanoni et al. (1983).
RESULTS
Physiology and global transcriptional profile of an sfp1Δ mutant and the isogenic CEN.PK 113-7D strain growing under glucose/ethanol limitation
The deletion of SFP1 in S. cerevisiae causes a slow growth phenotype during growth in shake flasks (Cipollina et al., 2005; Jorgensen et al., 2002). To eliminate growth rate difference between CEN.PK 113-7D and its isogenic sfp1Δ mutant, both strains were grown at the same growth rate of 0.10 h−1 (duplication time of 6.9 h) in chemostat cultures limited by a glucose–ethanol mixture (Cipollina et al., 2008; Novick & Szilard, 1950). As shown in Table 1⇓, under this growth regime both strains showed highly similar physiological characteristics. The very low residual concentrations of glucose and ethanol in the broth culture indicated that both substrates were being co-consumed.
Physiology of CEN.PK 113.7D and the isogenic sfp1Δ mutant during continuous growth at D=0.10 h−1, limited by glucose/ethanol mixed substrate
Data represent the mean and sd from four independent experiments. Res. Glu/Res. EtOH=residual concentration of glucose/ethanol in the culture broth; RQ, respiratory coefficient (
One of the most typical features of the sfp1Δ phenotype is a reduction in both cell size and RNA content, especially during batch growth on glucose (Cipollina et al., 2005; Jorgensen et al., 2002). The protein content, which provides a good estimate of the cell size (Alberghina & Porro, 1993; Vanoni et al., 1983), was analysed by flow cytometry. The mutant population had a clear small size (whi) phenotype (Fig. 1⇓, at t=0), with smaller cells both at birth and at division (Supplementary Fig. S1). The analysis of the volume distributions, performed with a Coulter counter analyser, showed results very similar to the ones obtained by the flow cytometric approach (data not shown). The total RNA content of the cells was also significantly lower in the sfp1Δ mutant compared with the reference strain (Fig. 2⇓, at t=0).
Protein content per cell at steady state (t=0, dashed line) and after the glucose pulse. Samples were collected at the indicated time points. The average protein content per cell [expressed as channel (ch.) number, relative units] was determined from the protein content distribution obtained by flow cytometry on FITC-stained cells. Data are shown for the reference (•) and the sfp1Δ mutant strain (○). The average and the mean deviation of two independent experiments are reported.
Changes in the RNA content per cell in response to the glucose pulse. Samples were collected at steady state (t=0, dashed line) and at the indicated time points after the pulse. The total RNA content was measured with the orcinol method. Data are reported for the reference (•) and the sfp1Δ mutant strain (○). The average and the mean deviation of two independent experiments are reported.
The transcriptional response of yeast to the deletion of SFP1 was investigated by microarray analysis. Statistical analysis of the array data for the sfp1Δ and reference strains was followed by an enrichment computation of specific functional categories among the significantly differentially expressed genes (see Methods). Only the use of lenient statistical criteria [i.e. significance analysis of microarrays (SAM) analysis, false discovery rate (FDR)=1 %, fold change=1.5; Supplementary Table S1] led to the identification of enriched categories. While the genes with a higher expression in the mutant were marginally enriched for the ‘pheromone response, mating-type determination, sex-specific proteins’ category, some of the genes displaying significant lower expression in the mutant could be more specifically related to ribosome biogenesis (RPA14, RPB5, RPA190, RPA135 and RRN3). These genes encoded subunits of the RNA polymerase I, the enzymic complex responsible for the transcription of rRNA (Moss, 2004). Since rRNA represents ∼80 % of the total cellular RNA (Warner, 1999), this finding was consistent with the lower RNA content typical of sfp1Δ mutant cells (Fig. 2⇑).
Adaptation of S. cerevisiae CEN.PK 113-7D and sfp1Δ strains to fast metabolic rates
To investigate the response of the sfp1Δ mutant to a sudden increase in the external glucose concentration, the carbon-limited chemostat cultures of both mutant and reference strains were pulsed with a highly concentrated glucose solution (final concentration ∼190 mM). Since the sfp1Δ phenotype is especially evident during anaerobic growth and hence under fully fermentative growth conditions (Cipollina et al., 2005, 2008), we decided to switch the cultures to anaerobiosis at the time of the glucose pulse. Moreover, we have recently reported an extensive analysis of the response of S. cerevisiae CEN.PK 113-7D to an anaerobic glucose pulse (van den Brink et al., 2008). Here, by choosing a similar experimental set-up, we could specifically focus on the phenotypic traits due to SFP1 deletion, leaving the detailed description of the behaviour of the reference strains to the previous manuscript.
After glucose addition, both strains, which displayed a fully respiratory metabolism in carbon-limited cultures, immediately switched to fermentative metabolism, as revealed by ethanol production (Fig. 3a⇓). Similarly, the specific glucose consumption rate (qglucose) immediately increased for both strains (Fig. 3b⇓). After the initial rapid induction, qglucose increased more slowly until it reached a plateau ∼2.5 h after the glucose pulse. The magnitude of qglucose throughout the pulse was substantially lower for the sfp1Δ mutant than for the reference strain (about twofold lower after 2.5 h). This resulted in a longer fermentation time for the mutant (Supplementary Fig. S2).
Time-course profile of ethanol specific production rate (a) and glucose specific consumption rate (b) at steady state and after the glucose pulse (t=0, dotted line). The q values are reported for the time window in which the glucose concentration remained above 50 mM (i.e. 3 h for the reference strain and 5 h for the mutant; see Supplementary Fig. S2). Data are reported for the reference (•) and the sfp1Δ mutant strain (○). Data represent the average and the mean deviation of four independent experiments.
Although sfp1Δ could increase its ethanol production rate, and therefore its glycolytic flux, the deletion of SFP1 clearly resulted in a strain with lower maximum glycolytic rates compared with the reference strain. One likely explanation of the inability of sfp1Δ to channel high carbon fluxes could be a lower glucose-uptake capacity. In yeast, the hexose transporters are mainly encoded by seven genes (HXT1–HXT7), whose expression is mostly (although not exclusively) regulated at the transcriptional level (Ozcan et al., 1996). The transcriptional response of HXT genes in the sfp1Δ mutant was very similar to that of the reference strain, apart from a slightly lower induction of HXT3 and HXT4 (data not shown), suggesting that glucose uptake was not responsible for the slower glycolytic flux in the sfp1Δ mutant. Alternatively, this lower glycolytic flux could result from a lower maximum capacity (Vmax) of the glycolytic enzymes. This Vmax can be quantified in vitro by assaying the specific activity of the glycolytic enzymes, which is expressed as enzyme units per milligram of protein. The Vmax of most glycolytic enzymes increased between 30 and 180 min following glucose addition (Supplementary Fig. S3), and both their profiles and the intensity of the change were very similar for both strains. The changes in the transcript levels of these enzymes were unaffected in the mutant strain, and mirrored rather well the changes of the activity profiles (data not shown).
Conversely, the capacities of the first and the last enzymes of the glycolytic pathway, i.e. hexokinase (Hxk) and alcohol dehydrogenase (Adh), were altered in the mutant phenotype compared with the reference strain (Fig. 4⇓). The specific activity of Hxk at steady state and throughout the pulse experiment was lower (∼1.5-fold) in the sfp1Δ strain. While HXK1 and GLK1 showed similar expression levels and regulation in response to glucose, HXK2 (encoding the predominant hexokinase isoform during growth on glucose; Johnston, 1999) showed, in good agreement with Hxk Vmax, a steady twofold lower expression in the mutant strain (Fig. 4⇓). This suggested that the difference in Hxk capacity between the reference and mutant strains was predominantly regulated at the transcriptional level.
Transcript levels and specific activity of hexokinase (a), alcohol dehydrogenase (b) and enolase (c) at steady state (t=0, dotted line) and after the glucose pulse. Transcript data are reported in the heat maps and derive from samples collected at steady state and at 5, 10, 30, 60 and 120 min after glucose addition. The transcript levels were normalized by dividing the expression level of each gene at a given time point by the mean expression level over the 12 conditions (six time points, two strains). The key to the heat maps is shown to the right of panel (c). Samples for enzyme activity determination were collected at steady state and 30, 60, 120 and 180 min after the addition of glucose. Data are reported for the reference (•) and the sfp1Δ strain (○), and represent the average and the mean deviation of two independent experiments.
The activity profile for Adh was similar to that for Hxk, with lower values in the mutant strain (∼1.7-fold) throughout the experiment. The correlation between activity and transcript levels for this enzyme was not straightforward, since in yeast six genes encode Adh. Fig. 4⇑ shows that while ADH2 was similarly downregulated in both strains after glucose addition, the regulation of the remaining ADH genes was affected in the mutant strain. This might correlate with the observed differences in enzyme activity.
Finally, enolase (Eno) activity, although increasing after glucose addition for both strains, was slightly lower in the sfp1Δ mutant after the pulse. Such a difference could not be predicted from the mRNA levels (Fig. 4⇑) and suggested a post-transcriptional regulation of the enzyme capacity (i.e. by protein synthesis and turnover or post-translational modification).
Effect of SFP1 deletion on the global transcriptional response to the glucose pulse
We have recently investigated the general transcriptional response of budding yeast to an anaerobic glucose pulse (van den Brink et al., 2008). Here we specifically studied the effect of SFP1 deletion on the transcriptional reprogramming triggered by this perturbation. Statistical analysis revealed that the expression profiles of 1524 genes differed significantly between the reference and the sfp1Δ strain after the anaerobic glucose pulse (Supplementary Table S2). These 1524 genes were spread across twelve clusters according to their expression profiles (Supplementary Fig. S4). Each cluster was subsequently searched for enrichment in specific MIPS functional categories and promoter elements.
As reported in Table 2⇓, significant enrichment was found for most of the 12 clusters. However, the biological meaning of such enrichment was in many cases unclear. For example, clusters 5, 9, 11 and 12 were enriched for categories related to carbon metabolism and cell cycle progression, and cluster 6 comprised a significant number of genes involved in the stress response. However, the number of genes represented in each cluster, although statistically significant, was too low compared with the total number of genes in the category to allow meaningful biological conclusions to be drawn. Neither clear patterns nor relationships were found between these genes and the observed phenotype, and no striking differences in expression profiles between the two strains were observed for these genes. Clusters 1 and 7 showed no enrichment for MIPS categories, while cluster 3 showed some enrichment for a small category related to vitamins/cofactors transport. Finally, clusters 2, 4, 8 and 10 showed the strongest differences in the expression profiles between the two strains and the most significant enrichments for functional categories.
Enrichment for MIPS functional categories and known transcription factor (TF) binding sites
Horizontal lines within the body of the table demarcate the MIPS category combinations that add up to the number of genes in each cluster.
Genes belonging to cluster 2 showed a comparable initial induction in the two strains until 30 min, but then diverged in the two strains during the remaining fermentation (Fig. 5⇓). This cluster was enriched for genes involved in amino acid metabolism and translation elongation, and for ribosomal protein genes, and accordingly promoter elements recognized by the transcription factors Fhl1 and Rap1 were over-represented (Table 2⇑; Knijnenburg et al., 2007). Genes grouped in clusters 4 and 10 showed a delay in induction in the sfp1Δ strain as compared with the reference strain. These clusters comprised genes involved in ribosome biogenesis, translation initiation and nucleotide metabolism, and showed accordingly a significant enrichment for RRPE and PAC promoter elements (Table 2⇑) (Hughes et al., 2000).
Z score profiles of clusters (Cl) 2, 4, 8 and 10 for the reference (solid line) and the sfp1Δ (dashed line) strain. The Z score indicates both the significance of the change at time tx relative to t=0 and the direction of the change (up-/downregulation). It was calculated as described in Methods.
Finally, cluster 8 displayed the opposite trend between the reference and the sfp1Δ strain, and was enriched for genes involved in proteasomal degradation. Thirteen out of 106 genes in this cluster encoded either structural or regulatory subunits of the 26S proteasome involved in the degradation of ubiquitinated substrates. The binding site recognized by Rpn4, a transcription factor specific for proteasome genes (Xie & Varshavsky, 2001), was accordingly enriched in cluster 8 (Table 2⇑). The proteasome 26S is involved, among other functions, in the post-transcriptional regulation of many important factors required for proper control of cell cycle progression (DeSalle & Pagano, 2001; Futcher, 1996). The altered regulation of these genes in the sfp1Δ strain might therefore be related to its small size phenotype. However, in the present work, this hypothesis was not investigated further.
Effect of SFP1 deletion on RiBi and RP gene expression in response to glucose pulse
We have previously shown that the deletion of SFP1 in S. cerevisiae resulted in decreased expression of the genes involved in ribosome biogenesis (the so called ‘RiBi cluster’), but did not affect expression of RP genes (Cipollina et al., 2008). These results suggested that Sfp1 was not involved in the transcriptional regulation of RP genes when S. cerevisiae was grown at a slow growth rate under limited carbon supply. Here, the addition of a glucose excess to carbon-limited cultures allowed us to analyse the involvement of Sfp1 in RiBi and RP gene transcriptional modulation when the cells express their full potential for sugar utilization and growth.
In agreement with our recently published data (van den Brink et al., 2008), RiBi (clusters 4 and 10) and RP (cluster 2) genes responded with different dynamics to sudden glucose excess (Fig. 6⇓). We focused on the 96 rRNA processing genes belonging to cluster 10 and the 47 RP genes belonging to cluster 2. While the 47 RP genes were similarly induced in both strains in the 30 min following glucose addition (Fig. 6a⇓), the induction of the RiBi genes showed a substantial delay in the mutant strain and remained weaker throughout the whole fermentation (Fig. 6b⇓). This finding indicated that while SFP1 was not essential for the early induction of RP genes, it was required for the initial transcriptional burst of the RiBi genes. After 30 min, RP genes were downregulated in the mutant strain in accordance with previous reports showing that during shake-flask growth the expression of these genes is reduced in a sfp1Δ strain. Interestingly, genes encoding factors required for nucleotide metabolism and translation initiation showed a response similar to the RiBi genes, while genes encoding factors required for amino acid biosynthesis and translation elongation showed a response similar to that of the RP genes (Table 2⇑).
Z score profiles of RiBi and RP genes in response to the glucose pulse. The Z score profiles relating to the 47 RP genes from cluster 2 (a) and the 96 ‘rRNA processing’ genes belonging to cluster 10 (b) are reported for the reference (solid line) and the sfp1Δ (dashed line) strain. The Z score indicates both the significance of the change at time tx relative to t=0 and the direction of the change (up-/downregulation). It was calculated as described in Methods.
The observed defect in the induction of RiBi genes was consistent with the smaller increase in the total RNA content in the sfp1Δ mutant compared with the reference strain (Fig. 2⇑).
Modulation of cell size and cell cycle progression
One of the typical traits of an sfp1Δ mutant is the small size (also called whi) phenotype (Cipollina et al., 2005; Jorgensen et al., 2002). Accordingly, during aerobic ethanol/glucose-limited growth, the sfp1Δ strain showed a whi phenotype (Fig. 1⇑). Approximately 1 h after the glucose pulse, the cell size of the reference strain increased to reach a value 30 % higher than the steady-state value at the end of the fermentation. Conversely, the size of the mutant cells remained steady for about 2 h after the glucose pulse, then finally slowly increased to reach a value as little as 10 % higher than the steady-state size (Fig. 1⇑). This resulted in an increase in the difference in cell size between the two strains caused by the addition of glucose excess, consistent with previous reports (Cipollina et al., 2005).
It has been shown elsewhere that after the addition of excess glucose to non-repressed cells, a yeast culture starts to grow (i.e. produce biomass) at its maximum specific growth rate (μmax) (Alberghina et al., 1998). Another typical effect induced by a glucose pulse is the increase in cell size, which is attained through a tight coordination between growth and cell cycle progression (Alberghina et al., 1998; Flick et al., 1998). During anaerobic batch growth on glucose, the reference and mutant strains grew with μmax values of 0.13 and 0.34 h−1, respectively. After the pulse, we observed a similar behaviour for both strains, with a transient decrease of the budding index (BI) from ∼40 to 30 % (1 h after the pulse) followed by an increase to about 60 % (Supplementary Fig. S5a). This could be explained by a temporary delay at the G1 to S transition (Alberghina et al., 1998; Flick et al., 1998). Accordingly, during the first 2 h after the pulse, no significant increase in cell number was observed for both strains, indicating that the cellular proliferation had temporarily slowed down (Supplementary Fig. S5b). The combined effect of this delay in the G1 to S transition and the change in the specific growth rate resulted in the observed adjustment to the new size setting, which was different for the two strains (Fig. 1⇑).
DISCUSSION
The involvement of Sfp1 in cell size modulation and transcriptional control is differentially modulated by the growth conditions
The use of chemostat cultures has recently been proven to be a very good approach to study the role of Sfp1 in cell-size modulation and transcriptional control (Cipollina et al., 2008). It allowed us to grow both a mutant deleted for SFP1 and its reference isogenic strain at the same growth rate under different carbon limitations. This would not have been possible in shake-flask cultures due to the slow growth phenotype of the sfp1Δ mutant (Cipollina et al., 2005; Jorgensen et al., 2002).
The present work concludes a series of chemostat experiments performed to assess the influence of the carbon source on sfp1Δ phenotype, independent of the growth rate-related phenotype. Here, we analysed the physiology of cultures growing aerobically under glucose/ethanol limitation. This growth regime created a metabolic condition intermediate between the previously analysed conditions (Cipollina et al., 2008). During anaerobic glucose-limited growth, the substrate was completely fermented to ethanol and the glycolytic flux was high [qglucose approx. −3.2 mmol (g biomass)−1 h−1] (Cipollina et al., 2008; Daran-Lapujade et al., 2007). During ethanol-limited growth, the cultures showed a purely respiratory metabolism and the glycolytic flux proceeded backwards (gluconeogenesis) for biosynthetic purposes (de Jong-Gubbels et al., 1995). In the present study, during aerobic glucose/ethanol limitation, with glucose supplying 80 % of the substrate carbon, metabolism was purely respiratory and glycolysis fully active, although the rate of glucose consumption was lower compared with the anaerobic growth previously analysed [qglucose approx. −0.90 mmol (g biomass)−1 h−1]. The finding that the sfp1Δ mutant showed a clear whi phenotype in all the growth regimes so far analysed indicates that the involvement of SFP1 in the modulation of cell size is independent of the growth rate, the supplied carbon source and the metabolism of the cell.
Conversely, the involvement of Sfp1 in RiBi transcriptional regulation appears to be strongly related to the carbon source provided and the growth regime. We observed that when no glucose was present in the medium (ethanol-limited growth), the transcriptional profile of an sfp1Δ mutant was similar to that of the reference strain. During anaerobic glucose-limited growth, serious transcriptional alterations involving the RiBi cluster and the RNA polymerase I coding genes were observed for the mutant strain (Cipollina et al., 2008). Here, under the intermediate metabolic regime typical of glucose/ethanol-limited growth, some transcriptional alteration was present, related to the small but important cluster of RNA Pol I coding genes. These findings suggest that a sustained glycolytic flux and/or a fermentative metabolism might generate a signal stimulating the Sfp1 transcriptional role. It has been reported previously that the presence of an active glycolytic flux contributes to the activation of important nutrient-signalling pathways that regulate growth and metabolism (Li et al., 2006; Newcomb et al., 2003; Rolland et al., 2001, 2002). Sfp1 might be involved in this regulatory network.
Sfp1 is required for efficient induction of RiBi, but not RP, genes in response to glucose excess
The involvement of Sfp1 in transcriptional regulation of RiBi and RP genes has so far been a controversial issue. Most of the data derived from shake-flask cultures indicate that during growth on glucose, Sfp1 is involved in transcriptional regulation of both RiBi and RP genes (Jorgensen et al., 2004; Marion et al., 2004). The former gene cluster is characterized by the presence of RRPE and PAC elements in the promoter regions. The latter is particularly enriched for the promoter elements recognized by the transcription factors Rap1 and Fhl1. All attempts to prove a physical interaction between Sfp1 and the RiBi gene promoters have failed, while some weak interactions have been observed between Sfp1 and the promoters of some RP genes (Fingerman et al., 2003; Harbison et al., 2004; Jorgensen et al., 2004; Marion et al., 2004). Recently it has been shown that the protein Stb3 binds the RRPE elements in vivo and that it is required for the proper regulation of RRPE-containing genes in response to glucose. It has been suggested that Sfp1 acts upstream of Stb3 in modulating the transcription of RiBi genes (Liko et al., 2007).
We recently proposed that the main target of Sfp1 is the RiBi cluster, and that the downregulation of RP genes observed in the sf1pΔ mutant during shake-flask growth might be a secondary effect due to the slow growth rate of the mutant strain. The glucose-pulse experiment presented here was successful in further supporting this hypothesis. Here we conclusively show that Sfp1 is required for the efficient and rapid increase in expression of RiBi genes in response to the glucose pulse, but not for the early transcriptional induction of RP genes (Fig. 6⇑). Since the transcriptional regulation of RP and RiBi genes has been shown to be tightly co-ordinated (Griffioen et al., 1996; Laferte et al., 2006), we suggest that the downregulation of RP genes observed in the mutant strain from 30 min after glucose addition is due to the adjustment of RP expression levels to the low expression of RiBi and RNA Pol I genes (Fig. 6⇑). This would in turn result in the slow growth-rate phenotype typical of sfp1Δ cells.
The observation that RiBi genes are still able to respond to glucose addition in the sfp1Δ strain, although with some delay and to a lesser extent compared with the reference strain, indicates that Sfp1 is not the only regulator of these genes. In the absence of SFP1, alternative mechanisms appear to operate to modulate the expression of RiBi genes depending on nutrient availability. Finally, while the observation of a weak binding of Sfp1 to RP promoter genes might suggest a role in the regulation of these genes (Fingerman et al., 2003; Jorgensen et al., 2004; Marion et al., 2004), our data indicate that this is not a primary function of Sfp1.
Sfp1 might play a role at the crossroads of glycolysis regulation and ribosome biogenesis
Most of the information so far available indicates that the phenotype of an sfp1Δ mutant is much more evident during fermentative growth on glucose than during respiratory growth on a non-fermentable carbon source (ethanol or glycerol). In the present work we chose an experimental set-up that allowed us to analyse the transition from a condition of respiratory growth in which the mutant phenotype was only partly present to a condition of anaerobic batch growth on glucose (fully fermentative). Therefore, we expected the mutant phenotype to dominate during the time interval that followed the perturbation. In particular, since at steady state both strains were growing at the same specific growth rate, the observed alterations in the early response to glucose addition could be ascribed to the absence of SFP1 rather than to growth rate-related defects.
Besides a different transcriptional phenotype, we observed significant defects in the induction of the glucose-specific consumption rate in the sfp1Δ mutant within minutes of glucose addition. This clearly indicated that deletion of SFP1 has a negative effect on the glycolytic activity of the cell. As reported in Results, the Vmax of most glycolytic enzymes at steady state and in response to glucose showed no difference between the two strains. The reported enzymic activities were normalized based on the protein content of the cell cultures used for the assays (see Methods). Since the sfp1Δ strain had a lower protein content than the reference strain at steady state and throughout the pulse experiment (Fig. 1⇑, Table 1⇑), the absolute capacity of glycolytic enzymes in the mutant cultures was lower than that of the reference strain. The effect of such lower activities was not evident at steady state, when the glucose consumption rate was relatively low (Table 1⇑). However, after the addition of excess glucose the cells were forced to exploit their maximum glycolytic potential. Under this condition the lower glycolytic potential of the sfp1Δ cultures could partly explain the observed reduction in the glycolytic flux. Nevertheless, it is noteworthy that the deletion of SFP1 did not affect the induction (relative to t=0) and thus the regulation of most glycolytic enzymes in response to glucose addition (Supplementary Fig. S3).
Conversely, the enzyme hexokinase showed differential expression between the two strains, with both lower transcriptional levels and lower activity in the sfp1Δ mutant throughout the experiment. This enzyme plays a key role in the regulation of central carbon metabolism for two reasons (Westergaard et al., 2007). It catalyses the first step of glycolysis, i.e. the phosphorylation of glucose, and it is involved in the regulation of protein kinase A and Mig1, which control growth and central metabolism in response to nutrient availability and environmental conditions (Ahuatzi et al., 2007; Moreno et al., 2005; Rolland et al., 2001). Given the crucial role of Hxk in the regulation of glycolytic flux, the reduced Vmax observed for this enzyme could be related to the lower glycolytic rate of the mutant strain. In addition, the observed defects in Hxk2 expression and activity may be linked to the carbon source-dependent phenotype of the sfp1Δ mutant. Whether Sfp1 affects Hxk activity via a direct mechanism or lower Hxk activity is caused by the low RNA or protein content typical of the mutant strain remains to be elucidated.
The finding that the deletion of SFP1 simultaneously affects the induction of both RiBi gene transcription and glycolytic activity in response to glucose suggests that Sfp1 is an upstream regulator that plays a role at the crossroads of ribosome biogenesis and central carbon metabolism.
Acknowledgments
This work was supported by the IOP-Genomics programme of Senter/Novem, and by the Netherlands Genomics Initiative (NGI) of the Netherlands Science Council (NWO). The authors wish to thank Marinka Almering and Erwin Suir for expert technical assistance, and Theo Knijnenburg for his expert contribution on microarray data mining. Marina Vai is specially thanked for thorough scientific discussions.