Synthetic microbial communities (SynComs) of the human gut: design, assembly, and applications.

The human gut harbors native microbial communities, forming a highly complex ecosystem. Synthetic microbial communities (SynComs) of the human gut are an assembly of microorganisms isolated from human mucosa or fecal samples. In recent decades, the ever-expanding culturing capacity and affordable sequencing, together with advanced computational modeling, started a ''golden age'' for harnessing the beneficial potential of SynComs to fight gastrointestinal disorders, such as infections and chronic inflammatory bowel diseases. As simplified and completely defined microbiota, SynComs offer a promising reductionist approach to understanding the multispecies and multikingdom interactions in the microbe-host-immune axis. However, there are still many challenges to overcome before we can precisely construct SynComs of designed function and efficacy that allow the translation of scientific findings to patients' treatments. Here, we discussed the strategies used to design, assemble, and test a SynCom, and address the significant challenges, which are of microbiological, engineering, and translational nature, that stand in the way of using SynComs as live bacterial therapeutics.

Evolutionary coexistence in a fluctuating environment by specialization on resource level.

Microbial communities in fluctuating environments, such as oceans or the human gut, contain a wealth of diversity. This diversity contributes to the stability of communities and the functions they have in their hosts and ecosystems. To improve stability and increase production of beneficial compounds, we need to understand the underlying mechanisms causing this diversity. When nutrient levels fluctuate over time, one possibly relevant mechanism is coexistence between specialists on low and specialists on high nutrient levels. The relevance of this process is supported by the observations of coexistence in the laboratory, and by simple models, which show that negative frequency dependence of two such specialists can stabilize coexistence. However, as microbial populations are often large and fast growing, they evolve rapidly. Our aim is to determine what happens when species can evolve; whether evolutionary branching can create diversity or whether evolution will destabilize coexistence. We derive an analytical expression of the invasion fitness in fluctuating environments and use adaptive dynamics techniques to find that evolutionarily stable coexistence requires a special type of trade-off between growth at low and high nutrients. We do not find support for the necessary evolutionary trade-off in data available for the bacterium Escherichia coli and the yeast Saccharomyces cerevisiae on glucose. However, this type of data is scarce and might exist for other species or in different conditions. Moreover, we do find evidence for evolutionarily stable coexistence of the two species together. Since we find this coexistence in the scarce data that are available, we predict that specialization on resource level is a relevant mechanism for species diversity in microbial communities in fluctuating environments in natural settings.

Towards evolutionary predictions: Current promises and challenges.

Evolution has traditionally been a historical and descriptive science, and predicting future evolutionary processes has long been considered impossible. However, evolutionary predictions are increasingly being developed and used in medicine, agriculture, biotechnology and conservation biology. Evolutionary predictions may be used for different purposes, such as to prepare for the future, to try and change the course of evolution or to determine how well we understand evolutionary processes. Similarly, the exact aspect of the evolved population that we want to predict may also differ. For example, we could try to predict which genotype will dominate, the fitness of the population or the extinction probability of a population. In addition, there are many uses of evolutionary predictions that may not always be recognized as such. The main goal of this review is to increase awareness of methods and data in different research fields by showing the breadth of situations in which evolutionary predictions are made. We describe how diverse evolutionary predictions share a common structure described by the predictive scope, time scale and precision. Then, by using examples ranging from SARS-CoV2 and influenza to CRISPR-based gene drives and sustainable product formation in biotechnology, we discuss the methods for predicting evolution, the factors that affect predictability and how predictions can be used to prevent evolution in undesirable directions or to promote beneficial evolution (i.e. evolutionary control). We hope that this review will stimulate collaboration between fields by establishing a common language for evolutionary predictions.

Population bottleneck has only marginal effect on fitness evolution and its repeatability in dioecious Caenorhabditis elegans.

The predictability of evolution is expected to depend on the relative contribution of deterministic and stochastic processes. This ratio is modulated by effective population size. Smaller effective populations harbor less genetic diversity and stochastic processes are generally expected to play a larger role, leading to less repeatable evolutionary trajectories. Empirical insight into the relationship between effective population size and repeatability is limited and focused mostly on asexual organisms. Here, we tested whether fitness evolution was less repeatable after a population bottleneck in obligately outcrossing populations of Caenorhabditis elegans. Replicated populations founded by 500, 50, or five individuals (no/moderate/strong bottleneck) were exposed to a novel environment with a different bacterial prey. As a proxy for fitness, population size was measured after one week of growth before and after 15 weeks of evolution. Surprisingly, we found no significant differences among treatments in their fitness evolution. Even though the strong bottleneck reduced the relative contribution of selection to fitness variation, this did not translate to a significant reduction in the repeatability of fitness evolution. Thus, although a bottleneck reduced the contribution of deterministic processes, we conclude that the predictability of evolution may not universally depend on effective population size, especially in sexual organisms.

Continual evolution through coupled fast and slow feedbacks.

Continual evolution describes the unceasing evolution of at least one trait involving at least one organism. The Red Queen Hypothesis is a specific case in which continual evolution results from coevolution of at least two species. While microevolutionary studies have described examples in which evolution does not cease, understanding which general conditions lead to continual evolution or to stasis remains a major challenge. In many cases, it is unclear which experimental features or model assumptions are necessary for the observed continual evolution to emerge, and whether the described behavior is robust to variations in the given setup. Here, we aim to find the minimal set of conditions under which continual evolution occurs. To this end, we present a theoretical framework that does not assume any specific functional form and, therefore, can be applied to a wide variety of systems. Our framework is also general enough to make predictions about both monomorphic and polymorphic populations. We show that the combination of a fast positive and a slow negative feedback between environment, population, and evolving traits causes continual evolution to emerge even from the evolution of a single evolving trait, provided that the ecological timescale is sufficiently faster than the timescales of mutation and the negative feedback. Our approach and results thus contribute to a deeper understanding of the evolutionary dynamics resulting from biotic interactions.

Low affinity uniporter carrier proteins can increase net substrate uptake rate by reducing efflux.

Many organisms have several similar transporters with different affinities for the same substrate. Typically, high-affinity transporters are expressed when substrate is scarce and low-affinity ones when it is abundant. The benefit of using low instead of high-affinity transporters remains unclear, especially when additional nutrient sensors are present. Here, we investigate two hypotheses. It was previously hypothesized that there is a trade-off between the affinity and the catalytic efficiency of transporters, and we find some but no definitive support for it. Additionally, we propose that for uptake by facilitated diffusion, at saturating substrate concentrations, lowering the affinity enhances the net uptake rate by reducing substrate efflux. As a consequence, there exists an optimal, external-substrate-concentration dependent transporter affinity. A computational model of Saccharomyces cerevisiae glycolysis shows that using the low affinity HXT3 transporter instead of the high affinity HXT6 enhances the steady-state flux by 36%. We tried to test this hypothesis with yeast strains expressing a single glucose transporter modified to have either a high or a low affinity. However, due to the intimate link between glucose perception and metabolism, direct experimental proof for this hypothesis remained inconclusive. Still, our theoretical results provide a novel reason for the presence of low-affinity transport systems.

Metabolic enzyme cost explains variable trade-offs between microbial growth rate and yield.

Microbes may maximize the number of daughter cells per time or per amount of nutrients consumed. These two strategies correspond, respectively, to the use of enzyme-efficient or substrate-efficient metabolic pathways. In reality, fast growth is often associated with wasteful, yield-inefficient metabolism, and a general thermodynamic trade-off between growth rate and biomass yield has been proposed to explain this. We studied growth rate/yield trade-offs by using a novel modeling framework, Enzyme-Flux Cost Minimization (EFCM) and by assuming that the growth rate depends directly on the enzyme investment per rate of biomass production. In a comprehensive mathematical model of core metabolism in E. coli, we screened all elementary flux modes leading to cell synthesis, characterized them by the growth rates and yields they provide, and studied the shape of the resulting rate/yield Pareto front. By varying the model parameters, we found that the rate/yield trade-off is not universal, but depends on metabolic kinetics and environmental conditions. A prominent trade-off emerges under oxygen-limited growth, where yield-inefficient pathways support a 2-to-3 times higher growth rate than yield-efficient pathways. EFCM can be widely used to predict optimal metabolic states and growth rates under varying nutrient levels, perturbations of enzyme parameters, and single or multiple gene knockouts.

Eco-evolutionary Red Queen dynamics regulate biodiversity in a metabolite-driven microbial system.

The Red Queen Hypothesis proposes that perpetual co-evolution among organisms can result from purely biotic drivers. After more than four decades, there is no satisfactory understanding as to which mechanisms trigger Red Queen dynamics or their implications for ecosystem features such as biodiversity. One reason for such a knowledge gap is that typical models are complicated theories where limit cycles represent an idealized Red Queen, and therefore cannot be used to devise experimental setups. Here, we bridge this gap by introducing a simple model for microbial systems able to show Red Queen dynamics. We explore diverse biotic sources that can drive the emergence of the Red Queen and that have the potential to be found in nature or to be replicated in the laboratory. Our model enables an analytical understanding of how Red Queen dynamics emerge in our setup, and the translation of model terms and phenomenology into general underlying mechanisms. We observe, for example, that in our system the Red Queen offers opportunities for the increase of biodiversity by facilitating challenging conditions for intraspecific dominance, whereas stasis tends to homogenize the system. Our results can be used to design and engineer experimental microbial systems showing Red Queen dynamics.

Evolutionary pressures on microbial metabolic strategies in the chemostat.

Protein expression is shaped by evolutionary processes that tune microbial fitness. The limited biosynthetic capacity of a cell constrains protein expression and forces the cell to carefully manage its protein economy. In a chemostat, the physiology of the cell feeds back on the growth conditions, hindering intuitive understanding of how changes in protein concentration affect fitness. Here, we aim to provide a theoretical framework that addresses the selective pressures and optimal evolutionary-strategies in the chemostat. We show that the optimal enzyme levels are the result of a trade-off between the cost of their production and the benefit of their catalytic function. We also show that deviations from optimal enzyme levels are directly related to selection coefficients. The maximal fitness strategy for an organism in the chemostat is to express a well-defined metabolic subsystem known as an elementary flux mode. Using a coarse-grained, kinetic model of Saccharomyces cerevisiae's metabolism and growth, we illustrate that the dynamics and outcome of evolution in a chemostat can be very counter-intuitive: Strictly-respiring and strictly-fermenting strains can evolve from a common ancestor. This work provides a theoretical framework that relates a kinetic, mechanistic view on metabolism with cellular physiology and evolutionary dynamics in the chemostat.

Interplay between constraints, objectives, and optimality for genome-scale stoichiometric models.

High-throughput data generation and genome-scale stoichiometric models have greatly facilitated the comprehensive study of metabolic networks. The computation of all feasible metabolic routes with these models, given stoichiometric, thermodynamic, and steady-state constraints, provides important insights into the metabolic capacities of a cell. How the feasible metabolic routes emerge from the interplay between flux constraints, optimality objectives, and the entire metabolic network of a cell is, however, only partially understood. We show how optimal metabolic routes, resulting from flux balance analysis computations, arise out of elementary flux modes, constraints, and optimization objectives. We illustrate our findings with a genome-scale stoichiometric model of Escherichia coli metabolism. In the case of one flux constraint, all feasible optimal flux routes can be derived from elementary flux modes alone. We found up to 120 million of such optimal elementary flux modes. We introduce a new computational method to compute the corner points of the optimal solution space fast and efficiently. Optimal flux routes no longer depend exclusively on elementary flux modes when we impose additional constraints; new optimal metabolic routes arise out of combinations of elementary flux modes. The solution space of feasible metabolic routes shrinks enormously when additional objectives---e.g. those related to pathway expression costs or pathway length---are introduced. In many cases, only a single metabolic route remains that is both feasible and optimal. This paper contributes to reaching a complete topological understanding of the metabolic capacity of organisms in terms of metabolic flux routes, one that is most natural to biochemists and biotechnologists studying and engineering metabolism.

Metabolic states with maximal specific rate carry flux through an elementary flux mode.

Specific product formation rates and cellular growth rates are important maximization targets in biotechnology and microbial evolution. Maximization of a specific rate (i.e. a rate expressed per unit biomass amount) requires the expression of particular metabolic pathways at optimal enzyme concentrations. In contrast to the prediction of maximal product yields, any prediction of optimal specific rates at the genome scale is currently computationally intractable, even if the kinetic properties of all enzymes are available. In the present study, we characterize maximal-specific-rate states of metabolic networks of arbitrary size and complexity, including genome-scale kinetic models. We report that optimal states are elementary flux modes, which are minimal metabolic networks operating at a thermodynamically-feasible steady state with one independent flux. Remarkably, elementary flux modes rely only on reaction stoichiometry, yet they function as the optimal states of mathematical models incorporating enzyme kinetics. Our results pave the way for the optimization of genome-scale kinetic models because they offer huge simplifications to overcome the concomitant computational problems.

Lost in transition: start-up of glycolysis yields subpopulations of nongrowing cells.

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.

Fatal attraction in glycolysis: how Saccharomyces cerevisiae manages sudden transitions to high glucose.

In the model eukaryote Saccharomyces cerevisiae, it has long been known that a functional trehalose pathway is indispensable for transitions to high glucose conditions. Upon addition of glucose, cells with a defect in trehalose 6-phosphate synthase (Tps1), the first committed step in the trehalose pathway, display what we have termed an imbalanced glycolytic state; in this state the flux through the upper part of glycolysis outpaces that through the lower part of glycolysis. As a consequence, the intermediate fructose 1,6-bisphosphate (FBP) accumulates at low concentrations of ATP and inorganic phosphate (Pi). Despite significant research efforts, a satisfactory understanding of the regulatory role that trehalose metabolism plays during such transitions has remained infamously unresolved. In a recent study, we demonstrate that the startup of glycolysis exhibits two dynamic fates: a proper, functional, steady state or the imbalanced state described above. Both states are stable, attracting states, and the probability distribution of initial states determines the fate of a yeast cell exposed to glucose. Trehalose metabolism steers the dynamics of glycolysis towards the proper functional state through its ATP hydrolysis activity; a mechanism that ensures that the demand and supply of ATP is balanced with Pi availability under dynamic conditions. [van Heerden et al. Science (2014), DOI: 10.1126/science.1245114.].

The timescale of phenotypic plasticity and its impact on competition in fluctuating environments.

Although phenotypic plasticity can be advantageous in fluctuating environments, it may come too late if the environment changes fast. Complementary chromatic adaptation is a colorful form of phenotypic plasticity, where cyanobacteria tune their pigmentation to the prevailing light spectrum. Here, we study the timescale of chromatic adaptation and its impact on competition among phytoplankton species exposed to fluctuating light colors. We parameterized a resource competition model using monoculture experiments with green and red picocyanobacteria and the cyanobacterium Pseudanabaena, which can change its color within approximately 7 days by chromatic adaptation. The model predictions were tested in competition experiments, where the incident light color switched between red and green at different frequencies (slow, intermediate, and fast). Pseudanabaena (the flexible phenotype) competitively excluded the green and red picocyanobacteria in all competition experiments. Strikingly, the rate of competitive exclusion was much faster when the flexible phenotype had sufficient time to fully adjust its pigmentation. Thus, the flexible phenotype benefited from its phenotypic plasticity if fluctuations in light color were relatively slow, corresponding to slow mixing processes or infrequent storms in their natural habitat. This shows that the timescale of phenotypic plasticity plays a key role during species interactions in fluctuating environments.