Genetics of skin color variation in Europeans: genome-wide association studies with functional follow-up

Prediction of Skin Color Using Forensic DNA Phenotyping in Asian Populations: A Focus on Thailand

Forensic DNA phenotyping (FDP) has emerged as an essential tool in criminal investigations, enabling the prediction of physical traits based on genetic information. This review explores the genetic factors influencing skin pigmentation, particularly within Asian populations, with a focus on Thailand. Key genes such as Oculocutaneous Albinism II (OCA2), Dopachrome Tautomerase (DCT), KIT Ligand (KITLG), and Solute Carrier Family 24 Member 2 (SLC24A2) are examined for their roles in melanin production and variations that lead to different skin tones. The OCA2 gene is highlighted for its role in transporting ions that help stabilize melanosomes, while specific variants in the DCT gene, including single nucleotide polymorphisms (SNPs) rs2031526 and rs3782974, are discussed for their potential effects on pigmentation in Asian groups. The KITLG gene, crucial for developing melanocytes, includes the SNP rs642742, which is linked to lighter skin in East Asians. Additionally, recent findings on the SLC24A2 gene are presented, emphasizing its connection to pigmentation through calcium regulation in melanin production. Finally, the review addresses the ethical considerations of using FDP in Thailand, where advances in genetic profiling raise concerns about privacy, consent, and discrimination. Establishing clear guidelines is vital to balancing the benefits of forensic DNA applications with the protection of individual rights.

Parents facing polygenic embryo scores: the 'best choice of a best life' and psychological counselling

Advances in reproductive medicine and genetic technologies now offer prospective parents the option to test IVF embryos for genetic predispositions to complex diseases, such as coronary heart disease and psychiatric disorders, through polygenic embryo screening (PES). However, limited clinical data on its real-world use leaves parents facing complex decisions based on probabilistic risk scores, requiring them to weigh uncertain benefits against potential harms. While clinical, ethical, and societal concerns regarding PES have been extensively discussed, the psychological considerations have received less attention. This paper highlights the importance of decision aids as part of psychological interventions, which are crucial for helping parents navigate these choices and make informed decisions based on individual perceptions and experiences. Additionally, determining how and when to disclose genetic risk information to children presents significant challenges for families. Early disclosure may lead to anxiety, while withholding information could undermine trust later in life. Psychological counseling is therefore an essential component in supporting families through these sensitive decisions. While PES offers opportunities to reduce genetic risks, it also introduces significant challenges that require thoughtful consideration and comprehensive support for both parents and children.

Comprehensive Profiling of Genetic and Nongenetic Factors that Influence Skin Traits in Asian Women from 4 Countries

Brightness, pore size, wrinkles, color, elasticity, moisture, transepidermal water loss, pH, and sebum are traits that directly affect skin health and aesthetics, which vary substantially among ethnic groups and individuals. Although a complex combination of genetic and environmental factors influences these traits, an accurate estimation of skin phenotypic determination in individuals from Asian countries remains to be determined. In this study, we established the ASTAGR (Asian Skin Traits and Genetics Research) cohort, comprising 2762 female participants from 4 Asian countries (Korea, Vietnam, Thailand, and Indonesia). We documented 26 skin traits and conducted extensive skin-related surveys among the participants. Estimation of the quantitative contribution from whole-genome sequencing-based genotyping, combined with the survey findings, led to a comprehensive documentation of skin trait determination and its variation by country. This approach also revealed major variables affecting each skin trait. Our transethnic GWAS revealed 10 significant associations, 5 of which, to our knowledge, have not been previously reported. This study provides a comprehensive profile of skin traits and their regulation in the 4 Asian countries.

Population Genomics Reveals Elevated Inbreeding and Accumulation of Deleterious Mutations in White Raccoon Dogs

The formation of animal breeds usually begins with a small subsample from their ancestral population. Deleterious mutations accumulate in the population under genetic drift, inbreeding, and artificial selection during the development and maintenance of traits desired by humans. White raccoon dogs are among the most popular breeds of farmed raccoon dogs, but white raccoon dogs are more susceptible to disease and have a lower reproductive ability. However, the accumulation of deleterious mutations in this white breed is largely unknown. By analyzing and comparing whole-genome sequencing data from 20 white raccoon dogs and 38 normal raccoon dogs, we detected an increased occurrence of loss-of-function (LoF) mutations in white raccoon dogs compared with normal raccoon dogs. With the finding of a significantly higher dosage of homozygous missense mutations in the white raccoon dog genome, we detected a greater fitness cost in white raccoon dogs. Although a much higher FROH level for ROH fragments longer than 1 Mb has been reported in white raccoon dogs, we did not detect a genetic signal of genetic purging in white raccoon dogs. This study provides valuable genomic resources and new insights into the accumulation of mutation loads in farmed raccoon dogs.

Sequential introgression of a carotenoid processing gene underlies sexual ornament diversity in a genus of manakins

In a hybrid zone between two tropical lekking birds, yellow male plumage of one species has introgressed asymmetrically replacing white plumage of another via sexual selection. Here, we present a detailed analysis of the plumage trait to uncover its physical and genetic bases and trace its evolutionary history. We determine that the carotenoid lutein underlies the yellow phenotype and describe microstructural feather features likely to enhance color appearance. These same features reduce predicted water shedding capacity of feathers, a potential liability in the tropics. Through genome-scale DNA sequencing of hybrids and each species in the genus, we identify BCO2 as the major gene responsible for the color polymorphism. The BCO2 gene tree and genome-wide allele frequency patterns suggest that carotenoid-pigmented collars initially arose in a third species and reached the hybrid zone through historical gene flow. Complex interplay between sexual selection and hybridization has thus shaped phenotypes of these species, where conspicuous sexual traits are key to male reproductive success.

Genome wide association study and meta-analysis identified multiple new risk loci for freckles in 4813 Chinese individuals

Freckle is a prevalent pigmentary dermatosis with an obvious hereditary component. Dozens of freckles risk loci have been discovered through research on multiple traits or other diseases, rather than as an independent trait. To discover novel variants associated with freckles, we performed GWAS and meta-analysis in 4813 Chinese individuals. We conducted GWAS and meta-analysis of two cohorts: 197 patients and 1603 controls (Cohort I), and 336 patients and 2677 controls (Cohort II), both from China. Then we performed linkage disequilibrium (LD) analysis, eQTL study, and enrichment analysis with association results for functional implications. Finally, we discovered 59 new SNPs and 13 novel susceptibility genes associated with freckles (Pmeta <5 × 10-8), which has enriched the genetic research on freckles.

Uncovering structural variants in Creole cattle from Guadeloupe and their impact on environmental adaptation through whole genome sequencing

Structural variants play an important role in evolutionary processes. Besides, they constitute a large source of inter individual genetic variation that might represent a major factor in the aetiology of complex, multifactorial traits. Their importance in adaptation is becoming increasingly evident in literature. Yet, the characterization of the genomic landscape of structural variants in local breeds remains scarce to date. Herein, we investigate patterns and gene annotation of structural variants in the Creole cattle from Guadeloupe breed using whole genome sequences from 23 bulls representative of the population. In total, we detected 32821 ascertained SV defining 15258 regions, representing ~ 17% of the Creole cattle genome. Among these, 6639 regions have not been previously reported in the Database of Genomic Variants archive. Average number of structural variants detected per individual in the studied population is in the same order of magnitude of that observed in indicine populations and higher than that reported in taurine breeds. We observe an important within-individual variability where approximately half of the detected structural variants have low frequency (MAF < 0.25). Most of the detected structural variants (55%) occurred in intergenic regions. Genic structural variants overlapped with 7793 genes and the predicted effect of most of them is ranked as "modifier". Among the structural variants that were predicted to have a high functional impact on the protein, a 5.5 Kb in length, highly frequent deletion on chromosome 2, affects ALPI, a gene associated with the interaction between gut microbiota and host immune system. The 6639 newly identified structural variants regions include three deletions and three duplications shared by more than 80% of individuals that are significantly enriched for genes related to tRNA threonylcarbamoyladenosine metabolic process, important for temperature adaptation in thermophilic organisms, therefore suggesting a potential role in the thermotolerance of Creole cattle from Guadeloupe cattle to tropical climate. Overall, highly frequent structural variants that are specific to the Creole cattle population encompass olfactory receptor and immunity genes as well as genes involved in muscle tone, muscle development and contraction. Beyond mapping and characterizing structural variants in the Creole cattle from Guadeloupe breed, this study provides valuable information for a better understanding of the potential role of chromosomal rearrangements in adaptive traits in cattle.

A sequence of SVA retrotransposon insertions in ASIP shaped human pigmentation

Retrotransposons comprise about 45% of the human genome1, but their contributions to human trait variation and evolution are only beginning to be explored2,3. Here, we find that a sequence of SVA retrotransposon insertions in an early intron of the ASIP (agouti signaling protein) gene has probably shaped human pigmentation several times. In the UK Biobank (n = 169,641), a recent 3.3-kb SVA insertion polymorphism associated strongly with lighter skin pigmentation (0.22 [0.21-0.23] s.d.; P = 2.8 × 10-351) and increased skin cancer risk (odds ratio = 1.23 [1.18-1.27]; P = 1.3 × 10-28), appearing to underlie one of the strongest common genetic influences on these phenotypes within European populations4-6. ASIP expression in skin displayed the same association pattern, with the SVA insertion allele exhibiting 2.2-fold (1.9-2.6) increased expression. This effect had an unusual apparent mechanism: an earlier, nonpolymorphic, human-specific SVA retrotransposon 3.9 kb upstream appeared to have caused ASIP hypofunction by nonproductive splicing, which the new (polymorphic) SVA insertion largely eliminated. Extended haplotype homozygosity indicated that the insertion allele has risen to allele frequencies up to 11% in European populations over the past several thousand years. These results indicate that a sequence of retrotransposon insertions contributed to a species-wide increase, then a local decrease, of human pigmentation.

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.

Comparisons between wrinkles and photo-ageing detected and self-reported by the participant or identified by trained assessors reveal insights from Chinese individuals in the Singapore/Malaysia Cross-sectional Genetics Epidemiology Study (SMCGES) cohort

BACKGROUND

Changes develop on the facial skin as a person ages. Other than chronological time, it has been discovered that gender, ethnicity, air pollution, smoking, nutrition, and sun exposure are notable risk factors that influence the development of skin ageing phenotypes such as wrinkles and photo-ageing. These risk factors can be quantified through epidemiological collection methods. We previously studied wrinkles and photo-ageing in detail using photo-numeric scales. The analysis was performed on the ethnic Chinese skin by three trained assessors. Recent studies have shown that it is possible to use self-reported data to identify skin-related changes including skin colour and skin cancer. In order to investigate the association between risk factors and skin ageing phenotypic outcomes in large-scale epidemiological studies, it would be useful to evaluate whether it is also possible for participants to self-report signs of ageing on their skin.

AIM

We have previously identified several validated photo-numeric scales for wrinkling and photo-ageing to use on ethnic Chinese skin. Using these scales, our trained assessors grade wrinkling and photo-ageing with moderately high inter-assessor concordance and agreement. The main objective of this study involves letting participants grade self-reported wrinkling and photo-ageing using these same scales. We aim to compare the concordance and agreement between signs of skin ageing by the participant and signs of ageing identified by our assessors.

METHOD

Three trained assessors studied facial photo-ageing on 1081 ethnic Chinese young adults from the Singapore/Malaysia Cross-sectional Genetics Epidemiology Study (SMCGES) cohort. Self-reported facial photo-ageing data by the same 1081 participants were also collated and the two sets of data are compared.

RESULTS

Here, we found that self-reported signs of photo-ageing are concordant with photo-ageing detected by our assessors. This finding is consistent whether photo-ageing is evaluated through studying wrinkle variations (Spearman's rank correlation (ρ) value: 0.246-0.329) or through studying dyspigmentation patterns (Spearman's rank correlation (ρ) value 0.203-0.278). When studying individual wrinkles, both participants and assessors often detect the presence of the same wrinkle (Spearman's rank correlation (ρ) value 0.249-0.366). A weak-to-fair level of agreement between both participants and assessors (Cohen's kappa (κ) values: 0.041-0.233) persists and is statistically significant after accounting for agreements due to chance. Both the participant and the assessor are largely consistent in evaluating the extent of photo-ageing (area under curve (AUC) values 0.689-0.769) and in discerning between the presence or absence of a given facial wrinkle (area under curve (AUC) values 0.601-0.856).

CONCLUSION

When we analyse the overall appearance of the face, our results show that signs of photo-ageing identified by the participant are concordant with signs of photo-ageing identified by our assessors. When we focused our analysis on specific areas of the face, we found that participants were more likely to identify and self-report the same wrinkles that our assessors have also detected. Here, we found that self-reported signs of skin ageing provide a satisfactory approximation to the signs of skin ageing identified by our assessors. The ability to use self-reported signs of skin ageing should also be evaluated on scales beyond the ones discussed in this study. Currently, there are not as many photo-numeric scales for quantifying dyspigmentation patterns as there are for quantifying wrinkle variations. As Chinese skin is known to become dyspigmented more easily with age, more photo-numeric scales need to be developed and properly validated.

Public Attitudes, Interests, and Concerns Regarding Polygenic Embryo Screening

Importance

Polygenic embryo screening (PES) is a novel technology that estimates the likelihood of developing future conditions (eg, diabetes or depression) and traits (eg, height or cognitive ability) in human embryos, with the goal of selecting which embryos to use. Given its commercial availability and concerns raised by researchers, clinicians, bioethicists, and professional organizations, it is essential to inform key stakeholders and relevant policymakers about the public's perspectives on this technology.

Objective

To survey US adults to examine general attitudes, interests, and concerns regarding PES use.

Design, Setting, and Participants

For this survey study, data were collected from 1 stratified sample and 1 nonprobability sample (samples 1 and 2, respectively) between March and July 2023. The surveys measured approval, interest, and concerns regarding various applications of PES. In the second sample, presentation of a list of potential concerns was randomized (presented at survey onset vs survey end). The survey was designed using Qualtrics and distributed to participants through Prolific, an online sampling firm. Sample 1 was nationally representative with respect to gender, age, and race and ethnicity; sample 2 was recruited without specific demographic criteria. Analyses were conducted between March 2023 and February 2024.

Main Outcomes and Measures

Participants reported their approval, interest, and concerns regarding various applications of PES and outcomes screened (eg, traits and conditions). Statistical analysis was conducted using independent samples t tests and repeated-measures analyses of variance.

Results

Of the 1435 respondents in sample 1, demographic data were available for 1427 (mean [SD] age, 45.8 [16.0] years; 724 women [50.7%]). Among these 1427 sample 1 respondents, 1027 (72.0%) expressed approval for PES and 1169 (81.9%) expressed some interest in using PES if already undergoing in vitro fertilization (IVF). Approval among these respondents for using PES for embryo selection was notably high for physical health conditions (1109 [77.7%]) and psychiatric health conditions (1028 [72.0%]). In contrast, there was minority approval for embryo selection based on PES for behavioral traits (514 [36.0%]) and physical traits (432 [30.3%]). Nevertheless, concerns about PES leading to false expectations and promoting eugenic practices were pronounced, with 787 of 1422 (55.3%) and 780 of 1423 (54.8%) respondents finding them very to extremely concerning, respectively. Sample 2 included 192 respondents (mean [SD] age 37.7 [12.2] years; 110 men [57.3%]). These respondents were presented concerns at survey onset (n = 95) vs survey end (n = 97), which was associated with less approval (28-percentage point decrease) and more uncertainty (24 percentage-point increase) but with only slightly higher disapproval (4 percentage-point increase).

Conclusions and Relevance

These findings suggest that it is critical for health care professionals and medical societies to consider and understand the perspectives of diverse stakeholders (eg, patients undergoing IVF, clinicians, and the general public), given the absence of regulation and the recent commercial availability of PES.

Multi-ancestry tandem repeat association study of hair colour using exome-wide sequencing

Hair colour variation is influenced by hundreds of positions across the human genome but this genetic contribution has only been narrowly explored. Genome-wide association studies identified single nucleotide polymorphisms (SNPs) influencing hair colour but the biology underlying these associations is challenging to interpret. We report 16 tandem repeats (TRs) with effects on different models of hair colour plus two TRs associated with hair colour in diverse ancestry groups. Several of these TRs expand or contract amino acid coding regions of their localized protein such that structure, and by extension function, may be altered. We also demonstrate that independent of SNP variation, these TRs can be used to great an additive polygenic score that predicts darker hair colour. This work adds to the growing body of evidence regarding TR influence on human traits with relatively large and independent effects relative to surrounding SNP variation.

Integrative functional genomic analyses identify genetic variants influencing skin pigmentation in Africans

Skin color is highly variable in Africans, yet little is known about the underlying molecular mechanism. Here we applied massively parallel reporter assays to screen 1,157 candidate variants influencing skin pigmentation in Africans and identified 165 single-nucleotide polymorphisms showing differential regulatory activities between alleles. We combine Hi-C, genome editing and melanin assays to identify regulatory elements for MFSD12, HMG20B, OCA2, MITF, LEF1, TRPS1, BLOC1S6 and CYB561A3 that impact melanin levels in vitro and modulate human skin color. We found that independent mutations in an OCA2 enhancer contribute to the evolution of human skin color diversity and detect signals of local adaptation at enhancers of MITF, LEF1 and TRPS1, which may contribute to the light skin color of Khoesan-speaking populations from Southern Africa. Additionally, we identified CYB561A3 as a novel pigmentation regulator that impacts genes involved in oxidative phosphorylation and melanogenesis. These results provide insights into the mechanisms underlying human skin color diversity and adaptive evolution.

The genomic architecture of continuous plumage colour variation in the European barn owl (Tyto alba)

The maintenance of colour variation in wild populations has long fascinated evolutionary biologists, although most studies have focused on discrete traits exhibiting rather simple inheritance patterns and genetic architectures. However, the study of continuous colour traits and their potentially oligo- or polygenic genetic bases remains rare in wild populations. We studied the genetics of the continuously varying white-to-rufous plumage coloration of the European barn owl (Tyto alba) using a genome-wide association approach on the whole-genome data of 75 individuals. We confirmed a mutation at the melanocortin-1-receptor gene (MC1R) is involved in the coloration and identified two new regions, located in super-scaffolds 9 and 42. The combination of the three regions explains most of the colour variation (80.37%, 95% credible interval 58.45-100%). One discovered region, located in the sex chromosome, differs between the most extreme colorations in owls sharing a specific MC1R genotype. This region may play a role in the colour sex dimorphism of this species, possibly in interaction with the autosomal MC1R. We thus provide insights into the genetic architecture of continuous colour variation, pointing to an oligogenic basis with potential epistatic effects among loci that should aid future studies understanding how continuous colour variation is maintained in nature.

Forensic DNA phenotyping and its politics of legitimation and contestation: Views of forensic geneticists in Europe

Forensic DNA Phenotyping (FDP) is a set of techniques that aim to infer externally visible characteristics in humans - such as eye, hair and skin color - and biogeographical ancestry of an unknown person, based on biological material. FDP has been applied in various jurisdictions in a limited number of high-profile cases to provide intelligence for criminal investigations. There are on-going controversies about the reliability and validity of FDP, which come together with debates about the ethical challenges emerging from the use of this technology in the criminal justice system. Our study explores how, in the context of complex politics of legitimation of and contestation over the use of FDP, forensic geneticists in Europe perceive this technology's potential applications, utility and risks. Forensic geneticists perform several forms of discursive boundary work, making distinctions between science and the criminal justice system, experts and non-experts, and good and bad science. Such forms of boundary work reconstruct the complex positioning vis-à-vis legal and scientific realities. In particular, while mobilizing interest in FDP, forensic geneticists simultaneously carve out notions of risk, accountability and scientific conduct that perform distance from FDP' implications in the criminal justice system.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

High-coverage genome of the Tyrolean Iceman reveals unusually high Anatolian farmer ancestry

The Tyrolean Iceman is known as one of the oldest human glacier mummies, directly dated to 3350-3120 calibrated BCE. A previously published low-coverage genome provided novel insights into European prehistory, despite high present-day DNA contamination. Here, we generate a high-coverage genome with low contamination (15.3×) to gain further insights into the genetic history and phenotype of this individual. Contrary to previous studies, we found no detectable Steppe-related ancestry in the Iceman. Instead, he retained the highest Anatolian-farmer-related ancestry among contemporaneous European populations, indicating a rather isolated Alpine population with limited gene flow from hunter-gatherer-ancestry-related populations. Phenotypic analysis revealed that the Iceman likely had darker skin than present-day Europeans and carried risk alleles associated with male-pattern baldness, type 2 diabetes, and obesity-related metabolic syndrome. These results corroborate phenotypic observations of the preserved mummified body, such as high pigmentation of his skin and the absence of hair on his head.

Evaluation of loci to predict ear morphology using two SNaPshot assays

Human ear morphology prediction with SNP-based genotypes is growing in forensic DNA phenotyping and is scarcely explored in Pakistan as a part of EVCs (externally visible characteristics). The ear morphology prediction assays with 21 SNPs were assessed for their potential utility in forensic identification of population. The SNaPshot™ multiplex chemistries, capillary electrophoresis methods and GeneMapper™ software were used for obtaining genotypic data. A total of 33 ear phenotypes were categorized with digital photographs of 300 volunteers. SHEsis software was applied to make LD plot. Ordinal and multinomial logistic regression was implemented for association testing. Multinomial logistic regression was executed to construct the prediction model in 90% training and 10% testing subjects. Several influential SNPs for ear phenotypic variation were found in association testing. The model based on genetic markers predicted ear phenotypes with moderate to good predictive accuracies demonstrated with the area under curve (AUC), sensitivity and specificity of predicted phenotypes. As an additional EVC, the estimated ear phenotypic profiles have the possibility of determining the human ear morphology differences in unknown biological samples found in crimes that do not result in a criminal database hit. Furthermore, this can help in facial reconstruction and act as an investigational lead.

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.

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.

Whole genome sequencing reveals signals of adaptive admixture in Creole cattle

The Creole cattle from Guadeloupe (GUA) are well adapted to the tropical environment. Its admixed genome likely played an important role in such adaptation. Here, we sought to detect genomic signatures of selection in the GUA genome. For this purpose, we sequenced 23 GUA individuals and combined our data with sequenced genomes of 99 animals representative of European, African and indicine groups. We detect 17,228,983 single nucleotide polymorphisms (SNPs) in the GUA genome, providing the most detailed exploration, to date, of patterns of genetic variation in this breed. We confirm the higher level of African and indicine ancestries, compared to the European ancestry and we highlight the African origin of indicine ancestry in the GUA genome. We identify five strong candidate regions showing an excess of indicine ancestry and consistently supported across the different detection methods. These regions encompass genes with adaptive roles in relation to immunity, thermotolerance and physical activity. We confirmed a previously identified horn-related gene, RXFP2, as a gene under strong selective pressure in the GUA population likely owing to human-driven (socio-cultural) pressure. Findings from this study provide insight into the genetic mechanisms associated with resilience traits in livestock.

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.

Fine-mapping of retinal vascular complexity loci identifies Notch regulation as a shared mechanism with myocardial infarction outcomes

There is increasing evidence that the complexity of the retinal vasculature measured as fractal dimension, Df, might offer earlier insights into the progression of coronary artery disease (CAD) before traditional biomarkers can be detected. This association could be partly explained by a common genetic basis; however, the genetic component of Df is poorly understood. We present a genome-wide association study (GWAS) of 38,000 individuals with white British ancestry from the UK Biobank aimed to comprehensively study the genetic component of Df and analyse its relationship with CAD. We replicated 5 Df loci and found 4 additional loci with suggestive significance (P < 1e-05) to contribute to Df variation, which previously were reported in retinal tortuosity and complexity, hypertension, and CAD studies. Significant negative genetic correlation estimates support the inverse relationship between Df and CAD, and between Df and myocardial infarction (MI), one of CAD's fatal outcomes. Fine-mapping of Df loci revealed Notch signalling regulatory variants supporting a shared mechanism with MI outcomes. We developed a predictive model for MI incident cases, recorded over a 10-year period following clinical and ophthalmic evaluation, combining clinical information, Df, and a CAD polygenic risk score. Internal cross-validation demonstrated a considerable improvement in the area under the curve (AUC) of our predictive model (AUC = 0.770 ± 0.001) when comparing with an established risk model, SCORE, (AUC = 0.741 ± 0.002) and extensions thereof leveraging the PRS (AUC = 0.728 ± 0.001). This evidences that Df provides risk information beyond demographic, lifestyle, and genetic risk factors. Our findings shed new light on the genetic basis of Df, unveiling a common control with MI, and highlighting the benefits of its application in individualised MI risk prediction.

Skin Phototype and Disease: A Comprehensive Genetic Approach to Pigmentary Traits Pleiotropy Using PRS in the GCAT Cohort

Human pigmentation has largely been associated with different disease prevalence among populations, but most of these studies are observational and inconclusive. Known to be genetically determined, pigmentary traits have largely been studied by Genome-Wide Association Study (GWAS), mostly in Caucasian ancestry cohorts from North Europe, identifying robustly, several loci involved in many of the pigmentary traits. Here, we conduct a detailed analysis by GWAS and Polygenic Risk Score (PRS) of 13 pigmentary-related traits in a South European cohort of Caucasian ancestry (n = 20,000). We observed fair phototype strongly associated with non-melanoma skin cancer and other dermatoses and confirmed by PRS-approach the shared genetic basis with skin and eye diseases, such as melanoma (OR = 0.95), non-melanoma skin cancer (OR = 0.93), basal cell carcinoma (OR = 0.97) and darker phototype with vitiligo (OR = 1.02), cataracts (OR = 1.04). Detailed genetic analyses revealed 37 risk loci associated with 10 out of 13 analyzed traits, and 16 genes significantly associated with at least two pigmentary traits. Some of them have been widely reported, such as MC1R, HERC2, OCA2, TYR, TYRP1, SLC45A2, and some novel candidate genes C1QTNF3, LINC02876, and C1QTNF3-AMACR have not been reported in the GWAS Catalog, with regulatory potential. These results highlight the importance of the assess phototype as a genetic proxy of skin functionality and disease when evaluating open mixed populations.

Phenome-wide association study of loci harboring de novo tandem repeat mutations in UK Biobank exomes

When present in coding regions, tandem repeats (TRs) may have large effects on protein structure and function contributing to health and disease. We use a family-based design to identify de novo TRs and assess their impact at the population level in 148,607 European ancestry participants from the UK Biobank. The 427 loci with de novo TR mutations are enriched for targets of microRNA-184 (21.1-fold, P = 4.30 × 10-5, FDR = 9.50 × 10-3). There are 123 TR-phenotype associations with posterior probabilities > 0.95. These relate to body structure, cognition, and cardiovascular, metabolic, psychiatric, and respiratory outcomes. We report several loci with large likely causal effects on tissue microstructure, including the FAN1-[TG]N and carotid intima-media thickness (mean thickness: beta = 5.22, P = 1.22 × 10-6, FDR = 0.004; maximum thickness: beta = 6.44, P = 1.12 × 10-6, FDR = 0.004). Two exonic repeats FNBP4-[GGT]N and BTN2A1-[CCT]N alter protein structure. In this work, we contribute clear and testable hypotheses of dose-dependent TR implications linking genetic variation and protein structure with health and disease outcomes.

Predicting Eye and Hair Color in a Turkish Population Using the HIrisPlex System

Forensic DNA Phenotyping (FDP) can reveal the appearance of an unknown individual by predicting the ancestry, phenotype (i.e., hair, eye, skin color), and age from DNA obtained at the crime scene. The HIrisPlex system has been developed to simultaneously predict eye and hair color. However, the prediction accuracy of the system needs to be assessed for the tested population before implementing FDP in casework. In this study, we evaluated the performance of the HIrisPlex system on 149 individuals from the Turkish population. We applied the single-based extension (SNaPshot chemistry) method and used the HIrisPlex online tool to test the prediction of the eye and hair colors. The accuracy of the HIrisPlex system was assessed through the calculation of the area under the receiver characteristic operating curves (AUC), sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). The results showed that the proposed method successfully predicted the eye and hair color, especially for blue (100%) and brown (95.60%) eye and black (95.23) and brown (98.94) hair colors. As observed in previous studies, the system failed to predict intermediate eye color, representing 25% in our cohort. The majority of incorrect predictions were observed for blond hair color (40.7%). Previous HIrisPlex studies have also noted difficulties with these phenotypes. Our study shows that the HIrisPlex system can be applied to forensic casework in Turkey with careful interpretation of the data, particularly intermediate eye color and blond hair color.

A systematic review of skin ageing genes: gene pleiotropy and genes on the chromosomal band 16q24.3 may drive skin ageing

Skin ageing is the result of intrinsic genetic and extrinsic lifestyle factors. However, there is no consensus on skin ageing phenotypes and ways to quantify them. In this systematic review, we first carefully identified 56 skin ageing phenotypes from multiple literature sources and sought the best photo-numeric grading scales to evaluate them. Next, we conducted a systematic review on all 44 Genome-wide Association Studies (GWAS) on skin ageing published to date and identified genetic risk factors (2349 SNPs and 366 genes) associated with skin ageing. We identified 19 promising SNPs found to be significantly (p-Value < 1E-05) associated with skin ageing phenotypes in two or more independent studies. Here we show, using enrichment analyses strategies and gene expression data, that (1) pleiotropy is a recurring theme among skin ageing genes, (2) SNPs associated with skin ageing phenotypes are mostly located in a small handful of 44 pleiotropic and hub genes (mostly on the chromosome band 16q24.3) and 32 skin colour genes. Since numerous genes on the chromosome band 16q24.3 and skin colour genes show pleiotropy, we propose that (1) genes traditionally identified to contribute to skin colour have more than just skin pigmentation roles, and (2) further progress towards understand the development of skin pigmentation requires understanding the contributions of genes on the chromosomal band 16q24.3. We anticipate our systematic review to serve as a hub to locate primary literature sources pertaining to the genetics of skin ageing and to be a starting point for more sophisticated work examining pleiotropic genes, hub genes, and skin ageing phenotypes.

Genome-Wide Association Study Identifies Two Common Loci Associated with Pigment Dispersion Syndrome/Pigmentary Glaucoma and Implicates Myopia in its Development

PURPOSE

To identify genetic variants associated with pigment dispersion syndrome (PDS) and pigmentary glaucoma (PG) in unrelated patients and to further understand the genetic and potentially causal relationships between PDS and associated risk factors.

DESIGN

A 2-stage genome-wide association meta-analysis with replication and subsequent in silico analyses including Mendelian randomization.

PARTICIPANTS

A total of 574 cases with PG or PDS and 52 627 controls of European descent.

METHODS

Genome-wide association analyses were performed in 4 cohorts and meta-analyzed in 3 stages: (1) a discovery meta-analysis was performed in 3 cohorts, (2) replication was performed in the fourth cohort, and (3) all 4 cohorts were meta-analyzed to increase statistical power. Two-sample Mendelian randomization was used to determine whether refractive error and intraocular pressure exert causal effects over PDS.

MAIN OUTCOME MEASURES

The association of genetic variants with PDS and whether myopia exerts causal effects over PDS.

RESULTS

Significant association was present at 2 novel loci for PDS/PG. These loci and follow-up analyses implicate the genes gamma secretase activator protein (GSAP) (lead single nucleotide polymorphism [SNP]: rs9641220, P = 6.0×10-10) and glutamate metabotropic receptor 5 (GRM5)/TYR (lead SNP: rs661177, P = 3.9×10-9) as important factors in disease risk. Mendelian randomization showed significant evidence that negative refractive error (myopia) exerts a direct causal effect over PDS (P = 8.86×10-7).

CONCLUSIONS

Common SNPs relating to the GSAP and GRM5/TYR genes are associated risk factors for the development of PDS and PG. Although myopia is a known risk factor, this study uses genetic data to demonstrate that myopia is, in part, a cause of PDS and PG.

A genome-wide epistatic network underlies the molecular architecture of continuous color variation of body extremities

Deciphering the molecular architecture of coat coloration for a better understanding of the biological mechanisms underlying pigmentation still remains a challenge. We took advantage of a rabbit French experimental population in which both a pattern and a gradient of coloration from white to brown segregated within the himalayan phenotype. The whole experimental design was genotyped using the high density Affymetrix® AxiomOrcun™ SNP Array and phenotyped into 6 different groups ordered from the lighter to the darker. Genome-wide association analyses pinpointed an oligogenic determinism, under recessive and additive inheritance, involving genes already known in melanogenesis (ASIP, KIT, MC1R, TYR), and likely processed pseudogenes linked to ribosomal function, RPS20 and RPS14. We also identified (i) gene-gene interactions through ASIP:MC1R affecting light cream/beige phenotypes while KIT:RPS responsible of dark chocolate/brown colors and (ii) a genome-wide epistatic network involving several others coloration genes such as POT1 or HPS5. Finally, we determined the recessive inheritance of the English spotting phenotype likely involving a copy number variation affecting at least the end of the coding sequence of the KIT gene. Our analyses of coloration as a continuous trait allowed us to go beyond much of the established knowledge through the detection of additional genes and gene-gene interactions that may contribute to the molecular architecture of the coloration phenotype.

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

Genome-wide association study of actinic keratosis identifies new susceptibility loci implicated in pigmentation and immune regulation pathways

Actinic keratosis (AK) is a common precancerous cutaneous neoplasm that arises on chronically sun-exposed skin. AK susceptibility has a moderate genetic component, and although a few susceptibility loci have been identified, including IRF4, TYR, and MC1R, additional loci have yet to be discovered. We conducted a genome-wide association study of AK in non-Hispanic white participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (n = 63,110, discovery cohort), with validation in the Mass-General Brigham (MGB) Biobank cohort (n = 29,130). We identified eleven loci (P < 5 × 10-8), including seven novel loci, of which four novel loci were validated. In a meta-analysis (GERA + MGB), one additional novel locus, TRPS1, was identified. Genes within the identified loci are implicated in pigmentation (SLC45A2, IRF4, BNC2, TYR, DEF8, RALY, HERC2, and TRPS1), immune regulation (FOXP1 and HLA-DQA1), and cell signaling and tissue remodeling (MMP24) pathways. Our findings provide novel insight into the genetics and pathogenesis of AK susceptibility.

Citrus-Gene interaction and melanoma risk in the UK Biobank

High citrus consumption may increase melanoma risk; however, little is known about the biological mechanisms of this association, or whether it is modified by genetic variants. We conducted a genome-wide analysis of gene-citrus consumption interactions on melanoma risk among 1563 melanoma cases and 193 296 controls from the UK Biobank. Both the 2-degrees-of-freedom (df) joint test of genetic main effect and gene-environment (G-E) interaction and the standard 1-df G-E interaction test were performed. Three index SNPs (lowest P-value SNP among highly correlated variants [r²  > .6]) were identified from among the 365 genome-wide significant 2-df test results (rs183783391 on chromosome 3 [MITF], rs869329 on chromosome 9 [MTAP] and rs11446223 on chromosome 16 [DEF8]). Although all three were statistically significant for the 2-df test (4.25e-08, 1.98e-10 and 4.93e-13, respectively), none showed evidence of interaction according to the 1-df test (P = .73, .24 and .12, respectively). Eight nonindex, 2-df test significant SNPs on chromosome 16 were significant (P < .05) according to the 1-df test, providing evidence of citrus-gene interaction. Seven of these SNPs were mapped to AFG3L1P (rs199600347, rs111822773, rs113178244, rs3803683, rs73283867, rs78800020, rs73283871), and one SNP was mapped to GAS8 (rs74583214). We identified several genetic loci that may elucidate the association between citrus consumption and melanoma risk. Further studies are needed to confirm these findings.

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

Bayesian model comparison for rare-variant association studies

Whole-genome sequencing studies applied to large populations or biobanks with extensive phenotyping raise new analytic challenges. The need to consider many variants at a locus or group of genes simultaneously and the potential to study many correlated phenotypes with shared genetic architecture provide opportunities for discovery not addressed by the traditional one variant, one phenotype association study. Here, we introduce a Bayesian model comparison approach called MRP (multiple rare variants and phenotypes) for rare-variant association studies that considers correlation, scale, and direction of genetic effects across a group of genetic variants, phenotypes, and studies, requiring only summary statistic data. We apply our method to exome sequencing data (n = 184,698) across 2,019 traits from the UK Biobank, aggregating signals in genes. MRP demonstrates an ability to recover signals such as associations between PCSK9 and LDL cholesterol levels. We additionally find MRP effective in conducting meta-analyses in exome data. Non-biomarker findings include associations between MC1R and red hair color and skin color, IL17RA and monocyte count, and IQGAP2 and mean platelet volume. Finally, we apply MRP in a multi-phenotype setting; after clustering the 35 biomarker phenotypes based on genetic correlation estimates, we find that joint analysis of these phenotypes results in substantial power gains for gene-trait associations, such as in TNFRSF13B in one of the clusters containing diabetes- and lipid-related traits. Overall, we show that the MRP model comparison approach improves upon useful features from widely used meta-analysis approaches for rare-variant association analyses and prioritizes protective modifiers of disease risk.

Prognostic prediction and gene regulation network of EIF2S2 in hepatocellular carcinoma based on data mining

BACKGROUND

Hepatocellular carcinoma (HCC) is a malignant tumor with a high fatality rate, predicting poor prognosis and therapeutic effect. Screening potential prognostic genes in HCC could be a creative way to advance clinical treatment. Eukaryotic translation initiation factor 2 subunit beta (EIF2S2) has reportedly been linked to several tumors, including liver cancer, but the prognostic predictions remain unknown. Therefore, we aimed to clarify the prognostic role and interaction network of EIF2S2 in HCC using bioinformatics data.

METHODS

We screened EIF2S2 using the Oncomine, Ualcan, and TCGA databases. R software was used to analyze the mRNA level and clinicopathological characteristics of hepatocellular carcinoma. Evaluation of the correlations between EIF2S2 and patients' survival was made using the Kaplan-Meier curves and Cox proportional hazards regression model. Then, the influence of EIF2S2 gene mutations on the prognosis of patients was explored by cBioPortal. The protein-protein interaction network of 50 similar genes related to EIF2S2 was implemented by GEPIA2 and Metascape. The LinkedOmics database allowed us to carry out Gene Set Enrichment Analysis. Finally, we constructed the EIF2S2 kinase, miRNA, and transcription factor target networks using GeneMANIA.

RESULTS

EIF2S2 mRNA was overexpressed in HCC and was closely associated with clinicopathological features, including gender, age, race, tumor grade, and stage. There was no correlation between EIF2S2 genetic mutations and prognostic survival. Combining Cox proportional hazards regression model analyses, high-expressed EIF2S2 predicted poor prognosis in HCC patients. Additionally, we screened the top three EIF2S2-related genes (PFDN4, HM13, and SNRPD1), the 50 similar genes, and then constructed a 50-similar-gene protein-protein interaction network identified by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using Metascape. EIF2S2 target networks in HCC were identified in kinase, miRNA, and transcription factor networks, including the mitogen-activated protein kinase 1 (MAPK1), miRNAs (Mir-144), and transcription factors (GGAANCGGAANY_UNKNOWN) using GeneMANIA.

CONCLUSIONS

EIF2S2 plays a crucial role in the gene-regulating network of HCC and may be a potential prognostic marker or therapeutic target for HCC patients.

Association of Melanoma-Risk Variants with Primary Melanoma Tumor Prognostic Characteristics and Melanoma-Specific Survival in the GEM Study

Genome-wide association studies (GWAS) and candidate pathway studies have identified low-penetrant genetic variants associated with cutaneous melanoma. We investigated the association of melanoma-risk variants with primary melanoma tumor prognostic characteristics and melanoma-specific survival. The Genes, Environment, and Melanoma Study enrolled 3285 European origin participants with incident invasive primary melanoma. For each of 47 melanoma-risk single nucleotide polymorphisms (SNPs), we used linear and logistic regression modeling to estimate, respectively, the per allele mean changes in log of Breslow thickness and odds ratios for presence of ulceration, mitoses, and tumor-infiltrating lymphocytes (TILs). We also used Cox proportional hazards regression modeling to estimate the per allele hazard ratios for melanoma-specific survival. Passing the false discovery threshold (p = 0.0026) were associations of IRF4 rs12203592 and CCND1 rs1485993 with log of Breslow thickness, and association of TERT rs2242652 with presence of mitoses. IRF4 rs12203592 also had nominal associations (p < 0.05) with presence of mitoses and melanoma-specific survival, as well as a borderline association (p = 0.07) with ulceration. CCND1 rs1485993 also had a borderline association with presence of mitoses (p = 0.06). MX2 rs45430 had nominal associations with log of Breslow thickness, presence of mitoses, and melanoma-specific survival. Our study indicates that further research investigating the associations of these genetic variants with underlying biologic pathways related to tumor progression is warranted.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

Quantitative Human Paleogenetics: What can Ancient DNA Tell us About Complex Trait Evolution?

Genetic association data from national biobanks and large-scale association studies have provided new prospects for understanding the genetic evolution of complex traits and diseases in humans. In turn, genomes from ancient human archaeological remains are now easier than ever to obtain, and provide a direct window into changes in frequencies of trait-associated alleles in the past. This has generated a new wave of studies aiming to analyse the genetic component of traits in historic and prehistoric times using ancient DNA, and to determine whether any such traits were subject to natural selection. In humans, however, issues about the portability and robustness of complex trait inference across different populations are particularly concerning when predictions are extended to individuals that died thousands of years ago, and for which little, if any, phenotypic validation is possible. In this review, we discuss the advantages of incorporating ancient genomes into studies of trait-associated variants, the need for models that can better accommodate ancient genomes into quantitative genetic frameworks, and the existing limits to inferences about complex trait evolution, particularly with respect to past populations.

A custom capture sequence approach for oculocutaneous albinism identifies structural variant alleles at the OCA2 locus

Oculocutaneous albinism (OCA) is a heritable disorder of pigment production that manifests as hypopigmentation and altered eye development. Exon sequencing of known OCA genes is unsuccessful in producing a complete molecular diagnosis for a significant number of affected individuals. We sequenced the DNA of individuals with OCA using short-read custom capture sequencing that targeted coding, intronic, and noncoding regulatory regions of known OCA genes, and genome-wide association study-associated pigmentation loci. We identified an OCA2 complex structural variant (CxSV), defined by a 143 kb inverted segment reintroduced in intron 1, upstream of the native location. The corresponding CxSV junctions were observed in 11/390 probands screened. The 143 kb CxSV presents in one family as a copy number variant duplication for the 143 kb region. In the remaining 10/11 families, the 143 kb CxSV acquired an additional 184 kb deletion across the same region, restoring exons 3-19 of OCA2 to a copy-number neutral state. Allele-associated haplotype analysis found rare SNVs rs374519281 and rs139696407 are linked with the 143 kb CxSV in both OCA2 alleles. For individuals in which customary molecular evaluation does not reveal a biallelic OCA diagnosis, we recommend preliminary screening for these haplotype-associated rare variants, followed by junction-specific validation for the OCA2 143 kb CxSV.

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.

Evolutionary genetics of skin pigmentation in African populations

Skin color is a highly heritable human trait, and global variation in skin pigmentation has been shaped by natural selection, migration and admixture. Ethnically diverse African populations harbor extremely high levels of genetic and phenotypic diversity, and skin pigmentation varies widely across Africa. Recent genome-wide genetic studies of skin pigmentation in African populations have advanced our understanding of pigmentation biology and human evolutionary history. For example, novel roles in skin pigmentation for loci near MFSD12 and DDB1 have recently been identified in African populations. However, due to an underrepresentation of Africans in human genetic studies, there is still much to learn about the evolutionary genetics of skin pigmentation. Here, we summarize recent progress in skin pigmentation genetics in Africans and discuss the importance of including more ethnically diverse African populations in future genetic studies. In addition, we discuss methods for functional validation of adaptive variants related to skin pigmentation.

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.

Skin pigmentation and its control: From ultraviolet radiation to stem cells

In the light of substantial discoveries in epithelial and hair pigmentation pathophysiology, this review summarizes the current understanding of skin pigmentation mechanisms. Melanocytes are pigment-producing cells, and their key regulating transcription factor is the melanocyte-specific microphthalmia-associated transcription factor (m-MITF). Ultraviolet (UV) radiation is a unique modulator of skin pigmentation influencing tanning pathways. The delayed tanning pathway occurs as UVB produces keratinocyte DNA damage, causing p53-mediated expression of the pro-opiomelanocortin (POMC) gene that is processed to release α-melanocyte-stimulating hormone (α-MSH). α-MSH stimulates the melanocortin 1 receptor (MC1R) on melanocytes, leading to m-MITF expression and melanogenesis. POMC cleavage also releases β-endorphin, which creates a neuroendocrine pathway that promotes UV-seeking behaviours. Mutations along the tanning pathway can affect pigmentation and increase the risk of skin malignancies. MC1R variants have received considerable attention, yet the allele is highly polymorphic with varied phenotypes. Vitiligo presents with depigmented skin lesions due to autoimmune destruction of melanocytes. UVB phototherapy stimulates melanocyte stem cells in the hair bulge to undergo differentiation and upwards migration resulting in perifollicular repigmentation of vitiliginous lesions, which is under sophisticated signalling control. Melanocyte stem cells, normally quiescent, undergo cyclic activation/differentiation and downward migration with the hair cycle, providing pigment to hair follicles. Physiological hair greying results from progressive loss of melanocyte stem cells and can be accelerated by acute stress-induced, sympathetic driven hyperproliferation of the melanocyte stem cells. Ultimately, by reviewing the pathways governing epithelial and follicular pigmentation, numerous areas of future research and potential points of intervention are highlighted.

Selection on a small genomic region underpins differentiation in multiple color traits between two warbler species

Speciation is one of the most important processes in biology, yet the study of the genomic changes underlying this process is in its infancy. North American warbler species Setophaga townsendi and Setophaga occidentalis hybridize in a stable hybrid zone, following a period of geographic separation. Genomic differentiation accumulated during geographic isolation can be homogenized by introgression at secondary contact, whereas genetic regions that cause low hybrid fitness can be shielded from such introgression. Here, we examined the genomic underpinning of speciation by investigating (1) the genetic basis of divergent pigmentation traits between species, (2) variation in differentiation across the genome, and (3) the evidence for selection maintaining differentiation in the pigmentation genes. Using tens of thousands of single nucleotide polymorphisms (SNPs) genotyped in hundreds of individuals within and near the hybrid zone, genome-wide association mapping revealed a single SNP associated with cheek, crown, breast coloration, and flank streaking, reflecting pleiotropy (one gene affecting multiple traits) or close physical linkage of different genes affecting different traits. This SNP is within an intron of the RALY gene, hence we refer to it as the RALY SNP. We then examined between-species genomic differentiation, using both genotyping-by-sequencing and whole genome sequencing. We found that the RALY SNP is within one of the highest peaks of differentiation, which contains three genes known to influence pigmentation: ASIP, EIF2S2, and RALY (the ASIP-RALY gene block). Heterozygotes at this gene block are likely of reduced fitness, as the geographic cline of the RALY SNP has been narrow over two decades. Together, these results reflect at least one barrier to gene flow within this narrow (∼200 kb) genomic region that modulates plumage difference between species. Despite extensive gene flow between species across the genome, this study provides evidence that selection on a phenotype-associated genomic region maintains a stable species boundary.

Overcoming Immune Evasion in Melanoma

Melanoma is the most aggressive and dangerous form of skin cancer that develops from transformed melanocytes. It is crucial to identify melanoma at its early stages, in situ, as it is "curable" at this stage. However, after metastasis, it is difficult to treat and the five-year survival is only 25%. In recent years, a better understanding of the etiology of melanoma and its progression has made it possible for the development of targeted therapeutics, such as vemurafenib and immunotherapies, to treat advanced melanomas. In this review, we focus on the molecular mechanisms that mediate melanoma development and progression, with a special focus on the immune evasion strategies utilized by melanomas, to evade host immune surveillances. The proposed mechanism of action and the roles of immunotherapeutic agents, ipilimumab, nivolumab, pembrolizumab, and atezolizumab, adoptive T- cell therapy plus T-VEC in the treatment of advanced melanoma are discussed. In this review, we implore that a better understanding of the steps that mediate melanoma onset and progression, immune evasion strategies exploited by these tumor cells, and the identification of biomarkers to predict treatment response are critical in the design of improved strategies to improve clinical outcomes for patients with this deadly disease.