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.

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.

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.

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.

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

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.

Does genetic risk modify the effect of skin screening on melanoma detection rates?

BACKGROUND: Skin screening is associated with higher melanoma detection rates, a potential indicator of overdiagnosis, but it remains possible that this effect is due to confounding by genetic risk. OBJECTIVES: To compare melanoma incidence among screened vs. unscreened participants within tertiles of genetic risk. METHODS: We investigated melanoma incidence in the QSkin study, a prospective cohort study which for this analysis comprised 15 283 participants aged 40-69 years with genotype data and no prior history of melanoma. We calculated a polygenic score (PGS) for melanoma. We first calculated the age-standardized rate (ASR) of melanoma within PGS tertiles, and then measured the association between skin examination and melanoma detection by calculating the hazard ratio (HR) and 95% confidence interval (95% CI), overall and within PGS tertiles. RESULTS: Melanoma incidence increased with PGS (ASR per 100 000 per year): tertile 1 = 442; tertile 2 = 519; tertile 3 = 871. We found that the HRs for all melanomas (i.e. in situ and invasive) associated with skin examination differed slightly across PGS tertiles [age- and sex-adjusted tertile 1 HR 1.88 (95% CI 1.26-2.81); tertile 2 HR 1.70 (95% CI 1.20-2.41); tertile 3 HR 1.96 (95% CI 1.43-2.70); fully adjusted tertile 1 HR 1.14 (95% CI 0.74-1.75); tertile 2 HR 1.21 (95% CI 0.82-1.78); tertile 3 HR 1.41 (95% CI 1.00-1.98)], but these differences were not statistically significant. HRs for in situ melanoma associated with skin examination were similar across PGS tertiles. For invasive melanomas, the point estimates appeared to be highest in PGS tertile 3 in both the minimally adjusted (age, sex) and fully adjusted models; however, these apparent differences were also not statistically significant. CONCLUSIONS: Genetic risk predicts subsequent melanoma incidence, and is weakly associated with screening behaviour, but it does not explain the higher rate of melanoma detection between screened and unscreened people.

The MC1R r allele does not increase melanoma risk in MITF E318K carriers.

BACKGROUND: Population-wide screening for melanoma is not cost-effective, but genetic characterization could facilitate risk stratification and targeted screening. Common Melanocortin-1 receptor (MC1R) red hair colour (RHC) variants and Microphthalmia-associated transcription factor (MITF) E318K separately confer moderate melanoma susceptibility, but their interactive effects are relatively unexplored. OBJECTIVES: To evaluate whether MC1R genotypes differentially affect melanoma risk in MITF E318K+ vs. E318K- individuals. MATERIALS AND METHODS: Melanoma status (affected or unaffected) and genotype data (MC1R and MITF E318K) were collated from research cohorts (five Australian and two European). In addition, RHC genotypes from E318K+ individuals with and without melanoma were extracted from databases (The Cancer Genome Atlas and Medical Genome Research Bank, respectively). χ2 and logistic regression were used to evaluate RHC allele and genotype frequencies within E318K+/- cohorts depending on melanoma status. Replication analysis was conducted on 200 000 general-population exomes (UK Biobank). RESULTS: The cohort comprised 1165 MITF E318K- and 322 E318K+ individuals. In E318K- cases MC1R R and r alleles increased melanoma risk relative to wild type (wt), P < 0.001 for both. Similarly, each MC1R RHC genotype (R/R, R/r, R/wt, r/r and r/wt) increased melanoma risk relative to wt/wt (P < 0.001 for all). In E318K+ cases, R alleles increased melanoma risk relative to the wt allele [odds ratio (OR) 2.04 (95% confidence interval 1.67-2.49); P = 0.01], while the r allele risk was comparable with the wt allele [OR 0.78 (0.54-1.14) vs. 1.00, respectively]. E318K+ cases with the r/r genotype had a lower but not significant melanoma risk relative to wt/wt [OR 0.52 (0.20-1.38)]. Within the E318K+ cohort, R genotypes (R/R, R/r and R/wt) conferred a significantly higher risk compared with non-R genotypes (r/r, r/wt and wt/wt) (P < 0.001). UK Biobank data supported our findings that r did not increase melanoma risk in E318K+ individuals. CONCLUSIONS: RHC alleles/genotypes modify melanoma risk differently in MITF E318K- and E318K+ individuals. Specifically, although all RHC alleles increase risk relative to wt in E318K- individuals, only MC1R R increases melanoma risk in E318K+ individuals. Importantly, in the E318K+ cohort the MC1R r allele risk is comparable with wt. These findings could inform counselling and management for MITF E318K+ individuals.

Multitrait genetic association analysis identifies 50 new risk loci for gastro-oesophageal reflux, seven new loci for Barrett's oesophagus and provides insights into clinical heterogeneity in reflux diagnosis.

OBJECTIVE: Gastro-oesophageal reflux disease (GERD) has heterogeneous aetiology primarily attributable to its symptom-based definitions. GERD genome-wide association studies (GWASs) have shown strong genetic overlaps with established risk factors such as obesity and depression. We hypothesised that the shared genetic architecture between GERD and these risk factors can be leveraged to (1) identify new GERD and Barrett's oesophagus (BE) risk loci and (2) explore potentially heterogeneous pathways leading to GERD and oesophageal complications. DESIGN: We applied multitrait GWAS models combining GERD (78 707 cases; 288 734 controls) and genetically correlated traits including education attainment, depression and body mass index. We also used multitrait analysis to identify BE risk loci. Top hits were replicated in 23andMe (462 753 GERD cases, 24 099 BE cases, 1 484 025 controls). We additionally dissected the GERD loci into obesity-driven and depression-driven subgroups. These subgroups were investigated to determine how they relate to tissue-specific gene expression and to risk of serious oesophageal disease (BE and/or oesophageal adenocarcinoma, EA). RESULTS: We identified 88 loci associated with GERD, with 59 replicating in 23andMe after multiple testing corrections. Our BE analysis identified seven novel loci. Additionally we showed that only the obesity-driven GERD loci (but not the depression-driven loci) were associated with genes enriched in oesophageal tissues and successfully predicted BE/EA. CONCLUSION: Our multitrait model identified many novel risk loci for GERD and BE. We present strong evidence for a genetic underpinning of disease heterogeneity in GERD and show that GERD loci associated with depressive symptoms are not strong predictors of BE/EA relative to obesity-driven GERD loci.

Genetically determined cutaneous nevi and risk of cancer.

Numerous epidemiologic studies have reported positive associations between higher nevus counts and internal cancers. Whether this association represents a true relationship or is due to bias or confounding by factors associated with both nevus counts and cancer remains unclear. We used germline genetic variants for nevus count to test whether this phenotypic trait is a risk-marker for cancer. We calculated polygenic risk scores (PRS) for nevus counts using individual-level data in the UK Biobank (n = 394 306) and QSkin cohort (n = 17 427). The association between the nevus PRS and each cancer site was assessed using logistic regression adjusted for the effects of age, sex and the first five principal components. In both cohorts, those in the highest nevus PRS quartile had higher risks of melanoma than those in the lowest quartile (UK Biobank odds ratio [OR] 1.42, 95% confidence interval [CI]: 1.29-1.55; QSkin OR 1.58, 95% CI: 1.29-1.94). We also observed increases in risk of basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) associated with higher nevus PRS quartiles (BCC UK Biobank OR 1.38, 95% CI: 1.33-1.44; QSkin OR 1.20, 95% CI: 1.05-1.38 and SCC UK Biobank OR 1.41, 95% CI: 1.28-1.55; QSkin OR 1.44, 95% CI: 1.19-1.77). We found no consistent evidence that nevus count PRS were associated with risks of developing internal cancers. We infer that associations between nevus counts and internal cancers reported in earlier observational studies arose because of unmeasured confounding or other biases.

The emerging field of polygenic risk scores and perspective for use in clinical care.

Genetic testing is used widely for diagnostic, carrier and predictive testing in monogenic diseases. Until recently, there were no genetic testing options available for multifactorial complex diseases like heart disease, diabetes and cancer. Genome-wide association studies (GWAS) have been invaluable in identifying single-nucleotide polymorphisms (SNPs) associated with increased or decreased risk for hundreds of complex disorders. For a given disease, SNPs can be combined to generate a cumulative estimation of risk known as a polygenic risk score (PRS). After years of research, PRSs are increasingly used in clinical settings. In this article, we will review the literature on how both genome-wide and restricted PRSs are developed and the relative merit of each. The validation and evaluation of PRSs will also be discussed, including the recognition that PRS validity is intrinsically linked to the methodological and analytical approach of the foundation GWAS together with the ethnic characteristics of that cohort. Specifically, population differences may affect imputation accuracy, risk magnitude and direction. Even as PRSs are being introduced into clinical practice, there is a push to combine them with clinical and demographic risk factors to develop a holistic disease risk. The existing evidence regarding the clinical utility of PRSs is considered across four different domains: informing population screening programs, guiding therapeutic interventions, refining risk for families at high risk, and facilitating diagnosis and predicting prognostic outcomes. The evidence for clinical utility in relation to five well-studied disorders is summarized. The potential ethical, legal and social implications are also highlighted.

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.

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.

The effect of screening on melanoma incidence and biopsy rates.

BACKGROUND: Cutaneous melanomas are common cancers in white-skinned populations, and early detection is promoted as a means of reducing morbidity and mortality. There is concern that increased skin screening is leading to overdiagnosis of indolent melanomas with low risk of lethality. The extent of melanoma overdiagnosis associated with screening is unknown. OBJECTIVES: To estimate possible overdiagnosis by comparing subsequent melanoma incidence and biopsy rates among people subjected to skin screening those who were not. METHODS: We recruited 43 762 residents of Queensland, Australia, aged 40-69 years, with no prior history of melanoma, selected at random from a population register in 2010. At baseline, participants completed a comprehensive melanoma risk factor survey and were asked if their skin had been examined by a doctor in the 3 years prior to baseline. We calculated incidence and relative risk of histologically confirmed melanoma (invasive and in situ) in years 2-7 of follow-up, obtained through linkage to the cancer registry. In secondary analyses, we measured biopsy rates in years 2-6 of follow-up. We used propensity score analysis to calculate adjusted hazard ratios (aHRs) and 95% confidence intervals (CIs). RESULTS: In total, 28 155 participants underwent skin screening prior to baseline. We observed 967 first-incident melanomas (381 invasive) during 197 191 person-years of follow-up. Those screened had higher rates of melanoma (aHR 1·29, 95% CI 1·02-1·63) and subsequent skin biopses (aHR 1·85, 95% CI 1·69-2·04) than unscreened participants. The higher risk associated with skin screening was evident for in situ melanoma (aHR 1·45, 95% CI 1·09-1·92) but not invasive melanoma (aHR 1·05, 95% CI 0·72-1·54). In secondary analyses, where screening was defined as having a skin biopsy in the first year after baseline, we observed significantly increased risks of melanoma (aHR 1·53, 95% CI 1·23-1·89) and subsequent biopsies (aHR 2·64, 95% CI 2·46-2·84) relative to those who did not have a biopsy. CONCLUSIONS: People who undergo skin screening subsequently experience higher rates of biopsies and melanoma (especially in situ melanoma), even after adjusting for all known risk factors, consistent with overdiagnosis. What is already known about this topic? Cutaneous melanomas are common cancers in white-skinned populations for which early detection is promoted as a means of reducing morbidity and mortality. There is concern that increased surveillance is leading to the overdiagnosis of indolent melanomas that are not destined to be lethal. The extent of melanoma overdiagnosis associated with surveillance is not known. What does this study add? People subjected to skin examinations by a doctor or who undergo skin biopsies subsequently have higher numbers of biopsies and higher rates of melanoma than people not subjected to either, even after adjusting for all known risk factors. These findings suggest that heightened surveillance leads to a proportion of melanomas being diagnosed that otherwise may not have come to clinical attention.

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.

Combined analysis of keratinocyte cancers identifies novel genome-wide loci.

The keratinocyte cancers (KC), basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are the most common cancers in fair-skinned people. KC treatment represents the second highest cancer healthcare expenditure in Australia. Increasing our understanding of the genetic architecture of KC may provide new avenues for prevention and treatment. We first conducted a series of genome-wide association studies (GWAS) of KC across three European ancestry datasets from Australia, Europe and USA, and used linkage disequilibrium (LD) Score regression (LDSC) to estimate their pairwise genetic correlations. We employed a multiple-trait approach to map genes across the combined set of KC GWAS (total N = 47 742 cases, 634 413 controls). We also performed meta-analyses of BCC and SCC separately to identify trait specific loci. We found substantial genetic correlations (generally 0.5-1) between BCC and SCC suggesting overlapping genetic risk variants. The multiple trait combined KC GWAS identified 63 independent genome-wide significant loci, 29 of which were novel. Individual separate meta-analyses of BCC and SCC identified an additional 13 novel loci not found in the combined KC analysis. Three new loci were implicated using gene-based tests. New loci included common variants in BRCA2 (distinct to known rare high penetrance cancer risk variants), and in CTLA4, a target of immunotherapy in melanoma. We found shared and trait specific genetic contributions to BCC and SCC. Considering both, we identified a total of 79 independent risk loci, 45 of which are novel.

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.

Impact of personal genomic risk information on melanoma prevention behaviors and psychological outcomes: a randomized controlled trial.

PURPOSE: We evaluated the impact of personal melanoma genomic risk information on sun-related behaviors and psychological outcomes. METHODS: In this parallel group, open, randomized controlled trial, 1,025 Australians of European ancestry without melanoma and aged 18-69 years were recruited via the Medicare database (3% consent). Participants were randomized to the intervention (n = 513; saliva sample for genetic testing, personalized melanoma risk booklet based on a 40-variant polygenic risk score, telephone-based genetic counseling, educational booklet) or control (n = 512; educational booklet). Wrist-worn ultraviolet (UV) radiation dosimeters (10-day wear) and questionnaires were administered at baseline, 1 month postintervention, and 12 months postbaseline. RESULTS: At 12 months, 948 (92%) participants completed dosimetry and 973 (95%) the questionnaire. For the primary outcome, there was no effect of the genomic risk intervention on objectively measured UV exposure at 12 months, irrespective of traditional risk factors. For secondary outcomes at 12 months, the intervention reduced sunburns (risk ratio: 0.72, 95% confidence interval: 0.54-0.96), and increased skin examinations among women. Melanoma-related worry was reduced. There was no overall impact on general psychological distress. CONCLUSION: Personalized genomic risk information did not influence sun exposure patterns but did improve some skin cancer prevention and early detection behaviors, suggesting it may be useful for precision prevention. There was no evidence of psychological harm.

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.

Vitamin D and overall cancer risk and cancer mortality: a Mendelian randomization study.

There is considerable debate regarding the role that 25-hydroxyvitamin D [25(OH)D] concentrations play in cancer risk or mortality, with earlier studies drawing mixed conclusions. Using data from the UK Biobank (UKB), we evaluate whether genetically predicted 25(OH)D concentrations are associated with overall cancer susceptibility and cancer mortality using five 25(OH)D genetic markers. Data comprised 438 870 white British UKB participants aged 37-73, including 46 155 cancer cases and 6998 cancer-specific deaths. Participants with keratinocyte cancers and/or benign tumors were excluded from the analysis. Odds ratios were calculated per 20 nmol/L increase in genetically predicted 25(OH)D for cancer risk and cancer mortality. For individual cancer risks, estimates were meta-analyzed with publicly available data using a fixed-effect inverse-variance-weighted model. We demonstrated that genetically low plasma 25(OH)D concentrations were not associated with increased cancer risk nor cancer mortality. Stratification by sex or cancer types did not reveal any meaningful differences albeit wider confidence intervals. Fixed-effect meta-analysis of our individual cancer risk estimates with those derived from publicly available cancer consortia data and previous studies further reinforced our null Mendelian randomization findings on prostate, lung, colorectal and breast cancers with tight confidence intervals; for ovarian and pancreatic cancers, our estimates were less precise despite being not statistically significant. Taken altogether, our results provide no genetic evidence for an association between vitamin D and overall cancer outcomes, with tight confidence intervals to exclude all but very small effect sizes.