Tumor-infiltrating neutrophils predict therapeutic benefit of tyrosine kinase inhibitors in metastatic renal cell carcinoma

SPP1 expression indicates outcome of immunotherapy plus tyrosine kinase inhibition in advanced renal cell carcinoma

Clinical guidelines have recently advised combination therapy involving immunotherapy (IO) and tyrosine kinase inhibitors (TKI) as the first-line therapy approach for advanced renal cell carcinoma (RCC). Nevertheless, there is currently no available biomarker that can effectively distinguish the progression-free survival (PFS). RNA-sequencing and immunohistochemistry were conducted on our cohort of metastatic RCC patients, namely ZS-MRCC, who received combination therapy consisting of IO and TKI. We further applied RNA-sequencing, immunohistochemistry, and flow cytometry to examine the immune cell infiltration and functionality inside the tumor microenvironment of high-risk localized RCC samples. SPP1 expression was significantly higher in non-responders to IO-TKI therapy. Elevated levels of SPP1 were associated with poor PFS in both the ZS-MRCC cohort (HR = 2.73, p = .018) and validated in the JAVELIN Renal 101 cohort (HR = 1.61, p = .004). By multivariate Cox analysis, SPP1 was identified as a significant independent prognosticator. Furthermore, there existed a negative correlation between elevated levels of SPP1 and the presence of GZMB+CD8+ T cells (Spearman's ρ= -0.48, p < .001). Conversely, SPP1 expression is associated with T cell exhaustion markers. A significant increase in the abundance of Tregs was observed in tumors with high levels of SPP1. Additionally, a machine-learning-based model was constructed to predict the benefit of IO-TKI treatment. High SPP1 is associated with therapeutic resistance and unfavorable PFS in IO-TKI therapy. SPP1 expression have also been observed to be indicative of malfunction and exhaustion in T cells. Increased SPP1 expression has the potential to serve as a potential biomarker for treatment selection of metastatic RCC.

Effect of Annexin A2 on prognosis and sensitivity to immune checkpoint plus tyrosine kinase inhibition in metastatic renal cell carcinoma

BACKGROUND

Immunotherapy (IO) plus tyrosine kinase inhibitor (TKI) therapy is the first-line recommendation for advanced renal cell carcinoma (RCC), but no biomarker has been approved for it. Annexin A2 (ANXA2) can induce immune escape in tumors.

METHODS

Two independent cohorts of advanced RCC treated by IO + TKI were utilized for survival analysis (ZS-MRCC, n = 45; Javelin-101, n = 726). ANXA2 expression was determined by RNA-sequencing. The impact of ANXA2 on the tumor microenvironment was assessed by RNA-sequencing, flow cytometry and immunohistochemistry in two localized RCC datasets (ZS-HRRCC, n = 40; TCGA-KIRC, n = 530).

RESULTS

ANXA2 was upregulated in non-responders of IO + TKI therapy (p = 0.027). High-ANXA2 group showed poor progression-free survival (PFS) in both the ZS-MRCC cohort (HR, 2.348; 95% CI 1.084-5.085; P = 0.025) and the Javelin-101 cohort (HR, 1.472; 95% CI 1.043-2.077; P = 0.027). Multivariate Cox regression determined ANXA2 as an independent prognostic factor (HR, 2.619; 95% CI 1.194-5.746; P = 0.016). High-ANXA2 was correlated with decreased proportion of granzyme B+ CD8+ T cells (Spearman's ρ = - 0.40, P = 0.01), and increased TIM-3+ (Spearman's ρ = 0.43, P < 0.001) and CTLA4+ (Spearman's ρ = 0.49, P < 0.001) tumor-infiltrating lymphocytes. A random forest (RF) score was further build by integrating ANXA2 and immune genes, which stratified patients who would benefit from IO + TKI therapy (low-RF score, IO + TKI vs TKI, HR = 0.453, 95% CI 0.328-0.626; high-RF score, IO + TKI vs TKI, HR = 0.877, 95% CI 0.661-1.165; interaction P = 0.003).

CONCLUSIONS

Upregulated ANXA2 was associated with poor PFS and therapeutic resistance in RCC treated by IO + TKI therapy, and related with T cell exhaustion. The integrated RF score could stratify patients who would benefit from IO + TKI therapy.

The relationship between tumor-infiltrating neutrophils and clinical outcomes in patients with resectable hepatocellular carcinoma

BACKGROUND

The impact of tumor-infiltrating neutrophils (TINs) on clinical outcomes has been reported in various cancer types, but their role in hepatocellular carcinoma (HCC) has not been fully evaluated. The aim of this study was to investigate the prognostic values for TINs in HCC patients undergoing curative resection.

METHODS

We assessed immune markers (CD3, CD4, CD8, CD66b) using immunohistochemistry in 115 patients who underwent curative resection for HCC. We analyzed the prognostic values for tumor-infiltrating immune cells, including neutrophils, and other clinicopathological factors.

RESULTS

In the Multivariate Cox analysis of overall survival (OS), alpha-fetoprotein (AFP) ≥ 100 ng/mL (hazard ratio (HR), 2.74, 95% confidence interval (CI), 1.17-6.44; P = 0.021) and Barcelona Clinic Liver Cancer (BCLC) B/C stage (HR, 3.98, 95% CI, 1.68-9.43; P = 0.020) were found to be independent poor prognostic factors in HCC patients undergoing resection. The presence of CD66b+TINs was observed in 66 (57.4%) patients. However, CD66b+TINs were not associated with recurrence-free survival and OS.

CONCLUSIONS

Our study identified low CD66b+TINs in resectable HCC, and CD66b+ TINs did not have a significant role for the clinical outcomes of patients undergoing curative resection. The results suggest that TINs may play a role in more advanced stages of HCC.

Functional enrichment analysis of LYSET and identification of related hub gene signatures as novel biomarkers to predict prognosis and immune infiltration status of clear cell renal cell carcinoma

PURPOSE

The latest research shows that the lysosomal enzyme trafficking factor (LYSET) encoded by TMEM251 is a key regulator of the amino acid metabolism reprogramming (AAMR) and related pathways significantly correlate with the progression of some tumors. The purpose of this study was to explore the potential pathways of the TMEM251 in clear cell renal cell carcinoma (ccRCC) and establish related predictive models based on the hub genes in these pathways for prognosis and tumor immune microenvironment (TIME).

METHODS

We obtained mRNA expression data and clinical information of ccRCC samples from The Cancer Genome Atlas (TCGA), E-MATE-1980, and immunotherapy cohorts. Single-cell sequencing data (GSE152938) were downloaded from the Gene Expression Omnibus (GEO) database. We explored biological pathways of the LYSET by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of TMEM251-coexpression genes. The correlation of LYSET-related pathways with the prognosis was conducted by Gene Set Variation Analysis (GSVA) and unsupervised cluster analysis. The least absolute shrinkage and selection operator (LASSO) and Cox regression were used to identify hub prognostic genes and construct the risk score. Immune infiltration analysis was conducted by CIBERSORTx and Tumor Immune Estimation Resource (TIMER) databases. The predictive value of the risk score and hub prognostic genes on immunotherapy responsiveness was analyzed through the tumor mutation burden (TMB) score, immune checkpoint expression, and survival analysis. Immunohistochemistry (IHC) was finally used to verify the expressions of hub prognostic genes.

RESULTS

The TMEM251 was found to be significantly correlated with some AAMR pathways. AAGAB, ENTR1, SCYL2, and WDR72 in LYSET-related pathways were finally identified to construct a risk score model. Immune infiltration analysis showed that LYSET-related gene signatures significantly influenced the infiltration of some vital immune cells such as CD4 + cells, NK cells, M2 macrophages, and so on. In addition, the constructed risk score was found to be positively correlated with TMB and some common immune checkpoint expressions. Different predictive values of these signatures for Nivolumab therapy responsiveness were also uncovered in immunotherapy cohorts. Finally, based on single-cell sequencing analysis, the TMEM251 and the hub gene signatures were found to be expressed in tumor cells and some immune cells. Interestingly, IHC verification showed a potential dual role of four hub genes in ccRCC progression.

CONCLUSION

The novel predictive biomarkers we built may benefit clinical decision-making for ccRCC. Our study may provide some evidence that LYSET-related gene signatures could be novel potential targets for treating ccRCC and improving immunotherapy efficacy. Our nomogram might be beneficial to clinical choices, but the results need more experimental verifications in the future.

Favorable Immunotherapy Plus Tyrosine Kinase Inhibition Outcome of Renal Cell Carcinoma Patients with Low CDK5 Expression

PURPOSE

Immunotherapy (IO) plus tyrosine kinase inhibitor (TKI) has become the first-line treatment for advanced renal cell carcinoma, despite the lack of prognostic biomarkers. Cyclin-dependent kinase 5 (CDK5) affects the tumor microenvironment, which may influence the efficacy of TKI+IO.

MATERIALS AND METHODS

Two cohorts from our center (Zhongshan Metastatic Renal Cell Carcinoma [ZS-MRCC] cohort, Zhongshan High-risk Localized Renal Cell Carcinoma [ZS-HRRCC] cohort) and one cohort from a clinical trial (JAVELIN-101) were enrolled. The expression of CDK5 of each sample was determined by RNA sequencing. Immune infiltration and T cell function were evaluated by flow cytometry and immunohistochemistry. Response and progression-free survival (PFS) were set as primary endpoints.

RESULTS

Patients of low CDK5 expression showed higher objective response rate (60.0% vs. 23.3%) and longer PFS in both cohorts (ZS-MRCC cohort, p=0.014; JAVELIN-101 cohort, p=0.040). CDK5 expression was enhanced in non-responders (p < 0.05). In the ZS-HRRCC cohort, CDK5 was associated with decreased tumor-infiltrating CD8+ T cells, which was proved by immunohistochemistry (p < 0.05) and flow cytometry (Spearman's ρ=-0.49, p < 0.001). In the high CDK5 subgroup, CD8+ T cells revealed a dysfunction phenotype with decreased granzyme B, and more regulatory T cells were identified. A predictive score was further constructed by random forest, involving CDK5 and T cell exhaustion features. The RFscore was also validated in both cohorts. By utilizing the model, more patients might be distinguished from the overall cohort. Additionally, only in the low RFscore did TKI+IO outperform TKI monotherapy.

CONCLUSION

High-CDK5 expression was associated with immunosuppression and TKI+IO resistance. RFscore based on CDK5 may be utilized as a biomarker to determine the optimal treatment strategy.

Effects of PYCR1 on prognosis and immunotherapy plus tyrosine kinase inhibition responsiveness in metastatic renal cell carcinoma patients

BACKGROUND

Immunotherapy plus tyrosine kinase inhibitor (IO-TKI) has become the first-line management for metastatic renal cell carcinoma (RCC), despite the absence of biomarkers. Recently, pyrroline-5-carboxylate reductase 1 (PYCR1) and proline metabolism have been reported regulatory roles in the anti-tumor response.

METHODS

There were three cohorts enrolled: two from our institution (ZS-MRCC and ZS-HRRCC) and one from a clinical trial (JAVELIN-101). The PYCR1expression in each sample was evaluated by RNA sequencing. Flow cytometry and immunohistochemistry were performed to assess immune infiltration. Single-cell RNA-seq (scRNA-seq) data was used for cluster analysis of T cells and macrophages. Primary endpoints were set as response and progression-free survival (PFS).

RESULTS

Patients in the low-PYCR1 group had greater objective response rate (52.2% vs 18.2%) and longer PFS in both cohorts (ZS-MRCC cohort, P=0.01, HR=2.80; JAVELIN-101 cohort, P<0.001, HR=1.85). In responders, PYCR1 expression was decreased (P<0.05). In the high PYCR1 group, CD8+ T cells exhibited an exhausted phenotype with decreased GZMB (Spearman's ρ=-0.36, P=0.02). scRNA-seq revealed tissue-resident memory T (Trm) (P<0.05) and tissue-resident macrophage (P<0.01) were decreased in samples with high PYCR1 expression. A machine learning score was further built by random forest, involving PYCR1 and Trm markers. Only in the subgroup with the lower RFscore did IO+TKI show a favorable outcome, compared to TKI monotherapy.

CONCLUSIONS

Immunosuppression and IO+TKI resistance were correlated with high PYCR1 expression. T cell exhaustion and dysfunction were also related with the expression of PYCR1. PYCR1 has the potential to be employed as a biomarker to discriminate between IO+TKI and TKI monotherapy as the optimal patient treatment strategy.

Unfavorable immunotherapy plus tyrosine kinase inhibition outcome of metastatic renal cell carcinoma after radical nephrectomy with increased ADAM9 expression

Immunotherapy plus tyrosine kinase inhibitor (IO-TKI) has become the standard first-line therapy for advanced renal cell carcinoma (RCC). However, the modest response rate of IO-TKI therapy and the absence of biomarkers limited the selection of treatment strategies for RCC patients. There were three cohorts enrolled: two from our facility (ZS-MRCC and ZS-HRRCC) and one from a clinical study (JAVELIN-101). By RNA sequencing, the expression of ADAM9 in each sample was measured. By flow cytometry and immunohistochemistry, immune infiltration and T cell function were examined. Primary outcomes were established as treatment response and progression-free survival (PFS). Patients with low-ADAM9 expression had a higher objective response rate (56.5% vs 13.6%, P = 0.01) and longer PFS in both cohorts. In the ZS-HRRCC cohort, the expression of ADAM9 was associated with increased tumor-infiltrating T cells, which was proved by immunohistochemistry (P < 0.05) and flow cytometry (Spearman's ρ = 0.42, P < 0.001). In the high-ADAM9 group, CD8⁺ and CD4⁺ T cells revealed an exhausted phenotype with decreased GZMB (Spearman's ρ =  - 0.31, P = 0.05, and Spearman's ρ =  - 0.49, P < 0.001, respectively), and fewer Macrophages were identified. A predictive RFscore was further constructed by random forest approach, involving ADAM9 and immunologic genes. Only in the subgroup with the lower RFscore did IO-TKI outperform TKI monotherapy. High-ADAM9 expression was associated with immunosuppression and IO-TKI resistance. Expression of ADAM9 was also associated with the exhaustion and dysfunction of T cells. ADAM9-based RFscore has the potential to be used as a biomarker to distinguish the optimal patient treatment methods between IO-TKI and TKI monotherapy.

Alternative Complement Pathway Signature Determines Immunosuppression and Resistance to Immunotherapy Plus Tyrosine Kinase Inhibitor Combinations in Renal Cell Carcinoma

BACKGROUND

Latest guidelines recommended immunotherapy (IO) plus tyrosine kinase inhibitor (TKI) combination as standard first-line therapy in renal cell carcinoma (RCC), with no predictive biomarker being applied. Complement system shapes tumor microenvironment, which may influence TKI+IO benefit.

METHODS

Two cohorts from our institute and 2 external cohorts were enrolled. RNA-sequencing was performed for each sample, and alternative complement pathway signature (ACPS) was defined by single sample gene set enrichment analysis. Immune infiltration and function were assessed by immunohistochemistry and flow cytometry.

RESULTS

Under TKI+IO therapy, ACPS was elevated in non-responders (P<0.01), and high-ACPS predicted lower response rate and shorter progression-free survival (P=0.040). Moreover, TKI+IO, rather than TKI monotherapy, may benefit patients of low-ACPS combined with SETD2-wild type (HR=0.55, P<0.001). In RCC, ACPS was associated with increased tumor-infiltrating T cells (Spearman's ρ=0.50, P=0.001). However, in high-ACPS samples, CD8+ T cells revealed an exhausted phenotype with decreased GZMB (P<0.001) and increased PD1 (P=0.008) expression. Elevated PD1 expression in high-ACPS samples was confirmed by immunohistochemistry (P=0.046). Besides, macrophage infiltration was increased in high-ACPS samples (P=0.045), along with suppressive cytokines.

CONCLUSIONS

Under TKI+IO, high-ACPS was linked to immunosuppression and treatment resistance. ACPS might be used as a biomarker for better treatment strategy between TKI+IO or TKI monotherapy.

Properties and functions of myeloperoxidase and its role in ovarian cancer

Background: Elevated levels of myeloperoxidase in body fluids are increasingly being used as an indicator for the diagnosis of cancer. Aim of the study: The aim of this study was to review the literature on the physical and chemical properties of myeloperoxidase, its role in carcinogenesis, the role of tumor-associated neutrophils in cancer, and the role of myeloperoxidase in ovarian cancer. Material and methods: The research literature published between January 1999 and December 2019 was reviewed. The properties and role of myeloperoxidase in the development of ovarian cancer were selected from publications available in selected online databases, including MEDLINE, PubMed, Scopus, and Web of Science. Searches were performed using the following word combinations: “myeloperoxidase”, “ovarian cancer”, “reactive oxygen species”, “expression”, “polymorphism”, and “tumor-associated neutrophils”. Results: Thirty-five scientific articles were included in the final review. Of the 35 articles, 11 discussed the role of myeloperoxidase in carcinogenesis, and five discussed its role in the development of ovarian cancer. Conclusions: Elevated myeloperoxidase levels are associated with many types of cancer, including ovarian cancer. In the studied group of invasive ovarian tumors, up to 65% exhibited elevated levels of myeloperoxidase. Continued research on myeloperoxidase expression in ovarian cancer cells is vital and warranted.

Macrophage Phenotype in Combination with Tumor Microbiome Composition Predicts RCC Patients’ Survival: A Pilot Study

The identification of new prognostic markers of renal cell carcinoma (RCC) is an urgent problem in oncourology. To investigate the potential prognostic significance of tumor microbiome and stromal inflammatory markers, we studied a cohort of 66 patients with RCC (23 clear cell RCC, 19 papillary RCC and 24 chromophobe RCC). The microbiome was analyzed in tumor and normal tissue by 16S rRNA amplicon sequencing. Characterization of the tumor stroma was performed using immunohistochemistry. A significant difference in alpha diversity was demonstrated between normal kidney tissue and all types of RCC. Further, we demonstrated that the bacterial burden was higher in adjacent normal tissue than in a tumor. For the first time, we demonstrated a significant correlation between bacterial burden and the content of PU.1+ macrophages and CD66b+ neutrophils in kidney tumors. Tumors with high content of PU.1+ cells and CD66b+ cells in the stroma were characterized by a lower bacterial burden. In the tumors with high bacterial burden, the number of PU.1+ cells and CD66b+ was associated with a poor prognosis. The identified associations indicate the great prognostic potential of a combined tumor microbiome and stromal cell analysis.

Prognostic Significance of NLR About NETosis and Lymphocytes Perturbations in Localized Renal Cell Carcinoma With Tumor Thrombus

Background

Non-metastatic renal cell carcinoma (RCC) with tumor thrombus showed a greater tendency for developing metastases after surgery. Early identification of patients with high risk of poor prognosis is especially important to explore adjuvant treatment of improving outcomes. Neutrophil-to-lymphocyte ratio (NLR) was a systemic inflammation marker and outcome predictor in RCC, reflecting the chaos in systemic immune status in cancer as myeloid cell expansion and lymphatic cell suppression. Neutrophil extracellular traps (NET) formation (NETosis) is the process of neutrophils generating an extracellular DNA net-like structure. NETosis in tumor was demonstrated to conduce to the subsequent metastases of tumor. However, the role of NLR for systemic immune status and tumor local immune infiltration, especially for neutrophil-associated NETs, in non-metastatic RCC with thrombus remains unclear.

Patients and Methods

In our clinical cohort, we enrolled the clinical, pathologic, and preoperative laboratory parameters of 214 RCC patients with tumor thrombus who were treated surgically. The clinical endpoint was defined as cancer-specific survival (CSS). In our basic research cohort, RNA-seq, TCR-seq, and scRNA-seq data were analyzed. Patients who reached the endpoint as recurrence-free survival (RFS) were defined as the “High-risk” group. Otherwise, they were separated into the “Low-risk” group.

Results

In the clinical cohort, NLR≥4 was an independent risk factor for 203 localized RCC with tumor thrombus. In the basic research cohort, tumor thrombi were separated into NETosis-thrombi belonging to the “High-risk” group and non-NETosis-thrombi to the “Low-risk” group. NETs induced by tumor-derived G-CSF in tumor thrombus has a mechanistic role in unfavorable prognosis. Besides, NETs-score from single sample GSEA (ssGSEA) algorithm was an independent prognostic factor validated in the TCGA data. Apart from the neutrophils-associated NETosis, systemic immune perturbations of lymphocytes occurred in the “High-risk” group, represented with decreased TCR diversity and increasingly high proportion of CD4-positive effector memory T (Tem) cells, which indirectly represented the state of lymphopenia.

Conclusions

Our findings firstly demonstrated that neutrophils-associated NETosis and systemic lymphocytes perturbations were considered as tumor progression in patients of localized RCC with tumor thrombus, which reflected NLR≥4 as an independent risk factor for patients.

Identification of a Novel Glycolysis-Related Gene Signature Correlates With the Prognosis and Therapeutic Responses in Patients With Clear Cell Renal Cell Carcinoma

Background

Accumulating evidences indicate significant alterations in the aerobic glycolysis in clear cell renal cell carcinoma (ccRCC). We aim to develop and validate a glycolysis-related genes signature for predicting the clinical outcomes of patients with ccRCC.

Methods

mRNA expression profiling of ccRCC was obtained from The Cancer Genome Atlas database. Univariate Cox regression analysis and lasso Cox regression model were performed to identify and construct the prognostic gene signature. The protein expression levels of the core genes were obtained from the Human Protein Atlas database. We used four external independent data sets to verify the predictive power of the model for prognosis, tyrosine kinase inhibitor (TKI) therapy, and immunotherapy responses, respectively. Finally, we explored the potential mechanism of this signature through gene set enrichment analysis (GSEA).

Results

Through the GSEA, glycolysis-related gene sets were significantly different between ccRCC tissues and normal tissues. Next, we identified and constructed a seven-mRNA signature (GALM, TGFA, RBCK1, CD44, HK3, KIF20A, and IDUA), which was significantly correlated with worse survival outcome and was an independent prognostic indicator for ccRCC patients. Furthermore, the expression levels of hub genes were validated based on the Human Protein Atlas databases. More importantly, the model can predict patients’ response to TKI therapy and immunotherapy. These findings were successfully validated in the external independent ccRCC cohorts. The mechanism exploration showed that the model may influence the prognosis by influencing tumor proliferation, base mismatch repair system and immune status of patients.

Conclusions

Our study has built up a robust glycolysis-based molecular signature that predicts the prognosis and TKI therapy and immunotherapy responses of patients with ccRCC with high accuracy, which might provide important guidance for clinical assessment. Also, clinical investigations in large ccRCC cohorts are greatly needed to validate our findings.

Intratumoral CXCL13+CD8+T cell infiltration determines poor clinical outcomes and immunoevasive contexture in patients with clear cell renal cell carcinoma

Background

Chemokine (C-X-C motif) ligand 13 (CXCL13) was known as a selective chemotaxis for B cells, a product of follicular helper CD4⁺T cells (T_FH) and a contributor to tertiary lymphoid structures (TLS). Although secretion and function of CXCL13 produced by T_FH have been deeply explored, the immune function and prognostic significance of CXCL13 secreted by CD8⁺T cells still remain unrevealed. This study aims to investigate the clinical merit of CXCL13⁺CD8⁺T cells in clear cell renal cell carcinoma (ccRCC).

Methods

We analyzed prognostic value and immune contexture that associated with CXCL13⁺CD8⁺T cells infiltration level in a total of 755 patients from Zhongshan Hospital cohort (n=223) and The Cancer Genome Atlas cohort (n=532). In vitro analyses were conducted on 42 samples of resected tumor tissue from Zhongshan Hospital in order to detect the immune status of CXCL13⁺CD8⁺T cells and total CD8⁺T cells. Immunohistochemistry (IHC) and flow cytometry were applied to characterize immune cells and portray the tumor microenvironment (TME) in ccRCC.

Results

Intratumoral CXCL13⁺CD8⁺T cells abundance was associated with inferior overall survival and disease-free survival. CXCL13⁺CD8⁺T cells possessed higher level of immune checkpoints like programmed cell-death protein 1 (PD-1), T-cell immunoglobulin mucin 3 (Tim-3), T cell immunoreceptor with Ig and ITIM domains (TIGIT) and cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), higher Ki-67 expression and lower tumor necrosis factor α (TNF-α), interferon γ (IFN-γ) expression. Total CD8⁺T cells in high-level CXCL13⁺CD8⁺T cells infiltration subgroup exhibited elevated exhausted markers (PD-1, Tim-3, TIGIT) and descended activated markers (TNF-α, IFN-γ) without quantity variance. Furthermore, the abundance of intratumoral CXCL13⁺CD8⁺T cell was correlated with immunoevasive TME accompanied by increased T helper 2 cells, tumor-associated macrophages, Foxp3⁺ regulatory T cells, TLS and decreased natural killer cells, GZMB⁺ cells.

Conclusions

Intratumoral CXCL13⁺CD8⁺T cells infiltration indicated inferior clinical outcome in patients with ccRCC. CXCL13⁺CD8⁺T cells possessed increased exhausted markers, decreased effector molecules and better proliferation ability. CXCL13⁺CD8⁺T cells abundance impaired total CD8⁺T cells’ immune function. Intratumoral CXCL13⁺CD8⁺T cells abundance was associated with immunoevasive contexture. The abundance of CXCL13⁺CD8⁺T cells was an independent prognosticator and a potential immunotherapeutic target marker for ccRCC treatment.