Turing pattern design principles and their robustness

Numerical investigation of new rete ridge formation in a multi-layer model of skin subjected to tissue expansion

The skin is a multilayered organ with microstructural and antomical heterogeneities that contribute to its unique mechanophysiology. Between the epidermis layer at the top and the dermis layer below, the basal keratinocytes form an interface with sinusoidal-like geometry termed rete ridges. In previous computational work we showed that the rete ridges contribute to lower delamination risk by increasing surface area and reducing the stress jump across the interface. Experimentally, we and others have shown that upon repeated tissue expansion and growth, physiological rete ridge frequency is preserved. Here we implement a 2D multilayered skin model where each layer is able to grow in response to applied loading toward recovering the layer-specific homeostatic stretch. Our simulations support the hypothesis that mechanics of growing tissue can explain secondary buckling and new rete ridge formation in tissue expansion. The process is robust with respect to parameters such as homeostatic stretch, layer thicknesses, and shear moduli of the different layers. Thicker epidermis suppresses higher frequency features, and so does a stiffer epidermis with respect to the basal layer. Interestingly, new rete ridge valleys are formed at locations that were originally peaks of the sine wave, whereas original valleys remain valleys. This pattern might have a connection to the localization of stem cell and transient amplifying cells in the epidermis. This study does not discard the role of cell-cell signaling dynamics, but rather emphasizes the possibility of achieving robust geometric patterns with simple rules of growing tissue, even in the absence of complex regulatory networks.

Automatic design of gene regulatory mechanisms for spatial pattern formation

Gene regulatory mechanisms (GRMs) control the formation of spatial and temporal expression patterns that can serve as regulatory signals for the development of complex shapes. Synthetic developmental biology aims to engineer such genetic circuits for understanding and producing desired multicellular spatial patterns. However, designing synthetic GRMs for complex, multi-dimensional spatial patterns is a current challenge due to the nonlinear interactions and feedback loops in genetic circuits. Here we present a methodology to automatically design GRMs that can produce any given two-dimensional spatial pattern. The proposed approach uses two orthogonal morphogen gradients acting as positional information signals in a multicellular tissue area or culture, which constitutes a continuous field of engineered cells implementing the same designed GRM. To efficiently design both the circuit network and the interaction mechanisms-including the number of genes necessary for the formation of the target spatial pattern-we developed an automated algorithm based on high-performance evolutionary computation. The tolerance of the algorithm can be configured to design GRMs that are either simple to produce approximate patterns or complex to produce precise patterns. We demonstrate the approach by automatically designing GRMs that can produce a diverse set of synthetic spatial expression patterns by interpreting just two orthogonal morphogen gradients. The proposed framework offers a versatile approach to systematically design and discover complex genetic circuits producing spatial patterns.

Robust cytoplasmic partitioning by solving an intrinsic cytoskeletal instability

Early development across vertebrates and insects critically relies on robustly reorganizing the cytoplasm of fertilized eggs into individualized cells. This intricate process is orchestrated by large microtubule structures that traverse the embryo, partitioning the cytoplasm into physically distinct and stable compartments. Despite the robustness of embryonic development, here we uncover an intrinsic instability in cytoplasmic partitioning driven by the microtubule cytoskeleton. We reveal that embryos circumvent this instability through two distinct mechanisms: either by matching the cell cycle duration to the time needed for the instability to unfold or by limiting microtubule nucleation. These regulatory mechanisms give rise to two possible strategies to fill the cytoplasm, which we experimentally demonstrate in zebrafish and Drosophila embryos, respectively. In zebrafish embryos, unstable microtubule waves fill the geometry of the entire embryo from the first division. Conversely, in Drosophila embryos, stable microtubule asters resulting from reduced microtubule nucleation gradually fill the cytoplasm throughout multiple divisions. Our results indicate that the temporal control of microtubule dynamics could have driven the evolutionary emergence of species-specific mechanisms for effective cytoplasmic organization. Furthermore, our study unveils a fundamental synergy between physical instabilities and biological clocks, uncovering universal strategies for rapid, robust, and efficient spatial ordering in biological systems.

Minimal reaction schemes for pattern formation

We link continuum models of reaction-diffusion systems that exhibit diffusion-driven instability to constraints on the particle-scale interactions underpinning this instability. While innumerable biological, chemical and physical patterns have been studied through the lens of Alan Turing's reaction-diffusion pattern-forming mechanism, the connections between models of pattern formation and the nature of the particle interactions generating them have been relatively understudied in comparison with the substantial efforts that have been focused on understanding proposed continuum systems. To derive the necessary reactant combinations for the most parsimonious reaction schemes, we analyse the emergent continuum models in terms of possible generating elementary reaction schemes. This analysis results in the complete list of such schemes containing the fewest reactions; these are the simplest possible hypothetical mass-action models for a pattern-forming system of two interacting species.

Turing Instabilities are Not Enough to Ensure Pattern Formation

Symmetry-breaking instabilities play an important role in understanding the mechanisms underlying the diversity of patterns observed in nature, such as in Turing's reaction-diffusion theory, which connects cellular signalling and transport with the development of growth and form. Extensive literature focuses on the linear stability analysis of homogeneous equilibria in these systems, culminating in a set of conditions for transport-driven instabilities that are commonly presumed to initiate self-organisation. We demonstrate that a selection of simple, canonical transport models with only mild multistable non-linearities can satisfy the Turing instability conditions while also robustly exhibiting only transient patterns. Hence, a Turing-like instability is insufficient for the existence of a patterned state. While it is known that linear theory can fail to predict the formation of patterns, we demonstrate that such failures can appear robustly in systems with multiple stable homogeneous equilibria. Given that biological systems such as gene regulatory networks and spatially distributed ecosystems often exhibit a high degree of multistability and nonlinearity, this raises important questions of how to analyse prospective mechanisms for self-organisation.

A multifunctional Wnt regulator underlies the evolution of rodent stripe patterns

Animal pigment patterns are excellent models to elucidate mechanisms of biological organization. Although theoretical simulations, such as Turing reaction-diffusion systems, recapitulate many animal patterns, they are insufficient to account for those showing a high degree of spatial organization and reproducibility. Here, we study the coat of the African striped mouse (Rhabdomys pumilio) to uncover how periodic stripes form. Combining transcriptomics, mathematical modelling and mouse transgenics, we show that the Wnt modulator Sfrp2 regulates the distribution of hair follicles and establishes an embryonic prepattern that foreshadows pigment stripes. Moreover, by developing in vivo gene editing in striped mice, we find that Sfrp2 knockout is sufficient to alter the stripe pattern. Strikingly, mutants exhibited changes in pigmentation, revealing that Sfrp2 also regulates hair colour. Lastly, through evolutionary analyses, we find that striped mice have evolved lineage-specific changes in regulatory elements surrounding Sfrp2, many of which may be implicated in modulating the expression of this gene. Altogether, our results show that a single factor controls coat pattern formation by acting both as an orienting signalling mechanism and a modulator of pigmentation. More broadly, our work provides insights into how spatial patterns are established in developing embryos and the mechanisms by which phenotypic novelty originates.

Diffusiophoresis-enhanced Turing patterns

Turing patterns are fundamental in biophysics, emerging from short-range activation and long-range inhibition processes. However, their paradigm is based on diffusive transport processes that yield patterns with shallower gradients than those observed in nature. A complete physical description of this discrepancy remains unknown. We propose a solution to this phenomenon by investigating the role of diffusiophoresis, which is the propulsion of colloids by a chemical gradient, in Turing patterns. Diffusiophoresis enables robust patterning of colloidal particles with substantially finer length scales than the accompanying chemical Turing patterns. A scaling analysis and a comparison to recent experiments indicate that chromatophores, ubiquitous in biological pattern formation, are likely diffusiophoretic and the colloidal Péclet number controls the pattern enhancement. This discovery suggests that important features of biological pattern formation can be explained with a universal mechanism that is quantified straightforwardly from the fundamental physics of colloids.

The diverse nature of intestinal fibroblasts in development, homeostasis, and disease

Only in recent years have we begun to appreciate the involvement of fibroblasts in intestinal development, tissue homeostasis, and disease. These insights followed the advent of single-cell transcriptomics that allowed researchers to explore the heterogeneity of intestinal fibroblasts in unprecedented detail. Since researchers often defined cell types and their associated function based on the biological process they studied, there are a plethora of partially overlapping markers for different intestinal fibroblast populations. This ambiguity complicates putting different research findings into context. Here, we provide a census on the function and identity of intestinal fibroblasts in mouse and human. We propose a simplified framework consisting of three colonic and four small intestinal fibroblast populations to aid navigating the diversity of intestinal fibroblasts.

The reaction-diffusion basis of animated patterns in eukaryotic flagella

The flagellar beat of bull spermatozoa and C. Reinhardtii are modelled by a minimal, geometrically exact, reaction-diffusion system. Spatio-temporal animated patterns describe flagellar waves, analogous to chemical-patterns from classical reaction-diffusion systems, with sliding-controlled molecular motor reaction-kinetics. The reaction-diffusion system is derived from first principles as a consequence of the high-internal dissipation by the flagellum relative to the external hydrodynamic dissipation. Quantitative comparison with nonlinear, large-amplitude simulations shows that animated reaction-diffusion patterns account for the experimental beating of both bull sperm and C. Reinhardtii. Our results suggest that a unified mechanism may exist for motors controlled by sliding, without requiring curvature-sensing, and uninfluenced by hydrodynamics. High-internal dissipation instigates autonomous travelling waves independently of the external fluid, enabling progressive swimming, otherwise not possible, in low viscosity environments, potentially critical for external fertilizers and aquatic microorganisms. The reaction-diffusion system may prove a powerful tool for studying pattern formation of movement on animated structures.

Physical interactions in non-ideal fluids promote Turing patterns

Turing's mechanism is often invoked to explain periodic patterns in nature, although direct experimental support is scarce. Turing patterns form in reaction-diffusion systems when the activating species diffuse much slower than the inhibiting species, and the involved reactions are highly nonlinear. Such reactions can originate from cooperativity, whose physical interactions should also affect diffusion. We here take direct interactions into account and show that they strongly affect Turing patterns. We find that weak repulsion between the activator and inhibitor can substantially lower the required differential diffusivity and reaction nonlinearity. By contrast, strong interactions can induce phase separation, but the resulting length scale is still typically governed by the fundamental reaction-diffusion length scale. Taken together, our theory connects traditional Turing patterns with chemically active phase separation, thus describing a wider range of systems. Moreover, we demonstrate that even weak interactions affect patterns substantially, so they should be incorporated when modelling realistic systems.

Learning system parameters from turing patterns

The Turing mechanism describes the emergence of spatial patterns due to spontaneous symmetry breaking in reaction-diffusion processes and underlies many developmental processes. Identifying Turing mechanisms in biological systems defines a challenging problem. This paper introduces an approach to the prediction of Turing parameter values from observed Turing patterns. The parameter values correspond to a parametrized system of reaction-diffusion equations that generate Turing patterns as steady state. The Gierer-Meinhardt model with four parameters is chosen as a case study. A novel invariant pattern representation based on resistance distance histograms is employed, along with Wasserstein kernels, in order to cope with the highly variable arrangement of local pattern structure that depends on the initial conditions which are assumed to be unknown. This enables us to compute physically plausible distances between patterns, to compute clusters of patterns and, above all, model parameter prediction based on training data that can be generated by numerical model evaluation with random initial data: for small training sets, classical state-of-the-art methods including operator-valued kernels outperform neural networks that are applied to raw pattern data, whereas for large training sets the latter are more accurate. A prominent property of our approach is that only a single pattern is required as input data for model parameter predicion. Excellent predictions are obtained for single parameter values and reasonably accurate results for jointly predicting all four parameter values.

Identifying patterning behavior in a plant infestation of insect pests

In this study, we developed a mechanistic model formulated as a system of reaction-diffusion equations (RDE) to explore the spatiotemporal dynamics of a theoretical pest with a tillering host plant in a controlled rectangular plant field. Local perturbation analysis, a recently developed method of analysis for wave propagation, was utilized to determine patterning regimes resulting from the local and global behaviors of the slow and fast diffusing components of the RDE system, respectively. Turing analysis was done to show that the RDE system does not exhibit Turing patterns. With bug mortality as the bifurcation parameter, regions with oscillations and stable coexistence of the pest and tillers were identified. Numerical simulations illustrate the patterning regimes in 1D and 2D settings. The oscillations suggest that recurrences in pest infestation is possible. Moreover, simulations showed that patterns produced in the model are strongly influenced by the pests' homogeneous dynamics inside the controlled environment.

Universal structures for adaptation in biochemical reaction networks

At the molecular level, the evolution of life is driven by the generation and diversification of adaptation mechanisms. A universal description of adaptation-capable chemical reaction network (CRN) structures has remained elusive until now, since currently-known criteria for adaptation apply only to a tiny subset of possible CRNs. Here we identify the definitive structural requirements that characterize all adaptation-capable collections of interacting molecules, however large or complex. We show that these network structures implement a form of integral control in which multiple independent integrals can collaborate to confer the capacity for adaptation on specific molecules. Using an algebraic algorithm informed by these findings, we demonstrate the existence of embedded integrals in a variety of biologically important CRNs that have eluded previous methods, and for which adaptation has been observed experimentally. This definitive picture of biological adaptation at the level of intermolecular interactions represents a blueprint for adaptation-capable signaling networks across all domains of life, and for the design of synthetic biosystems.

Universal dynamics of biological pattern formation in spatio-temporal morphogen variations

In biological systems, chemical signals termed morphogens self-organize into patterns that are vital for many physiological processes. As observed by Turing in 1952, these patterns are in a state of continual development, and are usually transitioning from one pattern into another. How do cells robustly decode these spatio-temporal patterns into signals in the presence of confounding effects caused by unpredictable or heterogeneous environments? Here, we answer this question by developing a general theory of pattern formation in spatio-temporal variations of ‘pre-pattern’ morphogens, which determine gene-regulatory network parameters. Through mathematical analysis, we identify universal dynamical regimes that apply to wide classes of biological systems. We apply our theory to two paradigmatic pattern-forming systems, and predict that they are robust with respect to non-physiological morphogen variations. More broadly, our theoretical framework provides a general approach to classify the emergent dynamics of pattern-forming systems based on how the bifurcations in their governing equations are traversed.

Appearance and suppression of Turing patterns under a periodically forced feed

Turing instability is a general and straightforward mechanism of pattern formation in reaction-diffusion systems, and its relevance has been demonstrated in different biological phenomena. Still, there are many open questions, especially on the robustness of the Turing mechanism. Robust patterns must survive some variation in the environmental conditions. Experiments on pattern formation using chemical systems have shown many reaction-diffusion patterns and serve as relatively simple test tools to study general aspects of these phenomena. Here, we present a study of sinusoidal variation of the input feed concentrations on chemical Turing patterns. Our experimental, numerical and theoretical analysis demonstrates that patterns may appear even at significant amplitude variation of the input feed concentrations. Furthermore, using time-dependent feeding opens a way to control pattern formation. The patterns settled at constant feed may disappear, or new patterns may appear from a homogeneous steady state due to the periodic forcing.

A computational modeling approach for predicting multicell spheroid patterns based on signaling-induced differential adhesion

Physiological and pathological processes including embryogenesis and tumorigenesis rely on the ability of individual cells to work collectively to form multicell patterns. In these heterogeneous multicell systems, cell-cell signaling induces differential adhesion between cells that leads to tissue-level patterning. However, the sensitivity of pattern formation to changes in the strengths of signaling or cell adhesion processes is not well understood. Prior work has explored these issues using synthetically engineered heterogeneous multicell spheroid systems, in which cell subpopulations engage in bidirectional intercellular signaling to regulate the expression of different cadherins. While engineered cell systems provide excellent experimental tools to observe pattern formation in cell populations, computational models of these systems may be leveraged to explore more systematically how specific combinations of signaling and adhesion parameters can drive the emergence of unique patterns. We developed and validated two- and three-dimensional agent-based models (ABMs) of spheroid patterning for previously described cells engineered with a bidirectional signaling circuit that regulates N- and P-cadherin expression. Systematic exploration of model predictions, some of which were experimentally validated, revealed how cell seeding parameters, the order of signaling events, probabilities of induced cadherin expression, and homotypic adhesion strengths affect pattern formation. Unsupervised clustering was also used to map combinations of signaling and adhesion parameters to these unique spheroid patterns predicted by the ABM. Finally, we demonstrated how the model may be deployed to design new synthetic cell signaling circuits based on a desired final multicell pattern.

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.

Introduction to 'Recent progress and open frontiers in Turing's theory of morphogenesis'

Elucidating pattern forming processes is an important problem in the physical, chemical and biological sciences. Turing's contribution, after being initially neglected, eventually catalysed a huge amount of work from mathematicians, physicists, chemists and biologists aimed towards understanding how steady spatial patterns can emerge from homogeneous chemical mixtures due to the reaction and diffusion of different chemical species. While this theory has been developed mathematically and investigated experimentally for over half a century, many questions still remain unresolved. This theme issue places Turing's theory of pattern formation in a modern context, discussing the current frontiers in foundational aspects of pattern formation in reaction-diffusion and related systems. It highlights ongoing work in chemical, synthetic and developmental settings which is helping to elucidate how important Turing's mechanism is for real morphogenesis, while highlighting gaps that remain in matching theory to reality. The theme issue also surveys a variety of recent mathematical research pushing the boundaries of Turing's original theory to more realistic and complicated settings, as well as discussing open theoretical challenges in the analysis of such models. It aims to consolidate current research frontiers and highlight some of the most promising future directions. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.