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.

An individual-based model to explore the impact of psychological stress on immune infiltration into tumour spheroids.

In recentexperiments on co-culture between breast tumour spheroids and activated immune cells, it was observed that the introduction of the stress hormone cortisol resulted in a decreased immune cell infiltration into the spheroids. Moreover, the presence of cortisol deregulated the normal levels of the pro- and anti-inflammatory cytokines IFN-and IL-10. We present an individual-based model to explore the interaction dynamics between tumour and immune cells under psychological stress conditions.

Evolutionary dynamics of glucose-deprived cancer cells: insights from experimentally informed mathematical modelling.

Glucose is a primary energy source for cancer cells. Several lines of evidence support the idea that monocarboxylate transporters, such as MCT1, elicit metabolic reprogramming of cancer cells in glucose-poor environments, allowing them to re-use lactate, a by-product of glucose metabolism, as an alternative energy source with serious consequences for disease progression. We employ a synergistic experimental and mathematical modelling approach to explore the evolutionary processes at the root of cancer cell adaptation to glucose deprivation, with particular focus on the mechanisms underlying the increase in MCT1 expression observed in glucose-deprived aggressive cancer cells.

Discrete and continuum models for the coevolutionary dynamics between CD8+ cytotoxic T lymphocytes and tumour cells.

We present an individual-based model for the coevolutionary dynamics between CD8+ cytotoxic T lymphocytes (CTLs) and tumour cells. In this model, every cell is viewed as an individual agent whose phenotypic state is modelled by a discrete variable. For tumour cells, this variable represents a parameterization of the antigen expression profiles, while for CTLs it represents a parameterization of the target antigens of T-cell receptors (TCRs).

The Impact of Phenotypic Heterogeneity on Chemotactic Self-Organisation.

The capacity to aggregate through chemosensitive movement forms a paradigm of self-organisation, with examples spanning cellular and animal systems. A basic mechanism assumes a phenotypically homogeneous population that secretes its own attractant, with the well known system introduced more than five decades ago by Keller and Segel proving resolutely popular in modelling studies. The typical assumption of population phenotypic homogeneity, however, often lies at odds with the heterogeneity of natural systems, where populations may comprise distinct phenotypes that vary according to their chemotactic ability, attractant secretion, etc.

A Hybrid Discrete-Continuum Modelling Approach to Explore the Impact of T-Cell Infiltration on Anti-tumour Immune Response.

We present a spatial hybrid discrete-continuum modelling framework for the interaction dynamics between tumour cells and cytotoxic T cells, which play a pivotal role in the immune response against tumours. In this framework, tumour cells and T cells are modelled as individual agents while chemokines that drive the chemotactic movement of T cells towards the tumour are modelled as a continuum. We formally derive the continuum counterpart of this model, which is given by a coupled system that comprises an integro-differential equation for the density of tumour cells, a partial differential equation for the density of T cells and a partial differential equation for the concentration of chemokines.

Travelling-wave analysis of a model of tumour invasion with degenerate, cross-dependent diffusion.

In this paper, we carry out a travelling-wave analysis of a model of tumour invasion with degenerate, cross-dependent diffusion. We consider two types of invasive fronts of tumour tissue into extracellular matrix (ECM), which represents healthy tissue. These types differ according to whether the density of ECM far ahead of the wave front is maximal or not.

A mathematical model to study the impact of intra-tumour heterogeneity on anti-tumour CD8 T cell immune response.

Intra-tumour heterogeneity (ITH) has a strong impact on the efficacy of the immune response against solid tumours. The number of sub-populations of cancer cells expressing different antigens and the percentage of immunogenic cells (i.e.

Evolutionary dynamics in an SI epidemic model with phenotype-structured susceptible compartment.

We present an SI epidemic model whereby a continuous structuring variable captures variability in proliferative potential and resistance to infection among susceptible individuals. The occurrence of heritable, spontaneous changes in these phenotypic characteristics and the presence of a fitness trade-off between resistance to infection and proliferative potential are explicitly incorporated into the model. The model comprises an ordinary differential equation for the number of infected individuals that is coupled with a partial integrodifferential equation for the population density function of susceptible individuals through an integral term.

A Mathematical Study of the Influence of Hypoxia and Acidity on the Evolutionary Dynamics of Cancer.

Hypoxia and acidity act as environmental stressors promoting selection for cancer cells with a more aggressive phenotype. As a result, a deeper theoretical understanding of the spatio-temporal processes that drive the adaptation of tumour cells to hypoxic and acidic microenvironments may open up new avenues of research in oncology and cancer treatment. We present a mathematical model to study the influence of hypoxia and acidity on the evolutionary dynamics of cancer cells in vascularised tumours.

Mechanical Models of Pattern and Form in Biological Tissues: The Role of Stress-Strain Constitutive Equations.

Mechanical and mechanochemical models of pattern formation in biological tissues have been used to study a variety of biomedical systems, particularly in developmental biology, and describe the physical interactions between cells and their local surroundings. These models in their original form consist of a balance equation for the cell density, a balance equation for the density of the extracellular matrix (ECM), and a force-balance equation describing the mechanical equilibrium of the cell-ECM system. Under the assumption that the cell-ECM system can be regarded as an isotropic linear viscoelastic material, the force-balance equation is often defined using the Kelvin-Voigt model of linear viscoelasticity to represent the stress-strain relation of the ECM.

A hybrid discrete-continuum approach to model Turing pattern formation.

Since its introduction in 1952, with a further refinement in 1972 by Gierer and Meinhardt, Turing's (pre-)pattern theory (the chemical basis of morphogenesis) has been widely applied to a number of areas in developmental biology, where evolving cell and tissue structures are naturally observed. The related pattern formation models normally comprise a system of reaction-diffusion equations for interacting chemical species (morphogens), whose heterogeneous distribution in some spatial domain acts as a template for cells to form some kind of pattern or structure through, for example, differentiation or proliferation induced by the chemical pre-pattern. Here we develop a hybrid discrete-continuum modelling framework for the formation of cellular patterns via the Turing mechanism.

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

Deterministic continuum models formulated as nonlocal 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 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 of two competing phenotype-structured cell populations that are exposed to time-varying nutrient levels and undergo spontaneous, heritable phenotypic changes with different probabilities.

A Mathematical Dissection of the Adaptation of Cell Populations to Fluctuating Oxygen Levels.

The disordered network of blood vessels that arises from tumour angiogenesis results in variations in the delivery of oxygen into the tumour tissue. This brings about regions of chronic hypoxia (i.e.

From a discrete model of chemotaxis with volume-filling to a generalized Patlak-Keller-Segel model.

We present a discrete model of chemotaxis whereby cells responding to a chemoattractant are seen as individual agents whose movement is described through a set of rules that result in a biased random walk. In order to take into account possible alterations in cellular motility observed at high cell densities (i.e.

Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments.

Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels.

A structured population model of clonal selection in acute leukemias with multiple maturation stages.

Recent progress in genetic techniques has shed light on the complex co-evolution of malignant cell clones in leukemias. However, several aspects of clonal selection still remain unclear. In this paper, we present a multi-compartmental continuously structured population model of selection dynamics in acute leukemias, which consists of a system of coupled integro-differential equations.

Bridging the gap between individual-based and continuum models of growing cell populations.

Continuum models for the spatial dynamics of growing cell populations have been widely used to investigate the mechanisms underpinning tissue development and tumour invasion. These models consist of nonlinear partial differential equations that describe the evolution of cellular densities in response to pressure gradients generated by population growth. Little prior work has explored the relation between such continuum models and related single-cell-based models.

A Mathematical Framework for Modelling the Metastatic Spread of Cancer.

Cancer is a complex disease that starts with mutations of key genes in one cell or a small group of cells at a primary site in the body. If these cancer cells continue to grow successfully and, at some later stage, invade the surrounding tissue and acquire a vascular network, they can spread to distant secondary sites in the body. This process, known as metastatic spread, is responsible for around 90% of deaths from cancer and is one of the so-called hallmarks of cancer.

The role of spatial variations of abiotic factors in mediating intratumour phenotypic heterogeneity.

We present here a space- and phenotype-structured model of selection dynamics between cancer cells within a solid tumour. In the framework of this model, we combine formal analyses with numerical simulations to investigate in silico the role played by the spatial distribution of abiotic components of the tumour microenvironment in mediating phenotypic selection of cancer cells. Numerical simulations are performed both on the 3D geometry of an in silico multicellular tumour spheroid and on the 3D geometry of an in vivo human hepatic tumour, which was imaged using computerised tomography.

Modelling the Immune Response to Cancer: An Individual-Based Approach Accounting for the Difference in Movement Between Inactive and Activated T Cells.

A growing body of experimental evidence indicates that immune cells move in an unrestricted search pattern if they are in the pre-activated state, whilst they tend to stay within a more restricted area upon activation induced by the presence of tumour antigens. This change in movement is not often considered in the existing mathematical models of the interactions between immune cells and cancer cells. With the aim to fill such a gap in the existing literature, in this work we present a spatially structured individual-based model of tumour-immune competition that takes explicitly into account the difference in movement between inactive and activated immune cells.

Examining the role of individual movement in promoting coexistence in a spatially explicit prisoner's dilemma.

The emergence of cooperation is a major conundrum of evolutionary biology. To unravel this evolutionary riddle, several models have been developed within the theoretical framework of spatial game theory, focussing on the interactions between two general classes of player, "cooperators" and "defectors". Generally, explicit movement in the spatial domain is not considered in these models, with strategies moving via imitation or through colonisation of neighbouring sites.

Emergence of spatial patterns in a mathematical model for the co-culture dynamics of epithelial-like and mesenchymal-like cells.

Accumulating evidence indicates that the interaction between epithelial and mesenchymal cells plays a pivotal role in cancer development and metastasis formation. Here we propose an integro-differential model for the co-culture dynamics of epithelial-like and mesenchymal-like cells. Our model takes into account the effects of chemotaxis, adhesive interactions between epithelial-like cells, proliferation and competition for nutrients.

Tracking the evolution of cancer cell populations through the mathematical lens of phenotype-structured equations.

Background: A thorough understanding of the ecological and evolutionary mechanisms that drive the phenotypic evolution of neoplastic cells is a timely and key challenge for the cancer research community. In this respect, mathematical modelling can complement experimental cancer research by offering alternative means of understanding the results of in vitro and in vivo experiments, and by allowing for a quick and easy exploration of a variety of biological scenarios through in silico studies.

Results: To elucidate the roles of phenotypic plasticity and selection pressures in tumour relapse, we present here a phenotype-structured model of evolutionary dynamics in a cancer cell population which is exposed to the action of a cytotoxic drug.