CRI-SPA: a high-throughput method for systematic genetic editing of yeast libraries.

Biological functions are orchestrated by intricate networks of interacting genetic elements. Predicting the interaction landscape remains a challenge for systems biology and new research tools allowing simple and rapid mapping of sequence to function are desirable. Here, we describe CRI-SPA, a method allowing the transfer of chromosomal genetic features from a CRI-SPA Donor strain to arrayed strains in large libraries of Saccharomyces cerevisiae. CRI-SPA is based on mating, CRISPR-Cas9-induced gene conversion, and Selective Ploidy Ablation. CRI-SPA can be massively parallelized with automation and can be executed within a week. We demonstrate the power of CRI-SPA by transferring four genes that enable betaxanthin production into each strain of the yeast knockout collection (≈4800 strains). Using this setup, we show that CRI-SPA is highly efficient and reproducible, and even allows marker-free transfer of genetic features. Moreover, we validate a set of CRI-SPA hits by showing that their phenotypes correlate strongly with the phenotypes of the corresponding mutant strains recreated by reverse genetic engineering. Hence, our results provide a genome-wide overview of the genetic requirements for betaxanthin production. We envision that the simplicity, speed, and reliability offered by CRI-SPA will make it a versatile tool to forward systems-level understanding of biological processes.

From genotype to phenotype: computational approaches for inferring microbial traits relevant to the food industry.

When selecting microbial strains for the production of fermented foods, various microbial phenotypes need to be taken into account to achieve target product characteristics, such as biosafety, flavor, texture, and health-promoting effects. Through continuous advances in sequencing technologies, microbial whole-genome sequences of increasing quality can now be obtained both cheaper and faster, which increases the relevance of genome-based characterization of microbial phenotypes. Prediction of microbial phenotypes from genome sequences makes it possible to quickly screen large strain collections in silico to identify candidates with desirable traits. Several microbial phenotypes relevant to the production of fermented foods can be predicted using knowledge-based approaches, leveraging our existing understanding of the genetic and molecular mechanisms underlying those phenotypes. In the absence of this knowledge, data-driven approaches can be applied to estimate genotype-phenotype relationships based on large experimental datasets. Here, we review computational methods that implement knowledge- and data-driven approaches for phenotype prediction, as well as methods that combine elements from both approaches. Furthermore, we provide examples of how these methods have been applied in industrial biotechnology, with special focus on the fermented food industry.

Genome-scale modeling drives 70-fold improvement of intracellular heme production in Saccharomyces cerevisiae.

Heme is an oxygen carrier and a cofactor of both industrial enzymes and food additives. The intracellular level of free heme is low, which limits the synthesis of heme proteins. Therefore, increasing heme synthesis allows an increased production of heme proteins. Using the genome-scale metabolic model (GEM) Yeast8 for the yeast Saccharomyces cerevisiae, we identified fluxes potentially important to heme synthesis. With this model, in silico simulations highlighted 84 gene targets for balancing biomass and increasing heme production. Of those identified, 76 genes were individually deleted or overexpressed in experiments. Empirically, 40 genes individually increased heme production (up to threefold). Heme was increased by modifying target genes, which not only included the genes involved in heme biosynthesis, but also those involved in glycolysis, pyruvate, Fe-S clusters, glycine, and succinyl-coenzyme A (CoA) metabolism. Next, we developed an algorithmic method for predicting an optimal combination of these genes by using the enzyme-constrained extension of the Yeast8 model, ecYeast8. The computationally identified combination for enhanced heme production was evaluated using the heme ligand-binding biosensor (Heme-LBB). The positive targets were combined using CRISPR-Cas9 in the yeast strain (IMX581-HEM15-HEM14-HEM3-Δshm1-HEM2-Δhmx1-FET4-Δgcv2-HEM1-Δgcv1-HEM13), which produces 70-fold-higher levels of intracellular heme.

Proteome allocations change linearly with the specific growth rate of Saccharomyces cerevisiae under glucose limitation.

Saccharomyces cerevisiae is a widely used cell factory; therefore, it is important to understand how it organizes key functional parts when cultured under different conditions. Here, we perform a multiomics analysis of S. cerevisiae by culturing the strain with a wide range of specific growth rates using glucose as the sole limiting nutrient. Under these different conditions, we measure the absolute transcriptome, the absolute proteome, the phosphoproteome, and the metabolome. Most functional protein groups show a linear dependence on the specific growth rate. Proteins engaged in translation show a perfect linear increase with the specific growth rate, while glycolysis and chaperone proteins show a linear decrease under respiratory conditions. Glycolytic enzymes and chaperones, however, show decreased phosphorylation with increasing specific growth rates; at the same time, an overall increased flux through these pathways is observed. Further analysis show that even though mRNA levels do not correlate with protein levels for all individual genes, the transcriptome level of functional groups correlates very well with its corresponding proteome. Finally, using enzyme-constrained genome-scale modeling, we find that enzyme usage plays an important role in controlling flux in amino acid biosynthesis.

Emergence of Phenotypically Distinct Subpopulations Is a Factor in Adaptation of Recombinant Saccharomyces cerevisiae under Glucose-Limited Conditions.

Cells cultured in a nutrient-limited environment can undergo adaptation, which confers improved fitness under long-term energy limitation. We have shown previously how a recombinant Saccharomyces cerevisiae strain, producing a heterologous insulin product, under glucose-limited conditions adapts over time at the average population level. Here, we investigated this adaptation at the single-cell level by application of fluorescence-activated cell sorting (FACS) and showed that the following three apparent phenotypes underlie the adaptive response observed at the bulk level: (i) cells that drastically reduced insulin production (23%), (ii) cells with reduced enzymatic capacity in central carbon metabolism (46%), and (iii) cells that exhibited pseudohyphal growth (31%). We speculate that the phenotypic heterogeneity is a result of different mechanisms to increase fitness. Cells with reduced insulin productivity have increased fitness by reducing the burden of the heterologous insulin production, and the populations with reduced enzymatic capacity of the central carbon metabolism and pseudohyphal growth have increased fitness toward the glucose-limited conditions. The results highlight the importance of considering population heterogeneity when studying adaptation and evolution. IMPORTANCE The yeast Saccharomyces cerevisiae is an attractive microbial host for industrial production and is used widely for manufacturing, e.g., pharmaceuticals. Chemostat cultivation mode is an efficient cultivation strategy for industrial production processes as it ensures a constant, well-controlled cultivation environment. Nevertheless, both the production of a heterologous product and the constant cultivation environment in the chemostat impose a selective pressure on the production organism, which may result in adaptation and loss of productivity. The exact mechanisms behind the observed adaptation and loss of performance are often unidentified. We used a recombinant S. cerevisiae strain producing heterologous insulin and investigated the adaptation occurring during chemostat growth at the single-cell level. We showed that three apparent phenotypes underlie the adaptive response observed at the bulk level in the chemostat. These findings highlight the importance of considering population heterogeneity when studying adaptation in industrial bioprocesses.

DebaryOmics: an integrative -omics study to understand the halophilic behaviour of Debaryomyces hansenii.

Debaryomyces hansenii is a non-conventional yeast considered to be a well-suited option for a number of different industrial bioprocesses. It exhibits a set of beneficial traits (halotolerant, oleaginous, xerotolerant, inhibitory compounds resistant) which translates to a number of advantages for industrial fermentation setups when compared to traditional hosts. Although D. hansenii has been highly studied during the last three decades, especially in regards to its salt-tolerant character, the molecular mechanisms underlying this natural tolerance should be further investigated in order to broadly use this yeast in biotechnological processes. In this work, we performed a series of chemostat cultivations in controlled bioreactors where D. hansenii (CBS 767) was grown in the presence of either 1M NaCl or KCl and studied the transcriptomic and (phospho)proteomic profiles. Our results show that sodium and potassium trigger different responses at both expression and regulation of protein activity levels and also complemented previous reports pointing to specific cellular processes as key players in halotolerance, moreover providing novel information about the specific genes involved in each process. The phosphoproteomic analysis, the first of this kind ever reported in D. hansenii, also implicated a novel and yet uncharacterized cation transporter in the response to high sodium concentrations.

Genome-scale metabolic modeling of P. thermoglucosidasius NCIMB 11955 reveals metabolic bottlenecks in anaerobic metabolism.

Parageobacillus thermoglucosidasius represents a thermophilic, facultative anaerobic bacterial chassis, with several desirable traits for metabolic engineering and industrial production. To further optimize strain productivity, a systems level understanding of its metabolism is needed, which can be facilitated by a genome-scale metabolic model. Here, we present p-thermo, the most complete, curated and validated genome-scale model (to date) of Parageobacillus thermoglucosidasius NCIMB 11955. It spans a total of 890 metabolites, 1175 reactions and 917 metabolic genes, forming an extensive knowledge base for P. thermoglucosidasius NCIMB 11955 metabolism. The model accurately predicts aerobic utilization of 22 carbon sources, and the predictive quality of internal fluxes was validated with previously published 13C-fluxomics data. In an application case, p-thermo was used to facilitate more in-depth analysis of reported metabolic engineering efforts, giving additional insight into fermentative metabolism. Finally, p-thermo was used to resolve a previously uncharacterised bottleneck in anaerobic metabolism, by identifying the minimal required supplemented nutrients (thiamin, biotin and iron(III)) needed to sustain anaerobic growth. This highlights the usefulness of p-thermo for guiding the generation of experimental hypotheses and for facilitating data-driven metabolic engineering, expanding the use of P. thermoglucosidasius as a high yield production platform.

Benchmarking accuracy and precision of intensity-based absolute quantification of protein abundances in Saccharomyces cerevisiae.

Protein quantification via label-free mass spectrometry (MS) has become an increasingly popular method for predicting genome-wide absolute protein abundances. A known caveat of this approach, however, is the poor technical reproducibility, that is, how consistent predictions are when the same sample is measured repeatedly. Here, we measured proteomics data for Saccharomyces cerevisiae with both biological and inter-batch technical triplicates, to analyze both accuracy and precision of protein quantification via MS. Moreover, we analyzed how these metrics vary when applying different methods for converting MS intensities to absolute protein abundances. We demonstrate that our simple normalization and rescaling approach can perform as accurately, yet more precisely, than methods which rely on external standards. Additionally, we show that inter-batch reproducibility is worse than biological reproducibility for all evaluated methods. These results offer a new benchmark for assessing MS data quality for protein quantification, while also underscoring current limitations in this approach.

Author Correction: A consensus S. cerevisiae metabolic model Yeast8 and its ecosystem for comprehensively probing cellular metabolism.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Combining mechanistic and machine learning models for predictive engineering and optimization of tryptophan metabolism.

Through advanced mechanistic modeling and the generation of large high-quality datasets, machine learning is becoming an integral part of understanding and engineering living systems. Here we show that mechanistic and machine learning models can be combined to enable accurate genotype-to-phenotype predictions. We use a genome-scale model to pinpoint engineering targets, efficient library construction of metabolic pathway designs, and high-throughput biosensor-enabled screening for training diverse machine learning algorithms. From a single data-generation cycle, this enables successful forward engineering of complex aromatic amino acid metabolism in yeast, with the best machine learning-guided design recommendations improving tryptophan titer and productivity by up to 74 and 43%, respectively, compared to the best designs used for algorithm training. Thus, this study highlights the power of combining mechanistic and machine learning models to effectively direct metabolic engineering efforts.

Model-Assisted Fine-Tuning of Central Carbon Metabolism in Yeast through dCas9-Based Regulation.

Engineering Saccharomyces cerevisiae for industrial-scale production of valuable chemicals involves extensive modulation of its metabolism. Here, we identified novel gene expression fine-tuning set-ups to enhance endogenous metabolic fluxes toward increasing levels of acetyl-CoA and malonyl-CoA. dCas9-based transcriptional regulation was combined together with a malonyl-CoA responsive intracellular biosensor to select for beneficial set-ups. The candidate genes for screening were predicted using a genome-scale metabolic model, and a gRNA library targeting a total of 168 selected genes was designed. After multiple rounds of fluorescence-activated cell sorting and library sequencing, the gRNAs that were functional and increased flux toward malonyl-CoA were assessed for their efficiency to enhance 3-hydroxypropionic acid (3-HP) production. 3-HP production was significantly improved upon fine-tuning genes involved in providing malonyl-CoA precursors, cofactor supply, as well as chromatin remodeling.

A consensus S. cerevisiae metabolic model Yeast8 and its ecosystem for comprehensively probing cellular metabolism.

Genome-scale metabolic models (GEMs) represent extensive knowledgebases that provide a platform for model simulations and integrative analysis of omics data. This study introduces Yeast8 and an associated ecosystem of models that represent a comprehensive computational resource for performing simulations of the metabolism of Saccharomyces cerevisiae--an important model organism and widely used cell-factory. Yeast8 tracks community development with version control, setting a standard for how GEMs can be continuously updated in a simple and reproducible way. We use Yeast8 to develop the derived models panYeast8 and coreYeast8, which in turn enable the reconstruction of GEMs for 1,011 different yeast strains. Through integration with enzyme constraints (ecYeast8) and protein 3D structures (proYeast8DB), Yeast8 further facilitates the exploration of yeast metabolism at a multi-scale level, enabling prediction of how single nucleotide variations translate to phenotypic traits.

SLIMEr: probing flexibility of lipid metabolism in yeast with an improved constraint-based modeling framework.

BACKGROUND: A recurrent problem in genome-scale metabolic models (GEMs) is to correctly represent lipids as biomass requirements, due to the numerous of possible combinations of individual lipid species and the corresponding lack of fully detailed data. In this study we present SLIMEr, a formalism for correctly representing lipid requirements in GEMs using commonly available experimental data. RESULTS: SLIMEr enhances a GEM with mathematical constructs where we Split Lipids Into Measurable Entities (SLIME reactions), in addition to constraints on both the lipid classes and the acyl chain distribution. By implementing SLIMEr on the consensus GEM of Saccharomyces cerevisiae, we can represent accurate amounts of lipid species, analyze the flexibility of the resulting distribution, and compute the energy costs of moving from one metabolic state to another. CONCLUSIONS: The approach shows potential for better understanding lipid metabolism in yeast under different conditions. SLIMEr is freely available at https://github.com/SysBioChalmers/SLIMEr .

RAVEN 2.0: A versatile toolbox for metabolic network reconstruction and a case study on Streptomyces coelicolor.

RAVEN is a commonly used MATLAB toolbox for genome-scale metabolic model (GEM) reconstruction, curation and constraint-based modelling and simulation. Here we present RAVEN Toolbox 2.0 with major enhancements, including: (i) de novo reconstruction of GEMs based on the MetaCyc pathway database; (ii) a redesigned KEGG-based reconstruction pipeline; (iii) convergence of reconstructions from various sources; (iv) improved performance, usability, and compatibility with the COBRA Toolbox. Capabilities of RAVEN 2.0 are here illustrated through de novo reconstruction of GEMs for the antibiotic-producing bacterium Streptomyces coelicolor. Comparison of the automated de novo reconstructions with the iMK1208 model, a previously published high-quality S. coelicolor GEM, exemplifies that RAVEN 2.0 can capture most of the manually curated model. The generated de novo reconstruction is subsequently used to curate iMK1208 resulting in Sco4, the most comprehensive GEM of S. coelicolor, with increased coverage of both primary and secondary metabolism. This increased coverage allows the use of Sco4 to predict novel genome editing targets for optimized secondary metabolites production. As such, we demonstrate that RAVEN 2.0 can be used not only for de novo GEM reconstruction, but also for curating existing models based on up-to-date databases. Both RAVEN 2.0 and Sco4 are distributed through GitHub to facilitate usage and further development by the community (https://github.com/SysBioChalmers/RAVEN and https://github.com/SysBioChalmers/Streptomyces_coelicolor-GEM).

Improving the phenotype predictions of a yeast genome-scale metabolic model by incorporating enzymatic constraints.

Genome-scale metabolic models (GEMs) are widely used to calculate metabolic phenotypes. They rely on defining a set of constraints, the most common of which is that the production of metabolites and/or growth are limited by the carbon source uptake rate. However, enzyme abundances and kinetics, which act as limitations on metabolic fluxes, are not taken into account. Here, we present GECKO, a method that enhances a GEM to account for enzymes as part of reactions, thereby ensuring that each metabolic flux does not exceed its maximum capacity, equal to the product of the enzyme's abundance and turnover number. We applied GECKO to a Saccharomyces cerevisiae GEM and demonstrated that the new model could correctly describe phenotypes that the previous model could not, particularly under high enzymatic pressure conditions, such as yeast growing on different carbon sources in excess, coping with stress, or overexpressing a specific pathway. GECKO also allows to directly integrate quantitative proteomics data; by doing so, we significantly reduced flux variability of the model, in over 60% of metabolic reactions. Additionally, the model gives insight into the distribution of enzyme usage between and within metabolic pathways. The developed method and model are expected to increase the use of model-based design in metabolic engineering.

Absolute Quantification of Protein and mRNA Abundances Demonstrate Variability in Gene-Specific Translation Efficiency in Yeast.

Protein synthesis is the most energy-consuming process in a proliferating cell, and understanding what controls protein abundances represents a key question in biology and biotechnology. We quantified absolute abundances of 5,354 mRNAs and 2,198 proteins in Saccharomyces cerevisiae under ten environmental conditions and protein turnover for 1,384 proteins under a reference condition. The overall correlation between mRNA and protein abundances across all conditions was low (0.46), but for differentially expressed proteins (n = 202), the median mRNA-protein correlation was 0.88. We used these data to model translation efficiencies and found that they vary more than 400-fold between genes. Non-linear regression analysis detected that mRNA abundance and translation elongation were the dominant factors controlling protein synthesis, explaining 61% and 15% of its variance. Metabolic flux balance analysis further showed that only mitochondrial fluxes were positively associated with changes at the transcript level. The present dataset represents a crucial expansion to the current resources for future studies on yeast physiology.

Genome scale models of yeast: towards standardized evaluation and consistent omic integration.

Genome scale models (GEMs) have enabled remarkable advances in systems biology, acting as functional databases of metabolism, and as scaffolds for the contextualization of high-throughput data. In the case of Saccharomyces cerevisiae (budding yeast), several GEMs have been published and are currently used for metabolic engineering and elucidating biological interactions. Here we review the history of yeast's GEMs, focusing on recent developments. We study how these models are typically evaluated, using both descriptive and predictive metrics. Additionally, we analyze the different ways in which all levels of omics data (from gene expression to flux) have been integrated in yeast GEMs. Relevant conclusions and current challenges for both GEM evaluation and omic integration are highlighted.

Construction of robust dynamic genome-scale metabolic model structures of Saccharomyces cerevisiae through iterative re-parameterization.

Dynamic flux balance analysis (dFBA) has been widely employed in metabolic engineering to predict the effect of genetic modifications and environmental conditions in the cell׳s metabolism during dynamic cultures. However, the importance of the model parameters used in these methodologies has not been properly addressed. Here, we present a novel and simple procedure to identify dFBA parameters that are relevant for model calibration. The procedure uses metaheuristic optimization and pre/post-regression diagnostics, fixing iteratively the model parameters that do not have a significant role. We evaluated this protocol in a Saccharomyces cerevisiae dFBA framework calibrated for aerobic fed-batch and anaerobic batch cultivations. The model structures achieved have only significant, sensitive and uncorrelated parameters and are able to calibrate different experimental data. We show that consumption, suboptimal growth and production rates are more useful for calibrating dynamic S. cerevisiae metabolic models than Boolean gene expression rules, biomass requirements and ATP maintenance.