Modelling chase-and-run migration in heterogeneous populations.

Cell migration is crucial for many physiological and pathological processes. During embryogenesis, neural crest cells undergo coordinated epithelial to mesenchymal transformations and migrate towards various forming organs. Here we develop a computational model to understand how mutual interactions between migrating neural crest cells (NCs) and the surrounding population of placode cells (PCs) generate coordinated migration. According to experimental findings, we implement a minimal set of hypotheses, based on a coupling between chemotactic movement of NCs in response to a placode-secreted chemoattractant (Sdf1) and repulsion induced from contact inhibition of locomotion (CIL), triggered by heterotypic NC-PC contacts. This basic set of assumptions is able to semi-quantitatively recapitulate experimental observations of the characteristic multispecies phenomenon of "chase-and-run", where the colony of NCs chases an evasive PC aggregate. The model further reproduces a number of in vitro manipulations, including full or partial disruption of NC chemotactic migration and selected mechanisms coordinating the CIL phenomenon. Finally, we provide various predictions based on altering other key components of the model mechanisms.

The impact of short- and long-range perception on population movements.

Navigation of cells and organisms is typically achieved by detecting and processing orienteering cues. Occasionally, a cue may be assessed over a much larger range than the individual's body size, as in visual scanning for landmarks. In this paper we formulate models that account for orientation in response to short- or long-range cue evaluation. Starting from an underlying random walk movement model, where a generic cue is evaluated locally or nonlocally to determine a preferred direction, we state corresponding macroscopic partial differential equations to describe population movements. Under certain approximations, these models reduce to well-known local and nonlocal biological transport equations, including those of Keller-Segel type. We consider a case-study application: "hilltopping" in Lepidoptera and other insects, a phenomenon in which populations accumulate at summits to improve encounter/mating rates. Nonlocal responses are shown to efficiently filter out the natural noisiness (or roughness) of typical landscapes and allow the population to preferentially accumulate at a subset of hilltopping locations, in line with field studies. Moreover, according to the timescale of movement, optimal responses may occur for different perceptual ranges.

A stochastic model of corneal epithelium maintenance and recovery following perturbation.

Various biological studies suggest that the corneal epithelium is maintained by active stem cells located in the limbus, the so-called limbal epithelial stem cell hypothesis. While numerous mathematical models have been developed to describe corneal epithelium wound healing, only a few have explored the process of corneal epithelium homeostasis. In this paper we present a purposefully simple stochastic mathematical model based on a chemical master equation approach, with the aim of clarifying the main factors involved in the maintenance process. Model analysis provides a set of constraints on the numbers of stem cells, division rates, and the number of division cycles required to maintain a healthy corneal epithelium. In addition, our stochastic analysis reveals noise reduction as the epithelium approaches its homeostatic state, indicating robustness to noise. Finally, recovery is analysed in the context of perturbation scenarios.

A Nonlocal Model for Contact Attraction and Repulsion in Heterogeneous Cell Populations.

Instructing others to move is fundamental for many populations, whether animal or cellular. In many instances, these commands are transmitted by contact, such that an instruction is relayed directly (e.g. by touch) from signaller to receiver: for cells, this can occur via receptor-ligand mediated interactions at their membranes, potentially at a distance if a cell extends long filopodia. Given that commands ranging from attractive to repelling can be transmitted over variable distances and between cells of the same (homotypic) or different (heterotypic) type, these mechanisms can clearly have a significant impact on the organisation of a tissue. In this paper, we extend a system of nonlocal partial differential equations (integrodifferential equations) to provide a general modelling framework to explore these processes, performing linear stability and numerical analyses to reveal its capacity to trigger the self-organisation of tissues. We demonstrate the potential of the framework via two illustrative applications: the contact-mediated dispersal of neural crest populations and the self-organisation of pigmentation patterns in zebrafish.

Mathematical modelling of glioma growth: the use of Diffusion Tensor Imaging (DTI) data to predict the anisotropic pathways of cancer invasion.

The nonuniform growth of certain forms of cancer can present significant complications for their treatment, a particularly acute problem in gliomas. A number of experimental results have suggested that invasion is facilitated by the directed movement of cells along the aligned neural fibre tracts that form a large component of the white matter. Diffusion tensor imaging (DTI) provides a window for visualising this anisotropy and gaining insight on the potential invasive pathways. In this paper we develop a mesoscopic model for glioma invasion based on the individual migration pathways of invading cells along the fibre tracts. Via scaling we obtain a macroscopic model that allows us to explore the overall growth of a tumour. To connect DTI data to parameters in the macroscopic model we assume that directional guidance along fibre tracts is described by a bimodal von Mises-Fisher distribution (a normal distribution on a unit sphere) and parametrised according to the directionality and degree of anisotropy in the diffusion tensors. We demonstrate the results in a simple model for glioma growth, exploiting both synthetic and genuine DTI datasets to reveal the potentially crucial role of anisotropic structure on invasion.

Intent matters: how flow and forms of information impact collective navigation.

The phenomenon of collective navigation has received considerable interest in recent years. A common line of thinking, backed by theoretical studies, is that collective navigation can improve navigation efficiency through the 'many-wrongs' principle, whereby individual error is reduced by comparing the headings of neighbours. When navigation takes place in a flowing environment, each individual's trajectory is influenced by drift. Consequently, a potential discrepancy emerges between an individual's intended heading and its actual heading. In this study, we develop a theoretical model to explore whether collective navigation benefits are altered according to the form of heading information transmitted between neighbours. Navigation based on each individual's intended heading is found to confer robust advantages across a wide spectrum of flows, via both a marked improvement in migration times and a capacity for a group to overcome flows unnavigable by solitary individuals. Navigation based on individual's actual headings is far less effective, only offering an improvement under highly favourable currents. For many currents, sharing actual heading information can even lead to journey times that exceed those of individual navigators.

How does cellular contact affect differentiation mediated pattern formation?

In this paper, we present a two-population continuous integro-differential model of cell differentiation, using a non-local term to describe the influence of the local environment on differentiation. We investigate three different versions of the model, with differentiation being cell autonomous, regulated via a community effect, or weakly dependent on the local cellular environment. We consider the spatial patterns that such different modes of differentiation produce, and investigate the formation of both stripes and spots by the model. We show that pattern formation only occurs when differentiation is regulated by a strong community effect. In this case, permanent spatial patterns only occur under a precise relationship between the parameters characterising cell dynamics, although transient patterns can persist for biologically relevant timescales when this condition is relaxed. In all cases, the long-lived patterns consist only of stripes, not spots.

The impact of rheotaxis and flow on the aggregation of organisms.

Dispersed populations often need to organize into groups. Chemical attractants provide one means of directing individuals into an aggregate, but whether these structures emerge can depend on various factors, such as there being a sufficiently large population or the response to the attractant being sufficiently sensitive. In an aquatic environment, fluid flow may heavily impact on population distribution and many aquatic organisms adopt a rheotaxis response when exposed to a current, orienting and swimming according to the flow field. Consequently, flow-induced transport could be substantially different for the population members and any aggregating signal they secrete. With the aim of investigating how flows and rheotaxis responses impact on an aquatic population's ability to form and maintain an aggregated profile, we develop and analyse a mathematical model that incorporates these factors. Through a systematic analysis into the effect of introducing rheotactic behaviour under various forms of environmental flow, we demonstrate that each of flow and rheotaxis can act beneficially or detrimentally on the ability to form and maintain a cluster. Synthesizing these findings, we test a hypothesis that density-dependent rheotaxis may be optimal for group formation and maintenance, in which individuals increase their rheotactic effort as they approach an aggregated state.

Modelling collective navigation via non-local communication.

Collective migration occurs throughout the animal kingdom, and demands both the interpretation of navigational cues and the perception of other individuals within the group. Navigational cues orient individuals towards a destination, while it has been demonstrated that communication between individuals enhances navigation through a reduction in orientation error. We develop a mathematical model of collective navigation that synthesizes navigational cues and perception of other individuals. Crucially, this approach incorporates uncertainty inherent to cue interpretation and perception in the decision making process, which can arise due to noisy environments. We demonstrate that collective navigation is more efficient than individual navigation, provided a threshold number of other individuals are perceptible. This benefit is even more pronounced in low navigation information environments. In navigation 'blindspots', where no information is available, navigation is enhanced through a relay that connects individuals in information-poor regions to individuals in information-rich regions. As an expository case study, we apply our framework to minke whale migration in the northeast Atlantic Ocean, and quantify the decrease in navigation ability due to anthropogenic noise pollution.

Towards an integrated experimental-theoretical approach for assessing the mechanistic basis of hair and feather morphogenesis.

In his seminal 1952 paper, 'The Chemical Basis of Morphogenesis', Alan Turing lays down a milestone in the application of theoretical approaches to understand complex biological processes. His deceptively simple demonstration that a system of reacting and diffusing chemicals could, under certain conditions, generate spatial patterning out of homogeneity provided an elegant solution to the problem of how one of nature's most intricate events occurs: the emergence of structure and form in the developing embryo. The molecular revolution that has taken place during the six decades following this landmark publication has now placed this generation of theoreticians and biologists in an excellent position to rigorously test the theory and, encouragingly, a number of systems have emerged that appear to conform to some of Turing's fundamental ideas. In this paper, we describe the history and more recent integration between experiment and theory in one of the key models for understanding pattern formation: the emergence of feathers and hair in the skins of birds and mammals.

Cellular automata and integrodifferential equation models for cell renewal in mosaic tissues.

Mosaic tissues are composed of two or more genetically distinct cell types. They occur naturally, and are also a useful experimental method for exploring tissue growth and maintenance. By marking the different cell types, one can study the patterns formed by proliferation, renewal and migration. Here, we present mathematical modelling suggesting that small changes in the type of interaction that cells have with their local cellular environment can lead to very different outcomes for the composition of mosaics. In cell renewal, proliferation of each cell type may depend linearly or nonlinearly on the local proportion of cells of that type, and these two possibilities produce very different patterns. We study two variations of a cellular automaton model based on simple rules for renewal. We then propose an integrodifferential equation model, and again consider two different forms of cellular interaction. The results of the continuous and cellular automata models are qualitatively the same, and we observe that changes in local environment interaction affect the dynamics for both. Furthermore, we demonstrate that the models reproduce some of the patterns seen in actual mosaic tissues. In particular, our results suggest that the differing patterns seen in organ parenchymas may be driven purely by the process of cell replacement under different interaction scenarios.

Modelling cell migration strategies in the extracellular matrix.

The extracellular matrix (ECM) is a highly organised structure with the capacity to direct cell migration through their tendency to follow matrix fibres, a process known as contact guidance. Amoeboid cell populations migrate in the ECM by making frequent shape changes and have minimal impact on its structure. Mesenchymal cells actively remodel the matrix to generate the space in which they can move. In this paper, these different types of movement are studied through simulation of a continuous transport model. It is shown that the process of contact guidance in a structured ECM can spatially organise cell populations. Furthermore, when combined with ECM remodelling, it can give rise to cellular pattern formation in the form of "cell-chains" or networks without additional environmental cues such as chemoattractants. These results are applied to a simple model for tumour invasion where it is shown that the interactions between invading cells and the ECM structure surrounding the tumour can have a profound impact on the pattern and rate of cell infiltration, including the formation of characteristic "fingering" patterns. The results are further discussed in the context of a variety of relevant processes during embryonic and adult stages.

A user's guide to PDE models for chemotaxis.

Mathematical modelling of chemotaxis (the movement of biological cells or organisms in response to chemical gradients) has developed into a large and diverse discipline, whose aspects include its mechanistic basis, the modelling of specific systems and the mathematical behaviour of the underlying equations. The Keller-Segel model of chemotaxis (Keller and Segel in J Theor Biol 26:399-415, 1970; 30:225-234, 1971) has provided a cornerstone for much of this work, its success being a consequence of its intuitive simplicity, analytical tractability and capacity to replicate key behaviour of chemotactic populations. One such property, the ability to display "auto-aggregation", has led to its prominence as a mechanism for self-organisation of biological systems. This phenomenon has been shown to lead to finite-time blow-up under certain formulations of the model, and a large body of work has been devoted to determining when blow-up occurs or whether globally existing solutions exist. In this paper, we explore in detail a number of variations of the original Keller-Segel model. We review their formulation from a biological perspective, contrast their patterning properties, summarise key results on their analytical properties and classify their solution form. We conclude with a brief discussion and expand on some of the outstanding issues revealed as a result of this work.

A chemotactic model for the advance and retreat of the primitive streak in avian development.

The formation of the primitive streak in early avian development marks the onset of gastrulation, during which large scale cell movement leads to a trilaminar blastoderm comprising prospective endodermal, mesodermal and ectodermal tissue. During streak formation a specialized group of cells first moves anteriorly as a coherent column, beginning from the posterior end of the prospective anterior-posterior axis (a process called progression), and then reverses course and returns to the most posterior point on the axis (a process called regression). To date little is known concerning the mechanisms controlling either progression or regression. Here we develop a model in which chemotaxis directs the cell movement and which is capable of reproducing the principal features connected with progression and regression of the primitive streak. We show that this model exhibits a number of experimentally-observed features of normal and abnormal streak development, and we propose a number of experimental tests which may serve to illuminate the mechanisms. This paper represents the first attempt to model the global features of primitive streak formation, and provides an initial stage in the development of a more biologically-realistic discrete cell model that will allow for variation of properties between cells and control over movement of individual cells.

Development and applications of a model for cellular response to multiple chemotactic cues.

The chemotactic response of a cell population to a single chemical species has been characterized experimentally for many cell types and has been extensively studied from a theoretical standpoint. However, cells frequently have multiple receptor types and can detect and respond chemotactically to more than one chemical. How these signals are integrated within the cell is not known. and we therefore adopt a macroscopic phenomenological approach to this problem. In this paper we derive and analyze chemotactic models based on partial differential (chemotaxis) equations for cell movement in response to multiple chemotactic cues. Our derivation generalizes the approach of Othmer and Stevens [29], who have recently developed a modeling framework for studying different chemotactic responses to a single chemical species. The importance of such a generalization is illustrated by the effect of multiple chemical cues on the chemotactic sensitivity and the spatial pattern of cell densities in several examples. We demonstrate that the model can generate the complex patterns observed on the skin of certain animal species and we indicate how the chemotactic response can be viewed as a form of positional indicator.

Pattern formation in a generalized chemotactic model.

Many models have been proposed for spatial pattern formation in embryology and analyzed for the standard case of zero-flux boundary conditions. However, relatively little attention has been paid to the role of boundary conditions on the form of the final pattern. Here we investigate numerically, the effect of nonstandard boundary conditions on a model pattern generator, which we choose to be of a cell-chemotactic type. We specifically focus on the role of boundary conditions and the effects of scale and aspect ratio, and study the spatiotemporal dynamics of pattern formation. We illustrate the properties of the model by application to the spatiotemporal sequence of skeletal development.

Stripe formation in juvenile Pomacanthus explained by a generalized turing mechanism with chemotaxis.

Current interest in pattern formation can be traced to a seminal paper by Turing, who demonstrated that a system of reacting and diffusing chemicals, called morphogens, can interact so as to produce stable nonuniform concentration patterns in space. Recently, a Turing model has been suggested to explain the development of pigmentation patterns on species of growing angelfish such as Pomacanthus semicirculatus, which exhibit readily observed changes in the number, size, and orientation of colored stripes during development of juvenile and adult stages, but the model fails to predict key features of the observations on stripe formation. Here we develop a generalized Turing model incorporating cell growth and movement, we analyze the effects of these processes on patterning, and we demonstrate that the model can explain important features of pattern formation in a growing system such as Pomacanthus. The applicability of classical Turing models to biological pattern formation is limited by virtue of the sensitivity of patterns to model parameters, but here we show that the incorporation of growth results in robustly generated patterns without strict parameter control. In the model, chemotaxis in response to gradients in a morphogen distribution leads to aggregation of one type of pigment cell into a striped spatial pattern.