Spots and stripes: pleomorphic patterning of stem cells via p-ERK-dependent cell chemotaxis shown by feather morphogenesis and mathematical simulation

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 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.

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.

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.

A shark-inspired general model of tooth morphogenesis unveils developmental asymmetries in phenotype transitions

Developmental complexity stemming from the dynamic interplay between genetic and biomechanic factors canalizes the ways genotypes and phenotypes can change in evolution. As a paradigmatic system, we explore how changes in developmental factors generate typical tooth shape transitions. Since tooth development has mainly been researched in mammals, we contribute to a more general understanding by studying the development of tooth diversity in sharks. To this end, we build a general, but realistic, mathematical model of odontogenesis. We show that it reproduces key shark-specific features of tooth development as well as real tooth shape variation in small-spotted catsharks Scyliorhinus canicula. We validate our model by comparison with experiments in vivo. Strikingly, we observe that developmental transitions between tooth shapes tend to be highly degenerate, even for complex phenotypes. We also discover that the sets of developmental parameters involved in tooth shape transitions tend to depend asymmetrically on the direction of that transition. Together, our findings provide a valuable base for furthering our understanding of how developmental changes can lead to both adaptive phenotypic change and trait convergence in complex, phenotypically highly diverse, structures.

The Impact of Phenotypic Heterogeneity on Chemotactic Self-Organisation

The capacity to aggregate through chemosensitive movement forms a paradigm of self-organisation, with examples spanning cellular and animal systems. A basic mechanism assumes a phenotypically homogeneous population that secretes its own attractant, with the well known system introduced more than five decades ago by Keller and Segel proving resolutely popular in modelling studies. The typical assumption of population phenotypic homogeneity, however, often lies at odds with the heterogeneity of natural systems, where populations may comprise distinct phenotypes that vary according to their chemotactic ability, attractant secretion, etc. To initiate an understanding into how this diversity can impact on autoaggregation, we propose a simple extension to the classical Keller and Segel model, in which the population is divided into two distinct phenotypes: those performing chemotaxis and those producing attractant. Using a combination of linear stability analysis and numerical simulations, we demonstrate that switching between these phenotypic states alters the capacity of a population to self-aggregate. Further, we show that switching based on the local environment (population density or chemoattractant level) leads to diverse patterning and provides a route through which a population can effectively curb the size and density of an aggregate. We discuss the results in the context of real world examples of chemotactic aggregation, as well as theoretical aspects of the model such as global existence and blow-up of solutions.

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.

Parameterized Computational Framework for the Description and Design of Genetic Circuits of Morphogenesis Based on Contact-Dependent Signaling and Changes in Cell-Cell Adhesion

Synthetic development is a nascent field of research that uses the tools of synthetic biology to design genetic programs directing cellular patterning and morphogenesis in higher eukaryotic cells, such as mammalian cells. One specific example of such synthetic genetic programs was based on cell-cell contact-dependent signaling using synthetic Notch pathways and was shown to drive the formation of multilayered spheroids by modulating cell-cell adhesion via differential expression of cadherin family proteins in a mouse fibroblast cell line (L929). The design method for these genetic programs relied on trial and error, which limited the number of possible circuits and parameter ranges that could be explored. Here, we build a parameterized computational framework that, given a cell-cell communication network driving changes in cell adhesion and initial conditions as inputs, predicts developmental trajectories. We first built a general computational framework where contact-dependent cell-cell signaling networks and changes in cell-cell adhesion could be designed in a modular fashion. We then used a set of available in vitro results (that we call the "training set" in analogy to similar pipelines in the machine learning field) to parameterize the computational model with values for adhesion and signaling. We then show that this parameterized model can qualitatively predict experimental results from a "testing set" of available in vitro data that varied the genetic network in terms of adhesion combinations, initial number of cells, and even changes to the network architecture. Finally, this parameterized model is used to recommend novel network implementation for the formation of a four-layered structure that has not been reported previously. The framework that we develop here could function as a testing ground to identify the reachable space of morphologies that can be obtained by controlling contact-dependent cell-cell communications and adhesion with these molecular tools and in this cellular system. Additionally, we discuss how the model could be expanded to include other forms of communication or effectors for the computational design of the next generation of synthetic developmental trajectories.

Coordination of local and long range signaling modulates developmental patterning

The development of multicellular organisms relies on correct patterns of cell fates to produce functional tissues in the mature organism. A commonly observed developmental pattern consists of alternating cell fates, where neighboring cells take on distinct cell fates characterized by contrasting gene and protein expression levels, and this cell fate pattern repeats over two or more cells. The patterns produced by these fate decisions are regulated by a small number of highly conserved signaling networks, some of which are mediated by long range diffusible signals and others mediated by local contact-dependent signals. However, it is not completely understood how local and long range signals associated with these networks interact to produce fate patterns that are both robust and flexible. Here we analyze mathematical models to investigate the patterning of cell fates in an array of cells, focusing on a two cell repeating pattern. Bifurcation analysis of a multicellular ODE model, where we consider the cells as discrete compartments, suggests that cells must balance sensitivity to external signals with robustness to perturbations. To focus on the patterning dynamics close to the bifurcation point, we derive a continuum PDE model that integrates local and long range signaling. For those cells with dynamics close to the bifurcation point, sensitivity to long range signals determines how far a pattern extends in space, while the number of local signaling connections determines the type of pattern produced. This investigation provides a general framework for understanding developmental patterning, and how both long range and local signals play a role in generating features observed across biology, such as species differences in nematode vulval development and insect bristle patterning, as well as medically relevant processes such as control of stem cell fate in the intestinal crypt.

Self-organizing hair peg-like structures from dissociated skin progenitor cells: New insights for human hair follicle organoid engineering and Turing patterning in an asymmetric morphogenetic field

Human skin progenitor cells will form new hair follicles, although at a low efficiency, when injected into nude mouse skin. To better study and improve upon this regenerative process, we developed an in vitro system to analyse the morphogenetic cell behaviour in detail and modulate physical-chemical parameters to more effectively generate hair primordia. In this three-dimensional culture, dissociated human neonatal foreskin keratinocytes self-assembled into a planar epidermal layer while fetal scalp dermal cells coalesced into stripes, then large clusters, and finally small clusters resembling dermal condensations. At sites of dermal clustering, subjacent epidermal cells protruded to form hair peg-like structures, molecularly resembling hair pegs within the sequence of follicular development. The hair peg-like structures emerged in a coordinated, formative wave, moving from periphery to centre, suggesting that the droplet culture constitutes a microcosm with an asymmetric morphogenetic field. In vivo, hair follicle populations also form in a progressive wave, implying the summation of local periodic patterning events with an asymmetric global influence. To further understand this global patterning process, we developed a mathematical simulation using Turing activator-inhibitor principles in an asymmetric morphogenetic field. Together, our culture system provides a suitable platform to (a) analyse the self-assembly behaviour of hair progenitor cells into periodically arranged hair primordia and (b) identify parameters that impact the formation of hair primordia in an asymmetric morphogenetic field. This understanding will enhance our future ability to successfully engineer human hair follicle organoids.

Quantitative gene expression dynamics of key placode signalling factors in the embryonic chicken scleral ossicle system

Development of the scleral ossicles, a ring of bony elements within the sclera, is directed by a series of papillae that arise from placodes in the conjunctival epithelium over a 1.5-day induction period in the chicken embryo. The regular spacing of the papillae around the corneal-scleral limbus suggests that their induction may be regulated by a reaction-diffusion mechanism, similar to other epithelial appendages. Some key placode signalling molecules, including β-catenin, are known to be expressed throughout the induction period. However, others have been studied only at certain stages or have not been successfully detected. Here we use qPCR to study the gene expression patterns of the wingless integration (WNT)/β-catenin, bone morphogenetic protein (BMP), ectodysplasin (EDA), fibroblast growth factor (FGF) and hedgehog (HH) signalling families in discrete regions of the eye throughout the complete conjunctival placode and papillae induction period. This comprehensive analysis revealed a variable level of gene expression within specific eye regions, with some genes exhibiting high, moderate or low changes. Most genes exhibited an initial increase in gene expression, followed by a decrease or plateau as development proceeded, suggesting that some genes are important for a brief initial period whilst the sustained elevated expression level of other genes is needed for developmental progression. The timing or magnitude of these changes, and/or the overall gene expression trend differed in the temporal, nasal and/or dorsal eye regions for some, but not all genes, demonstrating that gene expression may vary across different eye regions. Temporal and nasal EDA receptor (EDAR) had the greatest number of strong correlations (r > 0.700) with other genes and β-catenin had the greatest number of moderate correlations (r = 0.400-0.700), while EDA had the greatest range in correlation strengths. Among the strongly correlated genes, two distinct signalling modules were identified, connected by some intermediate genes. The dynamic gene expression patterns of the five signalling pathways studied here from conjunctival placode formation through to papillae development is consistent with other epithelial appendages and confirms the presence of a conserved induction and patterning signalling network. Two unique gene expression patterns and corresponding gene interaction modules suggest functionally distinct roles throughout placode development. Furthermore, spatial differences in gene expression patterns among the temporal, nasal and dorsal regions of the eye may indicate that the expression of certain genes is influenced by mechanical forces exerted throughout development. Therefore, this study identifies key placode signalling factors and their interactions, as well as some potential region-specific features of gene expression in the scleral ossicle system and provides a basis for further exploration of the spatial expression of these genes and the patterning mechanism(s) active throughout conjunctival placode and papillae formation.

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.

Synthetic development: building mammalian multicellular structures with artificial genetic programs

Synthetic biology efforts began in simple single-cell systems, which were relatively easy to manipulate genetically (Cameron et al., 2014). The field grew exponentially in the last two decades, and one of the latest frontiers are synthetic developmental programs for multicellular mammalian systems (Black et al., 2017; Wieland and Fussenegger, 2012) to genetically control features such as patterning or morphogenesis. These programs rely on engineered cell-cell communications, multicellular gene regulatory networks and effector genes. Here, we contextualize the first of these synthetic developmental programs, examine molecular and computational tools that can be used to generate next generation versions, and present the general logic that underpins these approaches. These advances are exciting as they represent a novel way to address both control and understanding in the field of developmental biology and tissue development (Elowitz and Lim, 2010; Velazquez et al., 2018; White et al., 2018; Morsut, 2017). This field is just at the beginning, and it promises to be of major interest in the upcoming years of biomedical research.

STOCHASTIC DYNAMICS OF CELL LINEAGE IN TISSUE HOMEOSTASIS

During epithelium tissue maintenance, lineages of cells differentiate and proliferate in a coordinated way to provide the desirable size and spatial organization of different types of cells. While mathematical models through deterministic description have been used to dissect role of feedback regulations on tissue layer size and stratification, how the stochastic effects influence tissue maintenance remains largely unknown. Here we present a stochastic continuum model for cell lineages to investigate how both layer thickness and layer stratification are affected by noise. We find that the cell-intrinsic noise often causes reduction and oscillation of layer size whereas the cell-extrinsic noise increases the thickness, and sometimes, leads to uncontrollable growth of the tissue layer. The layer stratification usually deteriorates as the noise level increases in the cell lineage systems. Interestingly, the morphogen noise, which mixes both cell-intrinsic noise and cell-extrinsic noise, can lead to larger size of layer with little impact on the layer stratification. By investigating different combinations of the three types of noise, we find the layer thickness variability is reduced when cell-extrinsic noise level is high or morphogen noise level is low. Interestingly, there exists a tradeoff between low thickness variability and strong layer stratification due to competition among the three types of noise, suggesting robust layer homeostasis requires balanced levels of different types of noise in the cell lineage systems.

Turing patterning with and without a global wave

Periodic patterning represents a fundamental process in tissue morphogenesis. In chicken dorsal skin, feather formation starts from the midline; then the morphogenetic wave propagates bilaterally, leaving a regular hexagonal array of feather germs. Yet, in vitro reconstitution showed feather germs appear simultaneously, leading to the hypothesis that the feather-forming wave results from the coupling of local Turing patterning processes with an unidentified global event. In this issue, Ho and colleagues showed such a global event in chicken feathers involves a spreading Ectodysplasin A (EDA) wave and Fibroblast Growth Factor 20 (FGF20)-cell aggregate-based mechanochemical coupling. In flightless birds, feather germs form periodically but without precise hexagonal patterning due to the lack of global wave.

Chicken feather buds form exquisite hexagonal patterns progressively from the midline to the flank, making them an excellent model of biological pattern formation. This Primer explores recent evidence that a local Turing periodic patterning, a global propagation mechanism, and mechano-chemical coupling work together to pattern the feather array progressively.

Calcium oscillations coordinate feather mesenchymal cell movement by SHH dependent modulation of gap junction networks

Collective cell migration mediates multiple tissue morphogenesis processes. Yet how multi-dimensional mesenchymal cell movements are coordinated remains mostly unknown. Here we report that coordinated mesenchymal cell migration during chicken feather elongation is accompanied by dynamic changes of bioelectric currents. Transcriptome profiling and functional assays implicate contributions from functional voltage-gated Ca²⁺ channels (VGCCs), Connexin-43 based gap junctions, and Ca²⁺ release activated Ca²⁺ (CRAC) channels. 4-Dimensional Ca²⁺ imaging reveals that the Sonic hedgehog-responsive mesenchymal cells display synchronized Ca²⁺ oscillations, which expand progressively in area during feather elongation. Inhibiting VGCCs, gap junctions, or Sonic hedgehog signaling alters the mesenchymal Ca²⁺ landscape, cell movement patterns and feather bud elongation. Ca²⁺ oscillations induced by cyclic activation of opto-cCRAC channels enhance feather bud elongation. Functional disruption experiments and promoter analysis implicate synergistic Hedgehog and WNT/β-Catenin signaling in activating Connexin-43 expression, establishing gap junction networks synchronizing the Ca²⁺ profile among cells, thereby coordinating cell movement patterns.

The molecular mechanisms regulating mesenchymal cell movements are unclear. Here, the authors show in chicken feather elongation that SHH/WNT signalling establishes gap-junctions, enabling synchronized Ca^2 + oscillations to emerge for and in turn coordinating directed cell migration.

Hair follicle dermal condensation forms via Fgf20 primed cell cycle exit, cell motility, and aggregation

Mesenchymal condensation is a critical step in organogenesis, yet the underlying molecular and cellular mechanisms remain poorly understood. The hair follicle dermal condensate is the precursor to the permanent mesenchymal unit of the hair follicle, the dermal papilla, which regulates hair cycling throughout life and bears hair inductive potential. Dermal condensate morphogenesis depends on epithelial Fibroblast Growth Factor 20 (Fgf20). Here, we combine mouse models with 3D and 4D microscopy to demonstrate that dermal condensates form de novo and via directional migration. We identify cell cycle exit and cell shape changes as early hallmarks of dermal condensate morphogenesis and find that Fgf20 primes these cellular behaviors and enhances cell motility and condensation. RNAseq profiling of immediate Fgf20 targets revealed induction of a subset of dermal condensate marker genes. Collectively, these data indicate that dermal condensation occurs via directed cell movement and that Fgf20 orchestrates the early cellular and molecular events.

All mammal hair springs from hair follicles under the skin. These follicles sit in the dermis, beneath the outermost skin layer, the epidermis. In the embryo, hair follicles develop from unspecialized cells in two tissues, the epithelium and the mesenchyme, which will later develop into the dermis and epidermis, respectively. As development progresses, the cells of these tissues begin to cluster, and signals passing back and forth between the epithelium and mesenchyme instruct the cells what to do. In the mesenchyme, cells called fibroblasts squeeze up against their neighbors, forming patches called dermal condensates. These mature into so-called dermal papillae, which supply specific molecules called growth factors that regulate hair formation throughout lifetime.

Fibroblasts in the developing skin respond to a signal from the epithelium called fibroblast growth factor 20 (Fgf20), but we do not yet understand its effects. It is possible that Fgf20 tells the cells to divide, forming clusters of daughter cells around their current location. Or, it could be that Fgf20 tells the cells to move, encouraging them to travel towards one another to form groups.

To address this question, Biggs, Mäkelä et al. examined developing mouse skin grown in the laboratory. They traced cells marked with fluorescent tags to analyze their behavior as the condensates formed. This revealed that the Fgf20 signal acts as a rallying call, triggering fibroblast movement. The cells changed shape and moved towards one another, rather than dividing to create their own clusters. In fact, they switched off their own cell cycle as the condensates formed, halting their ability to divide. A technique called RNA sequencing revealed that Fgf20 also promotes the use of genes known to be active in dermal condensates.

Dermal papillae control hair growth, and transplanting them under the skin can form new hair follicles. However, these cells lose this ability when grown in the laboratory. Understanding how they develop could be beneficial for future hair growth therapy. Further work could also address fundamental questions in embryology. Condensates of cells from the mesenchyme also precede the formation of limbs, bones, muscles and organs. Extending this work could help us to understand this critical developmental step.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

Hierarchical patterning modes orchestrate hair follicle morphogenesis

Two theories address the origin of repeating patterns, such as hair follicles, limb digits, and intestinal villi, during development. The Turing reaction-diffusion system posits that interacting diffusible signals produced by static cells first define a prepattern that then induces cell rearrangements to produce an anatomical structure. The second theory, that of mesenchymal self-organisation, proposes that mobile cells can form periodic patterns of cell aggregates directly, without reference to any prepattern. Early hair follicle development is characterised by the rapid appearance of periodic arrangements of altered gene expression in the epidermis and prominent clustering of the adjacent dermal mesenchymal cells. We assess the contributions and interplay between reaction-diffusion and mesenchymal self-organisation processes in hair follicle patterning, identifying a network of fibroblast growth factor (FGF), wingless-related integration site (WNT), and bone morphogenetic protein (BMP) signalling interactions capable of spontaneously producing a periodic pattern. Using time-lapse imaging, we find that mesenchymal cell condensation at hair follicles is locally directed by an epidermal prepattern. However, imposing this prepattern's condition of high FGF and low BMP activity across the entire skin reveals a latent dermal capacity to undergo spatially patterned self-organisation in the absence of epithelial direction. This mesenchymal self-organisation relies on restricted transforming growth factor (TGF) β signalling, which serves to drive chemotactic mesenchymal patterning when reaction-diffusion patterning is suppressed, but, in normal conditions, facilitates cell movement to locally prepatterned sources of FGF. This work illustrates a hierarchy of periodic patterning modes operating in organogenesis.

A spectrum of modularity in multi-functional gene circuits

A major challenge in systems biology is to understand the relationship between a circuit's structure and its function, but how is this relationship affected if the circuit must perform multiple distinct functions within the same organism? In particular, to what extent do multi‐functional circuits contain modules which reflect the different functions? Here, we computationally survey a range of bi‐functional circuits which show no simple structural modularity: They can switch between two qualitatively distinct functions, while both functions depend on all genes of the circuit. Our analysis reveals two distinct classes: hybrid circuits which overlay two simpler mono‐functional sub‐circuits within their circuitry, and emergent circuits, which do not. In this second class, the bi‐functionality emerges from more complex designs which are not fully decomposable into distinct modules and are consequently less intuitive to predict or understand. These non‐intuitive emergent circuits are just as robust as their hybrid counterparts, and we therefore suggest that the common bias toward studying modular systems may hinder our understanding of real biological circuits.

Genomic determinants of epidermal appendage patterning and structure in domestic birds

Variation in regional identity, patterning, and structure of epidermal appendages contributes to skin diversity among many vertebrate groups, and is perhaps most striking in birds. In pioneering work on epidermal appendage patterning, John Saunders and his contemporaries took advantage of epidermal appendage diversity within and among domestic chicken breeds to establish the importance of mesoderm-ectoderm signaling in determining skin patterning. Diversity in chickens and other domestic birds, including pigeons, is driving a new wave of research to dissect the molecular genetic basis of epidermal appendage patterning. Domestic birds are not only outstanding models for embryonic manipulations, as Saunders recognized, but they are also ideal genetic models for discovering the specific genes that control normal development and the mutations that contribute to skin diversity. Here, we review recent genetic and genomic approaches to uncover the basis of epidermal macropatterning, micropatterning, and structural variation. We also present new results that confirm expression changes in two limb identity genes in feather-footed pigeons, a case of variation in appendage structure and identity.

Pattern selection by dynamical biochemical signals

The development of multicellular organisms involves cells to decide their fate upon the action of biochemical signals. This decision is often spatiotemporally coordinated such that a spatial pattern arises. The dynamics that drive pattern formation usually involve genetic nonlinear interactions and positive feedback loops. These complex dynamics may enable multiple stable patterns for the same conditions. Under these circumstances, pattern formation in a developing tissue involves a selection process: why is a certain pattern formed and not another stable one? Herein we computationally address this issue in the context of the Notch signaling pathway. We characterize a dynamical mechanism for developmental selection of a specific pattern through spatiotemporal changes of the control parameters of the dynamics, in contrast to commonly studied situations in which initial conditions and noise determine which pattern is selected among multiple stable ones. This mechanism can be understood as a path along the parameter space driven by a sequence of biochemical signals. We characterize the selection process for three different scenarios of this dynamical mechanism that can take place during development: the signal either 1) acts in all the cells at the same time, 2) acts only within a cluster of cells, or 3) propagates along the tissue. We found that key elements for pattern selection are the destabilization of the initial pattern, the subsequent exploration of other patterns determined by the spatiotemporal symmetry of the parameter changes, and the speeds of the path compared to the timescales of the pattern formation process itself. Each scenario enables the selection of different types of patterns and creates these elements in distinct ways, resulting in different features. Our approach extends the concept of selection involved in cellular decision-making, usually applied to cell-autonomous decisions, to systems that collectively make decisions through cell-to-cell interactions.

Regulatory Differences in Natal Down Development between Altricial Zebra Finch and Precocial Chicken

Birds can be classified into altricial and precocial. The hatchlings of altricial birds are almost naked, whereas those of precocial birds are covered with natal down. This regulatory divergence is thought to reflect environmental adaptation, but the molecular basis of the divergence is unclear. To address this issue, we chose the altricial zebra finch and the precocial chicken as the model animals. We noted that zebra finch hatchlings show natal down growth suppressed anterior dorsal (AD) skin but partially down-covered posterior dorsal (PD) skin. Comparing the transcriptomes of AD and PD skins, we found that the feather growth promoter SHH (sonic hedgehog) was expressed higher in PD skin than in AD skin. Moreover, the data suggested that the FGF (fibroblast growth factor)/Mitogen-activated protein kinase (MAPK) signaling pathway is involved in natal down growth suppression and that FGF16 is a candidate upstream signaling suppressor. Ectopic expression of FGF16 on chicken leg skin showed downregulation of SHH, upregulation of the feather growth suppressor FGF10, and suppression of feather bud elongation, similar to the phenotype found in zebra finch embryonic AD skin. Therefore, we propose that FGF16-related signals suppress natal down elongation and cause the naked AD skin in zebra finch. Our study provides insights into the regulatory divergence in natal down formation between precocial and altricial birds.

Generation of bioengineered feather buds on a reconstructed chick skin from dissociated epithelial and mesenchymal cells

Various kinds of in vitro culture systems of tissues and organs have been developed, and applied to understand multicellular systems during embryonic organogenesis. In the research field of feather bud development, tissue recombination assays using an intact epithelial tissue and mesenchymal tissue/cells have contributed to our understanding the mechanisms of feather bud formation and development. However, there are few methods to generate a skin and its appendages from single cells of both epithelium and mesenchyme. In this study, we have developed a bioengineering method to reconstruct an embryonic dorsal skin after completely dissociating single epithelial and mesenchymal cells from chick skin. Multiple feather buds can form on the reconstructed skin in a single row in vitro. The bioengineered feather buds develop into long feather buds by transplantation onto a chorioallantoic membrane. The bioengineered bud sizes were similar to those of native embryo. The number of bioengineered buds was increased linearly with the initial contact length of epithelial and mesenchymal cell layers where the epithelial-mesenchymal interactions occur. In addition, the bioengineered bud formation was also disturbed by the inhibition of major signaling pathways including FGF (fibroblast growth factor), Wnt/β-catenin, Notch and BMP (bone morphogenetic protein). We expect that our bioengineering technique will motivate further extensive research on multicellular developmental systems, such as the formation and sizing of cutaneous appendages, and their regulatory mechanisms.

Cells or signals: which moves to drive skin pattern formation?

During its development, the skin produces an array of evenly spaced hair follicles. How the location of each follicle is determined to produce this pattern has been the subject of study and speculation for several decades. A central unresolved issue is the extent to which movement of scattered, precommitted follicle cells might play a role in this process. Xavier et al. now report the identification of subpopulations of dermal cells in developing sheep skin which are positive for Delta1 expression, suggesting that these cells may represent precommitted dermal papilla cells and that dermal Notch pathway signalling plays a role in hair follicle patterning.

Modelling from the experimental developmental biologists viewpoint

In this review we consider Reaction-Diffusion as the archetype of a model in developmental biology. We consider its history in relation to experimental work since it was first proposed in 1952 by Turing and revived in the 1970s by Meinhardt. We then discuss the most recent examples of experiments that address this model, including the challenges that remain in capturing the physico-chemical manifestation of the model mechanism in a real developmental system. Finally we discuss the model's current status and use in the experimental community.

Predicting the spatiotemporal dynamics of hair follicle patterns in the developing mouse

Reaction-diffusion models have been used as a paradigm for describing the de novo emergence of biological patterns such as stripes and spots. In many organisms, these initial patterns are typically refined and elaborated over the subsequent course of development. Here we study the formation of secondary hair follicle patterns in the skin of developing mouse embryos. We used the expression of sex-determining region Y box 2 to identify and distinguish the primary and secondary hair follicles and to infer the spatiotemporal dynamics of the follicle formation process. Quantitative analysis of the specific follicle patterns observed reveals a simple geometrical rule governing the formation of secondary follicles, and motivates an expansion-induction (EI) model in which new follicle formation is driven by the physical growth of the embryo. The EI model requires only one diffusible morphogen and provides quantitative, accurate predictions on the relative positions and timing of secondary follicle formation, using only the observed configuration of primary follicles as input. The same model accurately describes the positions of additional follicles that emerge from skin explants treated with an activator. Thus, the EI model provides a simple and robust mechanism for predicting secondary space-filling patterns in growing embryos.

Roles of EphB3/ephrin-B1 in feather morphogenesis

The ephrin receptor (Eph) tyrosine kinases and their ephrin ligands are involved in morphogenesis during organ formation. We studied their role in feather morphogenesis, focusing on ephrin-B1 and its receptor EphB3. Early in feather development, ephrin-B1 mRNA and protein were found to be expressed in the dermal condensation, but not in the inter-bud mesenchyme. Later, in feather buds, expression was found in both the epithelium and mesenchyme. In the feather follicle, ephrin-B1 protein expression was found to be enriched in the feather filament epithelium and in the marginal plate which sets the boundary between the barb ridges. EphB3 mRNA was also expressed in epithelia. In the feather bud, its expression was restricted to the posterior bud. In the follicle, its expression formed a circle at the bud base which may set the boundary between bud and inter-bud domains. Perturbation with ephrin-B1/Fc altered feather primordia segregation and feather bud elongation. Analyses revealed that ephrin-B1/Fc caused three types of changes: blurred placode boundaries with loose dermal condensations, incomplete follicle invagination with less compact dermal papillae, and aberrant barb ridge patterning in feather filament morphogenesis. Thus, while ephrin-B1 suppression does not inhibit the initial emergence of a new epithelial domain, Eph/ephrin-B1 interaction is required for its proper completion. Consequently, we propose that interaction between ephrin-B1 and its receptor is involved in boundary stabilization during feather morphogenesis.

Biology of tooth replacement in amniotes

Tooth replacement is a common trait to most vertebrates, including mammals. Mammals, however, have lost the capacity for continuous tooth renewal seen in most other vertebrates, and typically have only 1–2 generations of teeth. Here, we review the mechanisms of tooth replacement in reptiles and mammals, and discuss in detail the current and historical theories on control of timing and pattern of tooth replacement and development.

Multiscale mechanobiology: mechanics at the molecular, cellular, and tissue levels

Mechanical force is present in all aspects of living systems. It affects the conformation of molecules, the shape of cells, and the morphology of tissues. All of these are crucial in architecture-dependent biological functions. Nanoscience of advanced materials has provided knowledge and techniques that can be used to understand how mechanical force is involved in biological systems, as well as to open new avenues to tailor-made bio-mimetic materials with desirable properties.

In this article, we describe models and show examples of how force is involved in molecular functioning, cell shape patterning, and tissue morphology.

Fgf20 governs formation of primary and secondary dermal condensations in developing hair follicles

In hair follicle development, a placode-derived signal is believed to induce formation of the dermal condensation, an essential component of ectodermal organs. However, the identity of this signal is unknown. Furthermore, although induction and patterning of hair follicles are intimately linked, it is not known whether the mesenchymal condensation is necessary for inducing the initial epithelial pattern. Here, we show that fibroblast growth factor 20 (Fgf20) is expressed in hair placodes and is induced by and functions downstream from epithelial ectodysplasin (Eda)/Edar and Wnt/β-Catenin signaling to initiate formation of the underlying dermal condensation. Fgf20 governs formation of primary and secondary dermal condensations in developing hair follicles and subsequent formation of guard, awl, and auchene hairs. Although primary dermal condensations are absent in Fgf20 mutant mice, a regular array of hair placodes is formed, demonstrating that the epithelial patterning process is independent of known histological and molecular markers of underlying mesenchymal patterns during the initial stages of hair follicle development.

Module-based complexity formation: periodic patterning in feathers and hairs

Patterns describe order which emerges from homogeneity. Complex patterns on the integument are striking because of their visibility throughout an organism’s lifespan. Periodic patterning is an effective design because the ensemble of hair or feather follicles (modules) allows the generation of complexity, including regional variations and cyclic regeneration, giving the skin appendages a new lease on life. Spatial patterns include the arrangements of feathers and hairs in specific number, size, and spacing.We explorehowa field of equivalent progenitor cells can generate periodically arranged modules based on genetic information, physical–chemical rules and developmental timing. Reconstitution experiments suggest a competitive equilibrium regulated by activators/inhibitors involving Turing reaction-diffusion. Temporal patterns result from oscillating stem cell activities within each module (microenvironment regulation), reflected as growth (anagen) and resting (telogen) phases during the cycling of feather and hair follicles. Stimulating modules with activators initiates the spread of regenerative hair waves, while global inhibitors outside each module (macroenvironment) prevent this. Different wave patterns can be simulated by cellular automata principles. Hormonal status and seasonal changes can modulate appendage phenotypes, leading to ‘organ metamorphosis’, with multiple ectodermal organ phenotypes generated from the same precursors. We discuss potential novel evolutionary steps using this module-based complexity in several amniote integument organs, exemplified by the spectacular peacock feather pattern. We thus explore the application of the acquired knowledge of patterning in tissue engineering. New hair follicles can be generated after wounding. Hairs and feathers can be reconstituted through self-organization of dissociated progenitor cells.

Sprouty/FGF signaling regulates the proximal-distal feather morphology and the size of dermal papillae

In a feather, there are distinct morphologies along the proximal-distal axis. The proximal part is a cylindrical stalk (calamus), whereas the distal part has barb and barbule branches. Here we focus on what molecular signaling activity can modulate feather stem cells to generate these distinct morphologies. We demonstrate the drastic tissue remodeling during feather cycling which includes initiation, growth and resting phases. In the growth phase, epithelial components undergo progressive changes from the collar growth zone to the ramogenic zone, to maturing barb branches along the proximal-distal axis. Mesenchymal components also undergo progressive changes from the dermal papilla, to the collar mesenchyme, to the pulp along the proximal-distal axis. Over-expression of Spry4, a negative regulator of receptor tyrosine kinases, promotes barb branch formation at the expense of the epidermal collar. It even induces barb branches from the follicle sheath (equivalent to the outer root sheath in hair follicles). The results are feathers with expanded feather vane regions and small or missing proximal feather shafts (the calamus). Spry4 also expands the pulp region while reducing the size of dermal papillae, leading to a failure to regenerate. In contrast, over-expressing Fgf10 increases the size of the dermal papillae, expands collar epithelia and mesenchyme, but also prevents feather branch formation and feather keratin differentiation. These results suggest that coordinated Sprouty/FGF pathway activity at different stages is important to modulate feather epidermal stem cells to form distinct feather morphologies along the proximal-distal feather axis.

Focal high cell density generates a gradient of patterns in self-organizing vascular mesenchymal cells

In embryogenesis, structural patterns, such as vascular branching, may form via a reaction-diffusion mechanism in which activator and inhibitor morphogens guide cells into periodic aggregates. We previously found that vascular mesenchymal cells (VMCs) spontaneously aggregate into nodular structures and that morphogen pairs regulate the aggregation into patterns of spots and stripes. To test the effect of a focal change in activator morphogen on VMC pattern formation, we created a focal zone of high cell density by plating a second VMC layer within a cloning ring over a confluent monolayer. After 24 h, the ring was removed and pattern formation monitored by phase-contrast microscopy. At days 2-8, the patterns progressed from uniform distributions to swirl, labyrinthine and spot patterns. Within the focal high-density zone (HDZ) and a narrow halo zone, cells aggregated into spot patterns, whilst in the outermost zone of the plate, cells formed a labyrinthine pattern. The area occupied by aggregates was significantly greater in the outermost zone than in the HDZ or halo. The rate of pattern progression within the HDZ increased as a function of its plating density. Thus, focal differences in cell density may drive pattern formation gradients in tissue architecture, such as vascular branching.

Effects of intrinsic stochasticity on delayed reaction-diffusion patterning systems

Cellular gene expression is a complex process involving many steps, including the transcription of DNA and translation of mRNA; hence the synthesis of proteins requires a considerable amount of time, from ten minutes to several hours. Since diffusion-driven instability has been observed to be sensitive to perturbations in kinetic delays, the application of Turing patterning mechanisms to the problem of producing spatially heterogeneous differential gene expression has been questioned. In deterministic systems a small delay in the reactions can cause a large increase in the time it takes a system to pattern. Recently, it has been observed that in undelayed systems intrinsic stochasticity can cause pattern initiation to occur earlier than in the analogous deterministic simulations. Here we are interested in adding both stochasticity and delays to Turing systems in order to assess whether stochasticity can reduce the patterning time scale in delayed Turing systems. As analytical insights to this problem are difficult to attain and often limited in their use, we focus on stochastically simulating delayed systems. We consider four different Turing systems and two different forms of delay. Our results are mixed and lead to the conclusion that, although the sensitivity to delays in the Turing mechanism is not completely removed by the addition of intrinsic noise, the effects of the delays are clearly ameliorated in certain specific cases.

Towards an integrated experimental-theoretical approach for assessing the mechanistic basis of hair and feather morphogenesis

In his seminal 1952 paper, 'The Chemical Basis of Morphogenesis', Alan Turing lays down a milestone in the application of theoretical approaches to understand complex biological processes. His deceptively simple demonstration that a system of reacting and diffusing chemicals could, under certain conditions, generate spatial patterning out of homogeneity provided an elegant solution to the problem of how one of nature's most intricate events occurs: the emergence of structure and form in the developing embryo. The molecular revolution that has taken place during the six decades following this landmark publication has now placed this generation of theoreticians and biologists in an excellent position to rigorously test the theory and, encouragingly, a number of systems have emerged that appear to conform to some of Turing's fundamental ideas. In this paper, we describe the history and more recent integration between experiment and theory in one of the key models for understanding pattern formation: the emergence of feathers and hair in the skins of birds and mammals.

Genesis of two-dimensional patterns in cross-gradient fields

Tissue morphogenesis is controlled by the two-dimensional patterning of gene expression in epithelial layers, that determines cell fates. The mechanisms of pattern formation involve intracellular regulatory networks controlled by paracrine and autocrine signaling. We develop a general logical scheme to deduce the morphology of two-dimensional patterns in the field of two crossed morphogen gradients enriched by the action of autocrine signaling that may subdivide expression domains in nontrivial ways. A variety of persistent patterns, either stationary or oscillatory, are generated using the various interaction schemes, some of which have been generated by a special algorithm including random inputs and selected according to suitable criteria. We give the full classification of a variety of stationary and oscillatory expression patterns in the presence of a single autocrine signal based on logical arguments. These results are further confirmed by numerical computations based on reaction-diffusion equations for morphogens and ligands and the discrete (cell-level) description of intracellular dynamics. Model simulations also elucidate transient processes, in particular interaction schemes. Different internal schemes may lead to identical persistent patterns, although relaxation may proceed in distinct ways. As an illustration of the general method, we test a particular scheme suggested by genetic studies of dynamic gene expression patterns in the follicular epithelium of the Drosophila eggshell.

The influence of receptor-mediated interactions on reaction-diffusion mechanisms of cellular self-organisation

Understanding the mechanisms governing and regulating self-organisation in the developing embryo is a key challenge that has puzzled and fascinated scientists for decades. Since its conception in 1952 the Turing model has been a paradigm for pattern formation, motivating numerous theoretical and experimental studies, though its verification at the molecular level in biological systems has remained elusive. In this work, we consider the influence of receptor-mediated dynamics within the framework of Turing models, showing how non-diffusing species impact the conditions for the emergence of self-organisation. We illustrate our results within the framework of hair follicle pre-patterning, showing how receptor interaction structures can be constrained by the requirement for patterning, without the need for detailed knowledge of the network dynamics. Finally, in the light of our results, we discuss the ability of such systems to pattern outside the classical limits of the Turing model, and the inherent dangers involved in model reduction.

Stochastic reaction and diffusion on growing domains: understanding the breakdown of robust pattern formation

Many biological patterns, from population densities to animal coat markings, can be thought of as heterogeneous spatiotemporal distributions of mobile agents. Many mathematical models have been proposed to account for the emergence of this complexity, but, in general, they have consisted of deterministic systems of differential equations, which do not take into account the stochastic nature of population interactions. One particular, pertinent criticism of these deterministic systems is that the exhibited patterns can often be highly sensitive to changes in initial conditions, domain geometry, parameter values, etc. Due to this sensitivity, we seek to understand the effects of stochasticity and growth on paradigm biological patterning models. In this paper, we extend spatial Fourier analysis and growing domain mapping techniques to encompass stochastic Turing systems. Through this we find that the stochastic systems are able to realize much richer dynamics than their deterministic counterparts, in that patterns are able to exist outside the standard Turing parameter range. Further, it is seen that the inherent stochasticity in the reactions appears to be more important than the noise generated by growth, when considering which wave modes are excited. Finally, although growth is able to generate robust pattern sequences in the deterministic case, we see that stochastic effects destroy this mechanism for conferring robustness. However, through Fourier analysis we are able to suggest a reason behind this lack of robustness and identify possible mechanisms by which to reclaim it.

From buds to follicles: matrix metalloproteinases in developmental tissue remodeling during feather morphogenesis

Organogenesis involves a series of dynamic morphogenesis and remodeling processes. Since feathers exhibit complex forms, we have been using the feather as a model to analyze how molecular pathways and cellular events are used. While several major molecular pathways have been studied, the roles of matrix degrading proteases and inhibitors in feather morphogenesis are unknown. Here we addressed this knowledge gap by studying the temporal and spatial expression of proteases and inhibitors in developing feathers using mammalian antibodies that cross react with chicken proteins. We also investigated the effect of protease inhibitors on feather development employing an in vitro feather bud culture system. The results show that antibodies specific for mammalian MMP2 and TIMP2 stained positive in both feather epithelium and mesenchyme. The staining co-localized in structures of E10-E13 developing feathers. Interestingly, MMP2 and TIMP2 exhibited a complementary staining pattern in developing E15 and E20 feathers and in maturing feather filaments. Although they exhibited a slight delay in feather bud development, similar patterns of MMP2 and TIMP2 staining were observed in in vitro culture explants. The broad spectrum pharmacological inhibitors AG3340 and BB103 (MMP inhibitors) but not Aprotinin (a plasmin inhibitor) showed a reversible effect on epithelium invagination and feather bud elongation. TIMP2, a physiological inhibitor to MMPs, exhibited a similar effect. Markers of feather morphogenesis showed that MMP activity was required for both epithelium invagination and mesenchymal cell proliferation. Inhibition of MMP activity led to an overall delay in the expression of molecules that regulate either early feather bud growth and/or differentiation and thereby produced abnormal buds with incomplete follicle formation. This work demonstrates that MMPs and their inhibitors are not only important in injury repair, but also in development tissue remodeling as demonstrated here for the formation of feather follicles.

In search of the Golden Fleece: unraveling principles of morphogenesis by studying the integrative biology of skin appendages

The mythological story of the Golden Fleece symbolizes the magical regenerative power of skin appendages. Similar to the adventurous pursuit of the Golden Fleece by the multi-talented Argonauts, today we also need an integrated multi-disciplined approach to understand the cellular and molecular processes during development, regeneration and evolution of skin appendages. To this end, we have explored several aspects of skin appendage biology that contribute to the Turing activator/inhibitor model in feather pattern formation, the topo-biological arrangement of stem cells in organ shape determination, the macro-environmental regulation of stem cells in regenerative hair waves, and potential novel molecular pathways in the morphological evolution of feathers. Here we show our current integrative biology efforts to unravel the complex cellular behavior in patterning stem cells and the control of regional specificity in skin appendages. We use feather/scale tissue recombination to demonstrate the timing control of competence and inducibility. Feathers from different body regions are used to study skin regional specificity. Bioinformatic analyses of transcriptome microarrays show the potential involvement of candidate molecular pathways. We further show Hox genes exhibit some region specific expression patterns. To visualize real time events, we applied time-lapse movies, confocal microscopy and multiphoton microscopy to analyze the morphogenesis of cultured embryonic chicken skin explants. These modern imaging technologies reveal unexpectedly complex cellular flow and organization of extracellular matrix molecules in three dimensions. While these approaches are in preliminary stages, this perspective highlights the challenges we face and new integrative tools we will use. Future work will follow these leads to develop a systems biology view and understanding in the morphogenetic principles that govern the development and regeneration of ectodermal organs.

Cryptic patterning of avian skin confers a developmental facility for loss of neck feathering

Vertebrate skin is characterized by its patterned array of appendages, whether feathers, hairs, or scales. In avian skin the distribution of feathers occurs on two distinct spatial levels. Grouping of feathers within discrete tracts, with bare skin lying between the tracts, is termed the macropattern, while the smaller scale periodic spacing between individual feathers is referred to as the micropattern. The degree of integration between the patterning mechanisms that operate on these two scales during development and the mechanisms underlying the remarkable evolvability of skin macropatterns are unknown. A striking example of macropattern variation is the convergent loss of neck feathering in multiple species, a trait associated with heat tolerance in both wild and domestic birds. In chicken, a mutation called Naked neck is characterized by a reduction of body feathering and completely bare neck. Here we perform genetic fine mapping of the causative region and identify a large insertion associated with the Naked neck trait. A strong candidate gene in the critical interval, BMP12/GDF7, displays markedly elevated expression in Naked neck embryonic skin due to a cis-regulatory effect of the causative mutation. BMP family members inhibit embryonic feather formation by acting in a reaction-diffusion mechanism, and we find that selective production of retinoic acid by neck skin potentiates BMP signaling, making neck skin more sensitive than body skin to suppression of feather development. This selective production of retinoic acid by neck skin constitutes a cryptic pattern as its effects on feathering are not revealed until gross BMP levels are altered. This developmental modularity of neck and body skin allows simple quantitative changes in BMP levels to produce a sparsely feathered or bare neck while maintaining robust feather patterning on the body.

Waves and patterning in developmental biology: vertebrate segmentation and feather bud formation as case studies

In this article we will discuss the integration of developmental patterning mechanisms with waves of competency that control the ability of a homogeneous field of cells to react to pattern forming cues and generate spatially heterogeneous patterns. We base our discussion around two well known patterning events that take place in the early embryo: somitogenesis and feather bud formation. We outline mathematical models to describe each patterning mechanism, present the results of numerical simulations and discuss the validity of each model in relation to our example patterning processes.