Genome wide analysis of protein production load in Trichoderma reesei

BackgroundThe filamentous fungus Trichoderma reesei (teleomorph Hypocrea jecorina) is a widely used industrial host organism for protein production. In industrial cultivations, it can produce over 100 g/l of extracellular protein, mostly constituting of cellulases and hemicellulases. In order to improve protein production of T. reesei the transcriptional regulation of cellulases and secretory pathway factors have been extensively studied. However, the metabolism of T. reesei under protein production conditions has not received much attention.ResultsTo understand the physiology and metabolism of T. reesei under protein production conditions we carried out a well-controlled bioreactor experiment with extensive analysis. We used minimal media to make the data amenable for modelling and three strain pairs to cover different protein production levels. With RNA-sequencing transcriptomics we detected the concentration of the carbon source as the most important determinant of the transcriptome. As the major transcriptional response concomitant to protein production we detected the induction of selected genes that were putatively regulated by xyr1 and were related to protein transport, amino acid metabolism and transcriptional regulation. We found novel metabolic responses such as production of glycerol and a cellotriose-like compound. We then used this cultivation data for flux balance analysis of T. reesei metabolism and demonstrate for the first time the use of genome wide stoichiometric metabolic modelling for T. reesei. We show that our model can predict protein production rate and provides novel insight into the metabolism of protein production. We also provide this unprecedented cultivation and transcriptomics data set for future modelling efforts.ConclusionsThe use of stoichiometric modelling can open a novel path for the improvement of protein production in T. reesei. Based on this we propose sulphur assimilation as a major limiting factor of protein production. As an organism with exceptional protein production capabilities modelling of T. reesei can provide novel insight also to other less productive organisms.

[1]  C. Kubicek,et al.  Triggering of cellulase biosynthesis by cellulose in Trichoderma reesei. Involvement of a constitutive, sophorose-inducible, glucose-inhibited beta-diglucoside permease. , 1993, The Journal of biological chemistry.

[2]  M. Penttilä,et al.  The protein disulphide isomerase gene of the fungus Trichoderma reesei is induced by endoplasmic reticulum stress and regulated by the carbon source , 1999, Molecular and General Genetics MGG.

[3]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[4]  Susumu Goto,et al.  Data, information, knowledge and principle: back to metabolism in KEGG , 2013, Nucleic Acids Res..

[5]  Johannes P. van Dijken,et al.  Redox balances in the metabolism of sugars by yeasts (NAD(H); NADP(H); glucose metabolism; xylose fermentation; ethanol; Crabtree effect; Custers effect) , 1986 .

[6]  F. Henrique-Silva,et al.  Cellulase Induction in Trichoderma reesei by Cellulose Requires Its Own Basal Expression* , 1997, The Journal of Biological Chemistry.

[7]  R. Paules,et al.  An integrated stress response regulates amino acid metabolism and resistance to oxidative stress. , 2003, Molecular cell.

[8]  A. Lappalainen,et al.  Cellobiohydrolase from Trichoderma reesei. , 1983, The Biochemical journal.

[9]  G. F. Persinoti,et al.  Deciphering the Cis-Regulatory Elements for XYR1 and CRE1 Regulators in Trichoderma reesei , 2014, PloS one.

[10]  Hideaki Sugawara,et al.  The Sequence Read Archive , 2010, Nucleic Acids Res..

[11]  M. Penttilä,et al.  Isolation of the ace1 Gene Encoding a Cys2-His2 Transcription Factor Involved in Regulation of Activity of the Cellulase Promoter cbh1of Trichoderma reesei * , 2000, The Journal of Biological Chemistry.

[12]  Shuifang Zhu,et al.  Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads , 2014, BMC Bioinformatics.

[13]  M. Penttilä,et al.  Regulation of cellulase gene expression in the filamentous fungus Trichoderma reesei , 1997, Applied and environmental microbiology.

[14]  C. Cramer,et al.  Basic amino acids and inorganic polyphosphates in Neurospora crassa: independent regulation of vacuolar pools , 1980, Journal of bacteriology.

[15]  A. Amore,et al.  Send Orders of Reprints at Reprints@benthamscience.net Regulation of Cellulase and Hemicellulase Gene Expression in Fungi , 2022 .

[16]  Mikelina Gritzali,et al.  The Cellulase System ofTrichoderma: Relationships Between Purified Extracellular Enzymes from Induced or Cellulose-Grown Cells , 1979 .

[17]  Jesús Picó,et al.  Estimation of recombinant protein production in Pichia pastoris based on a constraint-based model , 2012 .

[18]  C. Kubicek,et al.  Systems Analysis of Lactose Metabolism in Trichoderma reesei Identifies a Lactose Permease That Is Essential for Cellulase Induction , 2013, PloS one.

[19]  V. Sahai,et al.  Comparison of growth and maintenance parameters for cellulase biosynthesis by Trichoderma reesei-C5 with some published data , 1994 .

[20]  B. Palsson,et al.  Metabolic Flux Balancing: Basic Concepts, Scientific and Practical Use , 1994, Bio/Technology.

[21]  D. Ussery,et al.  Comparison of protein coding gene contents of the fungal phyla Pezizomycotina and Saccharomycotina , 2007, BMC Genomics.

[22]  Arthur L. Koch,et al.  The Monod Model and Its Alternatives , 1998 .

[23]  Yang Xu,et al.  R/BHC: fast Bayesian hierarchical clustering for microarray data , 2009, BMC Bioinformatics.

[24]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[25]  Y. Anraku,et al.  Substrate specificities of active transport systems for amino acids in vacuolar-membrane vesicles of Saccharomyces cerevisiae. Evidence of seven independent proton/amino acid antiport systems. , 1984, The Journal of biological chemistry.

[26]  S. Baker,et al.  The putative protein methyltransferase LAE1 controls cellulase gene expression in Trichoderma reesei , 2012, Molecular microbiology.

[27]  J. Visser,et al.  A transcriptional activator, AoXlnR, controls the expression of genes encoding xylanolytic enzymes in Aspergillus oryzae. , 2002, Fungal genetics and biology : FG & B.

[28]  Merja Penttilä,et al.  Transcriptional regulation of plant cell wall degradation by filamentous fungi. , 2005, FEMS microbiology reviews.

[29]  L. Poughon,et al.  A new stoichiometric miniaturization strategy for screening of industrial microbial strains: application to cellulase hyper-producing Trichoderma reesei strains , 2012, Microbial Cell Factories.

[30]  J. Schwender,et al.  Transcript abundance on its own cannot be used to infer fluxes in central metabolism , 2014, Front. Plant Sci..

[31]  Masahiro Samejima,et al.  Cellotriose and Cellotetraose as Inducers of the Genes Encoding Cellobiohydrolases in the Basidiomycete Phanerochaete chrysosporium , 2010, Applied and Environmental Microbiology.

[32]  Tetsuo Kobayashi,et al.  Identification of specific binding sites for XYR1, a transcriptional activator of cellulolytic and xylanolytic genes in Trichoderma reesei. , 2009, Fungal genetics and biology : FG & B.

[33]  H. Blanch,et al.  Lactase production in continuous culture by Trichoderma reesei Rut-C30 , 1984, Biotechnology Letters.

[34]  G. Jan,et al.  Comparative secretome analyses of two Trichoderma reesei RUT-C30 and CL847 hypersecretory strains , 2008, Biotechnology for biofuels.

[35]  Christoph Wittmann,et al.  Integration of in vivo and in silico metabolic fluxes for improvement of recombinant protein production. , 2012, Metabolic engineering.

[36]  U. Sauer,et al.  Regulation and control of metabolic fluxes in microbes. , 2011, Current opinion in biotechnology.

[37]  M. Penttilä,et al.  Activation mechanisms of the HACI‐mediated unfolded protein response in filamentous fungi , 2003, Molecular microbiology.

[38]  J. Cherry,et al.  Directed evolution of industrial enzymes: an update. , 2003, Current opinion in biotechnology.

[39]  Juho Rousu,et al.  Comparative Genome-Scale Reconstruction of Gapless Metabolic Networks for Present and Ancestral Species , 2014, PLoS Comput. Biol..

[40]  R. Mahadevan,et al.  The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. , 2003, Metabolic engineering.

[41]  Monika Schmoll,et al.  Sulphur metabolism and cellulase gene expression are connected processes in the filamentous fungus Hypocrea jecorina (anamorph Trichoderma reesei) , 2008, BMC Microbiology.

[42]  M. Penttilä,et al.  The glucose repressor genecre1 ofTrichoderma: Isolation and expression of a full-length and a truncated mutant form , 1996, Molecular and General Genetics MGG.

[43]  Oskari Vinko.  Inferring Trichoderma reesei Gene Regulatory Network , 2013 .

[44]  Mikko Arvas,et al.  Screening of candidate regulators for cellulase and hemicellulase production in Trichoderma reesei and identification of a factor essential for cellulase production , 2014, Biotechnology for Biofuels.

[45]  C. Kubicek,et al.  Induction of cellulase formation in Trichoderma reesei by cellobiono-1,5-lacton , 1989, Archives of Microbiology.

[46]  Jens Nielsen,et al.  Imbalance of heterologous protein folding and disulfide bond formation rates yields runaway oxidative stress , 2012, BMC Biology.

[47]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[48]  Jens Nielsen,et al.  Quantitative Analysis of Glycerol Accumulation, Glycolysis and Growth under Hyper Osmotic Stress , 2013, PLoS Comput. Biol..

[49]  J. Visser,et al.  Isolation and analysis of xlnR, encoding a transcriptional activator co‐ordinating xylanolytic expression in Aspergillus niger , 1998, Molecular microbiology.

[50]  B. Palsson,et al.  Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. , 2003, Genome research.

[51]  N. Price,et al.  Genome-scale modeling for metabolic engineering , 2015, Journal of Industrial Microbiology & Biotechnology.

[52]  Matthias Heinemann,et al.  Functioning of a metabolic flux sensor in Escherichia coli , 2012, Proceedings of the National Academy of Sciences.

[53]  Sang Yup Lee,et al.  Model based engineering of Pichia pastoris central metabolism enhances recombinant protein production , 2014, Metabolic engineering.

[54]  M. Penttilä,et al.  Protein production and induction of the unfolded protein response in Trichoderma reesei strain Rut-C30 and its transformant expressing endoglucanase I with a hydrophobic tag. , 2005, Biotechnology and bioengineering.

[55]  Yongzhao Zhan,et al.  Maximum Neighborhood Margin Discriminant Projection for Classification , 2014, TheScientificWorldJournal.

[56]  Paula Jouhten,et al.  Metabolic flux profiling of recombinant protein secreting Pichia pastoris growing on glucose:methanol mixtures , 2012, Microbial Cell Factories.

[57]  S. Baker,et al.  The Genomes of Three Uneven Siblings: Footprints of the Lifestyles of Three Trichoderma Species , 2016, Microbiology and Molecular Reviews.

[58]  M. Penttilä,et al.  The Effects of Drugs Inhibiting Protein Secretion in the Filamentous Fungus Trichoderma reesei , 2003, Journal of Biological Chemistry.

[59]  N. Aro,et al.  The effects of extracellular pH and of the transcriptional regulator PACI on the transcriptome of Trichoderma reesei , 2015, Microbial Cell Factories.

[60]  Markus Heinonen,et al.  Detecting time periods of differential gene expression using Gaussian processes: an application to endothelial cells exposed to radiotherapy dose fraction , 2015, Bioinform..

[61]  Robert Gentleman,et al.  Using GOstats to test gene lists for GO term association , 2007, Bioinform..

[62]  M. Bailey,et al.  Efficient cellulase production by Trichoderma reesei in continuous cultivation on lactose medium with a computer-controlled feeding strategy , 2003, Applied Microbiology and Biotechnology.

[63]  John D. Storey,et al.  Strong control, conservative point estimation and simultaneous conservative consistency of false discovery rates: a unified approach , 2004 .

[64]  M. Penttilä,et al.  Transcriptional Regulation of xyn1, Encoding Xylanase I, in Hypocrea jecorina , 2006, Eukaryotic Cell.

[65]  S. Hohmann.  Osmotic Stress Signaling and Osmoadaptation in Yeasts , 2002, Microbiology and Molecular Biology Reviews.

[66]  S. Henry,et al.  Revising the Representation of Fatty Acid, Glycerolipid, and Glycerophospholipid Metabolism in the Consensus Model of Yeast Metabolism. , 2013, Industrial biotechnology.

[67]  George A. Milliken,et al.  Mathematical Modeling in Microbial Ecology , 1997, Chapman & Hall Microbiology Series.

[68]  Bernard Henrissat,et al.  Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei (syn. Hypocrea jecorina) , 2008, Nature Biotechnology.

[69]  Matthias G. Steiger,et al.  Transformation System for Hypocrea jecorina (Trichoderma reesei) That Favors Homologous Integration and Employs Reusable Bidirectionally Selectable Markers , 2010, Applied and Environmental Microbiology.

[70]  Merja Penttilä,et al.  The effect of specific growth rate on protein synthesis and secretion in the filamentous fungus Trichoderma reesei. , 2005, Microbiology.

[71]  Mikko Arvas,et al.  Common features and interesting differences in transcriptional responses to secretion stress in the fungi Trichoderma reesei and Saccharomyces cerevisiae , 2006, BMC Genomics.

[72]  N. Slonim,et al.  A universal framework for regulatory element discovery across all genomes and data types. , 2007, Molecular cell.

[73]  Bernard Henrissat,et al.  Corrigendum: Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei (syn. Hypocrea jecorina) , 2008, Nature Biotechnology.

[74]  R. Toledo,et al.  Cellulase production in continuous culture by Trichoderma reesei on xylose‐based media , 1992, Biotechnology and bioengineering.

[75]  Carl E. Rasmussen,et al.  Derivative Observations in Gaussian Process Models of Dynamic Systems , 2002, NIPS.

[76]  Dan M. Bolser,et al.  Ensembl Genomes 2016: more genomes, more complexity , 2015, Nucleic Acids Res..

[77]  Kara Dolinski,et al.  Saccharomyces Genome Database (SGD) provides biochemical and structural information for budding yeast proteins , 2003, Nucleic Acids Res..

[78]  E. de Nadal,et al.  Regulation of the Sko1 transcriptional repressor by the Hog1 MAP kinase in response to osmotic stress , 2001, The EMBO journal.

[79]  Matthias Heinemann,et al.  A flux-sensing mechanism could regulate the switch between respiration and fermentation. , 2012, FEMS yeast research.

[80]  Li-Jun Ma,et al.  Systematic discovery of regulatory motifs in Fusarium graminearum by comparing four Fusarium genomes , 2010, BMC Genomics.

[81]  Eric R. Ziegel,et al.  The Elements of Statistical Learning , 2003, Technometrics.

[82]  T. Pakula,et al.  The cargo and the transport system: secreted proteins and protein secretion in Trichoderma reesei (Hypocrea jecorina). , 2012, Microbiology.

[83]  S. Kuhara,et al.  The impact of a single-nucleotide mutation of bgl2 on cellulase induction in a Trichoderma reesei mutant , 2015, Biotechnology for Biofuels.

[84]  Mikko Arvas,et al.  Correlation of gene expression and protein production rate - a system wide study , 2011, BMC Genomics.

[85]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[86]  J. Buchert,et al.  Novel Coprinopsis cinerea Polyesterase That Hydrolyzes Cutin and Suberin , 2009, Applied and Environmental Microbiology.

[87]  Ben Taskar,et al.  Rich probabilistic models for gene expression , 2001, ISMB.

[88]  D. Botstein,et al.  Genomic expression programs in the response of yeast cells to environmental changes. , 2000, Molecular biology of the cell.

[89]  Intawat Nookaew,et al.  Genome-scale metabolic reconstructions of Pichia stipitis and Pichia pastoris and in silico evaluation of their potentials , 2012, BMC Systems Biology.

[90]  Martin J. Lercher,et al.  sybil – Efficient constraint-based modelling in R , 2013, BMC Systems Biology.

[91]  Jens Nielsen,et al.  A trispecies Aspergillus microarray: Comparative transcriptomics of three Aspergillus species , 2008, Proceedings of the National Academy of Sciences.

[92]  G. F. Persinoti,et al.  Comparative metabolism of cellulose, sophorose and glucose in Trichoderma reesei using high-throughput genomic and proteomic analyses , 2014, Biotechnology for Biofuels.

[93]  Robert Gentleman,et al.  Software for Computing and Annotating Genomic Ranges , 2013, PLoS Comput. Biol..

[94]  J. P. Craig,et al.  Conserved and essential transcription factors for cellulase gene expression in ascomycete fungi , 2012, Proceedings of the National Academy of Sciences.

[95]  S. Sze,et al.  Protein abundance in multiplexed samples (PAMUS) for quantitation of Trichoderma reesei secretome. , 2013, Journal of proteomics.

[96]  Eva Albers,et al.  Metabolic characteristics and importance of the universal methionine salvage pathway recycling methionine from 5′‐methylthioadenosine , 2009, IUBMB life.

[97]  Mikko Arvas,et al.  Re-annotation of the CAZy genes of Trichoderma reesei and transcription in the presence of lignocellulosic substrates , 2012, Microbial Cell Factories.

[98]  E. Heinzle,et al.  Engineering the supply chain for protein production/secretion in yeasts and mammalian cells , 2015, Journal of Industrial Microbiology & Biotechnology.

[99]  M. Penttilä,et al.  ACEII, a Novel Transcriptional Activator Involved in Regulation of Cellulase and Xylanase Genes of Trichoderma reesei * , 2001, The Journal of Biological Chemistry.

[100]  E. Reese,et al.  Enzymic production of cellotriose from cellulose. , 1959, Archives of Biochemistry and Biophysics.

[101]  Michael E. Himmel,et al.  Perspectives and New Directions for the Production of Bioethanol Using Consolidated Bioprocessing of Lignocellulose , 2009 .

[102]  Daniel Machado,et al.  Systematic Evaluation of Methods for Integration of Transcriptomic Data into Constraint-Based Models of Metabolism , 2014, PLoS Comput. Biol..

[103]  C. Blugeon,et al.  Kinetic transcriptome analysis reveals an essentially intact induction system in a cellulase hyper-producer Trichoderma reesei strain , 2014, Biotechnology for Biofuels.

[104]  M. Blythe,et al.  Genome-wide transcriptional response of Trichoderma reesei to lignocellulose using RNA sequencing and comparison with Aspergillus niger , 2013, BMC Genomics.