The DNA methylation landscape of human melanoma

DNA methylation variations of DNA damage response correlate survival and local immune status in melanomas

AIM

We aimed to explore the impact of DNA methylation alterations on the DNA damage response (DDR) in melanoma prognosis and immunity. MATERIAL & METHODS: Different melanoma cohorts with molecular and clinical data were included.

RESULTS

Hierarchical clustering utilizing different combinations of DDR-relevant CpGs yielded distinct melanoma subtypes, which were characteristic of different prognoses, transcriptional function profiles of DDR, and immunity and immunotherapy responses but were associated with similar tumor mutation burdens. We then constructed and validated a clinically applicable 4-CpG risk-score signature for predicting survival and immunotherapy response.

CONCLUSION

Our study describes the close interrelationship among DNA methylation, DDR machinery, local tumor immune status, melanoma prognosis, and immunotherapy response.

The relationship of MITF gene expression and promoter methylation with plumage colour in quail

1. This study focused on the relationship between MITF mRNA expression and plumage colour in quail and the effect of promoter methylation on the expression of MITF mRNA.2. The CDS region of MITF mRNA was cloned by RT-PCR, followed by DNA sequencing. The RT-qPCR method was used to analyse the expression levels of MITF mRNA in dorsal skin tissue in Korean quail and Beijing white quail. The promoter region of the MITF gene was cloned, and the CpG island was predicted by the CpGplot program. The methylation levels of the CpG island were analysed using BS-PCR technology.3. Quail MITF mRNA contains a 1,476 bp complete ORF, which encodes a 492 amino acid residue protein. The MITF protein has no signal peptide or transmembrane region. The expression of MITF mRNA in dorsal tissue of Korean quail was significantly higher than that in Beijing white quail (p < 0.01). Abundant cis-elements and a 346 bp CpG island were found in the promoter region of the MITF gene. The average methylation level of the CpG island was 22 (22%) in Korean quail, and 46 (30%) in Beijing white quail (p < 0.05).4. The hypermethylation of the MITF gene promoter region in Beijing white quail resulted in a decrease in expression level, which was related to white feather colour.

5-Methylcytosine immunohistochemistry for predicting cutaneous melanoma prognosis

There is a correlation between DNA methylation and the diseased stage and poor survival. 5-methylcytosine (5-mC) is one of the epigenetic modifications of bases that researchers focus on. Staining with 5-mC immunohistochemistry was used to examine pathological samples taken from individuals diagnosed with cutaneous melanoma. Between Breslow levels 2 and 4, there was a significant difference in the H-score of 5-mC expression (p = 0.046). A significant reduction in 5-mC expression H-scores was seen in patients who were diagnosed with ulcers (p = 0.039). It was shown that patients with low 5-mC had a significantly worse overall survival rate (p = 0.027).

RhoA activation promotes glucose uptake to elevate proliferation in MAPK inhibitor resistant melanoma cells

Cutaneous melanomas harboring a B-RafV600E mutation are treated with immune check point inhibitors or kinase inhibitor combination therapies relying on MAPK inhibitors (MAPKi) Dabrafenib and Trametinib (Curti and Faries, 2021). However, cells become resistant to treatments over the timespan of a few months. Resistance to MAPKi has been associated with adoption of an aggressive amoeboid phenotype characterized by elevated RhoA signaling, enhanced contractility and thick cortical filamentous actin (F-actin) structures (Kim et al., 2016; Misek et al., 2020). Targeting active RhoA through Rho-kinase (ROCK) inhibitors, either alone or in combination with immunotherapies, reverts MAPKi-resistance (Misek et al., 2020; Orgaz et al., 2020). Yet, the mechanisms for this behavior remain largely unknown. Given our recent findings of cytoskeleton's role in cancer cell proliferation (Mohan et al., 2019), survival (Weems et al., 2023), and metabolism (Park et al., 2020), we explored possibilities by which RhoA-driven changes in cytoskeleton structure may confer resistance. We confirmed elevated activation of RhoA in a panel of MAPKi-resistant melanoma cell lines, leading to a marked increase in the presence of contractile F-actin bundles. Moreover, these cells had increased glucose uptake and glycolysis, a phenotype disrupted by pharmacological perturbation of ROCK. However, glycolysis was unaffected by disruption of F-actin bundles, indicating that glycolytic stimulation in MAPKi-resistant melanoma is independent of F-actin organization. Instead, our findings highlight a mechanism in which elevated RhoA signaling activates ROCK, leading to the activation of insulin receptor substrate 1 (IRS1) and P85 of the PI3K pathway, which promotes cell surface expression of GLUT1 and elevated glucose uptake. Application of ROCK inhibitor GSK269962A results in reduced glucose uptake and glycolysis, thus impeding cell proliferation. Our study adds a mechanism to the proposed use of ROCK inhibitors for long-term treatments on MAPKi-resistant melanomas.

Diagnostic Algorithm to Subclassify Atypical Spitzoid Tumors in Low and High Risk According to Their Methylation Status

Current diagnostic algorithms are insufficient for the optimal clinical and therapeutic management of cutaneous spitzoid tumors, particularly atypical spitzoid tumors (AST). Therefore, it is crucial to identify new markers that allow for reliable and reproducible diagnostic assessment and can also be used as a predictive tool to anticipate the individual malignant potential of each patient, leading to tailored individual therapy. Using Reduced Representation Bisulfite Sequencing (RRBS), we studied genome-wide methylation profiles of a series of Spitz nevi (SN), spitzoid melanoma (SM), and AST. We established a diagnostic algorithm based on the methylation status of seven cg sites located in TETK4P2 (Tektin 4 Pseudogene 2), MYO1D (Myosin ID), and PMF1-BGLAP (PMF1-BGLAP Readthrough), which allows the distinction between SN and SM but is also capable of subclassifying AST according to their similarity to the methylation levels of Spitz nevi or spitzoid melanoma. Thus, our epigenetic algorithm can predict the risk level of AST and predict its potential clinical outcomes.

Functional analysis of recurrent CDC20 promoter variants in human melanoma

Small nucleotide variants in non-coding regions of the genome can alter transcriptional regulation, leading to changes in gene expression which can activate oncogenic gene regulatory networks. Melanoma is heavily burdened by non-coding variants, representing over 99% of total genetic variation, including the well-characterized TERT promoter mutation. However, the compendium of regulatory non-coding variants is likely still functionally under-characterized. We developed a pipeline to identify hotspots, i.e. recurrently mutated regions, in melanoma containing putatively functional non-coding somatic variants that are located within predicted melanoma-specific regulatory regions. We identified hundreds of statistically significant hotspots, including the hotspot containing the TERT promoter variants, and focused on a hotspot in the promoter of CDC20. We found that variants in the promoter of CDC20, which putatively disrupt an ETS motif, lead to lower transcriptional activity in reporter assays. Using CRISPR/Cas9, we generated an indel in the CDC20 promoter in human A375 melanoma cell lines and observed decreased expression of CDC20, changes in migration capabilities, increased growth of xenografts, and an altered transcriptional state previously associated with a more proliferative and less migratory state. Overall, our analysis prioritized several recurrent functional non-coding variants that, through downregulation of CDC20, led to perturbation of key melanoma phenotypes.

Mouse methylation profiles for leukocyte cell types, and estimation of leukocyte fractions in inflamed gastrointestinal DNA samples

Precise analysis of tissue DNA and RNA samples is often hampered by contaminating non-target cells whose amounts are highly variable. DNA methylation profiles are specific to cell types, and can be utilized for assessment of the fraction of such contaminating non-target cells. Here, we aimed 1) to identify methylation profiles specific to multiple types of mouse leukocytes, and 2) to estimate the fraction of leukocytes infiltrating inflamed tissues using DNA samples. First, genome-wide DNA methylation analysis was conducted for three myeloid-lineage cells and four lymphoid-lineage cells isolated by fluorescence-activated cell sorting after magnetic-activated cell sorting from leukocytes in the spleen. Clustering analysis using CpG sites within enhancers separated the three myeloid-lineage cells and four lymphoid-lineage cells while that using promoter CpG islands (TSS200CGIs) did not. Among the 266,108 CpG sites analyzed, one CpG site was specifically hypermethylated (β value ≥ 0.7) in B cells, and four, seven, 183, and 34 CpG sites were specifically hypomethylated (β value < 0.2) in CD4+ T cells, CD8+ T cells, B cells, and NK cells, respectively. Importantly, cell type-specific hypomethylated CpG sites were located at genes involved in cell type-specific biological functions. Then, marker CpG sites to estimate the leukocyte fraction in a tissue with leukocyte infiltration were selected, and an estimation algorithm was established. The fractions of infiltrating leukocytes were estimated to be 1.6-12.4% in the stomach (n = 10) with Helicobacter pylori-induced inflammation and 1.5-4.3% in the colon with dextran sulfate sodium-induced colitis (n = 4), and the fractions were highly correlated with those estimated histologically using Cd45-stained tissue sections [R = 0.811 (p = 0.004)]. These results showed that mouse methylation profiles at CpG sites within enhancers reflected leukocyte cell lineages, and the use of marker CpG sites successfully estimated the leukocyte fraction in inflamed gastric and colon tissues.

Cytosine epigenetic modifications and conformational changes in G-quadruplex DNA: An ultraviolet resonance Raman spectroscopy study

Epigenetic modifications of DNA are known to play important regulatory roles in biological systems, especially in regulation of gene expression, and are associated with many types of human diseases, including cancer. Alternative DNA secondary structures, such as G-quadruplexes, can also influence gene transcription, thus suggesting that such structures may represent a distinctive layer of epigenetic information. G-quadruplex structures and DNA epigenetic modifications often go side by side, and recent evidence reveals that cytosine modifications within loops of G-quadruplexes can play a role in modulating their stability and structural polymorphism. Therefore, the development and validation of experimental techniques that can easily and reliably analyse G-quadruplex structures are highly desirable. In the present study, we propose to exploit the advantages of UV resonance Raman (UVRR) spectroscopy to investigate cytosine epigenetic modifications along with conformational changes in G-quadruplex-forming DNA. Our findings show that clear and specific spectral changes occur when there is a change in a G-quadruplex structure. Moreover, UVRR spectral analysis can indirectly distinguish the spectral variations occurring because of modifications in the guanine glycosidic conformations, as well as detect changes in the loops induced by H-bond formation or hydration of nitrogenous bases. The results further underscore the utility of UVRR spectroscopy for G-quadruplex structure elucidation under biologically relevant solution conditions.

Lineage-coupled clonal capture identifies clonal evolution mechanisms and vulnerabilities of BRAFV600E inhibition resistance in melanoma

Targeted cancer therapies have revolutionized treatment but their efficacies are limited by the development of resistance driven by clonal evolution within tumors. We developed "CAPTURE", a single-cell barcoding approach to comprehensively trace clonal dynamics and capture live lineage-coupled resistant cells for in-depth multi-omics analysis and functional exploration. We demonstrate that heterogeneous clones, either preexisting or emerging from drug-tolerant persister cells, dominated resistance to vemurafenib in BRAFV600E melanoma. Further integrative studies uncovered diverse resistance mechanisms. This includes a previously unrecognized and clinically relevant mechanism, chromosome 18q21 gain, which leads to vulnerability of the cells to BCL2 inhibitor. We also identified targetable common dependencies of captured resistant clones, such as oxidative phosphorylation and E2F pathways. Our study provides new therapeutic insights into overcoming therapy resistance in BRAFV600E melanoma and presents a platform for exploring clonal evolution dynamics and vulnerabilities that can be applied to study treatment resistance in other cancers.

UVA Radiation, DNA Damage, and Melanoma

Melanoma is a lethal type of skin tumor that has been linked with sunlight exposure chiefly in fair-skinned human populations. Wavelengths from the sun that can reach the earth's surface include UVA radiation (320-400 nm) and UVB radiation (280-320 nm). UVB effectively induces the formation of dimeric DNA photoproducts, preferentially the cyclobutane pyrimidine dimers (CPDs). The characteristic UVB signature mutations in the form of C to T mutations at dipyrimidine sequences are prevalent in melanoma tumor genomes and have been ascribed to deamination of cytosines within CPDs before DNA polymerase bypass. However, evidence from epidemiological, animal, and other experimental studies also suggest that UVA radiation may participate in melanoma formation. The DNA damage relevant for UVA includes specific types of CPDs at TT sequences and perhaps oxidative DNA damage to guanine, both induced by direct or indirect, photosensitization-mediated chemical and biophysical processes. We summarize the evidence for a potential role of UVA in melanoma and discuss some of the mechanistic pathways of how UVA may induce mutagenesis in melanocytes.

Methylation Markers in Cutaneous Melanoma: Unravelling the Potential Utility of Their Tracking by Liquid Biopsy

Malignant melanoma is the most serious, life-threatening form of all dermatologic diseases, with a poor prognosis in the presence of metastases and advanced disease. Despite recent advances in targeted therapy and immunotherapy, there is still a critical need for a better understanding of the fundamental mechanisms behind melanoma progression and resistance onset. Recent advances in genome-wide methylation methods have revealed that aberrant changes in the pattern of DNA methylation play an important role in many aspects of cancer progression, including cell proliferation and migration, evasion of cell death, invasion, and metastasization. The purpose of the current review was to gather evidence regarding the usefulness of DNA methylation tracking in liquid biopsy as a potential biomarker in melanoma. We investigated the key genes and signal transduction pathways that have been found to be altered epigenetically in melanoma. We then highlighted the circulating tumor components present in blood, including circulating melanoma cells (CMC), circulating tumor DNA (ctDNA), and tumor-derived extracellular vesicles (EVs), as a valuable source for identifying relevant aberrations in DNA methylation. Finally, we focused on DNA methylation signatures as a marker for tracking response to therapy and resistance, thus facilitating personalized medicine and decision-making in the treatment of melanoma patients.

Machine Learning in Epigenomics: Insights into Cancer Biology and Medicine

The recent deluge of genome-wide technologies for the mapping of the epigenome and resulting data in cancer samples has provided the opportunity for gaining insights into and understanding the roles of epigenetic processes in cancer. However, the complexity, high-dimensionality, sparsity, and noise associated with these data pose challenges for extensive integrative analyses. Machine Learning (ML) algorithms are particularly suited for epigenomic data analyses due to their flexibility and ability to learn underlying hidden structures. We will discuss four overlapping but distinct major categories under ML: dimensionality reduction, unsupervised methods, supervised methods, and deep learning (DL). We review the preferred use cases of these algorithms in analyses of cancer epigenomics data with the hope to provide an overview of how ML approaches can be used to explore fundamental questions on the roles of epigenome in cancer biology and medicine.

Quantification of Global DNA Methylation in Canine Melanotic and Amelanotic Oral Mucosal Melanomas and Peripheral Blood Leukocytes From the Same Patients With OMM: First Study

Oral mucosal melanomas (OMMs) are aggressive and resistant cancers of high importance in veterinary oncology. Amelanotic OMM produces comparatively less melanin and is considered to be more aggressive than melanotic OMM. Global DNA methylation profiles with hypomethylated or hypermethylated patterns have both been associated with aggressive neoplasms; however, global DNA hypomethylation seems to correlate to higher aggressiveness. Accordingly, global DNA methylation in peripheral blood leukocytes has been investigated to understand the role of systemic or environmental factors in cancer development. This study aimed to quantify global DNA methylation in canine melanotic and amelanotic OMM samples and in the peripheral blood leukocytes of the same dogs. Tumor tissue samples were collected from 38 dogs, of which 19 were melanotic and 19 were amelanotic OMM. These were submitted to immunohistochemistry (IHC) with anti-5-methylcytosine (5mC) and anti-Ki67 primary antibodies. Ki67- and 5mC-positive nuclei were manually scored with the help of an image analysis system. Peripheral blood samples were collected from 18 among the 38 OMM-bearing dogs and from 7 additional healthy control dogs. Peripheral blood leukocytes were isolated from the 25 dogs, and DNA was extracted and analyzed by high-performance liquid chromatography (HPLC) for global DNA methylation. The pattern of global DNA methylation in both canine melanotic and amelanotic OMM indicated higher percentages of weakly or negatively stained nuclei in most of the OMM cells, presuming predominant global DNA hypomethylation. In addition, Ki67 counts in amelanotic OMM were significantly higher than those in melanotic OMM (p < 0.001). Global DNA methylation different immunostaining patterns (strong, weak or negative) correlated with Ki67 scores. Global DNA methylation in circulating leukocytes did not differ between the 9 melanotic and 9 amelanotic OMM or between the 18 OMM-bearing dogs and the 7 healthy dogs. This study provides new information on canine melanotic and amelanotic OMM based on global DNA methylation and cell proliferation.

AKT-mediated phosphorylation of Sox9 induces Sox10 transcription in a murine model of HER2-positive breast cancer

Background

Approximately 5–10% of HER2-positive breast cancers can be defined by low expression of the Ste20-like kinase, SLK, and high expression of SOX10. Our lab has observed that genetic deletion of SLK results in the induction of Sox10 and significantly accelerates tumor initiation in a HER2-induced mammary tumor model. However, the mechanism responsible for the induction of SOX10 gene expression in this context remains unknown.

Methods

Using tumor-derived cell lines from MMTV-Neu mice lacking SLK and biochemical approaches, we have characterized the signaling mechanisms and relevant DNA elements driving Sox10 expression.

Results

Biochemical and genetic analyses of the SOX10 regulatory region in SLK-deficient mammary tumor cells show that Sox10 expression is dependent on a novel −7kb enhancer that harbors three SoxE binding sites. ChIP analyses demonstrate that Sox9 is bound to those elements in vivo. Our data show that AKT can directly phosphorylate Sox9 in vitro at serine 181 and that AKT inhibition blocks Sox9 phosphorylation and Sox10 expression in SLK(-/-) tumor cells. AKT-mediated Sox9 phosphorylation increases its transcriptional activity on the Sox10 −7kb enhancer without altering its DNA-binding activity. Interestingly, analysis of murine and human mammary tumors reveals a direct correlation between the levels of active phospho-Sox9 S181 and Sox10 expression.

Conclusions

Our results have identified a novel Sox10 enhancer and validated Sox9 as a direct target for AKT. As Sox10 is a biomarker for triple-negative breast cancers (TNBC), these findings might have major implications in the targeting and treatment of those cancers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01435-6.

Quantum dynamics of the ππ*/nπ* decay of the epigenetic nucleobase 1,5-dimethyl-cytosine in the gas phase

We study the ultrafast dynamics of 1,5-dimethyl-cytosine, a model for 5-methyl-cytidine, after photoexcitation to the first two bright ππ* states, focusing on the possible population transfer to dark nπ* states. To that end we propagate the initial wave packets on the coupled potential energy surfaces of the seven lowest energy excited states modelled with a diabatic linear vibronic coupling (LVC) model, considering all the vibrational coordinates. Time-evolution is computed by the multilayer version of the multiconfigurational time dependent Hartree (ML-MCTDH) method. The LVC Hamiltonian is parametrized with time-dependent density functional theory (TD-DFT) calculations adopting PBE0 and CAM-B3LYP functionals, which provide a different energy gap between the lowest energy nπ* states and the spectroscopic ππ* state. Population of the lowest ππ* flows to a dark nπ* state which involves a lone pair (LP) of the carbonyl oxygen (nOπ*), but the extent of such transfer is much larger according to PBE0 than to CAM-B3LYP. Photoexcitation to the second bright state gives rise to much richer dynamics with an ultrafast (50 fs) complete decay to the lowest ππ*, to nOπ* and to another nπ* in which the excited electron comes from the LP of the ring nitrogen. We perform a detailed analysis of the vibronic dynamics both in terms of normal modes and valence coordinates (bond lengths and angles). The comparison with the analogous dynamics in 1-methyl-cytosine, a model for cytidine, provides insights into the effect of methylation at carbon 5 on the electronic and nuclear dynamics.

Stabilization of the Triplet Biradical Intermediate of 5-Methylcytosine Enhances Cyclobutane Pyrimidine Dimer (CPD) Formation in DNA

Cyclobutane pyrimidine dimer (CPD) is a photoproduct formed by two stacked pyrimidine bases through a cycloaddition reaction upon irradiation. Owing to its close association with skin cancer, the mechanism of CPD formation has been studied thoroughly. Among many aspects of CPD, its formation involving 5-methylcytosine (5mC) has been of special interest because the CPD yield is known to increase with C5-methylation of cytosine. In this work, high-level quantum mechanics/molecular mechanics (QM/MM) calculations are used to examine a previously experimentally detected pathway for CPD formation in hetero (thymine-cytosine and thymine-5mC) dipyrimidines, which is facilitated through intersystem crossing in thymine and formation of a triplet biradical intermediate. A DNA duplex model system containing a core sequence TmCG or TCG is used. The stabilization of a radical center in the biradical intermediate by the methyl group of 5mC can lead to increased CPD yield in TmCG compared with its non-methylated counterpart, TCG, thereby suggesting the existence of a new pathway of CPD formation enhanced by 5mC.

Epigenetic mechanisms and the hallmarks of cancer: an intimate affair

Epigenetic mechanisms comprising DNA methylation, histone modifications, and noncoding RNAs affect chromatin structure and regulate gene expression. These mechanisms control normal embryonic development and adult life and their deregulation contributes to several diseases including cancer. The process of tumorigenesis is complex and results from the evolution of different "hallmarks of cancer". Hanahan and Weinberg presented in 2000 and 2011 seminal contributions in the cancer field, first the six hallmarks of cancer and a decade later two additional hallmarks and two enabling characteristics were added. Here, we surmise that epigenetic mechanisms regulate and contribute to every single hallmark in cancer, and thus represent the hallmark of hallmarks in tumorigenesis. Focusing on epigenetics as a major hallmark in cancer formation has profound preventive, therapeutic, and clinical implications.

Targeted Therapy in Melanoma and Mechanisms of Resistance

The common mutation BRAFV600 in primary melanomas activates the mitogen-activated protein kinase/extracellular-signal-regulated kinase (MAPK/ERK) pathway and the introduction of proto-oncogene B-Raf (BRAF) and mitogen-activated protein kinase kinase (MEK) inhibitors (BRAFi and MEKi) was a breakthrough in the treatment of these cancers. However, 15–20% of tumors harbor primary resistance to this therapy, and moreover, patients develop acquired resistance to treatment. Understanding the molecular phenomena behind resistance to BRAFi/MEKis is indispensable in order to develop novel targeted therapies. Most often, resistance develops due to either the reactivation of the MAPK/ERK pathway or the activation of alternative kinase signaling pathways including phosphatase and tensin homolog (PTEN), neurofibromin 1 (NF-1) or RAS signaling. The hyperactivation of tyrosine kinase receptors, such as the receptor of the platelet-derived growth factor β (PDFRβ), insulin-like growth factor 1 receptor (IGF-1R) and the receptor for hepatocyte growth factor (HGF), lead to the induction of the AKT/3-phosphoinositol kinase (PI3K) pathway. Another pathway resulting in BRAFi/MEKi resistance is the hyperactivation of epidermal growth factor receptor (EGFR) signaling or the deregulation of microphthalmia-associated transcription factor (MITF).

Comprehensive analysis of tumor necrosis factor receptor TNFRSF9 (4-1BB) DNA methylation with regard to molecular and clinicopathological features, immune infiltrates, and response prediction to immunotherapy in melanoma

BACKGROUND

Immunotherapy, including checkpoint inhibition, has remarkably improved prognosis in advanced melanoma. Despite this success, acquired resistance is still a major challenge. The T cell costimulatory receptor TNFRSF9 (also known as 4-1BB and CD137) is a promising new target for immunotherapy and two agonistic antibodies are currently tested in clinical trials. However, little is known about epigenetic regulation of the encoding gene. In this study we investigate a possible correlation of TNFRSF9 DNA methylation with gene expression, clinicopathological parameters, molecular and immune correlates, and response to anti-PD-1 immunotherapy to assess the validity of TNFRSF9 methylation to serve as a biomarker.

METHODS

We performed a correlation analyses of methylation at twelve CpG sites within TNFRSF9 with regard to transcriptional activity, immune cell infiltration, mutation status, and survival in a cohort of N = 470 melanoma patients obtained from The Cancer Genome Atlas. Furthermore, we used quantitative methylation-specific PCR to confirm correlations in a cohort of N = 115 melanoma patients' samples (UHB validation cohort). Finally, we tested the ability of TNFRSF9 methylation and expression to predict progression-free survival (PFS) and response to anti-PD-1 immunotherapy in a cohort comprised of N = 121 patients (mRNA transcription), (mRNA ICB cohort) and a case-control study including N = 48 patients (DNA methylation, UHB ICB cohort).

FINDINGS

We found a significant inverse correlation between TNFRSF9 DNA methylation and mRNA expression levels at six of twelve analyzed CpG sites (P ≤ 0.005), predominately located in the promoter flank region. Consistent with its role as costimulatory receptor in immune cells, TNFRSF9 mRNA expression and hypomethylation positively correlated with immune cell infiltrates and an interferon-γ signature. Furthermore, elevated TNFRSF9 mRNA expression and TNFRSF9 hypomethylation correlated with superior overall survival. In patients receiving anti-PD-1 immunotherapy (mRNA ICB cohort), we found that TNFRSF9 hypermethylation and reduced mRNA expression correlated with poor PFS and response.

INTERPRETATION

Our study suggests that TNFRSF9 mRNA expression is regulated via DNA methylation. The observed correlations between TNFRSF9 DNA methylation or mRNA expression with known features of response to immune checkpoint blockage suggest TNFRSF9 methylation could serve as a biomarker in the context of immunotherapies. Concordantly, we identified a correlation between TNFRSF9 DNA methylation and mRNA expression with disease progression in patients under immunotherapy. Our study provides rationale for further investigating TNFRSF9 DNA methylation as a predictive biomarker for response to immunotherapy.

FUNDING

AF was partly funded by the Mildred Scheel Foundation. SF received funding from the University Hospital Bonn BONFOR program (O-105.0069). DN was funded in part by DFG Cluster of Excellence ImmunoSensation (EXC 1023). The funders had no role in study design, data collection and analysis, interpretation, decision to publish, or preparation of the manuscript; or any aspect pertinent to the study.

DNA hypermethylation is associated with invasive phenotype of malignant melanoma

Tumor cell invasion is one of the key processes during cancer progression, leading to life-threatening metastatic lesions in melanoma. As methylation of cancer-related genes plays a fundamental role during tumorigenesis and may lead to cellular plasticity which promotes invasion, our aim was to identify novel epigenetic markers on selected invasive melanoma cells. Using Illumina BeadChip assays and Affymetrix Human Gene 1.0 microarrays, we explored the DNA methylation landscape of selected invasive melanoma cells and examined the impact of DNA methylation on gene expression patterns. Our data revealed predominantly hypermethylated genes in the invasive cells affecting the neural crest differentiation pathway and regulation of the actin cytoskeleton. Integrative analysis of the methylation and gene expression profiles resulted in a cohort of hypermethylated genes (IL12RB2, LYPD6B, CHL1, SLC9A3, BAALC, FAM213A, SORCS1, GPR158, FBN1 and ADORA2B) with decreased expression. On the other hand, hypermethylation in the gene body of the EGFR and RBP4 genes was positively correlated with overexpression of the genes. We identified several methylation changes that can have role during melanoma progression, including hypermethylation of the promoter regions of the ARHGAP22 and NAV2 genes that are commonly altered in locally invasive primary melanomas as well as during metastasis. Interestingly, the down-regulation of the methylcytosine dioxygenase TET2 gene, which regulates DNA methylation, was associated with hypermethylated promoter region of the gene. This can probably lead to the observed global hypermethylation pattern of invasive cells and might be one of the key changes during the development of malignant melanoma cells.

Extracellular Vesicles and Epigenetic Modifications Are Hallmarks of Melanoma Progression

Cutaneous melanoma shows a high metastatic potential based on its ability to overcome the immune system's control. The mechanisms activated for these functions vary extremely and are also represented by the production of a number of extracellular vesicles including exosomes. Other vesicles showing a potential role in the melanoma progression include oncosomes and melanosomes and the majority of them mediate tumor processes including angiogenesis, immune regulation, and modifications of the micro-environment. Moreover, a number of epigenetic modifications have been described in melanoma and abundant production of altered microRNAs (mi-RNAs), non-coding RNAs, histones, and abnormal DNA methylation have been associated with different phases of melanoma progression. In addition, exosomes, miRNAs, and other molecular factors have been used as potential biomarkers reflecting disease evolution while others have been suggested to be potential druggable molecules for therapeutic application.

Epigenetic Mechanisms of Escape from BRAF Oncogene Dependency

About eight percent of all human tumors (including 50% of melanomas) carry gain-of-function mutations in the BRAF oncogene. Mutated BRAF and subsequent hyperactivation of the MAPK signaling pathway has motivated the use of MAPK-targeted therapies for these tumors. Despite great promise, however, MAPK-targeted therapies in BRAF-mutant tumors are limited by the emergence of drug resistance. Mechanisms of resistance include genetic, non-genetic and epigenetic alterations. Epigenetic plasticity, often modulated by histone-modifying enzymes and gene regulation, can influence a tumor cell's BRAF dependency and therefore, response to therapy. In this review, focusing primarily on class 1 BRAF-mutant cells, we will highlight recent work on the contribution of epigenetic mechanisms to inter- and intratumor cell heterogeneity in MAPK-targeted therapy response.

Gene Expression and Methylation Analysis in Melanomas and Melanocytes From the Same Patient: Loss of NPM2 Expression Is a Potential Immunohistochemical Marker for Melanoma

DNA methylation is considered the primary epigenetic mechanism underlying the development of malignant melanoma. Since DNA methylation can be influenced by environmental factors, it is preferable to compare cancer and normal cells from the same patient. In order to compare the methylation status in melanoma tissues and melanocytes from the same individuals, we employed a novel epidermal sheet cultivation technique to isolate normal melanocytes from unaffected sites of melanoma patients. We also analyzed primary and metastatic melanoma samples, three commercially available melanocytes, and four melanoma cell lines. Cluster analysis of DNA methylation data classified freshly isolated melanomas and melanocytes into the same group, whereas the four melanoma cell lines were clustered together in a distant clade. Moreover, our analysis discovered methylation at several novel loci (KRTCAP3, AGAP2, ZNF490), in addition to those identified in previous studies (COL1A2, GPX3); however, the latter two were not observed in fresh melanoma samples. Subsequent studies revealed that NPM2 was hypermethylated and downregulated in melanomas, which was consistent with previous reports. In many normal melanocytes, NPM2 showed distinct immunohistochemical staining, while its expression was lost in malignant melanoma cells. In particular, intraepithelial lesions of malignant melanoma, an important challenge in clinical practice, could be distinguished from benign nevi. The present findings indicate the importance of using fresh melanoma samples, not melanoma cell lines and melanocytes in epigenetic studies. In addition, NPM2 immunoreactivity could be used to differentiate melanomas from normal melanocytes or benign disease.

Demethylation by low-dose 5-aza-2'-deoxycytidine impairs 3D melanoma invasion partially through miR-199a-3p expression revealing the role of this miR in melanoma

BACKGROUND

Efficient treatments against metastatic melanoma dissemination are still lacking. Here, we report that low-cytotoxic concentrations of 5-aza-2'-deoxycytidine, a DNA demethylating agent, prevent in vitro 3D invasiveness of metastatic melanoma cells and reduce lung metastasis formation in vivo.

RESULTS

We unravelled that this beneficial effect is in part due to MIR-199A2 re-expression by promoter demethylation. Alone, this miR showed an anti-invasive and anti-metastatic effect. Throughout integration of micro-RNA target prediction databases with transcriptomic analysis after 5-aza-2'-deoxycytidine treatments, we found that miR-199a-3p downregulates set of genes significantly involved in invasion/migration processes. In addition, analysis of data from melanoma patients showed a stage- and tissue type-dependent modulation of MIR-199A2 expression by DNA methylation.

CONCLUSIONS

Thus, our data suggest that epigenetic- and/or miR-based therapeutic strategies can be relevant to limit metastatic dissemination of melanoma.

TET2-Dependent Hydroxymethylome Plasticity Reduces Melanoma Initiation and Progression

Although numerous epigenetic aberrancies accumulate in melanoma, their contribution to initiation and progression remain unclear. The epigenetic mark 5-hydroxymethylcytosine (5hmC), generated through TET-mediated DNA modification, is now referred to as the sixth base of DNA and has recently been reported as a potential biomarker for multiple types of cancer. Loss of 5hmC is an epigenetic hallmark of melanoma, but whether a decrease in 5hmc levels contributes directly to pathogenesis or whether it merely results from disease progression-associated epigenetic remodeling remains to be established. Here, we show that NRAS-driven melanomagenesis in mice is accompanied by an overall decrease in 5hmC and specific 5hmC gains in selected gene bodies. Strikingly, genetic ablation of Tet2 in mice cooperated with oncogenic NRAS^Q61K to promote melanoma initiation while suppressing specific gains in 5hmC. We conclude that TET2 acts as a barrier to melanoma initiation and progression, partly by promoting 5hmC gains in specific gene bodies. SIGNIFICANCE: This work emphasizes the importance of epigenome plasticity in cancer development and highlights the involvement of druggable epigenetic factors in cancer.

Identification of a Robust Methylation Classifier for Cutaneous Melanoma Diagnosis

Early diagnosis improves melanoma survival, yet the histopathological diagnosis of cutaneous primary melanoma can be challenging, even for expert dermatopathologists. Analysis of epigenetic alterations, such as DNA methylation, that occur in melanoma can aid in its early diagnosis. Using a genome-wide methylation screening, we assessed CpG methylation in a diverse set of 89 primary invasive melanomas, 73 nevi, and 41 melanocytic proliferations of uncertain malignant potential, classified based on interobserver review by dermatopathologists. Melanomas and nevi were split into training and validation sets. Predictive modeling in the training set using ElasticNet identified a 40-CpG classifier distinguishing 60 melanomas from 48 nevi. High diagnostic accuracy (area under the receiver operator characteristic curve = 0.996, sensitivity = 96.6%, and specificity = 100.0%) was independently confirmed in the validation set (29 melanomas, 25 nevi) and other published sample sets. The 40-CpG melanoma classifier included homeobox transcription factors and genes with roles in stem cell pluripotency or the nervous system. Application of the 40-CpG melanoma classifier to the diagnostically uncertain samples assigned melanoma or nevus status, potentially offering a diagnostic tool to assist dermatopathologists. In summary, the robust, accurate 40-CpG melanoma classifier offers a promising assay for improving primary melanoma diagnosis.

Defining Driver DNA Methylation Changes in Human Cancer

Human malignant tumors are characterized by pervasive changes in the patterns of DNA methylation. These changes include a globally hypomethylated tumor cell genome and the focal hypermethylation of numerous 5′-cytosine-phosphate-guanine-3′ (CpG) islands, many of them associated with gene promoters. It has been challenging to link specific DNA methylation changes with tumorigenesis in a cause-and-effect relationship. Some evidence suggests that cancer-associated DNA hypomethylation may increase genomic instability. Promoter hypermethylation events can lead to silencing of genes functioning in pathways reflecting hallmarks of cancer, including DNA repair, cell cycle regulation, promotion of apoptosis or control of key tumor-relevant signaling networks. A convincing argument for a tumor-driving role of DNA methylation can be made when the same genes are also frequently mutated in cancer. Many of the most commonly hypermethylated genes encode developmental transcription factors, the methylation of which may lead to permanent gene silencing. Inactivation of such genes will deprive the cells in which the tumor may initiate from the option of undergoing or maintaining lineage differentiation and will lock them into a perpetuated stem cell-like state thus providing an additional window for cell transformation.

Autophagy-dependent crosstalk between GILT and PAX-3 influences radiation sensitivity of human melanoma cells

Melanoma represents an ever-increasing problem in the western world as incidence rates continue to climb. Though manageable during early stages, late stage metastatic disease is highly resistant to current intervention. We have previously shown that gamma-interferon-inducible lysosomal thiol-reductase (GILT) enhances HLA class II antigen processing and immune detection of human melanoma cells. Here we report that GILT expression inhibits a potential target, paired box-3 (PAX-3) protein, in late stage human metastatic melanoma. We also show that GILT transfection or induction by IFN-γ, decreases PAX-3 protein expression while upregulating the expression of Daxx, which is also a repressor of PAX-3. Confocal microscopic analysis demonstrated that GILT co-localizes with PAX-3 protein, but not with Daxx within melanoma cells. Immunoprecipitation and immunoblotting studies suggest that GILT expression negatively regulates PAX-3 through the autophagy pathway, potentially resulting in increased susceptibility to conventional treatment in the form of chemotherapy or radiotherapy. While high-dose radiation is a common treatment for melanoma patients, our data suggest that GILT expression significantly increased the susceptibility of melanoma cells to low-dose radiation therapy via upregulation of tumor suppressor protein p53. Overall, these data suggest that GILT has multiple roles in inducing human melanoma cells as better targets for radiation and immunotherapy.

UV-induced damage to DNA: effect of cytosine methylation on pyrimidine dimerization

Methylation/demethylation of cytosine plays an important role in epigenetic signaling, the reversibility of epigenetic modifications offering important opportunities for targeted therapies. Actually, methylated sites have been correlated with mutational hotspots detected in skin cancers. The present brief review discusses the physicochemical parameters underlying the specific ultraviolet-induced reactivity of methylated cytosine. It focuses on dimerization reactions giving rise to cyclobutane pyrimidine dimers and pyrimidine (6–4) pyrimidone adducts. According to recent studies, four conformational and electronic factors that are affected by cytosine methylation may control these reactions: the red-shift of the absorption spectrum, the lengthening of the excited state lifetime, changes in the sugar puckering modifying the stacking between reactive pyrimidines and an increase in the rigidity of duplexes favoring excitation energy transfer toward methylated pyrimidines.

Higher Nevus Count Exhibits a Distinct DNA Methylation Signature in Healthy Human Skin: Implications for Melanoma

High nevus count is the strongest risk factor for melanoma, and although gene variants have been discovered for both traits, epigenetic variation is unexplored. We investigated 322 healthy human skin DNA methylomes associated with total body nevi count, incorporating genetic and transcriptomic variation. DNA methylation changes were identified at genes involved in melanocyte biology, such as RAF1 (P = 1.2 × 10^-6) and CTC1 (region: P = 6.3 × 10^-4), and other genes including ARRDC1 (P = 3.1 × 10^-7). A subset exhibited coordinated methylation and transcription changes within the same biopsy. The total analysis was also enriched for melanoma-associated DNA methylation variation (P = 6.33 × 10^-6). In addition, we show that skin DNA methylation is associated in cis with known genome-wide association study single nucleotide polymorphisms for nevus count, at PLA2G6 (P = 1.7 × 10^-49) and NID1 (P = 6.4 × 10^-14), as well as melanoma risk, including in or near MC1R, MX2, and TERT/CLPTM1L (P < 1 × 10^-10). Our analysis using a uniquely large dataset comprising healthy skin DNA methylomes identified known and additional regulatory loci and pathways in nevi and melanoma biology. This integrative study improves our understanding of predisposition to nevi and their potential contribution to melanoma pathogenesis.

Towards Precision Medicine in the Clinic: From Biomarker Discovery to Novel Therapeutics

Precision medicine continues to be the benchmark to which we strive in cancer research. Seeking out actionable aberrations that can be selectively targeted by drug compounds promises to optimize treatment efficacy and minimize toxicity. Utilizing these different targeted agents in combination or in sequence may further delay resistance to treatments and prolong antitumor responses. Remarkable progress in the field of immunotherapy adds another layer of complexity to the management of cancer patients. Corresponding advances in companion biomarker development, novel methods of serial tumor assessments, and innovative trial designs act synergistically to further precision medicine. Ongoing hurdles such as clonal evolution, intra- and intertumor heterogeneity, and varied mechanisms of drug resistance continue to be challenges to overcome. Large-scale data-sharing and collaborative networks using next-generation sequencing (NGS) platforms promise to take us further into the cancer 'ome' than ever before, with the goal of achieving successful precision medicine.

Precision Medicine for Molecularly Targeted Agents and Immunotherapies in Early-Phase Clinical Trials

Precision medicine in oncology promises the matching of genomic, molecular, and clinical data with underlying mechanisms of a range of novel anticancer therapeutics to develop more rational and effective antitumor strategies in a timely manner. However, despite the remarkable progress made in the understanding of novel drivers of different oncogenic processes, success rates for the approval of oncology drugs remain low with substantial fiscal consequences. In this article, we focus on how recent rapid innovations in technology have brought greater clarity to the biological and clinical complexities of different cancers and advanced the development of molecularly targeted agents and immunotherapies in clinical trials. We discuss the key challenges of identifying and validating predictive biomarkers of response and resistance using both tumor and surrogate tissues, as well as the hurdles associated with intratumor heterogeneity. Finally, we outline evolving strategies employed in early-phase trial designs that incorporate omics-based technologies.