A comparative transcriptomic, fluxomic and metabolomic analysis of the response of Saccharomyces cerevisiae to increases in NADPH oxidation
© Celton et al.; licensee BioMed Central Ltd. 2012
Received: 19 May 2012
Accepted: 28 June 2012
Published: 17 July 2012
Redox homeostasis is essential to sustain metabolism and growth. We recently reported that yeast cells meet a gradual increase in imposed NADPH demand by progressively increasing flux through the pentose phosphate (PP) and acetate pathways and by exchanging NADH for NADPH in the cytosol, via a transhydrogenase-like cycle. Here, we studied the mechanisms underlying this metabolic response, through a combination of gene expression profiling and analyses of extracellular and intracellular metabolites and 13 C-flux analysis.
NADPH oxidation was increased by reducing acetoin to 2,3-butanediol in a strain overexpressing an engineered NADPH-dependent butanediol dehydrogenase cultured in the presence of acetoin. An increase in NADPH demand to 22 times the anabolic requirement for NADPH was accompanied by the intracellular accumulation of PP pathway metabolites consistent with an increase in flux through this pathway. Increases in NADPH demand were accompanied by the successive induction of several genes of the PP pathway. NADPH-consuming pathways, such as amino-acid biosynthesis, were upregulated as an indirect effect of the decrease in NADPH availability. Metabolomic analysis showed that the most extreme modification of NADPH demand resulted in an energetic problem. Our results also highlight the influence of redox status on aroma production.
Combined 13 C-flux, intracellular metabolite levels and microarrays analyses revealed that NADPH homeostasis, in response to a progressive increase in NADPH demand, was achieved by the regulation, at several levels, of the PP pathway. This pathway is principally under metabolic control, but regulation of the transcription of PP pathway genes can exert a stronger effect, by redirecting larger amounts of carbon to this pathway to satisfy the demand for NADPH. No coordinated response of genes involved in NADPH metabolism was observed, suggesting that yeast has no system for sensing NADPH/NADP+ ratio. Instead, the induction of NADPH-consuming amino-acid pathways in conditions of NADPH limitation may indirectly trigger the transcription of a set of PP pathway genes.
Redox homeostasis is a fundamental requirement for the maintenance of metabolism. Intracellular redox potential is determined principally by the ratio of NADH/NAD+ and NADPH/NADP+ cofactors, which are involved in about 200 reactions in Saccharomyces cerevisiae. These metabolites are interconnected by many different pathways. Any change to intracellular redox balance therefore has far-reaching effects on the metabolic network. Redox engineering strategies thus often have wide-ranging effects on metabolism [2–4]. A better understanding of the way in which yeast efficiently adjusts its metabolic fluxes in response to genetic or environmental redox changes is therefore essential for the prediction and control of these effects.
Despite these significant advances in our understanding of NADPH metabolism, little is known about the mechanisms regulating NADPH homeostasis. It is generally thought that the pentose phosphate pathway is controlled principally at the enzymatic level, with NADPH and ATP competitively inhibiting both the glucose-6 phosphate dehydrogenase Zwf1p and the 6-phosphogluconate dehydrogenase Gnd1p . The coordinated regulation of genes involved in NADPH metabolism, including most of PP pathway genes, has been reported in conditions of oxidative stress. The activation of NADPH-dependent genes involves Stb5p, a zinc-binding factor , which also represses the expression of PGI1, encoding the phosphoglucose isomerase at the junction between glycolysis and the PP pathway. This transcription factor plays a key role in rerouting carbon flux to provide the additional NADPH required for the response to oxidative stress, as demonstrated by the greater susceptibility of the stb5 Δ mutant to several chemicals and oxidants [7, 8] and by the greater resistance to diamide of a strain overexpressing STB5. Stb5p was also recently shown to be essential for acetaldehyde tolerance  and, under anaerobic conditions, to play a major role in the maintenance of basal flux through the PP pathway .
High-throughput experiments and metabolic modeling are powerful approaches for studying the global response of cells to environmental changes and identifying the regulatory mechanisms involved in the rerouting of the metabolic network. We analyzed the transcriptome, extracellular and intracellular metabolites and metabolic fluxes, to identify the regulatory mechanisms involved in the response of yeast cells to various increases in NADPH oxidation. We found that redox homeostasis was maintained principally through metabolic control, with regulation of the transcription of certain genes of the PP pathway occurring at higher levels of NADPH demand. We also observed a general increase in the expression of genes involved in amino-acid synthesis, particularly those involved in the synthesis of sulfur-containing amino acids, as an indirect response to the lower availability of NADPH.
Transcriptomic response to changes in NADPH and NADH demands
Funspec classification into biological process categories of genes up- or downregulated in response to changes in the oxidation of NADPH and NADH
NADH-Bdh 200 mM
NADPH-Bdh 100 mM
NADPH-Bdh 200 mM
NADPH-Bdh 300 mM
alcoholic metabolic bioprocess GO:0006066
response to stress GO:0006950
PAU24 PAU14 PAU1
PAU7 PAU24 PAU3 PAU10 PAU2 PAU11 PAU12 PAU13 PAU14 PAU15 PAU1 PAU17 PAU18 PAU23 PAU4 PAU6 MET22 PAU20 FRT1
PAU7 PAU2 PAU14 PAU17 PAU4 MET22 ZEO1 PAU20
RPL31A RPS16B RPS17B RPL27B RPL23B RPL30 RPL24A RPL24B RPL40B RLP24 RPS31 RPS25B RPS22B RPS17A RPL15B RPL9B RPS3 RPS19A RPL25 RPS28A RPL20B RPL33A
methionine biosynthetic process GO:0009086
MET6 MET3 MET14 MET17 MET22 MET16
MET6 MET10 MET5 MET14 MET17 MET22 MET16
L-serine biosynthetic process GO:0006564
lysine biosynthetic process GO:0009085
hydrogen sulfide biosynthetic process GO:0070814
oxidation reduction GO:0055114
BDH1 ADH7 MET10 ARI1 ERG3 AHP1 AAD15
BDH1 ADH7 AAD3 YDR541C SER3 MET10 ARI1 GND1 SER33 LYS12 MET5 MDH1 LYS9 AAD15 IDH2 MET16
pentose-phosphate shunt GO:0006098
SOL3 GND1 TAL1 TKL1
cofactor metabolic process GO:0051186
glucose catabolic process to ethanol GO:0019655
vacuolar protein catabolic process GO:0007039
PEX29 HXK1 STF2 TFS1 PGM2 DDR48 GPH1
glucose catabolic process to ethanol GO:0019655
translational elongation GO:0006414
ZRT1 PRY1 TPO1 AGA1 WHI5 PDR12 SAM3
THI2 ZRT1 QDR2 PRY1 IZH2 WHI5 PDR12 SAM3
carbohydrate metabolic process GO:0005975
EMI2 HXK1 PGM2 GPH1
The repressed genes included ADH4, which encodes an NADH-dependent alcohol dehydrogenase. This gene displayed the highest level of repression, by a factor of 5 at an acetoin concentration of 200 mM and a factor of 10 at an acetoin concentration of 300 mM. ADH4 was also downregulated (by a factor of 3.7) in response to the modulation of NADH levels. These results suggest that the product of the ADH4 gene is involved in redox homeostasis. Adh4p is generally thought to make use of NAD+, but the role, distribution within the cell and cofactor specificity of this enzyme remain unclear .
At an acetoin concentration of 300 mM, several genes involved in carbohydrate metabolism and alcoholic fermentation, such as PGM2 EMI2 GPH1 HXK1 PDC1 and PDC5, were repressed. PGM2 is involved in the production of UDP-glucose, a precursor of trehalose and trehalose-6-phosphate (T6P) , whereas GPH1 and, probably, EMI2 are involved in glycogen catabolism. By contrast, PDC5 was upregulated in conditions in which NADH levels were modified.
Pathways involved in NADPH synthesis
Several of the downregulated genes are linked to the sulfate assimilation pathway: SAM3, which encodes a permease required for SAM utilization, and TPO1, which encodes an H+ transporter involved in exporting the spermidine produced from SAM. Another gene, QDR2 encodes a multidrug transporter activated in response to the limitation of amino-acid synthesis .
Genes encoding several other oxidoreductases, such as ADH7 ARI1, GND1, YDR541C (in the presence of 200 mM acetoin), MDH1, IDH2 and AAD15 (in the presence of 300 mM acetoin) were also upregulated in response to changes in NADPH demand. ADH7 showed the highest level of variation of these genes. It was overexpressed by a factor of 8 in the presence of 200 mM acetoin and by a factor of up to 14 in the presence of 300 mM acetoin. ADH7 encodes an NAD(P)H-dependent, medium-chain length alcohol dehydrogenase with broad substrate specificity thought to be involved in the accumulation of fusel alcohols, in lignin degradation and in the detoxification of furfural [17, 18]. ARI1 encoding an intermediate-subclass short-chain dehydrogenase/reductase that can use aromatic and aliphatic aldehydes as substrates and that contributes to the detoxification of furfural , was also strongly induced (by a factor of four) in the presence of 200 and 300 mM acetoin.
Despite the dependence of the corresponding enzymes on NADPH, ADH7 and ARI1 were also induced when NADH levels were modified, albeit to a lesser extent. One possible reason for this is that the changes in NADH levels also affected NADPH/NADP+ balance. Indeed, we previously suggested that a cycle exchanging NADPH for NADH operated in strain NADH-Bdh and that the activity of this cycle was increased by NADH oxidation, with the overall effect of decreasing NADPH availability . This would account for the induction of ADH7 and ARI1 by changes in NADH levels.
Rerouting metabolic flux in response to increases in NADPH demand
Flux distributions as a function of NADPH demand
Glucose uptake (mmol/gDW/h)
84.7 ± 0.9
78.8 ± 1.9
69.9 ± 0.6
65.2 ± 0.6
11.4 ± 0.9
18 ± 1.9
26 ± 0.9
33.0 ± 1.4
5.6 ± 1.5
9.0 ± 2.5
12.2 ± 2.7
18.6 ± 4.2
3.3 ± 0.1
3.4 ± 0.1
2.5 ± 0.1
2.1 ± 0.1
164.3 ± 0.1
164.8 ± 0.5
167.1 ± 0.1
167.9 ± 0.4
4.2 ± 0.8
2.7 ± 0.7
2.1 ± 0.6
2.5 ± 0.6
Import of oxaloacetate into the mitochondria
3.4 ± 0.8
2 ± 0.7
1.6 ± 0.6
2 ± 0.6
2.4 ± 0.1
3.8 ± 0.06
4.2 ± 0.1
5.6 ± 0.1
0.29 ± 0.01
0.17 ± 0.02
0.15 ± 0.02
0.09 ± 0.04
0.30 ± 0.03
0.20 ± 0.06
0.22 ± 0.1
0.13 ± 0.09
161.9 ± 0.1
161 ± 0.5
162.9 ± 0.5
162.3 ± 0.5
17.8 ± 0.1
17.9 ± 0.1
14.7 ± 0.1
13.7 ± 0.1
2.1 ± 0.1
3.6 ± 0. 1
4.0 ± 0. 1
5.5 ± 0.1
0.30 ± 0.01
0.26 ± 0. 01
0.26 ± 0. 01
0.28 ± 0. 01
Other changes in flux partitioning were observed at the pyruvate and acetaldehyde nodes. At these junctions, carbons were rerouted toward the formation of acetaldehyde and acetate, whereas flux through acetyl-CoA synthase decreased. These findings suggest an increase in the contribution of NADPH acetaldehyde dehydrogenase (Ald6p) to the provision of cytosolic NADPH (Table 2) and are consistent with the increase in acetate production observed (Figure2B).
Effects of an increase in NADPH oxidation on intracellular metabolites
We also determined the concentrations of nucleotides. The concentrations of the triphosphates decreased in response to increasing NADPH demand, by factors of 1.6 for ATP and 1.8 for UTP. By contrast, those of the monophosphate forms increased, by a factor of 2 for AMP and 1.6 for UMP. The glycolytic intermediates glucose 6-phosphate (G6P) and F6P accumulated, whereas the concentrations of G3P, 1,6-fructose-biphosphate (FBP) and fructose-2,6-bisphosphate (2,6 FBP) were similar at all levels of NADPH demand.
The UDP-glucose pyrophosphorylase (Ugp1p) is involved in the synthesis of UDP-glucose, a precursor of reserve carbohydrates, from glucose-1-phosphate (G1P). The accumulation of G1P observed in the presence of 300 mM acetoin may reflect a limitation of the flux through Ugp1p, due to a decrease in UTP availability.
An accumulation of 2-phosphoglycerate and 3-phosphoglycerate (measured as a single pool, denoted 2,3 PG) was also observed.
Impact of changes in NADPH and NADH availability on aroma compounds
Acetaldehyde, higher alcohols and esters produced by NADPH-Bdh and NADH-Bdh strains in response to various levels of redox disturbance
Acetoin initial (mM)
Acetoin consumed (mM)
Increases in the oxidation of NADH or NADPH (Table 3) decreased the overall production of ethyl esters, which are generated by the condensation of ethanol and a fatty acid derivative. However, the ethyl acetate production did increase (Table 3), consistent with the observed increase in acetate concentration (Table 2, Figure2B). By contrast, acetate esters, which are produced by the condensation of a higher alcohol and acetyl-CoA, did not change significantly in abundance, with the exception of isoamyl acetate in response to the greatest level of NADPH demand. The production of fusel alcohols was increased in response to altered NADH demand, but unaffected by the modification of NADPH demand. The level of isobutanol production was 2.4 times higher in cells in which NADH levels were modified than in control cells. In line with this, the most significant effect of NADH oxidation was a large increase in isobutyl acetate production. Similarly, the synthesis of propanol doubled in response to NADH oxidation but was unaffected by increases in NADPH oxidation (Table 3). Finally, NADPH oxidation tripled phenylethanol production. The precursor of this higher alcohol is erythrose-4-phosphate, an intermediate of the PP pathway. This result is therefore consistent with an increase in flux through the PP pathway (Table 2).
We investigated the molecular processes involved in the maintenance of intracellular NADPH homeostasis, by studying the fluxomic, metabolic and transcriptional responses to increases in NADPH oxidation in S. cerevisiae. The biological system used, which was based on the overproduction of an NADPH-dependent Bdh1p enzyme and the use of acetoin as an electron acceptor, made it possible to impose unprecedented levels of redox disturbance on yeast cells, corresponding to 22 times the anabolic demand for NADPH .
We previously showed using constraint based metabolic flux analysis (dMBA for dynamic Mass Balance Analysis, ) that (i) yeast responds to increases in NADPH demand by increasing the flux through the two main NADPH-producing pathways (the PP and acetate pathways) (ii) for NADPH demands of more than 22 times the anabolic demand, the PP pathway was saturated and additional NADPH was produced through the glycerol-DHAP cycle, involving exchanges of the NADH and NADPH redox cofactors. In this study, we characterized the mechanisms involved in flux rerouting for two moderate (100 and 200 mM acetoin) and one large (300 mM) increase in NADPH, corresponding to 8, 12 and 22 times the anabolic NADPH demand, respectively. In these conditions, little effect on growth was observed only for the two highest levels of NADPH demand, consistent with a decrease in biomass synthesis due to a lack of NADPH.
Regulatory mechanisms involved in the rerouting of carbon to NADPH-producing pathways
We combined 13 C-based metabolic flux, intracellular metabolite and global gene expression analyses, to investigate the way in which yeast adjusts its metabolism to respond to a high NADPH demand. We first confirmed experimentally that carbon was rerouted principally to the acetate and the PP pathways for the generation of NADPH. Flux through the PP pathway, which was estimated by 13 C-flux analysis, was increased by factors of 1.6, 2.4 and 2.9 in the presence of 100, 200 and 300 mM acetoin, respectively, consistent with the predictions of the DynamoYeast model .
We showed that the response of yeast to an increase in NADPH demand to eight times the anabolic demand was mediated mostly by metabolic control, as no transcriptional regulation of the genes involved in the acetate and PP pathways was observed. However, further increases in NADPH demand resulted in the induction of several PP pathway genes: SOL3 and GND1 in the presence of 200 mM and 300 mM acetoin and, in the presence of 300 mM acetoin, another two genes, TAL1 and TKL1. These genes were probably upregulated to increase the flux capacity of the PP pathway, which was limited in these conditions. Indeed, the accumulation of 6-phosphogluconate and of all intermediates of the nonoxidative part of the PP pathway (Figure6) suggests that the activities of 6-phosphogluconolactonase or 6-phosphogluconate dehydrogenase, catalyzed by Sol3p and Gnd1p, and of enzymes of the nonoxidative part of the pathway become limiting.
Both metabolic and genetic controls are therefore required to cover a demand for NADPH more than eight times greater than the anabolic demand for this compound.
GND1 upregulation has been observed in various situations in which the flux through the PP pathway increases [20–22]. This is not surprising, because this gene encodes a protein that catalyzes an irreversible reaction. By contrast, ZWF1, which encodes a protein that catalyzes the first irreversible NADPH-producing reaction of this pathway, was not transcriptionally regulated. These findings may be accounted for by differences in the cost of the proteins encoded by these genes. Indeed, a recent study in Escherichia coli suggested that the production of enzymes with a high cost (estimated from their abundance and molecular weight) is more tightly controlled than that of enzymes with a lower cost . If the same is true in yeast, it would account for the transcriptional regulation of GND1 rather than ZWF1, because the cost of Gnd1p is 10 times that of Zwf1p .
We initially hypothesized that the transcription factor Stb5p might be involved in the maintenance of redox homeostasis under anaerobic conditions. However, only eight of the 69 Stb5p target genes  (including 4 PP pathway genes) displayed differential regulation (Table 1, Additional file 1), consistent with Stb5p playing no role in the redirection of flux observed in response to increases in NADPH oxidation. This low proportion of Stb5p target genes responding to NADPH/NADP+ imbalance suggests that Stb5p does not sense NADPH changes directly during oxidative stress.
The repression of genes involved in the synthesis of reserve carbohydrates under conditions of the greatest NADPH demand might optimize NADPH synthesis by rerouting carbons towards the PP pathway. The same may also be true for the repression of genes involved in glycolysis and alcoholic fermentation. However, as pyruvate decarboxylase is also used to generate acetoin from acetaldehyde, the downregulation of PDC1 and PDC5 may result from disturbance of the acetaldehyde node. By contrast, PDC5 was upregulated in response to altered NADH demand. Whether this may have resulted from the accumulation of acetadehyde in this condition (Table 3) remains to be explored.
Indirect effects of changes in NADPH demand on amino-acid synthesis
The other main changes in transcription concerned the synthesis of amino acids, particularly the NADPH-consuming methionine and lysine pathways and serine synthesis, which is connected to the methionine pathway.
MET genes expression is normally activated by the Met4p transcription factor when the levels of these sulfur-containing compounds are low [25, 26]. Sulfate assimilation requires considerable amounts of NADPH (7 moles of NADPH are oxidized to NADP+ to generate 1 mole of methionine from 1 mole of SO42−). Our results therefore suggest that decreases in NADPH availability limit the flux through this pathway, resulting in low levels of sulfur-containing compound production and a global derepression of the sulfate assimilation pathway. These findings are consistent with the coordinated regulation of this pathway, which is controlled principally through feedback repression exerted primarily by cysteine . The relationship between the PP phosphate pathway and the sulfur pathway has long been known  and is illustrated by the sensitivity to osmotic stress of a zwf1 mutant. The connection between these two pathways does not seem to be based on transcriptional control. Indeed, no common regulator of the genes of these two pathways has been identified except Met4p, the main activator of the sulfur pathway, which can also bind GND2, the minor isoform of the 6-phosphogluconate dehydrogenase .
The upregulation of the LYS12 LYS9 and LYS2 genes suggests that the same mechanism may operate when lysine production decreases. On the other hand, SER3 and SER33 encoding the first enzyme of the serine synthesis pathway, were also upregulated despite the description of this pathway as NADH-dependent . However, as the biosynthetic reactions involving 3-phosphoglycerate, including those contributing to serine synthesis, constitute a highly branched network connecting the purine, thiamine, histidine and methionine biosynthesis pathways, the upregulation of these genes may be an indirect consequence of a decrease in the intracellular concentration of cysteine.
In addition, the repression of several genes linked to the sulfate assimilation pathway (SAM3, TPO1, QDR2) may reflect a need for the cells to conserve methionine, as the synthesis of this compound is probably compromised by the higher demand of NADPH.
Effects of redox alteration on the synthesis of aroma compounds
We observed distinct effects of redox alteration on the production of fermentative aroma compounds, depending on the cofactor (NADH or NADPH) and on the class of aroma compounds. Indeed, increases in the oxidation of NADH and NADPH decreased the production of ethyl esters but did not alter the production of acetate esters. Changes in esters profile could reflect a decrease in the availability of acetyl-CoA , a precursor of these molecules, as suggested by the lower level of flux through acetyl-CoA synthase (Table 2). This suggests that the alcohol-O-acetyltransferase responsible for catalyzing this reaction may be less sensitive to acetyl-CoA concentration than the esterases involved in ethyl ester synthesis.
On the other hand, although the synthesis of phenyl ethanol increased in response to altered NADPH demand, consistent with increased flux through the PP pathway, the synthesis of other fusel alcohols (isobutanol, propanol) was specifically increased in response to NADH alterations. Isobutanol is synthesized from valine, which is produced from α-acetolactate, a compound also involved in acetoin synthesis. We suggest that the synthesis of acetoin by the PDC route is favored in conditions of acetaldehyde accumulation (as described by Heux et al. ), and that α-acetolactate is rerouted towards the production of valine, isobutanol and isobutyl acetate. A key mechanism underlying this remodeling of flux may be the increase in PDC5 gene expression observed in response to NADH oxidation (Table 1). PDC5 encodes a protein involved in the decarboxylation step of the Erhlich pathway. The greater production of butanol and propanol may therefore be explained by increased PDC5 expression.
The energy problem and the glycerol-DHA cycle
Analysis of the metabolic response to an NADPH demand 22 times greater than anabolic demand revealed a metabolic shift that we attributed to saturation of the PP pathway . This conclusion is consistent with the intracellular metabolite profiles obtained in this study. Moreover, a previous constraint-based metabolic flux analysis suggested that a glycerol-DHA cycle exchanging NADH for NADPH was activated in these conditions . The intracellular metabolite profile revealed energetic inconsistencies (Figure6), consistent with the hypothetical operation of this cycle. Indeed, the ATP and UTP pools were markedly smaller in the presence of 300 mM acetoin than in the presence of the other two concentrations of this compound. AXP concentration is one of the major factors controlling the regulation of the glycolysis and carbohydrate pathways [14, 33]. The accumulation of glycolytic intermediates may result from a decrease in phosphofructokinase activity due to a decrease in ATP availability or from an increase in glycolytic flux. An increase in glycolytic flux was actually previously observed  in the presence of 300 mM acetoin (25.7 mmol/gDW/h instead of 19.5 mmol/gDW/h for lower −100 and 200 mM- acetoin concentrations).
Flux toward acetaldehyde, via pyruvate decarboxylase, increased, at the expense of the reactions supplying the TCA pathway. The NADPH imbalance results in the rerouting of the carbon atoms of acetaldehyde for the production of acetate via Ald6p. The marked accumulation of acetate when the cells were incubated in the presence of 300 mM acetoin (Figure2B) probably resulted from a lower level of flux through acetyl-CoA synthase, resulting from both ATP limitation and a smaller acetyl-CoA requirement for lipid biosynthesis. Consistent with the limited production of acetyl-CoA, the production of acetate esters and fatty acids from this precursor decreased in response to increases in NADPH demand.
This study highlights the remarkable robustness of yeast to changes in NADPH demand, based on efficient mechanisms for the reprogramming of metabolic flux to produce large amounts of NADPH. In line with previous studies [22, 34], we observed that fluxes are controlled by both metabolic and transcriptional regulation. However, the contribution of these mechanisms varies in function of the NADPH demand. Metabolic control was involved in the adaptation to increases in NADPH demand to up to eight times the anabolic demand. At higher levels of NADPH demand, transcriptional regulation of several genes also contributes to the rerouting of carbons towards the PP pathway.
Strains and culture conditions
The strains used in this study were derived from the S. cerevisiae haploid strain 59A . Strains NADH-Bdh and NADPH-Bdh were obtained by overexpressing the BDH1 gene encoding the NADH-dependent 2,3-butanediol dehydrogenase (Bdh1p) or a mutated version of this gene  encoding an NADPH-dependent Bdh1p . Overexpression was achieved by replacing the chromosomal BDH1 promoter with the TDH3 promoter and, for NADPH-Bdh, by replacing the native ORF with the mutated ORF.
S. cerevisiae strains were maintained and grown at 28 °C on YPD medium (1% Bacto yeast extract, 2% Bacto peptone, 2% glucose). An aliquot of each culture was mixed with glycerol to a final concentration of 20%, frozen and stored at −80 °C.
Fermentation medium and growth conditions
Batch fermentations were performed in one-liter bioreactors equipped with fermentation locks for the maintenance of anaerobiosis, at 28 °C, with continuous stirring (500 rpm). Fermentation experiments were carried out in 2X SD minimal medium (13.4% yeast nitrogen base without amino acids, 100% glucose). The SD minimal medium was supplemented with 1.25 mg/l ergosterol, 0.164 g/l Tween 80, and 0.35 mg/l oleic acid, to satisfy the lipid requirements of the yeast cells during anaerobic growth. Acetoin was added at various concentrations (from 0 to 300 mM) at the start of fermentation, as previously described by Celton et al.. Bioreactors were inoculated to a final OD600 of 0.05, with cells from a YPD preculture in 50 ml of medium in a 250 ml flask, which had been incubated at 28 °C, with shaking. Fermentations were performed at 28 °C. CO2 release was determined by automatic measurements of fermentor weight loss at 20-minute intervals . Fermentation experiments were performed in triplicate.
Yeast growth was monitored by measuring the optical density of the culture at 600 nm (OD600). Glucose and acetate were determined by high-performance liquid chromatography (HPLC) with an Aminex HPX-87 H ion exchange column. Volatile compounds were determined by GC-FID analysis, with a headspace autosampler and a BP20 (SGE).
Gene expression was analyzed from three independent cultures of each studied strain: 59A, 59A NADH-Bdh 200 mM, 59A NADPH-Bdh 100 mM, 200 mM and 300 mM. Aliquots of 1x109 cells were sampled from the fermentation medium when CO2 production reached 6 g/l/h. Cells were collected by centrifugation at 1000 g for 5 min at 4 °C and the cell pellets were washed with DEPC-treated water and frozen in methanol at −80 °C. Total RNA was extracted with Trizol reagent (Gibco BRL, Life Technologies)  and purified with the "RNA Cleanup kit" (Qiagen). We checked the quantity and quality of the extracted RNA by spectrometry (NanoDrop 1000, Thermo Scientific). The DNA microarrays were manufactured at the Biochip platform of Toulouse–Genopole, on dendrislides , with 70-mer oligonucleotides covering 99% of the yeast genome (Operon Inc.; a list of the corresponding genes is available from http://biopuce.insatoulouse.fr/oligosets/). Fluorescent cDNAs were synthesized from 5 μg of total RNA with the CleanUP System (Promega). Labeled cDNA was purified with the Pronto Purification kit (Promega). Microarrays were hybridized in static conditions for 16 h in 42 °C in a hybridization chamber (Corning), with the Pronto Universal Hybridization Quick kit.
Transcriptomic data acquisition and statistical analysis
Hybridization signals were detected with a GenePix 4000B laser scanner (Axon Instruments). Array images were acquired and quantified with integrated GenePix software version 3.0. Statistical analysis was done with R.2.11.1 software . The limma package [41–43] was used to normalize the microarray data (by the print-tip-loess method for normalization within the array followed by the quantile method for normalization between arrays). Differential expression between two sets of experimental conditions was detected by carrying out a modified t-test, with a significant threshold of p < 5 x 10-2, applying the Benjamini and Hochberg false discovery rate to correct for multiple testing . Only genes displaying a fold-change of at least 1.5 (positive or negative) in expression level were considered. Genes displaying differential expression were grouped according to Gene Ontology (GO) process terms, with the database for annotation, visualization and integrated discovery (DAVID), beta version 6.7 . The complete data set is available through the Gene Expression Omnibus database (accession number GSE34810).
Extraction and determination of intracellular metabolites
Sampling from two independent replicates of each studied strain (NADPH-Bdh 100 mM, 200 mM and 300 mM acetoin) and carried out when CO2 production reached 6 g/l/h. Cells (0.5 to 1.95 x 108) were added to 5 ml of methanol solution (40% vol/vol) to quench their metabolism. They were then collected by centrifugation (15000 g, 1 min, -10 °C) and stored at −80 °C. Extraction and quantitative analysis by isotope dilution mass spectrometry using uniformly 13 C-labeled cell extracts as internal standards  were carried out at the MetaToul Platform, Toulouse. Ethanol (5 μl of a 75% v/v solution) was added to the cell pellets and metabolites were extracted by incubating the mixture for three minutes at 80 °C. The samples were cooled on ice and centrifuged at 15000 g for 10 minutes. Ethanol was removed by evaporation and the metabolites were resuspended in 250 μl ultra-pure H2O and centrifuged at 20000 g for two minutes. Intracellular metabolites were analyzed by high-pressure anion exchange chromatography in a Dionex ICS 2000 machine (Dionex, USA) equipped with an AS11 analytical column (2 x 50 mm). Mass spectra were acquired with a QTrap 4000 electrochemical detector (Applied Biosystems) equipped with a Sciex Turbo V MSD (Toronto) for ESI electrospray ionization and MRM (Multiple Reaction Monitoring) in detector mode.
13 C-flux analysis
Cells were grown in penicillin vials filled with 10 ml 2X SD minimal medium supplemented with 1.25 mg/l ergosterol, 0.164 g/l Tween 80, 0.35 mg/l oleic acid and with 100 mg.l-1 of a mixture of labeled (40% 1-13 C-glucose) and unlabeled glucose and 1.4 mg.l-1 of NH4Cl as the sole nitrogen source. At an OD600 of 3.0, corresponding to mid-exponential growth phase, the concentration of extracellular metabolites (glucose, glycerol, organic acids) was measured in the supernatant (these data will be used to constrain the stoichiometric model, as described below). Cells were harvested and hydrolyzed overnight with 6 M HCl for the determination of amino acids labeling patterns. One glucose derivative (glucose pentacetate) and two amino-acid derivatives (ethyl chloroformate (ECF) and dimethyl formamide dimethyl acetal (DMFDMA)) were analyzed by GC-MS , as previously described by Gombert et al. . From the raw GC-MS data, the summed fractional labeling (SFL) of each fragment was calculated as follow: SFL = 100 × [(1.m 1 + 2.m 2 + …. + n.m n) × (m 0 + m 1 + m 2 + … + m n)-1, with m o the fractional abundance of the lowest corrected mass and mi > 0 the abundance of molecules with higher corrected masses. The labelling data are available in request.
The data set used for 13 C-flux analysis included the 25 SFLs calculated from labeling experiments and 22 measured fluxes, including the drain of metabolic intermediates to biomass and the formation of seven metabolites, giving a total of 47 items of experimental data. The distribution of carbon in central carbon metabolism (CCM) was estimated with the metabolic model described by Gombert et al., modified to take into account features specific to fermentation . The network used is described in additional file 2. Flux calculations (60 reactions) were carried out with Matlab 7, as previously described [47,48]. Differences between experimental and simulated SFLs and between experimental and simulated fluxes were minimized by an iterative procedure. Calculations were carried out 100 times and the fluxes reported are the means of the convergent solutions (Additional file 2).
We thank Lindsay Peyriga of the MetaSys (Functional Analysis of Metabolic Systems Platform), Toulouse, France, for the extraction and analysis of intracellular metabolites. We also thank Pierre Delobel for assistance with microarray analysis. This study was supported by the INRA AgroBi 2006–2008 program and by the BIOFLAVOUR Cost Action FA0907.
- Nielsen J: It is all about metabolic fluxes. J Bacteriol. 2003, 185: 7031-10.1128/JB.185.24.7031-7035.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Bro C, Regenberg B, Forster J, Nielsen J: In silico aided metabolic engineering of Saccharomyces cerevisiae for improved bioethanol production. Metab Eng. 2006, 8: 102-111. 10.1016/j.ymben.2005.09.007.View ArticlePubMedGoogle Scholar
- Heux S, Cadiere A, Dequin S: Glucose utilization of strains lacking PGI1 and expressing a transhydrogenase suggests differences in the pentose phosphate capacity among Saccharomyces cerevisiae strains. FEMS Yeast Res. 2008, 8: 217-224. 10.1111/j.1567-1364.2007.00330.x.View ArticlePubMedGoogle Scholar
- Michnick S, Roustan JL, Remize F, Barre P, Dequin S: Modulation of glycerol and ethanol yields during alcoholic fermentation in Saccharomyces cerevisiae strains overexpressed or disrupted for GPD1 encoding glycerol 3 phosphate dehydrogenase. Yeast. 1997, 13: 783-793. 10.1002/(SICI)1097-0061(199707)13:9<783::AID-YEA128>3.0.CO;2-W.View ArticlePubMedGoogle Scholar
- Celton M, Goelzer A, Camarasa C, Fromion V, Dequin S: A constraint-based model analysis of the metabolic consequences of increased NADPH oxidation in Saccharomyces cerevisiae. Metab Eng. 2012, 14: 366-79. 10.1016/j.ymben.2012.03.008.View ArticlePubMedGoogle Scholar
- Vaseghi S, Baumeister A, Rizzi M, Reuss M: In vivo dynamics of the pentose phosphate pathway in Saccharomyces cerevisiae. Metab Eng. 1999, 1: 128-140. 10.1006/mben.1998.0110.View ArticlePubMedGoogle Scholar
- Larochelle M, Drouin S, Robert F, Turcotte B: Oxidative stress-activated zinc cluster protein Stb5 has dual activator/repressor functions required for pentose phosphate pathway regulation and NADPH production. Mol Cell Biol. 2006, 26: 6690-6701. 10.1128/MCB.02450-05.PubMed CentralView ArticlePubMedGoogle Scholar
- Akache B, Wu K, Turcotte B: Phenotypic analysis of genes encoding yeast zinc cluster proteins. Nucl Acids Res. 2001, 29: 2181-10.1093/nar/29.10.2181.PubMed CentralView ArticlePubMedGoogle Scholar
- Hector RE, Bowman MJ, Skory CD, Cotta MA: The Saccharomyces cerevisiae YMR315W gene encodes an NADP(H)-specific oxidoreductase regulated by the transcription factor Stb5p in response to NADPH limitation. New biotechnol. 2009, 26: 171-180. 10.1016/j.nbt.2009.08.008.View ArticleGoogle Scholar
- Matsufuji Y, Nakagawa T, Fujimura S, Tani A, Nakagawa J: Transcription factor Stb5p is essential for acetaldehyde tolerance in Saccharomyces cerevisiae. J Basic Microbiol. 2010, 50: 494-498. 10.1002/jobm.200900391.View ArticlePubMedGoogle Scholar
- Cadière A, Galeote V, Dequin S: The Saccharomyces cerevisiae zinc factor protein Stb5p is required as a basal regulator of the pentose phosphate pathway. FEMS Yeast Res. 2010, 10: 819-827. 10.1111/j.1567-1364.2010.00672.x.View ArticlePubMedGoogle Scholar
- Regenberg B, Grotkjær T, Winther O, Fausbøll A, Åkesson M, Bro C, Hansen LK, Brunak S, Nielsen J: Growth-rate regulated genes have profound impact on interpretation of transcriptome profiling in Saccharomyces cerevisiae. Genome Biol. 2006, 7: R107-10.1186/gb-2006-7-11-r107.PubMed CentralView ArticlePubMedGoogle Scholar
- De Smidt O, Du Preez JC, Albertyn J: The alcohol dehydrogenases of Saccharomyces cerevisiae: a comprehensive review. FEMS Yeast Res. 2008, 8: 967-978. 10.1111/j.1567-1364.2008.00387.x.View ArticlePubMedGoogle Scholar
- Blomberg A: Metabolic surprises in Saccharomyces cerevisiae during adaptation to saline conditions: questions, some answers and a model. FEMS Microbiol Lett. 2000, 182: 1-8. 10.1111/j.1574-6968.2000.tb08864.x.View ArticlePubMedGoogle Scholar
- Clasquin MF, Melamud E, Singer A, Gooding JR, Xu X, Dong A, Cui H, Campagna SR, Savchenko A, Yakunin AF: Riboneogenesis in Yeast. Cell. 2011, 145: 969-980. 10.1016/j.cell.2011.05.022.PubMed CentralView ArticlePubMedGoogle Scholar
- Vargas RC, Garcia-Salcedo R, Tenreiro S, Teixeira MC, Fernandes AR, Ramos J, Sa-Correia I: Saccharomyces cerevisiae multidrug resistance transporter Qdr2 is implicated in potassium uptake, providing a physiological advantage to quinidine-stressed cells. Euk Cell. 2007, 6: 134-10.1128/EC.00290-06.View ArticleGoogle Scholar
- Larroy C, Parés X, Biosca JA: Characterization of a Saccharomyces cerevisiae NADP(H)‒dependent alcohol dehydrogenase (ADHVII), a member of the cinnamyl alcohol dehydrogenase family. Eur J Biochem. 2002, 269: 5738-5745. 10.1046/j.1432-1033.2002.03296.x.View ArticlePubMedGoogle Scholar
- Heer D, Heine D, Sauer U: Resistance of Saccharomyces cerevisiae to high concentrations of furfural is based on NADPH-dependent reduction by at least two oxireductases. Appl Environ Microbiol. 2009, 75: 7631-10.1128/AEM.01649-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu ZL, Moon J: A novel NADPH-dependent aldehyde reductase gene from Saccharomyces cerevisiae NRRL Y-12632 involved in the detoxification of aldehyde inhibitors derived from lignocellulosic biomass conversion. Gene. 2009, 446: 1-10. 10.1016/j.gene.2009.06.018.View ArticlePubMedGoogle Scholar
- Krüger A, Grüning NM, Wamelink MMC, Kerick M, Kirpy A, Parkhomchuk D, Bluemlein K, Schweiger MR, Soldatov A, Lehrach H: The pentose phosphate pathway is a metabolic redox sensor and regulates transcription during the anti-oxidant response. Antioxid Redox Signal. 2011, 15: 311-324. 10.1089/ars.2010.3797.View ArticlePubMedGoogle Scholar
- Sonderegger M, Jeppsson M, Hahn-Hagerdal B, Sauer U: Molecular basis for anaerobic growth of Saccharomyces cerevisiae on xylose, investigated by global gene expression and metabolic flux analysis. Appl Environ Microbiol. 2004, 70: 2307-10.1128/AEM.70.4.2307-2317.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Daran-Lapujade P, Jansen MLA, Daran JM, van Gulik W, de Winde JH, Pronk JT: Role of transcriptional regulation in controlling fluxes in central carbon metabolism of Saccharomyces cerevisiae. J Biol Chem. 2004, 279: 9125-9138. 10.1074/jbc.M309578200.View ArticlePubMedGoogle Scholar
- Wessely F, Bartl M, Guthke R, Li P, Schuster S, Kaleta C: Optimal regulatory strategies for metabolic pathways in Escherichia coli depending on protein costs. Mol Syst Biol. 2011, 7: 515-PubMed CentralView ArticlePubMedGoogle Scholar
- Lu P, Vogel C, Wang R, Yao X, Marcotte EM: Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol. 2006, 25: 117-124.View ArticlePubMedGoogle Scholar
- Kaiser P, Su NY, Yen J, Ouni I, Flick K: The yeast ubiquitin ligase SCFMet30: connecting environmental and intracellular conditions to cell division. Cell Div. 2006, 1: 16-10.1186/1747-1028-1-16.PubMed CentralView ArticlePubMedGoogle Scholar
- Thomas D, Cherest H, Surdin-Kerjan Y: Identification of the structural gene for glucose-6-phosphate dehydrogenase in yeast. Inactivation leads to a nutritional requirement for organic sulfur. EMBO J. 1991, 10: 547-PubMed CentralPubMedGoogle Scholar
- Hansen J, Johannesen PF: Cysteine is essential for transcriptional regulation of the sulfur assimilation genes in Saccharomyces cerevisiae. Mol Gen Genet. 2000, 263: 535-542. 10.1007/s004380051199.View ArticlePubMedGoogle Scholar
- Slekar KH, Kosman DJ, Culotta VC: The yeast copper/zinc superoxide dismutase and the pentose phosphate pathway play overlapping roles in oxidative stress protection. J Biol Chem. 1996, 271: 28831-10.1074/jbc.271.46.28831.View ArticlePubMedGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99-10.1038/nature02800.PubMed CentralView ArticlePubMedGoogle Scholar
- Albers E, Laizé V, Blomberg A, Hohmann S, Gustafsson L: Ser3p (Yer081wp) and Ser33p (Yil074cp) are phosphoglycerate dehydrogenases in Saccharomyces cerevisiae. J Biol Chem. 2003, 278: 10264-10.1074/jbc.M211692200.View ArticlePubMedGoogle Scholar
- Verstrepen KJ, Moonjai N, Bauer FF, Derdelinckx G, Dufour JP, Winderickx J, Thevelein JM, Pretorius IS, Delvaux FR: Genetic regulation of ester synthesis in yeast: new facts, insights and implications for the brewer. Brewing Yeast Fermentation Performance. Edited by: Smart K. 2003, Blackwell Science, Oxford, 234-248. 2Google Scholar
- Heux S, Cachon R, Dequin S: Cofactor engineering in Saccharomyces cerevisiae: Expression of a H2O-forming NADH oxidase and impact on redox metabolism. Met Eng. 2006, 8: 303-314.32. 10.1016/j.ymben.2005.12.003.View ArticleGoogle Scholar
- Larsson C, Nilsson A, Blomberg A, Gustafsson L: Glycolytic flux is conditionally correlated with ATP concentration in Saccharomyces cerevisiae: a chemostat study under carbon-or nitrogen-limiting conditions. J Bacteriol. 1997, 179: 7243-PubMed CentralPubMedGoogle Scholar
- Rossell S, Van Der Weijden CC, Lindenbergh A, Van Tuijl A, Francke C, Bakker BM, Westerhoff HV: Unraveling the complexity of flux regulation: a new method demonstrated for nutrient starvation in Saccharomyces cerevisiae. Proc Nat Acad Sci. 2006, 103: 2166-10.1073/pnas.0509831103.PubMed CentralView ArticlePubMedGoogle Scholar
- Ambroset C, Petit M, Brion C, Sanchez I, Delobel P, Guérin C, Chiapello H, Nicolas P, Bigey F, Dequin S: Deciphering the molecular basis of wine yeast fermentation traits using a combined genetic and genomic approach. G3. 2011, 1: 263-281. 2011.PubMed CentralView ArticlePubMedGoogle Scholar
- Ehsani M, Fernández MR, Biosca JA, Dequin S: Reversal of coenzyme specificity of 2, 3-butanediol dehydrogenase from Saccharomyces cerevisae and in vivo functional analysis. Biotech Bioeng. 2009, 104: 381-389. 10.1002/bit.22391.View ArticleGoogle Scholar
- Sablayrolles JM, Barre P, Grenier P: Design of a laboratory automatic system for studying alcoholic fermentations in anisothermal enological conditions. Biotechnol Tech. 1987, 1: 181-184.38.View ArticleGoogle Scholar
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159.View ArticlePubMedGoogle Scholar
- Le Berre V, Trévisiol E, Dagkessamanskaia A, Sokol S, Caminade AM, Majoral JP, Meunier B, François J: Dendrimeric coating of glass slides for sensitive DNA microarrays analysis. Nucl Acids Res. 2003, 31: e88-10.1093/nar/gng088.PubMed CentralView ArticlePubMedGoogle Scholar
- R: A language and environment for statistical computing. 2010, R Foundation for Statistical Computing, Vienna, Austria, URL R-project.org
- Smyth GK, Michaud J, Scott HS: Use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2067, 2005: 21-Google Scholar
- Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Appl Gen Mol Biol. 2004, 3: 3-Google Scholar
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate in behavior genetics research. Behavioural Brain Res. 2001, 125: 289-300.View ArticleGoogle Scholar
- Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC: DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucl Acids Res. 2007, 35 (suppl2): W169-PubMed CentralView ArticlePubMedGoogle Scholar
- Wu L, Mashego MR, van Dam JC, Proell AM, Vinke JL, Ras C, van Winden WA, van Gulik WM, Heijnen JJ: Quantitative analysis of the microbial metabolome by isotope dilution mass spectrometry using uniformly 13 C-labeled cell extracts as internal standards. Anal Biochem. 2005, 336: 164-171. 10.1016/j.ab.2004.09.001.View ArticlePubMedGoogle Scholar
- Christensen B, Nielsen J: Isotopomer analysis using GC-MS. Met Eng. 1999, 1: 282-290. 10.1006/mben.1999.0117.View ArticleGoogle Scholar
- Gombert AK, dos Moreira Santos M, Christensen B, Nielsen J: Network identification and flux quantification in the central metabolism of Saccharomyces cerevisiae under different conditions of glucose repression. J Bacteriol. 2001, 183: 1441-10.1128/JB.183.4.1441-1451.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Heux S: Ingénierie métabolique et analyse13C-flux du métabolisme central des levuresSaccharomyces cerevisiaeœnologiques.PhD Thesis. 2006, University of Burgundy, Sciences de la vie et de la santéGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.