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.

No abstract

1a,25-Dihydroxyvitamin D3 signals via the vitamin D receptor (VDR). Higher serum vitamin D is associated with thinner primary melanoma and better outcome, although a causal mechanism has not been established. As patients with melanoma commonly avoid sun exposure, and consequent vitamin D deficiency might worsen outcomes, we interrogated 703 primary melanoma transcriptomes to understand the role of vitamin D-VDR signaling and replicated the findings in The Cancer Genome Atlas metastases. VDR expression was independently protective for melanoma-related death in both primary and metastatic disease. High tumor VDR expression was associated with upregulation of pathways mediating antitumor immunity and corresponding with higher imputed immune cell scores and histologically detected tumor-infiltrating lymphocytes. High VDR-expressing tumors had downregulation of proliferative pathways, notably Wnt/b-catenin signaling. Deleterious low VDR levels resulted from promoter methylation and gene deletion in metastases. Vitamin D deficiency (<25 nmol/L $ 10 ng/mL) shortened survival in primary melanoma in a VDR-dependent manner. In vitro functional validation studies showed that elevated vitamin D-VDR signaling inhibited Wnt/b-catenin signaling genes. Murine melanoma cells overexpressing VDR produced fewer pulmonary metastases than controls in tail-vein metastasis assays. In summary, vitamin D-VDR signaling contributes to controlling pro-proliferative/immunosuppressive Wnt/ b-catenin signaling in melanoma and this is associated with less metastatic disease and stronger host immune responses. This is evidence of a causal relationship between vitamin D-VDR signaling and melanoma survival, which should be explored as a therapeutic target in primary resistance to checkpoint blockade.Significance: VDR expression could potentially be used as a biomarker to stratify patients with melanoma that may respond better to immunotherapy.

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×10−4), Breslow thickness (P=5×10−10), ulceration (P=9.x10−8) and mitotic rate (P=3×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.

Interferon regulatory factor 4 (IRF4) is central to the transcriptional network of activated B-cell-like diffuse large B-cell lymphoma (ABC-DLBCL), an aggressive lymphoma subgroup defined by gene expression profiling. Since cofactor association modifies transcriptional regulatory input by IRF4, we assessed genome occupancy by IRF4 and endogenous cofactors in ABC-DLBCL cell lines. IRF4 partners with SPIB, PU.1 and BATF genome-wide, but SPIB provides the dominant IRF4 partner in this context. Upon SPIB knockdown IRF4 occupancy is depleted and neither PU.1 nor BATF acutely compensates. Integration with ENCODE data from lymphoblastoid cell line GM12878, demonstrates that IRF4 adopts either SPIB- or BATF-centric genome-wide distributions in related states of post-germinal centre B-cell transformation. In primary DLBCL high-SPIB and low-BATF or the reciprocal low-SPIB and high-BATF mRNA expression links to differential gene expression profiles across nine data sets, identifying distinct associations with SPIB occupancy, signatures of B-cell differentiation stage and potential pathogenetic mechanisms. In a population-based patient cohort, SPIBhigh/BATFlow-ABC-DLBCL is enriched for mutation of MYD88, and SPIBhigh/BATFlow-ABC-DLBCL with MYD88-L265P mutation identifies a small subgroup of patients among this otherwise aggressive disease subgroup with distinct favourable outcome. We conclude that differential expression of IRF4 cofactors SPIB and BATF identifies biologically and clinically significant heterogeneity among ABC-DLBCL.

The immune response to melanoma improves survival in untreated patients and predicts response to immune checkpoint blockade. Here we report genetic and environmental predictors of the immune response in a large primary cutaneous melanoma cohort. Bioinformatic analysis of 703 tumor transcriptomes was used to infer immune cell infiltration and categorize tumors into immune subgroups, which were then investigated for association with biological pathways, clinicopathological factors, and copy number alterations. Three subgroups, with "Low”, “Intermediate", and "High" immune signals were identified in primary and replicated in metastatic tumors. Genes in the Low subgroup were enriched for cell cycle and metabolic pathways, whereas genes in the High subgroup were enriched for interferon and NF-κB signaling. We identified high MYC expression partially driven by amplification, HLA-B downregulation, and deletion of IFN-γ and NF-κB pathway genes as regulators of immune suppression. Furthermore, we showed that cigarette smoking, a globally detrimental environmental factor, modulates immunity, reducing survival primarily in patients with a strong immune response. Together these analyses identify a set of easily assessible factors that may serve as predictors of response to immunotherapy in melanoma patients.

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.

Purpose: Previously identified transcriptomic signatures have been based on primary and metastatic melanomas with relatively few American Joint Committee on Cancer (AJCC) stage I tumors, given difficulties in sampling small tumors. The advent of adjuvant therapies has highlighted the need for better prognostic and predictive biomarkers, especially for AJCC stage I and stage II disease.Experimental Design: A total of 687 primary melanoma transcriptomes were generated from the Leeds Melanoma Cohort (LMC). The prognostic value of existing signatures across all the AJCC stages was tested. Unsupervised clustering was performed, and the prognostic value of the resultant signature was compared with that of sentinel node biopsy (SNB) and tested as a biomarker in three published immunotherapy datasets.Results: Previous Lund and The Cancer Genome Atlas signatures predicted outcome in the LMC dataset (P ¼ 10 À8 to 10 À4 ) but showed a significant interaction with AJCC stage (P ¼ 0.04) and did not predict outcome in stage I tumors (P ¼ 0.3-0.7). Consensus-based classification of the LMC dataset identified six classes that predicted outcome, notably in stage I disease. LMC class was a similar indicator of prognosis when compared with SNB, and it added prognostic value to the genes reported by Gerami and colleagues. One particular LMC class consistently predicted poor outcome in patients receiving immunotherapy in two of three tested datasets. Biological characterization of this class revealed high JUN and AXL expression and evidence of epithelial-to-mesenchymal transition.Conclusions: A transcriptomic signature of primary melanoma was identified with prognostic value, including in stage I melanoma and in patients undergoing immunotherapy.