Somitic positional information guides self-organized patterning of snake scales

Cyclic renewal in three ectodermal appendage follicles: Hairs, feathers and teeth

Ectodermal appendages display a range of renewal mechanisms, with some undergoing continuous growth and others experiencing cyclic regeneration. The latter requires sustainable epithelial stem cells and mesenchymal niche essential for interacting with these stem cells. Furthermore, certain appendages dynamically adjust their mesenchymal niche in response to environmental factors, such as hormonal fluctuations, sex, and seasonal changes, enabling them to cyclically renew with different appendages phenotypes to adapt to different environments and to different life stages. Here we focus on amniotes, including reptiles, birds, and mammals, which exhibit integumentary adaptations that enable their survival across various ecological environments, from aquatic habitats and terrestrial landscapes to aerial domains. We highlight three representative integument appendage follicles: teeth, feathers, and hairs. Despite independent evolutionary origins, these structures share a fundamental architectural design characterized by the presence of stem cells and mesenchymal niches. They differ in the spatial arrangement and topology of these components. By examining the distinct architectural features of these follicles, we demonstrate the different strategies they use to orchestrate the physiological regenerative cycling, from growth initiation to cessation and molting, and regeneration after wounding. We delve into known molecular controls that govern these processes and unravel the evolutionary insights. We also identify new cell interactions that underlie the emergence of evolutionary novel follicle components. Various amniote scales have evolved independently with different configurations, but all lack follicle architecture and maintain homeostasis using a strategy similar to that of skin. The convergently evolved follicles in hairs, feathers, and teeth utilize different designs to achieve cyclic renewability, allowing them to produce spatially and temporally specific appendage phenotypes, thus enhancing the adaptability of the integumentary interface to external environmental pressures. This, in turn, enriches our understanding of evolutionary developmental biology (Evo-Devo) of the integument, shedding light on the intricate interplay between form and function across diverse taxa.

Organizational principles of integumentary organs: Maximizing variations for effective adaptation

The integument serves as the interface between an organism and its environment. It primarily comprises ectoderm-derived epithelium and mesenchyme derived from various embryonic sources. These integumentary organs serve as a barrier defining the physiological boundary between the internal and exterior environments and fulfill diverse functions. How does the integument generate such a large diversity? Here, we attempt to decipher the organizational principles. We focus on amniotes and use appendage follicles as the primary examples. The integument begins as a simple planar sheet of coupled epithelial and mesenchymal cells, then becomes more complex through the following patterning processes. 1) De novo Turing periodic patterning process: This process converts the integument into multiple skin appendage units. 2) Adaptive patterning process: Dermal muscle, blood vessels, adipose tissue, and other components are assembled and organized around appendage follicles when present. 3) Cyclic renewal: Skin appendage follicles contain stem cells and their niches, enabling physiological molting and regeneration in the adult animal. 4) Spatial variations: Multiple appendage units allow modulation of shape, size, keratin types, and color patterns of feathers and hairs across the animal's surface. 5) Temporal phenotypic plasticity: Cyclic renewal permits temporal transition of appendage phenotypes, i.e. regulatory patterning or integumentary metamorphosis, throughout an animal's lifetime. The diversities in (4) and (5) can be generated epigenetically within the same animal. Over the evolutionary timescale, different species can modulate the number, size, and distributions of existing ectodermal organs in the context of micro-evolution, allowing effective adaptation to new climates as seen in the variation of hair length among mammals. Novel ectodermal organs can also emerge in the context of macro-evolution, enabling animals to explore new ecological niches, as seen in the emergence of feathers on dinosaurs. These principles demonstrate how multi-scale organ adaption in the amniotes can maximize diverse and flexible integumentary organ phenotypes, producing a vast repertoire for natural selection and thereby providing effective adaptation and evolutionary advantages.

Exacerbated sonic hedgehog signalling promotes a transition from chemical pre-patterning of chicken reticulate scales to mechanical skin folding

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 in ovo 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 in situ 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 in ovo 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.

Regulatory and disruptive variants in the CLCN2 gene are associated with modified skin color pattern phenotypes in the corn snake

BACKGROUND

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.

RESULTS

We use mapping-by-sequencing to identify the genomic interval where the causal mutations reside. With our differential gene expression analyses, we find that CLCN2 (Chloride Voltage-Gated Channel 2), a gene within the genomic interval, is significantly downregulated in Motley embryonic skin. Furthermore, we identify the stripe allele as the insertion of an LTR-retrotransposon in CLCN2, resulting in a disruptive mutation of the protein. We confirm the involvement of CLCN2 in color pattern formation by producing knock-out snakes that present a phenotype similar to Stripe. In humans and mice, disruption of CLCN2 results in leukoencephalopathy, as well as retinal and testes degeneration. Our single-cell transcriptomic analyses in snakes reveal that CLCN2 is indeed expressed in chromatophores during embryogenesis and in the adult brain, but the behavior and fertility of Motley and Stripe corn snakes are not impacted.

CONCLUSIONS

Our genomic, transcriptomic, and functional analyses identify a plasma membrane anion channel to be involved in color pattern development in snakes and show that an active LTR-retrotransposon might be a key driver of trait diversification in corn snakes.

In vivo sonic hedgehog pathway antagonism temporarily results in ancestral proto-feather-like structures in the chicken

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.

Tending to the Facial Surfaces of a Mathematical Biology Head-Scratcher: Why Does the Head of the Sea Turtle Natator depressus Resemble a Convex Zygomorphic Dodecahedron?

Two convex polyhedra that markedly resemble the head of the flatback sea turtle hatchling are identified. The first example is a zygomorphic tetragonal dodecahedron, while the other, an even better matching structure, is a related tetradecahedron, herein speculated to arise from this particular dodecahedron via known mechanisms gleaned from studies of the behavior of foams. A segmented, biomorphic, convex polyhedral model to address cephalic topology is thus presented stemming from solid geometry, anatomical observations, and a recently computed densest local packing arrangement of fifteen slightly oblate spheroids in which fourteen oblate spheroids surround a central such spheroid. This particular array of oblate spheroids shares salient structural features with the aforementioned dodecahedron. Successful testing of the model has been achieved by converting this array of fifteen oblate spheroids constructed with putty to the cephaloid dodecahedron in a process involving ventral elongation induced by stretching in the anterior direction along the anteroposterior axis (convergent extension). During convergent extension, the two left most anterolateral oblate spheroids that are in direct contact with the ventral spheroid of the array merge into a single lateral facet of the incipient dodecahedron, while the corresponding two right such oblate spheroids do the same. Thus, the fourteen outer oblate spheroids of the array give rise to the twelve facets of the finalized dodecahedron, while the central oblate spheroid remnant assumes an interior dodecahedral position. The hypothetical dodecahedron to tetradecahedron transformation entails the collapse of a tetravalent vertex (which is known to occur in foams as part of a T1 transition) followed by bilateral facet splitting. Remarkably, a model stipulating that convexity is to be retained in connection with this sequence of steps necessitates that the starting dodecahedral template undergoes modification to become a tetradecahedron in possession of precisely the highly ordered feature found at the top of the head of numerous specimens of the flatback sea turtle hatchling, namely, a fused medial pentagon-heptagon pair in the form of a pentagonal frontal scute and heptagonal frontoparietal scute. Such a possible new instance of geometric biomorphy, taken together with the correct anticipation of the cephalic pentagon-heptagon pair, might serve to instill further confidence in renewed efforts to shed light on morphogenesis with foam embryo models.

Self-organized patterning of crocodile head scales by compressive folding

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.

A Ball Python Colour Morph Implicates MC1R in Melanophore-Xanthophore Distribution and Pattern Formation

Reptiles showcase an extensive array of skin colours and patterns, yet little is known about the genetics of reptile colouration. Here, we investigate the genetic basis of the Clown colour morph found in captive-bred ball pythons (Python regius) to study skin pigmentation and patterning in snakes. We obtained samples by crowdsourcing shed skin from commercial breeders and hobbyists. We applied a case-control design, whole-genome pool sequencing, variant annotation, histological analyses, and electron microscopy imaging. We identified a missense mutation in a transmembrane region of the melanocortin-1 receptor (MC1R) associated with the Clown phenotype. In classic avian and mammalian model species, MC1R is known for controlling the type and amount of melanin produced. In contrast, our results suggest that MC1R signalling might play a key role in pattern formation in ball pythons, affecting xanthophore-melanophore distribution. This work highlights the varied functions of MC1R across different vertebrate lineages and promotes a novel model system to study reptile colouration.

A new cell culture resource for investigations of reptilian gene function

The establishment of CRISPR/Cas9 gene editing in Anolis sagrei has positioned this species as a powerful model for studies of reptilian gene function. To enhance this model, we developed an immortalized lizard fibroblast cell line (ASEC-1) for the exploration of reptilian gene function in cellular processes. We demonstrate the use of this cell line by scrutinizing the role of primary cilia in lizard Hedgehog (Hh) signaling. Using CRISPR/Cas9 mutagenesis, we disrupted the ift88 gene, which is required for ciliogenesis in diverse organisms. We determined that loss of itf88 from lizard cells leads to an absence of primary cilia, a partial derepression of gli1 transcription, and an inability of the cells to respond to the Smoothened agonist, SAG. Through a cross-species analysis of SAG-induced transcriptional responses in cultured limb bud cells, we further determined that ∼46% of genes induced as a response to Hh pathway activation in A. sagrei are also SAG responsive in Mus musculus limb bud cells. Our results highlight conserved and diverged aspects of Hh signaling in anoles and establish a new resource for investigations of reptilian gene function.

Micropattern Fabricated by Acropetal Migration Controlled through Sequential Photo and Thermal Polymerization

Bottom-up patterning technology plays a significant role in both nature and synthetic materials, owing to its inherent advantages such as ease of implementation, spontaneity, and noncontact attributes, etc. However, constrained by the uncontrollability of molecular movement, energy interaction, and stress, obtained micropatterns tend to exhibit an inevitable arched outline, resulting in the limitation of applicability. Herein, inspired by auxin's action mode in apical dominance, a versatile strategy is proposed for fabricating precision self-organizing micropatterns with impressive height based on polymerization-induced acropetal migration. The copolymer containing fluorocarbon chains (low surface energy) and tertiary amine (coinitiator) is designed to self-assemble on the surface of the photo-curing system. The selective exposure under a photomask establishes a photocuring boundary and the radicals would be generated on the surface, which is pivotal in generating a vertical concentration difference of monomer. Subsequent heating treatment activates the material continuously transfers from the unexposed area to the exposed area and is accompanied by the obviously vertical upward mass transfer, resulting in the manufacture of a rectilinear profile micropattern. This strategy significantly broadens the applicability of self-organizing patterns, offering the potential to mitigate the complexity and time-consuming limitations associated with top-down methods.

PMEL is involved in snake colour pattern transition from blotches to stripes

Corn snakes are emerging models for animal colouration studies. Here, we focus on the Terrazzo morph, whose skin pattern is characterized by stripes rather than blotches. Using genome mapping, we discover a disruptive mutation in the coding region of the Premelanosome protein (PMEL) gene. Our transcriptomic analyses reveal that PMEL expression is significantly downregulated in Terrazzo embryonic tissues. We produce corn snake PMEL knockouts, which present a comparable colouration phenotype to Terrazzo and the subcellular structure of their melanosomes and xanthosomes is also similarly impacted. Our single-cell expression analyses of wild-type embryonic dorsal skin demonstrate that all chromatophore progenitors express PMEL at varying levels. Finally, we show that in wild-type embryos PMEL-expressing cells are initially uniformly spread before forming aggregates and eventually blotches, as seen in the adults. In Terrazzo embryos, the aggregates fail to form. Our results provide insights into the mechanisms governing colouration patterning in reptiles.

The dynamic behavior of chromatophores marks the transition from bands to spots in leopard geckos

Reptilian skin coloration is spectacular and diverse, yet little is known about the ontogenetic processes that govern its establishment and the molecular signaling pathways that determine it. Here, we focus on the development of the banded pattern of leopard gecko hatchlings and the transition to black spots in the adult. With our histological analyses, we show that iridophores are present in the white and yellow bands of the hatchling and they gradually perish in the adult skin. Furthermore, we demonstrate that melanophores can autonomously form spots in the absence of the other chromatophores both on the regenerated skin of the tail and on the dorsal skin of the Mack Super Snow (MSS) leopard geckos. This color morph is characterized by uniform black coloration in hatchlings and black spots in adulthood; we establish that their skin is devoid of xanthophores and iridophores at both stages. Our genetic analyses identified a 13-nucleotide deletion in the PAX7 transcription factor of MSS geckos, affecting its protein coding sequence. With our single-cell transcriptomics analysis of embryonic skin, we confirm that PAX7 is expressed in iridophores and xanthophores, suggesting that it plays a key role in the differentiation of both chromatophores. Our in situ hybridizations on whole-mount embryos document the dynamics of the skin pattern formation and how it is impacted in the PAX7 mutants. We hypothesize that the melanophores-iridophores interactions give rise to the banded pattern of the hatchlings and black spot formation is an intrinsic capacity of melanophores in the postembryonic skin.

Anatomy, development and regeneration of zebrafish elasmoid scales

Vertebrate skin appendages - particularly avian feathers and mammalian hairs, glands and teeth - are perennially useful systems for investigating fundamental mechanisms of development. The most common type of skin appendage in teleost fishes is the elasmoid scale, yet this structure has received much less attention than the skin appendages of tetrapods. Elasmoid scales are thin, overlapping plates of partially mineralized extracellular matrices, deposited in the skin in a hexagonal pattern by a specialized population of dermal cells in cooperation with the overlying epidermis. Recent years have seen rapid progress in our understanding of elasmoid scale development and regeneration, driven by the deployment of developmental genetics, live imaging and transcriptomics in larval and adult zebrafish. These findings are reviewed together with histological and ultrastructural approaches to understanding scale development and regeneration.

Periodic pattern formation during embryonic development

During embryonic development many organs and structures require the formation of series of repeating elements known as periodic patterns. Ranging from the digits of the limb to the feathers of the avian skin, the correct formation of these embryonic patterns is essential for the future form and function of these tissues. However, the mechanisms that produce these patterns are not fully understood due to the existence of several modes of pattern generation which often differ between organs and species. Here, we review the current state of the field and provide a perspective on future approaches to studying this fundamental process of embryonic development.

On the role of TFEC in reptilian coloration

Reptilian species, particularly snakes and lizards, are emerging models of animal coloration. Here, I focus on the role of the TFEC transcription factor in snake and lizard coloration based on a study on wild-type and piebald ball pythons. Genomic mapping previously identified a TFEC mutation linked to the piebald ball python phenotype. The association of TFEC with skin coloration was further supported by gene-editing experiments in the brown anole lizard. However, novel histological analyses presented here reveal discrepancies between the ball python and the anole TFEC mutants phenotype, cautioning against broad generalizations. Indeed, both wild-type and piebald ball pythons completely lack iridophores, whereas the TFEC anole lizard mutants lose their iridophores compared to the wild-type anole. Based on these findings, I discuss the potential role of the MiT/TFE family in skin pigmentation across vertebrate lineages and advocate the need for developmental analyses and additional gene-editing experiments to explore the reptilian coloration diversity.

Development shapes the evolutionary diversification of rodent stripe patterns

Vertebrate groups have evolved strikingly diverse color patterns. However, it remains unknown to what extent the diversification of such patterns has been shaped by the proximate, developmental mechanisms that regulate their formation. While these developmental mechanisms have long been inaccessible empirically, here we take advantage of recent insights into rodent pattern formation to investigate the role of development in shaping pattern diversification across rodents. Based on a broad survey of museum specimens, we first establish that various rodents have independently evolved diverse patterns consisting of longitudinal stripes, varying across species in number, color, and relative positioning. We then interrogate this diversity using a simple model that incorporates recent molecular and developmental insights into stripe formation in African striped mice. Our results suggest that, on the one hand, development has facilitated pattern diversification: The diversity of patterns seen across species can be generated by a single developmental process, and small changes in this process suffice to recapitulate observed evolutionary changes in pattern organization. On the other hand, development has constrained diversification: Constraints on stripe positioning limit the scope of evolvable patterns, and although pattern organization appears at first glance phylogenetically unconstrained, development turns out to impose a cryptic constraint. Altogether, this work reveals that pattern diversification in rodents can in part be explained by the underlying development and illustrates how pattern formation models can be leveraged to interpret pattern evolution.

Genetics and biology of coloration in reptiles: the curious case of the Lemon Frost geckos

Although there are more than 10,000 reptile species, and reptiles have historically contributed to our understanding of biology, genetics research into class Reptilia has lagged compared with other animals. Here, we summarize recent progress in genetics of coloration in reptiles, with a focus on the leopard gecko, Eublepharis macularius. We highlight genetic approaches that have been used to examine variation in color and pattern formation in this species as well as to provide insights into mechanisms underlying skin cancer. We propose that their long breeding history in captivity makes leopard geckos one of the most promising emerging reptilian models for genetic studies. More broadly, technological advances in genetics, genomics, and gene editing may herald a golden era for studies of reptile biology.

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.

An active traveling wave of Eda/NF-κB signaling controls the timing and hexagonal pattern of skin appendages in zebrafish

Periodic patterns drive the formation of a variety of tissues, including skin appendages such as feathers and scales. Skin appendages serve important and diverse functions across vertebrates, yet the mechanisms that regulate their patterning are not fully understood. Here, we have used live imaging to investigate dynamic signals regulating the ontogeny of zebrafish scales. Scales are bony skin appendages that develop sequentially along the anterior-posterior and dorsal-ventral axes to cover the fish in a hexagonal array. We have found that scale development requires cell-cell communication and is coordinated through an active wave mechanism. Using a live transcriptional reporter, we show that a wave of Eda/NF-κB activity precedes scale initiation and is required for scale formation. Experiments decoupling the propagation of the wave from dermal placode formation and osteoblast differentiation demonstrate that the Eda/NF-κB activity wavefront controls the timing of the sequential patterning of scales. Moreover, this decoupling resulted in defects in scale size and significant deviations in the hexagonal patterning of scales. Thus, our results demonstrate that a biochemical traveling wave coordinates scale initiation and proper hexagonal patterning across the fish body.