Common Variants in the Trichohyalin Gene Are Associated with Straight Hair in Europeans

Hair morphology is highly differentiated between populations and among people of European ancestry. Whereas hair morphology in East Asian populations has been studied extensively, relatively little is known about the genetics of this trait in Europeans. We performed a genome-wide association scan for hair morphology (straight, wavy, curly) in three Australian samples of European descent. All three samples showed evidence of association implicating the Trichohyalin gene (TCHH), which is expressed in the developing inner root sheath of the hair follicle, and explaining ∼6% of variance (p = 1.5 × 10⁻³¹). These variants are at their highest frequency in Northern Europeans, paralleling the distribution of the straight-hair EDAR variant in Asian populations.     

The adaptive variant EDARV370A is associated with straight hair in East Asians

Hair straightness/curliness is a highly heritable trait amongst human populations. Previous studies have reported European specific genetic variants influencing hair straightness, but those in East Asians remain unknown. One promising candidate is a derived coding variant of the ectodysplasin A receptor (EDAR), EDARV370A (370A), associated with several phenotypic changes of epidermal appendages. One of the strongest signals of natural selection in human genomes, 370A, has risen to high prevalence in East Asian and Native American populations, whilst being almost absent in Europeans and Africans. This striking frequency distribution and the pleiotropic nature of 370A led us to pursue if hair straightness, another epidermal appendage-related phenotype, is affected by this variant. By studying 1,718 individuals from four distinctive East Asian populations (Han, Tibetan, Mongolian, and Li), we found a significant association between 370A and the straight hair type in the Han (p = 2.90 × 10^?6), Tibetan (p = 3.07 × 10^?2), and Mongolian (p = 1.03 × 10^?5) populations. Combining all the samples, the association is even stronger (p = 5.18 × 10^?10). The effect of 370A on hair straightness is additive, with an odds ratio of 2.05. The results indicate very different biological mechanisms of straight hair in Europe and Asia, and also present a more comprehensive picture of the phenotypic consequences of 370A, providing important clues into the potential adaptive forces shaping the evolution of this extraordinary genetic variant.     

A Common Variation in EDAR Is a Genetic Determinant of Shovel-Shaped Incisors

Shovel shape of upper incisors is a common characteristic in Asian and Native American populations but is rare or absent in African and European populations. Like other common dental traits, genetic polymorphisms involved in the tooth shoveling have not yet been clarified. In ectodysplasin A receptor (EDAR), where dysfunctional mutations cause hypohidrotic ectodermal dysplasia, there is a nonsynonymous-derived variant, 1540C (rs3827760), that has a geographic distribution similar to that of the tooth shoveling. This allele has been recently reported to be associated with Asian-specific hair thickness. We aimed to clarify whether EDAR 1540C is also associated with dental morphology. For this purpose, we measured crown diameters and tooth-shoveling grades and analyzed the correlations between the dental traits and EDAR genotypes in two Japanese populations, inhabitants around Tokyo and in Sakishima Islands. The number of EDAR 1540C alleles in an individual was strongly correlated with the tooth-shoveling grade (p = 7.7 × 10⁻¹⁰). The effect of the allele was additive and explained 18.9% of the total variance in the shoveling grade, which corresponds to about one-fourth of the heritability of the trait reported previously. For data reduction of individual-level metric data, we applied a principal-component analysis, which yielded PC1-4, corresponding to four patterns of tooth size; this result implies that multiple factors are involved in dental morphology. The 1540C allele also significantly affected PC1 (p = 4.9 × 10⁻³), which denotes overall tooth size, and PC2 (p = 2.6 × 10⁻³), which denotes the ratio of mesiodistal diameter to buccolingual diameter.     

Ectodysplasin signalling genes and phenotypic evolution in sculpins (Cottus)

Despite their deeply conserved function among vertebrates, ectodysplasin (Eda) signalling genes are involved in microevolutionary change in humans and sticklebacks. If such a dual role is common, Eda signalling genes constitute hotspots for morphological evolution. Variation in sculpin (Cottus) skin prickling and body shape resembles patterns caused by variation in Eda signalling in sticklebacks. We mapped Eda signalling genes and performed quantitative trait locus mapping in crosses between Cottus rhenanus and Cottus perifretum. A genomic region containing the Eda receptor (Edar) was strongly associated with prickling and contributed to shape. The expression of Edar in developing prickles and skeletal elements in Cottus was confirmed by in situ hybridization. Coding sequence changes between Edar alleles in C. rhenanus and C. perifretum exceeded sequence differentiation in other vertebrates. However, it is likely that additional genetic elements besides coding changes affect the phenotypic variation. Although the phenotype in a natural hybrid lineage between C. rhenanus and C. perifretum resembles C. perifretum, the respective coding Edar alleles are not fully fixed (88.6%). Hence, our results support an involvement of Eda signalling in microevolutionary changes, but imply that the Edar gene is affected by multiple evolutionary processes that vary among freshwater sculpins.     

Identification of dkk4 as a target of Eda-A1/Edar pathway reveals an unexpected role of ectodysplasin as inhibitor of Wnt signalling in ectodermal placodes

The development of epithelial appendages, including hairs, glands and teeth starts from ectodermal placodes, and is regulated by interplay of stimulatory and inhibitory signals. Ectodysplasin-A1 (Eda-A1) and Wnts are high in hierarchy of placode activators. To identify direct targets of ectodysplasin pathway, we performed microarray profiling of genes differentially regulated by short exposure to recombinant Eda-A1 in embryonic eda^−/− skin explants. Surprisingly, there were only two genes with obvious involvement in Wnt pathway: dkk4 (most highly induced gene in the screen), and lrp4. Both genes colocalized with Eda-A1 receptor Edar in placodes of ectodermal organs. They were upregulated upon Edar activation while several other Wnt associated genes previously suggested as Edar targets were unaffected. However, low dkk4 and lrp4 expression was retained in the absence of NF-κB signalling in eda^−/− hair placodes. We provide evidence that this expression was dependent on Wnt activity present prior to Eda-A1/Edar signalling. Dkk4 was recently suggested as a key Wnt antagonist regulating lateral inhibition essential for correct patterning of hair follicles. Several pieces of evidence suggest Lrp4 as a Wnt inhibitor, as well. The finding that Eda-A1 induces placode inhibitors was unexpected, and underlines the importance of delicate fine-tuning of signalling during placode formation.     

Association of the OCA2 Polymorphism His615Arg with Melanin Content in East Asian Populations: Further Evidence of Convergent Evolution of Skin Pigmentation

The last decade has witnessed important advances in our understanding of the genetics of pigmentation in European populations, but very little is known about the genes involved in skin pigmentation variation in East Asian populations. Here, we present the results of a study evaluating the association of 10 Single Nucleotide Polymorphisms (SNPs) located within 5 pigmentation candidate genes (OCA2, DCT, ADAM17, ADAMTS20, and TYRP1) with skin pigmentation measured quantitatively in a sample of individuals of East Asian ancestry living in Canada. We show that the non-synonymous polymorphism rs1800414 (His615Arg) located within the OCA2 gene is significantly associated with skin pigmentation in this sample. We replicated this result in an independent sample of Chinese individuals of Han ancestry. This polymorphism is characterized by a derived allele that is present at a high frequency in East Asian populations, but is absent in other population groups. In both samples, individuals with the derived G allele, which codes for the amino acid arginine, show lower melanin levels than those with the ancestral A allele, which codes for the amino acid histidine. An analysis of this non-synonymous polymorphism using several programs to predict potential functional effects provides additional support for the role of this SNP in skin pigmentation variation in East Asian populations. Our results are consistent with previous research indicating that evolution to lightly-pigmented skin occurred, at least in part, independently in Europe and East Asia.     

A replication study confirmed the <Emphasis Type="Italic">EDAR</Emphasis> gene to be a major contributor to population differentiation regarding head hair thickness in Asia

Hair morphology is a highly divergent phenotype among human populations. We recently reported that a nonsynonymous SNP in the ectodysplasin A receptor (EDAR 1540T/C) is associated with head hair fiber thickness in an ethnic group in Thailand (Thai-Mai) and an Indonesian population. However, these Southeast Asian populations are genetically and geographically close, and thus the genetic contribution of EDAR to hair morphological variation in the other Asian populations has remained unclear. In this study, we examined the association of 1540T/C with hair morphology in a Japanese population (Northeast Asian). As observed in our previous study, 1540T/C showed a significant association with hair cross-sectional area (P = 2.7 × 10⁻⁶) in Japanese. When all populations (Thai-Mai, Indonesian, and Japanese) were combined, the association of 1540T/C was stronger (P = 3.8 × 10⁻¹⁰) than those of age, sex, and population. These results indicate that EDAR is the genetic determinant of hair thickness as well as a strong contributor to hair fiber thickness variation among Asian populations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ethnic related selection for an ADH Class I variant within East Asia

Background

The alcohol dehydrogenases (ADH) are widely studied enzymes and the evolution of the mammalian gene cluster encoding these enzymes is also well studied. Previous studies have shown that the ADH1B*47His allele at one of the seven genes in humans is associated with a decrease in the risk of alcoholism and the core molecular region with this allele has been selected for in some East Asian populations. As the frequency of ADH1B*47His is highest in East Asia, and very low in most of the rest of the world, we have undertaken more detailed investigation in this geographic region.

Methodology/Principal Findings

Here we report new data on 30 SNPs in the ADH7 and Class I ADH region in samples of 24 populations from China and Laos. These populations cover a wide geographic region and diverse ethnicities. Combined with our previously published East Asian data for these SNPs in 8 populations, we have typed populations from all of the 6 major linguistic phyla (Altaic including Korean-Japanese and inland Altaic, Sino-Tibetan, Hmong-Mien, Austro-Asiatic, Daic, and Austronesian). The ADH1B genotyping data are strongly related to ethnicity. Only some eastern ethnic phyla or subphyla (Korean-Japanese, Han Chinese, Hmong-Mien, Daic, and Austronesian) have a high frequency of ADH1B*47His. ADH1B haplotype data clustered the populations into linguistic subphyla, and divided the subphyla into eastern and western parts. In the Hmong-Mien and Altaic populations, the extended haplotype homozygosity (EHH) and relative EHH (REHH) tests for the ADH1B core were consistent with selection for the haplotype with derived SNP alleles. In the other ethnic phyla, the core showed only a weak signal of selection at best.

Conclusions/Significance

The selection distribution is more significantly correlated with the frequency of the derived ADH1B regulatory region polymorphism than the derived amino-acid altering allele ADH1B*47His. Thus, the real focus of selection may be the regulatory region. The obvious ethnicity-related distributions of ADH1B diversities suggest the existence of some culture-related selective forces that have acted on the ADH1B region.     

Ectodysplasin,Edar and TNFRSF19 are expressed in complementary and overlapping patterns during mouse embryogenesis

Ectodysplasin (Eda), a member of the tumor necrosis factor (TNF) superfamily, and its receptor Edar are necessary components of ectodermal organ development. Analysis of their expression patterns and mutant phenotypes has shown that during mouse hair and tooth development they may be involved in signalling between separate epithelial compartments. Here we have analysed ectodysplasin and Edar expression in other embryonic mouse tissues, and show that Edar mRNA is confined to the epithelium. Ectodysplasin and Edar are expressed in separate epithelial compartments in the developing brain and the lacrimal gland. In the salivary gland ectodysplasin is expressed in the mesenchyme and Edar in the epithelium. This is the first indication of ectodysplasin-Edar signalling between the epithelium and the mesenchyme. We also studied the expression pattern of a related TNF receptor, TNFRSF19, and show that it is expressed in an overlapping domain with Edar in the tooth, mammary gland, whiskers, and limb bud suggesting a potentially redundant role.     

Fast and robust detection of ancestral selective sweeps

There are few methods tailored for detecting signals of positive selection in populations directly ancestral to multiple descendent populations. We introduce the ancestral branch statistic (ABS), a four‐population summary statistic for identifying selective sweeps occurring in the direct ancestor of a pair of populations. Simulations show that ABS performs at least as well as, and often better under model violations, than the complementary likelihood approach of 3P‐CLR across diverse selection scenarios and parameter values. We first applied ABS to contemporary human genomic data to identify genes that may have been adaptive in ancestral East Asian populations, uncovering the well‐established candidate EDAR, as well as a novel candidate SLC35F3, which encodes a putative thiamine transporter that may have been involved in adaptation to eating polished grains. Next, we performed scans with ancient European genomic data to reexamine evidence of recent positive selection in ancestral Europeans. The MCM6/LCT cluster and the SLC45A2 and HERC2 genes are strong outliers, agreeing with previous studies. Novel candidates, such as SLC30A9 and CYP1A2, may have been involved in adaptation to local nutrient sufficiency and lifestyle changes. Finally, we provide open‐source software, CalcABS, which can perform genomic scans of ancestral sweeps with ABS from population allele frequency data.     

Involvement of and Interaction between WNT10A and EDA Mutations in Tooth Agenesis Cases in the Chinese Population

Background

Dental agenesis is the most common, often heritable, developmental anomaly in humans. Although WNT10A gene mutations are known to cause rare syndromes associated with tooth agenesis, including onycho-odontodermal dysplasia (OODD), Schöpf-Schulz-Passarge syndrome (SSPS), hypohidrotic ectodermal dysplasia (HED), and more than half of the cases of isolated oligodontia recently, the genotype-phenotype correlations and the mode of inheritance of WNT10A mutations remain unclear. The phenotypic expression with WNT10A mutations shows a high degree of variability, suggesting that other genes might function with WNT10A in regulating ectodermal organ development. Moreover, the involvement of mutations in other genes, such as EDA, which is also associated with HED and isolated tooth agenesis, is not clear. Therefore, we hypothesized that EDA mutations interact with WNT10A mutations to play a role in tooth agenesis. Additionally, EDA, EDAR, and EDARADD encode signaling molecules in the Eda/Edar/NF-κB signaling pathways, we also checked EDAR and EDARADD in this study.

Methods

WNT10A, EDA, EDAR and EDARADD were sequenced in 88 patients with isolated oligodontia and 26 patients with syndromic tooth agenesis. The structure of two mutated WNT10A and two mutated EDA proteins was analyzed.

Results

Digenic mutations of both WNT10A and EDA were identified in 2 of 88 (2.27%) isolated oligodontia cases and 4 of 26 (15.38%) syndromic tooth agenesis cases. No mutation in EDAR or EDARADD gene was found.

Conclusions

WNT10A and EDA digenic mutations could result in oligodontia and syndromic tooth agenesis in the Chinese population. Moreover, our results will greatly expand the genotypic spectrum of tooth agenesis.     

Edar/Eda interactions regulate enamel knot formation in tooth morphogenesis

tabby and downless mutant mice have apparently identical defects in teeth, hair and sweat glands. Recently, genes responsible for these spontaneous mutations have been identified. downless (Dl) encodes Edar, a novel member of the tumour necrosis factor (TNF) receptor family, containing the characteristic extracellular cysteine rich fold, a single transmembrane region and a death homology domain close to the C terminus. tabby (Ta) encodes ectodysplasin-A (Eda) a type II membrane protein of the TNF ligand family containing an internal collagen-like domain. As predicted by the similarity in adult mutant phenotype and the structure of the proteins, we demonstrate that Eda and Edar specifically interact in vitro. We have compared the expression pattern of Dl and Ta in mouse development, taking the tooth as our model system, and find that they are not expressed in adjacent cells as would have been expected. Teeth develop by a well recorded series of epithelial-mesenchymal interactions, similar to those in hair follicle and sweat gland development, the structures found to be defective in tabby and downless mice. We have analysed the downless mutant teeth in detail, and have traced the defect in cusp morphology back to initial defects in the structure of the tooth enamel knot at E13. Significantly, the defect is distinct from that of the tabby mutant. In the tabby mutant, there is a recognisable but small enamel knot, whereas in the downless mutant the knot is absent, but enamel knot cells are organised into a different shape, the enamel rope, showing altered expression of signalling factors (Shh, Fgf4, Bmp4 and Wnt10b). By adding a soluble form of Edar to tooth germs, we were able to mimic the tabby enamel knot phenotype, demonstrating the involvement of endogenous Eda in tooth development. We could not, however, reproduce the downless phenotype, suggesting the existence of yet another ligand or receptor, or of ligand-independent activation mechanisms for Edar. Changes in the structure of the enamel knot signalling centre in downless tooth germs provide functional data directly linking the enamel knot with tooth cusp morphogenesis. We also show that the Lef1 pathway, thought to be involved in these mutants, functions independently in a parallel pathway.     

Signaling and subcellular localization of the TNF receptor Edar

Tabby and downless mutant mice have identical phenotypes characterized by deficient development of several ectodermally derived organs such as teeth, hair, and sweat glands. Edar, encoded by the mouse downless gene and defective in human dominant and recessive forms of autosomal hypohidrotic ectodermal dysplasia (EDA) syndrome, is a new member of the tumor necrosis factor (TNF) receptor superfamily. The ligand of Edar is ectodysplasin, a TNF-like molecule mutated in the X-linked form of EDA and in the spontaneous mouse mutant Tabby. We have analyzed the response of Edar signaling in transfected cells and show that it activates nuclear factor-kappaB (NF-kappaB) in a dose-dependent manner. When Edar was expressed at low levels, the NF-kappaB response was enhanced by coexpression of ectodysplasin. The activation of NF-kappaB was greatly reduced in cells expressing mutant forms of Edar associated with the downless phenotype. Overexpression of Edar did not activate SAPK/JNK nor p38 kinase. Even though Edar harbors a death domain its overexpression did not induce apoptosis in any of the four cell lines analyzed, nor was there any difference in apoptosis in developing teeth of wild-type and Tabby mice. Additionally, we show that the subcellular localization of dominant negative alleles of downless is dramatically different from that of recessive or wild-type alleles. This together with differences in NF-kappaB responses suggests an explanation for the different mode of inheritance of the different downless alleles.     

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.     

GST M1-T1 null Allele Frequency Patterns in Geographically Assorted Human Populations: A Phylogenetic Approach

Genetic diversity in drug metabolism and disposition is mainly considered as the outcome of the inter-individual genetic variation in polymorphism of drug-xenobiotic metabolizing enzyme (XME). Among the XMEs, glutathione-S-transferases (GST) gene loci are an important candidate for the investigation of diversity in allele frequency, as the deletion mutations in GST M1 and T1 genotypes are associated with various cancers and genetic disorders of all major Population Affiliations (PAs). Therefore, the present population based phylogenetic study was focused to uncover the frequency distribution pattern in GST M1 and T1 null genotypes among 45 Geographically Assorted Human Populations (GAHPs). The frequency distribution pattern for GST M1 and T1 null alleles have been detected in this study using the data derived from literatures representing 44 populations affiliated to Africa, Asia, Europe, South America and the genome of PA from Gujarat, a region in western India. Allele frequency counting for Gujarat PA and scattered plot analysis for geographical distribution among the PAs were performed in SPSS-21. The GST M1 and GST T1 null allele frequencies patterns of the PAs were computed in Seqboot, Gendist program of Phylip software package (3.69 versions) and Unweighted Pair Group method with Arithmetic Mean in Mega-6 software. Allele frequencies from South African Xhosa tribe, East African Zimbabwe, East African Ethiopia, North African Egypt, Caucasian, South Asian Afghanistan and South Indian Andhra Pradesh have been identified as the probable seven patterns among the 45 GAHPs investigated in this study for GST M1-T1 null genotypes. The patternized null allele frequencies demonstrated in this study for the first time addresses the missing link in GST M1-T1 null allele frequencies among GAHPs.     

Regulation of hair follicle development by the TNF signal ectodysplasin and its receptor Edar

X-linked and autosomal forms of anhidrotic ectodermal dysplasia syndromes (HED) are characterized by deficient development of several ectodermal organs, including hair, teeth and exocrine glands. The recent cloning of the genes that underlie these syndromes, ectodysplasin (ED1) and the ectodysplasin A receptor (EDAR), and their identification as a novel TNF ligand-receptor pair suggested a role for TNF signaling in embryonic morphogenesis. In the mouse, the genes of the spontaneous mutations Tabby (Ta) and downless (dl) were identified as homologs of ED1 and EDAR, respectively. To gain insight into the function of this signaling pathway in development of skin and hair follicles, we analyzed the expression and regulation of Eda and Edar in wild type as well as Tabby and Lef1 mutant mouse embryos. We show that Eda and Edar expression is confined to the ectoderm and occurs in a pattern that suggests a role of ectodysplasin/Edar signaling in the interactions between the ectodermal compartments and the formation and function of hair placodes. By using skin explant cultures, we further show that this signaling pathway is intimately associated with interactions between the epithelial and mesenchymal tissues. We also find that Ta mutants lack completely the placodes of the first developing tylotrich hairs, and that they do not show patterned expression of placodal genes, including Bmp4, Lef1, Shh, Ptch and Edar, and the genes for beta-catenin and activin A. Finally, we identified activin as a mesenchymal signal that stimulates Edar expression and WNT as a signal that induces Eda expression, suggesting a hierarchy of distinct signaling pathways in the development of skin and hair follicles. In conclusion, we suggest that Eda and Edar are associated with the onset of ectodermal patterning and that ectodysplasin/edar signaling also regulates the morphogenesis of hair follicles.     

Development and Structure of Internal Glands and External Glandular Trichomes in Pogostemon cablin

Pogostemon cablin possesses two morphologically and ontogenetically different types of glandular trichomes, one type of bristle hair on the surfaces of leaves and stems and one type of internal gland inside the leaves and stems. The internal gland originates from elementary meristem and is associated with the biosynthesis of oils present inside the leaves and stems. However, there is little information on mechanism for the oil biosynthesis and secretion inside the leaves and stems. In this study, we identified three kinds of glandular trichome types and two kinds of internal gland in the Pogostemon cablin. The oil secretions from internal glands of stems and leaves contained lipids, flavones and terpenes. Our results indicated that endoplasmic reticulum and plastids and vacuoles are likely involved in the biosynthesis of oils in the internal glands and the synthesized oils are transported from endoplasmic reticulum to the cell wall via connecting endoplasmic reticulum membranes to the plasma membrane. And the comparative analysis of the development, distribution, histochemistry and ultrastructures of the internal and external glands in Pogostemon cablin leads us to propose that the internal gland may be a novel secretory structure which is different from external glands.     

An Indel Polymorphism in the MtnA 3' Untranslated Region Is Associated with Gene Expression Variation and Local Adaptation in Drosophila melanogaster

Insertions and deletions (indels) are a major source of genetic variation within species and may result in functional changes to coding or regulatory sequences. In this study we report that an indel polymorphism in the 3’ untranslated region (UTR) of the metallothionein gene MtnA is associated with gene expression variation in natural populations of Drosophila melanogaster. A derived allele of MtnA with a 49-bp deletion in the 3'' UTR segregates at high frequency in populations outside of sub-Saharan Africa. The frequency of the deletion increases with latitude across multiple continents and approaches 100% in northern Europe. Flies with the deletion have more than 4-fold higher MtnA expression than flies with the ancestral sequence. Using reporter gene constructs in transgenic flies, we show that the 3'' UTR deletion significantly contributes to the observed expression difference. Population genetic analyses uncovered signatures of a selective sweep in the MtnA region within populations from northern Europe. We also find that the 3’ UTR deletion is associated with increased oxidative stress tolerance. These results suggest that the 3'' UTR deletion has been a target of selection for its ability to confer increased levels of MtnA expression in northern European populations, likely due to a local adaptive advantage of increased oxidative stress tolerance.