Rare events in stochastic populations under bursty reproduction

Abstract

Recently, a first step was made by the authors towards a systematic investigation of the effect of reaction-step-size noise - uncertainty in the step size of the reaction - on the dynamics of stochastic populations. This was done by investigating the effect of bursty influx on the switching dynamics of stochastic populations. Here we extend this formalism to account for bursty reproduction processes, and improve the accuracy of the formalism to include subleading-order corrections. Bursty reproduction appears in various contexts, where notable examples include bursty viral production from infected cells, and reproduction of mammals involving varying number of offspring. The main question we quantitatively address is how bursty reproduction affects the overall fate of the population. We consider two complementary scenarios: population extinction and population survival; in the former a population gets extinct after maintaining a long-lived metastable state, whereas in the latter a population proliferates despite undergoing a deterministic drift towards extinction. In both models reproduction occurs in bursts, sampled from an arbitrary distribution. In the extinction problem, we show that bursty reproduction broadens the quasi-stationary distribution of population sizes in the metastable state, which results in an exponential decrease of the mean time to extinction. In the survival problem, bursty reproduction yields an exponential increase in survival probability of the population. Close to the bifurcation limit our analytical results simplify considerably and are shown to depend solely on the mean and variance of the burst-size distribution. Our formalism is demonstrated on several realistic distributions which all compare well with numerical Monte-Carlo simulations.

arXiv comment19 pages, 7 figures

References in corpus (19)

• Analytical distributions for stochastic gene expression
• Colored extrinsic fluctuations and stochastic gene expression
• Rare Events Statistics in Reaction--Diffusion Systems
• Extinction of metastable stochastic populations
• Determining the stability of genetic switches: explicitly accounting for mRNA noise
• Extinction Rates for Fluctuation-Induced Metastabilities : A Real-Space WKB Approach
• Population extinction in a fluctuating environment
• Spectral theory of metastability and extinction in birth-death systems
• Extrinsic noise driven phenotype switching in a self-regulating gene
• Switching rates of multi-step reactions
• Cooperation dilemma in finite populations under fluctuating environments
• Spectral theory of metastability and extinction in a branching-annihilation reaction
• The dynamics of simple gene-network motifs subject to extrinsic fluctuations
• Noise driven unlimited population growth
• Spectral formulation and WKB approximation for rare-event statistics in reaction systems
• Fixation in Evolutionary Games under Non-Vanishing Selection
• Switching barrier scaling near bifurcation points for non-Gaussian noise
• Immigration-extinction dynamics of stochastic populations
• Effect of reaction step-size noise on the switching dynamics of stochastic populations

Cited by in corpus (6)

• A Minimal Model of Burst-Noise Induced Bistability
• A framework of nonequilibrium statistical mechanics. II. Coarse-graining
• Efficient Low-Order Approximation of First-Passage Time Distributions
• Population Extinction under Bursty Reproduction in a Time Modulated Environment
• Reducing the extinction risk of stochastic populations via non-demographic noise
• Fast reactions with non-interacting species in stochastic reaction networks