No Arabic abstract
We study the dynamics of colonization of a territory by a stochastic population at low immigration pressure. We assume a sufficiently strong Allee effect that introduces, in deterministic theory, a large critical population size for colonization. At low immigration rates, the average pre-colonization population size is small thus invalidating the WKB approximation to the master equation. We circumvent this difficulty by deriving an exact zero-flux solution of the master equation and matching it with an approximate non-zero-flux solution of the pertinent Fokker-Planck equation in a small region around the critical population size. This procedure provides an accurate evaluation of the quasi-stationary probability distribution of population sizes in the pre-colonization state, and of the mean time to colonization, for a wide range of immigration rates. At sufficiently high immigration rates our results agree with WKB results obtained previously. At low immigration rates the results can be very different.
We study the extinction risk of a fragmented population residing on a network of patches coupled by migration, where the local patch dynamics include the Allee effect. We show that mixing between patches dramatically influences the populations viability. Slow migration is shown to always increase the populations global extinction risk compared to the isolated case. At fast migration, we demonstrate that synchrony between patches minimizes the populations extinction risk. Moreover, we discover a critical migration rate that maximizes the extinction risk of the population, and identify an early-warning signal when approaching this state. Our theoretical results are confirmed via the highly-efficient weighted ensemble method. Notably, our analysis can also be applied to studying switching in gene regulatory networks with multiple transcriptional states.
In recent years non-demographic variability has been shown to greatly affect dynamics of stochastic populations. For example, non-demographic noise in the form of a bursty reproduction process with an a-priori unknown burst size, or environmental variability in the form of time-varying reaction rates, have been separately found to dramatically impact the extinction risk of isolated populations. In this work we investigate the extinction risk of an isolated population under the combined influence of these two types of non-demographic variation. Using the so-called momentum-space WKB approach we arrive at a set of time-dependent Hamilton equations. In order to account for the explicit time dependence, we find the instanton of the time-perturbed Hamiltonian numerically, where analytical expressions are presented in particular limits using various perturbation techniques. We focus on two classes of time-varying environments: periodically-varying rates corresponding to seasonal effects, and a sudden decrease in the birth rate corresponding to a catastrophe. All our theoretical results are tested against numerical Monte Carlo simulations with time-dependent rates and also against a numerical solution of the corresponding time-dependent Hamilton equations.
Power-law-distributed species counts or clone counts arise in many biological settings such as multispecies cell populations, population genetics, and ecology. This empirical observation that the number of species $c_{k}$ represented by $k$ individuals scales as negative powers of $k$ is also supported by a series of theoretical birth-death-immigration (BDI) models that consistently predict many low-population species, a few intermediate-population species, and very high-population species. However, we show how a simple global population-dependent regulation in a neutral BDI model destroys the power law distributions. Simulation of the regulated BDI model shows a high probability of observing a high-population species that dominates the total population. Further analysis reveals that the origin of this breakdown is associated with the failure of a mean-field approximation for the expected species abundance distribution. We find an accurate estimate for the expected distribution $langle c_k rangle$ by mapping the problem to a lower-dimensional Moran process, allowing us to also straightforwardly calculate the covariances $langle c_k c_ell rangle$. Finally, we exploit the concepts associated with energy landscapes to explain the failure of the mean-field assumption by identifying a phase transition in the quasi-steady-state species counts triggered by a decreasing immigration rate.
Many unicellular organisms allocate their key proteins asymmetrically between the mother and daughter cells, especially in a stressed environment. A recent theoretical model is able to predict when the asymmetry in segregation of key proteins enhances the population fitness, extrapolating the solution at two limits where the segregation is perfectly asymmetric (asymmetry $a$ = 1) and when the asymmetry is small ($0 leq a ll 1$). We generalize the model by introducing stochasticity and use a transport equation to obtain a self-consistent equation for the population growth rate and the distribution of the amount of key proteins. We provide two ways of solving the self-consistent equation: numerically by updating the solution for the self-consistent equation iteratively and analytically by expanding moments of the distribution. With these more powerful tools, we can extend the previous model by Lin et al. to include stochasticity to the segregation asymmetry. We show the stochastic model is equivalent to the deterministic one with a modified effective asymmetry parameter ($a_{rm eff}$). We discuss the biological implication of our models and compare with other theoretical models.
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.