- Open Access
Large-scale 13C-flux analysis reveals mechanistic principles of metabolic network robustness to null mutations in yeast
Genome Biology volume 6, Article number: R49 (2005)
Quantification of intracellular metabolite fluxes by 13C-tracer experiments is maturing into a routine higher-throughput analysis. The question now arises as to which mutants should be analyzed. Here we identify key experiments in a systems biology approach with a genome-scale model of Saccharomyces cerevisiae metabolism, thereby reducing the workload for experimental network analyses and functional genomics.
Genome-scale 13C flux analysis revealed that about half of the 745 biochemical reactions were active during growth on glucose, but that alternative pathways exist for only 51 gene-encoded reactions with significant flux. These flexible reactions identified in silico are key targets for experimental flux analysis, and we present the first large-scale metabolic flux data for yeast, covering half of these mutants during growth on glucose. The metabolic lesions were often counteracted by flux rerouting, but knockout of cofactor-dependent reactions, as in the adh1, ald6, cox5A, fum1, mdh1, pda1, and zwf1 mutations, caused flux responses in more distant parts of the network. By integrating computational analyses, flux data, and physiological phenotypes of all mutants in active reactions, we quantified the relative importance of 'genetic buffering' through alternative pathways and network redundancy through duplicate genes for genetic robustness of the network.
The apparent dispensability of knockout mutants with metabolic function is explained by gene inactivity under a particular condition in about half of the cases. For the remaining 207 viable mutants of active reactions, network redundancy through duplicate genes was the major (75%) and alternative pathways the minor (25%) molecular mechanism of genetic network robustness in S. cerevisiae.
The availability of annotated genomes and accumulated biochemical evidence for individual enzymes triggered the reconstruction of stoichiometric reaction models for network-based pathway analysis [1, 2]. For many microbes, such network models are available at the genome scale, providing a largely comprehensive metabolic skeleton by interconnecting all known reactions in a given organism [3, 4]. Thus, network properties such as optimal performance, flexibility to cope with ever-changing environmental conditions, and enzyme dispensability (also referred to as robustness or genetic robustness [5, 6]) become mathematically tractable. These computational advances are matched with post-genomic advances in experimental methods that assess the cell's molecular make-up at the level of mRNA, protein, or metabolite concentrations. As the functional complement to these compositional data, quantification of intracellular in vivo reaction rates or molecular fluxes has been a focal point of method development in the realm of metabolism [7–9]. Recent progress in increasing the throughput of stable-isotope-based flux analyses [8, 10, 11] has allowed the quantification of flux responses to more than just a few intuitively chosen genetic or environmental perturbations [12–14]. Now that flux quantification in hundreds of null mutants under a particular condition is feasible in principle, the question arises of which mutants should be analyzed.
As perhaps the most widely used model eukaryote, the yeast Saccharomyces cerevisiae features a metabolic network of about 1,200 reactions that represent about 750 biochemically distinct reactions [3, 15]. Is it necessary to quantify flux responses to null mutations in all reactions for a comprehensive view of the metabolic capabilities under a given condition? To address this question, we used a recently modified version (iLL672; L Kuepfer, U Sauer and LM Blank, unpublished work) of the original iFF708 genome-scale model published by Förster et al. . On the basis of this model, we estimated the genome-scale flux distribution in wild-type S. cerevisiae from 13C-tracer experiments, to identify the 339 biochemical reactions that were active during growth on glucose. Yeast metabolism has the potential flexibility to use alternative pathways for 105 of these active reactions. For a major fraction of the potentially flexible reactions that catalyze significant flux, we then constructed prototrophic knockout mutants to elucidate whether or not the alternative pathway was used upon experimental knockout; that is, whether it contributes to the genetic robustness of the network [5, 6]. For the purpose of this work, robustness is defined as the ability to proliferate on glucose as the sole carbon source upon knockout of a single gene with metabolic function.
Identification of flexible reactions in yeast metabolism
To identify all potentially flexible reactions in yeast glucose metabolism that were active under a given condition, we used the recently reconciled metabolic network model iLL672 with 1,038 reactions (encoded by 672 genes) that represent 745 biochemically distinct reactions (L Kuepfer, U Sauer and LM Blank, unpublished work), which was based on the genome-scale S. cerevisiae model iFF708 . The main modifications to the original model include elimination of dead-end reactions and a new formulation of cell growth. It should be noted that none of the results below critically depended on the network model, but the reconciliation of iLL672 enabled a more accurate discrimination between lethal and viable reactions than iFF708, as was validated by large-scale growth experiments (L Kuepfer, U Sauer and LM Blank, unpublished work).
First, we identified all reactions active in wild-type glucose metabolism by genome-scale flux analysis. For this purpose, we determined the wild-type flux distribution in central metabolism from a stable isotope batch experiment with 20% [U-13C] and 80% unlabeled glucose. This flux solution was then mapped to the genome scale by using minimization of the Euclidean norm of fluxes as the objective function. In total, 339 of the 745 biochemical reactions were active during growth on glucose alone (Figure 1 and Additional data file 1), which agrees qualitatively with the estimate of Papp et al. . Most active reactions (234) were essential: 155 are encoded by singleton genes, 64 by two or more duplicate genes and 15 by yet unknown genes (Figure 1; Additional data file 1). In the entire network, only the remaining 105 reactions (30 encoded by yet unknown genes) were active and potentially flexible in the sense that they may be bypassed via alternative pathways (Figure 1). As fluxes in the peripheral reactions were typically below 0.1% of the glucose uptake rate (see Additional data file 1), we focused on the 51 gene-encoded flexible reactions that catalyzed a flux of at least 0.1%. These 51 reactions were encoded by 75 genes (43 duplicates, 23 singletons and 9 multiprotein complexes).
Physiological fitness of mutants deleted in flexible reactions
In 38 of these genes, which encoded 28 of the 51 potentially flexible and highly active reactions, we constructed prototrophic deletion mutants by homologous recombination  in the physiological model strain CEN.PK  (Figure 2). The prototrophic background was chosen to minimize potential problems of amino-acid supplementation for quantitative analysis . These 38 experimental knockouts were in the pentose phosphate (PP) pathway, tricarboxylic acid (TCA) cycle, glyoxylate cycle, polysaccharide synthesis, mitochondrial transporters, and by-product formation (Figure 2, Table 1). Genetically, the knockouts encompass 14 singleton and 24 duplicate genes, including six gene families of which all members were deleted.
With the exception of gnd1, all 38 mutants grew with glucose as the sole carbon source. The lethal phenotype of the gnd1 mutant is consistent with a previous report  and is similar to the gndA mutant in Bacillus subtilis . As in B. subtilis, we could select gnd1 suppressor mutants on glucose (data not shown). To assess the quantitative contribution of each gene to the organism's fitness, we determined maximum specific growth rates in minimal and complex medium using a well-aerated microtiter plate system . Mutant fitness was then expressed as the normalized growth rate, relative to the growth rate of the reference strain (Table 1). In contrast to the previously reported competitive fitness [20, 23, 24], the fitness determined here is a quantitative physiological value.
In complex YPD medium, physiological fitness in the 38 viable haploid mutants was generally in qualitative agreement with the competitive fitness . Quantitatively, however, our data seem to allow a better discrimination (Table 1), and significant differences between physiological and competitive fitness were seen in the adh1, fum1, and gpd1 mutants. Only threemutants - adh1, fum1, and gly1 - exhibited a fitness defect of 20% or greater (Table 2). gly1 lacks threonine aldolase, which catalyzes cleavage of threonine to glycine , hence its phenotype remains unexplained because glycine was present in the YPD medium.
In general, growth on the single substrate reduced the metabolic flexibility, as a much greater proportion of mutants exhibited significant fitness defects (Table 2). Major fitness defects were prominent in mutants of the PP pathway (gnd1, rpe1, sol3, and zwf1), which indicates an increased demand of NADPH for biosynthesis. Fitness of the fum1 mutant was clearly lower than that of other TCA-cycle mutants, for which duplicate genes exist. The strong phenotype of the fum1 mutant was somewhat unexpected because the flux through the TCA cycle is generally low or absent in glucose batch cultures of S. cerevisiae [13, 14, 26, 27].
Intracellular carbon flux redistribution in response to gene deletions
While physiological data quantify the fitness defect, they cannot differentiate between intracellular mechanisms that bring about robustness to the deletion. To identify how carbon flux was redistributed around a metabolic lesion, we used metabolic flux analysis based on 13C-glucose experiments [8, 9]. In contrast to in vitro enzyme activities and expression data, 13C-flux analysis provides direct evidence for such in vivo flux rerouting or its absence. The flux protocol consists of two distinct steps: first, analytical identification of seven independent metabolic flux ratios with probabilistic equations from the 13C distribution in proteinogenic amino acids [12, 28, 29]; and second, estimation of absolute fluxes (in vivo reaction rates) from physiological data and the flux ratios as constraints [10, 30]. The relative distribution of intracellular fluxes was rather invariant in the 37 mutants, with the fraction of mitochondrial oxaloacetate derived through the TCA cycle flux and the fraction of mitochondrial pyruvate originating from malate as prominent exceptions (Figure 3).
From the experimentally determined uptake/production rates and the flux ratios as constraints (Additional data file 3), absolute intracellular fluxes were calculated using a compartmentalized stoichiometric model that consists of 35 reactions and 30 metabolites (Additional data file 2). This flux model comprised mostly the reactions of central carbon metabolism that were most relevant to the genetic changes introduced. It should be noted that the deleted reactions, with the exception of pyruvate dehydrogenase (PDA1), were not omitted from the network model; thus the calculated absence of flux through a given reaction was independently verified from the 13C-labeling data. In contrast to the relative distribution of intracellular fluxes, absolute reaction rates varied significantly in the mutants. With the exception of the flux through the TCA cycle (Figure 4f) and the gluconeogenic PEP carboxykinase (Figure 4d), all other fluxes generally increased with increasing glucose uptake rates (Figure 4). Eleven of the 37 mutants, however, exhibited specific flux responses that deviated from this general trend (Table 2, Figure 4).
Specific flux responses in singleton and duplicate gene knockouts
Specific flux responses were more prominent among the singleton mutants (Table 2, Figure 4). Although the TCA cycle flux through the NAD+-dependent fumarase reaction from fumarate to malate was already very low in the reference strain (Figures 3, 4f), the fum1 mutant exhibited a pronounced phenotype with altered redox metabolism and significant glycerol production (Figure 5). Inactivation of the mitochondrial pyruvate dehydrogenase complex in the pda1 mutant was bypassed by the import of cytosolic acetyl-CoA into the mitochondria. Inactivation of the oxidative PP pathway branch in the zwf1 mutant was compensated by a reversed flux in the non-oxidative PP pathway to provide the biomass precursors pentose 5-phosphate and erythrose 4-phosphate (Figure 5). Because the primary role of the PP pathway on glucose is generation of NADPH, NADP+-dependent mitochondrial malic enzyme flux was significantly increased in the zwf1 mutant. This NADPH compensation by malic enzyme was also suggested recently from co-feeding experiments .
In contrast to singletons, deletion of flexible duplicate genes could be compensated by either alternative pathways or isoenzymes. In most cases, however, the isoenzymes were used because no flux alteration was detected, with the a dh1, ald6, cox5A, and mdh1 mutants as exceptions (Table 2). Deletion of the major acetate-producing acetaldehyde dehydrogenase, the cytoplasmic ALD6 , significantly reduced acetate formation. The primary effect of the deletion was the strongly reduced glucose-uptake rate (Figure 4). Although a major source of NADPH was inactivated in this mutant , the PP pathway flux was not increased, but was even lower than in the reference strain (Figure 6). This indicates that the strongly decreased fitness of the ald6 mutant (Table 1) could result from NADPH starvation - that is, a suboptimal rate of NADP+ reduction. Consistent with this, we estimated that the NADPH requirement exceeded the combined NADPH formation from the oxidative PP pathway and malic enzyme by 70%, indicating that an as-yet-unidentified reaction(s) substitutes for the remaining NADPH production. Candidates are the mitochondrial acetaldehyde dehydrogenase Ald4p , which can use either NAD+ or NADP+ as redox cofactors or the mitochondrial NADH kinase Pos5p . Deletion of the cytochrome c oxidase subunit Va COX5A in the mitochondrial respiratory chain increased glycerol production, which serves as means to reoxidize NADH (Figures 4b, 6). Because this mutant lacks functional mitochondria , glycerol production was driven by the limited NADH reoxidation through residual NADH oxidase activity in the electron-transport chain. Thus, robustness was brought about by using an alternative NADH sink. Considering that the flux through the mitochondrial malate dehydrogenase Mdh1 was already very low in the reference strain, the fitness defect of the mdh1 was surprising. Akin to the fum1 and ald6 mutants, the significantly reduced fitness of mdh1 may be explained by the imbalance between the TCA cycle and glucose catabolism (Figure 4f). Generally, the TCA cycle flux increases with decreasing glucose uptake rates , but remains non-proportionally low (absent) in the fum1, ald6, and mdh1 mutants (Figure 4f). The cytosolic and peroxisomal duplicate genes MDH2 and MDH3, respectively, did not compensate for the mitochondrial lesion, which is consistent with the observed lethal phenotype of the mdh1 mutant when grown on acetate .
Genetic network robustness
The above flux results reveal that knockouts of flexible reactions are bypassed through alternative pathways in about one third of the cases and through isoenzymes in the other two thirds. Does this reflect the relative contribution of alternative pathways and duplicate genes to genetic network robustness?  To address this question quantitatively for glucose metabolism, we grew the 196 duplicate (encoding 87 reactions) and 171 singleton (encoding 207 reactions) knockout mutants of all 294 gene-encoded active reactions on glucose plates.
In the 47 viable singleton knockouts, flux rerouting through an alternative pathway ensures survival, which was directly verified by flux data in 10 cases (Figure 4, Table 3 and Additional data file 3). Of the 196 experimental duplicate mutants, 180 grew on sole glucose, while 16 of the mutations were lethal. As these 16 duplicates obviously did not contribute to genetic robustness, their entire families (36 genes) were subtracted from the 150 duplicate-encoded essential reactions (Figure 1). For the remaining 114 duplicate genes we have strong evidence for network redundancy as the underlying mechanism of robustness, because they encode essential reactions (as determined in silico) and each of the experimental knockouts was viable (Figure 7). For the 46 duplicate genes that encode flexible reactions (Figure 1), both compensation by duplicates and/or alternative pathways might ensure proliferation. Where available, these mutants were classified according to their flux distribution; that is, of the 24 experimental duplicate mutants analyzed, four used alternative pathways and 20 an isoenzyme (Figure 4, Table 3 and Additional data file 3). In total we analyzed all 367 experimental mutants that encode the 294 active reactions of glucose metabolism, 140 of which were lethal and 227 viable. For the vast majority of the viable mutants, we identified the molecular mechanism that brought robustness to the knockout about: about 25% were alternative pathways and 75% duplicate genes (Figure 7).
Using an integrated computational and experimental approach, we show here that metabolic flexibility to knockout mutations is restricted to a relatively small set of biochemical reactions. About a third of all active reactions under the particular condition investigated may be bypassed by alternative pathways, of which about 30% support only negligible fluxes. The occurrence of flexible reactions might be even lower in prokaryotes, because several alternative pathways involved inter-compartmental transport. In general, the number of flexible reactions will differ substantially between species, with free-living yeast and fungi at the upper end of the scale, and intracellular pathogens with highly reduced genomes at the lower end.
Using flux balance (FBA) [1, 2], elementary flux mode [38, 39], or similar analyses , all in silico flexible reactions can be precisely identified. Hence, experimental analysis of intracellular flux responses to metabolic gene deletions can be limited to these potentially flexible mutants, rather than having to analyze the entire mutant collection. Using the systems biological approach described here, the true in vivo capability of metabolic network operation can be mapped with a reasonable workload. As the knowledge base on intracellular flux responses increases, a handful of flux experiments in computationally identified mutants will probably suffice to identify the in vivo network capability under a given condition. At the next level, such flux analyses will also include mutants affected in regulatory genes that modulate the network composition. Although not covered in the stoichiometric models employed for flux-balance analysis, several recently discussed computational approaches [39, 41, 42] may aid in identifying the most relevant regulatory mutants for in vivo flux quantification.
Consistent with the notion that metabolic networks undergo minimal flux redistributions with respect to the metabolic state of the parent , deletion of flexible singleton genes was mostly counteracted by local flux rerouting, for example, the lsc1,mae1, and oac1 mutants (see Additional data file 3). Deletions in redox cofactor-dependent singleton or duplicate reactions such as those mediated by adh1, ald6, cox5A, pda1, and zwf1, however, affected flux alterations in more distant reactions. While the relative flux distribution (in % values) was perturbed only very little in these mutants, the absolute magnitude of fluxes (in vivo reaction velocities) varied dramatically. In particular, knockout of fum1, whose encoded protein catalyzes only a rather small flux, led to an unexpectedly strong phenotype with about a 50% reduction in glucose-uptake rate. Although unexpected, this finding was qualitatively consistent with results from a recent genetic footprinting study , which also showed a significant fitness defect in this mutant. It was speculated that intramitochondrial shortage of amino acids such as aspartate and glutamate causes a lack of respiratory chain components, which leads to a petite-like phenotype . Another key mutant was pda1, whose knockout caused a substantial import of acetyl-CoA into the mitochondria; the mechanism for this remains elusive because the carnitine auxotrophic CEN.PK strain does not use the carnitine shuttle . As a consequence, a twofold overproduction of NADPH was estimated, which suggests that the NAD+-dependent acetaldehyde dehydrogenases instead of the NADP+-dependent ALD6 were active to balance NADPH formation/consumption. Consistent with this, the flux through the NADPH-producing PP pathway was significantly lower in this mutant. The strongly altered redox metabolism in pda1 is further evidenced by the substantial secretion of glycerol and succinate (Figure 5).
The metabolic flexibility to cope with metabolic lesions is generally known as genetic robustness , a concept that is used to explain the seemingly surprising number of phenotypically silent deletion mutations: only about 1,100 knockouts of the 5,700 genes are lethal in haploid S. cerevisiae [23, 46]. The causes and evolution of gene dispensability have been investigated in several theoretical analyses of pre-existing data, but the issue remains controversial [5, 6, 16, 47–49]. For metabolic networks, our flux data differentiate between the relative contributions of three mechanisms to the apparent genetic robustness: inactive, and thus dispensable, genes; 'genetic buffering' through alternative reactions; and functional complementation from duplicate genes ('redundancy').
In qualitative agreement with a recent estimate , genome-scale flux analysis revealed that about half of the available reactions (45% of the known metabolic genes) were not required for growth on glucose (Figure 1). Hence, deletion of these genes would not affect the growth phenotype on this substrate, making inactive reactions the primary reason for the apparent dispensability of genes with metabolic function. It should be noted that this apparent gene dispensability is a somewhat artificial classification that does not contribute to genetic robustness because most of these genes encode metabolic functions that are only relevant under conditions different from the one tested. The most important mechanism of true genetic robustness in yeast glucose metabolism was duplicate genes (Figure 7), the majority of which encoded essential reactions with no alternative pathway. Alternative pathways, contributed about 25% to genetic robustness by carbon flux rerouting. This leaves redundancy as the major, and modularity the minor, cause  of metabolic network robustness to single-gene deletions during growth on glucose.
Materials and methods
All prototrophic S. cerevisiae deletion mutants were constructed in the haploid, CEN.PK113-7D (Mat a MAL2-8cSUC2) background with the homolog flanking region approach  (Table 1). Briefly, genomic DNA was isolated from the corresponding amino-acid auxotrophic mutants . The kanMX4 cassettes of each mutant were amplified by PCR with primers located about 500 bp upstream and downstream of the deleted genes. The PCR reaction mixture was directly used for transformation and integrants were selected on YPD plates with 300 μg/ml geneticin. Correct cassette insertion was confirmed by overlapping PCR using either primer KanB (5'-CTGCAGCGAGGAGCCGTAAT-3') or KanC (5'-TGATTTTGATGACGAGCGTAA-3') primers in combination with one gene-specific primer. The reference strain was CEN.PK 113-7D with a deletion of the switching endonuclease, which was shown to be phenotypically neutral in chemostat competition experiments  and is commonly used as reference [52, 53].
Media and growth conditions
The composition of the yeast minimal medium (MM) was : 5 g (NH4)2SO4, 3 g KH2PO4, 0.5 g MgSO4·7H2O, 4.5 mg ZnSO4·7H2O, 0.3 mg CoCl2·6H2O, 1.0 mg MnCl2·4H2O, 0.3 mg CuSO4·5H2O, 4.5 mg CaCl2·2H2O, 3.0 mg FeSO4·7H2O, 0.4 mg NaMoO4·2H2O, 1.0 mg H3BO3, 0.1 mg KI, 15 mg EDTA, 0.05 mg biotin, 1.0 mg calcium pantothenate, 1.0 mg nicotinic acid, 25 mg inositol, 1.0 mg pyridoxine, 0.2 mg p-aminobenzoic acid, and 1.0 mg thiamine. The medium was buffered at pH 5.0 with 100 mM KH-phthalate to reduce pH changes throughout the growth experiments to 0.05. Filter-sterilized glucose and geneticin (300 μg/ml) were added freshly to the media. Batch growth experiments (1.2 ml) were carried out in deep-well plates (System Duetz, Kühner AG, Switzerland) using an orbital shaker with 5 cm amplitude at 300 rpm to allow optimal mixing .
Qualitative testing of mutant growth on glucose was done on agar plates. For this purpose, we used the haploid yeast mutant library in the BY4741 strain (MAT a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0) . The composition of the yeast minimal medium for the plate growth assay was as described above  with the exception of 20 g/l agar for solidification. Glucose was added to a final concentration of 20 g/l. Strain auxotrophies were complemented with 20 mg/l histidine, uracil, methionine, lysine and 60 mg/l leucine. The plates were incubated at 30°C for 3 days before scoring of the growth phenotype and further incubated for 1 week to score slow-growth phenotype mutants.
Analytical procedures and 13C-labeling experiments
Cell growth was monitored by following optical density changes at a wavelength of 600nm (OD600). Aliquots for extracellular metabolite analysis were centrifuged at 14,000 rpm in an Eppendorf tabletop centrifuge to remove cells. Glucose, acetate, ethanol and glycerol concentrations in the supernatant were determined with commercial enzymatic kits (Scil Diagnostics, Germany). Organic acids were quantified by high-pressure liquid chromatography (HPLC) using a Supelcogel C8 (4.6 by 250 mm) ion-exclusion column. The column was eluted at 30°C with 2% sulfuric acid at a flow rate of 0.3 ml/min. The organic acids were detected using a PerkinElmer UV detector (Series 2000) at a wavelength of 210 nm. The physiological parameters maximum specific growth rate, biomass yield on glucose, and specific glucose consumption rate were calculated during the exponential growth phase.
All labeling experiments were carried out in batch cultures assuming pseudosteady-state conditions during the exponential growth phase [12, 55]. 13C-labeling of proteinogenic amino acids was achieved either by growth on 5 g/l glucose as a mixture of 80% (w/w) unlabeled and 20% (w/w) uniformly labelled [U-13C]glucose (13C > 99%; Martek Biosciences, Columbia, MD) or 100% [1-13C]glucose (> 99%; Omicron Biochemicals, South Bend, IN). Cells from overnight cultures were harvested by centrifugation and washed using sugar-free MM to remove residual unlabeled carbon sources. Cultures were routinely inoculated to an maximum OD600 of 0.03 and harvested by centrifugation at an OD600 ≤ 1. Residual medium was removed by washing the pellet with water. Cell protein was hydrolyzed for 24 h at 105°C in 6 M HCL and dried in a heating block at 85°C for 6 h. The free amino acids were derivatized at 85°C for 1 h using 15 μl dimethylformamide and 15 μl N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide . Gas chromatography-mass spectrometry (GC-MS) analysis was carried out as reported  using a series 8000 GC in combination with an MD800 mass spectrometer (Fisons Instruments, Beverly, MA).
Metabolic flux ratio analysis
The recorded MS spectra include the distribution of mass isotopomers in 1-5 fragments of alanine, aspartate, glutamate, glycine, isoleucine, leucine, phenylalanine, proline, serine, threonine, tyrosine, and valine. For each amino-acid fragment α, a mass isotopomer distribution vector (MDV) was assigned:
with m0 being the fractional abundance of the lowest mass and mi>0 the abundances of molecules with higher masses. The MDV α values were corrected for naturally occurring stable isotopes  to obtain the exclusive mass isotopomer distributions of the carbon skeletons. The corrected MDV α were used to calculate the amino acids (MDVAA) and metabolites (MDVM) mass distribution vectors. Ratios of converging intracellular fluxes to a given metabolite were calculated from the MDVM as described previously [12, 29].
In addition, the relative contribution of the PP pathway was quantified from [1-13C]glucose experiments by tracking the positional 13C-labeling [10, 56]. The expected labeling pattern of triose phosphates or serine, which is derived exclusively through glycolysis, is 50% 13C-label in the C1 positions. Hence, the fraction of serine derived through the pentose phosphate (PP) pathway can be derived according to Equation 2 :
where GLU3unlabeled is an unlabeled three-carbon fragment from a source molecule of glucose. The remaining fraction of serine must then be derived through glycolysis. This flux ratio was not corrected for the potential withdrawal of 13C-label in dihydroxyacetone-phosphate-based biomass synthesis (such as phospholipids) and glycerol formation , because the influence was negligible under the condition used. The largest effect was found in the mutant with the highest specific glycerol formation rate (cox5A), where the estimated relative flux through the PP pathway would decrease from 12 ± 1% to 9 ± 1%.
13C-constrained flux analysis
Absolute values of intracellular fluxes were calculated with a flux model comprising all the major pathways of yeast central carbon metabolism (Additional data file 2). Deleted reactions were not omitted from the mutant models; thus the mutations were independently verified from the 13C data. The stoichiometric matrix of 34 linear equations and 30 metabolites has an infinite condition number ; it is thus underdetermined, and has a solution space with an infinite number of different flux vectors that fulfill the constraints from determined uptake and production rates. To uniquely solve the system for fluxes (ν), a set of linearly independent equations that quantify flux ratios (FlRs) were used to obtain eight constraints on the relative flux distribution from METAFoR analysis (see Additional data file 2).
The fraction of cytosolic oxaloacetate originating from cytosolic pyruvate is given by:
The fraction of mitochondrial oxaloacetate derived through anaplerosis is given by:
The fraction of PEP originating from cytosolic oxaloacetate is given by:
The fraction of serine derived through glycolysis is given by:
The upper and lower bounds for mitochondrial pyruvate derived through the malic enzyme (from mitochondrial malate) are given by:
The contribution of glycine to serine biosynthesis is given by:
and, finally, the contribution of serine to glycine biosynthesis is given by:
The stoichiometric matrix including Equations 3-10 has a condition number of 31, implying that the model is numerically robust . Error minimization was carried out as described by Fischer et al. . Balanced NADPH production and consumption were not added as additional constraints. In general, NADPH production was constrained by Equations 3 and 7/8, which estimate the relative use of the PP pathway and malic enzyme, respectively. As an additional source of NADPH, the flux through the NADPH-dependent acetaldehyde dehydrogenase  was estimated from the acetate production rate and the biomass requirement for cytosolic acetyl-CoA. Deviation of the NADPH production estimated in this way from the consumption for biosynthesis was generally below ± 20%, suggesting that the model assumptions and the experimental data are highly consistent. All extreme flux patterns were independently verified in 30-ml cultures (data not shown).
Genome-scale flux analysis
We used the experimentally determined in vivo flux data (νexp) to constrain the purely stoichiometric solution space of model iLL672 to obtain an experimentally validated genome-scale wild-type flux solution νWT. For glucose minimal medium, we constrained the model iLL672 with 30 fluxes that were derived from 13C-labeling experiments . In particular, we used 13C-constrained flux analysis  for GC-MS-detected mass isotope distributions in proteinogenic amino acids from a 20% [U-13C] glucose experiment and a compartmentalized yeast model . These experimental data were to be kept within an accuracy δ of ± 10% when mapping the determined central metabolic fluxes to the genome-scale reference flux solution. To overcome mathematical artifacts such as futile cycling (that is, a closed loop of fluxes that bring no net change), the original linear programming problem was modified. A minimization of the Euclidean norm of fluxes was chosen as the objective function such that (s.t.) the mass balance equations hold:
with j as the set of experimentally determined fluxes. Reactions were categorized as flexible when fulfilling the following criteria: the reactions carried a non-zero flux; and the reaction was not essential for growth.
In silicophenotyping of duplicate gene families
Phenotype predictions of deletion mutants were analyzed computationally with FBA [3, 59]. Assuming steady-state growth, mass balances were put up for each intracellular metabolite M i (1 × n) that have to be fulfilled, when multiplied with the overall flux vector ν (n × 1):
M i ·ν = 0. (12)
The entity of all m metabolite mass balances yields the stoichiometric matrix S (m × n), where:
S·ν = 0. (13)
To pick one solution out of the overall solution space formed by the stoichiometric constraints, FBA generally assumes maximization of biomass growth μ as the global cellular goal [3, 59]. Thus, the search for a single flux distribution ν results in a linear programming (LP) problem:
where i = 1,..,M and νlb,iand νub,icorrespond to upper and lower bounds of a specific reaction i. Gene knockout mutants can be simulated easily in silico by setting the deleted reactions to zero. All LPproblems were solved using the open-source GNU linear programming kit .
Additional data files
The following additional data are available with the online version of this paper. The classification of reactions according to Figure 1 is presented in Additional data file 1. The flux analysis model is defined in Additional data file 2. The physiological data, flux ratios and the calculated flux distributions are presented in Additional data file 3.
Papin JA, Stelling J, Price ND, Klamt S, Schuster S, Palsson BO: Comparison of network-based pathway analysis methods. Trends Biotechnol. 2004, 22: 400-405. 10.1016/j.tibtech.2004.06.010.
Price ND, Reed JL, Palsson BO: Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol. 2004, 2: 886-897. 10.1038/nrmicro1023.
Förster J, Famili I, Fu P, Palsson BO, Nielsen J: Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. Genome Res. 2003, 13: 244-253. 10.1101/gr.234503.
Reed JL, Palsson BO: Thirteen years of building constraint-based in silico models of Escherichia coli. J Bacteriol. 2003, 185: 2692-2699. 10.1128/JB.185.9.2692-2699.2003.
Gu X: Evolution of duplicate genes versus genetic robustness against null mutations. Trends Genet. 2003, 19: 354-356. 10.1016/S0168-9525(03)00139-2.
Hartman JLT, Garvik B, Hartwell L: Principles for the buffering of genetic variation. Science. 2001, 291: 1001-1004. 10.1126/science.291.5506.1001.
Hellerstein MK: In vivo measurement of fluxes through metabolic pathways: the missing link in functional genomics and pharmaceutical research. Annu Rev Nutr. 2003, 23: 379-402. 10.1146/annurev.nutr.23.011702.073045.
Sauer U: High-throughput phenomics: experimental methods for mapping fluxomes. Curr Opin Biotechnol. 2004, 15: 58-63. 10.1016/j.copbio.2003.11.001.
Wiechert W: 13C metabolic flux analysis. Metab Eng. 2001, 3: 195-206. 10.1006/mben.2001.0187.
Fischer E, Zamboni N, Sauer U: High-throughput metabolic flux analysis based on gas chromatography-mass spectrometry derived 13C constraints. Anal Biochem. 2004, 325: 308-316. 10.1016/j.ab.2003.10.036.
Zamboni N, Sauer U: Model-independent fluxome profiling from 2H and 13C experiments for metabolic variant discrimination. Genome Biol. 2004, 5: R99-10.1186/gb-2004-5-12-r99.
Fischer E, Sauer U: Metabolic flux profiling of Escherichia coli mutants in central carbon metabolism using GC-MS. Eur J Biochem. 2003, 270: 880-891. 10.1046/j.1432-1033.2003.03448.x.
Blank LM, Sauer U: TCA cycle activity in Saccharomyces cerevisiae is a function of the environmentally determined specific growth and glucose uptake rates. Microbiology. 2004, 150: 1085-1093. 10.1099/mic.0.26845-0.
Blank LM, Lehmbeck F, Sauer U: Metabolic flux and network analysis in fourteen hemiascomycetous yeasts. FEMS Yeast Res. 2005, 5: 545-558. 10.1016/j.femsyr.2004.09.008.
Duarte NC, Herrgard MJ, Palsson BO: Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model. Genome Res. 2004, 14: 1298-1309. 10.1101/gr.2250904.
Papp B, Pal C, Hurst LD: Metabolic network analysis of the causes and evolution of enzyme dispensability in yeast. Nature. 2004, 429: 661-664. 10.1038/nature02636.
Wach A, Brachat A, Pohlmann R, Philippsen P: New heterologous modules for classical or PCR-based gene disruptions in Saccharomyces cerevisiae. Yeast. 1994, 10: 1793-1808.
van Dijken JP, Bauer J, Brambilla L, Duboc P, Francois JM, Gancedo C, Giuseppin ML, Heijnen JJ, Hoare M, Lange HC, et al: An interlaboratory comparison of physiological and genetic properties of four Saccharomyces cerevisiae strains. Enzyme Microb Technol. 2000, 26: 706-714. 10.1016/S0141-0229(00)00162-9.
Pronk JT: Auxotrophic yeast strains in fundamental and applied research. Appl Environ Microbiol. 2002, 68: 2095-2100. 10.1128/AEM.68.5.2095-2100.2002.
Steinmetz LM, Scharfe C, Deutschbauer AM, Mokranjac D, Herman ZS, Jones T, Chu AM, Giaever G, Prokisch H, Oefner PJ, et al: Systematic screen for human disease genes in yeast. Nat Genet. 2002, 31: 400-404.
Zamboni N, Fischer E, Laudert D, Aymerich S, Hohmann HP, Sauer U: The Bacillus subtilis yqjI gene encodes the NADP+-dependent 6-P-gluconate dehydrogenase in the pentose phosphate pathway. J Bacteriol. 2004, 186: 4528-4534. 10.1128/JB.186.14.4528-4534.2004.
Duetz WA, Ruedi L, Hermann R, O'Connor K, Buchs J, Witholt B: Methods for intense aeration, growth, storage, and replication of bacterial strains in microtiter plates. Appl Environ Microbiol. 2000, 66: 2641-2646. 10.1128/AEM.66.6.2641-2646.2000.
Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, et al: Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002, 418: 387-391. 10.1038/nature00935.
Winzeler EA, Liang H, Shoemaker DD, Davis RW: Functional analysis of the yeast genome by precise deletion and parallel phenotypic characterization. Novartis Found Symp. 2000, 229: 105-109. discussion 109-111.
Monschau N, Stahmann KP, Sahm H, McNeil JB, Bognar AL: Identification of Saccharomyces cerevisiae GLY1 as a threonine aldolase: a key enzyme in glycine biosynthesis. FEMS Microbiol Lett. 1997, 150: 55-60. 10.1016/S0378-1097(97)00096-7.
Gombert AK, Moreira dos 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-1451. 10.1128/JB.183.4.1441-1451.2001.
Maaheimo H, Fiaux J, Cakar ZP, Bailey JE, Sauer U, Szyperski T: Central carbon metabolism of Saccharomyces cerevisiae explored by biosynthetic fractional 13C labeling of common amino acids. Eur J Biochem. 2001, 268: 2464-2479. 10.1046/j.1432-1327.2001.02126.x.
Szyperski T: Biosynthetically directed fractional 13C-labeling of proteinogenic amino acids. An efficient analytical tool to investigate intermediary metabolism. Eur J Biochem. 1995, 232: 433-448.
Blank LM, Sauer U: TCA cycle activity in Saccharomyces cerevisiae is a function of the environmentally determined specific growth and glucose uptake rates. Microbiology. 2004, 150: 1085-1093. 10.1099/mic.0.26845-0.
Sauer U, Hatzimanikatis V, Bailey JE, Hochuli M, Szyperski T, Wüthrich K: Metabolic fluxes in riboflavin-producing Bacillus subtilis. Nat Biotechnol. 1997, 15: 448-452. 10.1038/nbt0597-448.
Dos Santos MM, Gombert AK, Christensen B, Olsson L, Nielsen J: Identification of in vivo enzyme activities in the cometabolism of glucose and acetate by Saccharomyces cerevisiae by using (13)C-labeled substrates. Eukaryotic Cell. 2003, 2: 599-608. 10.1128/EC.2.3.599-608.2003.
Meaden PG, Dickinson FM, Mifsud A, Tessier W, Westwater J, Bussey H, Midgley M: The ALD6 gene of Saccharomyces cerevisiae encodes a cytosolic, Mg(2+)-activated acetaldehyde dehydrogenase. Yeast. 1997, 13: 1319-1327. 10.1002/(SICI)1097-0061(199711)13:14<1319::AID-YEA183>3.0.CO;2-T.
Grabowska D, Chelstowska A: The ALD6 gene product is indispensable for providing NADPH in yeast cells lacking glucose-6-phosphate dehydrogenase activity. J Biol Chem. 2003, 278: 13984-13988. 10.1074/jbc.M210076200.
Boubekeur S, Bunoust O, Camougrand N, Castroviejo M, Rigoulet M, Guerin B: A mitochondrial pyruvate dehydrogenase bypass in the yeast Saccharomyces cerevisiae. J Biol Chem. 1999, 274: 21044-21048. 10.1074/jbc.274.30.21044.
Outten CE, Culotta VC: A novel NADH kinase is the mitochondrial source of NADPH in Saccharomyces cerevisiae. EMBO J. 2003, 22: 2015-2024. 10.1093/emboj/cdg211.
McEwen JE, Cumsky MG, Ko C, Power SD, Poyton RO: Mitochondrial membrane biogenesis: characterization and use of pet mutants to clone the nuclear gene coding for subunit V of yeast cytochrome c oxidase. J Cell Biochem. 1984, 24: 229-242. 10.1002/jcb.240240305.
McAlister-Henn L, Thompson LM: Isolation and expression of the gene encoding yeast mitochondrial malate dehydrogenase. J Bacteriol. 1987, 169: 5157-5166.
Schuster S, Fell DA, Dandekar T: A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks. Nat Biotechnol. 2000, 18: 326-332. 10.1038/73786.
Stelling J, Klamt S, Bettenbrock K, Schuster S, Gilles ED: Metabolic network structure determines key aspects of functionality and regulation. Nature. 2002, 420: 190-193. 10.1038/nature01166.
Segre D, Vitkup D, Church GM: Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA. 2002, 99: 15112-15117. 10.1073/pnas.232349399.
Covert MW, Knight EM, Reed JL, Herrgard MJ, Palsson BO: Integrating high-throughput and computational data elucidates bacterial networks. Nature. 2004, 429: 92-96. 10.1038/nature02456.
Segre D: The regulatory software of cellular metabolism. Trends Biotechnol. 2004, 22: 261-265. 10.1016/j.tibtech.2004.04.013.
Dunn B, Ferea T, Spellman P, Schwarz J, Terraciano J, Troyanovich J, Walker S, Greene J, Shaw K, DiDomenico B, et al: Genetic footprinting: a functional analysis of the S. cerevisiae genome. [http://www.yeastgenome.org]
Wu M, Tzagoloff A: Mitochondrial and cytoplasmic fumarases in Saccharomyces cerevisiae are encoded by a single nuclear gene FUM1. J Biol Chem. 1987, 262: 12275-12282.
Van Maris AJ, Luttik MA, Winkler AA, Van Dijken JP, Pronk JT: Overproduction of threonine aldolase circumvents the biosynthetic role of pyruvate decarboxylase in glucose-limited chemostat cultures of Saccharomyces cerevisiae. Appl Environ Microbiol. 2003, 69: 2094-2099. 10.1128/AEM.69.4.2094-2099.2003.
Winzeler EA, Shoemaker DD, Astromoff A, Liang H, Anderson K, Andre B, Bangham R, Benito R, Boeke JD, Bussey H, et al: Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis. Science. 1999, 285: 901-906. 10.1126/science.285.5429.901.
Nowak MA, Boerlijst MC, Cooke J, Smith JM: Evolution of genetic redundancy. Nature. 1997, 388: 167-171. 10.1038/40618.
Wagner A: Robustness against mutations in genetic networks of yeast. Nat Genet. 2000, 24: 355-361. 10.1038/74174.
Gu Z, Steinmetz LM, Gu X, Scharfe C, Davis RW, Li WH: Role of duplicate genes in genetic robustness against null mutations. Nature. 2003, 421: 63-66. 10.1038/nature01198.
Stelling J, Sauer U, Szallasi Z, Doyle FJ, Doyle J: Robustness of cellular functions. Cell. 2004, 118: 675-685. 10.1016/j.cell.2004.09.008.
Baganz F, Hayes A, Marren D, Gardner DC, Oliver SG: Suitability of replacement markers for functional analysis studies in Saccharomyces cerevisiae. Yeast. 1997, 13: 1563-1573. 10.1002/(SICI)1097-0061(199712)13:16<1563::AID-YEA240>3.3.CO;2-Y.
Allen J, Davey HM, Broadhurst D, Heald JK, Rowland JJ, Oliver SG, Kell DB: High-throughput classification of yeast mutants for functional genomics using metabolic footprinting. Nat Biotechnol. 2003, 21: 692-696. 10.1038/nbt823.
Raamsdonk LM, Teusink B, Broadhurst D, Zhang N, Hayes A, Walsh MC, Berden JA, Brindle KM, Kell DB, Rowland JJ, et al: A functional genomics strategy that uses metabolome data to reveal the phenotype of silent mutations. Nat Biotechnol. 2001, 19: 45-50. 10.1038/83496.
Verduyn C, Postma E, Scheffers WA, Van Dijken JP: Effect of benzoic acid on metabolic fluxes in yeasts: a continuous-culture study on the regulation of respiration and alcoholic fermentation. Yeast. 1992, 8: 501-517.
Sauer U, Lasko DR, Fiaux J, Hochuli M, Glaser R, Szyperski T, Wüthrich K, Bailey JE: Metabolic flux ratio analysis of genetic and environmental modulations of Escherichia coli central carbon metabolism. J Bacteriol. 1999, 181: 6679-6688.
Christensen B, Christiansen T, Gombert AK, Thykaer J, Nielsen J: Simple and robust method for estimation of the split between the oxidative pentose phosphate pathway and the Embden-Meyerhof-Parnas pathway in microorganisms. Biotechnol Bioeng. 2001, 74: 517-523. 10.1002/bit.1143.
Nissen TL, Schulze U, Nielsen J, Villadsen J: Flux distributions in anaerobic, glucose-limited continuous cultures of Saccharomyces cerevisiae. Microbiology. 1997, 143: 203-218.
Fischer E, Zamboni N, Sauer U: High-throughput metabolic flux analysis based on gas chromatography-mass spectrometry derived 13C constraints. Anal Biochem. 2004, 325: 308-316. 10.1016/j.ab.2003.10.036.
Edwards JS, Palsson BO: The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities. Proc Natl Acad Sci USA. 2000, 97: 5528-5533. 10.1073/pnas.97.10.5528.
Makhorin A: GNU Linear Programming Kit. 2001, Moscow, Russia: Moscow Aviation Institute, [http://www.gnu.org/software/glpk/glpk.html]
Yeast Deletion Project and Proteomics of Mitochondria Database. [http://www-deletion.stanford.edu/YDPM]
We are grateful to Eckhard Boles for providing the mae1 mutant. LarsM. Blank gratefully acknowledges financial support by the Deutsche Akademie der Naturforscher Leopoldina (BMBF-LPD/8-78).
Electronic supplementary material
Additional File 3: Physiological data, flux ratios and the calculated flux distributions. Physiological data, flux ratios and the calculated flux distributions (XLS 323 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Blank, L.M., Kuepfer, L. & Sauer, U. Large-scale 13C-flux analysis reveals mechanistic principles of metabolic network robustness to null mutations in yeast. Genome Biol 6, R49 (2005). https://doi.org/10.1186/gb-2005-6-6-r49
- Additional Data File
- Pentose Phosphate Pathway
- Flux Analysis
- Glucose Uptake Rate
- Oxidative Pentose Phosphate Pathway