Derivation of marker gene signatures from human skin and their use in the interpretation trancriptional 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 specialised 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 transcriptomics 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 include 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 keratinisation 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

Distinct Impacts of Eda and Edar Loss of Function on the Mouse Dentition

The Eda-A1-Edar signaling pathway is involved in the development of organs with an ectodermal origin, including teeth. In mouse, mutants are known for both the ligand, Eda-A1 (Tabby), and the receptor, Edar (Downless). The adult dentitions of these two mutants have classically been considered to be similar. However, previous studies mentioned differences in embryonic dental development between EdaTa and Edardl-J mutants. A detailed study of tooth morphology in mutants bearing losses of functions of these two genes thus appears necessary to test the pattern variability induced by the developmental modifications. 3D-reconstructions of the cheek teeth have been performed at the ESRF (Grenoble, France) by X-ray synchrotron microtomography to assess dental morphology. The morphological variability observed in EdaTa and Edardl-J mutants have then been compared in detail. Despite patchy similarities, our detailed work on cheek teeth in EdaTa and Edardl-J mice show that all dental morphotypes defined in Edardl-J mice resolutely differ from those of EdaTa mice. This study reveals that losses of function of Eda and Edar have distinct impacts on the tooth size and morphology, contrary to what has previously been thought. The results indicate that unknown mechanisms of the Eda pathway are implicated in tooth morphogenesis. Three hypotheses could explain our results; an unexpected role of the Xedar pathway (which is influenced by the Eda gene product but not that of Edar), a more complex connection than has been appreciated between Edar and another protein, or a ligand-independent activity for Edar. Further work is necessary to test these hypotheses and improve our understanding of the mechanisms of development

Dkk4 and Eda Regulate Distinctive Developmental Mechanisms for Subtypes of Mouse Hair

The mouse hair coat comprises protective “primary” and thermo-regulatory “secondary” hairs. Primary hair formation is ectodysplasin (Eda) dependent, but it has been puzzling that Tabby (Eda-/y) mice still make secondary hair. We report that Dickkopf 4 (Dkk4), a Wnt antagonist, affects an auxiliary pathway for Eda-independent development of secondary hair. A Dkk4 transgene in wild-type mice had no effect on primary hair, but secondary hairs were severely malformed. Dkk4 action on secondary hair was further demonstrated when the transgene was introduced into Tabby mice: the usual secondary follicle induction was completely blocked. The Dkk4-regulated secondary hair pathway, like the Eda-dependent primary hair pathway, is further mediated by selective activation of Shh. The results thus reveal two complex molecular pathways that distinctly regulate subtype-based morphogenesis of hair follicles, and provide a resolution for the longstanding puzzle of hair formation in Tabby mice lacking Eda

A genome-wide association study identifies multiple loci for variation in human ear morphology

Here we report a genome-wide association study for non-pathological pinna morphology in over 5,000 Latin Americans. We find genome-wide significant association at seven genomic regions affecting: lobe size and attachment, folding of antihelix, helix rolling, ear protrusion and antitragus size (linear regression P values 2 × 10−8 to 3 × 10−14). Four traits are associated with a functional variant in the Ectodysplasin A receptor (EDAR) gene, a key regulator of embryonic skin appendage development. We confirm expression of Edar in the developing mouse ear and that Edar-deficient mice have an abnormally shaped pinna. Two traits are associated with SNPs in a region overlapping the T-Box Protein 15 (TBX15) gene, a major determinant of mouse skeletal development. Strongest association in this region is observed for SNP rs17023457 located in an evolutionarily conserved binding site for the transcription factor Cartilage paired-class homeoprotein 1 (CART1), and we confirm that rs17023457 alters in vitro binding of CART1

Positive Selection in East Asians for an EDAR Allele that Enhances NF-κB Activation

Genome-wide scans for positive selection in humans provide a promising approach to establish links between genetic variants and adaptive phenotypes. From this approach, lists of hundreds of candidate genomic regions for positive selection have been assembled. These candidate regions are expected to contain variants that contribute to adaptive phenotypes, but few of these regions have been associated with phenotypic effects. Here we present evidence that a derived nonsynonymous substitution (370A) in EDAR, a gene involved in ectodermal development, was driven to high frequency in East Asia by positive selection prior to 10,000 years ago. With an in vitro transfection assay, we demonstrate that 370A enhances NF-κB activity. Our results suggest that 370A is a positively selected functional genetic variant that underlies an adaptive human phenotype

A genome-wide association scan implicates <i>DCHS2, RUNX2, GLI3, PAX1</i> and <i>EDAR</i> in human facial variation

We report a genome-wide association scan for facial features in ∼6,000 Latin Americans. We evaluated 14 traits on an ordinal scale and found significant association (P values−8) at single-nucleotide polymorphisms (SNPs) in four genomic regions for three nose-related traits: columella inclination (4q31), nose bridge breadth (6p21) and nose wing breadth (7p13 and 20p11). In a subsample of ∼3,000 individuals we obtained quantitative traits related to 9 of the ordinal phenotypes and, also, a measure of nasion position. Quantitative analyses confirmed the ordinal-based associations, identified SNPs in 2q12 associated to chin protrusion, and replicated the reported association of nasion position with SNPs in PAX3. Strongest association in 2q12, 4q31, 6p21 and 7p13 was observed for SNPs in the EDAR, DCHS2, RUNX2 and GLI3 genes, respectively. Associated SNPs in 20p11 extend to PAX1. Consistent with the effect of EDAR on chin protrusion, we documented alterations of mandible length in mice with modified Edar function

Mutation screening of the Ectodysplasin-A receptor gene EDAR in hypohidrotic ectodermal dysplasia

Hypohidrotic ectodermal dysplasia (HED) can be caused by mutations in the X-linked ectodysplasin A (ED1) gene or the autosomal ectodysplasin A-receptor (EDAR) and EDAR-associated death domain (EDARADD) genes. X-linked and autosomal forms are sometimes clinically indistinguishable. For genetic counseling in families, it is therefore important to know the gene involved. In 24 of 42 unrelated patients with features of HED, we found a mutation in ED1. ED1-negative patients were screened for mutations in EDAR and EDARADD. We found mutations in EDAR in 5 of these 18 patients. One mutation, p.Glu354X, is novel. In EDARADD, a novel variant p.Ser93Phe, probably a neutral polymorphism, was also found. Clinically, there was a difference between autosomal dominant and autosomal recessive HED patients. The phenotype in patients with mutations in both EDAR alleles was comparable to males with X-linked HED. Patients with autosomal dominant HED had features comparable to those of female carriers of X-linked HED. The teeth of these patients were quite severely affected. Hypohidrosis and sparse hair were also evident, but less severe. This study confirms Chassaing et al's earlier finding that mutations in EDAR account for approximately 25% of non-ED1-related HED. Mutations leading to a premature stop codon have a recessive effect except when the stop codon is in the last exon. Heterozygous missense mutations in the functional domains of the gene may have a dominant-negative effect with much variation in expression. Patients with homozygous or compound heterozygous mutations in the EDAR gene have a more severe phenotype than those with a heterozygous missense, nonsense or frame-shift mutation