Non-equilibrium physics and evolution--adaptation, extinction, and ecology: a key issues review

Phase transitions in evolutionary dynamics

Sharp changes in state, such as transitions from survival to extinction, are hallmarks of evolutionary dynamics in biological systems. These transitions can be explored using the techniques of statistical physics and the physics of nonlinear and complex systems. For example, a survival-to-extinction transition can be characterized as a non-equilibrium phase transition to an absorbing state. Here, we review the literature on phase transitions in evolutionary dynamics. We discuss directed percolation transitions in cellular automata and evolutionary models, and models that diverge from the directed percolation universality class. We explore in detail an example of an absorbing phase transition in an agent-based model of evolutionary dynamics, including previously unpublished data demonstrating similarity to, but also divergence from, directed percolation, as well as evidence for phase transition behavior at multiple levels of the model system's evolutionary structure. We discuss phase transition models of the error catastrophe in RNA virus dynamics and phase transition models for transition from chemistry to biochemistry, i.e., the origin of life. We conclude with a review of phase transition dynamics in models of natural selection, discuss the possible role of phase transitions in unraveling fundamental unresolved questions regarding multilevel selection and the major evolutionary transitions, and assess the future outlook for phase transitions in the investigation of evolutionary dynamics.

Phenotypic-dependent variability and the emergence of tolerance in bacterial populations

Ecological and evolutionary dynamics have been historically regarded as unfolding at broadly separated timescales. However, these two types of processes are nowadays well-documented to intersperse much more tightly than traditionally assumed, especially in communities of microorganisms. Advancing the development of mathematical and computational approaches to shed novel light onto eco-evolutionary problems is a challenge of utmost relevance. With this motivation in mind, here we scrutinize recent experimental results showing evidence of rapid evolution of tolerance by lag in bacterial populations that are periodically exposed to antibiotic stress in laboratory conditions. In particular, the distribution of single-cell lag times-i.e., the times that individual bacteria from the community remain in a dormant state to cope with stress-evolves its average value to approximately fit the antibiotic-exposure time. Moreover, the distribution develops right-skewed heavy tails, revealing the presence of individuals with anomalously large lag times. Here, we develop a parsimonious individual-based model mimicking the actual demographic processes of the experimental setup. Individuals are characterized by a single phenotypic trait: their intrinsic lag time, which is transmitted with variation to the progeny. The model-in a version in which the amplitude of phenotypic variations grows with the parent's lag time-is able to reproduce quite well the key empirical observations. Furthermore, we develop a general mathematical framework allowing us to describe with good accuracy the properties of the stochastic model by means of a macroscopic equation, which generalizes the Crow-Kimura equation in population genetics. Even if the model does not account for all the biological mechanisms (e.g., genetic changes) in a detailed way-i.e., it is a phenomenological one-it sheds light onto the eco-evolutionary dynamics of the problem and can be helpful to design strategies to hinder the emergence of tolerance in bacterial communities. From a broader perspective, this work represents a benchmark for the mathematical framework designed to tackle much more general eco-evolutionary problems, thus paving the road to further research avenues.

Weak mixed phase in the mutator model

We consider the mutator model with unidirected transitions from the wild type to the mutator type, with different fitness functions for the wild types and mutator types. We calculate both the fraction of mutator types in the population and the surpluses, i.e., the mean number of mutations in the regular part of genomes for the wild type and mutator type, which have never been derived exactly. We identify the phase structure. Beside the mixed (ordinary evolution phase with finite fraction of wild types at large genome length) and the mutator phase (the absolute majority is mutators), we find another new phase as well-it has the mean fitness of the mixed phase but an exponentially small (in genome length) fraction of wild types. We identify the phase transition point and discuss its implications.

Tuning environmental timescales to evolve and maintain generalists

Natural environments can present diverse challenges, but some genotypes remain fit across many environments. Such "generalists" can be hard to evolve, outcompeted by specialists fitter in any particular environment. Here, inspired by the search for broadly neutralizing antibodies during B cell affinity maturation, we demonstrate that environmental changes on an intermediate timescale can reliably evolve generalists, even when faster or slower environmental changes are unable to do so. We find that changing environments on timescales comparable with evolutionary transients in a population enhance the rate of evolving generalists from specialists, without enhancing the reverse process. The yield of generalists is further increased in more complex dynamic environments, such as a "chirp" of increasing frequency. Our work offers design principles for how nonequilibrium fitness "seascapes" can dynamically funnel populations to genotypes unobtainable in static environments.

Evolving generalists in switching rugged landscapes

Evolving systems, be it an antibody repertoire in the face of mutating pathogens or a microbial population exposed to varied antibiotics, constantly search for adaptive solutions in time-varying fitness landscapes. Generalists refer to genotypes that remain fit across diverse selective pressures; while multi-drug resistant microbes are undesired yet prevalent, broadly-neutralizing antibodies are much wanted but rare. However, little is known about under what conditions such generalists with a high capacity to adapt can be efficiently discovered by evolution. In addition, can epistasis-the source of landscape ruggedness and path constraints-play a different role, if the environment varies in a non-random way? We present a generative model to estimate the propensity of evolving generalists in rugged landscapes that are tunably related and alternating relatively slowly. We find that environmental cycling can substantially facilitate the search for fit generalists by dynamically enlarging their effective basins of attraction. Importantly, these high performers are most likely to emerge at intermediate levels of ruggedness and environmental relatedness. Our approach allows one to estimate correlations across environments from the topography of experimental fitness landscapes. Our work provides a conceptual framework to study evolution in time-correlated complex environments, and offers statistical understanding that suggests general strategies for eliciting broadly neutralizing antibodies or preventing microbes from evolving multi-drug resistance.

Solution of the Crow-Kimura model with changing population size and Allee effect

The Crow-Kimura model is commonly used in the modeling of genetic evolution in the presence of mutations and associated selection pressures. We consider a modified version of the Crow-Kimura model, in which population sizes are not fixed and Allee saturation effects are present. We demonstrate the evolutionary dynamics in this system through an analytical approach, examining both symmetric and single-peak fitness landscape cases. Especially interesting are the dynamics of the populations near extinction. A special version of the model with saturation and degradation on the single-peak fitness landscape is investigated as a candidate of the Allee effect in evolution, revealing reduction tendencies of excessively large populations, and extinction tendencies for small populations. The analytical solutions for these dynamics are presented with accuracy O(1/N), where N is the number of nucleotides in the genome.

Sensitivity of asymmetric rate-dependent critical systems to initial conditions: Insights into cellular decision making

The work reported here aims to address the effects of time-dependent parameters and stochasticity on decision making in biological systems. We achieve this by extending previous studies that resorted to simple bifurcation normal forms, although in the present case we focus primarily on the issue of the system's sensitivity to initial conditions in the presence of two different noise distributions, Gaussian and Lévy. In addition, we also assess the impact of two-way sweeping at different rates through the critical region of a canonical Pitchfork bifurcation with a constant external asymmetry. The parallel with decision making in biocircuits is performed on this simple system since it is equivalent in its available states and dynamics to more complex genetic circuits published previously. Overall we verify that rate-dependent effects, previously reported as being important features of bifurcating systems, are specific to particular initial conditions. Processing of each starting state, which for the normal form underlying this study is akin to a classification task, is affected by the balance between sweeping speed through critical regions and the type of fluctuations added. For the heavy-tailed noise, two-way dynamic bifurcations are more efficient in processing the external signals, here understood to be jointly represented by the critical parameter profile and the external asymmetry amplitude, when compared to the system relying on escape dynamics. This is particular to the case when the system starts at an attractor not favored by the asymmetry and, in conjunction, when the sweeping amplitude is large.

Exploiting ecology in drug pulse sequences in favour of population reduction

A deterministic population dynamics model involving birth and death for a two-species system, comprising a wild-type and more resistant species competing via logistic growth, is subjected to two distinct stress environments designed to mimic those that would typically be induced by temporal variation in the concentration of a drug (antibiotic or chemotherapeutic) as it permeates through the population and is progressively degraded. Different treatment regimes, involving single or periodical doses, are evaluated in terms of the minimal population size (a measure of the extinction probability), and the population composition (a measure of the selection pressure for resistance or tolerance during the treatment). We show that there exist timescales over which the low-stress regime is as effective as the high-stress regime, due to the competition between the two species. For multiple periodic treatments, competition can ensure that the minimal population size is attained during the first pulse when the high-stress regime is short, which implies that a single short pulse can be more effective than a more protracted regime. Our results suggest that when the duration of the high-stress environment is restricted, a treatment with one or multiple shorter pulses can produce better outcomes than a single long treatment. If ecological competition is to be exploited for treatments, it is crucial to determine these timescales, and estimate for the minimal population threshold that suffices for extinction. These parameters can be quantified by experiment.

Spatial dispersal of bacterial colonies induces a dynamical transition from local to global quorum sensing

Bacteria communicate using external chemical signals called autoinducers (AI) in a process known as quorum sensing (QS). QS efficiency is reduced by both limitations of AI diffusion and potential interference from neighboring strains. There is thus a need for predictive theories of how spatial community structure shapes information processing in complex microbial ecosystems. As a step in this direction, we apply a reaction-diffusion model to study autoinducer signaling dynamics in a single-species community as a function of the spatial distribution of colonies in the system. We predict a dynamical transition between a local quorum sensing (LQS) regime, with the AI signaling dynamics primarily controlled by the local population densities of individual colonies, and a global quorum sensing (GQS) regime, with the dynamics being dependent on collective intercolony diffusive interactions. The crossover between LQS to GQS is intimately connected to a trade-off between the signaling network's latency, or speed of activation, and its throughput, or the total spatial range over which all the components of the system communicate.

The rich phase structure of a mutator model

We propose a modification of the Crow-Kimura and Eigen models of biological molecular evolution to include a mutator gene that causes both an increase in the mutation rate and a change in the fitness landscape. This mutator effect relates to a wide range of biomedical problems. There are three possible phases: mutator phase, mixed phase and non-selective phase. We calculate the phase structure, the mean fitness and the fraction of the mutator allele in the population, which can be applied to describe cancer development and RNA viruses. We find that depending on the genome length, either the normal or the mutator allele dominates in the mixed phase. We analytically solve the model for a general fitness function. We conclude that the random fitness landscape is an appropriate choice for describing the observed mutator phenomenon in the case of a small fraction of mutators. It is shown that the increase in the mutation rates in the regular and the mutator parts of the genome should be set independently; only some combinations of these increases can push the complex biomedical system to the non-selective phase, potentially related to the eradication of tumors.

Analytically tractable model for community ecology with many species

A fundamental problem in community ecology is understanding how ecological processes such as selection, drift, and immigration give rise to observed patterns in species composition and diversity. Here, we analyze a recently introduced, analytically tractable, presence-absence (PA) model for community assembly, and we use it to ask how ecological traits such as the strength of competition, the amount of diversity, and demographic and environmental stochasticity affect species composition in a community. In the PA model, species are treated as stochastic binary variables that can either be present or absent in a community: species can immigrate into the community from a regional species pool and can go extinct due to competition and stochasticity. Building upon previous work, we show that, despite its simplicity, the PA model reproduces the qualitative features of more complicated models of community assembly. In agreement with recent studies of large, competitive Lotka-Volterra systems, the PA model exhibits distinct ecological behaviors organized around a special ("critical") point corresponding to Hubbell's neutral theory of biodiversity. These results suggest that the concepts of ecological "phases" and phase diagrams can provide a powerful framework for thinking about community ecology, and that the PA model captures the essential ecological dynamics of community assembly.

Collective Fluctuations in the Dynamics of Adaptation and Other Traveling Waves

The dynamics of adaptation are difficult to predict because it is highly stochastic even in large populations. The uncertainty emerges from random genetic drift arising in a vanguard of particularly fit individuals of the population. Several approaches have been developed to analyze the crucial role of genetic drift on the expected dynamics of adaptation, including the mean fitness of the entire population, or the fate of newly arising beneficial deleterious mutations. However, little is known about how genetic drift causes fluctuations to emerge on the population level, where it becomes palpable as variations in the adaptation speed and the fitness distribution. Yet these phenomena control the decay of genetic diversity and variability in evolution experiments and are key to a truly predictive understanding of evolutionary processes. Here, we show that correlations induced by these emergent fluctuations can be computed at any arbitrary order by a suitable choice of a dynamical constraint. The resulting linear equations exhibit fluctuation-induced terms that amplify short-distance correlations and suppress long-distance ones. These terms, which are in general not small, control the decay of genetic diversity and, for wave-tip dominated ("pulled") waves, lead to anticorrelations between the tip of the wave and the lagging bulk of the population. While it is natural to consider the process of adaptation as a branching random walk in fitness space subject to a constraint (due to finite resources), we show that other traveling wave phenomena in ecology and evolution likewise fall into this class of constrained branching random walks. Our methods, therefore, provide a systematic approach toward analyzing fluctuations in a wide range of population biological processes, such as adaptation, genetic meltdown, species invasions, or epidemics.