Brandon J. Manley

Tumor immune microenvironment characterization in clear cell renal cell carcinoma identifies prognostic and immunotherapeutically relevant messenger RNA signatures.

BackgroundTumor-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.ResultsWe 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.ConclusionsOur 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.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.

Genomic landscape and evolution of metastatic chromophobe renal cell carcinoma

Chromophobe renal cell carcinoma (chRCC) typically shows ~7 chromosome losses (1, 2, 6, 10, 13, 17, and 21) and ~31 exonic somatic mutations, yet carries ~5%-10% metastatic incidence. Since extensive chromosomal losses can generate proteotoxic stress and compromise cellular proliferation, it is intriguing how chRCC, a tumor with extensive chromosome losses and a low number of somatic mutations, can develop lethal metastases. Genomic features distinguishing metastatic from nonmetastatic chRCC are unknown. An integrated approach, including whole-genome sequencing (WGS), targeted ultradeep cancer gene sequencing, and chromosome analyses (FACETS, OncoScan, and FISH), was performed on 79 chRCC patients including 38 metastatic (M-chRCC) cases. We demonstrate that TP53 mutations (58%), PTEN mutations (24%), and imbalanced chromosome duplication (ICD, duplication of ≥ 3 chromosomes) (25%) were enriched in M-chRCC. Reconstruction of the subclonal composition of paired primary-metastatic chRCC tumors supports the role of TP53, PTEN, and ICD in metastatic evolution. Finally, the presence of these 3 genomic features in primary tumors of both The Cancer Genome Atlas kidney chromophobe (KICH) (n = 64) and M-chRCC (n = 35) cohorts was associated with worse survival. In summary, our study provides genomic insights into the metastatic progression of chRCC and identifies TP53 mutations, PTEN mutations, and ICD as high-risk features.

Integration of Recurrent Somatic Mutations with Clinical Outcomes: A Pooled Analysis of 1049 Patients with Clear Cell Renal Cell Carcinoma

BACKGROUND Analyses of associations between clinicopathologic outcomes and recurrent somatic mutations in clear cell renal cell carcinoma (ccRCC) have been limited to individual cohorts. OBJECTIVE To define clinicopathologic associations between specific mutations and ccRCC disease characteristics. DESIGN, SETTING, AND PARTICIPANTS DNA sequencing data were pooled from three collaborative genomic cohorts (n=754) and our institutional database (n=295). All patients had clinical data and identification of somatic mutations from their primary tumors. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Analysis of gene mutations for associations with maximal tumor size (linear regression) and pathologic stage (logistic regression). Cancer-specific survival (CSS) and recurrence-free survival (RFS) were calculated using competing risks methods. Analyses were adjusted for cohort site, and results were adjusted for multiple testing (q value). Relevant genes were used in multivariable models that included confounding variables and the validated Mayo Clinic Stage, Size, Grade, and Necrosis (SSIGN) score. RESULTS AND LIMITATIONS Association with tumor size was found for mutations in BAP1 (q=0.013). No mutations were found to be associated with stage after adjusted analysis. Mutations in BAP1 (q=0.004) and TP53 (q=0.001) were associated with decreased CSS in a multivariable model; only TP53 (q=0.005) remained significant when SSIGN score was included. SETD2 mutations (q=0.047) were associated with decreased RFS in multivariable models, including models with SSIGN score. CONCLUSIONS In >1000 patients with ccRCC, pooled analysis and multivariable modeling demonstrated that three mutated genes have statistically significant associations with poor clinical outcomes. This included the more commonly mutated BAP1 and SETD2 and the less frequently mutated TP53. After adjustment for clinical confounders, mutations of TP53 and SETD2 were associated with decreased CSS and RFS, respectively. PATIENT SUMMARY Using rigorous statistical methods, this study affirmed that certain mutations in clear cell renal cell carcinoma may portend inferior survival and an increased risk of recurrence.

Characterizing recurrent and lethal small renal masses in clear cell renal cell carcinoma using recurrent somatic mutations

INTRODUCTION Small renal masses (SRMs) with evidence of clear cell renal cell carcinoma (ccRCC) are understudied. Current algorithms for the management of SRMs include surgical resection, ablation, and active surveillance. We sought to identify genomic biomarkers that could potentially refine the management of ccRCC in SRMs, especially in patients being evaluated for active surveillance. METHODS We identified patients who had SRMs (4cm or less) at time of surgery, had sequencing performed on their primary tumor and had a diagnosis of ccRCC. Patients were selected from 3 publicly available cohorts, The Cancer Genome Atlas (n = 110), University of Tokyo (n = 37), The International Cancer Genome Consortium (n = 31), and from our own institutional prospective database (n = 25). Among this cohort we analyzed mutations present in at least 5% of tumors, assessing for the enrichment of mutations and progression-free survival using the composite endpoint of recurrence or death of disease. Analysis was adjusted for multiple testing. A Cox regression model was used to assess clinical variables with significant mutations. RESULTS In total, 203 patients were available for analysis. Median follow-up was 43.1 months among survivors. Mutations in VHL, PBRM1, SETD2, BAP1, KDM5C, and MTOR were present in more than 5% of tumors. Twenty-three patients (11.3%) had recurrence or died of their disease. Mutations in KDM5C were associated with inferior survival from either recurrence or death from disease, adjusted P 0.033. CONCLUSIONS We identified mutations in SRMs in ccRCC that are associated with recurrence and lethality. The strongest association was seen in those with KDM5C mutations. Use of these genomic biomarkers may improve stratification of patients with SRMs and for those who may be appropriate for active surveillance. Prospective evaluation of these markers is needed.

The difficulty in selecting patients for cytoreductive nephrectomy: An evaluation of previously described predictive models ☆

PURPOSE To externally evaluate a preoperative points system and a preoperative nomogram, both created to assess time to death after cytoreductive nephrectomy (CN). MATERIALS AND METHODS We identified 298 patients who underwent CN at our institution, a tertiary cancer center, between 1989 and 2015. To validate the points system, we compared reported overall survival (OS) for each criterion to observed OS in our cohort. To evaluate the nomogram, we prognosticated risk of death at 6 months after surgery for 280 patients with sufficient follow-up in our cohort and evaluated discrimination using area under the curve (AUC) and calibration. Decision curve analysis was performed to assess clinical utility of the nomogram. RESULTS Significant differences in OS were observed between patients with and without 5 of 7 criteria on univariate analysis: low albumin (P<0.0001), high lactate dehydrogenase (P = 0.002), liver metastasis (P = 0.004), retroperitoneal lymphadenopathy (P = 0.002), and supradiaphragmatic lymphadenopathy (P = 0.019). Discrimination from the preoperative model, predicting death within 6 months of surgery was lower in our cohort (AUC = 0.65, 95% CI: 0.52-0.79) than the original publication (AUC = 0.76). Decision curve analysis demonstrated little benefit for applicability. CONCLUSIONS Five previously defined risk factors are predictive of decreased OS after CN in our cohort. We found lower discrimination using the preoperative model and minimal clinical utility according to decision analysis in our study cohort. These findings suggest the need for improved models to aid patient stratification and consequent treatment choice.

Genomic alterations as predictors of survival among patients within a combined cohort with clear cell renal cell carcinoma undergoing cytoreductive nephrectomy

PURPOSE To establish prognostic genomic biomarkers for patients with metastatic clear cell renal cell carcinoma (ccRCC). MATERIALS AND METHODS We identified 60 patients who presented with metastatic ccRCC at our institution between 2001 and 2015 and had genomic sequencing on their primary tumor. We pooled these patients with 107 other patients with the same inclusion criteria from three well-known public databases. Five commonly mutated genes were chosen for analysis: VHL, PBRM1, BAP1, SETD2, and KDM5C. Overall survival (OS) was estimated using the Kaplan-Meier method and the log-rank test was used for comparisons between groups. RESULTS Median OS in the cohort was 2.5 years. Higher Fuhrman grade was associated with decreased median OS (P<0.001). Mutations in SETD2 (P = 0.027) and KDM5C (P = 0.019) were associated with reduced risk of death (hazard ratio [HR] = 0.58 [95% CI: 0.35-0.94] and HR = 0.43 [95% CI: 0.22-0.85], respectively). BAP1 mutations (P = 0.008) were associated with increased risk of death (HR = 1.81 [95% CI: 1.16-2.83]). There were significantly more female patients with a BAP1 mutation than females in the overall cohort (P = 0.001). CONCLUSIONS Mutations in BAP1 negatively affected OS, whereas SETD2 and KDM5C mutations were associated with prolonged OS in our pooled cohort of 167 patients with metastatic ccRCC. Our results expand upon efforts at understanding genomic biomarkers in localized disease. Those efforts set the stage for our novel investigation examining associations of select recurrent somatic mutations in stage IV patients with ccRCC.

Comparative Genomic Profiling of Matched Primary and Metastatic Tumors in Renal Cell Carcinoma

BACKGROUND Next-generation sequencing (NGS) studies of matched pairs of primary and metastatic tumors in renal cell carcinoma (RCC) have been limited to small cohorts. OBJECTIVE To evaluate the discordance in somatic mutations between matched primary and metastatic RCC tumors. DESIGN, SETTING, AND PARTICIPANTS Primary tumor (P), metastasis (M), and germline DNA from 60 patients with RCC was subjected to NGS with a targeted exon capture-based assay of 341 cancer-associated genes. Somatic mutations were called using a validated pipeline. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Mutations were classified as shared (S) or private (Pr) in relation to each other within individual P-M pairs. The concordance score was calculated as (S-Pr)/(S+Pr). To calculate enrichment of Pr/S mutations for a particular gene, we calculated a two-sided p value from a binomial model for each gene with at least ten somatic mutation events, and also implemented a separate permutation test procedure. We adjusted p values for multiple hypothesis testing using the Benjamini-Hochberg procedure. The mutation discordance was calculated using Mann-Whitney U tests according to gene mutations or metastatic sites. RESULTS AND LIMITATIONS Twenty-one pairs (35%) showed Pr mutations in both P and M samples. Of the remaining 39 pairs (65%), 14 (23%) had Pr mutations specific to P samples, 12 (20%) had Pr mutations to M samples, and 13 (22%) had identical somatic mutations. No individual gene mutation was preferentially enriched in either P or M samples. P-M pairs with SETD2 mutations demonstrated higher discordance than pairs with wild-type SETD2. We observed that patients who received therapy before sampling of the P or M tissue had higher concordance of mutations for P-M pairs than patients who did not (Mann-Whitney p=0.088). CONCLUSIONS Our data show mutation discordance within matched P-M RCC tumor pairs. As most contemporary precision medicine trials do not differentiate mutations detected in P and M tumors, the prognostic and predictive value of mutations in P versus M tumors warrants further investigation. PATIENT SUMMARY In this study we evaluated the concordance of mutations between matched primary and metastatic tumors for 60 kidney cancer patients using a panel of 341 cancer genes. Forty-seven patients carried nonidentical cancer gene mutations within their matched primary-metastatic pair. The mutation profile of the primary tumor alone could compromise precision in selecting effective targeted therapies and result in suboptimal clinical outcomes.

Characterization and Impact of TERT Promoter Region Mutations on Clinical Outcome in Renal Cell Carcinoma

BACKGROUND Mutations in the promoter region of the TERT gene have been detected in a variety of cancers. These mutations can potentially lead to unlimited cell divisions and result in poor clinical prognosis. OBJECTIVE To determine the role and relevance of TERT promoter region mutations in both clear cell (ccRCC) and non-clear cell (nccRCC) renal cell carcinoma using ultra-deep and whole-genome sequencing methods on primary tumor samples. DESIGN, SETTING, AND PARTICIPANTS DNA from 281 kidney tumors (147 ccRCC and 134 nccRCC) was sequenced between 2013 and 2015, and clinical outcomes for these patients from a single institution were retrospectively analyzed. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Differences in patient characteristics and mutational status were tested using Fishers exact test for categorical variables and the Wilcoxon rank sum test for continuous variables. Survival times were estimated using the Kaplan-Meier method and differences were tested using the log-rank test. RESULTS AND LIMITATIONS TERT mutations occurred in 12.2% of ccRCC and 10.4% of nccRCC cases. In >80% of the cases, mutations were located at C228T and were found to co-occur only rarely with other relevant RCC driver genes. The median follow-up among survivors overall was 2.5 yr (range 0.1-18.3). TERT promoter mutations were significantly associated with cancer-specific survival in ccRCC (hazard ratio 2.68, 95% confidence interval 1.19-6.01; p=0.013). In nccRCC, TERT mutations were significantly associated with larger tumors and metastatic development. Assessment of further relevant clinical associations was precluded in the nccRCC group by the heterogeneous and small sample size. CONCLUSIONS Our data suggests that TERT mutational status reflects a distinct pathogenesis with an aggressive disease course in RCC. Stratifying patients with this unique tumorigenesis that leads to poor clinical outcomes could be a putative target for novel therapeutics. PATIENT SUMMARY We show a previously unrecognized frequency of TERT promoter mutations in both clear cell and non-clear cell renal cell carcinoma. TERT promoter mutations were associated with some worse outcomes in patients with clear cell renal cell carcinoma.

Cystic renal cell carcinoma: a report on outcomes of surgery and active surveillance in patients retrospectively identified on pretreatment imaging

Purpose: We evaluated the outcomes of surgical intervention and active surveillance in patients diagnosed with cystic renal cell carcinoma at our hypothesized radiological cutoff of greater than 50% cystic. Materials and Methods: We identified all 430 patients with a pathologically confirmed cystic renal mass that fit our criteria from 2000 to 2015. The 292 patients with a lack of computerized tomography, tumors less than 50% cystic on imaging, multifocal tumors and prior renal cell carcinoma were excluded from study. Patients were stratified into benign or malignant subgroups, and radiological, clinicopathological and oncologic features were determined. Univariate and multivariate associations between clinicoradiological parameters in each group were analyzed. We similarly reviewed the records of a separate cohort of patients treated with active surveillance for cystic renal cell carcinoma. Results: Of the 138 identified cases of cystic renal cell carcinoma 102 (73.9%) were renal cell carcinoma and 36 (26.1%) were benign masses. Of the tumors 77.5% were Fuhrman grade 1‐2, 83.4% were stage pT2 or less and 65.9% showed clear cell histology. On univariate analysis male gender, a solid component and increasing Bosniak classification were significant for malignancy. In a separate cohort we identified 38 patients on active surveillance. The growth rate was 1.0 mm per year overall and 2.3 mm per year for the solid component. At a median followup of more than 4 years in all cohorts there was no evidence of recurrence or metastasis of cystic renal cell carcinoma. Conclusions: Patients with unifocal cystic renal cell carcinoma evaluated using a standardized radiological threshold of greater than 50% cystic had an excellent prognosis on active surveillance and after surgical resection.

Genomic Alterations and Outcomes with VEGF-Targeted Therapy in Patients with Clear Cell Renal Cell Carcinoma

Background: Mutations in VHL, PBRM1, SETD2, BAP1, and KDM5C are common in clear cell renal cell carcinoma (ccRCC), and presence of certain mutations has been associated with outcomes in patients with non-metastatic disease. Limited information is available regarding the correlation between genomic alterations and outcomes in patients with metastatic disease, including response to VEGF-targeted therapy. Objective: To explore correlations between mutational profiles and cancer-specific outcomes, including response to standard VEGF-targeted agents, in patients with metastatic cc RCC. Methods: A retrospective review of 105 patients with metastatic ccRCC who had received systemic therapy and had targeted next-generation sequencing of tumors was conducted. Genomic alterations were correlated to outcomes, including overall survival and time to treatment failure to VEGF-targeted therapy. Results: The most frequent mutations were detected in VHL (83%), PBRM1 (51%), SETD2 (35%), BAP1 (24%), KDM5C (16%), and TERT (14%). Time to treatment failure with VEGF-targeted therapy differed significantly by PBRM1 mutation status (p = 0.01, median 12.0 months for MT versus 6.9 months for WT) and BAP1 mutation status (p = 0.01, median 6.4 months for MT versus 11.0 months for WT). Shorter overall survival was associated with TERT mutations (p = 0.03, median 29.6 months for MT versus 52.6 months for WT) or BAP1 mutations (p = 0.02, median 28.7 months for MT versus not reached for WT). Conclusions: Genomic alterations in ccRCC tumors have prognostic implications in patients with metastatic disease. BAP1 and TERT promoter mutations may be present in higher frequency than previously thought, and based on this data, deserve further study for their association with poor prognosis.