Interactions between pattern formation and domain growth

Bat Motion can be Described by Leap Frogging

We present models of bat motion derived from radio-tracking data collected over 14 nights. The data presents an initial dispersal period and a return to roost period. Although a simple diffusion model fits the initial dispersal motion we show that simple convection cannot provide a description of the bats returning to their roost. By extending our model to include non-autonomous parameters, or a leap frogging form of motion, where bats on the exterior move back first, we find we are able to accurately capture the bat's motion. We discuss ways of distinguishing between the two movement descriptions and, finally, consider how the different motion descriptions would impact a bat's hunting strategy.

A mathematical model for nutrient-limited uniaxial growth of a compressible tissue

We consider the uniaxial growth of a tissue or colony of cells, where a nutrient (or some other chemical) required for cell proliferation is supplied at one end, and is consumed by the cells. An example would be the growth of a cylindrical yeast colony in the experiments described by Vulin et al. (2014). We develop a reaction-diffusion model of this scenario which couples nutrient concentration and cell density on a growing domain. A novel element of our model is that the tissue is assumed to be compressible. We define replicative regions, where cells have sufficient nutrient to proliferate, and quiescent regions, where the nutrient level is insufficient for this to occur. We also define pathlines, which allow us to track individual cell paths within the tissue. We begin our investigation of the model by considering an incompressible tissue where cell density is constant before exploring the solution space of the full compressible model. In a large part of the parameter space, the incompressible and compressible models give qualitatively similar results for both the nutrient concentration and cell pathlines, with the key distinction being the variation in density in the compressible case. In particular, the replicative region is located at the base of the tissue, where nutrient is supplied, and nutrient concentration decreases monotonically with distance from the nutrient source. However, for a highly-compressible tissue with small nutrient consumption rate, we observe a counter-intuitive scenario where the nutrient concentration is not necessarily monotonically decreasing, and there can be two replicative regions. For parameter values given in the paper by Vulin et al. (2014), the incompressible model slightly overestimates the colony length compared to experimental observations; this suggests the colony may be somewhat compressible. Both incompressible and compressible models predict that, for these parameter values, cell proliferation is ultimately confined to a small region close to the colony base.

Concentration-Dependent Domain Evolution in Reaction-Diffusion Systems

Pattern formation has been extensively studied in the context of evolving (time-dependent) domains in recent years, with domain growth implicated in ameliorating problems of pattern robustness and selection, in addition to more realistic modelling in developmental biology. Most work to date has considered prescribed domains evolving as given functions of time, but not the scenario of concentration-dependent dynamics, which is also highly relevant in a developmental setting. Here, we study such concentration-dependent domain evolution for reaction-diffusion systems to elucidate fundamental aspects of these more complex models. We pose a general form of one-dimensional domain evolution and extend this to N-dimensional manifolds under mild constitutive assumptions in lieu of developing a full tissue-mechanical model. In the 1D case, we are able to extend linear stability analysis around homogeneous equilibria, though this is of limited utility in understanding complex pattern dynamics in fast growth regimes. We numerically demonstrate a variety of dynamical behaviours in 1D and 2D planar geometries, giving rise to several new phenomena, especially near regimes of critical bifurcation boundaries such as peak-splitting instabilities. For sufficiently fast growth and contraction, concentration-dependence can have an enormous impact on the nonlinear dynamics of the system both qualitatively and quantitatively. We highlight crucial differences between 1D evolution and higher-dimensional models, explaining obstructions for linear analysis and underscoring the importance of careful constitutive choices in defining domain evolution in higher dimensions. We raise important questions in the modelling and analysis of biological systems, in addition to numerous mathematical questions that appear tractable in the one-dimensional setting, but are vastly more difficult for higher-dimensional models.

Boundary Conditions Cause Different Generic Bifurcation Structures in Turing Systems

Turing’s theory of morphogenesis is a generic mechanism to produce spatial patterning from near homogeneity. Although widely studied, we are still able to generate new results by returning to common dogmas. One such widely reported belief is that the Turing bifurcation occurs through a pitchfork bifurcation, which is true under zero-flux boundary conditions. However, under fixed boundary conditions, the Turing bifurcation becomes generically transcritical. We derive these algebraic results through weakly nonlinear analysis and apply them to the Schnakenberg kinetics. We observe that the combination of kinetics and boundary conditions produce their own uncommon boundary complexities that we explore numerically. Overall, this work demonstrates that it is not enough to only consider parameter perturbations in a sensitivity analysis of a specific application. Variations in boundary conditions should also be considered.

Effect of obstructions on growing Turing patterns

We study how Turing pattern formation on a growing domain is affected by discrete domain discontinuities. We use the Lengyel-Epstein reaction-diffusion model to numerically simulate Turing pattern formation on radially expanding circular domains containing a variety of obstruction geometries, including obstructions spanning the length of the domain, such as walls and slits, and local obstructions, such as small blocks. The pattern formation is significantly affected by the obstructions, leading to novel pattern morphologies. We show that obstructions can induce growth mode switching and disrupt local pattern formation and that these effects depend on the shape and placement of the objects as well as the domain growth rate. This work provides a customizable framework to perform numerical simulations on different types of obstructions and other heterogeneous domains, which may guide future numerical and experimental studies. These results may also provide new insights into biological pattern growth and formation, especially in non-idealized domains containing noise or discontinuities.

Insights from chemical systems into Turing-type morphogenesis

In 1952, Alan Turing proposed a theory showing how morphogenesis could occur from a simple two morphogen reaction-diffusion system [Turing, A. M. (1952) Phil. Trans. R. Soc. Lond. A 237, 37-72. (doi:10.1098/rstb.1952.0012)]. While the model is simple, it has found diverse applications in fields such as biology, ecology, behavioural science, mathematics and chemistry. Chemistry in particular has made significant contributions to the study of Turing-type morphogenesis, providing multiple reproducible experimental methods to both predict and study new behaviours and dynamics generated in reaction-diffusion systems. In this review, we highlight the historical role chemistry has played in the study of the Turing mechanism, summarize the numerous insights chemical systems have yielded into both the dynamics and the morphological behaviour of Turing patterns, and suggest future directions for chemical studies into Turing-type morphogenesis. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Influence of the Turing instability on the motion of domain boundaries

Turing's theory of pattern formation has provided crucial insights into the behavior of various biological, geographical, and chemical systems over the last few decades. Existing studies have focused on moving-boundary Turing systems for which the motion of the boundary is prescribed by an external agent. In this paper, we present an extension of this theory to a class of systems in which the front motion is governed by the physical processes that occur within the domain. Biological systems exhibiting apically dominant growth and corrosion of metals and alloys highlight some of the noteworthy examples of such systems. In this study, we characterize the nature of interaction between the moving front and the Turing-instability for both an activator-inhibitor and an activator-substrate model. Behavioral regimes of periodic, as well as nonperiodic (nonconstant), growth rates are obtained. Furthermore, the trends in the first show striking similarities with the cyclic-boundary-kinetics observed in experimental systems. In general, a stationary, periodic structure is also left behind the moving front. If the periodicity of the boundary kinetics agrees with the allowed range of the stable-periodic solutions, the pattern formed tends to persist. Otherwise, it evolves to a nearby energy-minimum either by peak-splitting, peak-decay, or by settling down to a spatially homogeneous state.

Dynamics of a Diffusive Two-Prey-One-Predator Model with Nonlocal Intra-Specific Competition for Both the Prey Species

Investigation of interacting populations is an active area of research, and various modeling approaches have been adopted to describe their dynamics. Mathematical models of such interactions using differential equations are capable to mimic the stationary and oscillating (regular or irregular) population distributions. Recently, some researchers have paid their attention to explain the consequences of transient dynamics of population density (especially the long transients) and able to capture such behaviors with simple models. Existence of multiple stationary patches and settlement to a stable distribution after a long quasi-stable transient dynamics can be explained by spatiotemporal models with nonlocal interaction terms. However, the studies of such interesting phenomena for three interacting species are not abundant in literature. Motivated by these facts here we have considered a three species prey–predator model where the predator is generalist in nature as it survives on two prey species. Nonlocalities are introduced in the intra-specific competition terms for the two prey species in order to model the accessibility of nearby resources. Using linear analysis, we have derived the Turing instability conditions for both the spatiotemporal models with and without nonlocal interactions. Validation of such conditions indicates the possibility of existence of stationary spatially heterogeneous distributions for all the three species. Existence of long transient dynamics has been presented under certain parametric domain. Exhaustive numerical simulations reveal various scenarios of stabilization of population distribution due to the presence of nonlocal intra-specific competition for the two prey species. Chaotic oscillation exhibited by the temporal model is significantly suppressed when the populations are allowed to move over their habitat and prey species can access the nearby resources.

Pattern formation in multiphase models of chemotactic cell aggregation

We develop a continuum model for the aggregation of cells cultured in a nutrient-rich medium in a culture well. We consider a 2D geometry, representing a vertical slice through the culture well, and assume that the cell layer depth is small compared with the typical lengthscale of the culture well. We adopt a continuum mechanics approach, treating the cells and culture medium as a two-phase mixture. Specifically, the cells and culture medium are treated as fluids. Additionally, the cell phase can generate forces in response to environmental cues, which include the concentration of a chemoattractant that is produced by the cells within the culture medium. The model leads to a system of coupled nonlinear partial differential equations for the volume fraction and velocity of the cell phase, the culture medium pressure and the chemoattractant concentration, which must be solved subject to appropriate boundary and initial conditions. To gain insight into the system, we consider two model reductions, appropriate when the cell layer depth is thin compared to the typical length scale of the culture well: a (simple) 1D and a (more involved) thin-film extensional flow reduction. By investigating the resulting systems of equations analytically and numerically, we identify conditions under which small amplitude perturbations to a homogeneous steady state (corresponding to a spatially uniform cell distribution) can lead to a spatially varying steady state (pattern formation). Our analysis reveals that the simpler 1D reduction has the same qualitative features as the thin-film extensional flow reduction in the linear and weakly nonlinear regimes, motivating the use of the simpler 1D modelling approach when a qualitative understanding of the system is required. However, the thin-film extensional flow reduction may be more appropriate when detailed quantitative agreement between modelling predictions and experimental data is desired. Furthermore, full numerical simulations of the two model reductions in regions of parameter space when the system is not close to marginal stability reveal significant differences in the evolution of the volume fraction and velocity of the cell phase, and chemoattractant concentration.

Reaction-diffusion patterns in plant tip morphogenesis: bifurcations on spherical caps

We study a chemical reaction-diffusion model (the Brusselator) for pattern formation on developing plant tips. A family of spherical cap domains is used to represent tip flattening during development. Applied to conifer embryos, we model the chemical prepatterning underlying cotyledon ("seed leaf") formation, and demonstrate the dependence of patterns on tip flatness, radius, and precursor concentrations. Parameters for the Brusselator in spherical cap domains can be chosen to give supercritical pitchfork bifurcations of patterned solutions of the nonlinear reaction-diffusion system that correspond to the cotyledon patterns that appear on the flattening tips of conifer embryos.

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.

Diffusion-reaction model for Drosophila embryo development

During the early stages of gastrulation in Drosophila embryo, the epithelial cells composing the single tissue layer of the egg undergo large strains and displacements. These movements have been usually modelled by decomposing the total deformation gradient in an (imposed or strain/stress dependent) active part and a passive response. Although the influence of the chemical and genetic activity in the mechanical response of the cell has been experimentally observed, the effects of the mechanical deformation on the latter have been far less studied, and much less modelled. Here, we propose a model that couples morphogen transport and the cell mechanics during embryogenesis. A diffusion-reaction equation is introduced as an additional mechanical regulator of morphogenesis. Consequently, the active deformations are not directly imposed in the analytical formulation, but they rather depend on the morphogen concentration, which is introduced as a new variable. In this study, we show that strain patterns similar to those observed during biological experiments can be reproduced by properly combining the two phenomena. In addition, we use a novel technique to parameterise the embryo geometry by solving two Laplace problems with specific boundary conditions. We apply the method to two morphogenetic movements: ventral furrow invagination and germ band extension. The matching between our results and the observed experimental deformations confirms that diffusion-reaction of morphogens can actually be controlling large morphogenetic movements.

Stochastic reaction and diffusion on growing domains: understanding the breakdown of robust pattern formation

Many biological patterns, from population densities to animal coat markings, can be thought of as heterogeneous spatiotemporal distributions of mobile agents. Many mathematical models have been proposed to account for the emergence of this complexity, but, in general, they have consisted of deterministic systems of differential equations, which do not take into account the stochastic nature of population interactions. One particular, pertinent criticism of these deterministic systems is that the exhibited patterns can often be highly sensitive to changes in initial conditions, domain geometry, parameter values, etc. Due to this sensitivity, we seek to understand the effects of stochasticity and growth on paradigm biological patterning models. In this paper, we extend spatial Fourier analysis and growing domain mapping techniques to encompass stochastic Turing systems. Through this we find that the stochastic systems are able to realize much richer dynamics than their deterministic counterparts, in that patterns are able to exist outside the standard Turing parameter range. Further, it is seen that the inherent stochasticity in the reactions appears to be more important than the noise generated by growth, when considering which wave modes are excited. Finally, although growth is able to generate robust pattern sequences in the deterministic case, we see that stochastic effects destroy this mechanism for conferring robustness. However, through Fourier analysis we are able to suggest a reason behind this lack of robustness and identify possible mechanisms by which to reclaim it.

The dynamics of Turing patterns for morphogen-regulated growing domains with cellular response delays

Since its conception in 1952, the Turing paradigm for pattern formation has been the subject of numerous theoretical investigations. Experimentally, this mechanism was first demonstrated in chemical reactions over 20 years ago and, more recently, several examples of biological self-organisation have also been implicated as Turing systems. One criticism of the Turing model is its lack of robustness, not only with respect to fluctuations in the initial conditions, but also with respect to the inclusion of delays in critical feedback processes such as gene expression. In this work we investigate the possibilities for Turing patterns on growing domains where the morphogens additionally regulate domain growth, incorporating delays in the feedback between signalling and domain growth, as well as gene expression. We present results for the proto-typical Schnakenberg and Gierer-Meinhardt systems: exploring the dynamics of these systems suggests a reconsideration of the basic Turing mechanism for pattern formation on morphogen-regulated growing domains as well as highlighting when feedback delays on domain growth are important for pattern formation.

The role of chemical dynamics in plant morphogenesis(1)

In biological development, the generation of shape is preceded by the spatial localization of growth factors. Localization, and how it is maintained or changed during the process of growth, determines the shapes produced. Mathematical models have been developed to investigate the chemical, mechanical and transport properties involved in plant morphogenesis. These synthesize biochemical and biophysical data, revealing underlying principles, especially the importance of dynamics in generating form. Chemical kinetics has been used to understand the constraints on reaction and transport rates to produce localized concentration patterns. This approach is well developed for understanding de novo pattern formation, pattern spacing and transitions from one pattern to another. For plants, growth is continual, and a key use of the theory is in understanding the feedback between patterning and growth, especially for morphogenetic events which break symmetry, such as tip branching. Within the context of morphogenetic modelling in general, the present review gives a brief history of chemical patterning research and its particular application to shape generation in plant development.

Aberrant behaviours of reaction diffusion self-organisation models on growing domains in the presence of gene expression time delays

Turing's pattern formation mechanism exhibits sensitivity to the details of the initial conditions suggesting that, in isolation, it cannot robustly generate pattern within noisy biological environments. Nonetheless, secondary aspects of developmental self-organisation, such as a growing domain, have been shown to ameliorate this aberrant model behaviour. Furthermore, while in-situ hybridisation reveals the presence of gene expression in developmental processes, the influence of such dynamics on Turing's model has received limited attention. Here, we novelly focus on the Gierer-Meinhardt reaction diffusion system considering delays due the time taken for gene expression, while incorporating a number of different domain growth profiles to further explore the influence and interplay of domain growth and gene expression on Turing's mechanism. We find extensive pathological model behaviour, exhibiting one or more of the following: temporal oscillations with no spatial structure, a failure of the Turing instability and an extreme sensitivity to the initial conditions, the growth profile and the duration of gene expression. This deviant behaviour is even more severe than observed in previous studies of Schnakenberg kinetics on exponentially growing domains in the presence of gene expression (Gaffney and Monk in Bull. Math. Biol. 68:99-130, 2006). Our results emphasise that gene expression dynamics induce unrealistic behaviour in Turing's model for multiple choices of kinetics and thus such aberrant modelling predictions are likely to be generic. They also highlight that domain growth can no longer ameliorate the excessive sensitivity of Turing's mechanism in the presence of gene expression time delays. The above, extensive, pathologies suggest that, in the presence of gene expression, Turing's mechanism would generally require a novel and extensive secondary mechanism to control reaction diffusion patterning.

Nonlinear modelling of cancer: bridging the gap between cells and tumours

Despite major scientific, medical and technological advances over the last few decades, a cure for cancer remains elusive. The disease initiation is complex, and including initiation and avascular growth, onset of hypoxia and acidosis due to accumulation of cells beyond normal physiological conditions, inducement of angiogenesis from the surrounding vasculature, tumour vascularization and further growth, and invasion of surrounding tissue and metastasis. Although the focus historically has been to study these events through experimental and clinical observations, mathematical modelling and simulation that enable analysis at multiple time and spatial scales have also complemented these efforts. Here, we provide an overview of this multiscale modelling focusing on the growth phase of tumours and bypassing the initial stage of tumourigenesis. While we briefly review discrete modelling, our focus is on the continuum approach. We limit the scope further by considering models of tumour progression that do not distinguish tumour cells by their age. We also do not consider immune system interactions nor do we describe models of therapy. We do discuss hybrid-modelling frameworks, where the tumour tissue is modelled using both discrete (cell-scale) and continuum (tumour-scale) elements, thus connecting the micrometre to the centimetre tumour scale. We review recent examples that incorporate experimental data into model parameters. We show that recent mathematical modelling predicts that transport limitations of cell nutrients, oxygen and growth factors may result in cell death that leads to morphological instability, providing a mechanism for invasion via tumour fingering and fragmentation. These conditions induce selection pressure for cell survivability, and may lead to additional genetic mutations. Mathematical modelling further shows that parameters that control the tumour mass shape also control its ability to invade. Thus, tumour morphology may serve as a predictor of invasiveness and treatment prognosis.