Showing papers by "Bas Teusink published in 2019"
TL;DR: A computational approach is applied to find properties that distinguish a species from others in the community, as well as to follow individual species in a community to demonstrate how the distinguishing properties of an organism lead to experimentally testable hypotheses concerning the niche and the interactions with other species.
Abstract: Microbial life usually takes place in a community where individuals interact, by competition for nutrients, cross-feeding, inhibition by end-products, but also by their spatial distribution. Lactic acid bacteria are prominent members of microbial communities responsible for food fermentations. Their niche in a community depends on their own properties as well as those of the other species. Here, we apply a computational approach, which uses only genomic and metagenomic information and functional annotation of genes, to find properties that distinguish a species from others in the community, as well as to follow individual species in a community. We analyzed isolated and sequenced strains from a kefir community, and metagenomes from wine fermentations. We demonstrate how the distinguishing properties of an organism lead to experimentally testable hypotheses concerning the niche and the interactions with other species. We observe, for example, that L. kefiranofaciens, a dominant organism in kefir, stands out among the Lactobacilli because it potentially has more amino acid auxotrophies. Using metagenomic analysis of industrial wine fermentations we investigate the role of an inoculated L. plantarum in malolactic fermentation. We observed that L. plantarum thrives better on white than on red wine fermentations and has the largest number of phosphotransferase system among the bacteria observed in the wine communities. Also, L. plantarum together with Pantoea, Erwinia, Asaia, Gluconobacter, and Komagataeibacter genera had the highest number of genes involved in biosynthesis of amino acids.
199 citations
TL;DR: Model builders should carefully choose a tool (or combinations of tools) depending on the intended use of the metabolic model, as this benchmark study shows that none of the tools outperforms the others in all the defined features.
Abstract: Several genome-scale metabolic reconstruction software platforms have been developed and are being continuously updated. These tools have been widely applied to reconstruct metabolic models for hundreds of microorganisms ranging from important human pathogens to species of industrial relevance. However, these platforms, as yet, have not been systematically evaluated with respect to software quality, best potential uses and intrinsic capacity to generate high-quality, genome-scale metabolic models. It is therefore unclear for potential users which tool best fits the purpose of their research. In this work, we perform a systematic assessment of current genome-scale reconstruction software platforms. To meet our goal, we first define a list of features for assessing software quality related to genome-scale reconstruction. Subsequently, we use the feature list to evaluate the performance of each tool. To assess the similarity of the draft reconstructions to high-quality models, we compare each tool’s output networks with that of the high-quality, manually curated, models of Lactobacillus plantarum and Bordetella pertussis, representatives of gram-positive and gram-negative bacteria, respectively. We additionally compare draft reconstructions with a model of Pseudomonas putida to further confirm our findings. We show that none of the tools outperforms the others in all the defined features. Model builders should carefully choose a tool (or combinations of tools) depending on the intended use of the metabolic model. They can use this benchmark study as a guide to select the best tool for their research. Finally, developers can also benefit from this evaluation by getting feedback to improve their software.
137 citations
TL;DR: This work proves mathematically that the resulting optimal metabolic flux distribution is described by a limited number of subnetworks, known as Elementary Flux Modes (EFMs), and finds that the maximal number of flux-carrying EFMs is determined only by the number of imposed constraints on enzyme expression, not by the size, kinetics or topology of the network.
Abstract: Growth rate is a near-universal selective pressure across microbial species. High growth rates require hundreds of metabolic enzymes, each with different nonlinear kinetics, to be precisely tuned within the bounds set by physicochemical constraints. Yet, the metabolic behaviour of many species is characterized by simple relations between growth rate, enzyme expression levels and metabolic rates. We asked if this simplicity could be the outcome of optimisation by evolution. Indeed, when the growth rate is maximized-in a static environment under mass-conservation and enzyme expression constraints-we prove mathematically that the resulting optimal metabolic flux distribution is described by a limited number of subnetworks, known as Elementary Flux Modes (EFMs). We show that, because EFMs are the minimal subnetworks leading to growth, a small active number automatically leads to the simple relations that are measured. We find that the maximal number of flux-carrying EFMs is determined only by the number of imposed constraints on enzyme expression, not by the size, kinetics or topology of the network. This minimal-EFM extremum principle is illustrated in a graphical framework, which explains qualitative changes in microbial behaviours, such as overflow metabolism and co-consumption, and provides a method for identification of the enzyme expression constraints that limit growth under the prevalent conditions. The extremum principle applies to all microorganisms that are selected for maximal growth rates under protein concentration constraints, for example the solvent capacities of cytosol, membrane or periplasmic space.
82 citations
TL;DR: This work characterised 27 FPs in vivo using Saccharomyces cerevisiae as model organism, and codon optimized the best performing FPs for optimal expression in yeast, and found that codon-optimization alters FP characteristics.
Abstract: Fluorescent proteins (FPs) are widely used in many organisms, but are commonly characterised in vitro. However, the in vitro properties may poorly reflect in vivo performance. Therefore, we characterised 27 FPs in vivo using Saccharomyces cerevisiae as model organism. We linked the FPs via a T2A peptide to a control FP, producing equimolar expression of the 2 FPs from 1 plasmid. Using this strategy, we characterised the FPs for brightness, photostability, photochromicity and pH-sensitivity, achieving a comprehensive in vivo characterisation. Many FPs showed different in vivo properties compared to existing in vitro data. Additionally, various FPs were photochromic, which affects readouts due to complex bleaching kinetics. Finally, we codon optimized the best performing FPs for optimal expression in yeast, and found that codon-optimization alters FP characteristics. These FPs improve experimental signal readout, opening new experimental possibilities. Our results may guide future studies in yeast that employ fluorescent proteins.
62 citations
TL;DR: A genome-scale metabolic model of dairy-origin L. cremoris ATCC 19254 was reconstructed to explain the energetics and redox state mechanisms of the organism in full detail, and it was shown that flavor formation only occurs under carbon and ATP excess.
Abstract: Leuconostoc mesenteroides subsp. cremoris is an obligate heterolactic fermentative lactic acid bacterium that is mostly used in industrial dairy fermentations. The phosphoketolase pathway (PKP) is a unique feature of the obligate heterolactic fermentation, which leads to the production of lactate, ethanol, and/or acetate, and the final product profile of PKP highly depends on the energetics and redox state of the organism. Another characteristic of the L. mesenteroides subsp. cremoris is the production of aroma compounds in dairy fermentation, such as in cheese production, through the utilization of citrate. Considering its importance in dairy fermentation, a detailed metabolic characterization of the organism is necessary for its more efficient use in the industry. To this aim, a genome-scale metabolic model of dairy-origin L. mesenteroides subsp. cremoris ATCC 19254 (iLM.c559) was reconstructed to explain the energetics and redox state mechanisms of the organism in full detail. The model includes 559 genes governing 1088 reactions between 1129 metabolites, and the reactions cover citrate utilization and citrate-related flavor metabolism. The model was validated by simulating co-metabolism of glucose and citrate and comparing the in silico results to our experimental results. Model simulations further showed that, in co-metabolism of citrate and glucose, no flavor compounds were produced when citrate could stimulate the formation of biomass. Significant amounts of flavor metabolites (e.g., diacetyl and acetoin) were only produced when citrate could not enhance growth, which suggests that flavor formation only occurs under carbon and ATP excess. The effects of aerobic conditions and different carbon sources on product profiles and growth were also investigated using the reconstructed model. The analyses provided further insights for the growth stimulation and flavor formation mechanisms of the organism.
23 citations
TL;DR: The reconstruction of the A. pasteurianus 386B genome-scale metabolic model revealed knowledge gaps concerning the metabolism of this strain, especially related to the biosynthesis of its cell envelope and the presence or absence of metabolite transporters.
Abstract: Acetobacter pasteurianus 386B is a candidate functional starter culture for the cocoa bean fermentation process. To allow in silico simulations of its related metabolism in response to different environmental conditions, a genome-scale metabolic model for A. pasteurianus 386B was reconstructed. This is the first genome-scale metabolic model reconstruction for a member of the genus Acetobacter. The metabolic network reconstruction process was based on extensive genome re-annotation and comparative genomics analyses. The information content related to the functional annotation of metabolic enzymes and transporters was placed in a metabolic context by exploring and curating a Pathway/Genome Database of A. pasteurianus 386B using the Pathway Tools software. Metabolic reactions and curated gene-protein-reaction associations were bundled into a genome-scale metabolic model of A. pasteurianus 386B, named iAp386B454, containing 454 genes, 322 reactions, and 296 metabolites embedded in two cellular compartments. The reconstructed model was validated by performing growth experiments in a defined medium, which revealed that lactic acid as the sole carbon source could sustain growth of this strain. Further, the reconstruction of the A. pasteurianus 386B genome-scale metabolic model revealed knowledge gaps concerning the metabolism of this strain, especially related to the biosynthesis of its cell envelope and the presence or absence of metabolite transporters.
20 citations
TL;DR: A systematic assessment of the current genome-scale metabolic reconstruction software platforms shows that none of the tools outperforms the others in all the defined features and that model builders should carefully choose a tool (or combinations of tools) depending on the intended use of the metabolic model.
Abstract: Several genome-scale metabolic reconstruction software platforms have been developed and are being continuously updated. These tools have been widely applied to reconstruct metabolic models for hundreds of microorganisms ranging from important human pathogens to species of industrial relevance. However, these platforms, as yet, have not been systematically evaluated with respect to software quality, best potential uses and intrinsic capacity to generate high-quality, genome-scale metabolic models. It is therefore unclear for potential users which tool best fits the purpose of their research. In this work, we performed a systematic assessment of the current genome-scale reconstruction software platforms. To meet our goal, we first defined a list of features for assessing software quality related to genome-scale reconstruction, which we expect to be useful for the potential users of these tools. Subsequently, we used the feature list to evaluate the performance of each tool. In order to assess the similarity of the draft reconstructions to high-quality models, we compared each tool’s output networks with that of the high-quality, manually curated, models of Lactobacillus plantarum and Bordetella pertussis, representatives of gram-positive and gram-negative bacteria, respectively. We showed that none of the tools outperforms the others in all the defined features and that model builders should carefully choose a tool (or combinations of tools) depending on the intended use of the metabolic model. Author Summary Metabolic networks that comprise biochemical reactions at genome-scale have become very useful to study and predict the phenotype of important microorganisms. Several software platforms exist to build these metabolic networks. Based on different approaches and utilizing a variety of databases it is, unfortunately, unclear what are the best scenarios to use each of these tools. Hence, to understand the potential uses of these tools, we created a list of relevant features for metabolic reconstruction and we evaluated the tools in all these categories. Here, we show that none of the tools is better than the other in all the evaluated categories; instead, each tool is more suitable for particular purposes. Therefore, users should carefully select the tool(s) that best fit the purpose of their research. This is the first time these tools are systematically evaluated and this overview can be used as a guide for selecting the correct tool(s) for each case.
15 citations
TL;DR: This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display, and demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints.
Abstract: A central theme in (micro)biology is understanding the molecular basis of fitness i.e. which strategies are successful under which conditions; how do organisms implement such strategies at the molecular level; and which constraints shape the trade-offs between alternative strategies. Highly standardized microbial laboratory evolution experiments are ideally suited to approach these questions. For example, prolonged chemostats provide a constant environment in which the growth rate can be set, and the adaptive process of the organism to such environment can be subsequently characterized. We performed parallel laboratory evolution of Lactococcus lactis in chemostats varying the quantitative value of the selective pressure by imposing two different growth rates. A mutation in one specific amino acid residue of the global transcriptional regulator of carbon metabolism, CcpA, was selected in all of the evolution experiments performed. We subsequently showed that this mutation confers predictable fitness improvements at other glucose-limited growth rates as well. In silico protein structural analysis of wild type and evolved CcpA, as well as biochemical and phenotypic assays, provided the underpinning molecular mechanisms that resulted in the specific reprogramming favored in constant environments. This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display. It demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints, by fine-tuning the expression of a wide range of cellular processes.
14 citations
TL;DR: The theory starts from basic biochemical and evolutionary considerations, and describes unicellular life, both in growth-promoting and in stress-inducing environments, in terms of EGMs, the universal building blocks of self-fabrication and a cell’s phenotype.
Abstract: A major aim of biology is to predict phenotype from genotype. Here we ask if we can describe all possible molecular states (phenotypes) for a cell that fabricates itself at a constant rate, given its enzyme kinetics and the stoichiometry of all reactions (the genotype). For this, we must understand the autocatalytic process of cellular growth which is inherently nonlinear: steady-state self-fabrication requires a cell to synthesize all of its components, including metabolites, enzymes and ribosomes, in the proportions that exactly match its own composition -- the growth demand thus depends on the cellular composition. Simultaneously, the concentrations of these components should be tuned to accomplish this synthesis task -- the cellular composition thus depends on the growth demand. We here derive a theory that describes all phenotypes that solve this circular problem; the basic equations show how the concentrations of all cellular components and reaction rates must be balanced to get a constant self-fabrication rate. All phenotypes can be described as a combination of one or more minimal building blocks, which we call Elementary Growth Modes (EGMs). EGMs can be used as the theoretical basis for all models that explicitly model self-fabrication, such as the currently popular Metabolism and Expression models. We then used our theory to make concrete biological predictions: we find that natural selection for maximal growth rate drives microorganisms to states of minimal phenotypic complexity: only one EGM will be active when cellular growth rate is maximised. The phenotype of a cell is only extended with one more EGM whenever growth becomes limited by an additional biophysical constraint, such as a limited solvent capacity of a cellular compartment. Our theory starts from basic biochemical and evolutionary considerations, and describes unicellular life, both in growth-promoting and in stress-inducing environments, in terms of EGMs, the universal building blocks of self-fabrication and a cell9s phenotype.
13 citations
TL;DR: It is shown how kefir, a natural milk-fermenting community, realizes stable coexistence through spatio-temporal orchestration of species and metabolite dynamics, and how microbes poorly suited for milk survive in, and even dominate the community through metabolic cooperation and uneven partitioning between the grain and the liquid phase.
Abstract: Microbial communities in nature often feature complex compositional dynamics yet also stable coexistence of diverse species. The mechanistic underpinnings of such dynamic stability remain unclear as system-wide studies have been limited to small engineered communities or synthetic assemblies. Here we show how kefir, a natural milk-fermenting community, realizes stable coexistence through spatio-temporal orchestration of species and metabolite dynamics. During milk fermentation, kefir grains (a polysaccharide matrix synthesized by kefir microbes) grow in mass but remain unchanged in composition. In contrast, the milk is colonized in a dynamic fashion with early members opening metabolic niches for the followers. Through large-scale mapping of metabolic preferences and inter-species interactions, we show how microbes poorly suited for milk survive in, and even dominate, the community through metabolic cooperation and uneven partitioning between the grain and the liquid phase. Overall, our findings reveal how spatio-temporal dynamics promote stable coexistence and have implications for deciphering and modulating complex microbial ecosystems.
11 citations
TL;DR: All models predict overflow metabolism when two, model-specific, growth-limiting constraints are hit, consistent with recent theory.
Abstract: In each environment, living cells can express different metabolic pathways that support growth. The criteria that determine which pathways are selected remain unclear. One recurrent selection is overflow metabolism: the seemingly wasteful, simultaneous usage of an efficient and an inefficient pathway, shown for example in E. coli, S. cerevisiae and cancer cells. Many different models, based on different assumptions, can reproduce this observation. Therefore, they provide no conclusive evidence which mechanism is causing overflow metabolism. We compare the mathematical structure of these models. Although ranging from Flux Balance Analyses to self-fabricating Metabolism and Expression models, we could rewrite all models into one standard form. We conclude that all models predict overflow metabolism when two, model-specific, growth-limiting constraints are hit. This is consistent with recent theory. Thus, identifying these two constraints is essential for understanding overflow metabolism. We list all the imposed constraints by these models, so that they can hopefully be tested in future experiments.
TL;DR: It is shown that pregrowth conditions largely affect ATP dynamics during transitions, and single-cell analyses showed a large variety in ATP responses, which implies large differences of glycolytic startup between individual cells, confirming the need for new tools to study dynamic responses of individual cells in dynamic environments.
Abstract: Adenosine 5-triphosphate (ATP) is the main free energy carrier in metabolism. In budding yeast, shifts to glucose-rich conditions cause dynamic changes in ATP levels, but it is unclear how heterogeneous these dynamics are at the single-cell level. Furthermore, pH also changes and affects readout of fluorescence-based biosensors for single-cell measurements. To measure ATP changes reliably in single yeast cells, we developed yAT1.03, an adapted version of the AT1.03 ATP biosensor, that is pH-insensitive. We show that pregrowth conditions largely affect ATP dynamics during transitions. Moreover, single-cell analyses showed a large variety in ATP responses, which implies large differences of glycolytic startup between individual cells. We found three clusters of dynamic responses, and show that a small subpopulation of wild type cells reached an imbalanced state during glycolytic startup, characterised by low ATP levels. These results confirm the need for new tools to study dynamic responses of individual cells in dynamic environments.
TL;DR: YEPAC as mentioned in this paper is a FRET-based biosensor for cAMP measurements in yeast, which is used with flow cytometry for high-throughput single cell-level quantification during dynamic changes.
Abstract: The cAMP-PKA signalling cascade in budding yeast regulates adaptation to changing environments. We developed yEPAC, a FRET-based biosensor for cAMP measurements in yeast. We used this sensor with flow cytometry for high-throughput single cell-level quantification during dynamic changes in response to sudden nutrient transitions. We found that the characteristic cAMP peak differentiates between different carbon source transitions, and is rather homogenous among single-cells, especially for transitions to glucose. The peaks are mediated by a combination of extracellular sensing and intracellular metabolism. Moreover, the cAMP peak follows Webers law; its height scales with the relative, and not the absolute, change in glucose. Lastly, our results suggest that the cAMP peak height conveys information about prospective growth rates. In conclusion, our yEPAC-sensor makes possible new avenues for understanding yeast physiology, signalling and metabolic adaptation.