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.

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.

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.

Genome-wide association study of pigmentary traits (skin and iris color) in individuals of East Asian ancestry

Background

Currently, there is limited knowledge about the genetics underlying pigmentary traits in East Asian populations. Here, we report the results of the first genome-wide association study of pigmentary traits (skin and iris color) in individuals of East Asian ancestry.

Methods

We obtained quantitative skin pigmentation measures (M-index) in the inner upper arm of the participants using a portable reflectometer (N = 305). Quantitative measures of iris color (expressed as L*, a* and b* CIELab coordinates) were extracted from high-resolution iris pictures (N = 342). We also measured the color differences between the pupillary and ciliary regions of the iris (e.g., iris heterochromia). DNA samples were genotyped with Illumina’s Infinium Multi-Ethnic Global Array (MEGA) and imputed using the 1000 Genomes Phase 3 samples as reference haplotypes.

Results

For skin pigmentation, we did not observe any genome-wide significant signal. We followed-up in three independent Chinese samples the lead SNPs of five regions showing multiple common markers (minor allele frequency ≥ 5%) with good imputation scores and suggestive evidence of association (p-values < 10⁻⁵). One of these markers, rs2373391, which is located in an intron of the ZNF804B gene on chromosome 7, was replicated in one of the Chinese samples (p = 0.003). For iris color, we observed genome-wide signals in the OCA2 region on chromosome 15. This signal is driven by the non-synonymous rs1800414 variant, which explains 11.9%, 10.4% and 6% of the variation observed in the b*, a* and L* coordinates in our sample, respectively. However, the OCA2 region was not associated with iris heterochromia.

Discussion

Additional genome-wide association studies in East Asian samples will be necessary to further disentangle the genetic architecture of pigmentary traits in East Asian populations.

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.

Distribution of two OCA2 polymorphisms associated with pigmentation in East-Asian populations

Two OCA2 polymorphisms (rs1800414 and rs74653330) have been associated with pigmentation in East Asians. We explored the distribution of these markers in a panel of samples from populations around the world. The derived allele of rs1800414 has high frequencies in a broad East-Asian region, whereas the derived allele of rs74653330 is primarily restricted to northern East Asia. Our data suggest that these polymorphisms may have been selected independently in different regions of East Asia.

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.

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.

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.

Genetic evidence for the convergent evolution of light skin in Europeans and East Asians

Human skin pigmentation shows a strong positive correlation with ultraviolet radiation intensity, suggesting that variation in skin color is, at least partially, due to adaptation via natural selection. We investigated the evolution of pigmentation variation by testing for the presence of positive directional selection in 6 pigmentation genes using an empirical F(ST) approach, through an examination of global diversity patterns of these genes in the Centre d'Etude du Polymorphisme Humain (CEPH)-Diversity Panel, and by exploring signatures of selection in data from the International HapMap project. Additionally, we demonstrated a role for MATP in determining normal skin pigmentation variation using admixture mapping methods. Taken together (with the results of previous admixture mapping studies), these results point to the importance of several genes in shaping the pigmentation phenotype and a complex evolutionary history involving strong selection. Polymorphisms in 2 genes, ASIP and OCA2, may play a shared role in shaping light and dark pigmentation across the globe, whereas SLC24A5, MATP, and TYR have a predominant role in the evolution of light skin in Europeans but not in East Asians. These findings support a case for the recent convergent evolution of a lighter pigmentation phenotype in Europeans and East Asians.

Skin and hair pigmentation variation in Island Melanesia

Skin and hair pigmentation are two of the most easily visible examples of human phenotypic variation. Selection-based explanations for pigmentation variation in humans have focused on the relationship between melanin and ultraviolet radiation, which is largely dependent on latitude. In this study, skin and hair pigmentation were measured as the melanin (M) index, using narrow-band reflectance spectroscopy for 1,135 individuals from Island Melanesia. Overall, the results show remarkable pigmentation variation, given the small geographic region surveyed. This variation is discussed in terms of differences between males and females, among islands, and among neighborhoods within those islands. The relationship of pigmentation to age, latitude, and longitude is also examined. We found that male skin pigmentation was significantly darker than females in 5 of 6 islands examined. Hair pigmentation showed a negative, but weak, correlation with age, while skin pigmentation showed a positive, but also weak, correlation with age. Skin and hair pigmentation varied significantly between islands as well as between neighborhoods within those islands. Bougainvilleans showed significantly darker skin than individuals from any other island considered, and are darker than a previously described African-American population. These findings are discussed in relation to prevailing hypotheses about the role of natural selection in shaping pigmentation variation in the human species, as well as the role of demographic processes such as admixture and drift in Island Melanesia.

SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans

Lighter variations of pigmentation in humans are associated with diminished number, size, and density of melanosomes, the pigmented organelles of melanocytes. Here we show that zebrafish golden mutants share these melanosomal changes and that golden encodes a putative cation exchanger slc24a5 (nckx5) that localizes to an intracellular membrane, likely the melanosome or its precursor. The human ortholog is highly similar in sequence and functional in zebrafish. The evolutionarily conserved ancestral allele of a human coding polymorphism predominates in African and East Asian populations. In contrast, the variant allele is nearly fixed in European populations, is associated with a substantial reduction in regional heterozygosity, and correlates with lighter skin pigmentation in admixed populations, suggesting a key role for the SLC24A5 gene in human pigmentation.

The 8818G allele of the agouti signaling protein (ASIP) gene is ancestral and is associated with darker skin color in African Americans

Skin color, a predictor of social interactions and risk factor for several types of cancer, is due to two contrasting forms of melanin, the darker eumelanin and lighter phaeomelanin. The lighter pigment phaeomelanin is the product of the antagonistic function of the agouti signaling protein (ASIP) on the alpha-melanocyte stimulating hormone receptor (MC1R). Studies have shown that a single-nucleotide polymorphism (SNP) in the 3'UTR of the ASIP gene is associated with dark hair and eyes; however, little is known about its role in inter-individual variation in skin color. Here we examine the relationship between the ASIP g.8818A>G SNP and skin color (M index) as assessed by reflectometry in 234 African Americans. Analyses of variance (ANOVA) were performed to evaluate the effects of ASIP genotypes, age, individual ancestry, and sex on skin color variation. Significant effects on M index variation were observed for ASIP genotypes (F(2,236)=4.37, P=0.01), ancestry (F(1,243)=37.2, P<0.001), and sex (F(1,244)=4.08, P=0.05). Subsequent analyses revealed a strong effect on M index from ASIP genotypes in African American females (P<0.001). Our study suggests that the ASIP G>A polymorphism exhibits a dominant effect leading to lighter skin color and that variation in the ASIP gene may have been one of several factors contributing to reductions in pigmentation in some populations. Further study is needed to reveal how interactions between ASIP and several other genes, such as MC1R and P, predict human pigmentation.

Comparing quantitative measures of erythema, pigmentation and skin response using reflectometry

We measured a number of pigmentation and skin response phenotypes in a sample of volunteers (n=397) living in State College, PA. The majority of this sample was composed of four groups based on stated ancestry: African-American, European-American, Hispanic and East Asian. Several measures of melanin concentration (L*, melanin index and adjusted melanin index) were estimated by diffuse reflectance spectroscopy and compared. The efficacy of these measures for assessing constitutive pigmentation and melanogenic dose-response was evaluated. Similarly, several measures of erythema (a*, erythema index and adjusted erythema index) were compared and evaluated in their efficacy in measuring erythema and erythemal dose-response. We show a high correspondence among all of the measures for the assessment of constitutive pigmentation and baseline erythema. However, our results demonstrate that evaluating melanogenic dose-response is highly dependent on the summary statistic used: while L* is a valid measure of constitutive pigmentation it is not an effective measure of melanogenic dose-response. Our results also confirm the use of a*, as it is shown to be highly correlated with the adjusted erythema index, a more advanced measure of erythema based on the apparent absorbance. Diffuse reflectance spectroscopy can be used to quantify the constitutive pigmentation, melanogenic dose-response at 7 d and erythemal dose-response at both 24 h and 7 d postexposure.