Tumor-infiltrating immune cells have been linked to prognosis and response to immunotherapy; however, the levels of distinct immune cell subsets and the signals that draw them into a tumor, such as the expression of antigen presenting machinery genes, remain poorly characterized. Here, we employ a gene expression-based computational method to profile the infiltration levels of 24 immune cell populations in 19 cancer types. We compare cancer types using an immune infiltration score and a T cell infiltration score and find that clear cell renal cell carcinoma (ccRCC) is among the highest for both scores. Using immune infiltration profiles as well as transcriptomic and proteomic datasets, we characterize three groups of ccRCC tumors: T cell enriched, heterogeneously infiltrated, and non-infiltrated. We observe that the immunogenicity of ccRCC tumors cannot be explained by mutation load or neo-antigen load, but is highly correlated with MHC class I antigen presenting machinery expression (APM). We explore the prognostic value of distinct T cell subsets and show in two cohorts that Th17 cells and CD8+ T/Treg ratio are associated with improved survival, whereas Th2 cells and Tregs are associated with negative outcomes. Investigation of the association of immune infiltration patterns with the subclonal architecture of tumors shows that both APM and T cell levels are negatively associated with subclone number. Our analysis sheds light on the immune infiltration patterns of 19 human cancers and unravels mRNA signatures with prognostic utility and immunotherapeutic biomarker potential in ccRCC.

Purpose: To investigate the impact of newly identified chromosome 3p21 epigenetic tumor suppressors PBRM1, SETD2, and BAP1 on cancer-specific survival (CSS) of 609 patients with clear cell renal cell carcinoma (ccRCC) from 2 distinct cohorts. Experimental Design: Select sequencing on 3p tumor suppressors of 188 patients who underwent resection of primary ccRCC at the Memorial Sloan-Kettering Cancer Center (MSKCC) was conducted to interrogate the genotype–phenotype associations. These findings were compared with analyses of the genomic and clinical dataset from our nonoverlapping The Cancer Genome Atlas (TCGA) cohort of 421 patients with primary ccRCC. Results: 3p21 tumor suppressors are frequently mutated in both the MSKCC (PBRM1, 30.3%; SETD2, 7.4%; BAP1, 6.4%) and the TCGA (PBRM1, 33.5%; SETD2, 11.6%; BAP1, 9.7%) cohorts. BAP1 mutations are associated with worse CSS in both cohorts [MSKCC, P = 0.002; HR 7.71; 95% confidence interval (CI)2.08–28.6; TCGA, P = 0.002; HR 2.21; 95% CI 1.35–3.63]. SETD2 are associated with worse CSS in the TCGA cohort (P = 0.036; HR 1.68; 95% CI 1.04–2.73). On the contrary, PBRM1 mutations, the second most common gene mutations of ccRCC, have no impact on CSS. Conclusion: The chromosome 3p21 locus harbors 3 frequently mutated ccRCC tumor suppressor genes. BAP1 and SETD2 mutations (6%–12%) are associated with worse CSS, suggesting their roles in disease progression. PBRM1 mutations (30%–34%) do not impact CSS, implicating its principal role in the tumor initiation. Future efforts should focus on therapeutic interventions and further clinical, pathologic, and molecular interrogation of this novel class of tumor suppressors. Clin Cancer Res; 19(12); 3259–67. ©2013 AACR.

PBRM1 is the second most commonly mutated gene after VHL in clear cell renal cell carcinoma (ccRCC). However, the biological consequences of PBRM1 mutations for kidney tumorigenesis are unknown. Here, we find that kidney-specific deletion of Vhl and Pbrm1, but not either gene alone, results in bilateral, multifocal, transplantable clear cell kidney cancers. PBRM1 loss amplified the transcriptional outputs of HIF1 and STAT3 incurred by Vhl deficiency. Analysis of mouse and human ccRCC revealed convergence on mTOR activation, representing the third driver event after genetic inactivation of VHL and PBRM1. Our study reports a physiological preclinical ccRCC mouse model that recapitulates somatic mutations in human ccRCC and provides mechanistic and therapeutic insights into PBRM1 mutated subtypes of human ccRCC.

Purpose: Rapalogs are allosteric mTOR inhibitors and approved agents for advanced kidney cancer. Reports of clonal heterogeneity in this disease challenge the concept of targeted monotherapy, yet a small subset of patients derives extended benefit. Our aim was to analyze such outliers and explore the genomic background of extreme rapalog sensitivity in the context of intratumor heterogeneity. Experimental Design: We analyzed archived tumor tissue of 5 patients with renal cell carcinoma, who previously achieved durable disease control with rapalogs (median duration, 28 months). DNA was extracted from spatially separate areas of primary tumors and metastases. Custom target capture and ultradeep sequencing was used to identify alterations across 230 target genes. Whole-exome sequence analysis was added to investigate genes beyond this original target list. Results: Five long-term responders contributed 14 specimens to explore clonal heterogeneity. Genomic alterations with activating effect on mTOR signaling were detected in 11 of 14 specimens, offering plausible explanation for exceptional treatment response through alterations in two genes (TSC1 and MTOR). In two subjects, distinct yet functionally convergent alterations activated the mTOR pathway in spatially separate sites. In 1 patient, concurrent genomic events occurred in two separate pathway components across different tumor regions. Conclusions: Analysis of outlier cases can facilitate identification of potential biomarkers for targeted agents, and we implicate two genes as candidates for further study in this class of drugs. The previously reported phenomenon of clonal convergence can occur within a targetable pathway which might have implications for biomarker development beyond this disease and this class of agents. Clin Cancer Res; 20(7); 1955–64. ©2014 AACR.