Mapping condition-dependent regulation of metabolism in yeast through genome-scale modeling
© Österlund et al.; licensee BioMed Central Ltd. 2013
Received: 16 November 2012
Accepted: 23 April 2013
Published: 30 April 2013
The genome-scale metabolic model of Saccharomyces cerevisiae , first presented in 2003, was the first genome-scale network reconstruction for a eukaryotic organism. Since then continuous efforts have been made in order to improve and expand the yeast metabolic network.
Here we present iTO977, a comprehensive genome-scale metabolic model that contains more reactions, metabolites and genes than previous models. The model was constructed based on two earlier reconstructions, namely iIN800 and the consensus network, and then improved and expanded using gap-filling methods and by introducing new reactions and pathways based on studies of the literature and databases. The model was shown to perform well both for growth simulations in different media and gene essentiality analysis for single and double knock-outs. Further, the model was used as a scaffold for integrating transcriptomics, and flux data from four different conditions in order to identify transcriptionally controlled reactions, i.e. reactions that change both in flux and transcription between the compared conditions.
We present a new yeast model that represents a comprehensive up-to-date collection of knowledge on yeast metabolism. The model was used for simulating the yeast metabolism under four different growth conditions and experimental data from these four conditions was integrated to the model. The model together with experimental data is a useful tool to identify condition-dependent changes of metabolism between different environmental conditions.
KeywordsSaccharomyces cerevisiae Genome-scale metabolic model Integrated analysis Transcriptionally controlled reactions
The number of metabolic network reconstructions for microorganisms has increased rapidly during the past decade following the genome revolution and the rapid increase of genome sequencing projects. The yeast Saccharomyces cerevisiae was the first eukaryote to have its metabolic network reconstructed and its metabolism is very well-studied. During the last 10 years the yeast metabolic network has been further updated and further reconstructed leading to several genome-scale models[1–9]. Yeast genome-scale modeling have many different applications, which have been reviewed lately. Yeast has for example been used extensively in biotechnology as a workhorse for production of a wide range of chemicals, and metabolic modeling has helped in many cases to guide the strain construction by finding strategies for improved chemical production. Recently the S. cerevisiae metabolic network has also been used as the basis for constructing metabolic models for other yeasts, e.g. Schizosaccharomyces pombe, Yarrowia lipolytica, Pichia pastoris and Pichia stipitis.
A genome-scale metabolic model can be thought of 1) being a biochemically and genetically structured database of the metabolism linking enzymes, metabolites and reactions together, and 2) a predictive model that can estimate the fluxes of intracellular reactions as well as predicting the growth of the cell when constraining the extracellular (measured) fluxes to represent the growth conditions. The first part corresponds to the GENRE and the second part corresponds to the GEM concept. The genome-scale metabolic model can be applied as an in silico hypothesis testing tool to gain insight into the operations of metabolism, leading to identification of targets for superior strain development.
Genome-scale metabolic model reconstructions can also be used as a scaffold for integrated data analysis, where the network structure of the model is the basis for the integration of omics data in order to gain more mechanistic knowledge of the cell’s behavior in different environments or under different conditions. The reporter metabolite algorithm is a useful tool for integration of transcriptomics data into the metabolic network. The algorithm identifies key metabolites whose neighboring genes in the metabolic network are transcriptionally changed between two conditions. However, the regulation of the cellular response to different perturbations such as gene deletion or changes of environmental conditions may occur on many different levels in the cell. By using the genome-scale model to predict changes in metabolic flux between different conditions and correlate these changes with the change in expression of the genes involved, reactions that are directly controlled by transcription can be identified. The random sampling algorithm presented by Bordel et al. is a framework for identification of this kind of trancriptionally controlled reactions using the genome-scale metabolic model.
Here we present an updated genome-scale model (GEM) that is based on the consensus network and the iIN800 model. The goal of the model is to provide an up-to-date and comprehensive description of the yeast metabolism, so the new model contains more genes, reactions and metabolites than previous models. We show that the model performs well in simulating the growth under different environmental conditions. The model was also used as a tool for integrating a large amount of omics data from different conditions. It was used together with the random sampling algorithm to correlate the changes in fluxes between aerobic and anaerobic conditions as well as carbon limited and nitrogen limited conditions with changes in gene expression in order to identify condition-dependent transcriptionally controlled reactions.
Results and discussion
In order to construct a comprehensive genome-scale metabolic model, the consensus network was merged with the iIN800 model to construct a draft starting point reconstruction. The model was then improved by adding several new pathways, new reactions and new metabolites to the model based on database searches and literature evidence. The resulting model includes 977 metabolic genes and is called iTO977. Figure1A shows how the iTO977 model was created. The model is available in SBML and Excel format (Additional files1 and2 and through the BioMet Toolbox, http://www.sysbio.se/biomet ) and a list of newly added reactions and pathways is available in Additional file3.
iTO977 model properties
Summary of the characteristics of different yeast genome-scale metabolic models
iTO977 (Current model)
When comparing the number of metabolites and reactions between different models it is important to keep in mind that the same metabolite can be present several times in the model, but in different compartments. In genome-scale metabolic models the same metabolite is considered as two different species if it is localized to two different compartments. When measuring the levels of intracellular metabolites by metabolomics it is impossible to distinguish metabolites in different compartments, e.g. pyruvate in the cytoplasm and pyruvate in mitochondria will be considered as the same species in a metabolome experiment. Similarly a reaction taking place in two different compartments will also be counted twice in the model. Since the yeast 5 model has many more compartments than iTO977 it also has a higher number of metabolites and reactions when taking compartmentalization into account. However, iTO977 has a higher number of unique metabolites and unique reactions than all previous models, i.e. without taking compartmentalization into account when counting the number of metabolites and reactions.
The number of genes included in the iTO977 model is larger than in previous models. Figure1B shows the comparison of the iTO977 model and the Yeast 5 model in terms of number of open reading frames. A list of the ORFs included in the iTO977 model but not in Yeast 5 is presented in Additional file4. The 25 ORFs that are in yeast 5 but not in iTO977 were all considered to be out of scope of the iTO977 model, since they are mainly transporters over e.g. the Golgi membrane, the vacuole membrane and other organelles that are not included in the iTO977 model. Among the 104 ORFs that are included in iTO977, but not in Yeast 5, many of them belong to two pathways that are newly introduced in the iTO977 model. These two pathways represent the biosynthesis of lipid-linked oligosaccharides and Glycosylphosphaditylinositol, two molecules that are both used for protein modifications. These pathways provide a direct link between the iTO977 model and the recent model of the protein secretion machinery in yeast and is a step towards connecting the metabolic model with a model for protein secretion.
The predictive power of the model was evaluated in two ways: 1) by comparing the simulated growth rate to experimental growth rate under different conditions, and 2) by comparing the predicted viability of in silico single and double gene knock-outs with the experimental observations of single and double mutant growth.
Transcriptionally controlled reactions in glycolysis and ergosterol biosynthesis
Many of the steps in the ergosterol biosynthesis pathway were found to be down-regulated in both flux and expression in carbon limited conditions as compared with nitrogen limited conditions (Figure5B). Many of the reactions were also up-regulated in both flux and expression in anaerobic conditions compared to aerobic under nitrogen limitation. Only the reactions considered to be transcriptionally controlled are presented in Figure5.
Identification of responsible transcription factors
Over-represented transcription factor interactions among significantly changed reactions
Anaerobic - Aerobic (C-limited)
Anaerobic - Aerobic (N-limited)
C-limited - N-limited (Anaerobic)
C-limited - N-limited (Aerobic)
When comparing carbon limited growth with nitrogen limited growth during aerobic conditions we identify the transcription factors Opi1p, Yap7p, Dal80p and Dig1p as important regulators when considering both change in transcription and change in flux (Table 2). At glucose limited conditions yeast goes through full respiration and in glucose excess conditions the cell produces ethanol and the flux through the glycolysis is controlled by carbon catabolite repression. Similarly, at high nitrogen concentrations the nitrogen catabolite pathways are repressed. The transcription factor Dal80p is known to be involved in nitrogen catabolite repression. In this analysis we find that the majority of the genes regulated by Dal80p have a higher expression and higher flux in carbon limited (excess of nitrogen) conditions than in nitrogen limited conditions under aerobic cultivation, which indicates that these pathways are transcriptionally regulated by Dal80p when changing from nitrogen limitation to carbon limitation.
In this paper we present iTO977, a new updated genome-scale metabolic model for yeast. The model is bigger in scope than previous models and can be used both for simulations as well as a scaffold for integrated analysis. The model follows the standards that were introduced in the consensus reconstruction regarding annotation of metabolites and reactions and can therefore be seen as a well-annotated database representing the current knowledge of the yeast metabolism. By using the model to simulate different growth conditions and compare the simulated fluxes from different conditions using the random sampling algorithm we investigate the yeast metabolism in carbon limited, nitrogen limited, aerobic and anaerobic conditions. The change in predicted flux for each reaction in the model is compared with the changes in expression of the involved enzymes and reactions that are directly transcriptionally controlled are identified as well as transcription factors controlling these reactions. The genome-scale metabolic model together with experimental data is a valuable tool for data integration in order to find condition-dependent differences in yeast metabolism.
The first step in constructing the iTO977 model was to merge the iIN800 model and the consensus network to form a draft model. This was done by importing the SBML files for the two starting point reconstructions into Matlab format using the RAVEN toolbox. Additional pathways and biochemical reactions were added from Kyoto Encyclopedia of Genes and Genomes, Saccharomyces Genome Database and literature. A list of these reactions with references is available in Additional file3. Reactions from the consensus network taking place in other compartments than extracellular, cytoplasm, mitochondria or peroxisome was moved to the cytoplasm before merging the two models. Since the iIN800 model doesn’t have peroxisome as a model compartment some of the reactions from iIN800 involved in fatty acid synthesis and degradation (beta-oxidation) was moved from cytoplasm to peroxisome before merging. Water and hydroxyl ions were not balanced in the model, assuming that there are other processes in the cell that uses these molecules as well. The model was corrected for inconsistencies and network gaps, using gap filling methods in RAVEN. The reactions identified as inconsistent by Zomorrodi et al. were corrected to improve the ability to predict knock-out growth phenotypes. CHEBI identifiers and INCHI codes was added to all metabolites following the MIRIAM standard and the reactions in the model have links to EC-numbers or KEGG reactions. The biomass equation was taken from the iIN800 model. The final model was converted from Matlab format both to Excel format and SBML format.
In silico growth simulations
Experimental data were collected from literature for chemostat cultivations performed at various dilution rates under each of the four different conditions: carbon limited aerobic[20, 32–35], carbon limited anaerobic[36, 37], nitrogen limited aerobic[35, 36, 38, 39] and nitrogen limited anaerobic[40, 41] (Additional file5). For simulation purposes the uptake of glucose, NH3 and oxygen were constrained in the model according to the measured values for each experiment. FBA simulations were performed with maximization of biomass as the objective function and the deviation from the experimentally measured growth rate was identified.
Anaerobic conditions were simulated with the iTO977 model by constraining the oxygen uptake rate to zero, but also allowing for uptake of sterols to simulate that these compounds are supplied to the media for anaerobic growth.
Gene essentiality predictions
In silico simulations of single gene knock-out growth were carried out for each of the genes in the model by constraining the flux of all corresponding reactions to zero and then estimating the maximal biomass flux using FBA. A relative fitness value was calculated for each mutant as f = mutant growth rate/wild type growth rate. A mutation was considered as lethal if the relative fitness was smaller than a cutoff value. Three different cutoffs were considered, 0.8, 0.9 and 0.95. The simulations were carried out for growth on glucose, both for a simulated minimal media, allowing only uptake of glucose, ammonium, oxygen, sulfate and phosphate, and for a simulated rich media (YPD media), where also uptake of amino acids and nucleotides was allowed. The viability of the mutants was compared to experimental observations of growth phenotypes for 524 single gene knockouts in minimal media and YPD media[6, 30]. The model was also used to simulate the relative fitness for all pairs of double gene knock-outs, simulating growth on glucose in rich media (YPD) and compared to experimental observations for 234,719 double knock-outs.
Experimental data from chemostat cultivations from four different conditions (carbon limited, nitrogen limited, aerobic and anaerobic) was taken from Jewett et al.. The raw transcriptome data were normalized using the Probe Logarithmic Intensity Error (PLIER) method and the moderated t-statistic was used to identify significantly changed genes between two conditions. The upper and lower bounds of the exchange fluxes and lipid and biomass reactions in the iTO977 model were constrained according to experimentally measured fluxes (Additional file7) for each condition. A set of flux distributions compatible with the measured fluxes was generated for each of the considered conditions as it has been previously described. In order to avoid solutions involving loops we set the default bounds for the metabolic fluxes to Inf and –Inf instead of 1000 and −1000 as it is commonly done. For each pair of condition, we compared the first sample of the reference condition with all the samples of the second condition and computed how many times the reference had a higher (or lower) value for each reaction. The process was repeated for all the samples of the reference condition until all the possible combinations of samples at both conditions have been compared. The fraction of the comparisons in which the reference had a higher (or lower) value than the second condition was used as a probability score for the flux decreasing or increasing. The values of the moderated t-statistic from the transcriptome data were used to calculate a probability score P for each reaction representing the probability that the flux and the transcription is significantly changed in the same direction between two conditions. A hypergeometric enrichment test was performed to identify over-represented transcription factor interactions in the list of reactions from each comparison passing the cutoff P > 0.9 based on known transcription factor-gene interactions.
We acknowledge funding by the Novo Nordisk Foundation, the European Research Council (ERC grant no. 247013), Chalmers Foundation and Knut and Alice Wallenberg Foundation. We also thank the Swedish Bioinformatics Infrastructure for Life Sciences (BILS). The open access charge is funded by Chalmers Library.
- Forster J, Famili I, Fu P, Palsson B, Nielsen J: Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. Genome Res. 2003, 13 (2): 244-253. 10.1101/gr.234503.PubMedPubMed CentralView ArticleGoogle Scholar
- Österlund T, Nookaew I, Nielsen J: Fifteen years of large scale metabolic modeling of yeast: developments and impacts. Biotechnol Adv. 2012, 30 (5): 979-988. 10.1016/j.biotechadv.2011.07.021.PubMedView ArticleGoogle Scholar
- Herrgard M, Swainston N, Dobson P, Dunn W, Arga K, Arvas M, Bluthgen N, Borger S, Costenoble R, Heinemann M: A consensus yeast metabolic network reconstruction obtained from a community approach to systems biology. Nat Biotechnol. 2008, 26 (10): 1155-1160. 10.1038/nbt1492.PubMedPubMed CentralView ArticleGoogle Scholar
- Dobson P, Smallbone K, Jameson D, Simeonidis E, Lanthaler K, Pir P, Lu C, Swainston N, Dunn W, Fisher P: Further developments towards a genome-scale metabolic model of yeast. BMC Syst Biol. 2010, 4 (1): 145-10.1186/1752-0509-4-145.PubMedPubMed CentralView ArticleGoogle Scholar
- Heavner BD, Smallbone K, Barker B, Mendes P, Walker LP: Yeast 5–an expanded reconstruction of the Saccharomyces cerevisiae metabolic network. BMC Syst Biol. 2012, 6 (1): 55-10.1186/1752-0509-6-55.PubMedPubMed CentralView ArticleGoogle Scholar
- Mo ML, Palsson BØ, Herrgård MJ: Connecting extracellular metabolomic measurements to intracellular flux states in yeast. BMC Syst Biol. 2009, 3 (1): 37-10.1186/1752-0509-3-37.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuepfer L, Sauer U, Blank LM: Metabolic functions of duplicate genes in Saccharomyces cerevisiae. Genome Res. 2005, 15 (10): 1421-1430. 10.1101/gr.3992505.PubMedPubMed CentralView ArticleGoogle Scholar
- Nookaew I, Jewett M, Meechai A, Thammarongtham C, Laoteng K, Cheevadhanarak S, Nielsen J, Bhumiratana S: The genome-scale metabolic model iIN800 of Saccharomyces cerevisiae and its validation: a scaffold to query lipid metabolism. BMC Syst Biol. 2008, 2: 71-10.1186/1752-0509-2-71.PubMedPubMed CentralView ArticleGoogle Scholar
- Duarte N, Herrgard M, Palsson B: Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model. Genome Res. 2004, 14 (7): 1298-1309. 10.1101/gr.2250904.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim IK, Roldão A, Siewers V, Nielsen J: A systems‒level approach for metabolic engineering of yeast cell factories. FEMS yeast research. 2012, 12 (2): 228-248. 10.1111/j.1567-1364.2011.00779.x.PubMedView ArticleGoogle Scholar
- Sohn SB, Kim TY, Lee JH, Lee SY: Genome-scale metabolic model of the fission yeast Schizosaccharomyces pombe and the reconciliation of in silico/in vivo mutant growth. BMC Syst Biol. 2012, 6 (1): 49-10.1186/1752-0509-6-49.PubMedPubMed CentralView ArticleGoogle Scholar
- Loira N, Dulermo T, Nicaud JM, Sherman D: A genome-scale metabolic model of the lipid-accumulating yeast Yarrowia lipolytica. BMC Syst Biol. 2012, 6 (1): 35-10.1186/1752-0509-6-35.PubMedPubMed CentralView ArticleGoogle Scholar
- Caspeta L, Shoaie S, Agren R, Nookaew I, Nielsen J: Genome-scale metabolic reconstructions of Pichia stipitis and Pichia pastoris and in silico evaluation of their potentials. BMC Syst Biol. 2012, 6 (1): 24-10.1186/1752-0509-6-24.PubMedPubMed CentralView ArticleGoogle Scholar
- Patil K, Nielsen J: Uncovering transcriptional regulation of metabolism by using metabolic network topology. Proc Natl Acad Sci. 2005, 102 (8): 2685-2689. 10.1073/pnas.0406811102.PubMedPubMed CentralView ArticleGoogle Scholar
- Bordel S, Agren R, Nielsen J: Sampling the solution space in genome-scale metabolic networks reveals transcriptional regulation in key enzymes. PLoS Comput Biol. 2010, 6 (7): e1000859-10.1371/journal.pcbi.1000859.PubMedPubMed CentralView ArticleGoogle Scholar
- Cvijovic M, Olivares-Hernandez R, Agren R, Dahr N, Vongsangnak W, Nookaew I, Patil KR, Nielsen J: BioMet Toolbox: genome-wide analysis of metabolism. Nucleic Acids Res. 2010, 38 (suppl 2): W144-PubMedPubMed CentralView ArticleGoogle Scholar
- Burda P, Aebi M: The dolichol pathway of N-linked glycosylation. Biochim Biophys Acta. 1999, 1426 (2): 239-257. 10.1016/S0304-4165(98)00127-5.PubMedView ArticleGoogle Scholar
- Grimme SJ, Westfall BA, Wiedman JM, Taron CH, Orlean P: The essential Smp3 protein is required for addition of the side-branching fourth mannose during assembly of yeast glycosylphosphatidylinositols. J Biol Chem. 2001, 276 (29): 27731-27739. 10.1074/jbc.M101986200.PubMedView ArticleGoogle Scholar
- Feizi A, Österlund T, Petranovic D, Bordel S, Nielsen J: Genome-Scale Modeling of the Protein Secretory Machinery in Yeast. PLoS computational biology. 2013, : -10.1371/journal.pone.0063284.Google Scholar
- Jewett MC, Workman CT, Nookaew I, Pizarro FA, Agosin E, Hellgren LI, Nielsen J: Mapping Condition Dependent Regulation of Lipid Metabolism in Saccharomyces cerevisiae. Submitted for publication 2013
- Costanzo M, Baryshnikova A, Bellay J, Kim Y, Spear ED, Sevier CS, Ding H, Koh JLY, Toufighi K, Mostafavi S: The genetic landscape of a cell. Science Signalling. 2010, 327 (5964): 425-Google Scholar
- Förster J, Famili I, Palsson BØ, Nielsen J: Large-scale evaluation of in silico gene deletions in Saccharomyces cerevisiae. OMICS A Journal of Integrative Biology. 2003, 7 (2): 193-202. 10.1089/153623103322246584.PubMedView ArticleGoogle Scholar
- Smits HP, Hauf J, Müller S, Hobley TJ, Zimmermann FK, Hahn Hägerdal B, Nielsen J, Olsson L: Simultaneous overexpression of enzymes of the lower part of glycolysis can enhance the fermentative capacity of Saccharomyces cerevisiae. Yeast. 2000, 16 (14): 1325-1334. 10.1002/1097-0061(200010)16:14<1325::AID-YEA627>3.0.CO;2-E.View ArticleGoogle 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 (7004): 99-104. 10.1038/nature02800.PubMedPubMed CentralView ArticleGoogle Scholar
- Ambroziak J, Henry SA: INO2 and INO4 gene products, positive regulators of phospholipid biosynthesis in Saccharomyces cerevisiae, form a complex that binds to the INO1 promoter. J Biol Chem. 1994, 269 (21): 15344-15349.PubMedGoogle Scholar
- Wagner C, Dietz M, Wittmann J, Albrecht A, Schüller HJ: The negative regulator Opi1 of phospholipid biosynthesis in yeast contacts the pleiotropic repressor Sin3 and the transcriptional activator Ino2. Mol Microbiol. 2001, 41 (1): 155-166. 10.1046/j.1365-2958.2001.02495.x.PubMedView ArticleGoogle Scholar
- Hofman-Bang J: Nitrogen catabolite repression in Saccharomyces cerevisiae. Mol Biotechnol. 1999, 12 (1): 35-73. 10.1385/MB:12:1:35.PubMedView ArticleGoogle Scholar
- Agren R, Liu L, Shoaie S, Vongsangnak W, Nookaew I, Nielsen J: The RAVEN toolbox and its use for generating a genome-scale metabolic model for Penicillium chrysogenum. PLoS computational biology. 2013, 9 (3): e1002980-10.1371/journal.pcbi.1002980.PubMedPubMed CentralView ArticleGoogle Scholar
- Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita K, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 2006, 34: D354-D357. 10.1093/nar/gkj102.PubMedPubMed CentralView ArticleGoogle Scholar
- Cherry JM, Hong EL, Amundsen C, Balakrishnan R, Binkley G, Chan ET, Christie KR, Costanzo MC, Dwight SS, Engel SR: Saccharomyces Genome Database: the genomics resource of budding yeast. Nucleic Acids Res. 2012, 40 (D1): D700-D705. 10.1093/nar/gkr1029.PubMedPubMed CentralView ArticleGoogle Scholar
- Zomorrodi A, Maranas C: Improving the iMM904 S. cerevisiae metabolic model using essentiality and synthetic lethality data. BMC Syst Biol. 2010, 4 (1): 178-10.1186/1752-0509-4-178.PubMedPubMed CentralView ArticleGoogle Scholar
- Bakker BM, Bro C, Kötter P, Luttik MAH, Van Dijken JP, Pronk JT: The mitochondrial alcohol dehydrogenase Adh3p is involved in a redox shuttle in Saccharomyces cerevisiae. J Bacteriol. 2000, 182 (17): 4730-4737. 10.1128/JB.182.17.4730-4737.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Overkamp KM, Bakker BM, Kötter P, Van Tuijl A, De Vries S, Van Dijken JP, Pronk JT: In Vivo Analysis of the Mechanisms for Oxidation of Cytosolic NADH by Saccharomyces cerevisiaeMitochondria. J Bacteriol. 2000, 182 (10): 2823-2830. 10.1128/JB.182.10.2823-2830.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Gombert AK, Christensen B, Nielsen J, Moreira dos Santos M: Network identification and flux quantification in the central metabolism of Saccharomyces cerevisiae under different conditions of glucose repression. J Bacteriol. 2001, 183 (4): 1441-1451. 10.1128/JB.183.4.1441-1451.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Vemuri G, Eiteman M, McEwen J, Olsson L, Nielsen J: Increasing NADH oxidation reduces overflow metabolism in Saccharomyces cerevisiae. Proc Natl Acad Sci. 2007, 104 (7): 2402-2407. 10.1073/pnas.0607469104.PubMedPubMed CentralView ArticleGoogle Scholar
- Tai SL, Daran-Lapujade P, Luttik MAH, Walsh MC, Diderich JA, Krijger GC, van Gulik WM, Pronk JT, Daran JM: Control of the glycolytic flux in Saccharomyces cerevisiae grown at low temperature. J Biol Chem. 2007, 282 (14): 10243-10251. 10.1074/jbc.M610845200.PubMedView ArticleGoogle Scholar
- Nissen TL, Schulze U, Nielsen J, Villadsen J: Flux distributions in anaerobic, glucose-limited continuous cultures of Saccharomyces cerevisiae. Microbiology. 1997, 143 (1): 203-218. 10.1099/00221287-143-1-203.PubMedView ArticleGoogle Scholar
- Aon JC, Cortassa S: Involvement of nitrogen metabolism in the triggering of ethanol fermentation in aerobic chemostat cultures of Saccharomyces cerevisiae. Metab Eng. 2001, 3 (3): 250-264. 10.1006/mben.2001.0181.PubMedView ArticleGoogle Scholar
- Usaite R, Patil KR, Grotkjær T, Nielsen J, Regenberg B: Global transcriptional and physiological responses of Saccharomyces cerevisiae to ammonium, L-alanine, or L-glutamine limitation. Appl Environ Microbiol. 2006, 72 (9): 6194-6203. 10.1128/AEM.00548-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Tai SL, Boer VM, Daran-Lapujade P, Walsh MC, de Winde JH, Daran JM, Pronk JT: Two-dimensional Transcriptome Analysis in Chemostat Cultures - Combinatorial Effects of Oxygen Availability and Macronutrient Limitation in Saccharomyces cerevisiae. J Biol Chem. 2005, 280 (1): 437-447.PubMedView ArticleGoogle Scholar
- Lidén G, Persson A, Niklasson C, Gustafsson L: Energetics and product formation by Saccharomyces cerevisiae grown in anaerobic chemostats under nitrogen limitation. Appl Microbiol Biotechnol. 1995, 43 (6): 1034-1038. 10.1007/BF00166921.PubMedView ArticleGoogle Scholar