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

Conserved regulatory switches for the transition from natal down to juvenile feather in birds

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We report that extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial-mesenchymal interactions for branching morphogenesis. α-SMA (ACTA2) compartmentalizes dermal papilla stem cells for feather renewal cycling. LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We show that this primary feather transition is largely conserved in chicken (precocial) and zebra finch (altricial) and discuss the possibility that this evolutionary adaptation process started in feathered dinosaurs.

The avian ectodermal default competence to make feathers

Feathers originate as protofeathers before birds, in pterosaurs and basal dinosaurs. What characterizes a feather is not only its outgrowth, but its barb cells differentiation and a set of beta-corneous proteins. Reticula appear concomitantly with feathers, as small bumps on plantar skin, made only of keratins. Avian scales, with their own set of beta-corneous proteins, appear more recently than feathers on the shank, and only in some species. In the chick embryo, when feather placodes form, all the non-feather areas of the integument are already specified. Among them, midventral apterium, cornea, reticula, and scale morphogenesis appear to be driven by negative regulatory mechanisms, which modulate the inherited capacity of the avian ectoderm to form feathers. Successive dermal/epidermal interactions, initiated by the Wnt/β-catenin pathway, and involving principally Eda/Edar, BMP, FGF20 and Shh signaling, are responsible for the formation not only of feather, but also of scale placodes and reticula, with notable differences in the level of Shh, and probably FGF20 expressions. This sequence is a dynamic and labile process, the turning point being the FGF20 expression by the placode. This epidermal signal endows its associated dermis with the memory to aggregate and to stimulate the morphogenesis that follows, involving even a re-initiation of the placode.

Transition from natal downs to juvenile feathers: conserved regulatory switches in Neoaves

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We discovered that LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial -mesenchymal interactions for branching morphogenesis. ACTA2 compartments dermal papilla stem cells for feather cycling. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We found this primary feather transition largely conserved in chicken (precocious) and zebra finch (altricial) and discussed the possibility that this evolutionary adaptation process started in feathered dinosaurs.

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

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

Bat teeth illuminate the diversification of mammalian tooth classes

Tooth classes are an innovation that has contributed to the evolutionary success of mammals. However, our understanding of the mechanisms by which tooth classes diversified remain limited. We use the evolutionary radiation of noctilionoid bats to show how the tooth developmental program evolved during the adaptation to new diet types. Combining morphological, developmental and mathematical modeling approaches, we demonstrate that tooth classes develop through independent developmental cascades that deviate from classical models. We show that the diversification of tooth number and size is driven by jaw growth rate modulation, explaining the rapid gain/loss of teeth in this clade. Finally, we mathematically model the successive appearance of tooth buds, supporting the hypothesis that growth acts as a key driver of the evolution of tooth number and size. Our work reveal how growth, by tinkering with reaction/diffusion processes, drives the diversification of tooth classes and other repeated structure during adaptive radiations.

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.

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

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

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.

Cyclic growth of dermal papilla and regeneration of follicular mesenchymal components during feather cycling

How dermis maintains tissue homeostasis in cyclic growth and wounding is a fundamental unsolved question. Here, we study how dermal components of feather follicles undergo physiological (molting) and plucking injury-induced regeneration in chickens. Proliferation analyses reveal quiescent, transient-amplifying (TA) and long-term label-retaining dermal cell (LRDC) states. During the growth phase, LRDCs are activated to make new dermal components with distinct cellular flows. Dermal TA cells, enriched in the proximal follicle, generate both peripheral pulp, which extends distally to expand the epithelial-mesenchymal interactive interface for barb patterning, and central pulp, which provides nutrition. Entering the resting phase, LRDCs, accompanying collar bulge epidermal label-retaining cells, descend to the apical dermal papilla. In the next cycle, these apical dermal papilla LRDCs are re-activated to become new pulp progenitor TA cells. In the growth phase, lower dermal sheath can generate dermal papilla and pulp. Transcriptome analyses identify marker genes and highlight molecular signaling associated with dermal specification. We compare the cyclic topological changes with those of the hair follicle, a convergently evolved follicle configuration. This work presents a model for analyzing homeostasis and tissue remodeling of mesenchymal progenitors.

Symmetry breaking of tissue mechanics in wound induced hair follicle regeneration of laboratory and spiny mice

Tissue regeneration is a process that recapitulates and restores organ structure and function. Although previous studies have demonstrated wound-induced hair neogenesis (WIHN) in laboratory mice (Mus), the regeneration is limited to the center of the wound unlike those observed in African spiny (Acomys) mice. Tissue mechanics have been implicated as an integral part of tissue morphogenesis. Here, we use the WIHN model to investigate the mechanical and molecular responses of laboratory and African spiny mice, and report these models demonstrate opposing trends in spatiotemporal morphogenetic field formation with association to wound stiffness landscapes. Transcriptome analysis and K14-Cre-Twist1 transgenic mice show the Twist1 pathway acts as a mediator for both epidermal-dermal interactions and a competence factor for periodic patterning, differing from those used in development. We propose a Turing model based on tissue stiffness that supports a two-scale tissue mechanics process: (1) establishing a morphogenetic field within the wound bed (mm scale) and (2) symmetry breaking of the epidermis and forming periodically arranged hair primordia within the morphogenetic field (μm scale). Thus, we delineate distinct chemo-mechanical events in building a Turing morphogenesis-competent field during WIHN of laboratory and African spiny mice and identify its evo-devo advantages with perspectives for regenerative medicine.

Folding Keratin Gene Clusters during Skin Regional Specification

Regional specification is critical for skin development, regeneration, and evolution. The contribution of epigenetics in this process remains unknown. Here, using avian epidermis, we find two major strategies regulate β-keratin gene clusters. (1) Over the body, macro-regional specificities (scales, feathers, claws, etc.) established by typical enhancers control five subclusters located within the epidermal differentiation complex on chromosome 25; (2) within a feather, micro-regional specificities are orchestrated by temporospatial chromatin looping of the feather β-keratin gene cluster on chromosome 27. Analyses suggest a three-factor model for regional specification: competence factors (e.g., AP1) make chromatin accessible, regional specifiers (e.g., Zic1) target specific genome regions, and chromatin regulators (e.g., CTCF and SATBs) establish looping configurations. Gene perturbations disrupt morphogenesis and histo-differentiation. This chicken skin paradigm advances our understanding of how regulation of big gene clusters can set up a two-dimensional body surface map.

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.

Hair-bearing human skin generated entirely from pluripotent stem cells

The skin is a multi-layered organ equipped with appendages (i.e. follicles and glands) critical for regulating bodily fluid retention and temperature, guarding against external stresses, and mediating touch and pain sensation^1,2. Reconstruction of appendage-bearing skin in cultures and in bioengineered grafts remains an unmet biomedical challenge^3–9. Here, we report an organoid culture system that generates complex skin from human pluripotent stem cells. We use step-wise modulation of the TGFβ and FGF signalling pathways to co-induce cranial epithelial cells and neural crest cells within a spherical cell aggregate. During 4–5 months incubation, we observe the emergence of a cyst-like skin organoid composed of stratified epidermis, fat-rich dermis, and pigmented hair follicles equipped with sebaceous glands. A network of sensory neurons and Schwann cells form nerve-like bundles that target Merkel cells in organoid hair follicles, mimicking human touch circuitry. Single-cell RNA-sequencing and direct comparison to foetal specimens suggest that skin organoids are equivalent to human facial skin in the second-trimester of development. Moreover, we show that skin organoids form planar hair-bearing skin when grafted on nude mice. Together, our results demonstrate that nearly complete skin can self-assemble in vitro and be used to reconstitute skin in vivo. We anticipate skin organoids will be foundational to future studies of human skin development, disease modelling, or reconstructive surgery.

The Early Origin of Feathers

Feathers have long been regarded as the innovation that drove the success of birds. However, feathers have been reported from close dinosaurian relatives of birds, and now from ornithischian dinosaurs and pterosaurs, the cousins of dinosaurs. Incomplete preservation makes these reports controversial. If true, these findings shift the origin of feathers back 80 million years before the origin of birds. Gene regulatory networks show the deep homology of scales, feathers, and hairs. Hair and feathers likely evolved in the Early Triassic ancestors of mammals and birds, at a time when synapsids and archosaurs show independent evidence of higher metabolic rates (erect gait and endothermy), as part of a major resetting of terrestrial ecosystems following the devastating end-Permian mass extinction.

Self-assembly of biological networks via adaptive patterning revealed by avian intradermal muscle network formation

Networked structures integrate numerous elements into one functional unit, while providing a balance between efficiency, robustness, and flexibility. Understanding how biological networks self-assemble will provide insights into how these features arise. Here, we demonstrate how nature forms exquisite muscle networks that can repair, regenerate, and adapt to external perturbations using the feather muscle network in chicken embryos as a paradigm. The self-assembled muscle networks arise through the implementation of a few simple rules. Muscle fibers extend outward from feather buds in every direction, but only those muscle fibers able to connect to neighboring buds are eventually stabilized. After forming such a nearest-neighbor configuration, the network can be reconfigured, adapting to perturbed bud arrangement or mechanical cues. Our computational model provides a bioinspired algorithm for network self-assembly, with intrinsic or extrinsic cues necessary and sufficient to guide the formation of these regenerative networks. These robust principles may serve as a useful guide for assembling adaptive networks in other contexts.

How chemotherapy and radiotherapy damage the tissue: Comparative biology lessons from feather and hair models

Chemotherapy and radiotherapy are common modalities for cancer treatment. While targeting rapidly growing cancer cells, they also damage normal tissues and cause adverse effects. From the initial insult such as DNA double-strand break, production of reactive oxygen species (ROS) and a general stress response, there are complex regulatory mechanisms that control the actual tissue damage process. Besides apoptosis, a range of outcomes for the damaged cells are possible including cell cycle arrest, senescence, mitotic catastrophe, and inflammatory responses and fibrosis at the tissue level. Feather and hair are among the most actively proliferating (mini-)organs and are highly susceptible to both chemotherapy and radiotherapy damage, thus provide excellent, experimentally tractable model systems for dissecting how normal tissues respond to such injuries. Taking a comparative biology approach to investigate this has turned out to be particularly productive. Started in chicken feather and then extended to murine hair follicles, it was revealed that in addition to p53-mediated apoptosis, several other previously overlooked mechanisms are involved. Specifically, Shh, Wnt, mTOR, cytokine signalling and ROS-mediated degradation of adherens junctions have been implicated in the damage and/or reparative regeneration process. Moreover, we show here that inflammatory responses, which can be prominent upon histological examination of chemo- or radiotherapy-damaged hair follicle, may not be essential for the hair loss phenotype. These studies point to fundamental, evolutionarily conserved mechanisms in controlling tissue responses in vivo, and suggest novel strategies for the prevention and management of adverse effects that arise from chemo- or radiotherapy.

Morpho-regulation in diverse chicken feather formation: Integrating branching modules and sex hormone-dependent morpho-regulatory modules

Many animals can change the size, shape, texture and color of their regenerated coats in response to different ages, sexes, or seasonal environmental changes. Here, we propose that the feather core branching morphogenesis module can be regulated by sex hormones or other environmental factors to change feather forms, textures or colors, thus generating a large spectrum of complexity for adaptation. We use sexual dimorphisms of the chicken to explore the role of hormones. A long-standing question is whether the sex-dependent feather morphologies are autonomously controlled by the male or female cell types, or extrinsically controlled and reversible. We have recently identified core feather branching molecular modules which control the anterior-posterior (bone morphogenetic orotein [BMP], Wnt gradient), medio-lateral (Retinoic signaling, Gremlin), and proximo-distal (Sprouty, BMP) patterning of feathers. We hypothesize that morpho-regulation, through quantitative modulation of existing parameters, can act on core branching modules to topologically tune the dimension of each parameter during morphogenesis and regeneration. Here, we explore the involvement of hormones in generating sexual dimorphisms using exogenously delivered hormones. Our strategy is to mimic male androgen levels by applying exogenous dihydrotestosterone and aromatase inhibitors to adult females and to mimic female estradiol levels by injecting exogenous estradiol to adult males. We also examine differentially expressed genes in the feathers of wildtype male and female chickens to identify potential downstream modifiers of feather morphogenesis. The data show male and female feather morphology and their color patterns can be modified extrinsically through molting and resetting the stem cell niche during regeneration.

Hoxc-Dependent Mesenchymal Niche Heterogeneity Drives Regional Hair Follicle Regeneration

Mesenchymal niche cells instruct activity of tissue-resident stem and progenitor cell populations. Epithelial stem cells in hair follicles (HFs) have region-specific activity, which may arise from intrinsic cellular heterogeneity within mesenchymal dermal papilla (DP) cells. Here we show that expression of Hoxc genes is sufficient to reprogram mesenchymal DP cells and alter the regenerative potential of epithelial stem cells. Hoxc gene expression in adult skin dermis closely correlates with regional HF regeneration patterns. Disrupting the region-specific expression patterns of Hoxc genes, by either decreasing their epigenetic repression via Bmi1 loss or inducing ectopic interactions of the Hoxc locus with an active epigenetic region, leads to precocious HF regeneration. We further show that a single Hoxc gene is sufficient to activate dormant DP niches and promote regional HF regeneration through canonical Wnt signaling. Altogether, these results reveal that Hoxc genes bestow mesenchymal niches with tissue-level heterogeneity and plasticity.

Fishing for ancestry

The same genes and signalling pathways control the formation of skin appendages in both fish and land animals.

Wnt/β-catenin regulates an ancient signaling network during zebrafish scale development

Understanding how patterning influences cell behaviors to generate three dimensional morphologies is a central goal of developmental biology. Additionally, comparing these regulatory mechanisms among morphologically diverse tissues allows for rigorous testing of evolutionary hypotheses. Zebrafish skin is endowed with a coat of precisely patterned bony scales. We use in-toto live imaging during scale development and manipulations of cell signaling activity to elucidate core features of scale patterning and morphogenesis. These analyses show that scale development requires the concerted activity of Wnt/β-catenin, Ectodysplasin (Eda) and Fibroblast growth factor (Fgf) signaling. This regulatory module coordinates Hedgehog (HH) dependent collective cell migration during epidermal invagination, a cell behavior not previously implicated in skin appendage morphogenesis. Our analyses demonstrate the utility of zebrafish scale development as a tractable system in which to elucidate mechanisms of developmental patterning and morphogenesis, and suggest a single, ancient origin of skin appendage patterning mechanisms in vertebrates.

The hetZ gene indirectly regulates heterocyst development at the level of pattern formation in Anabaena sp. strain PCC 7120

Multicellular development requires the careful orchestration of gene expression to correctly create and position specialized cells. In the filamentous cyanobacterium Anabaena sp. strain PCC 7120, nitrogen-fixing heterocysts are differentiated from vegetative cells in a reproducibly periodic and physiologically relevant pattern. While many genetic factors required for heterocyst development have been identified, the role of HetZ has remained unclear. Here, we present evidence to clarify the requirement of hetZ for heterocyst production and support a model where HetZ functions in the patterning stage of differentiation. We show that a clean, nonpolar deletion of hetZ fails to express the developmental genes hetR, patS, hetP and hetZ correctly and fails to produce heterocysts. Complementation and overexpression of hetZ in a hetP mutant revealed that hetZ was incapable of bypassing hetP, suggesting that it acts upstream of hetP. Complementation and overexpression of hetZ in a hetR mutant, however, demonstrated bypass of hetR, suggesting that it acts downstream of hetR and is capable of bypassing the need for hetR for differentiation irrespective of nitrogen status. Finally, protein-protein interactions were observed between HetZ and HetR, Alr2902 and HetZ itself. Collectively, this work suggests a regulatory role for HetZ in the patterning phase of cellular differentiation in Anabaena.

Contraction of basal filopodia controls periodic feather branching via Notch and FGF signaling

Branching morphogenesis is a general mechanism that increases the surface area of an organ. In chicken feathers, the flat epithelial sheath at the base of the follicle is transformed into periodic branches. How exactly the keratinocytes are organized into this pattern remains unclear. Here we show that in the feather follicle, the pre-branch basal keratinocytes have extensive filopodia, which contract and smooth out after branching. Manipulating the filopodia via small GTPases RhoA/Cdc42 also regulates branch formation. These basal filopodia help interpret the proximal-distal FGF gradient in the follicle. Furthermore, the topological arrangement of cell adhesion via E-Cadherin re-distribution controls the branching process. Periodic activation of Notch signaling drives the differential cell adhesion and contraction of basal filopodia, which occurs only below an FGF signaling threshold. Our results suggest a coordinated adjustment of cell shape and adhesion orchestrates feather branching, which is regulated by Notch and FGF signaling.

Keratinocytes are organised into a periodic pattern in feather branching, but how this is regulated is unclear. Here, the authors show that there is a coordinated change in cell shape and adherence, mediated by Notch, FGF signalling and Rho GTPases, which in turn regulates feather branching.

Role of the boundary in feather bud formation on one-dimensional bioengineered skin

The role of a boundary in pattern formation from a homogenous state in Turing's reaction–diffusion equations is important, particularly when the domain size is comparable to the pattern scale. Such experimental conditions may be achieved for in vitro regeneration of ectodermal appendages such as feathers, via reconstruction of embryonic single cells. This procedure can eliminate a predefined genetic map, such as the midline of chick feather bud formation, leaving uniformly distributed identical cells as a bioengineered skin. Here, the self-organizing nature of multiple feather bud formation was examined in bioengineered 1D-skin samples. Primal formation of feather buds occurred at a fixed length from the skin edge. This formation was numerically recapitulated by a standard two-component reaction-diffusion model, suggesting that the boundary effect caused this observation. The proper boundary conditions were nonstandard, either mixed Dirichlet–Neumann or partial-flux. In addition, the model implies imperfect or hindered bud formation as well as nearly equal distances between buds. In contrast, experimental observations indicated that the skin curvature, which was not included in our model, also strongly affected bud formation. Thus, bioengineered skin may provide an ideal template for modeling a self-organized process from a homogenous state. This study will examine the possible diffusion activities of activator or inhibitor molecular candidates and mechanical activities during cell aggregation, which will advance our understanding of skin appendage regeneration from pluripotent or embryonic stem cells.

Regional differences

The skin is a complex landscape containing regions in which hair follicles exhibit different types of behavior.

A wide temporal window for conjunctival papillae development ensures the formation of a complete sclerotic ring

BACKGROUND

The conjunctival papillae are epithelial thickenings of the conjunctiva that are required for the induction of underlying bones (the scleral ossicles). These transient papillae develop and become inductively active over an extended temporal period (HH 30-36, 6.5-10 dpf). While their inductive capacity was discovered in the mid-1900s, little is known about their development.

RESULTS

Through a series of timed surgical ablations followed by in situ hybridization for Bmp2, we show that the ring of conjunctival papillae is not altered if the conjunctival epithelium is ablated either prior to or shortly after papillae induction (i.e., HH 29-30, 6.5-7 dpf). A conjunctival papilla ablated at or prior to HH 34 (8 dpf), when the complete ring is present, regenerates and quickly becomes inductively active, inducing an underlying scleral condensation with only a slight delay. This regenerative capacity extends until HH 35.5, a full 36 hours beyond the normal timeline of papillae induction. As such, the period of epithelial competency for papilla induction is longer than previously identified.

CONCLUSIONS

Papilla regeneration is a mechanism that ensures the formation of a complete sclerotic ring and provides another level of redundancy for the induction of a complete sclerotic ring during the normal inductive period. Developmental Dynamics 246:381-391, 2017. © 2017 Wiley Periodicals, Inc.

Identification and evolutionary analysis of long non-coding RNAs in zebra finch

Background

Long non-coding RNAs (lncRNAs) are important in various biological processes, but very few studies on lncRNA have been conducted in birds. To identify IncRNAs expressed during feather development, we analyzed single-stranded RNA-seq (ssRNA-seq) data from the anterior and posterior dorsal regions during zebra finch (Taeniopygia guttata) embryonic development. Using published transcriptomic data, we further analyzed the evolutionary conservation of IncRNAs in birds and amniotes.

Results

A total of 1,081 lncRNAs, including 965 intergenic lncRNAs (lincRNAs), 59 intronic lncRNAs, and 57 antisense lncRNAs (lncNATs), were identified using our newly developed pipeline. These avian IncRNAs share similar characteristics with lncRNAs in mammals, such as shorter transcript length, lower exon number, lower average expression level and less sequence conservation than mRNAs. However, the proportion of lncRNAs overlapping with transposable elements in birds is much lower than that in mammals. We predicted the functions of IncRNAs based on the enriched functions of co-expressed protein-coding genes. Clusters of lncRNAs associated with natal down development were identified. The sequences and expression levels of candidate lncRNAs that shared conserved sequences among birds were validated by qPCR in both zebra finch and chicken. Finally, we identified three highly conserved lncRNAs that may be associated with natal down development.

Conclusions

Our study provides the first systematical identification of avian lncRNAs using ssRNA-seq analysis and offers a resource of embryonically expressed lncRNAs in zebra finch. We also predicted the biological function of identified lncRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3506-z) contains supplementary material, which is available to authorized users.

Building and maintaining joints by exquisite local control of cell fate

We owe the flexibility of our bodies to sophisticated articulations between bones. Establishment of these joints requires the integration of multiple tissue types: permanent cartilage that cushions the articulating bones, synovial membranes that enclose a lubricating fluid-filled cavity, and a fibrous capsule and ligaments that provide structural support. Positioning the prospective joint region involves establishment of an "interzone" region of joint progenitor cells within a nascent cartilage condensation, which is achieved through the interplay of activators and inhibitors of multiple developmental signaling pathways. Within the interzone, tight regulation of BMP and TGFβ signaling prevents the hypertrophic maturation of joint chondrocytes, in part through downstream transcriptional repressors and epigenetic modulators. Synovial cells then acquire further specializations through expression of genes that promote lubrication, as well as the formation of complex structures such as cavities and entheses. Whereas genetic investigations in mice and humans have uncovered a number of regulators of joint development and homeostasis, recent work in zebrafish offers a complementary reductionist approach toward understanding joint positioning and the regulation of chondrocyte fate at joints. The complexity of building and maintaining joints may help explain why there are still few treatments for osteoarthritis, one of the most common diseases in the human population. A major challenge will be to understand how developmental abnormalities in joint structure, as well as postnatal roles for developmental genes in joint homeostasis, contribute to birth defects and degenerative diseases of joints. WIREs Dev Biol 2017, 6:e245. doi: 10.1002/wdev.245 For further resources related to this article, please visit the WIREs website.

The "tao" of integuments

Hair follicle and sweat gland fates can be switched by morphogens at specific skin regions or developmental stages

The heterocyst regulatory protein HetP and its homologs modulate heterocyst commitment in Anabaena sp. strain PCC 7120

The commitment of differentiating cells to a specialized fate is fundamental to the correct assembly of tissues within a multicellular organism. Because commitment is often irreversible, entry into and progression through this phase of development must be tightly regulated. Under nitrogen-limiting conditions, the multicellular cyanobacterium Anabaena sp. strain PCC 7120 terminally commits ∼10% of its cells to become specialized nitrogen-fixing heterocysts. Although commitment is known to occur 9-14 h after the induction of differentiation, the factors that regulate the initiation and duration of this phase have yet to be elucidated. Here, we report the identification of four genes that share a functional domain and modulate heterocyst commitment: hetP (alr2818), asl1930, alr2902, and alr3234 Epistatic relationships between all four genes relating to commitment were revealed by deleting them individually and in combination; asl1930 and alr3234 acted most upstream to delay commitment, alr2902 acted next in the pathway to inhibit development, and hetP acted most downstream to drive commitment forward. Possible protein-protein interactions between HetP, its homologs, and the heterocyst master regulator, HetR, were assessed, and interaction partners were defined. Finally, patterns of gene expression for each homolog, as determined by promoter fusions to gfp and reverse transcription-quantitative PCR, were distinct from that of hetP in both spatiotemporal organization and regulation. We posit that a dynamic succession of protein-protein interactions modulates the timing and efficiency of the commitment phase of development and note that this work highlights the utility of a multicellular cyanobacterium as a model for the study of developmental processes.

Quorum sensing and other collective regenerative behavior in organ populations

Stem cell and microenvironment molecular interactions have been studied in detail but regenerative behavior at the organ population level has remained unexplored. Organ renewal can occur continuously or in cyclic episodes. Progenitors may be distributed as one entity or compartmentalized into multiple units. Multiple units offer advantages as each unit can be regulated differently in different body regions or physiological stages, adapting animals to their niche with flexible functional forms. Using the hair paradigm, we show how periodic patterning can convert one morphogenetic field into many hair germs, how follicles can be renewed with different cycle times and phenotypes in a region-specific manner, and how new properties, such as regenerative waves and quorum sensing, emerge to coordinate collective regenerative behavior.

A Complex Structural Variation on Chromosome 27 Leads to the Ectopic Expression of HOXB8 and the Muffs and Beard Phenotype in Chickens

Muffs and beard (Mb) is a phenotype in chickens where groups of elongated feathers gather from both sides of the face (muffs) and below the beak (beard). It is an autosomal, incomplete dominant phenotype encoded by the Muffs and beard (Mb) locus. Here we use genome-wide association (GWA) analysis, linkage analysis, Identity-by-Descent (IBD) mapping, array-CGH, genome re-sequencing and expression analysis to show that the Mb allele causing the Mb phenotype is a derived allele where a complex structural variation (SV) on GGA27 leads to an altered expression of the gene HOXB8. This Mb allele was shown to be completely associated with the Mb phenotype in nine other independent Mb chicken breeds. The Mb allele differs from the wild-type mb allele by three duplications, one in tandem and two that are translocated to that of the tandem repeat around 1.70 Mb on GGA27. The duplications contain total seven annotated genes and their expression was tested during distinct stages of Mb morphogenesis. A continuous high ectopic expression of HOXB8 was found in the facial skin of Mb chickens, strongly suggesting that HOXB8 directs this regional feather-development. In conclusion, our results provide an interesting example of how genomic structural rearrangements alter the regulation of genes leading to novel phenotypes. Further, it again illustrates the value of utilizing derived phenotypes in domestic animals to dissect the genetic basis of developmental traits, herein providing novel insights into the likely role of HOXB8 in feather development and differentiation.

Genetic variation is a key part for the study of evolution, development and differentiation. In domestic animals, many breeds display striking phenotypes that differentiate them from their wild ancestors. Several of these have been related to structural variations, including Fibromelanosis and Rose-comb in chickens, Double-muscled and Osteopetrosis in cattle, Cone degeneration in dogs, and White coat color in pigs. The feather is a type of skin appendages that exists in multiple variants on different body parts, and the derived feathering phenotypes in domestic birds are perfect resources to decipher the mechanisms regulating feather development and differentiation. Here we study the genetics of the Muffs and beard trait, a variant that alters the feather development in the facial area of chickens. We show that this phenotype is associated with a genomic structural variant that leads to an ectopic expression of HOXB8 in the facial skin during feather development. This is thus another example of how structural variants in the genome lead to novel, derived phenotypic changes in domestic animals and suggests an important role for HOXB8 in feather development.

Deciphering principles of morphogenesis from temporal and spatial patterns on the integument

BACKGROUND

How tissue patterns form in development and regeneration is a fundamental issue remaining to be fully understood. The integument often forms repetitive units in space (periodic patterning) and time (cyclic renewal), such as feathers and hairs. Integument patterns are visible and experimentally manipulatable, helping us reveal pattern formative processes. Variability is seen in regional phenotypic specificities and temporal cycling at different physiological stages.

RESULTS

Here we show some cellular/molecular bases revealed by analyzing integument patterns. (1) Localized cellular activity (proliferation, rearrangement, apoptosis, differentiation) transforms prototypic organ primordia into specific shapes. Combinatorial positioning of different localized activity zones generates diverse and complex organ forms. (2) Competitive equilibrium between activators and inhibitors regulates stem cells through cyclic quiescence and activation.

CONCLUSIONS

Dynamic interactions between stem cells and their adjacent niche regulate regenerative behavior, modulated by multi-layers of macro-environmental factors (dermis, body hormone status, and external environment). Genomics studies may reveal how positional information of localized cellular activity is stored. In vivo skin imaging and lineage tracing unveils new insights into stem cell plasticity. Principles of self-assembly obtained from the integumentary organ model can be applied to help restore damaged patterns during regenerative wound healing and for tissue engineering to rebuild tissues. Developmental Dynamics 244:905-920, 2015. © 2015 Wiley Periodicals, Inc.

Evolution and functional significance of derived sternal ossification patterns in ornithothoracine birds

The midline pattern of sternal ossification characteristic of the Cretaceous enantiornithine birds is unique among the Ornithodira, the group containing birds, nonavian dinosaurs and pterosaurs. This has been suggested to indicate that Enantiornithes is not the sister group of Ornithuromorpha, the clade that includes living birds and their close relatives, which would imply rampant convergence in many nonsternal features between enantiornithines and ornithuromorphs. However, detailed comparisons reveal greater similarity between neornithine (i.e. crown group bird) and enantiornithine modes of sternal ossification than previously recognized. Furthermore, a new subadult enantiornithine specimen demonstrates that sternal ossification followed a more typically ornithodiran pattern in basal members of the clade. This new specimen, referable to the Pengornithidae, indicates that the unique ossification pattern observed in other juvenile enantiornithines is derived within Enantiornithes. A similar but clearly distinct pattern appears to have evolved in parallel in the ornithuromorph lineage. The atypical mode of sternal ossification in some derived enantiornithines should be regarded as an autapomorphic condition rather than an indication that enantiornithines are not close relatives of ornithuromorphs. Based on what is known about molecular mechanisms for morphogenesis and the possible selective advantages, the parallel shifts to midline ossification that took place in derived enantiornithines and living neognathous birds appear to have been related to the development of a large ventral keel, which is only present in ornithuromorphs and enantiornithines. Midline ossification can serve to medially reinforce the sternum at a relatively early ontogenetic stage, which would have been especially beneficial during the protracted development of the superprecocial Cretaceous enantiornithines.

Distinct developmental genetic mechanisms underlie convergently evolved tooth gain in sticklebacks

Teeth are a classic model system of organogenesis, as repeated and reciprocal epithelial and mesenchymal interactions pattern placode formation and outgrowth. Less is known about the developmental and genetic bases of tooth formation and replacement in polyphyodonts, which are vertebrates with continual tooth replacement. Here, we leverage natural variation in the threespine stickleback fish Gasterosteus aculeatus to investigate the genetic basis of tooth development and replacement. We find that two derived freshwater stickleback populations have both convergently evolved more ventral pharyngeal teeth through heritable genetic changes. In both populations, evolved tooth gain manifests late in development. Using pulse-chase vital dye labeling to mark newly forming teeth in adult fish, we find that both high-toothed freshwater populations have accelerated tooth replacement rates relative to low-toothed ancestral marine fish. Despite the similar evolved phenotype of more teeth and an accelerated adult replacement rate, the timing of tooth number divergence and the spatial patterns of newly formed adult teeth are different in the two populations, suggesting distinct developmental mechanisms. Using genome-wide linkage mapping in marine-freshwater F2 genetic crosses, we find that the genetic basis of evolved tooth gain in the two freshwater populations is largely distinct. Together, our results support a model whereby increased tooth number and an accelerated tooth replacement rate have evolved convergently in two independently derived freshwater stickleback populations using largely distinct developmental and genetic mechanisms.

Development, regeneration, and evolution of feathers

The feather is a complex ectodermal organ with hierarchical branching patterns. It provides functions in endothermy, communication, and flight. Studies of feather growth, cycling, and health are of fundamental importance to avian biology and poultry science. In addition, feathers are an excellent model for morphogenesis studies because of their accessibility, and their distinct patterns can be used to assay the roles of specific molecular pathways. Here we review the progress in aspects of development, regeneration, and evolution during the past three decades. We cover the development of feather buds in chicken embryos, regenerative cycling of feather follicle stem cells, formation of barb branching patterns, emergence of intrafeather pigmentation patterns, interplay of hormones and feather growth, and the genetic identification of several feather variants. The discovery of feathered dinosaurs redefines the relationship between feathers and birds. Inspiration from biomaterials and flight research further fuels biomimetic potential of feathers as a multidisciplinary research focal point.

The origin and loss of periodic patterning in the turtle shell

The origin of the turtle shell over 200 million years ago greatly modified the amniote body plan, and the morphological plasticity of the shell has promoted the adaptive radiation of turtles. The shell, comprising a dorsal carapace and a ventral plastron, is a layered structure formed by basal endochondral axial skeletal elements (ribs, vertebrae) and plates of bone, which are overlain by keratinous ectodermal scutes. Studies of turtle development have mostly focused on the bones of the shell; however, the genetic regulation of the epidermal scutes has not been investigated. Here, we show that scutes develop from an array of patterned placodes and that these placodes are absent from a soft-shelled turtle in which scutes were lost secondarily. Experimentally inhibiting Shh, Bmp or Fgf signaling results in the disruption of the placodal pattern. Finally, a computational model is used to show how two coupled reaction-diffusion systems reproduce both natural and abnormal variation in turtle scutes. Taken together, these placodal signaling centers are likely to represent developmental modules that are responsible for the evolution of scutes in turtles, and the regulation of these centers has allowed for the diversification of the turtle shell.

Parallel developmental genetic features underlie stickleback gill raker evolution

Background

Convergent evolution, the repeated evolution of similar phenotypes in independent lineages, provides natural replicates to study mechanisms of evolution. Cases of convergent evolution might have the same underlying developmental and genetic bases, implying that some evolutionary trajectories might be predictable. In a classic example of convergent evolution, most freshwater populations of threespine stickleback fish have independently evolved a reduction of gill raker number to adapt to novel diets. Gill rakers are a segmentally reiterated set of dermal bones important for fish feeding. A previous large quantitative trait locus (QTL) mapping study using a marine × freshwater F2 cross identified QTL on chromosomes 4 and 20 with large effects on evolved gill raker reduction.

Results

By examining skeletal morphology in adult and developing sticklebacks, we find heritable marine/freshwater differences in gill raker number and spacing that are specified early in development. Using the expression of the Ectodysplasin receptor (Edar) gene as a marker of raker primordia, we find that the differences are present before the budding of gill rakers occurs, suggesting an early change to a lateral inhibition process controlling raker primordia spacing. Through linkage mapping in F2 fish from crosses with three independently derived freshwater populations, we find in all three crosses QTL overlapping both previously identified QTL on chromosomes 4 and 20 that control raker number. These two QTL affect the early spacing of gill raker buds.

Conclusions

Collectively, these data demonstrate that parallel developmental genetic features underlie the convergent evolution of gill raker reduction in freshwater sticklebacks, suggesting that even highly polygenic adaptive traits can have a predictable developmental genetic basis.

Does osteoderm growth follow energy minimization principles?

Although the growth and development of tissues and organs of extinct species cannot be directly observed, their fossils can record and preserve evidence of these mechanisms. It is generally accepted that bone architecture is the result of genetically based biomechanical constraints, but what about osteoderms? In this article, the influence of physical constraints on cranial osteoderms growth is assessed. Comparisons among lepidosaurs, synapsids, and archosaurs are performed; according to these analyses, lepidosaur osteoderms growth is predicted to be less energy demanding than that of synapsids and archosaurs. Obtained results also show that, from an energetic viewpoint, ankylosaurid osteoderms growth resembles more that of mammals than the one of reptilians, adding evidence to debate whether dinosaurs were hot or cold blooded.

Macroenvironmental regulation of hair cycling and collective regenerative behavior

The hair follicle (HF) regeneration paradigm provides a unique opportunity for studying the collective behavior of stem cells in living animals. Activation of HF stem cells depends on the core inhibitory BMP and activating WNT signals operating within the HF microenvironment. Additionally, HFs receive multilayered signaling inputs from the extrafollicular macroenvironment, which includes dermis, adipocytes, neighboring HFs, hormones, and external stimuli. These activators/inhibitors are integrated across multiple stem-cell niches to produce dynamic hair growth patterns. Because of their pigmentation, these patterns can be easily studied on live shaved animals. Comparing to autonomous regeneration of one HF, populations of HFs display coupled decision making, allowing for more robust and adaptable regenerative behavior to occur collectively. The generic cellular automata model used to simulate coordinated HF cycling here can be extended to study population-level behavior of other complex biological systems made of cycling elements.