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

STAT3 promotes melanoma metastasis by CEBP-induced repression of the MITF pathway

Metastatic melanoma is hallmarked by its ability of phenotype switching to more slowly proliferating, but highly invasive cells. Here, we tested the impact of signal transducer and activator of transcription 3 (STAT3) on melanoma progression in association with melanocyte inducing transcription factor (MITF) expression levels. We established a mouse melanoma model for deleting Stat3 in melanocytes with specific expression of human hyperactive NRAS^Q61K in an Ink4a deficient background, two frequent driver mutations in human melanoma. Mice devoid of Stat3 showed early disease onset with higher proliferation in primary tumors, but displayed significantly diminished lung, brain and liver metastases. Whole genome expression profiling of tumor-derived cells also showed a reduced invasion phenotype, which was further corroborated by 3D melanoma model analysis. Notably, loss or knockdown of STAT3 in mouse or human cells resulted in up-regulation of MITF and induction of cell proliferation. Mechanistically we show that STAT3-induced CEBPa/b expression was sufficient to suppress MITF transcription. Epigenetic analysis by ATAC-seq confirmed that CEBPa/b binding to the MITF enhancer region silenced the MITF locus. Finally, by classification of patient-derived melanoma samples, we show that STAT3 and MITF act antagonistically and hence contribute differentially to melanoma progression. We conclude that STAT3 is a driver of the metastatic process in melanoma and able to antagonize MITF via direct induction of CEBP family member transcription.

Conserved Epigenetic Regulatory Logic Infers Genes Governing Cell Identity

Determining genes that orchestrate cell differentiation in development and disease remains a fundamental goal of cell biology. This study establishes a genome-wide metric based on the gene-repressive trimethylation of histone H3 at lysine 27 (H3K27me3) across hundreds of diverse cell types to identify genetic regulators of cell differentiation. We introduce a computational method, TRIAGE, which uses discordance between gene-repressive tendency and expression to identify genetic drivers of cell identity. We apply TRIAGE to millions of genome-wide single-cell transcriptomes, diverse omics platforms, and eukaryotic cells and tissue types. Using a wide range of data, we validate the performance of TRIAGE in identifying cell-type-specific regulatory factors across diverse species including human, mouse, boar, bird, fish, and tunicate. Using CRISPR gene editing, we use TRIAGE to experimentally validate RNF220 as a regulator of Ciona cardiopharyngeal development and SIX3 as required for differentiation of endoderm in human pluripotent stem cells. A record of this paper’s transparent peer review process is included in the Supplemental Information.

Perturbing genes controlling cell decisions have major implications in development or disease. However, identifying key regulatory genes from the thousands expressed in a cell is challenging. TRIAGE is a computational method that distills patterns of epigenetic repression across diverse cell types to infer regulatory genes using input gene expression data from any cell type. Demonstrating its utility, we combine single-cell RNA-seq and TRIAGE to identify and experimentally confirm novel regulators of heart development in evolutionarily distant species.

Cross-species analysis of enhancer logic using deep learning

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

Fra-2/AP-1 regulates melanoma cell metastasis by downregulating Fam212b

Metastatic melanoma remains a challenging disease. Understanding the molecular mechanisms how melanoma becomes metastatic is therefore of interest. Herein we show that downregulation of the AP-1 transcription factor member Fra-2 in melanoma cells is associated with an aggressive melanoma phenotype in vitro and in vivo. In vitro, Fra-2 knockdown in melanoma cells promoted cell migration and invasion associated with increased Snail-1, Twist-1/2, and matrix metalloproteinase-2 (MMP-2) expression. In vivo, Fra-2 knockdown in a melanoma cell line led to increased metastasis into the lungs and liver. The increased metastatic potential of Fra-2 knockdown melanoma cells was likely due to an accelerated cell cycle transition and increased tissue angiogenesis. Using Fra-2 knockdown cell lines microarray analysis, we identified the protein Fam212b (family with sequence similarity 212 member B) as a downstream target of Fra-2. By additional knockdown of Fam212b in Fra-2 mutant cells, we mitigated the cell migration, invasion, and cell cycle transition phenotype induced by Fra-2 knockdown. Furthermore, Fam212b overexpression enhanced β-catenin pathway. Finally, Fam212b expression is correlated with increased melanoma metastasis and poor clinical outcomes in human patients. In summary, these findings reveal the Fra-2-Fam212b axis as a new pathway of melanoma metastasis, which can be in the future used as potential marker of the metastatic properties of melanoma.

Human Endogenous Retrovirus K Rec forms a Regulatory Loop with MITF that Opposes the Progression of Melanoma to an Invasive Stage

The HML2 subfamily of HERV-K (henceforth HERV-K) represents the most recently endogenized retrovirus in the human genome. While the products of certain HERV-K genomic copies are expressed in normal tissues, they are upregulated in several pathological conditions, including various tumors. It remains unclear whether HERV-K(HML2)-encoded products overexpressed in cancer contribute to disease progression or are merely by-products of tumorigenesis. Here, we focus on the regulatory activities of the Long Terminal Repeats (LTR5_Hs) of HERV-K and the potential role of the HERV-K-encoded Rec in melanoma. Our regulatory genomics analysis of LTR5_Hs loci indicates that Melanocyte Inducing Transcription Factor (MITF) (also known as binds to a canonical E-box motif (CA(C/T)GTG) within these elements in proliferative type of melanoma, and that depletion of MITF results in reduced HERV-K expression. In turn, experimentally depleting Rec in a proliferative melanoma cell line leads to lower mRNA levels of MITF and its predicted target genes. Furthermore, Rec knockdown leads to an upregulation of epithelial-to-mesenchymal associated genes and an enhanced invasion phenotype of proliferative melanoma cells. Together these results suggest the existence of a regulatory loop between MITF and Rec that may modulate the transition from proliferative to invasive stages of melanoma. Because HERV-K(HML2) elements are restricted to hominoid primates, these findings might explain certain species-specific features of melanoma progression and point to some limitations of animal models in melanoma studies.

ShcD Binds DOCK4, Promotes Ameboid Motility and Metastasis Dissemination, Predicting Poor Prognosis in Melanoma

Simple Summary

Metastasis formation and dissemination is a complex process that relies on several steps. Even though highly inefficient, metastasis spreading is the primary cause of cancer morbidity and mortality in patients. The aim of our study was to investigate the molecular pathways leading to metastases making use of human-in-mouse melanoma models of patient-derived xenografts. We demonstrate that the modulation of the expression of an adaptor protein of the Shc family, ShcD, can change the phenotype and the invasive properties of melanoma cells when highly expressed. We also show that ShcD binds DOCK4 and confines it into the cytoplasm, blocking the Rac1 signaling pathways, thus leading to metastasis development. Moreover, our results indicate that melanoma cells are more sensitive to therapeutic treatments when the ShcD molecular pathway is inactivated, suggesting that new therapeutic strategies can be designed in melanomas.

Abstract

Metastases are the primary cause of cancer-related deaths. The underlying molecular and biological mechanisms remain, however, elusive, thus preventing the design of specific therapies. In melanomas, the metastatic process is influenced by the acquisition of metastasis-associated mutational and epigenetic traits and the activation of metastatic-specific signaling pathways in the primary melanoma. In the current study, we investigated the role of an adaptor protein of the Shc family (ShcD) in the acquisition of metastatic properties by melanoma cells, exploiting our cohort of patient-derived xenografts (PDXs). We provide evidence that the depletion of ShcD expression increases a spread cell shape and the capability of melanoma cells to attach to the extracellular matrix while its overexpression switches their morphology from elongated to rounded on 3D matrices, enhances cells’ invasive phenotype, as observed on collagen gel, and favors metastasis formation in vivo. ShcD overexpression sustains amoeboid movement in melanoma cells, by suppressing the Rac1 signaling pathway through the confinement of DOCK4 in the cytoplasm. Inactivation of the ShcD signaling pathway makes melanoma cells more sensitive to therapeutic treatments. Consistently, ShcD expression predicts poor outcome in a cohort of 183 primary melanoma patients.

Pan-cancer single-cell RNA-seq identifies recurring programs of cellular heterogeneity

Cultured cell lines are the workhorse of cancer research, but the extent to which they recapitulate the heterogeneity observed among malignant cells in tumors is unclear. Here we used multiplexed single-cell RNA-seq to profile 198 cancer cell lines from 22 cancer types. We identified 12 expression programs that are recurrently heterogeneous within multiple cancer cell lines. These programs are associated with diverse biological processes, including cell cycle, senescence, stress and interferon responses, epithelial-mesenchymal transition and protein metabolism. Most of these programs recapitulate those recently identified as heterogeneous within human tumors. We prioritized specific cell lines as models of cellular heterogeneity and used them to study subpopulations of senescence-related cells, demonstrating their dynamics, regulation and unique drug sensitivities, which were predictive of clinical response. Our work describes the landscape of heterogeneity within diverse cancer cell lines and identifies recurrent patterns of heterogeneity that are shared between tumors and specific cell lines.

WNT11-FZD7-DAAM1 signalling supports tumour initiating abilities and melanoma amoeboid invasion

Melanoma is a highly aggressive tumour that can metastasize very early in disease progression. Notably, melanoma can disseminate using amoeboid invasive strategies. We show here that high Myosin II activity, high levels of ki-67 and high tumour-initiating abilities are characteristic of invasive amoeboid melanoma cells. Mechanistically, we find that WNT11-FZD7-DAAM1 activates Rho-ROCK1/2-Myosin II and plays a crucial role in regulating tumour-initiating potential, local invasion and distant metastasis formation. Importantly, amoeboid melanoma cells express both proliferative and invasive gene signatures. As such, invasive fronts of human and mouse melanomas are enriched in amoeboid cells that are also ki-67 positive. This pattern is further enhanced in metastatic lesions. We propose eradication of amoeboid melanoma cells after surgical removal as a therapeutic strategy.

A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues

Background

Although considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues.

Results

Overall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions.

Conclusions

The lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

Enhancer Reprogramming Confers Dependence on Glycolysis and IGF Signaling in KMT2D Mutant Melanoma

Histone methyltransferase KMT2D harbors frequent loss-of-function somatic point mutations in several tumor types, including melanoma. Here, we identify KMT2D as a potent tumor suppressor in melanoma through an in vivo epigenome-focused pooled RNAi screen and confirm the finding by using a genetically engineered mouse model (GEMM) based on conditional and melanocyte-specific deletion of KMT2D. KMT2D-deficient tumors show substantial reprogramming of key metabolic pathways, including glycolysis. KMT2D deficiency aberrantly upregulates glycolysis enzymes, intermediate metabolites, and glucose consumption rates. Mechanistically, KMT2D loss causes genome-wide reduction of H3K4me1-marked active enhancer chromatin states. Enhancer loss and subsequent repression of IGFBP5 activates IGF1R-AKT to increase glycolysis in KMT2D-deficient cells. Pharmacological inhibition of glycolysis and insulin growth factor (IGF) signaling reduce proliferation and tumorigenesis preferentially in KMT2D-deficient cells. We conclude that KMT2D loss promotes tumorigenesis by facilitating an increased use of the glycolysis pathway for enhanced biomass needs via enhancer reprogramming, thus presenting an opportunity for therapeutic intervention through glycolysis or IGF pathway inhibitors.

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.

Insights into Differentiation of Melanocytes from Human Stem Cells and Their Relevance for Melanoma Treatment

Malignant melanoma represents a highly aggressive form of skin cancer. The metastatic process itself is mostly governed by the so-called epithelial mesenchymal transition (EMT), which confers cancer cells migrative, invasive and resistance abilities. Since EMT represents a conserved developmental process, it is worthwhile further examining the nature of early developmental steps fundamental for melanocyte differentiation. This can be done either in vivo by analyzing the physiologic embryo development in different species or by in vitro studies of melanocytic differentiation originating from embryonic human stem cells. Most importantly, external cues drive progenitor cell differentiation, which can be divided in stages favoring neural crest specification or melanocytic differentiation and proliferation. In this review, we describe ectopic factors which drive human pluripotent stem cell differentiation to melanocytes in 2D, as well as in organoid models. Furthermore, we compare developmental mechanisms with processes described to occur during melanoma development. Finally, we suggest differentiation factors as potential co-treatment options for metastatic melanoma patients.

The aryl hydrocarbon receptor promotes inflammation-induced dedifferentiation and systemic metastatic spread of melanoma cells

The aryl hydrocarbon receptor (AHR) is a ligand binding-transcription factor of the basic helix-loop-helix family regulating multiple cellular functions such as differentiation, cell cycle, apoptosis, and inflammatory reactions. In neoplastic diseases, the AHR has been described to modulate proliferation and differentiation in dichotomous ways, either inhibiting or augmenting the growth of tumors. The precise role of AHR in melanoma is mostly unknown. Here, we report a functional effect of AHR activation on inflammation-induced melanoma cell dedifferentiation and the development of lung metastases in a mouse model. Via in silico analyses of "The Cancer Genome Atlas" human melanoma cohort, we detected a correlation between AHR expression levels and a dedifferentiated melanoma cell phenotype with an invasive gene signature, which we were able to functionally recapitulate in a panel of human melanoma cell lines. Both human and mouse melanoma cell lines upregulated AHR expression after inflammatory stimulation with tumor necrosis factor-α (TNF-α). Activation of AHR in human and mouse melanoma cell lines with the endogenous ligand formylindolo(3,2-b)carbazole (FICZ) promoted inflammation-induced dedifferentiation in vitro. Importantly, mouse melanoma cells with CRISPR/Cas9-mediated disruption of the AHR gene showed impaired in vivo tumor growth after transplantation in the skin as well as decreased numbers of spontaneous lung metastases. Taken together, our results demonstrate a functional role for AHR expression in melanoma development and metastatic progression. This provides a scientific basis for future experiments that further dissect the underlying molecular mechanisms and assess the potential for AHR inhibition as part of multimodal melanoma treatment strategies.

New Insights into the Role of Sphingolipid Metabolism in Melanoma

Cutaneous melanoma is a deadly skin cancer whose aggressiveness is directly linked to its metastatic potency. Despite remarkable breakthroughs in term of treatments with the emergence of targeted therapy and immunotherapy, the prognosis for metastatic patients remains uncertain mainly because of resistances. Better understanding the mechanisms responsible for melanoma progression is therefore essential to uncover new therapeutic targets. Interestingly, the sphingolipid metabolism is dysregulated in melanoma and is associated with melanoma progression and resistance to treatment. This review summarises the impact of the sphingolipid metabolism on melanoma from the initiation to metastatic dissemination with emphasis on melanoma plasticity, immune responses and resistance to treatments.

Sometimes even apoptosis fails: implications for cancer

Apoptosis is vital for the correct morphogenesis of multi-cellular organisms. However, like most physiological programs, the cell's ability to commit suicide is hijacked by cancer in its own proliferative and invasive interest. We recently showed that inefficient execution of apoptosis (or failed apoptosis) is used by cancer to boost invasiveness.

Memory Sequencing Reveals Heritable Single-Cell Gene Expression Programs Associated with Distinct Cellular Behaviors

Non-genetic factors can cause individual cells to fluctuate substantially in gene expression levels over time. It remains unclear whether these fluctuations can persist for much longer than the time of one cell division. Current methods for measuring gene expression in single cells mostly rely on single time point measurements, making the duration of gene expression fluctuations or cellular memory difficult to measure. Here, we combined Luria and Delbrück's fluctuation analysis with population-based RNA sequencing (MemorySeq) for identifying genes transcriptome-wide whose fluctuations persist for several divisions. MemorySeq revealed multiple gene modules that expressed together in rare cells within otherwise homogeneous clonal populations. These rare cell subpopulations were associated with biologically distinct behaviors like proliferation in the face of anti-cancer therapeutics. The identification of non-genetic, multigenerational fluctuations can reveal new forms of biological memory in single cells and suggests that non-genetic heritability of cellular state may be a quantitative property.

Intrinsic Balance between ZEB Family Members Is Important for Melanocyte Homeostasis and Melanoma Progression

It has become clear that cellular plasticity is a main driver of cancer therapy resistance. Consequently, there is a need to mechanistically identify the factors driving this process. The transcription factors of the zinc-finger E-box-binding homeobox family, consisting of ZEB1 and ZEB2, are notorious for their roles in epithelial-to-mesenchymal transition (EMT). However, in melanoma, an intrinsic balance between ZEB1 and ZEB2 seems to determine the cellular state by modulating the expression of the master regulator of melanocyte homeostasis, microphthalmia-associated transcription factor (MITF). ZEB2 drives MITF expression and is associated with a differentiated/proliferative melanoma cell state. On the other hand, ZEB1 is correlated with low MITF expression and a more invasive, stem cell-like and therapy-resistant cell state. This intrinsic balance between ZEB1 and ZEB2 could prove to be a promising therapeutic target for melanoma patients. In this review, we will summarise what is known on the functional mechanisms of these transcription factors. Moreover, we will look specifically at their roles during melanocyte-lineage development and homeostasis. Finally, we will overview the current literature on ZEB1 and ZEB2 in the melanoma context and link this to the 'phenotype-switching' model of melanoma cellular plasticity.

Reversal of pre-existing NGFR-driven tumor and immune therapy resistance

Melanomas can switch to a dedifferentiated cell state upon exposure to cytotoxic T cells. However, it is unclear whether such tumor cells pre-exist in patients and whether they can be resensitized to immunotherapy. Here, we chronically expose (patient-derived) melanoma cell lines to differentiation antigen-specific cytotoxic T cells and observe strong enrichment of a pre-existing NGFR^hi population. These fractions are refractory also to T cells recognizing non-differentiation antigens, as well as to BRAF + MEK inhibitors. NGFR^hi cells induce the neurotrophic factor BDNF, which contributes to T cell resistance, as does NGFR. In melanoma patients, a tumor-intrinsic NGFR signature predicts anti-PD-1 therapy resistance, and NGFR^hi tumor fractions are associated with immune exclusion. Lastly, pharmacologic NGFR inhibition restores tumor sensitivity to T cell attack in vitro and in melanoma xenografts. These findings demonstrate the existence of a stable and pre-existing NGFR^hi multitherapy-refractory melanoma subpopulation, which ought to be eliminated to revert intrinsic resistance to immunotherapeutic intervention.

EMT-Inducing Transcription Factors, Drivers of Melanoma Phenotype Switching, and Resistance to Treatment

Transcription factors, extensively described for their role in epithelial-mesenchymal transition (EMT-TFs) in epithelial cells, also display essential functions in the melanocyte lineage. Recent evidence has shown specific expression patterns and functions of these EMT-TFs in neural crest-derived melanoma compared to carcinoma. Herein, we present an update of the specific roles of EMT-TFs in melanocyte differentiation and melanoma progression. As major regulators of phenotype switching between differentiated/proliferative and neural crest stem cell-like/invasive states, these factors appear as major drivers of intra-tumor heterogeneity and resistance to treatment in melanoma, which opens new avenues in terms of therapeutic targeting.

Robust gene expression programs underlie recurrent cell states and phenotype switching in melanoma

Melanoma cells can switch between a melanocytic and a mesenchymal-like state. Scattered evidence indicates that additional intermediate state(s) may exist. Here, to search for such states and decipher their underlying gene regulatory network (GRN), we studied 10 melanoma cultures using single-cell RNA sequencing (RNA-seq) as well as 26 additional cultures using bulk RNA-seq. Although each culture exhibited a unique transcriptome, we identified shared GRNs that underlie the extreme melanocytic and mesenchymal states and the intermediate state. This intermediate state is corroborated by a distinct chromatin landscape and is governed by the transcription factors SOX6, NFATC2, EGR3, ELF1 and ETV4. Single-cell migration assays confirmed the intermediate migratory phenotype of this state. Using time-series sampling of single cells after knockdown of SOX10, we unravelled the sequential and recurrent arrangement of GRNs during phenotype switching. Taken together, these analyses indicate that an intermediate state exists and is driven by a distinct and stable 'mixed' GRN rather than being a symbiotic heterogeneous mix of cells.

Genome-wide cooperation of EMT transcription factor ZEB1 with YAP and AP-1 in breast cancer

Invasion, metastasis and therapy resistance are the major cause of cancer‐associated deaths, and the EMT‐inducing transcription factor ZEB1 is a crucial stimulator of these processes. While work on ZEB1 has mainly focused on its role as a transcriptional repressor, it can also act as a transcriptional activator. To further understand these two modes of action, we performed a genome‐wide ZEB1 binding study in triple‐negative breast cancer cells. We identified ZEB1 as a novel interactor of the AP‐1 factors FOSL1 and JUN and show that, together with the Hippo pathway effector YAP, they form a transactivation complex, predominantly activating tumour‐promoting genes, thereby synergising with its function as a repressor of epithelial genes. High expression of ZEB1, YAP, FOSL1 and JUN marks the aggressive claudin‐low subtype of breast cancer, indicating the translational relevance of our findings. Thus, our results link critical tumour‐promoting transcription factors: ZEB1, AP‐1 and Hippo pathway factors. Disturbing their molecular interaction may provide a promising treatment option for aggressive cancer types.

Role of Flavonoids in the Prevention of AhR-Dependent Resistance During Treatment with BRAF Inhibitors

BRAF and MEK inhibitors (BRAFi and MEKi) are the standard of care for the treatment of metastatic melanoma in patients with BRAF^V600E mutations, greatly improving progression-free survival. However, the acquisition of resistance to BRAFi and MEKi remains a difficult clinical challenge, with limited therapeutic options available for these patients. Here, we investigated the therapeutic potential of natural flavonoids as specific AhR (Aryl hydrocarbon Receptor) transcription factor antagonists in combination with BRAFi. Experimental Design: Experiments were performed in vitro and in vivo with various human melanoma cell lines (mutated for BRAF^V600E) sensitive or resistant to BRAFi. We evaluated the role of various flavonoids on cell sensitivity to BRAFi and their ability to counteract resistance and the invasive phenotype of melanoma. Results: Flavonoids were highly effective in potentiating BRAFi therapy in human melanoma cell lines by increasing sensitivity and delaying the pool of resistant cells that arise during treatment. As AhR antagonists, flavonoids counteracted a gene expression program associated with the acquisition of resistance and phenotype switching that leads to an invasive and EMT-like phenotype. Conclusions: The use of natural flavonoids opens new therapeutic opportunities for the treatment of patients with BRAF-resistant disease.

Dysregulation of MITF Leads to Transformation in MC1R-Defective Melanocytes

The MC1R/cAMP/MITF pathway is a key determinant for growth, differentiation, and survival of melanocytes and melanoma. MITF-M is the melanocyte-specific isoform of Microphthalmia-associated Transcription Factor (MITF) in human melanoma. Here we use two melanocyte cell lines to show that forced expression of hemagglutinin (HA) -tagged MITF-M through lentiviral transduction represents an oncogenic insult leading to consistent cell transformation of the immortalized melanocyte cell line Hermes 4C, being a melanocortin-1 receptor (MC1R) compound heterozygote, while not causing transformation of the MC1R wild type cell line Hermes 3C. The transformed HA-tagged MITF-M transduced Hermes 4C cells form colonies in soft agar and tumors in mice. Further, Hermes 4C cells display increased MITF chromatin binding, and transcriptional reprogramming consistent with an invasive melanoma phenotype. Mechanistically, forced expression of MITF-M drives the upregulation of the AXL tyrosine receptor kinase (AXL), with concomitant downregulation of phosphatase and tensin homolog (PTEN), leading to increased activation of the PI3K/AKT pathway. Treatment with AXL inhibitors reduces growth of the transformed cells by reverting AKT activation. In conclusion, we present a model system of melanoma development, driven by MITF-M in the context of MC1R loss of function, and independent of UV exposure. This model provides a basis for further studies of critical changes in the melanocyte transformation process.

The Hippo pathway oncoprotein YAP promotes melanoma cell invasion and spontaneous metastasis

Melanoma is a deadly form of skin cancer that accounts for a disproportionally large proportion of cancer-related deaths in younger people. Compared with most other skin cancers, a feature of melanoma is its high metastatic capacity, although the mechanisms that confer this are not well understood. The Hippo pathway is a key regulator of organ growth and cell fate that is deregulated in many cancers. To analyse the Hippo pathway in cutaneous melanoma, we generated a transcriptional signature of melanoma cells that overexpressed YAP, the key downstream Hippo pathway oncoprotein. YAP-mediated transcriptional activity varied in melanoma cell lines but did not cluster with known genetic drivers of melanomagenesis such as BRAF and NRAS mutations. Instead, it correlated strongly with published gene expression profiles linked to melanoma cell invasiveness and varied throughout the metastatic cascade in melanoma patient tumours. Consistent with this, YAP was both necessary and sufficient for melanoma cell invasion in vitro. In vivo, YAP promoted spontaneous melanoma metastasis, whilst the growth of YAP-expressing primary tumours was impeded. Finally, we identified the YAP target genes AXL, THBS1 and CYR61 as key mediators of YAP-induced melanoma cell invasion. These data suggest that YAP is a critical regulator of melanoma metastasis.

Tuning Transcription Factor Availability through Acetylation-Mediated Genomic Redistribution

It is widely assumed that decreasing transcription factor DNA-binding affinity reduces transcription initiation by diminishing occupancy of sequence-specific regulatory elements. However, in vivo transcription factors find their binding sites while confronted with a large excess of low-affinity degenerate motifs. Here, using the melanoma lineage survival oncogene MITF as a model, we show that low-affinity binding sites act as a competitive reservoir in vivo from which transcription factors are released by mitogen-activated protein kinase (MAPK)-stimulated acetylation to promote increased occupancy of their regulatory elements. Consequently, a low-DNA-binding-affinity acetylation-mimetic MITF mutation supports melanocyte development and drives tumorigenesis, whereas a high-affinity non-acetylatable mutant does not. The results reveal a paradoxical acetylation-mediated molecular clutch that tunes transcription factor availability via genome-wide redistribution and couples BRAF to tumorigenesis. Our results further suggest that p300/CREB-binding protein-mediated transcription factor acetylation may represent a common mechanism to control transcription factor availability.

•

Reducing transcription factor DNA-binding affinity increases activity in vivo

•

Acetylation is triggered by MAPK signaling

•

Acetylation leads to genome-wide transcription factor redistribution

•

Acetylation of MITF drives tumorigenesis and melanocyte development

The EMT Transcription Factor ZEB2 Promotes Proliferation of Primary and Metastatic Melanoma While Suppressing an Invasive, Mesenchymal-Like Phenotype

Epithelial-to-mesenchymal transition (EMT)-inducing transcription factors (TF) are well known for their ability to induce mesenchymal states associated with increased migratory and invasive properties. Unexpectedly, nuclear expression of the EMT-TF ZEB2 in human primary melanoma has been shown to correlate with reduced invasion. We report here that ZEB2 is required for outgrowth for primary melanomas and metastases at secondary sites. Ablation of Zeb2 hampered outgrowth of primary melanomas in vivo, whereas ectopic expression enhanced proliferation and growth at both primary and secondary sites. Gain of Zeb2 expression in pulmonary-residing melanoma cells promoted the development of macroscopic lesions. In vivo fate mapping made clear that melanoma cells undergo a conversion in state where ZEB2 expression is replaced by ZEB1 expression associated with gain of an invasive phenotype. These findings suggest that reversible switching of the ZEB2/ZEB1 ratio enhances melanoma metastatic dissemination. SIGNIFICANCE: ZEB2 function exerts opposing behaviors in melanoma by promoting proliferation and expansion and conversely inhibiting invasiveness, which could be of future clinical relevance. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/14/2983/F1.large.jpg.

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma

Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAF^V600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

There are more than 20 known melanoma susceptibility genes. Here, using a massively parallel reporter assay, the authors identify risk-associated variants that alter gene transcription, and demonstrate that expression of one such gene, MX2, leads to the promotion of melanoma in a zebrafish model.

The Hippo Pathway as a Driver of Select Human Cancers

The Hippo pathway regulates myriad biological processes in diverse species and is a key cancer signaling network in humans. Although Hippo has been linked to multiple aspects of cancer, its role in this disease is incompletely understood. Large-scale pan-cancer analyses of core Hippo pathway genes reveal that the pathway is mutated at a high frequency only in select human cancers, including malignant mesothelioma and meningioma. Hippo pathway deregulation is also enriched in squamous epithelial cancers. We discuss cancer-related functions of the Hippo pathway and potential explanations for the cancer-restricted mutation profile of core Hippo pathway genes. Greater understanding of Hippo pathway deregulation in cancers will be essential to guide the imminent use of Hippo-targeted therapies.

Epigenetics and metabolism at the crossroads of stress-induced plasticity, stemness and therapeutic resistance in cancer

Despite the recent advances in the treatment of cancers, acquired drug resistance remains a major challenge in cancer management. While earlier studies suggest Darwinian factors driving acquired drug resistance, recent studies point to a more dynamic process involving phenotypic plasticity and tumor heterogeneity in the evolution of acquired drug resistance. Chronic stress after drug treatment induces intrinsic cellular reprogramming and cancer stemness through a slow-cycling persister state, which subsequently drives cancer progression. Both epigenetic and metabolic mechanisms play an important role in this dynamic process. In this review, we discuss how epigenetic and metabolic reprogramming leads to stress-induced phenotypic plasticity and acquired drug resistance, and how the two reprogramming mechanisms crosstalk with each other.

Differential responsiveness to BRAF inhibitors of melanoma cell lines BRAF V600E-mutated

Background

Most mutations in melanoma affect one critical amino acid on BRAF gene, resulting in the V600E substitution. Patient management is often based on the use of specific inhibitors targeting this mutation.

Methods

DNA and RNA mutation status was assessed in 15 melanoma cell lines by Sanger sequencing and RNA-seq. We tested the cell lines responsiveness to BRAF inhibitors (vemurafenib and PLX4720, BRAF-specific and sorafenib, BRAF non-specific). Cell proliferation was assessed by MTT colorimetric assay. BRAF V600E RNA expression was assessed by qPCR. Expression level of phosphorylated-ERK protein was assessed by Western Blotting as marker of BRAF activation.

Results

Three cell lines were discordant in the mutation detection (BRAF V600E at DNA level/Sanger sequencing and BRAF WT on RNA-seq). We initially postulated that those cell lines may express only the WT allele at the RNA level although mutated at the DNA level. A more careful analysis showed that they express low level of BRAF RNA and the expression may be in favor of the WT allele. We tested whether the discordant cell lines responded differently to BRAF-specific inhibitors. Their proliferation rate decreased after treatment with vemurafenib and PLX4720 but was not affected by sorafenib, suggesting a BRAF V600E biological behavior. Yet, responsiveness to the BRAF specific inhibitors was lower as compared to the control. Western Blot analysis revealed a decreased expression of p-ERK protein in the BRAF V600E control cell line and in the discordant cell lines upon treatment with BRAF-specific inhibitors. The discordant cell lines showed a lower responsiveness to BRAF inhibitors when compared to the BRAF V600E control cell line. The results obtained from the inhibition experiment and molecular analyses were also confirmed in three additional cell lines.

Conclusion

Cell lines carrying V600E mutation at the DNA level may respond differently to BRAF targeted treatment potentially due to a lower V600E RNA expression.

Tankyrase inhibition sensitizes melanoma to PD-1 immune checkpoint blockade in syngeneic mouse models

The development of immune checkpoint inhibitors represents a major breakthrough in cancer therapy. Nevertheless, a substantial number of patients fail to respond to checkpoint pathway blockade. Evidence for WNT/β-catenin signaling-mediated immune evasion is found in a subset of cancers including melanoma. Currently, there are no therapeutic strategies available for targeting WNT/β-catenin signaling. Here we show that a specific small-molecule tankyrase inhibitor, G007-LK, decreases WNT/β-catenin and YAP signaling in the syngeneic murine B16-F10 and Clone M-3 melanoma models and sensitizes the tumors to anti-PD-1 immune checkpoint therapy. Mechanistically, we demonstrate that the synergistic effect of tankyrase and checkpoint inhibitor treatment is dependent on loss of β-catenin in the tumor cells, anti-PD-1-stimulated infiltration of T cells into the tumor and induction of an IFNγ- and CD8⁺ T cell-mediated anti-tumor immune response. Our study uncovers a combinatorial therapeutical strategy using tankyrase inhibition to overcome β-catenin-mediated resistance to immune checkpoint blockade in melanoma.

Targeted chemotherapy overcomes drug resistance in melanoma

In this study, Yue et al. describe a therapeutic strategy termed “targeted chemotherapy” that involves depleting PP2A or inhibiting it using a small molecule inhibitor, phendione, in drug-resistant melanoma. The authors show phendione induces DNA damage response without causing DNA breaks or inducing cellular dormancy, therefore blocking tumor growth of BRAF mutant and NRAS mutant melanomas.

The emergence of drug resistance is a major obstacle for the success of targeted therapy in melanoma. Additionally, conventional chemotherapy has not been effective as drug-resistant cells escape lethal DNA damage effects by inducing growth arrest commonly referred to as cellular dormancy. We present a therapeutic strategy termed “targeted chemotherapy” by depleting protein phosphatase 2A (PP2A) or its inhibition using a small molecule inhibitor (1,10-phenanthroline-5,6-dione [phendione]) in drug-resistant melanoma. Targeted chemotherapy induces the DNA damage response without causing DNA breaks or allowing cellular dormancy. Phendione treatment reduces tumor growth of BRAF^V600E-driven melanoma patient-derived xenografts (PDX) and diminishes growth of NRAS^Q61R-driven melanoma, a cancer with no effective therapy. Remarkably, phendione treatment inhibits the acquisition of resistance to BRAF inhibition in BRAF^V600E PDX highlighting its effectiveness in combating the advent of drug resistance.

Clinical Implications of Primary Cilia in Skin Cancer

The primary cilium is a cell surface organelle that is an important component of cellular biology. While it was once believed to be a vestigial structure without biologic function, it is now known to have essential roles in critical cellular signaling pathways such as Hedgehog (HH) and Wnt. The HH and Wnt pathways are involved in pathogenesis of basal cell carcinoma and melanoma, respectively, and this knowledge is now beginning to inform therapeutic and diagnostic options for patients. The purpose of this review is to familiarize clinicians with primary cilia biology and how this complex cellular organelle has started to translate into clinical care.

BRAFi induced demethylation of miR-152-5p regulates phenotype switching by targeting TXNIP in cutaneous melanoma

Treatment of advanced BRAFV600-mutant melanoma using BRAF inhibitors (BRAFi) eventually leads to drug resistance and selects for highly metastatic tumor cells. We compared the most differentially dysregulated miRNA expression profiles of vemurafenib-resistant and highly-metastatic melanoma cell lines obtained from GEO DataSets. We discovered miR-152-5p was a potential regulator mediating melanoma drug resistance and metastasis. Functionally, knockdown of miR-152-5p significantly compromised the metastatic ability of BRAFi-resistant melanoma cells and overexpression of miR-152-5p promoted the formation of slow-cycling phenotype. Furthermore, we explored the cause of how and why miR-152-5p affected metastasis in depth. Mechanistically, miR-152-5p targeted TXNIP which affected metastasis and BRAFi altered the methylation status of MIR152 promoter. Our study highlights the crucial role of miR-152-5p on melanoma metastasis after BRAFi treatment and holds significant implying that discontinuous dosing strategy may improve the benefit of advanced BRAFV600-mutant melanoma patients.

Cooperative behaviour and phenotype plasticity evolve during melanoma progression

A major challenge for managing melanoma is its tumour heterogeneity based on individual co-existing melanoma cell phenotypes. These phenotypes display variable responses to standard therapies, and they drive individual steps of melanoma progression; hence, understanding their behaviour is imperative. Melanoma phenotypes are defined by distinct transcriptional states, which relate to different melanocyte lineage development phases, ranging from a mesenchymal, neural crest-like to a proliferative, melanocytic phenotype. It is thought that adaptive phenotype plasticity based on transcriptional reprogramming drives melanoma progression, but at which stage individual phenotypes dominate and moreover, how they interact is poorly understood. We monitored melanocytic and mesenchymal phenotypes throughout melanoma progression and detected transcriptional reprogramming at different stages, with a gain in mesenchymal traits in circulating melanoma cells (CTCs) and proliferative features in metastatic tumours. Intriguingly, we found that distinct phenotype populations interact in a cooperative manner, which generates tumours of greater "fitness," supports CTCs and expands organotropic cues in metastases. Fibronectin, expressed in mesenchymal cells, acts as key player in cooperativity and promotes survival of melanocytic cells. Our data reveal an important role for inter-phenotype communications at various stages of disease progression, suggesting these communications could act as therapeutic target.

A Feed-Forward Mechanosignaling Loop Confers Resistance to Therapies Targeting the MAPK Pathway in BRAF-Mutant Melanoma

Aberrant extracellular matrix (ECM) deposition and stiffening is a physical hallmark of several solid cancers and is associated with therapy failure. BRAF-mutant melanomas treated with BRAF and MEK inhibitors almost invariably develop resistance that is frequently associated with transcriptional reprogramming and a de-differentiated cell state. Melanoma cells secrete their own ECM proteins, an event that is promoted by oncogenic BRAF inhibition. Yet, the contribution of cancer cell-derived ECM and tumor mechanics to drug adaptation and therapy resistance remains poorly understood. Here, we show that melanoma cells can adapt to targeted therapies through a mechanosignaling loop involving the autocrine remodeling of a drug-protective ECM. Analyses revealed that therapy-resistant cells associated with a mesenchymal dedifferentiated state displayed elevated responsiveness to collagen stiffening and force-mediated ECM remodeling through activation of actin-dependent mechanosensors Yes-associated protein (YAP) and myocardin-related transcription factor (MRTF). Short-term inhibition of MAPK pathway also induced mechanosignaling associated with deposition and remodeling of an aligned fibrillar matrix. This provided a favored ECM reorganization that promoted tolerance to BRAF inhibition in a YAP- and MRTF-dependent manner. Matrix remodeling and tumor stiffening were also observed in vivo upon exposure of BRAF-mutant melanoma cell lines or patient-derived xenograft models to MAPK pathway inhibition. Importantly, pharmacologic targeting of YAP reversed treatment-induced excessive collagen deposition, leading to enhancement of BRAF inhibitor efficacy. We conclude that MAPK pathway targeting therapies mechanically reprogram melanoma cells to confer a drug-protective matrix environment. Preventing melanoma cell mechanical reprogramming might be a promising therapeutic strategy for patients on targeted therapies. SIGNIFICANCE: These findings reveal a biomechanical adaptation of melanoma cells to oncogenic BRAF pathway inhibition, which fuels a YAP/MRTF-dependent feed-forward loop associated with tumor stiffening, mechanosensing, and therapy resistance. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/10/1927/F1.large.jpg.

Nongenetic Mechanisms of Drug Resistance in Melanoma

Resistance to targeted and immune-based therapies limits cures in patients with metastatic melanoma. A growing number of reports have identified nongenetic primary resistance mechanisms including intrinsic microenvironment- and lineage plasticity–mediated processes serving critical functions in the persistence of disease throughout therapy. There is a temporally shifting spectrum of cellular identities fluidly occupied by therapy-persisting melanoma cells responsible for driving therapeutic resistance and metastasis. The key epigenetic, metabolic, and phenotypic reprogramming events requisite for the manifestation and maintenance of so-called persister melanoma populations remain poorly understood and underscore the need to comprehensively investigate actionable vulnerabilities. Here we attempt to integrate the field's observations on nongenetic mechanisms of drug resistance in melanoma. We postulate that the future design of therapeutic strategies specifically addressing therapy-persisting subpopulations of melanoma will improve the curative potential of therapy for patients with metastatic disease.

Misregulation of ELK1, AP1, and E12 Transcription Factor Networks Is Associated with Melanoma Progression

Melanoma is among the most malignant cutaneous cancers and when metastasized results in dramatically high mortality. Despite advances in high-throughput gene expression profiling in cancer transcriptomic studies, our understanding of mechanisms driving melanoma progression is still limited. We present here an in-depth bioinformatic analysis of the melanoma RNAseq, chromatin immunoprecipitation (ChIP)seq, and single-cell (sc)RNA seq data to understand cancer progression. Specifically, we have performed a consensus network analysis of RNA-seq data from clinically re-grouped melanoma samples to identify gene co-expression networks that are conserved in early (stage 1) and late (stage 4/invasive) stage melanoma. Overlaying the fold-change information on co-expression networks revealed several coordinately up or down-regulated subnetworks that may play a critical role in melanoma progression. Furthermore, by incorporating histone lysine-27 acetylation information and highly expressed genes identified from the single-cell RNA data from melanoma patient samples, we present a comprehensive list of pathways, putative protein-protein interactions (PPIs) and transcription factor (TF) networks that are driving cancer progression. From this analysis, we have identified Elk1, AP1 and E12 TF networks that coordinately change expression in late melanoma when compared to early melanoma, implicating these TFs in melanoma progression. Additionally, the sumoylation-associated interactome is upregulated in invasive melanoma. Together, this bioinformatic analysis potentially implicates a combination of TF networks and PPIs in melanoma progression, which if confirmed in the experimental systems, could be used as targets for drug intervention in melanoma.

Activated Hgf-Met Signaling Cooperates with Oncogenic BRAF to Drive Primary Cutaneous Melanomas and Angiotropic Lung Metastases in Mice

Oncogenic mutations in the BRAF kinase gene represent the most frequent genomic driver in acquired melanocytic nevi and in cutaneous melanomas. It is currently thought that oncogene-induced senescence and cell cycle arrest limit the ability of oncogenic BRAF to promote melanocyte proliferation in benign nevi. The molecular and cellular mechanisms that allow an oncogenic BRAF mutation to fully transform melanocytes into invasively growing melanoma cells that are able to metastasize systemically are only partially understood. In this study, we show in a genetic mouse model that constitutively enhanced Hgf-Met signaling cooperates with oncogenic BRAF to drive tumor development and metastatic spread. Activation of oncogenic BRAF in mice with transgenic Hgf overexpression and an oncogenic CDK4 germline mutation accelerated and increased the development of primary cutaneous melanomas. Primary melanomas showed considerable phenotypic heterogeneity with frequent signs of dedifferentiation. BRAF activation in Hgf-CDK4 mice also increased the number of lung metastases. Melanoma cells showed a pronounced angiotropic growth pattern both at the invasive front in primary tumors and in metastatic lesions of the lung. Taken together, our work supports the notion that activated Hgf-Met signaling and oncogenic BRAF can cooperate in melanoma pathogenesis.

Epigenetic Silencing of CDR1as Drives IGF2BP3-Mediated Melanoma Invasion and Metastasis

Metastasis is the primary cause of death of cancer patients. Dissecting mechanisms governing metastatic spread may uncover important tumor biology and/or yield promising therapeutic insights. Here, we investigated the role of circular RNAs (circRNA) in metastasis, using melanoma as a model aggressive tumor. We identified silencing of cerebellar degeneration-related 1 antisense (CDR1as), a regulator of miR-7, as a hallmark of melanoma progression. CDR1as depletion results from epigenetic silencing of LINC00632, its originating long non-coding RNA (lncRNA) and promotes invasion in vitro and metastasis in vivo through a miR-7-independent, IGF2BP3-mediated mechanism. Moreover, CDR1as levels reflect cellular states associated with distinct therapeutic responses. Our study reveals functional, prognostic, and predictive roles for CDR1as and expose circRNAs as key players in metastasis.

Dissecting Mechanisms of Melanoma Resistance to BRAF and MEK Inhibitors Revealed Genetic and Non-Genetic Patient- and Drug-Specific Alterations and Remarkable Phenotypic Plasticity

The clinical benefit of MAPK pathway inhibition in BRAF-mutant melanoma patients is limited by the development of acquired resistance. Using drug-naïve cell lines derived from tumor specimens, we established a preclinical model of melanoma resistance to vemurafenib or trametinib to provide insight into resistance mechanisms. Dissecting the mechanisms accompanying the development of resistance, we have shown that (i) most of genetic and non-genetic alterations are triggered in a cell line- and/or drug-specific manner; (ii) several changes previously assigned to the development of resistance are induced as the immediate response to the extent measurable at the bulk levels; (iii) reprogramming observed in cross-resistance experiments and growth factor-dependence restricted by the drug presence indicate that phenotypic plasticity of melanoma cells largely contributes to the sustained resistance. Whole-exome sequencing revealed novel genetic alterations, including a frameshift variant of RBMX found exclusively in phospho-AKT^high resistant cell lines. There was no similar pattern of phenotypic alterations among eleven resistant cell lines, including expression/activity of crucial regulators, such as MITF, AXL, SOX, and NGFR, which suggests that patient-to-patient variability is richer and more nuanced than previously described. This diversity should be considered during the development of new strategies to circumvent the acquired resistance to targeted therapies.

BORIS/CTCFL promotes a switch from a proliferative towards an invasive phenotype in melanoma cells

Melanoma is among the most aggressive cancers due to its tendency to metastasize early. Phenotype switching between a proliferative and an invasive state has been suggested as a critical process for metastasis, though the mechanisms that regulate state transitions are complex and remain poorly understood. Brother of Regulator of Imprinted Sites (BORIS), also known as CCCTC binding factor-Like (CTCFL), is a transcriptional modulator that becomes aberrantly expressed in melanoma. Yet, the role of BORIS in melanoma remains elusive. Here, we show that BORIS is involved in melanoma phenotype switching. Genetic modification of BORIS expression in melanoma cells combined with whole-transcriptome analysis indicated that BORIS expression contributes to an invasion-associated transcriptome. In line with these findings, inducible BORIS overexpression in melanoma cells reduced proliferation and increased migration and invasion, demonstrating that the transcriptional switch is accompanied by a phenotypic switch. Mechanistically, we reveal that BORIS binds near the promoter of transforming growth factor-beta 1 (TFGB1), a well-recognized factor involved in the transition towards an invasive state, which coincided with increased expression of TGFB1. Overall, our study indicates a pro-invasive role for BORIS in melanoma via transcriptional reprogramming.

DNA hypermethylation is associated with invasive phenotype of malignant melanoma

Tumor cell invasion is one of the key processes during cancer progression, leading to life-threatening metastatic lesions in melanoma. As methylation of cancer-related genes plays a fundamental role during tumorigenesis and may lead to cellular plasticity which promotes invasion, our aim was to identify novel epigenetic markers on selected invasive melanoma cells. Using Illumina BeadChip assays and Affymetrix Human Gene 1.0 microarrays, we explored the DNA methylation landscape of selected invasive melanoma cells and examined the impact of DNA methylation on gene expression patterns. Our data revealed predominantly hypermethylated genes in the invasive cells affecting the neural crest differentiation pathway and regulation of the actin cytoskeleton. Integrative analysis of the methylation and gene expression profiles resulted in a cohort of hypermethylated genes (IL12RB2, LYPD6B, CHL1, SLC9A3, BAALC, FAM213A, SORCS1, GPR158, FBN1 and ADORA2B) with decreased expression. On the other hand, hypermethylation in the gene body of the EGFR and RBP4 genes was positively correlated with overexpression of the genes. We identified several methylation changes that can have role during melanoma progression, including hypermethylation of the promoter regions of the ARHGAP22 and NAV2 genes that are commonly altered in locally invasive primary melanomas as well as during metastasis. Interestingly, the down-regulation of the methylcytosine dioxygenase TET2 gene, which regulates DNA methylation, was associated with hypermethylated promoter region of the gene. This can probably lead to the observed global hypermethylation pattern of invasive cells and might be one of the key changes during the development of malignant melanoma cells.

Extracellular Vesicles and Epigenetic Modifications Are Hallmarks of Melanoma Progression

Cutaneous melanoma shows a high metastatic potential based on its ability to overcome the immune system's control. The mechanisms activated for these functions vary extremely and are also represented by the production of a number of extracellular vesicles including exosomes. Other vesicles showing a potential role in the melanoma progression include oncosomes and melanosomes and the majority of them mediate tumor processes including angiogenesis, immune regulation, and modifications of the micro-environment. Moreover, a number of epigenetic modifications have been described in melanoma and abundant production of altered microRNAs (mi-RNAs), non-coding RNAs, histones, and abnormal DNA methylation have been associated with different phases of melanoma progression. In addition, exosomes, miRNAs, and other molecular factors have been used as potential biomarkers reflecting disease evolution while others have been suggested to be potential druggable molecules for therapeutic application.

Hippo-YAP/TAZ signalling in organ regeneration and regenerative medicine

The Hippo pathway and its downstream effectors, the transcriptional co-activators Yes-associated protein (YAP) and transcriptional co-activator with PDZ-binding motif (TAZ), regulate organ growth and cell plasticity during animal development and regeneration. Remarkably, experimental activation of YAP/TAZ in the mouse can promote regeneration in organs with poor or compromised regenerative capacity, such as the adult heart and the liver and intestine of old or diseased mice. However, therapeutic YAP/TAZ activation may cause serious side effects. Most notably, YAP/TAZ are hyperactivated in human cancers, and prolonged activation of YAP/TAZ triggers cancer development in mice. Thus, can the power of YAP/TAZ to promote regeneration be harnessed in a safe way? Here, we review the role of Hippo signalling in animal regeneration, examine the promises and risks of YAP/TAZ activation for regenerative medicine and discuss strategies to activate YAP/TAZ for regenerative therapy while minimizing adverse side effects.

The MITF-SOX10 regulated long non-coding RNA DIRC3 is a melanoma tumour suppressor

The MITF and SOX10 transcription factors regulate the expression of genes important for melanoma proliferation, invasion and metastasis. Despite growing evidence of the contribution of long noncoding RNAs (lncRNAs) in cancer, including melanoma, their functions within MITF-SOX10 transcriptional programmes remain poorly investigated. Here we identify 245 candidate melanoma associated lncRNAs whose loci are co-occupied by MITF-SOX10 and that are enriched at active enhancer-like regions. Our work suggests that one of these, Disrupted In Renal Carcinoma 3 (DIRC3), may be a clinically important MITF-SOX10 regulated tumour suppressor. DIRC3 depletion in human melanoma cells leads to increased anchorage-independent growth, a hallmark of malignant transformation, whilst melanoma patients classified by low DIRC3 expression have decreased survival. DIRC3 is a nuclear lncRNA that activates expression of its neighbouring IGFBP5 tumour suppressor through modulating chromatin structure and suppressing SOX10 binding to putative regulatory elements within the DIRC3 locus. In turn, DIRC3 dependent regulation of IGFBP5 impacts the expression of genes involved in cancer associated processes and is needed for DIRC3 control of anchorage-independent growth. Our work indicates that lncRNA components of MITF-SOX10 networks are an important new class of melanoma regulators and candidate therapeutic targets that can act not only as downstream mediators of MITF-SOX10 function but as feedback regulators of MITF-SOX10 activity.

Lineage-Restricted Regulation of SCD and Fatty Acid Saturation by MITF Controls Melanoma Phenotypic Plasticity

Phenotypic and metabolic heterogeneity within tumors is a major barrier to effective cancer therapy. How metabolism is implicated in specific phenotypes and whether lineage-restricted mechanisms control key metabolic vulnerabilities remain poorly understood. In melanoma, downregulation of the lineage addiction oncogene microphthalmia-associated transcription factor (MITF) is a hallmark of the proliferative-to-invasive phenotype switch, although how MITF promotes proliferation and suppresses invasion is poorly defined. Here, we show that MITF is a lineage-restricted activator of the key lipogenic enzyme stearoyl-CoA desaturase (SCD) and that SCD is required for MITF^High melanoma cell proliferation. By contrast MITF^Low cells are insensitive to SCD inhibition. Significantly, the MITF-SCD axis suppresses metastasis, inflammatory signaling, and an ATF4-mediated feedback loop that maintains de-differentiation. Our results reveal that MITF is a lineage-specific regulator of metabolic reprogramming, whereby fatty acid composition is a driver of melanoma phenotype switching, and highlight that cell phenotype dictates the response to drugs targeting lipid metabolism.