Ulcerated melanoma: Systems biology evidence of inflammatory imbalance towards pro-tumourigenicity.

Microscopic ulceration is an independent predictor of melanoma death. Here, we used systems biology to query the role of host and tumour-specific processes in defining the phenotype. Albumin level as a measure of systemic inflammation was predictive of fewer tumour-infiltrating lymphocytes and poorer survival in the Leeds Melanoma Cohort. Ulcerated melanomas were thicker and more mitotically active (with corresponding transcriptomic upregulated cell cycle pathways). Sequencing identified tumoural p53 and APC mutations, and TUBB2B amplification as associated with the phenotype. Ulcerated tumours had perturbed expression of cytokine genes, consistent with protumourigenic inflammation and histological and transcriptomic evidence for reduced adaptive immune cell infiltration. Pathway/network analysis of multiomic data using neural networks highlighted a role for the ß-catenin pathway in the ulceration, linking genomic changes in the tumour to immunosuppression and cell proliferation. In summary, the data suggest that ulceration is in part associated with genomic changes but that host factors also predict melanoma death with evidence of reduced immune responses to the tumour.

ß-Catenin-mediated immune evasion pathway frequently operates in primary cutaneous melanomas.

Immunotherapy prolongs survival in only a subset of melanoma patients, highlighting the need to better understand the driver tumor microenvironment. We conducted bioinformatic analyses of 703 transcriptomes to probe the immune landscape of primary cutaneous melanomas in a population-ascertained cohort. We identified and validated 6 immunologically distinct subgroups, with the largest having the lowest immune scores and the poorest survival. This poor-prognosis subgroup exhibited expression profiles consistent with ß-catenin-mediated failure to recruit CD141+ DCs. A second subgroup displayed an equally bad prognosis when histopathological factors were adjusted for, while 4 others maintained comparable survival profiles. The 6 subgroups were replicated in The Cancer Genome Atlas (TCGA) melanomas, where ß-catenin signaling was also associated with low immune scores predominantly related to hypomethylation. The survival benefit of high immune scores was strongest in patients with double-WT tumors for BRAF and NRAS, less strong in BRAF-V600 mutants, and absent in NRAS (codons 12, 13, 61) mutants. In summary, we report evidence for a ß-catenin-mediated immune evasion in 42% of melanoma primaries overall and in 73% of those with the worst outcome. We further report evidence for an interaction between oncogenic mutations and host response to melanoma, suggesting that patient stratification will improve immunotherapeutic outcomes.

Independent replication of a melanoma subtype gene signature and evaluation of its prognostic value and biological correlates in a population cohort.

Development and validation of robust molecular biomarkers has so far been limited in melanoma research. In this paper we used a large population-based cohort to replicate two published gene signatures for melanoma classification. We assessed the signatures prognostic value and explored their biological significance by correlating them with factors known to be associated with survival (vitamin D) or etiological routes (nevi, sun sensitivity and telomere length). Genomewide microarray gene expressions were profiled in 300 archived tumors (224 primaries, 76 secondaries). The two gene signatures classified up to 96% of our samples and showed strong correlation with melanoma specific survival (P=3 x 10(-4)), Breslow thickness (P=5 x 10(-10)), ulceration (P=9.x10-8) and mitotic rate (P=3 x 10(-7)), adding prognostic value over AJCC stage (adjusted hazard ratio 1.79, 95%CI 1.13-2.83), as previously reported. Furthermore, molecular subtypes were associated with season-adjusted serum vitamin D at diagnosis (P=0.04) and genetically predicted telomere length (P=0.03). Specifically, molecular high-grade tumors were more frequent in patients with lower vitamin D levels whereas high immune tumors came from patients with predicted shorter telomeres. Our data confirm the utility of molecular biomarkers in melanoma prognostic estimation using tiny archived specimens and shed light on biological mechanisms likely to impact on cancer initiation and progression.

The clinicopathological and gene expression patterns associated with ulceration of primary melanoma.

Ulceration of primary melanomas is associated with poor prognosis yet is reported to predict benefit from adjuvant interferon. To better understand the biological processes involved, clinicopathological factors associated with ulceration were determined in 1804 patients. From this cohort, 348 primary tumor blocks were sampled to generate gene expression data using a 502-gene cancer panel and 195 blocks were used for immunohistochemistry to detect macrophage infiltration and vessel density. Gene expression results were validated using a whole genome array in two independent sample sets. Ulceration of primary melanomas was associated with more proliferative tumors, tumor vessel invasion, and increased microvessel density. Infiltration of tumors with greater number of macrophages and gene expression pathways associated with wound healing and up-regulation of pro-inflammatory cytokines suggests that ulceration is associated with tumor-related inflammation. The relative benefit from interferon reported in patients with ulcerated tumors may reflect modification of signaling pathways involved in inflammation.

MGMT promoter methylation is associated with temozolomide response and prolonged progression-free survival in disseminated cutaneous melanoma.

To investigate the predictive and prognostic value of O(6) -methylguanine DNA methyltransferase (MGMT) inactivation by analyses of promoter methylation in pretreatment tumor biopsies from patients with cutaneous melanoma treated with dacarbazine (DTIC) or temozolomide (TMZ) were performed. The patient cohorts consisted of Belgian and Swedish disseminated melanoma patients. Patients were subdivided into those receiving single-agent treatment with DTIC/TMZ (cohort S, n = 74) and those treated with combination chemotherapy including DTIC/TMZ (cohort C, n = 79). Median follow-up was 248 and 336 days for cohort S and cohort C, respectively. MGMT promoter methylation was assessed by three methods. The methylation-related transcriptional silencing of MGMT mRNA expression was assessed by real-time RT-PCR. Response to chemotherapy and progression-free survival (PFS) and overall survival were correlated to MGMT promoter methylation status. MGMT promoter methylation was detected in tumor biopsies from 21.5 % of the patients. MGMT mRNA was found to be significantly lower in tumors positive for MGMT promoter methylation compared to tumors without methylation in both treatment cohorts (p < 0.005). DTIC/TMZ therapy response rate was found to be significantly associated with MGMT promoter methylation in cohort S (p = 0.0005), but did not reach significance in cohort C (p = 0.16). Significantly longer PFS was observed among patients with MGMT promoter-methylated tumors (p = 0.002). Multivariate Cox regression analysis identified presence of MGMT promoter methylation as an independent variable associated with longer PFS. Together, this implies that MGMT promoter methylation is associated with response to single-agent DTIC/TMZ and longer PFS in disseminated cutaneous melanoma.

Targeting human apurinic/apyrimidinic endonuclease 1 (APE1) in phosphatase and tensin homolog (PTEN) deficient melanoma cells for personalized therapy.

Phosphatase and tensin homolog (PTEN) loss is associated with genomic instability. APE1 is a key player in DNA base excision repair (BER) and an emerging drug target in cancer. We have developed small molecule inhibitors against APE1 repair nuclease activity. In the current study we explored a synthetic lethal relationship between PTEN and APE1 in melanoma. Clinicopathological significance of PTEN mRNA and APE1 mRNA expression was investigated in 191 human melanomas. Preclinically, PTEN-deficient BRAF-mutated (UACC62, HT144, and SKMel28), PTEN-proficient BRAF-wildtype (MeWo), and doxycycline-inducible PTEN-knockout BRAF-wildtype MeWo melanoma cells were DNA repair expression profiled and investigated for synthetic lethality using a panel of four prototypical APE1 inhibitors. In human tumours, low PTEN mRNA and high APE1 mRNA was significantly associated with reduced relapse free and overall survival. Pre-clinically, compared to PTEN-proficient cells, PTEN-deficient cells displayed impaired expression of genes involved in DNA double strand break (DSB) repair. Synthetic lethality in PTEN-deficient cells was evidenced by increased sensitivity, accumulation of DSBs and induction of apoptosis following treatment with APE1 inhibitors. We conclude that PTEN deficiency is not only a promising biomarker in melanoma, but can also be targeted by a synthetic lethality strategy using inhibitors of BER, such as those targeting APE1.

Inherited variation in the PARP1 gene and survival from melanoma.

We report the association of an inherited variant located upstream of the poly(adenosine diphosphate-ribose) polymerase 1 (PARP1) gene (rs2249844), with survival in 11 BioGenoMEL melanoma cohorts. The gene encodes a protein involved in a number of cellular processes including single-strand DNA repair. Survival analysis was conducted for each cohort using proportional hazards regression adjusting for factors known to be associated with survival. Survival was measured as overall survival (OS) and, where available, melanoma-specific survival (MSS). Results were combined using random effects meta-analysis. Evidence for a role of the PARP1 protein in melanoma ulceration and survival was investigated by testing gene expression levels taken from formalin-fixed paraffin-embedded tumors. A significant association was seen for inheritance of the rarer variant of PARP1, rs2249844 with OS (hazard ratio (HR) = 1.16 per allele, 95% confidence interval (CI) 1.04-1.28, p = 0.005, eleven cohorts) and MSS (HR = 1.20 per allele, 95% CI 1.01-1.39, p = 0.03, eight cohorts). We report bioinformatic data supportive of a functional effect for rs2249844. Higher levels of PARP1 gene expression in tumors were shown to be associated with tumor ulceration and poorer OS.