- Research article
- Open Access
Metabolic modelling in a dynamic evolutionary framework predicts adaptive diversification of bacteria in a long-term evolution experiment
© The Author(s). 2016
Received: 21 April 2016
Accepted: 4 August 2016
Published: 20 August 2016
Predicting adaptive trajectories is a major goal of evolutionary biology and useful for practical applications. Systems biology has enabled the development of genome-scale metabolic models. However, analysing these models via flux balance analysis (FBA) cannot predict many evolutionary outcomes including adaptive diversification, whereby an ancestral lineage diverges to fill multiple niches. Here we combine in silico evolution with FBA and apply this modelling framework, evoFBA, to a long-term evolution experiment with Escherichia coli.
Simulations predicted the adaptive diversification that occurred in one experimental population and generated hypotheses about the mechanisms that promoted coexistence of the diverged lineages. We experimentally tested and, on balance, verified these mechanisms, showing that diversification involved niche construction and character displacement through differential nutrient uptake and altered metabolic regulation.
The evoFBA framework represents a promising new way to model biochemical evolution, one that can generate testable predictions about evolutionary and ecosystem-level outcomes.
The ability to predict evolution would be valuable not only for understanding such processes as adaptation and speciation [1–3], but also for engineering robust industrial strains, anticipating ecosystem responses to climate change, and combatting antibiotic resistance [4–7]. Models that capture the relationship between genotypes and environments, the structure and state of regulatory and metabolic networks, and the resulting phenotypes are likely to be important for developing these predictive abilities [1, 3, 8]. Ultimately, models of the relationship between genotype and phenotype will need to be combined with models of evolutionary and ecological dynamics in integrated frameworks that can predict the trajectory of evolution [5, 9].
The dynamics of evolutionary change reflect multiple processes and varying selective pressures that are influenced by many ecological, physical, and cellular constraints that may conflict with one another. Understanding whether and how these dynamics lead to the splitting and divergence of lineages is of central interest, as these processes represent the initial steps towards speciation. To this end, several theoretical studies have shown that cellular tradeoffs can promote lineage divergence [10–16]. The importance of such tradeoffs can be readily understood in the context of metabolism and growth. For example, if there were no tradeoffs, then one would predict that cells should maximize their expression of transporters and their surface area to achieve the highest possible rate of substrate uptake . However, such cellular investments would impinge on other cellular processes owing to competing requirements for membrane and cytosol space [18, 19], ribosomes [15, 16], and redox carriers [20, 21]. Thus, cells may appear suboptimal for individual physiological parameters, but this might be merely a consequence of being optimal for the combined set of parameters and associated cellular tradeoffs.
Historically, the interplay between cellular tradeoffs and evolutionary and ecological dynamics has been analyzed using game theory and differential equation-based models that consider small or idealized metabolic systems [10, 11, 14, 22]. These studies have highlighted that tradeoffs in cellular metabolism can lead to incomplete degradation of a resource, resulting in the evolution of cross-feeding interactions [10, 11]. This phenomenon has been seen in several evolution experiments under both batch and chemostat conditions [23–27]. To increase predictive power in microbial ecology and evolution, it is now desirable to develop models that can take into account cellular metabolism at a larger scale and across different organisms. Stoichiometric models offer a promising approach because, in principle, they can capture all enzyme-mediated metabolic reactions of an organism in an unbiased and non-supervised way using genomic information .
Flux Balance Analysis (FBA) has been developed to determine the optimal metabolic state of an organism, given knowledge of its biochemical network, biomass composition, and uptake flux rates . This approach is based on the assumptions that evolution has optimized metabolism and that metabolic fluxes can be predicted by setting the growth rate for a given rate of substrate uptake (such that the ratio of the two rates represents a yield) as an optimization criterion that can be solved by linear programming [28–30]. Early applications of FBA ignored the essential role of tradeoffs in the computation of metabolic fluxes [28, 31, 32], but more recent applications have incorporated tradeoffs as constraints on total fluxes [18, 19, 33, 34] and thereby achieved better prediction of experimentally observed metabolic states, such as preferential substrate utilization  and acetate overflow . Experimentally measured reaction thermodynamics and gene expression levels have also been used to constrain optimal metabolic states that reflect tradeoffs [35–37], and there have been efforts to combine FBA with ecological interactions between multiple species in microbial communities [38–45]. These approaches use species-specific models in a shared environment to maximize a predefined, community-level objective [39, 41, 43, 44] or apply FBA within a dynamic framework . The latter approach enables prediction of ecological interactions such as competition and cross-feeding between different species making up the model community, given defined substrate uptake constraints for each model species [40, 42, 45]. However, none of these approaches can currently be used to predict the interplay between ecological and evolutionary dynamics.
Here, we begin to overcome these limitations by integrating a FBA model of multi-phenotype systems with both cellular constraints and evolutionary dynamics. We define an overall constraint on uptake rates to enforce tradeoffs while simulating multiple model organisms living in the same environment without the need to specify each organism’s uptake preferences a priori (for details on how evolution and mutations are simulated see Methods section). By limiting total uptake in the model, and including O2 “uptake” in that total, we seek to represent cellular limitations that can arise from many diverse processes, including redox cycling [20, 47], respiratory chain , enzyme expression [16, 48], and substrate uptake . Although O2 uptake per se might not be limiting, limitations in the electron transport chain can effectively limit O2 respiration. Accounting for all the different possible limitations arising from cellular processes in a mechanistic manner is beyond the scope of stoichiometric models; however, limiting total uptake provides a general constraint that allows us to implement the tradeoffs observed in different studies in a simple, consistent, albeit approximate manner [16–18, 20, 47].
This approach allows integration of evolutionary dynamics by mutations that change substrate uptake rates along with the optimization of each model organism in the context of other model organisms that are present and coevolving in the same environment. The combined framework, which we call evoFBA, thus aims to provide a more realistic way to model the interplay between ecological and evolutionary dynamics with global constraints arising from cellular tradeoffs. To the best of our knowledge, this is the first FBA modeling approach that captures the continuous adaptation of organisms to the interplay between ecological and evolutionary dynamics in systems with multiple strains or species.
To examine the ability of evoFBA to capture ecological and evolutionary dynamics, we used it to simulate the evolution of Escherichia coli populations in a defined glucose-limited environment with daily transfers. We then experimentally analyzed the predictions of evoFBA in the context of the long-term evolution experiment (LTEE) with E. coli, in which 12 populations started from a common ancestor have been propagated in a glucose-limited medium for more than 60,000 generations [2, 49]. We found that the evoFBA simulations predicted the emergence of cross-feeding model organisms as a stable end-point, which in fact has occurred in at least one of the LTEE populations [26, 50]. Moreover, we saw that key metabolic features of the model organisms were in qualitative agreement with the physiological properties we measured for the two biological lineages that emerged and subsequently coexisted for more than 50,000 generations.
Microbial communities and their underlying metabolic interactions reflect the ecological and evolutionary histories of the component species . To capture these interactions, we combine stoichiometric metabolic models with ecological and evolutionary dynamics in the multi-layered evoFBA framework (see Methods). To test the utility of this framework, we apply it to the LTEE in which E. coli populations evolve in a defined glucose-limited environment [2, 52].
To model the LTEE, we ran evoFBA simulations starting with a metabolic model of E. coli that accounts for 14 carbon sources including glucose and byproducts that can be scavenged from the environment to produce biomass and fuel associated core metabolic reactions. In each evoFBA simulation, we allowed the metabolic model to change by random mutations under global constraints that must be obeyed. Thus, each simulation produced mutant model organisms exhibiting different uptake rates, metabolic flux patterns, and resulting growth rates.
evoFBA predicts evolution of cross-feeding between lineages with different metabolic flux distributions
Adaptive diversification in one LTEE population, matching evoFBA predictions
Two distinct lineages had emerged in one of the LTEE populations, called Ara-2, by 6500 generations, and they have coexisted ever since [26, 50]. The lineages are called S (small) and L (large) after their colony sizes on agar plates. The maintenance of this polymorphism depends on a cross-feeding interaction in which the L type is a better competitor for the exogenously supplied glucose and the S type is better at using one or more secreted byproducts , although the precise ecological and metabolic mechanisms are still unknown. Therefore, we used predictions from the evoFBA simulations to generate hypotheses about these mechanisms.
We hypothesized that, first, L specializes on glucose and secretes acetate and, second, S specializes by improved acetate consumption. We tested this hypothesis by analyzing two evolved clones sampled at generation 6500 from the S and L lineages, named 6.5KS1 and 6.5KL4, respectively. HPLC analyses confirmed the presence of acetate in a 24-h supernatant of 6.5KL4 that was grown in the same medium as the LTEE (see Methods). Acetate was not detected after growing 6.5KS1 in that supernatant (Additional file 1: Figure S1). We then measured the acetate and glucose concentrations over time in cultures of the ancestor, 6.5KS1, and 6.5KL4 clones in DM250-glucose medium (Fig. 3b). Both the L and S clones consumed glucose faster than the ancestor, consistent with previous assays . Moreover, in agreement with the evoFBA results, 6.5KL4 secreted acetate, with its concentration remaining high for many hours in the monoculture, and 6.5KS1 drew down its own acetate secretion much faster than both 6.5KL4 and the ancestor. After exhausting the glucose by 6 h, 6.5KS1 showed diauxic growth and consumed acetate until it was depleted after 9 h, whereas 6.5KL4 had barely, if at all, begun to consume acetate at that time even as it had exhausted the glucose by 5 h (Fig. 3b). These results support the hypothesis that the stable coexistence of S and L depends on acetate cross-feeding, with acetate production by both the L and S lineages and more efficient acetate scavenging by the S lineage, which exhibits a faster metabolic switch from glucose to acetate (Additional file 2: Figure S2).
Physiology and fluxes in S and L clones agree qualitatively with evoFBA
We developed a modeling framework, called evoFBA, which combines metabolic models that are amenable to FBA with an evolutionary algorithm to simulate the interplay of evolutionary and ecological dynamics in systems with multiple strains or species. We applied evoFBA to the LTEE with E. coli and predicted the emergence of two stably coexisting lineages with distinct metabolic flux distributions that promote a cross-feeding interaction. These predictions fit with the polymorphism seen in the Ara–2 population, where two lineages emerged early in the LTEE that have now coexisted for tens of thousands of generations [26, 50]. The evoFBA simulations enabled us to hypothesize specific ecological and physiological mechanisms that generate and sustain this polymorphism, and we then tested these hypotheses using the bacteria from that LTEE population. On balance, we found that the ecological, physiological, and metabolic properties of these coexisting lineages agree reasonably well with the predictions of the evoFBA.
Only one of the 12 LTEE populations evolved a persistent polymorphism that has been studied in such detail [26, 50]. However, other LTEE populations show evidence of negative frequency-dependent fitness, deep phylogenetic divergence, or both, which are consistent with adaptive diversification [57–60]. These results suggest that other populations may have evolved cross-feeding interactions similar to the one studied here, even if they were not always so persistent . One possible explanation for why persistent polymorphisms did not evolve in the other populations is that the establishment of the S lineage appears to have involved epistatic interactions between multiple mutations , which may have limited its evolutionary accessibility . This possibility reflects one of the limitations of evoFBA, which cannot capture all of the intricacies of biological evolution but instead predicts optimal end states that emerge given the simplifying assumptions of this approach. For example, mutations in evoFBA affect the rates of resource uptake, but not the rates of internal reactions in the model. This limitation reflects the computational burden of simulating a multitude of mutant genotypes, the number of which would increase greatly if all reaction rates were subject to mutation. This limitation could be relieved by the development of more efficient algorithms (allowing mutations to affect all reactions in the model), but the final evolved model organisms might not differ functionally from those based on the current approach because changes in uptake rates can already affect downstream flux distributions. Another limitation of the evoFBA approach at this time is the assumption that constraints on the uptake fluxes can be changed only by mutation, while the optimization of fluxes within those constraints is immediate through FBA . In other words, metabolic fluxes change within physiological limits without delay, whereas changing the limits themselves requires mutations. These assumptions are reasonable starting points for incorporating evolutionary dynamics into an FBA framework, but physiological delays in metabolic adjustments are also sometimes important [61, 62]. Expanding the evoFBA framework to include the dynamics of physiological transitions could start by integrating previous work on incorporating gene regulation into FBA [63, 64].
Adaptive diversification is expected, and has occurred, in other evolution experiments besides the LTEE, such as when two exogenous carbon sources are provided [27, 65] and in high-glucose chemostats, where substantial acetate is produced [21, 23, 25]. However, the adaptive diversification observed in the LTEE was unexpected owing to the presence of a single carbon source, glucose, which was supplied at a low concentration . Using the evoFBA framework, we predicted that acetate secretion was the primary metabolic driver leading to the emergence of the polymorphism, and that prediction was supported by our experiments. The long duration of the LTEE—including several thousand generations to establish the S and L polymorphism  and its persistence for tens of thousands of generations —may have facilitated adaptive divergence under these more restrictive ecological conditions, in comparison with other studies of much shorter duration where glucose, acetate, or both were supplied exogenously at high concentrations [23, 25, 45]. In fact, low levels of acetate, as observed in our study, have previously been reported to favor generalists as opposed to divergence into coexisting specialists . Nonetheless, the L lineage evolved higher glucose uptake rates, which led to acetate overflow and the construction of a new niche that benefited the S lineage, as occurred in the evoFBA simulations. Thus, niche construction by the bacteria led to the emergence of this polymorphism, in contrast to experiments where both carbon sources were added to the medium . Despite the differences between the LTEE and previous evolution experiments [45, 65], similar metabolic processes emerged.
The combination of evoFBA and experimental evolution provides a useful approach that can give insights into general mechanisms involved in the emergence of bacterial diversity and community construction. This approach may stimulate the development of even more detailed and integrated studies aimed at predicting the outcomes of evolution experiments and dynamics in multi-species systems including synthetic microbial communities .
Evolutionary flux balance analysis (evoFBA)
In each evoFBA simulation, stoichiometric metabolic models were used to simulate clonal populations with distinct genotypes. Each genotype was represented by a metabolic model, which was simulated in a dynamical FBA formalism  to evaluate its growth and metabolic flux rates over time. At each time step of the dynamic FBA, the metabolic model was optimized using linear optimization and a pseudo-reaction representing biomass as the objective function . This optimization thus maximized growth rate given the constraints on uptake rates, i.e. it optimized biomass yield per substrate . Instead of defining specific uptake rates for a particular set of media components (as in standard FBA studies), we assumed a global constraint for all carbon and oxygen uptake reactions in each model organism. By limiting total uptake in the model (including O2 “uptake”), we represent cellular limitations that can arise from many different factors, including redox cycling [20, 47], respiratory chain , enzyme expression [16, 48], and substrate uptake . Similar implementations of global constraints in FBA models have been employed previously to study diauxic shift and substrate preference in E. coli [18, 19, 37]. The global uptake constraint implemented here favored a minimization of fluxes given the maximization of the FBA objective in order to achieve the most efficient use of cellular resources for growth .
New model organisms were generated by random mutation from existing ones. Mutations altered specific bounds on individual uptake rates, while maintaining an overall total flux constraint of carbon and oxygen into the model organism. Thus, mutations change how the overall uptake flux is distributed across different substrates, and they allow a second level of optimization to occur over evolutionary time in addition to the optimization that occurs by FBA over the physiological time scale. Focusing evolution on a subset of reactions made computation of the ecological and evolutionary dynamics feasible; even so, the simulations presented here took over 20 days on a dedicated high-performance computer to simulate over 90,000 different model organism genotypes and their associated population and metabolite dynamics. A more complete simulation might encompass genome-scale models with evolution of all reactions in the model and with global constraints on total protein biomass  and membrane space . The current implementation of evoFBA was unable to perform such simulations in a reasonable timeframe and with appropriate numbers of replicate simulations; efforts to run evoFBA with mutations allowed for all reactions caused a slowdown of over 10-fold relative to the current implementation.
Representing mutations in evoFBA
Simulating serial transfer and selection
To simulate the LTEE’s daily transfer cycles , we used dynamical FBA to compute growth over 24 h; selection is a direct consequence of the differential growth of the model genotypes. After 24 h, a dilution was performed by randomly drawing 1 % of the model organisms, which constituted the initial population for the next simulated day. The next day’s medium included 99 % of the initial medium and 1 % of the spent medium from the end of the previous day. The simulated growth and dilution ran for a total of 550 cycles. Results from replicate simulations (Fig. 2) are qualitatively similar to those in Fig. 1. The population dynamics arising from these simulations are expected to give rise eventually to one dominant clone in each stably coexisting lineage. However, similar model organisms may occur within a simulation as a result of independent mutations before any one of them has reached its population maximum (e.g., model genotypes 13437, 12364 and 12719 in Fig. 1). However, if the model organisms differ even slightly in their uptake rates, then one genotype should eventually prevail through competitive exclusion, unless the model organisms occupy distinct ecological niches (Fig. 1).
Simulations were performed using MATLAB (Math Works, Natick, Massachusetts) and dynamic FBA calculations using the COBRA toolbox . The MATLAB scripts used to run evoFBA and analyze the data are freely available at .
LTEE and bacterial strains
The LTEE consists of 12 populations founded from the same ancestral strain of E. coli, REL606 , that have been propagated since 1988 by daily 1:100 dilutions in Davis minimal medium  supplemented with glucose at 25 mg/l (DM25). Here, we focused on one population, called Ara–2, in which two lineages, S and L, diverged before 6500 generations and have co-existed ever since [26, 50, 60]. We studied the ancestor and one clone sampled from each lineage at 6500, 11,000, 18,000, 20,000, 30,000, 40,000 and 50,000 generations. Each evolved clone is named by its generation followed by S or L according to its lineage and an arbitrary numeral for a given clone. For example, 6.5KS1 is a clone from the S lineage that was sampled at 6500 generations.
Media and culture conditions
Bacteria were grown in the same medium as used in the LTEE , except that the carbon source was glucose at 250 mg/l (DM250-glucose), glucose at 1000 mg/l (DM1000-glucose), or acetate at 250 mg/l (DM250-acetate). These higher concentrations were used to increase cell density and thereby improve the accuracy of measurements of cell growth (e.g., Fig. 5) and concentrations of residual resources and secreted metabolites (e.g., Additional file 1: Figure S1). After overnight growth in DM1000-glucose, strains were inoculated by a 10,000-fold dilution into DM250-glucose, where they grew for 24 h at 37 °C with shaking at 120 rpm as an acclimation step. For each strain, three replicate acclimation cultures were then inoculated as duplicates, each at a 1:100 dilution, into DM250-glucose or DM250-acetate and incubated in 96-well microtiter plates at 37 °C for 24 h. Growth was monitored using an Infinite M200 microplate reader (Tecan, Lyon, France) by measuring the OD450nm every 10 min. Growth rates were computed from filtered OD data as dln(OD450)/dt over a sliding window of 1 h, using MATLAB. We report the mean of the three replicates. Filtering was performed by removing negative growth rates and mean growth rates that were more than 0.2 units above or below the immediately adjacent data points (outliers).
Measuring glucose and acetate concentrations
The ancestor, 6.5KS1, and 6.5KL4 clones were grown in DM250-glucose as described. Samples were taken at time 0 and every h for 9 h. After centrifugation to remove cells, we measured glucose and acetate concentrations in the supernatant using the Glucose Assay Kit (Merck Millipore, Lyon, France) and Acetic Acid Assay Kit (Megazyme, Pontcharra-sur-Turdine, France), respectively, following the manufacturers’ recommendations.
Analysis of flux patterns in individual model organisms
We simulated the growth of the evoFBA model organisms with IDs 44490 and 12364 (Fig. 1) to obtain the flux values for their biochemical reactions. Each model organism was simulated using dynamical FBA in medium containing 0.1389 mM glucose, the same concentration as in the LTEE. Each simulation ran for ten 24-h periods with daily 1:100 dilutions; the last day was used to record the flux values, in order to remove any effect of the initial conditions. The flux patterns for growth on glucose were taken 10 min after the onset of growth (Fig. 4a and c), and for growth on acetate at 388 min because glucose was exhausted while acetate was still present at a substantial level (Fig. 4b and d). From the flux patterns, we identified several reactions of interest that showed differences between the two evolved model organisms (highlighted in blue in Fig. 4).
Identification of metabolites in filtrates of spent cultures of 6.5KL4
We analyzed by HPLC and GC-MS the metabolic by-products secreted by clone 6.5KL4 using filtrates from 24-h spent cultures of that clone in DM25- and DM250-glucose, both before and after growth of clone 6.5KS1. For HPLC, 1 ml of filtrate was acidified with 5 μl 1 M H2SO4, incubated at room temperature for 5 min, and passed through a 0.45-μm regenerated cellulose syringe filter (PHENEX RC Membrane, Phenomenex, Le Pecq, France). Samples were then analyzed on an Agilent 1260 Infinity HPLC system equipped with a Rezex ROA-Organic Acid (8 %) 300 × 7.8-mm column (Phenomenex) and a diode array detector. The analytical conditions were as follows: mobile phase, 5 mM H2SO4; flow rate, 0.6 ml/min; column temperature, 35 °C; injection volume, 50 μl; wavelength scan range, 190–400 nm; detection wavelength, 210 nm; and run time, 35 min. Concentrations of acetate and fumarate in the L-clone filtrates were determined from linear standard curves over the ranges of 0–10 mM and 0–100 μM, respectively, and with lower detection limits of 0.1 mM and 0.3 μM, respectively. Succinate, lactate, formate, propionate, and butyrate can also be separated under these analytical conditions with detection limits similar to acetate, but they were not detected in any samples.
GC-MS analysis of volatile compounds was performed using an Agilent GC HP6890 gas chromatograph equipped with a Varian CP-WAX 58 column (length, 25 m; internal diameter, 0.25 mm; film thickness, 0.20 μm), and coupled to an MSD5973 mass sensitive detector. The sample (600 μl) was cooled on ice, acidified with 50 μl 4 M HCl, and extracted with 0.375 g NaCl and 650 μl ice-cold ether. After vortexing three times for 10 s each, with 30 s cooling intervals, the sample was centrifuged for 5 min at 10,000 rpm and placed on ice for 5 min. The upper organic layer (2.5 μl) was then injected manually into the GC using an ice-cold syringe (injection in split mode, split ratio = 10). The column was held at 40 °C for 1 min, ramped to 200 °C at a rate of 5 °C/min, and held for a further 3 min, giving a total run time of 36 min. The solvent delay for the MSD was 1.4 min and the mass scan range was set to 35–300 atomic mass units. The presence of acetate (retention time 12.2 min) in the L filtrate after growth in both DM25- and DM250-glucose was confirmed by this method. The estimated concentration from the DM250-glucose filtrate was 510 μM, which is close to the 480 μM detected by HPLC (Additional file 1: Figure S1). The concentration of acetate in the L filtrate from DM25-glucose was about one-tenth that detected in DM250-glucose. Ethanol (retention time 2.1 min) was also detected in the filtrates of all three strains tested (ancestor, 6.5KL4, and 6.5KS1). Other metabolites including isopropanol, butanol, acetoin, acetone, formic acid, propionic acid, butyric acid, isobutyric acid, and valeric acid were not detected (with lower detection limits around 50–100 μM in scan mode).
Analysis of promoter activities in LTEE clones
Genes used in the analysis of promoter activities
FBA model term
Malate synthase A
JC acknowledges the French ministry of research for a research fellowship. This work was supported by grants from the Biotechnology and Biological Sciences Research Council (BBSRC, BB/K003240/1 to OSS), Agence Nationale de la Recherche Programme Blanc (ANR-08-BLAN-0283-01 to DS), European Union programme FP7-ICT-2013-10 project EvoEvo (610427 to DS), Université Grenoble Alpes (to DS), Centre National de la Recherche Scientifique (CNRS, to DS), the U.S. National Science Foundation (DEB-1451740 to REL), and the BEACON Center for the Study of Evolution in Action (NSF Cooperative Agreement DBI-0939454 to REL).
JC, TG, JG performed the experiments, analysed and interpreted the data, and drafted the article; JW performed the HPLC experiments and drafted the article; JC, TG, DS, OSS designed the study, analysed and interpreted the data, drafted and revised the article; REL directs the LTEE and helped draft and revise the article. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Stern DL, Orgogozo V, Pierre U, Bâtiment A, Saint Q. Is genetic evolution predictable? Nonrandom distribution of evolutionarily relevant mutations. Science. 2009;323:746–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Wiser MJ, Ribeck N, Lenski RE. Long-term dynamics of adaptation in asexual populations. Science. 2013;342:1364–7.View ArticlePubMedGoogle Scholar
- de Visser JAGM, Krug J. Empirical fitness landscapes and the predictability of evolution. Nat Rev Genet. 2014;15:480–90.View ArticlePubMedGoogle Scholar
- Nesse RM, Stearns SC. The great opportunity: Evolutionary applications to medicine and public health. Evol Appl. 2008;1:28–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Soyer OS, O’Malley MA. Evolutionary systems biology: What it is and why it matters. BioEssays. 2013;35:696–705.View ArticlePubMedGoogle Scholar
- Lohbeck KT, Riebesell U, Reusch TBH. Adaptive evolution of a key phytoplankton species to ocean acidification. Nat Geosci. 2012;5:346–51.View ArticleGoogle Scholar
- Novais Â, Comas I, Baquero F, Cantón R, Coque TM, Moya A, et al. Evolutionary trajectories of beta-lactamase CTX-M-1 cluster enzymes: Predicting antibiotic resistance. PLoS Pathog. 2010;6:e1000735.View ArticlePubMedPubMed CentralGoogle Scholar
- Papp B, Notebaart RA, Pál C. Systems-biology approaches for predicting genomic evolution. Nat Rev Genet. 2011;12:591–602.View ArticlePubMedGoogle Scholar
- Prosser JI, Bohannan BJM, Curtis TP, Ellis RJ, Firestone MK, Freckleton RP, et al. The role of ecological theory in microbial ecology. Nat Rev Microbiol. 2007;5:384–92.View ArticlePubMedGoogle Scholar
- Doebeli M. A model for the evolutionary dynamics of cross-feeding polymorphisms in microorganisms. Popul Ecol. 2002;44:59–70.View ArticleGoogle Scholar
- Pfeiffer T, Bonhoeffer S. Evolution of cross-feeding in microbial populations. Am Nat. 2004;163:E126–35.View ArticlePubMedGoogle Scholar
- Gudelj I, Beardmore RE, Arkin SS, MacLean RC. Constraints on microbial metabolism drive evolutionary diversification in homogeneous environments. J Evol Biol. 2007;20:1882–9.View ArticlePubMedGoogle Scholar
- Schuster S, Pfeiffer T, Fell DA. Is maximization of molar yield in metabolic networks favoured by evolution? J Theor Biol. 2008;252:497–504.View ArticlePubMedGoogle Scholar
- Pfeiffer T, Schuster S, Bonhoeffer S. Cooperation and competition in the evolution of ATP-producing pathways. Science. 2001;292:504–7.View ArticlePubMedGoogle Scholar
- Molenaar D, van Berlo R, de Ridder D, Teusink B. Shifts in growth strategies reflect tradeoffs in cellular economics. Mol Syst Biol. 2009;5:323.View ArticlePubMedPubMed CentralGoogle Scholar
- Weiße AY, Oyarzún DA, Danos V, Swain PS. Mechanistic links between cellular trade-offs, gene expression, and growth. Proc Natl Acad Sci U S A. 2015;112:E1038–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Button D. Nutrient uptake by microorganisms according to kinetic parameters from theory as related to cytoarchitecture. Microbiol Mol Biol Rev. 1998;62:636–45.PubMedPubMed CentralGoogle Scholar
- Zhuang K, Vemuri GN, Mahadevan R. Economics of membrane occupancy and respiro-fermentation. Mol Syst Biol. 2011;7:500.View ArticlePubMedPubMed CentralGoogle Scholar
- Beg QK, Vazquez A, Ernst J, de Menezes MA, Bar-Joseph Z, Barabási A-L-L, et al. Intracellular crowding defines the mode and sequence of substrate uptake by Escherichia coli and constrains its metabolic activity. Proc Natl Acad Sci U S A. 2007;104:12663–8.View ArticlePubMedPubMed CentralGoogle Scholar
- van Hoek MJ, Merks RMH. Redox balance is key to explaining full vs. partial switching to low-yield metabolism. BMC Syst Biol. 2012;6:22.View ArticlePubMedPubMed CentralGoogle Scholar
- Vemuri GN, Altman E, Sangurdekar DP, Khodursky AB, Eiteman MA. Overflow metabolism in Escherichia coli during steady-state growth: Transcriptional regulation and effect of the redox ratio. Appl Environ Microbiol. 2006;72:3653–61.View ArticlePubMed CentralGoogle Scholar
- Schuster S, de Figueiredo LF, Schroeter A, Kaleta C. Combining metabolic pathway analysis with evolutionary game theory. Explaining the occurrence of low-yield pathways by an analytic optimization approach. Biosystems. 2011;105:147–53.View ArticlePubMedGoogle Scholar
- Rosenzweig RF, Sharp RR, Treves DS, Adams J. Microbial evolution in a simple unstructured environment: Genetic differentiation in Escherichia coli. Genetics. 1994;137:903–17.PubMedPubMed CentralGoogle Scholar
- Turner PE, Souza V, Lenski RE. Tests of ecological mechanisms promoting the stable coexistence of two bacterial genotypes. Ecology. 1996;77:2119–29.View ArticleGoogle Scholar
- Treves DS, Manning S, Adams J. Repeated evolution of an acetate-crossfeeding polymorphism in long-term populations of Escherichia coli. Mol Biol Evol. 1998;15:789–97.View ArticlePubMedGoogle Scholar
- Rozen DE, Lenski RE. Long-term experimental evolution in Escherichia coli. VIII. Dynamics of a balanced polymorphism. Am Nat. 2000;155:24–35.View ArticlePubMedGoogle Scholar
- Friesen ML, Saxer G, Travisano M, Doebeli M. Experimental evidence for sympatric ecological diversification due to frequency-dependent competition in Escherichia coli. Evolution. 2004;58:245–60.View ArticlePubMedGoogle Scholar
- Varma A, Palsson BO. Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol. 1994;60:3724–31.PubMedPubMed CentralGoogle Scholar
- Orth J, Thiele I, Palsson B. What is flux balance analysis? Nat Biotechnol. 2010;28:245–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Kauffman KJ, Prakash P, Edwards JS. Advances in flux balance analysis. Curr Opin Biotechnol. 2003;14:491–6.View ArticlePubMedGoogle Scholar
- Feist AM, Henry CS, Reed JL, Krummenacker M, Joyce AR, Karp PD, et al. A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol Syst Biol. 2007;3:121.View ArticlePubMedPubMed CentralGoogle Scholar
- Ibarra RU, Edwards JS, Palsson BO. Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth. Nature. 2002;420:20–3.View ArticleGoogle Scholar
- Heinrich R, Schuster S, Holzhütter H-GG. Mathematical analysis of enzymic reaction systems using optimization principles. Eur J Biochem. 1991;201:1–21.View ArticlePubMedGoogle Scholar
- Holzhütter HG. The principle of flux minimization and its application to estimate stationary fluxes in metabolic networks. Eur J Biochem. 2004;271:2905–22.View ArticlePubMedGoogle Scholar
- Hoppe A, Hoffmann S, Holzhütter H-G. Including metabolite concentrations into flux balance analysis: thermodynamic realizability as a constraint on flux distributions in metabolic networks. BMC Syst Biol. 2007;1:23.View ArticlePubMedPubMed CentralGoogle Scholar
- Thiele I, Jamshidi N, Fleming RMT, Palsson BØ. Genome-scale reconstruction of Escherichia coli’s transcriptional and translational machinery: A knowledge base, its mathematical formulation, and its functional characterization. PLoS Comput Biol. 2009;5:e1000312.View ArticlePubMedPubMed CentralGoogle Scholar
- OBrien EJ, Lerman JA, Chang RL, Hyduke DR, Palsson BO. Genome-scale models of metabolism and gene expression extend and refine growth phenotype prediction. Mol Syst Biol. 2013;9:693.Google Scholar
- Klitgord N, Segrè D. Ecosystems biology of microbial metabolism. Curr Opin Biotechnol. 2011;22:541–6.View ArticlePubMedGoogle Scholar
- Zomorrodi AR, Maranas CD. OptCom: a multi-level optimization framework for the metabolic modeling and analysis of microbial communities. PLoS Comput Biol. 2012;8:e1002363.View ArticlePubMedPubMed CentralGoogle Scholar
- Harcombe WR, Riehl WJ, Dukovski I, Granger BR, Betts A, Lang AH, et al. Metabolic resource allocation in individual microbes determines ecosystem interactions and spatial dynamics. Cell Rep. 2014;7:1104–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Klitgord N, Segrè D. Environments that induce synthetic microbial ecosystems. PLoS Comput Biol. 2010;6:e1001002.View ArticlePubMedPubMed CentralGoogle Scholar
- Freilich S, Zarecki R, Eilam O, Segal ES, Henry CS, Kupiec M, et al. Competitive and cooperative metabolic interactions in bacterial communities. Nat Commun. 2011;2:589.View ArticlePubMedGoogle Scholar
- Stolyar S, Van Dien S, Hillesland KL, Pinel N, Lie TJ, Leigh JA, et al. Metabolic modeling of a mutualistic microbial community. Mol Syst Biol. 2007;3:92.View ArticlePubMedPubMed CentralGoogle Scholar
- Khandelwal RA, Olivier BG, Röling WFM, Teusink B, Bruggeman FJ. Community flux balance analysis for microbial consortia at balanced growth. PLoS One. 2013;8:e64567.View ArticlePubMedPubMed CentralGoogle Scholar
- Louca S, Doebeli M. Calibration and analysis of genome-based models for microbial ecology. Elife. 2015;4:1–17.View ArticleGoogle Scholar
- Mahadevan R, Edwards JS, Doyle FJ. Dynamic flux balance analysis of diauxic growth in Escherichia coli. Biophys J. 2002;83:1331–40.View ArticlePubMedPubMed CentralGoogle Scholar
- el-Mansi EM, Holms WH. Control of carbon flux to acetate excretion during growth of Escherichia coli in batch and continuous cultures. Microbiology. 1989;135:2875–83.View ArticleGoogle Scholar
- Basan M, Hui S, Okano H, Zhang Z, Shen Y, Williamson JR, et al. Overflow metabolism in Escherichia coli results from efficient proteome allocation. Nature. 2015;528:99–104.View ArticlePubMedPubMed CentralGoogle Scholar
- Lenski RE, Wiser MJ, Ribeck N, Blount ZD, Maddamsetti R, Burmeister AR, et al. Sustained fitness gains and variability in fitness trajectories in the long-term evolution experiment with Escherichia coli. Proc R Soc B. 2015;282:20152292.View ArticlePubMedPubMed CentralGoogle Scholar
- Plucain J, Hindré T, Le Gac M, Tenaillon O, Cruveiller S, Médigue C, et al. Epistasis and allele specificity in the emergence of a stable polymorphism in Escherichia coli. Science. 2014;343:1366–9.View ArticlePubMedGoogle Scholar
- Großkopf T, Soyer OS. Synthetic microbial communities. Curr Opin Microbiol. 2014;18:72–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Lenski RE, Rose MR, Simpson SC, Tadler SC. Long-term experimental evolution in Escherichia coli. I. Adaptation and divergence during 2,000 generations. Am Nat. 1991;138:1315–41.View ArticleGoogle Scholar
- Leiby N, Marx CJ. Metabolic erosion primarily through mutation accumulation, and not tradeoffs, drives limited evolution of substrate specificity in Escherichia coli. PLoS Biol. 2014;12:e1001789.View ArticlePubMedPubMed CentralGoogle Scholar
- Harcombe WR, Delaney NF, Leiby N, Klitgord N, Marx CJ. The ability of flux balance analysis to predict evolution of central metabolism scales with the initial distance to the optimum. PLoS Comput Biol. 2013;9:e1003091. Stelling J, editor.View ArticlePubMedPubMed CentralGoogle Scholar
- Lewis NE, Hixson KK, Conrad TM, Lerman JA, Charusanti P, Polpitiya AD, et al. Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models. Mol Syst Biol. 2010;6:390.View ArticlePubMedPubMed CentralGoogle Scholar
- Rossell S, van der Weijden CC, Lindenbergh A, van Tuijl A, Francke C, Bakker BM, et al. Unraveling the complexity of flux regulation: a new method demonstrated for nutrient starvation in Saccharomyces cerevisiae. Proc Natl Acad Sci U S A. 2006;103:2166–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Elena SF, Lenski RE. Long-term experimental evolution in Escherichia coli: VII. Mechanisms maintaining genetic variability within populations. Evolution. 1997;51:1059–67.Google Scholar
- Blount ZD, Borland CZ, Lenski RE. Historical contingency and the evolution of a key innovation in an experimental population of Escherichia coli. Proc Natl Acad Sci U S A. 2008;105:7899–906.View ArticlePubMedPubMed CentralGoogle Scholar
- Maddamsetti R, Lenski RE, Barrick JE. Adaptation, clonal interference, and frequency dependent interactions in a long term evolution experiment with Escherichia coli. Genetics. 2015;200:619–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Tenaillon O, Barrick JE, Ribeck N, Deatherage DE, Blanchard JL, Dasgupta A, et al. Tempo and mode of genome evolution in a 50,000 - generation experiment. Nature. 2016; in press.Google Scholar
- Lenski RE, Souza V, Duong LP, Phan QC, Nguyen TNM, Bertrand KP. Epistatic effects of promoter and repressor functions of the Tn10 tetracycline-resistance operon on the fitness of Escherichia coli. Mol Ecol. 1994;3:127–35.View ArticlePubMedGoogle Scholar
- Le Gac M, Brazas MD, Bertrand M, Tyerman JG, Spencer CC, Hancock REW, et al. Metabolic changes associated with adaptive diversification in Escherichia coli. Genetics. 2008;178:1049–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Shlomi T, Eisenberg Y, Sharan R, Ruppin E. A genome-scale computational study of the interplay between transcriptional regulation and metabolism. Mol Syst Biol. 2007;3:101.View ArticlePubMedPubMed CentralGoogle Scholar
- Covert MW, Palsson B. Transcriptional regulation in constraints-based metabolic models of Escherichia coli. J Biol Chem. 2002;277:28058–64.View ArticlePubMedGoogle Scholar
- Spencer CC, Bertrand M, Travisano M, Doebeli M. Adaptive diversification in genes that regulate resource use in Escherichia coli. PLoS Genet. 2007;3:0083–8.View ArticleGoogle Scholar
- Feist A, Palsson B. The biomass objective function. Curr Opin Microbiol. 2010;13:344–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Edwards JS, Ibarra RU, Palsson BO. In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data. Nat Biotechnol. 2001;19:125–30.View ArticlePubMedGoogle Scholar
- Orth JD, Fleming RM, Palsson BØ (2010a) 10.2.1 -- Reconstruction and use of microbial metabolic networks: the core Escherichia coli metabolic model as an educational guide. In EcoSal -- Escherichia coli and Salmonella Cellular and Molecular Biology, Karp PD (ed), 10.2.1. Washington DC: ASM Press.Google Scholar
- Wielgoss S, Barrick JE, Tenaillon O, Cruveiller S, Chane-Woon-Ming B, Medigue C, et al. Mutation rate inferred from synonymous substitutions in a long-term evolution experiment with escherichia coli. G3. 2011;1:183–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Loferer-Krößbacher M, Klima J, Psenner R. Determination of bacterial cell dry mass by transmission electron microscopy and densitmetric image analysis. Appl Environ Microbiol. 1998;64:688–94.PubMedPubMed CentralGoogle Scholar
- Becker SA, Feist AM, Mo ML, Hannum G, Palsson BØ, Herrgard MJ. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protoc. 2007;2:727–38.View ArticlePubMedGoogle Scholar
- MATLAB scripts used to run evoFBA and analyze data. http://osslab.lifesci.warwick.ac.uk/?pid=resources.
- Jeong H, Barbe V, Lee CH, Vallenet D, Yu DS, Choi SH, et al. Genome sequences of Escherichia coli B strains REL606 and BL21(DE3). J Mol Biol. 2009;394:644–52.View ArticlePubMedGoogle Scholar
- Zaslaver A, Bren A, Ronen M, Itzkovitz S, Kikoin I, Shavit S, et al. A comprehensive library of fluorescent transcriptional reporters for Escherichia coli. Nat Methods. 2006;3:623–8.View ArticlePubMedGoogle Scholar