Tobias Österlund

Fifteen years of large scale metabolic modeling of yeast: Developments and impacts

Since the first large-scale reconstruction of the Saccharomyces cerevisiae metabolic network 15 years ago the development of yeast metabolic models has progressed rapidly, resulting in no less than nine different yeast genome-scale metabolic models. Here we review the historical development of large-scale mathematical modeling of yeast metabolism and the growing scope and impact of applications of these models in four different areas: as guide for metabolic engineering and strain improvement, as a tool for biological interpretation and discovery, applications of novel computational framework and for evolutionary studies.

Mapping condition-dependent regulation of metabolism in yeast through genome-scale modeling.

BackgroundThe 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.ResultsHere 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.ConclusionWe 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.

Genome-Scale Modeling of the Protein Secretory Machinery in Yeast

The protein secretory machinery in Eukarya is involved in post-translational modification (PTMs) and sorting of the secretory and many transmembrane proteins. While the secretory machinery has been well-studied using classic reductionist approaches, a holistic view of its complex nature is lacking. Here, we present the first genome-scale model for the yeast secretory machinery which captures the knowledge generated through more than 50 years of research. The model is based on the concept of a Protein Specific Information Matrix (PSIM: characterized by seven PTMs features). An algorithm was developed which mimics secretory machinery and assigns each secretory protein to a particular secretory class that determines the set of PTMs and transport steps specific to each protein. Protein abundances were integrated with the model in order to gain system level estimation of the metabolic demands associated with the processing of each specific protein as well as a quantitative estimation of the activity of each component of the secretory machinery.

Statistical evaluation of methods for identification of differentially abundant genes in comparative metagenomics.

BackgroundMetagenomics is the study of microbial communities by sequencing of genetic material directly from environmental or clinical samples. The genes present in the metagenomes are quantified by annotating and counting the generated DNA fragments. Identification of differentially abundant genes between metagenomes can provide important information about differences in community structure, diversity and biological function. Metagenomic data is however high-dimensional, contain high levels of biological and technical noise and have typically few biological replicates. The statistical analysis is therefore challenging and many approaches have been suggested to date.ResultsIn this article we perform a comprehensive evaluation of 14 methods for identification of differentially abundant genes between metagenomes. The methods are compared based on the power to detect differentially abundant genes and their ability to correctly estimate the type I error rate and the false discovery rate. We show that sample size, effect size, and gene abundance greatly affect the performance of all methods. Several of the methods also show non-optimal model assumptions and biased false discovery rate estimates, which can result in too large numbers of false positives. We also demonstrate that the performance of several of the methods differs substantially between metagenomic data sequenced by different technologies.ConclusionsTwo methods, primarily designed for the analysis of RNA sequencing data (edgeR and DESeq2) together with a generalized linear model based on an overdispersed Poisson distribution were found to have best overall performance. The results presented in this study may serve as a guide for selecting suitable statistical methods for identification of differentially abundant genes in metagenomes.

A novel method for cross-species gene expression analysis

BackgroundAnalysis of gene expression from different species is a powerful way to identify evolutionarily conserved transcriptional responses. However, due to evolutionary events such as gene duplication, there is no one-to-one correspondence between genes from different species which makes comparison of their expression profiles complex.ResultsIn this paper we describe a new method for cross-species meta-analysis of gene expression. The method takes the homology structure between compared species into account and can therefore compare expression data from genes with any number of orthologs and paralogs. A simulation study shows that the proposed method results in a substantial increase in statistical power compared to previously suggested procedures. As a proof of concept, we analyzed microarray data from heat stress experiments performed in eight species and identified several well-known evolutionarily conserved transcriptional responses. The method was also applied to gene expression profiles from five studies of estrogen exposed fish and both known and potentially novel responses were identified.ConclusionsThe method described in this paper will further increase the potential and reliability of meta-analysis of gene expression profiles from evolutionarily distant species. The method has been implemented in R and is freely available athttp://bioinformatics.math.chalmers.se/Xspecies/.

Anaerobic α-Amylase Production and Secretion with Fumarate as the Final Electron Acceptor in Saccharomyces cerevisiae

ABSTRACT In this study, we focus on production of heterologous α-amylase in the yeast Saccharomyces cerevisiae under anaerobic conditions. We compare the metabolic fluxes and transcriptional regulation under aerobic and anaerobic conditions, with the objective of identifying the final electron acceptor for protein folding under anaerobic conditions. We find that yeast produces more amylase under anaerobic conditions than under aerobic conditions, and we propose a model for electron transfer under anaerobic conditions. According to our model, during protein folding the electrons from the endoplasmic reticulum are transferred to fumarate as the final electron acceptor. This model is supported by findings that the addition of fumarate under anaerobic (but not aerobic) conditions improves cell growth, specifically in the α-amylase-producing strain, in which it is not used as a carbon source. Our results provide a model for the molecular mechanism of anaerobic protein secretion using fumarate as the final electron acceptor, which may allow for further engineering of yeast for improved protein secretion under anaerobic growth conditions.

Management of the endoplasmic reticulum stress by activation of the heat shock response in yeast

In yeast Saccharomyces cerevisiae, accumulation of misfolded proteins in the endoplasmic reticulum (ER) causes ER stress and activates the unfolded protein response (UPR), which is mediated by Hac1p. The heat shock response (HSR) mediated by Hsf1p, mainly regulates cytosolic processes and protects the cell from stresses. Here, we find that a constitutive activation of the HSR could increase ER stress resistance in both wild-type and UPR-deficient cells. Activation of HSR decreased UPR activation in the WT (as shown by the decreased HAC1 mRNA splicing). We analyzed the genome-wide transcriptional response in order to propose regulatory mechanisms that govern the interplay between UPR and HSR and followed up for the hypotheses by experiments in vivo and in vitro. Interestingly, we found that the regulation of ER stress response via HSR is (1) only partially dependent on over-expression of Kar2p (ER resident chaperone induced by ER stress); (2) does not involve the increase in protein turnover via the proteasome activity; (3) is related to the oxidative stress response. From the transcription data, we also propose that HSR enhances ER stress resistance mainly through facilitation of protein folding and secretion. We also find that HSR coordinates multiple stress-response pathways, including the repression of the overall transcription and translation.

Genome-scale analysis of the high-efficient protein secretion system of Aspergillus oryzae

BackgroundThe koji mold, Aspergillus oryzae is widely used for the production of industrial enzymes due to its particularly high protein secretion capacity and ability to perform post-translational modifications. However, systemic analysis of its secretion system is lacking, generally due to the poorly annotated proteome.ResultsHere we defined a functional protein secretory component list of A. oryzae using a previously reported secretory model of S. cerevisiae as scaffold. Additional secretory components were obtained by blast search with the functional components reported in other closely related fungal species such as Aspergillus nidulans and Aspergillus niger. To evaluate the defined component list, we performed transcriptome analysis on three α-amylase over-producing strains with varying levels of secretion capacities. Specifically, secretory components involved in the ER-associated processes (including components involved in the regulation of transport between ER and Golgi) were significantly up-regulated, with many of them never been identified for A. oryzae before. Furthermore, we defined a complete list of the putative A. oryzae secretome and monitored how it was affected by overproducing amylase.ConclusionIn combination with the transcriptome data, the most complete secretory component list and the putative secretome, we improved the systemic understanding of the secretory machinery of A. oryzae in response to high levels of protein secretion. The roles of many newly predicted secretory components were experimentally validated and the enriched component list provides a better platform for driving more mechanistic studies of the protein secretory pathway in this industrially important fungus.

Improved Production of a Heterologous Amylase in Saccharomyces cerevisiae by Inverse Metabolic Engineering

ABSTRACT The increasing demand for industrial enzymes and biopharmaceutical proteins relies on robust production hosts with high protein yield and productivity. Being one of the best-studied model organisms and capable of performing posttranslational modifications, the yeast Saccharomyces cerevisiae is widely used as a cell factory for recombinant protein production. However, many recombinant proteins are produced at only 1% (or less) of the theoretical capacity due to the complexity of the secretory pathway, which has not been fully exploited. In this study, we applied the concept of inverse metabolic engineering to identify novel targets for improving protein secretion. Screening that combined UV-random mutagenesis and selection for growth on starch was performed to find mutant strains producing heterologous amylase 5-fold above the level produced by the reference strain. Genomic mutations that could be associated with higher amylase secretion were identified through whole-genome sequencing. Several single-point mutations, including an S196I point mutation in the VTA1 gene coding for a protein involved in vacuolar sorting, were evaluated by introducing these to the starting strain. By applying this modification alone, the amylase secretion could be improved by 35%. As a complement to the identification of genomic variants, transcriptome analysis was also performed in order to understand on a global level the transcriptional changes associated with the improved amylase production caused by UV mutagenesis.

Controllability analysis of transcriptional regulatory networks reveals circular control patterns among transcription factors

Transcriptional regulation is the most committed type of regulation in living cells where transcription factors (TFs) control the expression of their target genes and TF expression is controlled by other TFs forming complex transcriptional regulatory networks that can be highly interconnected. Here we analyze the topology and organization of nine transcriptional regulatory networks for E. coli, yeast, mouse and human, and we evaluate how the structure of these networks influences two of their key properties, namely controllability and stability. We calculate the controllability for each network as a measure of the organization and interconnectivity of the network. We find that the number of driver nodes nD needed to control the whole network is 64% of the TFs in the E. coli transcriptional regulatory network in contrast to only 17% for the yeast network, 4% for the mouse network and 8% for the human network. The high controllability (low number of drivers needed to control the system) in yeast, mouse and human is due to the presence of internal loops in their regulatory networks where the TFs regulate each other in a circular fashion. We refer to these internal loops as circular control motifs (CCM). The E. coli transcriptional regulatory network, which does not have any CCMs, shows a hierarchical structure of the transcriptional regulatory network in contrast to the eukaryal networks. The presence of CCMs also has influence on the stability of these networks, as the presence of cycles can be associated with potential unstable steady-states where even small changes in binding affinities can cause dramatic rearrangements of the state of the network.