When transcriptome meets metabolome: fast cellular responses of yeast to sudden relief of glucose limitation

Within the first 5 min after a sudden relief from glucose limitation, Saccharomyces cerevisiae exhibited fast changes of intracellular metabolite levels and a major transcriptional reprogramming. Integration of transcriptome and metabolome data revealed tight relationships between the changes at these two levels. Transcriptome as well as metabolite changes reflected a major investment in two processes: adaptation from fully respiratory to respiro-fermentative metabolism and preparation for growth acceleration. At the metabolite level, a severe drop of the AXP pools directly after glucose addition was not accompanied by any of the other three NXP. To counterbalance this loss, purine biosynthesis and salvage pathways were transcriptionally upregulated in a concerted manner, reflecting a sudden increase of the purine demand. The short-term dynamics of the transcriptome revealed a remarkably fast decrease in the average half-life of downregulated genes. This acceleration of mRNA decay can be interpreted both as an additional nucleotide salvage pathway and an additional level of glucose-induced regulation of gene expression

Full humanization of the glycolytic pathway in Saccharomyces cerevisiae

Although transplantation of single genes in yeast plays a key role in elucidating gene functionality in metazoans, technical challenges hamper humanization of full pathways and processes. Empowered by advances in synthetic biology, this study demonstrates the feasibility and implementation of full humanization of glycolysis in yeast. Single gene and full pathway transplantation revealed the remarkable conservation of glycolytic and moonlighting functions and, combined with evolutionary strategies, brought to light context-dependent responses. Human hexokinase 1 and 2, but not 4, required mutations in their catalytic or allosteric sites for functionality in yeast, whereas hexokinase 3 was unable to complement its yeast ortholog. Comparison with human tissues cultures showed preservation of turnover numbers of human glycolytic enzymes in yeast and human cell cultures. This demonstration of transplantation of an entire essential pathway paves the way for establishment of species-, tissue-, and disease-specific metazoan models

Sampling the Solution Space in Genome-Scale Metabolic Networks Reveals Transcriptional Regulation in Key Enzymes

Genome-scale metabolic models are available for an increasing number of organisms and can be used to define the region of feasible metabolic flux distributions. In this work we use as constraints a small set of experimental metabolic fluxes, which reduces the region of feasible metabolic states. Once the region of feasible flux distributions has been defined, a set of possible flux distributions is obtained by random sampling and the averages and standard deviations for each of the metabolic fluxes in the genome-scale model are calculated. These values allow estimation of the significance of change for each reaction rate between different conditions and comparison of it with the significance of change in gene transcription for the corresponding enzymes. The comparison of flux change and gene expression allows identification of enzymes showing a significant correlation between flux change and expression change (transcriptional regulation) as well as reactions whose flux change is likely to be driven only by changes in the metabolite concentrations (metabolic regulation). The changes due to growth on four different carbon sources and as a consequence of five gene deletions were analyzed for Saccharomyces cerevisiae. The enzymes with transcriptional regulation showed enrichment in certain transcription factors. This has not been previously reported. The information provided by the presented method could guide the discovery of new metabolic engineering strategies or the identification of drug targets for treatment of metabolic diseases

De novo sequencing, assembly and analysis of the genome of the laboratory strain Saccharomyces cerevisiae CEN.PK113-7D, a model for modern industrial biotechnology

Saccharomyces cerevisiae CEN.PK 113-7D is widely used for metabolic engineering and systems biology research in industry and academia. We sequenced, assembled, annotated and analyzed its genome. Single-nucleotide variations (SNV), insertions/deletions (indels) and differences in genome organization compared to the reference strain S. cerevisiae S288C were analyzed. In addition to a few large deletions and duplications, nearly 3000 indels were identified in the CEN.PK113-7D genome relative to S288C. These differences were overrepresented in genes whose functions are related to transcriptional regulation and chromatin remodelling. Some of these variations were caused by unstable tandem repeats, suggesting an innate evolvability of the corresponding genes. Besides a previously characterized mutation in adenylate cyclase, the CEN.PK113-7D genome sequence revealed a significant enrichment of non-synonymous mutations in genes encoding for components of the cAMP signalling pathway. Some phenotypic characteristics of the CEN.PK113-7D strains were explained by the presence of additional specific metabolic genes relative to S288C. In particular, the presence of the BIO1 and BIO6 genes correlated with a biotin prototrophy of CEN.PK113-7D. Furthermore, the copy number, chromosomal location and sequences of the MAL loci were resolved. The assembled sequence reveals that CEN.PK113-7D has a mosaic genome that combines characteristics of laboratory strains and wild-industrial strains

Stress-induced expression is enriched for evolutionarily young genes in diverse budding yeasts

The Saccharomycotina subphylum (budding yeasts) spans 400 million years of evolution and includes species that thrive in diverse environments. To study niche-adaptation, we identify changes in gene expression in three divergent yeasts grown in the presence of various stressors. Duplicated and non-conserved genes are significantly more likely to respond to stress than genes that are conserved as single-copy orthologs. Next, we develop a sorting method that considers evolutionary origin and duplication timing to assign an evolutionary age to each gene. Subsequent analysis reveals that genes that emerged in recent evolutionary time are enriched amongst stress-responsive genes for each species. This gene expression pattern suggests that budding yeasts share a stress adaptation mechanism, whereby selective pressure leads to functionalization of young genes to improve growth in adverse conditions. Further characterization of young genes from species that thrive in harsh environments can inform the design of more robust strains for biotechnology

Transcriptional responses to glucose in Saccharomyces cerevisiae strains lacking a functional protein kinase A

Background The pattern of gene transcripts in the yeast Saccharomyces cerevisiae is strongly affected by the presence of glucose. An increased activity of protein kinase A (PKA), triggered by a rise in the intracellular concentration of cAMP, can account for many of the effects of glucose on transcription. In S. cerevisiae three genes, TPK1, TPK2, and TPK3, encode catalytic subunits of PKA. The lack of viability of tpk1 tpk2 tpk3 triple mutants may be suppressed by mutations such as yak1 or msn2/msn4. To investigate the requirement for PKA in glucose control of gene expression, we have compared the effects of glucose on global transcription in a wild-type strain and in two strains devoid of PKA activity, tpk1 tpk2 tpk3 yak1 and tpk1 tpk2 tpk3 msn2 msn4. Results We have identified different classes of genes that can be induced -or repressed- by glucose in the absence of PKA. Representative examples are genes required for glucose utilization and genes involved in the metabolism of other carbon sources, respectively. Among the genes responding to glucose in strains devoid of PKA some are also controlled by a redundant signalling pathway involving PKA activation, while others are not affected when PKA is activated through an increase in cAMP concentration. On the other hand, among genes that do not respond to glucose in the absence of PKA, some give a full response to increased cAMP levels, even in the absence of glucose, while others appear to require the cooperation of different signalling pathways. We show also that, for a number of genes controlled by glucose through a PKA-dependent pathway, the changes in mRNA levels are transient. We found that, in cells grown in gluconeogenic conditions, expression of a small number of genes, mainly connected with the response to stress, is reduced in the strains lacking PKA. Conclusions In S. cerevisiae, the transcriptional responses to glucose are triggered by a variety of pathways, alone or in combination, in which PKA is often involved. Redundant signalling pathways confer a greater robustness to the response to glucose, while cooperative pathways provide a greater flexibility.BT/BiotechnologyApplied Science

Correlation between transcript profiles and fitness of deletion mutants in anaerobic chemostat cultures of Saccharomyces cerevisiae

The applicability of transcriptomics for functional genome analysis rests on the assumption that global information on gene function can be inferred from transcriptional regulation patterns. This study investigated whether Saccharomyces cerevisiae genes that show a consistently higher transcript level under anaerobic than aerobic conditions do indeed contribute to fitness in the absence of oxygen. Tagged deletion mutants were constructed in 27 S. cerevisiae genes that showed a strong and consistent transcriptional upregulation under anaerobic conditions, irrespective of the nature of the growth-limiting nutrient (glucose, ammonia, sulfate or phosphate). Competitive anaerobic chemostat cultivation showed that only five out of the 27 mutants (eug1Δ, izh2Δ, plb2Δ, ylr413wΔ and yor012wΔ) conferred a significant disadvantage relative to a tagged reference strain. The implications of this study are that: (i) transcriptome analysis has a very limited predictive value for the contribution of individual genes to fitness under specific environmental conditions, and (ii) competitive chemostat cultivation of tagged deletion strains offers an efficient approach to select relevant leads for functional analysis studies

Oxygen dependence of metabolic fluxes and energy generation of Saccharomyces cerevisiae CEN.PK113-1A

<p>Abstract</p> <p>Background</p> <p>The yeast <it>Saccharomyces cerevisiae </it>is able to adjust to external oxygen availability by utilizing both respirative and fermentative metabolic modes. Adjusting the metabolic mode involves alteration of the intracellular metabolic fluxes that are determined by the cell's multilevel regulatory network. Oxygen is a major determinant of the physiology of <it>S. cerevisiae </it>but understanding of the oxygen dependence of intracellular flux distributions is still scarce.</p> <p>Results</p> <p>Metabolic flux distributions of <it>S. cerevisiae </it>CEN.PK113-1A growing in glucose-limited chemostat cultures at a dilution rate of 0.1 h^-1with 20.9%, 2.8%, 1.0%, 0.5% or 0.0% O₂in the inlet gas were quantified by ¹³C-MFA. Metabolic flux ratios from fractional [U-¹³C]glucose labelling experiments were used to solve the underdetermined MFA system of central carbon metabolism of <it>S. cerevisiae</it>.</p> <p>While ethanol production was observed already in 2.8% oxygen, only minor differences in the flux distribution were observed, compared to fully aerobic conditions. However, in 1.0% and 0.5% oxygen the respiratory rate was severely restricted, resulting in progressively reduced fluxes through the TCA cycle and the direction of major fluxes to the fermentative pathway. A redistribution of fluxes was observed in all branching points of central carbon metabolism. Yet only when oxygen provision was reduced to 0.5%, was the biomass yield exceeded by the yields of ethanol and CO₂. Respirative ATP generation provided 59% of the ATP demand in fully aerobic conditions and still a substantial 25% in 0.5% oxygenation. An extensive redistribution of fluxes was observed in anaerobic conditions compared to all the aerobic conditions. Positive correlation between the transcriptional levels of metabolic enzymes and the corresponding fluxes in the different oxygenation conditions was found only in the respirative pathway.</p> <p>Conclusion</p> <p>¹³C-constrained MFA enabled quantitative determination of intracellular fluxes in conditions of different redox challenges without including redox cofactors in metabolite mass balances. A redistribution of fluxes was observed not only for respirative, respiro-fermentative and fermentative metabolisms, but also for cells grown with 2.8%, 1.0% and 0.5% oxygen. Although the cellular metabolism was respiro-fermentative in each of these low oxygen conditions, the actual amount of oxygen available resulted in different contributions through respirative and fermentative pathways.</p

Stoichiometric representation of geneproteinreaction associations leverages constraint-based analysis from reaction to gene-level phenotype prediction

Genome-scale metabolic reconstructions are currently available for hundreds of organisms. Constraint-based modeling enables the analysis of the phenotypic landscape of these organisms, predicting the response to genetic and environmental perturbations. However, since constraint-based models can only describe the metabolic phenotype at the reaction level, understanding the mechanistic link between genotype and phenotype is still hampered by the complexity of gene-protein-reaction associations. We implement a model transformation that enables constraint-based methods to be applied at the gene level by explicitly accounting for the individual fluxes of enzymes (and subunits) encoded by each gene. We show how this can be applied to different kinds of constraint-based analysis: flux distribution prediction, gene essentiality analysis, random flux sampling, elementary mode analysis, transcriptomics data integration, and rational strain design. In each case we demonstrate how this approach can lead to improved phenotype predictions and a deeper understanding of the genotype-to-phenotype link. In particular, we show that a large fraction of reaction-based designs obtained by current strain design methods are not actually feasible, and show how our approach allows using the same methods to obtain feasible gene-based designs. We also show, by extensive comparison with experimental 13C-flux data, how simple reformulations of different simulation methods with gene-wise objective functions result in improved prediction accuracy. The model transformation proposed in this work enables existing constraint-based methods to be used at the gene level without modification. This automatically leverages phenotype analysis from reaction to gene level, improving the biological insight that can be obtained from genome-scale models.DM was supported by the Portuguese Foundationfor Science and Technologythrough a post-doc fellowship (ref: SFRH/BPD/111519/ 2015). This study was supported by the PortugueseFoundationfor Science and Technology (FCT) under the scope of the strategic fundingof UID/BIO/04469/2013 unitand COMPETE2020 (POCI-01-0145-FEDER-006684) and BioTecNorte operation (NORTE-01-0145FEDER-000004) fundedby EuropeanRegional Development Fund under the scope of Norte2020Programa Operacional Regional do Norte. This project has received fundingfrom the European Union’s Horizon 2020 research and innovation programme under grant agreementNo 686070. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript