BACKGROUND: It is unclear whether genetic variants affecting vitamin D metabolism are associated with melanoma prognosis. Two functional missense variants in the vitamin D-binding protein gene (GC), rs7041 and rs4588, determine 3 common haplotypes, Gc1s, Gc1f, and Gc2, of which Gc1f may be associated with decreased all-cause death among melanoma patients based on results of a prior study, but the association of Gc1f with melanoma-specific death is unclear. METHODS: We investigated the association of the Gc1s, Gc1f, and Gc2 haplotypes with melanoma-specific and all-cause death among 4490 individuals with incident, invasive primary melanoma in 2 population-based studies using multivariable Cox-proportional hazards regression. RESULTS: In the pooled analysis of both datasets, the patients with the Gc1f haplotype had a 37% lower risk of melanoma-specific death than the patients without Gc1f (hazard ratio [HR] = 0.63, 95% confidence interval [CI] = 0.47 to 0.83, P = .001), with adjustments for age, sex, study center, first- or higher-order primary melanoma, tumor site, pigmentary phenotypes, and Breslow thickness. Associations were similar in both studies. In pooled analyses stratified by Breslow thickness, the corresponding melanoma-specific death HRs for those patients with the Gc1f haplotype compared with those without Gc1f were 0.89 (95% CI = 0.63 to 1.27) among participants with tumor Breslow thickness equal to or less than 2.0 mm and 0.40 (95% CI = 0.25 to 0.63) among participants with tumor Breslow thickness greater than 2.0 mm (Pinteraction = .003). CONCLUSIONS: Our findings suggest that individuals with the GC haplotype Gc1f may have a lower risk of dying from melanoma-specifically from thicker, higher-risk melanoma-than individuals without this Gc1f haplotype.

MDM2-SNP309 (rs2279744), a common genetic modifier of cancer incidence in Li-Fraumeni syndrome, modifies risk, age of onset, or prognosis in a variety of cancers. Melanoma incidence and outcomes vary by sex, and although SNP309 exerts an effect on the estrogen receptor, no consensus exists on its effect on melanoma. MDM2 and MDM4 restrain p53-mediated tumor suppression, independently or together. We investigated SNP309, an a priori MDM4-rs4245739, and two coinherited variants, in a population-based cohort of 3663 primary incident melanomas. Per-allele and per-haplotype (MDM2_SNP309-SNP285; MDM4_rs4245739-rs1563828) odds ratios (OR) for multiple-melanoma were estimated with logistic regression models. Hazard ratios (HR) for melanoma death were estimated with Cox proportional hazards models. In analyses adjusted for covariates, females carrying MDM4-rs4245739*C were more likely to develop multiple melanomas (OR_per-allele = 1.25, 95% CI 1.03-1.51, and P_trend = 0.03), while MDM2-rs2279744*G was inversely associated with melanoma-death (HR_per-allele = 0.63, 95% CI 0.42-0.95, and P_trend = 0.03). We identified 16 coinherited expression quantitative loci that control the expression of MDM2, MDM4, and other genes in the skin, brain, and lungs. Our results suggest that MDM4/MDM2 variants are associated with the development of subsequent primaries and with the death of melanoma in a sex-dependent manner. Further investigations of the complex MDM2/MDM4 motif, and its contribution to the tumor microenvironment and observed associations, are warranted.

Genome-wide association studies (GWAS) and candidate pathway studies have identified low-penetrant genetic variants associated with cutaneous melanoma. We investigated the association of melanoma-risk variants with primary melanoma tumor prognostic characteristics and melanoma-specific survival. The Genes, Environment, and Melanoma Study enrolled 3285 European origin participants with incident invasive primary melanoma. For each of 47 melanoma-risk single nucleotide polymorphisms (SNPs), we used linear and logistic regression modeling to estimate, respectively, the per allele mean changes in log of Breslow thickness and odds ratios for presence of ulceration, mitoses, and tumor-infiltrating lymphocytes (TILs). We also used Cox proportional hazards regression modeling to estimate the per allele hazard ratios for melanoma-specific survival. Passing the false discovery threshold (p = 0.0026) were associations of IRF4 rs12203592 and CCND1 rs1485993 with log of Breslow thickness, and association of TERT rs2242652 with presence of mitoses. IRF4 rs12203592 also had nominal associations (p < 0.05) with presence of mitoses and melanoma-specific survival, as well as a borderline association (p = 0.07) with ulceration. CCND1 rs1485993 also had a borderline association with presence of mitoses (p = 0.06). MX2 rs45430 had nominal associations with log of Breslow thickness, presence of mitoses, and melanoma-specific survival. Our study indicates that further research investigating the associations of these genetic variants with underlying biologic pathways related to tumor progression is warranted.

Evidence on the relationship between the vitamin D pathway and outcomes in melanoma is growing, although it is not always clear. We investigated the impact of measured levels of sun exposure at diagnosis on associations of vitamin D receptor gene (VDR) polymorphisms and melanoma death in 3336 incident primary melanoma cases. Interactions between six SNPs and a common 3'-end haplotype were significant (p < .05). These SNPs, and a haplotype, had a statistically significant association with survival among subjects exposed to high UVB in multivariable regression models and exerted their effect in the opposite direction among those with low UVB. SNPs rs1544410/BsmI and rs731236/TaqI remained significant after adjustment for multiple testing. These results suggest that the association between VDR and melanoma-specific survival is modified by sun exposure around diagnosis, and require validation in an independent study. Whether the observed effects are dependent or independent of vitamin D activation remains to be determined.

Background

Genome-wide association studies have reported that genetic variation at ANRIL (CDKN2B-AS1) is associated with risk of several chronic diseases including coronary artery disease, coronary artery calcification, myocardial infarction, and type 2 diabetes mellitus. ANRIL is located at the CDKN2A/B locus, which encodes multiple melanoma tumor suppressors. We investigated the association of these variants with melanoma prognostic characteristics.

Methods

The Genes, Environment, and Melanoma Study enrolled 3,285 European origin participants with incident invasive primary melanoma. For each of ten disease-associated SNPs at or near ANRIL, we used linear and logistic regression modeling to estimate, respectively, the per allele mean changes in log of Breslow thickness and ORs for presence of ulceration and tumor-infiltrating lymphocytes (TIL). We also assessed effect modification by tumor NRAS/BRAF mutational status.

Results

Rs518394, rs10965215, and rs564398 passed false discovery and were each associated (P ≤ 0.005) with TILs, although only rs564398 was independently associated (P = 0.0005) with TILs. Stratified by NRAS/BRAF mutational status, rs564398*A was significantly positively associated with TILs among NRAS/BRAF mutant, but not wild-type, cases. We did not find SNP associations with Breslow thickness or ulceration.

Conclusions

ANRIL rs564398 was associated with TIL presence in primary melanomas, and this association may be limited to NRAS/BRAF-mutant cases.

Impact

Pathways related to ANRIL variants warrant exploration in relationship to TILs in melanoma, especially given the impact of TILs on immunotherapy and survival.

In an effort to understand the difference between melanomas diagnosed in Australia (New South Wales) and Canada, where the incidence in New South Wales is almost three times greater than in Canada, and mortality is twice as high although survival is slightly more favorable, we had one pathologist review 1,271 melanomas from British Columbia and Ontario, Canada, to compare these to melanomas in New South Wales, Australia. We hypothesized that histopathologic characteristics might provide insight into divergent pathways to melanoma development. We found a number of differences in risk factors and tumor characteristics between the two geographic areas. There were higher mole counts and darker phenotypes in the Canadian patients, while the Australian patients had greater solar elastosis, more lentigo maligna melanomas, and more tumor infiltrating lymphocytes. We hypothesize that the differences observed may illustrate different etiologies – the cumulative exposure pathway among Australian patients and the nevus pathway among Canadian patients. This is one of the largest studies investigating the divergent pathway hypothesis and is particularly robust due to the evaluation of all lesions by one dermatopathologist.

Importance

We previously reported that survival is poorer from histopathologically amelanotic than pigmented melanoma because of more advanced stage at diagnosis. Identifying patients at risk of amelanotic melanoma might enable earlier diagnosis and improved survival; however, the phenotypic characteristics and underlying genetics associated with amelanotic melanoma are unknown.

Objective

To determine whether phenotypic characteristics, carriage of MC1R variants, and history of amelanotic melanoma are associated with histopathologically amelanotic melanoma.

Design, setting, and participants

The Genes, Environment, and Melanoma (GEM) study is an international cohort study that enrolled patients with incident primary cutaneous melanomas from population-based and hospital-based cancer registries (1998 to 2003). The GEM participants included here were 2387 patients with data for phenotypes, MC1R genotype, and primary melanomas scored for histopathologic pigmentation. Of these 2387 patients with incident melanomas scored for pigmentation, 527 had prior primary melanomas also scored for pigmentation.

Main outcomes and measures

Associations of phenotypic characteristics (freckles, nevi, phenotypic index) and MC1R status with incident amelanotic melanomas were evaluated using logistic regression models adjusted for age, sex, study center, and primary status (single or multiple primary melanoma); odds ratios (ORs) and 95% CIs are reported. Association of histopathologic pigmentation between incident and prior melanomas was analyzed using an exact logistic regression model.

Results

This study included 2387 patients (1065 women, 1322 men; mean [SD] age at diagnosis, 58.3 [16.1] years) and 2917 primary melanomas. In a multivariable model including phenotypic characteristics, absence of back nevi, presence of many freckles, and a sun-sensitive phenotypic index were independently associated with amelanotic melanoma. Carriage of MC1R variants was associated with amelanotic melanoma but lost statistical significance in a model with phenotype. Further, patients with incident primary amelanotic melanomas were more likely to have had a prior primary amelanotic melanoma (OR, 4.62; 95% CI, 1.25-14.13) than those with incident primary pigmented melanomas.

Conclusions and relevance

Absence of back nevi, presence of many freckles, a sun-sensitive phenotypic index, and prior amelanotic melanoma increase odds for development of amelanotic melanoma. An increased index of suspicion for melanoma in presenting nonpigmented lesions and more careful examination for signs of amelanotic melanoma during periodic skin examination in patients at increased odds of amelanotic melanoma might lead to earlier diagnosis and improved survival.

Autophagy has been linked with melanoma risk and survival, but no polymorphisms in autophagy-related (ATG) genes have been investigated in relation to melanoma progression. We examined five single-nucleotide polymorphisms (SNPs) in three ATG genes (ATG5; ATG10; and ATG16L) with known or suspected impact on autophagic flux in an international population-based case-control study of melanoma. DNA from 911 melanoma patients was genotyped. An association was identified between (GG) (rs2241880) and earlier stage at diagnosis (OR 0.47; 95% Confidence Intervals (CI) = 0.27-0.81, P = 0.02) and a decrease in Breslow thickness (P = 0.03). The ATG16L heterozygous genotype (AG) (rs2241880) was associated with younger age at diagnosis (P = 0.02). Two SNPs in ATG5 were found to be associated with increased stage (rs2245214 CG, OR 1.47; 95% CI = 1.11-1.94, P = 0.03; rs510432 CC, OR 1.84; 95% CI = 1.12-3.02, P = 0.05). Finally, we identified inverse associations between ATG5 (GG rs2245214) and melanomas on the scalp or neck (OR 0.20, 95% CI = 0.05-0.86, P = 0.03); ATG10 (CC) (rs1864182) and brisk tumor infiltrating lymphocytes (TILs) (OR 0.42; 95% CI = 0.21-0.88, P = 0.02), and ATG5 (CC) (rs510432) with nonbrisk TILs (OR 0.55; 95% CI = 0.34-0.87, P = 0.01). Our data suggest that ATG SNPs might be differentially associated with specific host and tumor characteristics including age at diagnosis, TILs, and stage. These associations may be critical to understanding the role of autophagy in cancer, and further investigation will help characterize the contribution of these variants to melanoma progression.

Associations of MC1R with BRAF mutations in melanoma have been inconsistent between studies. We sought to determine for 1,227 participants in the international population-based Genes, Environment, and Melanoma (GEM) study whether MC1R and phenotypes were associated with melanoma BRAF/NRAS subtypes. We used logistic regression adjusted by age, sex, and study design features and examined effect modifications. BRAF⁺ were associated with younger age, blond/light brown hair, increased nevi, and less freckling, and NRAS⁺ with older age relative to the wild type (BRAF^-/NRAS^-) melanomas (all P < 0.05). Comparing specific BRAF subtypes to the wild type, BRAF V600E was associated with younger age, blond/light brown hair, and increased nevi and V600K with increased nevi and less freckling (all P < 0.05). MC1R was positively associated with BRAF V600E cases but only among individuals with darker phototypes or darker hair (P_interaction < 0.05) but inversely associated with BRAF V600K (P_trend = 0.006) with no significant effect modification by phenotypes. These results support distinct etiologies for BRAF V600E, BRAF V600K, NRAS⁺, and wild-type melanomas. MC1R's associations with BRAF V600E cases limited to individuals with darker phenotypes indicate that MC1R genotypes specifically provide information about BRAF V600E melanoma risk in those not considered high risk based on phenotype. Our results also suggest that melanin pathways deserve further study in BRAF V600E melanomagenesis.