Genetic loci associated with skin pigmentation in African Americans and their effects on vitamin D deficiency
Authors:
Ken Batai aff001; Zuxi Cui aff002; Amit Arora aff003; Ebony Shah-Williams aff004; Wenndy Hernandez aff005; Maria Ruden aff006; Courtney M. P. Hollowell aff006; Stanley E. Hooker aff007; Madhavi Bathina aff007; Adam B. Murphy aff008; Carolina Bonilla aff009; Rick A. Kittles aff007
Authors place of work:
Department of Urology, University of Arizona, Tucson, Arizona, United States of America
aff001; Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, Ohio, United States of America
aff002; Department of Epidemiology and Biostatistics, University of Arizona, Tucson, Arizona, United States of America
aff003; Department of Medical and Molecular Genetics, Indiana University, Indianapolis, Indiana United States of America
aff004; Department of Medicine, University of Chicago, Chicago, Illinois, United States of America
aff005; Department of Surgery, Cook County Health and Hospitals System, Chicago, Illinois, United States of America
aff006; Division of Health Equities, Department of Population Sciences, City of Hope Comprehensive Cancer Center, Duarte, California, United States of America
aff007; Department of Urology, Northwestern University, Chicago, Illinois, United States of America
aff008; Departamento de Medicina Preventiva, Faculdade de Medicina, Universidade de São Paulo, São Paulo, Brazil
aff009; Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, United Kingdom
aff010
Published in the journal:
Genetic loci associated with skin pigmentation in African Americans and their effects on vitamin D deficiency. PLoS Genet 17(2): e1009319. doi:10.1371/journal.pgen.1009319
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1009319
Summary
A recent genome-wide association study (GWAS) in African descent populations identified novel loci associated with skin pigmentation. However, how genomic variations affect skin pigmentation and how these skin pigmentation gene variants affect serum 25(OH) vitamin D variation has not been explored in African Americans (AAs). In order to further understand genetic factors that affect human skin pigmentation and serum 25(OH)D variation, we performed a GWAS for skin pigmentation with 395 AAs and a replication study with 681 AAs. Then, we tested if the identified variants are associated with serum 25(OH) D concentrations in a subset of AAs (n = 591). Skin pigmentation, Melanin Index (M-Index), was measured using a narrow-band reflectometer. Multiple regression analysis was performed to identify variants associated with M-Index and to assess their role in serum 25(OH)D variation adjusting for population stratification and relevant confounding variables. A variant near the SLC24A5 gene (rs2675345) showed the strongest signal of association with M-Index (P = 4.0 x 10−30 in the pooled dataset). Variants in SLC24A5, SLC45A2 and OCA2 together account for a large proportion of skin pigmentation variance (11%). The effects of these variants on M-Index was modified by sex (P for interaction = 0.009). However, West African Ancestry (WAA) also accounts for a large proportion of M-Index variance (23%). M-Index also varies among AAs with high WAA and high Genetic Score calculated from top variants associated with M-Index, suggesting that other unknown genomic factors related to WAA are likely contributing to skin pigmentation variation. M-Index was not associated with serum 25(OH)D concentrations, but the Genetic Score was significantly associated with vitamin D deficiency (serum 25(OH)D levels less than 12 ng/mL) (OR, 1.30; 95% CI, 1.04–1.64). The findings support the hypothesis suggesting that skin pigmentation evolved responding to increased demand for subcutaneous vitamin D synthesis in high latitude environments.
Keywords:
Genetic loci – Genetics – Genome-wide association studies – Genomics – Human genetics – Single nucleotide polymorphisms – Vitamin D deficiency
Introduction
A number of studies have identified single nucleotide polymorphisms (SNPs) associated with pigmentation traits, such as skin, eye, and hair color, sensitivity to sun or tanning, and freckles, using genome-wide approaches [1–7]. Recent studies focusing on skin pigmentation variation in African populations identified novel loci associated with skin pigmentation [8,9]. Another study in Latin American countries identified additional novel loci [7]. A GWAS meta-analysis of skin pigmentation in admixed populations did not identify a novel variant, but validated some of the major findings from previous studies in African and admixed populations [10]. These studies suggest that population-specific variants affecting skin pigmentation variation exist. However, there has been no genome-wide association study (GWAS) of skin pigmentation aiming to understand the genetic variation for skin pigmentation in African Americans (AAs). Pigmentation traits are complex phenotypes that show great variation across human populations. There is evidence suggesting that such variation has been shaped by natural selection at different latitudes, to prevent DNA damage by ultraviolet radiation to the skin and/or to guarantee enough synthesis of vitamin D, given that vitamin D synthesis is initiated in the skin [11]. Sexual selection has also been implicated as an important factor in the evolution of skin pigmentation and may explain the observed difference in skin color between males and females [12–15].
The importance of elucidating the genetic basis of pigmentation traits expands beyond a better understanding of the evolutionary mechanisms that shaped some of the most visible phenotypic traits. It will also provide a better understanding of genetic risk factors for skin cancer [16–19] and vitamin D deficiency in different populations [11,20]. While many studies have explored the relationship between skin pigmentation gene variants and skin cancer risk, there is paucity of studies aiming to understand the relationship between skin pigmentation gene variants and vitamin D deficiency in AAs who are disproportionately affected by vitamin D deficiency. Identifying genetic variants that affect pigmentation traits is also important for forensic science [21,22]. Furthermore, our understanding of the genetic basis of skin pigmentation may have social implications given its conspicuous nature [23–26].
In order to further understand genetic factors that affect human skin pigmentation, we performed a GWAS in AAs. Then, we assessed the effects of skin pigmentation variants on serum 25-hydroxyvitamin D [25(OH)D] variation and vitamin D deficiency [serum 25(OH)D levels <12 ng/mL] in male AAs who were part of a study on serum 25(OH)D. Instead of assessing relative homogeneous populations with respect to skin pigmentation, we studied AAs, as it is well documented that this population exhibits a vast range of variation with respect to skin pigmentation and genetic ancestral proportions [27,28].
Materials and methods
Ethics statement
Written consent was obtained from all study subjects. This study was approved by the Institutional Review Boards at Howard University (IRB-99-MED-20) and Northwestern University (STU00005398).
Study participants
Study participants in the discovery dataset were recruited in Washington, D.C. and Chicago, IL. Participants in the replication set were recruited in Washington, D.C., Cincinnati, OH, and Chicago, IL. The participants from Washington, D.C. were recruited at Howard University for one of two studies on human pigmentation [29] or serum vitamin D levels and prostate cancer in AAs [30]. The participants from Chicago were recruited for a vitamin D and prostate cancer study at five hospitals as well as from community health events [31,32]. The participants from Cincinnati were recruited as a part of a study to investigate the relationship between self-reported race, genetic ancestry, socio-economic status, and skin color [33]. All of the study participants were self-identified AAs. After receiving consent, study coordinators administered a questionnaire and obtained demographic information. Blood samples for DNA analysis and serum 25(OH)D assays were collected at the time of recruitment.
Skin pigmentation measurements
Constitutive skin pigmentation in the sun-protected area of skin on the inner upper arm was measured by trained study coordinators or research assistants using a computerized portable narrow-band reflectometer, called DermaSpectrometer (Cyberderm, PA) [34,35]. The DermaSpectrometer output is expressed in terms of erythema (E) and melanin (M) indices from 0–100%, where higher M values denote higher pigment content. Three separate measurements were recorded for one arm and the average of M was used in all the analyses.
Serum vitamin D measurements
Serum 25(OH)D were measured previously [31], and in the current study, we included 734 individuals with skin color measurements. The serum samples were stored at −20°C until the 25(OH)D assay was performed. Total 25(OH)D concentration was measured using the Diasorin chemiluminescence immunoassay method in the Department of Pathology NorthShore University HealthSystem. In this study, we considered individuals with serum 25(OH)D levels <12 ng/mL as severely vitamin D deficient based on the Institute of Medicine dietary reference for vitamin D intake [36] as well as an Endocrine Society Clinical Practice Guideline [37,38].
Genotyping
The discovery set (n = 395) was genotyped using two different genome-wide SNP arrays. Washington, D.C. samples were genotyped using the Illumina 1M array as part of a GWAS for prostate cancer in African descent populations (n = 215) [39]. Additional samples (n = 180) from the Washington D.C. pigmentation gene study and Chicago vitamin D study were genotyped using the Affymetrix PanAFR array for this study. Following standard GWAS quality control (QC) recommendations, SNPs with low genotyping call rate (<95%), Hardy-Weinberg Equilibrium (HWE) P< 0.001, and minor allele frequency (MAF) <1% were removed. Poorly genotyped samples (<95%) and related individuals (Identity-by-descent PI-HAT >0.2) were also removed. After QC, both GWAS datasets were merged and principal component analysis was performed to assess if there were systematic differences between them. We did not find evidence of systematic error. Genotyping to replicate the findings from the discovery cohort was performed using the Agena Bioscience MassARRAY iPLEX platform (n = 681). The pooled dataset consists of a total of 1,076 AAs.
We also selected 24 previously identified SNPs associated with pigmentation traits in European populations and admixed populations for replication. Some of these SNPs were not successfully imputed in the Affymetrix dataset, so genotyping was performed for the Affymetrix dataset using the MassARRAY platform, resulting in a total of 1,066 participants in this replication analysis after QC. Genotyping of 38 SNPs in 8 vitamin D metabolic and signaling pathway genes was performed in our previously study [31], and SNP selection criteria was described previously [40]. These include previously GWAS identified SNPs that were associated with serum vitamin D levels in European populations [41,42]. To estimate genetic ancestry, a validated set of ancestry informative markers to estimate continental ancestry information in admixed populations [43] were genotyped for all the samples including three parental population sets. The parental populations used to estimate admixture proportions included 243 Europeans (from England, Germany, Ireland and Spain), 279 West Africans (from Cameroon, Nigeria and Sierra Leone), and 214 Native Americans (Cheyenne, Maya, Pima, Pueblo, and Mayans).
Imputation
Genome lift-over was performed for the Illumina dataset to map the SNPs using the Human Genome version 19 (hg19). Shapeit.v2 was used for phasing genotype data [44]. Imputation was performed using Impute v2.3.2 and 1000 Genome Project as the reference panel [45,46]. After the imputation, variants with Infor Score <0.5 were removed. Before imputation, there were 1,013,952 SNPs in the Illumina dataset and 1,714,384 SNPs in the Affymetrix dataset. After imputation, there were 11,065,735 SNPs in the Illumina dataset and 9,218,475 SNPs in the Affymetrix dataset. We pooled both datasets and removed variants with a genotype missing rate > 5%, HWE P<1.0 x 10−5, and MAF < 1%. A total of 7,169,107 SNPs were used for subsequent analyses.
Statistical analysis
M-Index was log-transformed to normalize the distribution in the population. For the analysis of M-Index in the GWAS discovery dataset, a linear model was used adjusting for age, sex, and the first 3 principal components (PCs). The model building process included up to 10 PCs in the regression model. The final model includes minimum number of PCs necessary to correct for population structure and reduce genomic inflation. In the replication and pooled analysis, West African Ancestry (WAA) was used instead of principal components. Individual admixture proportion was estimated using STRUCTURE v2.3, a model-based clustering method [47,48]. STRUCTURE was run under the admixture model using K = 3 ancestral populations (West African, European, and Native American). We used a burn-in length of 100,000 for 100,000 repetitions. For the GWAS analyses we used P<5.0 x 10−8 as the genome-wide significant threshold and P<0.05 as a statistically significant cutoff for replication.
A weighted Genetic Score was calculated using the top 3 and 10 associated SNPs for skin pigmentation. The weighted Genetic Score is sum of the effects of each SNP weighted by its estimated effect size (β) from regression model, GeneticScore=∑j=1m(χijβj), where m is number of SNP included and χij is the genotype for the ith individual and jth SNP (coded as 0, 1, and 2 for increase number of allele associated with darker skin pigmentation) [49]. One SNP from a single genomic region was included for calculation. When there were more than 2 SNPs with P<0.05 in a same genomic region, the SNP with the lowest P-value was used after conditional analysis to test if the SNP with the second lowest P-value was independently associated with M-Index by including the lead SNP in the region in the regression model. The Genetic Score calculated from top SNPs from 3 and 10 loci were initially assessed to estimate variance in skin pigmentation. Because the top 3 and top 10 SNPs accounted for a similar amount of skin pigmentation variation, subsequently, we focused our analysis using the Genetic Score estimated from the top 3 SNPs. Sex-specific Genetic Scores were also calculated. First, linear regression analysis was performed separately for males and females for the top 3 SNPs associated with M-Index in the pooled dataset. Then, β coefficients obtained for each SNP in males and females were used for calculation of sex-specific Genetic Scores.
Associations with serum 25(OH)D levels were tested using linear regression models adjusting for age, UV season (season of blood draw), total vitamin D intake, recruitment site, and WAA as described previously [31]. Log-transformed serum 25(OH)D levels were used in the linear regression analysis. Genetic Score from vitamin D metabolic and signaling pathway gene variants were calculated using the same formula described above using the SNPs associated with serum 25(OH)D levels in our dataset. Binary logistic regression was used to examine the association between skin pigmentation gene variants and vitamin D deficiency adjusting for age, UV season, total vitamin D intake, recruitment site, WAA, and also Genetic Score calculated from vitamin D metabolic pathway gene variants associated with serum 25(OH)D levels. Statistical analysis was performed using PLINK 1.07 [50], SPSS (IBM Corp., Armonk, NY), and R.
Results
Study participants’ characteristics
Because the majority of study participants came from prostate cancer studies, our study participants tend to be older (mean age of 57.2 in the discovery and 51.4 in the replication dataset, Table 1). Men were also over-represented in both datasets (90.9% and 85.9% male in the discovery and replication dataset, respectively). The discovery dataset had older and more male study participants. Mean M-Index and WAA were similar in the discovery and replication datasets. Serum 25(OH)D measurements were available for 734 men in the discovery and replication datasets. Mean serum 25(OH)D concentration was 19.5 ng/mL, and over 50% of participants had 25(OH)D levels that were categorized as deficient to severely deficient.
Variants associated with skin pigmentation
One locus reached genome-wide significance in the discovery dataset (Fig 1A and 1B), and a SNP, rs2675345 near the SLC24A5 gene on chromosome 15, showed the strongest signal of association (P = 8.4 x 10−14). The previously identified SNP, rs1426654, was the second most significantly associated SNP (P = 9.9 x 10−14) [51,52]. These two SNPs were strongly linked with r2 = 0.97. We did not find evidence of genomic inflation (S1 Fig), and inclusion of additional PCs did not change the analysis results. The second signal of association was found on the TRHDE (thyrotropin-releasing hormone degrading enzyme) gene on chromosome 12 (S2 Fig). The SNP, rs11179301, showed the strongest association in the region (P = 6.28 x 10−7).
Associations between skin pigmentation and 4 previously identified loci in African populations [8,9] were explored. Index SNPs from these four loci (SNX13, SMARCA2/VLDLR, TMEM138/DDB1, and MFSD12) were not successfully genotyped or imputed in our study. SNP rs2093835 was the second most strongly associated SNP in the SMARCA2/VLDLR region in African populations [9], but this SNP was not significantly associated with M-Index (P = 0.59) in our AA samples. These 4 genomic regions were further explored to identify variants associated with M-Index. However, these 4 regions did not show strong evidence of association with skin pigmentation, and we did not find any variant with P< 1.0 x 10−3 (S3 Fig).
Seventeen SNPs with P-value less than 1.0 x 10−5 from 15 genomic regions were selected for replication analysis (Table 2). Two SNPs, rs2675345 near SLC24A5 on chromosome 15, and rs12644472 in the ATP8A1/ GRXCR1 region on chromosome 4, were replicated with P<0.05 and Beta coefficient in the same direction. The SLC24A5 SNP, rs2675345, was strongly associated with M-Index and in the pooled dataset, it had a P-value of 4.0 x 10−30. The ATP8A1/GRXCR1 region SNP, rs12644472, was replicated with P = 0.03. In the pooled dataset, the association between rs12644472 and M-Index did not reach genome-wide significance (P = 2.1x10-5). Three SNPs in the TRHDE gene were selected for the replication study. These SNPs were strongly linked to each other (r2>0.8), and they were not associated with M-Index in our replication dataset. Given strong suggestive associations with M-Index in the discovery cohort and a potential biological significance, we further explored the TRHDE gene by genotyping additional SNPs in our replication cohort, and rs76377291 was the most strongly associated with M-Index (P = 0.009; S4 Fig). Overall, the associations were heterogenous between the GWAS discovery and replication dataset.
We also performed a replication analysis of 24 previously identified variants associated with pigmentation traits in European and admixed populations with our combined GWAS and replication dataset. Fifteen SNPs in 10 regions were replicated with P<0.05 and SNPs from 3 regions reached genome-wide significance (Table 3). A SNP, rs2470102, on SLC24A5 showed the strongest association with M-Index (P = 9.86 x 10−30). The second strongest association was observed for a SLC45A2 variant, rs16891982 (1.93 x 10−13). A SNP, rs1800404 in the OCA2 gene was also strongly associated with M-Index (P = 4.94 x 10−8).
Association of Genetic Score with skin pigmentation and sex interaction
To examine how strongly the identified pigmentation gene variants combined were associated with skin pigmentation, a weighted Genetic Score was calculated with the top SNPs of the 10 most significantly associated loci as well as the top SNPs from 3 loci that reached genome-wide significance (S1 Table).
Table 4 shows the association between the Genetic Scores and M-Index. The Genetic Scores were positively associated with M-Index (P<0.001). Age and sex explained a very small proportion of variance in M-Index (0.3% and 1.7% respectively). WAA explained a large proportion of variance (23.2%). When the Genetic Score was added to the regression model with age, sex, and WAA, the Genetic Score from the top 3 SNPs explained additional 11% of skin pigmentation variance. Using the Genetic Score from the top 10 SNPs to the regression model instead of Genetic Score from 3 SNPs yielded a similar result, and the top 10 SNPs accounted for 12% of variance. In our subsequent analysis, we used the Genetic Score from the top 3 SNPs including rs2470102 (SLC24A5), rs16891982 (SLC45A2), and rs1800404 (OCA2).
Although both the Genetic Score and WAA were strongly associated with M-Index, high variance in M-Index was observed among the study participants with the highest Genetic Score and WAA (S5 Fig). When the association between the Genetic Score and M-Index was examined separately for male and female study participants, a stronger genetic effect was observed in male than female study participants (Table 4 and Fig 2). The interaction between sex and Genetic Score was statistically significant (P = 0.009).
Because we had an unbalanced representation of male and female study participants with a small sample size for females, sex-specific Genetic Scores were calculated to further assess the associations between Genetic Scores and M-Index. The top 3 SNPs included for calculation of sex-specific Genetic Scores were associated with M-Index in both males and females, when sex stratified analysis was performed (S2 Table). The results of sex-specific Genetic Scores and stratified analysis based on sex using overall Genetic Score were similar (S3 Table). The association between sex-specific Genetic Score and M-Index was stronger for males than for females, and males had a larger R2 and β estimate than females.
Association between skin pigmentation gene variants and vitamin D deficiency
Next, we assessed the role of skin pigmentation in serum 25(OH)D levels among male study participants who were part of vitamin D and prostate cancer studies. The M-Index was not associated with serum 25(OH)D levels, but the Genetic Score for skin pigmentation was significantly negatively associated with serum 25(OH)D levels (P = 0.01; S4 Table). The Genetic Score calculated with the top 3 skin pigmentation gene SNPs accounted for a small proportion of serum 25(OH)D variance (0.7%). WAA estimates were not significantly correlated with serum 25(OH)D levels (Spearman’s ρ = 0.015, P = 0.68), and the association was not significant in the linear regression model. When we examined the association between the top 10 SNPs associated with skin pigmentation (S1 Table) and serum 25(OH)D levels, two SNPs, rs2470102 on SLC24A5 (P = 0.004) and rs2733832 on TYRP1 (P = 0.04) were significantly associated with serum 25(OH)D levels (S5 Table).
Genotype data for 38 SNPs in 8 vitamin D metabolic and signaling pathway genes were also available for these individuals, and 3 SNPs were associated with serum 25(OH)D levels in this dataset (rs1155563 in GC, rs12800438 in DHCR7/NADSYN1, and rs11574143 in VDR, S6 Table). A Genetic Score for serum 25(OH)D levels was calculated using these 3 SNPs. The Genetic Score was strongly associated with serum 25(OH)D levels (P<0.001), but these 3 SNPs explained only 1.9% of variance. Genetic Scores using 3 skin pigmentation gene SNPs and 3 vitamin D metabolic and signaling gene SNPs together explained 2.9% of serum 25(OH)D variance, when they were included in the regression model together.
Finally, we examined how strongly the top 3 skin pigmentation gene SNPs are associated with vitamin D deficiency (Table 5). Increasing the number of alleles associated with darker pigmentation (when Genetic Score was treated as an ordinal variable) significantly increased the odds of having severe vitamin D deficiency (OR, 1.34; 95% CI, 1.02–1.86). The odds of severe vitamin D deficiency was even larger when Genetic Score Tertile 1 was compared to Tertile 3 (OR, 1.83; 95% CI, 1.02–3.29). Additional analysis was performed using the male-specific Genetic Score. Overall and male-specific Genetic Scores were highly correlated, and the analysis with male-specific Genetic Score produced identical results with overall Genetic Score.
Discussion
In this first GWAS of skin pigmentation in AAs, we demonstrated that variants near the SLC24A5 gene on chromosome 15 revealed the strongest signal of association in AAs, supporting findings from previous studies in African populations [8,9] and admixed populations [2,7,10]. A variant from this locus along with a variant from both SLC45A2 and OCA2 together accounted for most of the identified genetic variance in M-Index variance (Table 4). However, WAA accounted for a greater proportion of skin pigmentation variance, and we observed great variation in M-Index among AAs with high WAA and genetic score, suggesting that other unknown genomic factors related to WAA are likely contributing to skin pigmentation variation. We also demonstrated that skin pigmentation gene variants were associated with serum 25(OH)D levels in AAs providing evidence to support a role of skin pigmentation in serum 25(OH)D variation.
The top 3 regions that were associated with skin pigmentation were SLC24A5, SLC45A2, and OCA2, and variants in these genes explained a large proportion of skin pigmentation variation in our dataset. The variants in the SLC24A5 are the most strongly associated with skin pigmentation in African descent populations and the associations between them has been demonstrated in AAs [51,52], Cape Verde population [53], and African populations [8,9] as well as meta-analysis in admixed populations [10], but not in European populations [5,54]. The variants in SLC45A2 and OCA2 region also showed very strong associations with skin pigmentation in the African-European admixed population from Cape Verde [53]. Strong associations between variants in these two genes and skin pigmentation or other pigmentation traits have also been shown in European population [1,5,55]. However, African populations appear to have a more complex genetic architecture for skin pigmentation. Crawford and colleagues showed variants in OCA2 and HERC2 region as the fourth most significantly associated locus, but SLC45A2 region was not one of the most significantly associated loci. In the KhoeSan populations from southern Africa, variants in these two regions did not show strong association with skin pigmentation. Instead, these two studies in African populations identified novel loci in SNX13, SMARCA2/VLDLR, TMEM138/DDB1, and MFSD12 associated with skin pigmentation. The associations between MFSD12 variants and skin pigmentation were validated in Latin American populations [10] and meta-analysis of recently admixed populations [10]. These loci may explain additional skin pigmentation variance captured by WAA, but unexplained by the Genetic Score. Because the index SNP from each locus in these studies was not successfully imputed in our GWAS, we looked for SNPs associated with skin pigmentation in these regions. However, variants in these regions were not strongly associated with skin pigmentation. It is likely that variants in these regions contribute to skin pigmentation variation with moderate effect, and a larger sample size is necessary to replicate these findings. Furthermore, SNPs in IRF4, TPCN2, TYR, and MC1R associated with pigmentation traits in European populations were replicated [3–5,55], but these variants contribute a small proportion of skin pigmentation variation in AAs.
We identified a potentially novel locus on the TRHDE gene associated with skin pigmentation. The gene, TRHDE (thyrotropin releasing hormone degrading enzyme), encodes an enzyme that specifically cleaves and inactivates thyrotropin-releasing hormone (TRH). The TRH is produced in the hypothalamus and stimulates production of thyroid-stimulating hormone (TSH). A SNP, rs2044305, in TRHDE has been identified as a candidate variant influencing TSH levels [56]. Thyroid-stimulating hormone as well as melanocyte-stimulating hormone are produced in the anterior pituitary gland. TRH also stimulates growth hormone, prolactin, and α-melanocyte-stimulating hormone in fish, amphibians, and mammals, and these pituitary hormones play important roles in skin [57]. TRH also stimulate melanin synthesis in human hair follicle, potentially by binding to the melanocortin-1 receptor (MC1-R) [58,59]. It is interesting to note that horses with pituitary pars intermedia dysfunction have abnormal coat sometime with lighter color and elevated plasma α-melanocyte-stimulating hormone levels. Administering TRH to the healthy horses increases plasma α-melanocyte-stimulating hormone concentrations [60]. Polymorphisms in the TRHDE gene may alter TRH levels, and subsequently α-melanocyte-stimulating hormone production in the pituitary grand. We selected the top 3 SNPs in the TRHDE region for replication analysis. However, associations of these SNPs were not replicated in the independent dataset. Because of the strong suggestive association in GWAS dataset and a potential biological importance, additional TRHDE SNPs were genotyped to further explore the association with skin pigmentation. The associations of TRHDE SNPs with skin pigmentation were heterogeneous between discovery and replication dataset, but we were successful in replication at a gene level with a different SNP, rs76377291 showing a strong association with skin pigmentation in replication dataset. The SNP rs76377291 is located about 28kb from rs11179301, the lead SNP in the region in our GWAS and these two SNPs were not strongly correlated in 1000 Genome Project ASW population (r2 = 0.05, D’ = 1.0). It is possible that functional variants exist in a nearby location, and we did not capture that genomic variation in our replication data set. Further fine-mapping may help identify the functional variants within this region.
This study explored the associations between skin pigmentation gene variants and serum vitamin D levels in AAs. It has been hypothesized that human skin de-pigmentation evolved as our early ancestors migrated into low UV environment with increased needs for subcutaneous vitamin D synthesis [11] and numerous medical conditions associated with vitamin D deficiency, such as bone disorders, susceptibility to infections, autoimmune disorders, cancer, and reproductive health suggest importance of vitamin D in human health [61]. Despite the important role that skin pigmentation plays for vitamin D status, a few studies examined the relationship between skin pigmentation and vitamin D status in AAs. Studies have shown that individuals with darker skin pigmentation require longer and/or more intense UVR exposure to synthesize sufficient levels of vitamin D [62–64]. Genetic ancestry estimates were also associated with serum 25(OH)D levels in AAs in Southern Community Cohort Study and Black Women’s Health Study [65,66]. These observations led many to believe that disparities in vitamin D status between AAs and EAs is partly due to difference in skin pigmentation [65–67].
Contrary to these studies, genetic ancestry was not associated with serum 25(OH)D levels in our study of AA men. Location of residence, sex, and other factors affect serum 25(OH)D levels, and the differences between our study and previous studies may explain the inconsistent finding. Moreover, there has been no study which explored the association between skin color variation and serum vitamin D levels in AAs using objective measurements. In this study, skin pigmentation gene variants rather than skin pigmentation measured using a reflectometer were associated with serum vitamin D levels. Although skin pigmentation was measured in an area of the body unexposed to the sun, various factors, such as aging, outdoor activities, and consistent UV exposure over the years, may influence skin pigmentation and the association between skin pigmentation and serum vitamin D levels. Because genotype is assigned randomly at conception, the association between skin pigmentation gene variants and serum vitamin D levels is unlikely to be affected by confounding factors [68]. Pigmentation Genetic Scores were also associated with serum 25(OH)D levels in previous studies among children and adult men in the United Kingdom [69,70]. However, the contribution of skin pigmentation genomic variation to serum vitamin D variance was small in our study as well as previous studies. The relationship between skin pigmentation gene variants and serum vitamin D levels should be further examined.
There are some limitations of this study. First, this study had a small sample size for the initial GWAS discovery sample set (n = 395). Thus, we were only able to show association of variants with strong effects. It is likely there are many other yet unknown variants with small effects that influence skin pigmentation in African descent populations. SNPs with P-value less than 1.0 x 10−5 were included in the replication analysis to identify novel variants associated with skin pigmentation with smaller effect size in AAs. One SNP on chromosome 4, rs12644472, was replicated, but did not reach genome-wide significance in the pooled analysis. This SNP is located between two genes, GRXCR1 (glutaredoxin and cysteine rich domain containing 1) and ATP8A1 (ATPase phospholipid transporting 8A1). Mutations in the GRXCR1 genes are linked to hearing loss [71,72]. A GWAS of saggy eyelids showed suggestive associations of ATP8A1 SNPs in European populations, suggesting potential roles of ATP8A1 in skin [73].
Second, a less explored area of pigmentation research is genetic and sex effects on skin pigmentation [74,75]. A previous study found evidence of interaction between ASIP genotypes and sex in AAs [29]. Instead of using genotype, we calculated overall and sex-specific Genetic Scores and demonstrated that the strength of associations between the Genetic Scores and skin pigmentation were different between men and women. Sex hormones may modify genetic effects on skin pigmentation [76]. However, we had an over-representation of men in this study, because many study participants were recruited to prostate cancer risk studies. This observation should be validated incorporating much larger samples of AA women.
Finally, it is likely that the calculated Genetic Scores were over-fitted, because SNPs used for calculation of Genetics Score were selected based on our results instead of using a previously verified set of genetic markers for skin pigmentation and serum vitamin D levels in AA populations. Currently, there is no comprehensive genomic study of skin pigmentation and serum vitamin D levels in AA populations.
In conclusion, our results show strong associations between polymorphisms in 3 major pigmentation genes and skin pigmentation in AAs. The variants in these genes explained a large proportion of skin pigmentation variation and the effects of these variants on skin pigmentation was modified by sex. We also demonstrated that skin pigmentation gene variants were associated with serum vitamin D levels providing support for the vitamin D hypothesis of skin color evolution.
Supporting information
S1 Table [pdf]
10 SNPs used for calculation of Genetic Score in GWAS and replication combined dataset (n = 1,066).
S2 Table [pdf]
Top 3 SNPs associated with M-Index used for calculation of sex-specific Genetic Scores in GWAS and replication combined dataset (n = 1,066).
S3 Table [pdf]
Associations between sex-specific Genetic Scores and M-Index.
S4 Table [pdf]
Associations between Genetic Scores (from 3 pigmentation gene SNPs and 3 vitamin D pathway gene SNPs) and serum vitamin D levels.
S5 Table [pdf]
Association between 10 skin pigmentation SNPs and serum vitamin D levels.
S6 Table [pdf]
Variants in Vitamin D metabolic and signaling pathway gene variants associated with serum 25(OH)D levels (n = 606).
S1 Fig [pdf]
QQ plot indicating there was no evidence of genomic inflation (Inflation factor λ = 1.028) using first three PCs.
S2 Fig [pdf]
Association between SNPs on chromosome 12 and M-Index.
S3 Fig [pdf]
LocusZoom Plots of 4 Genomic Regions Identified in GWAS of African Populations.
S4 Fig [pdf]
Heterogeneous associations between variants and M-Index in replication and GWAS dataset.
Zdroje
1. Eriksson N, Macpherson JM, Tung JY, Hon LS, Naughton B, Saxonov S, et al. Web-based, participant-driven studies yield novel genetic associations for common traits. PLoS Genet. 2010;6:e1000993. doi: 10.1371/journal.pgen.1000993 20585627
2. Stokowski RP, Pant PVK, Dadd T, Fereday A, Hinds DA, Jarman C, et al. A genomewide association study of skin pigmentation in a South Asian population. Am J Hum Genet. 2007;81(6):1119–32. doi: 10.1086/522235 17999355
3. Sulem P, Gudbjartsson DF, Stacey SN, Helgason A, Rafnar T, Jakobsdottir M, et al. Two newly identified genetic determinants of pigmentation in Europeans. Nat Genet. 2008;40(7):835–7. http://www.nature.com/ng/journal/v40/n7/suppinfo/ng.160_S1.html doi: 10.1038/ng.160 18488028
4. Nan H, Kraft P, Qureshi AA, Guo Q, Chen C, Hankinson SE, et al. Genome-Wide Association Study of Tanning Phenotype in a Population of European Ancestry. J Invest Dermatol. 2009;129(9):2250–7. http://www.nature.com/jid/journal/v129/n9/suppinfo/jid200962s1.html doi: 10.1038/jid.2009.62 19340012
5. Han J, Kraft P, Nan H, Guo Q, Chen C, Qureshi A, et al. A genome-wide association study identifies novel alleles associated with hair color and skin pigmentation. PLoS Genet. 2008;4:e1000074. doi: 10.1371/journal.pgen.1000074 18483556
6. Morgan MD, Pairo-Castineira E, Rawlik K, Canela-Xandri O, Rees J, Sims D, et al. Genome-wide study of hair colour in UK Biobank explains most of the SNP heritability. Nat Commun. 2018;9(1):5271. doi: 10.1038/s41467-018-07691-z PubMed Central PMCID: PMC6288091 30531825
7. Adhikari K, Mendoza-Revilla J, Sohail A, Fuentes-Guajardo M, Lampert J, Chacón-Duque JC, et al. A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia. Nat Commun. 2019;10(1):358. doi: 10.1038/s41467-018-08147-0 30664655
8. Crawford NG, Kelly DE, Hansen MEB, Beltrame MH, Fan S, Bowman SL, et al. Loci associated with skin pigmentation identified in African populations. Science. 2017;358(6365). doi: 10.1126/science.aan8433 PubMed Central PMCID: PMC5759959. 29025994
9. Martin AR, Lin M, Granka JM, Myrick JW, Liu X, Sockell A, et al. An unexpectedly complex architecture for skin pigmentation in Africans. Cell. 2017;171(6):1340–53.e14. doi: 10.1016/j.cell.2017.11.015 29195075
10. Lona-Durazo F, Hernandez-Pacheco N, Fan S, Zhang T, Choi J, Kovacs MA, et al. Meta-analysis of GWA studies provides new insights on the genetic architecture of skin pigmentation in recently admixed populations. BMC Genet. 2019;20(1):59–. doi: 10.1186/s12863-019-0765-5 31315583.
11. Jablonski NG, Chaplin G. The colours of humanity: the evolution of pigmentation in the human lineage. Philos Trans R Soc Lond B Biol Sci. 2017;372(1724). doi: 10.1098/rstb.2016.0349 28533464
12. Frost P. Human skin color: a possible relationship between its sexual dimorphism and its social perception. Perspect Biol Med. 1988;32:38–58. doi: 10.1353/pbm.1988.0010 3059317
13. Ihara Y, Aoki K. Sexual selection by male choice in monogamous and polygynous human populations. Theor Popul Biol. 1999;55(1):77–93. doi: 10.1006/tpbi.1998.1388 9925810
14. Aoki K. Sexual selection as a cause of human skin colour variation: Darwin’s hypothesis revisited. Ann Hum Biol. 2002;29:589–608. doi: 10.1080/0301446021000019144 12573076
15. Madrigal L, Kelly W. Human skin-color sexual dimorphism: a test of the sexual selection hypothesis. Am J Phys Anthropol. 2007;132:470–82. doi: 10.1002/ajpa.20453 16685728
16. Bishop DT, Demenais F, Iles MM, Harland M, Taylor JC, Corda E, et al. Genome-wide association study identifies three loci associated with melanoma risk. Nat Genet. 2009;41(8):920–5. http://www.nature.com/ng/journal/v41/n8/suppinfo/ng.411_S1.html doi: 10.1038/ng.411 19578364
17. Duffy DL, Zhao ZZ, Sturm RA, Hayward NK, Martin NG, Montgomery GW. Multiple pigmentation gene polymorphisms account for a substantial proportion of risk of cutaneous malignant melanoma. J Invest Dermatol. 2010;130:520–8. http://www.nature.com/jid/journal/v130/n2/suppinfo/jid2009258s1.html doi: 10.1038/jid.2009.258 19710684
18. Gerstenblith MR, Shi J, Landi MT. Genome-wide association studies of pigmentation and skin cancer: a review and meta-analysis. Pigment Cell Melanoma Res. 2010;23:587–606. doi: 10.1111/j.1755-148X.2010.00730.x 20546537
19. Nan H, Xu M, Kraft P, Qureshi AA, Chen C, Guo Q, et al. Genome-wide association study identifies novel alleles associated with risk of cutaneous basal cell carcinoma and squamous cell carcinoma. Human Molecular Genetics. 2011. doi: 10.1093/hmg/ddr287 21700618
20. Vanchinathan V, Lim HW. A dermatologist’s perspective on vitamin D. Mayo Clin Proc. 2012;87:372–80. doi: 10.1016/j.mayocp.2011.12.010 22425213
21. Liu F, Wollstein A, Hysi PG, Ankra-Badu GA, Spector TD, Park D, et al. Digital Quantification of Human Eye Color Highlights Genetic Association of Three New Loci. PLoS Genet. 2010;6:e1000934. doi: 10.1371/journal.pgen.1000934 20463881
22. Spichenok O, Budimlija ZM, Mitchell AA, Jenny A, Kovacevic L, Marjanovic D, et al. Prediction of eye and skin color in diverse populations using seven SNPs. Forensic Sci Int Genet. 2011;5(5):472–8. doi: 10.1016/j.fsigen.2010.10.005 21050833
23. Goldsmith AH, Hamilton D, Jr WD. From dark to light: skin color and wages among African-Americans. J Hum Resour. 2007;42:701–38.
24. Blair IV, Steiner JF, Havranek EP. Unconscious (implicit) bias and health disparities: Where do we go from here? Perm J. 2011;15:71–8.
25. Krieger N, Sidney S, Coakley E. Racial discerimination and skin color in the CARDIA study: Implications for public health research. Am J Public Health. 1998;88:1308–13. doi: 10.2105/ajph.88.9.1308 9736868
26. Harrison MS, Thomas KM. The hidden prejudice in selection: a research investigation on skin color bias. J Appl Soc Psychol. 2009;39:134–68. doi: 10.1111/j.1559-1816.2008.00433.x
27. Kittles RA, Santos ER, Oji-Njideka N, Bonilla C. Race, skin color and genetic ancestry: Implications for biomedical research on health disparities. Californian Journal of Health Promotion. 2007;5:9–23.
28. Parra EJ, Kittles RA, Shriver MD. Implications of correlations between skin color and genetic ancestry for biomedical research. Nature Genetics. 2004;36:S54–S60. doi: 10.1038/ng1440 15508005
29. Bonilla C, Boxill L-A, Donald SAM, Williams T, Sylvester N, Parra EJ, et al. The 8818G allele of the agouti signaling protein (ASIP) gene is ancestral and is associated with darker skin color in African Americans. Hum Genet. 2005;116(5):402–6. doi: 10.1007/s00439-004-1251-2 15726415
30. Bonilla C, Hooker S, Mason T, Bock CH, Kittles RA. Prostate cancer susceptibility loci identified on chromosome 12 in African Americans. PLoS ONE. 2011;6(2):e16044. doi: 10.1371/journal.pone.0016044 21358824
31. Batai K, Murphy A, Shah E, Ruden M, Newsome J, Agate S, et al. Common vitamin D pathway gene variants reveal contrasting effects on serum vitamin D levels in African Americans and European Americans. Hum Genet. 2014;133(11):1395–405. doi: 10.1007/s00439-014-1472-y PubMed Central PMCID: PMC4185105. 25085266
32. Batai K, Murphy AB, Ruden M, Newsome J, Shah E, Dixon MA, et al. Race and BMI modify associations of calcium and vitamin D intake with prostate cancer. BMC Cancer. 2017;17(1):64. doi: 10.1186/s12885-017-3060-8 PubMed Central PMCID: PMC5248493. 28103838
33. Teteh DK, Dawkins-Moultin L, Hooker S, Hernandez W, Bonilla C, Galloway D, et al. Genetic ancestry, skin color and social attainment: The four cities study. PLoS One. 2020;15(8):e0237041. doi: 10.1371/journal.pone.0237041 32813691.
34. Shriver MD, Parra EJ. Comparison of narrow-band reflectance spectroscopy and tristimulus colorimetry for measurements of skin and hair color in persons of different biological ancestry. Am J Phys Anthropol. 2000;112(1):17–27. doi: 10.1002/(SICI)1096-8644(200005)112:1<17::AID-AJPA3>3.0.CO;2-D 10766940
35. Shriver MD, Parra EJ, Dios S, Bonilla C, Norton H, Jovel C, et al. Skin pigmentation, biogeographical ancestry and admixture mapping. Hum Genet. 2003;112(4):387–99. doi: 10.1007/s00439-002-0896-y 12579416
36. Institute of Medicine. Dietary reference intakes for calcium and vitamin D. Washington, DC: National Academy Press; 2011.
37. Holick MF. Vitamin D deficiency. N Engl J Med. 2007;357:266–81. doi: 10.1056/NEJMra070553 17634462
38. Holick MF, Binkley NC, Bischoff-Ferrari HA, Gordon CM, Hanley DA, Heaney RP, et al. Evaluation, treatment, and prevention of vitamin D deficiency: an Endocrine Society clinical practice guideline. The Journal of clinical endocrinology and metabolism. 2011;96:1911–30. doi: 10.1210/jc.2011-0385 21646368
39. Haiman CA, Chen GK, Blot WJ, Strom SS, Berndt SI, Kittles RA, et al. Genome-wide association study of prostate cancer in men of African ancestry identifies a susceptibility locus at 17q21. Nat Genet. 2011;43(6):570–3. doi: 10.1038/ng.839 PubMed Central PMCID: PMC3102788. 21602798
40. Pibiri F, Kittles RA, Sandler RS, Keku TO, Kupfer SS, Xicola RM, et al. Genetic variation in vitamin D-related genes and risk of colorectal cancer in African Americans. Cancer Causes Control. 2014;25:561–70. doi: 10.1007/s10552-014-0361-y 24562971
41. Ahn J, Yu K, Stolzenberg-Solomon R, Simon KC, McCullough ML, Gallicchio L, et al. Genome-wide association study of circulating vitamin D levels. Hum Mol Genet. 2010;19(13):2739–45. doi: 10.1093/hmg/ddq155 20418485
42. Wang TJ, Zhang F, Richards JB, Kestenbaum B, van Meurs JB, Berry D, et al. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet. 2010;376(9736):180–8. doi: 10.1016/S0140-6736(10)60588-0 20541252
43. Tian C, Hinds DA, Shigeta R, Kittles R, Ballinger DG, Seldin MF. A genomewide single-nucleotide polymorphism panel with high ancestry information for African American admixture mapping. Am J Hum Genet. 2006;79(4):640–9. doi: 10.1086/507954 16960800
44. Delaneau O, Marchini J, Zagury J-F. A linear complexity phasing method for thousands of genomes. Nat Methods. 2012;9(2):179–81. doi: 10.1038/nmeth.1785 22138821
45. Marchini J, Howie B. Genotype imputation for genome-wide association studies. Nat Rev Genet. 2010;11(7):499–511. http://www.nature.com/nrg/journal/v11/n7/suppinfo/nrg2796_S1.html doi: 10.1038/nrg2796 20517342
46. Howie B, Marchini J, Stephens M. Genotype imputation with thousands of genomes. G3 (Bethesda). 2011;1(6):457–70. doi: 10.1534/g3.111.001198 22384356
47. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164(4):1567–87. PubMed Central PMCID: PMC1462648. 12930761
48. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59. 10835412
49. Lambert SA, Abraham G, Inouye M. Towards clinical utility of polygenic risk scores. Hum Mol Genet. 2019;28(R2):R133–r42. Epub 2019/08/01. doi: 10.1093/hmg/ddz187 31363735.
50. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75. doi: 10.1086/519795 PubMed Central PMCID: PMC1950838. 17701901
51. Lamason RL, Mohideen M-APK, Mest JR, Wong AC, Norton HL, Aros MC, et al. SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans. Science. 2005;310(5755):1782–6. doi: 10.1126/science.1116238 16357253
52. Hernandez-Pacheco N, Flores C, Alonso S, Eng C, Mak ACY, Hunstman S, et al. Identification of a novel locus associated with skin colour in African-admixed populations. Sci Rep. 2017;7:44548. doi: 10.1038/srep44548 https://www.nature.com/articles/srep44548#supplementary-information 28300201
53. Beleza S, Johnson NA, Candille SI, Absher DM, Coram MA, Lopes J, et al. Genetic architecture of skin and eye color in an African-European admixed population. PLOS Genetics. 2013;9(3):e1003372. doi: 10.1371/journal.pgen.1003372 23555287
54. Candille SI, Absher DM, Beleza S, Bauchet M, McEvoy B, Garrison NA, et al. Genome-wide association studies of quantitatively measured skin, hair, and eye pigmentation in four European populations. PloS one. 2012;7(10):e48294–e. Epub 2012/10/31. doi: 10.1371/journal.pone.0048294 23118974.
55. Sulem P, Gudbjartsson DF, Stacey SN, Helgason A, Rafnar T, Magnusson KP, et al. Genetic determinants of hair, eye and skin pigmentation in Europeans. Nat Genet. 2007;39(12):1443–52. Epub 2007/10/21. doi: 10.1038/ng.2007.13 17952075.
56. Hwang SJ, Yang Q, Meigs JB, Pearce EN, Fox CS. A genome-wide association for kidney function and endocrine-related traits in the NHLBI’s Framingham Heart Study. BMC medical genetics. 2007;8 Suppl 1(Suppl 1):S10. Epub 2007/10/16. doi: 10.1186/1471-2350-8-S1-S10 17903292; PubMed Central PMCID: PMC1995611.
57. Galas L, Raoult E, Tonon M-C, Okada R, Jenks BG, Castaño JP, et al. TRH acts as a multifunctional hypophysiotropic factor in vertebrates. Gen Comp Endocrinol. 2009;164(1):40–50. doi: 10.1016/j.ygcen.2009.05.003 19435597
58. Gáspár E, Nguyen-Thi KT, Hardenbicker C, Tiede S, Plate C, Bodó E, et al. Thyrotropin-releasing hormone selectively stimulates human hair follicle pigmentation. J Invest Dermatol. 2011;131(12):2368–77. Epub 2011/10/01. doi: 10.1038/jid.2011.221 21956127.
59. van Beek N, Bodó E, Kromminga A, Gáspár E, Meyer K, Zmijewski MA, et al. Thyroid hormones directly alter human hair follicle functions: anagen prolongation and stimulation of both hair matrix keratinocyte proliferation and hair pigmentation. The Journal of clinical endocrinology and metabolism. 2008;93(11):4381–8. Epub 2008/08/30. doi: 10.1210/jc.2008-0283 18728176.
60. Funk RA, Stewart AJ, Wooldridge AA, Kwessi E, Kemppainen RJ, Behrend EN, et al. Seasonal changes in plasma adrenocorticotropic hormone and α-melanocyte-stimulating hormone in response to thyrotropin-releasing hormone in normal, aged horses. J Vet Intern Med. 2011;25(3):579–85. doi: 10.1111/j.1939-1676.2011.0712.x 21457320
61. Yuen AWC, Jablonski NG. Vitamin D: In the evolution of human skin colour. Medical Hypotheses. 2010;74(1):39–44. doi: 10.1016/j.mehy.2009.08.007 19717244
62. Clemens TL, Henderson SL, Adams JS, Holick MF. Increased skin pigment reduces the capacity of skin to synthesise vitamin D3. Lancet. 1982;319:74–6. doi: 10.1016/s0140-6736(82)90214-8 6119494
63. Armas LAG, Dowell S, Akhter M, Duthuluru S, Huerter C, Hollis BW, et al. Ultraviolet-B radiation increases serum 25-hydroxyvitamin D levels: the effect of UVB dose and skin color. J Am Acad Dermatol. 2007;57:588–93. doi: 10.1016/j.jaad.2007.03.004 17637484
64. Xiang F, Lucas R, de Gruijl F, Norval M. A systematic review of the influence of skin pigmentation on changes in the concentrations of vitamin D and 25-hydroxyvitamin D in plasma/serum following experimental UV irradiation. Photochem Photobiol Sci. 2015;14(12):2138–46. doi: 10.1039/c5pp00168d 26548800
65. Signorello LB, Williams SM, Zheng W, Smith JR, Long J, Cai Q, et al. Blood vitamin D levels in relation to genetic estimation of African ancestry. Cancer Epidemiol Biomarkers Prev. 2010;19(9):2325–31. doi: 10.1158/1055-9965.EPI-10-0482 20647395
66. Haddad SA, Ruiz-Narváez EA, Cozier YC, Gerlovin H, Rosenberg L, Palmer JR. Association of degree of European genetic ancestry with serum vitamin D levels in African Americans. Am J Epidemiol. 2018;187(7):1420–3. Epub 2018/02/02. doi: 10.1093/aje/kwy015 29390092; PubMed Central PMCID: PMC6030900.
67. Chan J, Jaceldo-Siegl K, Fraser G. Determinants of serum 25 hydroxyvitamin D levels in a nationwide cohort of blacks and non-Hispanic whites. Cancer Causes Control. 2010;21:501–11. doi: 10.1007/s10552-009-9481-1 PubMed Central PMCID: PMC3427006. 20012182
68. Sheehan NA, Didelez V, Burton PR, Tobin MD. Mendelian randomisation and causal inference in observational epidemiology. PLOS Medicine. 2008;5(8):e177. doi: 10.1371/journal.pmed.0050177 18752343
69. Bonilla C, Ness AR, Wills AK, Lawlor DA, Lewis SJ, Davey Smith G. Skin pigmentation, sun exposure and vitamin D levels in children of the Avon Longitudinal Study of Parents and Children. BMC Public Health. 2014;14:597–. doi: 10.1186/1471-2458-14-597 24924479.
70. Bonilla C, Gilbert R, Kemp JP, Timpson NJ, Evans DM, Donovan JL, et al. Using genetic proxies for lifecourse sun exposure to assess the causal relationship of sun exposure with circulating vitamin d and prostate cancer risk. Cancer Epidemiol Biomarkers Prev. 2013;22(4):597–606. Epub 2013/02/25. doi: 10.1158/1055-9965.EPI-12-1248 23441100.
71. Odeh H, Hunker KL, Belyantseva IA, Azaiez H, Avenarius MR, Zheng L, et al. Mutations in Grxcr1 are the basis for inner ear dysfunction in the pirouette mouse. Am J Hum Genet. 2010;86(2):148–60. Epub 2010/02/04. doi: 10.1016/j.ajhg.2010.01.016 20137774.
72. Schraders M, Lee K, Oostrik J, Huygen PLM, Ali G, Hoefsloot LH, et al. Homozygosity mapping reveals mutations of GRXCR1 as a cause of autosomal-recessive nonsyndromic hearing impairment. Am J Hum Genet. 2010;86(2):138–47. doi: 10.1016/j.ajhg.2009.12.017 20137778
73. Jacobs LC, Liu F, Bleyen I, et al. Intrinsic and extrinsic risk factors for sagging eyelids. JAMA Dermatol. 2014;150(8):836–43. doi: 10.1001/jamadermatol.2014.27 24869959
74. Hernando B, Ibarrola-Villava M, Fernandez LP, Peña-Chilet M, Llorca-Cardeñosa M, Oltra SS, et al. Sex-specific genetic effects associated with pigmentation, sensitivity to sunlight, and melanoma in a population of Spanish origin. Biol Sex Differ. 2016;7:17–. doi: 10.1186/s13293-016-0070-1 26998216.
75. Hernando B, Ibarrola-Villava M, Peña-Chilet M, Alonso S, Ribas G, Martinez-Cadenas C. Sex and MC1R variants in human pigmentation: Differences in tanning ability and sensitivity to sunlight between sexes. J Dermatol Sci. 2016;84(3):346–8. Epub 2016/09/13. doi: 10.1016/j.jdermsci.2016.09.004 27637409.
76. Natale CA, Duperret EK, Zhang J, Sadeghi R, Dahal A, O’Brien KT, et al. Sex steroids regulate skin pigmentation through nonclassical membrane-bound receptors. eLife. 2016;5:e15104. doi: 10.7554/eLife.15104 27115344
77. Graf J, Hodgson R, van Daal A. Single nucleotide polymorphisms in the MATP gene are associated with normal human pigmentation variation. Hum Mutat. 2005;25(3):278–84. doi: 10.1002/humu.20143 15714523.
78. Quillen EE, Bauchet M, Bigham AW, Delgado-Burbano ME, Faust FX, Klimentidis YC, et al. OPRM1 and EGFR contribute to skin pigmentation differences between Indigenous Americans and Europeans. Hum Genet. 2012;131(7):1073–80. Epub 2011/12/24. doi: 10.1007/s00439-011-1135-1 22198722.
79. Kenny EE, Timpson NJ, Sikora M, Yee M-C, Moreno-Estrada A, Eng C, et al. Melanesian blond hair is caused by an amino acid change in TYRP1. Science. 2012;336(6081):554–. doi: 10.1126/science.1217849 22556244.
80. Sturm RA. Molecular genetics of human pigmentation diversity. Hum Mol Genet. 2009;18(R1):R9–R17. doi: 10.1093/hmg/ddp003 19297406
81. Duffy DL, Montgomery GW, Chen W, Zhao ZZ, Le L, James MR, et al. A three-single-nucleotide polymorphism haplotype in intron 1 of OCA2 explains most human eye-color variation. Am J Hum Genet. 2007;80(2):241–52. Epub 2006/12/20. doi: 10.1086/510885 17236130.
82. Liu F, Wollstein A, Hysi PG, Ankra-Badu GA, Spector TD, Park D, et al. Digital quantification of human eye color highlights genetic association of three new loci. PLoS Genet. 2010;6(5):e1000934–e. doi: 10.1371/journal.pgen.1000934 20463881.
Článek vyšel v časopise
PLOS Genetics
2021 Číslo 2
- Může hubnutí souviset s vyšším rizikem nádorových onemocnění?
- Polibek, který mi „vzal nohy“ aneb vzácný výskyt EBV u 70leté ženy – kazuistika
- Zkoušku z bariatrické chirurgie nejlépe složil ChatGPT-4. Za ním zůstaly Bing a Bard
- Raději si zajděte na oční! Jak souvisí citlivost zraku s rozvojem demence?
- Metamizol jako analgetikum první volby: kdy, pro koho, jak a proč?
Nejčtenější v tomto čísle
- Glucocerebrosidase reduces the spread of protein aggregation in a Drosophila melanogaster model of neurodegeneration by regulating proteins trafficked by extracellular vesicles
- ATF3 downmodulates its new targets IFI6 and IFI27 to suppress the growth and migration of tongue squamous cell carcinoma cells
- Transcriptome-wide transmission disequilibrium analysis identifies novel risk genes for autism spectrum disorder
- Four families of folate-independent methionine synthases