Histologic-Based Tumor-Associated Immune Cells Status in Clear Cell Renal Cell Carcinoma Correlates with Gene Signatures Related to Cancer Immunity and Clinical Outcomes

Abstract

The three-tier immunophenotype (desert, excluded, and inflamed) and the four-tier immunophenotype (cold, immunosuppressed, excluded, and hot) have been linked to prognosis and immunotherapy response. This study aims to evaluate whether immunophenotypes of clear cell renal cell carcinoma, identified on hematoxylin and eosin-stained slides, correlate with gene expression signatures related to cancer immunity, and clinical outcomes. We evaluated tumor-associated immune cells (TAICs) status using three methodologies: three-tier immunophenotype based on the location of TAICs, four-tier immunophenotype considering both the location and degree of TAICs and inflammation score focusing only on the degree of TAICs, using a localized clear cell renal cell carcinoma cohort (n = 436) and The Cancer Genome Atlas (TCGA)-KIRC cohort (n = 162). We evaluated the association of the TAICs status assessed by three methodologies with CD8 and PD-L1 immunohistochemistry and immune gene expression signatures by TCGA RNA-sequencing data. All three methodologies correlated with immunohistochemical and immune gene expression signatures. The inflammation score and the four-tier immunophenotype showed similarly higher accuracy in predicting recurrence-free survival and overall survival compared to the three-tier immunophenotype. In conclusion, a simple histologic assessment of TIACs may predict clinical outcomes and immunotherapy responses.

1. Introduction

Recently, the prognosis of metastatic renal cell carcinoma (RCC) has improved by the efficacy of immune checkpoint inhibitors (ICIs) such as agents directed against cytotoxic T-lymphocyte antigen 4 (CTLA-4), programmed death-1 (PD-1), and programmed death ligand-1 (PD-L1) [1] in addition to the anti-vascular endothelial growth factor (VEGR) targeted therapy such as tyrosine kinase inhibitors (TKIs) [2]. Comprehensive genomic investigations have provided a biology-based tumor immune microenvironment for treatment selection using genomic and transcriptomic analysis [3,4,5]. Several clinical trials on clear cell RCC (ccRCC) to evaluate their potential predictive value of TKIs, ICIs alone, or in combination, for patients with metastatic RCC revealed the association of gene expression signatures related to angiogenesis, effector T-cell, and myeloid inflammation with clinical outcome [6,7].

Since multiregional molecular genetic studies have largely been performed without consideration of morphology, the investigation regarding the correlation between histologic immune status and underlying genomes has been limited [8]. Genetic analysis and multiplex immunofluorescence/immunohistochemical staining remain challenging in routine clinical practice. Thus, a simple histologic-based assessment of the tumor immune microenvironment is urgently needed. A standardized methodology to assess tumor-infiltrating lymphocytes (TILs) has been proposed in solid tumors to assess the immune response to tumors [9]. On the other hand, recent evidence suggests that tumor immunity has been categorized into three immunophenotypes based on the location of TILs as desert, non-inflamed; excluded, peritumoral immune infiltration; and inflamed, intratumoral immune infiltration [10,11]. However, the prognostic and predictive value of TILs for RCCs remains under investigation.

Regarding the histological features of ccRCC, although the most common architectural pattern was compact small nests with an extensive vascular network, various architectural patterns have been reported [12,13]. Our research group has recently established a vascularity-based architectural classification for ccRCC, which significantly correlated with angiogenesis and immune gene signatures [14]. Although we also confirmed that the immunophenotype, such as desert, excluded, and inflamed, was highly associated with the vascularity-based architectural classification, the degree of inflammation was variable in both the excluded and inflamed types. Recent advances of immune contexture, including the novel concept of a combination of immune variables, such as the nature, density, immune functional orientation, and distribution of immune cells within the tumor, has led to four classifications (cold, immunosuppressed, excluded, and hot) [15]. At present, there is little evidence regarding which histologic-based assessment of immunophenotype for ccRCC correlates with the underlying mechanism of cancer immunity.

The current study aims to evaluate whether immunophenotypes of ccRCC, identified on hematoxylin and eosin (H&E)-stained slides, correlated with the gene expression signature related to cancer immunity, and clinical outcomes. In the present study, we first validated the histologic-based tumor-associated immune cells (TAICs) status using three methodologies: three-tier immunophenotype based on the location of TAICs, four-tier immunophenotype considering both the location and degree of TAICs, and inflammation score focusing only on the TAIC degree. Then, we correlated findings with CD8 and PD-L1 immunohistochemical expression and gene expression signatures related to cancer immunity by RNA-sequencing data available from the Cancer Genome Atlas (TCGA) [16,17]. Furthermore, we compared the correlation and prognostic accuracy among three methods and evaluated the association of histologic-based TAICs with pathological prognostic factors.

2. Materials and Methods

2.1. Data Collection

In the present study, we used the same two cohorts as previously reported [14]; 436 cases with localized ccRCC (cT1-4N0-1M0) as the principal cohort and 162 ccRCC cases from the TCGA-KIRC cohort. Cases with presurgical treatment with TKIs or ICIs were not included in the principal cohort. The study of the principal cohort was approved by the Institutional Review Board (No. 2018109 (28/11/2018) and No. 2020222 (09/12/2020)) per the Declaration of Helsinki. TCGA whole slide images were accessed via the Cancer Digital Slide Archive [18]. We retrieved information regarding the tumor stage assessed by the 2017 TNM staging system [19], World Health Organization (WHO)/International Society of Urological Pathology (ISUP) nucleolar grade, and pathological prognostic factors, such as sarcomatoid/rhabdoid components and tumor-specific necrosis [20], from our institutional ccRCC database as previously described [14,21,22].

The study design is shown in Figure 1.

2.2. Histological Evaluation of Tumor Immune Microenvironment

To assess the tumor immune microenvironment, TAICs including both mononuclear cells and granulocytes were evaluated using whole H&E-stained slides. One representative slide including the highest-grade area (an average of five slides containing tumors per case) and one representative whole slide image was assessed in the principal cohort and TCGA cohort, respectively. The TAICs status was evaluated using three methods (shown in Table 1) by a genitourinary pathologist (C.O.) blinded to clinical outcomes. The three-tier immunophenotype was categorized based on the location of TAICs regardless of the degree of TAICs as follows: desert, no TAICs; excluded, peritumoral TAICs; and inflamed, intratumoral TAICs, as previously described [14]. The four-tier immunophenotype was categorized considering both the location and degree of TAICs as follows: cold, no TAICs; immunosuppressed, focal or low TAICs, regardless of the TAICs location; excluded, diffuse or high peritumoral TAICs; and hot, diffuse or high intratumoral TAICs, as previously described [15]. Inflammation score was categorized focusing only on the degree of TAICs, regardless of the location, as follows: score 0, no TAICs; score 1, focal or low TAICs; and score 2, diffuse or high TAICs. Figure 2 shows representative images of the categorization of TAICs status.

2.3. Immunohistochemical Analysis

We used the previously reported data assessed on tissue microarray sections from 2 mm-cores of formalin-fixed, paraffin-embedded (FFPE) blocks to evaluate the association of the histologic-based TAICs with CD8 and PD-L1 immunohistochemical expression for the principal cohort (n = 121 and 429, respectively) [21,22,23]. We performed immunohistochemistry (IHC) using Leica Bond-III (Leica Biosystems, Melbourne, Australia) and a Ventana Discovery Ultra autostainer (Roche Diagnostics K.K, Tokyo, Japan). Primary antibodies against CD8 (4B11, Prediluted; Leica Biosystems, Newcastle Upon Tyne, UK) and PD-L1 (28-8, 1:400; Abcam, Cambridge, MA, USA) were used to visualize with a BOND Polymer Refine Detection (Leica Biosystems) and OptiView DAB IHC Detection Kit (Ventana Medical System, Tucson, AZ, USA), respectively. We evaluated the density of CD8+ TILs (calculated as the number of cells/mm² per cores) and we semiquantitatively assessed the membranous staining pattern of PD-L1 in tumor cells using the H-score (0–300) as previously described [22,23].

2.4. Gene Expression Analysis

To validate the correlation between the histologic-based TAICs and gene expression signatures related to cancer immunity, the IMmotion 150 gene signatures [7], consisting of angiogenesis, immune and antigen presentation, and myeloid inflammation, were extracted from the RNA-sequencing data of TCGA according to previous reports [14,23]. Three gene signatures related to the underlying mechanisms of ICIs response were defined as follows: effector T-cell: CD8A, IFNG, GZMA, GZMB, PRF1, and EOMES; immune checkpoint: CD274 (PD-L1), CTLA4, and TIGIT; and myeloid: CXCL1, CXCL2, CXCL3, IL6, and PTGS2, as previously described [14]. To calculate gene signature scores, each gene score was normalized by the z-score across all patients and averaged to create signature scores for each patient [24]. The TCGA RNA-sequencing data were downloaded as described previously [25].

2.5. Statistical Analysis

All statistical analyses were performed using EZR version 1.54 (Saitama Medical Center, Jichi, Japan) [26]. All continuous data were shown as median-valued and interquartile ranges (IQR). A Chi-squared test for categorical variables and one-way ANOVA analysis, Mann–Whitney U test, or Kruskal–Wallis test for non-parametric variables were used to evaluate the statistical significance among three or four groups. The F-statistic in the linear regression analysis was calculated to determine the statistical significance among the three. The study outcome measure was recurrence-free survival (RFS), defined as the time from surgery to initial local or distant metastasis shown on imaging, and overall survival (OS) was defined as the time from surgery to any cause of death in the principal cohort and TCGA cohort, respectively. RFS or OS was assessed by the Kaplan–Meier method with the log-rank and the Cox proportional hazards models. Harrell’s concordance index (c-index) was used to compare the predictive accuracy of the Cox models. A two-sided p < 0.05 was considered statistically significant.

3. Results

3.1. Clinicopathological Characteristics

In the principal cohort, the median age at ccRCC diagnosis was 65 years (IQR, 56–73 years). In total, 103 (23.7%) of the tumors were in the high stage (TNM stage III or IV) and 142 (32.5%) were of a high WHO/ISUP grade (3 or 4). Of the 436 patients, 57 (13.1%) experienced a recurrence and 15 (3.4%) died of ccRCC during a median follow-up period of 61.7 months (IQR, 33.8–93.6), as previously reported [14].

In the TCGA cohort, of the 162 patients with ccRCC, 67 (41.4%) died during a median follow-up period of 1173 days (IQR, 563–1779). Of the tumors, 65 (40.1%) were in the high stage (TNM stage III or IV) and 69 (42.5%) were a high WHO/ISUP grade (3 or 4). The histologically identified TAICs status assessed by the three methods is summarized in Table 2. A statistically significant relationship was found between two cohorts in each methodology. The clinicopathological characteristics of three-tier and four-tier immunophenotype and inflammation scores in both the principal and TCGA cohorts are shown in Supplementary Tables S1–S6.

3.2. Comparison of Immunohistochemical Expression among Three-Tier and Four-Tier Immunophenotype and Inflammation Score

CD8+ TILs were significantly enriched in the inflamed compared to the desert (p < 0.001), followed by the excluded in the three-tier immunophenotype. CD8+ TILs were significantly distributed among all four types in the four-tier immunophenotype, whereas scores 1 and 2 made up a significantly higher density of CD8+ TILs compared to score 0 in the inflammation score (p < 0.01) (Figure 3A). Representative examples of CD8 IHC are shown in Figure 3B. Inflamed and excluded showed significantly higher PD-L1 expression compared to the desert in the three-tier immunophenotype (p < 0.01). In the four-tier immunophenotype, PD-L1 expression of hot and excluded was significantly higher compared to immunosuppressed and cold (p < 0.01). In the inflammation score, scores 2 and 1 had significantly higher expression compared to score 0 (p < 0.01) (Figure 3C). Representative examples of PD-L1 IHC are shown in Figure 3D.

3.3. Comparison of Gene Expression Signatures among Three-Tier and Four-Tier Immunophenotypes and Inflammation Scores

In the TCGA cohort, significant difference among the three scores were found in effector T-cell and myeloid gene signatures in the three-tier immunophenotype (both p < 0.05), whereas no significant difference was found in the immune checkpoint gene signature (Figure 4A). In the four-tier immunophenotype, effector T-cell and immune checkpoint gene signature was significantly enriched in hot and excluded compared to immunosuppressed and cold (p < 0.001). Furthermore, significant differences among the four types were found in the myeloid gene signature, with the highest expression of excluded (p < 0.05) (Figure 4B). In the inflammation score, the effector T-cell and immune checkpoint gene signature were significantly enriched in score 2 (p < 0.001), followed by score 1 (p < 0.001), and significant differences were noted among the three types in myeloid gene signature (p < 0.01) (Figure 4C). Correlation analysis among the three methodologies showed that the inflammation score was the highest in the effecter T-cell (f-statistic = 23.92, p < 0.001), immune checkpoint (f-statistic = 26.73, p < 0.001), and myeloid (f-statistic = 7.23, p = 0.008) gene signatures. Comparing the three-tier and four-tier immunophenotypes, the four-tier immunophenotype was higher in all gene signatures: effecter T-cell (f-statistic = 14.29, p < 0.001), immune checkpoint (f-statistic = 15.91, p < 0.001), and myeloid (f-statistic = 6.77, p < 0.001)), compared to the three-tier immunophenotypes (effecter T-cell (f-statistic = 3.48, p = 0.03), immune checkpoint (f-statistic = 4.08, p = 0.02), and myeloid (f-statistic = 1.59, p = 0.21) (Figure 5).

3.4. Comparison of Patient Outcome among Three-Tier and Four-Tier Immunophenotype and Inflammation Score

The Kaplan–Meier survival analysis showed a five-year RFS rate of 78.7% for inflamed and 74.5% for excluded versus 95% for desert in the three-tier immunophenotype; 79.8% for immunosuppressed, 78.6% for hot, and 67.3% for excluded versus 95% for cold in the four-tier immunophenotype; and 79.8% for score 1 and 71.7% for score 2 versus 95% for score 0 in the inflammation score (Figure 6A–C). The Kaplan–Meier survival analysis of the TCGA-cohort showed a similar trend compared to the principal cohort: a five-year OS rate of 45.9% for inflamed and 31.0% for excluded versus 74.9% for desert in the three-tier immunophenotype; 44.9% for immunosuppressed and 28.6% for hot versus 74.9% for cold in the four-tier immunophenotype (excluded could not be assessed due to the limited number of cases); and 44.9% for score 1 and 20.8% for score 2 versus 74.9% for score 0 in the inflammation score (Figure 6D–F). Moreover, the four-tier immunophenotype and inflammation score more accurately predicted RFS or OS than the three-tier immunophenotype in both cohorts (c-index = 0.692 vs. 0.672 for the principal cohort and 0.685 and 0.684 vs. 0.661 for the TCGA cohort) (Figure 7A,B). Univariate associations with recurrence after nephrectomy in patients with clear cell renal cell carcinoma are shown in Table 3.

3.5. Association of Inflammation Score with Pathological Prognostic Factors

Considering both correlation analysis and the prognostic accuracy among the three methodologies, the inflammation score was the most correlated with the gene expression signatures (Figure 5) and prognostic prediction (Figure 7). Regarding the association of the inflammation score with TNM stage pathological prognostic factors, WHO/ISUP grade, sarcomatoid/rhabdoid components, and necrosis factors, significant differences among scores were found in both the principal and TCGA cohort (all p < 0.001; Figure 8A,B).

4. Discussion

Using the localized ccRCC and TCGA-KIRC cohorts, we demonstrated an association between the three methodologies of histologic-based TAICs considering the location and/or degree of TAICs with immunohistochemical expression and gene expression signatures related to cancer immunity. Although all three methodologies correlated with immunohistochemical and immune gene expression signatures, we showed the inflammation score based on the degree of TAICs the most correlated with effector T-cell, immune checkpoint, and myeloid gene signatures among the three methodologies. Regarding the prognostic prediction, inflammation score and four-tier immunophenotype based on both the location and degree of TAICs showed similarly higher accuracy in predicting RFS or OS compared to the three-tier immunophenotype based on the TAICs location. Thus, a simple histologic assessment of TIACs may predict clinical outcomes and immunotherapy responses.

While the previously reported methodology focused on the TILs, we assessed the TAICs status including both mononuclear cells and granulocytes on only H&E-stained slides. Currently, the PD-L1 expression on immune or tumor cells, the extensive infiltration of CD8+ TILs, and the high tumor mutational burden are the most sensitive and specific biomarkers of clinical response to checkpoint blockade in several solid tumors [10]. Although many different approaches can assess and describe the immune response, a standardized methodology to assess TILs on H&E-stained sections for solid tumors is a valid and reproducible biomarker in routine clinical practice [9,27,28]. However, whether the assessment of only TILs reflects cancer immunity remains unclear.

Among various features of the tumor microenvironment of ccRCC, Cai and Christie et al. have recently shown that the intratumoral neutrophilic infiltration of ccRCC is one independent predictor for disease-free survival [13]. Although effective immunotherapy promotes the killing of cancer cells by cytotoxic T cells [10], CD4+ helper T cells, regulatory T cells (Tregs), effector T-cells, macrophages, and neutrophils are also important immune factors that can either promote or block tumor development [29]. Because host immunity, consisting of various types of TAICs and cellular matrix, plays an active role in controlling tumor growth and metastatic spreading [30,31], our simple and comprehensive assessment of TAICs is reasonable in this current study.

To the best of our knowledge, histologic-based TAICs status for ccRCC has not been extensively investigated. We first compared the association of the three TAICs status with immunohistochemical and gene expression and oncological outcomes, using three methodologies: the three-tier immunophenotypes inflamed, immune desert, or immune-excluded based on the spatial localization of immune cells [14,32]; the four-tier immunophenotypes referred to the novel concept of hot, excluded, immunosuppressed, and cold phenotypes, including the nature, density, immune functional orientation, and distribution of immune cell parameters [14]; and our proposed inflammation score focuses only on the degree of TAICs. Although these phenotypes were originally defined based on CD8+ TIL, we categorized these phenotypes by histologically identified TAICs. Nevertheless, TAICs detected on H&E-stained slides are highly correlated with immunohistochemical expression of CD8 and effecter T-cell gene signatures, including CD8A, which suggest histologically identified TAICs within ccRCC predominantly include CD8+ TILs.

We found that the four-tier immunophenotype considering both the location and degree of TAICs was more highly correlated with gene expression signatures related to cancer immunity compared to the three-tier immunophenotype based on the location of TAICs. In the present study, however, we demonstrated that the inflammation score considering only the degree of TAICs was the most correlated with both gene expression signatures and prognostic prediction among the three methodologies. Accordingly, our proposed simple inflammation score of TAICs by H&E-stained slides might be easier to apply than other methods in a routine clinical setting.

Concerning the prognostic prediction of infiltrating immune cells, evidence suggests that increased infiltrations with CD8+ and Th1-CD4+ T-cells are associated with a good prognosis in most solid tumors [33]. Our prior study of bladder cancer revealed that high TAICs were significantly associated with a favorable oncological outcome compared to low TAICs [34], which is supported by other studies [35]. Contrary to other solid tumors, RCC enriched in high-density CD8+ TILs are associated with poor prognosis [36,37,38]. Consist with previous studies of RCC, we found that degree of inflammation could stratify patient outcome, which implicated increased TAICs involvement in the formation of pre-metastatic niches.

Several studies have revealed that patients with the immune-excluded phenotype had a poor prognosis compared to patients with the inflamed or desert phenotype in various tumors because immune-excluded tumors were immunosuppressed by T-cells embedded in the tumor stromal microenvironment with upregulation of TGFβ, myeloid inflammation, and angiogenesis [10,14,39]. Consistent with previous studies, the prognosis of the excluded type was the worst in both the three-tier and four-tier immunophenotypes. Interestingly, our results showed that the myeloid gene signature of the excluded type was significantly high in the four-tier immunophenotype, although myeloid gene signatures were not statistically significant in the three-tier immunophenotype. Thus, when considering the myeloid-derived suppressor cells (MDSCs) involved in RCC progression and drug resistance, the four-tier immunophenotype would help determine patient selection for immunotherapy. Further studies are needed to determine which TAICs assessment methodology plays a predictive role in the selection of ICIs for ccRCC.

Regarding the association of pathological prognostic factors, we found that a dedifferentiation form of ccRCC, such as sarcomatoid/rhabdoid, was significantly correlated with inflammation score 3, followed by score 2, which corresponded with the enrichment of immune checkpoint signatures. Beuselinck et al. also demonstrated a correlation between poor tumor cell differentiation and overexpression of immunomodulatory molecules in ccRCCs using gene expression analysis [40]. These results corroborate that sarcomatoid/rhabdoid features could be potential surrogate markers of ICI response [41].

Our study has several limitations. First, this was a retrospective study including the excluded four-tier immunophenotype in the TCGA cohort with a small sample size. Because this study cohort was the same as that which was previously described [14], multiple institutional cohorts should be validated. Second, CD8 and PD-L1 immunohistochemical expression was evaluated with TMAs constructed with two representative cores and not using whole sections. Third, we correlated TAICs assessment and gene expression signatures using only the TCGA cohort. Fourth, the principal cohort candidates were localized ccRCC patients with only an endpoint of RFS. The candidate for the TCGA-KIRC cohort included metastatic ccRCC patients with only an available endpoint of OS, as previously described [14]. Although there was a statistically significant relationship between two cohorts in each methodology (Table 2), we have confirmed a similar trend in the association of TAICs status with clinicopathological factors (Supplementary Tables S1–S6). Despite these limitations, we comprehensively showed the association of histologic-based TAICs with gene signatures related to cancer immunity and clinical outcomes.

5. Conclusions

Our current work showed histologic-based TAICs status identified by H&E-stained slides correlated with effector T-cell, immune checkpoint, myeloid gene signatures, and clinical outcomes. Histologic assessment of TIACs may serve as a surrogate for gene expression signatures related to immunotherapy response and support prognostic prediction. Our approach using only H&E-stained slides can be performed in routine clinical practice at a low cost.

Author Contributions

Conceptualization, C.O. and T.Y.; Data curation, T.Y., J.I. and H.O.; Formal analysis, T.Y. and J.I.; Investigation, C.O. and R.Y.; Methodology, C.O., T.Y., T.T., R.O. and N.A.; Resources, R.O., Y.Y. and H.K.; Supervision, T.T., R.S., K.H., H.K. and K.T.; Visualization, C.O.; Writing—original draft, C.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Japan Society for the Promotion of Science KAKENHI fund (Grants No. 19K16875 to C.O., and Grants No. 20K16457 to H.O.).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Kansai Medical University Hospital (No. 2018109 (28/11/2018) and No. 2020222 (09/12/2020)).

Informed Consent Statement

Informed consent was obtained in the form of opt-out on the website of Kansai Medical University Hospital. No one expressed a refusal.

Data Availability Statement

The data are available upon reasonable request by contacting the corresponding author.

Acknowledgments

This work was supported by the research grant from Kansai Medical University (KMU) research consortium.

Conflicts of Interest

C.O. and T.T. received research funding from Chugai Pharmaceutical Co. Ltd. (Tokyo, Japan) outside the submitted work. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. The remaining authors have no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Table 1

Comparison of three evaluation methods of the tumor immune microenvironment.

Three-tier immunophenotype	Based on the location of TAICs, regardless of the TAICs degree
	Desert	Excluded	Inflamed
	No TAICs	Peritumoral TAICs	Intratumoral TAICs
Four-tier immunophenotype	Considering both the location and degree of TAICs
	Cold	Immunosuppressed	Excluded	Hot
	No TAICs	Focal or low TAICs, regardless of the TAICs location	Diffuse or highperitumoral TAICs	Diffuse or highintratumoral TAICs
Inflammation score	Based on the degree of TAICs, regardless of the TAICs location
	Score 0	Score 1	Score 2
	No TAICs	Focal or low TAICs	Diffuse or high TAICs

TAICs, tumor-associated immune cells.

Table 2

Tumor-associated immune cells status assessed by three methods.

Variables	Principal Cohort n = 436	TCGA Cohortn = 162	p Value
Three-tier immunophenotype, n (%)			0.029
Desert	243 (55.7)	75 (46.3)
Excluded	91 (20.9)	32 (19.8)
Inflamed	102 (23.4)	55 (34.0)
Four-tier immunophenotype, n (%)			0.003
Cold	243 (55.7)	75 (46.3)
Immunosuppressed	127 (29.1)	67 (41.4)
Excluded	35 (8.0)	4 (2.5)
Hot	31 (7.1)	16 (9.9)
Inflammation score, n (%)			0.018
Score 0	243 (55.7)	75 (46.3)
Score 1	127 (29.1)	67 (41.4)
Score 2	66 (15.1)	20 (12.3)

Table 3

Univariate associations with recurrence after nephrectomy in patients with clear cell renal cell carcinoma.

	Cohort	Principal Cohort		TCGA Cohort
	Variables	HR (95% CI)	p Value	HR (95% CI)	p Value
Three-tierimmunophenotype	Desert	1 (ref.)		1 (ref.)
	Excluded	5.23 (2.63–10.39)	<0.001	3.93 (2.05–7.54)	<0.001
	Inflamed	4.44 (2.23–8.81)	<0.001	3.10 (1.66–5.79)	<0.001
Four-tierimmunophenotype	Cold	1 (ref.)		1 (ref.)
	Immunosuppressed	3.87 (1.97–7.59)	<0.001	2.94 (1.62–5.36)	<0.001
	Excluded	7.60 (3.47–16.7)	<0.001	8.49 (2.79–25.9)	<0.001
	Hot	5.82 (2.41–14.1)	<0.001	5.20 (2.38–11.4)	<0.001
Inflammation score	Score 0	1 (ref.)		1 (ref.)
	Score 1	3.87 (1.97–7.59)	<0.001	2.94 (1.62–5.36)	<0.001
	Score 2	6.77 (3.37–13.6)	<0.001	5.79 (2.82–11.9)	<0.001

HR, hazard ratio; CI, confidence interval; ref., reference.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines10020323/s1, Table S1: Clinicopathological characteristics of three-tier immunophenotype in 436 cases with localized ccRCCs in the principal cohort; Table S2: Clinicopathological characteristics of three-tier immunophenotype in 162 cases with ccRCCs in the TCGA cohort; Table S3: Clinicopathological characteristics of four-tier immunophenotype in 436 cases with localized ccRCCs in the principal cohort; Table S4: Clinicopathological characteristics of four-tier immunophenotype in 162 cases with ccRCCs in the TCGA cohort; Table S5: Clinicopathological characteristics of inflammation score in 436 cases with localized ccRCCs in the principal cohort; Table S6: Clinicopathological characteristics of inflammation score in 162 cases with ccRCCs in the TCGA cohort. ccRCCs; clear cell renal cell carcinomas.

References

1. Vaddepally, R.K.; Kharel, P.; Pandey, R.; Garje, R.; Chandra, A.B. Review of Indications of FDA-Approved Immune Checkpoint Inhibitors per NCCN Guidelines with the Level of Evidence. Cancers; 2020; 12, 738. [DOI: https://dx.doi.org/10.3390/cancers12030738] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32245016]

2. Shuch, B.; Amin, A.; Armstrong, A.J.; Eble, J.N.; Ficarra, V.; Lopez-Beltran, A.; Martignoni, G.; Rini, B.I.; Kutikov, A. Understanding pathologic variants of renal cell carcinoma: Distilling therapeutic opportunities from biologic complexity. Eur. Urol.; 2015; 67, pp. 85-97. [DOI: https://dx.doi.org/10.1016/j.eururo.2014.04.029] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24857407]

3. Casuscelli, J.; Vano, Y.A.; Fridman, W.H.; Hsieh, J.J. Molecular classification of renal cell carcinoma and its implication in future clinical practice. Kidney Cancer; 2017; 1, pp. 3-13. [DOI: https://dx.doi.org/10.3233/KCA-170008]

4. Şenbabaoğlu, Y.; Gejman, R.S.; Winer, A.G.; Liu, M.; Van Allen, E.M.; de Velasco, G.; Miao, D.; Ostrovnaya, I.; Drill, E.; Luna, A. et al. Tumor immune microenvironment characterization in clear cell renal cell carcinoma identifies prognostic and immunotherapeutically relevant messenger RNA signatures. Genome. Biol.; 2016; 17, 231. [DOI: https://dx.doi.org/10.1186/s13059-016-1092-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27855702]

5. Ghatalia, P.; Gordetsky, J.; Kuo, F.; Dulaimi, E.; Cai, K.Q.; Devarajan, K.; Bae, S.; Naik, G.; Chan, T.A.; Uzzo, R. et al. Prognostic impact of immune gene expression signature and tumor infiltrating immune cells in localized clear cell renal cell carcinoma. J. Immunother. Cancer; 2019; 7, 139. [DOI: https://dx.doi.org/10.1186/s40425-019-0621-1]

6. Wallin, J.J.; Bendell, J.C.; Funke, R.; Sznol, M.; Korski, K.; Jones, S.; Hernandez, G.; Mier, J.; He, X.; Hodi, F.S. et al. Atezolizumab in combination with bevacizumab enhances antigen-specific T-cell migration in metastatic renal cell carcinoma. Nat. Commun.; 2016; 7, 12624. [DOI: https://dx.doi.org/10.1038/ncomms12624]

7. McDermott, D.F.; Huseni, M.A.; Atkins, M.B.; Motzer, R.J.; Rini, B.I.; Escudier, B.; Fong, L.; Joseph, R.W.; Pal, S.K.; Reeves, J.A. et al. Clinical activity and molecular correlates of response to atezolizumab alone or in combination with bevacizumab versus sunitinib in renal cell carcinoma. Nat. Med.; 2018; 24, pp. 749-757. [DOI: https://dx.doi.org/10.1038/s41591-018-0053-3]

8. Au, L.; Hatipoglu, E.; Robert de Massy, M.; Litchfield, K.; Beattie, G.; Rowan, A.; Schnidrig, D.; Thompson, R.; Byrne, F.; Horswell, S. et al. Determinants of anti-PD-1 response and resistance in clear cell renal cell carcinoma. Cancer Cell; 2021; 39, pp. 1497-1518.e11. [DOI: https://dx.doi.org/10.1016/j.ccell.2021.10.001]

9. Hendry, S.; Salgado, R.; Gevaert, T.; Russell, P.A.; John, T.; Thapa, B.; Christie, M.; van de Vijver, K.; Estrada, M.V.; Gonzalez-Ericsson, P.I. et al. Assessing Tumor-Infiltrating Lymphocytes in Solid Tumors: A Practical Review for Pathologists and Proposal for a Standardized Method from the International Immuno-Oncology Biomarkers Working Group: Part 2: TILs in Melanoma, Gastrointestinal Tract Carcinomas, Non-Small Cell Lung Carcinoma and Mesothelioma, Endometrial and Ovarian Carcinomas, Squamous Cell Carcinoma of the Head and Neck, Genitourinary Carcinomas, and Primary Brain Tumors. Adv. Anat. Pathol.; 2017; 24, pp. 311-335.

10. Hegde, P.S.; Karanikas, V.; Evers, S. The where, then when, and the how of immune monitoring for cancer immunotherapies in the era of checkpoint inhibition. Clin. Cancer Res.; 2016; 22, pp. 1865-1874. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-15-1507]

11. Hegde, P.S.; Chen, D.S. Top 10 Challenges in Cancer Immunotherapy. Immunity; 2020; 52, pp. 17-35. [DOI: https://dx.doi.org/10.1016/j.immuni.2019.12.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31940268]

12. Verine, J.; Colin, D.; Nheb, M.; Prapotnich, D.; Ploussard, G.; Cathelineau, X.; Desgrandchamps, F.; Mongiat-Artus, P.; Feugeas, J.P. Architectural Patterns are a Relevant Morphologic Grading System for Clear Cell Renal Cell Carcinoma Prognosis Assessment: Comparisons With WHO/ISUP Grade and Integrated Staging Systems. Am. J. Surg. Pathol.; 2018; 42, pp. 423-441. [DOI: https://dx.doi.org/10.1097/PAS.0000000000001025] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29356723]

13. Cai, Q.; Christie, A.; Rajaram, S.; Zhou, Q.; Araj, E.; Chintalapati, S.; Cadeddu, J.; Margulis, V.; Pedrosa, I.; Rakheja, D. et al. Ontological analyses reveal clinically-significant clear cell renal cell carcinoma subtypes with convergent evolutionary trajectories into an aggressive type. EBioMedicine; 2020; 51, 102526. [DOI: https://dx.doi.org/10.1016/j.ebiom.2019.10.052] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31859241]

14. Ohe, C.; Yoshida, T.; Amin, M.B.; Atsumi, N.; Ikeda, J.; Saiga, K.; Noda, Y.; Yasukochi, Y.; Ohashi, R.; Ohsugi, H. et al. Development and validation of a vascularity-based architectural classification for clear cell renal cell carcinoma: Correlation with conventional pathological prognostic factors, gene expression patterns, and clinical outcomes. Mod. Pathol.; 2021; [DOI: https://dx.doi.org/10.1038/s41379-021-00982-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34848832]

15. Galon, J.; Bruni, D. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies. Nat. Rev. Drug Discov.; 2019; 18, pp. 197-218. [DOI: https://dx.doi.org/10.1038/s41573-018-0007-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30610226]

16. Ricketts, C.J.; De Cubas, A.A.; Fan, H.; Smith, C.C.; Lang, M.; Reznik, E.; Bowlby, R.; Gibb, E.A.; Akbani, R.; Beroukhim, R. et al. The Cancer Genome Atlas comprehensive molecular characterization of renal cell carcinoma. Cell Rep.; 2018; 23, pp. 313-326.e5. [DOI: https://dx.doi.org/10.1016/j.celrep.2018.03.075] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29617669]

17. Broad GDAC FIREHOSE-Broad Institute. Available online: http://gdac.broadinstitute.org/ (accessed on 1 August 2021).

18. Digital Slide Archive (DSA). Available online: https://cancer.digitalslidearchive.org/ (accessed on 1 August 2021).

19. Brierley, J.D.; Gospodarowics, M.K.; Wittekind, C. Union for International Cancer Control. TNM Classification of Malignant Tumours; 8th ed. Wiley: New York, NY, USA, 2017.

20. Moch, H.; Humphrey, P.A.; Ulbright, T.M.; Reuter, V.E. WHO Classification of Tumours of the Urinary System and Male Genital Organs; 4th ed. IARC: Lyon, France, 2016.

21. Ohsugi, H.; Yoshida, T.; Ohe, C.; Ikeda, J.; Sugi, M.; Kinoshita, H.; Tsuta, K.; Matsuda, T. The SSPN Score, a novel scoring system incorporating PBRM1 expression, predicts postoperative recurrence for patients with non-metastatic clear cell renal cell carcinoma. Ann. Surg. Oncol.; 2021; 28, pp. 2359-2366. [DOI: https://dx.doi.org/10.1245/s10434-020-09075-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32940805]

22. Ikeda, J.; Ohe, C.; Yoshida, T.; Ohsugi, H.; Sugi, M.; Tsuta, K.; Kinoshita, H. PD-L1 Expression and Clinicopathological Factors in Renal Cell Carcinoma: A Comparison of Antibody Clone 73-10 With Clone 28-8. Anticancer Res.; 2021; 41, pp. 4577-4586. [DOI: https://dx.doi.org/10.21873/anticanres.15271]

23. Yoshida, T.; Ohe, C.; Ikeda, J.; Atsumi, N.; Ohsugi, H.; Sugi, M.; Higasa, K.; Saito, R.; Tsuta, K.; Matsuda, T. et al. Eosinophilic features in clear cell renal cell carcinoma correlate with outcomes of immune checkpoint and angiogenesis blockade. J. Immunother. Cancer; 2021; 9, e002922. [DOI: https://dx.doi.org/10.1136/jitc-2021-002922]

24. Bolen, C.R.; McCord, R.; Huet, S.; Frampton, G.M.; Bourgon, R.; Jardin, F.; Dartigues, P.; Punnoose, E.A.; Szafer-Glusman, E.; Xerri, L. et al. Mutation load and an effector T-cell gene signature may distinguish immunologically distinct and clinically relevant lymphoma subsets. Blood Adv.; 2017; 22, pp. 1884-1890. [DOI: https://dx.doi.org/10.1182/bloodadvances.2016000786]

25. Ohashi, R.; Angori, S.; Batavia, A.A.; Rupp, N.J.; Ajioka, Y.; Schraml, P.; Moch, H. Loss of CDKN1A mRNA and protein expression are independent predictors of poor outcome in chromophobe renal cell carcinoma patients. Cancers; 2020; 12, 465. [DOI: https://dx.doi.org/10.3390/cancers12020465] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32079343]

26. Kanda, Y. Investigation of the freely available easy-to-use software‘EZR’ for medical statistics. Bone Marrow Transpl.; 2013; 48, pp. 452-458. [DOI: https://dx.doi.org/10.1038/bmt.2012.244] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23208313]

27. Salgado, R.; Denkert, C.; Demaria, S.; Sirtaine, N.; Klauschen, F.; Pruneri, G.; Wienert, S.; Van den Eynden, G.; Baehner, F.L.; Penault-Llorca, F. et al. The evaluation of tumor-infiltrating lymphocytes (TILs) in breast cancer: Recommendations by an International TILs Working Group 2014. Ann. Oncol.; 2015; 26, pp. 259-271. [DOI: https://dx.doi.org/10.1093/annonc/mdu450] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25214542]

28. Denkert, C.; Wienert, S.; Poterie, A.; Loibl, S.; Budczies, J.; Badve, S.; Bago-Horvath, Z.; Bane, A.; Bedri, S.; Brock, J. et al. Standardized evaluation of tumor-infiltrating lymphocytes in breast cancer: Results of the ring studies of the international immuno-oncology biomarker working group. Mod. Pathol.; 2016; 29, pp. 1155-1164. [DOI: https://dx.doi.org/10.1038/modpathol.2016.109] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27363491]

29. Zhang, S.; Zhang, E.; Long, J.; Hu, Z.; Peng, J.; Liu, L.; Tang, F.; Li, L.; Ouyang, Y.; Zeng, Z. Immune infiltration in renal cell carcinoma. Cancer Sci.; 2019; 110, pp. 1564-1572. [DOI: https://dx.doi.org/10.1111/cas.13996]

30. Heidegger, I.; Pircher, A.; Pichler, R. Targeting the Tumor Microenvironment in Renal Cell Cancer Biology and Therapy. Front. Oncol.; 2019; 9, 490. [DOI: https://dx.doi.org/10.3389/fonc.2019.00490]

31. Grivennikov, S.I.; Greten, F.R.; Karin, M. Immunity, inflammation, and cancer. Cell; 2010; 140, pp. 883-899. [DOI: https://dx.doi.org/10.1016/j.cell.2010.01.025]

32. Brück, O.; Lee, M.H.; Turkki, R.; Uski, I.; Penttilä, P.; Paavolainen, L.; Kovanen, P.; Järvinen, P.; Bono, P.; Pellinen, T. et al. Spatial immunoprofiling of the intratumoral and peritumoral tissue of renal cell carcinoma patients. Mod. Pathol.; 2021; 34, pp. 2229-2241. [DOI: https://dx.doi.org/10.1038/s41379-021-00864-0]

33. Fridman, W.H.; Pagès, F.; Sautès-Fridman, C.; Galon, J. The immune contexture in human tumours: Impact on clinical outcome. Nat. Rev. Cancer; 2012; 12, pp. 298-306. [DOI: https://dx.doi.org/10.1038/nrc3245]

34. Ikeda, J.; Ohe, C.; Yoshida, T.; Kuroda, N.; Saito, R.; Kinoshita, H.; Tsuta, K.; Matsuda, T. Comprehensive pathological assessment of histological subtypes, molecular subtypes based on immunohistochemistry, and tumor-associated immune cell status in muscle-invasive bladder cancer. Pathol. Int.; 2021; 71, pp. 173-182. [DOI: https://dx.doi.org/10.1111/pin.13060]

35. Hodgson, A.; Xu, B.; Satkunasivam, R.; Downes, M.R. Tumour front inflammation and necrosis are independent prognostic predictors in high-grade urothelial carcinoma of the bladder. J. Clin. Pathol.; 2018; 71, pp. 154-160. [DOI: https://dx.doi.org/10.1136/jclinpath-2017-204562] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28768702]

36. Nakano, O.; Sato, M.; Naito, Y.; Suzuki, K.; Orikasa, S.; Aizawa, M.; Suzuki, Y.; Shintaku, I.; Nagura, H.; Ohtani, H. Proliferative activity of intratumoral CD8(+) T-lymphocytes as a prognostic factor in human renal cell carcinoma: Clinicopathologic demonstration of antitumor immunity. Cancer Res.; 2001; 61, pp. 5132-5136. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11431351]

37. Giraldo, N.A.; Becht, E.; Pagès, F.; Skliris, G.; Verkarre, V.; Vano, Y.; Mejean, A.; Saint-Aubert, N.; Lacroix, L.; Natario, I. et al. Orchestration and Prognostic Significance of Immune Checkpoints in the Microenvironment of Primary and Metastatic Renal Cell Cancer. Clin. Cancer Res.; 2015; 21, pp. 3031-3040. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-14-2926] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25688160]

38. Giraldo, N.A.; Becht, E.; Vano, Y.; Petitprez, F.; Lacroix, L.; Validire, P.; Sanchez-Salas, R.; Ingels, A.; Oudard, S.; Moatti, A. et al. Tumor-Infiltrating and Peripheral Blood T-cell Immunophenotypes Predict Early Relapse in Localized Clear Cell Renal Cell Carcinoma. Clin. Cancer Res.; 2017; 23, pp. 4416-4428. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-16-2848]

39. Desbois, M.; Udyavar, A.R.; Ryner, L.; Kozlowski, C.; Guan, Y.; Dürrbaum, M.; Lu, S.; Fortin, J.P.; Koeppen, H.; Ziai, J. et al. Integrated digital pathology and transcriptome analysis identifies molecular mediators of T-cell exclusion in ovarian cancer. Nat. Commun.; 2020; 11, 5583. [DOI: https://dx.doi.org/10.1038/s41467-020-19408-2]

40. Beuselinck, B.; Job, S.; Becht, E.; Karadimou, A.; Verkarre, V.; Couchy, G.; Giraldo, N.; Rioux-Leclercq, N.; Molinié, V.; Sibony, M. et al. Molecular subtypes of clear cell renal cell carcinoma are associated with sunitinib response in the metastatic setting. Clin. Cancer Res.; 2015; 21, pp. 1329-1339. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-14-1128]

41. Bakouny, Z.; Braun, D.A.; Shukla, S.A.; Pan, W.; Gao, X.; Hou, Y.; Flaifel, A.; Tang, S.; Bosma-Moody, A.; He, M.X. et al. Integrative molecular characterization of sarcomatoid and rhabdoid renal cell carcinoma. Nat. Commun.; 2021; 12, 808. [DOI: https://dx.doi.org/10.1038/s41467-021-21068-9]