Phenotypic switching mechanisms determine the structure of cell migration into extracellular matrix under the 'go-or-grow' hypothesis.

A fundamental feature of collective cell migration is phenotypic heterogeneity which, for example, influences tumour progression and relapse. While current mathematical models often consider discrete phenotypic structuring of the cell population, in-line with the 'go-or-grow' hypothesis (Hatzikirou et al., 2012; Stepien et al., 2018), they regularly overlook the role that the environment may play in determining the cells' phenotype during migration. Comparing a previously studied volume-filling model for a homogeneous population of generalist cells that can proliferate, move and degrade extracellular matrix (ECM) (Crossley et al., 2023) to a novel model for a heterogeneous population comprising two distinct sub-populations of specialist cells that can either move and degrade ECM or proliferate, this study explores how different hypothetical phenotypic switching mechanisms affect the speed and structure of the invading cell populations. Through a continuum model derived from its individual-based counterpart, insights into the influence of the ECM and the impact of phenotypic switching on migrating cell populations emerge. Notably, specialist cell populations that cannot switch phenotype show reduced invasiveness compared to generalist cell populations, while implementing different forms of switching significantly alters the structure of migrating cell fronts. This key result suggests that the structure of an invading cell population could be used to infer the underlying mechanisms governing phenotypic switching.

A biased random walk approach for modeling the collective chemotaxis of neural crest cells.

Collective cell migration is a multicellular phenomenon that arises in various biological contexts, including cancer and embryo development. 'Collectiveness' can be promoted by cell-cell interactions such as co-attraction and contact inhibition of locomotion. These mechanisms act on cell polarity, pivotal for directed cell motility, through influencing the intracellular dynamics of small GTPases such as Rac1. To model these dynamics we introduce a biased random walk model, where the bias depends on the internal state of Rac1, and the Rac1 state is influenced by cell-cell interactions and chemoattractive cues. In an extensive simulation study we demonstrate and explain the scope and applicability of the introduced model in various scenarios. The use of a biased random walk model allows for the derivation of a corresponding partial differential equation for the cell density while still maintaining a certain level of intracellular detail from the individual based setting.

Modelling microtube driven invasion of glioma.

Malignant gliomas are notoriously invasive, a major impediment against their successful treatment. This invasive growth has motivated the use of predictive partial differential equation models, formulated at varying levels of detail, and including (i) "proliferation-infiltration" models, (ii) "go-or-grow" models, and (iii) anisotropic diffusion models. Often, these models use macroscopic observations of a diffuse tumour interface to motivate a phenomenological description of invasion, rather than performing a detailed and mechanistic modelling of glioma cell invasion processes. Here we close this gap. Based on experiments that support an important role played by long cellular protrusions, termed tumour microtubes, we formulate a new model for microtube-driven glioma invasion. In particular, we model a population of tumour cells that extend tissue-infiltrating microtubes. Mitosis leads to new nuclei that migrate along the microtubes and settle elsewhere. A combination of steady state analysis and numerical simulation is employed to show that the model can predict an expanding tumour, with travelling wave solutions led by microtube dynamics. A sequence of scaling arguments allows us reduce the detailed model into simpler formulations, including models falling into each of the general classes (i), (ii), and (iii) above. This analysis allows us to clearly identify the assumptions under which these various models can be a posteriori justified in the context of microtube-driven glioma invasion. Numerical simulations are used to compare the various model classes and we discuss their advantages and disadvantages.

Feather arrays are patterned by interacting signalling and cell density waves.

Feathers are arranged in a precise pattern in avian skin. They first arise during development in a row along the dorsal midline, with rows of new feather buds added sequentially in a spreading wave. We show that the patterning of feathers relies on coupled fibroblast growth factor (FGF) and bone morphogenetic protein (BMP) signalling together with mesenchymal cell movement, acting in a coordinated reaction-diffusion-taxis system. This periodic patterning system is partly mechanochemical, with mechanical-chemical integration occurring through a positive feedback loop centred on FGF20, which induces cell aggregation, mechanically compressing the epidermis to rapidly intensify FGF20 expression. The travelling wave of feather formation is imposed by expanding expression of Ectodysplasin A (EDA), which initiates the expression of FGF20. The EDA wave spreads across a mesenchymal cell density gradient, triggering pattern formation by lowering the threshold of mesenchymal cells required to begin to form a feather bud. These waves, and the precise arrangement of feather primordia, are lost in the flightless emu and ostrich, though via different developmental routes. The ostrich retains the tract arrangement characteristic of birds in general but lays down feather primordia without a wave, akin to the process of hair follicle formation in mammalian embryos. The embryonic emu skin lacks sufficient cells to enact feather formation, causing failure of tract formation, and instead the entire skin gains feather primordia through a later process. This work shows that a reaction-diffusion-taxis system, integrated with mechanical processes, generates the feather array. In flighted birds, the key role of the EDA/Ectodysplasin A receptor (EDAR) pathway in vertebrate skin patterning has been recast to activate this process in a quasi-1-dimensional manner, imposing highly ordered pattern formation.

Mathematical models for chemotaxis and their applications in self-organisation phenomena.

Chemotaxis is a fundamental guidance mechanism of cells and organisms, responsible for attracting microbes to food, embryonic cells into developing tissues, immune cells to infection sites, animals towards potential mates, and mathematicians into biology. The Patlak-Keller-Segel (PKS) system forms part of the bedrock of mathematical biology, a go-to-choice for modellers and analysts alike. For the former it is simple yet recapitulates numerous phenomena; the latter are attracted to these rich dynamics. Here I review the adoption of PKS systems when explaining self-organisation processes. I consider their foundation, returning to the initial efforts of Patlak and Keller and Segel, and briefly describe their patterning properties. Applications of PKS systems are considered in their diverse areas, including microbiology, development, immunology, cancer, ecology and crime. In each case a historical perspective is provided on the evidence for chemotactic behaviour, followed by a review of modelling efforts; a compendium of the models is included as an Appendix. Finally, a half-serious/half-tongue-in-cheek model is developed to explain how cliques form in academia. Assumptions in which scholars alter their research line according to available problems leads to clustering of academics and the formation of "hot" research topics.

Aggregation of biological particles under radial directional guidance.

Many biological environments display an almost radially-symmetric structure, allowing proteins, cells or animals to move in an oriented fashion. Motivated by specific examples of cell movement in tissues, pigment protein movement in pigment cells and animal movement near watering holes, we consider a class of radially-symmetric anisotropic diffusion problems, which we call the star problem. The corresponding diffusion tensor D(x) is radially symmetric with isotropic diffusion at the origin. We show that the anisotropic geometry of the environment can lead to strong aggregations and blow-up at the origin. We classify the nature of aggregation and blow-up solutions and provide corresponding numerical simulations. A surprising element of this strong aggregation mechanism is that it is entirely based on geometry and does not derive from chemotaxis, adhesion or other well known aggregating mechanisms. We use these aggregate solutions to discuss the process of pigmentation changes in animals, cancer invasion in an oriented fibrous habitat (such as collagen fibres), and sheep distributions around watering holes.

Moments of von Mises and Fisher distributions and applications.

The von Mises and Fisher distributions are spherical analogues to the Normal distribution on the unit circle and unit sphere, respectively. The computation of their moments, and in particular the second moment, usually involves solving tedious trigonometric integrals. Here we present a new method to compute the moments of spherical distributions, based on the divergence theorem. This method allows a clear derivation of the second moments and can be easily generalized to higher dimensions. In particular we note that, to our knowledge, the variance-covariance matrix of the three dimensional Fisher distribution has not previously been explicitly computed. While the emphasis of this paper lies in calculating the moments of spherical distributions, their usefulness is motivated by their relationship to population statistics in animal/cell movement models and demonstrated in applications to the modelling of sea turtle navigation, wolf movement and brain tumour growth.

Spatio-temporal Models of Lymphangiogenesis in Wound Healing.

Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis), very few have been proposed for the regeneration of the lymphatic network. Lymphangiogenesis is a markedly different process from angiogenesis, occurring at different times and in response to different chemical stimuli. Two main hypotheses have been proposed: (1) lymphatic capillaries sprout from existing interrupted ones at the edge of the wound in analogy to the blood angiogenesis case and (2) lymphatic endothelial cells first pool in the wound region following the lymph flow and then, once sufficiently populated, start to form a network. Here, we present two PDE models describing lymphangiogenesis according to these two different hypotheses. Further, we include the effect of advection due to interstitial flow and lymph flow coming from open capillaries. The variables represent different cell densities and growth factor concentrations, and where possible the parameters are estimated from biological data. The models are then solved numerically and the results are compared with the available biological literature.

A mathematical model for lymphangiogenesis in normal and diabetic wounds.

Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis) very few have been proposed for the regeneration of the lymphatic network. Moreover, lymphangiogenesis is markedly distinct from angiogenesis, occurring at different times and in a different manner. Here a model of five ordinary differential equations is presented to describe the formation of lymphatic capillaries following a skin wound. The variables represent different cell densities and growth factor concentrations, and where possible the parameters are estimated from experimental and clinical data. The system is then solved numerically and the results are compared with the available biological literature. Finally, a parameter sensitivity analysis of the model is taken as a starting point for suggesting new therapeutic approaches targeting the enhancement of lymphangiogenesis in diabetic wounds. The work provides a deeper understanding of the phenomenon in question, clarifying the main factors involved. In particular, the balance between TGF-ß and VEGF levels, rather than their absolute values, is identified as crucial to effective lymphangiogenesis. In addition, the results indicate lowering the macrophage-mediated activation of TGF-ß and increasing the basal lymphatic endothelial cell growth rate, inter alia, as potential treatments. It is hoped the findings of this paper may be considered in the development of future experiments investigating novel lymphangiogenic therapies.

The impact of adhesion on cellular invasion processes in cancer and development.

In this paper we consider a simple continuous model to describe cell invasion, incorporating the effects of both cell-cell adhesion and cell-matrix adhesion, along with cell growth and proteolysis by cells of the surrounding extracellular matrix (ECM). We demonstrate that the model is capable of supporting both noninvasive and invasive tumour growth according to the relative strength of cell-cell to cell-matrix adhesion. Specifically, for sufficiently strong cell-matrix adhesion and/or sufficiently weak cell-cell adhesion, degradation of the surrounding ECM accompanied by cell-matrix adhesion pulls the cells into the surrounding ECM. We investigate the criticality of matrix heterogeneity on shaping invasion, demonstrating that a highly heterogeneous ECM can result in a "fingering" of the invasive front, echoing observations in real-life invasion processes ranging from malignant tumour growth to neural crest migration during embryonic development.

Continuous models for cell migration in tissues and applications to cell sorting via differential chemotaxis.

Chemotaxis, the guided migration of cells in response to chemical gradients, is vital to a wide variety of biological processes, including patterning of the slime mold Dictyostelium, embryonic morphogenesis, wound healing, and tumor invasion. Continuous models of chemotaxis have been developed to describe many such systems, yet few have considered the movements within a heterogeneous tissue composed of multiple subpopulations. In this paper, a partial differential equation (PDE) model is developed to describe a tissue formed from two distinct chemotactic populations. For a "crowded" (negligible extracellular space) tissue, it is demonstrated that the model reduces to a simpler one-species system while for an "uncrowded" tissue, it captures both movement of the entire tissue (via cells attaching to/migrating within an extracellular substrate) and the within-tissue rearrangements of the separate cellular subpopulations. The model is applied to explore the sorting of a heterogeneous tissue, where it is shown that differential-chemotaxis not only generates classical sorting patterns previously seen via differential-adhesion, but also demonstrates new classes of behavior. These new phenomena include temporal dynamics consisting of a traveling wave composed of spatially sorted subpopulations reminiscent of Dictyostelium slugs.

A continuum approach to modelling cell-cell adhesion.

Cells adhere to each other through the binding of cell adhesion molecules at the cell surface. This process, known as cell-cell adhesion, is fundamental in many areas of biology, including early embryo development, tissue homeostasis and tumour growth. In this paper we develop a new continuous mathematical model of this phenomenon by considering the movement of cells in response to the adhesive forces generated through binding. We demonstrate that our model predicts the aggregation behaviour of a disassociated adhesive cell population. Further, when the model is extended to represent the interactions between multiple populations, we demonstrate that it is capable of replicating the different types of cell sorting behaviour observed experimentally. The resulting pattern formation is a direct consequence of the relative strengths of self-population and cross-population adhesive bonds in the model. While cell sorting behaviour has been captured previously with discrete approaches, it has not, until now, been observed with a fully continuous model.

From a discrete to a continuous model of biological cell movement.

The process by which one may take a discrete model of a biophysical process and construct a continuous model based upon it is of mathematical interest as well as being of practical use. In this work, we take the extended Potts model applied to biological cell movement to its continuous limit. Beginning with a single cell moving in one dimension on a lattice and obeying Potts model rules of movement we develop an expression for the diffusion coefficient of a collection of noninteracting cells which depends explicitly on the Potts model parameters. We show how this coefficient varies when the Potts parameters for cell membrane elasticity and cell-medium adhesion are varied, and perform computer simulations which support our theoretical result. We explain the relationship between the probability of occupancy of lattice points and the density profile in the continuous limit, and extend our analysis by including interactions between the cells. In so doing we are able to develop a set of coupled ordinary differential equations showing the evolution of a density profile in the presence of significant cell-cell adhesion, and show how increases in the strength of this adhesion modulates diffusion. In so doing we develop some insights into how continuous models of physical systems can be based upon discrete models which describe the same system.