Bas Teusink

A functional genomics strategy that uses metabolome data to reveal the phenotype of silent mutations

A large proportion of the 6,000 genes present in the genome of Saccharomyces cerevisiae, and of those sequenced in other organisms, encode proteins of unknown function. Many of these genes are “silent,” that is, they show no overt phenotype, in terms of growth rate or other fluxes, when they are deleted from the genome. We demonstrate how the intracellular concentrations of metabolites can reveal phenotypes for proteins active in metabolic regulation. Quantification of the change of several metabolite concentrations relative to the concentration change of one selected metabolite can reveal the site of action, in the metabolic network, of a silent gene. In the same way, comprehensive analyses of metabolite concentrations in mutants, providing “metabolic snapshots,” can reveal functions when snapshots from strains deleted for unstudied genes are compared to those deleted for known genes. This approach to functional analysis, using comparative metabolomics, we call FANCY—an abbreviation for functional analysis by co-responses in yeast.

Shifts in growth strategies reflect tradeoffs in cellular economics.

The growth rate‐dependent regulation of cell size, ribosomal content, and metabolic efficiency follows a common pattern in unicellular organisms: with increasing growth rates, cell size and ribosomal content increase and a shift to energetically inefficient metabolism takes place. The latter two phenomena are also observed in fast growing tumour cells and cell lines. These patterns suggest a fundamental principle of design. In biology such designs can often be understood as the result of the optimization of fitness. Here we show that in basic models of self‐replicating systems these patterns are the consequence of maximizing the growth rate. Whereas most models of cellular growth consider a part of physiology, for instance only metabolism, the approach presented here integrates several subsystems to a complete self‐replicating system. Such models can yield fundamentally different optimal strategies. In particular, it is shown how the shift in metabolic efficiency originates from a tradeoff between investments in enzyme synthesis and metabolic yields for alternative catabolic pathways. The models elucidate how the optimization of growth by natural selection shapes growth strategies.

Analysis of growth of Lactobacillus plantarum WCFS1 on a complex medium using a genome-scale metabolic model

A genome-scale metabolic model of the lactic acid bacterium Lactobacillus plantarum WCFS1 was constructed based on genomic content and experimental data. The complete model includes 721 genes, 643 reactions, and 531 metabolites. Different stoichiometric modeling techniques were used for interpretation of complex fermentation data, as L. plantarum is adapted to nutrient-rich environments and only grows in media supplemented with vitamins and amino acids. (i) Based on experimental input and output fluxes, maximal ATP production was estimated and related to growth rate. (ii) Optimization of ATP production further identified amino acid catabolic pathways that were not previously associated with free-energy metabolism. (iii) Genome-scale elementary flux mode analysis identified 28 potential futile cycles. (iv) Flux variability analysis supplemented the elementary mode analysis in identifying parallel pathways, e.g. pathways with identical end products but different co-factor usage. Strongly increased flexibility in the metabolic network was observed when strict coupling between catabolic ATP production and anabolic consumption was relaxed. These results illustrate how a genome-scale metabolic model and associated constraint-based modeling techniques can be used to analyze the physiology of growth on a complex medium rather than a minimal salts medium. However, optimization of biomass formation using the Flux Balance Analysis approach, reported to successfully predict growth rate and by product formation in Escherichia coli and Saccharomyces cerevisiae, predicted too high biomass yields that were incompatible with the observed lactate production. The reason is that this approach assumes optimal efficiency of substrate to biomass conversion, and can therefore not predict the metabolically inefficient lactate formation.

In Silico Reconstruction of the Metabolic Pathways of Lactobacillus plantarum: Comparing Predictions of Nutrient Requirements with Those from Growth Experiments

ABSTRACT On the basis of the annotated genome we reconstructed the metabolic pathways of the lactic acid bacterium Lactobacillus plantarum WCFS1. After automatic reconstruction by the Pathologic tool of Pathway Tools (http://bioinformatics.ai.sri.com/ptools/ ), the resulting pathway-genome database, LacplantCyc, was manually curated extensively. The current database contains refinements to existing routes and new gram-positive bacterium-specific reactions that were not present in the MetaCyc database. These reactions include, for example, reactions related to cell wall biosynthesis, molybdopterin biosynthesis, and transport. At present, LacplantCyc includes 129 pathways and 704 predicted reactions involving some 670 chemical species and 710 enzymes. We tested vitamin and amino acid requirements of L. plantarum experimentally and compared the results with the pathways present in LacplantCyc. In the majority of cases (32 of 37 cases) the experimental results agreed with the final reconstruction. LacplantCyc is the most extensively curated pathway-genome database for gram-positive bacteria and is open to the microbiology community via the World Wide Web (www.lacplantcyc.nl ). It can be used as a reference pathway-genome database for gram-positive microbes in general and lactic acid bacteria in particular.

Modelling strategies for the industrial exploitation of lactic acid bacteria

Lactic acid bacteria (LAB) have a long tradition of use in the food industry, and the number and diversity of their applications has increased considerably over the years. Traditionally, process optimization for these applications involved both strain selection and trial and error. More recently, metabolic engineering has emerged as a discipline that focuses on the rational improvement of industrially useful strains. In the post-genomic era, metabolic engineering increasingly benefits from systems biology, an approach that combines mathematical modelling techniques with functional-genomics data to build models for biological interpretation and — ultimately — prediction. In this review, the industrial applications of LAB are mapped onto available global, genome-scale metabolic modelling techniques to evaluate the extent to which functional genomics and systems biology can live up to their industrial promise.

Accelerating the reconstruction of genome-scale metabolic networks

BackgroundThe genomic information of a species allows for the genome-scale reconstruction of its metabolic capacity. Such a metabolic reconstruction gives support to metabolic engineering, but also to integrative bioinformatics and visualization. Sequence-based automatic reconstructions require extensive manual curation, which can be very time-consuming. Therefore, we present a method to accelerate the time-consuming process of network reconstruction for a query species. The method exploits the availability of well-curated metabolic networks and uses high-resolution predictions of gene equivalency between species, allowing the transfer of gene-reaction associations from curated networks.ResultsWe have evaluated the method using Lactococcus lactis IL1403, for which a genome-scale metabolic network was published recently. We recovered most of the gene-reaction associations (i.e. 74 – 85%) which are incorporated in the published network. Moreover, we predicted over 200 additional genes to be associated to reactions, including genes with unknown function, genes for transporters and genes with specific metabolic reactions, which are good candidates for an extension to the previously published network. In a comparison of our developed method with the well-established approach Pathologic, we predicted 186 additional genes to be associated to reactions. We also predicted a relatively high number of complete conserved protein complexes, which are derived from curated metabolic networks, illustrating the potential predictive power of our method for protein complexes.ConclusionWe show that our methodology can be applied to accelerate the reconstruction of genome-scale metabolic networks by taking optimal advantage of existing, manually curated networks. As orthology detection is the first step in the method, only the translated open reading frames (ORFs) of a newly sequenced genome are necessary to reconstruct a metabolic network. When more manually curated metabolic networks will become available in the near future, the usefulness of our method in network prediction is likely to increase.

Lost in Transition: Start-Up of Glycolysis Yields Subpopulations of Nongrowing Cells

Introduction Cells use multilayered regulatory systems to respond adequately to changing environments or perturbations. Failure in regulation underlies cellular malfunctioning, loss of fitness, or disease. How molecular components dynamically interact to give rise to robust and adaptive responses is not well understood. Here, we studied how the model eukaryote Saccharomyces cerevisiae can cope with transition to high glucose levels, a failure of which results in metabolic malfunctioning and growth arrest. Initiation of glycolysis can have two outcomes. Upon glucose availability, glycolysis can end up in either a functional steady state or an unviable imbalanced state with imbalanced fluxes between ATP-consuming (Vupper) and ATP-producing steps (Vlower). In wild-type yeast, the transient activation of trehalose cycling pushes glycolysis toward the viable steady state. Failure to do so results in metabolic malfunctioning, as observed in mutants in trehalose biosynthesis (tps1Δ). Methods We combined experimental and modeling approaches to unravel the mechanisms used by yeast to cope with sudden glucose availability. We studied growth characteristics and metabolic state at population and single-cell levels (through flow cytometry and colony plating) of the wild type and of mutants unable to transit properly to excess glucose; such mutants are defective in trehalose synthesis, a disaccharide associated with stress resistance. Dynamic 13C tracer enrichment was used to estimate dynamic intracellular fluxes immediately after glucose addition. Mathematical modeling was used to interpret and generalize results and to suggest subsequent experiments. Results The failure to cope with glucose is caused by imbalanced reactions in glycolysis, the essential pathway in energy metabolism in most organisms. In the failure mode, the first steps of glycolysis carry more flux than the downstream steps, resulting in accumulating intermediates at constant low levels of adenosine triphosphate (ATP) and inorganic phosphate. We found that cells with such an unbalanced glycolysis coexist with vital cells with normal glycolytic function. Spontaneous, nongenetic metabolic variability among individual cells determines which state is reached and consequently which cells survive. In mutants of trehalose metabolism, only 0.01% of the cells started to grow on glucose; in the wild type, the success rate was still only 93% (i.e., 7% of wild-type yeast did not properly start up glycolysis). Mathematical models predicted that the dynamics of inorganic phosphate is a key determinant in successful transition to glucose, and that phosphate release through ATP hydrolysis reduces the probability of reaching an imbalanced state. 13C-labeling experiments confirmed the hypothesis that trehalose metabolism constitutes a futile cycle that would provide proper phosphate balance: Upon a glucose pulse, almost 30% of the glucose is transiently shuttled into trehalose metabolism. Discussion Our work reveals how cell fate can be determined by glycolytic dynamics combined with cell heterogeneity purely at the metabolic level. Specific regulatory mechanisms are required to initiate the glycolytic pathway; in yeast, trehalose cycling pushes glycolysis transiently into the right direction, after which cycling stops. The coexistence of two modes of glycolysis—an imbalanced state and the normal functional state—arises from the fundamental design of glycolysis. This makes the imbalanced state a generic risk for humans as well, extending our fundamental knowledge of this central pathway that is dysfunctional in diseases such as diabetes and cancer. Metabolic Heterogeneity We commonly think of genetic or epigenetic sources of variation in cells and individuals. However, biochemical regulatory pathways can potentially also exist in multiple stable states and confer variable phenotypes on cells in a population. Van Heerden et al. (10.1126/science.1245114, published online 16 January) demonstrate such a phenomenon in yeast cells. Two distinct types of cell were observed that differed in the state of glycolysis, the central pathway in energy metabolism for these cells. This allowed some members of a population of cells to survive changes in glucose concentrations, whereas most cells did not. One source of nongenetic variation in yeast can be traced to distinct steady-state levels of glycolysis. Cells need to adapt to dynamic environments. Yeast that fail to cope with dynamic changes in the abundance of glucose can undergo growth arrest. We show that this failure is caused by imbalanced reactions in glycolysis, the essential pathway in energy metabolism in most organisms. The imbalance arises largely from the fundamental design of glycolysis, making this state of glycolysis a generic risk. Cells with unbalanced glycolysis coexisted with vital cells. Spontaneous, nongenetic metabolic variability among individual cells determines which state is reached and, consequently, which cells survive. Transient ATP (adenosine 5′-triphosphate) hydrolysis through futile cycling reduces the probability of reaching the imbalanced state. Our results reveal dynamic behavior of glycolysis and indicate that cell fate can be determined by heterogeneity purely at the metabolic level.

Understanding the Adaptive Growth Strategy of Lactobacillus plantarum by In Silico Optimisation

In the study of metabolic networks, optimization techniques are often used to predict flux distributions, and hence, metabolic phenotype. Flux balance analysis in particular has been successful in predicting metabolic phenotypes. However, an inherent limitation of a stoichiometric approach such as flux balance analysis is that it can predict only flux distributions that result in maximal yields. Hence, previous attempts to use FBA to predict metabolic fluxes in Lactobacillus plantarum failed, as this lactic acid bacterium produces lactate, even under glucose-limited chemostat conditions, where FBA predicted mixed acid fermentation as an alternative pathway leading to a higher yield. In this study we tested, however, whether long-term adaptation on an unusual and poor carbon source (for this bacterium) would select for mutants with optimal biomass yields. We have therefore adapted Lactobacillus plantarum to grow well on glycerol as its main growth substrate. After prolonged serial dilutions, the growth yield and corresponding fluxes were compared to in silico predictions. Surprisingly, the organism still produced mainly lactate, which was corroborated by FBA to indeed be optimal. To understand these results, constraint-based elementary flux mode analysis was developed that predicted 3 out of 2669 possible flux modes to be optimal under the experimental conditions. These optimal pathways corresponded very closely to the experimentally observed fluxes and explained lactate formation as the result of competition for oxygen by the other flux modes. Hence, these results provide thorough understanding of adaptive evolution, allowing in silico predictions of the resulting flux states, provided that the selective growth conditions favor yield optimization as the winning strategy.

Exploring Metabolic Pathway Reconstruction and Genome-Wide Expression Profiling in Lactobacillus reuteri to Define Functional Probiotic Features

The genomes of four Lactobacillus reuteri strains isolated from human breast milk and the gastrointestinal tract have been recently sequenced as part of the Human Microbiome Project. Preliminary genome comparisons suggested that these strains belong to two different clades, previously shown to differ with respect to antimicrobial production, biofilm formation, and immunomodulation. To explain possible mechanisms of survival in the host and probiosis, we completed a detailed genomic comparison of two breast milk–derived isolates representative of each group: an established probiotic strain (L. reuteri ATCC 55730) and a strain with promising probiotic features (L. reuteri ATCC PTA 6475). Transcriptomes of L. reuteri strains in different growth phases were monitored using strain-specific microarrays, and compared using a pan-metabolic model representing all known metabolic reactions present in these strains. Both strains contained candidate genes involved in the survival and persistence in the gut such as mucus-binding proteins and enzymes scavenging reactive oxygen species. A large operon predicted to encode the synthesis of an exopolysaccharide was identified in strain 55730. Both strains were predicted to produce health-promoting factors, including antimicrobial agents and vitamins (folate, vitamin B12). Additionally, a complete pathway for thiamine biosynthesis was predicted in strain 55730 for the first time in this species. Candidate genes responsible for immunomodulatory properties of each strain were identified by transcriptomic comparisons. The production of bioactive metabolites by human-derived probiotics may be predicted using metabolic modeling and transcriptomics. Such strategies may facilitate selection and optimization of probiotics for health promotion, disease prevention and amelioration.

Bet-hedging during bacterial diauxic shift

Significance More than 70 years ago, Monod described the phenomenon of diauxic growth of bacteria, the observation that in the presence of two alternative sugars, cells first use one of them and then, after a short lag phase, switch to the other. Until now it had been assumed that all cells in a population engage in the outgrowth on the second sugar after major metabolic adaptation of enzymatic composition has occurred, which takes time (hence the lag phase in growth). Here, we show that actually only a subpopulation is fit enough to partake in the second growth phase and present an evolutionary model, suggesting that this phenomenon might entail a bet-hedging strategy that helps bacteria adapt to the unexpectedly changing environment. When bacteria grow in a medium with two sugars, they first use the preferred sugar and only then start metabolizing the second one. After the first exponential growth phase, a short lag phase of nongrowth is observed, a period called the diauxie lag phase. It is commonly seen as a phase in which the bacteria prepare themselves to use the second sugar. Here we reveal that, in contrast to the established concept of metabolic adaptation in the lag phase, two stable cell types with alternative metabolic strategies emerge and coexist in a culture of the bacterium Lactococcus lactis. Only one of them continues to grow. The fraction of each metabolic phenotype depends on the level of catabolite repression and the metabolic state-dependent induction of stringent response, as well as on epigenetic cues. Furthermore, we show that the production of alternative metabolic phenotypes potentially entails a bet-hedging strategy. This study sheds new light on phenotypic heterogeneity during various lag phases occurring in microbiology and biotechnology and adjusts the generally accepted explanation of enzymatic adaptation proposed by Monod and shared by scientists for more than half a century.