Temporal genomic analysis of melanoma rejection identifies regulators of tumor immune evasion

Decreased intra-tumor heterogeneity (ITH) correlates with increased patient survival and immunotherapy response. However, even highly homogenous tumors may display variability in their aggressiveness, and how immunologic-factors impinge on their aggressiveness remains understudied. Here we studied the mechanisms responsible for the immune-escape of murine tumors with low ITH. We compared the temporal growth of homogeneous, genetically-similar single-cell clones that are rejected vs. those that are not-rejected after transplantation in-vivo using single-cell RNA sequencing and immunophenotyping. Non-rejected clones showed high infiltration of tumor-associated-macrophages (TAMs), lower T-cell infiltration, and increased T-cell exhaustion compared to rejected clones. Comparative analysis of rejection-associated gene expression programs, combined with in-vivo CRISPR knockout screens of candidate mediators, identified Mif (macrophage migration inhibitory factor) as a regulator of immune rejection. Mif knockout led to smaller tumors and reversed non-rejection-associated immune composition, particularly, leading to the reduction of immunosuppressive macrophage infiltration. Finally, we validated these results in melanoma patient data.

The landscape of T cell antigens for cancer immunotherapy

The remarkable capacity of immunotherapies to induce durable regression in some patients with metastatic cancer relies heavily on T cell recognition of tumor-presented antigens. As checkpoint-blockade therapy has limited efficacy, tumor antigens have the potential to be exploited for complementary treatments, many of which are already in clinical trials. The surge of interest in this topic has led to the expansion of the tumor antigen landscape with the emergence of new antigen categories. Nonetheless, how different antigens compare in their ability to elicit efficient and safe clinical responses remains largely unknown. Here, we review known cancer peptide antigens, their attributes and the relevant clinical data and discuss future directions.

The proteasome regulator PSME4 modulates proteasome activity and antigen diversity to abrogate antitumor immunity in NSCLC

Immunotherapy revolutionized treatment options in cancer, yet the mechanisms underlying resistance in many patients remain poorly understood. Cellular proteasomes have been implicated in modulating antitumor immunity by regulating antigen processing, antigen presentation, inflammatory signaling and immune cell activation. However, whether and how proteasome complex heterogeneity may affect tumor progression and the response to immunotherapy has not been systematically examined. Here, we show that proteasome complex composition varies substantially across cancers and impacts tumor-immune interactions and the tumor microenvironment. Through profiling of the degradation landscape of patient-derived non-small-cell lung carcinoma samples, we find that the proteasome regulator PSME4 is upregulated in tumors, alters proteasome activity, attenuates presented antigenic diversity and associates with lack of response to immunotherapy. Collectively, our approach affords a paradigm by which proteasome composition heterogeneity and function should be examined across cancer types and targeted in the context of precision oncology.

Large-Scale Immuno-Peptidome Analysis Reveals Recurrent Post-Translational Splicing of Cancer and Immune-Associated Genes

Post-translational spliced peptides (PTSPs) are a unique class of peptides that have been found to be presented by HLA-class-I molecules in cancer (1). Thus far, no consensus has been reached on the proportion of PTSPs in the immunopeptidome, with estimates ranging from 2% to as high as 45% and stirring significant debate (2-8). Furthermore, the role of the HLA-class-II pathway in PTSP presentation has been studied only in diabetes (9). Here, we exploit our large-scale cancer peptidomics database and our newly devised pipeline for filtering spliced peptide predictions to identify recurring spliced peptides, both for HLA-class-I and -II complexes. Our results indicate that HLA-class-I spliced peptides account for a low percentage of the immunopeptidome (less than 3.1%), yet are larger in number relative to other types of identified aberrant peptides. Therefore, spliced peptides significantly contribute to the repertoire of presented peptides in cancer cells. In addition, we identified HLA-class-II-bound spliced peptides, but to a lower extent (less than 0.5%). The identified spliced peptides include cancer- and immune-associated genes, such as the MITF oncogene, DAPK1 tumor suppressor and HLA-E, which were validated using synthetic peptides. The potential immunogenicity of the DAPK1- and HLA-E-derived PTSPs was also confirmed. In addition, a reanalysis of our published mouse single-cell clone immunopeptidome dataset showed that most of the spliced peptides were found repeatedly in a large number of the single-cell clones. Establishing a novel search-scheme for the discovery and evaluation of recurring PTSPs among cancer patients may assist in identifying potential novel targets for immunotherapy.

NRas activity is regulated by dynamic interactions with nanoscale signaling clusters at the plasma membrane

NRas is a key mediator of the mitogenic pathway in normal cells and in cancer cells. Its dynamics and nanoscale organization at the plasma membrane (PM) facilitate its signaling. Here, we used two-color photoactivated localization microscopy to resolve the organization of individual NRas and associated signaling proteins in live melanoma cells, with resolution down to ∼20 nm. Upon EGF activation, a fraction of NRas and BRAF (dis)assembled synchronously at the PM in co-clusters. NRas and BRAF clusters associated with GPI-enriched domains, serving as possible nucleation sites for these clusters. NRas and BRAF association in mutual clusters was reduced by the NRas farnesylation inhibitor lonafarnib, yet enhanced by the BRAF inhibitor vemurafenib. Surprisingly, dispersed NRas molecules associated with the periphery of self-clusters of either Grb2 or NF1. Thus, NRas-mediated signaling, which is critical in health and disease, is regulated by dynamic interactions with functional clusters of BRAF or other related proteins at the PM.

Intratumor Heterogeneity and Antitumor Immunity Shape One Another Bidirectionally

Over the last decade, it has become clear that the genomic landscapes of tumors profoundly impact their immunogenicity and how tumor cells interact with immune cells. Whereas past discoveries mainly focused on the interplay between tumor immunogenicity and tumor mutational burden (TMB), under the assumption that a higher mutation load would give rise to a better patient response to immune checkpoint blockade therapies, we and others have underlined intratumor heterogeneity (ITH) as an important determinant of the magnitude of the antitumor response and the nature of the tumor microenvironment. In this review, we define TMB versus ITH and how the two factors are being inferred from data, examine key findings in the cancer immunogenomics literature deciphering the complex cross-talk between TMB, ITH, and antitumor immunity in human cancers and in vivo models, and discuss the mutual influence of ITH and immunity-how the antitumor response can give rise to tumors with higher ITH, and how higher ITH can put shackles on the antitumor response.

Mechanisms of immune activation and regulation: lessons from melanoma

Melanoma, a skin cancer that develops from pigment cells, has been studied intensively, particularly in terms of the immune response to tumours, and has been used as a model for the development of immunotherapy. This is due, in part, to the high mutational burden observed in melanomas, which increases both their immunogenicity and the infiltration of immune cells into the tumours, compared with other types of cancers. The immune response to melanomas involves a complex set of components and interactions. As the tumour evolves, it accumulates an increasing number of genetic and epigenetic alterations, some of which contribute to the immunogenicity of the tumour cells and the infiltration of immune cells. However, tumour evolution also enables the development of resistance mechanisms, which, in turn, lead to tumour immune escape. Understanding the interactions between melanoma tumour cells and the immune system, and the evolving changes within the melanoma tumour cells, the immune system and the microenvironment, is essential for the development of new cancer therapies. However, current research suggests that other extrinsic factors, such as the microbiome, may play a role in the immune response to melanomas. Here, we review the mechanisms underlying the immune response in the tumour and discuss recent advances as well as strategies for treatment development.

Dietary fiber and probiotics influence the gut microbiome and melanoma immunotherapy response

Gut bacteria modulate the response to immune checkpoint blockade (ICB) treatment in cancer, but the effect of diet and supplements on this interaction is not well studied. We assessed fecal microbiota profiles, dietary habits, and commercially available probiotic supplement use in melanoma patients and performed parallel preclinical studies. Higher dietary fiber was associated with significantly improved progression-free survival in 128 patients on ICB, with the most pronounced benefit observed in patients with sufficient dietary fiber intake and no probiotic use. Findings were recapitulated in preclinical models, which demonstrated impaired treatment response to anti–programmed cell death 1 (anti–PD-1)–based therapy in mice receiving a low-fiber diet or probiotics, with a lower frequency of interferon-γ–positive cytotoxic T cells in the tumor microenvironment. Together, these data have clinical implications for patients receiving ICB for cancer.

Combined presentation and immunogenicity analysis reveals a recurrent RAS.Q61K neoantigen in melanoma

Neoantigens are now recognized drivers of the antitumor immune response. Recurrent neoantigens, shared among groups of patients, have thus become increasingly coveted therapeutic targets. Here, we report on the data-driven identification of a robustly presented, immunogenic neoantigen that is derived from the combination of HLA-A*01:01 and RAS.Q61K. Analysis of large patient cohorts indicated that this combination applies to 3% of patients with melanoma. Using HLA peptidomics, we were able to demonstrate robust endogenous presentation of the neoantigen in 10 tumor samples. We detected specific reactivity to the mutated peptide within tumor-infiltrating lymphocytes (TILs) from 2 unrelated patients, thus confirming its natural immunogenicity. We further investigated the neoantigen-specific clones and their T cell receptors (TCRs) via a combination of TCR sequencing, TCR overexpression, functional assays, and single-cell transcriptomics. Our analysis revealed a diverse repertoire of neoantigen-specific clones with both intra- and interpatient TCR similarities. Moreover, 1 dominant clone proved to cross-react with the highly prevalent RAS.Q61R variant. Transcriptome analysis revealed a high association of TCR clones with specific T cell phenotypes in response to cognate melanoma, with neoantigen-specific cells showing an activated and dysfunctional phenotype. Identification of recurrent neoantigens and their reactive TCRs can promote "off-the-shelf" precision immunotherapies, alleviating limitations of personalized treatments.

Identification of presented SARS-CoV-2 HLA class I and HLA class II peptides using HLA peptidomics

The human leukocyte antigen (HLA)-bound viral antigens serve as an immunological signature that can be selectively recognized by T cells. As viruses evolve by acquiring mutations, it is essential to identify a range of presented viral antigens. Using HLA peptidomics, we are able to identify severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-derived peptides presented by highly prevalent HLA class I (HLA-I) molecules by using infected cells as well as overexpression of SARS-CoV-2 genes. We find 26 HLA-I peptides and 36 HLA class II (HLA-II) peptides. Among the identified peptides, some are shared between different cells and some are derived from out-of-frame open reading frames (ORFs). Seven of these peptides were previously shown to be immunogenic, and we identify two additional immunoreactive peptides by using HLA multimer staining. These results may aid the development of the next generation of SARS-CoV-2 vaccines based on presented viral-specific antigens that span several of the viral genes.

Melanoma models for the next generation of therapies

There is a lack of appropriate melanoma models that can be used to evaluate the efficacy of novel therapeutic modalities. Here, we discuss the current state of the art of melanoma models including genetically engineered mouse, patient-derived xenograft, zebrafish, and ex vivo and in vitro models. We also identify five major challenges that can be addressed using such models, including metastasis and tumor dormancy, drug resistance, the melanoma immune response, and the impact of aging and environmental exposures on melanoma progression and drug resistance. Additionally, we discuss the opportunity for building models for rare subtypes of melanomas, which represent an unmet critical need. Finally, we identify key recommendations for melanoma models that may improve accuracy of preclinical testing and predict efficacy in clinical trials, to help usher in the next generation of melanoma therapies.

Identification of bacteria-derived HLA-bound peptides in melanoma

A variety of species of bacteria are known to colonize human tumours1-11, proliferate within them and modulate immune function, which ultimately affects the survival of patients with cancer and their responses to treatment12-14. However, it is not known whether antigens derived from intracellular bacteria are presented by the human leukocyte antigen class I and II (HLA-I and HLA-II, respectively) molecules of tumour cells, or whether such antigens elicit a tumour-infiltrating T cell immune response. Here we used 16S rRNA gene sequencing and HLA peptidomics to identify a peptide repertoire derived from intracellular bacteria that was presented on HLA-I and HLA-II molecules in melanoma tumours. Our analysis of 17 melanoma metastases (derived from 9 patients) revealed 248 and 35 unique HLA-I and HLA-II peptides, respectively, that were derived from 41 species of bacteria. We identified recurrent bacterial peptides in tumours from different patients, as well as in different tumours from the same patient. Our study reveals that peptides derived from intracellular bacteria can be presented by tumour cells and elicit immune reactivity, and thus provides insight into a mechanism by which bacteria influence activation of the immune system and responses to therapy.

PhISCS-BnB: a fast branch and bound algorithm for the perfect tumor phylogeny reconstruction problem

MOTIVATION

Recent advances in single-cell sequencing (SCS) offer an unprecedented insight into tumor emergence and evolution. Principled approaches to tumor phylogeny reconstruction via SCS data are typically based on general computational methods for solving an integer linear program, or a constraint satisfaction program, which, although guaranteeing convergence to the most likely solution, are very slow. Others based on Monte Carlo Markov Chain or alternative heuristics not only offer no such guarantee, but also are not faster in practice. As a result, novel methods that can scale up to handle the size and noise characteristics of emerging SCS data are highly desirable to fully utilize this technology.

RESULTS

We introduce PhISCS-BnB (phylogeny inference using SCS via branch and bound), a branch and bound algorithm to compute the most likely perfect phylogeny on an input genotype matrix extracted from an SCS dataset. PhISCS-BnB not only offers an optimality guarantee, but is also 10-100 times faster than the best available methods on simulated tumor SCS data. We also applied PhISCS-BnB on a recently published large melanoma dataset derived from the sublineages of a cell line involving 20 clones with 2367 mutations, which returned the optimal tumor phylogeny in <4 h. The resulting phylogeny agrees with and extends the published results by providing a more detailed picture on the clonal evolution of the tumor.

AVAILABILITY AND IMPLEMENTATION

https://github.com/algo-cancer/PhISCS-BnB.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

MEK Inhibition Reverses Aberrant Signaling in Melanoma Cells through Reorganization of NRas and BRAF in Self Nanoclusters

Hotspot mutations of the oncogenes BRAF and NRas are the most common genetic alterations in cutaneous melanoma. Still, the nanoscale organization and signal coupling of these proteins remain incompletely understood, particularly upon expression of oncogenic NRas mutants. Here we employed single-molecule localization microscopy to study the nanoscale organization of NRas and BRAF at the plasma membrane (PM) of melanoma cells. NRas and BRAF resided in self-clusters that did not associate well in resting cells. In EGF-activated cells, NRas clusters became more diffused while overall protein levels at the PM increased; thus allowing enhanced association of NRas and BRAF and downstream signaling. In multiple melanoma cell lines, mutant NRas resided in more pronounced self-clusters relative to wild-type (WT) NRas yet associated more with the clustered and more abundant BRAF. In cells resistant to trametinib, a clinical MEK inhibitor (MEKi), a similar coclustering of NRas and BRAF was observed upon EGF activation. Strikingly, treatment of cells expressing mutant NRas with trametinib reversed the effect of mutant NRas expression by restoring the nonoverlapping self-clusters of NRas and BRAF and by reducing their PM levels and elevated pERK levels caused by mutant NRas. Our results indicate a new mechanism for signal regulation of NRas in melanoma through its nanoscale dynamic organization and a new mechanism for MEKi function in melanoma cells carrying NRas mutations but lacking MEK mutations. SIGNIFICANCE: Nanoscale dynamic organization of WT and mutant NRas relative to BRAF serves as a regulatory mechanism for NRas signaling and may be a viable therapeutic target for its sensitivity to MEKi.

Anti-tumour immunity induces aberrant peptide presentation in melanoma

Extensive tumour inflammation, which is reflected by high levels of infiltrating T cells and interferon-γ (IFNγ) signalling, improves the response of patients with melanoma to checkpoint immunotherapy^1,2. Many tumours, however, escape by activating cellular pathways that lead to immunosuppression. One such mechanism is the production of tryptophan metabolites along the kynurenine pathway by the enzyme indoleamine 2,3-dioxygenase 1 (IDO1), which is induced by IFNγ^3-5. However, clinical trials using inhibition of IDO1 in combination with blockade of the PD1 pathway in patients with melanoma did not improve the efficacy of treatment compared to PD1 pathway blockade alone^6,7, pointing to an incomplete understanding of the role of IDO1 and the consequent degradation of tryptophan in mRNA translation and cancer progression. Here we used ribosome profiling in melanoma cells to investigate the effects of prolonged IFNγ treatment on mRNA translation. Notably, we observed accumulations of ribosomes downstream of tryptophan codons, along with their expected stalling at the tryptophan codon. This suggested that ribosomes bypass tryptophan codons in the absence of tryptophan. A detailed examination of these tryptophan-associated accumulations of ribosomes-which we term 'W-bumps'-showed that they were characterized by ribosomal frameshifting events. Consistently, reporter assays combined with proteomic and immunopeptidomic analyses demonstrated the induction of ribosomal frameshifting, and the generation and presentation of aberrant trans-frame peptides at the cell surface after treatment with IFNγ. Priming of naive T cells from healthy donors with aberrant peptides induced peptide-specific T cells. Together, our results suggest that IDO1-mediated depletion of tryptophan, which is induced by IFNγ, has a role in the immune recognition of melanoma cells by contributing to diversification of the peptidome landscape.

Abstract IA07: UVB-induced tumor heterogeneity directs immune response in melanoma

Little is known regarding the relationship between intra-tumor heterogeneity (ITH) and immune response in melanoma. Here, we explored the role of ITH in tumor rejection by establishing a melanoma mouse model and inducing UVB-derived mutations that increase both ITH and mutational load. This induction gives rise to highly aggressive tumors and decreased cytotoxic activity of tumor infiltrating lymphocytes (TILs). Conversely, single cell-derived melanoma clones with reduced ITH are swiftly rejected. Tumor rejection is accompanied by increased TIL reactivity and increased infiltration into the tumor core. Using phylogenetic tree analyses and mixing experiments of 20 single cell clones that lie along the phylogenetic tee we show that tumor rejection is strongly affected by number of injected clones and genetic diversity. We have, thus set up a novel, highly controlled system that enables us to study the interphase between the immune system and different layers of intra-tumor heterogeneity. Finally, the analysis of melanoma patient data identifies parallel observations, supporting the importance of ITH in determining patient survival and response to checkpoint blockade.

Citation Format: Osnat Bartok, Sushant Patkar, Sapir Cohen, Kevin Litchfield, Hiren Karathia, Joo Sang Lee, Alejandro Jiménez-Sánchez, Chi-Ping Day, Lea Eisenbach, Martin Miller, Glenn Merlino, Eli Pikarsky, Arie Admon, Charles Swanton, Eytan Ruppin, Yardena Samuels, Yochai Wolf. UVB-induced tumor heterogeneity directs immune response in melanoma [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr IA07.

A genome-wide CRISPR screen identifies FBXO42 involvement in resistance toward MEK inhibition in NRAS-mutant melanoma

NRAS mutations are the most common alterations among RAS isoforms in cutaneous melanoma, with patients harboring these aggressive tumors having a poor prognosis and low survival rate. The main line of treatment for these patients is MAPK pathway‐targeted therapies, such as MEK inhibitors, but, unfortunately, the response to these inhibitors is variable due to tumor resistance. Identifying genetic modifiers involved in resistance toward MEK‐targeted therapy may assist in the development of new therapeutic strategies, enhancing treatment response and patient survival. Our whole‐genome CRISPR‐Cas9 knockout screen identified the target Kelch domain‐containing F‐Box protein 42 (FBXO42) as a factor involved in NRAS‐mutant melanoma‐acquired resistance to the MEK1/2 inhibitor trametinib. We further show that FBXO42, an E3 ubiquitin ligase, is involved in the TAK1 signaling pathway, possibly prompting an increase in active P38. In addition, we demonstrate that combining trametinib with the TAK1 inhibitor, takinib, is a far more efficient treatment than trametinib alone in NRAS‐mutant melanoma cells. Our findings thus show a new pathway involved in NRAS‐mutant melanoma resistance and provide new opportunities for novel therapeutic options.

UVB-Induced Tumor Heterogeneity Diminishes Immune Response in Melanoma

Although clonal neo-antigen burden is associated with improved response to immune therapy, the functional basis for this remains unclear. Here we study this question in a novel controlled mouse melanoma model that enables us to explore the effects of intra-tumor heterogeneity (ITH) on tumor aggressiveness and immunity independent of tumor mutational burden. Induction of UVB-derived mutations yields highly aggressive tumors with decreased anti-tumor activity. However, single-cell-derived tumors with reduced ITH are swiftly rejected. Their rejection is accompanied by increased T cell reactivity and a less suppressive microenvironment. Using phylogenetic analyses and mixing experiments of single-cell clones, we dissect two characteristics of ITH: the number of clones forming the tumor and their clonal diversity. Our analysis of melanoma patient tumor data recapitulates our results in terms of overall survival and response to immune checkpoint therapy. These findings highlight the importance of clonal mutations in robust immune surveillance and the need to quantify patient ITH to determine the response to checkpoint blockade.

Cross-Talk between Receptor Tyrosine Kinases AXL and ERBB3 Regulates Invadopodia Formation in Melanoma Cells

The invasive phenotype of metastatic cancer cells is accompanied by the formation of actin-rich invadopodia, which adhere to the extracellular matrix and degrade it. In this study, we explored the role of the tyrosine kinome in the formation of invadopodia in metastatic melanoma cells. Using a microscopy-based siRNA screen, we identified a series of regulators, the knockdown of which either suppresses (e.g., TYK2, IGFR1, ERBB3, TYRO3, FES, ALK, PTK7) or enhances (e.g., ABL2, AXL, CSK) invadopodia formation and function. Notably, the receptor tyrosine kinase AXL displayed a dual regulatory function, where both depletion or overexpression enhanced invadopodia formation and activity. This apparent contradiction was attributed to the capacity of AXL to directly stimulate invadopodia, yet its suppression upregulates the ERBB3 signaling pathway, which can also activate core invadopodia regulators and enhance invadopodia function. Bioinformatic analysis of multiple melanoma cell lines points to an inverse expression pattern of AXL and ERBB3. High expression of AXL in melanoma cells is associated with high expression of invadopodia components and an invasive phenotype. These results provide new insights into the complexity of metastasis-promoting mechanisms and suggest that targeting of multiple invadopodia signaling networks may serve as a potential anti-invasion therapy in melanoma. SIGNIFICANCE: These findings uncover a unique interplay between AXL and ERBB3 in invadopodia regulation that points to the need for combined therapy in order to prevent invadopodia-mediated metastasis in melanoma.

Exome Sequencing of ABCB5 Identifies Recurrent Melanoma Mutations that Result in Increased Proliferative and Invasive Capacities

ABCB5 is an ABC transporter that was shown to confer low-level multidrug resistance in cancer. In this study, we show that ABCB5 was mutated in 13.75% of the 640 melanoma samples analyzed. Besides nonsense mutations, two mutation hotspots were found in the ABCB5 protein, in the drug-binding pocket and the nucleotide-binding domains. Four mutations, which are representative of the mutation pattern, were selected. ATPase assays showed that these mutations resulted in a decrease in basal ATP hydrolysis by ABCB5. To select informative melanoma cell lines, mutational profiles of the clinical samples were further analyzed. This study showed mutations in the tumor suppressor CDKN2A gene and the NRAS oncogene in 62.5% and 75%, respectively of the samples that had mutations in the ABCB5 gene. No mutation was found in the tumor suppressor PTEN gene, whereas the activating V600E mutation in the BRAF oncogene was found in 25% of the samples with a mutated ABCB5 gene. Studies in four melanoma cell lines with various genetic backgrounds showed an increase in the proliferation and migration capacity of mutant ABCB5-expressing cells, suggesting that ABCB5 plays a role in the development of melanoma as a tumor suppressor gene.

Author Correction: RGS7 is recurrently mutated in melanoma and promotes migration and invasion of human cancer cells

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Correction: CAPN1 is a novel binding partner and regulator of the tumor suppressor NF1 in melanoma

[This corrects the article DOI: 10.18632/oncotarget.25805.].

Cancer Exome-Based Identification of Tumor Neo-Antigens Using Mass Spectrometry

Neo-antigens expressed on tumors are targets for development of cancer immunotherapy strategies. Use of prediction algorithms to identify neo-antigens yields a significant number of peptides that must be validated in laborious and time-consuming methods; many prove to be false-positive identifications. The use of HLA peptidomics allows the isolation of the HLA-peptide complexes directly from cells and can be done on fresh tumor, patient-derived xerographs, or cell lines when the tissue sample is limited. This method can be used to identify both HLA class I and HLA class II or any different MHC from different species. Here we describe the steps to create the immune-affinity columns used from the process, the immunoprecipitation procedure, and also the isolation of the peptides that will be analyzed by mass spectrometry.

Cancer research in the era of immunogenomics

The most meaningful advancement in cancer treatment in recent years has been the emergence of immunotherapy. Checkpoint inhibitor blockade and adoptive T cell therapy have shown remarkable clinical effects in a wide range of tumour types. Despite these advances, many tumours do not respond to these treatments, which raises the need to further investigate how patients can benefit from immunotherapy. This effort can now take advantage of the recent technological progress in single-cell, high-throughput sequencing and computational efforts. In this review, we will discuss advances in different immunotherapies and the principles of cancer immunogenomics, with an emphasis on the detection of cancer neoantigens with human leucocyte antigen peptidomics, and how these principles can be further used for more efficient clinical output.

Combined Analysis of Antigen Presentation and T-cell Recognition Reveals Restricted Immune Responses in Melanoma

The quest for tumor-associated antigens (TAA) and neoantigens is a major focus of cancer immunotherapy. Here, we combine a neoantigen prediction pipeline and human leukocyte antigen (HLA) peptidomics to identify TAAs and neoantigens in 16 tumors derived from seven patients with melanoma and characterize their interactions with their tumor-infiltrating lymphocytes (TIL). Our investigation of the antigenic and T-cell landscapes encompassing the TAA and neoantigen signatures, their immune reactivity, and their corresponding T-cell identities provides the first comprehensive analysis of cancer cell T-cell cosignatures, allowing us to discover remarkable antigenic and TIL similarities between metastases from the same patient. Furthermore, we reveal that two neoantigen-specific clonotypes killed 90% of autologous melanoma cells, both in vitro and in vivo, showing that a limited set of neoantigen-specific T cells may play a central role in melanoma tumor rejection. Our findings indicate that combining HLA peptidomics with neoantigen predictions allows robust identification of targetable neoantigens, which could successfully guide personalized cancer immunotherapies.Significance: As neoantigen targeting is becoming more established as a powerful therapeutic approach, investigating these molecules has taken center stage. Here, we show that a limited set of neoantigen-specific T cells mediates tumor rejection, suggesting that identifying just a few antigens and their corresponding T-cell clones could guide personalized immunotherapy. Cancer Discov; 8(11); 1366-75. ©2018 AACR. This article is highlighted in the In This Issue feature, p. 1333.

Urea Cycle Dysregulation Generates Clinically Relevant Genomic and Biochemical Signatures

The urea cycle (UC) is the main pathway by which mammals dispose of waste nitrogen. We find that specific alterations in the expression of most UC enzymes occur in many tumors, leading to a general metabolic hallmark termed "UC dysregulation" (UCD). UCD elicits nitrogen diversion toward carbamoyl-phosphate synthetase2, aspartate transcarbamylase, and dihydrooratase (CAD) activation and enhances pyrimidine synthesis, resulting in detectable changes in nitrogen metabolites in both patient tumors and their bio-fluids. The accompanying excess of pyrimidine versus purine nucleotides results in a genomic signature consisting of transversion mutations at the DNA, RNA, and protein levels. This mutational bias is associated with increased numbers of hydrophobic tumor antigens and a better response to immune checkpoint inhibitors independent of mutational load. Taken together, our findings demonstrate that UCD is a common feature of tumors that profoundly affects carcinogenesis, mutagenesis, and immunotherapy response.

CAPN1 is a novel binding partner and regulator of the tumor suppressor NF1 in melanoma

Neurofibromin 1 (NF1), a tumor suppressor that negatively regulates RAS through its GTPase activity, is highly mutated in various types of sporadic human cancers, including melanoma. However, the binding partners of NF1 and the pathways in which it is involved in melanoma have not been characterized in an in depth manner. Utilizing a mass spectrometry analysis of NF1 binding partners, we revealed Calpain1 (CAPN1), a calcium-dependent neutral cysteine protease, as a novel NF1 binding partner that regulates NF1 degradation in melanoma cells. ShRNA-mediated knockdown of CAPN1 or treatment with a CAPN1 inhibitor stabilizes NF1 protein levels, downregulates AKT signaling and melanoma cell growth. Combination treatment of Calpain inhibitor I with MEKi Trametinib in different melanoma cells is more effective in reducing melanoma cell growth compared to treatment with Trametinib alone, suggesting that this combination may have a therapeutic potential in melanoma. This novel mechanism for regulating NF1 in melanoma provides a molecular basis for targeting CAPN1 in order to stabilize NF1 levels and, in doing so, suppressing Ras activation; this mechanism can be exploited therapeutically in melanoma and other cancers.

Harnessing synthetic lethality to predict the response to cancer treatment

While synthetic lethality (SL) holds promise in developing effective cancer therapies, SL candidates found via experimental screens often have limited translational value. Here we present a data-driven approach, ISLE (identification of clinically relevant synthetic lethality), that mines TCGA cohort to identify the most likely clinically relevant SL interactions (cSLi) from a given candidate set of lab-screened SLi. We first validate ISLE via a benchmark of large-scale drug response screens and by predicting drug efficacy in mouse xenograft models. We then experimentally test a select set of predicted cSLi via new screening experiments, validating their predicted context-specific sensitivity in hypoxic vs normoxic conditions and demonstrating cSLi's utility in predicting synergistic drug combinations. We show that cSLi can successfully predict patients' drug treatment response and provide patient stratification signatures. ISLE thus complements existing actionable mutation-based methods for precision cancer therapy, offering an opportunity to expand its scope to the whole genome.

Refinement of the endogenous epitope tagging technology allows the identification of a novel NRAS binding partner in melanoma

The NRAS oncoprotein is highly mutated in melanoma. However, to date, no comprehensive proteomic study has been reported for NRAS. Here, we utilized the endogenous epitope tagging (EET) approach for the identification of novel NRAS binding partners. Using EET, an epitope tag is added to the endogenously expressed protein, via modification of its genomic coding sequence. Existing EET systems are not robust, suffer from high background, and are labor-intensive. To this end, we present a polyadenylation signal-trap construct for N'-tagging that generates a polycistronic mRNA with the gene of interest. This system requires the integration of the tagging cassette in frame with the target gene to be expressed. Using this design, we demonstrate, for the first time, endogenous tagging of NRAS in melanoma cells allowing the identification of the E3 ubiquitin ligase c-CBL as a novel NRAS binding partner. Thus, our developed EET technology allows the characterization of new RAS effectors, which could be beneficial for the design of future drugs that inhibit constitutive signaling of RAS oncogenic mutants.

Targeting the Senescence-Overriding Cooperative Activity of Structurally Unrelated H3K9 Demethylases in Melanoma

Oncogene-induced senescence, e.g., in melanocytic nevi, terminates the expansion of pre-malignant cells via transcriptional silencing of proliferation-related genes due to decoration of their promoters with repressive trimethylated histone H3 lysine 9 (H3K9) marks. We show here that structurally distinct H3K9-active demethylases-the lysine-specific demethylase-1 (LSD1) and several Jumonji C domain-containing moieties (such as JMJD2C)-disable senescence and permit Ras/Braf-evoked transformation. In mouse and zebrafish models, enforced LSD1 or JMJD2C expression promoted Braf-V600E-driven melanomagenesis. A large subset of established melanoma cell lines and primary human melanoma samples presented with a collective upregulation of related and unrelated H3K9 demethylase activities, whose targeted inhibition restored senescence, even in Braf inhibitor-resistant melanomas, evoked secondary immune effects and controlled tumor growth in vivo.

RGS7 is recurrently mutated in melanoma and promotes migration and invasion of human cancer cells

Analysis of 501 melanoma exomes revealed RGS7, which encodes a GTPase-accelerating protein (GAP), to be a tumor-suppressor gene. RGS7 was mutated in 11% of melanomas and was found to harbor three recurrent mutations (p.R44C, p.E383K and p.R416Q). Structural modeling of the most common recurrent mutation of the three (p.R44C) predicted that it destabilizes the protein due to the loss of an H-bond and salt bridge network between the mutated position and the serine and aspartic acid residues at positions 58 as 61, respectively. We experimentally confirmed this prediction showing that the p.R44C mutant protein is indeed destabilized. We further show RGS7 p.R44C has weaker catalytic activity for its substrate Gα_o, thus providing a dual mechanism for its loss of function. Both of these effects are expected to contribute to loss of function of RGS7 resulting in increased anchorage-independent growth, migration and invasion of melanoma cells. By mutating position 56 in the R44C mutant from valine to cysteine, thereby enabling the formation of a disulfide bridge between the two mutated positions, we slightly increased the catalytic activity and reinstated protein stability, leading to the rescue of RGS7′s function as a tumor suppressor. Our findings identify RGS7 as a novel melanoma driver and point to the clinical relevance of using strategies to stabilize the protein and, thereby, restore its function.

Analysis of Enzymatic Activity of Matrix Metalloproteinase (MMP) by Collagen Zymography in Melanoma

Protein zymography is the most commonly used technique to study the enzymatic activity of matrix metalloproteinases (MMPs) and their inhibitors. MMPs are proteolytic enzymes that promote extracellular matrix degradation. MMPs are frequently mutated in malignant melanomas as well as other cancers and are linked to increasing incidence of tumor metastasis. Substrate zymography characterizes MMP activity by their ability to degrade preferred substrates. Here we describe the collagen zymography technique to measure the active or latent form of MMPs using MMP-8 as an example, which is a frequently mutated MMP family member in malignant melanomas. The same technique can be used with the modification of substrate to detect metalloproteinase activity of other MMPs. Both wild-type and mutated forms of MMPs can be analyzed using a single gel using this method.

Abstract A188: Harnessing synthetic lethality to predict the response to cancer treatments

Synthetic lethality (SL) describes an interaction between a pair of genes whereby their double knockout is lethal, while their respective knockout is not. The identification of SL interactions (SLi) via large-scale genomic screens offers promising opportunities for developing selective therapies in cancer. However, our analysis of the TCGA cohort shows that many of the interactions do not carry predictive signal of patient survival or drug response. Here we present a data-driven approach termed ISLE (Identification of clinically relevant Synthetic LEthality) that mines the TCGA cohort to identify a subset of clinically relevant SL interactions (cSLi). ISLE consists of the following inference steps, analysis of tumor, cell line, and gene evolutionary data. We first create an initial pool of SL pairs identified through direct double knockout screens/isogenic cell line screens or inferred from large-scale shRNA/sgRNA single-gene knockout screens. Starting from this initial SL pool, ISLE first identifies putative SL gene pairs whose co-inactivation is under-represented in tumors, testifying that it is selected against. Second, it prioritizes candidate SL pairs whose co-inactivation is associated with improved patient’s prognosis, testifying that it may hamper tumor progression. Finally, it prioritizes SL-gene pairs with similar evolutionary phylogenetic profiles based on the notion that SL interactions are conserved across multiple species. We validate the identified SL pairs using an unseen large-scale in vitro drug response screen by showing the SL pairs marks a decent prediction accuracy (AUC~0.8). We compare ISLE’s performance to the standard supervised drug response prediction approaches in DREAM challenges, and our prediction based on generic pretreatment tumor samples (from TCGA) was within top 3 in prediction accuracy among the top predictors. ISLE-based approach also successfully distinguishes responders vs nonresponders to drug treatment (for &gt;70% of drugs) in mouse xenografts using the activity profile of the drug target’s SL-partners. We then experimentally show the utility of SL in predicting synergistic drug combinations in patient-derived cell lines based on the notion that the two drugs whose targets have SL interactions are synergistic. Most importantly, we demonstrate for the first time that an SL network can successfully predict the treatment outcome in cancer patients in multiple large-scale patient datasets including TCGA, where cSLi are successfully predict patients’ response for more than 70% of cancer drugs. ISLE is predictive of patients’ response for the majority of current cancer drugs without any drug-specific training. Of paramount importance, the predictions of ISLE are based on SLi between (potentially) all genes in the cancer genome, thus prioritizing treatments for patients whose tumors do not bear specific actionable mutations in cancer driver genes, offering a novel approach to precision-based cancer therapy.

Citation Format: Joo S. Lee, Avinash Das, Livnat Jerby-Arnon, Rand Arafeh, Matthew Davidson, Arnaud Amzallag, Seung Gu Park, Kuoyuan Cheng, Welles Robinson, Dikla Atias, Chani Stossel, Ella Buzhor, Gidi Stein, Joshua J. Waterfall, Paul S. Meltzer, Talia Golan, Sridhar Hannenhalli, Eyal Gottlieb, Cyril H. Benes, Yardena Samuels, Emma Shanks, Eytan Ruppin. Harnessing synthetic lethality to predict the response to cancer treatments [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2017 Oct 26-30; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Ther 2018;17(1 Suppl):Abstract nr A188.

Combined inhibition of MEK and nuclear ERK translocation has synergistic antitumor activity in melanoma cells

Genetic alterations in BRAF, NRAS and NF1 that activate the ERK cascade, account for over 80% of metastatic melanomas. However, ERK cascade inhibitors have been proven beneficial almost exclusively for BRAF mutant melanomas. One of the hallmarks of the ERK cascade is the nuclear translocation of ERK1/2, which is important mainly for the induction of proliferation. This translocation can be inhibited by the NTS-derived peptide (EPE) that blocks the ERK1/2-importin7 interaction, inhibits the nuclear translocation of ERK1/2, and arrests active ERK1/2 in the cytoplasm. In this study, we found that the EPE peptide significantly reduced the viability of not only BRAF, but also several NRAS and NF1 mutant melanomas. Importantly, combination of the EPE peptide and trametinib showed synergy in reducing the viability of some NRAS mutant melanomas, an effect driven by the partial preservation of negative feedback loops. The same combination significantly reduced the viability of other melanoma cells, including those resistant to mono-treatment with EPE peptide and ERK cascade inhibitors. Our study indicates that targeting the nuclear translocation of ERK1/2, in combination with MEK inhibitors can be used for the treatment of different mutant melanomas.

Abstract LB-031: Deciphering distinct roles of RASA2 in melanomagenesis

Melanoma is the deadliest form of human skin cancer. The incidence of melanoma continues to rise. Recent advances in knowledge of melanoma genetics, genomics and biology has led to an optimistic view of the therapeutic outlook for melanoma patients. We analyzed sequence data from &gt;500 melanoma genomes/exomes to identify novel tumor suppressor genes in melanoma. RASA2 was identified as the most highly somatically mutated novel tumor suppressor gene. RASA2 was mutated in 5% of melanomas and deleted in an additional 16.4% of cases. RASA2 is a GTPase Activating Protein (GAP) that regulates RAS; which is one of the most highly mutated oncogenes in melanoma but drugs targeting RAS have as yet shown poor efficacy. The role of RASA2 has not been investigated in melanoma. NF1, which encodes another RAS- specific GAP, was found to be frequently mutated in melanoma. Interestingly, mutations in RASA2 and NF1 co-occur in the same patients with high frequency. We plan to elucidate the roles of RASA2 in melanomagenesis and to understand why RASA2 and NF1 mutations co-occur despite the fact that both proteins are RasGAPs. Ras includes three isoforms: NRas, KRas and HRas. Our preliminary data show that RASA2 is more specific to NRAS and that NF1 is more specific to KRAS and HRAS. This finding highlights the existence of a paradigm of cooperativity in which combined loss of multiple negative regulators (RASA2 and NF1) of the RAS pathway is required for melanoma development. Therefore, this type of enhancement of RAS signaling is possibly selected for in some melanomas. We will apply a proteomic screen using BioID to identify RASA2 and NF1 binding partners to provide insights into the functional effects and consequences of alterations in RASA2 and NF1. We expect that these studies will not only identify the cellular components that contribute to the Ras signaling pathway but will also identify potential novel therapeutic targets.

Citation Format: Rand Arafeh, Nouar Qutob, Rafi Rafi Emmanuel, Jason Madore, Abdel Elkahloun, James S. James S. Wilmott, Jared J. Gartner, Antonella Di Pizio, Ron Rotkopf, Ken Dutton-Regester, Victoria Hill, Antonia Pritchard, Jimmy C. Lin, Steven A Rosenberg, Javed Khan, Shifra Ben-Dor, Masha Y. Masha Y. Niv, Igor Ulitsky, Graham J Mann, Richard A. Scolyer, Nicholas K. Hayward, Yardena Samuels. Deciphering distinct roles of RASA2 in melanomagenesis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-031. doi:10.1158/1538-7445.AM2017-LB-031

Synonymous Somatic Variants in Human Cancer Are Not Infamous: A Plea for Full Disclosure in Databases and Publications

Single-nucleotide variants (SNVs) are the most frequent genetic changes found in human cancer. Most driver alterations are missense and nonsense variants localized in the coding region of cancer genes. Unbiased cancer genome sequencing shows that synonymous SNVs (sSNVs) can be found clustered in the coding regions of several cancer oncogenes or tumor suppressor genes suggesting purifying selection. sSNVs are currently underestimated, as they are usually discarded during analysis. Furthermore, several public databases do not display sSNVs, which can lead to analytical bias and the false assumption that this mutational event is uncommon. Recent progress in our understanding of the deleterious consequences of these sSNVs for RNA stability and protein translation shows that they can act as strong drivers of cancer, as demonstrated for several cancer genes such as TP53 or BCL2L12. It is therefore essential that sSNVs be properly reported and analyzed in order to provide an accurate picture of the genetic landscape of the cancer genome.

The functional relevance of somatic synonymous mutations in melanoma and other cancers

Recent technological advances in sequencing have flooded the field of cancer research with knowledge about somatic mutations for many different cancer types. Most cancer genomics studies focus on mutations that alter the amino acid sequence, ignoring the potential impact of synonymous mutations. However, accumulating experimental evidence has demonstrated clear consequences for gene function, leading to a widespread recognition of the functional role of synonymous mutations and their causal connection to various diseases. Here, we review the evidence supporting the direct impact of synonymous mutations on gene function via gene splicing; mRNA stability, folding, and translation; protein folding; and miRNA-based regulation of expression. These results highlight the functional contribution of synonymous mutations to oncogenesis and the need to further investigate their detection and prioritization for experimental assessment.

Parkinson disease (PARK) genes are somatically mutated in cutaneous melanoma

OBJECTIVE

To assess whether Parkinson disease (PD) genes are somatically mutated in cutaneous melanoma (CM) tissue, because CM occurs in patients with PD at higher rates than in the general population and PD is more common than expected in CM cohorts.

METHODS

We cross-referenced somatic mutations in metastatic CM detected by whole-exome sequencing with the 15 known PD (PARK) genes. We computed the empirical distribution of the sum of mutations in each gene (Smut) and of the number of tissue samples in which a given gene was mutated at least once (SSampl) for each of the analyzable genes, determined the 90th and 95th percentiles of the empirical distributions of these sums, and verified the location of PARK genes in these distributions. Identical analyses were applied to adenocarcinoma of lung (ADENOCA-LUNG) and squamous cell carcinoma of lung (SQUAMCA-LUNG). We also analyzed the distribution of the number of mutated PARK genes in CM samples vs the 2 lung cancers.

RESULTS

Somatic CM mutation analysis (n = 246) detected 315,914 mutations in 18,758 genes. Somatic CM mutations were found in 14 of 15 PARK genes. Forty-eight percent of CM samples carried ≥1 PARK mutation and 25% carried multiple PARK mutations. PARK8 mutations occurred above the 95th percentile of the empirical distribution for SMut and SSampl. Significantly more CM samples harbored multiple PARK gene mutations compared with SQUAMCA-LUNG (p = 0.0026) and with ADENOCA-LUNG (p < 0.0001).

CONCLUSIONS

The overrepresentation of somatic PARK mutations in CM suggests shared dysregulated pathways for CM and PD.

The involvement of mutant Rac1 in the formation of invadopodia in cultured melanoma cells

In this article, we discuss the complex involvement of a Rho-family GTPase, Rac1, in cell migration and in invadopodia-mediated matrix degradation. We discuss the involvement of invadopodia in invasive cell migration, and their capacity to promote cancer metastasis. Considering the regulation of invadopodia formation, we describe studies that demonstrate the role of Rac1 in the metastatic process, and the suggestion that this effect is attributable to the capacity of Rac1 to promote invadopodia formation. This notion is demonstrated here by showing that knockdown of Rac1 in melanoma cells expressing a wild-type form of this GTPase, reduces invadopodia-dependent matrix degradation. Interestingly, we also show that excessive activity of Rac1, displayed by the P29S, hyperactive, "fast cycling" mutant of Rac1, which is present in 5-10% of melanoma tumors, inhibits invadopodia function. Moreover, knockdown of this hyperactive mutant enhanced matrix degradation, indicating that excessive Rac1 activity by this mutant can negatively regulate invadopodia formation and function.

Extreme Vulnerability of IDH1 Mutant Cancers to NAD+ Depletion

Heterozygous mutation of IDH1 in cancers modifies IDH1 enzymatic activity, reprogramming metabolite flux and markedly elevating 2-hydroxyglutarate (2-HG). Here, we found that 2-HG depletion did not inhibit growth of several IDH1 mutant solid cancer types. To identify other metabolic therapeutic targets, we systematically profiled metabolites in endogenous IDH1 mutant cancer cells after mutant IDH1 inhibition and discovered a profound vulnerability to depletion of the coenzyme NAD+. Mutant IDH1 lowered NAD+ levels by downregulating the NAD+ salvage pathway enzyme nicotinate phosphoribosyltransferase (Naprt1), sensitizing to NAD+ depletion via concomitant nicotinamide phosphoribosyltransferase (NAMPT) inhibition. NAD+ depletion activated the intracellular energy sensor AMPK, triggered autophagy, and resulted in cytotoxicity. Thus, we identify NAD+ depletion as a metabolic susceptibility of IDH1 mutant cancers.

Recurrent inactivating RASA2 mutations in melanoma

Analysis of 501 melanoma exomes identified RASA2, encoding a RasGAP, as a tumor-suppressor gene mutated in 5% of melanomas. Recurrent loss-of-function mutations in RASA2 were found to increase RAS activation, melanoma cell growth and migration. RASA2 expression was lost in ≥30% of human melanomas and was associated with reduced patient survival. These findings identify RASA2 inactivation as a melanoma driver and highlight the importance of RasGAPs in cancer.

The head and neck cancer cell oncogenome: a platform for the development of precision molecular therapies

The recent elucidation of the genomic landscape of head and neck squamous cell carcinoma (HNSCC) has provided a unique opportunity to develop selective cancer treatment options. These efforts will require the establishment of relevant HNSCC models for preclinical testing. Here, we performed full exome and transcriptome sequencing of a large panel of HNSCC-derived cells from different anatomical locations and human papillomavirus (HPV) infection status. These cells exhibit typical mutations in TP53, FAT1, CDK2NA, CASP8, and NOTCH1, and copy number variations (CNVs) and mutations in PIK3CA, HRAS, and PTEN that reflect the widespread activation of the PI3K-mTOR pathway. SMAD4 alterations were observed that may explain the decreased tumor suppressive effect of TGF-β in HNSCC. Surprisingly, we identified HPV⁺ HNSCC cells harboring TP53 mutations, and documented aberrant TP53 expression in a subset of HPV⁺ HNSCC cases. This analysis also revealed that most HNSCC cells harbor multiple mutations and CNVs in epigenetic modifiers (e.g., EP300, CREBP, MLL1, MLL2, MLL3, KDM6A, and KDM6B) that may contribute to HNSCC initiation and progression. These genetically-defined experimental HNSCC cellular systems, together with the identification of novel actionable molecular targets, may now facilitate the pre-clinical evaluation of emerging therapeutic agents in tumors exhibiting each precise genomic alteration.