Genome-wide association studies of quantitatively measured skin, hair, and eye pigmentation in four European populations

Predictive accuracy of genetic variants for eye color in a Kazakh population using the IrisPlex system

OBJECTIVE

This study assesses the accuracy of the IrisPlex system, a genetic eye color prediction tool for forensic analysis, in the Kazakh population. The study compares previously published genotypes of 515 Kazakh individuals from varied geographical and ethnohistorical contexts with phenotypic data on their eye color, introduced for the first time in this research.

RESULTS

The IrisPlex panel's effectiveness in predicting eye color in the Kazakh population was validated. It exhibited slightly lower accuracy than in Western European populations but was higher than in Siberian populations. The sensitivity was notably high for brown-eyed individuals (0.99), but further research is needed for blue and intermediate eye colors. This study establishes IrisPlex as a useful predictive tool in the Kazakh population and provides a basis for future investigations into the genetic basis of phenotypic variations in this diverse population.

Skin colour: A window into human phenotypic evolution and environmental adaptation

As modern humans ventured out of Africa and dispersed around the world, they faced novel environmental challenges that led to geographic adaptations including skin colour. Over the long history of human evolution, skin colour has changed dramatically, showing tremendous diversity across different geographical regions, for example, the majority of individuals from the expansive lands of Africa have darker skin, whereas the majority of people from Eurasia exhibit lighter skin. What adaptations did lighter skin confer upon modern humans as they migrated from Africa to Eurasia? What genetic mechanisms underlie the diversity of skin colour observed in different populations? In recent years, scientists have gradually gained a deeper understanding of the interactions between pigmentation gene and skin colour through population-based genomic studies of different groups around the world, particularly in East Asia and Africa. In this review, we summarize our current understanding of 26 skin colour-related pigmentation genes and 48 SNPs that influence skin colour. Important pigmentation genes across three major populations are described in detail: MFSD12, SLC24A5, PDPK1 and DDB1/CYB561A3/TMEM138 influence skin colour in African populations; OCA2, KITLG, SLC24A2, GNPAT and PAH are key to the evolution of skin pigmentation in East Asian populations; and SLC24A5, SLC45A2, TYR, TYRP1, ASIP, MC1R and IRF4 significantly contribute to the lightening of skin colour in European populations. We summarized recent findings in genomic studies of skin colour in populations that implicate diverse geographic environments, local adaptation among populations, gene flow and multi-gene interactions as factors influencing skin colour diversity.

Melanin-based plumage coloration and melanin content in organs in the barn owl

Although the evolutionary ecology of melanin pigments and melanin-based coloration has been studied in great details, particularly in birds, little is known about the function of melanin stored inside the body. In the barn owl Tyto alba, in which individuals vary in the degree of reddish pheomelanin-based coloration and in the size of black eumelanic feather spots, we measured the concentration in melanin pigments in seven organs. The eyes had by far the most melanin then the skin, pectoral muscle, heart, liver, trachea, and uropygial gland. The concentration in eumelanin was not necessarily correlated with the concentration in pheomelanin suggesting that their production can be regulated independently from each other. Redder barn owls had more pheomelanin in the skin and uropygial gland than white owls, while owls displaying larger black feather spots had more eumelanin in the skin than small-spotted owls. More data are required to evaluate whether melanin-based traits can evolve as an indirect response to selection exerted on melanin deposition in organs.

Forensic DNA Phenotyping: Genes and Genetic Variants for Eye Color Prediction

In recent decades, the use of genetic polymorphisms related to specific phenotypes, such as eye color, has greatly contributed to the development of the research field called forensic DNA phenotyping (FDP), enabling the investigators of crime cases to reduce the number of suspects, making their work faster and more precise. Eye color is a polygenic phenotype, and many genetic variants have been highlighted, with the major contributor being the HERC2-OCA2 locus, where many single nucleotide variations (SNPs) were identified. Interestingly, the HERC2-OCA2 locus, containing the intronic SNP rs12913832, the major eye color determinant, shows a high level of evolutionary conservation across many species of vertebrates. Currently, there are some genetic panels to predict eye color by genomic DNA analysis, even if the exact role of the SNP variants in the formation of eye color is still poorly understood, with a low level of predictivity in the so-called intermediate eye color. Many variants in OCA2, HERC2, and other genes lie in introns or correspond to synonymous variants, highlighting greater complexity in the mechanism of action of such genes than a simple missense variation. Here, we show the main genes involved in oculocutaneous pigmentation and their structural and functional features, as well as which genetic variants show the highest level of eye color predictivity in currently used FDP assays. Despite the great recent advances and impact of FDP in criminal cases, it is necessary to enhance scientific research to better understand the mechanism of action behind each genetic variant involved in eye color, with the goal of obtaining higher levels of prediction.

Exploring regional aspects of 3D facial variation within European individuals

Facial ancestry can be described as variation that exists in facial features that are shared amongst members of a population due to environmental and genetic effects. Even within Europe, faces vary among subregions and may lead to confounding in genetic association studies if unaccounted for. Genetic studies use genetic principal components (PCs) to describe facial ancestry to circumvent this issue. Yet the phenotypic effect of these genetic PCs on the face has yet to be described, and phenotype-based alternatives compared. In anthropological studies, consensus faces are utilized as they depict a phenotypic, not genetic, ancestry effect. In this study, we explored the effects of regional differences on facial ancestry in 744 Europeans using genetic and anthropological approaches. Both showed similar ancestry effects between subgroups, localized mainly to the forehead, nose, and chin. Consensus faces explained the variation seen in only the first three genetic PCs, differing more in magnitude than shape change. Here we show only minor differences between the two methods and discuss a combined approach as a possible alternative for facial scan correction that is less cohort dependent, more replicable, non-linear, and can be made open access for use across research groups, enhancing future studies in this field.

Population designations in biomedical research: Limitations and perspectives

In biomedical research, population differences are of central interest. Variations in the frequency and severity of diseases and in treatment effects among human subpopulation groups are common in many medical conditions. Unfortunately, the practices in terms of subpopulation labeling do not exhibit the level of rigor one would expect in biomedical research, especially when studying multifactorial diseases such as cancer or atherosclerosis. The reporting of population differences in clinical research is characterized by large disparities in practices, and fraught with methodological issues and inconsistencies. The actual designations such as "Black" or "Asian" refer to broad and heterogeneous groups, with a great discrepancy among countries. Moreover, the use of obsolete concepts such as "Caucasian" is unfortunate and imprecise. The use of adequate labeling to reflect the scientific hypothesis needs to be promoted. Furthermore, the use of "race/ethnicity" as a unique cause of human heterogeneity may distract from investigating other factors related to a medical condition, particularly if this label is employed as a proxy for cultural habits, diet, or environmental exposure. In addition, the wide range of opinions among researchers does not facilitate the attempts made for resolving this heterogeneity in labeling. "Race," "ethnicity," "ancestry," "geographical origin," and other similar concepts are saturated with meanings. Even if the feasibility of a global consensus on labeling seems difficult, geneticists, sociologists, anthropologists, and ethicists should help develop policies and practices for the biomedical field.

Evaluation of Personalized Skincare Through in-silico Gene Interactive Networks and Cellular Responses to UVR and Oxidative Stress

Purpose

Personalized approaches in dermatology are designed to match the specific requirements based on the individual genetic makeup. One major factor accounting for the differences in skin phenotypes is single nucleotide polymorphism (SNP) within several genes with diverse roles that extend beyond skin tone and pigmentation. Therefore, the cellular sensitivities to the environmental stress and damage linked to extrinsic aging could also underlie the individual characteristics of the skin and dictate the unique skin care requirements. This study aimed to identify the likely biomarkers and molecular signatures expressed in skin cells of different ethnic backgrounds, which could aid further the design of personalized skin products based on specific demands.

Methods

Using data mining and in-silico modeling, the association of SNP-affected genes with three major skin types of European, Asian and African origin was analyzed and compared within the structure-function gene interaction networks. Cultured dermal fibroblasts were subsequently subjected to ultraviolet radiation and oxidative stress and analyzed for DNA damage and senescent markers. The protective applications of two cosmetic ingredients, Resveratrol and Quercetin, were validated in both cellular and in-silico models.

Results

Each skin type was characterized by the presence of SNPs in the genes controlling facultative and constitutive pigmentation, which could also underlie the major differences in responses to photodamage, such as oxidative stress, inflammation, and barrier homeostasis. Skin-type-specific dermal fibroblasts cultured in-vitro demonstrated distinctive sensitivities to ultraviolet radiation and oxidative stress, which could be modulated further by the bioactive compounds with the predicted capacities to interact with some of the genes in the in-silico models.

Conclusion

Evaluation of the SNP-affected gene networks and likely sensitivities of skin cells, defined as low threshold levels to extrinsic stress factors, can provide a valuable tool for the design and formulation of personalized skin products that match more accurately diverse ethnic backgrounds.

Genetic adaptation of skin pigmentation in highland Tibetans

Strong ultraviolet (UV) radiation at high altitude imposes a serious selective pressure, which may induce skin pigmentation adaptation of indigenous populations. We conducted skin pigmentation phenotyping and genome-wide analysis of Tibetans in order to understand the underlying mechanism of adaptation to UV radiation. We observe that Tibetans have darker baseline skin color compared with lowland Han Chinese, as well as an improved tanning ability, suggesting a two-level adaptation to boost their melanin production. A genome-wide search for the responsible genes identifies GNPAT showing strong signals of positive selection in Tibetans. An enhancer mutation (rs75356281) located in GNPAT intron 2 is enriched in Tibetans (58%) but rare in other world populations (0 to 18%). The adaptive allele of rs75356281 is associated with darker skin in Tibetans and, under UVB treatment, it displays higher enhancer activities compared with the wild-type allele in in vitro luciferase assays. Transcriptome analyses of gene-edited cells clearly show that with UVB treatment, the adaptive variant of GNPAT promotes melanin synthesis, likely through the interactions of CAT and ACAA1 in peroxisomes with other pigmentation genes, and they act synergistically, leading to an improved tanning ability in Tibetans for UV protection.

Investigating the genetic architecture of eye colour in a Canadian cohort

Eye color is highly variable in populations with European ancestry, ranging from low to high quantities of melanin in the iris. Polymorphisms in the HERC2/OCA2 locus have the largest effect on eye color in these populations, although other genomic regions also influence eye color. We performed genome-wide association studies of eye color in a Canadian cohort of European ancestry (N = 5,641) and investigated candidate causal variants. We uncovered several candidate causal signals in the HERC2/OCA2 region, whereas other loci likely harbor a single causal signal. We observed colocalization of eye color signals with the expression or methylation profiles of cultured primary melanocytes. Genetic correlations of eye and hair color suggest high genome-wide pleiotropy, but locus-level differences in the genetic architecture of both traits. Overall, we provide a better picture of the polymorphisms underpinning eye color variation, which may be a consequence of specific molecular processes in the iris melanocytes.

•

Genome-wide association studies of eye color in 5,641 participants

•

Multiple independent candidate causal variants were identified across HERC2/OCA2

•

Single candidate causal variants observed on or near IRF4, SLC24A4, TYR, and TYRP1

•

Colocalization of eye color signals with expression and methylation profiles

Advances in integrative African genomics

There has been a rapid increase in human genome sequencing in the past two decades, resulting in the identification of millions of previously unknown genetic variants. However, African populations are under-represented in sequencing efforts. Additional sequencing from diverse African populations and the construction of African-specific reference genomes is needed to better characterize the full spectrum of variation in humans. However, sequencing alone is insufficient to address the molecular and cellular mechanisms underlying variable phenotypes and disease risks. Determining functional consequences of genetic variation using multi-omics approaches is a fundamental post-genomic challenge. We discuss approaches to close the knowledge gaps about African genomic diversity and review advances in African integrative genomic studies and their implications for precision medicine.

Evolutionary selection of alleles in the melanophilin gene that impacts on prostate organ function and cancer risk

Background and objectives

Several hundred inherited genetic variants or SNPs that alter the risk of cancer have been identified through genome-wide association studies. In populations of European ancestry, these variants are mostly present at relatively high frequencies. To gain insight into evolutionary origins, we screened a series of genes and SNPs linked to breast or prostate cancer for signatures of historical positive selection.

Methodology

We took advantage of the availability of the 1000 genome data and we performed genomic scans for positive selection in five different Caucasian populations as well as one African reference population. We then used prostate organoid cultures to provide a possible functional explanation for the interplay between the action of evolutionary forces and the disease risk association.

Results

Variants in only one gene showed genomic signatures of positive, evolutionary selection within Caucasian populations melanophilin (MLPH). Functional depletion of MLPH in prostate organoids, by CRISPR/Cas9 mutation, impacted lineage commitment of progenitor cells promoting luminal versus basal cell differentiation and on resistance to androgen deprivation.

Conclusions and implications

The MLPH variants influencing prostate cancer risk may have been historically selected for their adaptive benefit on skin pigmentation but MLPH is highly expressed in the prostate and the derivative, positively selected, alleles decrease the risk of prostate cancer. Our study suggests a potential functional mechanism via which MLPH and its genetic variants could influence risk of prostate cancer, as a serendipitous consequence of prior evolutionary benefits to another tissue.

Lay Summary

We screened a limited series of genomic variants associated with breast and prostate cancer risk for signatures of historical positive selection. Variants within the melanophilin (MLPH) gene fell into this category. Depletion of MLPH in prostate organoid cultures, suggested a potential functional mechanism for impacting on cancer risk, as a serendipitous consequence of prior evolutionary benefits to another tissue.

GenNet framework: interpretable deep learning for predicting phenotypes from genetic data

Applying deep learning in population genomics is challenging because of computational issues and lack of interpretable models. Here, we propose GenNet, a novel open-source deep learning framework for predicting phenotypes from genetic variants. In this framework, interpretable and memory-efficient neural network architectures are constructed by embedding biologically knowledge from public databases, resulting in neural networks that contain only biologically plausible connections. We applied the framework to seventeen phenotypes and found well-replicated genes such as HERC2 and OCA2 for hair and eye color, and novel genes such as ZNF773 and PCNT for schizophrenia. Additionally, the framework identified ubiquitin mediated proteolysis, endocrine system and viral infectious diseases as most predictive biological pathways for schizophrenia. GenNet is a freely available, end-to-end deep learning framework that allows researchers to develop and use interpretable neural networks to obtain novel insights into the genetic architecture of complex traits and diseases.

van Hilten and colleagues present GenNet, a deep-learning framework for predicting phenotype from genetic data. This framework generates interpretable neural networks that provide insight into the genetic basis of complex traits and diseases.

Association of pigmentation related-genes polymorphisms and geographic environmental variables in the Chinese population

Background

Human skin color is highly heritable and one of the most variable phenotypic traits. However, the genetic causes and environmental selective pressures underlying this phenotypic variation have remained largely unknown. To investigate whether the pigmentation related-genes polymorphisms are associated with the geographic environmental variables. We selected randomly 795 healthy individuals from eight ethnic groups in nine provinces in China. Six single nucleotide polymorphisms (SNPs) of SLC45A2 and TYR were genotyped using Agena MassARRAY. The Chi-square test and Spearman correlation analysis were used to compare the frequency distribution of genotypes among different ethnic groups and evaluate the relationship between SNP genetic diversity and environmental variables, respectively.

Results

The results indicated that rs28777 and rs183671 (SLC45A2) and rs1042602 (TYR) genotype frequency distributions were significantly different between the Xinjiang-Uighur and other ethnic groups (P < 0.05). Spearman correlation analysis found that rs28777-A (r = − 0.090, P = 0.011), rs183671-G (r = − 0.105, P = 0.003), rs1042602-A (r = − 0.108, P = 0.002), rs1126809-A (r = − 0.151, P < 0.001) allele frequencies were negatively correlated with the longitude; rs183671-G (r = 0.151), rs1042602-A (r = 0.157) and rs1126809-A (r = 0.138) allele frequencies were positively associated with the latitude (P < 0.001); rs183671-G (r = 0.116, P = 0.001), rs1042602-A (r = 0.105, P = 0.003) and rs1126809-A (r = 0.070, P = 0.048) allele frequencies were positively correlated with the sunshine hours; rs183671-G (r = − 0.076, P = 0.033), rs1042602-A (r = − 0.079, P = 0.027) and rs1126809-A (r = − 0.076, P = 0.031) were negatively correlated with the annual average temperature.

Conclusions

Our results confirmed the idea that environmental factors have been an important selective pressure upon pigmentation related gene polymorphisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-021-00189-7.

The evolution of skin pigmentation-associated variation in West Eurasia

Some of the genes responsible for the evolution of light skin pigmentation in Europeans show signals of positive selection in present-day populations. Recently, genome-wide association studies have highlighted the highly polygenic nature of skin pigmentation. It is unclear whether selection has operated on all of these genetic variants or just a subset. By studying variation in over a thousand ancient genomes from West Eurasia covering 40,000 y, we are able to study both the aggregate behavior of pigmentation-associated variants and the evolutionary history of individual variants. We find that the evolution of light skin pigmentation in Europeans was driven by frequency changes in a relatively small fraction of the genetic variants that are associated with variation in the trait today.

Skin pigmentation is a classic example of a polygenic trait that has experienced directional selection in humans. Genome-wide association studies have identified well over a hundred pigmentation-associated loci, and genomic scans in present-day and ancient populations have identified selective sweeps for a small number of light pigmentation-associated alleles in Europeans. It is unclear whether selection has operated on all of the genetic variation associated with skin pigmentation as opposed to just a small number of large-effect variants. Here, we address this question using ancient DNA from 1,158 individuals from West Eurasia covering a period of 40,000 y combined with genome-wide association summary statistics from the UK Biobank. We find a robust signal of directional selection in ancient West Eurasians on 170 skin pigmentation-associated variants ascertained in the UK Biobank. However, we also show that this signal is driven by a limited number of large-effect variants. Consistent with this observation, we find that a polygenic selection test in present-day populations fails to detect selection with the full set of variants. Our data allow us to disentangle the effects of admixture and selection. Most notably, a large-effect variant at SLC24A5 was introduced to Western Europe by migrations of Neolithic farming populations but continued to be under selection post-admixture. This study shows that the response to selection for light skin pigmentation in West Eurasia was driven by a relatively small proportion of the variants that are associated with present-day phenotypic variation.

Prediction of eye, hair and skin colour in Latin Americans

Here we evaluate the accuracy of prediction for eye, hair and skin pigmentation in a dataset of > 6500 individuals from Mexico, Colombia, Peru, Chile and Brazil (including genome-wide SNP data and quantitative/categorical pigmentation phenotypes - the CANDELA dataset CAN). We evaluated accuracy in relation to different analytical methods and various phenotypic predictors. As expected from statistical principles, we observe that quantitative traits are more sensitive to changes in the prediction models than categorical traits. We find that Random Forest or Linear Regression are generally the best performing methods. We also compare the prediction accuracy of SNP sets defined in the CAN dataset (including 56, 101 and 120 SNPs for eye, hair and skin colour prediction, respectively) to the well-established HIrisPlex-S SNP set (including 6, 22 and 36 SNPs for eye, hair and skin colour prediction respectively). When training prediction models on the CAN data, we observe remarkably similar performances for HIrisPlex-S and the larger CAN SNP sets for the prediction of hair (categorical) and eye (both categorical and quantitative), while the CAN sets outperform HIrisPlex-S for quantitative, but not for categorical skin pigmentation prediction. The performance of HIrisPlex-S, when models are trained in a world-wide sample (although consisting of 80% Europeans, https://hirisplex.erasmusmc.nl), is lower relative to training in the CAN data (particularly for hair and skin colour). Altogether, our observations are consistent with common variation of eye and hair colour having a relatively simple genetic architecture, which is well captured by HIrisPlex-S, even in admixed Latin Americans (with partial European ancestry). By contrast, since skin pigmentation is a more polygenic trait, accuracy is more sensitive to prediction SNP set size, although here this effect was only apparent for a quantitative measure of skin pigmentation. Our results support the use of HIrisPlex-S in the prediction of categorical pigmentation traits for forensic purposes in Latin America, while illustrating the impact of training datasets on its accuracy.

Evaluation of supervised machine-learning methods for predicting appearance traits from DNA

The prediction of human externally visible characteristics (EVCs) based solely on DNA information has become an established approach in forensic and anthropological genetics in recent years. While for a large set of EVCs, predictive models have already been established using multinomial logistic regression (MLR), the prediction performances of other possible classification methods have not been thoroughly investigated thus far. Motivated by the question to identify a potential classifier that outperforms these specific trait models, we conducted a systematic comparison between the widely used MLR and three popular machine learning (ML) classifiers, namely support vector machines (SVM), random forest (RF) and artificial neural networks (ANN), that have shown good performance outside EVC prediction. As examples, we used eye, hair and skin color categories as phenotypes and genotypes based on the previously established IrisPlex, HIrisPlex, and HIrisPlex-S DNA markers. We compared and assessed the performances of each of the four methods, complemented by detailed hyperparameter tuning that was applied to some of the methods in order to maximize their performance. Overall, we observed that all four classification methods showed rather similar performance, with no method being substantially superior to the others for any of the traits, although performances varied slightly across the different traits and more so across the trait categories. Hence, based on our findings, none of the ML methods applied here provide any advantage on appearance prediction, at least when it comes to the categorical pigmentation traits and the selected DNA markers used here.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

Genetic loci associated with skin pigmentation in African Americans and their effects on vitamin D deficiency

A recent genome-wide association study (GWAS) in African descent populations identified novel loci associated with skin pigmentation. However, how genomic variations affect skin pigmentation and how these skin pigmentation gene variants affect serum 25(OH) vitamin D variation has not been explored in African Americans (AAs). In order to further understand genetic factors that affect human skin pigmentation and serum 25(OH)D variation, we performed a GWAS for skin pigmentation with 395 AAs and a replication study with 681 AAs. Then, we tested if the identified variants are associated with serum 25(OH) D concentrations in a subset of AAs (n = 591). Skin pigmentation, Melanin Index (M-Index), was measured using a narrow-band reflectometer. Multiple regression analysis was performed to identify variants associated with M-Index and to assess their role in serum 25(OH)D variation adjusting for population stratification and relevant confounding variables. A variant near the SLC24A5 gene (rs2675345) showed the strongest signal of association with M-Index (P = 4.0 x 10-30 in the pooled dataset). Variants in SLC24A5, SLC45A2 and OCA2 together account for a large proportion of skin pigmentation variance (11%). The effects of these variants on M-Index was modified by sex (P for interaction = 0.009). However, West African Ancestry (WAA) also accounts for a large proportion of M-Index variance (23%). M-Index also varies among AAs with high WAA and high Genetic Score calculated from top variants associated with M-Index, suggesting that other unknown genomic factors related to WAA are likely contributing to skin pigmentation variation. M-Index was not associated with serum 25(OH)D concentrations, but the Genetic Score was significantly associated with vitamin D deficiency (serum 25(OH)D levels less than 12 ng/mL) (OR, 1.30; 95% CI, 1.04-1.64). The findings support the hypothesis suggesting that skin pigmentation evolved responding to increased demand for subcutaneous vitamin D synthesis in high latitude environments.

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.

Genetic ancestry, skin pigmentation, and the risk of cutaneous squamous cell carcinoma in Hispanic/Latino and non-Hispanic white populations

Although cutaneous squamous cell carcinoma (cSCC) is one of the most common malignancies in individuals of European ancestry, the incidence of cSCC in Hispanic/Latinos is also increasing. cSCC has both a genetic and environmental etiology. Here, we examine the role of genetic ancestry, skin pigmentation, and sun exposure in Hispanic/Latinos and non-Hispanic whites on cSCC risk. We observe an increased cSCC risk with greater European ancestry (P = 1.27 × 10⁻⁴²) within Hispanic/Latinos and with greater northern (P = 2.38 × 10⁻⁶⁵) and western (P = 2.28 × 10⁻⁴⁹) European ancestry within non-Hispanic whites. These associations are significantly, but not completely, attenuated after considering skin pigmentation-associated loci, history of actinic keratosis, and sun-protected versus sun-exposed anatomical sites. We also report an association of the well-known pigment variant Ala111Thr (rs1426654) at SLC24A5 with cSCC in Hispanic/Latinos. These findings demonstrate a strong correlation of northwestern European genetic ancestry with cSCC risk in both Hispanic/Latinos and non-Hispanic whites, largely but not entirely mediated through its impact on skin pigmentation.

Eric Jorgenson and Hélène Choquet et al. find that northwestern European genetic ancestry is associated with increased risk of cutaneous squamous cell carcinoma (cSCC) in non-Hispanic whites, and more so in Hispanic/Latinos of the US. The ancestry effect is largely, but not entirely explained by genetic determinants of skin pigmentation in both populations.

Optimizing the genetic prediction of the eye and hair color for North Eurasian populations

Background

Predicting the eye and hair color from genotype became an established and widely used tool in forensic genetics, as well as in studies of ancient human populations. However, the accuracy of this tool has been verified on the West and Central Europeans only, while populations from border regions between Europe and Asia (like Caucasus and Ural) also carry the light pigmentation phenotypes.

Results

We phenotyped 286 samples collected across North Eurasia, genotyped them by the standard HIrisPlex-S markers and found that predictive power in Caucasus/Ural/West Siberian populations is reasonable but lower than that in West Europeans. As these populations have genetic ancestries different from that of West Europeans, we hypothesized they may carry a somewhat different allele spectrum. Thus, for all samples we performed the exome sequencing additionally enriched with the 53 genes and intergenic regions known to be associated with the eye/hair color. Our association analysis replicated the importance of the key previously known SNPs but also identified five new markers whose eye color prediction power for the studied populations is compatible with the two major previously well-known SNPs. Four out of these five SNPs lie within the HERС2 gene and the fifth in the intergenic region. These SNPs are found at high frequencies in most studied populations. The released dataset of exomes from Russian populations can be further used for population genetic and medical genetic studies.

Conclusions

This study demonstrated that precision of the established systems for eye/hair color prediction from a genotype is slightly lower for the populations from the border regions between Europe and Asia that for the West Europeans. However, this precision can be improved if some newly revealed predictive SNPs are added into the panel. We discuss that the replication of these pigmentation-associated SNPs on the independent North Eurasian sample is needed in the future studies.

Insights on hair, skin and eye color of ancient and contemporary Native Americans

Over the past few years, tools capable of predicting pigmentation phenotypes have been developed aiming to contribute for criminal and anthropological investigations. In this study, we used eight genetic systems to infer eye, hair, and skin color of ancient and contemporary Native Americans. To achieve this goal, we retrieved 61 SNPs from 42 samples available in free online repositories of DNA sequences. We performed pigmentation predictions using two freely available tools, HIrisPlex-S and Snipper, in addition to two other published models. This workflow made possible to predict all three phenotypes with at least one tool for 29 out of the 42 samples. Considering these 29 individuals, predictions for eye, hair, and skin color were obtained with HIrisPlex-S for 27, 28 and 27 individuals, respectively, while 24, 25 and 25 individuals had such predictions with Snipper. In general, ancient and contemporary Native Americans were predicted to have intermediate/brown eyes, black hair, and intermediate/darker skin pigmentation.

A study in scarlet: MC1R as the main predictor of red hair and exemplar of the flip-flop effect

Genetic variation in melanocortin-1 receptor (MC1R) is a known contributor to disease-free red hair in humans. Three loss-of-function single-nucleotide variants (rs1805007, rs1805008 and rs1805009) have been established as strongly correlated with red hair. The contribution of other loss-of-function MC1R variants (in particular rs1805005, rs2228479 and rs885479) and the extent to which other genetic loci are involved in red hair colour is less well understood. Here, we used the UK Biobank cohort to capture a comprehensive list of MC1R variants contributing to red hair colour. We report a correlation with red hair for both strong-effect variants (rs1805007, rs1805008 and rs1805009) and weak-effect variants (rs1805005, rs2228479 and rs885479) and show that their coefficients differ by two orders of magnitude. On the haplotype level, both strong- and weak-effect variants contribute to the red hair phenotype, but when considered individually, weak-effect variants show a reverse, negative association with red hair. The reversal of association direction in the single-variant analysis is facilitated by a distinguishing structure of MC1R, in which loss-of-function variants are never found to co-occur on the same haplotype. The other previously reported hair colour genes' variants do not substantially improve the MC1R red hair colour predictive model. Our best model for predicting red versus other hair colours yields an unparalleled area under the receiver operating characteristic of 0.96 using only MC1R variants. In summary, we present a comprehensive statistically derived characterization of the role of MC1R variants in red hair colour and offer a powerful, economical and parsimonious model that achieves unsurpassed performance.

The accuracy of predicting eye and hair pigmentation based on genetic markers in Russian populations

Prediction of eye and hair color from DNA is being increasingly employed in forensics and the studies of ancient populations. HIrisPlex-S is a prediction tool trained on the Dutch dataset and verified using a few other European populations. The accuracy of its predictions for other regions of the world has not been studied yet. Russian populations pose a special interest because, unlike the majority of world populations, their representatives can have not only dark but also light color eyes and hair. The aim of this work was to evaluate the accuracy of eye and hair color prediction in Russian populations. We recruited 144 representatives of indigenous peoples of Russia (Avar, Aleut, Buryat, Itelmen, Karelian, Koryak, Mari, Nanai, Russian, Rutulian, Chuvash, Chukchi, Evenk, and Even populations). All study participants were photographed. Eye and hair colors were identified from the anthropological images by anthropologists. The SNP markers included in the HIrisPlex system were genotyped. Phenotypes were predicted from the obtained genotypes and subsequently compared to the actual phenotypes. Quality metrics were calculated for HIrisPlex prediction accuracy in the populations of European Russia and Siberia. On the whole, HIrisPlex prediction accuracy was satisfactory, although a bit lower than in Western European datasets. Further research could focus on identifying additional markers to improve the accuracy of predictions in Russian populations.

Contributions by MC1R Variants to Melanoma Risk in Males and Females

Importance

Recently, the red hair variants of MC1R were found to contribute differently to pigmentation phenotype in males and females.

Objective

To investigate the role of these variants in melanoma risk in males and females separately because carriers of the red hair variants of MC1R are at increased risk of melanoma.

Design, Setting, and Participants

In this hospital-based, case-control study, we evaluated the effect of MC1R and melanoma risk for males and females separately by performing multivariate logistic regression analyses.

Main Outcomes and Measures

Association of MC1R variants and melanoma risk in males and females.

Results

A total of 905 females (473 melanoma cases, 432 controls) and 886 males (518 melanoma cases, 368 controls) were included in the analyses. The mean (SD) age of the study population was 59.2 (15.6). In females, carrying any MC1R red hair variants remained an independent risk factor of melanoma in a multivariable analysis (adjusted odds ratio [OR], 2.19 [95% CI, 1.60-2.99]), whereas in males, only signs of actinic skin damage (lentigines on the back [OR, 2.56; 95% CI, 1.47-4.45; P = .001] and the hands [OR, 2.31; 95% CI, 1.24-4.29; P = .008] and wrinkling on the neck [OR, 2.17; 95% CI, 1.23-3.82; P = .007]) and sunburns (OR, 1.65; 95% CI, 1.12-2.42; P = .01) remained significant risk factors.

Conclusions and Relevance

MC1R variants contribute differently to melanoma risk in males and females. This could be helpful to better classify melanoma risk factors between the sexes.

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.

Genome-Wide Characterization of Arabian Peninsula Populations: Shedding Light on the History of a Fundamental Bridge between Continents

The Arabian Peninsula (AP) was an important crossroad between Africa, Asia, and Europe, being the cradle of the structure defining these main human population groups, and a continuing path for their admixture. The screening of 741,000 variants in 420 Arabians and 80 Iranians allowed us to quantify the dominant sub-Saharan African admixture in the west of the peninsula, whereas South Asian and Levantine/European influence was stronger in the east, leading to a rift between western and eastern sides of the Peninsula. Dating of the admixture events indicated that Indian Ocean slave trade and Islamization periods were important moments in the genetic makeup of the region. The western-eastern axis was also observable in terms of positive selection of diversity conferring lactose tolerance, with the West AP developing local adaptation and the East AP acquiring the derived allele selected in European populations and existing in South Asia. African selected malaria resistance through the DARC gene was enriched in all Arabian genomes, especially in the western part. Clear European influences associated with skin and eye color were equally frequent across the Peninsula.

Meta-analysis and prioritization of human skin pigmentation-associated GWAS-SNPs using ENCODE data-based web-tools

Skin pigmentation in human is a complex trait, which varies widely, both within and between human populations. The exact players governing the trait of skin pigmentation remain elusive till date. Various Genome Wide Association Studies (GWAS) have shown the association of different genomic variants with normal human skin pigmentation, often indicating genes with no direct implications in melanin biosynthesis or distribution. Little has been explained in terms of the functionality of the associated Single-Nucleotide Polymorphisms (SNPs) with respect to modulating the skin pigmentation phenotype. In the present study, which, to our knowledge, is the first of its kind, we tried to analyze and prioritize 519 non-coding SNPs and 24 3'UTR SNPs emerging from 14 different human skin pigmentation-related GWAS, primarily using several ENCODE-based web-tools like rSNPBase, RegulomeDB, HaploReg, etc., most of which incorporate experimentally validated evidences in their predictions. Using this comprehensive, in-silico, analytical approach, we successfully prioritized all the pigmentation-associated GWAS-SNPs and tried to annotate pigmentation-related functionality to them, which would pave the way for deeper understanding of the molecular basis of human skin pigmentation variations.

Pigment genes not skin pigmentation affect UVB-induced vitamin D

The long-term UVB induced serum 25(OH)D increase was influenced by pigment gene polymorphisms rather than measured facultative skin pigmentation.

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.

Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4 Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.

Clinical and Biological Characterization of Skin Pigmentation Diversity and Its Consequences on UV Impact

Skin color diversity is the most variable and noticeable phenotypic trait in humans resulting from constitutive pigmentation variability. This paper will review the characterization of skin pigmentation diversity with a focus on the most recent data on the genetic basis of skin pigmentation, and the various methodologies for skin color assessment. Then, melanocyte activity and amount, type and distribution of melanins, which are the main drivers for skin pigmentation, are described. Paracrine regulators of melanocyte microenvironment are also discussed. Skin response to sun exposure is also highly dependent on color diversity. Thus, sensitivity to solar wavelengths is examined in terms of acute effects such as sunburn/erythema or induced-pigmentation but also long-term consequences such as skin cancers, photoageing and pigmentary disorders. More pronounced sun-sensitivity in lighter or darker skin types depending on the detrimental effects and involved wavelengths is reviewed.

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.

Possible Role of Single Stranded DNA Binding Protein 3 on Skin Hydration by Regulating Epidermal Differentiation

Background

Skin hydration is a common problem both in elderly and young people as dry skin may cause irritation, dermatological disorders, and wrinkles. While both genetic and environmental factors seem to influence skin hydration, thorough genetic studies on skin hydration have not yet been conducted.

Objective

We used a genome-wide association study (GWAS) to explore the genetic elements underlying skin hydration by regulating epidermal differentiation and skin barrier function.

Methods

A GWAS was conducted to investigate the genetic factors influencing skin hydration in 100 Korean females along with molecular studies of genes in human epidermal keratinocytes for functional study in vitro.

Results

Among several single nucleotide polymorphisms identified in GWAS, we focused on Single Stranded DNA Binding Protein 3 (SSBP3) which is associated with DNA replication and DNA damage repair. To better understand the role of SSBP3 in skin cells, we introduced a calcium-induced differentiation keratinocyte culture system model and found that SSBP3 was upregulated in keratinocytes in a differentiation dependent manner. When SSBP3 was overexpressed using a recombinant adenovirus, the expression of differentiation-related genes such as loricrin and involucrin was markedly increased.

Conclusion

Taken together, our results suggest that genetic variants in the intronic region of SSBP3 could be determinants in skin hydration of Korean females. SSBP3 represents a new candidate gene to evaluate the molecular basis of the hydration ability in individuals.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

The iris signal: blue periphery, tan collaret and freckles pattern - strong indicators for epidermal skin cancer in South-Eastern Europe

BACKGROUND

Eye and skin share the embryological origin. Both are established risk factors in epidermal skin cancer. There are few reports using iris colour classification scales, most of them analyse colour in general or are too complex to use in daily practice.

OBJECTIVES

To investigate which iris colour pattern is associated with epidermal skin cancer in a S-E European Caucasian population.

METHODS

A case-control study was conducted on 480 patients: 229 skin cancers patients and 251 controls (dermatological patients free of skin cancers) admitted in two medical clinics of Dermatology in Bucharest, between October 2011 and May 2014. High-resolution iris photographs were taken for each patient. Three parameters of the iris were analysed individually and in association patterns for each patient: periphery, collaret and freckles.

RESULTS

The most frequent iris colour pattern associated with epidermal skin cancer was blue periphery with light brown collaret and freckles present. In terms of individual parameters, the strongest indicators for skin cancer patients were blue periphery and blue collaret.

CONCLUSIONS

The results of this study sustain the hypothesis that blue periphery with light brown collaret and freckles iris pattern is a reliable phenotypic marker for epidermal skin cancer. The results of this study differ from previous reports in which skin cancer risk was associated with a homogeneous blue iris. We account these differences in the characteristics of the recruited patients (S-E European, skin type II and III). The assessment of iris colour patterns is an easy and inexpensive detection tool in skin cancer risk assessment.

An Unexpectedly Complex Architecture for Skin Pigmentation in Africans

Approximately 15 genes have been directly associated with skin pigmentation variation in humans, leading to its characterization as a relatively simple trait. However, by assembling a global survey of quantitative skin pigmentation phenotypes, we demonstrate that pigmentation is more complex than previously assumed, with genetic architecture varying by latitude. We investigate polygenicity in the KhoeSan populations indigenous to southern Africa who have considerably lighter skin than equatorial Africans. We demonstrate that skin pigmentation is highly heritable, but known pigmentation loci explain only a small fraction of the variance. Rather, baseline skin pigmentation is a complex, polygenic trait in the KhoeSan. Despite this, we identify canonical and non-canonical skin pigmentation loci, including near SLC24A5, TYRP1, SMARCA2/VLDLR, and SNX13, using a genome-wide association approach complemented by targeted resequencing. By considering diverse, under-studied African populations, we show how the architecture of skin pigmentation can vary across humans subject to different local evolutionary pressures.

The Genetic Overlap Between Hair and Eye Color

We identified the genetic variants for eye color by Genome-Wide Association Study (GWAS) in a Dutch Caucasian family-based population sample and examined the genetic correlation between hair and eye color using data from unrelated participants from the Netherlands Twin Register. With the Genome-wide Complex Trait Analysis software package, we found strong genetic correlations between various combinations of hair and eye colors. The strongest positive correlations were found for blue eyes with blond hair (0.87) and brown eyes with dark hair (0.71), whereas blue eyes with dark hair and brown eyes with blond hair showed the strongest negative correlations (-0.64 and -0.94, respectively). Red hair with green/hazel eyes showed the weakest correlation (-0.14). All analyses were corrected for age and sex, and we explored the effects of correcting for principal components (PCs) that represent ancestry and describe the genetic stratification of the Netherlands. When including the first three PCs as covariates, the genetic correlations between the phenotypes disappeared. This is not unexpected since hair and eye colors strongly indicate the ancestry of an individual. This makes it difficult to separate the effects of population stratification and the true genetic effects of variants on these particular phenotypes.

Comparative Transcriptome Analysis of Mink (Neovison vison) Skin Reveals the Key Genes Involved in the Melanogenesis of Black and White Coat Colour

Farmed mink (Neovison vison) is one of the most important fur-bearing species worldwide, and coat colour is a crucial qualitative characteristic that contributes to the economic value of the fur. To identify additional genes that may play important roles in coat colour regulation, Illumina/Solexa high-throughput sequencing technology was used to catalogue the global gene expression profiles in mink skin with two different coat colours (black and white). RNA-seq analysis indicated that a total of 12,557 genes were differentially expressed in black versus white minks, with 3,530 genes up-regulated and 9,027 genes down-regulated in black minks. Significant differences were not observed in the expression of MC1R and TYR between the two different coat colours, and the expression of ASIP was not detected in the mink skin of either coat colour. The expression levels of KITLG, LEF1, DCT, TYRP1, PMEL, Myo5a, Rab27a and SLC7A11 were validated by qRT-PCR, and the results were consistent with RNA-seq analysis. This study provides several candidate genes that may be associated with the development of two coat colours in mink skin. These results will expand our understanding of the complex molecular mechanisms underlying skin physiology and melanogenesis in mink and will provide a foundation for future studies.

Variants at the OCA2/HERC2 locus affect time to first cutaneous squamous cell carcinoma in solid organ transplant recipients collected using two different study designs

BACKGROUND

Variants at the oculocutaneous albinism 2 (OCA2)/HECT and RLD domain containing E3 ubiquitin protein ligase 2 (HERC2) locus have been associated with pigmentation phenotypes and risk of developing several types of skin cancer.

OBJECTIVES

To evaluate OCA2/HERC2 locus variants for their impact on time to develop cutaneous squamous cell carcinoma (cSCC) in organ transplant recipients (OTRs) who are at elevated risk of developing cSCC.

METHODS

Participants were solid OTRs ascertained from two centres (n = 125 and 261) with an average of 13·1 years of follow-up post-transplant. DNA was available for genotyping for all participants, in addition to medical records and questionnaire data. The Ohio State University study had a case-control design with prospective follow-up, and the University of California San Francisco study was a national cross-sectional survey with retrospective chart review.

RESULTS

OCA2 variants rs12913832 and rs916977 were significantly associated with time to first cSCC post-transplant. OTRs homozygous for the brown-eye alleles of rs916977 (GG) and rs12913832 (AA) had significant delays of time to first cSCC post-transplant compared with individuals homozygous for the blue-eye alleles (hazard ratio 0·34, P < 0·001 and hazard ratio 0·54, P = 0·012, respectively). Both variants were highly associated with eye colour in the combined studies (P < 0·001).

CONCLUSIONS

This study is the first to show an association between OCA2/HERC2 variants and time to first cSCC post-transplant. This may impact dermatological screening recommendations for high-risk populations.

miR-128-1 is not required for hair pigmentation in mice

MicroRNAs are endogenous, regulatory RNAs implicated in many biological processes including pigmentation. Software algorithms and in vitro experiments predict that microRNAs can target pigmentation pathway genes, but few have been tested in vivo. MiR-128-1, a microRNA within the strongly selected lactase locus in the human genome, has predicted pigmentation targets. To test the role of miR-128-1 in pigment regulation, we created C57BL/6 agouti miR-128-1 knockout mice and quantified melanin deposition in hair. MiR-128-1 knockout mice have no detectable hair pigmentation phenotype. We conclude that miR-128-1 does not play a significant regulatory role in hair pigmentation in mice.

The colours of humanity: the evolution of pigmentation in the human lineage

Humans are a colourful species of primate, with human skin, hair and eye coloration having been influenced by a great variety of evolutionary forces throughout prehistory. Functionally naked skin has been the physical interface between the physical environment and the human body for most of the history of the genus Homo, and hence skin coloration has been under intense natural selection. From an original condition of protective, dark, eumelanin-enriched coloration in early tropical-dwelling Homo and Homo sapiens, loss of melanin pigmentation occurred under natural selection as Homo sapiens dispersed into non-tropical latitudes of Africa and Eurasia. Genes responsible for skin, hair and eye coloration appear to have been affected significantly by population bottlenecks in the course of Homo sapiens dispersals. Because specific skin colour phenotypes can be created by different combinations of skin colour-associated genetic markers, loss of genetic variability due to genetic drift appears to have had negligible effects on the highly redundant genetic 'palette' for the skin colour. This does not appear to have been the case for hair and eye coloration, however, and these traits appear to have been more strongly influenced by genetic drift and, possibly, sexual selection.This article is part of the themed issue 'Animal coloration: production, perception, function and application'.

Genetic 3'UTR variation is associated with human pigmentation characteristics and sensitivity to sunlight

Sunlight exposure induces signalling pathways leading to the activation of melanin synthesis and tanning response. MicroRNAs (miRNAs) can regulate the expression of genes involved in pigmentation pathways by binding to the complementary sequence in their 3'untranslated regions (3'UTRs). Therefore, 3'UTR SNPs are predicted to modify the ability of miRNAs to target genes, resulting in differential gene expression. In this study, we investigated the role in pigmentation and sun-sensitivity traits, as well as in melanoma susceptibility, of 38 different 3'UTR SNPs from 38 pigmentation-related genes. A total of 869 individuals of Spanish origin (526 melanoma cases and 343 controls) were analysed. The association of genotypic data with pigmentation traits was analysed via logistic regression. Web-based tools for predicting the effect of genetic variants in microRNA-binding sites in 3'UTR gene regions were also used. Seven 3'UTR SNPs showed a potential implication in melanoma risk phenotypes. This association is especially noticeable for two of them, rs2325813 in the MLPH gene and rs752107 in the WNT3A gene. These two SNPs were predicted to disrupt a miRNA-binding site and to impact on miRNA-mRNA interaction. To our knowledge, this is the first time that these two 3'UTR SNPs have been associated with sun-sensitivity traits. We state the potential implication of these SNPs in human pigmentation and sensitivity to sunlight, possibly as a result of changes in the level of gene expression through the disruption of putative miRNA-binding sites.

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.