|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 10, 9125-9138, March 5, 2004
Role of Transcriptional Regulation in Controlling Fluxes in Central Carbon Metabolism of Saccharomyces cerevisiaeA CHEMOSTAT CULTURE STUDY*![]() ![]() ![]() ![]() ![]() ¶![]()
From the
Received for publication, August 28, 2003 , and in revised form, November 20, 2003.
In contrast to batch cultivation, chemostat cultivation allows the identification of carbon source responses without interference by carbon-catabolite repression, accumulation of toxic products, and differences in specific growth rate. This study focuses on the yeast Saccharomyces cerevisiae, grown in aerobic, carbon-limited chemostat cultures. Genome-wide transcript levels and in vivo fluxes were compared for growth on two sugars, glucose and maltose, and for two C2-compounds, ethanol and acetate. In contrast to previous reports on batch cultures, few genes (180 genes) responded to changes of the carbon source by a changed transcript level. Very few transcript levels were changed when glucose as the growth-limiting nutrient was compared with maltose (33 transcripts), or when acetate was compared with ethanol (16 transcripts). Although metabolic flux analysis using a stoichiometric model revealed major changes in the central carbon metabolism, only 117 genes exhibited a significantly different transcript level when sugars and C2-compounds were provided as the growth-limiting nutrient. Despite the extensive knowledge on carbon source regulation in yeast, many of the carbon source-responsive genes encoded proteins with unknown or incompletely characterized biological functions. In silico promoter analysis of carbon source-responsive genes confirmed the involvement of several known transcriptional regulators and suggested the involvement of additional regulators. Transcripts involved in the glyoxylate cycle and gluconeogenesis showed a good correlation with in vivo fluxes. This correlation was, however, not observed for other important pathways, including the pentose-phosphate pathway, tricarboxylic acid cycle, and, in particular, glycolysis. These results indicate that in vivo fluxes in the central carbon metabolism of S. cerevisiae grown in steadystate, carbon-limited chemostat cultures are controlled to a large extent via post-transcriptional mechanisms.
The yeast Saccharomyces cerevisiae is widely used as a model organism to study carbon source-dependent metabolic regulation in eukaryotes. Wild-type S. cerevisiae strains have a narrow set of carbon sources that can support fast growth in synthetic media (1). The most widely studied of these are the hexoses glucose, fructose, galactose, and mannose, the disaccharides maltose and sucrose, and the C2-compounds ethanol and acetate. The metabolic networks employed for the metabolism of the hexoses and disaccharides are very similar and differ only in the initial steps of metabolism (Fig. 1). For example, glucose and maltose metabolism differ only with respect to two reactions. The first reaction is the sugar transport through the plasma membrane; maltose uptake is catalyzed by an energy-dependent maltose-proton symport mechanism (Fig. 1, step 30) (2), whereas glucose uptake is catalyzed exclusively by a facilitated diffusion mechanism (step 33) (3). The second reaction is the intracellular breakdown of maltose into glucose, which involves a specific -glucosidase ("maltase," step 29) (4). Similarly, the metabolism of the C2-compounds ethanol and acetate only differ by the initial substrate-uptake step (steps 31 and 32) (57) and by two sequential NAD(P)+-dependent oxidation reactions that convert ethanol into acetate (steps 22 and 23).
Drastic changes are observed in central metabolism when the metabolism of sugars is compared with that of the C2-compounds. During growth on sugars, all metabolic building blocks can be derived from glycolysis, tricarboxylic acid cycle, and pentose phosphate pathway. During growth on C2-compounds, gluconeogenesis and glyoxylate cycle are essential for the provision of some of these precursors. Furthermore, the higher ATP requirement for biosynthesis during growth on the C2-compounds (and in particular acetate (Ref. 8)) implies that, at a fixed specific growth rate, dissimilatory fluxes have to be higher with the C2-compounds than with a sugar as the sole carbon source. So far, most studies on regulation of central carbon metabolism in S. cerevisiae have been performed in batch mode in shake-flask or reactors. This cultivation method, however, has several drawbacks for quantitative analysis. First of all, the concentrations of substrates and products change throughout cultivation, which makes it difficult to accurately measure fluxes through specific pathways or to assess the influence of carbon sources on cellular regulation. Furthermore, batch cultivation by definition requires the use of excess concentrations of the carbon source. When different carbon sources are compared, this will lead to different specific growth rates, which, in itself, may already lead to changes in the make-up and activity of the metabolic network. Finally, the relatively high substrate concentrations in batch cultures lead to catabolite repression and inactivation phenomena (9, 10). Cultivation of microorganisms in chemostats offers numerous advantages for studying the structure and regulation of metabolic networks (11). In chemostat cultures, individual culture parameters can be changed, while keeping other relevant physical and chemical culture parameters (composition of synthetic medium, pH, temperature, aeration, etc.) constant. An especially important parameter in this respect is the specific growth rate, which, in a chemostat, is equal to the dilution rate, which can be accurately controlled. This allows the experimenter to investigate the effects of environmental changes or genetic interventions at a fixed specific growth rate, even if these changes result in different specific growth rates in batch cultures. In a chemostat, growth can be limited by a single, selected nutrient. The very low residual concentrations of this growth-limiting nutrient in chemostat cultures alleviate effects of catabolite repression and inactivation. Furthermore, these low residual substrate concentrations prevent substrate toxicity, which, for example, occurs when S. cerevisiae is grown on ethanol or acetate as the carbon source in batch cultures (12, 13). The central goal of the present study is to assess to what extent carbon source-dependent regulation of fluxes through central carbon metabolism in S. cerevisiae is regulated at the level of transcription. To this end, we compare the transcriptome of carbon-limited, aerobic chemostat cultures grown on four different carbon sources: glucose, maltose, ethanol, and acetate. Data from the transcriptome analysis are compared with flux distribution profiles calculated with a stoichiometric metabolic network model. Questions that will be addressed are as follows: (i) does glucose-limited aerobic cultivation lead to a complete alleviation of glucose-catabolite repression; (ii) how (in)complete is our understanding of the genes involved in the transcriptional response of S. cerevisiae to four of the most common carbon sources for this yeast; and (iii) to what extent do transcriptome analyses with microarrays provide a reliable indication of flux distribution in metabolic networks? The complete data set used in this study is available for download at www.bt.tudelft.nl/carbon-source.
Strain and Growth ConditionsWild-type S. cerevisiae strain CEN.PK1137D (MATa) (14) was grown at 30 °C in 2-liter chemostats (Applikon), with a working volume of 1.0 liter as described in Ref. 15. Cultures were fed with a defined mineral medium that limited growth by glucose, ethanol, acetate, or maltose with all other growth requirements in excess. The dilution rate was set at 0.10 h-1. The pH was measured on-line and kept constant at 5.0 by the automatic addition of 2 M KOH with the use of an Applikon ADI 1030 biocontroller. Stirrer speed was 800 rpm, and the airflow was 0.5 liters·min-1. Dissolved oxygen tension was measured online with an Ingold model 34-100-3002 probe, and was between 60 and 75% of air saturation. The off-gas was cooled by a condenser connected to a cryostat set at 2 °C and analyzed as previously described (16). Steady-state samples were taken after 1014 volume changes to avoid strain adaptation caused by long term cultivation (17). Dry weight, metabolite, dissolved oxygen, and off-gas profiles had to be constant over at least five volume changes prior to sampling for RNA extraction. MediaThe defined mineral medium composition was based on that described by Verduyn et al. (18). The carbon source was 256 ± 19 mmol of carbon/liter. Analytical MethodsCulture supernatants were obtained after centrifugation of samples from the chemostats. For the purpose of glucose, ethanol, acetate, and maltose determination and carbon recovery, culture supernatants and media were analyzed by HPLC,1 fitted with an Aminex HPX-87H ion exchange column using 5 mM H2SO4 as the mobile phase. Culture dry weights were determined via filtration as described by Postma et al. (19).
Metabolic Flux DistributionIntracellular metabolic fluxes were calculated through metabolic flux balancing using a compartmented stoichiometric model derived from the model developed by Lange (20). Because the intracellular localization of certain enzymes as well as the trafficking of certain metabolites in S. cerevisiae are still a matter of debate, assumptions had to be made on these aspects. The main assumptions concern acetyl-coenzyme A metabolism. It has been shown that transport of acetyl-CoA through the mitochondrial membrane cannot be performed in S. cerevisiae by the carnitine shuttle in the absence of exogenous carnitine (21, 22). It is, however, possible to conceive growth on sugar and gluconeogenic carbon sources without acetyl-CoA transport. During growth on sugar, mitochondrial acetyl-CoA is synthesized directly in this compartment via the pyruvate dehydrogenase complex, whereas the small amounts of necessary cytosolic acetyl-CoA are synthesized in the cytosol by the acetyl-CoA synthase (encoded by ACS1 and ACS2 (Refs. 15 and 23)). The gluconeogenic carbon sources acetate and ethanol can be converted by acetyl-CoA synthase into acetyl-CoA in the cytosol, acetyl-CoA being further converted in the cytosol to citrate by citrate synthase (CIT2 (Ref. 23)). Citrate can then be transported trough the mitochondrial membrane by the well described citrate transporter (CTP1 (Ref. 24)). The set-up of the stoichiometric models for growth of S. cerevisiae on glucose, maltose, ethanol, and acetate, as well as the flux balancing, was performed using dedicated software (SPAD it, Nijmegen, The Netherlands). The theory and practice of metabolic flux balancing has been described well in literature and will not be repeated here (2529). For each carbon source the specific rates of growth, substrate consumption, carbon dioxide production, and oxygen consumption during steady-state chemostat cultivation were calculated from the measured concentrations and flow rates from three independent experiments. The calculated specific conversion rates and their variances were used as input for the metabolic flux balancing procedure. In all cases the ATP balance was omitted as constraint in the flux balancing. This is a prerequisite for proper balancing in case the ATP stoichiometry of some reactions is insufficiently known (e.g. maintenance energy requirements, P/O ratio, etc.). However, without the ATP balance, the number of measurements was sufficient to result in an overdetermined system, thus making data reconciliation possible. In all cases the degree of redundancy was equal to 2. The complete list of reactions and components used to build the model can be found at www.bt.tudelft.nl/carbon-source. Microarray AnalysisSampling of cells from chemostats, probe preparation, and hybridization to Affymetrix GeneChip® microarrays were performed as described previously (30). The results for each growth condition were derived from three independently cultured replicates. Data Acquisition and AnalysisAcquisition and quantification of array images and data filtering were performed using the Affymetrix software packages Microarray Suite version 5.0, MicroDB version 3.0, and Data Mining Tool version 3.0. Before comparison, all arrays were globally scaled to a target value of 150 using the average signal from all gene features using Microarray Suite version 5.0. From the 9,335 transcript features on the YG-S98 arrays, a filter was applied to extract 6,383 yeast open reading frames, of which there were 6,084 different genes. This discrepancy was the result of several genes being represented more than once when suboptimal probe sets were used in the array design. To represent the variation in triplicate measurements, the coefficient of variation (standard deviation divided by the mean) was calculated as previously described by Boer et al. (31). For further statistical analyses, Microsoft Excel running the Significance Analysis of Microarrays (SAM, version 1.12) add-in was used (32) for all possible pairwise comparisons of the four data sets. Genes were considered as being changed in expression if they were called significantly changed using SAM (expected median false discovery rate of 1%) by at least 2-fold from each other's conditions. Hierarchical clustering of the obtained set of significantly changed expression levels was subsequently performed by GeneSpring (Silicon Genetics). Promoter analysis was performed using web-based softwares Regulatory Sequence Analysis Tools2 (RSA Tools (Ref. 33)) and AlignAce3 (34). The promoters (from -800 to -50) of each set of co-regulated genes were analyzed for over-represented motifs. When motifs shared largely overlapping sequences, they were aligned to form longer conserved elements. All the individual promoter sequences contributing to these elements were then aligned, and redundant elements were determined by counting the base representation at each position. The relative abundance of these redundant elements was then determined from a new enquiry of the co-regulated gene promoters and the entire set of yeast promoters in the genome using RSA Tools. The cluster coverage (Table V) is then expressed as the number of genes in the cluster containing the motif at least once, divided by the total number of genes in the cluster. Similarly, the genome coverage is expressed as the genome-wide number of genes containing the motif at least once, divided by the total number of genes used by RSA Tools (6,451 open reading frames).
Biomass Yields and Respiration Rates in Carbon-limited Chemostat CulturesCarbon-limited, aerobic chemostat cultures were grown on glucose, maltose, ethanol, and acetate as single growth-limiting nutrients. At a dilution rate of 0.1 h-1, the concentration of all four carbon sources in the reservoir medium was 250 mmol of carbon/liter (Table I), whereas their residual concentrations in steady-state cultures were below their respective detection limits (i.e. less than 0.5 mM). For the glucose- and maltose-grown cultures, the respiratory quotient (ratio of the specific rate of carbon dioxide production and oxygen consumption) was close to 1.0, indicating a fully respiratory metabolism of these sugars. In all cultures over 95% of the substrate carbon was recovered as either biomass or carbon dioxide (Table I), and HPLC analysis of culture supernatants did not reveal the production of any low molecular weight metabolites.
Biomass yields, as well as the specific rates of oxygen consumption and carbon dioxide production, were very similar for cultures grown on maltose and glucose. This was expected, as maltose metabolism is initiated by the uptake and hydrolysis of the disaccharide to two glucose molecules (Fig. 1). In contrast to glucose uptake, however, the uptake of maltose occurs via an energy-dependent proton-symport mechanism (2), which is likely to be responsible for the slightly lower (8%) biomass yield on maltose as compared with glucose (Table I). Consistent with earlier studies, biomass yields on ethanol and, in particular, acetate were lower than that on glucose (8, 35). These lower biomass yields and correspondingly higher respiration rates can be explained from the lower ATP yield from respiratory dissimilation of these substrates (which is largely the result of the investment of 2 ATP equivalents in the acetyl-coenzyme A synthetase reaction) and from the necessity to synthesize biosynthetic precursors via the glyoxylate cycle and gluconeogenesis. The difference in biomass yield between ethanol and acetate can be attributed to two factors: (i) energy-dependent uptake of acetate via a proton symport mechanism (6, 7) and (ii) the higher degree of reduction of ethanol, for which the ethanol and acetaldehyde dehydrogenases (Fig. 1, reactions 22 and 23, respectively) can yield NAD(P)H that can either be used for biosynthesis or yield ATP via oxidative phosphorylation. Flux Distribution in Central Carbon MetabolismIn vivo fluxes through central pathways in carbon metabolism were estimated by metabolic flux analysis, using a stoichiometric model of the S. cerevisiae metabolic network. As S. cerevisiae is a eukaryote, this model took into account metabolic compartmentation by discriminating between reactions that occur in the yeast cytosol and in the mitochondrial matrix (36, 37). To calculate intracellular fluxes, the model was fed with quantitative data on the biomass composition of S. cerevisiae (38) and with the substrate consumption and product formation rates observed in the carbon-limited chemostat cultures. The estimated fluxes of central carbon metabolism relevant for this study are summarized in Tables II and III.
With the exception of the sugar transport and maltose hydrolysis steps, there were only very few predicted changes in central carbon metabolism when either glucose or maltose was used as the carbon source (Tables II and III). The slight increase of catabolic fluxes in maltose-grown cultures was caused by the ATP requirement for maltose uptake (2). Larger changes were predicted between the C2 substrates ethanol and acetate. When ethanol is the carbon source, NADP-dependent acetaldehyde dehydrogenase can make an important contribution to fulfill the cellular demand for NADPH (Table II). In contrast, in acetate-grown cultures, this important reduced cofactor must be regenerated by NADP-dependent isocitrate dehydrogenase and possibly by the pentose-phosphate pathway (39). Furthermore, the lower degree of reduction of acetate and the associated lower yield of reducing equivalents during its dissimilation result in a lower ATP yield and necessitate higher fluxes through dissimilatory pathways. This was reflected by a substantially higher predicted in vivo activity of the tricarboxylic acid cycle in acetate-grown cells (Table III). Major metabolic rearrangements were predicted when ethanol or acetate replaced glucose or maltose as the sole carbon source for growth of S. cerevisiae. Predicted fluxes through acetyl-coenzyme A synthetase (Fig. 1, reaction 24) were low in glucose-grown cultures, where this enzyme is only required for the provision of relatively small amounts of cytosolic acetylcoenzyme A that are needed for lipid and lysine biosynthesis (40, 41). In contrast, predicted fluxes through this enzyme were high in ethanol- and acetate-grown cultures (Table II). Another important difference between growth on sugars and growth on C2-compounds is the involvement, for the latter substrates, of the glyoxylate cycle (Fig. 1, reactions 14, 15, 21, 25, and 26) and gluconeogenesis (steps 27 and 28). Furthermore, growth on C2-compounds led to a dramatic increase of predicted fluxes through the tricarboxylic acid cycle relative to sugar-grown cultures (Table III). Conversely, the flux through glycolysis was reversed (gluconeogenesis) and much lower in ethanol- and acetate-grown cells. Finally, the predicted flux through the oxidative pentose-phosphate pathway was reduced, as NADP-dependent acetaldehyde and/or isocitrate dehydrogenases provide alternative sources of NADPH during growth on C2-compounds. Global Transcriptome Changes in Chemostat Cultures Limited for Different Carbon SourcesIndependent triplicate chemostat cultures were run for each carbon limitation, followed by genome-wide transcriptional analysis with oligonucleotide DNA microarrays. Consistent with the excellent reproducibility reported in earlier studies in which DNA microarray analysis was applied to chemostat cultures (30, 32), the average coefficient of variation for the independent triplicate analyses did not exceed 0.18 (Table IV). Furthermore, the levels of ACT1 and PDA1 transcripts, which are commonly applied as loading standards for conventional Northern analysis, were not significantly different for cultures grown on different growth-limiting carbon sources (Table IV). The lowest measurable signal was 12 (arbitrary units) in all conditions (Table IV). 833 transcripts (13.7% of the genome) remained below this detection limit for all four growth-limiting carbon sources tested.
As few as 180 genes (only 2.8% of the genome) were carbon source-responsive, as defined by a significant change in the transcription profile (-fold change greater than 2 and false discovery rate of 1%; see "Experimental Procedures"). This number seems low, especially when considering the major rearrangements in metabolism that must occur when, for example, ethanol replaces glucose as the sole carbon source for growth (Table I). As shown on Fig. 2, many of the carbon source-responsive genes (72 genes, 40%) have not yet been assigned a biological function, as defined by MIPS4 (42). A significant number (50 genes, 28%) were found to be related to carbon metabolism-encoding enzymes (35 genes), transporters (6 genes), or proteins involved in regulation (9 genes). Finally, 11 genes (6%) were involved in nitrogen metabolism.
Specific Transcriptional Responses to Growth-limiting Carbon SourcesTo analyze the specific transcriptional response of S. cerevisiae to the four growth-limiting carbon sources, the transcriptome data were subjected to hierarchical cluster analysis. This resulted in six distinct clusters (Fig. 3). Most of the C2-responsive genes found in the present study (117 of 180) were present in clusters 2 and 5 (Fig. 3), which yielded different transcript levels for growth limited by sugar (glucose or maltose) and for growth limited by C2-compounds (ethanol or acetate), respectively. The other clusters consisted of genes that specifically responded to a single growth-limiting carbon source: low expression under acetate limitation (cluster 1, 16 genes), high expression under glucose limitation (cluster 3, 13 genes), high expression under maltose limitation (cluster 4, 16 genes), and, finally, low expression under maltose limitation (cluster 6, 18 genes). No genes were found with a specific response to growth under ethanol limitation, or with a specifically increased transcript level in acetate- or glucose-limited cultures. The six clusters that were identified will be briefly discussed below.
Different Transcript Levels in Sugar and C2-compound-limited Cultures (Clusters 2 and 5)Of the 117 genes that yielded different transcript levels in sugar-limited cultures and cultures limited by either ethanol or acetate, 79 were up-regulated (cluster 5) and 38 down-regulated (cluster 2) in cultures limited by the C2-compounds. Among the up-regulated transcripts, 21 encoded enzymes or regulatory proteins related to carbon metabolism and included the four structural genes for gluconeogenesis and glyoxylate cycle enzymes: PCK1 and FBP1 (encoding the gluconeogenic enzymes phosphoenolpyruvate carboxykinase and fructose-1,6-bisphosphatase) and ICL1 and MLS1 (encoding the glyoxylate cycle enzymes isocitrate lyase and malate synthase, respectively). An additional seven of these genes encoded enzymes or subunits active in the tricarboxylic acid cycle. This included five cytosolic enzymes (IDH1, IDH2, FUM1, SDH1, SDH3) and two mitochondrial enzymes (CIT2, MDH2) that were already known to be up-regulated during growth on non-fermentable carbon sources in batch cultures (43). Five additional genes encoded enzymes involved in acetyl-coenzyme A metabolism and its trafficking across intracellular membranes (ACS1, ACH1, CAT1, YAT1, YAT2), consistent with the key role of this intermediate in the metabolism of C2-compounds. Three transcripts that showed an increased level in cultures limited by C2-compounds are involved in transcriptional regulation of carbon metabolism: SIP4, a transcriptional activator of gluconeogenic genes SIP2 and REG2. Finally, RKI1 (encoding ribose-phosphate isomerase, involved in the pentose-phosphate pathway), INO1 (encoding inositol-1-phosphate synthase), and ICL2 (encoding a 2-methylisocitrate lyase involved in propionate metabolism) transcripts were also up-regulated in the presence of ethanol or acetate.
Eight of the genes that showed increased transcript levels in cultures limited by C2-compounds encoded proteins involved in transport across the cytosolic and the mitochondrial membranes. In addition to the CAT1, YAT1, and YAT2 genes mentioned above, two further up-regulated genes were involved in acetyl-CoA trafficking via the L-carnitine shuttle. AGP2 encodes a plasma membrane carnitine transporter, whereas the gene product of CRC1 transports acetylcarnitine across the mitochondrial inner membrane. Consistent with published results on batch cultures, SFC1, which encodes the mitochondrial succinate-fumarate exchanger, was expressed at increased levels in ethanol- and acetate-limited chemostat cultures (44). Surprisingly, STL1 also exhibited elevated transcript levels in ethanol- and acetate-limited cultures ( Of the 38 genes down-regulated in the presence of ethanol and acetate (cluster 2), more than half (20 genes) have not yet been assigned a clear biological role. Remarkably, of the remaining 18, 10 were linked to carbon metabolism: 2 glycolytic genes (HXK1, TDH1), 4 genes from the pentose-phosphate pathway (TKL2, GND1, GND2, and SOL3), VID24 encoding a protein involved in fructose-6-bisphosphate vacuolar transport and degradation, 2 members of the hexose transport family (HXT2 and HXT7), and finally MTH1, involved in glucose signaling and repression (46). Surprisingly, one gene (CYC7) encoding iso-2-cytochrome c and involved in respiration was repressed, whereas respiration rate was increased in cells grown with C2-compounds compared with sugars (Table I). Among the remaining seven genes were BAP2, a branched-chain amino acid permease; SPS100 and SWM1, both involved in sporulation; and FDH1 and FDH2, encoding formate dehydrogenases for which a precise role in yeast metabolism has not been clearly defined yet (47). Two other down-regulated genes respond to stress conditions, DOG2 responds to oxidative and osmotic stress and PDR12 is involved in weak organic acid resistance.
Up- and Down-regulation in Response to MaltoseA set of 34 genes specifically responded to growth with maltose. Of the 16 up-regulated transcripts, 6 had poorly described biological functions, 4 encoded ribosomal proteins (RPS10A, RPS26B, RPL15B, and RPL31B), and 6 coded for proteins involved in maltose utilization. Genes necessary for maltose degradation (i.e. maltose permeases MALx1, maltases MALx2 and maltose transcription activators MALx3) are clustered on MAL loci carried by subtelomeric regions (48). S. cerevisiae strains contain up to five highly homologous MAL loci, and one locus is enough to sustain growth on maltose (4). Two MAL loci have been sequenced in S. cerevisiae S288C, the strain used for the genome sequencing program, although this yeast is not able to grow on maltose as the sole carbon source (49, 50). Our laboratory strain, like all members from the CEN.PK family, can grow on maltose and contains four MAL loci.5 The Affymetrix microarrays contain four probe sets that are specific to known MAL genes (MAL11, MAL13, MAL33, and MAL23) and two that can hybridize to several homologues (MALx2 and MALx1). As expected, genes encoding for maltose permeases and maltases were strongly induced in the presence of maltose (610-fold), as did FSP2, YJL216C, and YGR287C, genes sharing high homology with maltases. Despite their recent identification as
Of the 18 genes down-regulated in the presence of maltose, most have not been assigned a function yet. The remaining genes cover a broad range of functional categories (budding, protein processing, DNA repair, cell wall maintenance, etc.) and could not be directly linked to maltose utilization. Down-regulation in Response to AcetateAs few as 16 transcriptional changes were observed between cells grown on ethanol or acetate, and all of them were down-regulations in the presence of acetate. Among these changes five were related to carbon metabolism: ALD5, PFK26, GPH1, YNL134C, and SOL4. ALD5, encoding a mitochondrial minor isoform of acetaldehyde dehydrogenase, has previously been shown to be induced in the presence of ethanol (53). In apparent contrast, genome-wide transcription analysis indicated a mild repression of its expression after the diauxic shift (43). The specific physiological role of this acetaldehyde dehydrogenase isoenzyme has not been elucidated (54), and a proposed role in maintenance of the electron transport chain (53) does not shed light on its down-regulation in the presence of acetate. PFK26 codes for a 6-phosphofructokinase catalyzing fructose-2,6-bisphosphate production. This metabolite has been shown to activate phosphofructokinase (Fig. 1, step 3) and to inactivate fructose bisphosphatase (Fig. 1, step 28), although its role in the switch between glycolysis and gluconeogenesis remains unclear (55, 56). This significant down-regulation of PFK26 in the presence of acetate (2.3-fold) and minor down-regulation (1.6-fold) in the presence of ethanol are, however, indicating a transcriptional control of PFK26 by carbon source in carbon-limited cultures. This may be relevant for controlling the in vivo fluxes through the antagonistic enzymes phosphofructokinase and fructose-1,6-biphosphatase. YNL134C is a member of the zinc-binding dehydrogenase family, which catalyzes the reversible oxidation of ethanol to acetaldehyde. This glucose-repressed gene (31, 57) is the only alcohol dehydrogenase-encoding ORF displaying a significant change in expression in response to carbon source identified in this study; however, its physiological role has not been clearly identified so far. Finally, GPH1 encodes a glycogen phosphorylase involved in glycogen degradation (58), and SOL4 expression product is a 6-phosphogluconolactonase involved in the pentose-phosphate pathway. No transcripts related to acidic stress response were up-regulated, confirming that acetate limitation efficiently resolved acid stress problems that are typically observed in batch cultures. Furthermore, three genes responding to acid and/or oxidative stress were down-regulated (MSC1, AHP1, GAD1). Up-regulation in Response to GlucoseA puzzling cluster is the small subset of 13 transcripts that were specifically up-regulated when S. cerevisiae was grown with glucose as the sole carbon source. Six genes have unknown or poorly described functions. Four genes (SUC2, SUC4, HXT4, and GIP2) are involved in carbon metabolism. Invertase, encoded by SUC homologues, is a well described target of catabolite repression by glucose, both at the level of transcription and mRNA stability (59). However, it has been shown that low concentrations of glucose (0.1%) are necessary for a maximum expression of SUC genes (60). These findings are in good agreement with the significantly higher expression level of SUC2 and SUC4 (34-fold) measured when cells were grown in glucose compared with maltose, ethanol, or acetate. HXT4, encoding a moderate to low affinity hexose transporter, was also responding to low extracellular glucose concentration in our chemostat cultivations. The up-regulation of GIP2 (protein phosphatase PP1-interacting protein), as well as GNP1 (glutamine permease), SKN1 (glucan synthase subunit), and MCH5 (member of the major facilitator superfamily), in the presence of glucose does not have a clear physiological relevance. In Silico Promoter Analysis of Carbon Source-regulated GenesCo-regulation of global transcription is generally controlled by the specific binding of common activating or repressing proteins (transcription factors) to short sequences contained in promoter regions of the regulated genes. Searching the promoter regions of co-regulated genes for over-represented short sequences can identify these binding sites. We analyzed the upstream sequences of the genes from the six clusters defined above using web-based tools (for further information see "Experimental Procedures"). No significantly over-represented sequences were recovered from three clusters, namely low expression on acetate, high expression with glucose, and low expression with maltose (Table V). 11 putative promoter elements were over-represented in the three remaining clusters, of which only 4 could be associated to known transcription factors. Among the promoter regions of the genes down-regulated in the presence of ethanol and acetate, three sequences were found to be over-represented compared with their genome coverage. One of them resembles the binding site targeted by the transcription factors Msn2p and Msn4p known as stress response element (STRE, Table V). However, the down-regulation of Msn2p/Msn4p-regulated genes seems unlikely as, in our experiments, cells were not exposed to any stress known to trigger the activity of these factors (reviewed by Estruch (Ref. 61)). Furthermore, a comparison of the genome-wide and cluster 2-wide coverage of the core sequence of the STRE (AGGGG) did not show any significant over-representation of this element among the genes down-regulated with C2-carbon sources. This renders it unlikely that the recovered sequence is a "true" STRE element. The induction of genes necessary for maltose utilization (MALx1 and MALx2) requires the activation by the transcription factor Malx3p (48). It was therefore expected that one of the binding sites targeted by Malx3p would be over-represented in the promoter region of genes up-regulated when S. cerevisiae is grown in the presence of maltose. Indeed, a 10-nucleotide sequence with significant homology to the Mal63p binding site (62) was retrieved from cluster 4 (Table V). Surprisingly, three other sequences were significantly over-represented (453-fold compared with genome coverage) in this relatively small cluster, but could not be associated to any known DNA-binding protein (Table V). Four rather well conserved sequences were recovered from the promoter regions of genes up-regulated in the presence of C2-carbon sources (cluster 5). Two of them could be identified as targets of known transcription factors: Mig1p and Cat8p/Sip4p (Table V). In glucose-grown batch cultures, many genes involved in the utilization of alternative carbon sources are repressed via a complex of signals known as carbon catabolite repression (10, 63). Glucose repression of transcription is mainly exerted via the Mig1p-binding protein. This is consistent with the finding that 24 genes (30%) up-regulated when S. cerevisiae is grown with acetate or ethanol contain a Mig1p binding site in their promoter region. From this set, most of the 18 genes with known function have previously been experimentally shown to contain a Mig1p binding site (REG2 (Ref. 64), FBP1 (Ref. 65), and ICL1 (Ref. 66)) or at least to be repressed by glucose (ACH1, MDH2, ICL1, NDE2, YAT1, CAT2, YAT2, CRC1, SFC1, and SDH1 (Refs. 31 and 6769)). MIG1 transcription itself did not respond to carbon source, which is consistent with previous reports of Mig1p activity control by post-translational modification and nuclear export (70, 71). An important cis-acting element for the transcriptional response to gluconeogenic carbon sources is Cat8p (for review, see Ref. 72). Cat8p binds to an 11-bp upstream activation sequence element named CSRE (carbon source-responsive element) and most recently defined as CCrTysrnCCG (r = AorG,y = CorT, s = C or G, n = A, T, G, or C) (73). Most of the genes that have been shown to contain functional Cat8p binding sites were indeed up-regulated in response to growth with C2-carbon sources (SIP4, FBP1, ICL1, MLS1, PCK1, MDH2, SFC1, and CAT2). Furthermore, screening the upstream untranslated regions of all genes up-regulated in the presence of C2-compounds allowed the identification of 20 additional genes containing one or more CSRE-related sequences (Table VI). In contrast to most of its targets, CAT8 itself was not affected by the carbon source (Figs. 3 and 4). This behavior is consistent with earlier reports of Cat8p activation at the post-translational level by phosphorylation (74).
No other binding sites corresponding to transcription factors known to be involved in carbon source adaptation (Adr1p, Hap complex, Rtg regulators) could be recovered from this set of C2-compound co-regulated genes. ACS1 and ADH2, the main targets of Adr1p (75, 76), were not significantly up-regulated in our culture conditions and were therefore not included in the promoter analysis. Similarly, few of the potential targets of the Hap complex and the Rtg regulators (SDH1, CIT2, IDH1, and IDH2) were up-regulated in this study and could therefore not result in a significant over-representation of their binding motifs in the cluster.
Transcription Analysis of Chemostat Cultures Versus Batch CulturesIn carbon-limited, aerobic chemostat cultures, remarkably few transcripts exhibited significant differences when the growth-limiting carbon source was changed. The change from a sugar (glucose or maltose) to a C2-compound (ethanol or acetate) resulted in only 117 genes with a significantly changed transcript level (Fig. 3). This robustness of the yeast transcriptome in response to changes of the carbon source is in contrast with previously reported data from batch cultures. In a transcriptome analysis of the diauxic shift in S. cerevisiae, which essentially represents the transition from growth on glucose to growth on ethanol, over 400 transcripts were found to change by more than 2-fold (43). Similarly, in an independent study on glucose- and ethanol-grown batch cultures, over 600 transcripts were found to change (77). When comparing the carbon source response in chemostat with these previous reports in batch cultivation, approximately half of the 117 genes that responded to C2-compounds in chemostat cultivations were also transcriptionally regulated during these two transcriptome studies.
Several factors may have contributed to these different carbon source-dependent transcriptional responses in batch and chemostat cultures. A large number of genes (225 based on work in our laboratory (Ref. 31)) are transcriptionally regulated by glucose, for instance via glucose catabolite repression. This phenomenon does not occur in glucose limited chemostat cultures, where the low residual glucose concentration ( Carbon Source-dependent Transcripts and Regulation MechanismsCarbon metabolism in S. cerevisiae is one of the most intensively studied metabolic systems. Indeed, in many cases, the carbon source-dependent transcription of genes could be teleologically explained from the known catalytic or regulatory functions of their gene products. Some obvious examples include the high transcript levels of glyoxylate cycle and gluconeogenic genes during growth on C2-compounds and the induction of MAL genes in maltose-limited chemostat cultures. However, in many other cases, the function of carbon source-responsive genes is either entirely unknown (the frequency of genes with unknown function was 40%, which is significantly higher than that of the entire yeast genome (27% according to YPDTM)) or difficult to interpret in terms of biochemical function of the gene product. Although there is no a priori proof that carbon source-dependent transcriptional regulation correlates with physiological function, our data provide an interesting lead for future functional analysis research.
As expected, very few transcriptional differences were measured between cultures grown on glucose and maltose (34 changes). The MAL structural genes (maltose permeases and maltases), but not the MAL activators, were up-regulated in the presence of maltose. Analysis of the upstream region of these co-regulated genes resulted in the identification of the MAL regulator binding site (85) (Table V). Cultivation on maltose, however, generated a puzzling set of 18 genes down-regulated compared with glucose cultivations. These changes suggest that maltose utilization, despite its closeness to glucose utilization, results in more metabolic perturbations than is generally recognized. Comparing the transcriptomes of ethanol- or acetate-grown cells resulted in even fewer differences. As few as 16 genes yielded different transcript levels for these two gluconeogenic compounds, all being down-regulated with acetate. Acetate uptake is mediated by a protein for which the encoding gene has not been identified so far (6, 7). Unfortunately no potential candidate for an acetate transporter could be identified from our data set. Finally, comparing C2-compound-limited cultures with sugar-limited cultures resulted in the identification of 117 carbon source-responsive genes. Among these genes, a high proportion have unknown or poorly defined biological functions (48%), but as many as half of the remainder (34 genes, 29%) are involved in carbon metabolism. The promoter analysis of co-regulated genes resulted in the definition of a new set of genes containing one or more sequences closely related to the carbon source-responsive element (Table VI). Using the totally different approach of integrating transcriptome and proteome comparison of a wild-type and a cat8 Transcript Levels Versus Metabolic FluxesFrom a combination of metabolic flux analysis and transcriptome analysis, we were able to compare metabolic fluxes in central carbon metabolism and transcript levels of key structural genes that encode enzymes of the central metabolic pathways (Tables II and III). Three distinct types of correlation between transcript levels and fluxes could be identified. For the pathways that are specific for maltose metabolism (maltase) or for the metabolism of gluconeogenic carbon sources (gluconeogenesis and glyoxylate cycle), there was a strong qualitative correspondence between transcriptional regulation and estimated in vivo metabolic fluxes. The genes involved in these pathways are known to be strongly transcriptionally regulated, but several of the corresponding enzymes are also subject to post-translational regulation (inactivation by phosphorylation and/or ubiquitin-catalyzed degradation (Refs. 88 and 89)). The time scale of these post-translational processes is much shorter than that of transcriptional regulation. Although one would intuitively assume that post-translational processes are predominantly important under dynamic conditions, we cannot presently exclude the possibility that they also contribute to the regulation of in vivo enzyme activity during carbon-limited steady-state cultivation in chemostat cultures. The changes in metabolic fluxes in the tricarboxylic acid cycle and in the pentose-phosphate pathway were only partially mirrored by changes in transcript levels. Little is known about the regulation of these two pathways. From our data it appears clearly that, to meet the new flux requirements when C2 carbon sources are used, S. cerevisiae only enhanced or repressed the transcription of a few genes. One can assume that only the rate controlling steps need to have their protein concentration optimized to adjust the flux to the new requirement. In such a case, transcriptome analysis would help to identify the potential rate-limiting steps, i.e. succinate dehydrogenase and fumarase for the tricarboxylic acid cycle (CIT2 and MDH2 induction mainly reflecting the need of the corresponding proteins in a different compartment rather than a real rate-limiting step) and 6P-gluconolactonase, 6P-gluconate dehydrogenase, ribose-P isomerase, and transketolase 2 for the pentose-phosphate pathway. However, the magnitude of change in transcript levels does not correlate well with the magnitude of change in metabolic fluxes, indicating that transcriptional control alone cannot explain the modifications in the flux distribution in the tricarboxylic acid cycle and in the pentose-phosphate pathway. The model-predicted fluxes in the glycolytic pathway and in the enzymic reactions surrounding pyruvate were strongly dependent on the studied carbon sources. However, the different in vivo activities of the key enzymes of these pathways were not at all mirrored by their transcript levels. In glycolysis, only HXK1 (encoding hexokinase I) and TDH1 (encoding a minor isoform of glyceraldehyde dehydrogenase) displayed a reduced transcription level when the glycolytic flux was decreased at least 3-fold. This is a clear indication that, during carbon-limited cultivation, fluxes through these central metabolic pathways in S. cerevisiae are not primarily controlled at the transcriptional level. Further research is required to assess the contribution of translational efficiency, post-translational modification, and regulation by intracellular concentrations of substrates, products, and effectors to the regulation of in vivo activity of these pathways. Such discrepancies between fluxes and/or enzyme activities and transcript levels have already been reported by other integrative approaches with S. cerevisiae (90) and several bacteria (9193). Our study underlines that DNA microarrays, however useful for studying transcriptional regulation, comparative genotyping, and purely correlation-based diagnostics, have limited value as indicators for in vivo activity of proteins. This limitation should be considered when applying DNA microarrays as a tool for activities such as metabolic engineering or identification of potential drug targets.
* This work was supported by a Technologiestichting STW grant and by a Marie-Curie fellowship (to P. D.-L.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 The abbreviations used are: HPLC, high performance liquid chromatography; STRE, stress response element; CSRE, carbon source-responsive element.
2 URL is rsat.ulb.ac.be/rsat/.
3 URL is atlas.med.harvard.edu//.
4 Munich Information Center for Protein Sequences URL is mips.gsf.de/.
5 P. Daran-Lapujade, J.-M. Daran, T. Petit, and J. T. Pronk, unpublished results.
The research group of Prof. J. T. Pronk is part of the Kluyver Centre for Genomics of Industrial Fermentation, which is supported by the Netherlands Genomics Initiative. We thank H. C. Lange for providing the compartmented stoichiometric model. We also thank L. M. Serrano, M. de Herdt, and J. M. A. Geertman for technical assistance and Prof. J. P. van Dijken for critically reading the manuscript.
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||