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

doi:10.1371/journal.pone.0002209

PLoS ONE. 2008; 3(5): e2209.

Published online 2008 May 21. doi: 10.1371/journal.pone.0002209

PMCID: PMC2374902

Copyright Bryk et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

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

Jarosław Bryk,¹ Emilie Hardouin,¹ Irina Pugach,¹ David Hughes,¹ Rainer Strotmann,² Mark Stoneking,¹ and Sean Myles¹^*^¤

¹Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany

²Institute of Biochemistry, Molecular Biochemistry, Medical Faculty, University of Leipzig, Leipzig, Germany

Jason E. Stajich, Editor

University of California, Berkeley, United States of America

* E-mail: smm367/at/cornell.edu

Conceived and designed the experiments: MS SM JB. Performed the experiments: SM RS EH IP JB. Analyzed the data: MS SM JB. Contributed reagents/materials/analysis tools: DH RS. Wrote the paper: MS SM JB.

^¤Current address: Institute for Genomic Diversity, Cornell University, Ithaca, New York, United States of America

Received March 3, 2008; Accepted April 7, 2008.

This article has been cited by other articles in PMC.

Abstract

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.

Introduction

Humans expanded within and out of Africa between 50,000 and 100,000 years ago and now inhabit radically different physical and cultural environments around the globe. A long-standing question in anthropology has been: How have humans genetically adapted to these environments? Numerous recent studies have attempted to address this issue for a review, see [1]. Most of these studies aim to identify genomic regions that have experienced local positive selection (i.e. geographically-restricted positive selection) from genome-wide polymorphism data e.g. [2], [3], [4]. A small number of these genomic regions are associated with phenotypes that have long been candidates for local positive selection, such as lactose tolerance [5], [6] and skin pigmentation [7]–[12]. For the large majority of these regions, however, no association to a locally adaptive phenotype has been established and the precise location of functional adaptive variants remains elusive. Moreover, sexual selection may be a plausible explanation for at least some of these candidate regions. Further investigation of candidate regions, including the possible functional consequences of the genetic variation observed in such regions, is required to fully understand the impact of local positive selection on our species. Here, we investigate one candidate gene, EDAR, in detail.

A strong signature of positive selection in East Asians has been found in the genomic region containing the EDAR gene (HGNC [ratio]

2895) according to tests based on the allele frequency spectrum [4], [13], [14], haplotype structure [2], [3] and population differentiation [15]–[19]. EDAR is known to be involved in the development of hair follicles, teeth and sweat glands [20] and harbors a nonsynonymous single nucleotide polymorphism (SNP) (rs3827760) that results in a valine to alanine substitution at position 370 of the amino acid sequence (V370A). Here we test the function of V370A in vitro and provide several lines of evidence suggesting that there was positive selection on the derived 370A allele in East Asians prior to 10,000 years ago.

Results

Fst analysis

We assessed worldwide population differentiation for V370A by use of the Fst statistic [21]. Alleles that have been targets of local positive selection tend to have unusually high Fst values [22]–[24]. We genotyped the V370A polymorphism in the 53 worldwide populations of the CEPH Human Genome Diversity Panel CEPH-HGDP; [25] and compared Fst values for V370A to an empirical Fst distribution derived from 2750 autosomal markers (2540 SNPs [26] and 210 indels [27]) previously typed in the same set of samples. Global Fst, the degree of differentiation among all 53 populations, for V370A is 0.760. This value is higher than all Fst values from the empirical global Fst distribution. To examine the patterns of population differentiation at a more refined geographical scale, we calculated Fst for every pairwise comparison among the 53 populations and 7 geographic regions to produce 53×53 and 7×7 Fst matrices, respectively. Each pairwise Fst value for V370A was then compared to the corresponding empirical distribution of pairwise Fst values to generate a P value. Worldwide allele frequencies and the P value matrices from the Fst analysis are depicted in Figure 1.

Figure 1

Worldwide allele frequencies and population differentiation for V370A.

Functional assay

EDAR (Swiss-Prot [ratio]

Q9UNE0) is a cell-surface receptor that, upon binding to its ligand, induces an intracellular cascade leading to the activation of NF-κB, a transcription factor [20]. To investigate the functional consequences of the V370A polymorphism, we measured NF-κB activation in vitro using a luciferase reporter assay, from a HEK293 cell line heterologously expressing the 370V and 370A variants of EDAR cDNA. As a positive control, we also performed the same transfections with 370V and 370A EDAR cDNAs that harbored an additional disease mutation (375H; Figure 2) that had previously been shown to severely reduce NF-κB activation in vitro [28]. In agreement with previous work, the cDNAs carrying 375H showed a ~6 fold reduction in activation of NF-κB compared to cDNAs without 375H (Figure 3). Moreover, the derived 370A allele results in increased activation of NF-κB compared to the ancestral 370V allele on both the normal (p

0.018, univariate ANOVA with experiment as a random factor and clone as a fixed factor) and the disease background (p<0.004, two-tailed t-test assuming equal variances; Figure 3).

Figure 2

Multiple species alignment of the death domain of EDAR.

Figure 3

The derived 370A allele of EDAR results in enhanced activation of NF-κB in vitro.

Estimation of time since fixation

Finally, we assessed the timing of the selection event on 370A from ~22 kb of EDAR sequence from 23 individuals of Chinese ancestry from the Seattle SNPs data set (http://pga.gs.washington.edu/). We employed a method that estimates the time since fixation of a beneficial allele [29]. Only one chromosome from the data set carried the ancestral 370V allele and we removed it from the analysis and assumed that the 370A allele had reached fixation. The resulting estimate for the time since fixation is 10,740 years (95% CI

1133–73996; Figure 4).

Figure 4

Density plot of the posterior distribution of estimates of the time since fixation of the 370A allele.

Discussion

Recent scans of genome-wide polymorphism data have generated long lists of genomic regions that are believed to have been targeted by local positive selection, but few of these regions have been shown to harbor functional variants or to be linked to a putatively adaptive phenotype. Detailed investigation of these candidate genomic regions is required for a comprehensive picture of local human adaptation. The present study provides an analysis of genetic variation and function at the EDAR gene, a candidate for positive selection in East Asians identified from recent genome scans [2]–[4], [13]–[16], [19].

We provide two lines of evidence supporting the hypothesis that the derived 370A allele in EDAR is functional and experienced positive selection in East Asians. First, worldwide population differentiation for 370A as measured by Fst is highly unusual. Figure 1 demonstrates that 370A has a highly unusual worldwide frequency distribution, supporting a scenario in which the 370A allele was driven to high frequency in East Asians and Native Americans by positive selection (Figure 1). Two Central/South Asia populations, the Uygur and the Hazara, have intermediate frequencies of 370A (0.44 and 0.5, respectively), in agreement with their close genetic relationship to East Asians [30], [31]. The 370A allele is also found at low frequency in Melanesia (0.12), and was likely introduced there via the recent Austronesian expansion [32]. Otherwise 370A is absent in Africans and Papuans and is observed at very low frequency in most of Central/South Asia, Europe and the Middle East (Figure 1).

The second line of evidence suggesting an adaptive functional role for 370A stems from its location in the amino acid sequence of EDAR. V370A is located in the death domain, a protein interaction module, of EDAR. The death domain of EDAR is highly conserved (Figure 2) and interacts with the death domain of EDARADD, an intracellular ligand to EDAR [33], [34]. This interaction initiates an intracellular signalling cascade that results in the activation of the transcription factor NF-κB [35]. Therefore, V370A may alter binding affinity with the death domain of EDARADD and thereby influence the activation of NF-κB. Moreover, seven nonsynonymous substitutions in the death domain of EDAR cause hypohidrotic ectodermal dysplasia (HED) in humans [36], a disease characterized by sparse and thin hair, missing teeth and the absence of sweat glands (OMIM [ratio]

604095; see Figure 2). In particular, the R375H substitution, only 5 amino acids upstream from V370A, results in a loss of affinity for EDARADD and reduced NF-κB activation [28]. We therefore tested the function of V370A in vitro and found that the 370A allele does differ from the ancestral 370V allele in that it results in enhanced NF-κB activation (Figure 3).

Our results contradict a recent report from Fujimoto et al. [37] in which the 370A allele was shown to reduce NF-κB activation in vitro [37]. We are confident that our results are correct for two reasons. First, in addition to observing enhanced NF-κB activation for the derived 370A allele on a normal genetic background, we measured NF-κB activation of the 370V and 370A alleles on the background of a disease mutation (375H) that was previously demonstrated to result in significantly reduced NF-κB activation [28]. In agreement with [28], we observed significantly reduced NF-κB activation in both clones carrying the 375H disease mutation. Moreover, NF-κB activation was significantly higher in the derived 370A+375H construct than in the ancestral 370V+375H construct (Figure 3). Thus, we observed enhanced NF-κB activation for the derived 370A allele in two independent constructs. Second, East Asians have thicker hair than Europeans and Africans [38] and an increase in NF-κB activation is arguably more likely to lead to the thicker East Asian hair phenotype. This is because a decrease in NF-κB activation, as observed for carriers of the 375H allele that causes hypohidrotic ectodermal dysplasia, is associated with thin hair. We suspect that the large doses of plasmid DNA (300 ng) and long post-transfection incubation period (48 h) could have induced cell death in the experiments of Fujimoto et al. [37]. Although EDAR-induced cell death is a matter of controversy [39], [40], several features of cell death (detachment, rounding and membrane permeation) are observed 36 hours after transfection with a high dosage of plasmid (500 ng) [39]. Control cells transfected with the same amount of an empty vector do not display features of cell death, suggesting that simply overdosing the cells with plasmid is not responsible for cell death [39]. Thus, we speculate that the reduced NF-κB activation from 370A observed by Fujimoto et al. [37] could be the result of induced cell death and that 370A in fact enhances NF-κB activation.

Our estimate of the time since fixation of 370A in a sample of 45 Chinese chromosomes is 10,740 years (Figure 4). This estimate involves several assumptions (see Materials and Methods) and should be interpreted with caution. Nevertheless the result suggests that 370A was likely at high frequency before the colonization of the Americas 10,500–14,000 years ago [41]–[43]. Thus, the high frequency of 370A in Native Americans (see Figure 1) is most likely due to positive selection prior to migrations from Asia to America.

In summary, we have demonstrated that the worldwide frequency distribution of 370A is highly unusual and that 370A was likely rising in frequency by positive selection in East Asia prior to 10,000 years ago. In addition, we have shown that the 370A allele results in enhanced NF-κB activation in vitro. What was the source of the selection pressure on 370A and what effect may 370A have on the phenotype? Since EDAR is involved in ectodermal development, 370A might be expected to affect teeth, hair, skin, nails and/or sweat glands. Fujimoto et al. [37] recently noted an association between 370A and hair thickness. Replication of this result is desirable since correction for population structure was inadequate: only a single SNP was used to correct for population structure. Nevertheless, the results of Fujimoto et al. [37] are suggestive, especially since East Asians have thicker hair than Europeans and Africans [38]. These observations lead us to question why thicker hair may have been advantageous in ancestral East Asian environments. Of course, thicker hair may not have been adaptive at all and may simply be the result of phenotypic hitchhiking: selection on 370A may have targeted a different phenotype (e.g. tooth morphology [16]) and hair thickness may have resulted as a by-product of this selection. Sexual selection also remains a possibility.

Regardless of the nature of the selective force, our results provide compelling evidence that positive selection has acted on the 370A allele in EDAR. In addition, our finding that 370A results in increased NF-κB activation suggests further lines for investigation: in particular, how does this increased NF-κB activation influence the expression of the target genes regulated by NF-κB? Such future studies will lead to a more complete understanding of the phenotypic effect of 370A, and will permit more explicit tests of hypotheses concerning the possible selection pressure(s) responsible for its rapid increase in frequency in East Asia.

Methods

Experimental Procedures

We genotyped the V370A polymorphism (rs3827760) in the HGDP-CEPH Human Genome Diversity Cell Line Panel [25] using the SNaPshot™ minisequencing kit (Applied Biosystems, Foster City, CA). Atypical and related individuals were removed [44], which resulted in 952 individuals from 53 populations with 10 missing genotypes. One population (Burusho) was out of Hardy-Weinberg equilibrium (p

0.008), but this value becomes non-significant after correction for multiple comparisons. We obtained an empirical Fst distribution from 2750 autosomal markers (2540 SNPs [31] and 210 indels [32]) previously typed in 927 individuals from the CEPH-HGDP panel. Fst values for V370A were calculated from the same set of 927 individuals to allow for an unbiased comparison to the empirical distribution. Each P value was calculated as the proportion of Fst values from the empirical distribution that were greater than or equal to the observed Fst value.

Conservation of EDAR [NP_071731.1] was assessed by aligning all orthologs of EDAR obtained from release 52 of Homologene (http://www.ncbi.nlm.nih.gov/entrez/query.fcgidbhomologene, Homologene ID

7699) using MUSCLE version 3.6 [45].

We analysed 22094 bp of resequencing data from the EDAR gene from 23 individuals of Asian origin (http://pga.mbt.washington.edu/) to investigate the age of the selective event. The method [29] assumes that the selective sweep has reached fixation, so the single non-derived haplotype was removed, resulting in 45 chromosomes for analysis. The method also assumes that the sequence data is contiguous, but the available data included 7 sequenced regions distributed across 98136 bp of the genome. We therefore concatenated the 7 regions to produce one contiguous sequence of 22094 bp. Simulations were performed to obtain a sample from the posterior distribution of the time since fixation (T), conditional on the number of segregating sites (17), Tajima's D (−2.13) and the number of haplotypes (13). The following parameters were used: μ

2.3×10⁻⁸ [46], [47]; ρ

0.001 [46], [48]; N_e

10,000 [49], and a generation time of 20 years. As the observed number of haplotypes is affected by recombination events across the ~98 kb of EDAR sequence, we increased ρ by a factor of 5 in our simulations of ~22 kb of sequence. We found, however, that point estimates and confidence intervals did not vary widely when different values of ρ, ranging from 0.0001 to 0.005, were used (data not shown). We used the mode of the distribution of T as the point estimate of the time since fixation as recommended in [29].

EDAR cDNA was generated from human fetal skin total RNA (Stratagene) using oligo-dT primers and SuperScript II reverse transcriptase (Invitrogen) and amplified with primers covering amino acid positions 31–448. The PCR product was cloned into the pCMV-Tag2c vector (Stratagene) and several inserts were sequenced to determine the allelic variant of EDAR corresponding to amino acid position 370. Clones resulting in valine or alanine at this position were termed p370V and p370A, respectively. The clones underwent site-directed mutagenesis using the QuikChange II Kit (Stratagene) to generate the disease variant with histidine at position 375 on both ancestral and derived backgrounds (p370V+375H and p370A+375H). Whole inserts in mutated clones were sequenced to ensure correctness.

For the in vitro functional study, HEK293 cells (RZPD, Germany) were seeded in a 96-well plate at a density of 5×10⁴ cells/well and cultured in DMEM supplemented with 10% (vol) FBS (Invitrogen). After approximately 48 hours, cells were transfected with 50 ng/well of total DNA, consisting of equal amounts (ng) of 370V or 370A (each with or without disease variant), pNF-κB-Luc vector (Clonetech) and pRL4.74 vector (Promega). Transfection was performed using FuGene HD reagent (Roche) diluted with OptiMEM medium (Invitrogen) in 100 µl total volume in 5 replicates for each version of the EDAR construct. Cells transfected with empty pCMV-Tag2 vector were used as a background signal control. After 18 h, luciferase signals were read using the Dual-Glo Luciferase Assay System (Promega) and a Victor-2 luminometer (PerkinElmer). The signal from pNF-κB-Luc was normalized to the signal from pRL4.74 and the background was subtracted. An 18 h timepoint was selected based on previous time-curve experiments (data not shown).

Acknowledgments

We thank Denis Headon, Wolfgang Enard, Kerstin Ullmann, Susan Ptak, Montgomery Slatkin and Hua Chen for useful discussions and Birgit Nickel and Richard E. Green for technical assistance.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was supported by the German Bundesministerium für Bildung und Forschung (BMBF: NGFN2) and the Max Planck Society.

References

1. Myles S, Hughes DA, Stoneking M. Santos C, Lima M, editors. Detecting local selection in humans. Recent Advances in Molecular Biology and Evolution: Applications to Biological Anthropology. Research Signpost. 2008.

2. Voight BF, Kudaravalli S, Wen X, Pritchard JK. A Map of Recent Positive Selection in the Human Genome. PLoS Biol. 2006;4:e72. [PMC free article] [PubMed]

3. Tang K, Thornton KR, Stoneking M. A New Approach for Using Genome Scans to Detect Recent Positive Selection in the Human Genome. PLoS Biol. 2007;5:e171. [PMC free article] [PubMed]

4. Williamson SH, Hubisz MJ, Clark AG, Payseur BA, Bustamante CD, et al. Localizing Recent Adaptive Evolution in the Human Genome. PLoS Genet. 2007;3:e90. [PubMed]

5. Bersaglieri T, Sabeti PC, Patterson N, Vanderploeg T, Schaffner SF, et al. Genetic signatures of strong recent positive selection at the lactase gene. Am J Hum Genet. 2004;74:1111–1120. [PubMed]

6. Tishkoff SA, Reed FA, Ranciaro A, Voight BF, Babbitt CC, et al. Convergent adaptation of human lactase persistence in Africa and Europe. Nat Genet. 2007;39:31–40. [PMC free article] [PubMed]

7. Lamason RL, Mohideen M-APK, Mest JR, Wong AC, Norton HL, et al. SLC24A5, a Putative Cation Exchanger, Affects Pigmentation in Zebrafish and Humans. Science. 2005;310:1782–1786. [PubMed]

8. Soejima M, Tachida H, Ishida T, Sano A, Koda Y. Evidence for Recent Positive Selection at the Human AIM1 Locus in a European Population. Mol Biol Evol. 2006;23:179–188. [PubMed]

9. Izagirre N, Garcia I, Junquera C, de la Rua C, Alonso S. A Scan for Signatures of Positive Selection in Candidate Loci for Skin Pigmentation in Humans. Mol Biol Evol. 2006;23:1697–1706. [PubMed]

10. Norton HL, Kittles RA, Parra E, McKeigue P, Mao X, et al. Genetic Evidence for the Convergent Evolution of Light Skin in Europeans and East Asians. Mol Biol Evol. 2007;24:710–722. [PubMed]

11. Lao O, de Gruijter JM, van Duijn K, Navarro A, Kayser M. Signatures of Positive Selection in Genes Associated with Human Skin Pigmentation as Revealed from Analyses of Single Nucleotide Polymorphisms. Ann Hum Genet. 2007;71:354–369. [PubMed]

12. Myles S, Somel M, Tang K, Kelso J, Stoneking M. Identifying genes underlying skin pigmentation differences among human populations. Hum Genet. 2007;120:613–621. [PubMed]

13. Carlson CS, Thomas DJ, Eberle MA, Swanson JE, Livingston RJ, et al. Genomic regions exhibiting positive selection identified from dense genotype data. Genome Res. 2005;15:1553–1565. [PubMed]

14. Kelley JL, Madeoy J, Calhoun JC, Swanson W, Akey JM. Genomic signatures of positive selection in humans and the limits of outlier approaches. Genome Res. 2006;16:980–989. [PubMed]

15. Akey JM, Eberle MA, Rieder MJ, Carlson CS, Shriver MD, et al. Population History and Natural Selection Shape Patterns of Genetic Variation in 132 Genes. PLoS Biol. 2004;2:e286. [PMC free article] [PubMed]

16. Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007;449:913–918. [PMC free article] [PubMed]

17. The International HapMap Consortium. A haplotype map of the human genome. Nature. 2005;437:1299. [PMC free article] [PubMed]

18. Barreiro LB, Laval G, Quach H, Patin E, Quintana-Murci L. Natural selection has driven population differentiation in modern humans. Nat Genet. 2008;40:340–345. [PubMed]

19. Myles S, Tang K, Somel M, Green RE, Kelso J, et al. Identification and analysis of high Fst regions from genome-wide SNP data from three human populations. Ann Hum Genet. 2008;72:99–110. [PubMed]

20. Botchkarev VA, Fessing MY. Edar Signaling in the Control of Hair Follicle Development. Journal of Investigative Dermatology Symposium Proceedings. 2005;10:247–251.

21. Wright S. The theory of gene frequencies. Chicago: University of Chicago Press; 1969.

22. Thornton KR, Jensen JD. Controlling the False-Positive Rate in Multilocus Genome Scans for Selection. Genetics. 2007;175:737–750. [PubMed]

23. Beaumont MA, Balding DJ. Identifying adaptive genetic divergence among populations from genome scans. Mol Ecol. 2004;13:969–980. [PubMed]

24. Pollinger JP, Bustamante CD, Fledel-Alon A, Schmutz S, Gray MM, et al. Selective sweep mapping of genes with large phenotypic effects. Genome Res. 2005;15:1809–1819. [PubMed]

25. Cann HM, de Toma C, Cazes L, Legrand MF, Morel V, et al. A human genome diversity cell line panel. Science. 2002;296:261–262. [PubMed]

26. Conrad DF, Jakobsson M, Coop G, Wen X, Wall JD, et al. A worldwide survey of haplotype variation and linkage disequilibrium in the human genome. Nat Genet. 2006;38:1251–1260. [PubMed]

27. Rosenberg NA, Mahajan S, Ramachandran S, Zhao C, Pritchard JK, et al. Clines, Clusters, and the Effect of Study Design on the Inference of Human Population Structure. PLoS Genet. 2005;1:e70. [PubMed]

28. Shimomura Y, Sato N, Miyashita A, Hashimoto T, Ito M, et al. A rare case of hypohidrotic ectodermal dysplasia caused by compound heterozygous mutations in the EDAR gene. J Invest Dermatol. 2004;123:649–655. [PubMed]

29. Przeworski M. Estimating the Time Since the Fixation of a Beneficial Allele. Genetics. 2003;164:1667–1676. [PubMed]

30. Rosenberg NA, Pritchard JK, Weber JL, Cann HM, Kidd KK, et al. Genetic Structure of Human Populations. Science. 2002;298:2381–2385. [PubMed]

31. Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, et al. Worldwide Human Relationships Inferred from Genome-Wide Patterns of Variation. Science. 2008;319:1100–1104. [PubMed]

32. Friedlaender JS, Friedlaender FR, Reed FA, Kidd KK, Kidd JR, et al. The genetic structure of Pacific Islanders. PLoS Genet. 2008;4:e19. [PubMed]

33. Headon DJ, Emmal SA, Ferguson BM, Tucker AS, Justice MJ, et al. Gene defect in ectodermal dysplasia implicates a death domain adapter in development. Nature. 2001;414:913–916. [PubMed]

34. Yan M, Zhang Z, Brady JR, Schilbach S, Fairbrother WJ, et al. Identification of a Novel Death Domain-Containing Adaptor Molecule for Ectodysplasin-A Receptor that Is Mutated in crinkled Mice. Curr Biol. 2002;12:409–413. [PubMed]

35. Thesleff I, Mikkola ML. Death receptor signaling giving life to ectodermal organs. Sci STKE. 2002;2002:PE22. [PubMed]

36. Chassaing N, Bourthoumieu S, Cossee M, Calvas P, Vincent M. Mutations in EDAR account for one-quarter of non-ED1-related hypohidrotic ectodermal dysplasia. Human Mutation. 2006;27:255–259. [PubMed]

37. Fujimoto A, Kimura R, Ohashi J, Omi K, Yuliwulandari R, et al. A scan for genetic determinants of human hair morphology: EDAR is associated with Asian hair thickness. Hum Mol Genet. 2008;17:835–843. [PubMed]

38. Franbourg A, Hallegot P, Baltenneck F, Toutain C, Leroy F. Current research on ethnic hair. J Am Acad Dermatol. 2003;48:S115–119. [PubMed]

39. Kumar A, Eby MT, Sinha S, Jasmin A, Chaudhary PM. The Ectodermal Dysplasia Receptor Activates the Nuclear Factor-kappa B, JNK, and Cell Death Pathways and Binds to Ectodysplasin A. J Biol Chem. 2001;276:2668–2677. [PubMed]

40. Koppinen P, Pispa J, Laurikkala J, Thesleff I, Mikkola ML. Signaling and Subcellular Localization of the TNF Receptor Edar. Experimental Cell Research. 2001;269:180. [PubMed]

41. Hey J. On the Number of New World Founders: A Population Genetic Portrait of the Peopling of the Americas. PLoS Biol. 2005;3:965–975.

42. Waters MR, Stafford TW., Jr Redefining the Age of Clovis: Implications for the Peopling of the Americas. Science. 2007;315:1122–1126. [PubMed]

43. Fagundes NJR, Ray N, Beaumont M, Neuenschwander S, Salzano FM, et al. Statistical evaluation of alternative models of human evolution. PNAS. 2007;104:17614–17619. [PubMed]

44. Rosenberg NA. Standardized Subsets of the HGDP-CEPH Human Genome Diversity Cell Line Panel, Accounting for Atypical and Duplicated Samples and Pairs of Close Relatives. Ann Hum Genet. 2006;70:841–847. [PubMed]

45. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res. 2004;32:1792–1797. [PMC free article] [PubMed]

46. Frisse L, Hudson RR, Bartoszewicz A, Wall JD, Donfack J, et al. Gene conversion and different population histories may explain the contrast between polymorphism and linkage disequilibrium levels. Am J Hum Genet. 2001;69:831–843. [PubMed]

47. Fischer A, Wiebe V, Paabo S, Przeworski M. Evidence for a Complex Demographic History of Chimpanzees. Mol Biol Evol. 2004;21:799–808. [PubMed]

48. Ptak SE, Hinds DA, Koehler K, Nickel B, Patil N, et al. Fine-scale recombination patterns differ between chimpanzees and humans. Nat Genet. 2005;37:429–434. [PubMed]

49. Voight BF, Adams AM, Frisse LA, Qian Y, Hudson RR, et al. Interrogating multiple aspects of variation in a full resequencing data set to infer human population size changes. PNAS. 2005;102:18508–18513. [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science