Fast-mode elimination in stochastic metapopulation models

Modeling spatial evolution of multi-drug resistance under drug environmental gradients

Multi-drug combinations to treat bacterial populations are at the forefront of approaches for infection control and prevention of antibiotic resistance. Although the evolution of antibiotic resistance has been theoretically studied with mathematical population dynamics models, extensions to spatial dynamics remain rare in the literature, including in particular spatial evolution of multi-drug resistance. In this study, we propose a reaction-diffusion system that describes the multi-drug evolution of bacteria based on a drug-concentration rescaling approach. We show how the resistance to drugs in space, and the consequent adaptation of growth rate, is governed by a Price equation with diffusion, integrating features of drug interactions and collateral resistances or sensitivities to the drugs. We study spatial versions of the model where the distribution of drugs is homogeneous across space, and where the drugs vary environmentally in a piecewise-constant, linear and nonlinear manner. Although in many evolution models, per capita growth rate is a natural surrogate for fitness, in spatially-extended, potentially heterogeneous habitats, fitness is an emergent property that potentially reflects additional complexities, from boundary conditions to the specific spatial variation of growth rates. Applying concepts from perturbation theory and reaction-diffusion equations, we propose an analytical metric for characterization of average mutant fitness in the spatial system based on the principal eigenvalue of our linear problem, λ1. This enables an accurate translation from drug spatial gradients and mutant antibiotic susceptibility traits to the relative advantage of each mutant across the environment. Our approach allows one to predict the precise outcomes of selection among mutants over space, ultimately from comparing their λ1 values, which encode a critical interplay between growth functions, movement traits, habitat size and boundary conditions. Such mathematical understanding opens new avenues for multi-drug therapeutic optimization.

Modeling spatial evolution of multi-drug resistance under drug environmental gradients

Multi-drug combinations to treat bacterial populations are at the forefront of approaches for infection control and prevention of antibiotic resistance. Although the evolution of antibiotic resistance has been theoretically studied with mathematical population dynamics models, extensions to spatial dynamics remain rare in the literature, including in particular spatial evolution of multi-drug resistance. In this study, we propose a reaction-diffusion system that describes the multi-drug evolution of bacteria, based on a rescaling approach (Gjini and Wood, 2021). We show how the resistance to drugs in space, and the consequent adaptation of growth rate is governed by a Price equation with diffusion. The covariance terms in this equation integrate features of drug interactions and collateral resistances or sensitivities to the drugs. We study spatial versions of the model where the distribution of drugs is homogeneous across space, and where the drugs vary environmentally in a piecewise-constant, linear and nonlinear manner. Applying concepts from perturbation theory and reaction-diffusion equations, we propose an analytical characterization of average mutant fitness in the spatial system based on the principal eigenvalue of our linear problem. This enables an accurate translation from drug spatial gradients and mutant antibiotic susceptibility traits, to the relative advantage of each mutant across the environment. Such a mathematical understanding allows to predict the precise outcomes of selection over space, ultimately from the fundamental balance between growth and movement traits, and their diversity in a population.

Agent-based neutral competition in two-community networks

Competition between alternative states is an essential process in social and biological networks. Neutral competition can be represented by an unbiased random drift process in which the states of vertices (e.g., opinions, genotypes, or species) in a network are updated by repeatedly selecting two connected vertices. One of these vertices copies the state of the selected neighbor. Such updates are repeated until all vertices are in the same "consensus" state. There is no unique rule for selecting the vertex pair to be updated. Real-world processes comprise three limiting factors that can influence the selected edge and the direction of spread: (1) the rate at which a vertex sends a state to its neighbors, (2) the rate at which a state is received by a neighbor, and (3) the rate at which a state can be exchanged through a connecting edge. We investigate how these three limitations influence neutral competition in networks with two communities generated by a stochastic block model. By using Monte Carlo simulations, we show how the community structure and update rule determine the states' success probabilities and the time until a consensus is reached. We present a heterogeneous mean-field theory that agrees well with the Monte Carlo simulations. The effectiveness of the heterogeneous mean-field theory implies that quantitative predictions about the consensus are possible even if empirical data (e.g., from ecological fieldwork or observations of social interactions) do not allow a complete reconstruction of all edges in the network.

Evolution in alternating environments with tunable interlandscape correlations

Natural populations are often exposed to temporally varying environments. Evolutionary dynamics in varying environments have been extensively studied, although understanding the effects of varying selection pressures remains challenging. Here, we investigate how cycling between a pair of statistically related fitness landscapes affects the evolved fitness of an asexually reproducing population. We construct pairs of fitness landscapes that share global fitness features but are correlated with one another in a tunable way, resulting in landscape pairs with specific correlations. We find that switching between these landscape pairs, depending on the ruggedness of the landscape and the interlandscape correlation, can either increase or decrease steady-state fitness relative to evolution in single environments. In addition, we show that switching between rugged landscapes often selects for increased fitness in both landscapes, even in situations where the landscapes themselves are anticorrelated. We demonstrate that positively correlated landscapes often possess a shared maximum in both landscapes that allows the population to step through sub-optimal local fitness maxima that often trap single landscape evolution trajectories. Finally, we demonstrate that switching between anticorrelated paired landscapes leads to ergodic-like dynamics where each genotype is populated with nonzero probability, dramatically lowering the steady-state fitness in comparison to single landscape evolution.

Elimination of fast variables in stochastic nonlinear kinetics

A reduced chemical scheme involving a small number of variables is often sufficient to account for the deterministic evolution of the concentration of the main species contributing to a reaction. However, its predictions are questionable in small systems used, for example in fluorescence correlation spectroscopy (FCS) or in explosive systems involving strong nonlinearities such as autocatalytic steps. We make precise dynamical criteria defining the validity domain of the quasi-steady-state approximation and the elimination of a fast concentration in deterministic dynamics. Designing two different three-variable models converging toward the same two-variable model, we show that the variances and covariance of the fluctuations of the slow variables are not correctly predicted using the two-variable model, even in the limit of a large system size. The more striking weaknesses of the reduced scheme are figured out in mesoscaled systems containing a small number of molecules. The results of two stochastic approaches are compared and the shortcomings of the Langevin equations with respect to the master equation are pointed out. We conclude that the description of the fluctuations and their coupling with nonlinearities of deterministic dynamics escape reduced chemical schemes.

Tuning Spatial Profiles of Selection Pressure to Modulate the Evolution of Drug Resistance

Spatial heterogeneity plays an important role in the evolution of drug resistance. While recent studies have indicated that spatial gradients of selection pressure can accelerate resistance evolution, much less is known about evolution in more complex spatial profiles. Here we use a stochastic toy model of drug resistance to investigate how different spatial profiles of selection pressure impact the time to fixation of a resistant allele. Using mean first passage time calculations, we show that spatial heterogeneity accelerates resistance evolution when the rate of spatial migration is sufficiently large relative to mutation but slows fixation for small migration rates. Interestingly, there exists an intermediate regime-characterized by comparable rates of migration and mutation-in which the rate of fixation can be either accelerated or decelerated depending on the spatial profile, even when spatially averaged selection pressure remains constant. Finally, we demonstrate that optimal tuning of the spatial profile can dramatically slow the spread and fixation of resistant subpopulations, even in the absence of a fitness cost for resistance. Our results may lay the groundwork for optimized, spatially resolved drug dosing strategies for mitigating the effects of drug resistance.

Clonal dominance and transplantation dynamics in hematopoietic stem cell compartments

Hematopoietic stem cells in mammals are known to reside mostly in the bone marrow, but also transitively passage in small numbers in the blood. Experimental findings have suggested that they exist in a dynamic equilibrium, continuously migrating between these two compartments. Here we construct an individual-based mathematical model of this process, which is parametrised using existing empirical findings from mice. This approach allows us to quantify the amount of migration between the bone marrow niches and the peripheral blood. We use this model to investigate clonal hematopoiesis, which is a significant risk factor for hematologic cancers. We also analyse the engraftment of donor stem cells into non-conditioned and conditioned hosts, quantifying the impact of different treatment scenarios. The simplicity of the model permits a thorough mathematical analysis, providing deeper insights into the dynamics of both the model and of the real-world system. We predict the time taken for mutant clones to expand within a host, as well as chimerism levels that can be expected following transplantation therapy, and the probability that a preconditioned host is reconstituted by donor cells.

Clonal hematopoiesis—where mature myeloid cells in the blood deriving from a single stem cell are over-represented—is a major risk factor for overt hematologic malignancies. To quantify how likely this phenomena is, we combine existing observations with a novel stochastic model and extensive mathematical analysis. This approach allows us to observe the hidden dynamics of the hematopoietic system. We conclude that for a clone to be detectable within the lifetime of a mouse, it requires a selective advantage. I.e. the clonal expansion cannot be explained by neutral drift alone. Furthermore, we use our model to describe the dynamics of hematopoiesis after stem cell transplantation. In agreement with earlier findings, we observe that niche-space saturation decreases engraftment efficiency. We further discuss the implications of our findings for human hematopoiesis where the quantity and role of stem cells is frequently debated.

Mapping of the stochastic Lotka-Volterra model to models of population genetics and game theory

The relationship between the M-species stochastic Lotka-Volterra competition (SLVC) model and the M-allele Moran model of population genetics is explored via timescale separation arguments. When selection for species is weak and the population size is large but finite, precise conditions are determined for the stochastic dynamics of the SLVC model to be mappable to the neutral Moran model, the Moran model with frequency-independent selection, and the Moran model with frequency-dependent selection (equivalently a game-theoretic formulation of the Moran model). We demonstrate how these mappings can be used to calculate extinction probabilities and the times until a species' extinction in the SLVC model.

Heterogeneous network epidemics: real-time growth, variance and extinction of infection

Recent years have seen a large amount of interest in epidemics on networks as a way of representing the complex structure of contacts capable of spreading infections through the modern human population. The configuration model is a popular choice in theoretical studies since it combines the ability to specify the distribution of the number of contacts (degree) with analytical tractability. Here we consider the early real-time behaviour of the Markovian SIR epidemic model on a configuration model network using a multitype branching process. We find closed-form analytic expressions for the mean and variance of the number of infectious individuals as a function of time and the degree of the initially infected individual(s), and write down a system of differential equations for the probability of extinction by time t that are numerically fast compared to Monte Carlo simulation. We show that these quantities are all sensitive to the degree distribution—in particular we confirm that the mean prevalence of infection depends on the first two moments of the degree distribution and the variance in prevalence depends on the first three moments of the degree distribution. In contrast to most existing analytic approaches, the accuracy of these results does not depend on having a large number of infectious individuals, meaning that in the large population limit they would be asymptotically exact even for one initial infectious individual.

Stochastic epidemic dynamics on extremely heterogeneous networks

Networks of contacts capable of spreading infectious diseases are often observed to be highly heterogeneous, with the majority of individuals having fewer contacts than the mean, and a significant minority having relatively very many contacts. We derive a two-dimensional diffusion model for the full temporal behavior of the stochastic susceptible-infectious-recovered (SIR) model on such a network, by making use of a time-scale separation in the deterministic limit of the dynamics. This low-dimensional process is an accurate approximation to the full model in the limit of large populations, even for cases when the time-scale separation is not too pronounced, provided the maximum degree is not of the order of the population size.

Demographic noise can reverse the direction of deterministic selection

Deterministic evolutionary theory robustly predicts that populations displaying altruistic behaviors will be driven to extinction by mutant cheats that absorb common benefits but do not themselves contribute. Here we show that when demographic stochasticity is accounted for, selection can in fact act in the reverse direction to that predicted deterministically, instead favoring cooperative behaviors that appreciably increase the carrying capacity of the population. Populations that exist in larger numbers experience a selective advantage by being more stochastically robust to invasions than smaller populations, and this advantage can persist even in the presence of reproductive costs. We investigate this general effect in the specific context of public goods production and find conditions for stochastic selection reversal leading to the success of public good producers. This insight, developed here analytically, is missed by the deterministic analysis as well as by standard game theoretic models that enforce a fixed population size. The effect is found to be amplified by space; in this scenario we find that selection reversal occurs within biologically reasonable parameter regimes for microbial populations. Beyond the public good problem, we formulate a general mathematical framework for models that may exhibit stochastic selection reversal. In this context, we describe a stochastic analog to [Formula: see text] theory, by which small populations can evolve to higher densities in the absence of disturbance.

Genetic Toggle Switch in the Absence of Cooperative Binding: Exact Results

We present an analytical treatment of a genetic switch model consisting of two mutually inhibiting genes operating without cooperative binding of the corresponding transcription factors. Previous studies have numerically shown that these systems can exhibit bimodal dynamics without possessing two stable fixed points at the deterministic level. We analytically show that bimodality is induced by the noise and find the critical repression strength that controls a transition between the bimodal and nonbimodal regimes. We also identify characteristic polynomial scaling laws of the mean switching time between bimodal states. These results, independent of the model under study, reveal essential differences between these systems and systems with cooperative binding, where there is no critical threshold for bimodality and the mean switching time scales exponentially with the system size.

Fixation properties of subdivided populations with balancing selection

In subdivided populations, migration acts together with selection and genetic drift and determines their evolution. Building upon a recently proposed method, which hinges on the emergence of a time scale separation between local and global dynamics, we study the fixation properties of subdivided populations in the presence of balancing selection. The approximation implied by the method is accurate when the effective selection strength is small and the number of subpopulations is large. In particular, it predicts a phase transition between species coexistence and biodiversity loss in the infinite-size limit and, in finite populations, a nonmonotonic dependence of the mean fixation time on the migration rate. In order to investigate the fixation properties of the subdivided population for stronger selection, we introduce an effective coarser description of the dynamics in terms of a voter model with intermediate states, which highlights the basic mechanisms driving the evolutionary process.

Stationary solutions for metapopulation Moran models with mutation and selection

We construct an individual-based metapopulation model of population genetics featuring migration, mutation, selection, and genetic drift. In the case of a single "island," the model reduces to the Moran model. Using the diffusion approximation and time-scale separation arguments, an effective one-variable description of the model is developed. The effective description bears similarities to the well-mixed Moran model with effective parameters that depend on the network structure and island sizes, and it is amenable to analysis. Predictions from the reduced theory match the results from stochastic simulations across a range of parameters. The nature of the fast-variable elimination technique we adopt is further studied by applying it to a linear system, where it provides a precise description of the slow dynamics in the limit of large time-scale separation.

Models of genetic drift as limiting forms of the Lotka-Volterra competition model

The relationship between the Moran model and stochastic Lotka-Volterra competition (SLVC) model is explored via time scale separation arguments. For neutral systems the two are found to be equivalent at long times. For systems with selective pressure, their behavior differs. It is argued that the SLVC is preferable to the Moran model since in the SLVC population size is regulated by competition, rather than arbitrarily fixed as in the Moran model. As a consequence, ambiguities found in the Moran model associated with the introduction of more complex processes, such as selection, are avoided.

Population genetics on islands connected by an arbitrary network: an analytic approach

We analyse a model consisting of a population of individuals which is subdivided into a finite set of demes, each of which has a fixed but differing number of individuals. The individuals can reproduce, die and migrate between the demes according to an arbitrary migration network. They are haploid, with two alleles present in the population; frequency-independent selection is also incorporated, where the strength and direction of selection can vary from deme to deme. The system is formulated as an individual-based model and the diffusion approximation systematically applied to express it as a set of nonlinear coupled stochastic differential equations. These can be made amenable to analysis through the elimination of fast-time variables. The resulting reduced model is analysed in a number of situations, including migration-selection balance leading to a polymorphic equilibrium of the two alleles and an illustration of how the subdivision of the population can lead to non-trivial behaviour in the case where the network is a simple hub. The method we develop is systematic, may be applied to any network, and agrees well with the results of simulations in all cases studied and across a wide range of parameter values.