Faculty Bibliography

As efforts to study the mechanisms of melanoma metastasis and novel therapeutic approaches multiply, researchers need accurate, high-throughput methods to evaluate the effects on tumor burden resulting from specific interventions. We show that automated quantification of tumor content from whole slide images is a compelling solution to assess in vivo experiments. In order to increase the outflow of data collection from preclinical studies, we assembled a large dataset with annotations and trained a deep neural network for the quantitative analysis of melanoma tumor content on histopathological sections of murine models. After assessing its performance in segmenting these images, the tool obtained consistent results with an orthogonal method (bioluminescence) of measuring metastasis in an experimental setting. This AI-based algorithm, made freely available to academic laboratories through a web-interface called MetFinder, promises to become an asset for melanoma researchers and pathologists interested in accurate, quantitative assessment of metastasis burden.

PURPOSE/OBJECTIVE:Patients with stage II and III cutaneous primary melanoma vary considerably in their risk of melanoma-related death. We explore the ability of methylation profiling to distinguish primary melanoma methylation classes and their associations with clinicopathologic characteristics and survival. MATERIALS AND METHODS/METHODS:InterMEL is a retrospective case-control study that assembled primary cutaneous melanomas from American Joint Committee on Cancer (AJCC) 8th edition stage II and III patients diagnosed between 1998 and 2015 in the United States and Australia. Cases are patients who died of melanoma within 5 years from original diagnosis. Controls survived longer than 5 years without evidence of melanoma recurrence or relapse. Methylation classes, distinguished by consensus clustering of 850K methylation data, were evaluated for their clinicopathologic characteristics, 5-year survival status, and differentially methylated gene sets. RESULTS:= .007). CIMP and IM had a 2-fold higher likelihood of 5-year death from melanoma than LM (CIMP odds ratio [OR], 2.16 [95% CI, 1.18 to 3.96]; IM OR, 2.00 [95% CI, 1.12 to 3.58]) in a multivariable model adjusted for age, sex, log Breslow thickness, ulceration, mitotic index, and N stage. Despite more extensive CpG island hypermethylation in CIMP, CIMP and IM shared similar patterns of differential methylation and gene set enrichment compared with LM. CONCLUSION/CONCLUSIONS:Melanoma MethylClasses may provide clinical value in predicting 5-year death from melanoma among patients with primary melanoma independent of other clinicopathologic factors.

Melanoma represents a critical clinical challenge owing to its unfavorable outcomes. This type of skin cancer exhibits unique adaptability to the brain microenvironment, but its underlying molecular mechanisms are poorly understood. Recent findings have suggested that melanoma brain metastases may share biological processes similar to those found in various neurodegenerative diseases. To further characterize melanoma brain metastasis development, we explore the relationship between the transcriptional profiles of melanoma brain metastases and the neurodegenerative diseases Alzheimer's disease, Parkinson's disease, and multiple sclerosis. We take an in silico approach to unveil a neurodegenerative signature of melanoma brain metastases compared with those of melanoma nonbrain metastasis (53 dysregulated genes were enriched in 11 functional terms, such as associated terms to the extracellular matrix and development) and with those of nontumor-bearing brain controls (195 dysregulated genes, mostly involved in development and cell differentiation, chromatin remodeling and nucleosome organization, and translation). Two genes, ITGA10 and DNAJC6, emerged as key potential markers being dysregulated in both scenarios. Finally, we developed an open-source, user-friendly web tool (https://bioinfo.cipf.es/metafun-mbm/) that allows interactive exploration of the complete results.

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

Glycosylation is a hallmark of cancer biology, and altered glycosylation influences multiple facets of melanoma growth and progression. To identify glycosyltransferases, glycans, and glycoproteins essential for melanoma maintenance, we conducted an in vivo growth screen with a pooled shRNA library of glycosyltransferases, lectin microarray profiling of benign nevi and melanoma patient samples, and mass spectrometry-based glycoproteomics. We found that α-2,3 sialyltransferases ST3GAL1 and ST3GAL2 and corresponding α-2,3-linked sialosides are upregulated in melanoma compared to nevi and are essential for melanoma growth in vivo and in vitro. Glycoproteomics revealed that glycoprotein targets of ST3GAL1 and ST3GAL2 are enriched in transmembrane proteins involved in growth signaling, including the amino acid transporter Solute Carrier Family 3 Member 2 (SLC3A2/CD98hc). CD98hc suppression mimicked the effect of ST3GAL1 and ST3GAL2 silencing, inhibiting melanoma cell proliferation. We found that both CD98hc protein stability and its pro-survival effect in melanoma are dependent upon α-2,3 sialylation mediated by ST3GAL1 and ST3GAL2. In summary, our studies reveal that α-2,3-sialosides functionally contribute to melanoma maintenance, supporting ST3GAL1 and ST3GAL2 as novel therapeutic targets in these tumors.

Super-enhancers (SEs) are stretches of enhancers ensuring a high level of expression of key genes associated with cell function. The identification of cancer-specific SE-driven genes is a powerful means for the development of innovative therapeutic strategies. Here, we identify a MITF/SOX10/TFIIH-dependent SE promoting the expression of BAHCC1 in a broad panel of melanoma cells. BAHCC1 is highly expressed in metastatic melanoma and is required for tumor engraftment, growth, and dissemination. Integrative genomics analyses reveal that BAHCC1 is a transcriptional regulator controlling expression of E2F/KLF-dependent cell-cycle and DNA-repair genes. BAHCC1 associates with BRG1-containing remodeling complexes at the promoters of these genes. BAHCC1 silencing leads to decreased cell proliferation and delayed DNA repair. Consequently, BAHCC1 deficiency cooperates with PARP inhibition to induce melanoma cell death. Our study identifies BAHCC1 as an SE-driven gene expressed in melanoma and demonstrates how its inhibition can be exploited as a therapeutic target.

Differential methylation plays an important role in melanoma development and is associated with survival, progression and response to treatment. However, the mechanisms by which methylation promotes melanoma development are poorly understood. The traditional explanation of selective advantage provided by differential methylation postulates that hypermethylation of regulatory 5'-cytosine-phosphate-guanine-3' dinucleotides (CpGs) downregulates the expression of tumor suppressor genes and therefore promotes tumorigenesis. We believe that other (not necessarily alternative) explanations of the selective advantages of methylation are also possible. Here, we hypothesize that melanoma cells use methylation to shut down transcription of nonessential genes - those not required for cell survival and proliferation. Suppression of nonessential genes allows tumor cells to be more efficient in terms of energy and resource usage, providing them with a selective advantage over the tumor cells that transcribe and subsequently translate genes they do not need. We named the hypothesis the Rule Out (RO) hypothesis. The RO hypothesis predicts higher methylation of CpGs located in regulatory regions (CpG islands) of nonessential genes. It also predicts the higher methylation of regulatory CpGs linked to nonessential genes in melanomas compared to nevi and lower expression of nonessential genes in malignant (derived from melanoma) versus normal (derived from nonaffected skin) melanocytes. The analyses conducted using in-house and publicly available data found that all predictions derived from the RO hypothesis hold, providing observational support for the hypothesis.

Metastatic melanoma develops once transformed melanocytic cells begin to de-differentiate into migratory and invasive melanoma cells with neural crest cell (NCC)-like and epithelial-to-mesenchymal transition (EMT)-like features. However, it is still unclear how transformed melanocytes assume a metastatic melanoma cell state. Here, we define DNA methylation changes that accompany metastatic progression in melanoma patients and discover Nuclear Receptor Subfamily 2 Group F, Member 2 - isoform 2 (NR2F2-Iso2) as an epigenetically regulated metastasis driver. NR2F2-Iso2 is transcribed from an alternative transcriptional start site (TSS) and it is truncated at the N-terminal end which encodes the NR2F2 DNA-binding domain. We find that NR2F2-Iso2 expression is turned off by DNA methylation when NCCs differentiate into melanocytes. Conversely, this process is reversed during metastatic melanoma progression, when NR2F2-Iso2 becomes increasingly hypomethylated and re-expressed. Our functional and molecular studies suggest that NR2F2-Iso2 drives metastatic melanoma progression by modulating the activity of full-length NR2F2 (Isoform 1) over EMT- and NCC-associated target genes. Our findings indicate that DNA methylation changes play a crucial role during metastatic melanoma progression, and their control of NR2F2 activity allows transformed melanocytes to acquire NCC-like and EMT-like features. This epigenetically regulated transcriptional plasticity facilitates cell state transitions and metastatic spread.

INTRODUCTION:We are conducting a multicenter study to identify classifiers predictive of disease-specific survival in patients with primary melanomas. Here we delineate the unique aspects, challenges, and best practices for optimizing a study of generally small-sized pigmented tumor samples including primary melanomas of at least 1.05mm from AJTCC TNM stage IIA-IIID patients. We also evaluated tissue-derived predictors of extracted nucleic acids' quality and success in downstream testing. This ongoing study will target 1,000 melanomas within the international InterMEL consortium. METHODS:Following a pre-established protocol, participating centers ship formalin-fixed paraffin embedded (FFPE) tissue sections to Memorial Sloan Kettering Cancer Center for the centralized handling, dermatopathology review and histology-guided coextraction of RNA and DNA. Samples are distributed for evaluation of somatic mutations using next gen sequencing (NGS) with the MSK-IMPACTTM assay, methylation-profiling (Infinium MethylationEPIC arrays), and miRNA expression (Nanostring nCounter Human v3 miRNA Expression Assay). RESULTS:Sufficient material was obtained for screening of miRNA expression in 683/685 (99%) eligible melanomas, methylation in 467 (68%), and somatic mutations in 560 (82%). In 446/685 (65%) cases, aliquots of RNA/DNA were sufficient for testing with all three platforms. Among samples evaluated by the time of this analysis, the mean NGS coverage was 249x, 59 (18.6%) samples had coverage below 100x, and 41/414 (10%) failed methylation QC due to low intensity probes or insufficient Meta-Mixed Interquartile (BMIQ)- and single sample (ss)- Noob normalizations. Six of 683 RNAs (1%) failed Nanostring QC due to the low proportion of probes above the minimum threshold. Age of the FFPE tissue blocks (p<0.001) and time elapsed from sectioning to co-extraction (p = 0.002) were associated with methylation screening failures. Melanin reduced the ability to amplify fragments of 200bp or greater (absent/lightly pigmented vs heavily pigmented, p<0.003). Conversely, heavily pigmented tumors rendered greater amounts of RNA (p<0.001), and of RNA above 200 nucleotides (p<0.001). CONCLUSION:Our experience with many archival tissues demonstrates that with careful management of tissue processing and quality control it is possible to conduct multi-omic studies in a complex multi-institutional setting for investigations involving minute quantities of FFPE tumors, as in studies of early-stage melanoma. The study describes, for the first time, the optimal strategy for obtaining archival and limited tumor tissue, the characteristics of the nucleic acids co-extracted from a unique cell lysate, and success rate in downstream applications. In addition, our findings provide an estimate of the anticipated attrition that will guide other large multicenter research and consortia.

Leiomyosarcoma (LMS) is an aggressive, often poorly differentiated cancer of the smooth muscle (SM) lineage for which the molecular drivers of transformation and progression are poorly understood. In microRNA (miRNA) profiling studies, miR-130b was previously found to be upregulated in LMS vs. normal SM, and down-regulated during the differentiation of mesenchymal stem cells (MSCs) into SM, suggesting a role in LMS tumor progression. In the present study, the effects of miR-130b on human LMS tumorigenesis were investigated. Stable miR-130b overexpression enhanced invasion of LMS cells in vitro, and led to the formation of undifferentiated, pleomorphic tumors in vivo, with increased growth and metastatic potential compared to control LMS cells. TSC1 was identified as a direct miR-130b target in luciferase-3'UTR assays, and shRNA-mediated knockdown of TSC1 replicated miR-130b effects. Loss-of-function and gain-of-function studies showed that miR-130b levels regulate cell morphology and motility. Following miR-130b suppression, LMS cells adopted a rounded morphology, amoeboid mode of cell movement and enhanced invasive capacity that was Rho/ROCK dependent. Conversely, miR-130b-overexpressing LMS cells exhibited Rho-independent invasion, accompanied by down-regulation of Rho-pathway effectors. In mesenchymal stem cells, both miR-130b overexpression and TSC1 silencing independently impaired SM differentiation in vitro. Together, the data reveal miR-130b as a pro-oncogenic miRNA in LMS and support a miR-130b-TSC1 regulatory network that enhances tumor progression via inhibition of SM differentiation.