The Morphostatic Limit for a Model of Skeletal Pattern Formation in the Vertebrate Limb

Abstract A recently proposed mathematical model of a “core” set of cellular and molecular interactions present in the developing vertebrate limb was shown to exhibit pattern-forming instabilities and limb skeleton-like patterns under certain restrictive conditions, suggesting that it may authentically represent the underlying embryonic process (Hentschel et al., Proc. R. Soc. B 271, 1713–1722, 2004). The model, an eight-equation system of partial differential equations, incorporates the behavior of mesenchymal cells as “reactors,” both participating in the generation of morphogen patterns and changing their state and position in response to them. The full system, which has smooth solutions that exist globally in time, is nonetheless highly complex and difficult to handle analytically or numerically. According to a recent classification of developmental mechanisms (Salazar-Ciudad et al., Development 130, 2027–2037, 2003), the limb model of Hentschel et al. is “morphodynamic,” since differentiation of new cell types occurs simultaneously with cell rearrangement. This contrasts with “morphostatic” mechanisms, in which cell identity becomes established independently of cell rearrangement. Under the hypothesis that development of some vertebrate limbs employs the core mechanism in a morphostatic fashion, we derive in an analytically rigorous fashion a pair of equations representing the spatiotemporal evolution of the morphogen fields under the assumption that cell differentiation relaxes faster than the evolution of the overall cell density (i.e., the morphostatic limit of the full system). This simple reaction–diffusion system is unique in having been derived analytically from a substantially more complex system involving multiple morphogens, extracellular matrix deposition, haptotaxis, and cell translocation. We identify regions in the parameter space of the reduced system where Turing-type pattern formation is possible, which we refer to as its “Turing space.” Obtained values of the parameters are used in numerical simulations of the reduced system, using a new Galerkin finite element method, in tissue domains with nonstandard geometry. The reduced system exhibits patterns of spots and stripes like those seen in developing limbs, indicating its potential utility in hybrid continuum-discrete stochastic modeling of limb development. Lastly, we discuss the possible role in limb evolution of selection for increasingly morphostatic developmental mechanisms.