Systematic classification of melanoma cells by phenotype-specific gene expression mapping

Peroxidasin is associated with a mesenchymal-like transcriptional phenotype and promotes invasion in metastatic melanoma

Cutaneous melanoma is a highly invasive, heterogeneous and treatment resistant cancer. It's ability to dynamically shift between transcriptional states or phenotypes results in an adaptive cell plasticity that may drive cancer cell invasion or the development of therapy resistance. The expression of peroxidasin (PXDN), an extracellular matrix peroxidase, has been proposed to be associated with the invasive metastatic melanoma phenotype. We have confirmed this association by analysing the transcriptomes of 70 metastatic melanoma cell lines with variable levels of PXDN expression. This analysis highlighted a strong association between high PXDN expression and the undifferentiated invasive melanoma phenotype. To assess the functional role of PXDN in melanoma invasion, we performed a knockout of PXDN in a highly invasive cell line (NZM40). PXDN knockout decreased the invasive potential by ∼50 % and decreased the expression of epithelial-mesenchymal transition and invasive marker genes as determined by RNAseq and substantiated by proteomics analysis. Bioinformatics analysis of differentially expressed genes following PXDN knockout highlighted decreases in genes linked to extracellular matrix formation, organization and degradation as well as signalling pathways such as the WNT pathway. This study provides compelling evidence that PXDN plays a functional role in melanoma invasion by promoting an invasive, mesenchymal-like transcriptional phenotype.

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

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

Crosstalk in Skin: Loss of Desmoglein 1 in Keratinocytes Inhibits BRAFV600E-Induced Cellular Senescence in Human Melanocytes

Melanoma arises from transformation of melanocytes in the basal layer of epidermis where they are surrounded by keratinocytes, with which they interact through cell contact and paracrine communication. Although research focuses on how the accumulation of oncogene and tumor suppressor gene mutations in melanocytes drive melanomagenesis, how alterations in keratinocytes serve as extrinsic drivers of melanoma initiation and progression is poorly understood. We recently identified keratinocyte desmoglein 1 (DSG1) as an mediator of keratinocyte:melanoma crosstalk. In this study, we address the extent to which DSG1 loss, which occurs in response to environmental stress such as UVR, affects early steps in melanomagenesis. RNA-sequencing analysis revealed that paracrine signals from DSG1-deficient keratinocytes mediate a transcriptional switch from a differentiated to undifferentiated cell state in melanocytes expressing BRAFV600E. Of 221 differentially expressed genes in BRAFV600E cells treated with conditioned media from DSG1-deficient keratinocytes, the laminin superfamily member Netrin-4 (NTN4), which inhibits senescence, stood out. Indeed, although BRAFV600E melanocytes treated with conditioned media from DSG1-deficient keratinocytes showed signs of senescence bypass, NTN4 knockdown reversed these effects, whereas ectopic Netrin-4 expression mimicked them. These results suggest that DSG1 loss in keratinocytes provides an extrinsic signal to push melanocytes toward oncogenic transformation once an initial mutation has been introduced.

Interferon regulatory factor 4 modulates epigenetic silencing and cancer-critical pathways in melanoma cells

Interferon regulatory factor 4 (IRF4) was initially identified as a key controller in lymphocyte differentiation and function, and subsequently as a dependency factor and therapy target in lymphocyte-derived cancers. In melanocytes, IRF4 takes part in pigmentation. Although genetic studies have implicated IRF4 in melanoma, how IRF4 functions in melanoma cells has remained largely elusive. Here, we confirmed prevalent IRF4 expression in melanoma and showed that high expression is linked to dependency in cells and mortality in patients. Analysis of genes activated by IRF4 uncovered, as a novel target category, epigenetic silencing factors involved in DNA methylation (DNMT1, DNMT3B, UHRF1) and histone H3K27 methylation (EZH2). Consequently, we show that IRF4 controls the expression of tumour suppressor genes known to be silenced by these epigenetic modifications, for instance cyclin-dependent kinase inhibitors CDKN1A and CDKN1B, the PI3-AKT pathway regulator PTEN, and primary cilium components. Furthermore, IRF4 modulates activity of key downstream oncogenic pathways, such as WNT/β-catenin and AKT, impacting cell proliferation and survival. Accordingly, IRF4 modifies the effectiveness of pertinent epigenetic drugs on melanoma cells, a finding that encourages further studies towards therapeutic targeting of IRF4 in melanoma.

Melanoma progression and prognostic models drawn from single-cell, spatial maps of benign and malignant tumors

Melanoma clinical outcomes emerge from incompletely understood genetic mechanisms operating within the tumor and its microenvironment. Here, we used single-cell RNA-based spatial molecular imaging (RNA-SMI) in patient-derived archival tumors to reveal clinically relevant markers of malignancy progression and prognosis. We examined spatial gene expression of 203,472 cells inside benign and malignant melanocytic neoplasms, including melanocytic nevi and primary invasive and metastatic melanomas. Algorithmic cell clustering paired with intratumoral comparative two-dimensional analyses visualized synergistic, spatial gene signatures linking cellular proliferation, metabolism, and malignancy, validated by protein expression. Metastatic niches included up-regulation of CDK2 and FABP5, which independently predicted poor clinical outcome in 473 patients with melanoma via Cox regression analysis. More generally, our work demonstrates a framework for applying single-cell RNA-SMI technology toward identifying gene regulatory landscapes pertinent to cancer progression and patient survival.

The role of yes activated protein (YAP) in melanoma metastasis

Hippo was first identified in a genetic screen as a protein that suppressed proliferation and cell growth. Subsequently, it was shown that hippo acted in a so-called canonical cascade to suppress Yorkie, the Drosophila equivalent of Yes-activated protein (YAP), a mechanosensitive transcriptional cofactor that enhances the activity of the TEAD family of transcription factors. YAP promotes fibrosis, activation of cancer-associated fibroblasts, angiogenesis and cancer cell invasion. YAP activates the expression of the matricellular proteins CCN1 (cyr61) and CCN2 (ctgf), themselves mediators of fibrogenesis and oncogenesis, and coordination of matrix deposition and angiogenesis. This review discusses how therapeutically targeting YAP through YAP inhibitors verteporfin and celastrol and its downstream mediators CCN1 and CCN2 might be useful in treating melanoma.

Thiol starvation triggers melanoma state switching in an ATF4 and NRF2-dependent manner

The cystine/glutamate antiporter xCT is an important source of cysteine for cancer cells. Once taken up, cystine is reduced to cysteine and serves as a building block for the synthesis of glutathione, which efficiently protects cells from oxidative damage and prevents ferroptosis. As melanomas are particularly exposed to several sources of oxidative stress, we investigated the biological role of cysteine and glutathione supply by xCT in melanoma. xCT activity was abolished by genetic depletion in the Tyr::CreER; BrafCA; Ptenlox/+ melanoma model and by acute cystine withdrawal in melanoma cell lines. Both interventions profoundly impacted melanoma glutathione levels, but they were surprisingly well tolerated by murine melanomas in vivo and by most human melanoma cell lines in vitro. RNA sequencing of human melanoma cells revealed a strong adaptive upregulation of NRF2 and ATF4 pathways, which orchestrated the compensatory upregulation of genes involved in antioxidant defence and de novo cysteine biosynthesis. In addition, the joint activation of ATF4 and NRF2 triggered a phenotypic switch characterized by a reduction of differentiation genes and induction of pro-invasive features, which was also observed after erastin treatment or the inhibition of glutathione synthesis. NRF2 alone was capable of inducing the phenotypic switch in a transient manner. Together, our data show that cystine or glutathione levels regulate the phenotypic plasticity of melanoma cells by elevating ATF4 and NRF2.

The lipid droplet protein DHRS3 is a regulator of melanoma cell state

Lipid droplets are fat storage organelles composed of a protein envelope and lipid rich core. Regulation of this protein envelope underlies differential lipid droplet formation and function. In melanoma, lipid droplet formation has been linked to tumor progression and metastasis, but it is unknown whether lipid droplet proteins play a role. To address this, we performed proteomic analysis of the lipid droplet envelope in melanoma. We found that lipid droplet proteins were differentially enriched in distinct melanoma states; from melanocytic to undifferentiated. DHRS3, which converts all-trans-retinal to all-trans-retinol, is upregulated in the MITFLO/undifferentiated/neural crest-like melanoma cell state and reduced in the MITFHI/melanocytic state. Increased DHRS3 expression is sufficient to drive MITFHI/melanocytic cells to a more undifferentiated/invasive state. These changes are due to retinoic acid mediated regulation of melanocytic genes. Our data demonstrate that melanoma cell state can be regulated by expression of lipid droplet proteins which affect downstream retinoid signaling.

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.

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.

Tracking the EMT-like phenotype switching during targeted therapy in melanoma by analyzing extracellular vesicle phenotypes

Melanoma continues to be a leading cause of mortality among skin cancers. Despite advancements in targeted therapy, patients frequently develop resistance, leading to disease progression within a year. This resistance may result from the epithelial-to-mesenchymal transition (EMT)-like phenotype switching of melanoma cells. Tracking EMT-related phenotypic changes on extracellular vesicles (EVs) has potential to inform early about response to targeted therapy and melanoma progression. However, the knowledge on protein biomarkers carried by melanoma EVs involved in the EMT-like process remains unexplored. Herein, we present a biosensor integrating surface-enhanced Raman scattering and alternating current electrohydrodynamics-induced nanomixing enhancement, for sensitive detection of EMT-associated biomarkers on EV surfaces during targeted therapy. This biosensor successfully tracks the EMT-like phenotype switching in melanoma cell lines treated with mitogen-activated protein kinase inhibitor (MAPKi). Longitudinal monitoring of patients who receive MAPKi therapy and develop resistance, our biosensor shows its ability to identify the EMT-like phenotype switching on circulating EVs. This ability potentially can be leveraged to predict the development of resistance to targeted therapy, allowing for timely intervention and personalized treatment strategies.

Slide-tags enables single-nucleus barcoding for multimodal spatial genomics

Recent technological innovations have enabled the high-throughput quantification of gene expression and epigenetic regulation within individual cells, transforming our understanding of how complex tissues are constructed1-6. However, missing from these measurements is the ability to routinely and easily spatially localize these profiled cells. We developed a strategy, Slide-tags, in which single nuclei within an intact tissue section are tagged with spatial barcode oligonucleotides derived from DNA-barcoded beads with known positions. These tagged nuclei can then be used as an input into a wide variety of single-nucleus profiling assays. Application of Slide-tags to the mouse hippocampus positioned nuclei at less than 10 μm spatial resolution and delivered whole-transcriptome data that are indistinguishable in quality from ordinary single-nucleus RNA-sequencing data. To demonstrate that Slide-tags can be applied to a wide variety of human tissues, we performed the assay on brain, tonsil and melanoma. We revealed cell-type-specific spatially varying gene expression across cortical layers and spatially contextualized receptor-ligand interactions driving B cell maturation in lymphoid tissue. A major benefit of Slide-tags is that it is easily adaptable to almost any single-cell measurement technology. As a proof of principle, we performed multiomic measurements of open chromatin, RNA and T cell receptor (TCR) sequences in the same cells from metastatic melanoma, identifying transcription factor motifs driving cancer cell state transitions in spatially distinct microenvironments. Slide-tags offers a universal platform for importing the compendium of established single-cell measurements into the spatial genomics repertoire.

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.

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.

Ircinia ramosa Sponge Extract (iSP) Induces Apoptosis in Human Melanoma Cells and Inhibits Melanoma Cell Migration and Invasiveness

Marine compounds represent a varied source of new drugs with potential anticancer effects. Among these, sponges, including those belonging to the Irciniidae family, have been demonstrated to exert cytotoxic effects on different human cancer cells. Here, we investigated, for the first time, the therapeutic effect of an extract (referred as iSP) from the sponge, Ircinia ramosa (Porifera, Dictyoceratida, and Irciniidae), on A375 human melanoma cells. We found that iSP impaired A375 melanoma cells proliferation, induced cell death through caspase-dependent apoptosis and arrested cells in the G1 phase of the cell cycle, as demonstrated via both flow cytometry and qPCR analysis. The proapoptotic effect of iSP is associated with increased ROS production and mitochondrial modulation, as observed by using DCF-DHA and mitochondrial probes. In addition, we performed wound healing, invasion and clonogenic assays and found that iSP was able to restrain A375 migration, invasion and clonogenicity. Importantly, we observed that an iSP treatment modulated the expression of the EMT-associated epithelial markers, E-CAD and N-CAD, unveiling the mechanism underlying the effect of iSP in modulating A375 migration and invasion. Collectively, this study provides the first evidence to support the role of Ircinia ramosa sponge extracts as a potential therapeutic resource for the treatment of human melanoma.

TARGETING OF PDGF-C/NRP-1 AUTOCRINE LOOP AS A NEW STRATEGY FOR COUNTERACTING THE INVASIVENESS OF MELANOMA RESISTANT TO BRAF INHIBITORS

Melanoma resistance to BRAF inhibitors (BRAFi) is often accompanied by a switch from a proliferative to an invasive phenotype. Therefore, the identification of signaling molecules involved in the development of metastatic properties by resistant melanoma cells is of primary importance. We have previously demonstrated that activation of neuropilin-1 (NRP-1) by platelet-derived growth factor (PDGF)-C confers melanoma cells with an invasive behavior similar to that of BRAFi resistant tumors. Aims of the present study were to evaluate the role of PDGF-C/NRP-1 autocrine loop in the acquisition of an invasive and BRAFi-resistant phenotype by melanoma cells and the effect of its inhibition on drug resistance and extracellular matrix (ECM) invasion. Furthermore, we investigated whether PDGF-C serum levels were differentially modulated by drug treatment in metastatic melanoma patients responsive or refractory to BRAFi as a single agent or in combination with MEK inhibitors (MEKi). The results indicated that human melanoma cells resistant to BRAFi express higher levels of PDGF-C and NRP-1 as compared to their susceptible counterparts. Overexpression occurs early during development of drug resistance and contributes to the invasive properties of resistant cells. Accordingly, silencing of NRP-1 or PDGF-C reduces tumor cell invasiveness. Analysis of PDGF-C in the serum collected from patients treated with BRAFi or BRAFi+MEKi, showed that in responders PDGF-C levels decrease after treatment and raise again at tumor progression. Conversely, in non-responders treatment does not affect PDGF-C serum levels. Thus, blockade of NRP-1 activation by PDGF-C might represent a new therapeutic approach to counteract the invasiveness of BRAFi-resistant melanoma.

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.

Slide-tags: scalable, single-nucleus barcoding for multi-modal spatial genomics

Recent technological innovations have enabled the high-throughput quantification of gene expression and epigenetic regulation within individual cells, transforming our understanding of how complex tissues are constructed. Missing from these measurements, however, is the ability to routinely and easily spatially localise these profiled cells. We developed a strategy, Slide-tags, in which single nuclei within an intact tissue section are 'tagged' with spatial barcode oligonucleotides derived from DNA-barcoded beads with known positions. These tagged nuclei can then be used as input into a wide variety of single-nucleus profiling assays. Application of Slide-tags to the mouse hippocampus positioned nuclei at less than 10 micron spatial resolution, and delivered whole-transcriptome data that was indistinguishable in quality from ordinary snRNA-seq. To demonstrate that Slide-tags can be applied to a wide variety of human tissues, we performed the assay on brain, tonsil, and melanoma. We revealed cell-type-specific spatially varying gene expression across cortical layers and spatially contextualised receptor-ligand interactions driving B-cell maturation in lymphoid tissue. A major benefit of Slide-tags is that it is easily adaptable to virtually any single-cell measurement technology. As proof of principle, we performed multiomic measurements of open chromatin, RNA, and T-cell receptor sequences in the same cells from metastatic melanoma. We identified spatially distinct tumour subpopulations to be differentially infiltrated by an expanded T-cell clone and undergoing cell state transition driven by spatially clustered accessible transcription factor motifs. Slide-tags offers a universal platform for importing the compendium of established single-cell measurements into the spatial genomics repertoire.

Crosstalk in skin: Loss of desmoglein 1 in keratinocytes inhibits BRAFV600E-induced cellular senescence in human melanocytes

Melanoma arises from transformation of melanocytes in the basal layer of the epidermis where they are surrounded by keratinocytes, with which they interact through cell contact and paracrine communication. Considerable effort has been devoted to determining how the accumulation of oncogene and tumor suppressor gene mutations in melanocytes drive melanoma development. However, the extent to which alterations in keratinocytes that occur in the developing tumor niche serve as extrinsic drivers of melanoma initiation and progression is poorly understood. We recently identified the keratinocyte-specific cadherin, desmoglein 1 (Dsg1), as an important mediator of keratinocyte:melanoma cell crosstalk, demonstrating that its chronic loss, which can occur through melanoma cell-dependent paracrine signaling, promotes behaviors that mimic a malignant phenotype. Here we address the extent to which Dsg1 loss affects early steps in melanomagenesis. RNA-Seq analysis revealed that paracrine signals from Dsg1-deficient keratinocytes mediate a transcriptional switch from a differentiated to undifferentiated cell state in melanocytes expressing BRAFV600E, a driver mutation commonly present in both melanoma and benign nevi and reported to cause growth arrest and oncogene-induced senescence (OIS). Of ~220 differentially expressed genes in BRAFV600E cells treated with Dsg1-deficient conditioned media (CM), the laminin superfamily member NTN4/Netrin-4, which inhibits senescence in endothelial cells, stood out. Indeed, while BRAFV600E melanocytes treated with Dsg1-deficient CM showed signs of senescence bypass as assessed by increased senescence-associated β-galactosidase activity and decreased p16, knockdown of NTN4 reversed these effects. These results suggest that Dsg1 loss in keratinocytes provides an extrinsic signal to push melanocytes towards oncogenic transformation once an initial mutation has been introduced.

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.

Antiproliferative and Proapoptotic Effects of Erucin, a Diet-Derived H2S Donor, on Human Melanoma Cells

Melanoma is the most dangerous form of skin cancer and is characterized by chemotherapy resistance and recurrence despite the new promising therapeutic approaches. In the last years, erucin (ERU), the major isothiocyanate present in Eruca sativa, commonly known as rocket salads, has demonstrated great efficacy as an anticancer agent in different in vitro and in vivo models. More recently, the chemopreventive effects of ERU have been associated with its property of being a H₂S donor in human pancreatic adenocarcinoma. Here, we investigated the effects of ERU in modulating proliferation and inducing human melanoma cell death by using multiple in vitro approaches. ERU significantly reduced the proliferation of different human melanoma cell lines. A flow cytometry analysis with annexin V/PI demonstrated that ERU was able to induce apoptosis and cell cycle arrest in A375 melanoma cells. The proapoptotic effect of ERU was associated with the modulation of the epithelial-to-mesenchymal transition (EMT)-related cadherins and transcription factors. Moreover, ERU thwarted the migration, invasiveness and clonogenic abilities of A375 melanoma cells. These effects were associated with melanogenesis impairment and mitochondrial fitness modulation. Therefore, we demonstrated that ERU plays an important role in inhibiting the progression of melanoma and could represent a novel add-on therapy for the treatment of human melanoma.

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 cellular hierarchy in melanoma uncouples growth and metastasis

Although melanoma is notorious for its high degree of heterogeneity and plasticity1,2, the origin and magnitude of cell-state diversity remains poorly understood. Equally, it is unclear whether growth and metastatic dissemination are supported by overlapping or distinct melanoma subpopulations. Here, by combining mouse genetics, single-cell and spatial transcriptomics, lineage tracing and quantitative modelling, we provide evidence of a hierarchical model of tumour growth that mirrors the cellular and molecular logic underlying the cell-fate specification and differentiation of the embryonic neural crest. We show that tumorigenic competence is associated with a spatially localized perivascular niche, a phenotype acquired through an intercellular communication pathway established by endothelial cells. Consistent with a model in which only a fraction of cells are fated to fuel growth, temporal single-cell tracing of a population of melanoma cells with a mesenchymal-like state revealed that these cells do not contribute to primary tumour growth but, instead, constitute a pool of metastatic initiating cells that switch cell identity while disseminating to secondary organs. Our data provide a spatially and temporally resolved map of the diversity and trajectories of melanoma cell states and suggest that the ability to support growth and metastasis are limited to distinct pools of cells. The observation that these phenotypic competencies can be dynamically acquired after exposure to specific niche signals warrant the development of therapeutic strategies that interfere with the cancer cell reprogramming activity of such microenvironmental cues.

Innate immune checkpoint inhibitor resistance is associated with melanoma sub-types exhibiting invasive and de-differentiated gene expression signatures

Melanoma is a highly aggressive skin cancer, which, although highly immunogenic, frequently escapes the body’s immune defences. Immune checkpoint inhibitors (ICI), such as anti-PD1, anti-PDL1, and anti-CTLA4 antibodies lead to reactivation of immune pathways, promoting rejection of melanoma. However, the benefits of ICI therapy remain limited to a relatively small proportion of patients who do not exhibit ICI resistance. Moreover, the precise mechanisms underlying innate and acquired ICI resistance remain unclear. Here, we have investigated differences in melanoma tissues in responder and non-responder patients to anti-PD1 therapy in terms of tumour and immune cell gene-associated signatures. We performed multi-omics investigations on melanoma tumour tissues, which were collected from patients before starting treatment with anti-PD1 immune checkpoint inhibitors. Patients were subsequently categorized into responders and non-responders to anti-PD1 therapy based on RECIST criteria. Multi-omics analyses included RNA-Seq and NanoString analysis. From RNA-Seq data we carried out HLA phenotyping as well as gene enrichment analysis, pathway enrichment analysis and immune cell deconvolution studies. Consistent with previous studies, our data showed that responders to anti-PD1 therapy had higher immune scores (median immune score for responders = 0.1335, median immune score for non-responders = 0.05426, p-value = 0.01, Mann-Whitney U two-tailed exact test) compared to the non-responders. Responder melanomas were more highly enriched with a combination of CD8+ T cells, dendritic cells (p-value = 0.03) and an M1 subtype of macrophages (p-value = 0.001). In addition, melanomas from responder patients exhibited a more differentiated gene expression pattern, with high proliferative- and low invasive-associated gene expression signatures, whereas tumours from non-responders exhibited high invasive- and frequently neural crest-like cell type gene expression signatures. Our findings suggest that non-responder melanomas to anti-PD1 therapy exhibit a de-differentiated gene expression signature, associated with poorer immune cell infiltration, which establishes a gene expression pattern characteristic of innate resistance to anti-PD1 therapy. Improved understanding of tumour-intrinsic gene expression patterns associated with response to anti-PD1 therapy will help to identify predictive biomarkers of ICI response and may help to identify new targets for anticancer treatment, especially with a capacity to function as adjuvants to improve ICI outcomes.

Bioinformatic analysis identifies epidermal development genes that contribute to melanoma progression

Several diagnostic and prognostic markers for melanoma have been identified in last few years. However, their actual contribution to melanoma progression have not been investigated in detail. This study was aimed to identify genes, biological processes, and signaling pathways implicated in melanoma progression by applying bioinformatics analysis. We identified nine differentially expressed genes (DEGs) (IL36RN, KRT6A, KRT6B, KRT16, S100A7, SPRR1A, SPRR1B, SPRR2B, and KLK7) that were upregulated in primary melanoma compared with metastatic melanoma in all five datasets analyzed. All these genes except IL36RN, both form a protein-protein interaction network and have cellular functions associated with constitutive processes of keratinocytes. Thus, they were generically termed Epidermal Development and Cornification (EDC) genes. The differential expression of these genes in primary and metastatic melanoma was confirmed in the TCGA-SKCM cohort. High expression of the EDC genes correlated with reduced tumor thickness in primary melanoma and shorter survival in metastatic melanoma. Analysis of DEGs from primary melanoma patients displaying high or low expression of all eight EDC revealed that the upregulated genes are enriched in biological process related to cell migration, extracellular matrix organization, invasion, and Epithelial-Mesenchymal Transition. Further analysis of enriched curated oncogenic genesets together with RPPA data of phosphorylated proteins revealed the activation of MEK, ATF2, and EGFR pathways in tumors displaying high expression of EDC genes. Thus, EDC genes may contribute to melanoma progression by promoting the activation of MEK, ATF2, and EGFR pathways together with biological processes associated with tumor aggressiveness.

Multi-modal molecular programs regulate melanoma cell state

Melanoma cells display distinct intrinsic phenotypic states. Here, we seek to characterize the molecular regulation of these states using multi-omic analyses of whole exome, transcriptome, microRNA, long non-coding RNA and DNA methylation data together with reverse-phase protein array data on a panel of 68 highly annotated early passage melanoma cell lines. We demonstrate that clearly defined cancer cell intrinsic transcriptomic programs are maintained in melanoma cells ex vivo and remain highly conserved within melanoma tumors, are associated with distinct immune features within tumors, and differentially correlate with checkpoint inhibitor and adoptive T cell therapy efficacy. Through integrative analyses we demonstrate highly complex multi-omic regulation of melanoma cell intrinsic programs that provide key insights into the molecular maintenance of phenotypic states. These findings have implications for cancer biology and the identification of new therapeutic strategies. Further, these deeply characterized cell lines will serve as an invaluable resource for future research in the field.

Scales of Cancer Evolution: Selfish Genome or Cooperating Cells?

Simple Summary

Cancer continuously evolves its ability to survive in time-varying microenvironment, which results, regarding the therapeutic context, in its therapeutic resistance. As it is accepted that the development of resistance is the direct consequence of intratumour heterogeneity, its evolutionary etiology is intensively studied. Models of carinogenesis are often assessed accordingly to how well they fit into the evolutionary scenario. In the paper, the relevant observations and concepts in cancer research, such as intratumour heterogeneity, cell plasticity, and Markov cell state dynamics, are reviewed and integrated into an evolutionary model. The possibility that the interaction between cancer cells can be interpreted as cooperation is proposed.

Abstract

The exploitation of the evolutionary modus operandi of cancer to steer its progression towards drug sensitive cancer cells is a challenging research topic. Integrating evolutionary principles into cancer therapy requires properly identified selection level, the relevant timescale, and the respective fitness of the principal selection unit on that timescale. Interpretation of some features of cancer progression, such as increased heterogeneity of isogenic cancer cells, is difficult from the most straightforward evolutionary view with the cancer cell as the principal selection unit. In the paper, the relation between the two levels of intratumour heterogeneity, genetic, due to genetic instability, and non-genetic, due to phenotypic plasticity, is reviewed and the evolutionary role of the latter is outlined. In analogy to the evolutionary optimization in a changing environment, the cell state dynamics in cancer clones are interpreted as the risk diversifying strategy bet hedging, optimizing the balance between the exploitation and exploration of the cell state space.

Heterogeneity in Melanoma

There is growing evidence that tumour heterogeneity has an imperative role in cancer development, evolution and resistance to therapy. Continuing advancements in biomedical research enable tumour heterogeneity to be observed and studied more critically. As one of the most heterogeneous human cancers, melanoma displays a high level of biological complexity during disease progression. However, much is still unknown regarding melanoma tumour heterogeneity, as well as the role it plays in disease progression and treatment response. This review aims to provide a concise summary of the importance of tumour heterogeneity in melanoma.

Loss of YY1, a Regulator of Metabolism in Melanoma, Drives Melanoma Cell Invasiveness and Metastasis Formation

Deregulation of cellular metabolism through metabolic rewiring and translational reprogramming are considered hallmark traits of tumor development and malignant progression. The transcription factor YY1 is a master regulator of metabolism that we have previously shown to orchestrate a metabolic program required for melanoma formation. In this study, we demonstrate that YY1, while being essential for primary melanoma formation, suppresses metastatic spreading. Its downregulation or loss resulted in the induction of an invasiveness gene program and sensitized melanoma cells for pro-invasive signaling molecules, such as TGF-β. In addition, NGFR, a key effector in melanoma invasion and phenotype switching, was among the most upregulated genes after YY1 knockdown. High levels of NGFR were also associated with other metabolic stress inducers, further indicating that YY1 knockdown mimics a metabolic stress program associated with an increased invasion potential in melanoma. Accordingly, while counteracting tumor growth, loss of YY1 strongly promoted melanoma cell invasiveness in vitro and metastasis formation in melanoma mouse models in vivo. Thus, our findings show that the metabolic regulator YY1 controls phenotype switching in melanoma.

Altered BAF occupancy and transcription factor dynamics in PBAF-deficient melanoma

ARID2 is the most recurrently mutated SWI/SNF complex member in melanoma; however, its tumor-suppressive mechanisms in the context of the chromatin landscape remain to be elucidated. Here, we model ARID2 deficiency in melanoma cells, which results in defective PBAF complex assembly with a concomitant genomic redistribution of the BAF complex. Upon ARID2 depletion, a subset of PBAF and shared BAF-PBAF-occupied regions displays diminished chromatin accessibility and associated gene expression, while BAF-occupied enhancers gain chromatin accessibility and expression of genes linked to the process of invasion. As a function of altered accessibility, the genomic occupancy of melanoma-relevant transcription factors is affected and significantly correlates with the observed transcriptional changes. We further demonstrate that ARID2-deficient cells acquire the ability to colonize distal organs in multiple animal models. Taken together, our results reveal a role for ARID2 in mediating BAF and PBAF subcomplex chromatin dynamics with consequences for melanoma metastasis.

Intermittent treatment of BRAFV600E melanoma cells delays resistance by adaptive resensitization to drug rechallenge

Preclinical studies of metastatic melanoma treated with targeted therapeutics have suggested that alternating periods of treatment and withdrawal might delay the onset of resistance. This has been attributed to drug addiction, where cells lose fitness upon drug removal due to the resulting hyperactivation of mitogen-activated protein (MAP) kinase signaling. This study presents evidence that the intermittent treatment response can also be explained by the resensitization of cells following drug removal and enhanced cell loss upon drug rechallenge. Resensitization is accompanied by adaptive transcriptomic switching and occurs despite the sustained expression of resistance genes throughout the intermittent treatment.

Patients with melanoma receiving drugs targeting BRAF^V600E and mitogen-activated protein (MAP) kinase kinases 1 and 2 (MEK1/2) invariably develop resistance and face continued progression. Based on preclinical studies, intermittent treatment involving alternating periods of drug withdrawal and rechallenge has been proposed as a method to delay the onset of resistance. The beneficial effect of intermittent treatment has been attributed to drug addiction, where drug withdrawal reduces the viability of resistant cells due to MAP kinase pathway hyperactivation. However, the mechanistic basis of the intermittent effect is incompletely understood. We show that intermittent treatment with the BRAF^V600E inhibitor, LGX818/encorafenib, suppresses growth compared with continuous treatment in human melanoma cells engineered to express BRAF^V600E, p61-BRAF^V600E, or MEK2^C125 oncogenes. Analysis of the BRAF^V600E-overexpressing cells shows that, while drug addiction clearly occurs, it fails to account for the advantageous effect of intermittent treatment. Instead, growth suppression is best explained by resensitization during periods of drug removal, followed by cell death after drug readdition. Continuous treatment leads to transcriptional responses prominently associated with chemoresistance in melanoma. By contrast, cells treated intermittently reveal a subset of transcripts that reverse expression between successive cycles of drug removal and rechallenge and include mediators of cell invasiveness and the epithelial-to-mesenchymal transition. These transcripts change during periods of drug removal by adaptive switching, rather than selection pressure. Resensitization occurs against a background of sustained expression of melanoma resistance genes, producing a transcriptome distinct from that of the initial drug-naive cell state. We conclude that phenotypic plasticity leading to drug resensitization can underlie the beneficial effect of intermittent treatment.

A DNA replication-independent function of pre-replication complex genes during cell invasion in C. elegans

Cell invasion is an initiating event during tumor cell metastasis and an essential process during development. A screen of C. elegans orthologs of genes overexpressed in invasive human melanoma cells has identified several components of the conserved DNA pre-replication complex (pre-RC) as positive regulators of anchor cell (AC) invasion. The pre-RC genes function cell-autonomously in the G1-arrested AC to promote invasion, independently of their role in licensing DNA replication origins in proliferating cells. While the helicase activity of the pre-RC is necessary for AC invasion, the downstream acting DNA replication initiation factors are not required. The pre-RC promotes the invasive fate by regulating the expression of extracellular matrix genes and components of the PI3K signaling pathway. Increasing PI3K pathway activity partially suppressed the AC invasion defects caused by pre-RC depletion, suggesting that the PI3K pathway is one critical pre-RC target. We propose that the pre-RC, or a part of it, acts in the postmitotic AC as a transcriptional regulator that facilitates the switch to an invasive phenotype.

YAP facilitates melanoma migration through regulation of actin-related protein 2/3 complex subunit 5 (ARPC5)

Yes-associated protein 1 (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ) are transcriptional coactivators that have been implicated in driving metastasis and progression in many cancers, mainly through their transcriptional regulation of downstream targets. Although YAP and TAZ have shown redundancy in many contexts, it is still unknown whether or not this is true in melanoma. Here, we show that while both YAP and TAZ are expressed in a panel of melanoma cell lines, depletion of YAP results in decreased cell numbers, focal adhesions, and the ability to invade matrigel. Using non-biased RNA-sequencing analysis, we find that melanoma cells depleted of YAP, TAZ, or YAP/TAZ exhibit drastically different transcriptomes. We further uncover the ARP2/3 subunit ARPC5 as a specific target of YAP but not TAZ and that ARPC5 is essential for YAP-dependent maintenance of melanoma cell focal adhesion numbers. Our findings suggest that in melanoma, YAP drives melanoma progression, survival, and invasion.

Six potential biomarkers for bladder cancer: key proteins in cell-cycle division and apoptosis pathways

BACKGROUND

The bladder cancer (BC) pathology is caused by both exogenous environmental and endogenous molecular factors. Several genes have been implicated, but the molecular pathogenesis of BC and its subtypes remains debatable. The bioinformatic analysis evaluates high numbers of proteins in a single study, increasing the opportunity to identify possible biomarkers for disorders.

METHODS

The aim of this study is to identify biomarkers for the identification of BC using several bioinformatic analytical tools and methods. BC and normal samples were compared for each probeset with T test in GSE13507 and GSE37817 datasets, and statistical probesets were verified with GSE52519 and E-MTAB-1940 datasets. Differential gene expression, hierarchical clustering, gene ontology enrichment analysis, and heuristic online phenotype prediction algorithm methods were utilized. Statistically significant proteins were assessed in the Human Protein Atlas database. GSE13507 (6271 probesets) and GSE37817 (3267 probesets) data were significant after the extraction of probesets without gene annotation information. Common probesets in both datasets (2888) were further narrowed by analyzing the first 100 upregulated and downregulated probesets in BC samples.

RESULTS

Among the total 400 probesets, 68 were significant for both datasets with similar fold-change values (Pearson r: 0.995). Protein-protein interaction networks demonstrated strong interactions between CCNB1, BUB1B, and AURKB. The HPA database revealed similar protein expression levels for CKAP2L, AURKB, APIP, and LGALS3 both for BC and control samples.

CONCLUSION

This study disclosed six candidate biomarkers for the early diagnosis of BC. It is suggested that these candidate proteins be investigated in a wet lab to identify their functions in BC pathology and possible treatment approaches.

Targeting Discoidin Domain Receptors DDR1 and DDR2 overcomes matrix-mediated tumor cell adaptation and tolerance to BRAF-targeted therapy in melanoma

Resistance to BRAF/MEK inhibitor therapy in BRAF^V600‐mutated advanced melanoma remains a major obstacle that limits patient benefit. Microenvironment components including the extracellular matrix (ECM) can support tumor cell adaptation and tolerance to targeted therapy; however, the underlying mechanisms remain poorly understood. Here, we investigated the process of matrix‐mediated drug resistance (MMDR) in response to BRAF^V600 pathway inhibition in melanoma. We demonstrate that physical and structural cues from fibroblast‐derived ECM abrogate anti‐proliferative responses to BRAF/MEK inhibition. MMDR is mediated by drug‐induced linear clustering of phosphorylated DDR1 and DDR2, two tyrosine kinase collagen receptors. Depletion and pharmacological targeting of DDR1 and DDR2 overcome ECM‐mediated resistance to BRAF‐targeted therapy. In xenografts, targeting DDR with imatinib enhances BRAF inhibitor efficacy, counteracts drug‐induced collagen remodeling, and delays tumor relapse. Mechanistically, DDR‐dependent MMDR fosters a targetable pro‐survival NIK/IKKα/NF‐κB2 pathway. These findings reveal a novel role for a collagen‐rich matrix and DDR in tumor cell adaptation and resistance. They also provide important insights into environment‐mediated drug resistance and a preclinical rationale for targeting DDR signaling in combination with targeted therapy in melanoma.

Hierarchy of TGFβ/SMAD, Hippo/YAP/TAZ, and Wnt/β-catenin signaling in melanoma phenotype switching

TGFβ, YAP/TAZ, and canonical Wnt/β-catenin signaling functionally interact in a hierarchical manner to induce the switching of melanoma cells from proliferative-to-invasive cell phenotype.

In melanoma, a switch from a proliferative melanocytic to an invasive mesenchymal phenotype is based on dramatic transcriptional reprogramming which involves complex interactions between a variety of signaling pathways and their downstream transcriptional regulators. TGFβ/SMAD, Hippo/YAP/TAZ, and Wnt/β-catenin signaling pathways are major inducers of transcriptional reprogramming and converge at several levels. Here, we report that TGFβ/SMAD, YAP/TAZ, and β-catenin are all required for a proliferative-to-invasive phenotype switch. Loss and gain of function experimentation, global gene expression analysis, and computational nested effects models revealed the hierarchy between these signaling pathways and identified shared target genes. SMAD-mediated transcription at the top of the hierarchy leads to the activation of YAP/TAZ and of β-catenin, with YAP/TAZ governing an essential subprogram of TGFβ-induced phenotype switching. Wnt/β-catenin signaling is situated further downstream and exerts a dual role: it promotes the proliferative, differentiated melanoma cell phenotype and it is essential but not sufficient for SMAD or YAP/TAZ–induced phenotype switching. The results identify epistatic interactions among the signaling pathways underlying melanoma phenotype switching and highlight the priorities in targets for melanoma therapy.

Cooperation between melanoma cell states promotes metastasis through heterotypic cluster formation

Melanomas can have multiple coexisting cell states, including proliferative (PRO) versus invasive (INV) subpopulations that represent a "go or grow" trade-off; however, how these populations interact is poorly understood. Using a combination of zebrafish modeling and analysis of patient samples, we show that INV and PRO cells form spatially structured heterotypic clusters and cooperate in the seeding of metastasis, maintaining cell state heterogeneity. INV cells adhere tightly to each other and form clusters with a rim of PRO cells. Intravital imaging demonstrated cooperation in which INV cells facilitate dissemination of less metastatic PRO cells. We identified the TFAP2 neural crest transcription factor as a master regulator of clustering and PRO/INV states. Isolation of clusters from patients with metastatic melanoma revealed a subset with heterotypic PRO-INV clusters. Our data suggest a framework for the co-existence of these two divergent cell populations, in which heterotypic clusters promote metastasis via cell-cell cooperation.

Melanoma Plasticity: Promoter of Metastasis and Resistance to Therapy

Melanoma is the deadliest form of skin cancer. Although targeted therapies and immunotherapies have revolutionized the treatment of metastatic melanoma, most patients are not cured. Therapy resistance remains a significant clinical challenge. Melanoma comprises phenotypically distinct subpopulations of cells, exhibiting distinct gene signatures leading to tumor heterogeneity and favoring therapeutic resistance. Cellular plasticity in melanoma is referred to as phenotype switching. Regardless of their genomic classification, melanomas switch from a proliferative and differentiated phenotype to an invasive, dedifferentiated and often therapy-resistant state. In this review we discuss potential mechanisms underpinning melanoma phenotype switching, how this cellular plasticity contributes to resistance to both targeted therapies and immunotherapies. Finally, we highlight novel strategies to target plasticity and their potential clinical impact in melanoma.

Human melanocyte development and melanoma dedifferentiation at single-cell resolution

In humans, epidermal melanocytes are responsible for skin pigmentation, defence against ultraviolet radiation and the deadliest common skin cancer, melanoma. Although there is substantial overlap in melanocyte development pathways between different model organisms, species-dependent differences are frequent and the conservation of these processes in human skin remains unresolved. Here, we used a single-cell enrichment and RNA-sequencing pipeline to study human epidermal melanocytes directly from the skin, capturing transcriptomes across different anatomical sites, developmental age, sexes and multiple skin tones. We uncovered subpopulations of melanocytes that exhibit anatomical site-specific enrichment that occurs during gestation and persists through adulthood. The transcriptional signature of the volar-enriched subpopulation is retained in acral melanomas. Furthermore, we identified human melanocyte differentiation transcriptional programs that are distinct from gene signatures generated from model systems. Finally, we used these programs to define patterns of dedifferentiation that are predictive of melanoma prognosis and response to immune checkpoint inhibitor therapy.

Transcriptional Reprogramming and Constitutive PD-L1 Expression in Melanoma Are Associated with Dedifferentiation and Activation of Interferon and Tumour Necrosis Factor Signalling Pathways

Simple Summary

Melanoma, an aggressive form of skin cancer, is frequently associated with drug resistance in the advanced stages. For instance, frequently resistance is observed in sequential treatment of melanoma with targeted therapy and immunotherapy. In this research, the authors investigated whether potential transcriptional mechanisms and pathways associated with PD-L1 protein expression could underlie targeted therapy drug resistance in melanoma. The authors found a PD-L1 expression transcriptional pattern underlies resistance to targeted therapy in a subgroup of melanomas. These melanomas were markedly dedifferentiated, as compared to melanomas that were not drug resistant. Understanding changes in transcription and molecular pathways that lead to drug resistance could allow researchers to develop interventions to prevent drug resistance from occurring in melanoma, which could also be relevant to other cancer types.

Abstract

Melanoma is the most aggressive type of skin cancer, with increasing incidence worldwide. Advances in targeted therapy and immunotherapy have improved the survival of melanoma patients experiencing recurrent disease, but unfortunately treatment resistance frequently reduces patient survival. Resistance to targeted therapy is associated with transcriptomic changes and has also been shown to be accompanied by increased expression of programmed death ligand 1 (PD-L1), a potent inhibitor of immune response. Intrinsic upregulation of PD-L1 is associated with genome-wide DNA hypomethylation and widespread alterations in gene expression in melanoma cell lines. However, an in-depth analysis of the transcriptomic landscape of melanoma cells with intrinsically upregulated PD-L1 expression is lacking. To determine the transcriptomic landscape of intrinsically upregulated PD-L1 expression in melanoma, we investigated transcriptomes in melanomas with constitutive versus inducible PD-L1 expression (referred to as PD-L1_CON and PD-L1_IND). RNA-Seq analysis was performed on seven PD-L1_CON melanoma cell lines and ten melanoma cell lines with low inducible PD-L1_IND expression. We observed that PD-L1_CON melanoma cells had a reprogrammed transcriptome with a characteristic pattern of dedifferentiated gene expression, together with active interferon (IFN) and tumour necrosis factor (TNF) signalling pathways. Furthermore, we identified key transcription factors that were also differentially expressed in PD-L1_CON versus PD-L1_IND melanoma cell lines. Overall, our studies describe transcriptomic reprogramming of melanomas with PD-L1_CON expression.

CDK7 and MITF repress a transcription program involved in survival and drug tolerance in melanoma

Melanoma cell phenotype switching between differentiated melanocytic and undifferentiated mesenchymal-like states drives metastasis and drug resistance. CDK7 is the serine/threonine kinase of the basal transcription factor TFIIH. We show that dedifferentiation of melanocytic-type melanoma cells into mesenchymal-like cells and acquisition of tolerance to targeted therapies is achieved through chronic inhibition of CDK7. In addition to emergence of a mesenchymal-type signature, we identify a GATA6-dependent gene expression program comprising genes such as AMIGO2 or ABCG2 involved in melanoma survival or targeted drug tolerance, respectively. Mechanistically, we show that CDK7 drives expression of the melanocyte lineage transcription factor MITF that in turn binds to an intronic region of GATA6 to repress its expression in melanocytic-type cells. We show that GATA6 expression is activated in MITF-low melanoma cells of patient-derived xenografts. Taken together, our data show how the poorly characterized repressive function of MITF in melanoma participates in a molecular cascade regulating activation of a transcriptional program involved in survival and drug resistance in melanoma.

How Neural Crest Transcription Factors Contribute to Melanoma Heterogeneity, Cellular Plasticity, and Treatment Resistance

Cutaneous melanoma represents one of the deadliest types of skin cancer. The prognosis strongly depends on the disease stage, thus early detection is crucial. New therapies, including BRAF and MEK inhibitors and immunotherapies, have significantly improved the survival of patients in the last decade. However, intrinsic and acquired resistance is still a challenge. In this review, we discuss two major aspects that contribute to the aggressiveness of melanoma, namely, the embryonic origin of melanocytes and melanoma cells and cellular plasticity. First, we summarize the physiological function of epidermal melanocytes and their development from precursor cells that originate from the neural crest (NC). Next, we discuss the concepts of intratumoral heterogeneity, cellular plasticity, and phenotype switching that enable melanoma to adapt to changes in the tumor microenvironment and promote disease progression and drug resistance. Finally, we further dissect the connection of these two aspects by focusing on the transcriptional regulators MSX1, MITF, SOX10, PAX3, and FOXD3. These factors play a key role in NC initiation, NC cell migration, and melanocyte formation, and we discuss how they contribute to cellular plasticity and drug resistance in melanoma.

RNA-seq Characterization of Melanoma Phenotype Switch in 3D Collagen after p38 MAPK Inhibitor Treatment

Melanoma phenotype plasticity underlies tumour dissemination and resistance to therapy, yet its regulation is incompletely understood. In vivo switching between a more differentiated, proliferative phenotype and a dedifferentiated, invasive phenotype is directed by the tumour microenvironment. We found that treatment of partially dedifferentiated, invasive A375M2 cells with two structurally unrelated p38 MAPK inhibitors, SB2021920 and BIRB796, induces a phenotype switch in 3D collagen, as documented by increased expression of melanocyte differentiation markers and a loss of invasive phenotype markers. The phenotype is accompanied by morphological change corresponding to amoeboid–mesenchymal transition. We performed RNA sequencing with an Illumina HiSeq platform to fully characterise transcriptome changes underlying the switch. Gene expression results obtained with RNA-seq were validated by comparing them with RT-qPCR. Transcriptomic data generated in the study will extend the present understanding of phenotype plasticity in melanoma and its contribution to invasion and metastasis.

IER2-induced senescence drives melanoma invasion through osteopontin

Expression of the immediate-early response gene IER2 has been associated with the progression of several types of cancer, but its functional role is poorly understood. We found that increased IER2 expression in human melanoma is associated with shorter overall survival, and subsequently investigated the mechanisms through which IER2 exerts this effect. In experimental melanoma models, sustained expression of IER2 induced senescence in a subset of melanoma cells in a p53/MAPK/AKT-dependent manner. The senescent cells produced a characteristic secretome that included high levels of the extracellular phosphoglycoprotein osteopontin. Nuclear localization of the IER2 protein was critical for both the induction of senescence and osteopontin secretion. Osteopontin secreted by IER2-expressing senescent cells strongly stimulated the migration and invasion of non-senescent melanoma cells. Consistently, we observed coordinate expression of IER2, p53/p21, and osteopontin in primary human melanomas and metastases, highlighting the pathophysiological relevance of IER2-mediated senescence in melanoma progression. Together, our study reveals that sustained IER2 expression drives melanoma invasion and progression through stimulating osteopontin secretion via the stochastic induction of senescence.

Prevention of Melanoma Extravasation as a New Treatment Option Exemplified by p38/MK2 Inhibition

Melanoma releases numerous tumor cells into the circulation; however, only a very small fraction of these cells is able to establish distant metastasis. Intravascular survival of circulating tumor cells is limited through hemodynamic forces and by the lack of matrix interactions. The extravasation step is, thus, of unique importance to establish metastasis. Similar to leukocyte extravasation, this process is under the control of adhesion molecule pairs expressed on melanoma and endothelial cells, and as for leukocytes, ligands need to be adequately presented on cell surfaces. Based on melanoma plasticity, there is considerable heterogeneity even within one tumor and one patient resulting in a mixture of invasive or proliferative cells. The molecular control for this switch is still ill-defined. Recently, the balance between two kinase pathways, p38 and JNK, has been shown to determine growth characteristics of melanoma. While an active JNK pathway induces a proliferative phenotype with reduced invasive features, an active p38/MK2 pathway results in an invasive phenotype and supports the extravasation step via the expression of molecules capable of binding to endothelial integrins. Therapeutic targeting of MK2 to prevent extravasation might reduce metastatic spread.

Cellular Phenotypic Plasticity of Cutaneous Melanoma: A Complex Puzzle

Wenzina et al. (2020) explore the potential role of E-cadherin (CDH1) as a marker for invasive behavior in melanoma. The authors show that CDH1 expression is modulated by p38 signaling, and that manipulation of this pathway can impede endothelial disruption and lung dissemination in vivo and in vitro. The downstream markers PODXL and DEL of the invasive phenotype are associated with a poor prognosis.

Inhibitory effects of cynaropicrin on human melanoma progression by targeting MAPK, NF-κB, and Nrf-2 signaling pathways in vitro

Malignant melanoma is the deadliest skin cancer, due to its propensity to metastasize. MAPKs and NF-κB pathways are constitutively activated in melanoma and promote cell proliferation, cell invasion, metastasis formation, and resistance to therapeutic regimens. Thus, they represent potential targets for melanoma prevention and treatment. Phytochemicals are gaining considerable attention for the management of melanoma because of their several cellular and molecular targets. A screening of a small library of sesquiterpenes lactones selected cynaropicrin, isolated from the aerial parts of Centaurea drabifolia subsp. detonsa, for its potential anticancer effect against melanoma cells. Treatment of human melanoma cells A375 with cynaropicrin resulted in inhibition of cell proliferation and induction of caspase-3-dependent apoptosis. Furthermore, cynaropicrin reduced several cellular malignant features such migration, invasion, and colonies formation through the inhibition of ERK1/2 and NF-κB activity. Cynaropicrin was able to reduce intracellular reactive oxygen species generation, which are involved in all the stages of carcinogenesis. Indeed, cynaropicrin increased the expression of several antioxidant genes, such as glutamate-cysteine ligase and heme oxygenase-1, by promoting the activation of the transcription factor Nrf-2. In conclusion, our results individuate cynaropicrin as a potential adjuvant chemotherapeutic agent for melanoma by targeting several protumorigenic signaling pathways.

Plasma Treatment Limits Human Melanoma Spheroid Growth and Metastasis Independent of the Ambient Gas Composition

Simple Summary

Despite recent advances in therapeutic options, melanoma remains a deadly disease with a poor prognosis. Physical gas plasma has been proposed as a promising technology for the treatment of melanoma. This study aimed to develop and investigate a convenient test system based on three-dimensional cell cultures (spheroids) of two melanoma cell lines in response to physical gas plasma. The experimental approach combined high-content imaging technology and different gas plasma treatment modalities (direct and indirect, gas compositions). Our results revealed that plasma treatment was toxic for both cell lines predominantly dependent on the treatment time. Furthermore, we addressed the question of safety and morphological changes in response to physical gas plasma exposure and found no support for metastatic progression. Treatment with physical gas plasma effectively limited the growth of human 3D melanoma spheroids and provided a versatile test system for more in vivo-like tumor tissue.

Abstract

Melanoma skin cancer is still a deadly disease despite recent advances in therapy. Previous studies have suggested medical plasma technology as a promising modality for melanoma treatment. However, the efficacy of plasmas operated under different ambient air conditions and the comparison of direct and indirect plasma treatments are mostly unexplored for this tumor entity. Moreover, exactly how plasma treatment affects melanoma metastasis has still not been explained. Using 3D tumor spheroid models and high-content imaging technology, we addressed these questions by utilizing one metastatic and one non-metastatic human melanoma cell line targeted with an argon plasma jet. Plasma treatment was toxic in both cell lines. Modulating the oxygen and nitrogen ambient air composition (100/0, 75/25, 50/50, 25/75, and 0/100) gave similar toxicity and reduced the spheroid growth for all conditions. This was the case for both direct and indirect treatments, with the former showing a treatment time-dependent response while the latter resulted in cytotoxicity with the longest treatment time investigated. Live-cell imaging of in-gel cultured spheroids indicated that plasma treatment did not enhance metastasis, and flow cytometry showed a significant modulation of S100A4 but not in any of the five other metastasis-related markers (β-catenin, E-cadherin, LEF1, SLUG, and ZEB1) investigated.

The Stress-Like Cancer Cell State Is a Consistent Component of Tumorigenesis

Transcriptional profiling of tumors has revealed a stress-like state among the cancer cells with the concerted expression of genes such as fos, jun, and heat-shock proteins, though this has been controversial given possible dissociation-effects associated with single-cell RNA sequencing. Here, we validate the existence of this state using a combination of zebrafish melanoma modeling, spatial transcriptomics, and human samples. We found that the stress-like subpopulation of cancer cells is present from the early stages of tumorigenesis. Comparing with previously reported single-cell RNA sequencing datasets from diverse cancer types, including triple-negative breast cancer, oligodendroglioma, and pancreatic adenocarcinoma, indicated the conservation of this state during tumorigenesis. We also provide evidence that this state has higher tumor-seeding capabilities and that its induction leads to increased growth under both MEK and BRAF inhibitors. Collectively, our study supports the stress-like cells as a cancer cell state expressing a coherent set of genes and exhibiting drug-resistance properties.