Comparative Analysis of Yeast Metabolic Network Models Highlights Progress, Opportunities for Metabolic Reconstruction

Emergence of Orchestrated and Dynamic Metabolism of Saccharomyces cerevisiae

Microbial metabolism is a fundamental cellular process that involves many biochemical events and is distinguished by its emergent properties. While the molecular details of individual reactions have been increasingly elucidated, it is not well understood how these reactions are quantitatively orchestrated to produce collective cellular behaviors. Here we developed a coarse-grained, systems, and dynamic mathematical framework, which integrates metabolic reactions with signal transduction and gene regulation to dissect the emergent metabolic traits of Saccharomyces cerevisiae. Our framework mechanistically captures a set of characteristic cellular behaviors, including the Crabtree effect, diauxic shift, diauxic lag time, and differential growth under nutrient-altered environments. It also allows modular expansion for zooming in on specific pathways for detailed metabolic profiles. This study provides a systems mathematical framework for yeast metabolic behaviors, providing insights into yeast physiology and metabolic engineering.

Reconstruction of the genome-scale metabolic network model of Sinorhizobium fredii CCBAU45436 for free-living and symbiotic states

Sinorhizobium fredii CCBAU45436 is an excellent rhizobium that plays an important role in agricultural production. However, there still needs more comprehensive understanding of the metabolic system of S. fredii CCBAU45436, which hinders its application in agriculture. Therefore, based on the first-generation metabolic model iCC541 we developed a new genome-scale metabolic model iAQY970, which contains 970 genes, 1,052 reactions, 942 metabolites and is scored 89% in the MEMOTE test. Cell growth phenotype predicted by iAQY970 is 81.7% consistent with the experimental data. The results of mapping the proteome data under free-living and symbiosis conditions to the model showed that the biomass production rate in the logarithmic phase was faster than that in the stable phase, and the nitrogen fixation efficiency of rhizobia parasitized in cultivated soybean was higher than that in wild-type soybean, which was consistent with the actual situation. In the symbiotic condition, there are 184 genes that would affect growth, of which 94 are essential; In the free-living condition, there are 143 genes that influence growth, of which 78 are essential. Among them, 86 of the 94 essential genes in the symbiotic condition were consistent with the prediction of iCC541, and 44 essential genes were confirmed by literature information; meanwhile, 30 genes were identified by DEG and 33 genes were identified by Geptop. In addition, we extracted four key nitrogen fixation modules from the model and predicted that sulfite reductase (EC 1.8.7.1) and nitrogenase (EC 1.18.6.1) as the target enzymes to enhance nitrogen fixation by MOMA, which provided a potential focus for strain optimization. Through the comprehensive metabolic model, we can better understand the metabolic capabilities of S. fredii CCBAU45436 and make full use of it in the future.

Integration of graph neural networks and genome-scale metabolic models for predicting gene essentiality

Genome-scale metabolic models are powerful tools for understanding cellular physiology. Flux balance analysis (FBA), in particular, is an optimization-based approach widely employed for predicting metabolic phenotypes. In model microbes such as Escherichia coli, FBA has been successful at predicting essential genes, i.e. those genes that impair survival when deleted. A central assumption in this approach is that both wild type and deletion strains optimize the same fitness objective. Although the optimality assumption may hold for the wild type metabolic network, deletion strains are not subject to the same evolutionary pressures and knock-out mutants may steer their metabolism to meet other objectives for survival. Here, we present FlowGAT, a hybrid FBA-machine learning strategy for predicting essentiality directly from wild type metabolic phenotypes. The approach is based on graph-structured representation of metabolic fluxes predicted by FBA, where nodes correspond to enzymatic reactions and edges quantify the propagation of metabolite mass flow between a reaction and its neighbours. We integrate this information into a graph neural network that can be trained on knock-out fitness assay data. Comparisons across different model architectures reveal that FlowGAT predictions for E. coli are close to those of FBA for several growth conditions. This suggests that essentiality of enzymatic genes can be predicted by exploiting the inherent network structure of metabolism. Our approach demonstrates the benefits of combining the mechanistic insights afforded by genome-scale models with the ability of deep learning to infer patterns from complex datasets.

The origins of the Guinness stout yeast

Beer is made via the fermentation of an aqueous extract predominantly composed of malted barley flavoured with hops. The transforming microorganism is typically a single strain of Saccharomyces cerevisiae, and for the majority of major beer brands the yeast strain is a unique component. The present yeast used to make Guinness stout brewed in Dublin, Ireland, can be traced back to 1903, but its origins are unknown. To that end, we used Illumina and Nanopore sequencing to generate whole-genome sequencing data for a total of 22 S. cerevisiae yeast strains: 16 from the Guinness collection and 6 other historical Irish brewing. The origins of the Guinness yeast were determined with a SNP-based analysis, demonstrating that the Guinness strains occupy a distinct group separate from other historical Irish brewing yeasts. Assessment of chromosome number, copy number variation and phenotypic evaluation of key brewing attributes established Guinness yeast-specific SNPs but no specific chromosomal amplifications. Our analysis also demonstrated the effects of yeast storage on phylogeny. Altogether, our results suggest that the Guinness yeast used today is related to the first deposited Guinness yeast; the 1903 Watling Laboratory Guinness yeast.

Building Yarrowia lipolytica Cell Factories for Advanced Biomanufacturing: Challenges and Solutions

Microbial cell factories have shown great potential for industrial production with the benefit of being environmentally friendly and sustainable. Yarrowia lipolytica is a promising and superior non-model host for biomanufacturing due to its cumulated advantages compared to model microorganisms, such as high fluxes of metabolic precursors (acetyl-CoA and malonyl-CoA) and its naturally hydrophobic microenvironment. However, although diverse compounds have been synthesized in Y. lipolytica cell factories, most of the relevant studies have not reached the level of industrialization and commercialization due to a number of remaining challenges, including unbalanced metabolic flux, conflict between cell growth and product synthesis, and cytotoxic effects. Here, various metabolic engineering strategies for solving the challenges are summarized, which is developing fast and extremely conducive to rational design and reconstruction of robust Y. lipolytica cell factories for advanced biomanufacturing. Finally, future engineering efforts for enhancing the production efficiency of this platform strain are highlighted.

Mitochondrial Metabolism in the Spotlight: Maintaining Balanced RNAP III Activity Ensures Cellular Homeostasis

RNA polymerase III (RNAP III) holoenzyme activity and the processing of its products have been linked to several metabolic dysfunctions in lower and higher eukaryotes. Alterations in the activity of RNAP III-driven synthesis of non-coding RNA cause extensive changes in glucose metabolism. Increased RNAP III activity in the S. cerevisiae maf1Δ strain is lethal when grown on a non-fermentable carbon source. This lethal phenotype is suppressed by reducing tRNA synthesis. Neither the cause of the lack of growth nor the underlying molecular mechanism have been deciphered, and this area has been awaiting scientific explanation for a decade. Our previous proteomics data suggested mitochondrial dysfunction in the strain. Using model mutant strains maf1Δ (with increased tRNA abundance) and rpc128-1007 (with reduced tRNA abundance), we collected data showing major changes in the TCA cycle metabolism of the mutants that explain the phenotypic observations. Based on 13C flux data and analysis of TCA enzyme activities, the present study identifies the flux constraints in the mitochondrial metabolic network. The lack of growth is associated with a decrease in TCA cycle activity and downregulation of the flux towards glutamate, aspartate and phosphoenolpyruvate (PEP), the metabolic intermediate feeding the gluconeogenic pathway. rpc128-1007, the strain that is unable to increase tRNA synthesis due to a mutation in the C128 subunit, has increased TCA cycle activity under non-fermentable conditions. To summarize, cells with non-optimal activity of RNAP III undergo substantial adaptation to a new metabolic state, which makes them vulnerable under specific growth conditions. Our results strongly suggest that balanced, non-coding RNA synthesis that is coupled to glucose signaling is a fundamental requirement to sustain a cell's intracellular homeostasis and flexibility under changing growth conditions. The presented results provide insight into the possible role of RNAP III in the mitochondrial metabolism of other cell types.

Microbial coexistence in the rhizosphere and the promotion of plant stress resistance: A review

Plants can recruit soil microorganisms into the rhizosphere when experiencing various environmental stresses, including biotic (e.g., insect pests) and abiotic (e.g., heavy metal pollution, droughts, floods, and salinity) stresses. However, species coexistence in plant resistance has not received sufficient attention. Current research on microbial coexistence is only at the community scale, and there is a limited understanding of the interaction patterns between species, especially microbe‒microbe interactions. The relevant interaction patterns are limited to a few model strains. The coexisting microbial communities form a stable system involving complex nutritional competition, metabolic exchange, and even interdependent interactions. This pattern of coexistence can ultimately enhance plant stress tolerance. Hence, a systematic understanding of the coexistence pattern of rhizosphere microorganisms under stress is essential for the precise development and utilization of synthetic microbial communities and the achievement of efficient ecological control. Here, we integrated current analytical methods and introduced several new experimental methods to elucidate rhizosphere microbial coexistence patterns. Some advancements (e.g., network analysis, coculture experiments, and synthetic communities) that can be applied to plant stress resistance are also updated. This review aims to summarize the key role and potential application prospects of microbial coexistence in the resistance of plants to environmental stresses. Our suggestions, enhancing plant resistance with coexisting microbes, would allow us to gain further knowledge on plant-microbial and microbial-microbial functions, and facilitate translation to more effective measures.

Pareto Optimal Metabolic Engineering for the Growth-coupled Overproduction of Sustainable Chemicals

Our research aims to help industrial biotechnology develop a sustainable economy using green technology based on microorganisms and synthetic biology through two case studies that improve metabolic capacity in yeast models Yarrowia lipolytica (Y. lipolytica) and Saccharomyces cerevisiae (S. cerevisiae). We aim to increase the production capacity of beta‐carotene (β‐carotene) and succinic acid, which are among the highest market demands due to their versatile use in numerous consumer products. We performed simulations to identify in silico ranking of strains based on multiple objectives: the growth rate of yeast microorganisms, the number of used chromosomes, and the production capability of β‐carotene (for Y. lipolytica) and succinate (for S. cerevisiae). Our multiobjective optimization methodology identified notable gene deletions by searching a vast solution space to highlight near‐optimal strains on Pareto Fronts, balancing the above‐cited three objectives. Moreover, preserving the metabolic constraints and the essential genes, this study produced robust results: seven significant strains of Y. lipolytica and seven strains of S. cerevisiae. We examined gene knockout to study the function of genes and pathways. In fact, by studying the frequently silenced genes, we found that when the GPH1 gene is knocked out in S. cerevisiae, the isocitrate lyase enzyme is activated, which converts the isocitrate into succinate. Our goals are to simplify and facilitate the in vitro processes. Hence, we present strains with the least possible number of knockout genes and solutions in which the genes are turned off on the same chromosome. Therefore, we present results where the constraints mentioned above are met, like the strains where only two genes are switched off and other strains where half of the knockout genes are on the same chromosome. This study offers solutions for developing an efficient in vitro mutagenesis for microorganisms and demonstrates the efficiency of multiobjective optimization in automatizing metabolic engineering processes.

Genome-scale modeling of yeast metabolism: retrospectives and perspectives

Yeasts have been widely used for production of bread, beer and wine, as well as for production of bioethanol, but they have also been designed as cell factories to produce various chemicals, advanced biofuels and recombinant proteins. To systematically understand and rationally engineer yeast metabolism, genome-scale metabolic models (GEMs) have been reconstructed for the model yeast Saccharomyces cerevisiae and nonconventional yeasts. Here, we review the historical development of yeast GEMs together with their recent applications, including metabolic flux prediction, cell factory design, culture condition optimization and multi-yeast comparative analysis. Furthermore, we present an emerging effort, namely the integration of proteome constraints into yeast GEMs, resulting in models with improved performance. At last, we discuss challenges and perspectives on the development of yeast GEMs and the integration of proteome constraints.

Rhizospheric microbiome: Bio-based emerging strategies for sustainable agriculture development and future perspectives

In the light of intensification of cropping practices and changing climatic conditions, nourishing a growing global population requires optimizing environmental sustainability and reducing ecosystem impacts of food production. The use of microbiological systems to ameliorate the agricultural production in a sustainable and eco-friendly way is widespread accepted as a future key-technology. However, the multitude of interaction possibilities between the numerous beneficial microbes and plants in their habitat calls for systematic analysis and management of the rhizospheric microbiome. This review exploits present and future strategies for rhizospheric microbiome management with the aim to generate a comprehensive understanding of the known tools and techniques. Significant information on the structure and dynamics of rhizospheric microbiota of isolated microbial communities is now available. These microbial communities have beneficial effects including increased plant growth, essential nutrient acquisition, pathogens tolerance, and increased abiotic as well as biotic stress tolerance such as drought, temperature, salinity and antagonistic activities against the phyto-pathogens. A better and comprehensive understanding of the various effects and microbial interactions can be gained by application of molecular approaches as extraction of DNA/RNA and other biochemical markers to analyze microbial soil diversity. Novel techniques like interactome network analysis and split-ubiquitin system framework will enable to gain more insight into communication and interactions between the proteins from microbes and plants. The aim of the analysis tasks leads to the novel approach of Rhizosphere microbiome engineering. The capability of forming the rhizospheric microbiome in a defined way will allow combining several microbes (e.g. bacteria and fungi) for a given environment (soil type and climatic zone) in order to exert beneficial influences on specific plants. This integration will require a large-scale effort among academic researchers, industry researchers and farmers to understand and manage interactions of plant-microbiomes within modern farming systems, and is clearly a multi-domain approach and can be mastered only jointly by microbiology, mathematics and information technology. These innovations will open up a new avenue for designing and implementing intensive farming microbiome management approaches to maximize resource productivity and stress tolerance of agro-ecosystems, which in return will create value to the increasing worldwide population, for both food production and consumption.

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome-scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence-level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.

Tolerance to NADH/NAD+ imbalance anticipates aging and anti-aging interventions

Redox couples coordinate cellular function, but the consequences of their imbalances are unclear. This is somewhat associated with the limitations of their experimental quantification. Here we circumvent these difficulties by presenting an approach that characterizes fitness-based tolerance profiles to redox couple imbalances using an in silico representation of metabolism. Focusing on the NADH/NAD⁺ redox couple in yeast, we demonstrate that reductive disequilibria generate metabolic syndromes comparable to those observed in cancer cells. The tolerance of yeast mutants to redox disequilibrium can also explain 30% of the variability in their experimentally measured chronological lifespan. Moreover, by predicting the significance of some metabolites to help stand imbalances, we correctly identify nutrients underlying mechanisms of pathology, lifespan-protecting molecules, or caloric restriction mimetics. Tolerance to redox imbalances becomes, in this way, a sound framework to recognize properties of the aging phenotype while providing a consistent biological rationale to assess anti-aging interventions.

•

We simulate how imbalances in NADH/NAD+ ratio modify cellular metabolic behavior

•

This reveals a mechanism to understand metabolic alterations at low growth rates

•

Tolerance to imbalance explains experimentally measured lifespan in yeast

•

We predict lifespan-protecting metabolites in yeast, animal, and human models

Evaluating accessibility, usability and interoperability of genome-scale metabolic models for diverse yeasts species

Metabolic network reconstructions have become an important tool for probing cellular metabolism in the field of systems biology. They are used as tools for quantitative prediction but also as scaffolds for further knowledge contextualization. The yeast Saccharomyces cerevisiae was one of the first organisms for which a genome-scale metabolic model (GEM) was reconstructed, in 2003, and since then 45 metabolic models have been developed for a wide variety of relevant yeasts species. A systematic evaluation of these models revealed that-despite this long modeling history-the sequential process of tracing model files, setting them up for basic simulation purposes and comparing them across species and even different versions, is still not a generalizable task. These findings call the yeast modeling community to comply to standard practices on model development and sharing in order to make GEMs accessible and useful for a wider public.

Rare genetic variants affecting urine metabolite levels link population variation to inborn errors of metabolism

Metabolite levels in urine may provide insights into genetic mechanisms shaping their related pathways. We therefore investigate the cumulative contribution of rare, exonic genetic variants on urine levels of 1487 metabolites and 53,714 metabolite ratios among 4864 GCKD study participants. Here we report the detection of 128 significant associations involving 30 unique genes, 16 of which are known to underlie inborn errors of metabolism. The 30 genes are strongly enriched for shared expression in liver and kidney (odds ratio = 65, p-FDR = 3e-7), with hepatocytes and proximal tubule cells as driving cell types. Use of UK Biobank whole-exome sequencing data links genes to diseases connected to the identified metabolites. In silico constraint-based modeling of gene knockouts in a virtual whole-body, organ-resolved metabolic human correctly predicts the observed direction of metabolite changes, highlighting the potential of linking population genetics to modeling. Our study implicates candidate variants and genes for inborn errors of metabolism.

Recent advances in systems and synthetic biology approaches for developing novel cell-factories in non-conventional yeasts

Microbial bioproduction of chemicals, proteins, and primary metabolites from cheap carbon sources is currently an advancing area in industrial research. The model yeast, Saccharomyces cerevisiae, is a well-established biorefinery host that has been used extensively for commercial manufacturing of bioethanol from myriad carbon sources. However, its Crabtree-positive nature often limits the use of this organism for the biosynthesis of commercial molecules that do not belong in the fermentative pathway. To avoid extensive strain engineering of S. cerevisiae for the production of metabolites other than ethanol, non-conventional yeasts can be selected as hosts based on their natural capacity to produce desired commodity chemicals. Non-conventional yeasts like Kluyveromyces marxianus, K. lactis, Yarrowia lipolytica, Pichia pastoris, Scheffersomyces stipitis, Hansenula polymorpha, and Rhodotorula toruloides have been considered as potential industrial eukaryotic hosts owing to their desirable phenotypes such as thermotolerance, assimilation of a wide range of carbon sources, as well as ability to secrete high titers of protein and lipid. However, the advanced metabolic engineering efforts in these organisms are still lacking due to the limited availability of systems and synthetic biology methods like in silico models, well-characterised genetic parts, and optimized genome engineering tools. This review provides an insight into the recent advances and challenges of systems and synthetic biology as well as metabolic engineering endeavours towards the commercial usage of non-conventional yeasts. Particularly, the approaches in emerging non-conventional yeasts for the production of enzymes, therapeutic proteins, lipids, and metabolites for commercial applications are extensively discussed here. Various attempts to address current limitations in designing novel cell factories have been highlighted that include the advances in the fields of genome-scale metabolic model reconstruction, flux balance analysis, 'omics'-data integration into models, genome-editing toolkit development, and rewiring of cellular metabolisms for desired chemical production. Additionally, the understanding of metabolic networks using ¹³C-labelling experiments as well as the utilization of metabolomics in deciphering intracellular fluxes and reactions have also been discussed here. Application of cutting-edge nuclease-based genome editing platforms like CRISPR/Cas9, and its optimization towards efficient strain engineering in non-conventional yeasts have also been described. Additionally, the impact of the advances in promising non-conventional yeasts for efficient commercial molecule synthesis has been meticulously reviewed. In the future, a cohesive approach involving systems and synthetic biology will help in widening the horizon of the use of unexplored non-conventional yeast species towards industrial biotechnology.

Advances in biotechnological production of santalenes and santalols

Sandalwood essential oil has been widely used not only as natural medicines but also in perfumery and food industries, with sesquiterpenoids as its major components including (Z)- α-santalol and (Z)-β-santalol and so on. The mature heartwoods of Santalum album, Santalum austrocaledonicum and Santalum spicatum are the major plant resources for extracting sandalwood essential oil, which have been overexploited. Synthetic biology approaches have been successfully applied to produce natural products on large scale. In this review, we summarize biosynthetic enzymes of santalenes and santalols, including various santalene synthases (STSs) and cytochrome P450 monooxygenases (CYPs), and then highlight the advances of biotechnological production of santalenes and santalols in heterologous hosts, especially metabolic engineering strategies for constructing santalene- and santalol-producing Saccharomyces cerevisiae.

CEG 2.0: an updated database of clusters of essential genes including eukaryotic organisms

Essential genes are key elements for organisms to maintain their living. Building databases that store essential genes in the form of homologous clusters, rather than storing them as a singleton, can provide more enlightening information such as the general essentiality of homologous genes in multiple organisms. In 2013, the first database to store prokaryotic essential genes in clusters, CEG (Clusters of Essential Genes), was constructed. Afterward, the amount of available data for essential genes increased by a factor >3 since the last revision. Herein, we updated CEG to version 2, including more prokaryotic essential genes (from 16 gene datasets to 29 gene datasets) and newly added eukaryotic essential genes (nine species), specifically the human essential genes of 12 cancer cell lines. For prokaryotes, information associated with drug targets, such as protein structure, ligand–protein interaction, virulence factor and matched drugs, is also provided. Finally, we provided the service of essential gene prediction for both prokaryotes and eukaryotes. We hope our updated database will benefit more researchers in drug targets and evolutionary genomics.

Database URL:http://cefg.uestc.cn/ceg

Curation and Analysis of a Saccharomyces cerevisiae Genome-Scale Metabolic Model for Predicting Production of Sensory Impact Molecules under Enological Conditions

One approach for elucidating strain-to-strain metabolic differences is the use of genome-scale metabolic models (GSMMs). To date GSMMs have not focused on the industrially important area of flavor production and, as such; do not cover all the pathways relevant to flavor formation in yeast. Moreover, current models for Saccharomyces cerevisiae generally focus on carbon-limited and/or aerobic systems, which is not pertinent to enological conditions. Here, we curate a GSMM (iWS902) to expand on the existing Ehrlich pathway and ester formation pathways central to aroma formation in industrial winemaking, in addition to the existing sulfur metabolism and medium-chain fatty acid (MCFA) pathways that also contribute to production of sensory impact molecules. After validating the model using experimental data, we predict key differences in metabolism for a strain (EC 1118) in two distinct growth conditions, including differences for aroma impact molecules such as acetic acid, tryptophol, and hydrogen sulfide. Additionally, we propose novel targets for metabolic engineering for aroma profile modifications employing flux variability analysis with the expanded GSMM. The model provides mechanistic insights into the key metabolic pathways underlying aroma formation during alcoholic fermentation and provides a potential framework to contribute to new strategies to optimize the aroma of wines.

Genetic buffering and potentiation in metabolism

Cells adjust their metabolism in response to mutations, but how this reprogramming depends on the genetic context is not well known. Specifically, the absence of individual enzymes can affect reprogramming, and thus the impact of mutations in cell growth. Here, we examine this issue with an in silico model of Saccharomyces cerevisiae's metabolism. By quantifying the variability in the growth rate of 10000 different mutant metabolisms that accumulated changes in their reaction fluxes, in the presence, or absence, of a specific enzyme, we distinguish a subset of modifier genes serving as buffers or potentiators of variability. We notice that the most potent modifiers refer to the glycolysis pathway and that, more broadly, they show strong pleiotropy and epistasis. Moreover, the evidence that this subset depends on the specific growing condition strengthens its systemic underpinning, a feature only observed before in a toy model of a gene-regulatory network. Some of these enzymes also modulate the effect that biochemical noise and environmental fluctuations produce in growth. Thus, the reorganization of metabolism induced by mutations has not only direct physiological implications but also transforms the influence that other mutations have on growth. This is a general result with implications in the development of cancer therapies based on metabolic inhibitors.

Genome-Scale Identification of Essential Metabolic Processes for Targeting the Plasmodium Liver Stage

Plasmodium gene functions in mosquito and liver stages remain poorly characterized due to limitations in the throughput of phenotyping at these stages. To fill this gap, we followed more than 1,300 barcoded P. berghei mutants through the life cycle. We discover 461 genes required for efficient parasite transmission to mosquitoes through the liver stage and back into the bloodstream of mice. We analyze the screen in the context of genomic, transcriptomic, and metabolomic data by building a thermodynamic model of P. berghei liver-stage metabolism, which shows a major reprogramming of parasite metabolism to achieve rapid growth in the liver. We identify seven metabolic subsystems that become essential at the liver stages compared with asexual blood stages: type II fatty acid synthesis and elongation (FAE), tricarboxylic acid, amino sugar, heme, lipoate, and shikimate metabolism. Selected predictions from the model are individually validated in single mutants to provide future targets for drug development.

•

1,342 barcoded P. berghei knockout (KO) mutants analyzed for stage-specific phenotypes

•

Life-stage-specific metabolic models reveal reprogramming of cellular function

•

High agreement between blood/liver stage metabolic models and genetic screening data

•

Essential metabolic pathways for parasite development and mechanistic origin revealed

Into the wild: new yeast genomes from natural environments and new tools for their analysis

Genomic studies of yeasts from the wild have increased considerably in the past few years. This revolution has been fueled by advances in high-throughput sequencing technologies and a better understanding of yeast ecology and phylogeography, especially for biotechnologically important species. The present review aims to first introduce new bioinformatic tools available for the generation and analysis of yeast genomes. We also assess the accumulated genomic data of wild isolates of industrially relevant species, such as Saccharomyces spp., which provide unique opportunities to further investigate the domestication processes associated with the fermentation industry and opportunistic pathogenesis. The availability of genome sequences of other less conventional yeasts obtained from the wild has also increased substantially, including representatives of the phyla Ascomycota (e.g. Hanseniaspora) and Basidiomycota (e.g. Phaffia). Here, we review salient examples of both fundamental and applied research that demonstrate the importance of continuing to sequence and analyze genomes of wild yeasts.

redLips: a comprehensive mechanistic model of the lipid metabolic network of yeast

Over the last decades, yeast has become a key model organism for the study of lipid biochemistry. Because the regulation of lipids has been closely linked to various physiopathologies, the study of these biomolecules could lead to new diagnostics and treatments. Before the field can reach this point, however, sufficient tools for integrating and analyzing the ever-growing availability of lipidomics data will need to be developed. To this end, genome-scale models (GEMs) of metabolic networks are useful tools, though their large size and complexity introduces too much uncertainty in the accuracy of predicted outcomes. Ideally, therefore, a model for studying lipids would contain only the pathways required for the proper analysis of these biomolecules, but would not be an ad hoc reduction. We hereby present a metabolic model that focuses on lipid metabolism constructed through the integration of detailed lipid pathways into an already existing GEM of Saccharomyces cerevisiae. Our model was then systematically reduced around the subsystems defined by these pathways to provide a more manageable model size for complex studies. We show that this model is as consistent and inclusive as other yeast GEMs regarding the focus and detail on the lipid metabolism, and can be used as a scaffold for integrating lipidomics data to improve predictions in studies of lipid-related biological functions.

Green production of silybin and isosilybin by merging metabolic engineering approaches and enzymatic catalysis

Silymarin extracted from milk thistle seeds, is used for treating hepatic diseases. Silybin and isosilybin are its main components, and synthesized from coupling of taxifolin and coniferyl alcohol. Here, the biosynthetic pathways of taxifolin and coniferyl alcohol were reconstructed in Saccharomyces cerevisiae for the first time. To alleviate substantial burden caused by a great deal of genetic manipulation, expression of the enzymes (e.g. ZWF1, TYR1 and ARO8) playing multiple roles in the relevant biosynthetic pathways was selectively optimized. The strain YT1035 overexpressing seven heterologous enzymes and five native enzymes and the strain YC1053 overexpressing seven heterologous enzymes and four native enzymes, respectively produce 336.8 mg/L taxifolin and 201.1 mg/L coniferyl alcohol. Silybin and isosilybin are synthesized from taxifolin and coniferyl alcohol under catalysis of APX1t (the truncated milk thistle peroxidase), with a yield of 62.5%. This study demonstrates an approach for producing silybin and isosilybin from glucose for the first time.

Functional and Computational Genomics Reveal Unprecedented Flexibility in Stage-Specific Toxoplasma Metabolism

To survive and proliferate in diverse host environments with varying nutrient availability, the obligate intracellular parasite Toxoplasma gondii reprograms its metabolism. We have generated and curated a genome-scale metabolic model (iTgo) for the fast-replicating tachyzoite stage, harmonized with experimentally observed phenotypes. To validate the importance of four metabolic pathways predicted by the model, we have performed in-depth in vitro and in vivo phenotyping of mutant parasites including targeted metabolomics and CRISPR-Cas9 fitness screening of all known metabolic genes. This led to unexpected insights into the remarkable flexibility of the parasite, addressing the dependency on biosynthesis or salvage of fatty acids (FAs), purine nucleotides (AMP and GMP), a vitamin (pyridoxal-5P), and a cofactor (heme) in both the acute and latent stages of infection. Taken together, our experimentally validated metabolic network leads to a deeper understanding of the parasite's biology, opening avenues for the development of therapeutic intervention against apicomplexans.

Effects of MCHM on yeast metabolism

On January 2014 approximately 10,000 gallons of crude 4-Methylcyclohexanemethanol (MCHM) and propylene glycol phenol ether (PPH) were accidentally released into the Elk River, West Virginia, contaminating the tap water of around 300,000 residents. Crude MCHM is an industrial chemical used as flotation reagent to clean coal. At the time of the spill, MCHM's toxicological data were limited, an issue that has been addressed by different studies focused on understanding the immediate and long-term effects of MCHM on human health and the environment. Using S. cerevisiae as a model organism we study the effect of acute exposure to crude MCHM on metabolism. Yeasts were treated with MCHM 550 ppm in YPD for 30 minutes. Polar and lipid metabolites were extracted from cells by a chloroform-methanol-water mixture. The extracts were then analyzed by direct injection ESI-MS and by GC-MS. The metabolomics analysis was complemented with flux balance analysis simulations done with genome-scale metabolic network models (GSMNM) of MCHM treated vs non-treated control. We integrated the effect of MCHM on yeast gene expression from RNA-Seq data within these GSMNM. A total of 215 and 73 metabolites were identified by the ESI-MS and GC-MS procedures, respectively. From these 26 and 23 relevant metabolites were selected from ESI-MS and GC-MS respectively, for 49 unique compounds. MCHM induced amino acid accumulation, via its effects on amino acid metabolism, as well as a potential impairment of ribosome biogenesis. MCHM affects phospholipid biosynthesis, with a potential impact on the biophysical properties of yeast cellular membranes. The FBA simulations were able to reproduce the deleterious effect of MCHM on cellular growth and suggest that the effect of MCHM on ubiquinol:ferricytochrome c reductase reaction, caused by the under-expression of CYT1 gene, could be the driven force behind the observed effect on yeast metabolism and growth.

Metabolic fluxes-oriented control of bioreactors: a novel approach to tune micro-aeration and substrate feeding in fermentations

Background

Fine-tuning the aeration for cultivations when oxygen-limited conditions are demanded (such as the production of vaccines, isobutanol, 2–3 butanediol, acetone, and bioethanol) is still a challenge in the area of bioreactor automation and advanced control. In this work, an innovative control strategy based on metabolic fluxes was implemented and evaluated in a case study: micro-aerated ethanol fermentation.

Results

The experiments were carried out in fed-batch mode, using commercial Saccharomyces cerevisiae, defined medium, and glucose as carbon source. Simulations of a genome-scale metabolic model for Saccharomyces cerevisiae were used to identify the range of oxygen and substrate fluxes that would maximize ethanol fluxes. Oxygen supply and feed flow rate were manipulated to control oxygen and substrate fluxes, as well as the respiratory quotient (RQ). The performance of the controlled cultivation was compared to two other fermentation strategies: a conventional “Brazilian fuel-ethanol plant” fermentation and a strictly anaerobic fermentation (with ultra-pure nitrogen used as the inlet gas). The cultivation carried out under the proposed control strategy showed the best average volumetric ethanol productivity (7.0 g L⁻¹ h⁻¹), with a final ethanol concentration of 87 g L⁻¹ and yield of 0.46 g_ethanol g_substrate⁻¹. The other fermentation strategies showed lower yields (close to 0.40 g_ethanol g_substrate⁻¹) and ethanol productivity around 4.0 g L⁻¹ h⁻¹.

Conclusion

The control system based on fluxes was successfully implemented. The proposed approach could also be adapted to control several bioprocesses that require restrict aeration.

Closed-loop cycles of experiment design, execution, and learning accelerate systems biology model development in yeast

Systems biology involves the development of large computational models of biological systems. The radical improvement of systems biology models will necessarily involve the automation of model improvement cycles. We present here a general approach to automating systems biology model improvement. Humans are eukaryotic organisms, and the yeast Saccharomyces cerevisiae is widely used in biology as a “model” for eukaryotic cells. The yeast diauxic shift is the most studied cellular transformation. We combined multiple software tools with integrated laboratory robotics to execute three semiautomated cycles of diauxic shift model improvement. All the experiments were formalized and communicated to a cloud laboratory automation system (Eve) for execution. The resulting improved model is relevant to understanding cancer, the immune system, and aging.

One of the most challenging tasks in modern science is the development of systems biology models: Existing models are often very complex but generally have low predictive performance. The construction of high-fidelity models will require hundreds/thousands of cycles of model improvement, yet few current systems biology research studies complete even a single cycle. We combined multiple software tools with integrated laboratory robotics to execute three cycles of model improvement of the prototypical eukaryotic cellular transformation, the yeast (Saccharomyces cerevisiae) diauxic shift. In the first cycle, a model outperforming the best previous diauxic shift model was developed using bioinformatic and systems biology tools. In the second cycle, the model was further improved using automatically planned experiments. In the third cycle, hypothesis-led experiments improved the model to a greater extent than achieved using high-throughput experiments. All of the experiments were formalized and communicated to a cloud laboratory automation system (Eve) for automatic execution, and the results stored on the semantic web for reuse. The final model adds a substantial amount of knowledge about the yeast diauxic shift: 92 genes (+45%), and 1,048 interactions (+147%). This knowledge is also relevant to understanding cancer, the immune system, and aging. We conclude that systems biology software tools can be combined and integrated with laboratory robots in closed-loop cycles.

Flux balance analysis with or without molecular crowding fails to predict two thirds of experimentally observed epistasis in yeast

Computational predictions of double gene knockout effects by flux balance analysis (FBA) have been used to characterized genome-wide patterns of epistasis in microorganisms. However, it is unclear how in silico predictions are related to in vivo epistasis, as FBA predicted only a minority of experimentally observed genetic interactions between non-essential metabolic genes in yeast. Here, we perform a detailed comparison of yeast experimental epistasis data to predictions generated with different constraint-based metabolic modeling algorithms. The tested methods comprise standard FBA; a variant of MOMA, which was specifically designed to predict fitness effects of non-essential gene knockouts; and two alternative implementations of FBA with macro-molecular crowding, which account approximately for enzyme kinetics. The number of interactions uniquely predicted by one method is typically larger than its overlap with any alternative method. Only 20% of negative and 10% of positive interactions jointly predicted by all methods are confirmed by the experimental data; almost all unique predictions appear to be false. More than two thirds of epistatic interactions are undetectable by any of the tested methods. The low prediction accuracies indicate that the physiology of yeast double metabolic gene knockouts is dominated by processes not captured by current constraint-based analysis methods.

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.

Systems Biology of Plant-Microbiome Interactions

In natural environments, plants are exposed to diverse microbiota that they interact with in complex ways. While plant-pathogen interactions have been intensely studied to understand defense mechanisms in plants, many microbes and microbial communities can have substantial beneficial effects on their plant host. Such beneficial effects include improved acquisition of nutrients, accelerated growth, resilience against pathogens, and improved resistance against abiotic stress conditions such as heat, drought, and salinity. However, the beneficial effects of bacterial strains or consortia on their host are often cultivar and species specific, posing an obstacle to their general application. Remarkably, many of the signals that trigger plant immune responses are molecularly highly similar and often identical in pathogenic and beneficial microbes. Thus, it is unclear what determines the outcome of a particular microbe-host interaction and which factors enable plants to distinguish beneficials from pathogens. To unravel the complex network of genetic, microbial, and metabolic interactions, including the signaling events mediating microbe-host interactions, comprehensive quantitative systems biology approaches will be needed.

Metabolism and Development during Conidial Germination in Response to a Carbon-Nitrogen-Rich Synthetic or a Natural Source of Nutrition in Neurospora crassa

Fungal spores germinate and undergo vegetative growth, leading to either asexual or sexual reproductive dispersal. Previous research has indicated that among developmental regulatory genes, expression is conserved across nutritional environments, whereas pathways for carbon and nitrogen metabolism appear highly responsive-perhaps to accommodate differential nutritive processing. To comprehensively investigate conidial germination and the adaptive life history decision-making underlying these two modes of reproduction, we profiled transcription of Neurospora crassa germinating on two media: synthetic Bird medium, designed to promote asexual reproduction; and a natural maple sap medium, on which both asexual reproduction and sexual reproduction manifest. A later start to germination but faster development was observed on synthetic medium. Metabolic genes exhibited altered expression in response to nutrients-at least 34% of the genes in the genome were significantly downregulated during the first two stages of conidial germination on synthetic medium. Knockouts of genes exhibiting differential expression across development altered germination and growth rates, as well as in one case causing abnormal germination. A consensus Bayesian network of these genes indicated especially tight integration of environmental sensing, asexual and sexual development, and nitrogen metabolism on a natural medium, suggesting that in natural environments, a more dynamic and tentative balance of asexual and sexual development may be typical of N. crassa colonies.IMPORTANCE One of the most remarkable successes of life is its ability to flourish in response to temporally and spatially varying environments. Fungi occupy diverse ecosystems, and their sensitivity to these environmental changes often drives major fungal life history decisions, including the major switch from vegetative growth to asexual or sexual reproduction. Spore germination comprises the first and simplest stage of vegetative growth. We examined the dependence of this early life history on the nutritional environment using genome-wide transcriptomics. We demonstrated that for developmental regulatory genes, expression was generally conserved across nutritional environments, whereas metabolic gene expression was highly labile. The level of activation of developmental genes did depend on current nutrient conditions, as did the modularity of metabolic and developmental response network interactions. This knowledge is critical to the development of future technologies that could manipulate fungal growth for medical, agricultural, or industrial purposes.

OptRAM: In-silico strain design via integrative regulatory-metabolic network modeling

The ultimate goal of metabolic engineering is to produce desired compounds on an industrial scale in a cost effective manner. To address challenges in metabolic engineering, computational strain optimization algorithms based on genome-scale metabolic models have increasingly been used to aid in overproducing products of interest. However, most of these strain optimization algorithms utilize a metabolic network alone, with few approaches providing strategies that also include transcriptional regulation. Moreover previous integrated approaches generally require a pre-existing regulatory network. In this study, we developed a novel strain design algorithm, named OptRAM (Optimization of Regulatory And Metabolic Networks), which can identify combinatorial optimization strategies including overexpression, knockdown or knockout of both metabolic genes and transcription factors. OptRAM is based on our previous IDREAM integrated network framework, which makes it able to deduce a regulatory network from data. OptRAM uses simulated annealing with a novel objective function, which can ensure a favorable coupling between desired chemical and cell growth. The other advance we propose is a systematic evaluation metric of multiple solutions, by considering the essential genes, flux variation, and engineering manipulation cost. We applied OptRAM to generate strain designs for succinate, 2,3-butanediol, and ethanol overproduction in yeast, which predicted high minimum predicted target production rate compared with other methods and previous literature values. Moreover, most of the genes and TFs proposed to be altered by OptRAM in these scenarios have been validated by modification of the exact genes or the target genes regulated by the TFs, for overproduction of these desired compounds by in vivo experiments cataloged in the LASER database. Particularly, we successfully validated the predicted strain optimization strategy for ethanol production by fermentation experiment. In conclusion, OptRAM can provide a useful approach that leverages an integrated transcriptional regulatory network and metabolic network to guide metabolic engineering applications.

Functional genetic discovery of enzymes using full-scan mass spectrometry metabolomics

Our understanding of metabolic networks is incomplete, and new enzymatic activities await discovery in well-studied organisms. Mass spectrometric measurement of cellular metabolites reveals compounds inside cells that are unexplained by current maps of metabolic reactions, and existing computational models are unable to account for all activities observed within cells. Additional large-scale genetic and biochemical approaches are required to elucidate metabolic gene function. We have used full-scan mass spectrometry metabolomics of polar small molecules to examine deletion mutants of candidate enzymes in the model yeast Saccharomyces cerevisiae. We report the identification of 25 genes whose deletion results in focal metabolic changes consistent with loss of enzymatic activity and describe the informatic approaches used to enrich for candidate enzymes from uncharacterized open reading frames. Triumphs and pitfalls of metabolic phenotyping screens are discussed, including estimates of the frequency of uncharacterized eukaryotic genes that affect metabolism and key issues to consider when searching for new enzymatic functions in other organisms.

Extension of the yeast metabolic model to include iron metabolism and its use to estimate global levels of iron-recruiting enzyme abundance from cofactor requirements

Metabolic networks adapt to changes in their environment by modulating the activity of their enzymes and transporters; often by changing their abundance. Understanding such quantitative changes can shed light onto how metabolic adaptation works, or how it can fail and lead to a metabolically dysfunctional state. We propose a strategy to quantify metabolic protein requirements for cofactor‐utilising enzymes and transporters through constraint‐based modelling. The first eukaryotic genome‐scale metabolic model to comprehensively represent iron metabolism was constructed, extending the most recent community model of the Saccharomyces cerevisiae metabolic network. Partial functional impairment of the genes involved in the maturation of iron‐sulphur (Fe‐S) proteins was investigated employing the model and the in silico analysis revealed extensive rewiring of the fluxes in response to this functional impairment, despite its marginal phenotypic effect. The optimal turnover rate of enzymes bearing ion cofactors can be determined via this novel approach; yeast metabolism, at steady state, was determined to employ a constant turnover of its iron‐recruiting enzyme at a rate of 3.02 × 10 ⁻¹¹ mmol·(g biomass) ⁻¹·h  ⁻¹.

Yeast Genome-Scale Metabolic Models for Simulating Genotype-Phenotype Relations

Understanding genotype-phenotype dependency is a universal aim for all life sciences. While the complete genotype-phenotype relations remain challenging to resolve, metabolic phenotypes are moving within the reach through genome-scale metabolic model simulations. Genome-scale metabolic models are available for commonly investigated yeasts, such as model eukaryote and domesticated fermentation species Saccharomyces cerevisiae, and automatic reconstruction methods facilitate obtaining models for any sequenced species. The models allow for investigating genotype-phenotype relations through simulations simultaneously considering the effects of nutrient availability, and redox and energy homeostasis in cells. Genome-scale models also offer frameworks for omics data integration to help to uncover how the translation of genotypes to the apparent phenotypes is regulated at different levels. In this chapter, we provide an overview of the yeast genome-scale metabolic models and the simulation approaches for using these models to interrogate genotype-phenotype relations. We review the methodological approaches according to the underlying biological reasoning in order to inspire formulating novel questions and applications that the genome-scale metabolic models could contribute to. Finally, we discuss current challenges and opportunities in the genome-scale metabolic model simulations.

Genome-Scale Metabolic Modeling from Yeast to Human Cell Models of Complex Diseases: Latest Advances and Challenges

Genome-scale metabolic models (GEMs) are mathematical models that enable systematic analysis of metabolism. This modeling concept has been applied to study the metabolism of many organisms including the eukaryal model organism, the yeast Saccharomyces cerevisiae, that also serves as an important cell factory for production of fuels and chemicals. With the application of yeast GEMs, our knowledge of metabolism is increasing. Therefore, GEMs have also been used for modeling human cells to study metabolic diseases. Here we introduce the concept of GEMs and provide a protocol for reconstructing GEMs. Besides, we show the historic development of yeast GEMs and their applications. Also, we review human GEMs as well as their uses in the studies of complex diseases.

Predicting growth rate from gene expression

Growth rate is one of the most important and most complex phenotypic characteristics of unicellular microorganisms, which determines the genetic mutations that dominate at the population level, and ultimately whether the population will survive. Translating changes at the genetic level to their growth-rate consequences remains a subject of intense interest, since such a mapping could rationally direct experiments to optimize antibiotic efficacy or bioreactor productivity. In this work, we directly map transcriptional profiles to growth rates by gathering published gene-expression data from Escherichia coli and Saccharomyces cerevisiae with corresponding growth-rate measurements. Using a machine-learning technique called k-nearest-neighbors regression, we build a model which predicts growth rate from gene expression. By exploiting the correlated nature of gene expression and sparsifying the model, we capture 81% of the variance in growth rate of the E. coli dataset, while reducing the number of features from >4,000 to 9. In S. cerevisiae, we account for 89% of the variance in growth rate, while reducing from >5,500 dimensions to 18. Such a model provides a basis for selecting successful strategies from among the combinatorial number of experimental possibilities when attempting to optimize complex phenotypic traits like growth rate.

CiliateGEM: an open-project and a tool for predictions of ciliate metabolic variations and experimental condition design

BACKGROUND

The study of cell metabolism is becoming central in several fields such as biotechnology, evolution/adaptation and human disease investigations. Here we present CiliateGEM, the first metabolic network reconstruction draft of the freshwater ciliate Tetrahymena thermophila. We also provide the tools and resources to simulate different growth conditions and to predict metabolic variations. CiliateGEM can be extended to other ciliates in order to set up a meta-model, i.e. a metabolic network reconstruction valid for all ciliates. Ciliates are complex unicellular eukaryotes of presumably monophyletic origin, with a phylogenetic position that is equal from plants and animals. These cells represent a new concept of unicellular system with a high degree of species, population biodiversity and cell complexity. Ciliates perform in a single cell all the functions of a pluricellular organism, including locomotion, feeding, digestion, and sexual processes.

RESULTS

After generating the model, we performed an in-silico simulation with the presence and absence of glucose. The lack of this nutrient caused a 32.1% reduction rate in biomass synthesis. Despite the glucose starvation, the growth did not stop due to the use of alternative carbon sources such as amino acids.

CONCLUSIONS

The future models obtained from CiliateGEM may represent a new approach to describe the metabolism of ciliates. This tool will be a useful resource for the ciliate research community in order to extend these species as model organisms in different research fields. An improved understanding of ciliate metabolism could be relevant to elucidate the basis of biological phenomena like genotype-phenotype relationships, population genetics, and cilia-related disease mechanisms.

Grape and Wine Metabolomics to Develop New Insights Using Untargeted and Targeted Approaches

Chemical analysis of grape juice and wine has been performed for over 50 years in a targeted manner to determine a limited number of compounds using Gas Chromatography, Mass-Spectrometry (GC-MS) and High Pressure Liquid Chromatography (HPLC). Therefore, it only allowed the determination of metabolites that are present in high concentration, including major sugars, amino acids and some important carboxylic acids. Thus, the roles of many significant but less concentrated metabolites during wine making process are still not known. This is where metabolomics shows its enormous potential, mainly because of its capability in analyzing over 1000 metabolites in a single run due to the recent advancements of high resolution and sensitive analytical instruments. Metabolomics has predominantly been adopted by many wine scientists as a hypothesis-generating tool in an unbiased and non-targeted way to address various issues, including characterization of geographical origin (terroir) and wine yeast metabolic traits, determination of biomarkers for aroma compounds, and the monitoring of growth developments of grape vines and grapes. The aim of this review is to explore the published literature that made use of both targeted and untargeted metabolomics to study grapes and wines and also the fermentation process. In addition, insights are also provided into many other possible avenues where metabolomics shows tremendous potential as a question-driven approach in grape and wine research.

Microbial Methylotrophic Metabolism: Recent Metabolic Modeling Efforts and Their Applications In Industrial Biotechnology

Developing methylotrophic bacteria into cell factories that meet the chemical demand of the future could be both economical and environmentally friendly. Methane is not only an abundant, low-cost resource but also a potent greenhouse gas, the capture of which could help to reduce greenhouse gas emissions. Rational strain design workflows rely on the availability of carefully combined knowledge often in the form of genome-scale metabolic models to construct high-producer organisms. In this review, the authors present the most recent genome-scale metabolic models in aerobic methylotrophy and their applications. Further, the authors present models for the study of anaerobic methanotrophy through reverse methanogenesis and suggest organisms that may be of interest for expanding one-carbon industrial biotechnology. Metabolic models of methylotrophs are scarce, yet they are important first steps toward rational strain-design in these organisms.

Population FBA predicts metabolic phenotypes in yeast

Using protein counts sampled from single cell proteomics distributions to constrain fluxes through a genome-scale model of metabolism, Population flux balance analysis (Population FBA) successfully described metabolic heterogeneity in a population of independent Escherichia coli cells growing in a defined medium. We extend the methodology to account for correlations in protein expression arising from the co-regulation of genes and apply it to study the growth of independent Saccharomyces cerevisiae cells in two different growth media. We find the partitioning of flux between fermentation and respiration predicted by our model agrees with recent ¹³C fluxomics experiments, and that our model largely recovers the Crabtree effect (the experimentally known bias among certain yeast species toward fermentation with the production of ethanol even in the presence of oxygen), while FBA without proteomics constraints predicts respirative metabolism almost exclusively. The comparisons to the ¹³C study showed improvement upon inclusion of the correlations and motivated a technique to systematically identify inconsistent kinetic parameters in the literature. The minor secretion fluxes for glycerol and acetate are underestimated by our method, which indicate a need for further refinements to the metabolic model. For yeast cells grown in synthetic defined (SD) medium, the calculated broad distribution of growth rates matches experimental observations from single cell studies, and we characterize several metabolic phenotypes within our modeled populations that make use of diverse pathways. Fast growing yeast cells are predicted to perform significant amount of respiration, use serine-glycine cycle and produce ethanol in mitochondria as opposed to slow growing cells. We use a genetic algorithm to determine the proteomics constraints necessary to reproduce the growth rate distributions seen experimentally. We find that a core set of 51 constraints are essential but that additional constraints are still necessary to recover the observed growth rate distribution in SD medium.

No two living cells are exactly the same. Even cells from a clonal population with identical genomes living in the same environment will express proteins in different numbers simply due to the random nature of the chemistry involved in gene expression. The consequences of this stochastic gene expression are complex and not well understood, especially at the level of large reaction networks like metabolism. Here we investigate how variability in the copy numbers of metabolic enzymes affects how individual cells extract nourishment from their environment and grow. We model 100,000 independent yeast cells, each with their own set of enzyme copy numbers sampled from experimental distributions, and use flux balance analysis (FBA) to compute the optimal way that each cell can use its metabolic pathways—an approach we dubbed Population FBA. We find that enzyme variability gives rise to a wide distribution of growth rates, and several metabolic phenotypes—subpopulations relying on diverse metabolic pathways. Most importantly, we compare the predicted fluxes through the different pathways to experimental values; we find that Population FBA is able to correctly predict Crabtree effect, while traditional FBA, which lacks the proteomics constraints our method imposes, differs both qualitatively and quantitatively from experiment.

Extracellular Microbial Metabolomics: The State of the Art

Microorganisms produce and secrete many primary and secondary metabolites to the surrounding environment during their growth. Therefore, extracellular metabolites provide important information about the changes in microbial metabolism due to different environmental cues. The determination of these metabolites is also comparatively easier than the extraction and analysis of intracellular metabolites as there is no need for cell rupture. Many analytical methods are already available and have been used for the analysis of extracellular metabolites from microorganisms over the last two decades. Here, we review the applications and benefits of extracellular metabolite analysis. We also discuss different sample preparation protocols available in the literature for both types (e.g., metabolites in solution and in gas) of extracellular microbial metabolites. Lastly, we evaluate the authenticity of using extracellular metabolomics data in the metabolic modelling of different industrially important microorganisms.

Genome-scale modeling of yeast: chronology, applications and critical perspectives

Over the last 15 years, several genome-scale metabolic models (GSMMs) were developed for different yeast species, aiding both the elucidation of new biological processes and the shift toward a bio-based economy, through the design of in silico inspired cell factories. Here, an historical perspective of the GSMMs built over time for several yeast species is presented and the main inheritance patterns among the metabolic reconstructions are highlighted. We additionally provide a critical perspective on the overall genome-scale modeling procedure, underlining incomplete model validation and evaluation approaches and the quest for the integration of regulatory and kinetic information into yeast GSMMs. A summary of experimentally validated model-based metabolic engineering applications of yeast species is further emphasized, while the main challenges and future perspectives for the field are finally addressed.