Genetic variants for smoking behaviour and risk of skin cancer

Observational studies have suggested that smoking may increase the risk of cutaneous squamous cell carcinoma (cSCC) while decreasing the risks of basal cell carcinoma (BCC), and melanoma. However, it remains possible that confounding by other factors may explain these associations. The aim of this investigation was to use Mendelian randomization (MR) to test whether smoking is associated with skin cancer, independently of other factors. Two-sample MR analyses were conducted to determine the causal effect of smoking measures on skin cancer risk using genome-wide association study (GWAS) summary statistics. We used the inverse-variance-weighted estimator to derive separate risk estimates across genetic instruments for all smoking measures. A genetic predisposition to smoking initiation was associated with lower risks of all skin cancer types, although none of the effect estimates reached statistical significance (OR 95% CI BCC 0.91, 0.82-1.01; cSCC 0.82, 0.66-1.01; melanoma 0.91, 0.82-1.01). Results for other measures were similar to smoking initiation with the exception of smoking intensity which was associated with a significantly reduced risk of melanoma (OR 0.67, 95% CI 0.51-0.89). Our findings support the findings of observational studies linking smoking to lower risks of melanoma and BCC. However, we found no evidence that smoking is associated with an elevated risk of cSCC; indeed, our results are most consistent with a decreased risk, similar to BCC and melanoma.

Uncovering the complex relationship between balding, testosterone and skin cancers in men

Male-pattern baldness (MPB) is related to dysregulation of androgens such as testosterone. A previously observed relationship between MPB and skin cancer may be due to greater exposure to ultraviolet radiation or indicate a role for androgenic pathways in the pathogenesis of skin cancers. We dissected this relationship via Mendelian randomization (MR) analyses, using genetic data from recent male-only meta-analyses of cutaneous melanoma (12,232 cases; 20,566 controls) and keratinocyte cancers (KCs) (up to 17,512 cases; >100,000 controls), followed by stratified MR analysis by body-sites. We found strong associations between MPB and the risk of KC, but not with androgens, and multivariable models revealed that this relationship was heavily confounded by MPB single nucleotide polymorphisms involved in pigmentation pathways. Site-stratified MR analyses revealed strong associations between MPB with head and neck squamous cell carcinoma and melanoma, suggesting that sun exposure on the scalp, rather than androgens, is the main driver. Men with less hair covering likely explains, at least in part, the higher incidence of melanoma in men residing in countries with high ambient UV.

Population-based analysis of POT1 variants in a cutaneous melanoma case-control cohort

Pathogenic germline variants in the protection of telomeres 1 gene (POT1) have been associated with predisposition to a range of tumour types, including melanoma, glioma, leukaemia and cardiac angiosarcoma. We sequenced all coding exons of the POT1 gene in 2928 European-descent melanoma cases and 3298 controls, identifying 43 protein-changing genetic variants. We performed POT1-telomere binding assays for all missense and stop-gained variants, finding nine variants that impair or disrupt protein-telomere complex formation, and we further define the role of variants in the regulation of telomere length and complex formation through molecular dynamics simulations. We determine that POT1 coding variants are a minor contributor to melanoma burden in the general population, with only about 0.5% of melanoma cases carrying germline pathogenic variants in this gene, but should be screened in individuals with a strong family history of melanoma and/or multiple malignancies.

Genome-wide analyses characterize shared heritability among cancers and identify novel cancer susceptibility regions

BACKGROUND

The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci.

METHODS

We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci.

RESULTS

We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci.

CONCLUSIONS

Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Metabolic drivers of dysglycemia in pregnancy: ethnic-specific GWAS of 146 metabolites and 1-sample Mendelian randomization analyses in a UK multi-ethnic birth cohort

Introduction

Gestational diabetes mellitus (GDM) is the most common pregnancy complication worldwide and is associated with short- and long-term health implications for both mother and child. Prevalence of GDM varies between ethnicities, with South Asians (SAs) experiencing up to three times the risk compared to white Europeans (WEs). Recent evidence suggests that underlying metabolic difference contribute to this disparity, but an investigation of causality is required.

Methods

To address this, we paired metabolite and genomic data to evaluate the causal effect of 146 distinct metabolic characteristics on gestational dysglycemia in SAs and WEs. First, we performed 292 GWASs to identify ethnic-specific genetic variants associated with each metabolite (P ≤ 1 x 10-5) in the Born and Bradford cohort (3688 SA and 3354 WE women). Following this, a one-sample Mendelian Randomisation (MR) approach was applied for each metabolite against fasting glucose and 2-hr post glucose at 26-28 weeks gestation. Additional GWAS and MR on 22 composite measures of metabolite classes were also conducted.

Results

This study identified 15 novel genome-wide significant (GWS) SNPs associated with tyrosine in the FOXN and SLC13A2 genes and 1 novel GWS SNP (currently in no known gene) associated with acetate in SAs. Using MR approach, 14 metabolites were found to be associated with postprandial glucose in WEs, while in SAs a distinct panel of 11 metabolites were identified. Interestingly, in WEs, cholesterols were the dominant metabolite class driving with dysglycemia, while in SAs saturated fatty acids and total fatty acids were most commonly associated with dysglycemia.

Discussion

In summary, we confirm and demonstrate the presence of ethnic-specific causal relationships between metabolites and dysglycemia in mid-pregnancy in a UK population of SA and WE pregnant women. Future work will aim to investigate their biological mechanisms on dysglycemia and translating this work towards ethnically tailored GDM prevention strategies.

The association between genetically elevated polyunsaturated fatty acids and risk of cancer

BACKGROUND

The causal relevance of polyunsaturated fatty acids (PUFAs) for risk of site-specific cancers remains uncertain.

METHODS

Using a Mendelian randomization (MR) framework, we assessed the causal relevance of PUFAs for risk of cancer in European and East Asian ancestry individuals. We defined the primary exposure as PUFA desaturase activity, proxied by rs174546 at the FADS locus. Secondary exposures were defined as omega 3 and omega 6 PUFAs that could be proxied by genetic polymorphisms outside the FADS region. Our study used summary genetic data on 10 PUFAs and 67 cancers, corresponding to 562,871 cases and 1,619,465 controls, collected by the Fatty Acids in Cancer Mendelian Randomization Collaboration. We estimated odds ratios (ORs) for cancer per standard deviation increase in genetically proxied PUFA exposures.

FINDINGS

Genetically elevated PUFA desaturase activity was associated (P < 0.0007) with higher risk (OR [95% confidence interval]) of colorectal cancer (1.09 [1.07-1.11]), esophageal squamous cell carcinoma (1.16 [1.06-1.26]), lung cancer (1.06 [1.03-1.08]) and basal cell carcinoma (1.05 [1.02-1.07]). There was little evidence for associations with reproductive cancers (OR = 1.00 [95% CI: 0.99-1.01]; Pheterogeneity = 0.25), urinary system cancers (1.03 [0.99-1.06], Pheterogeneity = 0.51), nervous system cancers (0.99 [0.95-1.03], Pheterogeneity = 0.92) or blood cancers (1.01 [0.98-1.04], Pheterogeneity = 0.09). Findings for colorectal cancer and esophageal squamous cell carcinoma remained compatible with causality in sensitivity analyses for violations of assumptions. Secondary MR analyses highlighted higher omega 6 PUFAs (arachidonic acid, gamma-linolenic acid and dihomo-gamma-linolenic acid) as potential mediators. PUFA biosynthesis is known to interact with aspirin, which increases risk of bleeding and inflammatory bowel disease. In a phenome-wide MR study of non-neoplastic diseases, we found that genetic lowering of PUFA desaturase activity, mimicking a hypothetical intervention to reduce cancer risk, was associated (P < 0.0006) with increased risk of inflammatory bowel disease but not bleeding.

INTERPRETATION

The PUFA biosynthesis pathway may be an intervention target for prevention of colorectal cancer and esophageal squamous cell carcinoma but with potential for increased risk of inflammatory bowel disease.

FUNDING

Cancer Resesrch UK (C52724/A20138, C18281/A19169). UK Medical Research Council (MR/P014054/1). National Institute for Health Research (NIHR202411). UK Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4). National Cancer Institute (R00 CA215360). National Institutes of Health (U01 CA164973, R01 CA60987, R01 CA72520, U01 CA74806, R01 CA55874, U01 CA164973 and U01 CA164973).

Investigation of Shared Genetic Risk Factors Between Parkinson's Disease and Cancers

BACKGROUND

Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties.

OBJECTIVE

We used results from genome-wide association studies (GWASs) to study the genetic correlation between PD and different cancers to identify common genetic risk factors.

METHODS

We used individual data for participants of European ancestry from the Courage-PD (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease; PD, N = 16,519) and EPITHYR (differentiated thyroid cancer, N = 3527) consortia and summary statistics of GWASs from iPDGC (International Parkinson Disease Genomics Consortium; PD, N = 482,730), Melanoma Meta-Analysis Consortium (MMAC), Breast Cancer Association Consortium (breast cancer), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (prostate cancer), International Lung Cancer Consortium (lung cancer), and Ovarian Cancer Association Consortium (ovarian cancer) (N comprised between 36,017 and 228,951 for cancer GWASs). We estimated the genetic correlation between PD and cancers using linkage disequilibrium score regression. We studied the association between PD and polymorphisms associated with cancers, and vice versa, using cross-phenotypes polygenic risk score (PRS) analyses.

RESULTS

We confirmed a previously reported positive genetic correlation of PD with melanoma (Gcorr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (Gcorr = 0.11 [0.03; 0.19]). There was a significant inverse association between the PRS for ovarian cancer and PD (odds ratio [OR] = 0.89 [0.84; 0.94]). Conversely, the PRS of PD was positively associated with breast cancer (OR = 1.08 [1.06; 1.10]) and inversely associated with ovarian cancer (OR = 0.95 [0.91; 0.99]). The association between PD and ovarian cancer was mostly driven by rs183211 located in an intron of the NSF gene (17q21.31).

CONCLUSIONS

We show evidence in favor of a contribution of pleiotropic genes to the association between PD and specific cancers. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Massively parallel reporter assays and variant scoring identified functional variants and target genes for melanoma loci and highlighted cell-type specificity

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.

Differences in the management of patients requiring an emergency resection for colonic cancer in two European populations

Background

Patients with colonic cancer who require emergency colonic cancer surgery often experience poorer outcomes compared with their elective counterparts. In this setting, several treatments approaches are available. In 2009, Danish guidelines recommended treatment with stent for obstruction in left-sided tumours as a bridge to surgery, if expertise is accessible. The aim of this study was to compare the use of elective and emergency resections for colonic cancer and postoperative mortality in two similar demographic populations.

Methods

All patients who underwent a major resection for colonic cancer, between 2005 and 2016 in Denmark and Yorkshire (UK) were identified. The proportion undergoing emergency surgery, the proportion receiving a stent procedure before their resection, and 30-day postoperative mortality were compared between the populations. Logistic regression was used to determine changes in the proportion of those undergoing emergency surgery and 30-day postoperative mortality.

Results

Out of 45 397 patients treated during the study interval, 41 880 were selected. Emergency surgery decreased in Denmark from 16.6 per cent in 2005–07 to 12.9 per cent in 2014–16, but increased in Yorkshire (13.5 per cent to 16.8 per cent). Danish patients with left-sided tumours were less likely to undergo emergency surgery (risk ratio 0.90, 95 per cent c.i. 0.82 to 0.99) and an increase in stent use coincided with a statistically significant decrease in emergency surgery in these patients. Thirty-day postoperative mortality in all resections (elective and emergency) decreased in both populations, but a larger decrease was observed in Denmark (7.7 per cent to 3.0 per cent in Denmark and 7.1 per cent to 3.3 per cent in Yorkshire).

Conclusion

Patients in Denmark experienced a reduction in the use of emergency resection and increase in stenting procedures, following the policy implemented in some departments of converting potential emergency resections into elective resections.

Interplay between demographic, clinical and polygenic risk factors for severe COVID-19

BACKGROUND

We aimed to identify clinical, socio-demographic and genetic risk factors for severe COVID-19 (hospitalization, critical care admission or death) in the general population.

METHODS

In this observational study, we identified 9560 UK Biobank participants diagnosed with COVID-19 during 2020. A polygenic risk score (PRS) for severe COVID-19 was derived and optimized using publicly available European and trans-ethnic COVID-19 genome-wide summary statistics. We estimated the risk of hospital or critical care admission within 28 days or death within 100 days following COVID-19 diagnosis, and assessed associations with socio-demographic factors, immunosuppressant use and morbidities reported at UK Biobank enrolment (2006-2010) and the PRS. To improve biological understanding, pathway analysis was performed using genetic variants comprising the PRS.

RESULTS

We included 9560 patients followed for a median of 61 (interquartile range = 34-88) days since COVID-19 diagnosis. The risk of severe COVID-19 increased with age and obesity, and was higher in men, current smokers, those living in socio-economically deprived areas, those with historic immunosuppressant use and individuals with morbidities and higher co-morbidity count. An optimized PRS, enriched for single-nucleotide polymorphisms in multiple immune-related pathways, including the 'oligoadenylate synthetase antiviral response' and 'interleukin-10 signalling' pathways, was associated with severe COVID-19 (adjusted odds ratio 1.32, 95% CI 1.11-1.58 for the highest compared with the lowest PRS quintile).

CONCLUSION

This study conducted in the pre-SARS-CoV-2-vaccination era, emphasizes the novel insights to be gained from using genetic data alongside commonly considered clinical and socio-demographic factors to develop greater biological understanding of severe COVID-19 outcomes.

Higher polygenic risk for melanoma is associated with improved survival in a high ultraviolet radiation setting

BACKGROUND

The role of germline genetic factors in determining survival from cutaneous melanoma (CM) is not well understood.

OBJECTIVE

To perform a genome-wide association study (GWAS) meta-analysis of melanoma-specific survival (MSS), and test whether a CM-susceptibility polygenic risk score (PRS) is associated with MSS.

METHODS

We conducted two Cox proportional-hazard GWAS of MSS using data from the Melanoma Institute Australia, a high ultraviolet (UV) radiation setting (MIA; 5,762 patients with melanoma; 800 melanoma deaths) and UK Biobank (UKB: 5,220 patients with melanoma; 241 melanoma deaths), and combined them in a fixed-effects meta-analysis. Significant (P < 5 × 10-8) results were investigated in the Leeds Melanoma Cohort (LMC; 1,947 patients with melanoma; 370 melanoma deaths). We also developed a CM-susceptibility PRS using a large independent GWAS meta-analysis (23,913 cases, 342,870 controls). The PRS was tested for an association with MSS in the MIA and UKB cohorts.

RESULTS

Two loci were significantly associated with MSS in the meta-analysis of MIA and UKB with lead SNPs rs41309643 (G allele frequency 1.6%, HR = 2.09, 95%CI = 1.61-2.71, P = 2.08 × 10-8) on chromosome 1, and rs75682113 (C allele frequency 1.8%, HR = 2.38, 95%CI = 1.77-3.21, P = 1.07 × 10-8) on chromosome 7. While neither SNP replicated in the LMC, rs75682113 was significantly associated in the combined discovery and replication sets. After adjusting for age at diagnosis, sex and the first ten principal components, a one standard deviation increase in the CM-susceptibility PRS was associated with improved MSS in the discovery meta-analysis (HR = 0.88, 95% CI = 0.83-0.94, P = 6.93 × 10-5; I2 = 88%). However, this was only driven by the high UV setting cohort (MIA HR = 0.84, 95% CI = 0.78-0.90).

CONCLUSION

We found two loci potentially associated with MSS. Increased genetic susceptibility to develop CM is associated with improved MSS in a high UV setting.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

Multi-Trait Genetic Analysis Identifies Autoimmune Loci Associated with Cutaneous Melanoma

Genome-wide association studies (GWAS) have identified a number of risk loci for cutaneous melanoma. Cutaneous melanoma shares overlapping genetic risk (genetic correlation) with a number of other traits, including its risk factors such as sunburn propensity. This genetic correlation can be exploited to identify additional cutaneous melanoma risk loci by multitrait analysis of GWAS (MTAG). We used bivariate linkage disequilibrium-score regression score regression to identify traits that are genetically correlated with clinically confirmed cutaneous melanoma and then used publicly available GWAS for these traits in a multitrait analysis of GWAS. Multitrait analysis of GWAS allows GWAS to be combined while accounting for sample overlap and incomplete genetic correlation. We identified a total of 74 genome-wide independent loci, 19 of them were not previously reported in the input cutaneous melanoma GWAS meta-analysis. Of these loci, 55 were replicated (P < 0.05/74, Bonferroni-corrected P-value in two independent cutaneous melanoma replication cohorts from Melanoma Institute Australia and 23andMe, Inc. Among the, to our knowledge, previously unreported cutaneous melanoma loci are ones that have also been associated with autoimmune traits including rs715199 near LPP and rs10858023 near AP4B1. Our analysis indicates genetic correlation between traits can be leveraged to identify new risk genes for cutaneous melanoma.

Independent evaluation of melanoma polygenic risk scores in UK and Australian prospective cohorts

BACKGROUND

Previous studies suggest that polygenic risk scores (PRSs) may improve melanoma risk stratification. However, there has been limited independent validation of PRS-based risk prediction, particularly assessment of calibration (comparing predicted to observed risks).

OBJECTIVES

To evaluate PRS-based melanoma risk prediction in prospective UK and Australian cohorts with European ancestry.

METHODS

We analysed invasive melanoma incidence in the UK Biobank (UKB; n = 395 647, 1651 cases) and a case-cohort nested within the Melbourne Collaborative Cohort Study (MCCS, Australia; n = 4765, 303 cases). Three PRSs were evaluated: 68 single-nucleotide polymorphisms (SNPs) at 54 loci from a 2020 meta-analysis (PRS68), 50 SNPs significant in the 2020 meta-analysis excluding UKB (PRS50) and 45 SNPs at 21 loci known in 2018 (PRS45). Ten-year melanoma risks were calculated from population-level cancer registry data by age group and sex, with and without PRS adjustment.

RESULTS

Predicted absolute melanoma risks based on age and sex alone underestimated melanoma incidence in the UKB [ratio of expected/observed cases: E/O = 0·65, 95% confidence interval (CI) 0·62-0·68] and MCCS (E/O = 0·63, 95% CI 0·56-0·72). For UKB, calibration was improved by PRS adjustment, with PRS50-adjusted risks E/O = 0·91, 95% CI 0·87-0·95. The discriminative ability for PRS68- and PRS50-adjusted absolute risks was higher than for risks based on age and sex alone (Δ area under the curve 0·07-0·10, P < 0·0001), and higher than for PRS45-adjusted risks (Δ area under the curve 0·02-0·04, P < 0·001).

CONCLUSIONS

A PRS derived from a larger, more diverse meta-analysis improves risk prediction compared with an earlier PRS, and might help tailor melanoma prevention and early detection strategies to different risk levels. Recalibration of absolute risks may be necessary for application to specific populations.

Functional annotation and investigation of the 10q24.33 melanoma risk locus identifies a common variant that influences transcriptional regulation of OBFC1

The 10q24.33 locus is known to be associated with susceptibility to cutaneous malignant melanoma (CMM), but the mechanisms underlying this association have been not extensively investigated. We carried out an integrative genomic analysis of 10q24.33 using epigenomic annotations and in vitro reporter gene assays to identify regulatory variants. We found two putative functional single nucleotide polymorphisms (SNPs) in an enhancer and in the promoter of OBFC1, respectively, in neural crest and CMM cells, one, rs2995264, altering enhancer activity. The minor allele G of rs2995264 correlated with lower OBFC1 expression in 470 CMM tumors and was confirmed to increase the CMM risk in a cohort of 484 CMM cases and 1801 controls of Italian origin. Hi-C and chromosome conformation capture (3C) experiments showed the interaction between the enhancer-SNP region and the promoter of OBFC1 and an isogenic model characterized by CRISPR-Cas9 deletion of the enhancer-SNP region confirmed the potential regulatory effect of rs2995264 on OBFC1 transcription. Moreover, the presence of G-rs2995264 risk allele reduced the binding affinity of the transcription factor MEOX2. Biologic investigations showed significant cell viability upon depletion of OBFC1, specifically in CMM cells that were homozygous for the protective allele. Clinically, high levels of OBFC1 expression associated with histologically favorable CMM tumors. Finally, preliminary results suggested the potential effect of decreased OBFC1 expression on telomerase activity in tumorigenic conditions. Our results support the hypothesis that reduced expression of OBFC1 gene through functional heritable DNA variation can contribute to malignant transformation of normal melanocytes.

Ethnic-specific associations between dietary consumption and gestational diabetes mellitus incidence: A meta-analysis

Globally, one in seven pregnant women are diagnosed with gestational diabetes mellitus (GDM), conferring short- and long-term health risks to both mother and child. While dietary prevention strategies are common in clinical practice, their effectiveness in different ethnicities is uncertain. To better inform prevention strategies, here the effects of unhealthy and healthy diets on GDM risk within distinct ethnic or cultural populations and geographic regions were evaluated and summarised. Pubmed, Scopus, Cochrane and OVID were systematically searched to identify randomised controlled trials (RCTs) and observational studies that investigated diet and GDM. A grouped analysis of common 'healthy' and 'unhealthy' diets was performed first, before analysing individual dietary patterns (e.g., prudent, Mediterranean). Random effect models and dose response analyses were performed where possible. PROSPERO (CRD42019140873). Thirty-eight publications provided information on 5 population groups: white European (WE), Asian, Iranian, Mediterranean and Australian. No associations were identified between healthy diets and GDM incidence in RCTs in any population. However, when synthesizing observational studies, healthy diets reduced odds of GDM by 23% (95% CI: 0.70-0.89, p<0.001, I2 = 75%), while unhealthy diets increased odds of GDM by 61% (95% CI: 1.41-1.81, p<0.0001, I2 = 0%) in WE women. No evidence of consistent effects in other populations were observed, even when adequately powered. Diet consistently associated with GDM risk in WEs but not in other populations. Heterogenous use and reporting of ethnically and culturally appropriate diets and dietary assessment tools, particularly in RCTs, raises uncertainty regarding the lack of association found in non-WE populations. Future studies require the use of culturally appropriate tools to confidently evaluate dietary and metabolic mediators of GDM and inform culturally-specific dietary prevention strategies.

Fine-scale population structure and demographic history of British Pakistanis

Previous genetic and public health research in the Pakistani population has focused on the role of consanguinity in increasing recessive disease risk, but little is known about its recent population history or the effects of endogamy. Here, we investigate fine-scale population structure, history and consanguinity patterns using genotype chip data from 2,200 British Pakistanis. We reveal strong recent population structure driven by the biraderi social stratification system. We find that all subgroups have had low recent effective population sizes (N_e), with some showing a decrease 15‒20 generations ago that has resulted in extensive identity-by-descent sharing and homozygosity, increasing the risk of recessive disorders. Our results from two orthogonal methods (one using machine learning and the other coalescent-based) suggest that the detailed reporting of parental relatedness for mothers in the cohort under-represents the true levels of consanguinity. These results demonstrate the impact of cultural practices on population structure and genomic diversity in Pakistanis, and have important implications for medical genetic studies.

An Overview of Strategies for Detecting Genotype-Phenotype Associations Across Ancestrally Diverse Populations

Genome-wide association studies (GWAS) have been very successful at identifying genetic variants influencing a large number of traits. Although the great majority of these studies have been performed in European-descent individuals, it has been recognised that including populations with differing ancestries enhances the potential for identifying causal SNPs due to their differing patterns of linkage disequilibrium. However, when individuals from distinct ethnicities are included in a GWAS, it is necessary to implement a number of control steps to ensure that the identified associations are real genotype-phenotype relationships. In this Review, we discuss the analyses that are required when performing multi-ethnic studies, including methods for determining ancestry at the global and local level for sample exclusion, controlling for ancestry in association testing, and post-GWAS interrogation methods such as genomic control and meta-analysis. We hope that this overview provides a primer for those researchers interested in including distinct populations in their studies.

Association between a 46-SNP Polygenic Risk Score and melanoma risk in Dutch patients with familial melanoma

BACKGROUND

Familial clustering of melanoma suggests a shared genetic predisposition among family members, but only 10%-40% of familial cases carry a pathogenic variant in a known high-risk melanoma susceptibility gene. We investigated whether a melanoma-specific Polygenic Risk Score (PRS) is associated with melanoma risk in patients with genetically unexplained familial melanoma.

METHODS

Dutch familial melanoma cases (n=418) were genotyped for 46 SNPs previously identified as independently associated with melanoma risk. The 46-SNP PRS was calculated and standardised to 3423 healthy controls (sPRS) and the association between PRS and melanoma risk was modelled using logistic regression. Within the case series, possible differences were further explored by investigating the PRS in relation to (1) the number of primary melanomas in a patient and (2) the extent of familial clustering of melanoma.

RESULTS

The PRS was significantly associated with melanoma risk, with a per-SD OR of 2.12 (95% CI 1.90 to 2.35, p<0.001), corresponding to a 5.70-fold increased risk (95% CI 3.93 to 8.28) when comparing the top 90th to the middle 40-60th PRS percentiles. The mean PRS was significantly higher in cases with multiple primary melanomas than in cases with a single melanoma (sPRS 1.17 vs 0.71, p=0.001). Conversely, cases from high-density melanoma families had a lower (but non-significant) mean PRS than cases from low-density families (sPRS 0.60 vs 0.94, p=0.204).

CONCLUSION

Our work underlines the significance of a PRS in determining melanoma susceptibility and encourages further exploration of the diagnostic value of a PRS in genetically unexplained melanoma families.

Influence of age on surgical treatment and postoperative outcomes of patients with colorectal cancer in Denmark and Yorkshire, England

AIM

Denmark and Yorkshire are demographically similar and both have undergone changes in their management of colorectal cancer to improve outcomes. The differential provision of surgical treatment, especially in the older age groups, may contribute to the magnitude of improved survival rates. This study aimed to identify differences in the management of colorectal cancer surgery and postoperative outcomes according to patient age between Denmark and Yorkshire.

METHOD

This was a retrospective population-based study of colorectal cancer patients diagnosed in Denmark and Yorkshire between 2005 and 2016. Proportions of patients undergoing major surgical resection, postoperative mortality and relative survival were compared between Denmark and Yorkshire across several age groups (18-59, 60-69, 70-79 and ≥80 years) and over time.

RESULTS

The use of major surgical resection was higher in Denmark than in Yorkshire, especially for patients aged ≥80 years (70.5% versus 50.5% for colon cancer, 49.3% versus 38.1% for rectal cancer). Thirty-day postoperative mortality for Danish patients aged ≥80 years was significantly higher than that for Yorkshire patients with colonic cancer [OR (95% CI) = 1.22 (1.07, 1.38)] but not for rectal cancer or for 1-year postoperative mortality. Relative survival significantly increased in all patients aged ≥80 years except for Yorkshire patients with colonic cancer.

CONCLUSION

This study suggests that there are major differences between the management of elderly patients with colorectal cancer between the two populations. Improved selection for surgery and better peri- and postoperative care in these patients appears to improve long-term outcomes, but may come at the cost of a higher 30-day mortality.

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

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

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Genetic Relationship Between Endometriosis and Melanoma

Epidemiological studies have observed that risk of endometriosis is associated with history of cutaneous melanoma and vice versa. Evidence for shared biological mechanisms between the two traits is limited. The aim of this study was to investigate the genetic correlation and causal relationship between endometriosis and melanoma. Summary statistics from genome-wide association meta-analyses (GWAS) for endometriosis and melanoma were used to estimate the genetic correlation between the traits and Mendelian randomization was used to test for a causal association. When using summary statistics from separate female and male melanoma cohorts we identified a significant positive genetic correlation between melanoma in females and endometriosis (rg = 0.144, se = 0.065, p = 0.025). However, we find no evidence of a correlation between endometriosis and melanoma in males or a combined melanoma dataset. Endometriosis was not genetically correlated with skin color, red hair, childhood sunburn occasions, ease of skin tanning, or nevus count suggesting that the correlation between endometriosis and melanoma in females is unlikely to be influenced by pigmentary traits. Mendelian Randomization analyses also provided evidence for a relationship between the genetic risk of melanoma in females and endometriosis. Colocalization analysis identified 27 genomic loci jointly associated with the two diseases regions that contain different causal variants influencing each trait independently. This study provides evidence of a small genetic correlation and relationship between the genetic risk of melanoma in females and endometriosis. Genetic risk does not equate to disease occurrence and differences in the pathogenesis and age of onset of both diseases means it is unlikely that occurrence of melanoma causes endometriosis. This study instead provides evidence that having an increased genetic risk for melanoma in females is related to increased risk of endometriosis. Larger GWAS studies with increased power will be required to further investigate these associations.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

Risk factors for melanoma by anatomical site: an evaluation of aetiological heterogeneity

BACKGROUND

Melanoma aetiology has been proposed to have two pathways, which are determined by naevi and type of sun exposure and related to the anatomical site where melanoma develops.

OBJECTIVES

We examined associations with melanoma by anatomical site for a comprehensive set of risk factors including pigmentary and naevus phenotypes, ultraviolet radiation exposure and polygenic risk.

METHODS

We analysed harmonized data from 2617 people with incident first invasive melanoma and 975 healthy controls recruited through two population-based case-control studies in Australia and the UK. Questionnaire data were collected by interview using a single protocol, and pathway-specific polygenic risk scores were derived from DNA samples. We estimated adjusted odds ratios using unconditional logistic regression that compared melanoma cases at each anatomical site with all controls.

RESULTS

When cases were compared with control participants, there were stronger associations for many naevi vs. no naevi for melanomas on the trunk, and upper and lower limbs than on the head and neck (P-heterogeneity < 0·001). Very fair skin (vs. olive/brown skin) was more weakly related to melanoma on the trunk than to melanomas at other sites (P-heterogeneity = 0·04). There was no significant difference by anatomical site for polygenic risk. Increased weekday sun exposure was positively associated with melanoma on the head and neck but not on other sites.

CONCLUSIONS

We found evidence of aetiological heterogeneity for melanoma, supporting the dual pathway hypothesis. These findings enhance understanding of risk factors for melanoma and can guide prevention and skin examination education and practices.

Body mass index and height and risk of cutaneous melanoma: Mendelian randomization analyses

BACKGROUND

Height and body mass index (BMI) have both been positively associated with melanoma risk, although findings for BMI have been less consistent than height. It remains unclear, however, whether these associations reflect causality or are due to residual confounding by environmental and lifestyle risk factors. We re-evaluated these associations using a two-sample Mendelian randomization (MR) approach.

METHODS

We identified single nucleotide polymorphisms (SNPs) for BMI and height from separate genome-wide association study (GWAS) meta-analyses. We obtained melanoma SNPs from the most recent melanoma GWAS meta-analysis comprising 12 874 cases and 23 203 controls. We used the inverse variance-weighted estimator to derive separate causal risk estimates across all SNP instruments for BMI and height.

RESULTS

Based on the combined estimate derived from 730 SNPs for BMI, we found no evidence of an association between genetically predicted BMI and melanoma [odds ratio (OR) per one standard deviation (1 SD) (4.6 kg/m2) increase in BMI 1.00, 95% confidence interval (CI): 0.91-1.11]. In contrast, we observed a positive association between genetically-predicted height (derived from a pooled estimate of 3290 SNPs) and melanoma risk [OR 1.08, 95% CI: 1.02-1.13, per 1 SD (9.27 cm) increase in height]. Sensitivity analyses using two alternative MR methods yielded similar results.

CONCLUSIONS

These findings provide no evidence for a causal association between higher BMI and melanoma, but support the notion that height is causally associated with melanoma risk. Mechanisms through which height influences melanoma risk remain unclear, and it remains possible that the effect could be mediated through diverse pathways including growth factors and even socioeconomic status.

A comprehensive re-assessment of the association between vitamin D and cancer susceptibility using Mendelian randomization

Previous Mendelian randomization (MR) studies on 25-hydroxyvitamin D (25(OH)D) and cancer have typically adopted a handful of variants and found no relationship between 25(OH)D and cancer; however, issues of horizontal pleiotropy cannot be reliably addressed. Using a larger set of variants associated with 25(OH)D (74 SNPs, up from 6 previously), we perform a unified MR analysis to re-evaluate the relationship between 25(OH)D and ten cancers. Our findings are broadly consistent with previous MR studies indicating no relationship, apart from ovarian cancers (OR 0.89; 95% C.I: 0.82 to 0.96 per 1 SD change in 25(OH)D concentration) and basal cell carcinoma (OR 1.16; 95% C.I.: 1.04 to 1.28). However, after adjustment for pigmentation related variables in a multivariable MR framework, the BCC findings were attenuated. Here we report that lower 25(OH)D is unlikely to be a causal risk factor for most cancers, with our study providing more precise confidence intervals than previously possible.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma.

Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness.

To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness.

The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry.

We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10⁻⁸) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Multiplex melanoma families are enriched for polygenic risk

Cancers, including cutaneous melanoma, can cluster in families. In addition to environmental etiological factors such as ultraviolet radiation, cutaneous melanoma has a strong genetic component. Genetic risks for cutaneous melanoma range from rare, high-penetrance mutations to common, low-penetrance variants. Known high-penetrance mutations account for only about half of all densely affected cutaneous melanoma families, and the causes of familial clustering in the remainder are unknown. We hypothesize that some clustering is due to the cumulative effect of a large number of variants of individually small effect. Common, low-penetrance genetic risk variants can be combined into polygenic risk scores. We used a polygenic risk score for cutaneous melanoma to compare families without known high-penetrance mutations with unrelated melanoma cases and melanoma-free controls. Family members had significantly higher mean polygenic load for cutaneous melanoma than unrelated cases or melanoma-free healthy controls (Bonferroni-corrected t-test P = 1.5 × 10-5 and 6.3 × 10-45, respectively). Whole genome sequencing of germline DNA from 51 members of 21 families with low polygenic risk for melanoma identified a CDKN2A p.G101W mutation in a single family but no other candidate high-penetrance melanoma susceptibility genes. This work provides further evidence that melanoma, like many other common complex disorders, can arise from the joint action of multiple predisposing factors, including rare high-penetrance mutations, as well as via a combination of large numbers of alleles of small effect.

Abstract 1194: Cross-cancer GWAS meta-analysis of more than 370,000 cases and 530,000 controls identifies multiple novel cancer risk regions

Genome-wide association studies (GWAS) have identified hundreds of common, low-penetrance alleles associated with cancer risk. However, known rare and common risk alleles only explain between 10% and 30% of the familial relative risk for different cancers and multiple lines of evidence indicate that many more risk alleles remain to be discovered. We have demonstrated genetic correlations between cancers, reflecting a shared genetic origin for solid tumors. These results suggest that jointly analyzing multiple cancer sites will lead to the discovery of novel risk regions.

We conducted a cross-cancer GWAS meta-analysis by leveraging GWAS summary statistics from 12 solid cancers (breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate and renal cancers) with a total of 373,818 cases and 532,382 controls of European ancestry. All studies had been imputed to either 1,000 Genomes or the Haplotype Reference Consortium panel. We conducted four meta-analysis using (1) fixed-effect, (2) random-effect, (3) one-sided subset (ASSET) and (4) two-sided subset (ASSET) models. The subset analysis were conducted assuming either the same direction of effects across cancers (one-sided ASSET) or allowed for opposite direction of effects across cancers (two-sided ASSET). In all analyses, we used tetrachoric correlations to account for sample overlap across cancer sites. In total, we tested 10,223,013 variants for association. We considered regions with a p-value&lt;1.25 × 10−8 in at least one of the four meta-analysis approaches and located at least 500kb away from known cancer risk SNPs as novel.

We identified eight novel regions that reached genome-wide significance. Of those eight regions, two were identified from fixed-effects meta-analysis, three from random effects meta-analysis, one from the two-sided subset analysis, and two regions (15.q15.3 and 21q22.3) were identified at p&lt;1.25 × 10−8 in three of the meta-analysis approaches. Among novel findings is a deleterious missense variant located in RREB1 previously associated with type 2 diabetes, a deleterious missense variant located in DSTYK previously associated with waist-to-hip ratio and triglycerides, and an intergenic variant in linkage disequilibrium (LD) with variants in TMEM18, previously associated with body mass index. Other potential target genes among the newly discovered regions include TP53BP1 and PCNT, both previously implicated in carcinogenesis.

Citation Format: Sara Lindström, Siddhartha Kar, Lu Wang, Constance Turman, James MacDonald, Theo Bammler, BCAC, OCAC, PRACTICAL, Jeroen Huyghe, Stephanie Schmit, Tracy A. O'Mara, Deborah J. Thompson, Puya Gharahkhani, Stuart MacGregor, Paul Brennan, Richard S. Houlston, Beatrice S. Melin, Christopher I. Amos, James McKay, Mark M. Iles, Matthew H. Law, Alison Klein, Laufey Amundadottir, Bogdan Pasaniuc, Paul Pharoah, Rayjean J. Hung, Peter Kraft. Cross-cancer GWAS meta-analysis of more than 370,000 cases and 530,000 controls identifies multiple novel cancer risk regions [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1194.

Abstract 30: Cross-cancer cross-tissue transcriptome-wide association study (TWAS) of 11 cancers identifies 56 novel genes

Although cancer is a heterogeneous disease, there are shared hallmark mechanisms across multiple tumor types. Because of this, identifying genes associated with multiple cancer types has the potential to shed light on general oncogenic mechanisms. Conversely, integrating evidence for genetic association across multiple cancers could identify novel genes missed in single-cancer analyses. Transcriptome Wide Association Studies (TWAS) have been successful in identifying genes associated with individual cancers. TWAS test whether genetically-predicted tissue-specific gene expression levels are associated with cancer risk. Although cross-cancer genome-wide association studies (GWAS) analyses have been performed previously, no cross-cancer TWAS has been conducted to date. Here, we implement a pipeline to perform cross-cancer, cross-tissue TWAS analysis. We use newly-developed multi-trait TWAS test statistics to integrate the TWAS results for association between 11 separated cancers and predicted gene expression in each of 43 GTEx tissues; these include a "sum" test and a "variance components" test, analogous to fixed- and random-effects meta-analyses. We then integrated the results across different tissues using the Aggregated Cauchy Association Test (ACAT), a novel powerful and robust test for combining association results under general correlation patterns. A total of 403 genes were significantly associated with at least one cancer type using predicted gene expression for at least one tissue (p&lt;0.05/899,338, Bonferroni adjusted for the total number of gene-cancer-tissue combinations tested); 96 additional genes were identified when combining test results across cancers (p&lt;0.05/81,758, adjusting for the total number of gene-tissue combinations tested); and 35 additional genes when further combining test results across tissue (p&lt;0.05/12,001, adjusting for the number of genes tested). Among these significant genes, 70 were not near previously-published GWAS index variants (&gt;250 kb distant). 14 of the 70 novel genes were identified from the single cancer single tissue test; an additional 43 were identified with the cross-cancer test; and another 13 were identified when further combining the results across tissues. The newly identified genes include RBBP8, which regulates cell proliferation and modulates BRCA1-mediated DNA repair, and TP53BP1, which is involved in double-strand break repair. Other newly identified genes are involved in chromatin structure, tumorigenesis, apoptosis, transcriptional regulation, DNA repair, immune system, oxidative damage and cell cycle, proliferation, progression, shape, structure, and migration.

Citation Format: Helian Feng, Arunabha Majumdar, Bogdan Pasaniuc, Hongjie Chen, Sara Lindström, BCAC, OCAC, PRACTICAL, Jeroen Huyghe, Stephanie L. Schmit, Tracy A. O'Mara, Deborah J. Thompson, Stuart MacGregor, Paul Brennan, James McKay, Richard S. Houlston, Beatrice S. Melin, Christopher Amos, Anne E. Cus, Mark M. Iles, Siddhartha Kar, Paul Pharoah, Rayjean J. Hung, Peter Kraft. Cross-cancer cross-tissue transcriptome-wide association study (TWAS) of 11 cancers identifies 56 novel genes [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 30.

Assessment of polygenic architecture and risk prediction based on common variants across fourteen cancers

Genome-wide association studies (GWAS) have led to the identification of hundreds of susceptibility loci across cancers, but the impact of further studies remains uncertain. Here we analyse summary-level data from GWAS of European ancestry across fourteen cancer sites to estimate the number of common susceptibility variants (polygenicity) and underlying effect-size distribution. All cancers show a high degree of polygenicity, involving at a minimum of thousands of loci. We project that sample sizes required to explain 80% of GWAS heritability vary from 60,000 cases for testicular to over 1,000,000 cases for lung cancer. The maximum relative risk achievable for subjects at the 99th risk percentile of underlying polygenic risk scores (PRS), compared to average risk, ranges from 12 for testicular to 2.5 for ovarian cancer. We show that PRS have potential for risk stratification for cancers of breast, colon and prostate, but less so for others because of modest heritability and lower incidence.

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma

Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAF^V600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

There are more than 20 known melanoma susceptibility genes. Here, using a massively parallel reporter assay, the authors identify risk-associated variants that alter gene transcription, and demonstrate that expression of one such gene, MX2, leads to the promotion of melanoma in a zebrafish model.

Neural crest-derived tumor neuroblastoma and melanoma share 1p13.2 as susceptibility locus that shows a long-range interaction with the SLC16A1 gene

Neuroblastoma (NB) and malignant cutaneous melanoma (CMM) are neural crest cells (NCC)-derived tumors and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association studies (GWAS). We took a three-staged approach to conduct cross-disease meta-analysis of GWAS for NB and CMM (2101 NB cases and 4202 controls; 12 874 CMM cases and 23 203 controls) to identify shared loci. Findings were replicated in 1403 NB cases and 1403 controls of European ancestry and in 636 NB, 508 CMM cases and 2066 controls of Italian origin. We found a cross-association at locus 1p13.2 (rs2153977, odds ratio = 0.91, P = 5.36 × 10-8). We also detected a suggestive (P < 10-7) NB-CMM cross-association at 2q37.1 with opposite effect on cancer risk. Pathway analysis of 110 NB-CMM risk loci with P < 10-4 demonstrated enrichment of biological processes such as cell migration, cell cycle, metabolism and immune response, which are essential of human NCC development, underlying both tumors. In vitro and in silico analyses indicated that the rs2153977-T protective allele, located in an NB and CMM enhancer, decreased expression of SLC16A1 via long-range loop formation and altered a T-box protein binding site. Upon depletion of SLC16A1, we observed a decrease of cellular proliferation and invasion in both NB and CMM cell lines, suggesting its role as oncogene. This is the largest study to date examining pleiotropy across two NC cell-derived tumors identifying 1p13.2 as common susceptibility locus for NB and CMM risk. We demonstrate that combining genome-wide association studies results across cancers with same origins can identify new loci common to neuroblastoma and melanoma arising from tissues which originate from neural crest cells. Our results also show 1p13.2 confer risk to neuroblastoma and melanoma by regulating SLC16A1.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Correction to: Overlapping genetic architecture between Parkinson disease and melanoma

The original version of this article unfortunately contained a mistake. Supplementary Tables 3 and 4 are not available with the rest of the supplementary material available online.

Overlapping genetic architecture between Parkinson disease and melanoma

Epidemiologic studies have reported inconsistent results regarding an association between Parkinson disease (PD) and cutaneous melanoma (melanoma). Identifying shared genetic architecture between these diseases can support epidemiologic findings and identify common risk genes and biological pathways. Here, we apply polygenic, linkage disequilibrium-informed methods to the largest available case-control, genome-wide association study summary statistic data for melanoma and PD. We identify positive and significant genetic correlation (correlation: 0.17, 95% CI 0.10-0.24; P = 4.09 × 10-06) between melanoma and PD. We further demonstrate melanoma and PD-inferred gene expression to overlap across tissues (correlation: 0.14, 95% CI 0.06 to 0.22; P = 7.87 × 10-04) and highlight seven genes including PIEZO1, TRAPPC2L, and SOX6 as potential mediators of the genetic correlation between melanoma and PD. These findings demonstrate specific, shared genetic architecture between PD and melanoma that manifests at the level of gene expression.

Does polygenic risk influence associations between sun exposure and melanoma? A prospective cohort analysis

BACKGROUND

Melanoma develops as the result of complex interactions between sun exposure and genetic factors. However, data on these interactions from prospective studies are scant.

OBJECTIVES

To quantify the association between ambient and personal ultraviolet exposure and incident melanoma in a large population-based prospective study of men and women residing in a setting of high ambient ultraviolet radiation, and to examine potential gene-environment interactions.

METHODS

Data were obtained from the QSkin Sun and Health Study, a prospective cohort study of men and women aged 40-69 years, randomly sampled from the Queensland population in 2011. Participants were genotyped and assessed for ultraviolet exposure.

RESULTS

Among participants with genetic data (n = 15 373), 420 (2·7%) developed cutaneous melanoma (173 invasive, 247 in situ) during a median follow-up time of 4·4 years. Country of birth, age at migration, having > 50 sunburns in childhood or adolescence, and a history of keratinocyte cancer or actinic lesions were significantly associated with melanoma risk.

CONCLUSIONS

An interaction with polygenic risk was suggested: among people at low polygenic risk, markers of cumulative sun exposure (as measured by actinic damage) were associated with melanoma. In contrast, among people at high polygenic risk, markers of high-level early-life ambient exposure (as measured by place of birth) were associated with melanoma (hazard ratio for born in Australia vs. overseas 3·16, 95% confidence interval 1·39-7·22). These findings suggest interactions between genotype and environment that are consistent with divergent pathways for melanoma development. What's already known about this topic? The relationship between sun exposure and melanoma is complex, and exposure effects are highly modified by host factors and behaviours. The role of genotype on the relationship between ultraviolet radiation exposure and melanoma risk is poorly understood. What does this study add? We found that country of birth, age at migration, sunburns in childhood or adolescence, and history of keratinocyte cancer or actinic lesions were significantly associated with melanoma risk, while other measures of continuous or more intermittent patterns of sun exposure were not. We found evidence for gene-environment interactions that are consistent with divergent pathways for melanoma development. Linked Comment: Cust. Br J Dermatol 2020; 183:205-206. Plain language summary available online.

Is there a causal relationship between vitamin D and melanoma risk? A Mendelian randomization study

BACKGROUND

Several preclinical studies have identified the antiproliferative effects of 25-hydroxyvitamin D [25(OH)D; vitamin D]. Ultraviolet radiation (UVR) is essential for vitamin D synthesis yet increases the risk of melanoma. Observational studies on the association of vitamin D levels with melanoma risk have reported inconclusive results, and are difficult to interpret owing to the potential confounding from the dual role of UVR.

OBJECTIVES

To determine whether there is a causal association between genetically predicted 25(OH)D concentrations and melanoma using a Mendelian randomization (MR) approach.

METHODS

We performed MR using summary data from a large genome-wide association study (GWAS) meta-analysis of melanoma risk, consisting of 12 874 cases and 23 203 controls. Five single nucleotide polymorphisms associated with 25(OH)D concentration - rs12785878, rs10741657, rs2282679, rs6013897 and rs116970203 - were selected as instrumental variables. An inverse variance weighted method was used to access the evidence for causality. MR results from the melanoma meta-analysis were combined with results from an MR study based on a melanoma risk GWAS using UK Biobank data.

RESULTS

A 20 nmol L^-1 decrease in 25(OH)D was not associated with melanoma risk [odds ratio (OR) 1·06, 95% confidence interval (CI) 0·95-1·19]. Results from the UK Biobank were concordant with this, with meta-analysis of our and UK Biobank-derived MR causal estimates showing no association (OR 1·02, 95% CI 0·92-1·13 for a 20 nmol L^-1 decrease).

CONCLUSIONS

The results suggest that vitamin D levels may not be causally associated with the risk of melanoma. What's already known about this topic? Antitumour activity of vitamin D has been identified in preclinical studies. Observational studies link vitamin D deficiency with an increased risk of a range of cancers. There is a growing public interest for vitamin D supplementation. Observational studies of melanoma are fraught with difficulties because while higher ultraviolet radiation levels increase vitamin D levels, such exposure is also associated with increased melanoma risk. Results from observational studies are inconclusive regarding the effect of vitamin D on melanoma risk. What does this study add? Using Mendelian randomization, an approach to causal inference, which is analogous to a natural randomized controlled trial, we found no causal association between vitamin D levels and melanoma.

Abstract 1592: Genome-wide meta-analysis of keratinocytic cancers identifies 26 novel risk loci

Background: Basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) (collectively referred to as keratinocytic skin cancers (KC)) are the most common forms of cancer in fair-skinned people. KCs account for the highest expenditure of any cancer type in the Australian healthcare system. Knowledge about the genetic architecture of skin cancer may help elucidate the biology of the disease.

Methods: We conducted a genome-wide meta-analysis of KC, including 28248 cases and 353855 controls from United Kingdom, United States and Australian populations. We performed LD score regression to determine the genetic correlation between KC, BCC, SCC, and melanoma.

Results: We identified 41 independent genome-wide significant SNPs associated with KC risk, 16 of which have not been reported previously. A further 10 loci were implicated using gene-based tests. New loci included common variants in BRCA2 (rare variants increase risk of various cancers) and CTLA4 (immune response; antibodies targeting CTLA-4 are used to treat melanoma). We found a strong genetic correlation between BCC and SCC, 0.93 (95% confidence interval [CI] 0.71 - 1.16). The genetic correlation between melanoma and BCC was 0.66 (95% CI: 0.41 - 0.91). Comparatively lower genetic correlations were identified between melanoma and SCC 0.43 (95% CI: 0.12 - 0.74), and melanoma and KC 0.52 (95% CI: 0.27 - 0.76).

Conclusions: We have identified novel genetic loci associated with KC risk which will lead to a greater understanding of the biology of KC. Future work will leverage the genetic correlations between BCC, SCC and melanoma to further advance our understanding of these cancers.

Citation Format: Upekha E. Liyanage, Matthew H. Law, Xikun Han, Jiyuan An, Jue-Sheng Ong, Puya Gharahkhani, Mark M. Iles, Rachel E. Neale, Catherine Olsen, Stuart Macgregor, David C. Whiteman. Genome-wide meta-analysis of keratinocytic cancers identifies 26 novel risk loci [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1592.

Publisher Correction: Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways

The original version of this Article contained errors in the spelling of the authors Fan Liu and M. Arfan Ikram, which were incorrectly given as Fan Lui and Arfan M. Ikram. In addition, the original version of this Article also contained errors in the author affiliations which are detailed in the associated Publisher Correction.

Assessing the Incremental Contribution of Common Genomic Variants to Melanoma Risk Prediction in Two Population-Based Studies

It is unclear to what degree genomic and traditional (phenotypic and environmental) risk factors overlap in their prediction of melanoma risk. We evaluated the incremental contribution of common genomic variants (in pigmentation, nevus, and other pathways) and their overlap with traditional risk factors, using data from two population-based case-control studies from Australia (n = 1,035) and the United Kingdom (n = 1,460) that used the same questionnaires. Polygenic risk scores were derived from 21 gene regions associated with melanoma and odds ratios from published meta-analyses. Logistic regression models were adjusted for age, sex, center, and ancestry. Adding the polygenic risk score to a model with traditional risk factors increased the area under the receiver operating characteristic curve (AUC) by 2.3% (P = 0.003) for Australia and by 2.8% (P = 0.002) for Leeds. Gene variants in the pigmentation pathway, particularly MC1R, were responsible for most of the incremental improvement. In a cross-tabulation of polygenic by traditional tertile risk scores, 59% (Australia) and 49% (Leeds) of participants were categorized in the same (concordant) tertile. Of participants with low traditional risk, 9% (Australia) and 21% (Leeds) had high polygenic risk. Testing of genomic variants can identify people who are susceptible to melanoma despite not having a traditional phenotypic risk profile.

Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways

The total number of acquired melanocytic nevi on the skin is strongly correlated with melanoma risk. Here we report a meta-analysis of 11 nevus GWAS from Australia, Netherlands, UK, and USA comprising 52,506 individuals. We confirm known loci including MTAP, PLA2G6, and IRF4, and detect novel SNPs in KITLG and a region of 9q32. In a bivariate analysis combining the nevus results with a recent melanoma GWAS meta-analysis (12,874 cases, 23,203 controls), SNPs near GPRC5A, CYP1B1, PPARGC1B, HDAC4, FAM208B, DOCK8, and SYNE2 reached global significance, and other loci, including MIR146A and OBFC1, reached a suggestive level. Overall, we conclude that most nevus genes affect melanoma risk (KITLG an exception), while many melanoma risk loci do not alter nevus count. For example, variants in TERC and OBFC1 affect both traits, but other telomere length maintenance genes seem to affect melanoma risk only. Our findings implicate multiple pathways in nevogenesis.

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.

Abstract 234: Understanding melanoma susceptibility through GWAS of risk phenotypes

Melanoma susceptibility SNPs can be categorized by their impact on known susceptibility phenotypes (pigmentation, nevus count, telomere length, or other). In an ongoing meta-analysis of melanoma, at this stage involving 17,800 cases but soon expanding to more than 30,000 cases, 26 SNPs have achieved genome-wide significance. UK Biobank is a population-based cohort of 500,000 UK persons aged at recruitment 40-69 years; participants reported pigmentation characteristics including ease of tanning, natural hair color, skin color (in the absence of tan), and number of childhood sunburns. To conduct analyses, we separately analyzed the red hair phenotype while ordering the other hair colors from black to blond. Among the 26 SNPs from the melanoma GWAS, 15 were associated with at least one pigmentation factor assessed in UK Biobank (each at p &lt;10-5) with 5 being associated with all pigmentation factors (e.g., SLC45A2, ASIP, MC1R) while others were associated with a subset of factors (e.g., RAD23B/TAL2 was strongly associated with hair color alone). 13 of the melanoma SNPs were associated with "ease of tanning." Detailed analysis of "ease of tanning" SNPs found 29 genome-wide significant in a subset of UK Biobank, of which 13 were associated with melanoma risk (p &lt;0.001) but notably 4 of the SNPs showed no evidence of melanoma risk even when being among the strongest effects for "ease of tanning" and the remainder showing weaker evidence of association. For all pigmentation factors (except having red hair), the effect sizes for a SNP on melanoma risk and pigmentation were highly correlated. We examined the 26 melanoma SNPs in recent GWAS of nevus count and telomere length, and similarly observed overlaps with these phenotypes. While the telomere, nevus and pigmentation instruments were largely distinct, we found evidence of some overlap in genetic susceptibility (for instance, the PLA2G6 SNP rs2092180 is strongly associated with nevus count, ease of tanning and skin color [all p &lt;10-12]). Combined analysis of the melanoma GWAS results with the GWAS of each of these phenotypes provides overwhelming evidence of further susceptibility loci, indicating that GWAS of intermediate phenotypes can benefit identification of the underlying genetic profile.

Supported by Cancer Research UK (C588/A19167) & NIH (CA083115). This research has been conducted using the UK Biobank Resource.

Citation Format: David T. Bishop, Mark M. Iles, Julia A. Newton-Bishop, Jennifer H. Barrett, Adam Trouwer, Matthew Law, Stuart MacGregor, David Duffy, Nicholas G. Martin, Nicholas Hayward, Anne Cust, Maria T. Landi, Alisa Goldstein, Jianxin Shi, Mitch Machiela, Florence Demenais, Peter Kanetsky, GenoMEL & MELANOSTRUM. Understanding melanoma susceptibility through GWAS of risk phenotypes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 234.

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.

Publisher correction: Functional characterization of a multi-cancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148

This corrects the article DOI: 10.1038/ncomms15034.

Polyunsaturated fatty acids and risk of melanoma: A Mendelian randomisation analysis

Melanoma is the deadliest form of skin cancer, mainly affecting populations of European ancestry. Some observational studies suggest that particular diets reduce melanoma risk, putatively through an increase in polyunsaturated fatty acid (PUFA) consumption. However, interpretation of these observational findings is difficult due to residual confounding or reverse causality. To date, a randomized controlled trial has not been carried out to examine the relationship between PUFAs and melanoma. Hence, we performed a Mendelian randomisation (MR) study to evaluate the link between PUFAs and melanoma. To perform MR, we used summary results from the largest risk genome-wide association study (GWAS) meta-analysis of melanoma, consisting of 12,874 cases and 23,203 controls. As instrumental variables we selected SNPs associated with PUFA levels from a GWAS meta-analysis of PUFA levels, from the CHARGE consortium. We used the inverse variance weighted method to estimate a causal odds ratio. To aid interpretation, we established a benchmark "large" predicted change in PUFAs in which, for example, an increase in docosahexaenoic acid (DPA) of 0.17 units (equal to 1 standard deviation) moves a person from the 17^th percentile to the median. Raising PUFA levels by a large amount (increasing DPA by 0.17 units) only negligibly changed melanoma risk: odds ratio [OR] = 1.03 (95% confidence interval [CI] = 0.96-1.10). Other PUFAs yielded similar results as DPA. Our MR analysis suggests that the effect of PUFA levels on melanoma risk is either zero or very small.

A common intronic variant of PARP1 confers melanoma risk and mediates melanocyte growth via regulation of MITF

Previous genome-wide association studies have identified a melanoma-associated locus at 1q42.1 that encompasses a ∼100-kb region spanning the PARP1 gene. Expression quantitative trait locus (eQTL) analysis in multiple cell types of the melanocytic lineage consistently demonstrated that the 1q42.1 melanoma risk allele (rs3219090[G]) is correlated with higher PARP1 levels. In silico fine-mapping and functional validation identified a common intronic indel, rs144361550 (-/GGGCCC; r2 = 0.947 with rs3219090), as displaying allele-specific transcriptional activity. A proteomic screen identified RECQL as binding to rs144361550 in an allele-preferential manner. In human primary melanocytes, PARP1 promoted cell proliferation and rescued BRAFV600E-induced senescence phenotypes in a PARylation-independent manner. PARP1 also transformed TERT-immortalized melanocytes expressing BRAFV600E. PARP1-mediated senescence rescue was accompanied by transcriptional activation of the melanocyte-lineage survival oncogene MITF, highlighting a new role for PARP1 in melanomagenesis.