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.

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

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.

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.