Modeling the morphodynamic galectin patterning network of the developing avian limb skeleton

We present a mathematical model for the morphogenesis and patterning of the mesenchymal condensations that serve as primordia of the avian limb skeleton. The model is based on the experimentally established dynamics of a multiscale regulatory network consisting of two glycan-binding proteins expressed early in limb development: CG (chicken galectin)-1A, CG-8 and their counterreceptors that determine the formation, size, number and spacing of the "protocondensations" that give rise to the condensations and subsequently the cartilaginous elements that serve as the templates of the bones. The model, a system of partial differential and integro-differential equations containing a flux term to represent local adhesion gradients, is simulated in a "full" and a "reduced" form to confirm that the system has pattern-forming capabilities and to explore the nature of the patterning instability. The full model recapitulates qualitatively and quantitatively the experimental results of network perturbation and leads to new predictions, which are verified by further experimentation. The reduced model is used to demonstrate that the patterning process is inherently morphodynamic, with cell motility being intrinsic to it. Furthermore, subtle relationships between cell movement and the positive and negative interactions between the morphogens produce regular patterns without the requirement for activators and inhibitors with widely separated diffusion coefficients. The described mechanism thus represents an extension of the category of activator-inhibitor processes capable of generating biological patterns with repetitive elements beyond the morphostatic mechanisms of the Turing/Gierer-Meinhardt type.