Melanoma plasticity and phenotypic diversity: therapeutic barriers and opportunities

α-Chitosan and β-Oligochitosan Mixtures-Based Formula for In Vitro Assessment of Melanocyte Cells Response

Chitosan is a natural polymer with numerous biomedical applications. The cellular activity of chitosan has been studied in various types of cancer, including melanoma, and indicates that these molecules can open new perspectives on antiproliferative action and anticancer therapy. This study analyzes how different chitosan conformations, such as α-chitosan (CH) or β-oligochitosan (CO), with various degrees of deacetylation (DDA) and molar mass (MM), both in different concentrations and in CH-CO mixtures, influence the cellular processes of SK-MEL-28 melanocytes, to estimate the reactivity of these cells to the applied treatments. The in vitro evaluation was carried out, aiming at the cellular metabolism (MTT assay), cellular morphology, and chitinase-like glycoprotein YKL-40 expression. The in vitro effect of the CH-CO mixture application on melanocytes is obvious at low concentrations of α-chitosan/β-oligochitosan (1:2 ratio), with the cell's response supporting the hypothesis that β-oligo-chitosan amplifies the effect. This oligochitosan mixture, favored by the β conformation and its small size, penetrates faster into the cells, being more reactive when interacting with some cellular components. Morphological effects expressed by the loss of cell adhesion and the depletion of YKL-40 synthesis are significant responses of melanocytes. β-oligochitosan (1.5 kDa) induces an extension of cytophysiological effects and limits the cell viability compared to α-chitosan (400-900 kDa). Statistical analysis using multivariate techniques showed differences between the CH samples and CH-CO mixtures.

AMBRA1 levels predict resistance to MAPK inhibitors in melanoma

Intrinsic and acquired resistance to mitogen-activated protein kinase inhibitors (MAPKi) in melanoma remains a major therapeutic challenge. Here, we show that the clinical development of resistance to MAPKi is associated with reduced tumor expression of the melanoma suppressor Autophagy and Beclin 1 Regulator 1 (AMBRA1) and that lower expression levels of AMBRA1 predict a poor response to MAPKi treatment. Functional analyses show that loss of AMBRA1 induces phenotype switching and orchestrates an extracellular signal-regulated kinase (ERK)-independent resistance mechanism by activating focal adhesion kinase 1 (FAK1). In both in vitro and in vivo settings, melanomas with low AMBRA1 expression exhibit intrinsic resistance to MAPKi therapy but higher sensitivity to FAK1 inhibition. Finally, we show that the rapid development of resistance in initially MAPKi-sensitive melanomas can be attributed to preexisting subclones characterized by low AMBRA1 expression and that cotreatment with MAPKi and FAK1 inhibitors (FAKi) effectively prevents the development of resistance in these tumors. In summary, our findings underscore the value of AMBRA1 expression for predicting melanoma response to MAPKi and supporting the therapeutic efficacy of FAKi to overcome MAPKi-induced resistance.

Valine aminoacyl-tRNA synthetase promotes therapy resistance in melanoma

Transfer RNA dynamics contribute to cancer development through regulation of codon-specific messenger RNA translation. Specific aminoacyl-tRNA synthetases can either promote or suppress tumourigenesis. Here we show that valine aminoacyl-tRNA synthetase (VARS) is a key player in the codon-biased translation reprogramming induced by resistance to targeted (MAPK) therapy in melanoma. The proteome rewiring in patient-derived MAPK therapy-resistant melanoma is biased towards the usage of valine and coincides with the upregulation of valine cognate tRNAs and of VARS expression and activity. Strikingly, VARS knockdown re-sensitizes MAPK-therapy-resistant patient-derived melanoma in vitro and in vivo. Mechanistically, VARS regulates the messenger RNA translation of valine-enriched transcripts, among which hydroxyacyl-CoA dehydrogenase mRNA encodes for a key enzyme in fatty acid oxidation. Resistant melanoma cultures rely on fatty acid oxidation and hydroxyacyl-CoA dehydrogenase for their survival upon MAPK treatment. Together, our data demonstrate that VARS may represent an attractive therapeutic target for the treatment of therapy-resistant melanoma.

A Spatial Transcriptomics Browser for Discovering Gene Expression Landscapes across Microscopic Tissue Sections

A crucial feature of life is its spatial organization and compartmentalization on the molecular, cellular, and tissue levels. Spatial transcriptomics (ST) technology has opened a new chapter of the sequencing revolution, emerging rapidly with transformative effects across biology. This technique produces extensive and complex sequencing data, raising the need for computational methods for their comprehensive analysis and interpretation. We developed the ST browser web tool for the interactive discovery of ST images, focusing on different functional aspects such as single gene expression, the expression of functional gene sets, as well as the inspection of the spatial patterns of cell-cell interactions. As a unique feature, our tool applies self-organizing map (SOM) machine learning to the ST data. Our SOM data portrayal method generates individual gene expression landscapes for each spot in the ST image, enabling its downstream analysis with high resolution. The performance of the spatial browser is demonstrated by disentangling the intra-tumoral heterogeneity of melanoma and the microarchitecture of the mouse brain. The integration of machine-learning-based SOM portrayal into an interactive ST analysis environment opens novel perspectives for the comprehensive knowledge mining of the organization and interactions of cellular ecosystems.

Targeting HIF-1α with Specific DNA Yokes for Effective Anticancer Therapy

Hypoxia, a ubiquitous hallmark in cancer, underscores the significance of targeting HIF-1α, the principal transcriptional factor of hypoxic responses, for effective cancer therapy. Herein, DNA yokes, a novel class of DNA nanomaterials harboring specific HIF-1α binding sequences (hypoxia response elements, HREs), are introduced as nanopharmaceuticals for cancer treatment. Comprising a basal tetrahedral DNA nanostructure and four HRE-bearing overhanging chains, DNA yokes exhibit exceptional stability and prolonged intracellular retention. The investigation reveals their capacity to bind HIF-1α, thereby disrupting its interaction with the downstream genomic DNAs and impeding transcriptional activity. Moreover, DNA yokes facilitate HIF-1α degradation via the ubiquitination pathway, thereby sequestering it from downstream targets and ultimately promoting its degradation. In addition, DNA yokes attenuate cancer cell proliferation, migration, and invasion under hypoxic conditions, while also displaying preferential accumulation within tumors, thereby inhibiting tumor growth and metastasis in vivo. This study pioneers a novel approach to cancer therapy through the development of DNA-based drugs characterized by high stability and low toxicity to normal cells, positioning DNA yokes as promising candidates for cancer treatment.

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.

Modeling Melanoma Heterogeneity In Vitro: Redox, Resistance and Pigmentation Profiles

Microenvironment and transcriptional plasticity generate subpopulations within the tumor, and the use of BRAF inhibitors (BRAFis) contributes to the rise and selection of resistant clones. We stochastically isolated subpopulations (C1, C2, and C3) from naïve melanoma and found that the clones demonstrated distinct morphology, phenotypic, and functional profiles: C1 was less proliferative, more migratory and invasive, less sensitive to BRAFis, less dependent on OXPHOS, more sensitive to oxidative stress, and less pigmented; C2 was more proliferative, less migratory and invasive, more sensitive to BRAFis, less sensitive to oxidative stress, and more pigmented; and C3 was less proliferative, more migratory and invasive, less sensitive to BRAFis, more dependent on OXPHOS, more sensitive to oxidative stress, and more pigmented. Hydrogen peroxide plays a central role in oxidative stress and cell signaling, and PRDXs are one of its main consumers. The intrinsically resistant C1 and C3 clones had lower MITF, PGC-1α, and PRDX1 expression, while C1 had higher AXL and decreased pigmentation markers, linking PRDX1 to clonal heterogeneity and resistance. PRDX2 is depleted in acquired BRAFi-resistant cells and acts as a redox sensor. Our results illustrate that decreased pigmentation markers are related to therapy resistance and decreased antioxidant defense.

Stochastic modelling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present the first formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Interestingly, sublines previously observed to be resistant to anti-CTLA-4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression.

Molecular patterns of resistance to immune checkpoint blockade in melanoma

Immune checkpoint blockade (ICB) has improved outcome for patients with metastatic melanoma but not all benefit from treatment. Several immune- and tumor intrinsic features are associated with clinical response at baseline. However, we need to further understand the molecular changes occurring during development of ICB resistance. Here, we collect biopsies from a cohort of 44 patients with melanoma after progression on anti-CTLA4 or anti-PD1 monotherapy. Genetic alterations of antigen presentation and interferon gamma signaling pathways are observed in approximately 25% of ICB resistant cases. Anti-CTLA4 resistant lesions have a sustained immune response, including immune-regulatory features, as suggested by multiplex spatial and T cell receptor (TCR) clonality analyses. One anti-PD1 resistant lesion harbors a distinct immune cell niche, however, anti-PD1 resistant tumors are generally immune poor with non-expanded TCR clones. Such immune poor microenvironments are associated with melanoma cells having a de-differentiated phenotype lacking expression of MHC-I molecules. In addition, anti-PD1 resistant tumors have reduced fractions of PD1+ CD8+ T cells as compared to ICB naïve metastases. Collectively, these data show the complexity of ICB resistance and highlight differences between anti-CTLA4 and anti-PD1 resistance that may underlie differential clinical outcomes of therapy sequence and combination.

Purinergic Ca2+ signaling as a novel mechanism of drug tolerance in BRAF mutant melanoma

Drug tolerance is a major cause of relapse after cancer treatment. In spite of intensive efforts1-9, its molecular basis remains poorly understood, hampering actionable intervention. We report a previously unrecognized signaling mechanism supporting drug tolerance in BRAF-mutant melanoma treated with BRAF inhibitors that could be of general relevance to other cancers. Its key features are cell-intrinsic intracellular Ca2+ signaling initiated by P2X7 receptors (purinergic ligand-gated cation channels), and an enhanced ability for these Ca2+ signals to reactivate ERK1/2 in the drug-tolerant state. Extracellular ATP, virtually ubiquitous in living systems, is the ligand that can initiate Ca2+ spikes via P2X7 channels. ATP is abundant in the tumor microenvironment and is released by dying cells, ironically implicating treatment-initiated cancer cell death as a source of trophic stimuli that leads to ERK reactivation and drug tolerance. Such a mechanism immediately offers an explanation of the inevitable relapse after BRAFi treatment in BRAF-mutant melanoma, and points to actionable strategies to overcome it.

BRAF - a tumour-agnostic drug target with lineage-specific dependencies

In June 2022, the FDA granted Accelerated Approval to the BRAF inhibitor dabrafenib in combination with the MEK inhibitor trametinib for the treatment of adult and paediatric patients (≥6 years of age) with unresectable or metastatic BRAFV600E-mutant solid tumours, except for BRAFV600E-mutant colorectal cancers. The histology-agnostic approval of dabrafenib plus trametinib marks the culmination of two decades of research into the landscape of BRAF mutations in human cancers, the biochemical mechanisms underlying BRAF-mediated tumorigenesis, and the clinical development of selective RAF and MEK inhibitors. Although the majority of patients with BRAFV600E-mutant tumours derive clinical benefit from BRAF inhibitor-based combinations, resistance to treatment develops in most. In this Review, we describe the biochemical basis for oncogenic BRAF-induced activation of MAPK signalling and pan-cancer and lineage-specific mechanisms of intrinsic, adaptive and acquired resistance to BRAF inhibitors. We also discuss novel RAF inhibitors and drug combinations designed to delay the emergence of treatment resistance and/or expand the population of patients with BRAF-mutant cancers who benefit from molecularly targeted therapies.

Deciphering cellular plasticity in pancreatic cancer for effective treatments

Cellular plasticity and therapy resistance are critical features of pancreatic cancer, a highly aggressive and fatal disease. The pancreas, a vital organ that produces digestive enzymes and hormones, is often affected by two main types of cancer: the pre-dominant ductal adenocarcinoma and the less common neuroendocrine tumors. These cancers are difficult to treat due to their complex biology characterized by cellular plasticity leading to therapy resistance. Cellular plasticity refers to the capability of cancer cells to change and adapt to different microenvironments within the body which includes acinar-ductal metaplasia, epithelial to mesenchymal/epigenetic/metabolic plasticity, as well as stemness. This plasticity allows heterogeneity of cancer cells, metastasis, and evasion of host's immune system and develops resistance to radiation, chemotherapy, and targeted therapy. To overcome this resistance, extensive research is ongoing exploring the intrinsic and extrinsic factors through cellular reprogramming, chemosensitization, targeting metabolic, key survival pathways, etc. In this review, we discussed the mechanisms of cellular plasticity involving cellular adaptation and tumor microenvironment and provided a comprehensive understanding of its role in therapy resistance and ways to overcome it.

Dormancy of cutaneous melanoma

Many cancer-related deaths including melanoma result from metastases that develop months or years after the initial cancer therapy. Even the most effective drugs and immune therapies rarely eradicate all tumor cells. Instead, they strongly reduce cancer burden, permitting dormant cancer cells to persist in niches, where they establish a cellular homeostasis with their host without causing clinical symptoms. Dormant cancers respond poorly to most drugs and therapies since they do not proliferate and hide in niches. It therefore remains a major challenge to develop novel therapies for dormant cancers. In this review we focus on the mechanisms regulating the initiation of cutaneous melanoma dormancy as well as those which are involved in reawakening of dormant cutaneous melanoma cells. In recent years the role of neutrophils and niche components in reawakening of melanoma cells came into focus and indicate possible future therapeutic applications. Sophisticated in vitro and in vivo melanoma dormancy models are needed to make progress in this field and are discussed.

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.

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.

Cancer Cell Biomechanical Properties Accompany Tspan8-Dependent Cutaneous Melanoma Invasion

The intrinsic biomechanical properties of cancer cells remain poorly understood. To decipher whether cell stiffness modulation could increase melanoma cells' invasive capacity, we performed both in vitro and in vivo experiments exploring cell stiffness by atomic force microscopy (AFM). We correlated stiffness properties with cell morphology adaptation and the molecular mechanisms underlying epithelial-to-mesenchymal (EMT)-like phenotype switching. We found that melanoma cell stiffness reduction was systematically associated with the acquisition of invasive properties in cutaneous melanoma cell lines, human skin reconstructs, and Medaka fish developing spontaneous MAP-kinase-induced melanomas. We observed a systematic correlation of stiffness modulation with cell morphological changes towards mesenchymal characteristic gains. We accordingly found that inducing melanoma EMT switching by overexpressing the ZEB1 transcription factor, a major regulator of melanoma cell plasticity, was sufficient to decrease cell stiffness and transcriptionally induce tetraspanin-8-mediated dermal invasion. Moreover, ZEB1 expression correlated with Tspan8 expression in patient melanoma lesions. Our data suggest that intrinsic cell stiffness could be a highly relevant marker for human cutaneous melanoma development.

CD133-Dependent Activation of Phosphoinositide 3-Kinase /AKT/Mammalian Target of Rapamycin Signaling in Melanoma Progression and Drug Resistance

Melanoma frequently harbors genetic alterations in key molecules leading to the aberrant activation of PI3K and its downstream pathways. Although the role of PI3K/AKT/mTOR in melanoma progression and drug resistance is well documented, targeting the PI3K/AKT/mTOR pathway showed less efficiency in clinical trials than might have been expected, since the suppression of the PI3K/mTOR signaling pathway-induced feedback loops is mostly associated with the activation of compensatory pathways such as MAPK/MEK/ERK. Consequently, the development of intrinsic and acquired resistance can occur. As a solid tumor, melanoma is notorious for its heterogeneity. This can be expressed in the form of genetically divergent subpopulations including a small fraction of cancer stem-like cells (CSCs) and non-cancer stem cells (non-CSCs) that make the most of the tumor mass. Like other CSCs, melanoma stem-like cells (MSCs) are characterized by their unique cell surface proteins/stemness markers and aberrant signaling pathways. In addition to its function as a robust marker for stemness properties, CD133 is crucial for the maintenance of stemness properties and drug resistance. Herein, the role of CD133-dependent activation of PI3K/mTOR in the regulation of melanoma progression, drug resistance, and recurrence is reviewed.

Mechanisms of Melanoma Progression and Treatment Resistance: Role of Cancer Stem-like Cells

Melanoma is the third most common type of skin cancer, characterized by its heterogeneity and propensity to metastasize to distant organs. Melanoma is a heterogeneous tumor, composed of genetically divergent subpopulations, including a small fraction of melanoma-initiating cancer stem-like cells (CSCs) and many non-cancer stem cells (non-CSCs). CSCs are characterized by their unique surface proteins associated with aberrant signaling pathways with a causal or consequential relationship with tumor progression, drug resistance, and recurrence. Melanomas also harbor significant alterations in functional genes (BRAF, CDKN2A, NRAS, TP53, and NF1). Of these, the most common are the BRAF and NRAS oncogenes, with 50% of melanomas demonstrating the BRAF mutation (BRAFV600E). While the successful targeting of BRAFV600E does improve overall survival, the long-term efficacy of available therapeutic options is limited due to adverse side effects and reduced clinical efficacy. Additionally, drug resistance develops rapidly via mechanisms involving fast feedback re-activation of MAPK signaling pathways. This article updates information relevant to the mechanisms of melanoma progression and resistance and particularly the mechanistic role of CSCs in melanoma progression, drug resistance, and recurrence.

Cell state dependent effects of Bmal1 on melanoma immunity and tumorigenicity

The circadian clock regulator Bmal1 modulates tumorigenesis, but its reported effects are inconsistent. Here, we show that Bmal1 has a context-dependent role in mouse melanoma tumor growth. Loss of Bmal1 in YUMM2.1 or B16-F10 melanoma cells eliminates clock function and diminishes hypoxic gene expression and tumorigenesis, which could be rescued by ectopic expression of HIF1α in YUMM2.1 cells. By contrast, over-expressed wild-type or a transcriptionally inactive mutant Bmal1 non-canonically sequester myosin heavy chain 9 (Myh9) to increase MRTF-SRF activity and AP-1 transcriptional signature, and shift YUMM2.1 cells from a Sox10high to a Sox9high immune resistant, mesenchymal cell state that is found in human melanomas. Our work describes a link between Bmal1, Myh9, mouse melanoma cell plasticity, and tumor immunity. This connection may underlie cancer therapeutic resistance and underpin the link between the circadian clock, MRTF-SRF and the cytoskeleton.

Differentiation States of Phenotypic Transition of Melanoma Cells Are Revealed by 3D Cell Cultures

Melanoma is characterized by high metastatic potential favored by the epithelial-to-mesenchymal transition (EMT), leading melanoma cells to exhibit a spectrum of typical EMT markers. This study aimed to analyze the expression of EMT markers in A375 and BLM melanoma cell lines cultured in 2D monolayers and 3D spheroids using morphological and molecular methods. The expression of EMT markers was strongly affected by 3D arrangement and revealed a hybrid phenotype for the two cell lines. Indeed, although E-cadherin was almost undetectable in both A375 and BLM cells, cortical actin was detected in A375 2D monolayers and 3D spheroids and was strongly expressed in BLM 3D spheroids. The mesenchymal marker N-cadherin was significantly up-regulated in A375 3D spheroids while undetectable in BLM cells, but vimentin was similarly expressed in both cell lines at the gene and protein levels. This pattern suggests that A375 cells exhibit a more undifferentiated/mesenchymal phenotype, while BLM cells have more melanocytic/differentiated characteristics. Accordingly, the Zeb1 and 2, Slug, Snail and Twist gene expression analyses showed that they were differentially expressed in 2D monolayers compared to 3D spheroids, supporting this view. Furthermore, A375 cells are characterized by a greater invasive potential, strongly influenced by 3D arrangement, compared to the BLM cell line, as evaluated by SDS-zymography and TIMPs gene expression analysis. Finally, TGF-β1, a master controller of EMT, and lysyl oxidase (LOX), involved in melanoma progression, were strongly up-regulated by 3D arrangement in the metastatic BLM cells alone, likely playing a role in the metastatic phases of melanoma progression. Overall, these findings suggest that A375 and BLM cells possess a hybrid/intermediate phenotype in relation to the expression of EMT markers. The former is characterized by a more mesenchymal/undifferentiated phenotype, while the latter shows a more melanocytic/differentiated phenotype. Our results contribute to the characterization of the role of EMT in melanoma cells and confirm that a 3D cell culture model could provide deeper insight into our understanding of the biology of melanoma.

Uncovering Minimal Pathways in Melanoma Initiation

Cutaneous melanomas are clinically and histologically heterogeneous. Most display activating mutations in Braf or Nras and complete loss of function of one or more tumor suppressor genes. Mouse models that replicate such mutations produce fast-growing, pigmented tumors. However, mice that combine Braf activation with only heterozygous loss of Pten also produce tumors and, as we show here, in an Albino background this occurs even with Braf activation alone. Such tumors arise rarely, grow slowly, and express low levels of pigmentation genes. The timing of their appearance was consistent with a single step stochastic event, but no evidence could be found that it required de novo mutation, suggesting instead the involvement of an epigenetic transition. Single-cell transcriptomic analysis revealed such tumors to be heterogeneous, including a minor cell type we term LNM ( L ow-pigment, N eural- and extracellular M atrix-signature) that displays gene expression resembling "neural crest"-like cell subsets detected in the fast-growing tumors of more heavily-mutated mice, as well as in human biopsy and xenograft samples. We provide evidence that LNM cells pre-exist in normal skin, are expanded by Braf activation, can transition into malignant cells, and persist with malignant cells through multiple rounds of transplantation. We discuss the possibility that LNM cells not only serve as a pre-malignant state in the production of some melanomas, but also as an important intermediate in the development of drug 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.

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.

Assessment of Autophagy in Tumor Cells of Human Skin Melanoma of Different Stages

High levels of autophagy can increase the viability of tumor cells as well as their resistance to chemotherapy. Evaluation of the dynamics of autophagy processes at different stages of carcinogenesis can extend our understanding of melanoma pathogenesis to develop new therapeutic approaches. We performed a comparative study of tumor cell autophagy in stages II and III human skin melanoma. Tumor cells were characterized by high content of structures associated with autophagy (autophagosomes and autolysosomes). In stage III melanoma characterized by the presence of regional metastases in the lymph nodes, tumor cells showed higher expression of the autophagy marker protein LC3beta in comparison with stage II melanoma cells, which can indicate the involvement of autophagy processes in tumor progression and the formation of metastases in the lymph nodes.

Master Transcription Factor Reprogramming Unleashes Selective Translation Promoting Castration Resistance and Immune Evasion in Lethal Prostate Cancer

Signaling rewiring allows tumors to survive therapy. Here we show that the decrease of the master regulator microphthalmia transcription factor (MITF) in lethal prostate cancer unleashes eukaryotic initiation factor 3B (eIF3B)-dependent translation reprogramming of key mRNAs conferring resistance to androgen deprivation therapy (ADT) and promoting immune evasion. Mechanistically, MITF represses through direct promoter binding eIF3B, which in turn regulates the translation of specific mRNAs. Genome-wide eIF3B enhanced cross-linking immunoprecipitation sequencing (eCLIP-seq) showed specialized binding to a UC-rich motif present in subsets of 5' untranslated regions. Indeed, translation of the androgen receptor and major histocompatibility complex I (MHC-I) through this motif is sensitive to eIF3B amount. Notably, pharmacologic targeting of eIF3B-dependent translation in preclinical models sensitizes prostate cancer to ADT and anti-PD-1 therapy. These findings uncover a hidden connection between transcriptional and translational rewiring promoting therapy-refractory lethal prostate cancer and provide a druggable mechanism that may transcend into effective combined therapeutic strategies.

SIGNIFICANCE

Our study shows that specialized eIF3B-dependent translation of specific mRNAs released upon downregulation of the master transcription factor MITF confers castration resistance and immune evasion in lethal prostate cancer. Pharmacologic targeting of this mechanism delays castration resistance and increases immune-checkpoint efficacy. This article is featured in Selected Articles from This Issue, p. 2489.

Has-miR-300-GADD45B promotes melanoma growth via cell cycle

Response to oncogenic factors like UV, GADD45 family in skin participates in scavenging ROS, DNA repair and cell cycle control. Because of this, the previous study of the chronic UVB injury model has found that hsa-miR-300 can conduct intercellular transport by exosomes and target regulation of GADD45B. Whether the hsa-miR-300-GADD45B still regulates tumor development by cell cycle pathway is unclear. Through transcriptomic analysis of primary (n=39) and metastatic (n=102) melanoma, it was confirmed that in metastatic samples, some of the 97 down-regulated genes participate in maintaining skin homeostasis while 42 up-regulated genes were enriched in cancer-related functions. Furthermore, CDKN1A, CDKN2A, CXCR4 and RAD51 in the melanoma pathway, were also differentially expressed between normal skin and melanoma. CDKN1A and CDKN2A were also found to be involved in TP53-dependent cell cycle regulation. In conclusion, it was speculated that CDKN1A, CDKN2A, TP53, GADD45B and hsa-miR-300 may have regulatory relationships. It was demonstrated that there is a bidirectional regulation between hsa-miR-300 and TP53. In addition, miR-300 can regulate CDKN1A by GADD45B/TP53 and promote melanoma growth by accelerating the cell cycle transition from G1/S to G2 phase.

Ex vivo mapping of enhancer networks that define the transcriptional program driving melanoma metastasis

Mortality from vmelanoma is associated with metastatic disease, but the mechanisms leading to spreading of the cancer cells remain obscure. Spatial profiling revealed that melanoma is characterized by a high degree of heterogeneity, which is established by the ability of melanoma cells to switch between different phenotypical stages. This plasticity, likely a heritage from embryonic pathways, accounts for a relevant part of the metastatic potential of these lesions, and requires the rapid and efficient reorganization of the transcriptional landscape of melanoma cells. A large part of the non-coding genome cooperates to control gene expression, specifically through the activity of enhancers (ENHs). In this study, we aimed to identify ex vivo the network of active ENHs and to outline their cooperative interactions in supporting transcriptional adaptation during melanoma metastatic progression. We conducted a genome-wide analysis to map active ENHs distribution in a retrospective cohort of 39 melanoma patients, comparing the profiles obtained in primary (N = 19) and metastatic (N = 20) melanoma lesions. Unsupervised clustering showed that the profile for acetylated histone H3 at lysine 27 (H3K27ac) efficiently segregates lesions into three different clusters corresponding to progressive stages of the disease. We reconstructed the map of super-ENHs (SEs) and cooperative ENHs that associate with metastatic progression in melanoma, which showed that cooperation among regulatory elements is a mandatory requirement for transcriptional plasticity. We also showed that these elements carry out specialized and non-redundant functions, and indicated the existence of a hierarchical organization, with SEs on top as masterminds of the entire transcriptional program and classical ENHs as executors. By providing an innovative vision of how the chromatin landscape of melanoma works during metastatic spreading, our data also point out the need to integrate functional profiling in the analysis of cancer lesions to increase definition and improve interpretation of tumor heterogeneity.

Ferroptosis: mechanisms and implications for cancer development and therapy response

As cancer cells develop resistance to apoptosis, non-apoptotic cell death modalities, such as ferroptosis, have emerged as promising strategies to combat therapy-resistant cancers. Cells that develop resistance to conventional therapies or metastatic cancer cells have been shown to have increased sensitivity to ferroptosis. Therefore, targeting the regulatory elements of ferroptosis in cancer could offer novel therapeutic opportunities. In this review, we first provide an overview of the known ferroptosis regulatory networks and discuss recent findings on how they contribute to cancer plasticity. We then expand into the critical role of selenium metabolism in regulating ferroptosis. Finally, we highlight specific cases where induction of ferroptosis could be used to sensitize cancer cells to this form of cell death.

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.

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.

Targeting ferroptosis in melanoma: cancer therapeutics

Melanoma is an aggressive kind of skin cancer; its rate has risen rapidly over the past few decades. Melanoma reports for only about 1% of skin cancers but leads to a high majority of skin cancer deaths. Thus, new useful therapeutic approaches are currently required, to state effective treatments to consistently enhance the overall survival rate of melanoma patients. Ferroptosis is a recently identified cell death process, which is different from autophagy, apoptosis, necrosis, and pyroptosis in terms of biochemistry, genetics, and morphology which plays an important role in cancer treatment. Ferroptosis happens mostly by accumulating iron and lipid peroxides in the cell. Recently, studies have revealed that ferroptosis has a key role in the tumor's progression. Especially, inducing ferroptosis in cells can inhibit the tumor cells' growth, leading to back warding tumorigenesis. Here, we outline the ferroptosis characteristics from its basic role in melanoma cancer and mention its possible applications in melanoma cancer treatment. Video Abstract.

Data-driven structural analysis of small cell lung cancer transcription factor network suggests potential subtype regulators and transition pathways

Small cell lung cancer (SCLC) is an aggressive disease and challenging to treat due to its mixture of transcriptional subtypes and subtype transitions. Transcription factor (TF) networks have been the focus of studies to identify SCLC subtype regulators via systems approaches. Yet, their structures, which can provide clues on subtype drivers and transitions, are barely investigated. Here, we analyze the structure of an SCLC TF network by using graph theory concepts and identify its structurally important components responsible for complex signal processing, called hubs. We show that the hubs of the network are regulators of different SCLC subtypes by analyzing first the unbiased network structure and then integrating RNA-seq data as weights assigned to each interaction. Data-driven analysis emphasizes MYC as a hub, consistent with recent reports. Furthermore, we hypothesize that the pathways connecting functionally distinct hubs may control subtype transitions and test this hypothesis via network simulations on a candidate pathway and observe subtype transition. Overall, structural analyses of complex networks can identify their functionally important components and pathways driving the network dynamics. Such analyses can be an initial step for generating hypotheses and can guide the discovery of target pathways whose perturbation may change the network dynamics phenotypically.

A human mucosal melanoma organoid platform for modeling tumor heterogeneity and exploring immunotherapy combination options

Mucosal melanoma (MM), an aggressive rare subtype of melanoma, is distinct from cutaneous melanoma and has poor prognoses. We addressed the lack of cell models for MM by establishing 30 organoids of human oral MM (OMM), which retained major histopathological and functional features of parental tumors. Organoid groups derived from chronologically or intratumorally distinct lesions within the same individual displayed heterogeneous genetics, expression profiles, and drug responses, indicating rapid tumor evolution and poor clinical response. Furthermore, transcriptome analysis revealed receptor tyrosine kinases (RTKs) signaling, particularly NGFR, a nerve growth factor receptor, was significantly up-regulated in OMMs and organoids from patients resistant to anti-programmed cell death protein 1 (anti-PD-1) therapy. Combining anti-PD-1 with anlotinib (a phase 2 multitarget RTK inhibitor for OMM) or NGFR knockdown enhanced the effective activity of immune cells in organoid-immune cell coculture systems. Together, our study suggested that OMM organoids serve as faithful models for exploring tumor evolution and immunotherapy combination strategies.

Application of single-cell sequencing to the research of tumor microenvironment

Single-cell sequencing is a technique for detecting and analyzing genomes, transcriptomes, and epigenomes at the single-cell level, which can detect cellular heterogeneity lost in conventional sequencing hybrid samples, and it has revolutionized our understanding of the genetic heterogeneity and complexity of tumor progression. Moreover, the tumor microenvironment (TME) plays a crucial role in the formation, development and response to treatment of tumors. The application of single-cell sequencing has ushered in a new age for the TME analysis, revealing not only the blueprint of the pan-cancer immune microenvironment, but also the heterogeneity and differentiation routes of immune cells, as well as predicting tumor prognosis. Thus, the combination of single-cell sequencing and the TME analysis provides a unique opportunity to unravel the molecular mechanisms underlying tumor development and progression. In this review, we summarize the recent advances in single-cell sequencing and the TME analysis, highlighting their potential applications in cancer research and clinical translation.

Multi-modal transcriptomic analysis unravels enrichment of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. It manifests along multiple phenotypic axes and decoding the interconnections among these different axes is crucial to understand its molecular origins and to develop novel therapeutic strategies to control it. Here, we use multi-modal transcriptomic data analysis - bulk, single-cell and spatial transcriptomics - from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity - two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. These patterns were inherent in methylation profiles, suggesting an epigenetic crosstalk between EMT and lineage plasticity in breast cancer. Mathematical modelling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes recapitulate and thus elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and to identify possible interventions to restrict it.

GREB1 isoform 4 is specifically transcribed by MITF and required for melanoma proliferation

Growth regulation by estrogen in breast cancer 1 (GREB1) is involved in hormone-dependent and -independent tumor development (e.g., hepatoblastoma). In this study, we found that a GREB1 splicing variant, isoform 4 (Is4), which encodes C-terminal half of full-length GREB1, is specifically expressed via microphthalmia-associated transcription factor (MITF) in melanocytic melanoma, and that two MITF-binding E-box CANNTG motifs at the 5'-upstream region of GREB1 exon 19 are necessary for GREB1 Is4 transcription. MITF and GREB1 Is4 were strongly co-expressed in approximately 20% of the melanoma specimens evaluated (17/89 cases) and their expression was associated with tumor thickness. GREB1 Is4 silencing reduced melanoma cell proliferation in association with altered expression of cell proliferation-related genes in vitro. In addition, GREB1 Is4 targeting by antisense oligonucleotide (ASO) decreased melanoma xenograft tumor formation and GREB1 Is4 expression in a BRAFV600E; PTENflox melanoma mouse model promoted melanoma formation, demonstrating the crucial role of GREB1 Is4 for melanoma proliferation in vivo. GREB1 Is4 bound to CAD, the rate-limiting enzyme of pyrimidine metabolism, and metabolic flux analysis revealed that GREBI Is4 is necessary for pyrimidine synthesis. These results suggest that MITF-dependent GREB1 Is4 expression leads to melanoma proliferation and GREB1 Is4 represents a new molecular target in melanoma.

Pro-Apoptotic Activity of MCL-1 Inhibitor in Trametinib-Resistant Melanoma Cells Depends on Their Phenotypes and Is Modulated by Reversible Alterations Induced by Trametinib Withdrawal

BACKGROUND

Although BRAFV600/MEK inhibitors improved the treatment of melanoma patients, resistance is acquired almost inevitably.

METHODS

Trametinib withdrawal/rechallenge and MCL-1 inhibition in trametinib-resistance models displaying distinct p-ERK1/2 levels were investigated.

RESULTS

Trametinib withdrawal/rechallenge caused reversible changes in ERK1/2 activity impacting the balance between pro-survival and pro-apoptotic proteins. Reversible alterations were found in MCL-1 levels and MCL-1 inhibitors, BIM and NOXA. Taking advantage of melanoma cell dependency on MCL-1 for survival, we used S63845. While it was designed to inhibit MCL-1 activity, we showed that it also significantly reduced NOXA levels. S63845-induced apoptosis was detected as the enhancement of Annexin V-positivity, caspase-3/7 activation and histone H2AX phosphorylation. Percentages of Annexin V-positive cells were increased most efficiently in trametinib-resistant melanoma cells displaying the p-ERK1/2low/MCL-1low/BIMhigh/NOXAlow phenotype with EC50 values at concentrations as low as 0.1 μM. Higher ERK1/2 activity associated with increased MCL-1 level and reduced BIM level limited pro-apoptotic activity of S63845 further influenced by a NOXA level.

CONCLUSIONS

Our study supports the notion that the efficiency of an agent designed to target a single protein can largely depend on the phenotype of cancer cells. Thus, it is important to define appropriate phenotype determinants to stratify the patients for the novel therapy.

Acetylation reprograms MITF target selectivity and residence time

The ability of transcription factors to discriminate between different classes of binding sites associated with specific biological functions underpins effective gene regulation in development and homeostasis. How this is achieved is poorly understood. The microphthalmia-associated transcription factor MITF is a lineage-survival oncogene that plays a crucial role in melanocyte development and melanoma. MITF suppresses invasion, reprograms metabolism and promotes both proliferation and differentiation. How MITF distinguishes between differentiation and proliferation-associated targets is unknown. Here we show that compared to many transcription factors MITF exhibits a very long residence time which is reduced by p300/CBP-mediated MITF acetylation at K206. While K206 acetylation also decreases genome-wide MITF DNA-binding affinity, it preferentially directs DNA binding away from differentiation-associated CATGTG motifs toward CACGTG elements. The results reveal an acetylation-mediated switch that suppresses differentiation and provides a mechanistic explanation of why a human K206Q MITF mutation is associated with Waardenburg syndrome.

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.

The CD73 is induced by TGF-β1 triggered by nutrient deprivation and highly expressed in dedifferentiated human melanoma

CD73 is the key enzyme in the generation of extracellular adenosine, a mediator involved in tumor progression, tumor immune escape and resistance to anti-cancer therapeutics. Microenvironmental conditions influence the expression of CD73 in tumor cells. However how CD73 expression and activity is regulated in a stress condition of lower nutrient availability are largely unknown. Our results indicate that serum starvation leads to a marked up-regulation of CD73 expression on A375 melanoma cells in a time-dependent manner. The cell-surface expression of CD73 is associated with an increased release of TGF-β1 by starved cells. Blockade of TGF-β1 receptors or TGFβ/SMAD3 signaling pathway significantly reduce the expression of CD73 induced by starvation. Treatment of cells with rTGF-β1 up-regulates the expression of CD73 in a concentration-dependent manner, confirming the role of this pathway in regulating CD73 in melanoma A375 cells. The increased expression of CD73 is associated with enhanced AMPase activity, which is selectively reduced by inhibitors of CD73 activity, APCP and PSB-12489. Pharmacological blockade of CD73 significantly inhibits invasion of melanoma cells in a transwell system. Furthermore, using multiplex immunofluorescence imaging we found that, within human melanoma metastases, tumor cells at the dedifferentiated stage show the highest CD73 protein expression. In summary, our data provide new insights into the mechanism regulating the expression/activity of CD73 in melanoma cells in a condition of lower availability of nutrients, which is a common feature of the tumor microenvironment. Within human metastatic melanoma tissues elevated protein expression of CD73 is associated with an invasive-like phenotype.

BRAF Mutations in Melanoma: Biological Aspects, Therapeutic Implications, and Circulating Biomarkers

Melanoma is an aggressive form of skin cancer resulting from the malignant transformation of melanocytes. Recent therapeutic approaches, including targeted therapy and immunotherapy, have improved the prognosis and outcome of melanoma patients. BRAF is one of the most frequently mutated oncogenes recognised in melanoma. The most frequent oncogenic BRAF mutations consist of a single point mutation at codon 600 (mostly V600E) that leads to constitutive activation of the BRAF/MEK/ERK (MAPK) signalling pathway. Therefore, mutated BRAF has become a useful target for molecular therapy and the use of BRAF kinase inhibitors has shown promising results. However, several resistance mechanisms invariably develop leading to therapeutic failure. The aim of this manuscript is to review the role of BRAF mutational status in the pathogenesis of melanoma and its impact on differentiation and inflammation. Moreover, this review focuses on the mechanisms responsible for resistance to targeted therapies in BRAF-mutated melanoma and provides an overview of circulating biomarkers including circulating tumour cells, circulating tumour DNA, and non-coding RNAs.

Novel mechanisms of MITF regulation and melanoma predisposition identified in a mouse suppressor screen

MITF, a basic-Helix-Loop-Helix Zipper (bHLHZip) transcription factor, plays vital roles in melanocyte development and functions as an oncogene. To explore MITF regulation and its role in melanoma, we conducted a genetic screen for suppressors of the Mitf-associated pigmentation phenotype. An intragenic Mitf mutation was identified, leading to termination of MITF at the K316 SUMOylation site and loss of the C-end intrinsically disordered region (IDR). The resulting protein is more nuclear but less stable than wild-type MITF and retains DNA-binding ability. Interestingly, as a dimer, it can translocate wild-type and mutant MITF partners into the nucleus, improving its own stability and ensuring an active nuclear MITF supply. Interactions between K316 SUMOylation and S409 phosphorylation sites across monomers largely explain the observed effects. Notably, the recurrent melanoma-associated E318K mutation in MITF, which affects K316 SUMOylation, also alters protein regulation in concert with S409, unraveling a novel regulatory mechanism with unexpected disease insights.

Accurate prediction of pan-cancer types using machine learning with minimal number of DNA methylation sites

DNA methylation analysis has been applied to determine the primary site of cancer; however, robust and accurate prediction of cancer types with a minimum number of sites is still a significant scientific challenge. To build an accurate and robust cancer type prediction tool with a minimum number of DNA methylation sites, we internally benchmarked different DNA methylation site selection and ranking procedures, as well as different classification models. We used The Cancer Genome Atlas dataset (26 cancer types with 8296 samples) to train and test models and used an independent dataset (17 cancer types with 2738 samples) for model validation. A deep neural network model using a combined feature selection procedure (named MethyDeep) can predict 26 cancer types using 30 methylation sites with superior performance compared with the known methods for both primary and metastatic cancers in independent validation datasets. In conclusion, MethyDeep is an accurate and robust cancer type predictor with the minimum number of DNA methylation sites; it could help the cost-effective clarification of cancer of unknown primary patients and the liquid biopsy-based early screening of cancers.

Epigenetic suppression of PGC1α (PPARGC1A) causes collateral sensitivity to HMGCR-inhibitors within BRAF-treatment resistant melanomas

While targeted treatment against BRAF(V600E) improve survival for melanoma patients, many will see their cancer recur. Here we provide data indicating that epigenetic suppression of PGC1α defines an aggressive subset of chronic BRAF-inhibitor treated melanomas. A metabolism-centered pharmacological screen further identifies statins (HMGCR inhibitors) as a collateral vulnerability within PGC1α-suppressed BRAF-inhibitor resistant melanomas. Lower PGC1α levels mechanistically causes reduced RAB6B and RAB27A expression, whereby their combined re-expression reverses statin vulnerability. BRAF-inhibitor resistant cells with reduced PGC1α have increased integrin-FAK signaling and improved extracellular matrix detached survival cues that helps explain their increased metastatic ability. Statin treatment blocks cell growth by lowering RAB6B and RAB27A prenylation that reduces their membrane association and affects integrin localization and downstream signaling required for growth. These results suggest that chronic adaptation to BRAF-targeted treatments drive novel collateral metabolic vulnerabilities, and that HMGCR inhibitors may offer a strategy to treat melanomas recurring with suppressed PGC1α expression.

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.

In vitro cancer models as an approach to identify targetable developmental phenotypes in cancer stem cells

Cancer therapeutics are often highly toxic to the patient, and they often elicit rapid resistance in the tumour. Recent advances have suggested a potential new way in which we may improve on this, through two important concepts: (1) that multitudinous pathway alterations converge on a limited number of cancer cellular phenotypes, and (2) that these cancer cellular phenotypes depend on reactivation of developmental processes that are only minimally active in adult tissues. This provides a rationale for pursuing an approach of 'drugging the phenotype' focussed on targeting reactivated cellular processes from embryonic development. In this concepts paper, we cover these recent developments and their implications for the development of new cancer therapeutics that can avoid patient toxicity and acquired resistance. We then propose that in vitro tumour and developmental models can provide an experimental approach to identify and target the specific developmental processes at play, with a focus on the reactivation of developmental processes in the cancer stem cells that drive tumour progression and spread. Ultimately, the aim is to identify cellular processes that are specific to developmental phenotypes, are reactivated in cancer stem cells, and are essential to tumour progression. Therapeutically targeting these cellular processes could represent a new approach of 'drugging the phenotype' that treats the tumour whilst avoiding patient toxicity or the acquisition of therapeutic resistance.

Therapeutic targeting of anoikis resistance in cutaneous melanoma metastasis

The acquisition of resistance to anoikis, the cell death induced by loss of adhesion to the extracellular matrix, is an absolute requirement for the survival of disseminating and circulating tumour cells (CTCs), and for the seeding of metastatic lesions. In melanoma, a range of intracellular signalling cascades have been identified as potential drivers of anoikis resistance, however a full understanding of the process is yet to be attained. Mechanisms of anoikis resistance pose an attractive target for the therapeutic treatment of disseminating and circulating melanoma cells. This review explores the range of small molecule, peptide and antibody inhibitors targeting molecules involved in anoikis resistance in melanoma, and may be repurposed to prevent metastatic melanoma prior to its initiation, potentially improving the prognosis for patients.

Localized immune surveillance of primary melanoma in the skin deciphered through executable modeling

While skin is a site of active immune surveillance, primary melanomas often escape detection. Here, we have developed an in silico model to determine the local cross-talk between melanomas and Langerhans cells (LCs), the primary antigen-presenting cells at the site of melanoma development. The model predicts that melanomas fail to activate LC migration to lymph nodes until tumors reach a critical size, which is determined by a positive TNF-α feedback loop within melanomas, in line with our observations of murine tumors. In silico drug screening, supported by subsequent experimental testing, shows that treatment of primary tumors with MAPK pathway inhibitors may further prevent LC migration. In addition, our in silico model predicts treatment combinations that bypass LC dysfunction. In conclusion, our combined approach of in silico and in vivo studies suggests a molecular mechanism that explains how early melanomas develop under the radar of immune surveillance by LC.

Data-driven structural analysis of Small Cell Lung Cancer transcription factor network suggests potential subtype regulators and transition pathways

Small Cell Lung Cancer (SCLC) is an aggressive disease and challenging to treat due to its mixture of transcriptional subtypes and subtype transitions. Transcription factor (TF) networks have been the focus of studies to identify SCLC subtype regulators via systems approaches. Yet, their structures, which can provide clues on subtype drivers and transitions, are barely investigated. Here, we analyze the structure of an SCLC TF network by using graph theory concepts and identify its structurally important components responsible for complex signal processing, called hubs. We show that the hubs of the network are regulators of different SCLC subtypes by analyzing first the unbiased network structure and then integrating RNA-seq data as weights assigned to each interaction. Data-driven analysis emphasizes MYC as a hub, consistent with recent reports. Furthermore, we hypothesize that the pathways connecting functionally distinct hubs may control subtype transitions and test this hypothesis via network simulations on a candidate pathway and observe subtype transition. Overall, structural analyses of complex networks can identify their functionally important components and pathways driving the network dynamics. Such analyses can be an initial step for generating hypotheses and can guide the discovery of target pathways whose perturbation may change the network dynamics phenotypically.