The developmental basis of fingerprint pattern formation and variation.

Fingerprints are complex and individually unique patterns in the skin. Established prenatally, the molecular and cellular mechanisms that guide fingerprint ridge formation and their intricate arrangements are unknown. Here we show that fingerprint ridges are epithelial structures that undergo a truncated hair follicle developmental program and fail to recruit a mesenchymal condensate. Their spatial pattern is established by a Turing reaction-diffusion system, based on signaling between EDAR, WNT, and antagonistic BMP pathways. These signals resolve epithelial growth into bands of focalized proliferation under a precociously differentiated suprabasal layer. Ridge formation occurs as a set of waves spreading from variable initiation sites defined by the local signaling environments and anatomical intricacies of the digit, with the propagation and meeting of these waves determining the type of pattern that forms. Relying on a dynamic patterning system triggered at spatially distinct sites generates the characteristic types and unending variation of human fingerprint patterns.

Limb development genes underlie variation in human fingerprint patterns.

Fingerprints are of long-standing practical and cultural interest, but little is known about the mechanisms that underlie their variation. Using genome-wide scans in Han Chinese cohorts, we identified 18 loci associated with fingerprint type across the digits, including a genetic basis for the long-recognized "pattern-block" correlations among the middle three digits. In particular, we identified a variant near EVI1 that alters regulatory activity and established a role for EVI1 in dermatoglyph patterning in mice. Dynamic EVI1 expression during human development supports its role in shaping the limbs and digits, rather than influencing skin patterning directly. Trans-ethnic meta-analysis identified 43 fingerprint-associated loci, with nearby genes being strongly enriched for general limb development pathways. We also found that fingerprint patterns were genetically correlated with hand proportions. Taken together, these findings support the key role of limb development genes in influencing the outcome of fingerprint patterning.

Rapid mechanosensitive migration and dispersal of newly divided mesenchymal cells aid their recruitment into dermal condensates.

Embryonic mesenchymal cells are dispersed within an extracellular matrix but can coalesce to form condensates with key developmental roles. Cells within condensates undergo fate and morphological changes and induce cell fate changes in nearby epithelia to produce structures including hair follicles, feathers, or intestinal villi. Here, by imaging mouse and chicken embryonic skin, we find that mesenchymal cells undergo much of their dispersal in early interphase, in a stereotyped process of displacement driven by 3 hours of rapid and persistent migration followed by a long period of low motility. The cell division plane and the elevated migration speed and persistence of newly born mesenchymal cells are mechanosensitive, aligning with tissue tension, and are reliant on active WNT secretion. This behaviour disperses mesenchymal cells and allows daughters of recent divisions to travel long distances to enter dermal condensates, demonstrating an unanticipated effect of cell cycle subphase on core mesenchymal behaviour.

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

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

Feather arrays are patterned by interacting signalling and cell density waves.

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

The Enrichment of Specific Hair Follicle-Associated Cell Populations in Plucked Hairs Offers an Opportunity to Study Gene Expression Underlying Hair Traits.

Gene expression differences can assist in characterizing important underlying genetic mechanisms between different phenotypic traits. However, when population-dense tissues are studied, the signals from scarce populations are diluted. Therefore, appropriately choosing a sample collection method that enriches a particular type of effector cells might yield more specific results. To address this issue, we performed a polyA-selected RNA-seq experiment of domestic horse (Equus ferus caballus) plucked-hair samples and skin biopsies. Then, we layered the horse gene abundance results against cell type-specific marker genes generated from a scRNA-seq supported with spatial mapping of laboratory mouse (Mus musculus) skin to identify the captured populations. The hair-plucking and skin-biopsy sample-collection methods yielded comparable quality and quantity of RNA-seq results. Keratin-related genes, such as KRT84 and KRT75, were among the genes that showed higher abundance in plucked hairs, while genes involved in cellular processes and enzymatic activities, such as MGST1, had higher abundance in skin biopsies. We found an enrichment of hair-follicle keratinocytes in plucked hairs, but detected an enrichment of other populations, including epidermis keratinocytes, in skin biopsies. In mammalian models, biopsies are often the method of choice for a plethora of gene expression studies and to our knowledge, this is a novel study that compares the cell-type enrichment between the non-invasive hair-plucking and the invasive skin-biopsy sample-collection methods. Here, we show that the non-invasive and ethically uncontroversial plucked-hair method is recommended depending on the research question. In conclusion, our study will allow downstream -omics approaches to better understand integumentary conditions in both health and disease in horses as well as other mammals.

A genome-wide association analysis for body weight at 35 days measured on 137,343 broiler chickens.

BACKGROUND: Body weight (BW) is an economically important trait in the broiler (meat-type chickens) industry. Under the assumption of polygenicity, a "large" number of genes with "small" effects is expected to control BW. To detect such effects, a large sample size is required in genome-wide association studies (GWAS). Our objective was to conduct a GWAS for BW measured at 35 days of age with a large sample size. METHODS: The GWAS included 137,343 broilers spanning 15 pedigree generations and 392,295 imputed single nucleotide polymorphisms (SNPs). A false discovery rate of 1% was adopted to account for multiple testing when declaring significant SNPs. A Bayesian ridge regression model was implemented, using AlphaBayes, to estimate the contribution to the total genetic variance of each region harbouring significant SNPs (1 Mb up/downstream) and the combined regions harbouring non-significant SNPs. RESULTS: GWAS revealed 25 genomic regions harbouring 96 significant SNPs on 13 Gallus gallus autosomes (GGA1 to 4, 8, 10 to 15, 19 and 27), with the strongest associations on GGA4 at 65.67-66.31 Mb (Galgal4 assembly). The association of these regions points to several strong candidate genes including: (i) growth factors (GGA1, 4, 8, 13 and 14); (ii) leptin receptor overlapping transcript (LEPROT)/leptin receptor (LEPR) locus (GGA8), and the STAT3/STAT5B locus (GGA27), in connection with the JAK/STAT signalling pathway; (iii) T-box gene (TBX3/TBX5) on GGA15 and CHST11 (GGA1), which are both related to heart/skeleton development); and (iv) PLAG1 (GGA2). Combined together, these 25 genomic regions explained ~ 30% of the total genetic variance. The region harbouring significant SNPs that explained the largest portion of the total genetic variance (4.37%) was on GGA4 (~ 65.67-66.31 Mb). CONCLUSIONS: To the best of our knowledge, this is the largest GWAS that has been conducted for BW in chicken to date. In spite of the identified regions, which showed a strong association with BW, the high proportion of genetic variance attributed to regions harbouring non-significant SNPs supports the hypothesis that the genetic architecture of BW35 is polygenic and complex. Our results also suggest that a large sample size will be required for future GWAS of BW35.

Systems for intricate patterning of the vertebrate anatomy.

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

Isolating Patterns in Open Reaction-Diffusion Systems.

Realistic examples of reaction-diffusion phenomena governing spatial and spatiotemporal pattern formation are rarely isolated systems, either chemically or thermodynamically. However, even formulations of 'open' reaction-diffusion systems often neglect the role of domain boundaries. Most idealizations of closed reaction-diffusion systems employ no-flux boundary conditions, and often patterns will form up to, or along, these boundaries. Motivated by boundaries of patterning fields related to the emergence of spatial form in embryonic development, we propose a set of mixed boundary conditions for a two-species reaction-diffusion system which forms inhomogeneous solutions away from the boundary of the domain for a variety of different reaction kinetics, with a prescribed uniform state near the boundary. We show that these boundary conditions can be derived from a larger heterogeneous field, indicating that these conditions can arise naturally if cell signalling or other properties of the medium vary in space. We explain the basic mechanisms behind this pattern localization and demonstrate that it can capture a large range of localized patterning in one, two, and three dimensions and that this framework can be applied to systems involving more than two species. Furthermore, the boundary conditions proposed lead to more symmetrical patterns on the interior of the domain and plausibly capture more realistic boundaries in developmental systems. Finally, we show that these isolated patterns are more robust to fluctuations in initial conditions and that they allow intriguing possibilities of pattern selection via geometry, distinct from known selection mechanisms.

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.

Ectodysplasin signalling deficiency in mouse models of hypohidrotic ectodermal dysplasia leads to middle ear and nasal pathology.

Hypohidrotic ectodermal dysplasia (HED) results from mutation of the EDA, EDAR or EDARADD genes and is characterized by reduced or absent eccrine sweat glands, hair follicles and teeth, and defective formation of salivary, mammary and craniofacial glands. Mouse models with HED also carry Eda, Edar or Edaradd mutations and have defects that map to the same structures. Patients with HED have ear, nose and throat disease, but this has not been investigated in mice bearing comparable genetic mutations. We report that otitis media, rhinitis and nasopharyngitis occur at high frequency in Eda and Edar mutant mice and explore the pathogenic mechanisms related to glandular function, microbial and immune parameters in these lines. Nasopharynx auditory tube glands fail to develop in HED mutant mice and the functional implications include loss of lysozyme secretion, reduced mucociliary clearance and overgrowth of nasal commensal bacteria accompanied by neutrophil exudation. Heavy nasopharynx foreign body load and loss of gland protection alters the auditory tube gating function and the auditory tubes can become pathologically dilated. Accumulation of large foreign body particles in the bulla stimulates granuloma formation. Analysis of immune cell populations and myeloid cell function shows no evidence of overt immune deficiency in HED mutant mice. Our findings using HED mutant mice as a model for the human condition support the idea that ear and nose pathology in HED patients arises as a result of nasal and nasopharyngeal gland deficits, reduced mucociliary clearance and impaired auditory tube gating function underlies the pathological sequelae in the bulla.

A chemotaxis model of feather primordia pattern formation during avian development.

The orderly formation of the avian feather array is a classic example of periodic pattern formation during embryonic development. Various mathematical models have been developed to describe this process, including Turing/activator-inhibitor type reaction-diffusion systems and chemotaxis/mechanical-based models based on cell movement and tissue interactions. In this paper we formulate a mathematical model founded on experimental findings, a set of interactions between the key cellular (dermal and epidermal cell populations) and molecular (fibroblast growth factor, FGF, and bone morphogenetic protein, BMP) players and a medially progressing priming wave that acts as the trigger to initiate patterning. Linear stability analysis is used to show that FGF-mediated chemotaxis of dermal cells is the crucial driver of pattern formation, while perturbations in the form of ubiquitous high BMP expression suppress patterning, consistent with experiments. Numerical simulations demonstrate the capacity of the model to pattern the skin in a spatial-temporal manner analogous to avian feather development. Further, experimental perturbations in the form of bead-displacement experiments are recapitulated and predictions are proposed in the form of blocking mesenchymal cell proliferation.

Derivation of marker gene signatures from human skin and their use in the interpretation of the transcriptional changes associated with dermatological disorders.

Numerous studies have explored the altered transcriptional landscape associated with skin diseases to understand the nature of these disorders. However, data interpretation represents a significant challenge due to a lack of good maker sets for many of the specialized cell types that make up this tissue, whose composition may fundamentally alter during disease. Here we have sought to derive expression signatures that define the various cell types and structures that make up human skin, and demonstrate how they can be used to aid the interpretation of transcriptomic data derived from this organ. Two large normal skin transcriptomic datasets were identified, one RNA-seq (n = 578), the other microarray (n = 165), quality controlled and subjected separately to network-based analyses to identify clusters of robustly co-expressed genes. The biological significance of these clusters was then assigned using a combination of bioinformatics analyses, literature, and expert review. After cross comparison between analyses, 20 gene signatures were defined. These included expression signatures for hair follicles, glands (sebaceous, sweat, apocrine), keratinocytes, melanocytes, endothelia, muscle, adipocytes, immune cells, and a number of pathway systems. Collectively, we have named this resource SkinSig. SkinSig was then used in the analysis of transcriptomic datasets for 18 skin conditions, providing in-context interpretation of these data. For instance, conventional analysis has shown there to be a decrease in keratinization and fatty metabolism with age; we more accurately define these changes to be due to loss of hair follicles and sebaceous glands. SkinSig also highlighted the over-/under-representation of various cell types in skin diseases, reflecting an influx in immune cells in inflammatory disorders and a relative reduction in other cell types. Overall, our analyses demonstrate the value of this new resource in defining the functional profile of skin cell types and appendages, and in improving the interpretation of disease data. © 2016 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

FGF and EDA pathways control initiation and branching of distinct subsets of developing nasal glands.

Hypertrophy, hyperplasia and altered mucus secretion from the respiratory submucosal glands (SMG) are characteristics of airway diseases such as cystic fibrosis, asthma and chronic bronchitis. More commonly, hyper-secretion of the nasal SMGs contributes to allergic rhinitis and upper airway infection. Considering the role of these glands in disease states, there is a significant dearth in understanding the molecular signals that regulate SMG development and patterning. Due to the imperative role of FGF signalling during the development of other branched structures, we investigated the role of Fgf10 during initiation and branching morphogenesis of murine nasal SMGs. Fgf10 is expressed in the mesenchyme around developing SMGs while expression of its receptor Fgfr2 is seen within glandular epithelial cells. In the Fgf10 null embryo, Steno's gland and the maxillary sinus gland were completely absent while other neighbouring nasal glands showed normal duct elongation but defective branching. Interestingly, the medial nasal glands were present in Fgf10 homozygotes but missing in Fgfr2b mutants, with expression of Fgf7 specifically expressed around these developing glands, indicating that Fgf7 might compensate for loss of Fgf10 in this group of glands. Intriguingly the lateral nasal glands were only mildly affected by loss of FGF signalling, while these glands were missing in Eda mutant mice, where the Steno's and maxillary sinus gland developed as normal. This analysis reveals that regulation of nasal gland development is complex with different subsets of glands being regulated by different signalling pathways. This analysis helps shed light on the nasal gland defects observed in patients with hypohidrotic ectodermal dysplasia (HED) (defect EDA pathway) and LADD syndrome (defect FGFR2b pathway).

Hair Follicle Bulge Stem Cells Appear Dispensable for the Acute Phase of Wound Re-epithelialization.

The cutaneous healing response has evolved to occur rapidly, in order to minimize infection and to re-establish epithelial homeostasis. Rapid healing is achieved through complex coordination of multiple cell types, which importantly includes specific cell populations within the hair follicle (HF). Under physiological conditions, the epithelial compartments of HF and interfollicular epidermis remain discrete, with K15(+ve) bulge stem cells contributing progeny for HF reconstruction during the hair cycle and as a basis for hair shaft production during anagen. Only upon wounding do HF cells migrate from the follicle to contribute to the neo-epidermis. However, the identity of the first-responding cells, and in particular whether this process involves a direct contribution of K15(+ve) bulge cells to the early stage of epidermal wound repair remains unclear. Here we demonstrate that epidermal injury in murine skin does not induce bulge activation during early epidermal wound repair. Specifically, bulge cells of uninjured HFs neither proliferate nor appear to migrate out of the bulge niche upon epidermal wounding. In support of these observations, Diphtheria toxin-mediated partial ablation of K15(+ve) bulge cells fails to delay wound healing. Our data suggest that bulge cells only respond to epidermal wounding during later stages of repair. We discuss that this response may have evolved as a protective safeguarding mechanism against bulge stem cell exhaust and tumorigenesis. Stem Cells 2016;34:1377-1385.

Regionalisation of the skin.

The skin displays marked anatomical variation in thickness, colour and in the appendages that it carries. These regional distinctions arise in the embryo, likely founded on a combinatorial positional code of transcription factor expression. Throughout adult life, the skin's distinct anatomy is maintained through both cell autonomous epigenetic processes and by mesenchymal-epithelial induction. Despite the readily apparent anatomical differences in skin characteristics across the body, several fundamental questions regarding how such regional differences first arise and then persist are unresolved. However, it is clear that the skin's positional code is at the molecular level far more detailed than that discernible at the phenotypic level. This provides a latent reservoir of anatomical complexity ready to surface if perturbed by mutation, hormonal changes, ageing or experiment.

Generation and characterization of function-blocking anti-ectodysplasin A (EDA) monoclonal antibodies that induce ectodermal dysplasia.

Development of ectodermal appendages, such as hair, teeth, sweat glands, sebaceous glands, and mammary glands, requires the action of the TNF family ligand ectodysplasin A (EDA). Mutations of the X-linked EDA gene cause reduction or absence of many ectodermal appendages and have been identified as a cause of ectodermal dysplasia in humans, mice, dogs, and cattle. We have generated blocking antibodies, raised in Eda-deficient mice, against the conserved, receptor-binding domain of EDA. These antibodies recognize epitopes overlapping the receptor-binding site and prevent EDA from binding and activating EDAR at close to stoichiometric ratios in in vitro binding and activity assays. The antibodies block EDA1 and EDA2 of both mammalian and avian origin and, in vivo, suppress the ability of recombinant Fc-EDA1 to rescue ectodermal dysplasia in Eda-deficient Tabby mice. Moreover, administration of EDA blocking antibodies to pregnant wild type mice induced in developing wild type fetuses a marked and permanent ectodermal dysplasia. These function-blocking anti-EDA antibodies with wide cross-species reactivity will enable study of the developmental and postdevelopmental roles of EDA in a variety of organisms and open the route to therapeutic intervention in conditions in which EDA may be implicated.

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.

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.