1
|
Krause AL, Gaffney EA, Jewell TJ, Klika V, Walker BJ. Turing Instabilities are Not Enough to Ensure Pattern Formation. Bull Math Biol 2024; 86:21. [PMID: 38253936 PMCID: PMC10803432 DOI: 10.1007/s11538-023-01250-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Accepted: 12/22/2023] [Indexed: 01/24/2024]
Abstract
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.
Collapse
Affiliation(s)
- Andrew L Krause
- Department of Mathematical Sciences, Durham University, Upper Mountjoy Campus, Stockton Road, Durham, DH1 3LE, UK.
| | - Eamonn A Gaffney
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Oxford, OX2 6GG, UK
| | - Thomas Jun Jewell
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Oxford, OX2 6GG, UK
| | - Václav Klika
- Department of Mathematics, FNSPE, Czech Technical University in Prague, Trojanova 13, 120 00, Prague, Czech Republic
| | - Benjamin J Walker
- Department of Mathematical Sciences, University of Bath, Bath, BA2 7AY, UK
- Department of Mathematics, University College London, London, WC1E 6BT, UK
| |
Collapse
|
2
|
Diez A, Krause AL, Maini PK, Gaffney EA, Seirin-Lee S. Turing Pattern Formation in Reaction-Cross-Diffusion Systems with a Bilayer Geometry. Bull Math Biol 2024; 86:13. [PMID: 38170298 PMCID: PMC10764571 DOI: 10.1007/s11538-023-01237-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Accepted: 11/06/2023] [Indexed: 01/05/2024]
Abstract
Conditions for self-organisation via Turing's mechanism in biological systems represented by reaction-diffusion or reaction-cross-diffusion models have been extensively studied. Nonetheless, the impact of tissue stratification in such systems is under-explored, despite its ubiquity in the context of a thin epithelium overlying connective tissue, for instance the epidermis and underlying dermal mesenchyme of embryonic skin. In particular, each layer can be subject to extensively different biochemical reactions and transport processes, with chemotaxis - a special case of cross-diffusion - often present in the mesenchyme, contrasting the solely molecular transport typically found in the epidermal layer. We study Turing patterning conditions for a class of reaction-cross-diffusion systems in bilayered regions, with a thin upper layer and coupled by a linear transport law. In particular, the role of differential transport through the interface is explored together with the presence of asymmetry between the homogeneous equilibria of the two layers. A linear stability analysis is carried out around a spatially homogeneous equilibrium state in the asymptotic limit of weak and strong coupling strengths, where quantitative approximations of the bifurcation curve can be computed. Our theoretical findings, for an arbitrary number of reacting species, reveal quantitative Turing conditions, highlighting when the coupling mechanism between the layered regions can either trigger patterning or stabilize a spatially homogeneous equilibrium regardless of the independent patterning state of each layer. We support our theoretical results through direct numerical simulations, and provide an open source code to explore such systems further.
Collapse
Affiliation(s)
- Antoine Diez
- Institute for the Advanced Study of Human Biology (ASHBi), Kyoto University Institute for Advanced Study (KUIAS), Kyoto University, Yoshida-Konoe-cho, Sakyo-ku, Kyoto, 606-8501, Japan
| | - Andrew L Krause
- Mathematical Sciences Department, Durham University, Upper Mountjoy Campus, Stockton Rd, Durham, DH1 3LE, UK
| | - Philip K Maini
- Wolfson Centre For Mathematical Biology, Mathematical Institute, Andrew Wiles Building, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG, UK
| | - Eamonn A Gaffney
- Wolfson Centre For Mathematical Biology, Mathematical Institute, Andrew Wiles Building, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG, UK.
| | - Sungrim Seirin-Lee
- Institute for the Advanced Study of Human Biology (ASHBi), Kyoto University Institute for Advanced Study (KUIAS), Kyoto University, Yoshida-Konoe-cho, Sakyo-ku, Kyoto, 606-8501, Japan.
- Department of Mathematical Medicine, Graduate School of Medicine, Kyoto University, Yoshida-Konoe-cho, Sakyo-ku, Kyoto, 606-8501, Japan.
| |
Collapse
|
3
|
Cisneros-Ake LA, Gonzalez-Rodriguez JC, González-Ramírez LR. Turing Instabilities and Rotating Spiral Waves in Glycolytic Processes. Bull Math Biol 2022; 84:100. [PMID: 35951127 DOI: 10.1007/s11538-022-01060-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 07/25/2022] [Indexed: 11/29/2022]
Abstract
We study single-frequency oscillations and pattern formation in the glycolytic process modeled by a reduction in the well-known Sel'kov's equations (Sel'kov in Eur J Biochem 4:79, 1968), which describe, in the whole cell, the phosphofructokinase enzyme reaction. By using averaging theory, we establish the existence conditions for limit cycles and their limiting average radius in the kinetic reaction equations. We analytically establish conditions on the model parameters for the appearance of unstable nonlinear modes seeding the formation of two-dimensional patterns in the form of classical spots and stripes. We also establish the existence of a Hopf bifurcation, which characterizes the reaction dynamics, producing glycolytic rotating spiral waves. We numerically establish parameter regions for the existence of these spiral waves and address their linear stability. We show that as the model tends toward a suppression of the relative source rate, the spiral wave solution loses stability. All our findings are validated by full numerical simulations of the model equations. Finally, we discuss in vitro evidence of spatiotemporal activity patterns found in glycolytic experiments, and propose plausible biological implications of our model results.
Collapse
Affiliation(s)
- Luis A Cisneros-Ake
- Departamento de Matemáticas, ESFM, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos Edificio 9, 07738, Mexico, Mexico.
| | - Juan C Gonzalez-Rodriguez
- Posgrado en Ciencias Fisicomatemáticas, ESFM, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos Edificio 9, 07738, Mexico, Mexico
| | - Laura R González-Ramírez
- Departamento de Matemáticas, ESFM, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos Edificio 9, 07738, Mexico, Mexico
| |
Collapse
|
4
|
Krause AL, Gaffney EA, Maini PK, Klika V. Introduction to 'Recent progress and open frontiers in Turing's theory of morphogenesis'. Philos Trans A Math Phys Eng Sci 2021; 379:20200280. [PMID: 34743606 PMCID: PMC8580473 DOI: 10.1098/rsta.2020.0280] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
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'.
Collapse
Affiliation(s)
- Andrew L. Krause
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
- Department of Mathematical Sciences, Durham University, Upper Mountjoy Campus, Stockton Rd, Durham DH1 3LE, UK
| | - Eamonn A Gaffney
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Philip K. Maini
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Václav Klika
- Department of Mathematics, FNSPE, Czech Technical University in Prague, Trojanova, 13, 120 00 Praha, Czech Republic
| |
Collapse
|
5
|
Abstract
Reaction-diffusion processes across layered media arise in several scientific domains such as pattern-forming E. coli on agar substrates, epidermal-mesenchymal coupling in development, and symmetry-breaking in cell polarization. We develop a modeling framework for bilayer reaction-diffusion systems and relate it to a range of existing models. We derive conditions for diffusion-driven instability of a spatially homogeneous equilibrium analogous to the classical conditions for a Turing instability in the simplest nontrivial setting where one domain has a standard reaction-diffusion system, and the other permits only diffusion. Due to the transverse coupling between these two regions, standard techniques for computing eigenfunctions of the Laplacian cannot be applied, and so we propose an alternative method to compute the dispersion relation directly. We compare instability conditions with full numerical simulations to demonstrate impacts of the geometry and coupling parameters on patterning, and explore various experimentally relevant asymptotic regimes. In the regime where the first domain is suitably thin, we recover a simple modulation of the standard Turing conditions, and find that often the broad impact of the diffusion-only domain is to reduce the ability of the system to form patterns. We also demonstrate complex impacts of this coupling on pattern formation. For instance, we exhibit non-monotonicity of pattern-forming instabilities with respect to geometric and coupling parameters, and highlight an instability from a nontrivial interaction between kinetics in one domain and diffusion in the other. These results are valuable for informing design choices in applications such as synthetic engineering of Turing patterns, but also for understanding the role of stratified media in modulating pattern-forming processes in developmental biology and beyond.
Collapse
Affiliation(s)
- Andrew L Krause
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG, UK.
| | - Václav Klika
- Department of Mathematics, FNSPE, Czech Technical University in Prague, Trojanova 13, 120 00, Prague, Czech Republic
| | - Jacob Halatek
- Microsoft Research, 21 Station Rd, Cambridge, CB1 2FB, UK
| | - Paul K Grant
- Microsoft Research, 21 Station Rd, Cambridge, CB1 2FB, UK
| | - Thomas E Woolley
- Cardiff School of Mathematics, Cardiff University, Senghennydd Road, Cardiff, CF24 4AG, UK
| | - Neil Dalchau
- Microsoft Research, 21 Station Rd, Cambridge, CB1 2FB, UK
| | - Eamonn A Gaffney
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG, UK
| |
Collapse
|
6
|
Krause AL, Klika V, Woolley TE, Gaffney EA. From one pattern into another: analysis of Turing patterns in heterogeneous domains via WKBJ. J R Soc Interface 2020; 17:20190621. [PMID: 31937231 DOI: 10.1098/rsif.2019.0621] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
Pattern formation from homogeneity is well studied, but less is known concerning symmetry-breaking instabilities in heterogeneous media. It is non-trivial to separate observed spatial patterning due to inherent spatial heterogeneity from emergent patterning due to nonlinear instability. We employ WKBJ asymptotics to investigate Turing instabilities for a spatially heterogeneous reaction-diffusion system, and derive conditions for instability which are local versions of the classical Turing conditions. We find that the structure of unstable modes differs substantially from the typical trigonometric functions seen in the spatially homogeneous setting. Modes of different growth rates are localized to different spatial regions. This localization helps explain common amplitude modulations observed in simulations of Turing systems in heterogeneous settings. We numerically demonstrate this theory, giving an illustrative example of the emergent instabilities and the striking complexity arising from spatially heterogeneous reaction-diffusion systems. Our results give insight both into systems driven by exogenous heterogeneity, as well as successive pattern forming processes, noting that most scenarios in biology do not involve symmetry breaking from homogeneity, but instead consist of sequential evolutions of heterogeneous states. The instability mechanism reported here precisely captures such evolution, and extends Turing's original thesis to a far wider and more realistic class of systems.
Collapse
Affiliation(s)
- Andrew L Krause
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Václav Klika
- Department of Mathematics, FNSPE, Czech Technical University in Prague, Trojanova 13, 120 00 Prague, Czech Republic
| | - Thomas E Woolley
- Cardiff School of Mathematics, Cardiff University, Senghennydd Road, Cardiff CF24 4AG, UK
| | - Eamonn A Gaffney
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| |
Collapse
|