A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia

Polygenic Risk Score Analysis of 37 SNPs Associated with Melanoma Risk in Colombian Population

Melanoma incidence is increasing, with distinct genetic and clinical patterns observed in the Latin American population. This study aimed to evaluate melanoma risk in a Colombian cohort through polygenic risk analysis using 37 variants across nine genes previously associated with melanoma. We performed polygenic risk score (PRS) analysis on 85 melanoma patients and 165 controls. Genotyping was performed for 37 melanoma-associated SNPs, and on the basis of previous GWAS reports, individual PRSs were calculated for each participant. The participants were then stratified into quartiles to examine risk gradients. In addition, phenotypic features such as eye and hair color were evaluated, and genetic models and haplotype analyses were performed, adjusting for sex and family history of cancer. PRS quartile stratification revealed a clear risk gradient. Notably, 31.8% of the melanoma cases were clustered in the highest-risk quartile (Q4), with a maximum PRS of 1.04. Variants in TYR, TYRP1, CDKN2A, and HERC2 significantly contributed to risk, and light brown eye and hair colors were strongly associated with increased melanoma risk. Moreover, a protective haplotype in the OCA2-HERC2 region was identified among males. The integration of the PRS with clinical and phenotypic factors has potential for improving melanoma risk stratification in the Colombian population, warranting further investigation in larger, diverse cohorts.

Choice between DNA primer sets (A or B) of the ForenSeq kit: forensic evaluation in a Mexican admixed population sample

Massively parallel sequencing (MPS) overcomes many PCR-CE limitations to analyze STRs and allow simultaneous inclusion of SNPs in forensic cases. By MPS, the ForenSeq™ DNA Signature Prep kit analyzes 27 aSTRs, 7 X-STRs, 24Y-STRs, and 94 identity-informative SNPs (iiSNPs) with the DNA Primer Set-A (DPS-A). Optionally, the DNA Primer Set-B (DPS-B) adds to the analysis 56 ancestry-informative SNPs (aiSNPs) and 24 phenotype-informative SNPs (piSNPs), but diminishes from 96 to 32 the number of samples per sequencing run. We assessed the forensic informativity provided by the loci analyzed by these two DPS in admixed individuals from Mexico City (Center, Mexico). For STRs, we report length-based (LB) and sequence-based (SB) allele frequencies and forensic parameters of the 152 identity informative markers (DPS-A). For aSTRs, the combined PD of SB genotypes (PD ~ 100%) was ~ 2949 times larger than that from LB. Conversely, the observed phenotype distribution offered low PD levels (PD = 6.6% and 10.4%), whereas piSNPs predicted accurately only the modal brown eye and dark hair colors, respectively. Similarly, aiSNPs detected a large prevalence of admixed individuals (97.3%; PD = 5.4%). Although few individuals were inferred as Europeans and Native Americans (1.37% each), they were self-declared as admixed, which result confusing for HID purposes. In brief, SB genotypes increased significantly the informativity of STRs to solve complex cases (DPS-A), whereas aiSNPs and piSNPs added mostly irrelevant information (DPS-B). We provide useful cost-benefit criteria in one Latin American population to choose DPS-A (96 samples) instead of DPS-B (32 samples) of the Forenseq kit.

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.

Genetics of Latin American Diversity Project: Insights into population genetics and association studies in admixed groups in the Americas

Latin Americans are underrepresented in genetic studies, increasing disparities in personalized genomic medicine. Despite available genetic data from thousands of Latin Americans, accessing and navigating the bureaucratic hurdles for consent or access remains challenging. To address this, we introduce the Genetics of Latin American Diversity (GLAD) Project, compiling genome-wide information from 53,738 Latin Americans across 39 studies representing 46 geographical regions. Through GLAD, we identified heterogeneous ancestry composition and recent gene flow across the Americas. Additionally, we developed GLAD-match, a simulated annealing-based algorithm, to match the genetic background of external samples to our database, sharing summary statistics (i.e., allele and haplotype frequencies) without transferring individual-level genotypes. Finally, we demonstrate the potential of GLAD as a critical resource for evaluating statistical genetic software in the presence of admixture. By providing this resource, we promote genomic research in Latin Americans and contribute to the promises of personalized medicine to more people.

Selection scan in Native Americans of Mexico identifies FADS2 rs174616: Evidence of gene-diet interactions affecting lipid levels and Delta-6-desaturase activity

Searching for positive selection signals across genomes has identified functional genetic variants responding to environmental change. In Native Americans of Mexico, we used the fixation index (Fst) and population branch statistic (PBS) to identify SNPs suggesting positive selection. The 103 most differentiated SNPs were tested for associations with metabolic traits, the most significant association was FADS2/rs174616 with body mass index (BMI). This variant lies within a linkage disequilibrium (LD) block independent of previously reported FADS selection signals and has not been clearly associated with metabolic phenotypes. We tested this variant in two independent cohorts with cardiometabolic data. In the Genetics of Atherosclerotic Disease (GEA) cohort, the derived allele (T) was associated with increased BMI, lower LDL-C levels and a decreased risk of subclinical atherosclerosis in women. Significant gene-diet interactions affected lipid, apolipoprotein and adiponectin levels with differences according to sex, involving mainly total and complex dietary carbohydrate%. In the Genotype-related Effects of PUFA trial, the derived allele was associated with lower Δ-6 desaturase activity and erythrocyte membrane dihomo-gamma-linolenic acid (DGLA) levels, and with increased Δ-5 desaturase activity and eicosapentaenoic acid levels. This variant interacted with dietary carbohydrate% affecting Δ-6 desaturase activity. Notably, the relationship of DGLA and other erythrocyte membrane LC-PUFA indices with HOMA-IR differed according to rs174616 genotype, which has implications regarding how these indices should be interpreted. In conclusion, this observational study identified rs174616 as a signal suggesting selection in an independent linkage disequilibrium block, was associated with cardiometabolic and erythrocyte measurements of LC-PUFA in two independent Mexican cohorts and showed significant gene-diet interactions.

Polygenic Indices (a.k.a. Polygenic Scores) in Social Science: A Guide for Interpretation and Evaluation

Polygenic indices (PGI)-the new recommended label for polygenic scores (PGS) in social science-are genetic summary scales often used to represent an individual's liability for a disease, trait, or behavior based on the additive effects of measured genetic variants. Enthusiasm for linking genetic data with social outcomes and the inclusion of premade PGIs in social science datasets have facilitated increased uptake of PGIs in social science research-a trend that will likely continue. Yet, most social scientists lack the expertise to interpret and evaluate PGIs in social science research. Here, we provide a primer on PGIs for social scientists focusing on key concepts, unique statistical genetic considerations, and best practices in calculation, estimation, reporting, and interpretation. We summarize our recommended best practices as a checklist to aid social scientists in evaluating and interpreting studies with PGIs. We conclude by discussing the similarities between PGIs and standard social science scales and unique interpretative considerations.

Skin pigmentation related variants in Mexican population and interaction effects on serum 25(OH)D concentration and vitamin D deficiency

Skin pigmentation is negatively associated with circulating vitamin D (VD) concentration. Therefore, genetic factors involved in skin pigmentation could influence the risk of vitamin D deficiency (VDD). We evaluated the impact genetic variants related to skin pigmentation on VD in Mexican population. This cross-sectional analysis included 848 individuals from the Health Worker Cohort Study (ratio males to females ~ 1:3). Eight genetic variants: rs16891982 (SLC45A2), rs12203592 (IRF4), rs1042602 and rs1126809 (TYR), rs1800404 (OCA2), rs12913832 (HERC2), rs1426654 (SLC24A5), and rs2240751 (MFSD12); involved in skin pigmentation were genotyped. Skin pigmentation was assessed by self-report. Linear and logistic regression were used to assess the association between the variants of interest and VD and VDD, as appropriate. In our study, eight genetic variants were associated with skin pigmentation. A genetic risk score built with the variants rs1426654 and rs224075 was associated with lower VD levels (β = - 1.38, 95% CI - 2.59, - 0.17, p = 0.025). Nevertheless, when examining gene-gene interactions, we observed that rs2240751 × rs12203592 were associated with VD levels (P interaction = 0.021). Whereas rs2240751 × rs12913832 (P interaction = 0.0001) were associated with VDD. Our results suggest that skin pigmentation-related gene variants are associated with lower VD levels in Mexican population. These results underscore the importance of considering genetic interactions when assessing the impact of genetic polymorphisms on VD levels.

Proteomic Biomarkers Associated with Low Bone Mineral Density: A Systematic Review

Osteoporosis is a globally relevant public health issue. Our study aimed to summarize the knowledge on the proteomic biomarkers for low bone mineral density over the last years. We conducted a systematic review following the PRISMA guidelines; the scoured databases were PubMed, Web of Sciences, Scopus, and EBSCO, from inception to 2 June 2023. A total of 610 relevant studies were identified and 33 were assessed for eligibility. Finally, 29 studies met the criteria for this systematic review. The risk of bias was evaluated using the Joanna Briggs Institute Critical Appraisal Checklist tool. From the studies selected, 154 proteins were associated with changes of bone mineral density, from which only 10 were reported in at least two articles. The protein-protein network analysis indicated potential biomarkers involved in the skeletal system, immune system process, regulation of protein metabolic process, regulation of signaling, transport, cellular component assembly, cell differentiation, hemostasis, and extracellular matrix organization. Mass spectrometry-based proteomic profiling has allowed the discovery of new biomarkers with diagnostic potential. However, it is necessary to compare and validate the potential biomarkers in different populations to determine their association with bone metabolism and evaluate their translation to the clinical management of osteoporosis.

Weakened tanning ability is an important mechanism for evolutionary skin lightening in East Asians

The evolution of light-skin pigmentation among Eurasians is considered as an adaptation to the high-latitude environments. East Asians are ideal populations for studying skin color evolution because of the complex environment of East Asia. Here, we report a strong selection signal for the pigmentation gene phenylalanine hydroxylase (PAH) in light-skinned Han Chinese individuals. The intron mutation rs10778203 in PAH is enriched in East Asians and is significantly associated with skin color of the back of the hand in Han Chinese males (P < 0.05). In vitro luciferase and transcription factor binding assays show that the ancestral allele of rs10778203 could bind to SMAD2 and has a significant enhancer activity for PAH. However, the derived T allele (the major allele in East Asians) of rs10778203 decreases the binding activity of transcription factors and enhancer activity. Meanwhile, the derived T allele of rs10778203 shows a weaker ultraviolet radiation response in A375 cells and zebrafish embryos. Furthermore, rs10778203 decreases melanin production in transgenic zebrafish embryos after ultraviolet B (UVB) treatment. Collectively, PAH is a potential pigmentation gene that regulates skin tanning ability. Natural selection has enriched the adaptive allele, resulting in weakened tanning ability in East Asians, suggesting a unique genetic mechanism for evolutionary skin lightening in East Asians.

Mapping and annotating genomic loci to prioritize genes and implicate distinct polygenic adaptations for skin color

Evidence for adaptation of human skin color to regional ultraviolet radiation suggests shared and distinct genetic variants across populations. However, skin color evolution and genetics in East Asians are understudied. We quantified skin color in 48,433 East Asians using image analysis and identified associated genetic variants and potential causal genes for skin color as well as their polygenic interplay with sun exposure. This genome-wide association study (GWAS) identified 12 known and 11 previously unreported loci and SNP-based heritability was 23-24%. Potential causal genes were determined through the identification of nonsynonymous variants, colocalization with gene expression in skin tissues, and expression levels in melanocytes. Genomic loci associated with pigmentation in East Asians substantially diverged from European populations, and we detected signatures of polygenic adaptation. This large GWAS for objectively quantified skin color in an East Asian population improves understanding of the genetic architecture and polygenic adaptation of skin color and prioritizes potential causal genes.

Giant pandas in captivity undergo short-term adaptation in nerve-related pathways

BACKGROUND

Behaviors in captive animals, including changes in appetite, activity level, and social interaction, are often seen as adaptive responses. However, these behaviors may become progressively maladaptive, leading to stress, anxiety, depression, and other negative reactions in animals.

RESULTS

In this study, we investigated the whole-genome sequencing data of 39 giant panda individuals, including 11 in captivity and 28 in the wild. To eliminate the mountain range effect and focus on the factor of captivity only, we first performed a principal component analysis. We then enumerated the 21,474,180 combinations of wild giant pandas (11 chosen from 28) and calculated their distances from the 11 captive individuals. The 11 wild individuals with the closest distances were used for the subsequent analysis. The linkage disequilibrium (LD) patterns demonstrated that the population was almost eliminated. We identified 505 robust selected genomic regions harboring at least one SNP, and the absolute frequency difference was greater than 0.6 between the two populations. GO analysis revealed that genes in these regions were mainly involved in nerve-related pathways. Furthermore, we identified 22 GO terms for which the selection strength significantly differed between the two populations, and there were 10 nerve-related pathways among them. Genes in the differentially abundant regions were involved in nerve-related pathways, indicating that giant pandas in captivity underwent minor genomic selection. Additionally, we investigated the relationship between genetic variation and chromatin conformation structures. We found that nucleotide diversity (θπ) in the captive population was correlated with chromatin conformation structures, which included A/B compartments, topologically associated domains (TADs) and TAD-cliques. For each GO term, we then compared the expression level of genes regulated by the above four factors (AB index, TAD intactness, TAD clique and PEI) with the corresponding genomic background. The retained 10 GO terms were all coordinately regulated by the four factors, and three of them were associated with nerve-related pathways.

CONCLUSIONS

This study revealed that giant pandas in captivity undergo short-term adaptation in nerve-related pathways. Furthermore, it provides new insights into the molecular mechanism of gene expression regulation under short-term adaptation to environmental change.

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.

Genetic landscape of forensic DNA phenotyping markers among Mediterranean populations

Forensic DNA Phenotyping can reveal the appearance of an unknown individual by predicting the External Visible Characteristics (EVC) from DNA obtained at the crime scene. Our aim is to characterize the genetic landscape of Human identification markers responsible for EVC among Mediterranean populations compared to other worldwide groups. We conducted an exhaustive search for genes involved in EVC variation. Then, variants located on these genes were extracted from public genotypic data of Mediterranean, American, African and East Asiatic populations. The genetic landscape of these Human identification markers, their allelic distribution and admixture analyses, were determined using plink, R and ADMIXTURE softwares. Our results showed that the Mediterranean populations appear close to the Mexican populations and distinguished from sub Saharan African populations living in the USA and from East Asiatic populations. We highlighted a total of 103454 common variants shared between the studied populations and among them, 25 common variants associated with EVC. Interestingly, genotype frequencies results showed that the rs17646946, rs13016869, rs977588, rs1805008 and rs2240751 variants located respectively in the TCHH, PRKCE, OCA2, MC1R and MFSD12 genes are significantly different between the Mediterranean and Asiatic populations. The genotype frequencies of the variants rs977589 and rs7179994 located in the OCA2 gene, and of rs12913832 and rs2240751 located respectively in HERC2 and MFSD12 genes are significantly different between the Mediterranean and American populations. Our work generates a large number of EVC variants that could be a valuable resource for future studies in the forensic field.

Effects of facial skin pigmentation on social judgments in a Mexican population

People quickly and involuntarily form impressions of others based on their facial physical attributes, which can modulate critical social interactions. Skin pigmentation is one of the most variable and conspicuous facial traits among human populations. Empirical evidence suggests that these variations reflect ancestral ecological selective pressures balancing cutaneous vitamin D synthesis with the protection of the dermis from ultraviolet radiation. Nevertheless, skin pigmentation may currently be subject to additional selective pressures. For instance, the colonial era in Central and South America developed a highly stratified society based on ethnic origins, and light skin pigmentation became associated with higher social status and deference. This association could have originated through historical social learning that promoted favorable social perceptions towards individuals with lighter skin color and unfavorable perceptions towards individuals with darker skin color, which could still be present in the perception of current populations. Facial skin pigmentation is also sexually dimorphic, with males tending to exhibit darker skin than females, a difference that could be driven by sexual selection. To explore whether social learning and sexual selection represent additional selective pressures on skin pigmentation, we tested how this facial trait influences fundamental social perceptions in a Mexican population (N = 700, 489 female). We sampled facial images of eight European American males with natural lighter facial skin and eight males from an indigenous pre-Columbian community from Mexico, the Me'Phaa, with natural darker facial skin. We produced stimuli from these images by varying the skin pigmentation while preserving the facial shape. Stimuli were rated on attractiveness, trustworthiness, perceived health, dominance, aggressiveness, and femininity/masculinity. We found that the natural light-skinned faces were perceived as more attractive, trustworthy, and healthy but less dominant than the natural dark faces. Furthermore, by varying the facial skin color in these original groups, we altered the perceptions of them, mainly their attractiveness. These results partially support the hypothesis that dark facial skin color may help males compete for mates. Also, the results strongly support the view that lighter facial skin color became associated with social benefits through social learning in this Mexican population. Our findings, when viewed through the lens of cultural evolution, align with previous research in social psychology and anthropology. They hold the potential to offer a comprehensive understanding of the origin of this social phenomenon of cultural transmission, which currently plays a role in the formation of racial attitudes, stereotyping, and racial inequality in Mexican and other Latin American populations.

Neanderthal introgression in SCN9A impacts mechanical pain sensitivity

The Nav1.7 voltage-gated sodium channel plays a key role in nociception. Three functional variants in the SCN9A gene (encoding M932L, V991L, and D1908G in Nav1.7), have recently been identified as stemming from Neanderthal introgression and to associate with pain symptomatology in UK BioBank data. In 1000 genomes data, these variants are absent in Europeans but common in Latin Americans. Analysing high-density genotype data from 7594 Latin Americans, we characterized Neanderthal introgression in SCN9A. We find that tracts of introgression occur on a Native American genomic background, have an average length of ~123 kb and overlap the M932L, V991L, and D1908G coding positions. Furthermore, we measured experimentally six pain thresholds in 1623 healthy Colombians. We found that Neanderthal ancestry in SCN9A is significantly associated with a lower mechanical pain threshold after sensitization with mustard oil and evidence of additivity of effects across Nav1.7 variants. Our findings support the reported association of Neanderthal Nav1.7 variants with clinical pain, define a specific sensory modality affected by archaic introgression in SCN9A and are consistent with independent effects of the Neanderthal variants on Nav1.7 function.

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.

The role of Beringia in human adaptation to Arctic conditions based on results of genomic studies of modern and ancient populations

The results of studies in Quaternary geology, archeology, paleoanthropology and human genetics demonstrate that the ancestors of Native Americans arrived in mid-latitude North America mainly along the Pacific Northwest Coast, but had previously inhabited the Arctic and during the last glacial maximum were in a refugium in Beringia, a land bridge connecting Eurasia and North America. The gene pool of Native Americans is represented by unique haplogroups of mitochondrial DNA and the Y chromosome, the evolutionary age of which ranges from 13 to 22 thousand years. The results of a paleogenomic analysis also show that during the last glacial maximum Beringia was populated by human groups that had arisen as a result of interaction between the most ancient Upper Paleolithic populations of Northern Eurasia and newcomer groups from East Asia. Approximately 20 thousand years ago the Beringian populations began to form, and the duration of their existence in relative isolation is estimated at about 5 thousand years. Thus, the adaptation of the Beringians to the Arctic conditions could have taken several millennia. The adaptation of Amerindian ancestors to high latitudes and cold climates is supported by genomic data showing that adaptive genetic variants in Native Americans are associated with various metabolic pathways: melanin production processes in the skin, hair and eyes, the functioning of the cardiovascular system, energy metabolism and immune response characteristics. Meanwhile, the analysis of the existing hypotheses about the selection of some genetic variants in the Beringian ancestors of the Amerindians in connection with adaptation to the Arctic conditions (for example, in the FADS, ACTN3, EDAR genes) shows the ambiguity of the testing results, which may be due to the loss of some traces of the "Beringian" adaptation in the gene pools of modern Native Americans. The most optimal strategy for further research seems to be the search for adaptive variant.

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.

Genetic Variants in Telomerase Reverse Transcriptase Contribute to Solar Lentigines

Solar lentigines (SLs) are a hallmark of human skin aging. They result from chronic exposure to sunlight and other environmental stressors. Recent studies also imply genetic factors, but findings are partially conflicting and lack of replication. Through a multi-trait based analysis strategy, we discovered that genetic variants in telomerase reverse transcriptase were significantly associated with non-facial SL in two East Asian (Taizhou longitudinal cohort, n = 2,964 and National Survey of Physical Traits, n = 2,954) and one Caucasian population (SALIA, n = 462), top SNP rs2853672 (P-value for Taizhou longitudinal cohort = 1.32 × 10‒28 and P-value for National Survey of Physical Traits = 3.66 × 10‒17 and P-value for SALIA = 0.0007 and Pmeta = 4.93 × 10‒44). The same variants were nominally associated with facial SL but not with other skin aging or skin pigmentation traits. The SL-enhanced allele/haplotype upregulated the transcription of the telomerase reverse transcriptase gene. Of note, well-known telomerase reverse transcriptase‒related aging markers such as leukocyte telomere length and intrinsic epigenetic age acceleration were not associated with SL. Our results indicate a previously unrecognized role of telomerase reverse transcriptase in skin aging‒related lentigines formation.

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.

Automatic landmarking identifies new loci associated with face morphology and implicates Neanderthal introgression in human nasal shape

We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value <5 × 10-8) at 42 genome regions, nine of which have been previously reported. In follow-up analyses, 26 of the 33 novel regions replicate in East Asians, Europeans, or Africans, and one mouse homologous region influences craniofacial morphology in mice. The novel region in 1q32.3 shows introgression from Neanderthals and we find that the introgressed tract increases nasal height (consistent with the differentiation between Neanderthals and modern humans). Novel regions include candidate genes and genome regulatory elements previously implicated in craniofacial development, and show preferential transcription in cranial neural crest cells. The automated approach used here should simplify the collection of large study samples from across the world, facilitating a cosmopolitan characterization of the genetics of facial features.

DNA methylation abnormalities induced by advanced maternal age in villi prime a high-risk state for spontaneous abortion

BACKGROUND

Advanced maternal age (AMA) has increased in many high-income countries in recent decades. AMA is generally associated with a higher risk of various pregnancy complications, and the underlying molecular mechanisms are largely unknown. In the current study, we profiled the DNA methylome of 24 human chorionic villi samples (CVSs) from early pregnancies in AMA and young maternal age (YMA), 11 CVSs from early spontaneous abortion (SA) cases using reduced representation bisulfite sequencing (RRBS), and the transcriptome of 10 CVSs from AMA and YMA pregnancies with mRNA sequencing(mRNA-seq). Single-cell villous transcriptional atlas presented expression patterns of targeted AMA-/SA-related genes. Trophoblast cellular impairment was investigated through the knockdown of GNE expression in HTR8-S/Vneo cells.

RESULTS

AMA-induced local DNA methylation changes, defined as AMA-related differentially methylated regions (DMRs), may be derived from the abnormal expression of genes involved in DNA demethylation, such as GADD45B. These DNA methylation changes were significantly enriched in the processes involved in NOTCH signaling and extracellular matrix organization and were reflected in the transcriptional alterations in the corresponding biological processes and specific genes. Furthermore, the DNA methylation level of special AMA-related DMRs not only significantly changed in AMA but also showed more excessive defects in CVS from spontaneous abortion (SA), including four AMA-related DMRs whose nearby genes overlapped with AMA-related differentially expressed genes (DEGs) (CDK11A, C19orf71, COL5A1, and GNE). The decreased DNA methylation level of DMR near GNE was positively correlated with the downregulated expression of GNE in AMA. Single-cell atlas further revealed comparatively high expression of GNE in the trophoblast lineage, and knockdown of GNE in HTR8-S/Vneo cells significantly impaired cellular proliferation and migration.

CONCLUSION

Our study provides valuable resources for investigating AMA-induced epigenetic abnormalities and provides new insights for explaining the increased risks of pregnancy complications in AMA pregnancies.

Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

Piebaldism and chromatophore development in reptiles are linked to the tfec gene

Reptiles display great diversity in color and pattern, yet much of what we know about vertebrate coloration comes from classic model species such as the mouse and zebrafish.1,2,3,4 Captive-bred ball pythons (Python regius) exhibit a remarkable degree of color and pattern variation. Despite the wide range of Mendelian color phenotypes available in the pet trade, ball pythons remain an overlooked species in pigmentation research. Here, we investigate the genetic basis of the recessive piebald phenotype, a pattern defect characterized by patches of unpigmented skin (leucoderma). We performed whole-genome sequencing and used a case-control approach to discover a nonsense mutation in the gene encoding the transcription factor tfec, implicating this gene in the leucodermic patches in ball pythons. We functionally validated tfec in a lizard model (Anolis sagrei) using the gene editing CRISPR/Cas9 system and TEM imaging of skin. Our findings show that reading frame mutations in tfec affect coloration and lead to a loss of iridophores in Anolis, indicating that tfec is required for chromatophore development. This study highlights the value of captive-bred ball pythons as a model species for accelerating discoveries on the genetic basis of vertebrate coloration.

LncRNA evolution and DNA methylation variation participate in photosynthesis pathways of distinct lineages of Populus

During the independent process of evolution in plants, photosynthesis appears to have been under convergent evolution to adapt to specific selection pressure in their geographical regions. However, it is unclear how lncRNA regulation and DNA methylation are involved in the phenotypic convergence in distinct lineages. Here, we present a large-scale comparative study of lncRNA transcription profile and whole-genome bisulfite sequencing (WGBS) data in two unrelated Populus species, selected from three relatively overlapping geographical regions. The results indicated that 39.75% lncRNAs of Populus tomentosa were shown to have homologous sequences in the 46.99% lncRNA of Populus simonii. Evolutionary analysis revealed that lncRNAs showed a rapid gain rate in the Populus lineage. Furthermore, co-expression networks in two Populus species identified eight lncRNAs that have the potential to simultaneously cis- or trans-regulate eight photosynthetic-related genes. These photosynthetic lncRNAs and genes were predominantly expressed in accessions from the southern region, indicating a conserved spatial expression in photosynthetic pathways in Populus. We also detected that most lncRNA targeted photosynthetic genes hypomethylated in promoter regions of Southern accessions compared with Northern accessions. Geographical DMRs correlated with genetic SNP variations in photosynthetic genes among Populus from the three geographic regions, indicating that DNA methylation coordinated with lncRNAs in convergent evolution of photosynthesis in Populus. Our results shed light on the evolutionary forces acting on patterns of lncRNA and DNA methylation, and provide a better understanding of the genetic and epigenetic mechanism in photosynthetic convergence evolution.

Electrically conductive carbon-based (bio)-nanomaterials for cardiac tissue engineering

A proper self-regenerating capability is lacking in human cardiac tissue which along with the alarming rate of deaths associated with cardiovascular disorders makes tissue engineering critical. Novel approaches are now being investigated in order to speedily overcome the challenges in this path. Tissue engineering has been revolutionized by the advent of nanomaterials, and later by the application of carbon-based nanomaterials because of their exceptional variable functionality, conductivity, and mechanical properties. Electrically conductive biomaterials used as cell bearers provide the tissue with an appropriate microenvironment for the specific seeded cells as substrates for the sake of protecting cells in biological media against attacking mechanisms. Nevertheless, their advantages and shortcoming in view of cellular behavior, toxicity, and targeted delivery depend on the tissue in which they are implanted or being used as a scaffold. This review seeks to address, summarize, classify, conceptualize, and discuss the use of carbon-based nanoparticles in cardiac tissue engineering emphasizing their conductivity. We considered electrical conductivity as a key affecting the regeneration of cells. Correspondingly, we reviewed conductive polymers used in tissue engineering and specifically in cardiac repair as key biomaterials with high efficiency. We comprehensively classified and discussed the advantages of using conductive biomaterials in cardiac tissue engineering. An overall review of the open literature on electroactive substrates including carbon-based biomaterials over the last decade was provided, tabulated, and thoroughly discussed. The most commonly used conductive substrates comprising graphene, graphene oxide, carbon nanotubes, and carbon nanofibers in cardiac repair were studied.

Association of TYR SNP rs1042602 with Melanoma Risk and Prognosis

Cutaneous melanoma is the most aggressive of skin tumors. In order to discover new biomarkers that could help us improve prognostic prediction in melanoma patients, we have searched for germline DNA variants associated with melanoma progression. Thus, after exome sequencing of a set of melanoma patients and healthy control individuals, we identified rs1042602, an SNP within TYR, as a good candidate. After genotyping rs1042602 in 1025 patients and 773 healthy donors, we found that the rs1042602-A allele was significantly associated with susceptibility to melanoma (CATT test: p = 0.0035). Interestingly, we also observed significant differences between patients with good and bad prognosis (5 years of follow-up) (n = 664) (CATT test for all samples p = 0.0384 and for men alone p = 0.0054). Disease-free-survival (DFS) analyses also showed that patients with the A allele had shorter DFS periods. In men, the association remained significant even in a multivariate Cox Proportional-hazards model, which was adjusted for age at diagnosis, Breslow thickness, ulceration and melanoma subtype (HR 0.4; 95% confidence interval (CI) 0.20-0.83; p = 0.0139). Based on our results, we propose that rs1042602-A is a risk allele for melanoma, which also seems to be responsible for a poorer prognosis of the disease, particularly in men.

The potential impact of melanosomal pH and metabolism on melanoma

Melanin is synthesized in melanocytes and is transferred into keratinocytes to block the effects of ultraviolet (UV) radiation and is important for preventing skin cancers including melanoma. However, it is known that after melanomagenesis and melanoma invasion or metastases, melanin synthesis still occurs. Since melanoma cells are no longer involved in the sun tanning process, it is unclear why melanocytes would maintain melanin synthesis after melanomagenesis has occurred. Aside from blocking UV-induced DNA mutation, melanin may provide other metabolic functions that could benefit melanoma. In addition, studies have suggested that there may be a selective advantage to melanin synthesis in melanoma; however, mechanisms regulating melanin synthesis outside the epidermis or hair follicle is unknown. We will discuss how melanosomal pH controls melanin synthesis in melanocytes and how melanosomal pH control of melanin synthesis might function in melanoma. We will also discuss potential reasons why melanin synthesis might be beneficial for melanoma cellular metabolism and provide a rationale for why melanin synthesis is not limited to benign melanocytes.

Candidate pigmentation genes related to feather color variation in an indigenous chicken breed revealed by whole genome data

Chicken plumage color is an inheritable phenotype that was naturally and artificially selected for during domestication. The Baicheng You chicken is an indigenous Chinese chicken breed presenting three main feather colors, lavender, black, and yellow plumages. To explore the genetic mechanisms underlying the pigmentation in Baicheng You chickens, we re-sequenced the whole genome of Baicheng You chicken with the three plumage colors. By analyzing the divergent regions of the genome among the chickens with different feather colors, we identified some candidate genomic regions associated with the feather colors in Baicheng You chickens. We found that EGR1, MLPH, RAB17, SOX5, and GRM5 genes were the potential genes for black, lavender, and yellow feathers. MLPH, GRM5, and SOX5 genes have been found to be related to plumage colors in birds. Our results showed that EGR1 is a most plausible candidate gene for black plumage, RAB17, MLPH, and SOX5 for lavender plumage, and GRM5 for yellow plumage in Baicheng You chicken.

Genetic adaptation of skin pigmentation in highland Tibetans

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

Applications of massively parallel sequencing in forensic genetics

Massively parallel sequencing, also referred to as next-generation sequencing, has positively changed DNA analysis, allowing further advances in genetics. Its capability of dealing with low quantity/damaged samples makes it an interesting instrument for forensics. The main advantage of MPS is the possibility of analyzing simultaneously thousands of genetic markers, generating high-resolution data. Its detailed sequence information allowed the discovery of variations in core forensic short tandem repeat loci, as well as the identification of previous unknown polymorphisms. Furthermore, different types of markers can be sequenced in a single run, enabling the emergence of DIP-STRs, SNP-STR haplotypes, and microhaplotypes, which can be very useful in mixture deconvolution cases. In addition, the multiplex analysis of different single nucleotide polymorphisms can provide valuable information about identity, biogeographic ancestry, paternity, or phenotype. DNA methylation patterns, mitochondrial DNA, mRNA, and microRNA profiling can also be analyzed for different purposes, such as age inference, maternal lineage analysis, body-fluid identification, and monozygotic twin discrimination. MPS technology also empowers the study of metagenomics, which analyzes genetic material from a microbial community to obtain information about individual identification, post-mortem interval estimation, geolocation inference, and substrate analysis. This review aims to discuss the main applications of MPS in forensic genetics.

Sodium-calcium exchanger-3 regulates pain "wind-up": From human psychophysics to spinal mechanisms

Repeated application of noxious stimuli leads to a progressively increased pain perception; this temporal summation is enhanced in and predictive of clinical pain disorders. Its electrophysiological correlate is "wind-up," in which dorsal horn spinal neurons increase their response to repeated nociceptor stimulation. To understand the genetic basis of temporal summation, we undertook a GWAS of wind-up in healthy human volunteers and found significant association with SLC8A3 encoding sodium-calcium exchanger type 3 (NCX3). NCX3 was expressed in mouse dorsal horn neurons, and mice lacking NCX3 showed normal, acute pain but hypersensitivity to the second phase of the formalin test and chronic constriction injury. Dorsal horn neurons lacking NCX3 showed increased intracellular calcium following repetitive stimulation, slowed calcium clearance, and increased wind-up. Moreover, virally mediated enhanced spinal expression of NCX3 reduced central sensitization. Our study highlights Ca2+ efflux as a pathway underlying temporal summation and persistent pain, which may be amenable to therapeutic targeting.

Novel genetic associations with five aesthetic facial traits: A genome-wide association study in the Chinese population

Background: The aesthetic facial traits are closely related to life quality and strongly influenced by genetic factors, but the genetic predispositions in the Chinese population remain poorly understood.

Methods: A genome-wide association studies (GWAS) and subsequent validations were performed in 26,806 Chinese on five facial traits: widow’s peak, unibrow, double eyelid, earlobe attachment, and freckles. Functional annotation was performed based on the expression quantitative trait loci (eQTL) variants, genome-wide polygenic scores (GPSs) were developed to represent the combined polygenic effects, and single nucleotide polymorphism (SNP) heritability was presented to evaluate the contributions of the variants.

Results: In total, 21 genetic associations were identified, of which ten were novel: GMDS-AS1 (rs4959669, p = 1.29 × 10⁻⁴⁹) and SPRED2 (rs13423753, p = 2.99 × 10⁻¹⁴) for widow’s peak, a previously unreported trait; FARSB (rs36015125, p = 1.96 × 10⁻²¹) for unibrow; KIF26B (rs7549180, p = 2.41 × 10⁻¹⁵), CASC2 (rs79852633, p = 4.78 × 10⁻¹¹⁾, RPGRIP1L (rs6499632, p = 9.15 × 10⁻¹¹), and PAX1 (rs147581439, p = 3.07 × 10⁻⁸) for double eyelid; ZFHX3 (rs74030209, p = 9.77 × 10⁻¹⁴) and LINC01107 (rs10211400, p = 6.25 × 10⁻¹⁰) for earlobe attachment; and SPATA33 (rs35415928, p = 1.08 × 10⁻⁸) for freckles. Functionally, seven identified SNPs tag the missense variants and six may function as eQTLs. The combined polygenic effect of the associations was represented by GPSs and contributions of the variants were evaluated using SNP heritability.

Conclusion: These identifications may facilitate a better understanding of the genetic basis of features in the Chinese population and hopefully inspire further genetic research on facial development.

Population differentiation of polygenic score predictions under stabilizing selection

Given the many small-effect loci uncovered by genome-wide association studies (GWAS), polygenic scores have become central to genomic medicine, and have found application in diverse settings including evolutionary studies of adaptation. Despite their promise, polygenic scores have been found to suffer from limited portability across human populations. This at first seems in conflict with the observation that most common genetic variation is shared among populations. We investigate one potential cause of this discrepancy: stabilizing selection on complex traits. Counterintuitively, while stabilizing selection constrains phenotypic evolution, it accelerates the loss and fixation of alleles underlying trait variation within populations (GWAS loci). Thus even when populations share an optimum phenotype, stabilizing selection erodes the variance contributed by their shared GWAS loci, such that predictions from GWAS in one population explain less of the phenotypic variation in another. We develop theory to quantify how stabilizing selection is expected to reduce the prediction accuracy of polygenic scores in populations not represented in GWAS samples. In addition, we find that polygenic scores can substantially overstate average genetic differences of phenotypes among populations. We emphasize stabilizing selection around a common optimum as a useful null model to connect patterns of allele frequency and polygenic score differentiation. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

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.

Genome-wide association studies in Japanese women identified genetic loci associated with wrinkles and sagging

Wrinkles and sagging are caused by various factors, such as ultraviolet rays; however, recent findings demonstrated that some individuals are genetically predisposed to these phenotypes of skin aging. The contribution of single nucleotide polymorphisms (SNPs) to the development of wrinkles and sagging has been demonstrated in genome-wide association studies (GWAS). However, these findings were mainly obtained from European and Chinese populations. Limited information is currently available on the involvement of SNPs in the development of wrinkles and sagging in a Japanese population. Therefore, we herein performed GWAS on wrinkles at the outer corners of the eyes and nasolabial folds in 1041 Japanese women. The results obtained revealed that 5 SNPs (19p13.2: rs2303098 (p = 3.39×10^-8 ), rs56391955 (p = 3.39×10^-8 ), rs67560822 (p = 3.50×10^-8 ), rs889126 (p = 3.78×10^-8 ), rs57490083 (p = 3.99×10^-8 )) located within the COL5A3 gene associated with wrinkles at the outer corners of the eyes. Regarding nasolabial folds, 8q24.11 (rs4876369; p = 1.05×10^-7 , rs6980503; p = 1.25×10^-7 , rs61027543; p = 1.25×10^-7 , rs16889363; p = 1.38×10^-7 ) was suggested to be associated with RAD21 gene expression. These SNPs have not been reported in other populations, and were first found in Japanese women population. These SNPs may be used as markers to examine the genetic predisposition of individuals to wrinkles and sagging.

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

Genetic diversity of the melanocortin-1 receptor in an admixed population of Rio de Janeiro: Structural and functional impacts of Cys35Tyr variant

The melanocortin-1 receptor (MC1R) is one of the key proteins involved in the regulation of melanin production and several polymorphisms have been associated with different phenotypes of skin and hair color in human and nonhuman species. Most of the knowledge is centered on more homogeneous populations and studies involving an admixed group of people should be encouraged due to the great importance of understanding the human color variation. This work evaluates the MC1R diversity and the possible impacts of MC1R variants in an admixed sample population of Rio de Janeiro, Brazil, which is a product of Native American, African, and European miscegenation. Sequencing of complete coding region and part of the 3´UTR of MC1R gene identified 31 variants including one insertion and three novel synonymous substitutions in sample population grouped according to skin, hair and eye pigmentation levels. In nonmetric multidimensional scaling analysis (NMDS), three main clusters were identified, in which the Brazilian dark skin group remained in the African cluster whereas the intermediate and the light skin color phenotype in the European one. None gathered with Asians since their immigration to Brazil was a recent event. In silico analyses demonstrated that Cys35Tyr, Ile155Thr and Pro256Ser, found in our population, have a negative effect on receptor function probably due to changes on the receptor structure. Notably, Cys35Tyr mutation could potentially impair agonist binding. Altogether, this work contributes to the understanding of the genetic background of color variation on an admixed population and gives insights into the damaging effects of MC1R variants.

From musk to body odor: Decoding olfaction through genetic variation

The olfactory system combines input from multiple receptor types to represent odor information, but there are few explicit examples relating olfactory receptor (OR) activity patterns to odor perception. To uncover these relationships, we performed genome-wide scans on odor-perception phenotypes for ten odors in 1000 Han Chinese and validated results for six of these odors in an ethnically diverse population (n = 364). In both populations, consistent with previous studies, we replicated three previously reported associations (β-ionone/OR5A1, androstenone/OR7D4, cis-3-hexen-1-ol/OR2J3 LD-band), but not for odors containing aldehydes, suggesting that olfactory phenotype/genotype studies are robust across populations. Two novel associations between an OR and odor perception contribute to our understanding of olfactory coding. First, we found a SNP in OR51B2 that associated with trans-3-methyl-2-hexenoic acid, a key component of human underarm odor. Second, we found two linked SNPs associated with the musk Galaxolide in a novel musk receptor, OR4D6, which is also the first human OR shown to drive specific anosmia to a musk compound. We noticed that SNPs detected for odor intensity were enriched with amino acid substitutions, implying functional changes of odor receptors. Furthermore, we also found that the derived alleles of the SNPs tend to be associated with reduced odor intensity, supporting the hypothesis that the primate olfactory gene repertoire has degenerated over time. This study provides information about coding for human body odor, and gives us insight into broader mechanisms of olfactory coding, such as how differential OR activation can converge on a similar percept.

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.

Canine coat pigmentation genetics: a review

Our understanding of canine coat colour genetics and the associated health implications is developing rapidly. To date, there are 15 genes with known roles in canine coat colour phenotypes. Many coat phenotypes result from complex and/or epistatic genetic interactions among variants within and between loci, some of which remain unidentified. Some genes involved in canine pigmentation have been linked to aural, visual and neurological impairments. Consequently, coat pigmentation in the domestic dog retains considerable ethical and economic interest. In this paper we discuss coat colour phenotypes in the domestic dog, the genes and variants responsible for these phenotypes and any proven coat colour-associated health effects.

A large Canadian cohort provides insights into the genetic architecture of human hair colour

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

Melanosome Biogenesis in the Pigmentation of Mammalian Skin

Melanins, the main pigments of the skin and hair in mammals, are synthesized within membrane-bound organelles of melanocytes called melanosomes. Melanosome structure and function are determined by a cohort of resident transmembrane proteins, many of which are expressed only in pigment cells and localize specifically to melanosomes. Defects in the genes that encode melanosome-specific proteins or components of the machinery required for their transport in and out of melanosomes underlie various forms of ocular or oculocutaneous albinism, characterized by hypopigmentation of the hair, skin, and eyes and by visual impairment. We review major components of melanosomes, including the enzymes that catalyze steps in melanin synthesis from tyrosine precursors, solute transporters that allow these enzymes to function, and structural proteins that underlie melanosome shape and melanin deposition. We then review the molecular mechanisms by which these components are biosynthetically delivered to newly forming melanosomes-many of which are shared by other cell types that generate cell type-specific lysosome-related organelles. We also highlight unanswered questions that need to be addressed by future investigation.

GWAS Identifies Multiple Genetic Loci for Skin Color in Korean Women

Human skin color is largely determined by genetic factors. Recent GWASs have reported several genetic variants associated with skin color, mostly in European and African populations. In this study, we performed GWAS in 17,019 Korean women to identify genetic variants associated with facial skin color, quantitatively measured as CIELAB color index. We identified variants in three, one, and six genomic loci associated with facial skin color index L∗, a∗, and b∗ values, respectively, and replicated the associations (combined analysis P-value < 5.0 × 10^-8). The significant loci included variants in known genes (OCA2 rs74653330, BNC2 rs16935073, rs72620727 near KITLG, and SLC6A17 rs6689641) and to our knowledge previously unreported genes (SCARB1 rs10846744, SYN2 rs12629034, and LINC00486 rs6543678). This is GWAS to elucidate genetic variants of facial skin color in a Korean female population. Further functional characterizations of the investigated genes are warranted to elucidate their contribution to skin pigmentation-related traits.

Genome-wide association study of occupational attainment as a proxy for cognitive reserve

Occupational attainment, which represents middle-age cognitive activities, is a known proxy marker of cognitive reserve for Alzheimer's disease. Previous genome-wide association studies (GWAS) have identified numerous genetic variants and revealed the genetic architecture of educational attainment, another marker of cognitive reserve. However, the genetic architecture and heritability for occupational attainment remain elusive. We performed a large-scale GWAS of occupational attainment with 248,847 European individuals from the UK Biobank using the proportional odds logistic mixed model method. In this analysis, we defined occupational attainment using the classified job levels formulated in the UK Standard Occupational Classification system considering the individual professional skill and academic level. We identified 30 significant loci (P < 5 × 10-8); 12 were novel variants, unassociated with other traits. Among them, four lead variants were associated with genes expressed in brain tissues by expression quantitative trait loci mapping from 10 brain regions: rs13002946, rs3741368, rs11654986, and rs1627527. The single nucleotide polymorphism (SNP)-based heritability was estimated to be 8.5% (s.e. = 0.004) and partitioned heritability was enriched in the central nervous system and brain tissues. Genetic correlation analysis showed shared genetic backgrounds between occupational attainment and multiple traits, including education, intelligence, leisure activities, life satisfaction, and neuropsychiatric disorders. In two-sample Mendelian randomization (MR) analysis, we demonstrated that high occupation levels were associated with reduced risk for Alzheimer's disease (OR = 0.78, 95% CI = 0.65-0.92 in inverse variance weighted (IVW) method; OR = 0.73, 95% CI = 0.57-0.92 in the weighted median (WM) method). This causal relationship between occupational attainment and Alzheimer's disease was robust in additional sensitivity analysis that excluded potentially pleiotropic SNPs (OR = 0.72, 95% CI = 0.57-0.91 in the IVW method; OR = 0.72, 95% CI = 0.53-0.97 in the WM method). Multivariable MR confirmed that occupational attainment had an independent effect on the risk for Alzheimer's disease even after taking educational attainment into account (OR = 0.72, 95% CI = 0.54-0.95 in the IVW method; OR = 0.68, 95% CI = 0.48-0.97 in the WM method). Overall, our analyses provide insights into the genetic architecture of occupational attainment and demonstrate that occupational attainment is a potential causal protective factor for Alzheimer's disease as a proxy marker of cognitive reserve.