A comparative study between discrete and continuum models for the evolution of competing phenotype-structured cell populations in dynamical environments

Abstract

Deterministic continuum models formulated in terms of non-local partial differential equations for the evolutionary dynamics of populations structured by phenotypic traits have been used recently to address open questions concerning the adaptation of asexual species to periodically fluctuating environmental conditions. These deterministic continuum models are usually defined on the basis of population-scale phenomenological assumptions and cannot capture adaptive phenomena that are driven by stochastic variability in the evolutionary paths of single individuals. In light of these considerations, in this paper we develop a stochastic individual-based model for the coevolution between two competing phenotype-structured cell populations that are exposed to time-varying nutrient levels and undergo spontaneous, heritable phenotypic variations with different probabilities.Here, the evolution of every cell is described by a set of rules that result in a discrete-time branching random walk on the space of phenotypic states, and nutrient levels are governed by a difference equation in which a sink term models nutrient consumption by the cells. We formally show that the deterministic continuum counterpart of this model comprises a system of non-local partial differential equations for the cell population density functions coupled with an ordinary differential equation for the nutrient concentration. We compare the individual-based model and its continuum analogue, focussing on scenarios whereby the predictions of the two models differ. The results obtained clarify the conditions under which significant differences between the two models can emerge due to stochastic effects associated with small population levels. In particular, these differences arise in the presence of low probabilities of phenotypic variation, and become more apparent when the two populations are characterised by less fit initial mean phenotypes and smaller initial levels of phenotypic heterogeneity. The agreement between the two modelling approaches is also dependent on the initial proportions of the two populations.