Reconstruction and analysis of a Kluyveromyces marxianus genome-scale metabolic model

BackgroundKluyveromyces marxianus is a thermotolerant yeast with multiple biotechnological potentials for industrial applications, which can metabolize a broad range of carbon sources, including less conventional sugars like lactose, xylose, arabinose and inulin. These phenotypic traits are sustained even up to 45 °C, what makes it a relevant candidate for industrial biotechnology applications, such as ethanol production. It is therefore of much interest to get more insight into the metabolism of this yeast. Recent studies suggested, that thermotolerance is achieved by reducing the number of growth-determining proteins or suppressing oxidative phosphorylation. Here we aimed to find related factors contributing to the thermotolerance of K. marxianus.ResultsHere, we reported the first genome-scale metabolic model of Kluyveromyces marxianus, iSM996, using a publicly available Kluyveromyces lactis model as template. The model was manually curated and refined to include the missing species-specific metabolic capabilities. The iSM996 model includes 1913 reactions, associated with 996 genes and 1531 metabolites. It performed well to predict the carbon source utilization and growth rates under different growth conditions. Moreover, the model was coupled with transcriptomics data and used to perform simulations at various growth temperatures.ConclusionsK. marxianus iSM996 represents a well-annotated metabolic model of thermotolerant yeast, which provides a new insight into theoretical metabolic profiles at different temperatures of K. marxianus. This could accelerate the integrative analysis of multi-omics data, leading to model-driven strain design and improvement.

[1]  Paul Horton,et al.  Nucleic Acids Research Advance Access published May 21, 2007 WoLF PSORT: protein localization predictor , 2007 .

[2]  Anne Richelle,et al.  Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0 , 2019, Nature Protocols.

[3]  Amos Bairoch,et al.  The ENZYME database in 2000 , 2000, Nucleic Acids Res..

[4]  Christoph Wittmann,et al.  Physiology of the yeast Kluyveromyces marxianus during batch and chemostat cultures with glucose as the sole carbon source. , 2007, FEMS yeast research.

[5]  Ian T. Paulsen,et al.  TransportDB 2.0: a database for exploring membrane transporters in sequenced genomes from all domains of life , 2016, Nucleic Acids Res..

[6]  C. Wittmann,et al.  The yeast Kluyveromyces marxianus and its biotechnological potential , 2008, Applied Microbiology and Biotechnology.

[7]  S. Placzek,et al.  The BRENDA enzyme information system-From a database to an expert system. , 2017, Journal of biotechnology.

[8]  Suprayogi,et al.  Genetic basis of the highly efficient yeast Kluyveromyces marxianus: complete genome sequence and transcriptome analyses , 2015, Biotechnology for Biofuels.

[9]  J. Morrissey,et al.  Kluyveromyces marxianus: A yeast emerging from its sister's shadow , 2010 .

[10]  M. Uhlén,et al.  Genome-scale metabolic modelling of hepatocytes reveals serine deficiency in patients with non-alcoholic fatty liver disease , 2014, Nature Communications.

[11]  Helen Treichel,et al.  Inulinase production by Kluyveromyces marxianus NRRL Y-7571 using solid state fermentation , 2006, Applied biochemistry and biotechnology.

[12]  Sunil Bansal,et al.  Production of β-galactosidase by Kluyveromyces marxianus MTCC 1388 using whey and effect of four different methods of enzyme extraction on β-galactosidase activity , 2008, Indian Journal of Microbiology.

[13]  C. Maranas,et al.  Recent advances in the reconstruction of metabolic models and integration of omics data. , 2014, Current opinion in biotechnology.

[14]  Mamoru Yamada,et al.  Genomic analyses of thermotolerant microorganisms used for high-temperature fermentations , 2016, Bioscience, biotechnology, and biochemistry.

[15]  Wenqing Shui,et al.  Understanding the Mechanism of Thermotolerance Distinct From Heat Shock Response Through Proteomic Analysis of Industrial Strains of Saccharomyces cerevisiae* , 2015, Molecular & Cellular Proteomics.

[16]  Isabel Rocha,et al.  iOD907, the first genome-scale metabolic model for the milk yeast Kluyveromyces lactis. , 2014, Biotechnology journal.

[17]  H. Hoshida,et al.  High-temperature fermentation: how can processes for ethanol production at high temperatures become superior to the traditional process using mesophilic yeast? , 2009, Applied Microbiology and Biotechnology.

[18]  Carlos Gancedo,et al.  Construction and characterization of a Saccharomyces cerevisiae strain able to grow on glucosamine as sole carbon and nitrogen source , 2018, Scientific Reports.

[19]  Rinji Akada,et al.  Identification of auxotrophic mutants of the yeast Kluyveromyces marxianus by non‐homologous end joining‐mediated integrative transformation with genes from Saccharomyces cerevisiae , 2013, Yeast.

[20]  G. H. Fleet,et al.  Composition of the cell walls of several yeast species , 1998, Applied Microbiology and Biotechnology.

[21]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[22]  Markus Krummenacker,et al.  The MetaCyc database of metabolic pathways and enzymes , 2017, Nucleic acids research.

[23]  Steffen Klamt,et al.  Memote: A community driven effort towards a standardized genome-scale metabolic model test suite , 2018, bioRxiv.

[24]  Kenta Nakai,et al.  Genome-wide characterization of transcriptional start sites in humans by integrative transcriptome analysis. , 2011, Genome research.

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

[26]  Peter Ruhdal Jensen,et al.  Oxidative Stress at High Temperatures in Lactococcus lactis Due to an Insufficient Supply of Riboflavin , 2013, Applied and Environmental Microbiology.

[27]  M. Gorwa-Grauslund,et al.  Adaptation of Scheffersomyces stipitis to hardwood spent sulfite liquor by evolutionary engineering , 2015, Biotechnology for Biofuels.

[28]  Hilde van der Togt,et al.  Publisher's Note , 2003, J. Netw. Comput. Appl..

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

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

[31]  W. A. Scheffers,et al.  Effect of benzoic acid on metabolic fluxes in yeasts: A continuous‐culture study on the regulation of respiration and alcoholic fermentation , 1992, Yeast.

[32]  Ayca Cankorur-Cetinkaya,et al.  Genome-Wide Transcriptional Response of Saccharomyces cerevisiae to Stress-Induced Perturbations , 2016, Front. Bioeng. Biotechnol..

[33]  Hao Wang,et al.  RAVEN 2.0: A versatile toolbox for metabolic network reconstruction and a case study on Streptomyces coelicolor , 2018, bioRxiv.

[34]  Marc-André Lachance,et al.  Kluyveromyces van der Walt (1971) , 2011 .

[35]  Susumu Goto,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 2000, Nucleic Acids Res..

[36]  Muhammad Ibrahim Rajoka,et al.  Hyper-production of a thermotolerant β-xylosidase by a deoxy-D-glucose and cycloheximide resistant mutant derivative of Kluyveromyces marxianus PPY 125 , 2005 .

[37]  Dirk Walther,et al.  Dynamic transcriptional and metabolic responses in yeast adapting to temperature stress. , 2010, Omics : a journal of integrative biology.

[38]  R. Marchant,et al.  Partial characterization of β-glucosidase activity produced by Kluyveromyces marxianus IMB3 during growth on cellobiose-containing media at 45°C , 1995, Biotechnology Letters.

[39]  B. Palsson,et al.  A protocol for generating a high-quality genome-scale metabolic reconstruction , 2010 .

[40]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[41]  Laura O'Hara,et al.  Characterization of the Yeast DGK1-encoded CTP-dependent Diacylglycerol Kinase*♦ , 2008, Journal of Biological Chemistry.