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

Metabolic buffering suppresses phenotype switching in cancer

The impact of the microenvironment on epigenetically plastic cancer cells underpins phenotypic heterogeneity, a major cause of metastatic dissemination and therapy resistance that together represent the primary cause of cancer-related death. Nutrient limitation is a key microenvironmental stress that can cause a phenotypic transition from proliferation to invasion via activation of the integrated stress response. However, whether and how the capacity to store and mobilize nutrients impacts phenotype-switching through metabolic buffering remains unknown. Here, using melanoma as a model, we reveal that the ability to accumulate and mobilize glycogen, that buffers glucose availability, plays a key role in phenotypic transitions in melanoma. While proliferative phenotype cells exhibit high levels of glycogen, invasion is marked by low glycogen levels. Significantly, an inability to store and metabolize glycogen leads to phenotype instability and a switch to invasion. Accordingly, glycogen levels inversely correlate with Clark levels in primary melanomas, with low expression of the glycogen phosphorylases PYGB/L and phosphoglucomutase 1 (PGM1) being associated with worse overall survival. The importance of metabolic buffering in suppressing phenotypic transitions likely extrapolates to other cancer types.

Highlights

Melanoma phenotypes are distinguished by their ability to store and mobilize glycogen. Proliferative MITF High melanoma cells store glycogen to improve survival under stressful conditions. Inhibition of glycogen degradation impairs proliferation in MITF High melanoma cells. Lack of PGM1 drives invasion and metastatic dissemination.

Therapeutic Strategies for Drug-resistant Melanoma and Their Clinical Implications

Melanoma is a malignant tumor originating from melanocytes, characterized by its high invasiveness and metastasis, leading to poor prognosis and high mortality. Early-stage melanoma is primarily treated with surgery; however, due to its metastatic nature, surgery becomes challenging in advanced stages. Treatment strategies for advanced or metastatic melanoma include chemotherapy, radiation therapy, and targeted therapy. However, melanoma's propensity for rapid drug resistance remains a significant clinical challenge. This review summarizes the developments in the treatment of drug-resistant melanoma over the past decade and discusses the advantages and disadvantages of various therapeutic approaches and their clinical significance implications.

Loss of fatty acid-binding protein 7 promotes B16F10 melanoma metastasis

Melanoma possesses the characteristic phenotypic plasticity, enhancing its metastatic formation and drug resistance. Lipid and fatty acid metabolism are usually altered to support melanoma progression and can be targeted for therapeutic development. Fatty acid binding protein 7 (FABP7) is highly expressed in melanomas and is shown to support its proliferation, migration, and invasion, but the mechanisms remain unclear. Our study aimed to link FABP7 to lipid metabolism and phenotypic shift in melanomas. We established the Fabp7-knockout (KO) B16F10 melanoma cells, which showed an enhanced invasion through matrix-coated membrane, without significant change in proliferation. Similar outcomes were obtained when using RNA interference targeting FABP7. Fabp7-KO cells injected into mice exhibited slower primary tumor growth, but formed higher metastatic foci count in the lungs. We also discovered a higher saturation in overall lipids, phosphatidylcholines, and triacylglycerols. We observed transcriptional shifts toward the invasive MITFLow/AXLHigh phenotype, with upregulation of transforming growth factor-beta (TGF-β) receptor mRNAs. In conclusion, FABP7 may help balancing lipid saturation and maintain the proliferative state of melanomas, mitigating invasiveness and metastatic formation.

Schwann cells in regeneration and cancer: an epithelial-mesenchymal transition perspective

In the peripheral nervous system, glial cells, known as Schwann cells (SCs), are responsible for supporting and maintaining nerves. One of the most important characteristics of SCs is their remarkable plasticity. In various injury contexts, SCs undergo a reprogramming process that generates specialized cells to promote tissue regeneration and repair. However, in pathological conditions, this same plasticity and regenerative potential can be hijacked. Different studies highlight the activation of the epithelial-mesenchymal transition (EMT) as a driver of SC phenotypic plasticity. Although SCs are not epithelial, their neural crest origin makes EMT activation crucial for their ability to adopt repair phenotypes, mirroring the plasticity observed during development. These adaptive processes are essential for regeneration. However, EMT activation in SCs-derived tumours enhances cancer progression and aggressiveness. Furthermore, in the tumour microenvironment (TME), SCs also acquire activated phenotypes that contribute to tumour migration and invasion by activating EMT in cancer cells. In this review, we will discuss how EMT impacts SC plasticity and function from development and tissue regeneration to pathological conditions, such as cancer.

The long non-coding RNA ROSALIND protects the mitochondrial translational machinery from oxidative damage

Upregulation of mitochondrial respiration coupled with high ROS-scavenging capacity is a characteristic shared by drug-tolerant cells in several cancers. As translational fidelity is essential for cell fitness, protection of the mitochondrial and cytosolic ribosomes from oxidative damage is pivotal. While mechanisms for recognising and repairing such damage exist in the cytoplasm, the corresponding process in the mitochondria remains unclear.By performing Ascorbate PEroXidase (APEX)-proximity ligation assay directed to the mitochondrial matrix followed by isolation and sequencing of RNA associated to mitochondrial proteins, we identified the nuclear-encoded lncRNA ROSALIND as an interacting partner of ribosomes. ROSALIND is upregulated in recurrent tumours and its expression can discriminate between responders and non-responders to immune checkpoint blockade in a melanoma cohort of patients. Featuring an unusually high G content, ROSALIND serves as a substrate for oxidation. Consequently, inhibiting ROSALIND leads to an increase in ROS and protein oxidation, resulting in severe mitochondrial respiration defects. This, in turn, impairs melanoma cell viability and increases immunogenicity both in vitro and ex vivo in preclinical humanised cancer models. These findings underscore the role of ROSALIND as a novel ROS buffering system, safeguarding mitochondrial translation from oxidative stress, and shed light on potential therapeutic strategies for overcoming cancer therapy resistance.

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.

PIEZO1 Affects Cell Growth and Migration via Microfilament-Mediated YAP trans-Latitudinal Regulation

Environmental mechanical forces, such as cell membrane stress, cell extrusion, and stretch, have been proven to affect cell growth and migration. Piezo1, a mechanosensitive channel protein, responds directly to endogenous or exogenous mechanical stimuli. Here, we explored the Piezo1 distribution and microfilament morphological changes induced by mechanical forces in the tumor and normal cells. In addition, Piezo1 activation in tumor cells resulted in the nuclear accumulation of YAP, whereas nuclear export of YAP and microfilament depolymerization occurred with the prolonged activation, while a removal stimulation restored the YAP distribution and microfilament polymerization. Combining the morphological changes of the microfilament under Piezo1 activation and the function of YAP in regulating cell growth and development, we suggest that Piezo1 senses changes in environmental mechanical forces and regulates YAP distribution through the microfilament cytoskeleton network, which in turn affects the growth and migration more obviously in tumor cells rather than normal cells. Our results are essential for understanding the trans-latitudinal transmission of mechanical forces and exploring the role of environmental mechanical forces in tumor therapy.

Pan-inhibition of super-enhancer-driven oncogenic transcription by next-generation synthetic ecteinascidins yields potent anti-cancer activity

The plasticity of cancer cells facilitates their ability to adopt heterogeneous differentiation states, posing a significant challenge to therapeutic interventions. Specific gene expression programs, driven in part by super-enhancers (SEs), underlie cancer cell states. Here we successfully inhibit SE-driven transcription in phenotypically distinct metastatic melanoma cells using next-generation synthetic ecteinascidins. Through functional genomic methodologies, we demonstrate that these compounds inhibit the expression of genes encoding lineage-specific or ubiquitous transcription factors/coactivators by selectively targeting the CpG-rich sequences within their promoters and/or enhancers. This prevents the formation of transcription factor/coactivator condensates necessary for SE-dependent gene expression. Consequently, these compounds exhibit cytotoxic activity across distinct subpopulations of metastatic melanoma cells and inhibit tumor proliferation, including those resistant to current therapies. These findings extend to other cancers, like small cell lung cancer, recently approved for ecteinascidin-based treatment. Overall, our study provides preclinical proof that pan-inhibition of SE-dependent genes with synthetic ecteinascidins is a promising therapeutic approach for tumors with heterogeneous transcriptional landscapes.

TGFβ signaling sensitizes MEKi-resistant human melanoma to targeted therapy-induced apoptosis

The TGFβ signaling pathway is known for its pleiotropic functions in a plethora of biological processes. In melanoma, TGFβ signaling promotes invasiveness and metastasis formation. However, its involvement in the response to therapy is controversial. While several studies have linked TGFβ signaling to elevated resistance to targeted therapy in melanoma, separate findings have indicated a favorable treatment response through TGFβ-mediated increase of cell death. We now found that the outcome of TGFβ signaling in the context of targeted therapy is dose dependent. Unlike low doses, high levels of TGFβ signal activation induce apoptosis upon simultaneous MAPK pathway inhibition, even in targeted therapy resistant melanoma cell lines. Using transcriptomic analyses, combined with genomic target identification of the critical TGFβ signaling effector SMAD4, we demonstrate that parallel activation of TGFβ signaling and MAPK pathway inhibition causes a complete switch of TGFβ target genes from promoting pro-invasive processes to fueling pro-apoptotic pathways. Investigations of underlying mechanisms identified a novel apoptosis-inducing gene signature. Functional validation of signature members highlighted a central role of the pro-apoptotic BCL2 family member BCL2L11 (BIM) in mediating apoptosis in this condition. Using a modified, synthetic version of the TGFB1 mRNA for intra-tumoral injections, we additionally showcase a potential therapeutic application of this treatment combination.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

The classification of melanocytic gene signatures

Gene expression profiling technologies have revolutionized cell biology, enabling researchers to identify gene signatures linked to various biological attributes of melanomas, such as pigmentation status, differentiation state, proliferative versus invasive capacity, and disease progression. Although the discovery of gene signatures has significantly enhanced our understanding of melanocytic phenotypes, reconciling the numerous signatures reported across independent studies and different profiling platforms remains a challenge. Current methods for classifying melanocytic gene signatures depend on exact gene overlap and comparison with unstandardized baseline transcriptomes. In this study, we aimed to categorize published gene signatures into clusters based on their similar patterns of expression across clinical cutaneous melanoma specimens. We analyzed nearly 800 melanoma samples from six gene expression repositories and developed a classification framework for gene signatures that is resilient against biases in gene identification across profiling platforms and inconsistencies in baseline standards. Using 39 frequently cited published gene signatures, our analysis revealed seven principal classes of gene signatures that correlate with previously identified phenotypes: Differentiated, Mitotic/MYC, AXL, Amelanotic, Neuro, Hypometabolic, and Invasive. Each class is consistent with the phenotypes that the constituent gene signatures represent, and our classification method does not rely on overlapping genes between signatures. To facilitate broader application, we created WIMMS (what is my melanocytic signature, available at https://wimms.tanlab.org/), a user-friendly web application. WIMMS allows users to categorize any gene signature, determining its relationship to predominantly cited signatures and its representation within the seven principal classes.

Melanoma-specific mutation hotspots in distal, non-coding, promoter-interacting regions implicate novel candidate driver genes

BACKGROUND

To develop targeted treatments, it is crucial to identify the full spectrum of genetic drivers in melanoma, including those in non-coding regions. However, recent efforts to explore non-coding regions have primarily focused on gene-adjacent elements such as promoters and non-coding RNAs, leaving intergenic distal regulatory elements largely unexplored.

METHODS

We used Hi-C chromatin contact data from melanoma cells to map distal, non-coding, promoter-interacting regulatory elements genome-wide in melanoma. Using this "promoter-interaction network", alongside whole-genome sequence and gene expression data from the Pan Cancer Analysis of Whole Genomes, we developed multivariate linear regression models to identify distal somatic mutation hotspots that affect promoter activity.

RESULTS

We identified eight recurrently mutated hotspots that are novel, melanoma-specific, located in promoter-interacting distal regulatory elements, alter transcription factor binding motifs, and affect the expression of genes (e.g., HSPB7, CLDN1, ADCY9 and FDXR) previously implicated as tumour suppressors/oncogenes in various cancers.

CONCLUSIONS

Our study suggests additional non-coding drivers beyond the well-characterised TERT promoter in melanoma, offering new insights into the disruption of complex regulatory networks by non-coding mutations that may contribute to melanoma development. Furthermore, our study provides a framework for integrating multiple levels of biological data to uncover cancer-specific non-coding drivers.

Single-Nuclei Transcriptome Profiling Reveals Intra-Tumoral Heterogeneity and Characterizes Tumor Microenvironment Architecture in a Murine Melanoma Model

Malignant melanoma is an aggressive cancer, with a high risk of metastasis and mortality rates, characterized by cancer cell heterogeneity and complex tumor microenvironment (TME). Single cell biology is an ideal and powerful tool to address these features at a molecular level. However, this approach requires enzymatic cell dissociation that can influence cellular coverage. By contrast, single nucleus RNA sequencing (snRNA-seq) has substantial advantages including compatibility with frozen samples and the elimination of a dissociation-induced, transcriptional stress response. To better profile and understand the functional diversity of different cellular components in melanoma progression, we performed snRNA-seq of 16,839 nuclei obtained from tumor samples along the growth of murine syngeneic melanoma model carrying a BRAFV600E mutation and collected 9 days or 23 days after subcutaneous cell injection. We defined 11 different subtypes of functional cell clusters among malignant cells and 5 different subsets of myeloid cells that display distinct global transcriptional program and different enrichment in early or advanced stage of tumor growth, confirming that this approach was useful to accurately identify intratumor heterogeneity and dynamics during tumor evolution. The current study offers a deep insight into the biology of melanoma highlighting TME reprogramming through tumor initiation and progression, underlying further discovery of new TME biomarkers which may be potentially druggable.

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.

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.

Leveraging gene correlations in single cell transcriptomic data

BACKGROUND

Many approaches have been developed to overcome technical noise in single cell RNA-sequencing (scRNAseq). As researchers dig deeper into data-looking for rare cell types, subtleties of cell states, and details of gene regulatory networks-there is a growing need for algorithms with controllable accuracy and fewer ad hoc parameters and thresholds. Impeding this goal is the fact that an appropriate null distribution for scRNAseq cannot simply be extracted from data in which ground truth about biological variation is unknown (i.e., usually).

RESULTS

We approach this problem analytically, assuming that scRNAseq data reflect only cell heterogeneity (what we seek to characterize), transcriptional noise (temporal fluctuations randomly distributed across cells), and sampling error (i.e., Poisson noise). We analyze scRNAseq data without normalization-a step that skews distributions, particularly for sparse data-and calculate p values associated with key statistics. We develop an improved method for selecting features for cell clustering and identifying gene-gene correlations, both positive and negative. Using simulated data, we show that this method, which we call BigSur (Basic Informatics and Gene Statistics from Unnormalized Reads), captures even weak yet significant correlation structures in scRNAseq data. Applying BigSur to data from a clonal human melanoma cell line, we identify thousands of correlations that, when clustered without supervision into gene communities, align with known cellular components and biological processes, and highlight potentially novel cell biological relationships.

CONCLUSIONS

New insights into functionally relevant gene regulatory networks can be obtained using a statistically grounded approach to the identification of gene-gene correlations.

The JNK and Hippo pathways control epithelial integrity and prevent tumor initiation by regulating an overlapping transcriptome

Epithelial organs maintain their integrity and prevent tumor initiation by actively removing defective cells, such as those that have lost apicobasal polarity. Here, we identify how transcription factors of two key signaling pathways-Jun-N-terminal kinase (JNK) and Hippo-regulate epithelial integrity by controlling transcription of an overlapping set of target genes. Targeted DamID experiments reveal that, in proliferating cells of the Drosophila melanogaster eye, the AP-1 transcription factor Jun and the Hippo pathway transcription regulators Yorkie and Scalloped bind to a common suite of target genes that promote organ growth. In defective neoplastic cells, AP-1 transcription factors repress transcription of growth genes together with the C-terminal binding protein (CtBP) co-repressor. If gene repression by AP-1/CtBP fails, neoplastic tumor growth ensues, driven by Yorkie/Scalloped. Thus, AP-1/CtBP eliminates defective cells and prevents tumor initiation by acting in parallel to Yorkie/Scalloped to repress expression of a shared transcriptome. These findings shed new light on the maintenance of epithelial integrity and tumor suppression.

Efficacy of PARP inhibitor therapy after targeted BRAF/MEK failure in advanced melanoma

Modern advancements in targeted therapy and immunotherapy have significantly improved survival outcomes for advanced melanoma; however, there remains a need for novel approaches to overcome disease progression and treatment resistance. In recent years, PARPi therapy has shown great promise both as a single regimen and in combination with other therapeutics in melanoma. Here, we describe three unique cases of advanced BRAF V600 mutated melanoma that progressed on targeted BRAF/MEK agents that subsequently exhibited partial to near-complete responses to combinatory PARPi and BRAF/MEK inhibitors. This highlights both a potential synergy underlying this combinatory approach and its efficacy as a treatment option for patients with advanced melanoma refractory to targeted and/or immunotherapies. Prospective clinical trials are needed to explore this synergic effect in larger melanoma cohorts to investigate this combination for treating refractory advanced melanoma.

Netrin-1 and UNC5B Cooperate with Integrins to Mediate YAP-Driven Cytostasis

Opposite expression and pro- or anti-cancer function of YAP and its paralog TAZ/WWTR1 stratify cancers into binary YAPon and YAPoff classes. These transcriptional coactivators are oncogenic in YAPon cancers. In contrast, YAP/TAZ are silenced epigenetically along with their integrin and extracellular matrix adhesion target genes in neural and neuroendocrine YAPoff cancers (e.g., small cell lung cancer, retinoblastoma). Forced YAP/TAZ expression induces these targets, causing cytostasis in part through Integrin-αV/β5, independent of the integrin-binding RGD ligand. Other effectors of this anticancer YAP function are unknown. Here, using clustered regularly interspaced short palindromic repeats (CRISPR) screens, we link the Netrin receptor UNC5B to YAP-induced cytostasis in YAPoff cancers. Forced YAP expression induces UNC5B through TEAD DNA-binding partners, as either TEAD1/4-loss or a YAP mutation that disrupts TEAD-binding (S94A) blocks, whereas a TEAD-activator fusion (TEAD(DBD)-VP64) promotes UNC5B induction. Ectopic YAP expression also upregulates UNC5B relatives and their netrin ligands in YAPoff cancers. Netrins are considered protumorigenic, but knockout and peptide/decoy receptor blocking assays reveal that in YAPoff cancers, UNC5B and Netrin-1 can cooperate with integrin-αV/β5 to mediate YAP-induced cytostasis. These data pinpoint an unsuspected Netrin-1/UNC5B/integrin-αV/β5 axis as a critical effector of YAP tumor suppressor activity.

SIGNIFICANCE

Netrins are widely perceived as procancer proteins; however, we uncover an anticancer function for Netrin-1 and its receptor UNC5B.

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.

Malignant Melanoma With Neural Marker-Positive Distant Organ Cancers

Malignant melanoma is a melanocyte-derived tumor known for its aggressive clinical behavior. Melanocytes originate from the neural crest, which also gives rise to neural tissues. Malignant melanoma can occasionally exhibit neural differentiation. We report a case of a 70-year-old male with malignant melanoma exhibiting neural marker positivity in the absence of typical melanoma markers. The patient initially presented with a dark nodule on his left heel, which was confirmed as malignant melanoma through cytology. Surgical resection and lymph node dissection were performed, revealing atypical melanocytes. Despite postoperative nivolumab treatment, metastases in the brain and lungs were observed. Histological examination of the brain tumor showed neural differentiation markers (thyroid transcription factor 1 (TTF-1), cytokeratin 7 (CK7), AE1/AE3, and epidermal growth factor receptor (EGFR)) with negative melanoma markers. The patient eventually succumbed to the disease despite multiple treatments. An autopsy revealed multiple organ tumors (brain, duodenum, stomach, liver, and bile duct) negative for melanoma markers but positive for neuroendocrine markers (CD56, synaptophysin, and chromogranin A). This case suggests two possibilities: the coexistence of malignant melanoma with neuroendocrine tumors or a transformation of melanoma into a neuroendocrine phenotype. This case highlights the need for clinicians to consider the potential for melanoma to lose typical markers and transform into neuroendocrine cancer.

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.

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.

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.