Metastatic potential of melanomas defined by specific gene expression profiles with no BRAF signature

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.

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.

A developmental constraint model of cancer cell states and tumor heterogeneity

Cancer is a disease that stems from a fundamental liability inherent to multicellular life forms in which an individual cell is capable of reneging on the interests of the collective organism. Although cancer is commonly described as an evolutionary process, a less appreciated aspect of tumorigenesis may be the constraints imposed by the organism's developmental programs. Recent work from single-cell transcriptomic analyses across a range of cancer types has revealed the recurrence, plasticity, and co-option of distinct cellular states among cancer cell populations. Here, we note that across diverse cancer types, the observed cell states are proximate within the developmental hierarchy of the cell of origin. We thus posit a model by which cancer cell states are directly constrained by the organism's "developmental map." According to this model, a population of cancer cells traverses the developmental map, thereby generating a heterogeneous set of states whose interactions underpin emergent tumor behavior.

Calcifying Spindle Cell Soft Tissue Tumor With SOX10::PLAG1 Fusion: A Case Report of a Morphologically Distinctive and Potentially Novel Soft Tissue Tumor

The widespread use of advanced molecular techniques has led to the identification of several tumor types with PLAG1 gene fusions some of which also affect the skin and soft tissues. Herein, we present a 38-year-old female with a subcutaneous tumor affecting her forearm, which does not seem to fit into any currently recognized entity. It was a well-circumscribed tumor measuring 6 × 4,5 × 4 cm. It had a thick capsule composed of bland spindle cells forming palisades and Verocay body-like structures within a myxocollagenous background. Scattered calcifications were dispersed throughout the lesion. No cytological atypia, mitotic activity, or necrosis were present. Targeted NGS revealed a SOX10::PLAG1 fusion and fluorescent in situ hybridization confirmed the presence of PLAG1 gene rearrangement. The neoplastic cells showed a diffuse immunohistochemical expression of S100, SOX10, and PLAG1, as well as patchy desmin and CD34 positivity. The methylation profile of this tumor did not match any other entity covered by the DKFZ sarcoma classifier and apart from the gain of chromosome 12, the copy number profile was normal. The tumor was completely excised, and the patient has been free of disease for 4 years since the excision. While more cases are needed to confirm this tumor as a distinct entity, we propose a provisional name "SOX10::PLAG1-rearranged calcifying spindle cell tumor."

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.

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.

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.

Single-cell RNA sequencing in melanoma: what have we learned so far?

Over the past decade, there have been remarkable improvements in the treatment and survival rates of melanoma patients. Treatment resistance remains a persistent challenge, however, and is partly attributable to intratumoural heterogeneity. Melanoma cells can transition through a series of phenotypic and transcriptional cell states that vary in invasiveness and treatment responsiveness. The diverse stromal and immune contexture of the tumour microenvironment also contributes to intratumoural heterogeneity and disparities in treatment response in melanoma patients. Recent advances in single-cell sequencing technologies have enabled a more detailed understanding of melanoma heterogeneity and the underlying transcriptional programs that regulate melanoma cell diversity and behaviour. In this review, we examine the concept of intratumoural heterogeneity and the challenges it poses to achieving long-lasting treatment responses. We focus on the significance of next generation single-cell sequencing in advancing our understanding of melanoma diversity and the unique insights gained from single-cell studies.

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.

Targeting the RAS/RAF/MAPK pathway for cancer therapy: from mechanism to clinical studies

Metastatic dissemination of solid tumors, a leading cause of cancer-related mortality, underscores the urgent need for enhanced insights into the molecular and cellular mechanisms underlying metastasis, chemoresistance, and the mechanistic backgrounds of individuals whose cancers are prone to migration. The most prevalent signaling cascade governed by multi-kinase inhibitors is the mitogen-activated protein kinase (MAPK) pathway, encompassing the RAS-RAF-MAPK kinase (MEK)-extracellular signal-related kinase (ERK) pathway. RAF kinase is a primary mediator of the MAPK pathway, responsible for the sequential activation of downstream targets, such as MEK and the transcription factor ERK, which control numerous cellular and physiological processes, including organism development, cell cycle control, cell proliferation and differentiation, cell survival, and death. Defects in this signaling cascade are associated with diseases such as cancer. RAF inhibitors (RAFi) combined with MEK blockers represent an FDA-approved therapeutic strategy for numerous RAF-mutant cancers, including melanoma, non-small cell lung carcinoma, and thyroid cancer. However, the development of therapy resistance by cancer cells remains an important barrier. Autophagy, an intracellular lysosome-dependent catabolic recycling process, plays a critical role in the development of RAFi resistance in cancer. Thus, targeting RAF and autophagy could be novel treatment strategies for RAF-mutant cancers. In this review, we delve deeper into the mechanistic insights surrounding RAF kinase signaling in tumorigenesis and RAFi-resistance. Furthermore, we explore and discuss the ongoing development of next-generation RAF inhibitors with enhanced therapeutic profiles. Additionally, this review sheds light on the functional interplay between RAF-targeted therapies and autophagy in cancer.

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.

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.

LADON, a Natural Antisense Transcript of NODAL, Promotes Tumour Progression and Metastasis in Melanoma

The TGFβ family member NODAL, repeatedly required during embryonic development, has also been associated with tumour progression. Our aim was to clarify the controversy surrounding its involvement in melanoma tumour progression. We found that the deletion of the NODAL exon 2 in a metastatic melanoma cell line impairs its ability to form tumours and colonize distant tissues. However, we show that this phenotype does not result from the absence of NODAL, but from a defect in the expression of a natural antisense transcript of NODAL, here called LADON. We show that LADON expression is specifically activated in metastatic melanoma cell lines, that its transcript is packaged in exosomes secreted by melanoma cells, and that, via its differential impact on the expression of oncogenes and tumour suppressors, it promotes the mesenchymal to amoeboid transition that is critical for melanoma cell invasiveness. LADON is, therefore, a new player in the regulatory network governing tumour progression in melanoma and possibly in other types of cancer.

MITF regulates IDH1 and NNT and drives a transcriptional program protecting cutaneous melanoma from reactive oxygen species

Microphthalmia-associated transcription factor (MITF) plays pivotal roles in melanocyte development, function, and melanoma pathogenesis. MITF amplification occurs in melanoma and has been associated with resistance to targeted therapies. Here, we show that MITF regulates a global antioxidant program that increases survival of melanoma cell lines by protecting the cells from reactive oxygen species (ROS)-induced damage. In addition, this redox program is correlated with MITF expression in human melanoma cell lines and patient-derived melanoma samples. Using a zebrafish melanoma model, we show that MITF decreases ROS-mediated DNA damage in vivo . Some of the MITF target genes involved, such as IDH1 and NNT , are regulated through direct MITF binding to canonical enhancer box (E-BOX) sequences proximal to their promoters. Utilizing functional experiments, we demonstrate the role of MITF and its target genes in reducing cytosolic and mitochondrial ROS. Collectively, our data identify MITF as a significant driver of the cellular antioxidant state.

One Sentence Summary

MITF promote melanoma survival via increasing ROS tolerance.

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.

Transglutaminase Type 2-MITF axis regulates phenotype switching in skin cutaneous melanoma

Skin cutaneous melanoma (SKCM) is the deadliest form of skin cancer due to its high heterogeneity that drives tumor aggressiveness. Melanoma plasticity consists of two distinct phenotypic states that co-exist in the tumor niche, the proliferative and the invasive, respectively associated with a high and low expression of MITF, the master regulator of melanocyte lineage. However, despite efforts, melanoma research is still far from exhaustively dissecting this phenomenon. Here, we discovered a key function of Transglutaminase Type-2 (TG2) in regulating melanogenesis by modulating MITF transcription factor expression and its transcriptional activity. Importantly, we demonstrated that TG2 expression affects melanoma invasiveness, highlighting its positive value in SKCM. These results suggest that TG2 may have implications in the regulation of the phenotype switching by promoting melanoma differentiation and impairing its metastatic potential. Our findings offer potential perspectives to unravel melanoma vulnerabilities via tuning intra-tumor heterogeneity.

GABA Regulates Electrical Activity and Tumor Initiation in Melanoma

Oncogenes can initiate tumors only in certain cellular contexts, which is referred to as oncogenic competence. In melanoma, whether cells in the microenvironment can endow such competence remains unclear. Using a combination of zebrafish transgenesis coupled with human tissues, we demonstrate that GABAergic signaling between keratinocytes and melanocytes promotes melanoma initiation by BRAFV600E. GABA is synthesized in melanoma cells, which then acts on GABA-A receptors in keratinocytes. Electron microscopy demonstrates specialized cell-cell junctions between keratinocytes and melanoma cells, and multielectrode array analysis shows that GABA acts to inhibit electrical activity in melanoma/keratinocyte cocultures. Genetic and pharmacologic perturbation of GABA synthesis abrogates melanoma initiation in vivo. These data suggest that GABAergic signaling across the skin microenvironment regulates the ability of oncogenes to initiate melanoma.

SIGNIFICANCE

This study shows evidence of GABA-mediated regulation of electrical activity between melanoma cells and keratinocytes, providing a new mechanism by which the microenvironment promotes tumor initiation. This provides insights into the role of the skin microenvironment in early melanomas while identifying GABA as a potential therapeutic target in melanoma. See related commentary by Ceol, p. 2128. This article is featured in Selected Articles from This Issue, p. 2109.

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.

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

The Genomic Landscape of Melanoma and Its Therapeutic Implications

Melanoma is one of the most aggressive malignancies of the skin. The genetic composition of melanoma is complex and varies among different subtypes. With the aid of recent technologies such as next generation sequencing and single-cell sequencing, our understanding of the genomic landscape of melanoma and its tumor microenvironment has become increasingly clear. These advances may provide explanation to the heterogenic treatment outcomes of melanoma patients under current therapeutic guidelines and provide further insights to the development of potential new therapeutic targets. Here, we provide a comprehensive review on the genetics related to melanoma tumorigenesis, metastasis, and prognosis. We also review the genetics affecting the melanoma tumor microenvironment and its relation to tumor progression and treatment.

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.

Trametinib-Resistant Melanoma Cells Displaying MITFhigh/NGFRlow/IL-8low Phenotype Are Highly Responsive to Alternating Periods of Drug Withdrawal and Drug Rechallenge

Despite significant advances in targeted therapies against the hyperactivated BRAF^V600/MEK pathway for patients with unresectable metastatic melanoma, acquired resistance remains an unsolved clinical problem. In this study, we focused on melanoma cells resistant to trametinib, an agent broadly used in combination therapies. Molecular and cellular changes were assessed during alternating periods of trametinib withdrawal and rechallenge in trametinib-resistant cell lines displaying either a differentiation phenotype (MITF^high/NGFR^low) or neural crest stem-like dedifferentiation phenotype (NGFR^high/MITF^low). Neither drug withdrawal nor drug rechallenge induced cell death, and instead of loss of fitness, trametinib-resistant melanoma cells adapted to altered conditions by phenotype switching. In resistant cells displaying a differentiation phenotype, trametinib withdrawal markedly decreased MITF level and activity, which was associated with reduced cell proliferation capacity, and induced stemness assessed as NGFR-positive cells and senescence features, including IL-8 expression and secretion. All these changes could be reversed by trametinib re-exposure, which emphasizes melanoma cell plasticity. Trametinib-resistant cells displaying a dedifferentiation phenotype were less responsive presumably due to the already low level of MITF, a master regulator of the melanoma phenotype. Considering new directions of the development of anti-melanoma treatment, our study suggests that the phenotype of melanomas resistant to targeted therapy might be a crucial determinant of the selection of second-line therapy for melanoma patients.

The journey from melanocytes to melanoma

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

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

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

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.

Dissecting the roles and clinical potential of YY1 in the tumor microenvironment

Yin-Yang 1 (YY1) is a member of the GLI-Kruppel family of zinc finger proteins and plays a vital dual biological role in cancer as an oncogene or a tumor suppressor during tumorigenesis and tumor progression. The tumor microenvironment (TME) is identified as the “soil” of tumor that has a critical role in both tumor growth and metastasis. Many studies have found that YY1 is closely related to the remodeling and regulation of the TME. Herein, we reviewed the expression pattern of YY1 in tumors and summarized the function and mechanism of YY1 in regulating tumor angiogenesis, immune and metabolism. In addition, we discussed the potential value of YY1 in tumor diagnosis and treatment and provided a novel molecular strategy for the clinical diagnosis and treatment of tumors.

The molecular and functional landscape of resistance to immune checkpoint blockade in melanoma

Resistance to immune checkpoint inhibitor therapies in melanoma is common and remains an intractable clinical challenge. In this study, we comprehensively profile immune checkpoint inhibitor resistance mechanisms in short-term tumor cell lines and matched tumor samples from melanoma patients progressing on immune checkpoint inhibitors. Combining genome, transcriptome, and high dimensional flow cytometric profiling with functional analysis, we identify three distinct programs of immunotherapy resistance. Here we show that resistance programs include (1) the loss of wild-type antigen expression, resulting from tumor-intrinsic IFNγ signaling and melanoma de-differentiation, (2) the disruption of antigen presentation via multiple independent mechanisms affecting MHC expression, and (3) immune cell exclusion associated with PTEN loss. The dominant role of compromised antigen production and presentation in melanoma resistance to immune checkpoint inhibition highlights the importance of treatment salvage strategies aimed at the restoration of MHC expression, stimulation of innate immunity, and re-expression of wild-type differentiation antigens.

Ferroptosis in life: To be or not to be

Ferroptosis is a novel type of programmed cell death, characterized by a dysregulated iron metabolism and accumulation of lipid peroxides. It features the alteration of mitochondria and aberrant accumulation of excessive iron as well as loss of the cysteine-glutathione-GPX4 axis. Eventually, the accumulated lipid peroxides result in lethal damage to the cells. Ferroptosis is induced by the overloading of iron and the accumulation of ROS and can be inhibited by the activation of the GPX4 pathway, FS1-CoQ10 pathway, GCH1-BH4 pathway, and the DHODH pathway, it is also regulated by the oncogenes and tumor suppressors. Ferroptosis involves various physiological and pathological processes, and increasing evidence indicates that ferroptosis play a critical role in cancers and other diseases. It inhibits the proliferation of malignant cells in various types of cancers and inducing ferroptosis may become a new method of cancer treatment. Many inhibitors targeting the key factors of ferroptosis such as SLC7A11, GPX4, and iron overload have been developed. The application of ferroptosis is mainly divided into two directions, i.e. to avoid ferroptosis in healthy cells and selectively induce ferroptosis in cancers. In this review, we provide a critical analysis of the concept, and regulation pathways of ferroptosis and explored its roles in various diseases, we also summarized the compounds targeting ferroptosis, aiming to promote the speed of clinical use of ferroptosis induction in cancer treatment.

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

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

PPARs and the Kynurenine Pathway in Melanoma-Potential Biological Interactions

Peroxisome proliferator-activated receptors (PPARs) are ligand-activated transcription factors involved in various physiological and pathological processes within the skin. PPARs regulate several processes in one of the most aggressive skin cancers, melanoma, including proliferation, cell cycle, metabolic homeostasis, cell death, and metastasis. In this review, we focused not only on the biological activity of PPAR isoforms in melanoma initiation, progression, and metastasis but also on potential biological interactions between the PPAR signaling and the kynurenine pathways. The kynurenine pathway is a major pathway of tryptophan metabolism leading to nicotinamide adenine dinucleotide (NAD+) production. Importantly, various tryptophan metabolites exert biological activity toward cancer cells, including melanoma. Previous studies confirmed the functional relationship between PPAR and the kynurenine pathway in skeletal muscles. Despite the fact this interaction has not been reported in melanoma to date, some bioinformatics data and biological activity of PPAR ligands and tryptophan metabolites may suggest a potential involvement of these metabolic and signaling pathways in melanoma initiation, progression, and metastasis. Importantly, the possible relationship between the PPAR signaling pathway and the kynurenine pathway may relate not only to the direct biological effect on melanoma cells but also to the tumor microenvironment and the immune system.

S-adenosylmethionine blocks tumorigenesis and with immune checkpoint inhibitor enhances anti-cancer efficacy against BRAF mutant and wildtype melanomas

Despite marked success in treatment with immune checkpoint inhibitor (CPI), only a third of patients are responsive. Thus, melanoma still has one of the highest prevalence and mortality rates; which has led to a search for novel combination therapies that might complement CPI. Aberrant methylomes are one of the mechanisms of resistance to CPI therapy. S-adenosylmethionine (SAM), methyl donor of important epigenetic processes, has significant anti-cancer effects in several malignancies; however, SAM's effect has never been extensively investigated in melanoma. We demonstrate that SAM modulates phenotype switching of melanoma cells and directs the cells towards differentiation indicated by increased melanogenesis (melanin and melanosome synthesis), melanocyte-like morphology, elevated Mitf and Mitf activators' expression, increased antigen expression, reduced proliferation, and reduced stemness genes' expression. Consistently, providing SAM orally, reduced tumor growth and progression, and metastasis of syngeneic BRAF mutant and wild-type (WT) melanoma mouse models. Of note, SAM and anti-PD-1 antibody combination treatment had enhanced anti-cancer efficacy compared to monotherapies, showed significant reduction in tumor growth and progression, and increased survival. Furthermore, SAM and anti-PD-1 antibody combination triggered significantly higher immune cell infiltration, higher CD8+ T cells infiltration and effector functions, and polyfunctionality of CD8+ T cells in YUMMER1.7 tumors. Therefore, SAM combined with CPI provides a novel therapeutic strategy against BRAF mutant and WT melanomas and provides potential to be translated into clinic.

Mechanisms of resistance to immune checkpoint inhibitors in melanoma: What we have to overcome?

Marching into the second decade after the approval of ipilimumab, it is clear that immune checkpoint inhibitors (ICIs) have dramatically improved the prognosis of melanoma. Although the current edge is already high, with a 4-year OS% of 77.9% for adjuvant nivolumab and a 6.5-year OS% of 49% for nivolumab/ipilimumab combination in the metastatic setting, a high proportion of patients with advanced melanoma have no benefit from immunotherapy, or experience an early disease relapse/progression in the first few months of treatment, surviving much less. Reasonably, the primary and acquired resistance to ICIs has entered into the focus of clinical research with positive (e.g., nivolumab and relatlimab combination) and negative feedbacks (e.g., nivolumab with pegylated-IL2, pembrolizumab with T-VEC, nivolumab with epacadostat, and combinatorial triplets of BRAF/MEK inhibitors with immunotherapy). Many intrinsic (intracellular or intra-tumoral) but also extrinsic (systematic) events are considered to be involved in the development of this resistance to ICIs: i) melanoma cell immunogenicity (e.g., tumor mutational burden, antigen-processing machinery and immunogenic cell death, neoantigen affinity and heterogeneity, genomic instability, melanoma dedifferentiation and phenotypic plasticity), ii) immune cell trafficking, T-cell priming, and cell death evasion, iii) melanoma neovascularization, cellular TME components(e.g., T_regs, CAFs) and extracellular matrix modulation, iv) metabolic antagonism in the TME(highly glycolytic status, upregulated CD39/CD73/adenosine pathway, iDO-dependent tryptophan catabolism), v) T-cell exhaustion and negative immune checkpoints, and vi) gut microbiota. In the present overview, we discuss how these parameters compromise the efficacy of ICIs, with an emphasis on the lessons learned by the latest melanoma studies; and in parallel, we describe the main ongoing approaches to overcome the resistance to immunotherapy. Summarizing this information will improve the understanding of how these complicated dynamics contribute to immune escape and will help to develop more effective strategies on how anti-tumor immunity can surpass existing barriers of ICI-refractory melanoma.

The Regulators of Peroxisomal Acyl-Carnitine Shuttle CROT and CRAT Promote Metastasis in Melanoma

Circulating tumor cells are the key link between a primary tumor and distant metastases, but once in the bloodstream, loss of adhesion induces cell death. To identify the mechanisms relevant for melanoma circulating tumor cell survival, we performed RNA sequencing and discovered that detached melanoma cells and isolated melanoma circulating tumor cells rewire lipid metabolism by upregulating fatty acid (FA) transport and FA beta-oxidation‒related genes. In patients with melanoma, high expression of FA transporters and FA beta-oxidation enzymes significantly correlates with reduced progression-free and overall survival. Among the highest expressed regulators in melanoma circulating tumor cells were the carnitine transferases carnitine O-octanoyltransferase and carnitine acetyltransferase, which control the shuttle of peroxisome-derived medium-chain FAs toward mitochondria to fuel mitochondrial FA beta-oxidation. Knockdown of carnitine O-octanoyltransferase or carnitine acetyltransferase and short-term treatment with peroxisomal or mitochondrial FA beta-oxidation inhibitors thioridazine or ranolazine suppressed melanoma metastasis in mice. Carnitine O-octanoyltransferase and carnitine acetyltransferase depletion could be rescued by medium-chain FA supplementation, indicating that the peroxisomal supply of FAs is crucial for the survival of nonadherent melanoma cells. Our study identifies targeting the FA-based cross-talk between peroxisomes and mitochondria as a potential therapeutic opportunity to challenge melanoma progression. Moreover, the discovery of the antimetastatic activity of the Food and Drug Administration‒approved drug ranolazine carries translational potential.

A Comprehensive Analysis of Cutaneous Melanoma Patients in Greece Based on Multi-Omic Data

Cutaneous melanoma (CM) is the most aggressive type of skin cancer, and it is characterised by high mutational load and heterogeneity. In this study, we aimed to analyse the genomic and transcriptomic profile of primary melanomas from forty-six Formalin-Fixed, Paraffin-Embedded (FFPE) tissues from Greek patients. Molecular analysis for both germline and somatic variations was performed in genomic DNA from peripheral blood and melanoma samples, respectively, exploiting whole exome and targeted sequencing, and transcriptomic analysis. Detailed clinicopathological data were also included in our analyses and previously reported associations with specific mutations were recognised. Most analysed samples (43/46) were found to harbour at least one clinically actionable somatic variant. A subset of samples was profiled at the transcriptomic level, and it was shown that specific melanoma phenotypic states could be inferred from bulk RNA isolated from FFPE primary melanoma tissue. Integrative bioinformatics analyses, including variant prioritisation, differential gene expression analysis, and functional and gene set enrichment analysis by group and per sample, were conducted and molecular circuits that are implicated in melanoma cell programmes were highlighted. Integration of mutational and transcriptomic data in CM characterisation could shed light on genes and pathways that support the maintenance of phenotypic states encrypted into heterogeneous primary tumours.

New panel of biomarkers to discriminate between amelanotic and melanotic metastatic melanoma

Melanoma is a form of skin cancer that can rapidly invade distant organs. A distinctive feature of melanomas is their pigmentation status, as melanin is present in most skin melanomas, whilst many metastatic tumors could become amelanotic. Besides the obvious malfunction of the key genes of the melanin pathway, the amelanotic tumors could bear a characteristic molecular signature accounting for their aggressivity. Using mass spectrometry-based proteomics we report here a distinctive panel of biomarkers for amelanotic aggressive melanoma that differ from the less invasive pigmented cells. The developed method allows the label-free quantification of proteins identified by LC-MS/MS analysis. We found a set of proteins comprising AHNAK, MYOF, ANXA1, CAPN2, ASPH, EPHA2, THBS1, TGM2, ACTN4 along with proteins involved in cell adhesion/migration (integrins, PLEC, FSCN1, FN1) that are highly expressed in amelanotic melanoma. Accompanying the down regulation of pigmentation specific proteins such as tyrosinase and TYRP1, these biomarkers are highly specific for a type of highly invasive melanoma. Interestingly, the LC-MS/MS proteomics analysis in hypoxia revealed that the abundance of this specific set of proteins found in normoxia was rather unaltered in these conditions. These biomarkers could therefore predict a metastatic behaviour for the amelanotic cells in the early stages of the tumor development and thus serve in melanoma prognostic. Applying this algorithm to related databases including melanoma samples published by independent laboratories/public databases we confirm the specificity of the newly found signatures. Overall, we begin to unravel the molecular alterations in the amelanotic melanoma and how basic proteomics offers insights into how to assess the clinical, pathological and misdiagnosis differences between the main subtypes of melanoma.

Phenotype Switching and the Melanoma Microenvironment; Impact on Immunotherapy and Drug Resistance

Melanoma, a highly heterogeneous tumor, is comprised of a functionally diverse spectrum of cell phenotypes and subpopulations, including stromal cells in the tumor microenvironment (TME). Melanoma has been shown to dynamically shift between different transcriptional states or phenotypes. This is referred to as phenotype switching in melanoma, and it involves switching between quiescent and proliferative cell cycle states, and dramatic shifts in invasiveness, as well as changes in signaling pathways in the melanoma cells, and immune cell composition in the TME. Melanoma cell plasticity is associated with altered gene expression in immune cells and cancer-associated fibroblasts, as well as changes in extracellular matrix, which drive the metastatic cascade and therapeutic resistance. Therefore, resistance to therapy in melanoma is not only dependent on genetic evolution, but it has also been suggested to be driven by gene expression changes and adaptive phenotypic cell plasticity. This review discusses recent findings in melanoma phenotype switching, immunotherapy resistance, and the balancing of the homeostatic TME between the different melanoma cell subpopulations. We also discuss future perspectives of the biology of neural crest-like state(s) in melanoma.

Growth and adaptation mechanisms of tumour spheroids with time-dependent oxygen availability

Tumours are subject to external environmental variability. However, in vitro tumour spheroid experiments, used to understand cancer progression and develop cancer therapies, have been routinely performed for the past fifty years in constant external environments. Furthermore, spheroids are typically grown in ambient atmospheric oxygen (normoxia), whereas most in vivo tumours exist in hypoxic environments. Therefore, there are clear discrepancies between in vitro and in vivo conditions. We explore these discrepancies by combining tools from experimental biology, mathematical modelling, and statistical uncertainty quantification. Focusing on oxygen variability to develop our framework, we reveal key biological mechanisms governing tumour spheroid growth. Growing spheroids in time-dependent conditions, we identify and quantify novel biological adaptation mechanisms, including unexpected necrotic core removal, and transient reversal of the tumour spheroid growth phases.

A unique hyperdynamic dimer interface permits small molecule perturbation of the melanoma oncoprotein MITF for melanoma therapy

Microphthalmia transcription factor (MITF) regulates melanocyte development and is the "lineage-specific survival" oncogene of melanoma. MITF is essential for melanoma initiation, progression, and relapse and has been considered an important therapeutic target; however, direct inhibition of MITF through small molecules is considered impossible, due to the absence of a ligand-binding pocket for drug design. Here, our structural analyses show that the structure of MITF is hyperdynamic because of its out-of-register leucine zipper with a 3-residue insertion. The dynamic MITF is highly vulnerable to dimer-disrupting mutations, as we observed that MITF loss-of-function mutations in human Waardenburg syndrome type 2 A are frequently located on the dimer interface and disrupt the dimer forming ability accordingly. These observations suggest a unique opportunity to inhibit MITF with small molecules capable of disrupting the MITF dimer. From a high throughput screening against 654,650 compounds, we discovered compound TT-012, which specifically binds to dynamic MITF and destroys the latter's dimer formation and DNA-binding ability. Using chromatin immunoprecipitation assay and RNA sequencing, we showed that TT-012 inhibits the transcriptional activity of MITF in B16F10 melanoma cells. In addition, TT-012 inhibits the growth of high-MITF melanoma cells, and inhibits the tumor growth and metastasis with tolerable toxicity to liver and immune cells in animal models. Together, this study demonstrates a unique hyperdynamic dimer interface in melanoma oncoprotein MITF, and reveals a novel approach to therapeutically suppress MITF activity.

LAP1 supports nuclear adaptability during constrained melanoma cell migration and invasion

Metastasis involves dissemination of cancer cells away from a primary tumour and colonization at distal sites. During this process, the mechanical properties of the nucleus must be tuned since they pose a challenge to the negotiation of physical constraints imposed by the microenvironment and tissue structure. We discovered increased expression of the inner nuclear membrane protein LAP1 in metastatic melanoma cells, at the invasive front of human primary melanoma tumours and in metastases. Human cells express two LAP1 isoforms (LAP1B and LAP1C), which differ in their amino terminus. Here, using in vitro and in vivo models that recapitulate human melanoma progression, we found that expression of the shorter isoform, LAP1C, supports nuclear envelope blebbing, constrained migration and invasion by allowing a weaker coupling between the nuclear envelope and the nuclear lamina. We propose that LAP1 renders the nucleus highly adaptable and contributes to melanoma aggressiveness.

Frizzled receptors in melanomagenesis: From molecular interactions to target identification

Frizzled (FZD) proteins are receptors for the WNT family ligands. Inherited human diseases and genetic experiments using knockout mice have revealed a central role of FZDs in multiple aspects of embryonic development and tissue homeostasis. Misregulated FZD signaling has also been found in many cancers. Recent studies on three out of the ten mammalian FZDs in melanoma have shown that they promote tumor cell proliferation and invasion, via the activation of the canonical WNT/β-catenin or non-canonical PCP signaling pathway. In this concise review, we summarize our current knowledge of individual FZDs in melanoma, discuss the involvement of both the canonical and non-canonical pathways, and describe ongoing efforts to target the FZD receptors for melanoma treatment.

PI3K/AKT signaling allows for MAPK/ERK pathway independency mediating dedifferentiation-driven treatment resistance in melanoma

BACKGROUND

Current therapeutic management of advanced melanoma patients largely depends on their BRAF mutation status. However, the vast heterogeneity of the tumors hampers the success of therapies targeting the MAPK/ERK pathway alone. Dissecting this heterogeneity will contribute to identifying key players in the oncogenic progression to tailor more effective therapies.

METHODS

We performed a comprehensive molecular and phenotypic characterization of a panel of patient-derived BRAF^V600E-positive melanoma cell lines. Transcriptional profiling was used to identify groups of coregulated genes whose expression relates to an increased migratory potential and a higher resistance.

RESULTS

A decrease in sensitivity to MAPK/ERK pathway inhibition with vemurafenib or trametinib corresponded with an increasing quiescence and migratory properties of the cells. This was accompanied by the loss of transcriptional signatures of melanocytic differentiation, and the gain of stem cell features that conferred highly-resistant/mesenchymal-like cells with increased xenobiotic efflux capacity. Nevertheless, targeting of the implicated ABC transporters did not improve the response to vemurafenib, indicating that incomplete BRAF inhibition due to reduced drug uptake is not a main driver of resistance. Rather, indifference to MAPK/ERK pathway inhibition arose from the activation of compensatory signaling cascades. The PI3K/AKT pathway in particular showed a higher activity in mesenchymal-like cells, conferring a lower dependency on MAPK/ERK signaling and supporting stem-like properties that could be reverted by dual PI3K/mTOR inhibition with dactolisib.

CONCLUSIONS

In case of MAPK/ERK independency, therapeutic focus may be shifted to the PI3K/AKT pathway to overcome late-stage resistance in melanoma tumors that have acquired a mesenchymal phenotype. Video Abstract.

Targeting Key Signaling Pathways in Glioblastoma Stem Cells for the Development of Efficient Chemo- and Immunotherapy

Glioblastoma multiforme (GBM) is the most aggressive and most common malignant brain tumor with poor patient survival despite therapeutic intervention. On the cellular level, GBM comprises a rare population of glioblastoma stem cells (GSCs), driving therapeutic resistance, invasion, and recurrence. GSCs have thus come into the focus of therapeutic strategies, although their targeting remains challenging. In the present study, we took advantage of three GSCs-populations recently established in our lab to investigate key signaling pathways and subsequent therapeutic strategies targeting GSCs. We observed that NF-κB, a crucial transcription factor in GBM progression, was expressed in all CD44+/CD133+/Nestin+-GSC-populations. Exposure to TNFα led to activation of NF-κB-RELA and/or NF-κB-c-REL, depending on the GBM type. GSCs further expressed the proto-oncogene MYC family, with MYChigh GSCs being predominantly located in the tumor spheres ("GROW"-state) while NF-κB-RELAhigh GSCs were migrating out of the sphere ("GO"-state). We efficiently targeted GSCs by the pharmacologic inhibition of NF-κB using PTDC/Bortezomib or inhibition of MYC by KJ-Pyr-9, which significantly reduced GSC-viability, even in comparison to the standard chemotherapeutic drug temozolomide. As an additional cell-therapeutic strategy, we showed that NK cells could kill GSCs. Our findings offer new perspectives for developing efficient patient-specific chemo- and immunotherapy against GBM.

SRC-RAC1 signaling drives drug resistance to BRAF inhibition in de-differentiated cutaneous melanomas

Rare gain-of-function mutations in RAC1 drive drug resistance to targeted BRAF inhibition in cutaneous melanoma. Here, we show that wildtype RAC1 is a critical driver of growth and drug resistance, but only in a subset of melanomas with elevated markers of de-differentiation. Similarly, SRC inhibition also selectively sensitized de-differentiated melanomas to BRAF inhibition. One possible mechanism may be the suppression of the de-differentiated state, as SRC and RAC1 maintained markers of de-differentiation in human melanoma cells. The functional differences between melanoma subtypes suggest that the clinical management of cutaneous melanoma can be enhanced by the knowledge of differentiation status. To simplify the task of classification, we developed a binary classification strategy based on a small set of ten genes. Using this gene set, we reliably determined the differentiation status previously defined by hundreds of genes. Overall, our study informs strategies that enhance the precision of BRAFi by discovering unique vulnerabilities of the de-differentiated cutaneous melanoma subtype and creating a practical method to resolve differentiation status.

A nuclear cAMP microdomain suppresses tumor growth by Hippo pathway inactivation

Cyclic AMP (cAMP) signaling is localized to multiple spatially distinct microdomains, but the role of cAMP microdomains in cancer cell biology is poorly understood. Here, we present a tunable genetic system that allows us to activate cAMP signaling in specific microdomains. We uncover a nuclear cAMP microdomain that activates a tumor-suppressive pathway in a broad range of cancers by inhibiting YAP, a key effector protein of the Hippo pathway, inside the nucleus. We show that nuclear cAMP induces a LATS-dependent pathway leading to phosphorylation of nuclear YAP solely at serine 397 and export of YAP from the nucleus with no change in YAP protein stability. Thus, nuclear cAMP inhibition of nuclear YAP is distinct from other known mechanisms of Hippo regulation. Pharmacologic targeting of specific cAMP microdomains remains an untapped therapeutic approach for cancer; thus, drugs directed at the nuclear cAMP microdomain may provide avenues for the treatment of cancer.

Role of extracellular matrix architecture and signaling in melanoma therapeutic resistance

The extracellular matrix (ECM) is critical for maintaining tissue homeostasis therefore its production, assembly and mechanical stiffness are highly regulated in normal tissues. However, in solid tumors, increased stiffness resulting from abnormal ECM structural changes is associated with disease progression, an increased risk of metastasis and poor survival. As a dynamic and key component of the tumor microenvironment, the ECM is becoming increasingly recognized as an important feature of tumors, as it has been shown to promote several hallmarks of cancer via biochemical and biomechanical signaling. In this regard, melanoma cells are highly sensitive to ECM composition, stiffness and fiber alignment because they interact directly with the ECM in the tumor microenvironment via cell surface receptors, secreted factors or enzymes. Importantly, seeing as the ECM is predominantly deposited and remodeled by myofibroblastic stromal fibroblasts, it is a key avenue facilitating their paracrine interactions with melanoma cells. This review gives an overview of melanoma and further describes the critical roles that ECM properties such as ECM remodeling, ECM-related proteins and stiffness play in cutaneous melanoma progression, tumor cell plasticity and therapeutic resistance. Finally, given the emerging importance of ECM dynamics in melanoma, future perspectives on therapeutic strategies to normalize the ECM in tumors are discussed.

AP-1 transcription factor network explains diverse patterns of cellular plasticity in melanoma cells

Cellular plasticity associated with fluctuations in transcriptional programs allows individual cells in a tumor to adopt heterogeneous differentiation states and switch phenotype during their adaptive responses to therapies. Despite increasing knowledge of such transcriptional programs, the molecular basis of cellular plasticity remains poorly understood. Here, we combine multiplexed transcriptional and protein measurements at population and single-cell levels with multivariate statistical modeling to show that the state of AP-1 transcription factor network plays a unifying role in explaining diverse patterns of plasticity in melanoma. We find that a regulated balance among AP-1 factors cJUN, JUND, FRA2, FRA1, and cFOS determines the intrinsic diversity of differentiation states and adaptive responses to MAPK inhibitors in melanoma cells. Perturbing this balance through genetic depletion of specific AP-1 proteins, or by MAPK inhibitors, shifts cellular heterogeneity in a predictable fashion. Thus, AP-1 may serve as a critical node for manipulating cellular plasticity with potential therapeutic implications.