Modeling Uniaxial Nonuniform Cell Proliferation

A Model for Cell Proliferation in a Developing Organism

In mathematical biology, there is a great deal of interest in producing continuum models by scaling discrete agent-based models governed by local stochastic rules. We discuss a particular example of this approach: a model for the proliferation of neural crest cells that can help us understand the development of Hirschprung's disease, a potentially-fatal condition in which the enteric nervous system of a new-born child does not extend all the way through the intestine and colon. Our starting point is a discrete-state, continuous-time Markov chain model proposed by Hywood et al. (2013a) for the location of the neural crest cells that make up the enteric nervous system. Hywood et al. (2013a) scaled their model to derive an approximate second order partial differential equation describing how the limiting expected number of neural crest cells evolve in space and time. In contrast, we exploit the relationship between the above-mentioned Markov chain model and the well-known Yule-Furry process to derive the exact form of the scaled version of the process. Furthermore, we provide expressions for other features of the domain agent occupancy process, such as the variance of the marginal occupancy at a particular site, the distribution of the number of agents that are yet to reach a given site and a stochastic description of the process itself.

Modelling uniaxial non-uniform yeast colony growth: Comparing an agent-based model and continuum approximations

Biological experiments have shown that yeast can be restricted to grow in a uniaxial direction, vertically upwards from an agar plate to form a colony. The growth occurs as a consequence of cell proliferation driven by a nutrient supply at the base of the colony, and the height of the colony has been observed to increase linearly with time. Within the colony the nutrient concentration is non-constant and yeast cells throughout the colony will therefore not have equal access to nutrient, resulting in non-uniform growth. In this work, an agent based model is developed to predict the microscopic spatial distribution of labelled cells within the colony when the probability of cell proliferation can vary in space and time. We also describe a method for determining the average trajectories or pathlines of labelled cells within a colony growing in a uniaxial direction, enabling us to connect the microscopic and macroscopic behaviours of the system. We present results for six cases, which involve different assumptions for the presence or absence of a quiescent region (where no cell proliferation occurs), the size of the proliferative region, and the spatial variation of proliferation rates within the proliferative region. These six cases are designed to provide qualitative insight into likely growth scenarios whilst remaining amenable to analysis. We compare our macroscopic results to experimental observations of uniaxial colony growth for two cases where only a fixed number of cells at the base of the colony can proliferate. The model predicts that the height of the colony will increase linearly with time in both these cases, which is consistent with experimental observations. However, our model shows how different functional forms for the spatial dependence of the proliferation rate can be distinguished by tracking the pathlines of cells at different positions in the colony. More generally, our methodology can be applied to other biological systems exhibiting uniaxial growth, providing a framework for classifying or determining regions of uniform and non-uniform growth.