Phenotypic switching of populations of cells in a stochastic environment

Peter G. Hufton,Yen Ting Lin,Tobias Galla

DOI: https://doi.org/10.1088/1742-5468/aaa78e

2018-01-06

Abstract: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.

Populations and Evolution,Biological Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how cell populations optimize their survival and growth through phenotypic switching strategies in randomly changing environments. Specifically, the authors focus on how cells can maximize the growth rate of the population by adjusting their own phenotypic switching rates when the environmental state randomly switches between two discrete states. A mathematical framework is proposed in the paper to describe and analyze this dynamic process, and to explore how the optimal phenotypic switching strategy changes under different environmental switching rates. In addition, the paper also compares the optimal phenotypic responses in periodic and random environments, and finds that the optimal response in random environments is more likely to be heterogeneous, and the net growth rate is often higher under random environmental dynamics.

### Summary of Main Research Questions:
1. **Phenotypic Switching in Random Environments**: Study how cell populations adapt to environmental changes through phenotypic switching under the condition of random environmental state switching.
2. **Optimal Phenotypic Response**: Explore what the best phenotypic switching strategy for cell populations is under different environmental switching rates.
3. **Comparison between Periodic and Random Environments**: Compare the optimal phenotypic responses of cell populations in periodic and random environments, and the impact of these responses on the population growth rate.
4. **Theoretical Model and Application**: Propose a mathematical framework to describe and analyze the phenotypic switching dynamics of cell populations in random environments, and discuss the application of this model in ecology and treatment strategy design.

### Key Formulas:
- **Instantaneous Growth Rate**:
\[
\mu_t = \mu_A^\sigma \varphi_t + \mu_B^\sigma (1 - \varphi_t)
\]
where \(\mu_A^\sigma\) and \(\mu_B^\sigma\) represent the net growth rates of phenotype A and phenotype B respectively in environmental state \(\sigma\), and \(\varphi_t\) represents the proportion of phenotype A cells.

- **Average Growth Rate**:
\[
E(\mu) = P_0^* E(\mu | 0) + P_1^* E(\mu | 1)
\]
where \(P_0^*\) and \(P_1^*\) represent the steady - state probabilities of the environment being in state 0 and state 1 respectively, and \(E(\mu | \sigma)\) represents the average growth rate in environmental state \(\sigma\).

- **Fixed Point**:
\[
\varphi_\pm^\sigma = \frac{\Delta_\sigma - p - q \pm \sqrt{(\Delta_\sigma - p - q)^2 + 4q \Delta_\sigma}}{2 \Delta_\sigma}
\]
where \(\Delta_\sigma = \mu_A^\sigma - \mu_B^\sigma\), \(\varphi_+^\sigma\) is the stable fixed point, and \(\varphi_-^\sigma\) is the unstable fixed point.

Through these formulas and models, the authors can systematically analyze and predict cell population dynamics under different environmental conditions, providing a theoretical basis for understanding the survival strategies of microorganisms in complex environments.