Skin pigmentation variation among populations of West Maharashtra, India

A Cross-Sectional Assessment of Quantitative Epidermal Melanin and Erythema Indices among North Indians

Background

Remarkable diversity of skin tones among Indians ranging from pale pinkish to dark brown appears to be an ideal choice for the assessment of skin pigment variation.

Aim

The present study was designed to assess the variation observed in melanin and erythema indices among north Indians.

Material and Methods

Skin reflectance data (n = 574) was collected from six diverse populations of north India using DermaSpectrometer (DSM II ColorMeter) followed by statistical analysis to investigate the impact of geographical location and gender on constitutive skin pigmentation.

Results

The melanin index (MI) varied between 17.93 and 56.92 (Mean (M) = 35.80 ± 6.26) whereas the erythema index (EI) varied between 4.92 and 18.82 (M = 10.48 ± 2.68). MI and EI of females were found to be significantly lower than males (P < 0.001). Geographical location exhibited a significant association with MI and EI (P < 0.001). Furthermore, we have noted a positive correlation between MI and EI (P < 0.001).

Conclusion

The study has refined our understanding of skin pigmentation variation among north Indians in terms of significant association with geographical location {MI: F _(5,568) = 31.07, P < 0.001; EI: F _(5,568) = 73.37, P < 0.001} and gender {MI: t ₍₃₈₆₎ = -4.06, P < 0.001; EI: t ₍₃₈₆₎ = -11.96, P < 0.001} and rendered opportunities for further studies.

Identifying signatures of natural selection in Indian populations

In this study, we present the results of a genome-wide scan for signatures of positive selection using data from four tribal groups (Kokana, Warli, Bhil, and Pawara) and two caste groups (Deshastha Brahmin and Kunbi Maratha) from West of the Maharashtra State In India, as well as two samples of South Asian ancestry from the 1KG project (Gujarati Indian from Houston, Texas and Indian Telugu from UK). We used an outlier approach based on different statistics, including PBS, xpEHH, iHS, CLR, Tajima’s D, as well as two recently developed methods: Graph-aware Retrieval of Selective Sweeps (GRoSS) and Ascertained Sequentially Markovian Coalescent (ASMC). In order to minimize the risk of false positives, we selected regions that are outliers in all the samples included in the study using more than one method. We identified putative selection signals in 107 regions encompassing 434 genes. Many of the regions overlap with only one gene. The signals observed using microarray-based data are very consistent with our analyses using high-coverage sequencing data, as well as those identified with a novel coalescence-based method (ASMC). Importantly, at least 24 of these genomic regions have been identified in previous selection scans in South Asian populations or in other population groups. Our study highlights genomic regions that may have played a role in the adaptation of anatomically modern humans to novel environmental conditions after the out of Africa migration.

MC1R diversity and its role in skin pigmentation variation in West Maharashtra, India

OBJECTIVES

MC1R polymorphisms have been reported to be under a selective constraint in populations inhabiting high UVR regions such as Africans; however, these patterns are not consistent. Here we analyze the MC1R gene in West Maharashtra, India to see if sequence diversity corresponds to their diverse pigmentary profiles and if MC1R is constrained in dark skinned tribal as compared to lighter skinned caste populations.

METHODS

A 2648 bp region of this gene was sequenced in 102 individuals and the data was compared for π, ϴ diversity indices. Tajima's D was assessed for signatures of purifying selection and MC1R variants were associated with MI measures using the additive, dominant, and recessive models. Pairwise F_ST was tested among study populations and between study populations and 1000 Genomes regional samples.

RESULTS

MC1R diversity was not uniquely patterned among castes and tribes. Non-synonymous variants rs2228479A, rs1805007_T, and rs885479_A showed low variability in these populations. Selection tests did not indicate any constraint on MC1R and pairwise F_ST were also low among the study populations (-0.0163 to 0.06112). The SNP rs3212359 was significantly associated with MI measures when tested using different association models.

CONCLUSIONS

We do not find evidence of a selective constraint on MC1R. The presence of a large number of unique haplotypes and low F_ST values at this locus suggests that MC1R polymorphisms may not be influencing pigmentation variation among castes and tribes in this region. Observed associations between rs3212359 and MI measures need to be validated through studies on larger samples and in-vitro functional studies.

A study of genomic diversity in populations of Maharashtra, India, inferred from 20 autosomal STR markers

OBJECTIVE

This study was planned to evaluate the genetic diversity in the admixed and Teli (a Hindu caste) populations of Maharashtra, India using 20 autosomal Short Tandem Repeat (STR) genetic markers. We further investigated the genetic relatedness of the studied populations with other Indian populations.

RESULTS

The studied populations showed a wide range of observed heterozygosity viz. 0.690 to 0.918 for the admixed population and 0.696 to 0.942 for the Teli population. This might be due to the multi-directional gene flow. The admixed and Teli populations also showed a high degree polymorphism which ranged from 0.652 to 0.903 and 0.644 to 0.902, respectively. Their combined value of matching probability for all the studied loci was 4.29 × 10-25 and 5.01 × 10-24, respectively. The results of Neighbor-Joining tree and Principal Component Analysis showed that the studied populations clustered with the general populations of Jharkhand, UttarPradesh, Rajasthan and Central Indian States, as well as with the specific populations of Maharashtra (Konkanastha Brahmins) and Tamil Nadu (Kurmans). Overall, the obtained data showed a high degree of forensic efficacy and would be useful for forensic applications as well as genealogical studies.

Novel insights on demographic history of tribal and caste groups from West Maharashtra (India) using genome-wide data

The South Asian subcontinent is characterized by a complex history of human migrations and population interactions. In this study, we used genome-wide data to provide novel insights on the demographic history and population relationships of six Indo-European populations from the Indian State of West Maharashtra. The samples correspond to two castes (Deshastha Brahmins and Kunbi Marathas) and four tribal groups (Kokana, Warli, Bhil and Pawara). We show that tribal groups have had much smaller effective population sizes than castes, and that genetic drift has had a higher impact in tribal populations. We also show clear affinities between the Bhil and Pawara tribes, and to a lesser extent, between the Warli and Kokana tribes. Our comparisons with available modern and ancient DNA datasets from South Asia indicate that the Brahmin caste has higher Ancient Iranian and Steppe pastoralist contributions than the Kunbi Marathas caste. Additionally, in contrast to the two castes, tribal groups have very high Ancient Ancestral South Indian (AASI) contributions. Indo-European tribal groups tend to have higher Steppe contributions than Dravidian tribal groups, providing further support for the hypothesis that Steppe pastoralists were the source of Indo-European languages in South Asia, as well as Europe.

A Genome-Wide Association Study of Skin and Iris Pigmentation among Individuals of South Asian Ancestry

South Asia has a complex history of migrations and is characterized by substantial pigmentary and genetic diversity. For this reason, it is an ideal region to study the genetic architecture of normal pigmentation variation. Here, we present a meta-analysis of two genome-wide association studies (GWASs) of skin pigmentation using skin reflectance (M-index) as a quantitative phenotype. The meta-analysis includes a sample of individuals of South Asian descent living in Canada (N = 348), and a sample of individuals from two caste and four tribal groups from West Maharashtra, India (N = 480). We also present the first GWAS of iris color in South Asian populations. This GWAS was based on quantitative measures of iris color obtained from high-resolution iris pictures. We identified genome-wide significant associations of variants within the well-known gene SLC24A5, including the nonsynonymous rs1426654 polymorphism, with both skin pigmentation and iris color, highlighting the pleiotropic effects of this gene on pigmentation. Variants in the HERC2 gene (e.g., rs12913832) were also associated with iris color and iris heterochromia. Our study emphasizes the usefulness of quantitative methods to study iris color variation. We also identified novel genome-wide significant associations with skin pigmentation and iris color, but we could not replicate these associations due to the lack of independent samples. It will be critical to expand the number of studies in South Asian populations in order to better understand the genetic variation driving the diversity of skin pigmentation and iris color observed in this region.

Shades of complexity: New perspectives on the evolution and genetic architecture of human skin

Like many highly variable human traits, more than a dozen genes are known to contribute to the full range of skin color. However, the historical bias in favor of genetic studies in European and European-derived populations has blinded us to the magnitude of pigmentation's complexity. As deliberate efforts are being made to better characterize diverse global populations and new sequencing technologies, better measurement tools, functional assessments, predictive modeling, and ancient DNA analyses become more widely accessible, we are beginning to appreciate how limited our understanding of the genetic bases of human skin color have been. Novel variants in genes not previously linked to pigmentation have been identified and evidence is mounting that there are hundreds more variants yet to be found. Even for genes that have been exhaustively characterized in European populations like MC1R, OCA2, and SLC24A5, research in previously understudied groups is leading to a new appreciation of the degree to which genetic diversity, epistatic interactions, pleiotropy, admixture, global and local adaptation, and cultural practices operate in population-specific ways to shape the genetic architecture of skin color. Furthermore, we are coming to terms with how factors like tanning response and barrier function may also have influenced selection on skin throughout human history. By examining how our knowledge of pigmentation genetics has shifted in the last decade, we can better appreciate how far we have come in understanding human diversity and the still long road ahead for understanding many complex human traits.

The influences of genes, the environment, and social factors on the evolution of skin color diversity in India

OBJECTIVES

Skin color is a highly visible and variable trait across human populations. It is not yet clear how evolutionary forces interact to generate phenotypic diversity. Here, we sought to unravel through an integrative framework the role played by three factors-demography and migration, sexual selection, and natural selection-in driving skin color diversity in India.

METHODS

Skin reflectance data were collected from 10 diverse socio-cultural populations along the latitudinal expanse of India, including both sexes. We first looked at how skin color varies within and between these populations. Second, we compared patterns of sexual dimorphism in skin color. Third, we studied the influence of ultraviolet radiation on skin color throughout India. Finally, we attempted to disentangle the interactions between these factors in the context of available genetic data.

RESULTS

We found that the relative importance of these forces varied between populations. Social factors and population structure have played a stronger role than natural selection in shaping skin color diversity across India. Phenotypic overprinting resulted from additional genetic mutations overriding the skin lightening effect of variants such as the SLC24A5 rs1426654-A allele in some populations, in the context of the variable influence of sexual selection. Furthermore, specific genotypes are not associated reliably with specific skin color phenotypes. This result has relevance for DNA forensics and ancient DNA research.

CONCLUSIONS

India is a crucible of macro- and micro-evolutionary forces, and the complex interactions of physical and social forces are visible in the patterns of skin color seen today in the country.

Loci associated with skin pigmentation identified in African populations

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

Identifying signatures of positive selection in pigmentation genes in two South Asian populations

OBJECTIVES

Skin pigmentation is a polygenic trait showing wide phenotypic variations among global populations. While numerous pigmentation genes have been identified to be under positive selection among European and East populations, genes contributing to phenotypic variation in skin pigmentation within and among South Asian populations are still poorly understood. The present study uses data from the Phase 3 of the 1000 genomes project focusing on two South Asian populations-GIH (Gujarati Indian from Houston, Texas) and ITU (Indian Telugu from UK), so as to decode the genetic architecture involved in adaptation to ultraviolet radiation in South Asian populations.

METHODS

Statistical tests included were (1) tests to identify deviations of the Site Frequency Spectrum (SFS) from neutral expectations (Tajima's D, Fay and Wu's H and Fu and Li's D* and F*), (2) tests focused on the identification of high-frequency haplotypes with extended linkage disequilibrium (iHS and Rsb), and (3) tests based on genetic differentiation between populations (LSBL).

RESULTS

Twenty-two pigmentation genes fall in the top 1% for at least one statistic in the GIH population, 5 of which (LYST, OCA2, SLC24A5, SLC45A2, and TYR) have been previously associated with normal variation in skin, hair, or eye color. In comparison, 17 genes fall in the top 1% for at least one statistic in the ITU population. Twelve loci which are identified as outliers in the ITU scan were also identified in the GIH population.

CONCLUSIONS

These results suggest that selection may have affected these loci broadly across the region.

Genotype-Phenotype Study of the Middle Gangetic Plain in India Shows Association of rs2470102 with Skin Pigmentation

Our understanding of the genetics of skin pigmentation has been largely skewed towards populations of European ancestry, imparting less attention to South Asian populations, who behold huge pigmentation diversity. Here, we investigate skin pigmentation variation in a cohort of 1,167 individuals in the Middle Gangetic Plain of the Indian subcontinent. Our data confirm the association of rs1426654 with skin pigmentation among South Asians, consistent with previous studies, and also show association for rs2470102 single nucleotide polymorphism. Our haplotype analyses further help us delineate the haplotype distribution across social categories and skin color. Taken together, our findings suggest that the social structure defined by the caste system in India has a profound influence on the skin pigmentation patterns of the subcontinent. In particular, social category and associated single nucleotide polymorphisms explain about 32% and 6.4%, respectively, of the total phenotypic variance. Phylogeography of the associated single nucleotide polymorphisms studied across 52 diverse populations of the Indian subcontinent shows wide presence of the derived alleles, although their frequencies vary across populations. Our results show that both polymorphisms (rs1426654 and rs2470102) play an important role in the skin pigmentation diversity of South Asians.

Association of genetic variants with skin pigmentation phenotype among populations of west Maharashtra, India

OBJECTIVES

South Asians exhibit extensive variation in skin melanin index (MI) which is observed across the broader region of South Asia as well as within restricted geographic regions. However, the genetic variants associated with variation in the skin pigmentation phenotype are poorly understood in these populations. The present study examines the association between MI measures and genetic variants from 5 candidate pigmentation genes among 533 individuals representing 6 populations of West Maharashtra.

METHODS

Associations between five single nucleotide polymorphisms (SNPs) known to play a role in pigmentation (rs1426654-SLC24A5, rs1042602-TYR, rs16891982-SLC45A2, rs6058017-ASIP, and rs642742-KITLG) and MI measures were tested using standard one-way analysis of variance (ANOVA) within each population. Multiple linear regression was used to test the effects of these SNPs in the full West Maharashtra sample using sex, age, and population or social group as covariates.

RESULTS

rs1426654 showed significant association with MI in all six study populations (P < 0.01). Association tests using sex, age, and population as covariates showed rs1426654 and rs1042602 to be significantly (P < 0.01) associated with lighter skin pigmentation in West Maharashtra as a whole. By contrast, when social group was added as a covariate instead of population, rs1426654, rs1042602, and rs16891982 were significantly (P < 0.01) associated with lighter skin pigmentation.

CONCLUSIONS

Only rs1426654 is significantly associated with MI in each individual population; however, rs1426654, rs1042602, and rs16891982 are significantly associated with pigmentation in the broader West Maharashtra region after controlling for population and social group, with rs1426654 (SLC24A5) explaining the majority of the observed variation. Am. J. Hum. Biol. 28:610-618, 2016. © 2016 Wiley Periodicals, Inc.

Quantitative assessment of skin, hair, and iris variation in a diverse sample of individuals and associated genetic variation

OBJECTIVES

The main goals of this study are to 1) quantitatively measure skin, hair, and iris pigmentation in a diverse sample of individuals, 2) describe variation within and between these samples, and 3) demonstrate how quantitative measures can facilitate genotype-phenotype association tests.

MATERIALS AND METHODS

We quantitatively characterize skin, hair, and iris pigmentation using the Melanin (M) Index (skin) and CIELab values (hair) in 1,450 individuals who self-identify as African American, East Asian, European, Hispanic, or South Asian. We also quantify iris pigmentation in a subset of these individuals using CIELab values from high-resolution iris photographs. We compare mean skin M index and hair and iris CIELab values among populations using ANOVA and MANOVA respectively and test for genotype-phenotype associations in the European sample.

RESULTS

All five populations are significantly different for skin (P <2 × 10(-16) ) and hair color (P <2 × 10(-16) ). Our quantitative analysis of iris and hair pigmentation reinforces the continuous, rather than discrete, nature of these traits. We confirm the association of three loci (rs16891982, rs12203592, and rs12913832) with skin pigmentation and four loci (rs12913832, rs12203592, rs12896399, and rs16891982) with hair pigmentation. Interestingly, the derived rs12203592 T allele located within the IRF4 gene is associated with lighter skin but darker hair color.

DISCUSSION

The quantitative methods used here provide a fine-scale assessment of pigmentation phenotype and facilitate genotype-phenotype associations, even with relatively small sample sizes. This represents an important expansion of current investigations into pigmentation phenotype and associated genetic variation by including non-European and admixed populations. Am J Phys Anthropol 160:570-581, 2016. © 2015 Wiley Periodicals, Inc.