We study individual-based dynamics in finite populations, subject to randomly switching environmental conditions. These are inspired by models in which genes transition between on and off states, regulating underlying protein dynamics. Similarly, switches between environmental states are relevant in bacterial populations and in models of epidemic spread. Existing piecewise-deterministic Markov process approaches focus on the deterministic limit of the population dynamics while retaining the randomness of the switching. Here we go beyond this approximation and explicitly include effects of intrinsic stochasticity at the level of the linear-noise approximation. Specifically, we derive the stationary distributions of a number of model systems, in good agreement with simulations. This improves existing approaches which are limited to the regimes of fast and slow switching.

In biology phenotypic switching is a common bet-hedging strategy in the face of uncertain environmental conditions. Existing mathematical models often focus on periodically changing environments to determine the optimal phenotypic response. We focus on the case in which the environment switches randomly between discrete states. Starting from an individual-based model we derive stochastic differential equations to describe the dynamics, and obtain analytical expressions for the mean instantaneous growth rates based on the theory of piecewise-deterministic Markov processes. We show that optimal phenotypic responses are non-trivial for slow and intermediate environmental processes, and systematically compare the cases of periodic and random environments. The best response to random switching is more likely to be heterogeneity than in the case of deterministic periodic environments, net growth rates tend to be higher under stochastic environmental dynamics.The combined system of environment and population of cells can be interpreted as host-pathogen interaction, in which the host tries to choose environmental switching so as to minimise growth of the pathogen, and in which the pathogen employs a phenotypic switching optimised to increase its growth rate. We discuss the existence of Nash-like mutual best-response scenarios for such host-pathogen games.

Pluripotent embryonic stem cells are of paramount importance for biomedical sciences because of their innate ability for self-renewal and differentiation into all major cell lines. The fateful decision to exit or remain in the pluripotent state is regulated by complex genetic regulatory networks. The rapid growth of single-cell sequencing data has greatly stimulated applications of statistical and machine learning methods for inferring topologies of pluripotency regulating genetic networks. The inferred network topologies, however, often only encode Boolean information while remaining silent about the roles of dynamics and molecular stochasticity inherent in gene expression. Herein we develop a framework for systematically extending Boolean-level network topologies into higher resolution models of networks which explicitly account for the promoter architectures and gene state switching dynamics. We show the framework to be useful for disentangling the various contributions that gene switching, external signaling, and network topology make to the global heterogeneity and dynamics of transcription factor populations. We find the pluripotent state of the network to be a steady state which is robust to global variations of gene switching rates which we argue are a good proxy for epigenetic states of individual promoters. The temporal dynamics of exiting the pluripotent state, on the other hand, is significantly influenced by the rates of genetic switching which makes cells more responsive to changes in extracellular signals.

We study classical stochastic systems with discrete states, coupled to switching external environments. For fast environmental processes we derive reduced dynamics for the system itself, focusing on corrections to the adiabatic limit of infinite time scale separation. In some cases, this leads to master equations with negative transition 'rates' or bursting events. We devise a simulation algorithm in discrete time to unravel these master equations into sample paths, and provide an interpretation of bursting events. Focusing on stochastic population dynamics coupled to external environments, we discuss a series of approximation schemes combining expansions in the inverse switching rate of the environment, and a Kramers-Moyal expansion in the inverse size of the population. This places the different approximations in relation to existing work on piecewise-deterministic and piecewisediffusive Markov processes. We apply the model reduction methods to different examples including systems in biology and a model of crack propagation.

No abstract

The time of a stochastic process first passing through a boundary is important to many diverse applications. However, we can rarely compute the analytical distribution of these first-passage times. We develop an approximation to the first and second moments of a general first-passage time problem in the limit of large, but finite, populations using Kramers–Moyal expansion techniques. We demonstrate these results by application to a stochastic birth-death model for a population of cells in order to develop several approximations to the normal tissue complication probability (NTCP): a problem arising in the radiation treatment of cancers. We specifically allow for interaction between cells, via a nonlinear logistic growth model, and our approximations capture the effects of intrinsic noise on NTCP. We consider examples of NTCP in both a simple model of normal cells and in a model of normal and damaged cells. Our analytical approximation of NTCP could help optimise radiotherapy planning, for example by estimating the probability of complication-free tumour under different treatment protocols.

No abstract