Next Article in Journal
Modifier Genes in Microcephaly: A Report on WDR62, CEP63, RAD50 and PCNT Variants Exacerbating Disease Caused by Biallelic Mutations of ASPM and CENPJ
Next Article in Special Issue
Prediction of Eye Colour in Scandinavians Using the EyeColour 11 (EC11) SNP Set
Previous Article in Journal
Can Knowledge of Genetic Distances, Genome Sizes and Chromosome Numbers Support Breeding Programs in Hardy Geraniums?
Previous Article in Special Issue
Animal Forensic Genetics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 5
10.3390/genes12050733
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch

by
Muhammad Adnan Shan
1,2,*,†,
Olivia Strunge Meyer
1,†,
Mie Refn
1,
Niels Morling
1,
Jeppe Dyrberg Andersen
1 and
Claus Børsting
1
1
Section of Forensic Genetics, Department of Forensic Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, 2100 Copenhagen, Denmark
2
Centre for Applied Molecular Biology (CAMB), University of the Punjab, Lahore 54590, Pakistan
*
Author to whom correspondence should be addressed.
†
Authors contributed equally.
Genes 2021, 12(5), 733; https://0-doi-org.brum.beds.ac.uk/10.3390/genes12050733
Submission received: 10 April 2021 / Revised: 29 April 2021 / Accepted: 11 May 2021 / Published: 13 May 2021
(This article belongs to the Special Issue Advances in Forensic Genetics)
Browse Figures
Versions Notes

Abstract

:
Skin pigmentation is one of the most prominent and variable phenotypes in humans. We compared the alleles of 163 SNPs and indels from the Human Pigmentation (HuPi) AmpliSeq™ Custom panel, and biogeographic ancestry with the quantitative skin pigmentation levels on the upper arm, lower arm, and forehead of 299 Pakistani individuals from three subpopulations: Baloch, Pashtun, and Punjabi. The biogeographic ancestry of each individual was estimated using the Precision ID Ancestry Panel. All individuals were mainly of mixed South-Central Asian and European ancestry. However, the Baloch individuals also had an average proportion of Sub-Saharan African ancestry of approximately 10%, whereas it was <1% in the Punjabi and Pashtun individuals. The pairwise genetic distances between the Pashtun, Punjabi, and Baloch subpopulations based on the ancestry markers were statistically significantly different. Individuals from the Pashtun subpopulation had statistically significantly lower skin pigmentation than individuals from the Punjabi and Baloch subpopulations (p < 0.05). The proportions of European and Sub-Saharan African ancestry and five SNPs (rs1042602, rs10831496, rs1426654, rs16891982, and rs12913832) were statistically significantly associated with skin pigmentation at either the upper arm, lower arm or forehead in the Pakistani population after correction for multiple testing (p < 10−3). A model based on four of these SNPs (rs1426654, rs1042602, rs16891982, and rs12913832) explained 33% of the upper arm skin pigmentation. The four SNPs and the proportions of European and Sub-Saharan African ancestry explained 37% of the upper arm skin pigmentation. Our results indicate that the four likely causative SNPs, rs1426654, rs1042602, rs16891982, and rs12913832 located in SLC24A5, TYR, SLC45A2, and HERC2, respectively, are essential for skin color variation in the admixed Pakistani subpopulations.

1. Introduction

Human pigmentation is one of the most variable externally visible characteristics (EVCs). Prediction of EVCs may result in helpful leads early in a police investigation by providing a ‘genetic eye witness’ of a possible perpetrator and may allow the police investigators to focus their attention on specific groups of individuals and decrease the number of potential suspects [1,2,3]. The genetics of human pigmentation is a field of great interest in forensic genetics because it offers a prediction of three prominent EVCs: hair, eye, and skin color [4,5]. Skin pigmentation is strongly associated with biogeographic ancestry [6] and influenced by environmental factors (e.g., UV radiation) [7], which complicates the identification of causative variants and the reliability of the predictions. Prediction of the biogeographic ancestry may indicate skin pigmentation of the individual. However, skin color prediction based on the assumption that biogeographic ancestry and skin color are well correlated only offers an indirect prediction, and the assumption is not necessarily true in admixed populations [8,9]. Studies of admixed individuals from North and South American populations showed that the correlation between ancestry and skin pigmentation varied considerably and explained 20–65% of the skin pigmentation [3,6,10,11]. Nevertheless, admixed populations are ideal for studying the effect of genetic variants associated with skin pigmentation because they have a wide range of skin colors, and the ancestry proportions of each individual may be estimated using standard population genetic methods.
The Pakistani population is an admixed population with several subpopulations and ethnic groups. Pakistan is situated at the junction of the Middle East, Central Asia, and Southeast Asia. This area was the home of various ancient cultures, including the Bronze Age Indus Valley Civilization [12]. The region was ruled by dynasties and empires of different religions and cultures, including Muslims, Turco-Mongols, Hindus, Indo-Greeks, Sikhs, and Afghans [13]. Hence, the Pakistani population is well suited for studying biogeographic ancestry and skin pigmentation genetics.
We studied the associations among pigmentary variants, biogeographic ancestry, and quantitative skin pigmentation in three admixed subpopulations from Pakistan (n = 299). We evaluated the individuals’ ancestry with the Precision ID Ancestry Panel (Thermo Fisher Scientific, Waltham, MA, USA) [14,15,16] and investigated the association between skin pigmentation and genetic pigmentary variants in the HuPi AmpliSeq™ Custom panel [17]. We also analyzed the correlation between skin pigmentation and the 36 HIrisPlex-S SNPs [18,19] as well as nine SNPs found to be associated with skin pigmentation variation in an admixed Brazilian population [3]. Lastly, we analyzed the efficiency of the skin color prediction using the HIrisPlex-S prediction model [18,19].

2. Materials and Methods

2.1. Samples, DNA Extraction, and DNA Quantification

Buccal swabs were collected from 299 unrelated healthy Pakistani individuals and stored on FTA cards. The individuals belonged to three major Pakistani subpopulations: 107 Baloch from Baluchistan, 103 Pashtuns from Khyber Pakhtunkhwa, and 89 Punjabis from northern Punjab. Signed informed consent was obtained from all participants, and the samples were anonymized. The study was approved by the Review Board/Ethical Committee of the University of the Punjab, Pakistan (D/No. 019/DFEMS). Genomic DNA was extracted using the Qiagen BioRobot® EZ1 Workstation and the EZ1 DNA investigator kit (Qiagen, Hilden, Germany). Purified DNA was quantified using the Qubit™ dsDNA HS Assay Kit and a Qubit® 4.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s recommendations.

2.2. Measurement of Quantitative Skin Pigmentation

The UV-Optimize Scientific 555 (Chromo Light APS, Vedbæk, Denmark) calibrated with a white standard (ISO 2469) was used to measure quantitative skin pigmentation. Quantitative skin pigmentation was determined as the pigment protection factor, PPF, which is a measure for skin pigmentation levels [20,21]. The skin pigmentation of each individual was measured in triplicates from skin areas on the lower inner forearm (lower arm), upper inner arm (upper arm), and forehead. The median PPF-value for each measuring site was used for further analysis. All skin measurements were taken from areas without hair, freckles, nevi, or tattoos.

2.3. Typing with the Precision ID Ancestry Panel and the Human Pigmentation (HuPi) AmpliSeq™ Custom Panel

DNA was amplified with the Precision ID Ancestry Panel (Thermo Fisher Scientific, Waltham, MA, USA) and the Human Pigmentation (HuPi) AmpliSeq™ Custom panel [17] in two separate reactions. The Precision ID Ancestry Panel targets 165 ancestry-informative markers (AIMs). The HuPi AmpliSeq™ Custom panel targets 183 SNPs and indels, which were previously found to be associated with human pigmentary variation [17]. Sequencing libraries were prepared with the Ion AmpliSeq™ Library Kit 2.0 (Thermo Fisher Scientific, Waltham, MA, USA) using half volume of all reagents but otherwise following the manufacturer’s instructions. Amplification of targets was carried out using approximately 1 ng DNA (1–3 µL) and 24 and 26 cycles for the Precision ID Ancestry Panel and the HuPi AmpliSeq™ Custom panel, respectively. Sequencing libraries were purified using the Biomek 3000 pipetting robot [22]. Purified libraries were quantified using the Qubit™ 4.0 Fluorometer with the Qubit™ dsDNA HS Assay kit (Thermo Fisher Scientific, Waltham, MA, USA) and pooled in equimolar concentrations to a final volume of 28 pM for the Precision ID Ancestry panel and 35–60 pM for the HuPi AmpliSeq™ Custom Panel. Template preparation was performed with the Ion Chef™ instrument (Thermo Fisher Scientific, Waltham, MA, USA) using the Ion S5™ Precision ID Chef Kit (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer’s recommendations (Thermo Fisher Scientific, Waltham, MA, USA). Sequencing was carried out with Ion 530™ Chips and Ion S5™ Precision ID Sequencing Reagents on the Ion S5™ (Thermo Fisher Scientific, Waltham, MA, USA) with up to 80 samples per chip.

2.4. Analysis of Sequencing Data

The data were initially analyzed with the Torrent Suite Software v.5.10.1 (Thermo Fisher Scientific, Waltham, MA, USA). For the Precision ID Ancestry panel, data were analyzed with the HID_SNP_Genotyper_v5_2_2 plugin using GRCh37 (hg19) as the reference genome and default parameters (Thermo Fisher Scientific, Waltham, MA, USA). The data obtained with the HuPi AmpliSeq™ Custom Panel were analyzed with the VariantCaller v5.10.1.20 plugin and the GRCh38.p2 (hg38) reference genome. The “Down-sample to coverage” was changed from 400 to 10,000 reads. Otherwise, the analyses were performed with the default parameters (Thermo Fisher Scientific, Waltham, MA, USA). The resulting Excel-files were processed using R v.4.0.2 (R core team, https://www.R-project.org/, accessed on 1 November 2020) and an in-house developed script. The data quality control comprised evaluation of the locus balance, heterozygote balance (Hb), and noise levels for each target. Hb was calculated as the number of reads for one nucleotide divided by the number of reads for the other nucleotide in the order A, C, G, and T. Genotypes were accepted when a locus had ≥45 reads and the Hb was between 0.3 and 3. SNPs with 30–44 reads were inspected and accepted if the Hb was between 0.7 and 1.5 or if the noise did not exceed one read for homozygous genotype calls. Genotypes that did not meet the criteria were annotated as NN. For both panels, loci with NNs in more than 15% of the samples were excluded. Moreover, samples with more than 15% missing data (NN) were excluded from further analysis. Hardy–Weinberg Equilibrium (HWE) and pairwise linkage disequilibrium (LD) were calculated with Haploview v. 4.2 [23], and pairwise FST-values were calculated using Arlequin v3.5.2.2 software [24].

2.5. Classification of Biogeographic Ancestry

The alleles of the AIM-SNPs in each sample were investigated using STRUCTURE v.2.3.4.21 [25,26]. The analyses were carried out using 100,000 steps of burn-in followed by 100,000 repetitions for the MCMC. The ‘admixture’ and the ‘correlated allele frequencies’ models were used [25,26]. After testing the number of clusters (K) from K = 3 to K = 6, K was set to 4, corresponding to the Sub-Saharan African, South-Central Asian, European, and East Asian populations. The results were visualized using CLUMPP v.1.1.222 [27] and Distruct v.1.1.23 [28]. Principal component analysis (PCA) was carried out using an in-house written Python script. Reference population data (Supplementary Table S1) were collected as previously described [14,29]. Four SNPs, rs1800414, rs12913832, rs1426654, and rs16891982, were included in both the Precision ID Ancestry Panel and the HuPi AmpliSeq™ Custom panel. One SNP, rs10954737, lacked genotype data in all reference populations. Hence, these five SNPs were not considered in the PCA and STRUCTURE analyses and were not used to classify ancestry.

2.6. Correlations between Biogeographic Ancestry, Pigmentary Variants, and Skin Pigmentation

We excluded monomorphic genetic variants (i.e., variants with only one observed genotype) and variants in LD (r2 > 0.8) with other variants in the HuPi AmpliSeq™ Custom Panel. Hence, the dataset was limited to 102 genetic variants and four metapopulations (Sub-Saharan Africa, South-Central Asia, Europe, and East Asia). Differences in skin pigmentation levels of the lower arm, upper arm, and forehead, as well as differences in skin pigmentation of the three subpopulations, were investigated using Student’s t-test for paired data or the Welch two-sample t-test for unpaired data. The correlations among skin pigmentation levels, genetic variants in the HuPi AmpliSeq™ Custom Panel, and the estimation of the ancestry proportions were carried out using multiple linear regression (MLR). The correlation was evaluated using the adjusted R2. To select a subset of variants, backwards model selection was performed using MLR with Akaike Information Criterion (AIC) using the stepAIC function in the MASS R-package [30].
We also evaluated the correlation between skin pigmentation and the genotypes of the 36 HIrisPlex-S SNPs [18,19], and between skin pigmentation and nine SNPs (rs1426654, rs1448484, rs16891982, rs4424881, rs10831496, rs6119471, rs12913832, rs10424065, and rs1408799) previously found to explain up to 65% of the skin pigmentation variation in a Brazilian population of primarily European ancestry [3].

2.7. Prediction of Skin Colour Using the HIrisPlex-S

Prediction of skin color using the HIrisPlex-S model was carried out with the online web-tool (https://hirisplex.erasmusmc.nl/, accessed on 1 Novmber 2020). The HIrisPlex-S model predicts skin color in five categories: Very pale, Pale, Intermediate, Dark, and Dark to Black [18,19]. The skin color category with the highest predictive value was used as the skin pigmentation prediction. The HIrisPlex-S skin color prediction and the quantitative skin pigmentation measurements were compared.

2.8. Statistical Methods

The statistical methods used are described in each section. When multiple comparisons were performed, the statistical significance was corrected with the Bonferroni method.

3. Results

3.1. Skin Pigmentation Measurements

Quantitative skin pigmentation measurements were performed for the 299 Pakistani individuals, Baloch (n = 107), Pashtun (n = 103), and Punjabi (n = 89). The lowest pigmentation levels were observed on the upper arm (mean PPF of 11.08), followed by the lower arm (mean PPF of 12.30) and the forehead (mean PPF of 12.43). We observed statistically significant differences between the pigmentation levels on the upper arm and lower arm (p < 0.05) and the upper arm and forehead (p < 0.05) among all subpopulations. Statistically significant difference between pigmentation levels of the lower arm and forehead was only observed in the Baloch subpopulation (p < 0.05).
Individuals from the Pashtun subpopulation showed statistically significant (p < 0.05) lower pigmentation levels than individuals in the Baloch and Punjabi subpopulations on all measured areas (Figure 1).

3.2. The Precision ID Ancestry Panel

The allele frequencies for the 165 loci in the Punjabi, Pashtun, and Baloch populations are presented in Supplementary Table S2. The median number of reads per target was 493 (range: 45–12,105), the median Hb was 1.0 (range: 0.3–3.0) for heterozygous SNPs, and the median level of noise was 0.0% (range: 0–14.3%). One locus, rs2070586, deviated from HWE (p < 10−3) in the Baloch population. The loci rs3811801 and rs671 were monomorphic in all three populations, rs1871534 and rs2814778 were monomorphic in Pashtuns and Punjabis, and rs1800414 was monomorphic in the Baloch. No pairwise LD (pairwise r2 > 0.8) was detected in any of the three populations. The observed values of pairwise FST genetic distances were statistically significantly different among all subpopulations (p < 10−5) (Supplementary Table S3).

3.3. Proportion of Ancestry Components

The biogeographic ancestry of each individual was investigated with STRUCTURE and PCA using reference data from 34 populations grouped into six metapopulations: Sub-Saharan African, North African, European, Middle Eastern, South-Central Asian, and East Asian (Supplementary Table S1). The results of the STRUCTURE analysis with K = 4 are shown in Figure 2. The Baloch, Punjabi, and Pashtun subpopulations were admixed populations with major contributions from South-Central Asian and European populations. The Baloch differed from the other two subpopulations due to approximately 10% Sub-Saharan African genetic contribution, whereas the Sub-Saharan African contribution was <1% in the Punjabis and Pashtuns. The Baloch had the lowest proportion of South-Central Asian ancestry (44.7%), followed by the Pashtuns (56.6%) and the Punjabis (69.3%).
The first two principal components (PC1 and PC2) in the PCA analysis separated the Sub-Saharan African, European, and East Asian populations from each other (Supplementary Figure S1). The North Africans clustered closest to the European and Middle Eastern populations, while the South-Central Asian populations clustered between the East Asian and European populations. The positions of the three Pakistani subpopulations, Punjabis, Pashtuns, and the Baloch, overlapped with each other and the South-Central Asian and Middle Eastern populations. The Baloch individuals clustered closer to the African populations than the Pashtuns and Punjabis. In contrast, the Punjabis clustered closer to the South-Central Asian populations than the Pashtuns and Baloch, supporting the STRUCTURE analysis conclusions.

3.4. The HuPi AmpliSeq™ Custom Panel

A total of 163 variants were successfully typed with the HuPi AmpliSeq™ Custom Panel (Supplementary Table S4). The median number of reads was 597 (range: 60–1532) per target. The median Hb for heterozygous SNP calls was 1.0 (range: 0.35–2.86). The median level of noise was 0.0% (range: 0–12.7%). The allele frequencies of rs4778241 in OCA2 deviated statistically significantly from HWE (p < 10−3). Twelve variants were monomorphic (Supplementary Table S4), and these data were removed from the subsequent analysis. For the correlation measurements between skin pigmentation and genetic variants, only independent genetic variants were considered (pairwise r2 < 0.8), reducing the number of genetic variants to 102 (Supplementary Table S4).

3.5. Correlation between Skin Pigmentation and Biogeographic Ancestry

The proportions of Sub-Saharan African, European, South-Central Asian, and East Asian ancestry were compared with the skin pigmentation at each measured site using MLR. We found statistically significant (p < 10−3) adjusted R2-values for the upper arm: 0.10, lower arm: 0.13, and forehead: 0.12 when comparing the full ancestry profiles with the skin pigmentation (Supplementary Table S5). We found statistically significant (p < 0.05) negative correlations between the proportion of European ancestry and skin pigmentation at all three measured sites (adjusted R2: 0.04–0.09) and statistically significant positive correlation (p < 0.05) between Sub-Saharan African ancestry and skin pigmentation at the upper arm and forehead (adjusted R2 = 0.06 and 0.07, respectively). In contrast, the proportions of East Asian and South-Central Asian ancestry showed no statistically significant correlation with skin pigmentation (Supplementary Table S5). We analyzed the correlation between the proportions of only European and Sub-Saharan African ancestry with skin pigmentation at all three measuring sites using MLR. The resulting adjusted R2-values were 0.10 (p < 10−3), 0.13 (p < 10−6), and 0.12 (p < 10−6) for the upper arm, lower arm, and forehead pigmentation, respectively (Supplementary Table S5).

3.6. Correlations among Skin Pigmentation, Pigmentary Variants, and Biogeographic Ancestry

Each of the independent variants typed with the HuPi AmpliSeq™ Custom Panel was compared with the skin pigmentation at the upper arm, lower arm, and forehead using linear regression. We observed statistically significant (p < 0.05) correlations between skin pigmentation on either the upper arm, lower arm, or forehead and the SNPs rs1042602 in TYR and rs10831496 in GRM5, rs1426654 in SLC24A5, rs16891982 in SLC45A2, and rs12913832 in HERC2 (Table 1). Only rs1426654 and rs1042602 were statistically significantly correlated with skin pigmentation at all three sites (p < 10−3) (Table 1).
The five SNPs were compared with skin pigmentation using MLR, and the adjusted R2 ranged from 0.24 to 0.33 (p < 10−6) (Table 2). Subsequently, model selection was performed using stepwise AIC, giving the best model with rs1426654, rs1042602, rs16891982, and rs12913832 for all three sites (upper arm, lower arm, and forehead). When the proportions of European and Sub-Saharan African ancestry were included, the adjusted R2-values increased slightly (Table 2, Figure 3). The cumulative correlations between (1) the four SNPs and ancestry and (2) the skin pigmentation are shown in Figure 3. The four SNPs, ancestry, and the upper arm pigmentation showed the highest adjusted R2-values (Table 2 and Figure 3).

3.7. Correlation of Skin Pigmentation with SNPs from Existing Skin Colour Models

The skin pigmentation in the Pakistani individuals was compared with nine SNPs previously found to be associated with skin pigmentation in a Brazilian population of primarily European ancestry [3]. The nine SNPs included rs10831496, rs12913832, rs1426654, and rs16891982, which were associated with skin pigmentation in the Pakistani population (Table 1), as well as rs1448484 in OCA2, rs4424881 in APBA2, rs6119471 in ASIP, rs1408799 in TYRP1, and rs10424065 in MFSD12, which were not associated with skin pigmentation in our study. In the Pakistani population, these nine SNPs explained up to 26% of the skin pigmentation (upper arm). The adjusted R2-values were 0.26 (p < 0.05), 0.15 (p < 0.05), and 0.19 (p < 0.05) for the upper arm, lower arm, and forehead pigmentation, respectively.
We also correlated the HIrisPlex-S SNPs with the skin pigmentation on the upper arm, lower arm, and forehead using MLR. Four SNPs in the MC1R-gene, rs3212355, rs1805006, rs11547464, and rs1110400 were monomorphic and excluded from the analysis. Thus, the correlation was based on 32 of the 36 skin color predictive SNPs in the HIrisPlex-S. The adjusted R2-values were 0.36 (p < 10−6), 0.29 (p < 10−6), and 0.27 (p < 10−6), for pigmentation on the upper arm, lower arm, and forehead, respectively.
Lastly, we evaluated the performance of the HIrisPlex-S model that is a forensically validated skin color prediction model [18,19]. The predicted skin colors of the 299 Pakistani individuals were: Intermediate: 102, Dark: 167, and Black: 30. No individual was predicted to have very pale or pale skin colors. Comparisons between the predicted skin color categories and the quantitative skin pigmentation measurements are shown in Figure 4.

4. Discussion

The pairwise genetic distances calculated among the three subpopulations were statistically significantly different (Supplementary Table S3). The Baloch were more distant from Punjabis and Pashtuns than Punjabis were from Pashtuns. The results agreed with our previous results obtained with short tandem repeats that depicted genetic differences between the Baloch and other subpopulations from Pakistan [31]. The STRUCTURE analysis (Figure 2) showed that all three subpopulations from Pakistan were admixed, mainly with South-Central Asian and European genetic contributions. However, the Baloch population also had a genetic contribution from Sub-Saharan African populations, which may be remnants of African individuals settling in the Indian subcontinent [32,33].
The proportion of European and Sub-Saharan African ancestry affected the skin pigmentation, whereas the proportion of South-Central Asian ancestry did not (Supplementary Table S5). However, the proportions of European and Sub-Saharan African ancestry only explained approximately 10% of the variation (Supplementary Table S5) and only influenced the adjusted R2-values of the MLR model marginally (Table 2). This showed that ancestry was a relatively poor predictor of skin pigmentation in the highly admixed Pakistani populations. It also indicated that any association between a locus and skin pigmentation levels would most likely come from causative SNPs or indels, and not be a consequence of population genetic differences.
Of the 163 variants that were previously shown to be associated with pigmentary traits, five SNPs were statistically significantly correlated with the skin pigmentation levels of the investigated individuals from the Pakistani subpopulations (Table 1). These five SNPs are located in five genes, TYR, GRM5, SLC24A5, SLC45A2, and HERC2. Four of the five SNPs, rs1042602 in TYR, rs1426654 in SLC24A5, rs16891982 in SLC45A2, and rs12913832 in HERC2, explained 33% of the skin pigmentation variation on the upper arm area and 23–24% of the skin pigmentation variation on the lower arm and forehead that are less protected from UV irradiation (Table 2). Three of the SNPs, rs1426654 (ranked 1), rs1042602 (ranked 2), and rs16891982 (ranked 3), were previously found to be associated with skin pigmentation in individuals of South Asian descent [7,34,35]. In agreement with our results, the SLC24A5 SNP rs1426654 had the most pronounced effect on skin pigmentation levels [7,34,36]. The HERC2 SNP rs12913832 (ranked 4) has not previously reported to be associated with skin pigmentation in South Asian populations. However, it was reported to be associated with skin color in other populations [11,19]. The fifth SNP, GRM5 rs10831496 (ranked 5), did not increase the cumulative adjusted R2-value of the MLR model at any of the measured pigmentation areas. This SNP is highly polymorphic and located in intron 3 of the GRM5 gene, 340 kbp upstream of the TYR gene on chromosome 11. GRM5 rs10831496 was previously reported to be associated with tanning response after exposure to sunlight [37]. In our study, it was associated with skin pigmentation on the lower arm.
The three highest ranking SNPs, rs1426654, rs1042602, and rs16891982, are non-synonymous variants: Thr111Ala in SLC24A5, Ser192Tyr in TYR, and Phe374Leu in SLC45A2, respectively, and they all influence melanin production. In human primary melanocytes that were homozygous for the SLC24A5 111Ala variant, the TYR activity, melanin content, and the number of mature melanosomes were decreased [38]. The 192Tyr TYR variant reduced the TYR activity by 40% [39] and the number of primary melanocytes, and human skin tissue samples that were homozygous for the 192Tyr variant had reduced TYR activity and expression [40]. The SLC45A2 374Leu variant was found in increased frequency in non-European populations, and 374Leu homozygous melanocytes had higher TYR activity and melanin content than melanocytes that were homozygous for the 374Phe variant [38]. The fourth-ranking SNP, HERC2 rs12913832, is the most important variant for eye color determination [41,42,43]. It is positioned in a key enhancer element of the OCA2 pigmentary gene [44]. Human primary melanocytes that were homozygous for the derived rs12913832 G allele had reduced TYR activity, melanin content, and numbers of mature melanosomes [38]. When we analyzed five additional SNPs that were previously found to be associated with pigmentation in an admixed Brazilian population with primarily European and Sub-Saharan African genetic contributions [3], the correlation with skin pigmentation of the upper arm was reduced from 33% to 26%. This indicates that the biogeographic background of an individual is important for the selection of markers that eventually will give the best prediction of skin color, which is well in line with the knowledge that some genetic markers are associated with skin pigmentation in certain populations and not in others [45]. It also indicates that the estimation of the biogeographic background should be the first step in the estimation of skin pigmentation levels from DNA samples of unknown individuals. The second step should be to investigate causative variants found to be associated with skin pigmentation in the most likely biogeographic background of the collected trace sample and these variants should form the basis for the estimation of skin pigmentation.
In forensics genetic casework, ancestry inference can provide helpful information about an unidentified perpetrator who left a biological trace sample at a crime scene. Several well-established AIM panels that successfully differentiate the major human populations, including Europeans, East Asians, South Central Asians, Native Americans, Sub-Saharan Africans, and Oceanians, have been developed [14,46,47,48,49]. Ancestry information may also provide indications on EVCs that are typical for a given population. However, caution must be taken when using ancestry to infer the phenotypes of an individual. An individual of mixed ancestry may have AIM alleles inherited from one parental population, but phenotypic characteristics, e.g., skin color, of the other parental population [3,9]. Thus, skin color and other EVCs should be estimated using causative genetic variants. We identified four likely causative variants that influence skin pigmentation in Pakistanis and are included in the leading skin color prediction model, HIrisPlex-S [18,19]. However, the variants only explained approximately one-third of the skin pigmentation variation and the remaining SNPs of the HIrisPlex-S only slightly increased the adjusted R2-values underlining the need for additional genetic studies of human skin pigmentation.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes12050733/s1, Figure S1: PCA plot of the studied populations and selected reference populations, Table S1: Reference data used in the population comparison analysis, Table S2: Allele frequencies of 165 variants typed with the Precision ID Ancestry panel in 299 Pakistani individuals from three subpopulations (Baloch, Pashtun, and Punjabi), Table S3: Pairwise FST-values among the three Pakistani subpopulations Baloch, Pashtun, and Punjabi based on the allele frequencies of 165 AIMs typed with the Precision ID Ancestry panel, Table S4: Allele frequencies of 163 variants typed with the Human Pigmentation (HuPi) AmpliSeq™ Custom panel in 299 Pakistani individuals from three subpopulations (Baloch, Pashtun, and Punjabi), Table S5: Estimated ancestry proportions associated with skin pigmentation on the upper arm, lower arm, or forehead in 299 Pakistani individuals from three subpopulations (Baloch, Pashtun, and Punjabi).

Author Contributions

Conceptualization, C.B., J.D.A., M.A.S., N.M., and O.S.M.; analysis and interpretation, M.A.S., M.R., and O.S.M.; drafting the work, M.A.S. and O.S.M.; critically revising the work for important intellectual content, C.B., J.D.A., M.R., and N.M. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by an overseas research grant to Muhammad Adnan Shan from the University of the Punjab, Pakistan (No. D-1829-Est-I/2017).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Review Board/Ethical Committee of the University of the Punjab, Pakistan (D/No. 019/DFEMS). The project was notified to the Ethics Committee for the Capital Region of Denmark (Journal no. H-20024636). According to the Danish Act on Research Ethics Review of Health Research Projects, the work did not require approval by the Ethics Committee.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data generated in the present study are included within the manuscript and its supplementary file.

Acknowledgments

The authors would like to thank Nadia Jochumsen for laboratory assistance and Vania Pereira for help with the population genetic analyses. We thank Kenneth Kidd and his staff for providing the genotype SNP data of the reference populations assembled from publicly available data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kayser, M.; Schneider, P.M. DNA-based prediction of human externally visible characteristics in forensics: Motivations, scientific challenges, and ethical considerations. Forensic Sci. Int. Genet. 2009, 3, 154–161. [Google Scholar] [CrossRef] [PubMed]
  2. Samuel, G.; Prainsack, B. Forensic DNA phenotyping in Europe: Views “on the ground” from those who have a professional stake in the technology. New Genet. Soc. 2019, 38, 119–141. [Google Scholar] [CrossRef] [Green Version]
  3. Andersen, J.D.; Meyer, O.S.; Simão, F.; Jannuzzi, J.; Carvalho, E.; Andersen, M.M.; Pereira, V.; Børsting, C.; Morling, N.; Gusmão, L. Skin pigmentation and genetic variants in an admixed Brazilian population of primarily European ancestry. Int. J. Legal Med. 2020, 134, 1569–1579. [Google Scholar] [CrossRef]
  4. Kayser, M.; De Knijff, P. Improving human forensics through advances in genetics, genomics and molecular biology. Nat. Rev. Genet. 2011, 12, 179–192. [Google Scholar] [CrossRef] [PubMed]
  5. Maroñas, O.; Phillips, C.; Söchtig, J.; Gomez-Tato, A.; Cruz, R.; Alvarez-Dios, J.; de Cal, M.C.; Ruiz, Y.; Fondevila, M.; Carracedo, Á.; et al. Development of a forensic skin colour predictive test. Forensic Sci. Int. Genet. 2014, 13, 34–44. [Google Scholar] [CrossRef]
  6. Parra, E.J.; Kittles, R.A.; Shriver, M.D. Implications of correlations between skin color and genetic ancestry for biomedical research. Nat. Genet. 2004, 36, S54–S60. [Google Scholar] [CrossRef] [PubMed]
  7. Stokowski, R.P.; Pant, P.V.K.; Dadd, T.; Fereday, A.; Hinds, D.A.; Jarman, C.; Filsell, W.; Ginger, R.S.; Green, M.R.; Van Der Ouderaa, F.J.; et al. A genomewide association study of skin pigmentation in a South Asian population. Am. J. Hum. Genet. 2007, 81, 1119–1132. [Google Scholar] [CrossRef] [Green Version]
  8. Shriver, M.D.; Parra, E.J.; Dios, S.; Bonilla, C.; Norton, H.; Jovel, C.; Pfaff, C.; Jones, C.; Massac, A.; Cameron, N.; et al. Skin pigmentation, biogeographical ancestry and admixture mapping. Hum. Genet. 2003, 112, 387–399. [Google Scholar] [CrossRef] [PubMed]
  9. Valenzuela, R.K.; Henderson, M.S.; Walsh, M.H.; Garrison, N.; Kelch, J.T.; Cohen-Barak, O.; Erickson, D.T.; John Meaney, F.; Bruce Walsh, J.; Cheng, K.C. Predicting phenotype from genotype: Normal pigmentation. J. Forensic Sci. 2010, 55, 315–322. [Google Scholar] [CrossRef] [Green Version]
  10. Beleza, S.; Johnson, N.A.; Candille, S.I.; Absher, D.M.; Coram, M.A.; Lopes, J.; Campos, J.; Araújo, I.I.; Anderson, T.M.; Vilhjálmsson, B.J.; et al. Genetic Architecture of Skin and Eye Color in an African-European Admixed Population. PLOS Genet. 2013, 9, e1003372. [Google Scholar] [CrossRef] [Green Version]
  11. Lona-Durazo, F.; Hernandez-Pacheco, N.; Fan, S.; Zhang, T.; Choi, J.; Kovacs, M.A.; Loftus, S.K.; Le, P.; Edwards, M.; Fortes-Lima, C.A.; 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, 59. [Google Scholar] [CrossRef] [Green Version]
  12. Wright, R.P. The Ancient Indus: Urbanism, Economy, and Society; Case Studies in Early Societies; Cambridge University Press: Cambridge, UK, 2010; ISBN 9780521576529. [Google Scholar]
  13. Wynbrandt, J. A Brief History of Pakistan; Brief history; Infobase Publishing: New York, NY, USA, 2009; ISBN 9780816061846. [Google Scholar]
  14. Pereira, V.; Mogensen, H.S.; Børsting, C.; Morling, N. Evaluation of the Precision ID Ancestry Panel for crime case work: A SNP typing assay developed for typing of 165 ancestral informative markers. Forensic Sci. Int. Genet. 2017, 28, 138–145. [Google Scholar] [CrossRef] [PubMed]
  15. Mogensen, H.S.; Tvedebrink, T.; Børsting, C.; Pereira, V.; Morling, N. Ancestry prediction efficiency of the software GenoGeographer using a z-score method and the ancestry informative markers in the Precision ID Ancestry Panel. Forensic Sci. Int. Genet. 2020, 44, 102154. [Google Scholar] [CrossRef] [PubMed]
  16. Al-Asfi, M.; McNevin, D.; Mehta, B.; Power, D.; Gahan, M.E.; Daniel, R. Assessment of the Precision ID Ancestry panel. Int. J. Legal Med. 2018, 132, 1581–1594. [Google Scholar] [CrossRef] [PubMed]
  17. Meyer, O.S.; Andersen, J.D.; Børsting, C. Presentation of the Human Pigmentation (HuPi) AmpliSeqTM custom panel. Forensic Sci. Int. Genet. Suppl. Ser. 2019, 7, 478–479. [Google Scholar] [CrossRef]
  18. Walsh, S.; Chaitanya, L.; Breslin, K.; Muralidharan, C.; Bronikowska, A.; Pospiech, E.; Koller, J.; Kovatsi, L.; Wollstein, A.; Branicki, W.; et al. Global skin colour prediction from DNA. Hum. Genet. 2017, 136, 847–863. [Google Scholar] [CrossRef] [Green Version]
  19. Chaitanya, L.; Breslin, K.; Zuñiga, S.; Wirken, L.; Pośpiech, E.; Kukla-Bartoszek, M.; Sijen, T.; de Knijff, P.; Liu, F.; Branicki, W.; et al. The HIrisPlex-S system for eye, hair and skin colour prediction from DNA: Introduction and forensic developmental validation. Forensic Sci. Int. Genet. 2018, 35, 123–135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kongshoj, B.; Thorleifsson, A.; Wulf, H.C. Pheomelanin and eumelanin in human skin determined by high-performance liquid chromatography and its relation to in vivo reflectance measurements. Photodermatol. Photoimmunol. Photomed. 2006, 22, 141–147. [Google Scholar] [CrossRef] [PubMed]
  21. Ravnbak, M.H. Objective determination of Fitzpatrick skin type. Dan. Med. Bull. 2010, 57, B4153. [Google Scholar]
  22. Farzad, M.S.; Pedersen, B.M.; Mogensen, H.S.; Børsting, C. Development of an automated AmpliSeqTM library building workflow for biological stain samples on the Biomek(®) 3000. Biotechniques 2020, 68, 342–344. [Google Scholar] [CrossRef] [Green Version]
  23. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [Green Version]
  24. Excoffier, L.; Lischer, H.E.L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef] [PubMed]
  25. Falush, D.; Stephens, M.; Pritchard, J.K. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 2003, 164, 1567–1587. [Google Scholar] [CrossRef]
  26. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of Population Structure Using Multilocus Genotype Data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  27. Jakobsson, M.; Rosenberg, N.A. CLUMPP: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 2007, 23, 1801–1806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Rosenberg, N.A. DISTRUCT: A program for the graphical display of population structure. Mol. Ecol. Notes 2004, 4, 137–138. [Google Scholar] [CrossRef]
  29. Simayijiang, H.; Børsting, C.; Tvedebrink, T.; Morling, N. Analysis of Uyghur and Kazakh populations using the Precision ID Ancestry Panel. Forensic Sci. Int. Genet. 2019, 43, 102144. [Google Scholar] [CrossRef] [PubMed]
  30. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S; Springer: New York, NY, USA, 2002. [Google Scholar]
  31. Shan, M.A.; Børsting, C.; Morling, N. Forensic application and genetic diversity of 21 autosomal STR loci in five major population groups of Pakistan. Int. J. Legal Med. 2021, 135, 775–777. [Google Scholar] [CrossRef]
  32. Korn, A.; Nourzaei, M. “Those were the hungry years”: A glimpse of Coastal Afro-Balochi. J. R. Asiat. Soc. 2018, 28, 661–695. [Google Scholar] [CrossRef] [Green Version]
  33. Siddiqi, M.H.; Akhtar, T.; Rakha, A.; Abbas, G.; Ali, A.; Haider, N.; Ali, A.; Hayat, S.; Masooma, S.; Ahmad, J.; et al. Genetic characterization of the Makrani people of Pakistan from mitochondrial DNA control-region data. Leg. Med. 2015, 17, 134–139. [Google Scholar] [CrossRef] [PubMed]
  34. Mallick, C.B.; Iliescu, F.M.; Möls, M.; Hill, S.; Tamang, R.; Chaubey, G.; Goto, R.; Ho, S.Y.W.; Gallego Romero, I.; Crivellaro, F.; et al. The light skin allele of SLC24A5 in South Asians and Europeans shares identity by descent. PLoS Genet. 2013, 9, e1003912. [Google Scholar] [CrossRef] [Green Version]
  35. Jonnalagadda, M.; Faizan, M.A.; Ozarkar, S.; Ashma, R.; Kulkarni, S.; Norton, H.L.; Parra, E. A Genome-Wide Association Study of Skin and Iris Pigmentation among Individuals of South Asian Ancestry. Genome Biol. Evol. 2019, 11, 1066–1076. [Google Scholar] [CrossRef] [Green Version]
  36. Lamason, R.L.; Mohideen, M.-A.P.K.; Mest, J.R.; Wong, A.C.; Norton, H.L.; Aros, M.C.; Jurynec, M.J.; Mao, X.; Humphreville, V.R.; Humbert, J.E.; et al. SLC24A5, a Putative Cation Exchanger, Affects Pigmentation in Zebrafish and Humans. Science 2005, 310, 1782–1786. [Google Scholar] [CrossRef] [Green Version]
  37. Nan, H.; Kraft, P.; Qureshi, A.A.; Guo, Q.; Chen, C.; Hankinson, S.E.; Hu, F.B.; Thomas, G.; Hoover, R.N.; Chanock, S.; et al. Genome-wide association study of tanning phenotype in a population of European ancestry. J. Investig. Dermatol. 2009, 129, 2250–2257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Cook, A.L.; Chen, W.; Thurber, A.E.; Smit, D.J.; Smith, A.G.; Bladen, T.G.; Brown, D.L.; Duffy, D.L.; Pastorino, L.; Bianchi-Scarra, G. Analysis of cultured human melanocytes based on polymorphisms within the SLC45A2/MATP, SLC24A5/NCKX5, and OCA2/P loci. J. Investig. Dermatol. 2009, 129, 392–405. [Google Scholar] [CrossRef] [Green Version]
  39. Chaki, M.; Sengupta, M.; Mondal, M.; Bhattacharya, A.; Mallick, S. Molecular and functional studies of tyrosinase variants among Indian oculocutaneous albinism type 1 patients. J. Investig. Dermatol. 2011, 131, 260–262. [Google Scholar] [CrossRef]
  40. Jagirdar, K.; Smit, D.J.; Ainger, S.A.; Lee, K.J.; Brown, D.L.; Chapman, B.; Zhen Zhao, Z.; Montgomery, G.W.; Martin, N.G.; Stow, J.L. Molecular analysis of common polymorphisms within the human Tyrosinase locus and genetic association with pigmentation traits. Pigment Cell Melanoma Res. 2014, 27, 552–564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Eiberg, H.; Troelsen, J.; Nielsen, M.; Mikkelsen, A.; Mengel-From, J.; Kjaer, K.W.; Hansen, L. Blue eye color in humans may be caused by a perfectly associated founder mutation in a regulatory element located within the HERC2 gene inhibiting OCA2 expression. Hum. Genet. 2008, 123, 177–187. [Google Scholar] [CrossRef] [PubMed]
  42. Sturm, R.A.; Duffy, D.L.; Zhao, Z.Z.; Leite, F.P.N.; Stark, M.S.; Hayward, N.K.; Martin, N.G.; Montgomery, G.W. A single SNP in an evolutionary conserved region within intron 86 of the HERC2 gene determines human blue-brown eye color. Am. J. Hum. Genet. 2008, 82, 424–431. [Google Scholar] [CrossRef] [Green Version]
  43. Liu, F.; van Duijn, K.; Vingerling, J.R.; Hofman, A.; Uitterlinden, A.G.; Janssens, A.C.J.W.; Kayser, M. Eye color and the prediction of complex phenotypes from genotypes. Curr. Biol. 2009, 19, R192–R193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Visser, M.; Kayser, M.; Palstra, R.-J. HERC2 rs12913832 modulates human pigmentation by attenuating chromatin-loop formation between a long-range enhancer and the OCA2 promoter. Genome Res. 2012, 22, 446–455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Pavan, W.J.; Sturm, R.A. The genetics of human skin and hair pigmentation. Annu. Rev. Genomics Hum. Genet. 2019, 20, 41–72. [Google Scholar] [CrossRef]
  46. Kosoy, R.; Nassir, R.; Tian, C.; White, P.A.; Butler, L.M.; Silva, G.; Kittles, R.; Alarcon-Riquelme, M.E.; Gregersen, P.K.; Belmont, J.W.; et al. Ancestry informative marker sets for determining continental origin and admixture proportions in common populations in America. Hum. Mutat. 2009, 30, 69–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Kidd, K.K.; Speed, W.C.; Pakstis, A.J.; Furtado, M.R.; Fang, R.; Madbouly, A.; Maiers, M.; Middha, M.; Friedlaender, F.R.; Kidd, J.R. Progress toward an efficient panel of SNPs for ancestry inference. Forensic Sci. Int. Genet. 2014, 10, 23–32. [Google Scholar] [CrossRef] [Green Version]
  48. Phillips, C.; Parson, W.; Lundsberg, B.; Santos, C.; Freire-Aradas, A.; Torres, M.; Eduardoff, M.; Børsting, C.; Johansen, P.; Fondevila, M.; et al. Building a forensic ancestry panel from the ground up: The EUROFORGEN Global AIM-SNP set. Forensic Sci. Int. Genet. 2014, 11, 13–25. [Google Scholar] [CrossRef] [Green Version]
  49. Phillips, C.; McNevin, D.; Kidd, K.K.; Lagacé, R.; Wootton, S.; de la Puente, M.; Freire-Aradas, A.; Mosquera-Miguel, A.; Eduardoff, M.; Gross, T.; et al. MAPlex—A massively parallel sequencing ancestry analysis multiplex for Asia-Pacific populations. Forensic Sci. Int. Genet. 2019, 42, 213–226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 12 00733 g001 550
Figure 1. Boxplots of skin pigmentation measurements (pigment protection factor) of the upper arm, lower arm, and forehead in three Pakistani subpopulations. The pigment protection factor was measured in triplicate at each measured skin area.
Figure 1. Boxplots of skin pigmentation measurements (pigment protection factor) of the upper arm, lower arm, and forehead in three Pakistani subpopulations. The pigment protection factor was measured in triplicate at each measured skin area.
Genes 12 00733 g001
Genes 12 00733 g002 550
Figure 2. STRUCTURE analysis with K = 4 using 160 ancestry informative markers (AIMs) and six meta-populations. Each cluster (K) is represented by a color. Population abbreviations used: SS Africa: Sub-Saharan Africa; N Africa: North Africa; M East: Middle East; SC Asia: South-Central Asia; E Asia: East Asia.
Figure 2. STRUCTURE analysis with K = 4 using 160 ancestry informative markers (AIMs) and six meta-populations. Each cluster (K) is represented by a color. Population abbreviations used: SS Africa: Sub-Saharan Africa; N Africa: North Africa; M East: Middle East; SC Asia: South-Central Asia; E Asia: East Asia.
Genes 12 00733 g002
Genes 12 00733 g003 550
Figure 3. Cumulative correlations between (1) four SNPs and ancestry and (2) upper arm, lower arm, and forehead pigmentation. The SNPs were ranked based on their correlations with skin color (adjusted R2). Ancestry is the proportion of Sub-Saharan African and European ancestry.
Figure 3. Cumulative correlations between (1) four SNPs and ancestry and (2) upper arm, lower arm, and forehead pigmentation. The SNPs were ranked based on their correlations with skin color (adjusted R2). Ancestry is the proportion of Sub-Saharan African and European ancestry.
Genes 12 00733 g003
Genes 12 00733 g004 550
Figure 4. Boxplots showing the skin color category predicted with the HIrisPlex-S prediction model and the pigmentation protection factor on the upper arm, lower arm, and forehead.
Figure 4. Boxplots showing the skin color category predicted with the HIrisPlex-S prediction model and the pigmentation protection factor on the upper arm, lower arm, and forehead.
Genes 12 00733 g004
Table
Table 1. Five SNPs with statistically significant association with skin pigmentation on the upper arm, lower arm, or forehead in 299 Pakistani individuals from three subpopulations (Baloch, Pashtun, and Punjabi).
Table 1. Five SNPs with statistically significant association with skin pigmentation on the upper arm, lower arm, or forehead in 299 Pakistani individuals from three subpopulations (Baloch, Pashtun, and Punjabi).
Upper ArmLower ArmForehead
GeneSNPAdjusted R2Sign.1Adjusted R2Sign. 1Adjusted R2Sign. 1
TYRrs10426020.13***0.11**0.071**
GRM5rs108314960.029–0.047*0.014–
HERC2rs129138320.025–0.016–0.064*
SLC24A5rs14266540.16***0.10**0.082**
SLC45A2rs168919820.075**0.028–0.058**
1 Statistical significance (Sign.) after Bonferroni correction: * p < 0.05, ** p < 10−3, *** p < 10−6.
Table
Table 2. Multiple loci correlations with skin pigmentation on the upper arm, lower arm, or forehead in 299 Pakistani individuals from three subpopulations (Baloch, Pashtun, and Punjabi).
Table 2. Multiple loci correlations with skin pigmentation on the upper arm, lower arm, or forehead in 299 Pakistani individuals from three subpopulations (Baloch, Pashtun, and Punjabi).
Upper ArmLower ArmForehead
Adjusted R2Sign. 1Adjusted R2Sign. 1Adjusted R2Sign. 1
Five SNPs 20.33***0.24***0.24***
Best model 30.33***0.23***0.24***
Best model + ancestry 40.37***0.28***0.30***
1 Statistical significance (Sign.): *** p < 10−6. 2 rs1042602, rs10831496, rs12913832, rs1426654, and rs16891982. 3 rs1042602, rs12913832, rs1426654, and rs16891982. 4 Ancestry is the estimated proportion of European and Sub-Saharan African ancestry.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shan, M.A.; Meyer, O.S.; Refn, M.; Morling, N.; Andersen, J.D.; Børsting, C. Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch. Genes 2021, 12, 733. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12050733

AMA Style

Shan MA, Meyer OS, Refn M, Morling N, Andersen JD, Børsting C. Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch. Genes. 2021; 12(5):733. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12050733

Chicago/Turabian Style

Shan, Muhammad Adnan, Olivia Strunge Meyer, Mie Refn, Niels Morling, Jeppe Dyrberg Andersen, and Claus Børsting. 2021. "Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch" Genes 12, no. 5: 733. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12050733

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop