A chemotaxis model of feather primordia pattern formation during avian development

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.

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.

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.

Hydrocortisone treatment as a tool to study conjunctival placode induction

BACKGROUND

Conjunctival placodes are a series of placodes that develop into the conjunctival (scleral) papillae and ultimately induce a series of scleral ossicles in the eyes of many vertebrates. This study establishes a hydrocortisone injection procedure (incl. dosage) that consistently inhibits all conjunctival papillae in the embryonic chicken eye. The effects of this hydrocortisone treatment on apoptosis, vasculature, and placode-related gene expression were assessed.

RESULTS

Hydrocortisone treatment does not increase apoptotic cell death or have a major effect on the ciliary artery or vascular plexus in the eye. β-catenin and Eda expression levels were not significantly altered following hydrocortisone treatment, despite the absence of conjunctival papillae. Notably, Fgf20 expression was significantly reduced following hydrocortisone treatment, and the distribution of β-catenin was altered.

CONCLUSIONS

Our study showed that conjunctival papillae induction begins as early as HH27.5 (E5.5). Hydrocortisone treatment reduces Fgf20 expression independently of β-catenin and Eda and may instead affect other members of the Wnt/β-catenin or Eda/Edar pathways, or it may affect the ability of morphogens to diffuse through the extracellular matrix. This study contributes to a growing profile of gene expression data during placode development and enhances our understanding of how some vertebrate eyes develop these fascinating bones.

Structure, function and formation of the amniote skin pattern

From feather and hair dotted arrays to pigmented stripes and spots, the spatial distribution of skin appendages and colouration often forms visible ornaments crucial for fitness in the coat of birds and mammals. These geometrical motifs are extremely diverse in nature. Yet, phenotypic surveys evidenced common themes in variation: the orientation, appendage-specificity or pigmentation of a given region may be conserved across groups or species. Here, we review naturalist observations of natural variation in the anatomy and ecological function of the skin pattern in amniotes. We then describe several decades of genetics, mathematical modelling and experimental embryology work aiming at understanding the molecular and morphogenetic mechanisms responsible for pattern formation. We discuss how these studies provided evidence that the morphological trends and differences representative of the phenotypic landscape of skin patterns in wild amniote species is rooted in the mechanisms controlling the production of distinct compartments in the embryonic skin.

Mechanical induction in metazoan development and evolution: from earliest multi-cellular organisms to modern animal embryos

The development and origin of animal body forms have long been intensely explored, from the analysis of morphological traits during antiquity to Newtonian mechanical conceptions of morphogenesis. Advent of molecular biology then focused most interests on the biochemical patterning and genetic regulation of embryonic development. Today, a view is arising of development of multicellular living forms as a phenomenon emerging from non-hierarchical, reciprocal mechanical and mechanotransductive interactions between biochemical patterning and biomechanical morphogenesis. Here we discuss the nature of these processes and put forward findings on how early biochemical and biomechanical patterning of metazoans may have emerged from a primitive behavioural mechanotransducive feeding response to marine environment which might have initiated the development of first animal multicellular organisms.

Gap junctions in Turing-type periodic feather pattern formation

Periodic patterning requires coordinated cell-cell interactions at the tissue level. Turing showed, using mathematical modeling, how spatial patterns could arise from the reactions of a diffusive activator-inhibitor pair in an initially homogeneous 2D field. Most activators and inhibitors studied in biological systems are proteins, and the roles of cell-cell interaction, ions, bioelectricity, etc. are only now being identified. Gap junctions (GJs) mediate direct exchanges of ions or small molecules between cells, enabling rapid long-distance communications in a cell collective. They are therefore good candidates for propagating nonprotein-based patterning signals that may act according to the Turing principles. Here, we explore the possible roles of GJs in Turing-type patterning using feather pattern formation as a model. We found 7 of the 12 investigated GJ isoforms are highly dynamically expressed in the developing chicken skin. In ovo functional perturbations of the GJ isoform, connexin 30, by siRNA and the dominant-negative mutant applied before placode development led to disrupted primary feather bud formation. Interestingly, inhibition of gap junctional intercellular communication (GJIC) in the ex vivo skin explant culture allowed the sequential emergence of new feather buds at specific spatial locations relative to the existing primary buds. The results suggest that GJIC may facilitate the propagation of long-distance inhibitory signals. Thus, inhibition of GJs may stimulate Turing-type periodic feather pattern formation during chick skin development, and the removal of GJ activity would enable the emergence of new feather buds if the local environment were competent and the threshold to form buds was reached. We further propose Turing-based computational simulations that can predict the sequential appearance of these ectopic buds. Our models demonstrate how a Turing activator-inhibitor system can continue to generate patterns in the competent morphogenetic field when the level of intercellular communication at the tissue scale is modulated.

Turing Pattern Formation in Reaction-Cross-Diffusion Systems with a Bilayer Geometry

Conditions for self-organisation via Turing's mechanism in biological systems represented by reaction-diffusion or reaction-cross-diffusion models have been extensively studied. Nonetheless, the impact of tissue stratification in such systems is under-explored, despite its ubiquity in the context of a thin epithelium overlying connective tissue, for instance the epidermis and underlying dermal mesenchyme of embryonic skin. In particular, each layer can be subject to extensively different biochemical reactions and transport processes, with chemotaxis - a special case of cross-diffusion - often present in the mesenchyme, contrasting the solely molecular transport typically found in the epidermal layer. We study Turing patterning conditions for a class of reaction-cross-diffusion systems in bilayered regions, with a thin upper layer and coupled by a linear transport law. In particular, the role of differential transport through the interface is explored together with the presence of asymmetry between the homogeneous equilibria of the two layers. A linear stability analysis is carried out around a spatially homogeneous equilibrium state in the asymptotic limit of weak and strong coupling strengths, where quantitative approximations of the bifurcation curve can be computed. Our theoretical findings, for an arbitrary number of reacting species, reveal quantitative Turing conditions, highlighting when the coupling mechanism between the layered regions can either trigger patterning or stabilize a spatially homogeneous equilibrium regardless of the independent patterning state of each layer. We support our theoretical results through direct numerical simulations, and provide an open source code to explore such systems further.

Novel Aspects in Pattern Formation Arise from Coupling Turing Reaction-Diffusion and Chemotaxis

Recent experimental studies on primary hair follicle formation and feather bud morphogenesis indicate a coupling between Turing-type diffusion driven instability and chemotactic patterning. Inspired by these findings we develop and analyse a mathematical model that couples chemotaxis to a reaction-diffusion system exhibiting diffusion-driven (Turing) instability. While both systems, reaction-diffusion systems and chemotaxis, can independently generate spatial patterns, we were interested in how the coupling impacts the stability of the system, parameter region for patterning, pattern geometry, as well as the dynamics of pattern formation. We conduct a classical linear stability analysis for different model structures, and confirm our results by numerical analysis of the system. Our results show that the coupling generally increases the robustness of the patterning process by enlarging the pattern region in the parameter space. Concerning time scale and pattern regularity, we find that an increase in the chemosensitivity can speed up the patterning process for parameters inside and outside of the Turing space, but generally reduces spatial regularity of the pattern. Interestingly, our analysis indicates that pattern formation can also occur when neither the Turing nor the chemotaxis system can independently generate pattern. On the other hand, for some parameter settings, the coupling of the two processes can extinguish the pattern formation, rather than reinforce it. These theoretical findings can be used to corroborate the biological findings on morphogenesis and guide future experimental studies. From a mathematical point of view, this work sheds a light on coupling classical pattern formation systems from the parameter space perspective.

Morphogens enable interacting supracellular phases that generate organ architecture

During vertebrate organogenesis, increases in morphological complexity are tightly coupled to morphogen expression. In this work, we studied how morphogens influence self-organizing processes at the collective or "supra"-cellular scale in avian skin. We made physical measurements across length scales, which revealed morphogen-enabled material property differences that were amplified at supracellular scales in comparison to cellular scales. At the supracellular scale, we found that fibroblast growth factor (FGF) promoted "solidification" of tissues, whereas bone morphogenetic protein (BMP) promoted fluidity and enhanced mechanical activity. Together, these effects created basement membrane-less compartments within mesenchymal tissue that were mechanically primed to drive avian skin tissue budding. Understanding this multiscale process requires the ability to distinguish between proximal effects of morphogens that occur at the cellular scale and their functional effects, which emerge at the supracellular scale.

Somitic positional information guides self-organized patterning of snake scales

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.

Gap junctions in Turing-type periodic feather pattern formation

Periodic patterning requires coordinated cell-cell interactions at the tissue level. Turing showed, using mathematical modeling, how spatial patterns could arise from the reactions of a diffusive activator-inhibitor pair in an initially homogenous two-dimensional field. Most activators and inhibitors studied in biological systems are proteins, and the roles of cell-cell interaction, ions, bioelectricity, etc. are only now being identified. Gap junctions (GJs) mediate direct exchanges of ions or small molecules between cells, enabling rapid long-distance communications in a cell collective. They are therefore good candidates for propagating non-protein-based patterning signals that may act according to the Turing principles. Here, we explore the possible roles of GJs in Turing-type patterning using feather pattern formation as a model. We found seven of the twelve investigated GJ isoforms are highly dynamically expressed in the developing chicken skin. In ovo functional perturbations of the GJ isoform, connexin 30, by siRNA and the dominant-negative mutant applied before placode development led to disrupted primary feather bud formation, including patches of smooth skin and buds of irregular sizes. Later, after the primary feather arrays were laid out, inhibition of gap junctional intercellular communication in the ex vivo skin explant culture allowed the emergence of new feather buds in temporal waves at specific spatial locations relative to the existing primary buds. The results suggest that gap junctional communication may facilitate the propagation of long-distance inhibitory signals. Thus, the removal of GJ activity would enable the emergence of new feather buds if the local environment is competent and the threshold to form buds is reached. We propose Turing-based computational simulations that can predict the appearance of these ectopic bud waves. Our models demonstrate how a Turing activator-inhibitor system can continue to generate patterns in the competent morphogenetic field when the level of intercellular communication at the tissue scale is modulated.

Extensive embryonic patterning without cellular differentiation primes the plant epidermis for efficient post-embryonic stomatal activities

Plant leaves feature epidermal stomata that are organized in stereotyped patterns. How does the pattern originate? We provide transcriptomic, imaging, and genetic evidence that Arabidopsis embryos engage known stomatal fate and patterning factors to create regularly spaced stomatal precursor cells. Analysis of embryos from 36 plant species indicates that this trait is widespread among angiosperms. Embryonic stomatal patterning in Arabidopsis is established in three stages: first, broad SPEECHLESS (SPCH) expression; second, coalescence of SPCH and its targets into discrete domains; and third, one round of asymmetric division to create stomatal precursors. Lineage progression is then halted until after germination. We show that the embryonic stomatal pattern enables fast stomatal differentiation and photosynthetic activity upon germination, but it also guides the formation of additional stomata as the leaf expands. In addition, key stomatal regulators are prevented from driving the fate transitions they can induce after germination, identifying stage-specific layers of regulation that control lineage progression during embryogenesis.

Cell shape anisotropy contributes to self-organized feather pattern fidelity in birds

Developing tissues can self-organize into a variety of patterned structures through the stabilization of stochastic fluctuations in their molecular and cellular properties. While molecular factors and cell dynamics contributing to self-organization have been identified in vivo, events channeling self-organized systems such that they achieve stable pattern outcomes remain unknown. Here, we described natural variation in the fidelity of self-organized arrays formed by feather follicle precursors in bird embryos. By surveying skin cells prior to and during tissue self-organization and performing species-specific ex vivo drug treatments and mechanical stress tests, we demonstrated that pattern fidelity depends on the initial amplitude of cell anisotropy in regions of the developing dermis competent to produce a pattern. Using live imaging, we showed that cell shape anisotropy is associated with a limited increase in cell motility for sharp and precisely located primordia formation, and thus, proper pattern geometry. These results evidence a mechanism through which initial tissue properties ensure stability in self-organization and thus, reproducible pattern production.

A study of natural variation in feather pattern geometry and a combination of pharmacological and mechanical manipulations ex vivo provides evidence for a mechanism by which initial tissue properties ensure stability in self-organization and species-specific pattern production.

Emergence of a geometric pattern of cell fates from tissue-scale mechanics in the Drosophila eye

Pattern formation of biological structures involves the arrangement of different types of cells in an ordered spatial configuration. In this study, we investigate the mechanism of patterning the Drosophila eye epithelium into a precise triangular grid of photoreceptor clusters called ommatidia. Previous studies had led to a long-standing biochemical model whereby a reaction-diffusion process is templated by recently formed ommatidia to propagate a molecular prepattern across the eye. Here, we find that the templating mechanism is instead, mechanochemical in origin; newly born columns of differentiating ommatidia serve as a template to spatially pattern flows that move epithelial cells into position to form each new column of ommatidia. Cell flow is generated by a source and sink, corresponding to narrow zones of cell dilation and contraction respectively, that straddle the growing wavefront of ommatidia. The newly formed lattice grid of ommatidia cells are immobile, deflecting, and focusing the flow of other cells. Thus, the self-organization of a regular pattern of cell fates in an epithelium is mechanically driven.

Systems for intricate patterning of the vertebrate anatomy

Periodic patterns form intricate arrays in the vertebrate anatomy, notably the hair and feather follicles of the skin, but also internally the villi of the gut and the many branches of the lung, kidney, mammary and salivary glands. These tissues are composite structures, being composed of adjoined epithelium and mesenchyme, and the patterns that arise within them require interaction between these two tissue layers. In embryonic development, cells change both their distribution and state in a periodic manner, defining the size and relative positions of these specialized structures. Their placement is determined by simple spacing mechanisms, with substantial evidence pointing to a variety of local enhancement/lateral inhibition systems underlying the breaking of symmetry. The nature of the cellular processes involved, however, has been less clear. While much attention has focused on intercellular soluble signals, such as protein growth factors, experimental evidence has grown for contributions of cell movement or mechanical forces to symmetry breaking. In the mesenchyme, unlike the epithelium, cells may move freely and can self-organize into aggregates by chemotaxis, or through generation and response to mechanical strain on their surrounding matrix. Different modes of self-organization may coexist, either coordinated into a single system or with hierarchical relationships. Consideration of a broad range of distinct biological processes is required to advance understanding of biological pattern formation. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Sweat Gland Organoids Originating from Reprogrammed Epidermal Keratinocytes Functionally Recapitulated Damaged Skin

Restoration of sweat glands (SwGs) represents a great issue in patients with extensive skin defects. Recent methods combining organoid technology with cell fate reprogramming hold promise for developing new regenerative methods for SwG regeneration. Here, a practical strategy for engineering functional human SwGs in vitro and in vivo is provided. First, by forced expression of the ectodysplasin-A in human epidermal keratinocytes (HEKs) combined with specific SwG culture medium, HEKs are efficiently converted into SwG cells (iSwGCs). The iSwGCs show typical morphology, gene expression pattern, and functions resembling human primary SwG cells. Second, by culturing the iSwGCs in a special 3D culturing system, SwG organoids (iSwGOs) that exhibit structural and biological features characteristic of native SwGs are obtained. Finally, these iSwGOs are successfully transplanted into a mouse skin damage model and they develop into fully functioning SwGs in vivo. Regeneration of functional SwG organoids from reprogrammed HEKs highlights the great translational potential for personalized SwG regeneration in patients with large skin defects.

Trends and variation in vertebrate patterns as outcomes of self-organization

In extant vertebrates, natural motifs such as coat markings, spongy bone structures, neuronal arborescence or collective swarms correspond to diverse pattern types, from fractals to periodic stripes or tessellations, and so on. In this subphylum, evolution produced an apparent paradox: a given pattern may vary tremendously in its extent, periodicity or regularity, but follows general geometrical trends and is produced with meticulous precision. In this review, we discuss the role of self-organization, a patterning strategy in which spontaneous order arises through local interactions without gradient formation, in shaping both natural pattern differences and common themes. Mathematical models evidenced a wide high adaptability of self-organizing dynamics, long-advocating for their contribution to natural pattern diversity. Recent empirical and theoretical approaches taking into account network topologies and natural variation also replaced outcomes of self-organization in more constrained biological contexts, shedding light on mechanisms ensuring pattern fidelity.

A "Numerical Evo-Devo" Synthesis for the Identification of Pattern-Forming Factors

Animals display extensive diversity in motifs adorning their coat, yet these patterns have reproducible orientation and periodicity within species or groups. Morphological variation has been traditionally used to dissect the genetic basis of evolutionary change, while pattern conservation and stability in both mathematical and organismal models has served to identify core developmental events. Two patterning theories, namely instruction and self-organisation, emerged from this work. Combined, they provide an appealing explanation for how natural patterns form and evolve, but in vivo factors underlying these mechanisms remain elusive. By bridging developmental biology and mathematics, novel frameworks recently allowed breakthroughs in our understanding of pattern establishment, unveiling how patterning strategies combine in space and time, or the importance of tissue morphogenesis in generating positional information. Adding results from surveys of natural variation to these empirical-modelling dialogues improves model inference, analysis, and in vivo testing. In this evo-devo-numerical synthesis, mathematical models have to reproduce not only given stable patterns but also the dynamics of their emergence, and the extent of inter-species variation in these dynamics through minimal parameter change. This integrative approach can help in disentangling molecular, cellular and mechanical interaction during pattern establishment.

From a discrete model of chemotaxis with volume-filling to a generalized Patlak-Keller-Segel model

We present a discrete model of chemotaxis whereby cells responding to a chemoattractant are seen as individual agents whose movement is described through a set of rules that result in a biased random walk. In order to take into account possible alterations in cellular motility observed at high cell densities (i.e. volume-filling), we let the probabilities of cell movement be modulated by a decaying function of the cell density. We formally show that a general form of the celebrated Patlak-Keller-Segel (PKS) model of chemotaxis can be formally derived as the appropriate continuum limit of this discrete model. The family of steady-state solutions of such a generalized PKS model are characterized and the conditions for the emergence of spatial patterns are studied via linear stability analysis. Moreover, we carry out a systematic quantitative comparison between numerical simulations of the discrete model and numerical solutions of the corresponding PKS model, both in one and in two spatial dimensions. The results obtained indicate that there is excellent quantitative agreement between the spatial patterns produced by the two models. Finally, we numerically show that the outcomes of the two models faithfully replicate those of the classical PKS model in a suitable asymptotic regime.

A dot-stripe Turing model of joint patterning in the tetrapod limb

Iterative joints are a hallmark of the tetrapod limb, and their positioning is a key step during limb development. Although the molecular regulation of joint formation is well studied, it remains unclear what controls the location, number and orientation (i.e. the pattern) of joints within each digit. Here, we propose the dot-stripe mechanism for joint patterning, comprising two coupled Turing systems inspired by published gene expression patterns. Our model can explain normal joint morphology in wild-type limbs, hyperphalangy in cetacean flippers, mutant phenotypes with misoriented joints and suggests a reinterpretation of the polydactylous Ichthyosaur fins as a polygonal joint lattice. By formulating a generic dot-stripe model, describing joint patterns rather than molecular joint markers, we demonstrate that the insights from the model should apply regardless of the biological specifics of the underlying mechanism, thus providing a unifying framework to interrogate joint patterning in the tetrapod limb.

Self-Organization in Pattern Formation

Self-organization is pervasive in development, from symmetry breaking in the early embryo to tissue patterning and morphogenesis. For a few model systems, the underlying molecular and cellular processes are now sufficiently characterized that mathematical models can be confronted with experiments, to explore the dynamics of pattern formation. Here, we review selected systems, ranging from cyanobacteria to mammals, where different forms of cell-cell communication, acting alone or together with positional cues, drive the patterning of cell fates, highlighting the insights that even very simple models can provide as well as the challenges on the path to a predictive understanding of development.

Exploration of key regulators driving primary feather follicle induction in goose skin

The primary feather follicles are universal skin appendages widely distributed in the skin of feathered birds. The morphogenesis and development of the primary feather follicles in goose skin remain largely unknown. Here, the induction of primary feather follicles in goose embryonic skin (pre-induction vs induction) was investigated by de novo transcriptome analyses to reveal 409 differentially expressed genes (DEGs). The DEGs were characterized to potentially regulate the de novo formation of feather follicle primordia consisting of placode (4 genes) and dermal condensate (12 genes), and the thickening of epidermis (5 genes) and dermal fibroblasts (17 genes), respectively. Further analyses enriched DEGs into GO terms represented as cell adhesion and KEGG pathways including Wnt and Hedgehog signaling pathways that are highly correlated with cell communication and molecular regulation. Six selected Wnt pathway genes were detected by qPCR with up-regulation in goose skin during the induction of primary feather follicles. The localization of WNT16, SFRP1 and FRZB by in situ hybridization showed weak expression in the primary feather primordia, whereas FZD1, LEF1 and DKK1 were expressed initially in the inter-follicular skin and feather follicle primordia, then mainly restricted in the feather primordia. The spatial-temporal expression patterns indicate that Wnt pathway genes DKK1, FZD1 and LEF1 are the important regulators functioned in the induction of primary feather follicle in goose skin. The dynamic molecular changes and specific gene expression patterns revealed in this report provide the general knowledge of primary feather follicle and skin development in waterfowl, and contribute to further understand the diversity of hair and feather development beyond the mouse and chicken models.

Symmetry breaking in the embryonic skin triggers directional and sequential plumage patterning

The development of an organism involves the formation of patterns from initially homogeneous surfaces in a reproducible manner. Simulations of various theoretical models recapitulate final states of natural patterns, yet drawing testable hypotheses from those often remains difficult. Consequently, little is known about pattern-forming events. Here, we surveyed plumage patterns and their emergence in Galliformes, ratites, passerines, and penguins, together representing the three major taxa of the avian phylogeny, and built a unified model that not only reproduces final patterns but also intrinsically generates shared and varying directionality, sequence, and duration of patterning. We used in vivo and ex vivo experiments to test its parameter-based predictions. We showed that directional and sequential pattern progression depends on a species-specific prepattern: an initial break in surface symmetry launches a travelling front of sharply defined, oriented domains with self-organising capacity. This front propagates through the timely transfer of increased cell density mediated by cell proliferation, which controls overall patterning duration. These results show that universal mechanisms combining prepatterning and self-organisation govern the timely emergence of the plumage pattern in birds.

A survey of plumage patterns and their emergence in Galliformes, ratites, passerines, and penguins shows that their formation depends on a species-specific prepattern in the embryo and demonstrates that universal mechanisms govern the timely emergence of natural patterns in birds.

Non-uniform curvature and anisotropic deformation control wrinkling patterns on tori

We investigate wrinkling patterns in a tri-layer torus consisting of an expanding thin outer layer, an intermediate soft layer and an inner core with a tunable shear modulus, inspired by pattern formation in developmental biology, such as follicle pattern formation during the development of chicken embryos. We show from large-scale finite element simulations that hexagonal wrinkling patterns form for stiff cores whereas stripe wrinkling patterns develop for soft cores. Hexagons and stripes co-exist to form hybrid patterns for cores with intermediate stiffness. The governing mechanism for the pattern transition is that the stiffness of the inner core controls the degree to which the major radius of the torus expands - this has a greater effect on deformation in the long direction as compared to the short direction of the torus. This anisotropic deformation alters stress states in the outer layer which change from biaxial (preferred hexagons) to uniaxial (preferred stripes) compression as the core stiffness is reduced. As the outer layer continues to expand, stripe and hexagon patterns will evolve into zigzags and segmented labyrinths, respectively. Stripe wrinkles are observed to initiate at the inner surface of the torus while hexagon wrinkles start from the outer surface as a result of curvature-dependent stresses in the torus. We further discuss the effects of elasticities and geometries of the torus on the wrinkling patterns.

Inherency of Form and Function in Animal Development and Evolution

I discuss recent work on the origins of morphology and cell-type diversification in Metazoa – collectively the animals – and propose a scenario for how these two properties became integrated, with the help of a third set of processes, cellular pattern formation, into the developmental programs seen in present-day metazoans. Inherent propensities to generate familiar forms and cell types, in essence a parts kit for the animals, are exhibited by present-day organisms and were likely more prominent in primitive ones. The structural motifs of animal bodies and organs, e.g., multilayered, hollow, elongated and segmented tissues, internal and external appendages, branched tubes, and modular endoskeletons, can be accounted for by the properties of mesoscale masses of metazoan cells. These material properties, in turn, resulted from the recruitment of “generic” physical forces and mechanisms – adhesion, contraction, polarity, chemical oscillation, diffusion – by toolkit molecules that were partly conserved from unicellular holozoan antecedents and partly novel, distributed in the different metazoan phyla in a fashion correlated with morphological complexity. The specialized functions of the terminally differentiated cell types in animals, e.g., contraction, excitability, barrier function, detoxification, excretion, were already present in ancestral unicellular organisms. These functions were implemented in metazoan differentiation in some cases using the same transcription factors as in single-celled ancestors, although controlled by regulatory mechanisms that were hybrids between earlier-evolved processes and regulatory innovations, such as enhancers. Cellular pattern formation, mediated by released morphogens interacting with biochemically responsive and excitable tissues, drew on inherent self-organizing processes in proto-metazoans to transform clusters of holozoan cells into animal embryos and organs.

Modeling Edar expression reveals the hidden dynamics of tooth signaling center patterning

When patterns are set during embryogenesis, it is expected that they are straightly established rather than subsequently modified. The patterning of the three mouse molars is, however, far from straight, likely as a result of mouse evolutionary history. The first-formed tooth signaling centers, called MS and R2, disappear before driving tooth formation and are thought to be vestiges of the premolars found in mouse ancestors. Moreover, the mature signaling center of the first molar (M1) is formed from the fusion of two signaling centers (R2 and early M1). Here, we report that broad activation of Edar expression precedes its spatial restriction to tooth signaling centers. This reveals a hidden two-step patterning process for tooth signaling centers, which was modeled with a single activator-inhibitor pair subject to reaction-diffusion (RD). The study of Edar expression also unveiled successive phases of signaling center formation, erasing, recovering, and fusion. Our model, in which R2 signaling center is not intrinsically defective but erased by the broad activation preceding M1 signaling center formation, predicted the surprising rescue of R2 in Edar mutant mice, where activation is reduced. The importance of this R2-M1 interaction was confirmed by ex vivo cultures showing that R2 is capable of forming a tooth. Finally, by introducing chemotaxis as a secondary process to RD, we recapitulated in silico different conditions in which R2 and M1 centers fuse or not. In conclusion, pattern formation in the mouse molar field relies on basic mechanisms whose dynamics produce embryonic patterns that are plastic objects rather than fixed end points.

Feather arrays are patterned by interacting signalling and cell density waves

Feathers are arranged in a precise pattern in avian skin. They first arise during development in a row along the dorsal midline, with rows of new feather buds added sequentially in a spreading wave. We show that the patterning of feathers relies on coupled fibroblast growth factor (FGF) and bone morphogenetic protein (BMP) signalling together with mesenchymal cell movement, acting in a coordinated reaction-diffusion-taxis system. This periodic patterning system is partly mechanochemical, with mechanical-chemical integration occurring through a positive feedback loop centred on FGF20, which induces cell aggregation, mechanically compressing the epidermis to rapidly intensify FGF20 expression. The travelling wave of feather formation is imposed by expanding expression of Ectodysplasin A (EDA), which initiates the expression of FGF20. The EDA wave spreads across a mesenchymal cell density gradient, triggering pattern formation by lowering the threshold of mesenchymal cells required to begin to form a feather bud. These waves, and the precise arrangement of feather primordia, are lost in the flightless emu and ostrich, though via different developmental routes. The ostrich retains the tract arrangement characteristic of birds in general but lays down feather primordia without a wave, akin to the process of hair follicle formation in mammalian embryos. The embryonic emu skin lacks sufficient cells to enact feather formation, causing failure of tract formation, and instead the entire skin gains feather primordia through a later process. This work shows that a reaction-diffusion-taxis system, integrated with mechanical processes, generates the feather array. In flighted birds, the key role of the EDA/Ectodysplasin A receptor (EDAR) pathway in vertebrate skin patterning has been recast to activate this process in a quasi-1-dimensional manner, imposing highly ordered pattern formation.