Highlights d Novel mouse system to uncouple tumor mutational load and tumor heterogeneity d Lower tumor heterogeneity leads to decreased tumor growth because of immune rejection d Both clone numbers and their genetic diversity mediate tumor growth and rejection d Tumor heterogeneity is linked to patient survival and checkpoint blockade response

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.

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 and, 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. 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. .

The tumor microenvironment (TME) is a complex mixture of cell types whose interactions affect tumor growth and clinical outcome. To discover such interactions, we developed CODEFACS (COnfident DEconvolution For All Cell Subsets), a tool deconvolving cell type–specific gene expression in each sample from bulk expression, and LIRICS (Ligand–Receptor Interactions between Cell Subsets), a statistical framework prioritizing clinically relevant ligand–receptor interactions between cell types from the deconvolved data. We first demonstrate the superiority of CODEFACS versus the state-of-the-art deconvolution method CIBERSORTx. Second, analyzing The Cancer Genome Atlas, we uncover cell type–specific ligand–receptor interactions uniquely associated with mismatch-repair deficiency across different cancer types, providing additional insights into their enhanced sensitivity to anti–programmed cell death protein 1 (PD-1) therapy compared with other tumors with high neoantigen burden. Finally, we identify a subset of cell type–specific ligand–receptor interactions in the melanoma TME that stratify survival of patients receiving anti–PD-1 therapy better than some recently published bulk transcriptomics-based methods. Significance: This work presents two new computational methods that can deconvolve a large collection of bulk tumor gene expression profiles into their respective cell type–specific gene expression profiles and identify cell type–specific ligand–receptor interactions predictive of response to immune-checkpoint blockade therapy. This article is highlighted in the In This Issue feature, p. 873