Evolution of a Fluctuating Population in a Randomly Switching Environment

Role of Noise-Modulated Self-Propulsion in Driving Spatiotemporal Orders in Active Systems

Fluctuations play a pivotal role in driving spatiotemporal order in active matter systems. In this study, we employ a novel analytical framework to investigate the impact of dichotomous noise on the self-propelling velocity of active particle systems such as polymerizing actin filaments or reproducing elongated bacteria. By incorporating dichotomous fluctuations with Ornstein-Zernike correlations into a continuum-based model, we derive a bifurcation condition in the noise parameter space, revealing a noise-induced instability that drives the emergence of traveling waves. This approach demonstrates how specific noise strengths and correlation times expand the instability region by introducing effective new degrees of freedom that alter the system's stability matrix. Advance numerical simulations, meticulously designed to handle the properties of dichotomous noise, validate these theoretical predictions and reveal excellent agreement. A key finding is the observation of wave-reversal behavior, driven by the sign alternation of the noise-modulated advection term and modulated by the relaxation time. Remarkably, we identify a finite parameter range where this reversal is suppressed, offering new insights into noise-induced bifurcations and spatiotemporal dynamics. Our combined analytical and numerical approach provides a deeper understanding of the role of noise in shaping self-organization and pattern formation in biological and synthetic active systems.

Population dynamics in a time-varying environment with fat-tailed correlations

Temporal environmental noise (EN) is a prevalent natural phenomenon that controls population and community dynamics, shaping the destiny of biological species and genetic types. Conventional theoretical models often depict EN as a Markovian process with an exponential distribution of correlation times, resulting in two distinct qualitative dynamical categories: quenched (long environmental timescales) and annealed (short environmental timescales). However, numerous empirical studies demonstrate a fat-tailed decay of correlation times. Here we study the consequences of power-law correlated EN on the dynamics of isolated and competing populations. We analyze the intermediate region that lies between the quenched and annealed regimes and show that it can widen indefinitely. Within this region, dynamics is primarily driven by rare, yet not exceedingly rare, long periods of almost-steady environmental conditions. For an isolated population, the time to extinction in this region exhibits a power-law scaling with the logarithm of the abundance and also a nonmonotonic dependence on the spectral exponent.

Hypochaos prevents tragedy of the commons in discrete-time eco-evolutionary game dynamics

While quite a few recent papers have explored game-resource feedback using the framework of evolutionary game theory, almost all the studies are confined to using time-continuous dynamical equations. Moreover, in such literature, the effect of ubiquitous chaos in the resulting eco-evolutionary dynamics is rather missing. Here, we present a deterministic eco-evolutionary discrete-time dynamics in generation-wise non-overlapping population of two types of harvesters-one harvesting at a faster rate than the other-consuming a self-renewing resource capable of showing chaotic dynamics. In the light of our finding that sometimes chaos is confined exclusively to either the dynamics of the resource or that of the consumer fractions, an interesting scenario is realized: The resource state can keep oscillating chaotically, and hence, it does not vanish to result in the tragedy of the commons-extinction of the resource due to selfish indiscriminate exploitation-and yet the consumer population, whose dynamics depends directly on the state of the resource, may end up being composed exclusively of defectors, i.e., high harvesters. This appears non-intuitive because it is well known that prevention of tragedy of the commons usually requires substantial cooperation to be present.

Coupled environmental and demographic fluctuations shape the evolution of cooperative antimicrobial resistance

There is a pressing need to better understand how microbial populations respond to antimicrobial drugs, and to find mechanisms to possibly eradicate antimicrobial-resistant cells. The inactivation of antimicrobials by resistant microbes can often be viewed as a cooperative behaviour leading to the coexistence of resistant and sensitive cells in large populations and static environments. This picture is, however, greatly altered by the fluctuations arising in volatile environments, in which microbial communities commonly evolve. Here, we study the eco-evolutionary dynamics of a population consisting of an antimicrobial-resistant strain and microbes sensitive to antimicrobial drugs in a time-fluctuating environment, modelled by a carrying capacity randomly switching between states of abundance and scarcity. We assume that antimicrobial resistance (AMR) is a shared public good when the number of resistant cells exceeds a certain threshold. Eco-evolutionary dynamics is thus characterised by demographic noise (birth and death events) coupled to environmental fluctuations which can cause population bottlenecks. By combining analytical and computational means, we determine the environmental conditions for the long-lived coexistence and fixation of both strains, and characterise a fluctuation-driven AMR eradication mechanism, where resistant microbes experience bottlenecks leading to extinction. We also discuss the possible applications of our findings to laboratory-controlled experiments.

Periodic temporal environmental variations induce coexistence in resource competition models

Natural ecosystems, in particular on the microbial scale, are inhabited by a large number of species. The population size of each species is affected by interactions of individuals with each other and by spatial and temporal changes in environmental conditions, such as resource abundance. Here, we use a generic population dynamics model to study how, and under what conditions, a periodic temporal environmental variation can alter an ecosystem's composition and biodiversity. We demonstrate that using timescale separation allows one to qualitatively predict the long-term population dynamics of interacting species in varying environments. We show that the notion of Tilman's R* rule, a well-known principle that applies for constant environments, can be extended to periodically varying environments if the timescale of environmental changes (e.g., seasonal variations) is much faster than the timescale of population growth (doubling time in bacteria). When these timescales are similar, our analysis shows that a varying environment deters the system from reaching a steady state, and stable coexistence between multiple species becomes possible. Our results posit that biodiversity can in part be attributed to natural environmental variations.

Lotka-Volterra predator-prey model with periodically varying carrying capacity

We study the stochastic spatial Lotka-Volterra model for predator-prey interaction subject to a periodically varying carrying capacity. The Lotka-Volterra model with on-site lattice occupation restrictions (i.e., finite local carrying capacity) that represent finite food resources for the prey population exhibits a continuous active-to-absorbing phase transition. The active phase is sustained by the existence of spatiotemporal patterns in the form of pursuit and evasion waves. Monte Carlo simulations on a two-dimensional lattice are utilized to investigate the effect of seasonal variations of the environment on species coexistence. The results of our simulations are also compared to a mean-field analysis in order to specifically delineate the impact of stochastic fluctuations and spatial correlations. We find that the parameter region of predator and prey coexistence is enlarged relative to the stationary situation when the carrying capacity varies periodically. The (quasi-)stationary regime of our periodically varying Lotka-Volterra predator-prey system shows qualitative agreement between the stochastic model and the mean-field approximation. However, under periodic carrying capacity-switching environments, the mean-field rate equations predict period-doubling scenarios that are washed out by internal reaction noise in the stochastic lattice model. Utilizing visual representations of the lattice simulations and dynamical correlation functions, we study how the pursuit and evasion waves are affected by ensuing resonance effects. Correlation function measurements indicate a time delay in the response of the system to sudden changes in the environment. Resonance features are observed in our simulations that cause prolonged persistent spatial correlations. Different effective static environments are explored in the extreme limits of fast and slow periodic switching. The analysis of the mean-field equations in the fast-switching regime enables a semiquantitative description of the (quasi-)stationary state.

Controlling the Mean Time to Extinction in Populations of Bacteria

Populations of ecological systems generally have demographic fluctuations due to birth and death processes. At the same time, they are exposed to changing environments. We studied populations composed of two phenotypes of bacteria and analyzed the impact that both types of fluctuations have on the mean time to extinction of the entire population if extinction is the final fate. Our results are based on Gillespie simulations and on the WKB approach applied to classical stochastic systems, here in certain limiting cases. As a function of the frequency of environmental changes, we observe a non-monotonic dependence of the mean time to extinction. Its dependencies on other system parameters are also explored. This allows the control of the mean time to extinction to be as large or as small as possible, depending on whether extinction should be avoided or is desired from the perspective of bacteria or the perspective of hosts to which the bacteria are deleterious.

Evolution in fluctuating environments: A generic modular approach

Evolutionary processes take place in fluctuating environments, where carrying capacities and selective forces vary over time. The fate of a mutant type and the persistence time of polymorphic states were studied in some specific cases of varying environments, but a generic methodology is still lacking. Here, we present such a general analytic framework. We first identify a set of elementary building blocks, a few basic demographic processes like logistic or exponential growth, competition at equilibrium, sudden decline, and so on. For each of these elementary blocks, we evaluate the mean and the variance of the changes in the frequency of the mutant population. Finally, we show how to find the relevant terms of the diffusion equation for each arbitrary combination of these blocks. Armed with this technique one may calculate easily the quantities that govern the evolutionary dynamics, like the chance of ultimate fixation, the time to absorption, and the time to fixation.

Mutual information in changing environments: Nonlinear interactions, out-of-equilibrium systems, and continuously varying diffusivities

Biochemistry, ecology, and neuroscience are examples of prominent fields aiming at describing interacting systems that exhibit nontrivial couplings to complex, ever-changing environments. We have recently shown that linear interactions and a switching environment are encoded separately in the mutual information of the overall system. Here we first generalize these findings to a broad class of nonlinear interacting models. We find that a new term in the mutual information appears, quantifying the interplay between nonlinear interactions and environmental changes, and leading to either constructive or destructive information interference. Furthermore, we show that a higher mutual information emerges in out-of-equilibrium environments with respect to an equilibrium scenario. Finally, we generalize our framework to the case of continuously varying environments. We find that environmental changes can be mapped exactly into an effective spatially varying diffusion coefficient, shedding light on modeling of biophysical systems in inhomogeneous media.

Mutual Information Disentangles Interactions from Changing Environments

Real-world systems are characterized by complex interactions of their internal degrees of freedom, while living in ever-changing environments whose net effect is to act as additional couplings. Here, we introduce a paradigmatic interacting model in a switching, but unobserved, environment. We show that the limiting properties of the mutual information of the system allow for a disentangling of these two sources of couplings. Further, our approach might stand as a general method to discriminate complex internal interactions from equally complex changing environments.

Game-environment feedback dynamics in growing population: Effect of finite carrying capacity

The tragedy of the commons (TOC) is an unfortunate situation where a shared resource is exhausted due to uncontrolled exploitation by the selfish individuals of a population. Recently, the paradigmatic replicator equation has been used in conjunction with a phenomenological equation for the state of the shared resource to gain insight into the influence of the games on the TOC. The replicator equation, by construction, models a fixed infinite population undergoing microevolution. Thus, it is unable to capture any effect of the population growth and the carrying capacity of the population although the TOC is expected to be dependent on the size of the population. Therefore, in this paper, we present a mathematical framework that incorporates the density dependent payoffs and the logistic growth of the population in the eco-evolutionary dynamics modeling the game-resource feedback. We discover a bistability in the dynamics: a finite carrying capacity can either avert or cause the TOC depending on the initial states of the resource and the initial fraction of cooperators. In fact, depending on the type of strategic game-theoretic interaction, a finite carrying capacity can either avert or cause the TOC when it is exactly the opposite for the corresponding case with infinite carrying capacity.

Exclusion of the fittest predicts microbial community diversity in fluctuating environments

Microorganisms live in environments that inevitably fluctuate between mild and harsh conditions. As harsh conditions may cause extinctions, the rate at which fluctuations occur can shape microbial communities and their diversity, but we still lack an intuition on how. Here, we build a mathematical model describing two microbial species living in an environment where substrate supplies randomly switch between abundant and scarce. We then vary the rate of switching as well as different properties of the interacting species, and measure the probability of the weaker species driving the stronger one extinct. We find that this probability increases with the strength of demographic noise under harsh conditions and peaks at either low, high, or intermediate switching rates depending on both species' ability to withstand the harsh environment. This complex relationship shows why finding patterns between environmental fluctuations and diversity has historically been difficult. In parameter ranges where the fittest species was most likely to be excluded, however, the beta diversity in larger communities also peaked. In sum, how environmental fluctuations affect interactions between a few species pairs predicts their effect on the beta diversity of the whole community.

Beyond the adiabatic limit in systems with fast environments: A τ-leaping algorithm

We propose a τ-leaping simulation algorithm for stochastic systems subject to fast environmental changes. Similar to conventional τ-leaping the algorithm proceeds in discrete time steps, but as a principal addition it captures environmental noise beyond the adiabatic limit. The key idea is to treat the input rates for the τ-leaping as (clipped) Gaussian random variables with first and second moments constructed from the environmental process. In this way, each step of the algorithm retains environmental stochasticity to subleading order in the timescale separation between system and environment. We test the algorithm on several toy examples with discrete and continuous environmental states and find good performance in the regime of fast environmental dynamics. At the same time, the algorithm requires significantly less computing time than full simulations of the combined system and environment. In this context we also discuss several methods for the simulation of stochastic population dynamics in time-varying environments with continuous states.

Aggregative cycles evolve as a solution to conflicts in social investment

Multicellular organization is particularly vulnerable to conflicts between different cell types when the body forms from initially isolated cells, as in aggregative multicellular microbes. Like other functions of the multicellular phase, coordinated collective movement can be undermined by conflicts between cells that spend energy in fuelling motion and ‘cheaters’ that get carried along. The evolutionary stability of collective behaviours against such conflicts is typically addressed in populations that undergo extrinsically imposed phases of aggregation and dispersal. Here, via a shift in perspective, we propose that aggregative multicellular cycles may have emerged as a way to temporally compartmentalize social conflicts. Through an eco-evolutionary mathematical model that accounts for individual and collective strategies of resource acquisition, we address regimes where different motility types coexist. Particularly interesting is the oscillatory regime that, similarly to life cycles of aggregative multicellular organisms, alternates on the timescale of several cell generations phases of prevalent solitary living and starvation-triggered aggregation. Crucially, such self-organized oscillations emerge as a result of evolution of cell traits associated to conflict escalation within multicellular aggregates.

In aggregative multicellular life cycles, cells come together in heterogenous aggregates, whose collective function benefits all the constituent cells. Current explanations for the evolutionary stability of such organization presume that alternating phases of aggregation and dispersal are already in place. Here we propose that, instead of being externally driven, the temporal arrangement of aggregative life cycles may emerge from the interplay between ecology and evolution in populations with differential motility. In our model, cell motility underpins group formation and allows cells to forage individually and collectively. Notably, slower cells can exploit the propulsion by faster cells within multicellular groups. When the level of such exploitation is let evolve, increasing social conflicts are associated to the evolutionary emergence of self-sustained oscillations. Akin to aggregative life cycles, resource exhaustion triggers group formation, whereas conflicts within multicellular groups restrain resource consumption, thus paving the way for the subsequent unicellular phase. The evolutionary transition from equilibrium coexistence to life cycles solves conflicts among heterogenous cell types by integrating them on a timescale longer than cell division, that comes to be associated to multicellular organization.

Stochastic control of single-species population dynamics model subject to jump ambiguity

A logistic type stochastic control model for cost-effective single-species population management subject to an ambiguous jump intensity is presented based on the modern multiplier-robust formulation. The Hamilton-Jacobi-Bellman-Isaacs (HJBI) equation for finding the optimal control is then derived. Mathematical analysis of the HJBI equation from the viewpoint of viscosity solutions is carried out with an emphasis on the non-linear and non-local term, which is a key term arising due to the jump ambiguity. We show that this term can be efficiently handled in the framework of viscosity solutions by utilizing its monotonicity property. A numerical scheme to discretize the HJBI equation is presented as well. Our model is finally applied to management of algae population in river environment. Optimal management policies ranging from the short-term to long-term viewpoints are numerically investigated.

Taming the diffusion approximation through a controlling-factor WKB method

The diffusion approximation (DA) is widely used in the analysis of stochastic population dynamics, from population genetics to ecology and evolution. The DA is an uncontrolled approximation that assumes the smoothness of the calculated quantity over the relevant state space and fails when this property is not satisfied. This failure becomes severe in situations where the direction of selection switches sign. Here we employ the WKB (Wentzel-Kramers-Brillouin) large-deviations method, which requires only the logarithm of a given quantity to be smooth over its state space. Combining the WKB scheme with asymptotic matching techniques, we show how to derive the diffusion approximation in a controlled manner and how to produce better approximations, applicable for much wider regimes of parameters. We also introduce a scalable (independent of population size) WKB-based numerical technique. The method is applied to a central problem in population genetics and evolution, finding the chance of ultimate fixation in a zero-sum, two-types competition.

Adapt or Perish: Evolutionary Rescue in a Gradually Deteriorating Environment

We investigate the evolutionary rescue of a microbial population in a gradually deteriorating environment, through a combination of analytical calculations and stochastic simulations. We consider a population destined for extinction in the absence of mutants, which can survive only if mutants sufficiently adapted to the new environment arise and fix. We show that mutants that appear later during the environment deterioration have a higher probability to fix. The rescue probability of the population increases with a sigmoidal shape when the product of the carrying capacity and of the mutation probability increases. Furthermore, we find that rescue becomes more likely for smaller population sizes and/or mutation probabilities if the environment degradation is slower, which illustrates the key impact of the rapidity of environment degradation on the fate of a population. We also show that our main conclusions are robust across various types of adaptive mutants, including specialist and generalist ones, as well as mutants modeling antimicrobial resistance evolution. We further express the average time of appearance of the mutants that do rescue the population and the average extinction time of those that do not. Our methods can be applied to other situations with continuously variable fitnesses and population sizes, and our analytical predictions are valid in the weak-to-moderate mutation regime.

The impact of random microenvironmental fluctuations on tumor control probability

The tumor control probability (TCP) is a metric used to calculate the probability of controlling or eradicating tumors through radiotherapy. Cancer cells vary in their response to radiation, and although many factors are involved, the tumor microenvironment is a crucial one that determines radiation efficacy. The tumor microenvironment plays a significant role in cancer initiation and propagation, as well as in treatment outcome. We have developed stochastic formulations to study the impact of arbitrary microenvironmental fluctuations on TCP and extinction probability (EP), which is defined as the probability of cancer cells removal in the absence of treatment. Since the derivation of analytical solutions are not possible for complicated cases, we employ a modified Gillespie algorithm to analyze TCP and EP, considering the random variations in cellular proliferation and death rates. Our results show that increasing the standard deviation in kinetic rates initially enhances the probability of tumor eradication. However, if the EP does not reach a probability of 1, the increase in the standard deviation subsequently has a negative impact on probability of cancer cells removal, decreasing the EP over time. The greatest effect on EP has been observed when both birth and death rates are being randomly modified and are anticorrelated. In addition, similar results are observed for TCP, where radiotherapy is included, indicating that increasing the standard deviation in kinetic rates at first enhances the probability of tumor eradication. But, it has a negative impact on treatment effectiveness if the TCP does not reach a probability of 1.

Population Dynamics in a Changing Environment: Random versus Periodic Switching

Environmental changes greatly influence the evolution of populations. Here, we study the dynamics of a population of two strains, one growing slightly faster than the other, competing for resources in a time-varying binary environment modeled by a carrying capacity switching either randomly or periodically between states of abundance and scarcity. The population dynamics is characterized by demographic noise (birth and death events) coupled to a varying environment. We elucidate the similarities and differences of the evolution subject to a stochastically and periodically varying environment. Importantly, the population size distribution is generally found to be broader under intermediate and fast random switching than under periodic variations, which results in markedly different asymptotic behaviors between the fixation probability of random and periodic switching. We also determine the detailed conditions under which the fixation probability of the slow strain is maximal.

Resist or perish: Fate of a microbial population subjected to a periodic presence of antimicrobial

The evolution of antimicrobial resistance can be strongly affected by variations of antimicrobial concentration. Here, we study the impact of periodic alternations of absence and presence of antimicrobial on resistance evolution in a microbial population, using a stochastic model that includes variations of both population composition and size, and fully incorporates stochastic population extinctions. We show that fast alternations of presence and absence of antimicrobial are inefficient to eradicate the microbial population and strongly favor the establishment of resistance, unless the antimicrobial increases enough the death rate. We further demonstrate that if the period of alternations is longer than a threshold value, the microbial population goes extinct upon the first addition of antimicrobial, if it is not rescued by resistance. We express the probability that the population is eradicated upon the first addition of antimicrobial, assuming rare mutations. Rescue by resistance can happen either if resistant mutants preexist, or if they appear after antimicrobial is added to the environment. Importantly, the latter case is fully prevented by perfect biostatic antimicrobials that completely stop division of sensitive microorganisms. By contrast, we show that the parameter regime where treatment is efficient is larger for biocidal drugs than for biostatic drugs. This sheds light on the respective merits of different antimicrobial modes of action.

Neutral competition in a deterministically changing environment: Revisiting continuum approaches

Environmental variation can play an important role in ecological competition by influencing the relative advantage between competing species. Here, we consider such effects by extending a classical, competitive Moran model to incorporate an environment that fluctuates periodically in time. We adapt methods from work on these classical models to investigate the effects of the magnitude and frequency of environmental fluctuations on two important population statistics: the probability of fixation and the mean time to fixation. In particular, we find that for small frequencies, the system behaves similar to a system with a constant fitness difference between the two species, and for large frequencies, the system behaves similar to a neutrally competitive model. Most interestingly, the system exhibits nontrivial behavior for intermediate frequencies. We conclude by showing that our results agree quite well with recent theoretical work on competitive models with a stochastically changing environment, and discuss how the methods we develop ease the mathematical analysis required to study such models.

Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments

Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels. The two populations undergo heritable, spontaneous phenotypic variations at different rates. The phenotypic state of each individual is represented by a continuous variable, and the phenotypic landscape of the populations evolves in time due to variations in the nutrient level. Exploiting the analytical tractability of our model, we study the long-time behaviour of the solutions to obtain a detailed mathematical depiction of the evolutionary dynamics. The results suggest that when nutrient levels undergo small and slow oscillations, it is evolutionarily more convenient to rarely undergo spontaneous phenotypic variations. Conversely, under relatively large and fast periodic oscillations in the nutrient levels, which bring about alternating cycles of starvation and nutrient abundance, higher rates of spontaneous phenotypic variations confer a competitive advantage. We discuss the implications of our results in the context of cancer metabolism.

Cooperation in Microbial Populations: Theory and Experimental Model Systems

Cooperative behavior, the costly provision of benefits to others, is common across all domains of life. This review article discusses cooperative behavior in the microbial world, mediated by the exchange of extracellular products called public goods. We focus on model species for which the production of a public good and the related growth disadvantage for the producing cells are well described. To unveil the biological and ecological factors promoting the emergence and stability of cooperative traits we take an interdisciplinary perspective and review insights gained from both mathematical models and well-controlled experimental model systems. Ecologically, we include crucial aspects of the microbial life cycle into our analysis and particularly consider population structures where ensembles of local communities (subpopulations) continuously emerge, grow, and disappear again. Biologically, we explicitly consider the synthesis and regulation of public good production. The discussion of the theoretical approaches includes general evolutionary concepts, population dynamics, and evolutionary game theory. As a specific but generic biological example, we consider populations of Pseudomonas putida and its regulation and use of pyoverdines, iron scavenging molecules, as public goods. The review closes with an overview on cooperation in spatially extended systems and also provides a critical assessment of the insights gained from the experimental and theoretical studies discussed. Current challenges and important new research opportunities are discussed, including the biochemical regulation of public goods, more realistic ecological scenarios resembling native environments, cell-to-cell signaling, and multispecies communities.

Comprehensive phase diagram for logistic populations in fluctuating environment

Population dynamics reflects an underlying birth-death process, where the rates associated with different events may depend on external environmental conditions and on the population density. A whole family of simple and popular deterministic models (such as logistic growth) supports a transcritical bifurcation point between an extinction phase and an active phase. Here we provide a comprehensive analysis of the phases of that system, taking into account both the endogenous demographic noise (random birth and death events) and the effect of environmental stochasticity that causes variations in birth and death rates. Three phases are identified: in the inactive phase the mean time to extinction T is independent of the carrying capacity N and scales logarithmically with the initial population size. In the power-law phase T∼N^{q}, and in the exponential phase T∼exp(αN). All three phases and the transitions between them are studied in detail. The breakdown of the continuum approximation is identified inside the power-law phase, and the accompanying changes in decline modes are analyzed. The applicability of the emerging picture to the analysis of ecological time series and to the management of conservation efforts is briefly discussed.

Classical stochastic systems with fast-switching environments: Reduced master equations, their interpretation, and limits of validity

We study classical Markovian stochastic systems with discrete states, coupled to randomly switching external environments. For fast environmental processes we derive reduced dynamics for the system itself, focusing on corrections to the adiabatic limit of infinite timescale separation. We show that this can lead to master equations with bursting events. Negative transition rates can result in the reduced master equation, leading to unphysical short-time behavior. However, the reduced master equation can describe stationary states better than a leading-order adiabatic calculation, similar to what is known for Kramers-Moyal expansions in the context of the Pawula theorem [R. F. Pawula, Phys. Rev. 162, 186 (1967)PHRVAO0031-899X10.1103/PhysRev.162.186; H. Risken and H. Vollmer, Z. Phys. B 35, 313 (1979)ZPBBDJ0340-224X10.1007/BF01319854]. We provide an interpretation of the reduced dynamics in discrete time and a criterion for the occurrence of negative rates for systems with two environmental states.

Model reduction methods for population dynamics with fast-switching environments: Reduced master equations, stochastic differential equations, and applications

We study stochastic population dynamics coupled to fast external environments and combine expansions in the inverse switching rate of the environment and a Kramers-Moyal expansion in the inverse size of the population. This leads to a series of approximation schemes, capturing both intrinsic and environmental noise. These methods provide a means of efficient simulation and we show how they can be used to obtain analytical results for the fluctuations of population dynamics in switching environments. We place the approximations in relation to existing work on piecewise-deterministic and piecewise-diffusive Markov processes. Finally, we demonstrate the accuracy and efficiency of these model-reduction methods in different research fields, including systems in biology and a model of crack propagation.

Phase Diagram for Logistic Systems under Bounded Stochasticity

Extinction is the ultimate absorbing state of any stochastic birth-death process; hence, the time to extinction is an important characteristic of any natural population. Here we consider logistic and logisticlike systems under the combined effect of demographic and bounded environmental stochasticity. Three phases are identified: an inactive phase where the mean time to extinction T increases logarithmically with the initial population size, an active phase where T grows exponentially with the carrying capacity N, and a temporal Griffiths phase, with a power-law relationship between T and N. The system supports an exponential phase only when the noise is bounded, in which case the continuum (diffusion) approximation breaks down within the Griffiths phase. This breakdown is associated with a crossover between qualitatively different survival statistics and decline modes. To study the power-law phase we present a new WKB scheme, which is applicable both in the diffusive and in the nondiffusive regime.

Exploration-exploitation tradeoffs dictate the optimal distributions of phenotypes for populations subject to fitness fluctuations

We study a minimal model for the growth of a phenotypically heterogeneous population of cells subject to a fluctuating environment in which they can replicate (by exploiting available resources) and modify their phenotype within a given landscape (thereby exploring novel configurations). The model displays an exploration-exploitation trade-off whose specifics depend on the statistics of the environment. Most notably, the phenotypic distribution corresponding to maximum population fitness (i.e., growth rate) requires a nonzero exploration rate when the magnitude of environmental fluctuations changes randomly over time, while a purely exploitative strategy turns out to be optimal in two-state environments, independently of the statistics of switching times. We obtain analytical insight into the limiting cases of very fast and very slow exploration rates by directly linking population growth to the features of the environment.

Quantifying the impact of a periodic presence of antimicrobial on resistance evolution in a homogeneous microbial population of fixed size

The evolution of antimicrobial resistance generally occurs in an environment where antimicrobial concentration is variable, which has dramatic consequences on the microorganisms' fitness landscape, and thus on the evolution of resistance. We investigate the effect of these time-varying patterns of selection within a stochastic model. We consider a homogeneous microbial population of fixed size subjected to periodic alternations of phases of absence and presence of an antimicrobial that stops growth. Combining analytical approaches and stochastic simulations, we quantify how the time necessary for fit resistant bacteria to take over the microbial population depends on the alternation period. We demonstrate that fast alternations strongly accelerate the evolution of resistance, reaching a plateau for sufficiently small periods. Furthermore, this acceleration is stronger in larger populations. For asymmetric alternations, featuring a different duration of the phases with and without antimicrobial, we shed light on the existence of a minimum for the time taken by the population to fully evolve resistance. The corresponding dramatic acceleration of the evolution of antimicrobial resistance likely occurs in realistic situations, and may have an important impact both in clinical and experimental situations.

Eco-evolutionary dynamics of a population with randomly switching carrying capacity

Environmental variability greatly influences the eco-evolutionary dynamics of a population, i.e. it affects how its size and composition evolve. Here, we study a well-mixed population of finite and fluctuating size whose growth is limited by a randomly switching carrying capacity. This models the environmental fluctuations between states of resources abundance and scarcity. The population consists of two strains, one growing slightly faster than the other, competing under two scenarios: one in which competition is solely for resources, and one in which the slow (cooperating) strain produces a public good (PG) that benefits also the fast (free-riding) strain. We investigate how the coupling of demographic and environmental (external) noise affects the population's eco-evolutionary dynamics. By analytical and computational means, we study the correlations between the population size and its composition, and discuss the social-dilemma-like 'eco-evolutionary game' characterizing the PG production. We determine in what conditions it is best to produce a PG; when cooperating is beneficial but outcompeted by free riding, and when the PG production is detrimental for cooperators. Within a linear noise approximation to populations of varying size, we also accurately analyse the coupled effects of demographic and environmental noise on the size distribution.

Theory of time-averaged neutral dynamics with environmental stochasticity

Competition is the main driver of population dynamics, which shapes the genetic composition of populations and the assembly of ecological communities. Neutral models assume that all the individuals are equivalent and that the dynamics is governed by demographic (shot) noise, with a steady state species abundance distribution (SAD) that reflects a mutation-extinction equilibrium. Recently, many empirical and theoretical studies emphasized the importance of environmental variations that affect coherently the relative fitness of entire populations. Here we consider two generic time-averaged neutral models; in both the relative fitness of each species fluctuates independently in time but its mean is zero. The first (model A) describes a system with local competition and linear fitness dependence of the birth-death rates, while in the second (model B) the competition is global and the fitness dependence is nonlinear. Due to this nonlinearity, model B admits a noise-induced stabilization mechanism that facilitates the invasion of new mutants. A self-consistent mean-field approach is used to reduce the multispecies problem to two-species dynamics, and the large-N asymptotics of the emerging set of Fokker-Planck equations is presented and solved. Our analytic expressions are shown to fit the SADs obtained from extensive Monte Carlo simulations and from numerical solutions of the corresponding master equations.

Survival behavior in the cyclic Lotka-Volterra model with a randomly switching reaction rate

We study the influence of a randomly switching reproduction-predation rate on the survival behavior of the nonspatial cyclic Lotka-Volterra model, also known as the zero-sum rock-paper-scissors game, used to metaphorically describe the cyclic competition between three species. In large and finite populations, demographic fluctuations (internal noise) drive two species to extinction in a finite time, while the species with the smallest reproduction-predation rate is the most likely to be the surviving one (law of the weakest). Here we model environmental (external) noise by assuming that the reproduction-predation rate of the strongest species (the fastest to reproduce and predate) in a given static environment randomly switches between two values corresponding to more and less favorable external conditions. We study the joint effect of environmental and demographic noise on the species survival probabilities and on the mean extinction time. In particular, we investigate whether the survival probabilities follow the law of the weakest and analyze their dependence on the external noise intensity and switching rate. Remarkably, when, on average, there is a finite number of switches prior to extinction, the survival probability of the predator of the species whose reaction rate switches typically varies nonmonotonically with the external noise intensity (with optimal survival about a critical noise strength). We also outline the relationship with the case where all reaction rates switch on markedly different time scales.