Epigenetic alterations are increasingly implicated in metastasis, whereas very few genetic mutations have been identified as authentic drivers of cancer metastasis. Yet, to date, few studies have identified metastasis-related epigenetic drivers, in part because a framework for identifying driver epigenetic changes in metastasis has not been established. Using reduced representation bisulfite sequencing (RRBS), we mapped genome-wide DNA methylation patterns in three cutaneous primary and metastatic melanoma cell line pairs to identify metastasis-related epigenetic drivers. Globally, metastatic melanoma cell lines were hypomethylated compared to the matched primary melanoma cell lines. Using whole genome RRBS we identified 75 shared (10 hyper- and 65 hypomethylated) differentially methylated fragments (DMFs), which were associated with 68 genes showing significant methylation differences. One gene, Early B Cell Factor 3 (EBF3), exhibited promoter hypermethylation in metastatic cell lines, and was validated with bisulfite sequencing and in two publicly available independent melanoma cohorts (n = 40 and 458 melanomas, respectively). We found that hypermethylation of the EBF3 promoter was associated with increased EBF3 mRNA levels in metastatic melanomas and subsequent inhibition of DNA methylation reduced EBF3 expression. RNAi-mediated knockdown of EBF3 mRNA levels decreased proliferation, migration and invasion in primary and metastatic melanoma cell lines. Overall, we have identified numerous epigenetic changes characterising metastatic melanoma cell lines, including EBF3-induced aggressive phenotypic behaviour with elevated EBF3 expression in metastatic melanoma, suggesting that EBF3 promoter hypermethylation may be a candidate epigenetic driver of metastasis.

SummaryConstitutive expression of the immune checkpoint, PD-L1, inhibits anti-tumor immune responses in cancer, although the factors involved in PD-L1 regulation are poorly understood. Here we show that loss of global DNA methylation, particularly in intergenic regions and repeat elements, is associated with constitutive (PD-L1CON), versus inducible (PD-L1IND), PD-L1 expression in melanoma cell lines. We further show this is accompanied by transcriptomic up-regulation. De novo epigenetic regulators (e.g., DNMT3A) are strongly correlated with PD-L1 expression and methylome status. Accordingly, decitabine-mediated inhibition of global methylation in melanoma cells leads to increased PD-L1 expression. Moreover, viral mimicry and immune response genes are highly expressed in lymphocyte-negative plus PD-L1-positive melanomas, versus PD-L1-negative melanomas in The Cancer Genome Atlas (TCGA). In summary, using integrated genomic analysis we identified that global DNA methylation influences PD-L1 expression in melanoma, and hence melanoma's ability to evade anti-tumor immune responses. These results have implications for combining epigenetic therapy with immunotherapy.

Prognostic markers for glioblastoma multiforme (GBM) are important for patient management. Recent advances have identified prognostic markers for GBMs that use telomerase or the alternative lengthening of telomeres (ALT) mechanism for telomere maintenance. Approximately 40% of GBMs have no defined telomere maintenance mechanism (NDTMM), with a mixed survival for affected individuals. This study examined genetic variants in the cyclin-dependent kinase inhibitor 2A (CDKN2A) gene that encodes the p16INK4a and p14ARF tumor suppressors, and the isocitrate dehydrogenase 1 (IDH1) gene as potential markers of survival for 40 individuals with NDTMM GBMs (telomerase negative and ALT negative by standard assays), 50 individuals with telomerase, and 17 individuals with ALT positive tumors. The analysis of CDKN2A showed NDTMM GBMs had an increased minor allele frequency for the C500G (rs11515) polymorphism compared to those with telomerase and ALT positive GBMs (p = 0.002). Patients with the G500 allele had reduced survival that was independent of age, extent of surgery, and treatment. In the NDTMM group G500 allele carriers had increased loss of CDKN2A gene dosage compared to C500 homozygotes. An analysis of IDH1 mutations showed the R132H mutation was associated with ALT positive tumors, and was largely absent in NDTMM and telomerase positive tumors. In the ALT positive tumors cohort, IDH1 mutations were associated with a younger age for the affected individual. In conclusion, the G500 CDKN2A allele was associated with NDTMM GBMs from older individuals with poorer survival. Mutations in IDH1 were not associated with NDTMM GBMs, and instead were a marker for ALT positive tumors in younger individuals.

Treatment resistance in metastatic melanoma is a longstanding issue. Current targeted therapy regimes in melanoma largely target the proliferating cancer population, leaving slow-cycling cancer cells undamaged. Consequently, slow-cycling cells are enriched upon drug therapy and can remain in the body for years until acquiring proliferative potential that triggers cancer relapse. Here we overview the molecular mechanisms of slow-cycling cells that underlie treatment resistance in melanoma. Three main areas of molecular reprogramming are discussed that mediate slow cycling and treatment resistance. First, a low microphthalmia-associated transcription factor (MITF) dedifferentiated state activates various signaling pathways. This includes WNT5A, EGFR, as well as other signaling activators, such as AXL and NF-κB. Second, the chromatin-remodeling factor Jumonji/ARID domain-containing protein 1B (JARID1B, ) orchestrates and maintains slow cycling and treatment resistance in a small subpopulation of melanoma cells. Finally, a shift in metabolic state toward oxidative phosphorylation has been demonstrated to regulate treatment resistance in slow-cycling cells. Elucidation of the underlying processes of slow cycling and its utilization by melanoma cells may reveal new vulnerable characteristics as therapeutic targets. Moreover, combining current therapies with targeting slow-cycling subpopulations of melanoma cells may allow for more durable and greater treatment responses..

Melanoma is a very aggressive neoplasm with a propensity to undergo progression and invasion early in its evolution. The molecular pathways underpinning invasion in melanoma are now just beginning to be elucidated, but a clear understanding of the transition from non-invasive to invasive melanoma cells remains elusive. Microphthalmia-associated transcription factor (MITF), is thought to be a central player in melanoma biology, and it controls many aspects of the phenotypic expression of the melanocytic lineage. However, recently the paired box transcription factor PAX3 was shown to transcriptionally activate POU3F2/BRN2, leading to direct repression of MITF expression. Here we present a theory to explain melanoma phenotype switching and discuss the predictions that this theory makes. One prediction is that independent and opposing roles for MITF and PAX3 in melanoma would be expected, and we present empirical evidence supporting this: in melanoma tissues PAX3 expression occurs independently of MITF, and PAX3 does not play a key role in melanoma cell proliferation. Furthermore, we show that knockdown of PAX3 inhibits cell migration in a group of “lower MITF” melanoma cell lines, while knockdown of MITF promotes cell migration in a complementary “higher MITF” group of melanoma cell lines. Moreover, the morphological effects of knocking down PAX3 versus MITF in melanoma cells were found to differ. While these data support the notion of independent roles for MITF and PAX3, additional experiments are required to provide robust examination of the proposed genetic switch theory. Only upon clear delineation of the mechanisms associated with progression and invasion of melanoma cells will successful treatments for invasive melanoma be developed.