Geospatial characterization of immune cell distributions and dynamics across the microenvironment in clear cell renal cell carcinoma.

INTRODUCTION: In clear cell renal cell carcinoma (ccRCC), tumor-associated macrophage (TAM) induction of CD8+T cells into a terminally exhausted state has been implicated as a major mechanism of immunotherapy resistance, but a deeper biological understanding is necessary. METHODS: Primary ccRCC tumor samples were obtained from 97 patients between 2004 and 2018. Multiplex immunofluorescence using lymphoid and myeloid markers was performed in seven regions of interest per patient across three predefined zones, and geospatial analysis was performed using Ripley's K analysis, a methodology adapted from ecology. RESULTS: Clustering of CD163+M2 like TAMs into the stromal compartment at the tumor-stroma interface was associated with worse clinical stage (tumor/CD163+nK(75): stage I/II: 4.4 (IQR -0.5 to 5.1); stage III: 1.4 (IQR -0.3 to 3.5); stage IV: 0.6 (IQR -2.1 to 2.1); p=0.04 between stage I/II and stage IV), and worse overall survival (OS) and cancer-specific survival (CSS) (tumor/CD163+nK(75): median OS-hi=149 months, lo=86 months, false-discovery rate (FDR)-adj. Cox p<0.001; median CSS-hi=174 months, lo=85 months; FDR-adj. Cox p<0.001). An RNA-seq differential gene expression score was developed using this geospatial metric, and was externally validated in multiple independent cohorts of patients with ccRCC including: TCGA KIRC, and the IMmotion151, IMmotion150, and JAVELIN Renal 101 clinical trials. In addition, this CD163+ geospatial pattern was found to be associated with a higher TIM-3+ proportion of CD8+T cells, indicative of terminal exhaustion (tumor-core: 0.07 (IQR 0.04-0.14) vs 0.40 (IQR 0.15-0.66), p=0.05). CONCLUSIONS: Geospatial clustering of CD163+M2 like TAMs into the stromal compartment at the tumor-stromal interface was associated with poor clinical outcomes and CD8+T cell terminal exhaustion.

Correlating Immune Cell Infiltration Patterns with Recurrent Somatic Mutations in Advanced Clear Cell Renal Cell Carcinoma.

BACKGROUND: Clear cell renal cell carcinoma (ccRCC) tumors have low frequencies of genetic alterations compared with other malignancies, but very high levels of immune cell infiltration and favorable response rates to immunotherapy. Currently, the interplay between specific ccRCC somatic mutations and immune infiltration pattern is unclear. OBJECTIVE: To analyze the associations between common ccRCC somatic mutations and immune cell infiltration patterns within the tumor immune microenvironment (TIME). DESIGN, SETTING, AND PARTICIPANTS: The study included tumor samples (24 primary and 24 metastatic) from 48 patients with stage IV ccRCC. Targeted sequencing was performed for well-characterized recurrent somatic mutations in ccRCC, with the analysis focusing on the six most common ones: VHL, BAP1, PBRM1, SETD2, TP53, and KDM5C. For each sample, multiplex immunofluorescence (IF) was performed in lymphoid and myeloid panels, for seven regions of interest in three zones (tumor core, stroma, and tumor-stroma interface). IF-derived cellular densities were compared across patients, stratified by their somatic mutation status, using a linear mixed-model analysis. External validation was pursued using RNA-seq enrichment scoring from three large external data sources. RESULTS AND LIMITATIONS: Tumors with SETD2 mutations demonstrated significantly decreased levels of FOXP3+ T cells in the tumor core, stroma, and tumor-stroma interface. PBRM1 mutations were associated with decreased FOXP3+ T cells in the tumor core. Primary KDM5C mutations were associated with significantly increased CD206+ macrophage tumor infiltration in the tumor core. A computational method estimating immune cell types in the TIME using bulk RNA-seq data, xCell scoring, failed to validate associations from the IF analysis in large external data sets. A major limitation of the study is the relatively small patient population studied. CONCLUSIONS: This study provides evidence that common somatic mutations in ccRCC, such as SETD2, PBRM1, and KDM5C, are associated with distinct immune infiltration patterns within the TIME. PATIENT SUMMARY: In this study, we analyzed tumor samples from patients with metastatic kidney cancer to determine whether common genetic mutations that arise from the cancer cells are associated with the density of immune cells found within those tumors. We found several distinct immune cell patterns that were associated with specific genetic mutations. These findings provide insight into the interaction between cancer genetics and the immune system in kidney cancer.

Spatial clustering of CD68+ tumor associated macrophages with tumor cells is associated with worse overall survival in metastatic clear cell renal cell carcinoma.

Immune infiltration is typically quantified using cellular density, not accounting for cellular clustering. Tumor-associated macrophages (TAM) activate oncogenic signaling through paracrine interactions with tumor cells, which may be better reflected by local cellular clustering than global density metrics. Using multiplex immunohistochemistry and digital pathologic analysis we quantified cellular density and cellular clustering for myeloid cell markers in 129 regions of interest from 55 samples from 35 patients with metastatic ccRCC. CD68+ cells were found to be clustered with tumor cells and dispersed from stromal cells, while CD163+ and CD206+ cells were found to be clustered with stromal cells and dispersed from tumor cells. CD68+ density was not associated with OS, while high tumor/CD68+ cell clustering was associated with significantly worse OS. These novel findings would not have been identified if immune infiltrate was assessed using cellular density alone, highlighting the importance of including spatial analysis in studies of immune cell infiltration of tumors. Significance: Increased clustering of CD68+ TAMs and tumor cells was associated with worse overall survival for patients with metastatic ccRCC. This effect would not have been identified if immune infiltrate was assessed using cell density alone, highlighting the importance of including spatial analysis in studies of immune cell infiltration of tumors.