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 Immunopeptidome Analysis Reveals Recurrent Posttranslational Splicing of Cancer- and Immune-Associated Genes.

Posttranslational spliced peptides (PTSPs) are a unique class of peptides that have been found to be presented by HLA class-I molecules in cancer. 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. Furthermore, the role of the HLA class-II pathway in PTSP presentation has been studied only in diabetes. 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 class-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.

Post-translational modifications reshape the antigenic landscape of the MHC I immunopeptidome in tumors.

Post-translational modification (PTM) of antigens provides an additional source of specificities targeted by immune responses to tumors or pathogens, but identifying antigen PTMs and assessing their role in shaping the immunopeptidome is challenging. Here we describe the Protein Modification Integrated Search Engine (PROMISE), an antigen discovery pipeline that enables the analysis of 29 different PTM combinations from multiple clinical cohorts and cell lines. We expanded the antigen landscape, uncovering human leukocyte antigen class I binding motifs defined by specific PTMs with haplotype-specific binding preferences and revealing disease-specific modified targets, including thousands of new cancer-specific antigens that can be shared between patients and across cancer types. Furthermore, we uncovered a subset of modified peptides that are specific to cancer tissue and driven by post-translational changes that occurred in the tumor proteome. Our findings highlight principles of PTM-driven antigenicity, which may have broad implications for T cell-mediated therapies in cancer and beyond.

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.

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.

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 immunotherapy1,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 alone6,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.

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.

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.

Immunoproteasome expression is associated with better prognosis and response to checkpoint therapies in melanoma.

Predicting the outcome of immunotherapy treatment in melanoma patients is challenging. Alterations in genes involved in antigen presentation and the interferon gamma (IFNγ) pathway play an important role in the immune response to tumors. We describe here that the overexpression of PSMB8 and PSMB9, two major components of the immunoproteasome, is predictive of better survival and improved response to immune-checkpoint inhibitors of melanoma patients. We study the mechanism underlying this connection by analyzing the antigenic peptide repertoire of cells that overexpress these subunits using HLA peptidomics. We find a higher response of patient-matched tumor infiltrating lymphocytes against antigens diferentially presented after immunoproteasome overexpression. Importantly, we find that PSMB8 and PSMB9 expression levels are much stronger predictors of melanoma patients' immune response to checkpoint inhibitors than the tumors' mutational burden. These results suggest that PSMB8 and PSMB9 expression levels can serve as important biomarkers for stratifying melanoma patients for immune-checkpoint treatment.

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.

Author Correction: Combined inhibition of MEK and nuclear ERK translocation has synergistic antitumor activity in melanoma 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.