Laboratory of Artificial & Natural Evolution

Vertebrate skin appendages are diverse micro-organs such as scales, feathers, and hair. These units typically develop from placodes, whose spatial patterning involves conserved chemical reaction-diffusion dynamics. Crocodile head scales are a spectacular exception to this paradigm, as they instead arise from a mechanically dominated process of compressive folding driven by constrained skin growth. Here, we reveal that chemical versus mechanical processes pattern tortoise scales in different regions of their head. Indeed, we show that placode-derived scales emerge across the peripheral head surfaces while remaining absent from the central dorsal region where scales subsequently form through a mechanical folding process. Using light sheet microscopy, we build a three-dimensional mechanical model that qualitatively recapitulates the diversity of scale patterns observed in this central head region in different tortoise species. Overall, our analyses indicate that mechanical head-scale patterning likely arose before the divergence between Testudinata and Archosauria, and was subsequently lost in birds.

Many examples of self-organized embryonic patterning can be attributed to chemically mediated systems comprising interacting morphogens. However, mechanical patterning also contributes to the emergence of biological forms. For example, various studies have demonstrated that diverse patterns arise from elastic instabilities associated with the constrained growth of soft materials, which generate wrinkles, creases and folds. Here, we show that between days 12 and 13 of development, transient experimentally increased activity of the sonic hedgehog pathway in the chicken embryo, through a single intravenous injection of smoothened agonist (SAG), abolishes the Turing-like chemical patterning of reticulate scales on the ventral footpad and promotes a transition to mechanical labyrinthine skin folding. Using hybridization, nanoindentation and labelling of proliferating cells, we confirm that skin surface folding is associated with the loss of signalling placode pre-patterning as well as increased epidermal growth and stiffness. Using additional hydrocortisone treatments, we also demonstrate that experimentally induced hyper-keratinization of the skin mechanically restricts SAG-induced folding. Finally, we verify our experimental findings with mechanical growth simulations built from volumetric light sheet fluorescence microscopy data. Overall, we reveal that pharmacological perturbation of the underlying gene regulatory network can abolish chemical skin appendage patterning and replace it with self-organized mechanical folding.

Snakes exhibit a broad variety of adaptive colors and color patterns, generated by the spatial arrangement of chromatophores, but little is known of the mechanisms responsible for these spectacular traits. Here, we investigate a mono-locus trait with two recessive alleles, motley and stripe, that both cause pattern aberrations in the corn snake.

The morphological intricacies of avian feathers make them an ideal model for investigating embryonic patterning and morphogenesis. In particular, the sonic hedgehog (Shh) pathway is an important mediator of feather outgrowth and branching. However, functional in vivo evidence regarding its role during feather development remains limited. Here, we demonstrate that an intravenous injection of sonidegib, a potent Shh pathway inhibitor, at embryonic day 9 (E9) temporarily produces striped domains (instead of spots) of Shh expression in the skin, arrests morphogenesis, and results in unbranched and non-invaginated feather buds-akin to proto-feathers-in embryos until E14. Although feather morphogenesis partially recovers, hatched treated chickens exhibit naked skin regions with perturbed follicles. Remarkably, these follicles are subsequently reactivated by seven weeks post-hatching. Our RNA-sequencing data and rescue experiment using Shh-agonism confirm that sonidegib specifically down-regulates Shh pathway activity. Overall, we provide functional evidence for the role of the Shh pathway in mediating feather morphogenesis and confirm its role in the evolutionary emergence and diversification of feathers.

Amniote integumentary appendages constitute a diverse group of micro-organs, including feathers, hair and scales. These structures typically develop as genetically controlled units1, the spatial patterning of which emerges from a self-organized chemical Turing system2,3 with integrated mechanical feedback4,5. The seemingly purely mechanical patterning of polygonal crocodile head scales provides an exception to this paradigm6. However, the nature and origin of the mechanical stress field driving this patterning remain unclear. Here, using precise in ovo intravenous injections of epidermal growth factor protein, we generate Nile crocodile embryos with substantially convoluted head skin, as well as hatchlings with smaller polygonal head scales resembling those of caimans. We then use light-sheet fluorescence microscopy to quantify embryonic tissue-layer geometry, collagen architecture and the spatial distribution of proliferating cells. Using these data, we build a phenomenological three-dimensional mechanical growth model that recapitulates both normal and experimentally modified patterning of crocodile head scales. Our experiments and numerical simulations demonstrate that crocodile head scales self-organize through compressive folding, originating from near-homogeneous skin growth with differential stiffness of the dermis versus the epidermis. Our experiments and theoretical morphospace analyses indicate that variation in embryonic growth and material properties of skin layers provides a simple evolutionary mechanism that produces a diversity of head-scale patterns among crocodilian species.

Amniote integumentary appendages constitute a diverse group of micro-organs, including feathers, hair and scales. These structures typically develop as genetically controlled units, the spatial patterning of which emerges from a self-organized chemical Turing system with integrated mechanical feedback. The seemingly purely mechanical patterning of polygonal crocodile head scales provides an exception to this paradigm. However, the nature and origin of the mechanical stress field driving this patterning remain unclear. Here, using precise in ovo intravenous injections of epidermal growth factor protein, we generate Nile crocodile embryos with substantially convoluted head skin, as well as hatchlings with smaller polygonal head scales resembling those of caimans. We then use light-sheet fluorescence microscopy to quantify embryonic tissue-layer geometry, collagen architecture and the spatial distribution of proliferating cells. Using these data, we build a phenomenological three-dimensional mechanical growth model that recapitulates both normal and experimentally modified patterning of crocodile head scales. Our experiments and numerical simulations demonstrate that crocodile head scales self-organize through compressive folding, originating from near-homogeneous skin growth with differential stiffness of the dermis versus the epidermis. Our experiments and theoretical morphospace analyses indicate that variation in embryonic growth and material properties of skin layers provides a simple evolutionary mechanism that produces a diversity of head-scale patterns among crocodilian species.

The glabrous skin of the rhinarium (naked nose) of many mammalian species exhibits a polygonal pattern of grooves that retain physiological fluid, thereby keeping their nose wet and, among other effects, facilitating the collection of chemosensory molecules. Here, we perform volumetric imaging of whole-mount rhinaria from sequences of embryonic and juvenile cows, dogs, and ferrets. We demonstrate that rhinarial polygonal domains are not placode-derived skin appendages but arise through a self-organized mechanical process consisting of the constrained growth and buckling of the epidermal basal layer, followed by the formation of sharp epidermal creases exactly facing an underlying network of stiff blood vessels. Our numerical simulations show that the mechanical stress generated by excessive epidermal growth concentrates at the positions of vessels that form rigid base points, causing the epidermal layers to move outward and shape domes-akin to arches rising against stiff pillars. Remarkably, this gives rise to a larger length scale (the distance between the vessels) in the surface folding pattern than would otherwise occur in the absence of vessels. These results hint at a concept of "mechanical positional information" by which material properties of anatomical elements can impose local constraints on an otherwise globally self-organized mechanical pattern. In addition, our analyses of the rhinarial patterns in cow clones highlight a substantial level of stochasticity in the pre-pattern of vessels, while our numerical simulations also recapitulate the disruption of the folding pattern in cows affected by a hereditary disorder that causes hyperextensibility of the skin.

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 () recapitulate the diversity of actual skin patterns, () capture the underlying dynamics of cellular interactions, () interact with tissue size and shape, () can lead to ordered sequential patterning, () generate cellular automaton dynamics in lizards and snakes, () predict actual patterns beyond their statistical features, and () 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.

Two influential concepts in tissue patterning are Wolpert's positional information and Turing's self-organized reaction-diffusion (RD). The latter establishes the patterning of hair and feathers. Here, our morphological, genetic, and functional-by CRISPR-Cas9-mediated gene disruption-characterization of wild-type versus "scaleless" snakes reveals that the near-perfect hexagonal pattern of snake scales is established through interactions between RD in the skin and somitic positional information. First, we show that ventral scale development is guided by hypaxial somites and, second, that ventral scales and epaxial somites guide the sequential RD patterning of the dorsolateral scales. The RD intrinsic length scale evolved to match somite periodicity, ensuring the alignment of ribs and scales, both of which play a critical role in snake locomotion.

We present a technique for precise drug delivery into the vascular system of developing amniote embryos via injection into chorioallantoic veins underlying the eggshell membrane. We describe steps for incubating and candling eggs, removing the shell to expose underlying veins, and precise intravenous injection. In addition to chicken embryos, this protocol is applicable to other amniote species that lay hard-shell eggs, including crocodiles and tortoises. This technique is rapid, is reproducible, is of low cost, and will provide an important resource for developmental biologists. For complete details on the use and execution of this protocol, please refer to Cooper & Milinkovitch..

Vertebrate skin appendage early development is mediated by conserved molecular signaling composing a dynamical reaction-diffusion-like system. Variations to such systems contribute to the remarkable diversity of skin appendage forms within and among species. Here, we demonstrate that stage-specific transient agonism of sonic hedgehog (Shh) pathway signaling in chicken triggers a complete and permanent transition from reticulate scales to feathers on the ventral surfaces of the foot and digits. Resulting ectopic feathers are developmentally comparable to feathers adorning the body, with down-type feathers transitioning into regenerative, bilaterally symmetric contour feathers in adult chickens. Crucially, this spectacular transition of skin appendage fate (from nodular reticulate scales to bona fide adult feathers) does not require sustained treatment. Our RNA sequencing analyses confirm that smoothened agonist treatment specifically promotes the expression of key Shh pathway-associated genes. These results indicate that variations in Shh pathway signaling likely contribute to the natural diversity and regionalization of avian integumentary appendages.

Although notoriously difficult, imaging collagen network architecture, a key element affecting tissue mechanical properties, is of paramount importance in developmental and cancer biology. Here, we introduce a simple and robust method of whole-mount collagen staining with the 'Fast Green' dye that provides unmatched visualization of collagen 3D network architecture, via confocal or light-sheet microscopy, compatible with solvent-based tissue clearing and immunostaining.

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.

In this paper we propose a new mathematical model for describing the complex interplay between skin cell populations with fibroblast growth factor and bone morphogenetic protein, occurring within deformable porous media describing feather primordia patterning. Tissue growth, in turn, modifies the transport of morphogens (described by reaction-diffusion equations) through diverse mechanisms such as advection from the solid velocity generated by mechanical stress, and mass supply. By performing an asymptotic linear stability analysis on the coupled poromechanical-chemotaxis system (assuming rheological properties of the skin cell aggregates that reside in the regime of infinitesimal strains and where the porous structure is fully saturated with interstitial fluid and encoding the coupling mechanisms through active stress) we obtain the conditions on the parameters—especially those encoding coupling mechanisms—under which the system will give rise to spatially heterogeneous solutions. We also extend the mechanical model to the case of incompressible poro-hyperelasticity and include the mechanisms of anisotropic solid growth and feedback by means of standard Lee decompositions of the tensor gradient of deformation. Because the model in question involves the coupling of several nonlinear PDEs, we cannot straightforwardly obtain closed-form solutions. We therefore design a suitable numerical method that employs backward Euler time discretisation, linearisation of the semidiscrete problem through Newton–Raphson’s method, a seven-field finite element formulation for the spatial discretisation, and we also advocate the construction and efficient implementation of tailored robust solvers. We present a few illustrative computational examples in 2D and 3D, briefly discussing different spatio-temporal patterns of growth factors as well as the associated solid response scenario depending on the specific poromechanical regime. Our findings confirm the theoretically predicted behaviour of spatio-temporal patterns, and the produced results reveal a qualitative agreement with respect to the expected experimental behaviour. We stress that the present study provides insight on several biomechanical properties of primordia patterning.

The ocellated lizard (Timon lepidus) exhibits an intricate skin color pattern made of monochromatic black and green skin scales, whose dynamics of color flipping are known to be well modeled by a stochastic cellular automaton. We show that the late-time probability distribution of the pattern corresponds to the canonical probability distribution of the antiferromagnetic Ising model and can be generated by dynamics different from the commonly-used Glauber. We comment on skin scale patterns generated by the Ising model on the triangular lattice in the low-temperature limit.

The elephant proboscis (trunk), which functions as a muscular hydrostat with a virtually infinite number of degrees of freedom, is a spectacular organ for delicate to heavy object manipulation as well as social and sensory functions. Using high-resolution motion capture and functional morphology analyses, we show here that elephants evolved strategies that reduce the biomechanical complexity of their trunk. Indeed, our behavioral experiments with objects of various shapes, sizes, and weights indicate that (1) complex behaviors emerge from the combination of a finite set of basic movements; (2) curvature, torsion, and strain provide an appropriate kinematic representation, allowing us to extract motion primitives from the trunk trajectories; (3) transport of objects involves the proximal propagation of an inward curvature front initiated at the tip; (4) the trunk can also form pseudo-joints for point-to-point motion; and (5) the trunk tip velocity obeys a power law with its path curvature, similar to human hand drawing movements. We also reveal with unprecedented precision the functional anatomy of the African and Asian elephant trunks using medical imaging and macro-scale serial sectioning, thus drawing strong connections between motion primitives and muscular synergies. Our study is the first combined quantitative analysis of the mechanical performance, kinematic strategies, and functional morphology of the largest animal muscular hydrostat on Earth. It provides data for developing innovative "soft-robotic" manipulators devoid of articulations, replicating the high compliance, flexibility, and strength of the elephant trunk. VIDEO ABSTRACT.

We previously showed that the adult ocellated lizard skin colour pattern is effectively generated by a stochastic cellular automaton (CA) of skin scales. We additionally suggested that the canonical continuous 2D reaction-diffusion (RD) process of colour pattern development is transformed into this discrete CA by reduced diffusion coefficients at the borders of scales (justified by the corresponding thinning of the skin). Here, we use RD numerical simulations in 3D on realistic lizard skin geometries and demonstrate that skin thickness variation on its own is sufficient to cause scale-by-scale coloration and CA dynamics during RD patterning. In addition, we show that this phenomenon is robust to RD model variation. Finally, using dimensionality-reduction approaches on large networks of skin scales, we show that animal growth affects the scale-colour flipping dynamics by causing a substantial decrease of the relative length scale of the labyrinthine colour pattern of the lizard skin.

In this issue of Developmental Cell, Chuyen et al. (2021) suggest that the Scf/Kit pathway controls mutual repulsion of multiciliated cells and their affinity for epidermal cell junctions through soluble and membrane-associated Scf ligands, respectively. Effective self-organizational patterning emerges at the mesoscopic scale as a small set of effective behaviors.

Reptiles exhibit a spectacular diversity of skin colors and patterns brought about by the interactions among three chromatophore types: black melanophores with melanin-packed melanosomes, red and yellow xanthophores with pteridine- and/or carotenoid-containing vesicles, and iridophores filled with light-reflecting platelets generating structural colors. Whereas the melanosome, the only color-producing endosome in mammals and birds, has been documented as a lysosome-related organelle, the maturation paths of xanthosomes and iridosomes are unknown. Here, we first use 10x Genomics linked-reads and optical mapping to assemble and annotate a nearly chromosome-quality genome of the corn snake The assembly is 1.71 Gb long, with an N50 of 16.8 Mb and L50 of 24. Second, we perform mapping-by-sequencing analyses and identify a 3.9-Mb genomic interval where the variant resides. The lavender color morph in corn snakes is characterized by gray, rather than red, blotches on a pink, instead of orange, background. Third, our sequencing analyses reveal a single nucleotide polymorphism introducing a premature stop codon in the lysosomal trafficking regulator gene () that shortens the corresponding protein by 603 amino acids and removes evolutionary-conserved domains. Fourth, we use light and transmission electron microscopy comparative analyses of wild type versus lavender corn snakes and show that the color-producing endosomes of all chromatophores are substantially affected in the mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

Females of some Asian salamanders of the genus Hynobius deposit in streams their eggs embedded in a translucent envelope called an 'egg sac'. The edges of the envelope exhibit a spectacular blue-to-yellow iridescent glow, which instantaneously disappears when the sac is removed from water. First, our scanning electron microscopy analyses reveal that the inner surface of the 100 μm-thick envelope displays striations (length scale of about 3 μm), which are themselves covered by much smaller (190 ± 30 nm) and quasi-periodic corrugations. The latter could constitute a surface diffraction grating generating iridescence by light interference. Second, our transmission electron microscopy and focused-ion-beam scanning electron microscopy analyses show that the bulk of the egg sac wall is composed of meandering fibres with a quasi-periodic modulation of 190 ± 60 nm along the thickness of the envelope, generating a photonic crystal. Third, Fourier power analyses of 450 electron microscopy images with varying incident angles indicate that changing the surrounding medium from water to air shifts most of the backscattered power spectrum to the ultraviolet range, hence, explaining that the egg sac loses visible iridescence when removed out of the water. Fourth, the results of our photography and optical spectroscopy experiments of submerged and emerged egg sacs rule out the possibility that the iridescence is due to a thin film or a multilayer, whereas the observed non-specular response is compatible with the backscattering expected from surface diffraction gratings and volumetric photonic crystals with spatial 1D modulation. Finally, although we mention several potential biological functions of the egg sac structural colours and iridescence, we emphasise that these optical properties might be the by-products of the envelope material internal structure selected during evolution for its mechanical properties.