Search

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.

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 with its risk factors such as sunburn propensity. This genetic correlation can be exploited to identify additional cutaneous melanoma risk loci by multi-trait analysis of GWAS (MTAG). We used bivariate LD-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 MTAG. MTAG 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. 55 of these loci 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 new 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.

Background Previous studies suggest polygenic risk scores (PRS) 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 UK Biobank (UKB; n=395,647; 1,651 cases) and a case-cohort nested within the Melbourne Collaborative Cohort Study (MCCS, Australia; n=4,765; 303 cases). Three PRS were evaluated: 68 SNPs at 54 loci from a 2020 meta-analysis (PRS68); 50 SNPs significant in the 2020 meta-analysis excluding UKB (PRS50); 45 SNPs at 21 loci known in 2018 (PRS45). 10-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 UKB (ratio expected/observed cases E/O=0.65, 95% confidence interval (95%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, e.g. PRS50-adjusted risks E/O=0.91, 95%CI 0.87-0.95. Discriminative ability for PRS68- and PRS50-adjusted absolute risks was higher than for risks based on age and sex alone (ΔAUC 0.07-0.10, p

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.

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.

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.

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.

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.

BackgroundMelanoma 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.ObjectivesWe 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.MethodsWe 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.ResultsWhen 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&nbsp

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

Melanoma develops as the result of complex interactions between sun exposure and genetic factors. Data on these interactions from prospective studies are scant however. We aimed to quantify the association between ambient and personal ultraviolet (UV) exposure and incident melanoma in a large population‐based prospective study of men and women residing in a high ambient UV setting, and to examine potential gene‐environment interactions. 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 greater than 50 sunburns in childhood/adolescence and a history of keratinocyte cancer/actinic lesions were significantly associated with melanoma risk. 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 (HR for born in Australia vs. overseas 3.16, 95% CI 1.39‐7.22). These findings suggest interactions between genotype and environment that are consistent with divergent pathways for melanoma development.

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.

Cancers, including cutaneous melanoma, can cluster in families. In addition to environmental aetiological 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 is unknown. We hypothesise 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 new 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.

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. These are included here.

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.

Meta-analysis of 36,760 cases and 375,188 controls identifies 54 loci associated with susceptibility to cutaneous melanoma. Further analysis combining nevus count and hair color GWAS results provide insights into the genetic architecture of melanoma. 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 x 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.

BACKGROUND: Several pre-clinical studies have identified the anti-proliferative 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 due to the potential confounding from the dual role of UVR. OBJECTIVES: Our objective was 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 SNPs that are associated with 25(OH) D concentration rs12785878, rs10741657, rs2282679, rs6013897 and rs116970203 were selected as instrumental variables (IVs). 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 decrease in 25(OH) D was not associated with melanoma risk (OR = 1.06, 95% CI = 0.94 - 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 20nmol/L decrease). CONCLUSIONS: Our study results suggest that the genetically vitamin D levels may not be causally associated with the risk of melanoma. This article is protected by copyright. All rights reserved.

Genome-wide association studies (GWAS) have identified ∼20 melanoma susceptibility loci. To identify susceptibility genes and variants simultaneously from multiple GWAS loci, we integrated massively-parallel reporter assays (MPRA) with cell type-specific epigenomic data as well as melanocyte-specific expression quantitative trait loci (eQTL) profiling. Starting from 16 melanoma loci, we selected 832 variants overlapping active regions of chromatin in cells of melanocytic lineage and identified 39 candidate functional variants displaying allelic transcriptional activity by MPRA. For four of these loci, we further identified four colocalizing melanocyte cis-eQTL genes (CTSS, CASP8, MX2, and MAFF) matching the allelic activity of MPRA functional variants. Among these, we further characterized the locus encompassing the HIV-1 restriction gene, MX2, on chromosome band Chr21q22.3 and validated a functional variant, rs398206, among multiple high LD variants. rs398206 mediates allelic transcriptional activity via binding of the transcription factor, YY1. This allelic transcriptional regulation is consistent with a significant cis-eQTL of MX2 in primary human melanocytes, where the melanoma risk-associated A allele of rs398206 is correlated with higher MX2 levels. Melanocyte-specific transgenic expression of human MX2 in a zebrafish model demonstrated accelerated melanoma formation in a BRAFV600E background. Thus, using an efficient scalable approach to streamline GWAS follow-up functional studies, we identified multiple candidate melanoma susceptibility genes and variants, and uncovered a pleiotropic function of MX2 in melanoma susceptibility.

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 UK (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, centre and ancestry. Adding the polygenic risk score to a model with traditional risk factors increased the area under 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.

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

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

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.

This corrects the article DOI: 10.1038/ncomms15034

Copyright 2017 American Medical Association. All rights reserved. IMPORTANCE: The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation. OBJECTIVE: To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases. DATA SOURCES: Genomewide association studies (GWAS) published up to January 15, 2015. STUDY SELECTION: GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available. DATA EXTRACTION AND SYNTHESIS: Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population. MAIN OUTCOMES AND MEASURES: Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation. RESULTS: Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]). CONCLUSIONS AND RELEVANCE: It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.

Genome wide association studies (GWAS) have mapped multiple independent cancer susceptibility loci to chr5p15.33. Here, we show that fine-mapping of pancreatic and testicular cancer GWAS within one of these loci (Region 2 in CLPTM1L) focuses the signal to nine highly correlated SNPs. Of these, rs36115365-C associated with increased pancreatic and testicular but decreased lung cancer and melanoma risk, and exhibited preferred protein-binding and enhanced regulatory activity. Transcriptional gene silencing of this regulatory element repressed TERT expression in an allele-specific manner. Proteomic analysis identifies allele-preferred binding of Zinc finger protein 148 (ZNF148) to rs36115365-C, further supported by binding of purified recombinant ZNF148. Knockdown of ZNF148 results in reduced TERT expression, telomerase activity and telomere length. Our results indicate that the association with chr5p15.33-Region 2 may be explained by rs36115365, a variant influencing TERT expression via ZNF148 in a manner consistent with elevated TERT in carriers of the C allele., Genetic variants at multiple loci of chr5p15.33 have been associated with susceptibility to numerous cancers. Here the authors show that the association of one of these loci may be explained by a variant, rs36115365, influencing telomerase reverse transcriptase (TERT) expression via ZNF148.

A shared genetic susceptibility between cutaneous malignant melanoma (CMM) and Parkinson's disease (PD) has been suggested. We investigated this by assessing the contribution of rare variants in genes involved in CMM to PD risk. We studied rare variation across 29 CMM risk genes using high-quality genotype data in 6875 PD cases and 6065 controls and sought to replicate findings using whole-exome sequencing data from a second independent cohort totaling 1255 PD cases and 473 controls. No statistically significant enrichment of rare variants across all genes, per gene, or for any individual variant was detected in either cohort. There were nonsignificant trends toward different carrier frequencies between PD cases and controls, under different inheritance models, in the following CMM risk genes: BAP1, DCC, ERBB4, KIT, MAPK2, MITF, PTEN, and TP53. The very rare TYR p.V275F variant, which is a pathogenic allele for recessive albinism, was more common in PD cases than controls in 3 independent cohorts. Tyrosinase, encoded by TYR, is the rate-limiting enzyme for the production of neuromelanin, and has a role in the production of dopamine. These results suggest a possible role for another gene in the dopamine-biosynthetic pathway in susceptibility to neurodegenerative Parkinsonism, but further studies in larger PD cohorts are needed to accurately determine the role of these genes/variants in disease pathogenesis.

We have genotyped 14,436 nonsynonymous SNPs (nsSNPs) and 897 major histocompatibility complex (MHC) tag SNPs from 1,000 independent cases of ankylosing spondylitis (AS), autoimmune thyroid disease (AITD), multiple sclerosis (MS) and breast cancer (BC). Comparing these data against a common control dataset derived from 1,500 randomly selected healthy British individuals, we report initial association and independent replication in a North American sample of two new loci related to ankylosing spondylitis, ARTS1 and IL23R, and confirmation of the previously reported association of AITD with TSHR and FCRL3. These findings, enabled in part by increased statistical power resulting from the expansion of the control reference group to include individuals from the other disease groups, highlight notable new possibilities for autoimmune regulation and suggest that IL23R may be a common susceptibility factor for the major 'seronegative' diseases.

We performed a multistage genome-wide association study of melanoma. In a discovery cohort of 1804 melanoma cases and 1026 controls, we identified loci at chromosomes 15q13.1 (HERC2/OCA2 region) and 16q24.3 (MC1R) regions that reached genome-wide significance within this study and also found strong evidence for genetic effects on susceptibility to melanoma from markers on chromosome 9p21.3 in the p16/ARF region and on chromosome 1q21.3 (ARNT/LASS2/ANXA9 region). The most significant single-nucleotide polymorphisms (SNPs) in the 15q13.1 locus (rs1129038 and rs12913832) lie within a genomic region that has profound effects on eye and skin color; notably, 50% of variability in eye color is associated with variation in the SNP rs12913832. Because eye and skin colors vary across European populations, we further evaluated the associations of the significant SNPs after carefully adjusting for European substructure. We also evaluated the top 10 most significant SNPs by using data from three other genome-wide scans. Additional in silico data provided replication of the findings from the most significant region on chromosome 1q21.3 rs7412746 (P = 6 × 10(-10)). Together, these data identified several candidate genes for additional studies to identify causal variants predisposing to increased risk for developing melanoma.

© 2015 Nature America, Inc. All rights reserved.Thirteen common susceptibility loci have been reproducibly associated with cutaneous malignant melanoma (CMM). We report the results of an international 2-stage meta-analysis of CMM genome-wide association studies (GWAS). This meta-analysis combines 11 GWAS (5 previously unpublished) and a further three stage 2 data sets, totaling 15,990 CMM cases and 26,409 controls. Five loci not previously associated with CMM risk reached genome-wide significance (P < 5 × 10-8), as did 2 previously reported but unreplicated loci and all 13 established loci. Newly associated SNPs fall within putative melanocyte regulatory elements, and bioinformatic and expression quantitative trait locus (eQTL) data highlight candidate genes in the associated regions, including one involved in telomere biology.

The genetic loci that have been found by genome-wide association studies to modulate risk of coronary heart disease explain only a fraction of its total variance, and gene-gene interactions have been proposed as a potential source of the remaining heritability. Given the potentially large testing burden, we sought to enrich our search space with real interactions by analyzing variants that may be more likely to interact on the basis of two distinct hypotheses: a biological hypothesis, under which MI risk is modulated by interactions between variants that are known to be relevant for its risk factors; and a statistical hypothesis, under which interacting variants individually show weak marginal association with MI. In a discovery sample of 2,967 cases of early-onset myocardial infarction (MI) and 3,075 controls from the MIGen study, we performed pair-wise SNP interaction testing using a logistic regression framework. Despite having reasonable power to detect interaction effects of plausible magnitudes, we observed no statistically significant evidence of interaction under these hypotheses, and no clear consistency between the top results in our discovery sample and those in a large validation sample of 1,766 cases of coronary heart disease and 2,938 controls from the Wellcome Trust Case-Control Consortium. Our results do not support the existence of strong interaction effects as a common risk factor for MI. Within the scope of the hypotheses we have explored, this study places a modest upper limit on the magnitude that epistatic risk effects are likely to have at the population level (odds ratio for MI risk 1.3–2.0, depending on allele frequency and interaction model).

BACKGROUND: We tested whether genetic factors distinctly contribute to either development of coronary atherosclerosis or, specifically, to myocardial infarction in existing coronary atherosclerosis. METHODS: We did two genome-wide association studies (GWAS) with coronary angiographic phenotyping in participants of European ancestry. To identify loci that predispose to angiographic coronary artery disease (CAD), we compared individuals who had this disorder (n=12,393) with those who did not (controls, n=7383). To identify loci that predispose to myocardial infarction, we compared patients who had angiographic CAD and myocardial infarction (n=5783) with those who had angiographic CAD but no myocardial infarction (n=3644). FINDINGS: In the comparison of patients with angiographic CAD versus controls, we identified a novel locus, ADAMTS7 (p=4·98×10(-13)). In the comparison of patients with angiographic CAD who had myocardial infarction versus those with angiographic CAD but no myocardial infarction, we identified a novel association at the ABO locus (p=7·62×10(-9)). The ABO association was attributable to the glycotransferase-deficient enzyme that encodes the ABO blood group O phenotype previously proposed to protect against myocardial infarction. INTERPRETATION: Our findings indicate that specific genetic predispositions promote the development of coronary atherosclerosis whereas others lead to myocardial infarction in the presence of coronary atherosclerosis. The relation to specific CAD phenotypes might modify how novel loci are applied in personalised risk assessment and used in the development of novel therapies for CAD. FUNDING: The PennCath and MedStar studies were supported by the Cardiovascular Institute of the University of Pennsylvania, by the MedStar Health Research Institute at Washington Hospital Center and by a research grant from GlaxoSmithKline. The funding and support for the other cohorts contributing to the paper are described in the webappendix.