OPRM1 and EGFR contribute to skin pigmentation differences between Indigenous Americans and Europeans

Global and Local Ancestry and its Importance: A Review

The fastest way to significantly change the composition of a population is through admixture, an evolutionary mechanism. In animal breeding history, genetic admixture has provided both short-term and long-term advantages by utilizing the phenomenon of complementarity and heterosis in several traits and genetic diversity, respectively. The traditional method of admixture analysis by pedigree records has now been replaced greatly by genome-wide marker data that enables more precise estimations. Among these markers, SNPs have been the popular choice since they are cost-effective, not so laborious, and automation of genotyping is easy. Certain markers can suggest the possibility of a population's origin from a sample of DNA where the source individual is unknown or unwilling to disclose their lineage, which are called Ancestry-Informative Markers (AIMs). Revealing admixture level at the locus-specific level is termed as local ancestry and can be exploited to identify signs of recent selective response and can account for genetic drift. Considering the importance of genetic admixture and local ancestry, in this mini-review, both concepts are illustrated, encompassing basics, their estimation/identification methods, tools/software used and their applications.

MC1R and age heteroclassification of face phenotypes in the Rio Grande do Sul population

BACKGROUND

Forensic DNA phenotyping (FDP) consists of the use of methodologies for predicting externally visible characteristics (EVCs) from the genetic material of biological samples found in crime scenes and has proven to be a promising tool in aiding human identification in police activities. Currently, methods based on multiplex assays and statistical models of prediction of EVCs related to hair, skin, and iris pigmentation using panels of SNP and INDEL biomarkers have already been developed and validated by the forensic scientific community. As well as traces of pigmentation, an individual's perceived age (PA) can also be considered an EVC and its estimation in unknown individuals can be useful for the progress of investigations. Liu and colleagues (2016) were pioneers in evidencing that, in addition to lifestyle and environmental factors, the presence of SNP and INDEL variants in the MC1R gene - which encodes a transmembrane receptor responsible for regulating melanin production - seems to contribute to an individual's PA. The group highlighted the association between these MC1R gene polymorphisms and the PA in the European population, where carriers of risk haplotypes appeared to be up to 2 years older in comparison to their chronological age (CA).

PURPOSE

Understanding that genotype-phenotype relationships cannot be extrapolated between different population groups, this study aimed to test this hypothesis and verify the applicability of this variant panel in the Rio Grande do Sul admixed population.

METHODS

Based on genomic data from a sample of 261 volunteers representative of gaucho population and using a multiple linear regression (MLR) model, our group was able to verify a significant association among nine intronic variants in loci adjacent to MC1R (e.g., AFG3L1P, TUBB3, FANCA) and facial age appearance, whose PA was defined after age heteroclassification of standard frontal face images through 11 assessors.

RESULTS

Different from that observed in European populations, our results show that the presence of effect alleles (R) of the selected variants in our sample influenced both younger and older face phenotypes. The influence of each variant on PA is expressed as β values.

CONCLUSIONS

There are important molecular mechanisms behind the effects of MC1R locus on PA, and the genomic background of each population seems to be crucial to determine this influence.

Deconstructing Eurocentrism in skin pigmentation research via the incorporation of diverse populations and theoretical perspectives

The evolution of skin pigmentation has been shaped by numerous biological and cultural shifts throughout human history. Vitamin D is considered a driver of depigmentation evolution in humans, given the deleterious health effects associated with vitamin D deficiency, which is often shaped by cultural factors. New advancements in genomics and epigenomics have opened the door to a deeper exploration of skin pigmentation evolution in both contemporary and ancient populations. Data from ancient Europeans has offered great context to the spread of depigmentation alleles via the evaluation of migration events and cultural shifts that occurred during the Neolithic. However, novel insights can further be gained via the inclusion of diverse ancient and contemporary populations. Here we present on how potential biases and limitations in skin pigmentation research can be overcome with the integration of interdisciplinary data that includes both cultural and biological elements, which have shaped the evolutionary history of skin pigmentation in humans.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Native American genetic ancestry and pigmentation allele contributions to skin color in a Caribbean population

Our interest in the genetic basis of skin color variation between populations led us to seek a Native American population with genetically African admixture but low frequency of European light skin alleles. Analysis of 458 genomes from individuals residing in the Kalinago Territory of the Commonwealth of Dominica showed approximately 55% Native American, 32% African, and 12% European genetic ancestry, the highest Native American genetic ancestry among Caribbean populations to date. Skin pigmentation ranged from 20 to 80 melanin units, averaging 46. Three albino individuals were determined to be homozygous for a causative multi-nucleotide polymorphism OCA2NW273KV contained within a haplotype of African origin; its allele frequency was 0.03 and single allele effect size was -8 melanin units. Derived allele frequencies of SLC24A5A111T and SLC45A2L374F were 0.14 and 0.06, with single allele effect sizes of -6 and -4, respectively. Native American genetic ancestry by itself reduced pigmentation by more than 20 melanin units (range 24-29). The responsible hypopigmenting genetic variants remain to be identified, since none of the published polymorphisms predicted in prior literature to affect skin color in Native Americans caused detectable hypopigmentation in the Kalinago.

Comparative Study of Combination of Oral Tranexamic Acid With Modified Kligman's Formula Versus Oral Tranexamic Acid With Azelaic Acid 15% in the Treatment of Melasma

Background Melasma is a persistent skin condition affecting individuals of Asian, African, and Hispanic backgrounds. It causes dark patches on sun-exposed areas of the face. The exact causes are unclear, but UV light and hormonal factors play a role. Melasma significantly impacts physical appearance and quality of life, causing emotional and social distress. Objective The objective was to compare the efficacy of a combination of oral tranexamic acid and modified Kligman's formula vs. oral tranexamic acid and 15% azelaic acid. Material and methods This two-year interventional study occurred at the Outpatient Department of Dermatology, Venereology, and Leprosy in Sawangi, Maharashtra. It included male and female patients aged 18-50 with melasma seeking treatment. Ethical approval was obtained, and data collection involved medical histories, skin examinations, and calculating the Melasma Area and Severity Index (MASI). Results The study found no significant association between age groups and subject distribution in Groups A and B. Significant differences were observed in MASI scores within each group over time. There was a significant difference in mean MASI scores between Group A and Group B at the eight-week mark. A burning sensation was significantly associated with the groups, while no significant association was found for erythema. Conclusion This study concludes that combining oral tranexamic acid with a modified Kligman's formula is more effective in treating melasma than combining oral tranexamic acid with azelaic acid 15%.

The evolution of human skin pigmentation: A changing medley of vitamins, genetic variability, and UV radiation during human expansion

This review examines putative, yet likely critical evolutionary pressures contributing to human skin pigmentation and subsequently, depigmentation phenotypes. To achieve this, it provides a synthesis of ideas that frame contemporary thinking, without limiting the narrative to pigmentation genes alone. It examines how geography and hence the quality and quantity of UV exposure, pigmentation genes, diet-related genes, vitamins, anti-oxidant nutrients, and cultural practices intersect and interact to facilitate the evolution of human skin color. The article has a strong focus on the vitamin D-folate evolutionary model, with updates on the latest biophysical research findings to support this paradigm. This model is examined within a broad canvas that takes human expansion out of Africa and genetic architecture into account. A thorough discourse on the biology of melanization is provided (includes relationship to BH4 and DNA damage repair), with the relevance of this to the UV sensitivity of folate and UV photosynthesis of vitamin D explained in detail, including the relevance of these vitamins to reproductive success. It explores whether we might be able to predict vitamin-related gene polymorphisms that pivot metabolism to the prevailing UVR exposome within the vitamin D-folate evolutionary hypothesis context. This is discussed in terms of a primary adaptive phenotype (pigmentation/depigmentation), a secondary adaptive phenotype (flexible metabolic phenotype based on vitamin-related gene polymorphism profile), and a tertiary adaptive strategy (dietary anti-oxidants to support the secondary adaptive phenotype). Finally, alternative evolutionary models for pigmentation are discussed, as are challenges to future research in this area.

The evolution of human skin pigmentation involved the interactions of genetic, environmental, and cultural variables

The primary biological role of human skin pigmentation is as a mediator of penetration of ultraviolet radiation (UVR) into the deep layers of skin and the cutaneous circulation. Since the origin of Homo sapiens, dark, protective constitutive pigmentation and strong tanning abilities have been favored under conditions of high UVR and represent the baseline condition for modern humans. The evolution of partly depigmented skin and variable tanning abilities has occurred multiple times in prehistory, as populations have dispersed into environments with lower and more seasonal UVR regimes, with unique complements of genes and cultural practices. The evolution of extremes of dark pigmentation and depigmentation has been rare and occurred only under conditions of extremely high or low environmental UVR, promoted by positive selection on variant pigmentation genes followed by limited gene flow. Over time, the evolution of human skin pigmentation has been influenced by the nature and course of human dispersals and modifications of cultural practices, which have modified the nature and actions of skin pigmentation genes. Throughout most of prehistory and history, the evolution of human skin pigmentation has been a contingent and non-deterministic process.

Prediction of eye and hair pigmentation phenotypes using the HIrisPlex system in a Brazilian admixed population sample

Human pigmentation is a complex trait, probably involving more than 100 genes. Predicting phenotypes using SNPs present in those genes is important for forensic purpose. For this, the HIrisPlex tool was developed for eye and hair color prediction, with both models achieving high accuracy among Europeans. Its evaluation in admixed populations is important, since they present a higher frequency of intermediate phenotypes, and HIrisPlex has demonstrated limitations in such predictions; therefore, the performance of this tool may be impaired in such populations. Here, we evaluate the set of 24 markers from the HIrisPlex system in 328 individuals from Ribeirão Preto (SP) region, predicting eye and hair color and comparing the predictions with their real phenotypes. We used the HaloPlex Target Enrichment System and MiSeq Personal Sequencer platform for massively parallel sequencing. The prediction of eye and hair color was accomplished by the HIrisPlex online tool, using the default prediction settings. Ancestry was estimated using the SNPforID 34-plex to observe if and how an individual's ancestry background would affect predictions in this admixed sample. Our sample presented major European ancestry (70.5%), followed by African (21.1%) and Native American/East Asian (8.4%). HIrisPlex presented an overall sensitivity of 0.691 for hair color prediction, with sensitivities ranging from 0.547 to 0.782. The lowest sensitivity was observed for individuals with black hair, who present a reduced European contribution (48.4%). For eye color prediction, the overall sensitivity was 0.741, with sensitivities higher than 0.85 for blue and brown eyes, although it failed in predicting intermediate eye color. Such struggle in predicting this phenotype category is in accordance with what has been seen in previous studies involving HIrisPlex. Individuals with brown eye color are more admixed, with European ancestry decreasing to 62.6%; notwithstanding that, sensitivity for brown eyes was almost 100%. Overall sensitivity increases to 0.791 when a 0.7 threshold is set, though 12.5% of the individuals become undefined. When combining eye and hair prediction, hit rates between 51.3 and 68.9% were achieved. Despite the difficulties with intermediate phenotypes, we have shown that HIrisPlex results can be very helpful when interpreted with caution.

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.

Adaptation and co-adaptation of skin pigmentation and vitamin D genes in native Americans

We carried out an exhaustive review regarding human skin color variation and how much it may be related to vitamin D metabolism and other photosensitive molecules. We discuss evolutionary contexts that modulate this variability and hypotheses postulated to explain them; for example, a small amount of melanin in the skin facilitates vitamin D production, making it advantageous to have fair skin in an environment with little radiation incidence. In contrast, more melanin protects folate from degradation in an environment with a high incidence of radiation. Some Native American populations have a skin color at odds with what would be expected for the amount of radiation in the environment in which they live, a finding challenging the so-called "vitamin D-folate hypothesis." Since food is also a source of vitamin D, dietary habits should also be considered. Here we argue that a gene network approach provides tools to explain this phenomenon since it indicates potential alleles co-evolving in a compensatory way. We identified alleles of the vitamin D metabolism and pigmentation pathways segregated together, but in different proportions, in agriculturalists and hunter-gatherers. Finally, we highlight how an evolutionary approach can be useful to understand current topics of medical interest.

Identifying adaptive alleles in the human genome: from selection mapping to functional validation

The suite of phenotypic diversity across geographically distributed human populations is the outcome of genetic drift, gene flow, and natural selection throughout human evolution. Human genetic variation underlying local biological adaptations to selective pressures is incompletely characterized. With the emergence of population genetics modeling of large-scale genomic data derived from diverse populations, scientists are able to map signatures of natural selection in the genome in a process known as selection mapping. Inferred selection signals further can be used to identify candidate functional alleles that underlie putative adaptive phenotypes. Phenotypic association, fine mapping, and functional experiments facilitate the identification of candidate adaptive alleles. Functional investigation of candidate adaptive variation using novel techniques in molecular biology is slowly beginning to unravel how selection signals translate to changes in biology that underlie the phenotypic spectrum of our species. In addition to informing evolutionary hypotheses of adaptation, the discovery and functional annotation of adaptive alleles also may be of clinical significance. While selection mapping efforts in non-European populations are growing, there remains a stark under-representation of diverse human populations in current public genomic databases, of both clinical and non-clinical cohorts. This lack of inclusion limits the study of human biological variation. Identifying and functionally validating candidate adaptive alleles in more global populations is necessary for understanding basic human biology and human disease.

Skin pigmentation and genetic variants in an admixed Brazilian population of primarily European ancestry

Although many genes have been shown to be associated with human pigmentary traits and forensic prediction assays exist (e.g. HIrisPlex-S), the genetic knowledge about skin colour remains incomplete. The highly admixed Brazilian population is an interesting study population for investigation of the complex genotype-phenotype architecture of human skin colour because of its large variation. Here, we compared variants in 22 pigmentary genes with quantitative skin pigmentation levels on the buttock, arm, and forehead areas of 266 genetically admixed Brazilian individuals. The genetic ancestry of each individual was estimated by typing 46 AIM-InDels. The mean proportion of genetic ancestry was 68.8% European, 20.8% Sub-Saharan African, and 10.4% Native American. A high correlation (adjusted R² = 0.65, p < 0.05) was observed between nine SNPs and quantitative skin pigmentation using multiple linear regression analysis. The correlations were notably smaller between skin pigmentation and biogeographic ancestry (adjusted R² = 0.45, p < 0.05), or markers in the leading forensic skin colour prediction system, the HIrisPlex-S (adjusted R² = 0.54, p < 0.05). Four of the nine SNPs, OCA2 rs1448484 (rank 2), APBA2 rs4424881 (rank 4), MFSD12 rs10424065 (rank 8), and TYRP1 1408799 (rank 9) were not investigated as part of the HIrisPlex-S selection process, and therefore not included in the HIrisPlex-S model. Our results indicate that these SNPs account for a substantial part of the skin colour variation in individuals of admixed ancestry. Hence, we suggest that these SNPs are considered when developing future skin colour prediction models.

MC1R Variation in a New Mexico Population

BACKGROUND

The Melanocortin 1 Receptor (MC1R) contributes to pigmentation, an important risk factor for developing melanoma. Evaluating SNPs in MC1R and association with race/ethnicity, skin type, and perceived cancer risk in a New Mexico (NM) population will elucidate the role of MC1R in a multicultural population.

METHODS

We genotyped MC1R in 191 NMs attending a primary care clinic in Albuquerque. We obtained individuals' self-identified race/ethnicity, skin type, and perceived cancer risk. We defined genetic risk as carriage of any one or more of the nine most common SNPs in MC1R.

RESULTS

We found that one MC1R SNP, R163Q (rs885479), was identified in 47.6% of self-identified Hispanics and 12.9% of non-Hispanic whites (NHW), making Hispanics at higher "genetic risk" (as defined by carrying one of the MC1R common variants). When we deleted R163Q from analyses, Hispanics were no longer at higher genetic risk (33.3%) compared with NHW (48.3%), consistent with melanoma rates, tanning ability, and lower perceived risk. Hispanics had a perceived risk significantly lower than NHW and a nonsignificant better tanning ability than NHW.

CONCLUSIONS

The R163Q variant in MC1R may not be a risk factor for melanoma among NM Hispanics. This suggestion points to the need to carefully interpret genetic risk factors among specific populations.

IMPACT

Genetic risk cannot be extrapolated from Northern European populations directly to non-European populations.

The Evolutionary History of Human Skin Pigmentation

Skin pigmentation is a complex, conspicuous, highly variable human trait that exhibits a remarkable correlation with latitude. The evolutionary history and genetic basis of skin color variation has been the subject of intense research in the last years. This article reviews the major hypotheses explaining skin color diversity and explores the implications of recent findings about the genes associated with skin pigmentation for understanding the evolutionary forces that have shaped the current patterns of skin color variation. A major aspect of these findings is that the genetic basis of skin color is less simple than previously thought and that geographic variation in skin pigmentation was influenced by the concerted action of different types of natural selection, rather than just by selective sweeps in a few key genes.

Rapid evolution of a skin-lightening allele in southern African KhoeSan

Skin pigmentation is under strong directional selection in northern European and Asian populations. The indigenous KhoeSan populations of far southern Africa have lighter skin than other sub-Saharan African populations, potentially reflecting local adaptation to a region of Africa with reduced UV radiation. Here, we demonstrate that a canonical Eurasian skin pigmentation gene, SLC24A5, was introduced to southern Africa via recent migration and experienced strong adaptive evolution in the KhoeSan. To reconstruct the evolution of skin pigmentation, we collected phenotypes from over 400 ≠Khomani San and Nama individuals and high-throughput sequenced candidate pigmentation genes. The derived causal allele in SLC24A5, p.Ala111Thr, significantly lightens basal skin pigmentation in the KhoeSan and explains 8 to 15% of phenotypic variance in these populations. The frequency of this allele (33 to 53%) is far greater than expected from colonial period European gene flow; however, the most common derived haplotype is identical among European, eastern African, and KhoeSan individuals. Using four-population demographic simulations with selection, we show that the allele was introduced into the KhoeSan only 2,000 y ago via a back-to-Africa migration and then experienced a selective sweep (s = 0.04 to 0.05 in ≠Khomani and Nama). The SLC24A5 locus is both a rare example of intense, ongoing adaptation in very recent human history, as well as an adaptive gene flow at a pigmentation locus in humans.

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.

Influence of maternal and own genotype at tanning dependence-related SNPs on sun exposure in childhood

Background

Research suggests there may be a genetic influence on the likelihood of becoming tanning dependent (TD). The way in which mothers regulate their children’s sun exposure may be affected by being TD. We investigated the associations between single nucleotide polymorphisms (SNPs) related to being TD and early sun exposure.

Methods

Data from the Avon Longitudinal Study of Parents and Children (ALSPAC) were used. Associations between 17 TD related SNPs in children and their mothers and 10 sun exposure variables in children (assessed via questionnaire at age 8) were analyzed in logistic and ordinal logistic regressions. Analyses were adjusted for principal components of population structure and age (at time of questionnaire response). Models with additional adjustment for maternal or offspring genotypes were also tested. Secondary analyses included adjustment for sex and skin pigmentation.

Results

Among ALSPAC children, the rs29132 SNP in the Vesicle-associated membrane protein-associated protein A (VAPA) gene was associated with five sun exposure variables whilst the rs650662 SNP in the Opioid Receptor Mu 1 (OPRM1) gene was associated with three. The remaining SNPs did not show associations beyond what was expected by chance. After Bonferroni correction one SNP in the children was associated with an increased likelihood of using sun cream whilst in the sun at 8 years old (rs60050811 in the Spermatogenesis and Centriole Associated 1 (SPATC1) gene, OR per C allele = 1.34, 95% CI 1.11–1.62, p = .003). In the mothers, rs650662 in OPRM1 was associated with the use of a lower factor of sun cream in their children, (OR per A allele = 0.89, 95% CI 0.82–0.96, p = .002). Whilst rs2073478 in the Aldehyde Dehydrogenase 1 Family Member B1 (ALDH1B1) gene was associated with a reduced odds of their child using a sun block or cream with a 4 star rating (OR per T allele = 0.68, 95% CI 0.53–0.88, p = .003). Similar but weaker associations were observed for the main findings in the secondary analyses.

Conclusions

We found weak evidence to suggest that genes previously associated with TD are associated with sun exposure in children of European ancestry from southwest England.

Electronic supplementary material

The online version of this article (10.1186/s12881-018-0575-z) contains supplementary material, which is available to authorized users.

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'.

Global skin colour prediction from DNA

Human skin colour is highly heritable and externally visible with relevance in medical, forensic, and anthropological genetics. Although eye and hair colour can already be predicted with high accuracies from small sets of carefully selected DNA markers, knowledge about the genetic predictability of skin colour is limited. Here, we investigate the skin colour predictive value of 77 single-nucleotide polymorphisms (SNPs) from 37 genetic loci previously associated with human pigmentation using 2025 individuals from 31 global populations. We identified a minimal set of 36 highly informative skin colour predictive SNPs and developed a statistical prediction model capable of skin colour prediction on a global scale. Average cross-validated prediction accuracies expressed as area under the receiver-operating characteristic curve (AUC) ± standard deviation were 0.97 ± 0.02 for Light, 0.83 ± 0.11 for Dark, and 0.96 ± 0.03 for Dark-Black. When using a 5-category, this resulted in 0.74 ± 0.05 for Very Pale, 0.72 ± 0.03 for Pale, 0.73 ± 0.03 for Intermediate, 0.87±0.1 for Dark, and 0.97 ± 0.03 for Dark-Black. A comparative analysis in 194 independent samples from 17 populations demonstrated that our model outperformed a previously proposed 10-SNP-classifier approach with AUCs rising from 0.79 to 0.82 for White, comparable at the intermediate level of 0.63 and 0.62, respectively, and a large increase from 0.64 to 0.92 for Black. Overall, this study demonstrates that the chosen DNA markers and prediction model, particularly the 5-category level; allow skin colour predictions within and between continental regions for the first time, which will serve as a valuable resource for future applications in forensic and anthropologic genetics.

African ancestry is associated with facial melasma in women: a cross-sectional study

Background

Melasma is a chronic acquired focal hypermelanosis affecting photoexposed areas, especially for women during fertile age. Several factors contribute to its development: sun exposure, sex steroids, medicines, and family history. Melanic pigmentation pathway discloses several SNPs in different populations. Here, we evaluated the association between genetic ancestry and facial melasma.

Methods

A cross-sectional study involving women with melasma and an age-matched control group from outpatients at FMB-Unesp, Botucatu-SP, Brazil was performed. DNA was extracted from oral mucosa swabs and ancestry determined by studying 61 INDELs. The genetic ancestry components were adjusted by other known risk factors by multiple logistic regression.

Results

We evaluated 119 women with facial melasma and 119 controls. Mean age was 39 ± 9 years. Mean age at beginning of disease was 27 ± 8 years. Pregnancy (40%), sun exposure (37%), and hormonal oral contraception (22%) were the most frequently reported melasma triggers. All subjects presented admixed ancestry, African and European genetic contributions were significantly different between cases and controls (respectively 10% vs 6%; 77% vs 82%; p < 0.05). African ancestry (OR = 1.04; 95% CI 1.01 to 1.07), first generation family history (OR = 3.04; 95% CI 1.56 to 5.94), low education level (OR = 4.04; 95% CI 1.56 to 5.94), and use of antidepressants by individuals with affected family members (OR = 6.15; 95% CI 1.13 to 33.37) were associated with melasma, independently of other known risk factors.

Conclusions

Facial melasma was independently associated with African ancestry in a highly admixed population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-017-0378-7) contains supplementary material, which is available to authorized users.

Hyperpigmentation after surgery for a deep dermal burn of the dorsum of the hand: partial-thickness debridement followed by medium split-thickness skin grafting vs full-thickness debridement followed by thick split-thickness skin grafting

BACKGROUND

Early excision and skin grafting are commonly used to treat deep dermal burns (DDBs) of the dorsum of the hand. Partial-thickness debridement (PTD) is one of the most commonly used procedures for the excision of burned tissue of the dorsum of the hand. In contrast, full-thickness debridement (FTD) has also been reported. However, it is unclear whether PTD or FTD is better.

METHODS

In this hospital-based retrospective study, we compared the outcomes of PTD followed by a medium split-thickness skin graft (STSG) with FTD followed by a thick STSG to treat a DDB of the dorsum of the hand in Japanese patients. To evaluate postoperative pigmentation of the skin graft, quantitative analyses were performed using the red, green, and blue (RGB) and the hue, saturation, and brightness (HSB) color spaces. We have organized the manuscript in a manner compliant with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement.

RESULTS

Data from 11 patients were analyzed. Six hands (five patients) received grafts in the PTD group and eight hands (six patients) received grafts in the FTD group. Graft take was significantly better in the FTD group (median 98 %, interquartile range 95-99) than in the PTD group (median 90 %, interquartile range 85-90) (P < 0.01). Quantitative skin color analyses in both the RGB and HSB color spaces showed that postoperative grafted skin was significantly darker than the adjacent control area in the PTD group, but not in the FTD group.

CONCLUSIONS

There is a possibility that FTD followed by a thick STSG is an option that can reduce the risk of hyperpigmentation after surgery for DDB of the dorsum of the hand in Japanese patients. Further investigation is needed to clarify whether the FTD or the thick STSG or both are the factor for the control of hyperpigmentation.

Bioinformatics and evolution of vertebrate nociceptin and opioid receptors

G protein-coupled receptors (GPCRs) are ancestrally related membrane proteins on cells that mediate the pharmacological effect of most drugs and neurotransmitters. GPCRs are the largest group of membrane receptor proteins encoded in the human genome. One of the most famous types of GPCRs is the opioid receptors. Opioid family receptors consist of four closely related proteins expressed in all vertebrate brains and spinal cords examined to date. The three classical types of opioid receptors shown unequivocally to mediate analgesia in animal models and in humans are the mu- (MOR), delta- (DOR), and kappa-(KOR) opioid receptor proteins. The fourth and most recent member of the opioid receptor family discovered is the nociceptin or orphanin FQ receptor (ORL). The role of ORL and its ligands in producing analgesia is not as clear, with both analgesic and hyperalgesic effects reported. All four opioid family receptor genes were cloned from expressed mRNA in a number of vertebrate species, and there are enough sequences presently available to carry out bioinformatic analysis. This chapter presents the results of a comparative analysis of vertebrate opioid receptors using pharmacological studies, bioinformatics, and the latest data from human whole-genome studies. Results confirm our initial hypotheses that the four opioid receptor genes most likely arose by whole-genome duplication, that there is an evolutionary vector of opioid receptor type divergence in sequence and function, and that the hMOR gene shows evidence of positive selection or adaptive evolution in Homo sapiens.

Implications of the admixture process in skin color molecular assessment

The understanding of the complex genotype-phenotype architecture of human pigmentation has clear implications for the evolutionary history of humans, as well as for medical and forensic practices. Although dozens of genes have previously been associated with human skin color, knowledge about this trait remains incomplete. In particular, studies focusing on populations outside the European-North American axis are rare, and, until now, admixed populations have seldom been considered. The present study was designed to help fill this gap. Our objective was to evaluate possible associations of 18 single nucleotide polymorphisms (SNPs), located within nine genes, and one pseudogene with the Melanin Index (MI) in two admixed Brazilian populations (Gaucho, N = 352; Baiano, N = 148) with different histories of geographic and ethnic colonization. Of the total sample, four markers were found to be significantly associated with skin color, but only two (SLC24A5 rs1426654, and SLC45A2 rs16891982) were consistently associated with MI in both samples (Gaucho and Baiano). Therefore, only these 2 SNPs should be preliminarily considered to have forensic significance because they consistently showed the association independently of the admixture level of the populations studied. We do not discard that the other two markers (HERC2 rs1129038 and TYR rs1126809) might be also relevant to admixed samples, but additional studies are necessary to confirm the real importance of these markers for skin pigmentation. Finally, our study shows associations of some SNPs with MI in a modern Brazilian admixed sample, with possible applications in forensic genetics. Some classical genetic markers in Euro-North American populations are not associated with MI in our sample. Our results point out the relevance of considering population differences in selecting an appropriate set of SNPs as phenotype predictors in forensic practice.

The evolution of skin pigmentation and hair texture in people of African ancestry

Our species, Homo sapiens, evolved in Africa, and humanity's highest levels of genetic diversity are maintained there today. Underlying genetic diversity combined with the great range of solar regimes and climatic conditions found in Africa has contributed to a wide range of human integumentary phenotypes within the continent. Millions of Africans have moved, voluntarily and involuntarily, to other continents in the past 2000 years, and the range of integumentary phenotypes among admixed African diaspora populations is enormous. In this contribution, we do not catalog this variation, but provide basic evolutionary background as to how it developed in the first place.

Lactase persistence and augmented salivary alpha-amylase gene copy numbers might have been selected by the combined toxic effects of gluten and (food born) pathogens

Various positively selected adaptations to new nutrients have been identified. Lactase persistence is among the best known, conferring the ability for drinking milk at post weaning age. An augmented number of amylase gene (AMY1) copies, giving rise to higher salivary amylase activity, has been implicated in the consumption of starch-rich foods. Higher AMY1 copy numbers have been demonstrated in populations with recent histories of starchy-rich diets. It is however questionable whether the resulting polymorphisms have exerted positive selection only by providing easily available sources of macro and micronutrients. Humans have explored new environments more than any other animal. Novel environments challenge the host, but especially its immune system with new climatic conditions, food and especially pathogens. With the advent of the agricultural revolution and the concurrent domestication of cattle came new pathogens. We contend that specific new food ingredients (e.g., gluten) and novel pathogens drove selection for lactase persistence and higher AMY gene copy numbers. Both adaptations provide ample glucose for activating the sodium glucose-dependent co-transporter 1 (SGLT1), which is the principal glucose, sodium and water transporter in the gastro-intestinal tract. Their rapid uptake confers protection against potentially lethal dehydration, hyponatremia and ultimately multiple organ failure. Oral rehydration therapy aims at SGLT1 activity and is the current treatment of choice for chronic diarrhoea and vomiting. We hypothesize that lifelong lactase activity and rapid starch digestion should be looked at as the evolutionary covalent of oral rehydration therapy.

Endogenous opiates and behavior: 2012

This paper is the thirty-fifth consecutive installment of the annual review of research concerning the endogenous opioid system. It summarizes papers published during 2012 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides, opioid receptors, opioid agonists and opioid antagonists. The particular topics that continue to be covered include the molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors related to behavior (Section 2), and the roles of these opioid peptides and receptors in pain and analgesia (Section 3); stress and social status (Section 4); tolerance and dependence (Section 5); learning and memory (Section 6); eating and drinking (Section 7); alcohol and drugs of abuse (Section 8); sexual activity and hormones, pregnancy, development and endocrinology (Section 9); mental illness and mood (Section 10); seizures and neurologic disorders (Section 11); electrical-related activity and neurophysiology (Section 12); general activity and locomotion (Section 13); gastrointestinal, renal and hepatic functions (Section 14); cardiovascular responses (Section 15); respiration and thermoregulation (Section 16); and immunological responses (Section 17).

The light skin allele of SLC24A5 in South Asians and Europeans shares identity by descent

Skin pigmentation is one of the most variable phenotypic traits in humans. A non-synonymous substitution (rs1426654) in the third exon of SLC24A5 accounts for lighter skin in Europeans but not in East Asians. A previous genome-wide association study carried out in a heterogeneous sample of UK immigrants of South Asian descent suggested that this gene also contributes significantly to skin pigmentation variation among South Asians. In the present study, we have quantitatively assessed skin pigmentation for a largely homogeneous cohort of 1228 individuals from the Southern region of the Indian subcontinent. Our data confirm significant association of rs1426654 SNP with skin pigmentation, explaining about 27% of total phenotypic variation in the cohort studied. Our extensive survey of the polymorphism in 1573 individuals from 54 ethnic populations across the Indian subcontinent reveals wide presence of the derived-A allele, although the frequencies vary substantially among populations. We also show that the geospatial pattern of this allele is complex, but most importantly, reflects strong influence of language, geography and demographic history of the populations. Sequencing 11.74 kb of SLC24A5 in 95 individuals worldwide reveals that the rs1426654-A alleles in South Asian and West Eurasian populations are monophyletic and occur on the background of a common haplotype that is characterized by low genetic diversity. We date the coalescence of the light skin associated allele at 22–28 KYA. Both our sequence and genome-wide genotype data confirm that this gene has been a target for positive selection among Europeans. However, the latter also shows additional evidence of selection in populations of the Middle East, Central Asia, Pakistan and North India but not in South India.

Human skin color is one of the most visible aspects of human diversity. The genetic basis of pigmentation in Europeans has been understood to some extent, but our knowledge about South Asians has been restricted to a handful of studies. It has been suggested that a single nucleotide difference in SLC24A5 accounts for 25–38% European-African pigmentation differences and correlates with lighter skin. This genetic variant has also been associated with skin color variation among South Asians living in the UK. Here, we report a study based on a homogenous cohort of South India. Our results confirm that SLC24A5 plays a key role in pigmentation diversity of South Asians. Country-wide screening of the variant reveals that the light skin associated allele is widespread in the Indian subcontinent and its complex patterning is shaped by a combination of processes involving selection and demographic history of the populations. By studying the variation of SLC24A5 sequences among a diverse set of individuals, we show that the light skin associated allele in South Asians is identical by descent to that found in Europeans. Our study also provides new insights into positive selection acting on the gene and the evolutionary history of light skin in humans.

Exploring signatures of positive selection in pigmentation candidate genes in populations of East Asian ancestry

Background

Currently, there is very limited knowledge about the genes involved in normal pigmentation variation in East Asian populations. We carried out a genome-wide scan of signatures of positive selection using the 1000 Genomes Phase I dataset, in order to identify pigmentation genes showing putative signatures of selective sweeps in East Asia. We applied a broad range of methods to detect signatures of selection including: 1) Tests designed 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). Based on the results obtained from a genome wide analysis of 25 kb windows, we constructed an empirical distribution for each statistic across all windows, and identified pigmentation genes that are outliers in the distribution.

Results

Our tests identified twenty genes that are relevant for pigmentation biology. Of these, eight genes (ATRN, EDAR, KLHL7, MITF, OCA2, TH, TMEM33 and TRPM1,) were extreme outliers (top 0.1% of the empirical distribution) for at least one statistic, and twelve genes (ADAM17, BNC2, CTSD, DCT, EGFR, LYST, MC1R, MLPH, OPRM1, PDIA6, PMEL (SILV) and TYRP1) were in the top 1% of the empirical distribution for at least one statistic. Additionally, eight of these genes (BNC2, EGFR, LYST, MC1R, OCA2, OPRM1, PMEL (SILV) and TYRP1) have been associated with pigmentary traits in association studies.

Conclusions

We identified a number of putative pigmentation genes showing extremely unusual patterns of genetic variation in East Asia. Most of these genes are outliers for different tests and/or different populations, and have already been described in previous scans for positive selection, providing strong support to the hypothesis that recent selective sweeps left a signature in these regions. However, it will be necessary to carry out association and functional studies to demonstrate the implication of these genes in normal pigmentation variation.

The probability of genetic parallelism and convergence in natural populations

Genomic and genetic methods allow investigation of how frequently the same genes are used by different populations during adaptive evolution, yielding insights into the predictability of evolution at the genetic level. We estimated the probability of gene reuse in parallel and convergent phenotypic evolution in nature using data from published studies. The estimates are surprisingly high, with mean probabilities of 0.32 for genetic mapping studies and 0.55 for candidate gene studies. The probability declines with increasing age of the common ancestor of compared taxa, from about 0.8 for young nodes to 0.1–0.4 for the oldest nodes in our study. Probability of gene reuse is higher when populations begin from the same ancestor (genetic parallelism) than when they begin from divergent ancestors (genetic convergence). Our estimates are broadly consistent with genomic estimates of gene reuse during repeated adaptation to similar environments, but most genomic studies lack data on phenotypic traits affected. Frequent reuse of the same genes during repeated phenotypic evolution suggests that strong biases and constraints affect adaptive evolution, resulting in changes at a relatively small subset of available genes. Declines in the probability of gene reuse with increasing age suggest that these biases diverge with time.

Comprehensive candidate gene study highlights UGT1A and BNC2 as new genes determining continuous skin color variation in Europeans

Natural variation in human skin pigmentation is primarily due to genetic causes rooted in recent evolutionary history. Genetic variants associated with human skin pigmentation confer risk of skin cancer and may provide useful information in forensic investigations. Almost all previous gene-mapping studies of human skin pigmentation were based on categorical skin color information known to oversimplify the continuous nature of human skin coloration. We digitally quantified skin color into hue and saturation dimensions for 5,860 Dutch Europeans based on high-resolution skin photographs. We then tested an extensive list of 14,185 single nucleotide polymorphisms in 281 candidate genes potentially involved in human skin pigmentation for association with quantitative skin color phenotypes. Confirmatory association was revealed for several known skin color genes including HERC2, MC1R, IRF4, TYR, OCA2, and ASIP. We identified two new skin color genes: genetic variants in UGT1A were significantly associated with hue and variants in BNC2 were significantly associated with saturation. Overall, digital quantification of human skin color allowed detecting new skin color genes. The variants identified in this study may also contribute to the risk of skin cancer. Our findings are also important for predicting skin color in forensic investigations.

Human pigmentation genes under environmental selection

Genome-wide association studies and comparative genomics have established major loci and specific polymorphisms affecting human skin, hair and eye color. Environmental changes have had an impact on selected pigmentation genes as populations have expanded into different regions of the globe.

A high-coverage genome sequence from an archaic Denisovan individual

We present a DNA library preparation method that has allowed us to reconstruct a high-coverage (30×) genome sequence of a Denisovan, an extinct relative of Neandertals. The quality of this genome allows a direct estimation of Denisovan heterozygosity indicating that genetic diversity in these archaic hominins was extremely low. It also allows tentative dating of the specimen on the basis of "missing evolution" in its genome, detailed measurements of Denisovan and Neandertal admixture into present-day human populations, and the generation of a near-complete catalog of genetic changes that swept to high frequency in modern humans since their divergence from Denisovans.