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.

A pilot study to evaluate tissue- and plasma-based DNA driver mutations in a cohort of patients with pancreatic intraductal papillary mucinous neoplasms.

Intraductal papillary mucinous neoplasms (IPMNs) are precursor lesions to pancreatic ductal adenocarcinoma that are challenging to manage due to limited imaging, cytologic, and molecular markers that accurately classify lesions, grade of dysplasia, or focus of invasion preoperatively. The objective of this pilot study was to determine the frequency and type of DNA mutations in a cohort of surgically resected, pathologically confirmed IPMN, and to determine if concordant mutations are detectable in paired pretreatment plasma samples. Formalin-fixed paraffin-embedded (FFPE) tissue from 46 surgically resected IPMNs (31 low-grade, 15 high-grade) and paired plasma from a subset of 15 IPMN cases (10 low-grade, 5 high-grade) were subjected to targeted mutation analysis using a QIAseq Targeted DNA Custom Panel. Common driver mutations were detected in FFPE from 44 of 46 (95.6%) IPMN cases spanning all grades; the most common DNA mutations included: KRAS (80%), RNF43 (24%), and GNAS (43%). Of note, we observed a significant increase in the frequency of RNF43 mutations from low-grade to high-grade IPMNs associated or concomitant with invasive carcinoma (trend test, P = 0.01). Among the subset of cases with paired plasma, driver mutations identified in the IPMNs were not detected in circulation. Overall, our results indicate that mutational burden for IPMNs is a common occurrence, even in low-grade IPMNs. Furthermore, although blood-based biopsies are an attractive, noninvasive method for detecting somatic DNA mutations, the QIAseq panel was not sensitive enough to detect driver mutations that existed in IPMN tissue using paired plasma in the volume we were able to retrieve for this retrospective study.

Quantification of T- and B-cell Immune Receptor Distribution Diversity Characterizes Immune Cell Infiltration and Lymphocyte Heterogeneity in Clear Cell Renal Cell Carcinoma.

Immune-modulating systemic therapies are often used to treat advanced cancer such as metastatic clear cell renal cell carcinoma (ccRCC). Used alone, sequence-based biomarkers neither accurately capture patient dynamics nor the tumor immune microenvironment. To better understand the tumor ecology of this immune microenvironment, we quantified tumor infiltration across three distinct ccRCC patient tumor cohorts using complementarity determining region-3 (CDR3) sequence recovery counts in tumor-infiltrating lymphocytes and a generalized diversity index (GDI) for CDR3 sequence distributions. GDI can be understood as a curve over a continuum of diversity scales that allows sensitive characterization of distributions to capture sample richness, evenness, and subsampling uncertainty, along with other important metrics that characterize tumor heterogeneity. For example, richness quantified the total unique sequence count, while evenness quantified similarities across sequence frequencies. Significant differences in receptor sequence diversity across gender and race revealed that patients with larger and more clinically aggressive tumors had increased richness of recovered tumoral CDR3 sequences, specifically in those from T-cell receptor alpha and B-cell immunoglobulin lambda light chain. The GDI inflection point (IP) allowed for a novel and robust measure of distribution evenness. High IP values were associated with improved overall survival, suggesting that normal-like sequence distributions lead to better outcomes. These results propose a new quantitative tool that can be used to better characterize patient-specific differences related to immune cell infiltration, and to identify unique characteristics of tumor-infiltrating lymphocyte heterogeneity in ccRCC and other malignancies. SIGNIFICANCE: Assessment of tumor-infiltrating T-cell and B-cell diversity in renal cell carcinoma advances the understanding of tumor-immune system interactions, linking tumor immune ecology with tumor burden, aggressiveness, and patient survival. See related commentary by Krishna and Hakimi, p. 764.

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.

Geospatial Cellular Distribution of Cancer-Associated Fibroblasts Significantly Impacts Clinical Outcomes in Metastatic Clear Cell Renal Cell Carcinoma.

Cancer-associated fibroblasts (CAF) are highly prevalent cells in the tumor microenvironment in clear cell renal cell carcinoma (ccRCC). CAFs exhibit a pro-tumor effect in vitro and have been implicated in tumor cell proliferation, metastasis, and treatment resistance. Our objective is to analyze the geospatial distribution of CAFs with proliferating and apoptotic tumor cells in the ccRCC tumor microenvironment and determine associations with survival and systemic treatment. Pre-treatment primary tumor samples were collected from 96 patients with metastatic ccRCC. Three adjacent slices were obtained from 2 tumor-core regions of interest (ROI) per patient, and immunohistochemistry (IHC) staining was performed for αSMA, Ki-67, and caspase-3 to detect CAFs, proliferating cells, and apoptotic cells, respectively. H-scores and cellular density were generated for each marker. ROIs were aligned, and spatial point patterns were generated, which were then used to perform spatial analyses using a normalized Ripley's K function at a radius of 25 µm (nK(25)). The survival analyses used an optimal cut-point method, maximizing the log-rank statistic, to stratify the IHC-derived metrics into high and low groups. Multivariable Cox regression analyses were performed accounting for age and International Metastatic RCC Database Consortium (IMDC) risk category. Survival outcomes included overall survival (OS) from the date of diagnosis, OS from the date of immunotherapy initiation (OS-IT), and OS from the date of targeted therapy initiation (OS-TT). Therapy resistance was defined as progression-free survival (PFS) <6 months, and therapy response was defined as PFS >9 months. CAFs exhibited higher cellular clustering with Ki-67+ cells than with caspase-3+ cells (nK(25): Ki-67 1.19; caspase-3 1.05; p = 0.04). The median nearest neighbor (NN) distance from CAFs to Ki-67+ cells was shorter compared to caspase-3+ cells (15 µm vs. 37 µm, respectively; p < 0.001). Multivariable Cox regression analyses demonstrated that both high Ki-67+ density and H-score were associated with worse OS, OS-IT, and OS-TT. Regarding αSMA+CAFs, only a high H-score was associated with worse OS, OS-IT, and OS-TT. For caspase-3+, high H-score and density were associated with worse OS and OS-TT. Patients whose tumors were resistant to targeted therapy (TT) had higher Ki-67 density and H-scores than those who had TT responses. Overall, this ex vivo geospatial analysis of CAF distribution suggests that close proximity clustering of tumor cells and CAFs potentiates tumor cell proliferation, resulting in worse OS and resistance to TT in metastatic ccRCC.

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.