Spatial pattern formation in reaction-diffusion models: a computational approach. (English) Zbl 1432.92009
Summary: Reaction-diffusion equations have been widely used to describe biological pattern formation. Nonuniform steady states of reaction-diffusion models correspond to stationary spatial patterns supported by these models. Frequently these steady states are not unique and correspond to various spatial patterns observed in biology. Traditionally, time-marching methods or steady state solvers based on Newton’s method were used to compute such solutions. However, the solutions that these methods converge to highly depend on the initial conditions or guesses. In this paper, we present a systematic method to compute multiple nonuniform steady states for reaction-diffusion models and determine their dependence on model parameters. The method is based on homotopy continuation techniques and involves mesh refinement, which significantly reduces computational cost. The method generates one-parameter steady state bifurcation diagrams that may contain multiple unconnected components, as well as two-parameter solution maps that divide the parameter space into different regions according to the number of steady states. We applied the method to two classic reaction-diffusion models and compared our results with available theoretical analysis in the literature. The first is the Schnakenberg model which has been used to describe biological pattern formation due to diffusion-driven instability. The second is the Gray-Scott model which was proposed in the 1980s to describe autocatalytic glycolysis reactions. In each case, the method uncovers many, if not all, nonuniform steady states and their stabilities.
MSC:
| 92C15 | Developmental biology, pattern formation |
| 35Q92 | PDEs in connection with biology, chemistry and other natural sciences |
| 35K57 | Reaction-diffusion equations |
References:
| [1] | Baker, Re; Schnell, S.; Maini, Pk, Waves and patterning in developmental biology: vertebrate segmentation and feather bud formation as case studies, Int J Dev Biol, 53, 783 (2009) |
| [2] | Bates, Dj; Hauenstein, Jd; Sommese, Aj, Efficient path tracking methods, Numer Algorithms, 58, 4, 451-459 (2011) · Zbl 1230.65059 |
| [3] | Bates, Dj; Hauenstein, Jd; Sommese, Aj; Wampler, Ii; Charles, W., Adaptive multiprecision path tracking, SIAM J Numer Anal, 46, 2, 722-746 (2008) · Zbl 1162.65026 |
| [4] | Ben-Jacob, E.; Cohen, I.; Levine, H., Cooperative self-organization of microorganisms, Adv Phys, 49, 4, 395-554 (2000) |
| [5] | Briggs, Wl; Mccormick, Sf, A multigrid tutorial (2000), University City: SIAM, University City |
| [6] | Cross, M.; Hohenberg, P., Pattern formation outside of equilibrium, Rev Mod Phys, 65, 3, 851 (1993) · Zbl 1371.37001 |
| [7] | Ermentrout, B., Simulating, analyzing, and animating dynamical systems: a guide to XPPAUT for researchers and students (2002), University City: SIAM, University City · Zbl 1003.68738 |
| [8] | Farrell, Pe; Birkisson, A.; Funke, Sw, Deflation techniques for finding distinct solutions of nonlinear partial differential equations, SIAM J Sci Comput, 37, 4, A2026-A2045 (2015) · Zbl 1327.65237 |
| [9] | Freiling, G.; Yurko, V., Inverse Sturm-Liouville problems and their applications (2001), New York: NOVA Science Publishers, New York · Zbl 1037.34005 |
| [10] | Gierer, A.; Meinhardt, H., A theory of biological pattern formation, Kybernetik, 12, 1, 30-39 (1972) · Zbl 1434.92013 |
| [11] | Gray, P.; Scott, Sk, Autocatalytic reactions in the isothermal continuous stirred tank reactor: isolas and other forms of multistability, Chem Eng Sci, 38, 1, 29-43 (1983) |
| [12] | Gray, P.; Scott, Sk, Autocatalytic reactions in the isothermal continuous stirred tank reactor: oscillations and instabilities in the system \({A}+2{B} \rightarrow 3{B};\, {B}\rightarrow{C} \), Chem Eng Sci, 39, 6, 1087-1097 (1984) |
| [13] | Gray, P.; Scott, Sk, Sustained oscillations and other exotic patterns of behavior in isothermal reactions, J Phys Chem, 89, 22-32 (1985) |
| [14] | Hao, W.; Hauenstein, J.; Hu, B.; Sommese, A., A bootstrapping approach for computing multiple solutions of differential equations, J Comput Appl Math, 258, 181-190 (2014) · Zbl 1294.65085 |
| [15] | Hillen, T.; Painter, Kj, A user’s guide to PDE models for chemotaxis, J Math Biol, 58, 1-2, 183-217 (2009) · Zbl 1161.92003 |
| [16] | Iron, D.; Wei, J.; Winter, M., Stability analysis of Turing patterns generated by the Schnakenberg model, J Math Biol, 49, 4, 358-390 (2004) · Zbl 1057.92011 |
| [17] | Jilkine, A.; Edelstein-Keshet, L., A comparison of mathematical models for polarization of single eukaryotic cells in response to guided cues, PLoS Comput Biol, 7, 4, e1001121 (2011) |
| [18] | Kaper, Hg; Wang, S.; Yari, M., Dynamical transitions of turing patterns, Nonlinearity, 22, 3, 601 (2009) · Zbl 1157.92003 |
| [19] | Kelley, Ct, Iterative methods for linear and nonlinear equations, Front Appl Math, 16, 575-601 (1995) · Zbl 0832.65046 |
| [20] | Koch, Aj; Meinhardt, H., Biological pattern formation: from basic mechanisms to complex structures, Rev Mod Phys, 66, 4, 1481 (1994) |
| [21] | Kondo, S.; Miura, T., Reaction-diffusion model as a framework for understanding biological pattern formation, Science, 329, 5999, 1616-1620 (2010) · Zbl 1226.35077 |
| [22] | Kong, Q.; Zettl, A., Eigenvalues of regular Sturm-Liouville problems, J Differ Equ, 131, 1, 1-19 (1996) · Zbl 0862.34020 |
| [23] | Lee, Kj; Mccormick, Wd; Ouyang, Q.; Swinney, Hl, Pattern formation by interacting chemical fronts, Science, 261, 5118, 192-194 (1993) |
| [24] | Liu, C.; Fu, X.; Liu, L.; Ren, X.; Chau, Ckl; Li, S.; Xiang, L.; Zeng, H.; Chen, G.; Tang, L., Sequential establishment of stripe patterns in an expanding cell population, Science, 334, 6053, 238-241 (2011) |
| [25] | Lo, W-C; Chen, L.; Wang, M.; Nie, Q., A robust and efficient method for steady state patterns in reaction-diffusion systems, J Comput Phys, 231, 15, 5062-5077 (2012) · Zbl 1250.35017 |
| [26] | Maini, Pk; Woolley, Te; Baker, Re; Gaffney, Ea; Lee, Ss, Turing’s model for biological pattern formation and the robustness problem, Interface Focus, 2, 4, 487-496 (2012) |
| [27] | Murray, Jd, Mathematical biology (2002), Berlin: Springer, Berlin · Zbl 1006.92001 |
| [28] | Othmer, Hg; Painter, Kj; Umulis, D.; Xue, C., The intersection of theory and application in elucidating pattern formation in developmental biology, Math Model Nat Phenom, 4, 4, 3-82 (2009) · Zbl 1181.35297 |
| [29] | Othmer, Hg; Stevens, A., Aggregation, blowup, and collapse: the ABC’s of taxis in reinforced random walks, SIAM J Appl Math, 57, 4, 1044-1081 (1997) · Zbl 0990.35128 |
| [30] | Painter, Kj; Maini, Pk; Othmer, Hg, Stripe formation in juvenile Pomacanthus explained by a generalized Turing mechanism with chemotaxis, Proc Natl Acad Sci USA, 96, 10, 5549-5554 (1999) |
| [31] | Painter, Kj; Hillen, T., Spatio-temporal chaos in a chemotaxis model, Physica D, 240, 4-5, 363-375 (2011) · Zbl 1255.37026 |
| [32] | Pearson, Je, Complex patterns in a simple system, Science, 261, 5118, 189-07 (1993) |
| [33] | Robinson, M.; Luo, C.; Farrell, Pe; Erban, R.; Majumdar, A., From molecular to continuum modelling of bistable liquid crystal devices, Liq Cryst, 44, 14-15, 2267-2284 (2017) |
| [34] | Sommese, A.; Wampler, C., The numerical solution of systems of polynomials arising in engineering and science (2005), Singapore: World Scientific, Singapore · Zbl 1091.65049 |
| [35] | Stoodley, P.; Sauer, K.; Davies, Dg; Costerton, Jw, Biofilms as complex differentiated communities, Ann Rev Microbiol, 56, 1, 187-209 (2002) |
| [36] | Sun, W.; Ward, Mj; Russell, R., The slow dynamics of two-spike solutions for the Gray-Scott and Gierer-Meinhardt systems: competition and oscillatory instabilities, SIAM J Appl Dyn Syst, 4, 4, 904-953 (2005) · Zbl 1145.35404 |
| [37] | Thomas, Jw, Numerical partial differential equations: finite difference methods (2013), Berlin: Springer, Berlin |
| [38] | Turing, Am, The chemical basis of morphogenesis, Philos Trans R Soc Lond B, 237, 641, 37-72 (1952) · Zbl 1403.92034 |
| [39] | Uecker, H.; Wetzel, D.; Rademacher, J., pde2path-A Matlab package for continuation and bifurcation in 2D elliptic systems, Numer Math Theory Methods Appl, 7, 1, 58-106 (2014) · Zbl 1313.65311 |
| [40] | Volkening, A.; Sandstede, B., Modelling stripe formation in zebrafish: an agent-based approach, J R Soc Interface, 12, 112, 20150812 (2015) |
| [41] | Wang Q, Oh JW, Lee H-L, Dhar A, Peng T, Ramos R, Guerrero-Juarez CF, Wang X, Zhao R, Cao X et al (2017) A multi-scale model for hair follicles reveals heterogeneous domains driving rapid spatiotemporal hair growth patterning. eLife 6:e22772 |
| [42] | Wei, J., Existence and stability of spikes for the Gierer-Meinhardt system, Handb Differ Equ Station Par Differ Equ, 5, 487-585 (2008) · Zbl 1223.35007 |
| [43] | Wei, J.; Winter, M., Stationary multiple spots for reaction-diffusion systems, J Math Biol, 57, 1, 53-89 (2008) · Zbl 1141.92007 |
| [44] | Wilkinson, Jh, Rounding errors in algebraic processes (1994), Chelmsford: Courier Corporation, Chelmsford · Zbl 0868.65027 |
| [45] | Wollkind, Dj; Manoranjan, Vs; Zhang, L., Weakly nonlinear stability analyses of prototype reaction-diffusion model equations, SIAM Rev, 36, 2, 176-214 (1994) · Zbl 0808.35056 |
| [46] | Xu, J., A novel two-grid method for semilinear elliptic equations, SIAM J Sci Comput, 15, 1, 231-237 (1994) · Zbl 0795.65077 |
| [47] | Xue, C., Macroscopic equations for bacterial chemotaxis: integration of detailed biochemistry of cell signaling, J Math Biol, 70, 1-2, 1-44 (2015) · Zbl 1375.92012 |
| [48] | Xue, C.; Budrene, Eo; Othmer, Hg, Radial and spiral stream formation in Proteus mirabilis colonies, PLoS Comput Biol, 7, 12, 12 (2011) |
| [49] | Xue, C.; Othmer, Hg, Multiscale models of taxis-driven patterning in bacterial populations, SIAM J Appl Math, 70, 1, 133-167 (2009) · Zbl 1184.35308 |
| [50] | Xue, C.; Shtylla, B.; Brown, A., A stochastic multiscale model that explains the segregation of axonal microtubules and neurofilaments in toxic neuropathies, PLoS Comput Biol, 11, e1004406 (2015) |
| [51] | Xue, X.; Xue, C.; Tang, M., The role of intracellular signaling in the stripe formation in engineered escherichia coli populations, PLoS Comput Biol, 14, 6, 1-23 (2018) |
| [52] | Zhao, S.; Ovadia, J.; Liu, X.; Zhang, Y-T; Nie, Q., Operator splitting implicit integration factor methods for stiff reaction-diffusion-advection systems, J Comput Phys, 230, 15, 5996-6009 (2011) · Zbl 1220.65120 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
