Influence of stochastic domain growth on pattern nucleation for diffusive systems with internal noise

A simple and fast method for estimating bat roost locations

Bats play a pivotal role in pest control, pollination and seed dispersal. Despite their ecological significance, locating bat roosts remains a challenging task for ecologists. Traditional field surveys are time-consuming, expensive and may disturb sensitive bat populations. In this article, we combine data from static audio detectors with a bat movement model to facilitate the detection of bat roosts. Crucially, our technique not only provides a point prediction for the most likely location of a bat roost, but because of the algorithm's speed, it can be applied over an entire landscape, resulting in a likelihood map, which provides optimal searching regions. To illustrate the success of the algorithm and highlight limitations, we apply our technique to greater horseshoe bat (Rhinolophus ferrumequinum) acoustic data acquired from six surveys from four different UK locations and over six different times in the year. Furthermore, we investigate what happens to the accuracy of our predictions in the case that the roost is not contained within the area spanned by the detectors. This innovative approach to searching rural environments holds the potential to greatly reduce the labour required for roost finding, and, hence, enhance the conservation efforts of bat populations and their habitats.

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.

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.

Accounting for dimensional differences in stochastic domain invasion with applications to precancerous cell removal

We consider a specific form of domain invasion that is an abstraction of pancreatic tissue eliminating precancerous mutant cells through juxtacrine signalling. The model is explored discretely, continuously, stochastically and deterministically, highlighting unforeseen nonlinear dependencies on the dimension of the solution domain. Specifically, stochastically simulated populations invade with a dimension dependent wave speed that can be over twice as fast as their deterministic analogues. Although the wave speed can be analytically derived in the cases of small domains, the probabilistic state space grows exponentially and, thus, we use numeric simulation and curve fitting to predict limiting dynamics.

Turing pattern design principles and their robustness

Turing patterns have morphed from mathematical curiosities into highly desirable targets for synthetic biology. For a long time, their biological significance was sometimes disputed but there is now ample evidence for their involvement in processes ranging from skin pigmentation to digit and limb formation. While their role in developmental biology is now firmly established, their synthetic design has so far proved challenging. Here, we review recent large-scale mathematical analyses that have attempted to narrow down potential design principles. We consider different aspects of robustness of these models and outline why this perspective will be helpful in the search for synthetic Turing-patterning systems. We conclude by considering robustness in the context of developmental modelling more generally. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

The design principles of discrete turing patterning systems

The formation of spatial structures lies at the heart of developmental processes. However, many of the underlying gene regulatory and biochemical processes remain poorly understood. Turing patterns constitute a main candidate to explain such processes, but they appear sensitive to fluctuations and variations in kinetic parameters, raising the question of how they may be adopted and realised in naturally evolved systems. The vast majority of mathematical studies of Turing patterns have used continuous models specified in terms of partial differential equations. Here, we complement this work by studying Turing patterns using discrete cellular automata models. We perform a large-scale study on all possible two-species networks and find the same Turing pattern producing networks as in the continuous framework. In contrast to continuous models, however, we find these Turing pattern topologies to be substantially more robust to changes in the parameters of the model. We also find that diffusion-driven instabilities are substantially weaker predictors for Turing patterns in our discrete modelling framework in comparison to the continuous case, in the sense that the presence of an instability does not guarantee a pattern emerging in simulations. We show that a more refined criterion constitutes a stronger predictor. The similarity of the results for the two modelling frameworks suggests a deeper underlying principle of Turing mechanisms in nature. Together with the larger robustness in the discrete case this suggests that Turing patterns may be more robust than previously thought.

EPHA2-dependent outcompetition of KRASG12D mutant cells by wild-type neighbors in the adult pancreas

As we age, our tissues are repeatedly challenged by mutational insult, yet cancer occurrence is a relatively rare event. Cells carrying cancer-causing genetic mutations compete with normal neighbors for space and survival in tissues. However, the mechanisms underlying mutant-normal competition in adult tissues and the relevance of this process to cancer remain incompletely understood. Here, we investigate how the adult pancreas maintains tissue health in vivo following sporadic expression of oncogenic Kras (KrasG12D), the key driver mutation in human pancreatic cancer. We find that when present in tissues in low numbers, KrasG12D mutant cells are outcompeted and cleared from exocrine and endocrine compartments in vivo. Using quantitative 3D tissue imaging, we show that before being cleared, KrasG12D cells lose cell volume, pack into round clusters, and E-cadherin-based cell-cell adhesions decrease at boundaries with normal neighbors. We identify EphA2 receptor as an essential signal in the clearance of KrasG12D cells from exocrine and endocrine tissues in vivo. In the absence of functional EphA2, KrasG12D cells do not alter cell volume or shape, E-cadherin-based cell-cell adhesions increase and KrasG12D cells are retained in tissues. The retention of KRasG12D cells leads to the early appearance of premalignant pancreatic intraepithelial neoplasia (PanINs) in tissues. Our data show that adult pancreas tissues remodel to clear KrasG12D cells and maintain tissue health. This study provides evidence to support a conserved functional role of EphA2 in Ras-driven cell competition in epithelial tissues and suggests that EphA2 is a novel tumor suppressor in pancreatic cancer.

Coloured Noise from Stochastic Inflows in Reaction-Diffusion Systems

In this paper, we present a framework for investigating coloured noise in reaction-diffusion systems. We start by considering a deterministic reaction-diffusion equation and show how external forcing can cause temporally correlated or coloured noise. Here, the main source of external noise is considered to be fluctuations in the parameter values representing the inflow of particles to the system. First, we determine which reaction systems, driven by extrinsic noise, can admit only one steady state, so that effects, such as stochastic switching, are precluded from our analysis. To analyse the steady-state behaviour of reaction systems, even if the parameter values are changing, necessitates a parameter-free approach, which has been central to algebraic analysis in chemical reaction network theory. To identify suitable models, we use tools from real algebraic geometry that link the network structure to its dynamical properties. We then make a connection to internal noise models and show how power spectral methods can be used to predict stochastically driven patterns in systems with coloured noise. In simple cases, we show that the power spectrum of the coloured noise process and the power spectrum of the reaction-diffusion system modelled with white noise multiply to give the power spectrum of the coloured noise reaction-diffusion system.

Pattern production through a chiral chasing mechanism

Recent experiments on zebrafish pigmentation suggests that their typical black and white striped skin pattern is made up of a number of interacting chromatophore families. Specifically, two of these cell families have been shown to interact through a nonlocal chasing mechanism, which has previously been modeled using integro-differential equations. We extend this framework to include the experimentally observed fact that the cells often exhibit chiral movement, in that the cells chase, and run away, at angles different to the line connecting their centers. This framework is simplified through the use of multiple small limits leading to a coupled set of partial differential equations which are amenable to Fourier analysis. This analysis results in the production of dispersion relations and necessary conditions for a patterning instability to occur. Beyond the theoretical development and the production of new pattern planiforms we are able to corroborate the experimental hypothesis that the global pigmentation patterns can be dependent on the chirality of the chromatophores.

Noise-induced temporal dynamics in Turing systems

We examine the ability of intrinsic noise to produce complex temporal dynamics in Turing pattern formation systems, with particular emphasis on the Schnakenberg kinetics. Using power spectral methods, we characterize the behavior of the system using stochastic simulations at a wide range of points in parameter space and compare with analytical approximations. Specifically, we investigate whether polarity switching of stochastic patterns occurs at a defined frequency. We find that it can do so in individual realizations of a stochastic simulation, but that the frequency is not defined consistently across realizations in our samples of parameter space. Further, we examine the effect of noise on deterministically predicted traveling waves and find them increased in amplitude and decreased in speed.

Turing's model for biological pattern formation and the robustness problem

One of the fundamental questions in developmental biology is how the vast range of pattern and structure we observe in nature emerges from an almost uniformly homogeneous fertilized egg. In particular, the mechanisms by which biological systems maintain robustness, despite being subject to numerous sources of noise, are shrouded in mystery. Postulating plausible theoretical models of biological heterogeneity is not only difficult, but it is also further complicated by the problem of generating robustness, i.e. once we can generate a pattern, how do we ensure that this pattern is consistently reproducible in the face of perturbations to the domain, reaction time scale, boundary conditions and so forth. In this paper, not only do we review the basic properties of Turing's theory, we highlight the successes and pitfalls of using it as a model for biological systems, and discuss emerging developments in the area.

Effects of intrinsic stochasticity on delayed reaction-diffusion patterning systems

Cellular gene expression is a complex process involving many steps, including the transcription of DNA and translation of mRNA; hence the synthesis of proteins requires a considerable amount of time, from ten minutes to several hours. Since diffusion-driven instability has been observed to be sensitive to perturbations in kinetic delays, the application of Turing patterning mechanisms to the problem of producing spatially heterogeneous differential gene expression has been questioned. In deterministic systems a small delay in the reactions can cause a large increase in the time it takes a system to pattern. Recently, it has been observed that in undelayed systems intrinsic stochasticity can cause pattern initiation to occur earlier than in the analogous deterministic simulations. Here we are interested in adding both stochasticity and delays to Turing systems in order to assess whether stochasticity can reduce the patterning time scale in delayed Turing systems. As analytical insights to this problem are difficult to attain and often limited in their use, we focus on stochastically simulating delayed systems. We consider four different Turing systems and two different forms of delay. Our results are mixed and lead to the conclusion that, although the sensitivity to delays in the Turing mechanism is not completely removed by the addition of intrinsic noise, the effects of the delays are clearly ameliorated in certain specific cases.