Population switching under a time-varying environment

Escape from a metastable state in non-Markovian population dynamics

We study the long-time dynamics in non-Markovian single-population stochastic models, where one or more reactions are modeled as a stochastic process with a fat-tailed nonexponential distribution of waiting times, mimicking long-term memory. We focus on three prototypical examples: genetic switching, population establishment, and population extinction, all with nonexponential production rates. The system is studied in two regimes. In the first, the distribution of waiting times has a finite mean. Here, the system approaches a (quasi)stationary steady state at long times, and we develop a general Wentzel-Kramers-Brillouin approach for these non-Markovian systems. We derive explicit results for the mean population size and mean escape time from the metastable state of the stochastic dynamics. In this realm, we reveal that for sufficiently strong memory, a memory-induced (meta)stable state can emerge in the system. In the second regime, the waiting time distribution is assumed to have an infinite mean. Here, for bistable systems we find two distinct scaling regimes, separated by an exponentially long time which may strongly depend on the initial conditions of the system.

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.

First detection of threshold crossing events under intermittent sensing

The time taken by a random variable to cross a threshold for the first time, known as the first passage time, is of interest in many areas of sciences and engineering. Conventionally, there is an implicit assumption that the notional "sensor" monitoring the threshold crossing event is always active. In many realistic scenarios, the sensor monitoring the stochastic process works intermittently. Then, the relevant quantity of interest is the first detection time, which denotes the time when the sensor detects the random variable to be above the threshold for the first time. In this Letter, a birth-death process monitored by a random intermittent sensor is studied for which the first detection time distribution is obtained. In general, it is shown that the first detection time is related to and is obtainable from the first passage time distribution. Our analytical results display an excellent agreement with simulations. Furthermore, this framework is demonstrated in several applications-the susceptible infected susceptible compartmental and logistic models and birth-death processes with resetting. Finally, we solve the practically relevant problem of inferring the first passage time distribution from the first detection time.

Reconstructing an epigenetic landscape using a genetic pulling approach

Cells use genetic switches to shift between alternate stable gene expression states, e.g., to adapt to new environments or to follow a developmental pathway. Conceptually, these stable phenotypes can be considered as attractive states on an epigenetic landscape with phenotypic changes being transitions between states. Measuring these transitions is challenging because they are both very rare in the absence of appropriate signals and very fast. As such, it has proved difficult to experimentally map the epigenetic landscapes that are widely believed to underly developmental networks. Here, we introduce a nonequilibrium perturbation method to help reconstruct a regulatory network's epigenetic landscape. We derive the mathematical theory needed and then use the method on simulated data to reconstruct the landscapes. Our results show that with a relatively small number of perturbation experiments it is possible to recover an accurate representation of the true epigenetic landscape. We propose that our theory provides a general method by which epigenetic landscapes can be studied. Finally, our theory suggests that the total perturbation impulse required to induce a switch between metastable states is a fundamental quantity in developmental dynamics.