Figures
Abstract
We conducted a multi-stage genome-wide association study of natural hair color in more than 10,000 men and women of European ancestry from the United States and Australia. An initial analysis of 528,173 single nucleotide polymorphisms (SNPs) genotyped on 2,287 women identified IRF4 and SLC24A4 as loci highly associated with hair color, along with three other regions encompassing known pigmentation genes. We confirmed these associations in 7,028 individuals from three additional studies. Across these four studies, SLC24A4 rs12896399 and IRF4 rs12203592 showed strong associations with hair color, with p = 6.0×10−62 and p = 7.46×10−127, respectively. The IRF4 SNP was also associated with skin color (p = 6.2×10−14), eye color (p = 6.1×10−13), and skin tanning response to sunlight (p = 3.9×10−89). A multivariable analysis pooling data from the initial GWAS and an additional 1,440 individuals suggested that the association between rs12203592 and hair color was independent of rs1540771, a SNP between the IRF4 and EXOC2 genes previously found to be associated with hair color. After adjustment for rs12203592, the association between rs1540771 and hair color was not significant (p = 0.52). One variant in the MATP gene was associated with hair color. A variant in the HERC2 gene upstream of the OCA2 gene showed the strongest and independent association with hair color compared with other SNPs in this region, including three previously reported SNPs. The signals detected in a region around the MC1R gene were explained by MC1R red hair color alleles. Our results suggest that the IRF4 and SLC24A4 loci are associated with human hair color and skin pigmentation.
Author Summary
It has been a longstanding hypothesis that human pigmentation is tightly regulated by genetic variation. However, very few genes have been identified that contain common genetic variants associated with human pigmentation. We scanned the genome for genetic variants associated with natural hair color and other pigmentary characteristics in a multi-stage study of more than 10,000 men and women of European ancestry from the United States and Australia. We identified IRF4 and SLC24A4 as loci highly associated with hair color, along with three other regions encompassing known pigmentation genes. Further work is needed to identify the causal variants at these loci. Improved understanding of the genetic determinants of human pigmentation may help identify the molecular mechanisms of pigmentation-associated conditions such as the tanning response and skin cancers.
Citation: Han J, Kraft P, Nan H, Guo Q, Chen C, Qureshi A, et al. (2008) A Genome-Wide Association Study Identifies Novel Alleles Associated with Hair Color and Skin Pigmentation. PLoS Genet 4(5): e1000074. https://doi.org/10.1371/journal.pgen.1000074
Editor: Gonçalo Abecasis, University of Michigan, United States of America
Received: December 24, 2007; Accepted: April 14, 2008; Published: May 16, 2008
Copyright: © 2008 Han 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.
Funding: This work was funded in part by NCI grants CA128080, CA132175, CA087969, CA049449, CA88363, and CA098233.
Competing interests: The authors have declared that no competing interests exist.
Introduction
There is substantial variation in human pigmentation within and across populations. Ultraviolet radiation (UV) exposure is the most important environmental factor influencing evolutionary selection pressure on pigmentation. In addition to UV-induced DNA damage, UVA can break down folic acid, and the major source of circulating vitamin D is synthesized in UVB-exposed skin. Because both nutrients are essential for human reproduction, it has been proposed that human pigmentation is selected, at least in part, to optimize levels of these two UV-related nutrients [1]. UV light is also the major environmental risk factor for skin cancer in humans. Red and blonde hair color, light skin pigmentation, and blue eye color are major host susceptibility factors for skin cancer [2].
Human pigmentation is a polygenic quantitative trait with high heritability [3]–[5]. A handful of genes underlying rare, extreme pigmentation phenotypes have been discovered [6], although until recently, only six genes were known to contain common genetic variants associated with human pigmentation in the normal range (MC1R, TYR, OCA2, SLC24A5, MATP, and ASIP). The proteins that these genes encode contribute to the control of melanin production and the maturation of melanosomes in melanogenesis, which determines human pigmentation. With new technologies that enable genotyping of hundreds of thousands of single nucleotide polymorphisms (SNPs), together with new insights into the structure of variation in the human genome [7], it is now possible to scan the genome in an agnostic manner in search of common genetic variants associated with human pigmentation. To identify common genetic variants associated with variation in natural diversity of human pigmentation, we performed a genome-wide association study (GWAS) of natural hair color in 2,287 U.S. women of European ancestry using data on 528,173 SNPs genotyped as part of the Cancer Genetic Markers of Susceptibility breast cancer GWAS [8]. Promising SNPs were examined in four additional studies with data on hair color and other pigmentation phenotypes: 870 U.S. women controls free of diagnosed skin cancer from a skin-cancer case-control study; 3,750 U.S. women from a diabetes case-control study; 2,405 U.S. men from a diabetes case-control study; and 1,440 parents of twins from an ongoing Australian family-based study of genetic and environmental factors contributing to the development of pigmented nevi [9].
Results
The frequencies of pigmentary phenotypes collected in the 5 component studies are presented in Table 1. The samples were broadly similar.
We compared the distribution of observed p-values from each of the 528,173 SNPs in the GWAS with those expected under the global null hypothesis that none of the tested SNPs is associated with natural hair color (Figure 1). The distribution of the observed p-values for the crude analyses that restricted analysis to women of self-reported European ancestry but did not further adjust for potential population stratification shows evidence for systematic bias: the genomic control inflation factor for the crude analyses (the ratio of the median observed test statistic to the theoretical median) is λGC = 1.24. This systematic bias is most likely due to confounding by latent population stratification. Hair color varies along a light-dark gradient from northern to southern Europe, so it will be associated with any SNP marker whose minor allele frequency also varies along a North-South gradient, even if that marker is not in linkage disequilibrium (LD) with a causal hair-color locus [10]. Adjusting for the top four principal components of genetic variance [11] eliminated most of the apparent residual confounding due to population stratification (λGC = 1.02 for the adjusted analyses); further control for up to 50 principal components did not alter the λGC. All of the association results from the initial GWAS reported below are from analyses that adjusted for the top four principal components of genetic variation.
p-values smaller than 10−8 are not plotted.
The GWAS identified several genomic locations as potentially associated with hair color (Figure 2). Of 528,173 SNPs tested, the 38 SNPs with the most extreme p-values associated with hair color are listed in Table 2.
Only p-values smaller than 0.05 are plotted.
We selected 31 of these 38 SNPs for further study in an independent sample. The remaining seven SNPs were in strong LD (r2>0.8) with one of these 31 SNPs (Table 2). The sample consisted of 870 controls of European ancestry from a nested case-control study of skin cancer within the Nurses' Health Study (NHS). Thirty of the 31 attempted SNPs were genotyped successfully.
Twenty-two of these 30 SNPs showed very strong evidence for association with natural hair color (p<9.5×10−8 = 0.05/528,173) in a pooled analysis of the initial GWAS and the validation sample (Table 3). Of the remaining eight SNPs, three showed very strong evidence for association with hair color either after excluding women with red hair or when comparing women with red hair to those without (Table 3). The associations between these 30 SNPs and with tanning ability and skin color are presented in Table S1 and Table S2.
The SNP rs12203592 in intron 4 of the IRF4 gene was strongly associated with hair color in the initial GWAS and validation study (black to red, pooled p value for trend = 8.5×10−28; black to blonde, pooled p value for trend = 7.1×10−49). The percentage of residual variation in hair color from black to blonde explained by this SNP after controlling for the top four principal components of genetic variation was 7.0%. This SNP is within 69.7 kb of two SNPs (rs4959270 and rs1540771) that were identified by a recent GWAS of natural hair color in women of European ancestry resident in Iceland [12]. However, neither of these variants, which lie between EXOC2 (SEC5L1) and IRF4, was as strongly associated with natural hair color in our initial GWAS as the IRF4 SNP rs12203592 (Figure 3). In our GWAS, the p values for association between hair color (black to blonde) and rs4959270 and rs1540771 were 2.9×10−4 and 0.007, respectively, and those for tanning ability were 0.002 and 0.001, respectively. In fact, the p-value for association between rs12203592 and natural hair color was more than 13 orders of magnitude smaller than the p-value for any other SNP on chromosome 6. This should not be taken as evidence that the loci that influence hair color in Iceland are different from those for the rest of Europe; rather, the previous GWAS may have failed to identify rs12203592 because this SNP is not on the Illumina HumanHap300 array used in that study, while it is on the Illumina HumanHap550 array used here.
-log10 p-values from the primary test of association with hair color in the initial GWAS (top panel), from the test of association with hair color excluding individuals with red hair (middle), and from the test of association with red hair versus non-red hair (bottom). The region plotted spans Chr6:95272.789899 (NCBI build 35).
We genotyped rs12203592 in an additional 6,155 individuals of predominantly European ancestry from the United States, including 3,750 women from the NHS and 2,405 men from the Health Professionals Follow-up Study (HPFS), and in an additional 1440 individuals of European ancestry from Australia (Queensland Institute of Medical Research (QIMR)). The association with hair color showed strong reproducibility; the independent p-values for association with hair color (black to blonde) in these three follow-up studies was 3.2×10−40, 2.8×10−35, and 4.6×10−23 respectively, and the pooled p-value across all five studies was 1.5×10−137 (Table 4). There was no statistical evidence for heterogeneity in the magnitude or direction of the hair color-minor allele correlation across the studies. The rs12203592 SNP was also highly associated with skin color (6.2×10−14), eye color (6.1×10−13), and tanning ability (3.9×10−89) in subsets of these individuals for which this information was available (Table 4). This variant allele was associated with lighter skin color, less tanning ability, and blue/light eye color (Table 5).
On the same chromosome, 145.8 kb centromeric from the IRF4 rs12203592, the SNP rs6918152 in the EXOC2 gene was associated with hair color (black to blond) in the initial GWAS, the NHS skin cancer controls, and the Australian samples (Tables 3 and 6). These two SNPs are in very weak LD (r2 = 0.04). Genotypes for the SNP rs1540771 previously reported by Sulem et al. [12] were available in the initial scan and the Australian samples. In a mutually adjusted multivariable regression of rs12203592, rs6918152, rs1540771, the strength of the association between the first SNP with hair color was attenuated but remained significant (p = 2.7×10−17 for IRF4 rs12203592 in the Australian samples) (Table 6). The previously reported SNP rs1540771 was not significant (p>0.05) after adjustment for the other SNPs, and association between rs6918152 and hair color was much weaker in the GWAS and no longer significant in the Australian samples after adjustment. While the IRF4 SNP rs12203592 was also associated with skin color, eye color and tanning ability, the EXOC2 SNP rs6918152 was not associated with these phenotypes. These results suggest that the IRF4 SNP rs12203592 is most likely to be in strong LD with the causal variant in this region.
The rs12896399 SNP 15.5 kb upstream of the SLC24A4 gene was highly associated with light hair color, and relatively weakly associated with less tanning ability in the pooled analysis of four studies (p = 6.0×10−62 for hair color, and p = 0.01 for tanning ability). The percentage of residual variation in hair color from black to blonde explained by this SNP after controlling for the top four principal components of genetic variation was 2.6%. This variant was also associated with blue/light eye color (p = 2.9×10−6 in the HPFS set). The SLC24A4 gene belongs to a family of potassium-dependent sodium/calcium exchangers. At least two other members of this family are associated with skin pigmentation. The SLC24A5 gene was recently shown to be involved in skin pigmentation in both zebrafish and humans [13].
Another member of this family, MATP (SLC45A2), is a pigmentation gene transcriptionally regulated by MITF [14],[15]. We identified the SNP rs28777 in the MATP gene from the GWAS, and the association with hair color was replicated in the controls of the skin cancer study (pooled P value = 8.9×10−14). This SNP was also associated with skin color (pooled P value = 9.5×10−4) and tanning ability (pooled P value = 2.2×10−10). Three SNPs in the MATP gene have been associated with human pigmentation: rs16891982 (Phe374Leu), rs26722 (Glu272Lys), and rs13289 C/G (-1721 in the promoter region) [16],[17]. We genotyped these three SNPs in the controls of the skin cancer study. None of the three previously reported SNPs were in LD with rs28777 (r2≤0.01), which is an intronic SNP. A multivariable analysis mutually adjusting for rs28777, rs16891982, rs26722, and rs13289 simultaneously showed that only rs16891982 remained significant in the model (P = 0.036 for hair color (black to blonde), P = 0.016 for tanning ability, and P = 0.0009 for skin color) and other SNPs became non-significant (p>0.05). These data suggested that rs16891982 is most likely to be the causal variant or in strong LD with the causal variant in the MATP gene.
Eleven SNPs spanning 1 Mb on chromosome 15 were strongly associated with hair color in the initial GWAS. These SNPs were located in the OCA2 5′ regulatory region and the HERC2 gene region and included the 3 SNPs reported previously with eye color: rs7495174, p = 7.1×10−7; rs6497268, p = 8.5×10−15; and rs11855019, p = 2.1×10−12) [18]. In an analysis mutually adjusting for all 11 SNPs simultaneously, only the HERC2 SNP rs12913832 (not on the HumanHap 300 version used in Sulem et al. [12]) remained significantly associated with hair color (p = 2.73×10−32) and tanning ability (p = 3.03×10−9). The associations between all other SNPs and hair color became non-significant (p>0.05). This suggested that the SNP rs12913832 was most likely to be in strong linkage disequilibrium with the causal variant in this region. The percentage of residual variation in hair color from black to blonde explained by this SNP after controlling for the top four principal components of genetic variation was 10.7%.
We observed 12 SNPs on chromosome 16 associated with hair color in the GWAS, spanning >756 kb. The MC1R gene, well established to be associated with red hair color, is located within this region. We had previously genotyped 7 common MC1R variants among the NHS skin cancer controls [19]. The analysis mutually adjusting for all 19 SNPs in the controls of the skin cancer study indicates that the signals that we detected in this region were mainly due to the three MC1R red hair color alleles (Arg151Cys, Arg160Trp, and Asp294His) (Table S3). The pairwise LD among these 19 SNPs was very low (the pattern of LD across these 19 SNPs is shown in Figure S1).
Discussion
It has been a longstanding hypothesis that human pigmentation is tightly regulated by genetic variation. However, very few genes have been identified that contain common genetic variants associated with human pigmentation. We conducted a genome-wide association study of hair color and identified several new variants associated with variation in hair color, skin color, eye color, and tanning ability among individuals of European ancestry.
Among the loci identified from our GWAS, the IRF4 and SLC24A4 loci had not been linked to human pigmentation before we began our study. Recently, Sulem et al. [12] reported a pigmentation GWAS using 316,515 SNPs in the Icelandic population. The associations of hair color in our GWAS with the 60 SNPs reported by Sulem et al. are listed in Table S4. These authors identified two SNPs (rs4959270 and rs1540771) between the EXOC2 and IRF4 genes in relation to freckles, hair color, and skin sensitivity to sun [12]. We identified a SNP in intron 4 of the IRF4 gene, not genotyped by Sulem et al. [12], with a much stronger association than they observed with hair color and tanning ability. The IRF4 gene product is a member of the interferon regulatory factor family of transcription factors [20]–[23], which are involved in the regulation of gene expression in response to interferon and other cytokines. The IRF4 gene encodes a B-cell proliferation/differentiation protein, which has been proposed as a sensitive and specific marker for conventional primary and metastatic melanomas and benign melanocytic nevi [24].
The SNP rs12896399 upstream of the SLC24A4 gene showed strong association with hair color in our study and that of Sulem et al. [12]. In addition, we identified three chromosomal regions adjacent to the previously known pigmentation genes: MC1R, OCA2, and MATP. The MC1R gene encodes a 317-amino acid 7-pass-transmembrane G protein coupled receptor and has been shown as the rate-limiting step in the activation of the cAMP pathway in terms of melanin production. Although the LD between the MC1R variants and surrounding highly significant SNPs was relatively low, the multivariable models mutually adjusting for all surrounding SNPs suggests that the signals that we identified on chromosome 16 were explained by the functional variants in the MC1R gene.
A previous report showed that three SNPs in intron 1 of the OCA2 gene were associated with eye, skin, and hair color [18]. We identified a SNP (rs12913832) in the upstream HERC2 gene with a much stronger association with hair color and tanning ability. Because the HERC2 gene has not been linked to human pigmentation to our knowledge, the SNP may be involved in the regulation of the expression or the function of the OCA2 gene. Similarly, Sulem et al. identified rs1667394 (∼165 kb upstream from rs12913832, r2 = 0.58) as the strongest hit in this region in their study [12]. However, we found a stronger association of rs12913832 with hair color than that of rs1667394 in our study.
Similar to our findings, Sulem et al. reported associations in the regions encompassing MC1R, OCA2, and SLC24A4. In addition, they reported loci in two established pigmentation genes, TYR and KITLG, which were not among the 38 SNPs with the strongest associations with hair color (black to red) that we sought to genotype in additional samples. The p-value for the association between KITLG rs12821256 and hair color (black to red) in the initial GWAS was 0.0002; the p-value for the association with hair color excluding women with red hair (black to blonde) was 1.28×10−8. The p-values for association between TYR rs1393350 and hair color coded as black to red or black to blonde were 0.05 and 0.02, respectively. We additionally identified another previously reported pigmentation gene from our GWAS, the MATP gene that was not reported by Sulem et al. [12]. In our analysis of testing the trend across hair color from black to blonde, there were no other loci reaching genome-wide significance level.
Two of the four regions we found to be associated with variation in hair color among Europeans without red hair (MATP and HERC2/OCA2) show strong evidence of recent positive selection, based on a comparison of allele frequencies across samples from three continental populations (Africa, Asia, and Europe) [25]. Both of the markers we identified in the two remaining regions showed significant differences in allele frequency across the HapMap CEU, CHB, JPT and YRI panels: the IRF4 SNP rs12203592 was monomorphic in CHB, JPT and YRI panels (the minor allele in Europeans was absent from these samples); the minor allele among Europeans for SLC24A4 SNP rs12896399 (G) was the major allele for the CHB and YRI panels, with the G allele frequency in the YRI sample being above 99%. Moreover, all of the markers with strongest association with hair color in these four regions were significantly associated with one or more of the top four principal components of genetic variation (Table S5), suggesting that allele frequencies for these markers also vary among European populations. Because we adjusted for latent population structure using these four principal components—and there are multiple lines of evidence suggesting these regions influence hair color among Europeans—we believe it unlikely that the strong associations we see between these markers and hair color are solely due to population stratification bias. Rather it is likely that differences in the distribution of hair color across Europe are due in part to differences in allele frequencies at these loci and other as-yet-unknown loci. Taken together, these four regions explain approximately 21.9% of the residual variation in hair color (black-blond) after adjusting for the top four principal components of genetic variation. (Conversely, after adjusting for these four regions, the top four principal components of genetic variation explain 2.6% of the residual variation in hair color.)
In our study of men and women of European ancestry we focused on the most statistically significant associations from our GWAS among women, identifying the IRF4 variant as reproducibly associated with human pigmentation. Further work is needed to identify the causal variant at this locus. Because a subset of true associations would be weakly associated with outcome in any given GWAS, large-scale replication is necessary for confirmation, and some true associations may be missed if they are not carried forward into replication studies. In this regard, the precomputed rankings and P values for all the SNPs included in the GWAS conducted in the NHS are freely available (http://www.channing.harvard.edu/nhs/publications/index.shtml) for others to use in subsequent studies.
Materials and Methods
Nurses' Health Study (NHS)
The NHS was established in 1976, when 121,700 female U.S. registered nurses between the ages of 30 and 55, residing in 11 larger U.S. states, completed and returned the initial self-administered questionnaire on their medical histories and baseline health-related exposures, forming the basis for the NHS cohort. Biennial questionnaires with collection of exposure information on risk factors and (every 4 years since 1980) nutritional assessments have been collected prospectively. Along with exposures every 2 years, outcome data with appropriate follow-up of reported disease events, including melanoma and non-melanoma skin cancers, are collected. Overall, follow-up has been very high; after more than 20 years approximately 90% of participants continue to complete questionnaires. From May 1989 through September 1990, we collected blood samples from 32,826 participants in the NHS cohort. Subsequent follow-up has been greater than 99% for this subcohort. The information on natural hair color at age of 20 and childhood and adolescence tanning ability were collected in the 1982 questionnaire.
Initial GWAS
We initially performed genotyping in a nested case-control study of postmenopausal invasive breast cancer within the Nurses' Health Study (NHS) cohort [26] using the Illumina HumanHap550 array, as part of the National Cancer Institute's Cancer Genetic Markers of Susceptibility (CGEMS) Project [8]. We performed our initial genome-wide analysis on 528,173 SNPs in 2,287 women [8]. All cases and controls were self-described as being of European ancestry. Four samples were excluded because of evidence of intercontinental admixture. Controlling for breast cancer case-control status made no material difference to the GWAS results. Information on natural hair color at age 20 was collected in the NHS main questionnaire and grouped into five categories (black, dark brown, light brown, blonde, and red).
Detailed methods related to the initial GWAS were published previously [8], including genotyping and quality control, initial assessment of sample completion rates, assessment of SNP call rates, concordance rate, deviation from Hardy–Weinberg proportions in control DNA, and final sample selection and exclusion for association analysis.
The Controls in the Skin Cancer Nested Case-Control Study within the NHS
The promising SNPs from the initial GWAS were further genotyped among 870 controls in the skin cancer nested case-control study within the NHS. The distribution of risk factors for skin cancer in the subcohort of those who donated blood samples was very similar to that in the overall cohort [2]. A common control series was randomly selected from participants who gave a blood sample and were free of diagnosed skin cancer up to and including the questionnaire cycle in which the corresponding case was diagnosed.
Health Professionals Followup Study (HPFS)
In 1986, 51,529 men from all 50 U.S. States in health professions (dentists, pharmacists, optometrists, osteopath physicians, podiatrists, and veterinarians) aged 40–75 answered a detailed mailed questionnaire, forming the basis of the study. Between 1993 and 1994, 18,159 study participants provided blood samples by overnight courier. The information on natural hair color and eye color was collected in the 1988 questionnaire, and the information on tanning ability was asked in the 1992 questionnaire.
The Diabetes Nested Case-Control Studies within the NHS and HPFS
Two additional studies were used to genotype novel pigmentation loci: 3,750 samples from the nested case-control study of diabetes in the NHS and 2,405 samples from the nested case-control study of diabetes in the HPFS [27]. All samples that we used were cases and controls from these two studies. Cases were incident cases of diabetes after blood collection, and controls were matched on age. Controlling for case-control status made no material difference to the results.
There was no sample overlap among the initial GWAS, the skin cancer case-control study, and the two diabetes case-control studies. The study protocol was approved by the Institutional Review Board of Brigham and Women's Hospital and Harvard School of Public Health. Informed consent was obtained from all patients.
Australian Study from the Queensland Institute of Medical Research (QIMR)
The Australian sample comprised 1,442 parents of twins taking part in a long-running study of melanoma risk factors [9],[18],[28],[29]. Participants rated their own hair color (at age 20 years) on a five-point classification (blonde, light brown, dark brown, black, red), eye color (blue/grey, green/hazel, brown), and skin color (light, medium, or dark).
Statistical Analysis
For the primary analysis of hair color we regressed an ordinal coding for hair color (1 = black; 2 = dark brown; 3 = light brown; 4 = blonde; and 5 = red) on an ordinal coding for genotype (0, 1 or 2 copies of SNP minor allele) separately for each SNP that passed quality control filters [8]. Crude analyses that did not adjust for any other variables showed evidence of systematic bias (see results); as this bias was greatly reduced by adjusting for the four largest principal components of genetic variation, all subsequent association analyses in the initial GWAS included these four components in the regression model. These principal components were calculated for all individuals on the basis of ca. 10,000 unlinked markers using the EIGENSTRAT software [8],[11]. The top four eigenvectors were chosen on the basis of significant (p<0.05) Tracy-Wisdom tests [30]. Adjusting for up to the top 50 principal components did not further reduce the genomic control inflation factor λGC. We chose markers for genotyping in subsequent validation studies based on the p-values for association from the primary analysis. Partial correlation coefficients (i.e., adjusted r2 or the percent of residual variance explained by the SNP marker) were calculated from the linear regressions adjusted for the top four principal components of genetic variation.
There is some evidence that determinants of hair color may act along two phenotypic axes: red hair color versus non-red color and light-dark variation among those without red hair. For example, alleles at the MC1R locus primarily determine presence or absence of red hair [31]. Hence, we conducted further analyses among individuals without red hair and comparing those with red hair to those without to evaluate whether discovered loci act on one or both phenotypic axes. We regressed the ordinal coding for hair color on minor allele counts excluding the individuals with red hair and used logistic regression to test the association between the ordinal genotype coding and a binary red-hair phenotype (red vs. non-red hair color). The regression parameter beta refers to the mean change in hair color scoring (or change in log odds of red hair for red hair analyses) per copy of the SNP minor allele.
We also used linear regression to test association between minor allele counts and self-reported tanning in response to sunlight (1 = deep tan, 2 = average tan, 3 = light tan, 4 = no tan), eye color (1 = brown/dark, 2 = hazel/green/medium, and 3 = blue/light), and skin color (1 = black, 2 = medium, and 3 = fair). Pooled analyses of multiple studies were conducted by merging data sets and including separate baseline parameters for each study.
Genotyping in Followup Studies
The TaqMan/BioTrove assays on the 31 SNPs in the skin cancer controls were performed at the Dana Farber/Harvard Cancer Center Polymorphism Detection Core (primers and probe sequences are available on request). Two loci (IRF4 rs12203592 and SLC24A4 rs12896399) were further genotyped in diabetes samples in the NHS and HPFS studies using the Taqman assay. Laboratory personnel were blinded to the case-control status, and 10% blinded quality control samples were inserted to validate genotyping procedures; concordance for the blinded samples was 100%. Primers, probes, and conditions for genotyping assays are available upon request.
For the Australian study, genotyping was performed as a single multiplex reaction on the Sequenom high-throughput genotyping platform on IRF4 SNP rs12203592 and EXOC2 rs6918152 and the best pigmentation-associated SNPs in the region of 6p25 (rs1540771) reported by Sulem et al. [12].
Supporting Information
Figure S1.
The LD pattern of SNPs around the MC1R locus on Chromosome 16.
https://doi.org/10.1371/journal.pgen.1000074.s001
(0.06 MB PDF)
Table S1.
Associations between the 30 most significant SNPs, which were identified in the GWAS of hair color, and tanning ability in the GWAS, the skin cancer controls, and the pooled data.
https://doi.org/10.1371/journal.pgen.1000074.s002
(0.02 MB XLS)
Table S2.
Associations between the 30 most significant SNPs, which were identified in the GWAS of hair color, and skin color in the skin cancer controls.
https://doi.org/10.1371/journal.pgen.1000074.s003
(0.02 MB XLS)
Table S3.
Associations between SNPs around the MC1R locus on chromosome 16 and hair color and tanning ability.
https://doi.org/10.1371/journal.pgen.1000074.s004
(0.02 MB XLS)
Table S4.
The associations between hair color and the 60 SNPs in Supplementary Table 1 of Sulem et al. [12].
https://doi.org/10.1371/journal.pgen.1000074.s005
(0.03 MB XLS)
Table S5.
The associations between SNPs in the MATP, IRF4, SLC24A4, and HERC2 genes and EVs (from 1 to 10).
https://doi.org/10.1371/journal.pgen.1000074.s006
(0.02 MB XLS)
Acknowledgments
We thank the investigators involved in the CGEMS project, including Kevin B. Jacobs at Bioinformed Consulting Services, Gaithersburg, Maryland; Sholom Wacholder, Nick Orr, Kai Yu, Nilanjan Chatterjee, Joseph F. Fraumeni, Jr., Richard B. Hayes, Margaret Tucker, and Daniela S. Gerhard at the Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Department of Health and Human Services, Bethesda, Maryland, United States; and Meredith Yeager, Robert Welch, Amy Hutchinson, Andrew Crenshaw, and Zhaoming Wang at the Core Genotyping Facility, Advanced Technology Program, SAIC Frederick, NCI-Frederick, Frederick, Maryland, United States. We thank Pati Soule and Dr. Hardeep Ranu of the Dana Farber/Harvard Cancer Center High-Throughput Polymorphism Detection Core for sample handling and genotyping of the NHS samples. We are also indebted to the participants in all of these studies.
Author Contributions
Conceived and designed the experiments: JH PK HN AQ SH DD ZZ NM GM NH GT RH DH. Performed the experiments: DD ZZ NM GM NH GT RH DH. Analyzed the data: JH PK HN QG CC DD ZZ NM GM NH GT RH DH. Contributed reagents/materials/analysis tools: PK HN QG CC SH FH DD ZZ NM GM NH GT RH DH. Wrote the paper: JH PK HN DD ZZ NM GM NH GT RH DH.
References
- 1. Jablonski NG, Chaplin G (2000) The evolution of human skin coloration. J Hum Evol 39: 57–106.
- 2. Han J, Colditz GA, Hunter DJ (2006) Risk factors for skin cancers: a nested case-control study within the Nurses' Health Study. Int J Epidemiol 35: 1514–1521.
- 3. Clark P, Stark AE, Walsh RJ, Jardine R, Martin NG (1981) A twin study of skin reflectance. Ann Hum Biol 8: 529–541.
- 4. Frisancho AR, Wainwright R, Way A (1981) Heritability and components of phenotypic expression in skin reflectance of Mestizos from the Peruvian lowlands. Am J Phys Anthropol 55: 203–208.
- 5. Harrison GA, Owen JJ (1964) Studies on the inheritance of human skin colour. Ann Hum Genet 28: 27–37.
- 6. Sturm RA (2006) A golden age of human pigmentation genetics. Trends Genet 22: 464–468.
- 7. Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, et al. (2007) A second generation human haplotype map of over 3.1 million SNPs. Nature 449: 851–861.
- 8. Hunter DJ, Kraft P, Jacobs KB, Cox DG, Yeager M, et al. (2007) A genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancer. Nat Genet 39: 870–874.
- 9. Duffy DL, Box NF, Chen W, Palmer JS, Montgomery GW, et al. (2004) Interactive effects of MC1R and OCA2 on melanoma risk phenotypes. Hum Mol Genet 13: 447–461.
- 10. Campbell CD, Ogburn EL, Lunetta KL, Lyon HN, Freedman ML, et al. (2005) Demonstrating stratification in a European American population. Nat Genet 37: 868–872.
- 11. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, et al. (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38: 904–909.
- 12. Sulem P, Gudbjartsson DF, Stacey SN, Helgason A, Rafnar T, et al. (2007) Genetic determinants of hair, eye and skin pigmentation in Europeans. Nat Genet 39: 1443–1452.
- 13. Lamason RL, Mohideen MA, Mest JR, Wong AC, Norton HL, et al. (2005) SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans. Science 310: 1782–1786.
- 14. Newton JM, Cohen-Barak O, Hagiwara N, Gardner JM, Davisson MT, et al. (2001) Mutations in the human orthologue of the mouse underwhite gene (uw) underlie a new form of oculocutaneous albinism, OCA4. Am J Hum Genet 69: 981–988.
- 15. Du J, Fisher DE (2002) Identification of Aim-1 as the underwhite mouse mutant and its transcriptional regulation by MITF. J Biol Chem 277: 402–406.
- 16. Graf J, Hodgson R, van Daal A (2005) Single nucleotide polymorphisms in the MATP gene are associated with normal human pigmentation variation. Hum Mutat 25: 278–284.
- 17. Graf J, Voisey J, Hughes I, van Daal A (2007) Promoter polymorphisms in the MATP (SLC45A2) gene are associated with normal human skin color variation. Hum Mutat 28: 710–717.
- 18. Duffy DL, Montgomery GW, Chen W, Zhao ZZ, Le L, et al. (2007) A three-single-nucleotide polymorphism haplotype in intron 1 of OCA2 explains most human eye-color variation. Am J Hum Genet 80: 241–252.
- 19. Han J, Kraft P, Colditz GA, Wong J, Hunter DJ (2006) Melanocortin 1 receptor variants and skin cancer risk. Int J Cancer 119: 1976–1984.
- 20. Grossman A, Mittrucker HW, Nicholl J, Suzuki A, Chung S, et al. (1996) Cloning of human lymphocyte-specific interferon regulatory factor (hLSIRF/hIRF4) and mapping of the gene to 6p23-p25. Genomics 37: 229–233.
- 21. Eisenbeis CF, Singh H, Storb U (1995) Pip, a novel IRF family member, is a lymphoid-specific, PU.1-dependent transcriptional activator. Genes Dev 9: 1377–1387.
- 22. Matsuyama T, Grossman A, Mittrucker HW, Siderovski DP, Kiefer F, et al. (1995) Molecular cloning of LSIRF, a lymphoid-specific member of the interferon regulatory factor family that binds the interferon-stimulated response element (ISRE). Nucleic Acids Res 23: 2127–2136.
- 23. Yamagata T, Nishida J, Tanaka S, Sakai R, Mitani K, et al. (1996) A novel interferon regulatory factor family transcription factor, ICSAT/Pip/LSIRF, that negatively regulates the activity of interferon-regulated genes. Mol Cell Biol 16: 1283–1294.
- 24. Sundram U, Harvell JD, Rouse RV, Natkunam Y (2003) Expression of the B-cell proliferation marker MUM1 by melanocytic lesions and comparison with S100, gp100 (HMB45), and MelanA. Mod Pathol 16: 802–810.
- 25. Barreiro LB, Laval G, Quach H, Patin E, Quintana-Murci L (2008) Natural selection has driven population differentiation in modern humans. Nat Genet 40: 340–345.
- 26. Tworoger SS, Eliassen AH, Sluss P, Hankinson SE (2007) A prospective study of plasma prolactin concentrations and risk of premenopausal and postmenopausal breast cancer. J Clin Oncol 25: 1482–1488.
- 27. Qi L, van Dam RM, Meigs JB, Manson JE, Hunter D, et al. (2006) Genetic variation in IL6 gene and type 2 diabetes: tagging-SNP haplotype analysis in large-scale case-control study and meta-analysis. Hum Mol Genet 15: 1914–1920.
- 28. Palmer JS, Duffy DL, Box NF, Aitken JF, O'Gorman LE, et al. (2000) Melanocortin-1 receptor polymorphisms and risk of melanoma: is the association explained solely by pigmentation phenotype? Am J Hum Genet 66: 176–186.
- 29. Zhu G, Montgomery GW, James MR, Trent JM, Hayward NK, et al. (2007) A genome-wide scan for naevus count: linkage to CDKN2A and to other chromosome regions. Eur J Hum Genet 15: 94–102.
- 30. Patterson N, Price AL, Reich D (2006) Population structure and eigenanalysis. PLoS Genet 2(12): e190.
- 31. Rees JL (2004) The genetics of sun sensitivity in humans. Am J Hum Genet 75: 739–751.