Decoding the regulatory landscape of melanoma reveals TEADS as regulators of the invasive cell state

Monitoring melanoma patients on treatment reveals a distinct macrophage population driving targeted therapy resistance

Resistance to targeted therapy remains a major clinical challenge in melanoma. To uncover resistance mechanisms, we perform single-cell RNA sequencing on fine-needle aspirates from resistant and responding tumors of patients undergoing BRAFi/MEKi treatment. Among the genes most prominently expressed in resistant tumors is POSTN, predicted to signal to a macrophage population associated with targeted therapy resistance (TTR). Accordingly, tumors from patients with fast disease progression after therapy exhibit high POSTN expression levels and high numbers of TTR macrophages. POSTN polarizes human macrophages toward a TTR phenotype and promotes resistance to targeted therapy in a melanoma mouse model, which is associated with a phenotype change in intratumoral macrophages. Finally, polarized TTR macrophages directly protect human melanoma cells from MEKi-induced killing via CD44 receptor expression on melanoma cells. Thus, interfering with the protective activity of TTR macrophages may offer a strategy to overcome resistance to targeted therapy in melanoma.

PRRX1 silencing is required for metastatic outgrowth in melanoma and is an independent prognostic of reduced survival in patients

Paired related homeobox 1 (PRRX1) is an inducer of epithelial-to-mesenchymal transition (EMT) in different types of cancer cells. We detected low PRRX1 expression in nevus but increased levels in primary human melanoma and cell lines carrying the BRAFV600E mutation. High expression of PRRX1 correlates with invasiveness and enrichment of genes belonging to the EMT programme. Conversely, we found that loss of PRRX1 in metastatic samples is an independent prognostic predictor of poor survival for melanoma patients. Here, we show that stable depletion of PRRX1 improves the growth of melanoma xenografts and increases the number of distant spontaneous metastases, compared to controls. We provide evidence that loss of PRRX1 counteracts the EMT phenotype, impairing the expression of other EMT-related transcription factors, causing dysregulation of the ERK and signal transducer and activator of transcription 3 (STAT3) signaling pathways, and abrogating the invasive and migratory properties of melanoma cells while triggering the up-regulation of proliferative/melanocytic genes and the expression of the neural-crest-like markers nerve growth factor receptor (NGFR; also known as neurotrophin receptor p75NTR) and neural cell adhesion molecule L1 (L1CAM). Overall, our results indicate that loss of PRRX1 triggers a switch in the invasive programme, and cells de-differentiate towards a neural crest stem cell (NCSC)-like phenotype that accounts for the metastatic aggressiveness.

Sirtuin 5-mediated deacetylation of TAZ at K54 promotes melanoma development

PURPOSE

Nuclear accumulation of YAP/TAZ promotes tumorigenesis in several cancers, including melanoma. Although the mechanisms underlying the nuclear retention of YAP are known, those underlying the retention of TAZ remain unclear. Our study investigates a novel acetylation/deacetylation switch in TAZ, governing its subcellular localization in melanoma tumorigenesis.

METHODS

Immunoprecipitation/Western blot assessed TAZ protein interactions and acetylation. SIRT5 activity was quantified with enzyme-linked immunosorbent assay. Immunofluorescence indicated TAZ nuclear localization. TEAD transcriptional activity was measured through luciferase reporter assays. ChIP detected TAZ binding to the CTGF promoter. Transwell and wound healing assays quantified melanoma cell invasiveness and migration. Metastasis was evaluated using a mouse model via tail vein injections. Clinical relevance was explored via immunohistochemical staining of patient tumors.

RESULTS

CBP facilitated TAZ acetylation at K54 in response to epidermal growth factor stimulation, while SIRT5 mediated deacetylation. Acetylation correlated with phosphorylation, regulating TAZ's binding with LATS2 or TEAD. TAZ K54 acetylation enhanced its S89 phosphorylation, promoting cytosolic retention via LATS2 interaction. SIRT5-mediated deacetylation enhanced TAZ-TEAD interaction and nuclear retention. Chromatin IP showed SIRT5-deacetylated TAZ recruited to CTGF promoter, boosting transcriptional activity. In a mouse model, SIRT5 overexpression induced melanoma metastasis to lung tissue following the injection of B16F10 melanocytes via the tail vein, and this effect was prevented by verteporfin treatment.

CONCLUSIONS

Our study revealed a novel mechanism of TAZ nuclear retention regulated by SIRT5-mediated K54 deacetylation and demonstrated the significance of TAZ deacetylation in CTGF expression. This study highlights the potential implications of the SIRT5/TAZ axis for treating metastatic melanoma.

Genetic dependencies associated with transcription factor activities in human cancer cell lines

Transcription factors (TFs) are important mediators of aberrant transcriptional programs in cancer cells. In this study, we focus on TF activity (TFa) as a biomarker for cell-line-selective anti-proliferative effects, in that high TFa predicts sensitivity to loss of function of a given gene (i.e., genetic dependencies [GDs]). Our linear-regression-based framework identifies 3,047 pan-cancer and 3,952 cancer-type-specific candidate TFa-GD associations from cell line data, which are then cross-examined for impact on survival in patient cohorts. One of the most prominent biomarkers is TEAD1 activity, whose associations with its predicted GDs are validated through experimental evidence as proof of concept. Overall, these TFa-GD associations represent an attractive resource for identifying innovative, biomarker-driven hypotheses for drug discovery programs in oncology.

Nuclear Localization of Yes-Associated Protein Is Associated With Tumor Progression in Cutaneous Melanoma

Yes-associated protein (YAP), an effector molecule of the Hippo signaling pathway, is expressed at high levels in cutaneous melanoma. However, the role of YAP in melanoma progression according to cellular localization is poorly understood. Tissues from 140 patients with invasive melanoma were evaluated by immunohistochemistry. Flow cytometry, western blotting, viability assays, wound healing assays, verteporfin treatment, and xenograft assays were conducted using melanoma cell lines B16F1 and B16F10 subjected to YapS127A transfection and siYap knockdown. Nuclear YAP localization was identified in 63 tumors (45.0%) and was more frequent than cytoplasmic YAP in acral lentiginous and nodular subtypes (P = .007). Compared with cytoplasmic YAP melanomas, melanomas with nuclear YAP had higher mitotic activity (P = .016), deeper invasion (P < .001), and more frequently metastasized to lymph nodes (P < .001) and distant organs (P < .001). Patients with nuclear YAP melanomas had poorer disease-free survival (P < .001) and overall survival (P < .001). Nuclear YAP was an independent risk factor for distant metastasis (hazard ratio: 3.206; 95% CI, 1.032-9.961; P = .044). Proliferative ability was decreased in siYapB16F1 (P < .001) and siYapB16F10 (P = .001) cells and increased in YapS127AB16F1 (P = .003) and YapS127AB16F10 (P = .002) cells. Cell cycle analysis demonstrated relative G1 retention in siYapB16F1 (P < .001) and siYapB16F10 (P < .001) cells and S retention in YapS127AB16F1 cells (P = .008). Wound healing assays showed that Yap knockdown inhibited cell invasion (siYapB16F1, P = .001; siYapB16F10, P < .001), whereas nuclear YAP promoted it (YapS127AB16F, P < .001; YapS127AB16F1, P = .017). Verteporfin, a direct YAP inhibitor, reduced cellular proliferation in B16F1 (P = .003) and B16F10 (P < .001) cells. Proliferative effects of nuclear YAP were confirmed in xenograft mice (P < .001). In conclusion, nuclear YAP in human melanomas showed subtype specificity and correlated with proliferative activity and proinvasiveness. It is expected that YAP becomes a useful prognostic marker, and its inhibition may be a potential therapy for melanoma patients.

ZEB1 controls a lineage-specific transcriptional program essential for melanoma cell state transitions

Cell plasticity sustains intra-tumor heterogeneity and treatment resistance in melanoma. Deciphering the transcriptional mechanisms governing reversible phenotypic transitions between proliferative/differentiated and invasive/stem-like states is required. Expression of the ZEB1 transcription factor is frequently activated in melanoma, where it fosters adaptive resistance to targeted therapies. Here, we performed a genome-wide characterization of ZEB1 transcriptional targets, by combining ChIP-sequencing and RNA-sequencing, upon phenotype switching in melanoma models. We identified and validated ZEB1 binding peaks in the promoter of key lineage-specific genes crucial for melanoma cell identity. Mechanistically, ZEB1 negatively regulates SOX10-MITF dependent proliferative/melanocytic programs and positively regulates AP-1 driven invasive and stem-like programs. Comparative analyses with breast carcinoma cells revealed lineage-specific ZEB1 binding, leading to the design of a more reliable melanoma-specific ZEB1 regulon. We then developed single-cell spatial multiplexed analyses to characterize melanoma cell states intra-tumoral heterogeneity in human melanoma samples. Combined with scRNA-Seq analyses, our findings confirmed increased ZEB1 expression in Neural-Crest-like cells and mesenchymal cells, underscoring its significance in vivo in both populations. Overall, our results define ZEB1 as a major transcriptional regulator of cell states transitions and provide a better understanding of lineage-specific transcriptional programs sustaining intra-tumor heterogeneity in melanoma.

Links between melanoma germline risk loci, driver genes and comorbidities: insight from a tissue-specific multi-omic analysis

Genome-wide association studies (GWAS) have associated 76 loci with the risk of developing melanoma. However, understanding the molecular basis of such associations has remained a challenge because most of these loci are in non-coding regions of the genome. Here, we integrated data on epigenomic markers, three-dimensional (3D) genome organization, and expression quantitative trait loci (eQTL) from melanoma-relevant tissues and cell types to gain novel insights into the mechanisms underlying melanoma risk. This integrative approach revealed a total of 151 target genes, both near and far away from the risk loci in linear sequence, with known and novel roles in the etiology of melanoma. Using protein-protein interaction networks, we identified proteins that interact-directly or indirectly-with the products of the target genes. The interacting proteins were enriched for known melanoma driver genes. Further integration of these target genes into tissue-specific gene regulatory networks revealed patterns of gene regulation that connect melanoma to its comorbidities. Our study provides novel insights into the biological implications of genetic variants associated with melanoma risk.

Identification of a Gene Signature That Predicts Dependence upon YAP/TAZ-TEAD

Targeted therapies are effective cancer treatments when accompanied by accurate diagnostic tests that can help identify patients that will respond to those therapies. The YAP/TAZ-TEAD axis is activated and plays a causal role in several cancer types, and TEAD inhibitors are currently in early-phase clinical trials in cancer patients. However, a lack of a reliable way to identify tumors with YAP/TAZ-TEAD activation for most cancer types makes it difficult to determine which tumors will be susceptible to TEAD inhibitors. Here, we used a combination of RNA-seq and bioinformatic analysis of metastatic melanoma cells to develop a YAP/TAZ gene signature. We found that the genes in this signature are TEAD-dependent in several melanoma cell lines, and that their expression strongly correlates with YAP/TAZ activation in human melanomas. Using DepMap dependency data, we found that this YAP/TAZ signature was predictive of melanoma cell dependence upon YAP/TAZ or TEADs. Importantly, this was not limited to melanoma because this signature was also predictive when tested on a panel of over 1000 cancer cell lines representing numerous distinct cancer types. Our results suggest that YAP/TAZ gene signatures like ours may be effective tools to predict tumor cell dependence upon YAP/TAZ-TEAD, and thus potentially provide a means to identify patients likely to benefit from TEAD inhibitors.

Melanoma and sex hormones: Pathogenesis, progressive disease and response to treatments

Cutaneous melanoma represents the fifth tumor in terms of incidence in young adults, with a major involvement of males than females. Despite the significant changes in available effective treatments for cutaneous melanoma, there is still a proportion of patients that do not benefit long-term disease control with immune checkpoint inhibitors and/or BRAF/MEK inhibitors, and eventually develop progressive disease. In addition to the emerging biomarkers under investigation to understand resistance to treatments, recent studies resumed the role of sex hormones (estrogens, progesterone and androgens) in melanoma patients. In the last decades, the impact of sex hormones has been considered controversial in melanoma patients, but actual growing preclinical and clinical evidence underline the potential influence on melanoma cells' growth, tumor microenvironment, the immune system and consequently on the course of disease.This review will provide available insights on the role of sex hormones in melanoma pathogenesis, disease progression and response/resistance to systemic treatments. We will also offer an overview on the recent studies on the theme, describing the hormonal contribution to disease response and the interaction with targeted therapies and immune-checkpoint inhibitors in cutaneous melanoma patients, illustrating an insight into future research in this field.

DNA damage remodels the MITF interactome to increase melanoma genomic instability

Since genome instability can drive cancer initiation and progression, cells have evolved highly effective and ubiquitous DNA damage response (DDR) programs. However, some cells (for example, in skin) are normally exposed to high levels of DNA-damaging agents. Whether such high-risk cells possess lineage-specific mechanisms that tailor DNA repair to the tissue remains largely unknown. Using melanoma as a model, we show here that the microphthalmia-associated transcription factor MITF, a lineage addition oncogene that coordinates many aspects of melanocyte and melanoma biology, plays a nontranscriptional role in shaping the DDR. On exposure to DNA-damaging agents, MITF is phosphorylated at S325, and its interactome is dramatically remodeled; most transcription cofactors dissociate, and instead MITF interacts with the MRE11-RAD50-NBS1 (MRN) complex. Consequently, cells with high MITF levels accumulate stalled replication forks and display defects in homologous recombination-mediated repair associated with impaired MRN recruitment to DNA damage. In agreement with this, high MITF levels are associated with increased single-nucleotide and copy number variant burdens in melanoma. Significantly, the SUMOylation-defective MITF-E318K melanoma predisposition mutation recapitulates the effects of DNA-PKcs-phosphorylated MITF. Our data suggest that a nontranscriptional function of a lineage-restricted transcription factor contributes to a tissue-specialized modulation of the DDR that can impact cancer initiation.

Differential histone acetylation and super-enhancer regulation underlie melanoma cell dedifferentiation

Dedifferentiation or phenotype switching refers to the transition from a proliferative to an invasive cellular state. We previously identified a 122-gene epigenetic gene signature that classifies primary melanomas as low versus high risk (denoted as Epgn1 or Epgn3). We found that the transcriptomes of the Epgn1 low-risk and Epgn3 high-risk cells are similar to the proliferative and invasive cellular states, respectively. These signatures were further validated in melanoma tumor samples. Examination of the chromatin landscape revealed differential H3K27 acetylation in the Epgn1 low-risk versus Epgn3 high-risk cell lines that corroborated with a differential super-enhancer and enhancer landscape. Melanocytic lineage genes (MITF, its targets and regulators) were associated with super-enhancers in the Epgn1 low-risk state, whereas invasiveness genes were linked with Epgn3 high-risk status. We identified the ITGA3 gene as marked by a super-enhancer element in the Epgn3 invasive cells. Silencing of ITGA3 enhanced invasiveness in both in vitro and in vivo systems, suggesting it as a negative regulator of invasion. In conclusion, we define chromatin landscape changes associated with Epgn1/Epgn3 and phenotype switching during early steps of melanoma progression that regulate transcriptional reprogramming. This super-enhancer and enhancer-driven epigenetic regulatory mechanism resulting in major changes in the transcriptome could be important in future therapeutic targeting efforts.

The CoREST repressor complex mediates phenotype switching and therapy resistance in melanoma

Virtually all patients with BRAF-mutant melanoma develop resistance to MAPK inhibitors largely through nonmutational events. Although the epigenetic landscape is shown to be altered in therapy-resistant melanomas and other cancers, a specific targetable epigenetic mechanism has not been validated. Here, we evaluated the corepressor for element 1-silencing transcription factor (CoREST) epigenetic repressor complex and the recently developed bivalent inhibitor corin within the context of melanoma phenotype plasticity and therapeutic resistance. We found that CoREST was a critical mediator of the major distinct melanoma phenotypes and that corin treatment of melanoma cells led to phenotype reprogramming. Global assessment of transcript and chromatin changes conferred by corin revealed specific effects on histone marks connected to epithelial-mesenchymal transition-associated (EMT-associated) transcription factors and the dual-specificity phosphatases (DUSPs). Remarkably, treatment of BRAF inhibitor-resistant (BRAFi-R) melanomas with corin promoted resensitization to BRAFi therapy. DUSP1 was consistently downregulated in BRAFi-R melanomas, which was reversed by corin treatment and associated with inhibition of p38 MAPK activity and resensitization to BRAFi therapies. Moreover, this activity was recapitulated by the p38 MAPK inhibitor BIRB 796. These findings identify the CoREST repressor complex as a central mediator of melanoma phenotype plasticity and resistance to targeted therapy and suggest that CoREST inhibitors may prove beneficial for patients with BRAFi-resistant melanoma.

Cell-type-directed design of synthetic enhancers

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

Systematic analysis of the transcriptional landscape of melanoma reveals drug-target expression plasticity

Metastatic melanoma originates from melanocytes of the skin. Melanoma metastasis results in poor treatment prognosis for patients and is associated with epigenetic and transcriptional changes that reflect the developmental program of melanocyte differentiation from neural crest stem cells. Several studies have explored melanoma transcriptional heterogeneity using microarray, bulk and single-cell RNA-sequencing technologies to derive data-driven models of the transcriptional-state change which occurs during melanoma progression. No study has systematically examined how different models of melanoma progression derived from different data types, technologies and biological conditions compare. Here, we perform a cross-sectional study to identify averaging effects of bulk-based studies that mask and distort apparent melanoma transcriptional heterogeneity; we describe new transcriptionally distinct melanoma cell states, identify differential co-expression of genes between studies and examine the effects of predicted drug susceptibilities of different cell states between studies. Importantly, we observe considerable variability in drug-target gene expression between studies, indicating potential transcriptional plasticity of melanoma to down-regulate these drug targets and thereby circumvent treatment. Overall, observed differences in gene co-expression and predicted drug susceptibility between studies suggest bulk-based transcriptional measurements do not reliably gauge heterogeneity and that melanoma transcriptional plasticity is greater than described when studies are considered in isolation.

ATP1A1 is a promising new target for melanoma treatment and can be inhibited by its physiological ligand bufalin to restore targeted therapy efficacy

Despite advancements in treating metastatic melanoma, many patients exhibit resistance to targeted therapies. Our study focuses on ATP1A1, a sodium pump subunit associated with cancer development. We aimed to assess ATP1A1 prognostic value in melanoma patients and examine the impact of its ligand, bufalin, on melanoma cell lines in vitro and in vivo. High ATP1A1 expression (IHC) correlated with reduced overall survival in melanoma patients. Resistance to BRAF inhibitor was linked to elevated ATP1A1 levels in patient biopsies (IHC, qPCR) and cell lines (Western blot, qPCR). Additionally, high ATP1A1 mRNA expression positively correlated with differentiation/pigmentation markers based on data from The Cancer Genome Atlas (TCGA) databases and Verfaillie proliferative gene signature analysis. Bufalin specifically targeted ATP1A1 in caveolae, (proximity ligation assay) and influenced Src phosphorylation (Western blot), thereby disrupting multiple signaling pathways (phosphokinase array). In vitro, bufalin induced apoptosis in melanoma cell lines by acting on ATP1A1 (siRNA experiments) and, in vivo, significantly impeded melanoma growth using a nude mouse xenograft model with continuous bufalin delivery via an osmotic pump. In conclusion, our study demonstrates that ATP1A1 could serve as a prognostic marker for patient survival and a predictive marker for response to BRAF inhibitor therapy. By targeting ATP1A1, bufalin inhibited cell proliferation, induced apoptosis in vitro, and effectively suppressed tumor development in mice. Thus, our findings strongly support ATP1A1 as a promising therapeutic target, with bufalin as a potential agent to disrupt its tumor-promoting activity.

A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Expression pattern of osteopontin isoforms in malignant melanoma cell lines

Osteopontin (OPN) is a secreted integrin-binding protein that plays a role in inflammation, cellular viability, cell adhesion and migration, cancer development, and diabetes through different mechanisms. The splice variants of OPN can play essential roles in cancer development, progression, and metastasis formation; however, limited data are available about the role of OPN isoforms in human malignant melanoma. Our goal was to define the gene expression patterns of five OPN variants (OPN4, OPN5, OPNa, OPNb, and OPNc), integrin, and CD44 receptor genes in primary and metastatic melanoma-originated cell lines (n = 19), and to explore the association of the expression patterns with clinicopathological parameters. We evaluated the invasive property of the cell lines and investigated the potential association between the invasion and gene expression of OPN isoforms. We found a significant rise in the expression of OPNc in the invasive cell lines compared to the noninvasive cells and detected significantly higher expression of the OPN splice variants in melanoma cell lines originating from more advanced stages tumors than cell lines originating from early-stage melanomas. The correlation analysis revealed that all five OPN variants positively correlated with ITGB3 and ITGA9, whereas OPN5 positively correlated with ITGB1, ITGAV, ITGA6, and CD44. OPN can activate extracellular signal-regulated kinase signaling through binding to α9β1 integrin, promoting melanoma tumor cell migration. It is possible that such associations between OPN splice variants and integrin receptors may play a role in melanoma progression. In conclusion, our findings suggest that high expression of OPNc correlates with the invasive behavior of melanoma cells.

Sox10-Deficient Drug-Resistant Melanoma Cells Are Refractory to Oncolytic RNA Viruses

Targeted therapy resistance frequently develops in melanoma due to intratumor heterogeneity and epigenetic reprogramming. This also typically induces cross-resistance to immunotherapies. Whether this includes additional modes of therapy has not been fully assessed. We show that co-treatments of MAPKi with VSV-based oncolytics do not function in a synergistic fashion; rather, the MAPKis block infection. Melanoma resistance to vemurafenib further perturbs the cells' ability to be infected by oncolytic viruses. Resistance to vemurafenib can be induced by the loss of SOX10, a common proliferative marker in melanoma. The loss of SOX10 promotes a cross-resistant state by further inhibiting viral infection and replication. Analysis of RNA-seq datasets revealed an upregulation of interferon-stimulated genes (ISGs) in SOX10 knockout populations and targeted therapy-resistant cells. Interestingly, the induction of ISGs appears to be independent of type I IFN production. Overall, our data suggest that the pathway mediating oncolytic resistance is due to the loss of SOX10 during acquired drug resistance in melanoma.

Clinical Significance of SOX10 Expression in Human Pathology

The embryonic development of neural crest cells and subsequent tissue differentiation are intricately regulated by specific transcription factors. Among these, SOX10, a member of the SOX gene family, stands out. Located on chromosome 22q13, the SOX10 gene encodes a transcription factor crucial for the differentiation, migration, and maintenance of tissues derived from neural crest cells. It plays a pivotal role in developing various tissues, including the central and peripheral nervous systems, melanocytes, chondrocytes, and odontoblasts. Mutations in SOX10 have been associated with congenital disorders such as Waardenburg-Shah Syndrome, PCWH syndrome, and Kallman syndrome, underscoring its clinical significance. Furthermore, SOX10 is implicated in neural and neuroectodermal tumors, such as melanoma, malignant peripheral nerve sheath tumors (MPNSTs), and schwannomas, influencing processes like proliferation, migration, and differentiation. In mesenchymal tumors, SOX10 expression serves as a valuable marker for distinguishing between different tumor types. Additionally, SOX10 has been identified in various epithelial neoplasms, including breast, ovarian, salivary gland, nasopharyngeal, and bladder cancers, presenting itself as a potential diagnostic and prognostic marker. However, despite these associations, further research is imperative to elucidate its precise role in these malignancies.

Exploring the role of epigenetic alterations and non-coding RNAs in melanoma pathogenesis and therapeutic strategies

Melanoma is a rare but highly lethal type of skin cancer whose incidence is increasing globally. Melanoma is characterized by high resistance to therapy and relapse. Despite significant advances in the treatment of metastatic melanoma, many patients experience progression due to resistance mechanisms. Epigenetic changes, including alterations in chromatin remodeling, DNA methylation, histone modifications, and non-coding RNA rearrangements, contribute to neoplastic transformation, metastasis, and drug resistance in melanoma. This review summarizes current research on epigenetic mechanisms in melanoma and their therapeutic potential. Specifically, we discuss the role of histone acetylation and methylation in gene expression regulation and melanoma pathobiology, as well as the promising results of HDAC inhibitors and DNMT inhibitors in clinical trials. We also examine the dysregulation of non-coding RNA, particularly miRNAs, and their potential as targets for melanoma therapy. Finally, we highlight the challenges of epigenetic therapies, such as the complexity of epigenetic mechanisms combined with immunotherapies and the need for combination therapies to overcome drug resistance. In conclusion, epigenetic changes may be reversible, and the use of combination therapy between traditional therapies and epigenetically targeted drugs could be a viable solution to reverse the increasing number of patients who develop treatment resistance or even prevent it. While several clinical trials are underway, the complexity of these mechanisms presents a significant challenge to the development of effective therapies. Further research is needed to fully understand the role of epigenetic mechanisms in melanoma and to develop more effective and targeted therapies.

Combinatorial Treatment with PARP and MAPK Inhibitors Overcomes Phenotype Switch-Driven Drug Resistance in Advanced Melanoma

Metastatic melanoma is either intrinsically resistant or rapidly acquires resistance to targeted therapy treatments, such as MAPK inhibitors (MAPKi). A leading cause of resistance to targeted therapy is a dynamic transition of melanoma cells from a proliferative to a highly invasive state, a phenomenon called phenotype switching. Mechanisms regulating phenotype switching represent potential targets for improving treatment of patients with melanoma. Using a drug screen targeting chromatin regulators in patient-derived three-dimensional MAPKi-resistant melanoma cell cultures, we discovered that PARP inhibitors (PARPi) restore sensitivity to MAPKis, independent of DNA damage repair pathways. Integrated transcriptomic, proteomic, and epigenomic analyses demonstrated that PARPis induce lysosomal autophagic cell death, accompanied by enhanced mitochondrial lipid metabolism that ultimately increases antigen presentation and sensitivity to T-cell cytotoxicity. Moreover, transcriptomic and epigenetic rearrangements induced by PARP inhibition reversed epithelial-mesenchymal transition-like phenotype switching, which redirected melanoma cells toward a proliferative and MAPKi-sensitive state. The combination of PARP and MAPKis synergistically induced cancer cell death both in vitro and in vivo in patient-derived xenograft models. Therefore, this study provides a scientific rationale for treating patients with melanoma with PARPis in combination with MAPKis to abrogate acquired therapy resistance.

SIGNIFICANCE

PARP inhibitors can overcome resistance to MAPK inhibitors by activating autophagic cell death and reversing phenotype switching, suggesting that this synergistic combination could help improve the prognosis of patients with melanoma.

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.

HDAC8-mediated inhibition of EP300 drives a transcriptional state that increases melanoma brain metastasis

Melanomas can adopt multiple transcriptional states. Little is known about the epigenetic drivers of these cell states, limiting our ability to regulate melanoma heterogeneity. Here, we identify stress-induced HDAC8 activity as driving melanoma brain metastasis development. Exposure of melanocytes and melanoma cells to multiple stresses increases HDAC8 activation leading to a neural crest-stem cell transcriptional state and an amoeboid, invasive phenotype that increases seeding to the brain. Using ATAC-Seq and ChIP-Seq we show that increased HDAC8 activity alters chromatin structure by increasing H3K27ac and enhancing accessibility at c-Jun binding sites. Functionally, HDAC8 deacetylates the histone acetyltransferase EP300, causing its enzymatic inactivation. This, in turn, increases binding of EP300 to Jun-transcriptional sites and decreases binding to MITF-transcriptional sites. Inhibition of EP300 increases melanoma cell invasion, resistance to stress and increases melanoma brain metastasis development. HDAC8 is identified as a mediator of transcriptional co-factor inactivation and chromatin accessibility that drives brain metastasis.

Super-enhancer-driven expression of BAHCC1 promotes melanoma cell proliferation and genome stability

Super-enhancers (SEs) are stretches of enhancers ensuring a high level of expression of key genes associated with cell function. The identification of cancer-specific SE-driven genes is a powerful means for the development of innovative therapeutic strategies. Here, we identify a MITF/SOX10/TFIIH-dependent SE promoting the expression of BAHCC1 in a broad panel of melanoma cells. BAHCC1 is highly expressed in metastatic melanoma and is required for tumor engraftment, growth, and dissemination. Integrative genomics analyses reveal that BAHCC1 is a transcriptional regulator controlling expression of E2F/KLF-dependent cell-cycle and DNA-repair genes. BAHCC1 associates with BRG1-containing remodeling complexes at the promoters of these genes. BAHCC1 silencing leads to decreased cell proliferation and delayed DNA repair. Consequently, BAHCC1 deficiency cooperates with PARP inhibition to induce melanoma cell death. Our study identifies BAHCC1 as an SE-driven gene expressed in melanoma and demonstrates how its inhibition can be exploited as a therapeutic target.

Disrupting cellular memory to overcome drug resistance

Gene expression states persist for varying lengths of time at the single-cell level, a phenomenon known as gene expression memory. When cells switch states, losing memory of their prior state, this transition can occur in the absence of genetic changes. However, we lack robust methods to find regulators of memory or track state switching. Here, we develop a lineage tracing-based technique to quantify memory and identify cells that switch states. Applied to melanoma cells without therapy, we quantify long-lived fluctuations in gene expression that are predictive of later resistance to targeted therapy. We also identify the PI3K and TGF-β pathways as state switching modulators. We propose a pretreatment model, first applying a PI3K inhibitor to modulate gene expression states, then applying targeted therapy, which leads to less resistance than targeted therapy alone. Together, we present a method for finding modulators of gene expression memory and their associated cell fates.

Addressing erroneous scale assumptions in microbe and gene set enrichment analysis

By applying Differential Set Analysis (DSA) to sequence count data, researchers can determine whether groups of microbes or genes are differentially enriched. Yet sequence count data suffer from a scale limitation: these data lack information about the scale (i.e., size) of the biological system under study, leading some authors to call these data compositional (i.e., proportional). In this article, we show that commonly used DSA methods that rely on normalization make strong, implicit assumptions about the unmeasured system scale. We show that even small errors in these scale assumptions can lead to positive predictive values as low as 9%. To address this problem, we take three novel approaches. First, we introduce a sensitivity analysis framework to identify when modeling results are robust to such errors and when they are suspect. Unlike standard benchmarking studies, this framework does not require ground-truth knowledge and can therefore be applied to both simulated and real data. Second, we introduce a statistical test that provably controls Type-I error at a nominal rate despite errors in scale assumptions. Finally, we discuss how the impact of scale limitations depends on a researcher's scientific goals and provide tools that researchers can use to evaluate whether their goals are at risk from erroneous scale assumptions. Overall, the goal of this article is to catalyze future research into the impact of scale limitations in analyses of sequence count data; to illustrate that scale limitations can lead to inferential errors in practice; yet to also show that rigorous and reproducible scale reliant inference is possible if done carefully.

Podoplanin Drives Amoeboid Invasion in Canine and Human Mucosal Melanoma

Mucosal melanoma metastasizes at an early stage of the disease in human and dog. We revealed that overexpression of podoplanin in tumor invasion fronts (IF) was related to poor prognosis of dogs with mucosal melanoma. Moreover, podoplanin expressed in canine mucosal melanoma cells promotes proliferation and aggressive amoeboid invasion by activating Rho-associated kinase (ROCK)-myosin light chain 2 (MLC2) signaling. PDPN-ROCK-MLC2 signaling plays a role in cell-cycle arrest and cellular senescence escape as a mechanism for regulating proliferation. Podoplanin induces amoeboid invasion in the IFs of mouse xenografted tumor tissues, similar to canine mucosal melanoma clinical samples. We further identified that podoplanin expression was related to poor prognosis of human patients with mucosal melanoma, and human mucosal melanoma with podoplanin-high expression enriched gene signatures related to amoeboid invasion, similar to canine mucosal melanoma. Overall, we propose that podoplanin promotes canine and human mucosal melanoma metastasis by inducing aggressive amoeboid invasion and naturally occurring canine mucosal melanoma can be a novel research model for podoplanin expressing human mucosal melanoma.

IMPLICATIONS

Podoplanin could be a new therapeutic target to restrict the metastatic dissemination of canine and human mucosal melanoma.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Landscape of enhancer disruption and functional screen in melanoma cells

BACKGROUND

The high mutation rate throughout the entire melanoma genome presents a major challenge in stratifying true driver events from the background mutations. Numerous recurrent non-coding alterations, such as those in enhancers, can shape tumor evolution, thereby emphasizing the importance in systematically deciphering enhancer disruptions in melanoma.

RESULTS

Here, we leveraged 297 melanoma whole-genome sequencing samples to prioritize highly recurrent regions. By performing a genome-scale CRISPR interference (CRISPRi) screen on highly recurrent region-associated enhancers in melanoma cells, we identified 66 significant hits which could have tumor-suppressive roles. These functional enhancers show unique mutational patterns independent of classical significantly mutated genes in melanoma. Target gene analysis for the essential enhancers reveal many known and hidden mechanisms underlying melanoma growth. Utilizing extensive functional validation experiments, we demonstrate that a super enhancer element could modulate melanoma cell proliferation by targeting MEF2A, and another distal enhancer is able to sustain PTEN tumor-suppressive potential via long-range interactions.

CONCLUSIONS

Our study establishes a catalogue of crucial enhancers and their target genes in melanoma growth and progression, and illuminates the identification of novel mechanisms of dysregulation for melanoma driver genes and new therapeutic targeting strategies.

Understanding cancer drug resistance with Sleeping Beauty functional genomic screens: Application to MAPK inhibition in cutaneous melanoma

Combined BRAF and MEK inhibition is an effective treatment for BRAF-mutant cutaneous melanoma. However, most patients progress on this treatment due to drug resistance. Here, we applied the Sleeping Beauty transposon system to understand how melanoma evades MAPK inhibition. We found that the specific drug resistance mechanisms differed across melanomas in our genetic screens of five cutaneous melanoma cell lines. While drivers that reactivated MAPK were highly conserved, many others were cell-line specific. One such driver, VAV1, activated a de-differentiated transcriptional program like that of hyperactive RAC1, RAC1P29S. To target this mechanism, we showed that an inhibitor of SRC, saracatinib, blunts the VAV1-induced transcriptional reprogramming. Overall, we highlighted the importance of accounting for melanoma heterogeneity in treating cutaneous melanoma with MAPK inhibitors. Moreover, we demonstrated the utility of the Sleeping Beauty transposon system in understanding cancer drug resistance.

AKT2 Loss Impairs BRAF-Mutant Melanoma Metastasis

Despite recent advances in treatment, melanoma remains the deadliest form of skin cancer due to its highly metastatic nature. Melanomas harboring oncogenic BRAFV600E mutations combined with PTEN loss exhibit unrestrained PI3K/AKT signaling and increased invasiveness. However, the contribution of different AKT isoforms to melanoma initiation, progression, and metastasis has not been comprehensively explored, and questions remain about whether individual isoforms play distinct or redundant roles in each step. We investigate the contribution of individual AKT isoforms to melanoma initiation using a novel mouse model of AKT isoform-specific loss in a murine melanoma model, and we investigate tumor progression, maintenance, and metastasis among a panel of human metastatic melanoma cell lines using AKT isoform-specific knockdown studies. We elucidate that AKT2 is dispensable for primary tumor formation but promotes migration and invasion in vitro and metastatic seeding in vivo, whereas AKT1 is uniquely important for melanoma initiation and cell proliferation. We propose a mechanism whereby the inhibition of AKT2 impairs glycolysis and reduces an EMT-related gene expression signature in PTEN-null BRAF-mutant human melanoma cells to limit metastatic spread. Our data suggest that the elucidation of AKT2-specific functions in metastasis might inform therapeutic strategies to improve treatment options for melanoma patients.

Exploring the Complex and Multifaceted Interplay between Melanoma Cells and the Tumor Microenvironment

Malignant melanoma is a very aggressive skin cancer, characterized by a heterogeneous nature and high metastatic potential. The incidence of melanoma is continuously increasing worldwide, and it is one of the most common cancers in young adults. In the past twenty years, our understanding of melanoma biology has increased profoundly, and disease management for patients with disseminated disease has improved due to the emergence of immunotherapy and targeted therapy. However, a significant fraction of patients relapse or do not respond adequately to treatment. This can partly be explained by the complex signaling between the tumor and its microenvironment, giving rise to melanoma phenotypes with different patterns of disease progression. This review focuses on the key aspects and complex relationship between pathogenesis, genetic abnormalities, tumor microenvironment, cellular plasticity, and metabolic reprogramming in melanoma. By acquiring a deeper understanding of the multifaceted features of melanomagenesis, we can reach a point of more individualized and patient-centered disease management and reduced costs of ineffective treatments.

Gene signature reveals decreased SOX10-dependent transcripts in malignant cells from immune checkpoint inhibitor-resistant cutaneous melanomas

Evidence is mounting for cross-resistance between immune checkpoint and targeted kinase inhibitor therapies in cutaneous melanoma patients. Since the loss of the transcription factor, SOX10, causes tolerance to MAPK pathway inhibitors, we used bioinformatic techniques to determine if reduced SOX10 expression/activity is associated with immune checkpoint inhibitor resistance. We integrated SOX10 ChIP-seq, knockout RNA-seq, and knockdown ATAC-seq data from melanoma cell models to develop a robust SOX10 gene signature. We used computational methods to validate this signature as a measure of SOX10-dependent activity in independent single-cell and bulk RNA-seq SOX10 knockdown, cell line panel, and MAPK inhibitor drug-resistant datasets. Evaluation of patient single-cell RNA-seq data revealed lower levels of SOX10-dependent transcripts in immune checkpoint inhibitor-resistant tumors. Our results suggest that SOX10-deficient melanoma cells are associated with cross-resistance between targeted and immune checkpoint inhibitors and highlight the need to identify therapeutic strategies that target this subpopulation.

SCENIC+: single-cell multiomic inference of enhancers and gene regulatory networks

Joint profiling of chromatin accessibility and gene expression in individual cells provides an opportunity to decipher enhancer-driven gene regulatory networks (GRNs). Here we present a method for the inference of enhancer-driven GRNs, called SCENIC+. SCENIC+ predicts genomic enhancers along with candidate upstream transcription factors (TFs) and links these enhancers to candidate target genes. To improve both recall and precision of TF identification, we curated and clustered a motif collection with more than 30,000 motifs. We benchmarked SCENIC+ on diverse datasets from different species, including human peripheral blood mononuclear cells, ENCODE cell lines, melanoma cell states and Drosophila retinal development. Next, we exploit SCENIC+ predictions to study conserved TFs, enhancers and GRNs between human and mouse cell types in the cerebral cortex. Finally, we use SCENIC+ to study the dynamics of gene regulation along differentiation trajectories and the effect of TF perturbations on cell state. SCENIC+ is available at scenicplus.readthedocs.io .

Dynamical modeling of proliferative-invasive plasticity and IFNγ signaling in melanoma reveals mechanisms of PD-L1 expression heterogeneity

BACKGROUND

Phenotypic heterogeneity of melanoma cells contributes to drug tolerance, increased metastasis, and immune evasion in patients with progressive disease. Diverse mechanisms have been individually reported to shape extensive intra-tumor and inter-tumor phenotypic heterogeneity, such as IFNγ signaling and proliferative to invasive transition, but how their crosstalk impacts tumor progression remains largely elusive.

METHODS

Here, we integrate dynamical systems modeling with transcriptomic data analysis at bulk and single-cell levels to investigate underlying mechanisms behind phenotypic heterogeneity in melanoma and its impact on adaptation to targeted therapy and immune checkpoint inhibitors. We construct a minimal core regulatory network involving transcription factors implicated in this process and identify the multiple 'attractors' in the phenotypic landscape enabled by this network. Our model predictions about synergistic control of PD-L1 by IFNγ signaling and proliferative to invasive transition were validated experimentally in three melanoma cell lines-MALME3, SK-MEL-5 and A375.

RESULTS

We demonstrate that the emergent dynamics of our regulatory network comprising MITF, SOX10, SOX9, JUN and ZEB1 can recapitulate experimental observations about the co-existence of diverse phenotypes (proliferative, neural crest-like, invasive) and reversible cell-state transitions among them, including in response to targeted therapy and immune checkpoint inhibitors. These phenotypes have varied levels of PD-L1, driving heterogeneity in immunosuppression. This heterogeneity in PD-L1 can be aggravated by combinatorial dynamics of these regulators with IFNγ signaling. Our model predictions about changes in proliferative to invasive transition and PD-L1 levels as melanoma cells evade targeted therapy and immune checkpoint inhibitors were validated in multiple RNA-seq data sets from in vitro and in vivo experiments.

CONCLUSION

Our calibrated dynamical model offers a platform to test combinatorial therapies and provide rational avenues for the treatment of metastatic melanoma. This improved understanding of crosstalk among PD-L1 expression, proliferative to invasive transition and IFNγ signaling can be leveraged to improve the clinical management of therapy-resistant and metastatic melanoma.

AKT2 Loss Impairs BRAF-Mutant Melanoma Metastasis

Despite recent advances in treatment, melanoma remains the deadliest form of skin cancer, due to its highly metastatic nature. Melanomas harboring oncogenic BRAF V600E mutations combined with PTEN loss exhibit unrestrained PI3K/AKT signaling and increased invasiveness. However, the contribution of different AKT isoforms to melanoma initiation, progression, and metastasis has not been comprehensively explored, and questions remain whether individual isoforms play distinct or redundant roles in each step. We investigate the contribution of individual AKT isoforms to melanoma initiation using a novel mouse model of AKT isoform-specific loss in a murine melanoma model, and investigate tumor progression, maintenance, and metastasis among a panel of human metastatic melanoma cell lines using AKT-isoform specific knockdown studies. We elucidate that AKT2 is dispensable for primary tumor formation but promotes migration and invasion in vitro and metastatic seeding in vivo , while AKT1 is uniquely important for melanoma initiation and cell proliferation. We propose a mechanism whereby inhibition of AKT2 impairs glycolysis and reduces an EMT-related gene expression signature in PTEN-null BRAF-mutant human melanoma cells to limit metastatic spread. Our data suggest that elucidation of AKT2-specific functions in metastasis could inform therapeutic strategies to improve treatment options for melanoma patients.

Multiple Fra-1-bound enhancers showing different molecular and functional features can cooperate to repress gene transcription

BACKGROUND

How transcription factors (TFs) down-regulate gene expression remains ill-understood, especially when they bind to multiple enhancers contacting the same gene promoter. In particular, it is not known whether they exert similar or significantly different molecular effects at these enhancers.

RESULTS

To address this issue, we used a particularly well-suited study model consisting of the down-regulation of the TGFB2 gene by the TF Fra-1 in Fra-1-overexpressing cancer cells, as Fra-1 binds to multiple enhancers interacting with the TGFB2 promoter. We show that Fra-1 does not repress TGFB2 transcription via reducing RNA Pol II recruitment at the gene promoter but by decreasing the formation of its transcription-initiating form. This is associated with complex long-range chromatin interactions implicating multiple molecularly and functionally heterogeneous Fra-1-bound transcriptional enhancers distal to the TGFB2 transcriptional start site. In particular, the latter display differential requirements upon the presence and the activity of the lysine acetyltransferase p300/CBP. Furthermore, the final transcriptional output of the TGFB2 gene seems to depend on a balance between the positive and negative effects of Fra-1 at these enhancers.

CONCLUSION

Our work unveils complex molecular mechanisms underlying the repressive actions of Fra-1 on TGFB2 gene expression. This has consequences for our general understanding of the functioning of the ubiquitous transcriptional complex AP-1, of which Fra-1 is the most documented component for prooncogenic activities. In addition, it raises the general question of the heterogeneity of the molecular functions of TFs binding to different enhancers regulating the same gene.

Dynamical modelling of proliferative-invasive plasticity and IFNγ signaling in melanoma reveals mechanisms of PD-L1 expression heterogeneity

Phenotypic heterogeneity of melanoma cells contributes to drug tolerance, increased metastasis, and immune evasion in patients with progressive disease. Diverse mechanisms have been individually reported to shape extensive intra- and inter-tumoral phenotypic heterogeneity, such as IFNγ signaling and proliferative to invasive transition, but how their crosstalk impacts tumor progression remains largely elusive. Here, we integrate dynamical systems modeling with transcriptomic data analysis at bulk and single-cell levels to investigate underlying mechanisms behind phenotypic heterogeneity in melanoma and its impact on adaptation to targeted therapy and immune checkpoint inhibitors. We construct a minimal core regulatory network involving transcription factors implicated in this process and identify the multiple "attractors" in the phenotypic landscape enabled by this network. Our model predictions about synergistic control of PD-L1 by IFNγ signaling and proliferative to invasive transition were validated experimentally in three melanoma cell lines - MALME3, SK-MEL-5 and A375. We demonstrate that the emergent dynamics of our regulatory network comprising MITF, SOX10, SOX9, JUN and ZEB1 can recapitulate experimental observations about the co-existence of diverse phenotypes (proliferative, neural crest-like, invasive) and reversible cell-state transitions among them, including in response to targeted therapy and immune checkpoint inhibitors. These phenotypes have varied levels of PD-L1, driving heterogeneity in immune-suppression. This heterogeneity in PD-L1 can be aggravated by combinatorial dynamics of these regulators with IFNγ signaling. Our model predictions about changes in proliferative to invasive transition and PD-L1 levels as melanoma cells evade targeted therapy and immune checkpoint inhibitors were validated in multiple data sets from in vitro and in vivo experiments. Our calibrated dynamical model offers a platform to test combinatorial therapies and provide rational avenues for the treatment of metastatic melanoma. This improved understanding of crosstalk among PD-L1 expression, proliferative to invasive transition and IFNγ signaling can be leveraged to improve the clinical management of therapy-resistant and metastatic melanoma.

Lipid droplets are a metabolic vulnerability in melanoma

Melanoma exhibits numerous transcriptional cell states including neural crest-like cells as well as pigmented melanocytic cells. How these different cell states relate to distinct tumorigenic phenotypes remains unclear. Here, we use a zebrafish melanoma model to identify a transcriptional program linking the melanocytic cell state to a dependence on lipid droplets, the specialized organelle responsible for lipid storage. Single-cell RNA-sequencing of these tumors show a concordance between genes regulating pigmentation and those involved in lipid and oxidative metabolism. This state is conserved across human melanoma cell lines and patient tumors. This melanocytic state demonstrates increased fatty acid uptake, an increased number of lipid droplets, and dependence upon fatty acid oxidative metabolism. Genetic and pharmacologic suppression of lipid droplet production is sufficient to disrupt cell cycle progression and slow melanoma growth in vivo. Because the melanocytic cell state is linked to poor outcomes in patients, these data indicate a metabolic vulnerability in melanoma that depends on the lipid droplet organelle.

The journey from melanocytes to melanoma

Over the past decade, melanoma has led the field in new cancer treatments, with impressive gains in on-treatment survival but more modest improvements in overall survival. Melanoma presents heterogeneity and transcriptional plasticity that recapitulates distinct melanocyte developmental states and phenotypes, allowing it to adapt to and eventually escape even the most advanced treatments. Despite remarkable advances in our understanding of melanoma biology and genetics, the melanoma cell of origin is still fiercely debated because both melanocyte stem cells and mature melanocytes can be transformed. Animal models and high-throughput single-cell sequencing approaches have opened new opportunities to address this question. Here, we discuss the melanocytic journey from the neural crest, where they emerge as melanoblasts, to the fully mature pigmented melanocytes resident in several tissues. We describe a new understanding of melanocyte biology and the different melanocyte subpopulations and microenvironments they inhabit, and how this provides unique insights into melanoma initiation and progression. We highlight recent findings on melanoma heterogeneity and transcriptional plasticity and their implications for exciting new research areas and treatment opportunities. The lessons from melanocyte biology reveal how cells that are present to protect us from the damaging effects of ultraviolet radiation reach back to their origins to become a potentially deadly cancer.

DNA damage-induced interaction between a lineage addiction oncogenic transcription factor and the MRN complex shapes a tissue-specific DNA Damage Response and cancer predisposition

Since genome instability can drive cancer initiation and progression, cells have evolved highly effective and ubiquitous DNA Damage Response (DDR) programs. However, some cells, in skin for example, are normally exposed to high levels of DNA damaging agents. Whether such high-risk cells possess lineage-specific mechanisms that tailor DNA repair to the tissue remains largely unknown. Here we show, using melanoma as a model, that the microphthalmia-associated transcription factor MITF, a lineage addition oncogene that coordinates many aspects of melanocyte and melanoma biology, plays a non-transcriptional role in shaping the DDR. On exposure to DNA damaging agents, MITF is phosphorylated by ATM/DNA-PKcs, and unexpectedly its interactome is dramatically remodelled; most transcription (co)factors dissociate, and instead MITF interacts with the MRE11-RAD50-NBS1 (MRN) complex. Consequently, cells with high MITF levels accumulate stalled replication forks, and display defects in homologous recombination-mediated repair associated with impaired MRN recruitment to DNA damage. In agreement, high MITF levels are associated with increased SNV burden in melanoma. Significantly, the SUMOylation-defective MITF-E318K melanoma predisposition mutation recapitulates the effects of ATM/DNA-PKcs-phosphorylated MITF. Our data suggest that a non-transcriptional function of a lineage-restricted transcription factor contributes to a tissue-specialised modulation of the DDR that can impact cancer initiation.

BORIS/CTCFL-mediated chromatin accessibility alterations promote a pro-invasive transcriptional signature in melanoma cells

Melanoma is the deadliest form of skin cancer, due to its tendency to metastasize early. Brother of regulator of imprinted sites (BORIS), also known as CCCTC binding factor-like (CTCFL), is a transcription regulator that becomes ectopically expressed in melanoma. We recently showed that BORIS contributes to melanoma phenotype switching by altering the gene expression program of melanoma cells from an intermediate melanocytic state toward a more mesenchymal-like state. However, the mechanism underlying this transcriptional switch remains unclear. Here, ATAC-seq was used to study BORIS-mediated chromatin accessibility alterations in melanoma cells harboring an intermediate melanocytic state. The gene set that gained promoter accessibility, following ectopic BORIS expression, showed enrichment for biological processes associated with melanoma invasion, while promoters of genes associated with proliferation showed reduced accessibility. Integration of ATAC-seq and RNA-seq data demonstrated that increased chromatin accessibility was associated with transcriptional upregulation of genes involved in tumor progression processes, and the aberrant activation of oncogenic transcription factors, while reduced chromatin accessibility and downregulated genes were associated with repressed activity of tumor suppressors and proliferation factors. Together, these findings indicate that BORIS mediates transcriptional reprogramming in melanoma cells by altering chromatin accessibility and gene expression, shifting the cellular transcription landscape of melanoma cells toward a mesenchymal-like genetic signature.

An in vivo avian model of human melanoma to perform rapid and robust preclinical studies

Metastatic melanoma patients carrying a BRAF^V600 mutation can be treated with a combination of BRAF and MEK inhibitors (BRAFi/MEKi), but innate and acquired resistance invariably occurs. Predicting patient response to targeted therapies is crucial to guide clinical decision. We describe here the development of a highly efficient patient-derived xenograft model adapted to patient melanoma biopsies, using the avian embryo as a host (AVI-PDX^TM ). In this in vivo paradigm, we depict a fast and reproducible tumor engraftment of patient samples within the embryonic skin, preserving key molecular and phenotypic features. We show that sensitivity and resistance to BRAFi/MEKi can be reliably modeled in these AVI-PDX^TM , as well as synergies with other drugs. We further provide proof-of-concept that the AVI-PDX^TM models the diversity of responses of melanoma patients to BRAFi/MEKi, within days, hence positioning it as a valuable tool for the design of personalized medicine assays and for the evaluation of novel combination strategies.

The mechanical phenotypic plasticity of melanoma cell: an emerging driver of therapy cross-resistance

Advanced cutaneous melanoma is the deadliest form of skin cancer and one of the most aggressive human cancers. Targeted therapies (TT) against BRAF mutated melanoma and immune checkpoints blockade therapies (ICB) have been a breakthrough in the treatment of metastatic melanoma. However, therapy-driven resistance remains a major hurdle in the clinical management of the metastatic disease. Besides shaping the tumor microenvironment, current treatments impact transition states to promote melanoma cell phenotypic plasticity and intratumor heterogeneity, which compromise treatment efficacy and clinical outcomes. In this context, mesenchymal-like dedifferentiated melanoma cells exhibit a remarkable ability to autonomously assemble their own extracellular matrix (ECM) and to biomechanically adapt in response to therapeutic insults, thereby fueling tumor relapse. Here, we review recent studies that highlight mechanical phenotypic plasticity of melanoma cells as a hallmark of adaptive and non-genetic resistance to treatment and emerging driver in cross-resistance to TT and ICB. We also discuss how targeting BRAF-mutant dedifferentiated cells and ECM-based mechanotransduction pathways may overcome melanoma cross-resistance.

Cooperative induction of receptor tyrosine kinases contributes to adaptive MAPK drug resistance in melanoma through the PI3K pathway

Vemurafenib-induced drug resistance in melanoma has been linked to receptor tyrosine kinase (RTK) upregulation. The MITF and SOX10 genes play roles as master regulators of melanocyte and melanoma development. Here, we aimed to explore the complex mechanisms behind the MITF/SOX10-controlled RTK-induced drug resistance in melanoma. To achieve this, we used a number of molecular techniques, including melanoma patient data from TCGA, vemurafenib-resistant melanoma cell lines, and knock-down studies. The melanoma cell lines were classified as proliferative or invasive based upon their MITF/AXL expression activity. We measured the change of expression activity for MITF/SOX10 and their receptor (AXL/ERBB3) and ligand (NRG1/GAS6) targets known to be involved in RTK-induced drug resistance after vemurafenib treatment. We find that melanoma cell lines characterized as proliferative (high MITF low AXL), transform into an invasive (low MITF, high AXL) cell state after vemurafenib resistance, indicating novel feedback loops and advanced compensatory regulation mechanisms between the master regulators, receptors, and ligands involved in vemurafenib-induced resistance. Together, our data disclose fine-tuned mechanisms involved in RTK-facilitated vemurafenib resistance that will be challenging to overcome by using single drug targeting strategies against melanoma.

Escape from NK cell tumor surveillance by NGFR-induced lipid remodeling in melanoma

Metastatic disease is a major cause of death for patients with melanoma. Melanoma cells can become metastatic not only due to cell-intrinsic plasticity but also due to cancer-induced protumorigenic remodeling of the immune microenvironment. Here, we report that innate immune surveillance by natural killer (NK) cells is bypassed by human melanoma cells expressing the stem cell marker NGFR. Using in vitro and in vivo cytotoxic assays, we show that NGFR protects melanoma cells from NK cell-mediated killing and, furthermore, boosts metastasis formation in a mouse model with adoptively transferred human NK cells. Mechanistically, NGFR leads to down-regulation of NK cell activating ligands and simultaneous up-regulation of the fatty acid stearoyl-coenzyme A desaturase (SCD) in melanoma cells. Notably, pharmacological and small interfering RNA-mediated inhibition of SCD reverted NGFR-induced NK cell evasion in vitro and in vivo. Hence, NGFR orchestrates immune control antagonizing pathways to protect melanoma cells from NK cell clearance, which ultimately favors metastatic disease.

Identifying the Transcriptional Drivers of Metastasis Embedded within Localized Melanoma

In melanoma, predicting which tumors will ultimately metastasize guides treatment decisions. Transcriptional signatures of primary tumors have been utilized to predict metastasis, but which among these are driver or passenger events remains unclear. We used data from the adjuvant AVAST-M trial to identify a predictive gene signature in localized tumors that ultimately metastasized. Using a zebrafish model of primary melanoma, we interrogated the top genes from the AVAST-M signature in vivo. This identified GRAMD1B, a cholesterol transfer protein, as a bona fide metastasis suppressor, with a majority of knockout animals rapidly developing metastasis. Mechanistically, excess free cholesterol or its metabolite 27-hydroxycholesterol promotes invasiveness via activation of an AP-1 program, which is associated with increased metastasis in humans. Our data demonstrate that the transcriptional seeds of metastasis are embedded within localized tumors, suggesting that early targeting of these programs can be used to prevent metastatic relapse.

SIGNIFICANCE

We analyzed human melanoma transcriptomics data to identify a gene signature predictive of metastasis. To rapidly test clinical signatures, we built a genetic metastasis platform in adult zebrafish and identified GRAMD1B as a suppressor of melanoma metastasis. GRAMD1B-associated cholesterol overload activates an AP-1 program to promote melanoma invasion. This article is highlighted in the In This Issue feature, p. 1.

Immunoediting instructs tumor metabolic reprogramming to support immune evasion

Immunoediting sculpts immunogenicity and thwarts host anti-tumor responses in tumor cells during tumorigenesis; however, it remains unknown whether metabolic programming of tumor cells can be guided by immunosurveillance. Here, we report that T cell-mediated immunosurveillance in early-stage tumorigenesis instructs c-Myc upregulation and metabolic reprogramming in tumor cells. This previously unexplored tumor-immune interaction is controlled by non-canonical interferon gamma (IFNγ)-STAT3 signaling and supports tumor immune evasion. Our findings uncover that immunoediting instructs deregulated bioenergetic programs in tumor cells to empower them to disarm the T cell-mediated immunosurveillance by imposing metabolic tug-of-war between tumor and infiltrating T cells and forming the suppressive tumor microenvironment.

Genomics and Epigenomics in the Molecular Biology of Melanoma-A Prerequisite for Biomarkers Studies

Melanoma is a common and aggressive tumor originating from melanocytes. The increasing incidence of cutaneous melanoma in recent last decades highlights the need for predictive biomarkers studies. Melanoma development is a complex process, involving the interplay of genetic, epigenetic, and environmental factors. Genetic aberrations include BRAF, NRAS, NF1, MAP2K1/MAP2K2, KIT, GNAQ, GNA11, CDKN2A, TERT mutations, and translocations of kinases. Epigenetic alterations involve microRNAs, non-coding RNAs, histones modifications, and abnormal DNA methylations. Genetic aberrations and epigenetic marks are important as biomarkers for the diagnosis, prognosis, and prediction of disease recurrence, and for therapeutic targets. This review summarizes our current knowledge of the genomic and epigenetic changes in melanoma and discusses the latest scientific information.

A novel lncRNA MDHDH suppresses glioblastoma multiforme by acting as a scaffold for MDH2 and PSMA1 to regulate NAD+ metabolism and autophagy

BACKGROUND

To identify potential targets related to nicotinamide adenine dinucleotide (NAD+) metabolism in gliomas, we used RNA immunoprecipitation to identify a novel long noncoding RNA renamed malate dehydrogenase degradation helper (MDHDH) (NONCODE annotation ID: NONHSAT138800.2, NCBI Reference Sequence: NR_028345), which bound to MDH2 (malate dehydrogenase 2), that is downregulated in glioblastoma multiforme (GBM) and associated with metabolic regulation. However, its underlying mechanisms in the progression of GBM have not been well studied.

METHODS

To investigate the clinical significance of MDHDH, we analyzed its expression levels in publicly available datasets and collected clinical samples from Shandong Provincial Hospital, affiliated with Shandong University. Functional assays, including FISH/CISH, CCK8, EdU, wound healing, and transwell assays, were used to determine the cellular/subcellular localization, tissue expression profile and anti-oncogenic role of MDHDH. Furthermore, RNA pulldown, mass spectrometry RNA immunoprecipitation, coimmunoprecipitation, JC-1 probe, and cell energy-production assays were used to determine the mechanisms of MDHDH in the development of GBM. Animal experiments were conducted to determine the antitumorigenic role of MDHDH in GBM in vivo.

RESULTS

In public datasets, MDHDH expression was significantly downregulated in GBM and LGG compared with GTEx normal brain tissues. The results of the tissue microarray showed that the MDHDH expression level negatively correlated with the tumor grade. Altered MDHDH expression led to significant changes in the proliferation, migration and invasion of GBM cells both in vitro and in vivo. Mechanistically, we found that MDHDH directly bound to MDH2 and PSMA1 (20S proteasomal core subunit alpha-type 1) as a molecular scaffold and accelerated the degradation of MDH2 by promoting the binding of ubiquitinated MDH2 to the proteasome. The degradation of MDH2 subsequently led to changes in the mitochondrial membrane potential and NAD+/NADH ratio, which impeded glycolysis in glioma cells.

CONCLUSIONS

In conclusion, this study broadened our understanding of the functions of lncRNAs in GBM. We demonstrated that the tumor suppressor MDHDH might act as a clinical biomarker and that the overexpression of MDHDH might be a novel synergistic strategy for enhancing metabolism-based, epigenetic-based, and autophagy regulation-based therapies with clinical benefits for glioblastoma multiforme patients.