Lizard Skin Patterns and the Ising Model

The Unreasonable Effectiveness of Reaction Diffusion in Vertebrate Skin Color Patterning

In 1952, Alan Turing published the reaction-diffusion (RD) mathematical framework, laying the foundations of morphogenesis as a self-organized process emerging from physicochemical first principles. Regrettably, this approach has been widely doubted in the field of developmental biology. First, we summarize Turing's line of thoughts to alleviate the misconception that RD is an artificial mathematical construct. Second, we discuss why phenomenological RD models are particularly effective for understanding skin color patterning at the meso/macroscopic scales, without the need to parameterize the profusion of variables at lower scales. More specifically, we discuss how RD models (a) recapitulate the diversity of actual skin patterns, (b) capture the underlying dynamics of cellular interactions, (c) interact with tissue size and shape, (d) can lead to ordered sequential patterning, (e) generate cellular automaton dynamics in lizards and snakes, (f) predict actual patterns beyond their statistical features, and (g) are robust to model variations. Third, we discuss the utility of linear stability analysis and perform numerical simulations to demonstrate how deterministic RD emerges from the underlying chaotic microscopic agents.

Colour patterns: Predicting patterns without knowing the details

A bewildering array of colour patterns occurs across animals. New work studying five lizard species reveals that reaction-diffusion models can be remarkably predictive of future adult skin patterning, even though molecular details are unknown. This has implications for understanding how complex patterns evolve.

Modeling convergent scale-by-scale skin color patterning in multiple species of lizards

Skin color patterning in vertebrates emerges at the macroscale from microscopic cell-cell interactions among chromatophores. Taking advantage of the convergent scale-by-scale skin color patterning dynamics in five divergent species of lizards, we quantify the respective efficiencies of stochastic (Lenz-Ising and cellular automata, sCA) and deterministic reaction-diffusion (RD) models to predict individual patterns and their statistical attributes. First, we show that all models capture the underlying microscopic system well enough to predict, with similar efficiencies, neighborhood statistics of adult patterns. Second, we show that RD robustly generates, in all species, a substantial gain in scale-by-scale predictability of individual adult patterns without the need to parametrize the system down to its many cellular and molecular variables. Third, using 3D numerical simulations and Lyapunov spectrum analyses, we quantitatively demonstrate that, given the non-linearity of the dynamical system, uncertainties in color measurements at the juvenile stage and in skin geometry variation explain most, if not all, of the residual unpredictability of adult individual scale-by-scale patterns. We suggest that the efficiency of RD is due to its intrinsic ability to exploit mesoscopic information such as continuous scale colors and the relations among growth, scales geometries, and the pattern length scale. Our results indicate that convergent evolution of CA patterning dynamics, leading to dissimilar macroscopic patterns in different species, is facilitated by their spontaneous emergence under a large range of RD parameters, as long as a Turing instability occurs in a skin domain with periodic thickness. VIDEO ABSTRACT.

Boundary fluctuation dynamics of a phase-separated domain in planar geometry

Using phase-ordering kinetics and of renormalization group theories, we derive analytically the relaxation times of the long wavelength fluctuations of a phase-separated domain boundary in the vicinity of (and below) the critical temperature in the planar Ising universality class. For a conserved order parameter, the relaxation time grows like Λ^{3} at wavelength Λ and can be expressed in terms of parameters relevant at the microscopic scale: lattice spacing, bulk diffusion coefficient of the minority phase, and temperature. These results are supported by numerical simulations of 2D Ising models, enabling in addition calculating the nonuniversal numerical prefactor. We discuss the applications of these findings to the determination of the real timescale associated with elementary Monte Carlo moves from the measurement of long wavelength relaxation times on experimental systems or molecular dynamics simulations.