A gene expression signature of invasive potential in metastatic melanoma cells

Deciphering pre-existing and induced 3D genome architecture changes involved in constricted melanoma migration

Metastatic cancer cells traverse constricted spaces that exert forces on their nucleus and the genomic contents within. Cancerous tumors are highly heterogeneous and not all cells within them can achieve such a feat. Here, we investigated what initial genome architecture characteristics favor the constricted migratory ability of cancer cells and which arise only after passage through multiple constrictions. We identified a cell surface protein (ITGB4) whose expression correlates with increased initial constricted migration ability in human melanoma A375 cells. Sorting out this subpopulation allowed us to identify cellular and nuclear features that pre-exist and favor migration, as well as alterations that only appear after cells have passed through constrictions. We identified specific genomic regions that experienced altered genome spatial compartment profiles only after constricted migration. Our study reveals 3D genome structure contributions to both selection and induction mechanisms of cell fate change during cancer metastasis.

The classification of melanocytic gene signatures

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

Deciphering Pre-existing and Induced 3D Genome Architecture Changes involved in Constricted Melanoma Migration

Metastatic cancer cells traverse constricted spaces that exert forces on their nucleus and the genomic contents within. Cancerous tumors are highly heterogeneous and not all cells within them can achieve such a feat. Here, we investigated what initial genome architecture characteristics favor the constricted migratory ability of cancer cells and which arise only after passage through multiple constrictions. We identified a cell surface protein (ITGB4) whose expression correlates with increased initial constricted migration ability in human melanoma A375 cells. Sorting out this subpopulation allowed us to identify cellular and nuclear features that pre-exist and favor migration, as well as alterations that only appear after cells have passed through constrictions. We identified specific genomic regions that experienced altered genome spatial compartment profiles only after constricted migration. Our study reveals 3D genome structure contributions to both selection and induction mechanisms of cell fate change during cancer metastasis.

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.

Proteomic Profiling of Advanced Melanoma Patients to Predict Therapeutic Response to Anti-PD-1 Therapy

PURPOSE

Despite high clinical need, there are no biomarkers that accurately predict the response of patients with metastatic melanoma to anti-PD-1 therapy.

EXPERIMENTAL DESIGN

In this multicenter study, we applied protein depletion and enrichment methods prior to various proteomic techniques to analyze a serum discovery cohort (n = 56) and three independent serum validation cohorts (n = 80, n = 12, n = 17). Further validation analyses by literature and survival analysis followed.

RESULTS

We identified several significantly regulated proteins as well as biological processes such as neutrophil degranulation, cell-substrate adhesion, and extracellular matrix organization. Analysis of the three independent serum validation cohorts confirmed the significant differences between responders (R) and nonresponders (NR) observed in the initial discovery cohort. In addition, literature-based validation highlighted 30 markers overlapping with previously published signatures. Survival analysis using the TCGA database showed that overexpression of 17 of the markers we identified correlated with lower overall survival in patients with melanoma.

CONCLUSIONS

Ultimately, this multilayered serum analysis led to a potential marker signature with 10 key markers significantly altered in at least two independent serum cohorts: CRP, LYVE1, SAA2, C1RL, CFHR3, LBP, LDHB, S100A8, S100A9, and SAA1, which will serve as the basis for further investigation. In addition to patient serum, we analyzed primary melanoma tumor cells from NR and found a potential marker signature with four key markers: LAMC1, PXDN, SERPINE1, and VCAN.

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.

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.

BG34-200 Immunotherapy of Advanced Melanoma

High levels of myeloid-derived cells are characteristic of the tumor microenvironment (TME) of advanced melanoma. These cells interact with tumor cells to suppress the development of antitumor immune responses, regulate tumor metastasis, and drive cancer's resistance to virtually all types of therapy. Therefore, methods to disrupt tumor-associated myeloid cell function are actively being sought to find a cure. Our team has recently developed a plant-derived carbohydrate molecule, BG34-200, that modulates tumor-associated myeloid cells by targeting the cell surface receptor CD11b. In this study, we found that BG34-200 IV administration could significantly inhibit tumor growth and improve survival in B16F10 mice with advanced melanoma. Our data supported a model that the entry of BG34-200 into circulating melanoma tumor-associated inflammatory monocytes (TAIMs) could trigger a sequential immune activation: the BG34-200⁺ TAIM subsets migrated to tumor and differentiated into monocyte-derived dendritic cells (mo-DCs); then, the BG34-200⁺ mo-DCs migrated to tumor draining lymph nodes, where they triggered the generation of tumor-antigen-specific T cells. Based upon these results, we combined BG34-200 treatment with adoptive transfer of TdLN-derived T cells to treat advanced melanoma, which significantly improved animal survival and helped tumor-free survivors be resistant to a second tumor-cell challenge. The scientific findings from this study will allow us to develop new technology and apply BG34-200-based immunotherapy to patients with advanced melanoma who have not responded to current standard of care therapies with and without immunotherapy.

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.

Peroxidasin protein expression and enzymatic activity in metastatic melanoma cell lines are associated with invasive potential

Peroxidasin, a heme peroxidase, has been shown to play a role in cancer progression. mRNA expression has been reported to be upregulated in metastatic melanoma cell lines and connected to the invasive phenotype, but little is known about how peroxidasin acts in cancer cells. We have analyzed peroxidasin protein expression and activity in eight metastatic melanoma cell lines using an ELISA developed with an in-house peroxidasin binding protein. RNAseq data analysis confirmed high peroxidasin mRNA expression in the five cell lines classified as invasive and low expression in the three non-invasive cell lines. Protein levels of peroxidasin were higher in the cell lines with an invasive phenotype. Active peroxidasin was secreted to the cell culture medium, where it accumulated over time, and peroxidasin protein levels in the medium were also much higher in invasive than non-invasive cell lines. The only well-established physiological role of peroxidasin is in the formation of a sulfilimine bond, which cross-links collagen IV in basement membranes via catalyzed oxidation of bromide to hypobromous acid. We found that peroxidasin secreted from melanoma cells formed sulfilimine bonds in uncross-linked collagen IV, confirming peroxidasin activity and hypobromous acid formation. Moreover, 3-bromotyrosine, a stable product of hypobromous acid reacting with tyrosine residues, was detected in invasive melanoma cells, substantiating that their expression of peroxidasin generates hypobromous acid, and showing that it does not exclusively react with collagen IV, but also with other biomolecules.

Genetic and Genomic Pathways of Melanoma Development, Invasion and Metastasis

Melanoma is a serious form of skin cancer that accounts for 80% of skin cancer deaths. Recent studies have suggested that melanoma invasiveness is attributed to phenotype switching, which is a reversible type of cell behaviour with similarities to epithelial to mesenchymal transition. Phenotype switching in melanoma is reported to be independent of genetic alterations, whereas changes in gene transcription, and epigenetic alterations have been associated with invasiveness in melanoma cell lines. Here, we review mutational, transcriptional, and epigenomic alterations that contribute to tumour heterogeneity in melanoma, and their potential to drive melanoma invasion and metastasis. We also discuss three models that are hypothesized to contribute towards aspects of tumour heterogeneity and tumour progression in melanoma, namely the clonal evolution model, the cancer stem cell model, and the phenotype switching model. We discuss the merits and disadvantages of each model in explaining tumour heterogeneity in melanoma, as a precursor to invasion and metastasis.

Transcriptomic profiling in canines and humans reveals cancer specific gene modules and biological mechanisms common to both species

Understanding relationships between spontaneous cancer in companion (pet) canines and humans can facilitate biomarker and drug development in both species. Towards this end we developed an experimental-bioinformatic protocol that analyzes canine transcriptomics data in the context of existing human data to evaluate comparative relevance of canine to human cancer. We used this protocol to characterize five canine cancers: melanoma, osteosarcoma, pulmonary carcinoma, B- and T-cell lymphoma, in 60 dogs. We applied an unsupervised, iterative clustering method that yielded five co-expression modules and found that each cancer exhibited a unique module expression profile. We constructed cancer models based on the co-expression modules and used the models to successfully classify the canine data. These canine-derived models also successfully classified human tumors representing the same cancers, indicating shared cancer biology between canines and humans. Annotation of the module genes identified cancer specific pathways relevant to cells-of-origin and tumor biology. For example, annotations associated with melanin production (PMEL, GPNMB, and BACE2), synthesis of bone material (COL5A2, COL6A3, and COL12A1), synthesis of pulmonary surfactant (CTSH, LPCAT1, and NAPSA), ribosomal proteins (RPL8, RPS7, and RPLP0), and epigenetic regulation (EDEM1, PTK2B, and JAK1) were unique to melanoma, osteosarcoma, pulmonary carcinoma, B- and T-cell lymphoma, respectively. In total, 152 biomarker candidates were selected from highly expressing modules for each cancer type. Many of these biomarker candidates are under-explored as drug discovery targets and warrant further study. The demonstrated transferability of classification models from canines to humans enforces the idea that tumor biology, biomarker targets, and associated therapeutics, discovered in canines, may translate to human medicine.

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.

Comparison of Global DNA Methylation Patterns in Human Melanoma Tissues and Their Derivative Cell Lines

Simple Summary

Cancer cell lines are a defined population of cells, originally sourced from tumour tissue, that can be maintained in culture for an extended period of time. They are a critical laboratory-based model, and are frequently used to unravel mechanisms of cancer cell biology. In all cells, gene activity is in part regulated by DNA methylation, the epigenetic process by which methyl groups are added to DNA. In this study, we demonstrate that at a global level, DNA methylation profiles are globally well conserved, but we identify specific sites that are consistently more methylated in tumour-derived cell lines compared to the original tumour tissue. The genes associated with these common differentially methylated regions are involved in important cellular processes and are strongly enriched for epigenetic mechanisms associated with suppression of gene activity. This study provides a valuable resource for identifying false positives due to cell culture and for better interpretation of future cancer epigenetics studies.

Abstract

DNA methylation is a heritable epigenetic mark that is fundamental to mammalian development. Aberrant DNA methylation is an epigenetic hallmark of cancer cells. Cell lines are a critical in vitro model and very widely used to unravel mechanisms of cancer cell biology. However, limited data are available to assess whether DNA methylation patterns in tissues are retained when cell lines are established. Here, we provide the first genome-scale sequencing-based methylation map of metastatic melanoma tumour tissues and their derivative cell lines. We show that DNA methylation profiles are globally conserved in vitro compared to the tumour tissue of origin. However, we identify sites that are consistently hypermethylated in cell lines compared to their tumour tissue of origin. The genes associated with these common differentially methylated regions are involved in cell metabolism, cell cycle and apoptosis and are also strongly enriched for the H3K27me3 histone mark and PRC2 complex-related genes. Our data indicate that although global methylation patterns are similar between tissues and cell lines, there are site-specific epigenomic differences that could potentially impact gene expression. Our work provides a valuable resource for identifying false positives due to cell culture and for better interpretation of cancer epigenetics studies in the future.

Genome-wide DNA methylation and RNA expression differences correlate with invasiveness in melanoma cell lines

Aims & objectives: The aim of this study was to investigate the role of DNA methylation in invasiveness in melanoma cells. Materials & methods: The authors carried out genome-wide transcriptome (RNA sequencing) and reduced representation bisulfite sequencing methylome profiling between noninvasive (n = 4) and invasive melanoma cell lines (n = 5). Results: The integration of differentially expressed genes and differentially methylated fragments (DMFs) identified 12 DMFs (two in AVPI1, one in HMG20B, two in BCL3, one in NTSR1, one in SYNJ2, one in ROBO2 and four in HORMAD2) that overlapped with either differentially expressed genes (eight DMFs and six genes) or cis-targets of lncRNAs (five DMFs associated with cis-targets and four differentially expressed lncRNAs). Conclusions: DNA methylation changes are associated with a number of transcriptional differences observed in noninvasive and invasive phenotypes 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.

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.

Targeting PSMD14 inhibits melanoma growth through SMAD3 stabilization

Although melanoma therapy is improved by novel molecular targeted reagents, including vemurafenib, aberrant proliferation and early metastasis remain obstacles for melanoma; therefore, novel target molecules for melanoma need to be identified. In this study, we focused on deubiquitinating enzymes, which regulate protein stability through ubiquitin–proteasome systems, and identified 26S proteasome non-ATPase regulatory subunit 14 (PSMD14) as a molecule related to melanoma growth using siRNA library screening. Similar to a previous report, PSMD14 knockdown strongly induced p21 expression and inhibited RB phosphorylation in melanoma. After in silico analysis, TGF-β signaling was identified as a negatively correlated gene set with PSMD14 expression. Although TGF-β signaling is also related to the invasive phenotype of melanoma, PSMD14 knockdown suppressed melanoma migration and reduced SLUG expression, suggesting that targeting PSMD14 suppresses both growth and migration. Furthermore, SMAD3 expression increased in nucleus and SMAD3 degradation was delayed after PSMD14 knockdown. Thus, our present study suggests that targeting PSMD14 can inhibit melanoma growth and migration through either SMAD3 accumulation or SLUG reduction, respectively.

Screening and identification of key genes and pathways in metastatic uveal melanoma based on gene expression using bioinformatic analysis

The current study aimed to elucidate the molecular mechanisms and identify the potential key genes and pathways for metastatic uveal melanoma (UM) using bioinformatics analysis.Gene expression microarray data from GSE39717 included 39 primary UM tissue samples and 2 metastatic UM tissue samples. Differentially expressed genes (DEGs) were generated using Gene Expression Omnibus 2R. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the online Database for Annotation, Visualization and Integrated Discovery (DAVID) tool. The web-based STRING tool was adopted to construct a protein--protein interaction (PPI) network. The MCODE tool in Cytoscape was used to generate significant modules of the PPI network.A total of 213 DEGs were identified. GO and KEGG analyses revealed that the upregulated genes were mainly enriched in extracellular matrix organization and blood coagulation cascades, while the downregulated DEGs were mainly related to protein binding, negative regulation of ERK cascade, nucleus and chromatin modification, and lung and renal cell carcinoma. The most significant module was extracted from the PPI network. GO and KEGG enrichment analyses of the module revealed that the genes were mainly enriched in the extracellular region and space organization, blood coagulation process, and PI3K-Akt signaling pathway. Hub genes, including FN1, APOB, F2, SERPINC1, SERPINA1, APOA1, FGG, PROC, ITIH2, VCAN, TFPI, CXCL8, CDH2, and HP, were identified from DEGs. Survival analysis and hierarchical clustering results revealed that most of the hub genes were associated with prognosis and clinical progression.Results of this bioinformatics analysis may provide predictive biomarkers and potential candidate therapeutic targets for individuals with metastatic UM.

Genomic and signalling pathway characterization of the NZM panel of melanoma cell lines: A valuable model for studying the impact of genetic diversity in melanoma

Melanoma is a disease associated with a very high mutation burden and thus the possibility of a diverse range of oncogenic mechanisms that allow it to evade therapeutic interventions and the immune system. Here, we describe the characterization of a panel of 102 cell lines from metastatic melanomas (the NZM lines), including using whole‐exome and RNA sequencing to analyse genetic variants and gene expression changes in a subset of this panel. Lines possessing all major melanoma genotypes were identified, and hierarchical clustering of gene expression profiles revealed four broad subgroups of cell lines. Immunogenotyping identified a range of HLA haplotypes as well as expression of neoantigens and cancer–testis antigens in the lines. Together, these characteristics make the NZM panel a valuable resource for cell‐based, immunological and xenograft studies to better understand the diversity of melanoma biology and the responses of melanoma to therapeutic interventions.

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.

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.

Hypoxic Gene Signature of Primary and Metastatic Melanoma Cell Lines: Focusing on HIF-1β and NDRG-1

Background:

Hypoxia is an important microenvironmental factor significantly affecting tumor proliferation and progression. The importance of hypoxia is, however, not well known in oncogenesis of malignant melanoma.

Aims:

To evaluate the difference of hypoxic gene expression signatures in primary melanoma cell lines and metastatic melanoma cell lines and to find the expression changes of hypoxia-related genes in primary melanoma cell lines at experimental hypoxic conditions.

Study Design:

Cell study.

Methods:

The mRNA expression levels of hypoxia-related genes in primary melanoma cell lines and metastatic melanoma cell lines and at experimental hypoxic conditions in primary melanoma cell lines were evaluated by using real-time polymerase chain reaction. Depending on the experimental data, we focused on two genes/proteins, the hypoxia-inducible factor-1 beta and the N-myc downstream regulated gene-1. The expression levels of the two proteins were investigated by immunohistochemistry methods in 16 primary and metastatic melanomas, 10 intradermal nevi, and a commercial tissue array comprised of 208 cores including 192 primary and metastatic malignant melanomas.

Results:

The real-time polymerase chain reaction study showed that hypoxic gene expression signature was different between metastatic melanoma cell lines and primary melanoma cell lines. Hypoxic experimental conditions significantly affected the hypoxic gene expression signature. In immunohistochemical study, N-myc downstream regulated gene-1 expression was found to be lower in primary cutaneous melanoma compared to in intradermal nevi (p=0.001). In contrast, the cytoplasmic expression of hypoxia-inducible factor-1 beta was higher in primary cutaneous melanoma than in intradermal nevi (p=0.001). We also detected medium/strong significant correlations between the two proteins studied in the study groups.

Conclusion:

Hypoxic response consists of closely related proteins in more complex pathways. These findings will shed light on hypoxic processes in melanoma and unlock a Pandora’s box for development of new therapeutic strategies.

HVEM has a broader expression than PD-L1 and constitutes a negative prognostic marker and potential treatment target for melanoma

HVEM (Herpes Virus Entry Mediator) engagement of BTLA (B and T Lymphocyte Attenuator) triggers inhibitory signals in T cells and could play a role in evading antitumor immunity. Here, HVEM expression levels in melanoma metastases were analyzed by immunohistochemistry, correlated with overall survival (OS) in 116 patients, and validated by TCGA transcriptomic data. Coincident expression of HVEM and its ligand BTLA was studied in tumor cells and tumor-infiltrating lymphocytes (TILs) by flow cytometry (n = 21) and immunofluorescence (n = 5). Candidate genes controlling HVEM expression in melanoma were defined by bioinformatics studies and validated by siRNA gene silencing. We found that in patients with AJCC stage III and IV melanoma, OS was poorer in those with high HVEM expression on melanoma cells, than in those with a low expression, by immunohistochemistry (p = .0160) or TCGA transcriptomics (p = .0282). We showed a coincident expression of HVEM at the surface of melanoma cells and of BTLA on TILs. HVEM was more widely expressed than PD-L1 in melanoma cells. From a mechanistic perspective, in contrast to PDL1, HVEM expression did not correlate with an IFNγ signature but with an aggressive gene signature. Interestingly, this signature contained MITF, a key player in melanoma biology, whose expression correlated strongly with HVEM. Finally, siRNA gene silencing validated MITF control of HVEM expression. In conclusion, HVEM expression seems to be a prognosis marker and targeting this axis by checkpoint-inhibitors may be of interest in metastatic melanoma.

Distinct Transcriptional Programming Drive Response to MAPK Inhibition in BRAF V600-Mutant Melanoma Patient-Derived Xenografts

Inhibitors targeting BRAF and its downstream kinase MEK produce robust response in patients with advanced BRAF ^V600-mutant melanoma. However, the duration and depth of response vary significantly between patients; therefore, predicting response a priori remains a significant challenge. Here, we utilized the Novartis collection of patient-derived xenografts to characterize transcriptional alterations elicited by BRAF and MEK inhibitors in vivo, in an effort to identify mechanisms governing differential response to MAPK inhibition. We show that the expression of an MITF-high, "epithelial-like" transcriptional program is associated with reduced sensitivity and adaptive response to BRAF and MEK inhibitor treatment. On the other hand, xenograft models that express an MAPK-driven "mesenchymal-like" transcriptional program are preferentially sensitive to MAPK inhibition. These gene-expression programs are somewhat similar to the MITF-high and -low phenotypes described in cancer cell lines, but demonstrate an inverse relationship with drug response. This suggests a discrepancy between in vitro and in vivo experimental systems that warrants future investigations. Finally, BRAF ^V600-mutant melanoma relies on either MAPK or alternative pathways for survival under BRAF and MEK inhibition in vivo, which in turn predicts their response to further pathway suppression using a combination of BRAF, MEK, and ERK inhibitors. Our findings highlight the intertumor heterogeneity in BRAF ^V600-mutant melanoma, and the need for precision medicine strategies to target this aggressive cancer.

Phenotype plasticity as enabler of melanoma progression and therapy resistance

Malignant melanoma is notorious for its inter- and intratumour heterogeneity, based on transcriptionally distinct melanoma cell phenotypes. It is thought that these distinct phenotypes are plastic in nature and that their transcriptional reprogramming enables heterogeneous tumours both to undergo different stages of melanoma progression and to adjust to drug exposure during treatment. Recent advances in genomic technologies and the rapidly expanding availability of large gene expression datasets have allowed for a refined definition of the gene signatures that characterize these phenotypes and have revealed that phenotype plasticity plays a major role in the resistance to both targeted therapy and immunotherapy. In this Review we discuss the definition of melanoma phenotypes through particular transcriptional states and reveal the prognostic relevance of the related gene expression signatures. We review how the establishment of phenotypes is controlled and which roles phenotype plasticity plays in melanoma development and therapy. Because phenotype plasticity in melanoma bears a great resemblance to epithelial-mesenchymal transition, the lessons learned from melanoma will also benefit our understanding of other cancer types.

Reprogramming human A375 amelanotic melanoma cells by catalase overexpression: Upregulation of antioxidant genes correlates with regression of melanoma malignancy and with malignant progression when downregulated

Reactive oxygen species (ROS) are implicated in tumor transformation. The antioxidant system (AOS) protects cells from ROS damage. However, it is also hijacked by cancers cells to proliferate within the tumor. Thus, identifying proteins altered by redox imbalance in cancer cells is an attractive prognostic and therapeutic tool. Gene expression microarrays in A375 melanoma cells with different ROS levels after overexpressing catalase were performed. Dissimilar phenotypes by differential compensation to hydrogen peroxide scavenging were generated. The melanotic A375-A7 (A7) upregulated TYRP1, CNTN1 and UCHL1 promoting melanogenesis. The metastatic A375-G10 (G10) downregulated MTSS1 and TIAM1, proteins absent in metastasis. Moreover, differential coexpression of AOS genes (EPHX2, GSTM3, MGST1, MSRA, TXNRD3, MGST3 and GSR) was found in A7 and G10. Their increase in A7 improved its AOS ability and therefore, oxidative stress response, resembling less aggressive tumor cells. Meanwhile, their decrease in G10 revealed a disruption in the AOS and therefore, enhanced its metastatic capacity.

These gene signatures, not only bring new insights into the physiopathology of melanoma, but also could be relevant in clinical prognostic to classify between non aggressive and metastatic melanomas.

PAX3d mRNA over 2.76 copies/µL in the bloodstream predicts cutaneous malignant melanoma relapse

OBJECTIVE

The aim of this study was to evaluate if our molecular algorithm, based on tumor circulating transcripts, may predict relapse risk in cutaneous malignant melanoma (CMM).

RESULTS

The multi-marker panel was able to differentiate patients with CMM from HC with high diagnostic sensitivity and specificity, especially for MITF-m and TGFB2 (91-100%) whose levels decreased during follow-up of recurrence-free patients, and remained stable in the case of relapse. PAX3d higher than 2.76 copies/µL emerged as a promising biomarker [specificity = 75-93% and negative predictive value = 75-98%] to stratify subjects at high risk of CMM recurrence independently of age, gender and AJCC staging [OD = 9.5(3.2-28.0), p < 0.001]. The survival analysis confirmed PAX3d performance in relapse prediction with significant differences in recurrence risk 12 months after the basal time-point (p = 0.008).

MATERIALS AND METHODS

Peripheral blood was collected from 111 CMM patients and from 87 healthy controls (HC) randomly selected. Each specimen was examined by qRT-PCR analysis for the expression of 3 tumor-related transcripts (PAX3d, MITF-m and TGFB2) at diagnosis, and at the following 6 and 12 months during clinical monitoring.

CONCLUSIONS

We demonstrated the usefulness of our molecular algorithm to indirectly detect circulating melanoma cells in blood, along with PAX3d capability to assess patients' progression and relapse prediction.

Genome-wide methylation sequencing of paired primary and metastatic cell lines identifies common DNA methylation changes and a role for EBF3 as a candidate epigenetic driver of melanoma metastasis

Epigenetic alterations are increasingly implicated in metastasis, whereas very few genetic mutations have been identified as authentic drivers of cancer metastasis. Yet, to date, few studies have identified metastasis-related epigenetic drivers, in part because a framework for identifying driver epigenetic changes in metastasis has not been established. Using reduced representation bisulfite sequencing (RRBS), we mapped genome-wide DNA methylation patterns in three cutaneous primary and metastatic melanoma cell line pairs to identify metastasis-related epigenetic drivers. Globally, metastatic melanoma cell lines were hypomethylated compared to the matched primary melanoma cell lines. Using whole genome RRBS we identified 75 shared (10 hyper- and 65 hypomethylated) differentially methylated fragments (DMFs), which were associated with 68 genes showing significant methylation differences. One gene, Early B Cell Factor 3 (EBF3), exhibited promoter hypermethylation in metastatic cell lines, and was validated with bisulfite sequencing and in two publicly available independent melanoma cohorts (n = 40 and 458 melanomas, respectively). We found that hypermethylation of the EBF3 promoter was associated with increased EBF3 mRNA levels in metastatic melanomas and subsequent inhibition of DNA methylation reduced EBF3 expression. RNAi-mediated knockdown of EBF3 mRNA levels decreased proliferation, migration and invasion in primary and metastatic melanoma cell lines. Overall, we have identified numerous epigenetic changes characterising metastatic melanoma cell lines, including EBF3-induced aggressive phenotypic behaviour with elevated EBF3 expression in metastatic melanoma, suggesting that EBF3 promoter hypermethylation may be a candidate epigenetic driver of metastasis.

Phenotype-based variation as a biomarker of sensitivity to molecularly targeted therapy in melanoma

Transcriptomic phenotypes defined for melanoma have been reported to correlate with sensitivity to various drugs. In this study, we aimed to define a minimal signature that could be used to distinguish melanoma sub-types in vitro, and to determine suitable drugs by which these sub-types can be targeted. By using primary melanoma cell lines, as well as commercially available melanoma cell lines, we find that the evaluation of MLANA and INHBA expression is as capable as one based on a combined analysis performed with genes for stemness, EMT and invasion/proliferation, in identifying melanoma subtypes that differ in their sensitivity to molecularly targeted drugs. Using this approach, we find that 75% of melanoma cell lines can be treated with either the MEK inhibitor AZD6244 or the HSP90 inhibitor 17AAG.

The importance of microenvironment: the role of CCL8 in metastasis formation of melanoma

We have attempted to characterize the changes occurring on the host side during the progression of human melanoma. To investigate the role of tumor microenvironment, we set up such an animal model, which was able to isolate the host related factors playing central role in metastasis formation. One of these ‘factors’, CCL12, was consequently selected and its behavior was examined alongside its human homologue (CCL8). In our animal model, metastasis forming primary melanoma in the host exhibited increased level of CCL12 mRNA expression. In clinical samples, when examining the tumor and the host together, the cumulative (tumor and host) CCL8 expression was lower in the group in which human primary melanoma formed lung metastasis compared to non-metastatic primary tumors. We could not detect significant difference in CCL8 receptor (CCR1) expression between the two groups. Increased migration of the examined tumor cell lines was observed when CCL8 was applied as a chemoattractant. The tumor cells and their interactions can be influenced the expression of CCL8 by dermal fibroblasts, as a significant change in the metastatic microenvironment. Furthermore, we examined changes in miRNA profile resulted by CCL8 and miR146a appears to be a promising prognostic marker for following this process.

Fibroblast-induced switching to the mesenchymal-like phenotype and PI3K/mTOR signaling protects melanoma cells from BRAF inhibitors

The knowledge on how tumor-associated stroma influences efficacy of anti-cancer therapy just started to emerge. Here we show that lung fibroblasts reduce melanoma sensitivity to the BRAF inhibitor (BRAFi) vemurafenib only if the two cell types are in close proximity. In the presence of fibroblasts, the adjacent melanoma cells acquire de-differentiated mesenchymal-like phenotype. Upon treatment with BRAFi, such melanoma cells maintain high levels of phospho ribosomal protein S6 (pS6), i.e. active mTOR signaling, which is suppressed in the BRAFi sensitive cells without stromal contacts. Inhibitors of PI3K/mTOR in combination with BRAFi eradicate pS6high cell subpopulations and potentiate anti-cancer effects in melanoma protected by the fibroblasts. mTOR and BRAF co-inhibition also delayed the development of early-stage lung metastases in vivo. In conclusion, we demonstrate that upon influence from fibroblasts, melanoma cells undergo a phenotype switch to the mesenchymal state, which can support PI3K/mTOR signaling. The lost sensitivity to BRAFi in such cells can be overcome by co-targeting PI3K/mTOR. This knowledge could be explored for designing BRAFi combination therapies aiming to eliminate both stroma-protected and non-protected counterparts of metastases.

A melanoma subtype with intrinsic resistance to BRAF inhibition identified by receptor tyrosine kinases gene-driven classification

Dysregulation of receptor tyrosine kinases (RTKs) contributes to several aspects of oncogenesis including drug resistance. In melanoma, distinct RTKs have been involved in BRAF inhibitors (BRAFi) resistance, yet the utility of RTKs expression pattern to identify intrinsically resistant tumors has not been assessed. Transcriptional profiling of RTKs and integration with a previous classification, reveals three robust subtypes in two independent datasets of melanoma cell lines and one cohort of melanoma samples. This classification was validated by Western blot in a panel of patient-derived melanoma cell lines. One of the subtypes identified here for the first time displayed the highest and lowest expression of EGFR and ERBB3, respectively, and included BRAF-mutant tumors all intrinsically resistant to BRAFi PLX4720, as assessed by analysis of the Cancer Cell Line Encyclopedia pharmacogenomic study and by in vitro growth inhibition assays. High levels of EGFR were detected, even before therapy, in tumor cells of one of three melanoma patients unresponsive to BRAFi. Use of different pharmacological inhibitors highlighted the relevance of PI3K/mTOR signaling for growth of this PLX4720-resistant subtype. Our results identify a specific molecular profile of melanomas intrinsically resistant to BRAFi and suggest the PI3K/mTOR pathway as a potential therapeutic target for these tumors.

Metabolic reprogramming supports the invasive phenotype in malignant melanoma

Invasiveness is a hallmark of aggressive cancer like malignant melanoma, and factors involved in acquisition or maintenance of an invasive phenotype are attractive targets for therapy. We investigated melanoma phenotype modulation induced by the metastasis-promoting microenvironmental protein S100A4, focusing on the relationship between enhanced cellular motility, dedifferentiation and metabolic changes. In poorly motile, well-differentiated Melmet 5 cells, S100A4 stimulated migration, invasion and simultaneously down-regulated differentiation genes and modulated expression of metabolism genes. Metabolic studies confirmed suppressed mitochondrial respiration and activated glycolytic flux in the S100A4 stimulated cells, indicating a metabolic switch toward aerobic glycolysis, known as the Warburg effect. Reversal of the glycolytic switch by dichloracetate induced apoptosis and reduced cell growth, particularly in the S100A4 stimulated cells. This implies that cells with stimulated invasiveness get survival benefit from the glycolytic switch and, therefore, become more vulnerable to glycolysis inhibition. In conclusion, our data indicate that transition to the invasive phenotype in melanoma involves dedifferentiation and metabolic reprogramming from mitochondrial oxidation to glycolysis, which facilitates survival of the invasive cancer cells. Therapeutic strategies targeting the metabolic reprogramming may therefore be effective against the invasive phenotype.

Melanoma cells revive an embryonic transcriptional network to dictate phenotypic heterogeneity

Compared to the overwhelming amount of literature describing how epithelial-to-mesenchymal transition (EMT)-inducing transcription factors orchestrate cellular plasticity in embryogenesis and epithelial cells, the functions of these factors in non-epithelial contexts, such as melanoma, are less clear. Melanoma is an aggressive tumor arising from melanocytes, endowed with unique features of cellular plasticity. The reversible phenotype-switching between differentiated and invasive phenotypes is increasingly appreciated as a mechanism accounting for heterogeneity in melanoma and is driven by oncogenic signaling and environmental cues. This phenotypic switch is coupled with an intriguing and somewhat counterintuitive signaling switch of EMT-inducing transcription factors. In contrast to carcinomas, different EMT-inducing transcription factors have antagonizing effects in melanoma. Balancing between these different EMT transcription factors is likely the key to successful metastatic spread of melanoma.

The neural crest: a versatile organ system

The neural crest is the name given to the strip of cells at the junction between neural and epidermal ectoderm in neurula-stage vertebrate embryos, which is later brought to the dorsal neural tube as the neural folds elevate. The neural crest is a heterogeneous and multipotent progenitor cell population whose cells undergo EMT then extensively and accurately migrate throughout the embryo. Neural crest cells contribute to nearly every organ system in the body, with derivatives of neuronal, glial, neuroendocrine, pigment, and also mesodermal lineages. This breadth of developmental capacity has led to the neural crest being termed the fourth germ layer. The neural crest has occupied a prominent place in developmental biology, due to its exaggerated migratory morphogenesis and its remarkably wide developmental potential. As such, neural crest cells have become an attractive model for developmental biologists for studying these processes. Problems in neural crest development cause a number of human syndromes and birth defects known collectively as neurocristopathies; these include Treacher Collins syndrome, Hirschsprung disease, and 22q11.2 deletion syndromes. Tumors in the neural crest lineage are also of clinical importance, including the aggressive melanoma and neuroblastoma types. These clinical aspects have drawn attention to the selection or creation of neural crest progenitor cells, particularly of human origin, for studying pathologies of the neural crest at the cellular level, and also for possible cell therapeutics. The versatility of the neural crest lends itself to interlinked research, spanning basic developmental biology, birth defect research, oncology, and stem/progenitor cell biology and therapy.

The motile breast cancer phenotype roles of proteoglycans/glycosaminoglycans

The consecutive stages of cancer growth and dissemination are obligatorily perpetrated through specific interactions of the tumor cells with their microenvironment. Importantly, cell-associated and tumor microenvironment glycosaminoglycans (GAGs)/proteoglycan (PG) content and distribution are markedly altered during tumor pathogenesis and progression. GAGs and PGs perform multiple functions in specific stages of the metastatic cascade due to their defined structure and ability to interact with both ligands and receptors regulating cancer pathogenesis. Thus, GAGs/PGs may modulate downstream signaling of key cellular mediators including insulin growth factor receptor (IGFR), epidermal growth factor receptor (EGFR), estrogen receptors (ERs), or Wnt members. In the present review we will focus on breast cancer motility in correlation with their GAG/PG content and critically discuss mechanisms involved. Furthermore, new approaches involving GAGs/PGs as potential prognostic/diagnostic markers or as therapeutic agents for cancer-related pathologies are being proposed.

Distinct microRNA expression signatures are associated with melanoma subtypes and are regulated by HIF1A

The complex genetic changes underlying metastatic melanoma need to be deciphered to develop new and effective therapeutics. Previously, genome-wide microarray analyses of human melanoma identified two reciprocal gene expression programs, including transcripts regulated by either transforming growth factor, beta 1 (TGFβ1) pathways, or microphthalmia-associated transcription factor (MITF)/SRY-box containing gene 10 (SOX10) pathways. We extended this knowledge by discovering that melanoma cell lines with these two expression programs exhibit distinctive microRNA (miRNA) expression patterns. We also demonstrated that hypoxia-inducible factor 1 alpha (HIF1A) is increased in TGFβ1 pathway-expressing melanoma cells and that HIF1A upregulates miR-210, miR-218, miR-224, and miR-452. Reduced expression of these four miRNAs in TGFβ1 pathway-expressing melanoma cells arrests the cell cycle, while their overexpression in mouse melanoma cells increases the expression of the hypoxic response gene Bnip3. Taken together, these data suggest that HIF1A may regulate some of the gene expression and biological behavior of TGFβ1 pathway-expressing melanoma cells, in part via alterations in these four miRNAs.

Retrotransposon hypomethylation in melanoma and expression of a placenta-specific gene

In the human placenta, DNA hypomethylation permits the expression of retrotransposon-derived genes that are normally silenced by methylation in somatic tissues. We previously identified hypomethylation of a retrotransposon-derived transcript of the voltage-gated potassium channel gene KCNH5 that is expressed only in human placenta. However, an RNA sequence from this placental-specific transcript has been reported in melanoma. This study examined the promoter methylation and expression of the retrotransposon-derived KCNH5 transcript in 25 melanoma cell lines to determine whether the acquisition of 'placental' epigenetic marks is a feature of melanoma. Methylation and gene expression analysis revealed hypomethylation of this retrotransposon in melanoma cell lines, particularly in those samples that express the placental KCNH5 transcript. Therefore we propose that hypomethylation of the placental-specific KCNH5 promoter is frequently associated with KCNH5 expression in melanoma cells. Our findings show that melanoma can develop hypomethylation of a retrotransposon-derived gene; a characteristic notably shared with the normal placenta.

Isolation and molecular characterization of circulating melanoma cells

Melanoma is an invasive malignancy with a high frequency of blood-borne metastases, but circulating tumor cells (CTCs) have not been readily isolated. We adapted microfluidic CTC capture to a tamoxifen-driven B-RAF/PTEN mouse melanoma model. CTCs were detected in all tumor-bearing mice and rapidly declined after B-RAF inhibitor treatment. CTCs were shed early from localized tumors, and a short course of B-RAF inhibition following surgical resection was sufficient to dramatically suppress distant metastases. The large number of CTCs in melanoma-bearing mice enabled a comparison of RNA-sequencing profiles with matched primary tumors. A mouse melanoma CTC-derived signature correlated with invasiveness and cellular motility in human melanoma. CTCs were detected in smaller numbers in patients with metastatic melanoma and declined with successful B-RAF-targeted therapy. Together, the capture and molecular characterization of CTCs provide insight into the hematogenous spread of melanoma.

Gene expression profiling identifies microphthalmia-associated transcription factor (MITF) and Dickkopf-1 (DKK1) as regulators of microenvironment-driven alterations in melanoma phenotype

BACKGROUND

The diversity of functional phenotypes observed within a tumor does not exclusively result from intratumoral genetic heterogeneity but also from the response of cancer cells to the microenvironment. We have previously demonstrated that the morphological and functional phenotypes of melanoma can be dynamically altered upon external stimuli.

FINDINGS

In the present study, transcriptome profiles were generated to explore the molecules governing phenotypes of melanospheres grown in the bFGF(+)EGF(+) serum-free cultures and monolayers maintained in the serum-containing medium. Higher expression levels of MITF-dependent genes that are responsible for differentiation, e.g., TYR and MLANA, and stemness-related genes, e.g., ALDH1A1, were detected in melanospheres. These results were supported by the observation that the melanospheres contained more pigmented cells and cells exerting the self-renewal capacity than the monolayers. In addition, the expression of the anti-apoptotic, MITF-dependent genes e.g., BCL2A1 was also higher in the melanospheres. The enhanced activity of MITF in melanospheres, as illustrated by the increased expression of 74 MITF-dependent genes, identified MITF as a central transcriptional regulator in melanospheres. Importantly, several genes including MITF-dependent ones were expressed in melanospheres and original tumors at similar levels. The reduced MITF level in monolayers might be partially explained by suppression of the Wnt/β-catenin pathway, and DKK1, a secreted inhibitor of this pathway, was highly up-regulated in monolayers in comparison to melanospheres and original tumors. Furthermore, the silencing of DKK1 in monolayers increased the percentage of cells with self-renewing capacity.

CONCLUSIONS

Our study indicates that melanospheres can be used to unravel the molecular pathways that sustain intratumoral phenotypic heterogeneity. Melanospheres directly derived from tumor specimens more accurately mirrored the morphology and gene expression profiles of the original tumors compared to monolayers. Therefore, melanospheres represent a relevant preclinical tool to study new anticancer treatment strategies.

MITF and PAX3 Play Distinct Roles in Melanoma Cell Migration; Outline of a "Genetic Switch" Theory Involving MITF and PAX3 in Proliferative and Invasive Phenotypes of Melanoma

Melanoma is a very aggressive neoplasm with a propensity to undergo progression and invasion early in its evolution. The molecular pathways underpinning invasion in melanoma are now just beginning to be elucidated, but a clear understanding of the transition from non-invasive to invasive melanoma cells remains elusive. Microphthalmia-associated transcription factor (MITF), is thought to be a central player in melanoma biology, and it controls many aspects of the phenotypic expression of the melanocytic lineage. However, recently the paired box transcription factor PAX3 was shown to transcriptionally activate POU3F2/BRN2, leading to direct repression of MITF expression. Here we present a theory to explain melanoma phenotype switching and discuss the predictions that this theory makes. One prediction is that independent and opposing roles for MITF and PAX3 in melanoma would be expected, and we present empirical evidence supporting this: in melanoma tissues PAX3 expression occurs independently of MITF, and PAX3 does not play a key role in melanoma cell proliferation. Furthermore, we show that knockdown of PAX3 inhibits cell migration in a group of "lower MITF" melanoma cell lines, while knockdown of MITF promotes cell migration in a complementary "higher MITF" group of melanoma cell lines. Moreover, the morphological effects of knocking down PAX3 versus MITF in melanoma cells were found to differ. While these data support the notion of independent roles for MITF and PAX3, additional experiments are required to provide robust examination of the proposed genetic switch theory. Only upon clear delineation of the mechanisms associated with progression and invasion of melanoma cells will successful treatments for invasive melanoma be developed.

Human melanoma metastasis in NSG mice correlates with clinical outcome in patients

Studies of human cancer metastasis have been limited by a lack of experimental assays in which cancer cells from patients metastasize in vivo in a way that correlates with clinical outcome. This makes it impossible to study intrinsic differences in the metastatic properties of cancers from different patients. We recently developed an assay in which human melanomas readily engraft in nonobese diabetic/severe combined immunodeficient interleukin-2 receptor-γ chain null (NSG) mice. We show that melanomas from 25 patients exhibited reproducible differences in the rate of spontaneous metastasis after transplantation into NSG mice and that these differences correlated with clinical outcome in the patients. Stage IIIB/C melanomas that formed distant metastases within 22 months in patients also formed tumors that metastasized widely in NSG mice, whereas stage IIIB/C melanomas that did not form distant metastases within 22 to 50 months in patients metastasized more slowly in NSG mice. These differences in the efficiency of metastasis correlated with the presence of circulating melanoma cells in the blood of NSG mice, suggesting that the rate of entry into the blood is one factor that limits the rate of metastasis. The study of NSG mice can therefore yield information about the metastasis of human melanomas in vivo, in this case revealing intrinsic differences among stage III melanomas in their ability to circulate/survive in the blood and to metastasize.

Efficient gene expression profiling of laser-microdissected melanoma metastases

Comparing the transcriptomes of primary and metastatic tumour tissues is a useful strategy for studying tumour progression. One factor limiting the interpretation of tissue-based transcriptomic data is the lack of cell-type purity. Laser capture microdissection (LCM) has been shown to be useful for overcoming this limitation. We established an efficient protocol for gene expression profiling of LCM and matched metastatic melanomas using a transgenic mouse model. This optimized workflow combines microsurgical recovery of mouse lungs, appropriate tissue freezing, laser microdissection of homogeneous tumour cell populations from cryosections, isolation of high-quality RNA and gene expression analysis. The RNA isolated from laser-microdissected material was not contaminated by stroma cells, was of excellent quality, and the synthesis of cDNAs was homogeneous and highly reproducible. Subsequent custom-based Taqman-low-density-array (TLDA)-based gene expression profiling identified stronger expression of five genes (M-MITF, TYR, STAT3, CCND1 and PAX3) in primary than metastatic melanoma. We detected only minor transcriptomic differences between primary and metastatic melanoma tissue. This optimized workflow could be very valuable for various studies requiring cell type-specific transcriptomic analysis.

Variable Expression of GLIPR1 Correlates with Invasive Potential in Melanoma Cells

GLI pathogenesis-related 1 (GLIPR1) was previously identified as an epigenetically regulated tumor suppressor in prostate cancer and, conversely, an oncoprotein in glioma. More recently, GLIPR1 was shown to be differentially expressed in other cancers including ovarian, acute myeloid leukemia, and Wilms' tumor. Here we investigated GLIPR1 expression in metastatic melanoma cell lines and tissue. GLIPR1 was variably expressed in metastatic melanoma cells, and transcript levels correlated with degree of GLIPR1 promoter methylation in vitro. Elevated GLIPR1 levels were correlated with increased invasive potential, and siRNA-mediated knockdown of GLIPR1 expression resulted in reduced cell migration and proliferation in vitro. Immunohistochemical studies of melanoma tissue microarrays showed moderate to high staining for GLIPR1 in 50% of specimens analyzed. GLIPR1 staining was observed in normal skin in merocrine sweat glands, sebaceous glands, and hair follicles within the dermis.

Heterogeneity of expression of epithelial-mesenchymal transition markers in melanocytes and melanoma cell lines

The epithelial–mesenchymal transition (EMT) describes a reversible switch from an epithelial-like to a mesenchymal-like phenotype. It is essential for the development of the normal epithelium and also contributes to the invasive properties of carcinomas. At the molecular level, the EMT transition is characterized by a series of coordinated changes including downregulation of the junctional protein E-cadherin (CDH1), up-regulation of transcriptional repressors of E-cadherin such as Snail (SNAI1) and Slug (SNAI2), and up-regulation of N-cadherin. We wished to determine whether cultured normal melanocytes and melanoma cell lines, which are derived from the neural crest, showed signs of a similarly coordinated phenotypic switch. We investigated normal melanocytes and 25 cell lines derived from New Zealand patients with metastatic melanoma. Most lines had been previously genotyped for common mutations such as BRAF, NRAS, PIK3CA (phosphatidylinositol-3-kinase), TP53 (p53), and CDKN2A (p16). Expression of E-cadherin, N-cadherin, microphthalmia-associated transcription factor (MITF), Snail, Slug, Axl, p53, and Hdm2 was compared by western blotting. Normal melanocytes expressed each of these proteins except for Snail, while normal melanocytes and almost every melanoma line expressed Slug. Expression of individual markers among different melanoma lines varied from high to low or undetectable. Quantitation of western blots showed that expression of MITF-M, the melanocyte-specific isoform of MITF, was positively related to that of E-cadherin but inversely related to that of N-cadherin and Axl. There was also no apparent relationship between expression of any particular marker and the presence of BRAF, NRAS, PIK3CA, TP53, or CDKN2A mutations. The results suggest that melanomas do not show the classical epithelial and mesenchymal phenotypes but rather display either high E-cadherin/high MITF-M expression on one hand, or high N-cadherin/high Axl expression on the other. These may correspond to differentiated and invasive phenotypes in vivo.

Gene therapy for advanced melanoma: selective targeting and therapeutic nucleic acids

Despite recent advances, the treatment of malignant melanoma still results in the relapse of the disease, and second line treatment mostly fails due to the occurrence of resistance. A wide range of mutations are known to prevent effective treatment with chemotherapeutic drugs. Hence, approaches with biopharmaceuticals including proteins, like antibodies or cytokines, are applied. As an alternative, regimens with therapeutically active nucleic acids offer the possibility for highly selective cancer treatment whilst avoiding unwanted and toxic side effects. This paper gives a brief introduction into the mechanism of this devastating disease, discusses the shortcoming of current therapy approaches, and pinpoints anchor points which could be harnessed for therapeutic intervention with nucleic acids. We bring the delivery of nucleic acid nanopharmaceutics into perspective as a novel antimelanoma therapeutic approach and discuss the possibilities for melanoma specific targeting. The latest reports on preclinical and already clinical application of nucleic acids in melanoma are discussed.

HOXA1 drives melanoma tumor growth and metastasis and elicits an invasion gene expression signature that prognosticates clinical outcome

Metastatic melanoma is a highly lethal disease notorious for its aggressive clinical course and eventual resistance to existing therapies. Currently we possess a limited understanding of the genetic events driving melanoma progression, and much effort is focused on identifying pro-metastatic aberrations or perturbed signaling networks that constitute new therapeutic targets. In this study, we validate and assess the mechanism by which homeobox transcription factor A1 (HOXA1), a pro-invasion oncogene previously identified in a metastasis screen by our group, contributes to melanoma progression. Transcriptome and pathway profiling analyses of cells expressing HOXA1 reveals up-regulation of factors involved in diverse cytokine pathways that include the TGFβ signaling axis, which we further demonstrate to be required for HOXA1-mediated cell invasion in melanoma cells. Transcriptome profiling also shows HOXA1’s ability to potently down-regulate expression of microphthalmia-associated transcription factor (MITF) and other genes required for melanocyte differentiation, suggesting a mechanism by which HOXA1 expression de-differentiates cells into a pro-invasive cell state concomitant with TGFβ activation. Our analysis of publicly available datasets indicate that the HOXA1-induced gene signature successfully categorizes melanoma specimens based on their metastatic potential and, importantly, is capable of stratifying melanoma patient risk for metastasis based on expression in primary tumors. Together, these validation data and mechanistic insights suggest that patients whose primary tumors express HOXA1 are among a high-risk metastasis subgroup that should be considered for anti-TGFβ therapy in adjuvant settings. Moreover, further analysis of HOXA1 target genes in melanoma may reveal new pathways or targets amenable to therapeutic intervention.

Integrative genomics identifies gene signature associated with melanoma ulceration

Background

Despite the extensive research approaches applied to characterise malignant melanoma, no specific molecular markers are available that are clearly related to the progression of this disease. In this study, our aims were to define a gene expression signature associated with the clinical outcome of melanoma patients and to provide an integrative interpretation of the gene expression -, copy number alterations -, and promoter methylation patterns that contribute to clinically relevant molecular functional alterations.

Methods

Gene expression profiles were determined using the Affymetrix U133 Plus2.0 array. The NimbleGen Human CGH Whole-Genome Tiling array was used to define CNAs, and the Illumina GoldenGate Methylation platform was applied to characterise the methylation patterns of overlapping genes.

Results

We identified two subclasses of primary melanoma: one representing patients with better prognoses and the other being characteristic of patients with unfavourable outcomes. We assigned 1,080 genes as being significantly correlated with ulceration, 987 genes were downregulated and significantly enriched in the p53, Nf-kappaB, and WNT/beta-catenin pathways. Through integrated genome analysis, we defined 150 downregulated genes whose expression correlated with copy number losses in ulcerated samples. These genes were significantly enriched on chromosome 6q and 10q, which contained a total of 36 genes. Ten of these genes were downregulated and involved in cell-cell and cell-matrix adhesion or apoptosis. The expression and methylation patterns of additional genes exhibited an inverse correlation, suggesting that transcriptional silencing of these genes is driven by epigenetic events.

Conclusion

Using an integrative genomic approach, we were able to identify functionally relevant molecular hotspots characterised by copy number losses and promoter hypermethylation in distinct molecular subtypes of melanoma that contribute to specific transcriptomic silencing and might indicate a poor clinical outcome of melanoma.