The color of normal: How a Eurocentric focus erases pigmentation complexity.

OBJECTIVES: Skin pigmentation is both a highly variable and highly visible human phenotypic trait. Investigations into the biology and origins of this variation have been the focus of research in the fields of dermatology, anthropology, and forensic science, among others. This manuscript explores how much of what we know about the biology, genetics, and evolutionary origins of pigmentation has been strongly influenced by investigations and applications that focus on lighter skin. METHODS: I reviewed literature from the fields of dermatology, anthropology and evolutionary genetics, and forensic science to assess how perceptions of lighter skin as the "normal" state in humans can shape the ways that knowledge is gathered and applied in these fields. RESULTS: This normalization of lighter skin has impacted common tools used in dermatology and shaped the framework of dermatological education. A strong Eurocentric bias has shaped our understanding of the genetic architecture of pigmentary traits, which influences the ways in we understand the evolutionary processes leading to modern pigmentation diversity. Finally, I discuss how these biases in pigmentation genetics work in combination with phenotypic systems that privilege predicting lighter pigmentation variation to impede accurate prediction of intermediate phenotypes, particularly in individuals with ancestry from multiple populations. This can lead to a disproportionate targeting of already over-policed populations with darker skin. CONCLUSIONS: Potential changes to how we conceptualize clinical and basic pigmentation research may help to reduce existing health disparities and improve understanding of pigmentation genetic architecture and how this knowledge is applied in forensic contexts.

Meta-analysis of GWA studies provides new insights on the genetic architecture of skin pigmentation in recently admixed populations.

BACKGROUND: Association studies in recently admixed populations are extremely useful to identify the genetic architecture of pigmentation, due to their high genotypic and phenotypic variation. However, to date only four Genome-Wide Association Studies (GWAS) have been carried out in these populations. RESULTS: We present a GWAS of skin pigmentation in an admixed sample from Cuba (N = 762). Additionally, we conducted a meta-analysis including the Cuban sample, and admixed samples from Cape Verde, Puerto Rico and African-Americans from San Francisco. This meta-analysis is one of the largest efforts so far to characterize the genetic basis of skin pigmentation in admixed populations (N = 2,104). We identified five genome-wide significant regions in the meta-analysis, and explored if the markers observed in these regions are associated with the expression of relevant pigmentary genes in human melanocyte cultures. In three of the regions identified in the meta-analysis (SLC24A5, SLC45A2, and GRM5/TYR), the association seems to be driven by non-synonymous variants (rs1426654, rs16891982, and rs1042602, respectively). The rs16891982 polymorphism is strongly associated with the expression of the SLC45A2 gene. In the GRM5/TYR region, in addition to the rs1042602 non-synonymous SNP located on the TYR gene, variants located in the nearby GRM5 gene have an independent effect on pigmentation, possibly through regulation of gene expression of the TYR gene. We also replicated an association recently described near the MFSD12 gene on chromosome 19 (lead variant rs112332856). Additionally, our analyses support the presence of multiple signals in the OCA2/HERC2/APBA2 region on chromosome 15. A clear causal candidate is the HERC2 intronic variant rs12913832, which has a profound influence on OCA2 expression. This variant has pleiotropic effects on eye, hair, and skin pigmentation. However, conditional and haplotype-based analyses indicate the presence of other variants with independent effects on melanin levels in OCA2 and APBA2. Finally, a follow-up of genome-wide signals identified in a recent GWAS for tanning response indicates that there is a substantial overlap in the genetic factors influencing skin pigmentation and tanning response. CONCLUSIONS: Our meta-analysis of skin pigmentation GWAS in recently admixed populations provides new insights about the genetic architecture of this complex trait.

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.

Understanding influences of culture and history on mtDNA variation and population structure in three populations from Assam, Northeast India.

OBJECTIVES: Positioned at the nexus of India, China, and Southeast Asia, Northeast India is presumed to have served as a channel for land-based human migration since the Upper Pleistocene. Assam is the largest state in the Northeast. We characterized the genetic background of three populations and examined the ways in which their population histories and cultural practices have influenced levels of intrasample and intersample variation. METHODS: We examined sequence data from the mtDNA hypervariable control region and selected diagnostic mutations from the coding region in 128 individuals from three ethnic groups currently living in Assam: two Scheduled tribes (Sonowal Kachari and Rabha), and the non-Scheduled Tai Ahom. RESULTS: The populations of Assam sampled here express mtDNA lineages indicative of South Asian, Southeast Asian, and East Asian ancestry. We discovered two completely novel haplogroups in Assam that accounted for 6.2% of the lineages in our sample. We also identified a new subhaplogroup of M9a that is prevalent in the Sonowal Kachari of Assam (19.1%), but not present in neighboring Arunachal Pradesh, indicating substantial regional population structuring. Employing a large comparative dataset into a series of multidimensional scaling (MDS) analyses, we saw the Rabha cluster with populations sampled from Yunnan Province, indicating that the historical matrilineality of the Rabha has maintained lineages from Southern China. CONCLUSION: Assam has undergone multiple colonization events in the time since the initial peopling event, with populations from Southern China and Southeast Asia having the greatest influence on maternal lineages in the region.

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.

The rs387907171 SNP in TYRP1 is not associated with blond hair color on the Island of Bougainville.

OBJECTIVES: The aims of this study are to characterize the frequency of the derived allele at rs387907171 in populations from the islands of New Britain and Bougainville in Northern Island Melanesia, to confirm its association with lighter hair color, and to refine hypotheses regarding its evolutionary history. METHODS: rs387907171 was genotyped in 93 individuals from New Britain and 101 from Bougainville for whom quantitative assessments of skin and hair pigmentation were available. Combining these with existing data from other Melanesian islands we tested for differences in allele frequencies between islands and for associations with skin and hair pigmentation using ANOVA, including sex, age, and island affiliations as covariates. RESULTS: The derived allele at rs387907171 was observed in a single copy in the New Britain and Bougainville populations genotyped here. Its frequency differs significantly among islands in the region (χ(2) = 206.5, df = 3, P < 0.001). rs387907171 remains significantly, although weakly, associated with lighter hair pigmentation (F = 10.28, R(2) = 0.0125, P = 0.0014). This association increases when sex and age (F = 20.68, R(2) = 0.074, P < 7.92 × 10(-13) ) are included as covariates. CONCLUSIONS: The rs387907171 SNP exhibits strong allele frequency differences among islands in Northern Island Melanesia. Its absence from Bougainville, as well as the weak association with decreased hair color, indicates that additional alleles contribute to the blondism phenotype. Its geographic distribution suggests that a Lapita-mediated model for the dispersal of the derived allele at rs387907171 remains a viable evolutionary scenario. Am. J. Hum. Biol. 28:431-435, 2016. © 2015 Wiley Periodicals, Inc.

MC1R diversity in Northern Island Melanesia has not been constrained by strong purifying selection and cannot explain pigmentation phenotype variation in the region.

BACKGROUND: Variation in human skin pigmentation evolved in response to the selective pressure of ultra-violet radiation (UVR). Selection to maintain darker skin in high UVR environments is expected to constrain pigmentation phenotype and variation in pigmentation loci. Consistent with this hypothesis, the gene MC1R exhibits reduced diversity in African populations from high UVR regions compared to low-UVR non-African populations. However, MC1R diversity in non-African populations that have evolved under high-UVR conditions is not well characterized. METHODS: In order to test the hypothesis that MC1R variation has been constrained in Melanesians the coding region of the MC1R gene was sequenced in 188 individuals from Northern Island Melanesia. The role of purifying selection was assessed using a modified McDonald Kreitman's test. Pairwise FST was calculated between Melanesian populations and populations from the 1000 Genomes Project. The SNP rs2228479 was genotyped in a larger sample (n = 635) of Melanesians and tested for associations with skin and hair pigmentation. RESULTS: We observe three nonsynonymous and two synonymous mutations. A modified McDonald Kreitman's test failed to detect a significant signal of purifying selection. Pairwise FST values calculated between the four islands sampled here indicate little regional substructure in MC1R. When compared to African, European, East and South Asian populations, Melanesians do not exhibit reduced population divergence (measured as FST) or a high proportion of haplotype sharing with Africans, as one might expect if ancestral haplotypes were conserved across high UVR populations in and out of Africa. The only common nonsynonymous polymorphism observed, rs2228479, is not significantly associated with skin or hair pigmentation in a larger sample of Melanesians. CONCLUSIONS: The pattern of sequence diversity here does not support a model of strong selective constraint on MC1R in Northern Island Melanesia This absence of strong constraint, as well as the recent population history of the region, may explain the observed frequencies of the derived rs2228479 allele. These results emphasize the complex genetic architecture of pigmentation phenotypes, which are controlled by multiple, possibly interacting loci. They also highlight the role that population history can play in influencing phenotypic diversity in the absence of strong natural selection.

Distribution of an allele associated with blond hair color across Northern Island Melanesia.

Pigmentation of the skin, hair, and eyes is a complex trait controlled by multiple genetic loci. Recently a non-synonymous mutation in the pigmentation candidate gene TYRP1 was shown to be significantly associated with a blond-hair phenotype in populations from the Solomon Islands. The distribution of this mutation in the islands of Northern Island Melanesia, where the blondism phenotype is also prevalent, was unknown. Here, we present data describing the distribution of this allele in 550 individuals sampled from across this region, and test for associations between genotype at this locus and quantitatively measured skin and hair pigmentation phenotype. We report that the frequency of the 93C allele is notably lower than observed in the Solomons (0.12 vs. 0.26). The allele exhibits significant geographic heterogeneity across the islands sampled (χ(2) = 108.4, P < 0.0001). It is observed at its highest frequencies on the islands of New Ireland and New Hanover, while being almost completely absent from the large island of New Britain. Using linear regression with age, sex, and island as covariates we report that, as in the Solomons, the 93C allele is significantly associated with a decrease in hair pigmentation but not skin pigmentation. We discuss the distribution of the 93C allele across the Southwest Pacific in light of its possible place of origin and dispersal.

Insights into the genetic architecture of the human face.

The human face is complex and multipartite, and characterization of its genetic architecture remains challenging. Using a multivariate genome-wide association study meta-analysis of 8,246 European individuals, we identified 203 genome-wide-significant signals (120 also study-wide significant) associated with normal-range facial variation. Follow-up analyses indicate that the regions surrounding these signals are enriched for enhancer activity in cranial neural crest cells and craniofacial tissues, several regions harbor multiple signals with associations to different facial phenotypes, and there is evidence for potential coordinated actions of variants. In summary, our analyses provide insights into the understanding of how complex morphological traits are shaped by both individual and coordinated genetic actions.

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.