Genome-wide association studies identify several new loci associated with pigmentation traits and skin cancer risk in European Americans

Design and quality control of large-scale two-sample Mendelian randomization studies

Background

Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors.

Methods

We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set.

Results

We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls.

Conclusions

In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

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.

Risk of Skin Cancer and Actinic Keratosis in Patients with Rosacea: A Nationwide Population-based Cohort Study

The association between rosacea and skin cancer remains inconclusive, with conflicting reports. The aim of this nationwide population-based cohort study was to determine the risk of skin cancer in patients with rosacea. A rosacea cohort (n = 11,420) was formulated and evaluated from 2010 to 2019. The incidence rate ratios of actinic keratosis, cutaneous melanoma, keratinocyte carcinoma and gastric, colorectal, and liver cancer were analysed in comparison with a matched control group, and multivariable stratified Cox proportional hazards model analysis was performed. The risk of actinic keratosis and keratinocyte carcinoma was increased in the rosacea group compared with the control group, with adjusted hazard ratios of 6.05 (95% confidence interval 3.63–10.09) and 2.66 (1.53–4.61), respectively. The risk of cutaneous melanoma and gastric, colorectal and liver cancer was not increased, with adjusted hazard ratios of 1.69 (0.25–11.37), 0.81 (0.59–1.10), 0.91 (0.69–1.18) and 1.32 (0.89–1.95), respectively. These results reveal an increased risk of actinic keratosis and keratinocyte carcinoma in patients with rosacea.

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

Purpose

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

Methods

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

Results

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

Conclusion

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

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.

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

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.

Epigenetic regulation during melanocyte development and homeostasis

Melanocytes originate in the neural crest as precursor cells which then migrate and proliferate to reach their destination where they differentiate into pigment-producing cells. Melanocytes not only determine the colour of hair, skin and eyes but also protect against the harmful effects of UV irradiation. The establishment of the melanocyte lineage is regulated by a defined set of transcription factors and signalling pathways that direct the specific gene expression programmes underpinning melanoblast specification, survival, migration, proliferation and differentiation. In addition, epigenetic modifiers and replacement histones play key roles in regulating gene expression and its timing during the different steps of this process. Here, we discuss the evidence for the role of epigenetic regulators in melanocyte development and function and how they interact with transcription factors and signalling pathways to establish and maintain this important cell lineage.

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

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

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.

The color of normal: How a Eurocentric focus erases pigmentation complexity

OBJECTIVES

Skin pigmentation is both a highly variable and highly visible human phenotypic trait. Investigations into the biology and origins of this variation have been the focus of research in the fields of dermatology, anthropology, and forensic science, among others. This manuscript explores how much of what we know about the biology, genetics, and evolutionary origins of pigmentation has been strongly influenced by investigations and applications that focus on lighter skin.

METHODS

I reviewed literature from the fields of dermatology, anthropology and evolutionary genetics, and forensic science to assess how perceptions of lighter skin as the "normal" state in humans can shape the ways that knowledge is gathered and applied in these fields.

RESULTS

This normalization of lighter skin has impacted common tools used in dermatology and shaped the framework of dermatological education. A strong Eurocentric bias has shaped our understanding of the genetic architecture of pigmentary traits, which influences the ways in we understand the evolutionary processes leading to modern pigmentation diversity. Finally, I discuss how these biases in pigmentation genetics work in combination with phenotypic systems that privilege predicting lighter pigmentation variation to impede accurate prediction of intermediate phenotypes, particularly in individuals with ancestry from multiple populations. This can lead to a disproportionate targeting of already over-policed populations with darker skin.

CONCLUSIONS

Potential changes to how we conceptualize clinical and basic pigmentation research may help to reduce existing health disparities and improve understanding of pigmentation genetic architecture and how this knowledge is applied in forensic contexts.

VWF maturation and release are controlled by 2 regulators of Weibel-Palade body biogenesis: exocyst and BLOC-2

von Willebrand factor (VWF) is an essential hemostatic protein that is synthesized in endothelial cells and stored in Weibel-Palade bodies (WPBs). Understanding the mechanisms underlying WPB biogenesis and exocytosis could enable therapeutic modulation of endogenous VWF, yet optimal targets for modulating VWF release have not been established. Because biogenesis of lysosomal related organelle-2 (BLOC-2) functions in the biogenesis of platelet dense granules and melanosomes, which like WPBs are lysosome-related organelles, we hypothesized that BLOC-2-dependent endolysosomal trafficking is essential for WPB biogenesis and sought to identify BLOC-2-interacting proteins. Depletion of BLOC-2 caused misdirection of cargo-carrying transport tubules from endosomes, resulting in immature WPBs that lack endosomal input. Immunoprecipitation of BLOC-2 identified the exocyst complex as a binding partner. Depletion of the exocyst complex phenocopied BLOC-2 depletion, resulting in immature WPBs. Furthermore, releasates of immature WPBs from either BLOC-2 or exocyst-depleted endothelial cells lacked high-molecular weight (HMW) forms of VWF, demonstrating the importance of BLOC-2/exocyst-mediated endosomal input during VWF maturation. However, BLOC-2 and exocyst showed very different effects on VWF release. Although BLOC-2 depletion impaired exocytosis, exocyst depletion augmented WPB exocytosis, indicating that it acts as a clamp. Exposure of endothelial cells to a small molecule inhibitor of exocyst, Endosidin2, reversibly augmented secretion of mature WPBs containing HMW forms of VWF. These studies show that, although BLOC-2 and exocyst cooperate in WPB formation, only exocyst serves to clamp WPB release. Exocyst function in VWF maturation and release are separable, a feature that can be exploited to enhance VWF release.

Genome-Wide Association Study Identifies Genetic Associations with Perceived Age

Failure of dermal protection or repair mechanisms might lead to visibly aged skin. The study aimed to identify genetic associations with perceived age. A genome-wide association study was undertaken in 423,992 adult participants of UK Biobank, using questionnaire data on perceived age and genetic data imputed to the Haplotype Reference Consortium imputation panel. The study identified 74 independently associated genetic loci, to our knowledge previously unreported (P < 5 × 10-8), which were enriched for cell signaling pathways, including the NEK6 and SMAD2 subnetworks. Common genetic variation was estimated to account for 14% of variation in perceived age, and the heritability of perceived age was partially shared with that of 75 other traits, including multiple traits representing adiposity, suggesting that perceived age may be a useful proxy trait in genetic association studies.

Germline and somatic albinism variants in amelanotic/hypomelanotic melanoma: Increased carriage of TYR and OCA2 variants

Amelanotic/hypomelanotic melanoma is a clinicopathologic subtype with absent or minimal melanin. This study assessed previously reported coding variants in albinism genes (TYR, OCA2, TYRP1, SLC45A2, SLC24A5, LRMDA) and common intronic, regulatory variants of OCA2 in individuals with amelanotic/hypomelanotic melanoma, pigmented melanoma cases and controls. Exome sequencing was available for 28 individuals with amelanotic/hypomelanotic melanoma and 303 individuals with pigmented melanoma, which were compared to whole exome data from 1144 Australian controls. Microarray genotyping was available for a further 17 amelanotic/hypomelanotic melanoma, 86 pigmented melanoma, 147 melanoma cases (pigmentation unknown) and 652 unaffected controls. Rare deleterious variants in TYR/OCA1 were more common in amelanotic/hypomelanotic melanoma cases than pigmented melanoma cases (set mixed model association tests P = 0.0088). The OCA2 hypomorphic allele p.V443I was more common in melanoma cases (1.8%) than controls (1.0%, X2 P = 0.02), and more so in amelanotic/hypomelanotic melanoma (4.4%, X2 P = 0.007). No amelanotic/hypomelanotic melanoma cases carried an eye and skin darkening haplotype of OCA2 (including rs7174027), present in 7.1% of pigmented melanoma cases (P = 0.0005) and 9.4% controls. Variants in TYR and OCA2 may play a role in amelanotic/hypomelanotic melanoma susceptibility. We suggest that somatic loss of function at these loci could contribute to the loss of tumor pigmentation, consistent with this we found a higher rate of somatic mutation in TYR/OCA2 in amelanotic/hypomelanotic melanoma vs pigmented melanoma samples (28.6% vs 3.0%; P = 0.021) from The Cancer Genome Atlas Skin Cutaneous Melanoma collection.

An Agouti-Signaling-Protein Mutation is Strongly Associated with Melanism in European Roe Deer (Capreolus capreolus)

Although the European roe deer (Capreolus capreolus) population of North-West Germany has a remarkable number of melanistic specimens between 10% and 25%, the underlying genetic mutation-causing melanism is still unknown. We used a gene targeting approach focusing on MC1R and ASIP as important genes of coat coloration. Overall, 1384 bp of MC1R and 2039 bp of ASIP were sequenced in 24 specimens and several SNPs were detected. But only the ASIP-SNP c.33G>T completely segregated both phenotypes leading to the amino acid substitution p.Leu11Phe. The SNP was further evaluated in additional 471 samples. Generally, all black specimens (n = 33) were homozygous TT, whereas chestnut individuals were either homozygote GG (n = 436) or heterozygote GT (n = 26). Considering the fact that all melanistic animals shared two mutated alleles of the strongly associated SNP, we concluded that melanism is inherited in a recessive mode in European roe deer.

Recognizable neonatal clinical features of aplasia cutis congenita

Background

Aplasia cutis congenita (ACC), classified in nine groups, is likely to be underreported, since milder isolated lesions in wellbeing newborns could often be undetected, and solitary lesions in the context of polymalformative syndromes could not always be reported. Regardless of form and cause, therapeutic options have in common the aim to restore the deficient mechanical and immunological cutaneous protection and to limit the risk of fluid leakage or rupture of the exposed organs. We aimed to review our institutional prevalence, comorbidities, treatment and outcome of newborns with ACC.

Methods

We conducted a retrospective study including all newborns affected by ACC and admitted at the University Mother-Child Department from October 2010 to October 2019. Anthropometric and clinical characteristics of ACC1 versus a non-isolated ACC group were analyzed.

Results

We encountered 37 newborns, 16 with ACC1 versus 21 with non-isolated ACC. The incidence rate of 0.1% in ACC1 was higher than expected, while 19% of cases showed intrafamilial autosomal dominant transmission. Higher birth weight centile, though lower than reference population, being adequate for gestational age, normal Apgar score and euglycemia characterizing ACC1 resulted associated to a rapid tissue regeneration. Non-isolated ACC, in relation to concomitant congenital anomalies and higher prematurity rate, showed more surgical and medical complications along with the risk of neonatal death. Specifically, newborns with ACC4 were characterized by the frequent necessity of abdominal wall defect repair, responsible for the occurrence of an abdominal compartment syndrome.

Conclusion

Prompt carefully assessment of the newborn with ACC in order to exclude concomitant other congenital malformations, provides clues to the underlying pathophysiology, and to the short-term prognosis. Family should be oriented toward identification of other family members affected by similar pathology, while obstetric history should exclude initial multiple pregnancy with death of a co-twin, placental anomalies and drug assumption. Molecular-genetic diagnosis and genetic counseling are integrative in individualized disease approach.

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

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

The value of endothelin receptor type B promoter methylation as a biomarker for the risk assessment and diagnosis of prostate cancer: A meta-analysis

Previous researches have demonstrated that the methylation status of the EDNRB promoter was associated with the prostate cancer (PCa), but these conclusions remained controversial. Thus, the aim of this meta-analysis was to evaluate the association between EDNRB promoter methylation and the PCa. According to the PRISMA statement, the Web of Science, PubMed, EMBASE, and Cochrane Library databases were retrieved. The ORs and 95 % CIs were analyzed to evaluate the associations between EDNRB promoter methylation and the risk and clinical features of PCa. Heterogeneity among the included studies was estimated by I² statistic and Q test. Publication bias and sensitivity analysis were utilized to test the robustness of our outcomes. In addition, the pooled sensitivity and specificity were calculated to assess the diagnostic value of EDNRB methylation for PCa. Ultimately, 11 eligible studies were included. Under the random-effects model, the pooled OR shown that the frequency of EDNRB methylation was substantially higher in cases compared with controls (OR = 5.42, 95 % CI = 1.98-14.88, P = 0.001). The similar results were also found by the data from TCGA database. Subgroup analysis according to the methylation detection method showed that the heterogeneity in quantitative methylation-specific polymerase chain reaction (qMSP) group was insignificant (I² = 0.0 %, P = 0.669). Moreover, the pooled sensitivity for all-inclusive studies was 0.55 (95 % CI: 0.26-0.81), and the pooled specificity was 0.93 (95 % CI: 0.55-0.99). The methylation of EDNRB promoter might increase the risk of PCa. Meanwhile, EDNRB promoter methylation test combined with PSA testing and/or other biomarkers could be promising diagnostic biomarkers for more accurate detection of PCa.

MEK inhibition remodels the active chromatin landscape and induces SOX10 genomic recruitment in BRAF(V600E) mutant melanoma cells

Background

The MAPK/ERK signaling pathway is an essential regulator of numerous cell processes that are crucial for normal development as well as cancer progression. While much is known regarding MAPK/ERK signal conveyance from the cell membrane to the nucleus, the transcriptional and epigenetic mechanisms that govern gene expression downstream of MAPK signaling are not fully elucidated.

Results

This study employed an integrated epigenome analysis approach to interrogate the effects of MAPK/ERK pathway inhibition on the global transcriptome, the active chromatin landscape, and protein–DNA interactions in 501mel melanoma cells. Treatment of these cells with the small-molecule MEK inhibitor AZD6244 induces hyperpigmentation, widespread gene expression changes including alteration of genes linked to pigmentation, and extensive epigenomic reprogramming of transcriptionally distinct regulatory regions associated with the active chromatin mark H3K27ac. Regulatory regions with differentially acetylated H3K27ac regions following AZD6244 treatment are enriched in transcription factor binding motifs of ETV/ETS and ATF family members as well as the lineage-determining factors MITF and SOX10. H3K27ac-dense enhancer clusters known as super-enhancers show similar transcription factor motif enrichment, and furthermore, these super-enhancers are associated with genes encoding MITF, SOX10, and ETV/ETS proteins. Along with genome-wide resetting of the active enhancer landscape, MEK inhibition also results in widespread SOX10 recruitment throughout the genome, including increased SOX10 binding density at H3K27ac-marked enhancers. Importantly, these MEK inhibitor-responsive enhancers marked by H3K27ac and occupied by SOX10 are located near melanocyte lineage-specific and pigmentation genes and overlap numerous human SNPs associated with pigmentation and melanoma phenotypes, highlighting the variants located within these regions for prioritization in future studies.

Conclusions

These results reveal the epigenetic reprogramming underlying the re-activation of melanocyte pigmentation and developmental transcriptional programs in 501mel cells in response to MEK inhibition and suggest extensive involvement of a MEK-SOX10 axis in the regulation of these processes. The dynamic chromatin changes identified here provide a rich genomic resource for further analyses of the molecular mechanisms governing the MAPK pathway in pigmentation- and melanocyte-associated diseases.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0297-2) contains supplementary material, which is available to authorized users.

Immunomodulatory germline variation associated with the development of multiple primary melanoma (MPM)

Multiple primary melanoma (MPM) has been associated with a higher 10-year mortality risk compared to patients with single primary melanoma (SPM). Given that 3–8% of patients with SPM develop additional primary melanomas, new markers predictive of MPM risk are needed. Based on the evidence that the immune system may regulate melanoma progression, we explored whether germline genetic variants controlling the expression of 41 immunomodulatory genes modulate the risk of MPM compared to patients with SPM or healthy controls. By genotyping these 41 variants in 977 melanoma patients, we found that rs2071304, linked to the expression of SPI1, was strongly associated with MPM risk reduction (OR = 0.60; 95% CI = 0.45–0.81; p = 0.0007) when compared to patients with SPM. Furthermore, we showed that rs6695772, a variant affecting expression of BATF3, is also associated with MPM-specific survival (HR = 3.42; 95% CI = 1.57–7.42; p = 0.0019). These findings provide evidence that the genetic variation in immunomodulatory pathways may contribute to the development of secondary primary melanomas and also associates with MPM survival. The study suggests that inherited host immunity may play an important role in MPM development.

Exploring various polygenic risk scores for skin cancer in the phenomes of the Michigan genomics initiative and the UK Biobank with a visual catalog: PRSWeb

Polygenic risk scores (PRS) are designed to serve as single summary measures that are easy to construct, condensing information from a large number of genetic variants associated with a disease. They have been used for stratification and prediction of disease risk. The primary focus of this paper is to demonstrate how we can combine PRS and electronic health records data to better understand the shared and unique genetic architecture and etiology of disease subtypes that may be both related and heterogeneous. PRS construction strategies often depend on the purpose of the study, the available data/summary estimates, and the underlying genetic architecture of a disease. We consider several choices for constructing a PRS using data obtained from various publicly-available sources including the UK Biobank and evaluate their abilities to predict not just the primary phenotype but also secondary phenotypes derived from electronic health records (EHR). This study was conducted using data from 30,702 unrelated, genotyped patients of recent European descent from the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort within Michigan Medicine. We examine the three most common skin cancer subtypes in the USA: basal cell carcinoma, cutaneous squamous cell carcinoma, and melanoma. Using these PRS for various skin cancer subtypes, we conduct a phenome-wide association study (PheWAS) within the MGI data to evaluate PRS associations with secondary traits. PheWAS results are then replicated using population-based UK Biobank data and compared across various PRS construction methods. We develop an accompanying visual catalog called PRSweb that provides detailed PheWAS results and allows users to directly compare different PRS construction methods.

3'UTR variants of TNS3, PHLDB1, NTN4, and GNG2 genes are associated with IgA nephropathy risk in Chinese Han population

BACKGROUND

Immunoglobulin A nephropathy (IgAN) is the most common primary glomerulonephritis and is characterized by mesangial cell proliferation and agglomeration of the mesangial matrix.

METHODS

In this study, we aimed to explore the role of TNS3, PHLDB1, NTN4, and GNG2 3'untranslated region (3'UTR) polymorphisms with the risk of IgAN in a Chinese Han cohort. A logistic recession model was used to calculate the effects of candidate single nucleotide polymorphism (SNP) on IgAN risk after adjusting age and gender difference. In silico prediction was conducted to identify potential functions of SNPs.

RESULTS

The analysis revealed a significant relationship between the homozygotic genotype for NTN4 rs1362970 A/A and higher risk of IgAN (p = 0.003). Statistically significant associations were found when the sample was stratified by gender and Lee's grade. As a result, NTN4 rs1362970 A/A and GNG2 rs3204008 G/G genotypes were associated with enhanced IgAN risk in males (p = 0.006, p = 0.023, respectively), and the association between the PHLDB1 rs7389 G/T genotype and higher IgAN risk was found in females (p = 0.008). In the Lee's grade III-V subgroup, the rs1369270 in NTN4 was significantly correlated with the risk of IgAN (p = 0.004). Bioinformatics prediction suggested that rs1362970 within NTN4 3'UTR was located in the potential target sequence of miR-483-5p.

CONCLUSIONS

Our research confirmed that NTN4, GNG2, and PHLDB1 gene polymorphisms were implicated in IgAN susceptibility in Chinese Han population. Further research should be conducted to investigate and validate the mechanism by which the above-mentioned polymorphisms affect IgAN.

Unexpected phenotype in a patient with two chromosomal deletions involving 6pter and 22q11

The 6p terminal deletions are rare and usually early diagnosed because of their association with eye and cranio-facial anomalies, particularly as part of Axenfeld-Rieger syndrome in relation with the haploinsufficiency of FOXC1 gene. Deletions in the 22q11 region are frequent, highly correlated with DiGeorge syndrome also named CATCH22, and may be associated with many clinical features of various severities. We report a 31-year-old man with an unbalanced 45,XY,der(6)t(6;22)(p25;q11.2),-22 karyotype leading to monosomies in both 6p25 and 22q11 regions, confirmed by FISH and array-CGH. The length of the deletions was respectively 770 Kb for 6pter and 2.9 Mb for 22q11. This karyotype was discovered at adult age following problems of fertility. The chromosomal formula was unexpected, regarding the patient's medical history and clinical features. This case makes a great example of the difficulties to correlate genotype and phenotype, and furthermore demonstrates the complexity of genetic counselling even in a case with two different chromosomal unbalances.

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

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

Identification of 12 genetic loci associated with human healthspan

Aging populations face diminishing quality of life due to increased disease and morbidity. These challenges call for longevity research to focus on understanding the pathways controlling healthspan. We use the data from the UK Biobank (UKB) cohort and observe that the risks of major chronic diseases increased exponentially and double every eight years, i.e., at a rate compatible with the Gompertz mortality law. Assuming that aging drives the acceleration in morbidity rates, we build a risk model to predict the age at the end of healthspan depending on age, gender, and genetic background. Using the sub-population of 300,447 British individuals as a discovery cohort, we identify 12 loci associated with healthspan at the whole-genome significance level. We find strong genetic correlations between healthspan and all-cause mortality, life-history, and lifestyle traits. We thereby conclude that the healthspan offers a promising new way to interrogate the genetics of human longevity.

Host Characteristics and Risk of Incident Melanoma by Breslow Thickness

BACKGROUND

Several host characteristics, including pigmentary traits (hair color, sunburn susceptibility and tanning ability), number of common nevi (moles), and family history of melanoma, have been associated with risk of melanoma.

METHODS

We prospectively examined the associations between host characteristics and risk of incident melanoma by Breslow thickness (≤1 mm, thin melanoma; or >1 mm, "thicker melanoma") based on the Nurses' Health Study (NHS, n = 86,380 women), NHS II (n = 104,100 women), and Health Professionals Follow-up Study (HPFS, n = 46,934 men).

RESULTS

During 22-30 years' follow-up, a total of 1,813 incident melanoma cases were identified with information on Breslow thickness, 1,392 (76.8%) of which had thin melanoma. No significant differences were observed for thin and thicker melanoma in associations with hair color, sunburn susceptibility, and tanning ability. However, we found significant differences for the association with family history of melanoma, with a higher risk estimate for thicker melanoma [HR = 2.55; 95% confidence interval (CI): 1.91-3.42] than thin melanoma (HR = 1.59; 95% CI: 1.21-2.08; P heterogeneity = 0.02). Interestingly, women and men displayed differential associations between nevi count and risk of melanoma by Breslow thickness, with the association appearing stronger for thicker melanoma than thin melanoma in men (P heterogeneity = 0.01), but not in women.

CONCLUSIONS

Individuals with family history of melanoma may be more likely to develop thicker melanoma. Men with high number of common nevi may tend to develop thicker melanoma, which was not found for women.

IMPACT

The findings further stress the risk of thicker melanoma for individuals with a family history of melanoma and men with a high nevi count.

Genome-wide study of hair colour in UK Biobank explains most of the SNP heritability

Natural hair colour within European populations is a complex genetic trait. Previous work has established that MC1R variants are the principal genetic cause of red hair colour, but with variable penetrance. Here, we have extensively mapped the genes responsible for hair colour in the white, British ancestry, participants in UK Biobank. MC1R only explains 73% of the SNP heritability for red hair in UK Biobank, and in fact most individuals with two MC1R variants have blonde or light brown hair. We identify other genes contributing to red hair, the combined effect of which accounts for ~90% of the SNP heritability. Blonde hair is associated with over 200 genetic variants and we find a continuum from black through dark and light brown to blonde and account for 73% of the SNP heritability of blonde hair. Many of the associated genes are involved in hair growth or texture, emphasising the cellular connections between keratinocytes and melanocytes in the determination of hair colour.

Natural hair colour in Europeans is a complex genetic trait. Here, the authors carry out a genome-wide association study using UK BioBank data, suggesting that in combination with pigmentation genes, variants with roles in hair texture and growth can affect hair colouration or our perception of it.

Shades of complexity: New perspectives on the evolution and genetic architecture of human skin

Like many highly variable human traits, more than a dozen genes are known to contribute to the full range of skin color. However, the historical bias in favor of genetic studies in European and European-derived populations has blinded us to the magnitude of pigmentation's complexity. As deliberate efforts are being made to better characterize diverse global populations and new sequencing technologies, better measurement tools, functional assessments, predictive modeling, and ancient DNA analyses become more widely accessible, we are beginning to appreciate how limited our understanding of the genetic bases of human skin color have been. Novel variants in genes not previously linked to pigmentation have been identified and evidence is mounting that there are hundreds more variants yet to be found. Even for genes that have been exhaustively characterized in European populations like MC1R, OCA2, and SLC24A5, research in previously understudied groups is leading to a new appreciation of the degree to which genetic diversity, epistatic interactions, pleiotropy, admixture, global and local adaptation, and cultural practices operate in population-specific ways to shape the genetic architecture of skin color. Furthermore, we are coming to terms with how factors like tanning response and barrier function may also have influenced selection on skin throughout human history. By examining how our knowledge of pigmentation genetics has shifted in the last decade, we can better appreciate how far we have come in understanding human diversity and the still long road ahead for understanding many complex human traits.

Inherited Genetic Variants Associated with Melanoma BRAF/NRAS Subtypes

BRAF and NRAS mutations arise early in melanoma development, but their associations with low-penetrance melanoma susceptibility loci remain unknown. In the Genes, Environment and Melanoma Study, 1,223 European-origin participants had their incident invasive primary melanomas screened for BRAF/NRAS mutations and germline DNA genotyped for 47 single-nucleotide polymorphisms identified as low-penetrant melanoma-risk variants. We used multinomial logistic regression to simultaneously examine each single-nucleotide polymorphism's relationship to BRAF V600E, BRAF V600K, BRAF other, and NRAS+ relative to BRAF-/NRAS- melanoma adjusted for study features. IRF4 rs12203592*T was associated with BRAF V600E (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.43-0.79) and V600K (OR = 0.65, 95% CI = 0.41-1.03), but not BRAF other or NRAS+ melanoma. A global test of etiologic heterogeneity (Pglobal = 0.001) passed false discovery (Pglobal = 0.0026). PLA2G6 rs132985*T was associated with BRAF V600E (OR = 1.32, 95% CI = 1.05-1.67) and BRAF other (OR = 1.82, 95% CI = 1.11-2.98), but not BRAF V600K or NRAS+ melanoma. The test for etiologic heterogeneity (Pglobal) was 0.005. The IRF4 rs12203592 associations were slightly attenuated after adjustment for melanoma-risk phenotypes. The PLA2G6 rs132985 associations were independent of phenotypes. IRF4 and PLA2G6 inherited genotypes may influence melanoma BRAF/NRAS subtype development.

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

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

Off-target phenotypes in forensic DNA phenotyping and biogeographic ancestry inference: A resource

With recent advances in DNA sequencing technologies it has become feasible and cost effective to genotype larger marker sets for forensic purposes. Two technologies that make use of the larger marker sets have come into focus in forensic research and applications; inference of biogeographic ancestry (BGA) and forensic DNA phenotyping (FDP). These methods hold the promise to reveal information about a yet unknown perpetrator from a DNA sample. In contrast, DNA-profiling, that is a standard practice in case work, relies on matching DNA-profiles between crime scene material and suspects on a database of DNA-profiles. Markers for DNA-profiling were developed under the premise to reveal as little additional information about the human source of the profile as possible, the rationale being that personal privacy rights have to be balanced against the public interest in solving a crime. The same argument holds for markers used in BGA and FDP; these markers might also reveal information on off-target phenotypes (OTPs), that go beyond BGA and the phenotypes targeted in FDP. In particular, health related OTPs might shift the balance between privacy protection and public interest. However, to our knowledge, there is currently no convenient resource available to incorporate knowledge on OTPs in BGA and FDP assay design and application. In order to provide such a resource, we performed a systematic search for OTPs associated with a comprehensive set of markers (1766 SNPs) used or suggested to be used for BGA inference and FDP. In this set, we identified a relatively small number of 27 SNPs (1.53%) that convey information on diverse health related OTPs such as cancer risk, induced asthma, or risk of alcoholism. Some of these SNPs are commonly used for FDP and BGA across different marker sets. We conclude that the effects of SNP markers used in FDP and BGA on OTPs are currently limited, with few exceptions that should be considered in a balanced decision on assay design and application.

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

Background

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

Objective

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

Methods

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

Results

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

Conclusion

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

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

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

Prioritization and functional assessment of noncoding variants associated with complex diseases

Unraveling functional noncoding variants associated with complex diseases is still a great challenge. We present a novel algorithm, Prioritization And Functional Assessment (PAFA), that prioritizes and assesses the functionality of genetic variants by introducing population differentiation measures and recalibrating training variants. Comprehensive evaluations demonstrate that PAFA exhibits much higher sensitivity and specificity in prioritizing noncoding risk variants than existing methods. PAFA achieves improved performance in distinguishing both common and rare recurrent variants from non-recurrent variants by integrating multiple annotations and metrics. An integrated platform was developed, providing comprehensive functional annotations for noncoding variants by integrating functional genomic data, which can be accessed at http://159.226.67.237:8080/pafa .

Genome-wide association study in 176,678 Europeans reveals genetic loci for tanning response to sun exposure

The skin’s tendency to sunburn rather than tan is a major risk factor for skin cancer. Here we report a large genome-wide association study of ease of skin tanning in 176,678 subjects of European ancestry. We identify significant association with tanning ability at 20 loci. We confirm previously identified associations at six of these loci, and report 14 novel loci, of which ten have never been associated with pigmentation-related phenotypes. Our results also suggest that variants at the AHR/AGR3 locus, previously associated with cutaneous malignant melanoma the underlying mechanism of which is poorly understood, might act on disease risk through modulation of tanning ability.

The skin’s tanning response to sun exposure shows great interindividual variability. Here, Visconti et al. perform a genome-wide association study for ease of skin tanning and identify 20 genetic loci, ten of which had not previously been associated with pigmentation-related traits.

Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability

Hair color is one of the most recognizable visual traits in European populations and is under strong genetic control. Here we report the results of a genome-wide association study meta-analysis of almost 300,000 participants of European descent. We identified 123 autosomal and one X-chromosome loci significantly associated with hair color; all but 13 are novel. Collectively, single-nucleotide polymorphisms associated with hair color within these loci explain 34.6% of red hair, 24.8% of blond hair, and 26.1% of black hair heritability in the study populations. These results confirm the polygenic nature of complex phenotypes and improve our understanding of melanin pigment metabolism in humans.

The Ensembl Genome Browser: Strategies for Accessing Eukaryotic Genome Data

The Ensembl Genome Browser provides a wealth of freely available genomic data that can be accessed for many purposes by genetics, genomics, and molecular biology researchers. Herein we present two protocols for exploring different aspects of these data: a phenotype and its associated variants and genes, and a promoter and the epigenetic marks and protein-binding activity associated with it. These workflows illustrate a subset of the data types available through the Ensembl Browser, and can be considered a springboard for further exploration.

The Genetic Overlap Between Hair and Eye Color

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

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

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

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

BACKGROUND

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

A Fast and Accurate Algorithm to Test for Binary Phenotypes and Its Application to PheWAS

The availability of electronic health record (EHR)-based phenotypes allows for genome-wide association analyses in thousands of traits and has great potential to enable identification of genetic variants associated with clinical phenotypes. We can interpret the phenome-wide association study (PheWAS) result for a single genetic variant by observing its association across a landscape of phenotypes. Because a PheWAS can test thousands of binary phenotypes, and most of them have unbalanced or often extremely unbalanced case-control ratios (1:10 or 1:600, respectively), existing methods cannot provide an accurate and scalable way to test for associations. Here, we propose a computationally fast score-test-based method that estimates the distribution of the test statistic by using the saddlepoint approximation. Our method is much (∼100 times) faster than the state-of-the-art Firth's test. It can also adjust for covariates and control type I error rates even when the case-control ratio is extremely unbalanced. Through application to PheWAS data from the Michigan Genomics Initiative, we show that the proposed method can control type I error rates while replicating previously known association signals even for traits with a very small number of cases and a large number of controls.

Adaptation of human skin color in various populations

BACKGROUND

Skin color is a well-recognized adaptive trait and has been studied extensively in humans. Understanding the genetic basis of adaptation of skin color in various populations has many implications in human evolution and medicine.

DISCUSSION

Impressive progress has been made recently to identify genes associated with skin color variation in a wide range of geographical and temporal populations. In this review, we discuss what is currently known about the genetics of skin color variation. We enumerated several cases of skin color adaptation in global modern humans and archaic hominins, and illustrated why, when, and how skin color adaptation occurred in different populations. Finally, we provided a summary of the candidate loci associated with pigmentation, which could be a valuable reference for further evolutionary and medical studies.

CONCLUSION

Previous studies generally indicated a complex genetic mechanism underlying the skin color variation, expanding our understanding of the role of population demographic history and natural selection in shaping genetic and phenotypic diversity in humans. Future work is needed to dissect the genetic architecture of skin color adaptation in numerous ethnic minority groups around the world, which remains relatively obscure compared with that of major continental groups, and to unravel the exact genetic basis of skin color adaptation.

Association study of genetic variation in DNA repair pathway genes and risk of basal cell carcinoma

DNA repair plays a critical role in protecting the genome from ultraviolet radiation and maintaining the genomic integrity of cells. Genetic variants in DNA repair-related genes can influence an individual's DNA repair capacity, which may be related to the risk of developing basal cell carcinoma (BCC). We comprehensively assessed the associations of 2,965 independent single-nucleotide polymorphisms (SNPs) across 165 DNA repair pathway genes with BCC risk in a genome-wide association meta-analysis totaling 17,187 BCC cases and 287,054 controls from two data sets. After multiple testing corrections, we identified three SNPs (rs2805831 upstream of XPA: OR = 0.93, P = 1.35 × 10^-6 ; rs659857 in exon of MUS81: OR = 1.06, P = 3.09 × 10^-6 and rs57343616 in 3' UTR of NABP2: OR = 1.11, P = 6.47 × 10^-6 ) as significantly associated with BCC risk in meta-analysis, and all of them were nominally significant in both data sets. Furthermore, rs659857 [T] was significantly associated with decreased expression of MUS81 mRNA in the expression quantitative trait locus (eQTL) analysis. Our findings suggest that the inherited common variation in three DNA repair genes-XPA, MUS81 and NABP2-may be involved in the development of BCC. To our knowledge, our study is the first report thoroughly examining the effects of SNPs across DNA repair pathway genes on BCC risk based on a genome-wide association meta-analysis.