Trade-offs between microbial growth phases lead to frequency-dependent and non-transitive selection

Dynamic coexistence driven by physiological transitions in microbial communities

Microbial ecosystems are commonly modeled by fixed interactions between species in steady exponential growth states. However, microbes often modify their environments so strongly that they are forced out of the exponential state into stressed or non-growing states. Such dynamics are typical of ecological succession in nature and serial-dilution cycles in the laboratory. Here, we introduce a phenomenological model, the Community State model, to gain insight into the dynamic coexistence of microbes due to changes in their physiological states. Our model bypasses specific interactions (e.g., nutrient starvation, stress, aggregation) that lead to different combinations of physiological states, referred to collectively as "community states", and modeled by specifying the growth preference of each species along a global ecological coordinate, taken here to be the total community biomass density. We identify three key features of such dynamical communities that contrast starkly with steady-state communities: increased tolerance of community diversity to fast growth rates of species dominating different community states, enhanced community stability through staggered dominance of different species in different community states, and increased requirement on growth dominance for the inclusion of late-growing species. These features, derived explicitly for simplified models, are proposed here to be principles aiding the understanding of complex dynamical communities. Our model shifts the focus of ecosystem dynamics from bottom-up studies based on idealized inter-species interaction to top-down studies based on accessible macroscopic observables such as growth rates and total biomass density, enabling quantitative examination of community-wide characteristics.

Dynamic coexistence driven by physiological transitions in microbial communities

Microbial ecosystems are commonly modeled by fixed interactions between species in steady exponential growth states. However, microbes often modify their environments so strongly that they are forced out of the exponential state into stressed or non-growing states. Such dynamics are typical of ecological succession in nature and serial-dilution cycles in the laboratory. Here, we introduce a phenomenological model, the Community State model, to gain insight into the dynamic coexistence of microbes due to changes in their physiological states. Our model bypasses specific interactions (e.g., nutrient starvation, stress, aggregation) that lead to different combinations of physiological states, referred to collectively as "community states", and modeled by specifying the growth preference of each species along a global ecological coordinate, taken here to be the total community biomass density. We identify three key features of such dynamical communities that contrast starkly with steady-state communities: increased tolerance of community diversity to fast growth rates of species dominating different community states, enhanced community stability through staggered dominance of different species in different community states, and increased requirement on growth dominance for the inclusion of late-growing species. These features, derived explicitly for simplified models, are proposed here to be principles aiding the understanding of complex dynamical communities. Our model shifts the focus of ecosystem dynamics from bottom-up studies based on idealized inter-species interaction to top-down studies based on accessible macroscopic observables such as growth rates and total biomass density, enabling quantitative examination of community-wide characteristics.

Frequent asymmetric migrations suppress natural selection in spatially structured populations

Natural microbial populations often have complex spatial structures. This can impact their evolution, in particular the ability of mutants to take over. While mutant fixation probabilities are known to be unaffected by sufficiently symmetric structures, evolutionary graph theory has shown that some graphs can amplify or suppress natural selection, in a way that depends on microscopic update rules. We propose a model of spatially structured populations on graphs directly inspired by batch culture experiments, alternating within-deme growth on nodes and migration-dilution steps, and yielding successive bottlenecks. This setting bridges models from evolutionary graph theory with Wright-Fisher models. Using a branching process approach, we show that spatial structure with frequent migrations can only yield suppression of natural selection. More precisely, in this regime, circulation graphs, where the total incoming migration flow equals the total outgoing one in each deme, do not impact fixation probability, while all other graphs strictly suppress selection. Suppression becomes stronger as the asymmetry between incoming and outgoing migrations grows. Amplification of natural selection can nevertheless exist in a restricted regime of rare migrations and very small fitness advantages, where we recover the predictions of evolutionary graph theory for the star graph.

Microbial population dynamics decouple growth response from environmental nutrient concentration

How the growth rate of a microbial population responds to the environmental availability of chemical nutrients and other resources is a fundamental question in microbiology. Models of this response, such as the widely used Monod model, are generally characterized by a maximum growth rate and a half-saturation concentration of the resource. What values should we expect for these half-saturation concentrations, and how should they depend on the environmental concentration of the resource? We survey growth response data across a wide range of organisms and resources. We find that the half-saturation concentrations vary across orders of magnitude, even for the same organism and resource. To explain this variation, we develop an evolutionary model to show that demographic fluctuations (genetic drift) can constrain the adaptation of half-saturation concentrations. We find that this effect fundamentally differs depending on the type of population dynamics: Populations undergoing periodic bottlenecks of fixed size will adapt their half-saturation concentrations in proportion to the environmental resource concentrations, but populations undergoing periodic dilutions of fixed size will evolve half-saturation concentrations that are largely decoupled from the environmental concentrations. Our model not only provides testable predictions for laboratory evolution experiments, but it also reveals how an evolved half-saturation concentration may not reflect the organism's environment. In particular, this explains how organisms in resource-rich environments can still evolve fast growth at low resource concentrations. Altogether, our results demonstrate the critical role of population dynamics in shaping fundamental ecological traits.

Exploring the expanse between theoretical questions and experimental approaches in the modern study of evolvability

Despite several decades of computational and experimental work across many systems, evolvability remains on the periphery with regards to its status as a widely accepted and regularly applied theoretical concept. Here we propose that its marginal status is partly a result of large gaps between the diverse but disconnected theoretical treatments of evolvability and the relatively narrower range of studies that have tested it empirically. To make this case, we draw on a range of examples-from experimental evolution in microbes, to molecular evolution in proteins-where attempts have been made to mend this disconnect. We highlight some examples of progress that has been made and point to areas where synthesis and translation of existing theory can lead to further progress in the still-new field of empirical measurements of evolvability.

Application of Electrolyzed Water in the Food Industry: A Review

Electrolyzed water is a novel disinfectant and cleaner that has been widely utilized in the food sector for several years to ensure that surfaces are sterilized, and that food is safe. It is produced by the electrolysis of a dilute salt solution, and the reaction products include sodium hydroxide (NaOH) and hypochlorous acid. In comparison to conventional cleaning agents, electrolyzed water is economical and eco-friendly, easy to use, and strongly effective. Electrolyzed water is also used in its acidic form, but it is non-corrosive to the human epithelium and other organic matter. The electrolyzed water can be utilized in a diverse range of foods; thus, it is an appropriate choice for synergistic microbial control in the food industry to ensure food safety and quality without damaging the organoleptic parameters of the food. The present review article highlights the latest information on the factors responsible for food spoilage and the antimicrobial potential of electrolyzed water in fresh or processed plant and animal products.

Diauxic lags explain unexpected coexistence in multi-resource environments

How the coexistence of species is affected by the presence of multiple resources is a major question in microbial ecology. We experimentally demonstrate that differences in diauxic lags, which occur as species deplete their own environments and adapt their metabolisms, allow slow-growing microbes to stably coexist with faster-growing species in multi-resource environments despite being excluded in single-resource environments. In our focal example, an Acinetobacter species (Aci2) competitively excludes Pseudomonas aurantiaca (Pa) on alanine and on glutamate. However, they coexist on the combination of both resources. Experiments reveal that Aci2 grows faster but Pa has shorter diauxic lags. We establish a tradeoff between Aci2's fast growth and Pa's short lags as their mechanism for coexistence. We model this tradeoff to accurately predict how environmental changes affect community composition. We extend our work by surveying a large set of competitions and observe coexistence nearly four times as frequently when the slow-grower is the fast-switcher. Our work illustrates a simple mechanism, based entirely on supplied-resource growth dynamics, for the emergence of multi-resource coexistence.

Community lifespan, niche expansion and the evolution of interspecific cooperation

Microbes live in dense and diverse communities where they deploy many traits that promote the growth and survival of neighbouring species, all the while also competing for shared resources. Because microbial communities are highly dynamic, the costs and benefits of species interactions change over the growth cycle of a community. How mutualistic interactions evolve under such demographic and ecological conditions is still poorly understood. Here, we develop an eco-evolutionary model to explore how different forms of helping with distinct fitness effects (rate-enhancing and yield-enhancing) affect the multiple phases of community growth, and its consequences for the evolution of mutualisms. We specifically focus on a form of yield-enhancing trait in which cooperation augments the common pool of resources, termed niche expansion. We show that although mutualisms in which cooperation increases partners growth rate are generally favoured at early stages of community growth, niche expansion can evolve at later stages where densities are high. Further, we find that niche expansion can promote the evolution of reproductive restraint, in which a focal species adaptively reduces its own growth rate to increase the density of partner species. Our findings suggest that yield-enhancing mutualisms are more prevalent in stable habitats with a constant supply of resources, and where populations typically live at high densities. In general, our findings highlight the need to integrate different components of population growth in the analysis of mutualisms to understand the composition and function of microbial communities.

Evolutionary stalling and a limit on the power of natural selection to improve a cellular module

Cells consist of molecular modules which perform vital biological functions. Cellular modules are key units of adaptive evolution because organismal fitness depends on their performance. Theory shows that in rapidly evolving populations, such as those of many microbes, adaptation is driven primarily by common beneficial mutations with large effects, while other mutations behave as if they are effectively neutral. As a consequence, if a module can be improved only by rare and/or weak beneficial mutations, its adaptive evolution would stall. However, such evolutionary stalling has not been empirically demonstrated, and it is unclear to what extent stalling may limit the power of natural selection to improve modules. Here we empirically characterize how natural selection improves the translation machinery (TM), an essential cellular module. We experimentally evolved populations of Escherichia coli with genetically perturbed TMs for 1,000 generations. Populations with severe TM defects initially adapted via mutations in the TM, but TM adaptation stalled within about 300 generations. We estimate that the genetic load in our populations incurred by residual TM defects ranges from 0.5 to 19%. Finally, we found evidence that both epistasis and the depletion of the pool of beneficial mutations contributed to evolutionary stalling. Our results suggest that cellular modules may not be fully optimized by natural selection despite the availability of adaptive mutations.

Collateral fitness effects of mutations

The distribution of fitness effects of mutation plays a central role in constraining protein evolution. The underlying mechanisms by which mutations lead to fitness effects are typically attributed to changes in protein specific activity or abundance. Here, we reveal the importance of a mutation's collateral fitness effects, which we define as effects that do not derive from changes in the protein's ability to perform its physiological function. We comprehensively measured the collateral fitness effects of missense mutations in the Escherichia coli TEM-1 β-lactamase antibiotic resistance gene using growth competition experiments in the absence of antibiotic. At least 42% of missense mutations in TEM-1 were deleterious, indicating that for some proteins collateral fitness effects occur as frequently as effects on protein activity and abundance. Deleterious mutations caused improper posttranslational processing, incorrect disulfide-bond formation, protein aggregation, changes in gene expression, and pleiotropic effects on cell phenotype. Deleterious collateral fitness effects occurred more frequently in TEM-1 than deleterious effects on antibiotic resistance in environments with low concentrations of the antibiotic. The surprising prevalence of deleterious collateral fitness effects suggests they may play a role in constraining protein evolution, particularly for highly expressed proteins, for proteins under intermittent selection for their physiological function, and for proteins whose contribution to fitness is buffered against deleterious effects on protein activity and protein abundance.

Evolution of Microbial Growth Traits Under Serial Dilution

Selection of mutants in a microbial population depends on multiple cellular traits. In serial-dilution evolution experiments, three key traits are the lag time when transitioning from starvation to growth, the exponential growth rate, and the yield (number of cells per unit resource). Here, we investigate how these traits evolve in laboratory evolution experiments using a minimal model of population dynamics, where the only interaction between cells is competition for a single limiting resource. We find that the fixation probability of a beneficial mutation depends on a linear combination of its growth rate and lag time relative to its immediate ancestor, even under clonal interference. The relative selective pressure on growth rate and lag time is set by the dilution factor; a larger dilution factor favors the adaptation of growth rate over the adaptation of lag time. The model shows that yield, however, is under no direct selection. We also show how the adaptation speeds of growth and lag depend on experimental parameters and the underlying supply of mutations. Finally, we investigate the evolution of covariation between these traits across populations, which reveals that the population growth rate and lag time can evolve a nonzero correlation even if mutations have uncorrelated effects on the two traits. Altogether these results provide useful guidance to future experiments on microbial evolution.

Dynamic metabolic adaptation can promote species coexistence in competitive microbial communities

Microbes are capable of physiologically adapting to diverse environmental conditions by differentially varying the rates at which they uptake different nutrients. In particular, microbes can switch hierarchically between different energy sources, consuming first those that ensure the highest growth rate. Experimentally, this can result in biphasic growth curves called "diauxic shifts" that typically arise when microbes are grown in media containing several nutrients. Despite these observations are well known in microbiology and molecular biology, the mathematical models generally used to describe the population dynamics of microbial communities do not account for dynamic metabolic adaptation, thus implicitly assuming that microbes cannot switch dynamically from one resource to another. Here, we introduce dynamic metabolic adaptation in the framework of consumer-resource models, which are commonly used to describe competitive microbial communities, allowing each species to temporally change its preferred energy source to maximize its own relative fitness. We show that dynamic metabolic adaptation enables the community to self-organize, allowing several species to coexist even in the presence of few resources, and to respond optimally to a time-dependent environment, thus showing that dynamic metabolic adaptation could be an important mechanism for maintaining high levels of diversity even in environments with few energy sources. We show that introducing dynamic metabolic strategies in consumer-resource models is necessary for reproducing experimental growth curves of the baker's yeast Saccharomyces cerevisiae growing in the presence of two carbon sources. Even though diauxic shifts emerge naturally from the model when two resources are qualitatively very different, the model predicts that the existence of such shifts is not a prerequisite for species coexistence in competitive communities.

Coexistence and cooperation in structured habitats

Background

Natural habitats are typically structured, imposing constraints on inhabiting populations and their interactions. Which conditions are important for coexistence of diverse communities, and how cooperative interaction stabilizes in such populations, have been important ecological and evolutionary questions.

Results

We investigate a minimal ecological framework of microbial population dynamics that exhibits crucial features to show coexistence: Populations repeatedly undergo cycles of separation into compartmentalized habitats and mixing with new resources. The characteristic time-scale is longer than that typical of individual growth. Using analytic approximations, averaging techniques and phase-plane methods of dynamical systems, we provide a framework for analyzing various types of microbial interactions. Population composition and population size are both dynamic variables of the model; they are found to be decoupled both in terms of time-scale and parameter dependence. We present specific results for two examples of cooperative interaction by public goods: collective antibiotics resistance, and enhanced iron-availability by pyoverdine. We find stable coexistence to be a likely outcome.

Conclusions

The two simple features of a long mixing time-scale and spatial compartmentalization are enough to enable coexisting strains. In particular, costly social traits are often stabilized in such an environment—and thus cooperation established.

Trusting the hand that feeds: microbes evolve to anticipate a serial transfer protocol as individuals or collectives

BACKGROUND

Experimental evolution of microbes often involves a serial transfer protocol, where microbes are repeatedly diluted by transfer to a fresh medium, starting a new growth cycle. This has revealed that evolution can be remarkably reproducible, where microbes show parallel adaptations both on the level of the phenotype as well as the genotype. However, these studies also reveal a strong potential for divergent evolution, leading to diversity both between and within replicate populations. We here study how in silico evolved Virtual Microbe "wild types" (WTs) adapt to a serial transfer protocol to investigate generic evolutionary adaptations, and how these adaptations can be manifested by a variety of different mechanisms.

RESULTS

We show that all WTs evolve to anticipate the regularity of the serial transfer protocol by adopting a fine-tuned balance of growth and survival. This anticipation is done by evolving either a high yield mode, or a high growth rate mode. We find that both modes of anticipation can be achieved by individual lineages and by collectives of microbes. Moreover, these different outcomes can be achieved with or without regulation, although the individual-based anticipation without regulation is less well adapted in the high growth rate mode.

CONCLUSIONS

All our in silico WTs evolve to trust the hand that feeds by evolving to anticipate the periodicity of a serial transfer protocol, but can do so by evolving two distinct growth strategies. Furthermore, both these growth strategies can be accomplished by gene regulation, a variety of different polymorphisms, and combinations thereof. Our work reveals that, even under controlled conditions like those in the lab, it may not be possible to predict individual evolutionary trajectories, but repeated experiments may well result in only a limited number of possible outcomes.

Time Series Genome-Centric Analysis Unveils Bacterial Response to Operational Disturbance in Activated Sludge

Understanding ecosystem response to disturbances and identifying the most critical traits for the maintenance of ecosystem functioning are important goals for microbial community ecology. In this study, we used 16S rRNA amplicon sequencing and metagenomics to investigate the assembly of bacterial populations in a full-scale municipal activated sludge wastewater treatment plant over a period of 3 years, including a 9-month period of disturbance characterized by short-term plant shutdowns. Following the reconstruction of 173 metagenome-assembled genomes, we assessed the functional potential, the number of rRNA gene operons, and the in situ growth rate of microorganisms present throughout the time series. Operational disturbances caused a significant decrease in bacteria with a single copy of the rRNA (rrn) operon. Despite moderate differences in resource availability, replication rates were distributed uniformly throughout time, with no differences between disturbed and stable periods. We suggest that the length of the growth lag phase, rather than the growth rate, is the primary driver of selection under disturbed conditions. Thus, the system could maintain its function in the face of disturbance by recruiting bacteria with the capacity to rapidly resume growth under unsteady operating conditions.IMPORTANCE Disturbance is a key determinant of community assembly and dynamics in natural and engineered ecosystems. Microbiome response to disturbance is thought to be influenced by bacterial growth traits and life history strategies. In this time series observational study, the response to disturbance of microbial communities in a full-scale activated sludge wastewater treatment plant was assessed by computing specific cellular traits of genomes retrieved from metagenomes. It was found that the genomes observed in disturbed periods have more copies of the rRNA operon than genomes observed in stable periods, whereas the in situ mean relative growth rates of bacteria present during stable and disturbed periods were indistinguishable. From these intriguing observations, we infer that the length of the lag phase might be a growth trait that affects the microbial response to disturbance. Further exploration of this hypothesis could contribute to better understanding of the adaptive response of microbiomes to unsteady environmental conditions.

Disentangling bacterial invasiveness from lethality in an experimental host-pathogen system

Quantifying virulence remains a central problem in human health, pest control, disease ecology, and evolutionary biology. Bacterial virulence is typically quantified by the LT50 (i.e., the time taken to kill 50% of infected hosts); however, such an indicator cannot account for the full complexity of the infection process, such as distinguishing between the pathogen's ability to colonize versus kill the hosts. Indeed, the pathogen needs to breach the primary defenses in order to colonize, find a suitable environment to replicate, and finally express the virulence factors that cause disease. Here, we show that two virulence attributes, namely pathogen lethality and invasiveness, can be disentangled from the survival curves of a laboratory population of Caenorhabditis elegans nematodes exposed to three bacterial pathogens: Pseudomonas aeruginosa, Serratia marcescens, and Salmonella enterica We first show that the host population eventually experiences a constant mortality rate, which quantifies the lethality of the pathogen. We then show that the time necessary to reach this constant mortality rate regime depends on the pathogen growth rate and colonization rate, and thus determines the pathogen invasiveness. Our framework reveals that Serratia marcescens is particularly good at the initial colonization of the host, whereas Salmonella enterica is a poor colonizer yet just as lethal once established. Pseudomonas aeruginosa, on the other hand, is both a good colonizer and highly lethal after becoming established. The ability to quantitatively characterize the ability of different pathogens to perform each of these steps has implications for treatment and prevention of disease and for the evolution and ecology of pathogens.

Mortality causes universal changes in microbial community composition

All organisms are sensitive to the abiotic environment, and a deteriorating environment can cause extinction. However, survival in a multispecies community depends upon interactions, and some species may even be favored by a harsh environment that impairs others, leading to potentially surprising community transitions as environments deteriorate. Here we combine theory and laboratory microcosms to predict how simple microbial communities will change under added mortality, controlled by varying dilution. We find that in a two-species coculture, increasing mortality favors the faster grower, confirming a theoretical prediction. Furthermore, if the slower grower dominates under low mortality, the outcome can reverse as mortality increases. We find that this tradeoff between growth and competitive ability is prevalent at low dilution, causing outcomes to shift dramatically as dilution increases, and that these two-species shifts propagate to simple multispecies communities. Our results argue that a bottom-up approach can provide insight into how communities change under stress.

Environmental stress can affect the outcome of ecological competition. Here, the authors use theory and experiments with a synthetic microbial community to show that a tradeoff between growth rate and competitive ability determines which species prevails when the population faces variable mortality rates.

Growth tradeoffs produce complex microbial communities on a single limiting resource

The relationship between the dynamics of a community and its constituent pairwise interactions is a fundamental problem in ecology. Higher-order ecological effects beyond pairwise interactions may be key to complex ecosystems, but mechanisms to produce these effects remain poorly understood. Here we model microbial growth and competition to show that higher-order effects can arise from variation in multiple microbial growth traits, such as lag times and growth rates, on a single limiting resource with no other interactions. These effects produce a range of ecological phenomena: an unlimited number of strains can exhibit multistability and neutral coexistence, potentially with a single keystone strain; strains that coexist in pairs do not coexist all together; and a strain that wins all pairwise competitions can go extinct in a mixed competition. Since variation in multiple growth traits is ubiquitous in microbial populations, our results indicate these higher-order effects may also be widespread, especially in laboratory ecology and evolution experiments.

Higher-order interactions occur when one species mediates the interaction between two others. Here, the authors model microbial growth and competition to show that higher-order interactions can arise from tradeoffs in growth traits, leading to neutral coexistence and other complex dynamics.

Optimization of lag phase shapes the evolution of a bacterial enzyme

Mutations provide the variation that drives evolution, yet their effects on fitness remain poorly understood. Here we explore how mutations in the essential enzyme adenylate kinase (Adk) of Escherichia coli affect multiple phases of population growth. We introduce a biophysical fitness landscape for these phases, showing how they depend on molecular and cellular properties of Adk. We find that Adk catalytic capacity in the cell (the product of activity and abundance) is the major determinant of mutational fitness effects. We show that bacterial lag times are at a well-defined optimum with respect to Adk's catalytic capacity, while exponential growth rates are only weakly affected by variation in Adk. Direct pairwise competitions between strains show how environmental conditions modulate the outcome of a competition where growth rates and lag times have a tradeoff, shedding light on the multidimensional nature of fitness and its importance in the evolutionary optimization of enzymes.