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.

Adaptive patterning of vascular network during avian skin development: Mesenchymal plasticity and dermal vasculogenesis.

A vasculature network supplies blood to feather buds in the developing skin. Does the vasculature network during early skin development form by sequential sprouting from the central vasculature or does local vasculogenesis occur first that then connect with the central vascular tree? Using transgenic Japanese quail Tg(TIE1p.H2B-eYFP), we observe that vascular progenitor cells appear after feather primordia formation. The vasculature then radiates out from each bud and connects with primordial vessels from neighboring buds. Later they connect with the central vasculature. Epithelial-mesenchymal recombination shows local vasculature is patterned by the epithelium, which expresses FGF2 and VEGF. Perturbing noggin expression leads to abnormal vascularization. To study endothelial origin, we compare transcriptomes of TIE1p.H2B-eYFP+ cells collected from the skin and aorta. Endothelial cells from the skin more closely resemble skin dermal cells than those from the aorta. The results show developing chicken skin vasculature is assembled by (1) physiological vasculogenesis from the peripheral tissue, and (2) subsequently connects with the central vasculature. The work implies mesenchymal plasticity and convergent differentiation play significant roles in development, and such processes may be re-activated during adult regeneration. SUMMARY STATEMENT: We show the vasculature network in the chicken skin is assembled using existing feather buds as the template, and endothelia are derived from local bud dermis and central vasculature.

Epidermal-dermal coupled spheroids are important for tissue pattern regeneration in reconstituted skin explant cultures.

Tissue patterning is critical for the development and regeneration of organs. To advance the use of engineered reconstituted skin organs, we study cardinal features important for tissue patterning and hair regeneration. We find they spontaneously form spheroid configurations, with polarized epidermal cells coupled with dermal cells through a newly formed basement membrane. Functionally, the spheroid becomes competent morphogenetic units (CMU) that promote regeneration of tissue patterns. The emergence of new cell types and molecular interactions during CMU formation was analyzed using scRNA-sequencing. Surprisingly, in newborn skin explants, IFNr signaling can induce apical-basal polarity in epidermal cell aggregates. Dermal-Tgfb induces basement membrane formation. Meanwhile, VEGF signaling mediates dermal cell attachment to the epidermal cyst shell, thus forming a CMU. Adult mouse and human fetal scalp cells fail to form a CMU but can be restored by adding IFNr or VEGF to achieve hair regeneration. We find different multi-cellular configurations and molecular pathways are used to achieve morphogenetic competence in developing skin, wound-induced hair neogenesis, and reconstituted explant cultures. Thus, multiple paths can be used to achieve tissue patterning. These insights encourage more studies of "in vitro morphogenesis" which may provide novel strategies to enhance regeneration.

Using Avian Skin Explants to Study Tissue Patterning and Organogenesis.

The developing avian skin during embryogenesis is a unique model that can provide valuable insights into tissue patterning. Here three variations on skin explant cultures to examine different aspects of skin development are described. First, ex vivo organ cultures and manipulations offer researchers opportunities to observe and study the development of feather buds directly. Skin explant culture can grow for 7 days enabling direct analysis of cellular behavior and 4D imaging at intervals during this growth period. This also allows for physical and molecular manipulations of culture conditions to visualize tissue response. For example, growth factor-coated beads can be applied locally to induce changes in feather patterning in a limited area. Alternatively, viral transduction can be delivered globally in the culture media to up or downregulate gene expression. Second, the skin recombination protocol allows researchers to investigate tissue interactions between the epidermis and mesenchyme that are derived from different skin regions, different life stages, or different species. This affords an opportunity to test the time window in which the epithelium is competent to respond to signals and its ability to form different skin appendages in response to signals from different mesenchymal sources. Third, skin reconstitution using dissociated dermal cells overlaid with intact epithelium resets skin development and enables the study of the initial processes of periodic patterning. This approach also enhances our ability to manipulate gene expression among the dissociated cells before creating the reconstituted skin explant. This paper provides the three culture protocols and exemplary experiments to demonstrate their utility.

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.

Regional Specific Differentiation of Integumentary Organs: Regulation of Gene Clusters within the Avian Epidermal Differentiation Complex and Impacts of SATB2 Overexpression.

The epidermal differentiation complex (EDC) encodes a group of unique proteins expressed in late epidermal differentiation. The EDC gave integuments new physicochemical properties and is critical in evolution. Recently, we showed ß-keratins, members of the EDC, undergo gene cluster switching with overexpression of SATB2 (Special AT-rich binding protein-2), considered a chromatin regulator. We wondered whether this unique regulatory mechanism is specific to ß-keratins or may be derived from and common to EDC members. Here we explore (1) the systematic expression patterns of non-ß-keratin EDC genes and their preferential expression in different skin appendages during development, (2) whether the expression of non-ß-keratin EDC sub-clusters are also regulated in clusters by SATB2. We analyzed bulk RNA-seq and ChIP-seq data and also evaluated the disrupted expression patterns caused by overexpressing SATB2. The results show that the expression of whole EDDA and EDQM sub-clusters are possibly mediated by enhancers in E14-feathers. Overexpressing SATB2 down-regulates the enriched EDCRP sub-cluster in feathers and the EDCH sub-cluster in beaks. These results reveal the potential of complex epigenetic regulation activities within the avian EDC, implying transcriptional regulation of EDC members acting at the gene and/or gene cluster level in a temporal and skin regional-specific fashion, which may contribute to the evolution of diverse avian integuments.

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.

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.

Integrating Bioelectrical Currents and Ca2+ Signaling with Biochemical Signaling in Development and Pathogenesis.

Roles of bioelectrical signals are increasingly recognized in excitable and nonexcitable non-neural tissues. Diverse ion-selective channels, pumps, and gap junctions participate in bioelectrical signaling, including those transporting calcium ions (Ca2+). Ca2+ is the most versatile transported ion, because it serves as an electrical charge carrier and a biochemical regulator for multiple molecular binding, enzyme, and transcription activities. We aspire to learn how bioelectrical signals crosstalk to biochemical/biomechanical signals. In this study, we review four recent studies showing how bioelectrical currents and Ca2+ signaling affect collective dermal cell migration during feather bud elongation, affect chondrogenic differentiation in limb development, couple with mechanical tension in aligning gut smooth muscle, and affect mitochondrial function and skeletal muscle atrophy. We observe bioelectrical signals involved in several developmental and pathological conditions in chickens and mice at multiple spatial scales: cellular, cellular collective, and subcellular. These examples inspire novel concept and approaches for future basic and translational studies.

The Making of a Flight Feather: Bio-architectural Principles and Adaptation.

The evolution of flight in feathered dinosaurs and early birds over millions of years required flight feathers whose architecture features hierarchical branches. While barb-based feather forms were investigated, feather shafts and vanes are understudied. Here, we take a multi-disciplinary approach to study their molecular control and bio-architectural organizations. In rachidial ridges, epidermal progenitors generate cortex and medullary keratinocytes, guided by Bmp and transforming growth factor ß (TGF-ß) signaling that convert rachides into adaptable bilayer composite beams. In barb ridges, epidermal progenitors generate cylindrical, plate-, or hooklet-shaped barbule cells that form fluffy branches or pennaceous vanes, mediated by asymmetric cell junction and keratin expression. Transcriptome analyses and functional studies show anterior-posterior Wnt2b signaling within the dermal papilla controls barbule cell fates with spatiotemporal collinearity. Quantitative bio-physical analyses of feathers from birds with different flight characteristics and feathers in Burmese amber reveal how multi-dimensional functionality can be achieved and may inspire future composite material designs. VIDEO ABSTRACT.

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.

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.

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 Ca2+ channels (VGCCs), Connexin-43 based gap junctions, and Ca2+ release activated Ca2+ (CRAC) channels. 4-Dimensional Ca2+ imaging reveals that the Sonic hedgehog-responsive mesenchymal cells display synchronized Ca2+ oscillations, which expand progressively in area during feather elongation. Inhibiting VGCCs, gap junctions, or Sonic hedgehog signaling alters the mesenchymal Ca2+ landscape, cell movement patterns and feather bud elongation. Ca2+ 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 Ca2+ profile among cells, thereby coordinating cell movement patterns.

Transcriptome analyses of reprogrammed feather / scale chimeric explants revealed co-expressed epithelial gene networks during organ specification.

BACKGROUND: The molecular mechanism controlling regional specific skin appendage phenotypes is a fundamental question that remains unresolved. We recently identified feather and scale primordium associated genes and with functional studies, proposed five major modules are involved in scale-to-feather conversion and their integration is essential to form today's feathers. Yet, how the molecular networks are wired and integrated at the genomic level is still unknown. RESULTS: Here, we combine classical recombination experiments and systems biology technology to explore the molecular mechanism controlling cell fate specification. In the chimeric explant, dermal fate is more stable, while epidermal fate is reprogrammed to be similar to the original appendage type of the mesenchyme. We analyze transcriptome changes in both scale-to-feather and feather-to-scale transition in the epidermis. We found a highly interconnected regulatory gene network controlling skin appendage types. These gene networks are organized around two molecular hubs, ß-catenin and retinoic acid (RA), which can bind to regulatory elements controlling downstream gene expression, leading to scale or feather fates. ATAC sequencing analyses revealed about 1000 altered widely distributed chromatin open sites. We find that perturbation of a key gene alters the expression of many other co-expressed genes in the same module. CONCLUSIONS: Our findings suggest that these feather / scale fate specification genes form an interconnected network and rewiring of the gene network can lead to changes of appendage phenotypes, acting similarly to endogenous reprogramming at the tissue level. This work shows that key hub molecules, ß-catenin and retinoic acid, regulate scale / feather fate specification gene networks, opening up new possibilities to understand the switches controlling organ phenotypes in a two component (epithelial and mesenchyme) system.

Comparative genomics and transcriptomics of Chrysolophus provide insights into the evolution of complex plumage coloration.

Background: As one of the most recognizable characteristics in birds, plumage color has a high impact on understanding the evolution and mechanisms of coloration. Feather and skin are ideal tissues to explore the genomics and complexity of color patterns in vertebrates. Two species of the genus Chrysolophus, golden pheasant (Chrysolophus pictus) and Lady Amherst's pheasant (Chrysolophus amherstiae), exhibit brilliant colors in their plumage, but with extreme phenotypic differences, making these two species great models to investigate plumage coloration mechanisms in birds. Results: We sequenced and assembled a genome of golden pheasant with high coverage and annotated 15,552 protein-coding genes. The genome of Lady Amherst's pheasant is sequenced with low coverage. Based on the feather pigment identification, a series of genomic and transcriptomic comparisons were conducted to investigate the complex features of plumage coloration. By identifying the lineage-specific sequence variations in Chrysolophus and golden pheasant against different backgrounds, we found that four melanogenesis biosynthesis genes and some lipid-related genes might be candidate genomic factors for the evolution of melanin and carotenoid pigmentation, respectively. In addition, a study among 47 birds showed some candidate genes related to carotenoid coloration in a broad range of birds. The transcriptome data further reveal important regulators of the two colorations, particularly one splicing transcript of the microphthalmia-associated transcription factor gene for pheomelanin synthesis. Conclusions: Analysis of the golden pheasant and its sister pheasant genomes, as well as comparison with other avian genomes, are helpful to reveal the underlying regulation of their plumage coloration. The present study provides important genomic information and insights for further studies of avian plumage evolution and diversity.

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.

Msx2 Supports Epidermal Competency during Wound-Induced Hair Follicle Neogenesis.

Cutaneous wounds in adult mammals typically heal by scarring. However, large full-thickness wounds undergo wound-induced hair follicle neogenesis (WIHN), a form of regeneration. Here, we show that WIHN requires transient expression of epidermal Msx2 in two phases: the wound margin early and the wound center late. Msx2 expression is present in the migrating epithelium during early wound healing and then presents in the epithelium and mesenchyme later in the wound center. WIHN is abrogated in germline and epithelial-specific Msx2 mutant mice. Unlike the full-length Msx2 promoter, a minimal Msx2 promoter fails activation in the wound center, suggesting complex regulation of Msx2 expression. The Msx2 promoter binding sites include Tcf/Lef, Jun/Creb, Pax3, and three SMAD sites. However, basal epithelial-induced BMP suppression by noggin overexpression did not affect WIHN. We propose that Msx2 signaling is required for the epidermis to acquire spatiotemporal competence during WIHN. Topologically, hair regeneration dominates in the wound center, coinciding with late Msx2 expression. Together, these results suggest that intrinsic Msx2 expression supports epithelial competency during hair follicle neogenesis. This work provides insight into endogenous mechanisms modulating competency of adult epidermal progenitors for mammalian ectodermal appendage neogenesis, and offers the target Msx2 for future regeneration-promoting therapies.

Systems Biology Analyses in Chicken: Workflow for Transcriptome and ChIP-Seq Analyses Using the Chicken Skin Paradigm.

With advances in molecular biology, various biological phenomena can now be explored at higher resolution using mRNA sequencing (RNA-Seq) and chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq), two powerful high-throughput next-generation sequencing (NGS) technologies. While methods are used widely in mouse, human, etc., less information is available in other animals, such as the chicken. Here we assemble a workflow of the RNA-Seq and ChIP-Seq analyses for the chicken studies using chicken skin appendage tissue as an example. We present guidelines for RNA-Seq quality control, alignment, quantification, normalization, and differentially expressed gene analysis. In the meantime, we outline a bioinformatics pipeline for ChIP-Seq quality control, alignment, peak calling, super-enhancer identification, and differential enrichment analysis.

MicroCT Imaging on Living Alligator Teeth Reveals Natural Tooth Cycling.

To study tooth cycling in polyphyodont animals, we chose to work on alligators. Alligators have teeth in three phases of development at each tooth location. This assembly of three teeth is called a tooth family unit. As part of the study, in order to study tooth cycling in alligators, we wanted to know the configuration of the tooth family unit in every tooth position. From the surface of the mouth, this is difficult to assess. Therefore, we decided to use MicroCT which can image X-ray dense materials providing a three-dimensional view. MicroCT provided us with valuable information for this study. The method described below can be applied to study tooth cycling in other vertebrate species.

The Modulatable Stem Cell Niche: Tissue Interactions during Hair and Feather Follicle Regeneration.

Hair and feathers are unique because (1) their stem cells are contained within a follicle structure, (2) they undergo cyclic regeneration repetitively throughout life, (3) regeneration occurs physiologically in healthy individuals and (4) regeneration is also induced in response to injury. Precise control of this cyclic regeneration process is essential for maintaining the homeostasis of living organisms. While stem cells are regulated by the intra-follicle-adjacent micro-environmental niche, this niche is also modulated dynamically by extra-follicular macro-environmental signals, allowing stem cells to adapt to a larger changing environment and physiological needs. Here we review several examples of macro-environments that communicate with the follicles: intradermal adipose tissue, innate immune system, sex hormones, aging, circadian rhythm and seasonal rhythms. Related diseases are also discussed. Unveiling the mechanisms of how stem cell niches are modulated provides clues for regenerative medicine. Given that stem cells are hard to manipulate, focusing translational therapeutic applications at the environments appears to be a more practical approach.