PTP4A3 Is a Prognostic Biomarker Correlated With Immune Infiltrates in Papillary Renal Cell Carcinoma

PTP4A3 plays an important role in the tumorigenesis and metastasis of multiple tumors, but its prognostic role in renal cancer is not well understood. We utilized the Oncomine and Tumor Immunoassay Resource databases to examine the differential expression of PTP4A3 in tumor tissues and normal tissues in breast, urinary tract, gastrointestinal tract and skin. Using the GEPIA and PrognoScan databases, the independent prognostic role of PTP4A3 was confirmed in clear cell renal cell cancer and papillary renal cell cancer. Expression of PTP4A3 were obviously higher in tumor tissue compare with normal tissues (P=0.028). We haven’t found the associations of PTP4A3 and clinicopathological features in our IHC cohort. Ectopic expression of PTP4A3 promotes proliferation, migration and invasion and increased the mRNA level of TGFB1 in RCC cell lines. Immunohistochemical staining indicated that the expression of PTP4A3 associates with CD3+ (P =0.037)/CD8+ (P =0.037) intratumor TILs, not with invasive margins in renal cancer. Comprehensive analysis of immune infiltration in the TIMER database correlated PTP4A3 expression with the infiltration of B cells, CD8+ T cells, CD4+ T cells and neutrophils in both clear cell renal cell carcinoma and papillary renal cell carcinoma. PTP4A3 expression was associated with the infiltration of dendritic cells in papillary renal cell carcinoma. We further confirmed that the infiltration of B cells and CD8+ T cells was associated with poor prognosis in papillary renal cell carcinoma patients, consistent with the prognostic role of PTP4A3 in papillary renal cell carcinoma. PTP4A3 expression correlated genes involved in B cells, monocytes, M1 macrophages, Th2 and Treg cells in papillary renal cell carcinoma. These results suggest PTP4A3 as a prognostic factor with a role in regulating immune cell infiltration in papillary renal cell carcinoma.


INTRODUCTION
Renal cell carcinoma (RCC) is a common malignant disease expected to account for approximately 76,080 cases and 13,780 deaths in 2021 in the United States (1). RCC is divided into three major histologic subtypes: clear cell RCC (ccRCC or KIRC), papillary cell RCC (pRCC or KIRP), and chromophobe RCC (chRCC). According to The Cancer Genome Atlas (TCGA), a comprehensive database of cancer genomic profiles, these three major RCC subtypes have distinctive genetic alterations including genomic mutations, methylation status, RNA expression signatures, and immune signatures (2)(3)(4). At diagnosis, 30% of RCC patients were confirmed as distant metastatic (5). While localized RCC could be cured by surgical resection, about 40% of surgical patients eventually experienced recurrence. Recently, adjuvant therapy, and especially immune therapy, has received significant attention as an option for treating metastatic RCC (mRCC). Immune therapy combined with a tyrosine kinase inhibitor; bevacizumab, a monoclonal antibody against vascular endothelial growth factor; or ipilimumab, a monoclonal antibody against cytotoxic T-lymphocyte antigen-4 are now considered optimal treatment strategies for mRCC. The biomarkers Programmed Cell Death Protein 1 (PD1) and Programmed Death-Ligand 1 (PD-L1) exhibited higher expression in ccRCC than in chRCC or pRCC, and high expression was correlated with decreased survival in all RCC patients (P=0.02) (6,7). Immune cell genes signatures, including those of B cells, dendritic cells, CD68, CD8, Treg cells, Th1/2 and immunoglobin G, were examined in ccRCC, pRCC and chRCC. In general, expression was higher in ccRCC than pRCC and chRCC, except for Th17, IL-8, and CD56 bright natural killer (NK) cell signatures. Th2 gene signatures were universally correlated with poor overall survival (OS) in all major histological subtypes (2). IL-8 inflammatory pathwayrelated genes and CD56 bright NK cell genes had higher expression in pRCC than in ccRCC and chRCC, whereas IL-17 producing T helper cells (T H 17 cells) had lower expression in pRCC than in ccRCC and chRCC. T H 17 gene signatures positively correlated with OS in ccRCC and chRCC (P=0.0021 and P=0.0362, respectively). CpG island methylator phenotypeassociated RCC had increased expression of T H 2 cells and Treg cells compared with pRCC and chRCC. Therefore, immune biomarkers for different subtypes of RCC need to be explored.
Due to the important role of PTP4A3 in cancers, a PRLzumab was generated by the Zeng group and was shown to specifically inhibit PTP4A3-positive liver cancer cell in vivo, and recruit B cells, NK cells and macrophages to the PTP4A3positive tumor microenvironment (27,28). A DNA vaccine targeting PTP4A3 triggered high expression of interferon-g and TNF-a in breast cancer through the CTL and T helper type 1 cells immune response and also stimulated the accumulation of PTP4A3 antibody in immunized mice (29). IL-6/8 secreted by tumor-associated macrophages (TAMs) facilitated the metastasis of colon cancer in a PTP4A3-KCNN4-dependent manner (30). Few reports have examined the role of PTP4A3 in kidney cancer, and this research has mostly focused on the correlation of clinicopathological characteristics and prognosis. The mechanism of how PTP4A3 regulates immune cells and its impact on renal cancer prognosis remains unclear.
Herein, we developed a comprehensive analysis across multiple databases to elucidate the role of PTP4A3 in ccRCC and pRCC. In addition, we analyzed the relationship of PTP4A3 and immune cell tumor infiltration (TIMER) (31,32). This study offers new insights into the function of PTP4A3 in renal cancer and proposes a mechanism for PTP4A3 regulating tumor infiltration of immune cells.

Oncomine Database Analysis
The Oncomine database comprises a wide range of gene expression datasets covering various cancer types. PTP4A3 expression in diverse tumor types was investigated using this database (https://www.oncomine.org/resource/login.html) (33)  with the screening thresholds set as: P-value of 0.001, fold change of 1.5 and gene rank of top 5%.

PrognoScan Database Analysis
The PrognoScan database complies accessible tumor microarray datasets, facilitating the correlation of gene expression and survival including OS and disease-free survival (DFS). We evaluated the relationship between PTP4A3 expression and prognosis in diverse tumor types using the PrognoScan database (http://www.abren.net/PrognoScan/) (34) with a threshold P-value of 0.05.

Lentivirus Infection and Western Blot
To stable expression of ectopic PTP4A3 in renal cancer cell, 786-O and Caki-2 cells were infected with 50 MOI control or myc-PTP4A3 expression lentivirus for 72 hours. Cells were directly homogenized in 2 × loading buffer to get the whole cell extracts. Protein samples were separated by the SDS-polyacrylamide gel electrophoresis and transferred onto the nitrocellulose membranes, and blocked in 5% milk in PBS buffer for 1 hour. The lentivirus infection were detected by the myc-tag antibody (1mg/ml) at 4°C overnight and washed in 0.1% Tween in PBS for 3 times. After incubating with HRP-conjugated anti-mouse secondary antibody, the protein bands were visualized by the enhanced chemiluminescence detection system (Thermo Fisher Scientific, Pittsburgh, PA, USA).

Quantitative Real-Time PCR
Total RNA was extracted from cells with Trizol reagent (Invitrogen, Carlsbad, CA, USA) based on the manufacturer's instruction. 1 mg of RNA was utilized to synthesize cDNA with GoScriptTM Reverse Transcription system (A5001, Promega). qRT-PCR was conducted with SYBR Green PCR master mix reagents (TOYOBO) and a StepOne Real-time PCR system (Applied Biosystems). Expression data of indicative gene was normalized to that of GAPDH. Primers used are listed in Supplementary Table 4.

Immunohistochemical Analysis
The kidney cancer tissue array (DC-Kid11051) were bought from Avilabio company (Shanxi). DC-Kid11051 were including 33 KIRC, 9 KIRP, 5 KICH and 6 normal kidney tissues. The sections were deparaffinized in xylene and hydrate in alcohol. Endogenous peroxidase activity was then blocked by incubation in 3% hydrogen peroxide-methanol for 10 min. For the antigen recovery, the section was heated in a citrate buffer (pH 6.0) for 15mins. Incubation with PTP4A3 antibody (11) (5mg/ml, monoclonal antibody 3B6 previous validated), CD3 antibody ((Monoclonal Rabbit Anti-Human CD3, Clone 2GV6, prediluted, Roche) and CD8 antibody (Monoclonal Mouse Anti-Human CD8, Clone C8/144B, Dako) at 4°C overnight. EnVision + TM (Dako, Carpinteria, CA, USA) was used as the secondary antibody. Antibody binding was visualized by a standard streptavidin immunoperoxidase reaction. Immunoreactivity in the cytoplasma and cytoplasmic membrane of PTP4A3 was evaluated.
Evaluation immunoreactivity was carried out independently by three experienced pathologists without any knowledge of the clinical data. The IHC score was assessed according to the percentage of positive cells and the intensity of cytoplasmic reactivity. These averaged immunoreactivity values were stratified into five scoring groups: -, not detected; ±, <10% positive cells; +, 10-20% weakly to moderately positive cells; ++, 10-20%intensely positive cells or 20-50% weakly positive cells; and +++, 20-50% positive cells with moderate to marked reactivity or >50% positive cells. In the statistical analysis, ± were considered negative, + and above were considered positive.

Quantification of CD3+ and CD8+TILs
Each entire section was evaluated for CD3+ and CD8+ TILs by two pathologists blind to PTP4A3 expression. The intratumor and invasive margins were identified by the immunohistochemistry. Each positive cell was manually counted under 40 ×magnification field in 5 independent tumor fields or invasive margins with most abundant TILs. The consecutive sections were stained for CD3 and CD8.

TIMER Database Analysis
Tumor Immune Estimation Resource (TIMER, https://cistrome. shinyapps.io/timer/) was developed to evaluate the abundance of immune infiltrates using 10,897 cancer samples from TCGA (31,32). TIMER gene modules were utilized to investigate PTP4A3 expression in diverse tumors and the relationship between PTP4A3 expression and the abundance of immune infiltrates. In addition, correlation modules were used to evaluate the relationship of PTP4A3 expression with gene biomarkers of tumor-infiltrating immune cells. These gene biomarkers have been previously studied (35)(36)(37).

GEPIA Database Analysis
We confirmed the correlation of PTP4A3 expression and tumor infiltration of immune cells with the Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/index. html) (38). We also assessed the prognostic value, including OS and DFS, of gene expression in diverse tumors with this database. Meanwhile, gene expression correlation was investigated based on TCGA expression data. Correlation analyses were conducted in both cancer and paired normal tissues. The correlation coefficient was calculated by the Spearman method. PTP4A3 was represented on the X-axis, and related biomarker genes were utilized for the Y-axis.

Statistical Analysis
Results from the Oncomine database included fold change, P-value, and gene rank. The PrognoScan and GEPIA databases produced HR and P-values or Cox P-values according to a logrank test. The correlation between PTP4A3 expression and other related immune genes was investigated using the Spearman method and the correlation strength was categorized according to r values. R values of 0.80-1.00 were deemed very strong, r values of 0.60-0.79 were strong, r values of 0.40-0.59 were moderate, r values of 0.20-0.39 were weak, and r values of 0.00-0.19 were categorized as very weak. P-values less than 0.05 were considered statistically significant. Chi-square test and nonparametric test were utilized to analyze the Categorical data by SPSS 19.0 software (SPSS Inc., Chicago, IL, USA).

PTP4A3 mRNA Expression Level in Diverse Tumors and Adjacent Normal Tissues
We investigated differential PTP4A3 mRNA levels in tumor and normal tissues across various cancer types in the Oncomine database. Compared with adjacent normal tissues, PTP4A3 mRNA levels were elevated in most cancer types, including breast, colorectal, esophageal, head and neck, kidney, leukemia, liver, melanoma, myeloma, pancreatic, prostate, and sarcoma. However, PTP4A3 mRNA levels were lower in kidney and sarcoma cancer in other datasets ( Figure 1A). PTP4A3 mRNA levels in diverse tumors are detailed in Supplementary Table 1.
PTP4A3 expression in diverse tumor types was further assessed using the TCGA and TIMER databases. Differential PTP4A3 expression in tumors compared with adjacent normal tissues in the TCGA datasets is presented in Figure 1B  thyroid carcinoma (THCA) relative to adjacent normal tissues (P<0.001). PTP4A3 expression was also slightly elevated in chRCC and stomach adenocarcinoma (STAD) compared with adjacent normal tissues (P<0.01). However, PTP4A3 expression was significantly lower in urothelial carcinoma (BLCA), and squamous cell lung carcinoma (LUSC) compared with normal tissues. Notably, PTP4A3 expression was significantly higher in head and neck squamous cell carcinoma (HNSC) HPV-positive cancer tissues compared with HNSC HPV-negative cancer tissues ( Figure 1B). Single cell RNA sequencing has shown that HNSC HPV-positive and HNSC HPV-negative patients have distinct immune gene signatures, indicating that PTP4A3 may be correlated with immune-related genes (39,40).

Prognostic Value of PTP4A3 in Multiple Cancer Types
To better understand the impact of PTP4A3 expression on clinical outcomes, the association between PTP4A3 expression and prognosis in 33 cancer types was analyzed using the GEPIA database. OS and DFS in urothelial cancer are presented in Figure 2.    Supplementary  Figures 2A-L). These results indicated that high PTP4A3 expression is significantly associated with poorer clinical outcomes in many tumors types, including KIRC and KIRP.

PTP4A3 Plays as an Oncogene and Regulates the Cytokines in RCC
To further explore the mechanism of PTP4A3 in RCC, we stained the PTP4A3 in the kidney cancer tissue array which containing 5 different tumor histological subtypes and 6 kidney normal tissues. As shown in Figures 3A, B, PTP4A3 expression is obviously higher in the kidney tumor tissues (n=49) compare with the normal tissues (n=6) (P=0.028). The representative high expression of PTP4A3 in KIRC, KIRP and KICH, kidney squamous cancer and kidney normal tissues (Figures 3Aa-e). And there was no significantly difference in the various histological subtype, such as KIRC, KIRP and other types. (P=0.538, Table 1). We also calculated the association of PTP4A3 in general clinicopathological features, there was no significant difference found in gender, age, Fuhrman grade, TNM stage and lymph node metastasis ( Table 1). Notably, PTP4A3 expression was slightly higher in the Fuhrman grade 3 compare with the expression in Grade 1 and 2(77.8% VS 40%, P=0.139). Although there was not rigorous statistic significant, PTP4A3 expression probably indicated the high grade, high recurrence and poor survival.
Results from the GEPIA and Prognoscan database also shown that PTP4A3 may play as an oncogene in the RCC. We Overexpression of PTP4A3 also triggered the metastasis phenotype of 786-O or Caki-2 cell lines ( Figures 3E, F). Stable expression of PTP4A3 elevated the secretion level and mRNA level of IL-1alpha in colon cancer cell lines. And IL-1RN (IL-1 alpha inhibitor) decreased the migration and invasion of PTP4A3-stable expression cells (22). Therefore, we detected the cytokines (CD115, CCR8, IL-1 alpha, MMP-11 and TGFB1) in 786-O and Caki-2 cell lines with PTP4A3 stable expression. As shown in Figure 3G

PTP4A3 Expression Is Associated With CD3+ and CD8+ Intratumor T Cells in Renal Cancer
TGFB1 expressing and FOXP3 expressing on Treg cell is important in the maintenance of immune homeostasis. Antitumor effect of PTP4A3-target monoclonal antibody was described in gastric cancer (27). Immunotherapy combined with target therapy began the new era of malignant tumor such as lung cancer, gastric cancer (41)(42)(43). Considering tissue array is not suitable for the research of immune cell infiltration, we stained the CD3 (general T cell markers) and CD8+ T cells in KIRC surgery sections (n=10) to further explore the role of PTP4A3 in immnuo-target combined therapy. The representative IHC stains were shown in Figure 4A. Each count of CD3+/CD8+ intratumor and CD3/ CD8+ invasive margins were shown in Figure 4B Figure 4C and Table 2). Therefore, PTP4A3 associated with low CD3+/CD8+TILs and indicated poor prognosis in in renal cancers. Figure 4 showed PTP4A3 may associate with T cell infiltration in renal cancer. Renal cancer samples were not easy to collect, especially the relative low incidence subtype. Therefore, we utilized the TIMER database to comprehensive explore the relation of PTP4A3 and immune infiltration. Previous research has revealed that tumorinfiltrating lymphocytes could serve as an independent prognostic predictor in various cancers (44,45). Using the databases (PrognoScan and GEPIA) and cell model, oncogene PTP4A3 expression was found to be correlated with decreased survival in KIRC and KIRP and promoted cell proliferation, migration and invasion (Figures 2 and 3).

Preliminary experiments in
Consistent with the unique immune gene profiles of renal cancer subtypes, the correlation of PTP4A3 expression with immune cell infiltration was distinct in KIRP and KIRC. In KIRC, significant correlations were observed between PTP4A3  Kaplan-Meier curves showing immune cell infiltration and PTP4A3 expression in KIRC and KIRP were generated using the TIMER database. B cell infiltration (P=0.035), CD8+T cell infiltration (P=0.024), and PTP4A3 expression (P=0.017) were significantly associated with survival in KIRP ( Figure 5B). However, there was no significant relationship between immune cell infiltration (including B cells, CD8+T cells, CD4+ T cells, neutrophils, macrophages, and dendritic cells) and survival in KIRC. Taken together, these results suggest that PTP4A3 is involved with regulating immune cell infiltration in KIRP and KIRC and, along with B cell infiltration and CD8+ T cell infiltration, influences clinical outcomes in KIRP.

Correlation Between PTP4A3 Expression and Immune Marker Sets
To further evaluate the association between PTP4A3 expression and immune cell infiltration, the TIMER and GEPIA databases were utilized to investigate the relationship between PTP4A3 and immune marker sets in KIRP and KIRC.

DISCUSSION
PTP4A3 is a 22kD tyrosine phosphatase and is highly homologous with PTP4A1 and PTP4A2. While PTP4A1 is expressed ubiquitously, PTP4A2 and PTP4A3 are specifically expressed in heart and skeletal muscles. PTP4A3 is a dual phosphatase and is known to act on ezrin, cytokeratin and integrin b1 (46)(47)(48). PTP4A3 is involved in the regulation of cell proliferation, migration and invasion through the PI3K-AKT-ERK pathway and destabilizes telomeres through the NFкB-RAP1-TRF2 axis (18,23,24,49). Moreover, PTP4A3 contributes to the secretion of inflammatory factors, which recruit more immune cells to the tumor microenvironment. Activation of PTP4A3 expression increases IL-1 alpha secretion through the NF-kB and   (27,50). Therefore, PTP4A3 could be a potential target for cancer immunotherapy. Clinical trials have shown that immunotherapy in renal cell carcinoma is superior to tyrosine kinase inhibitors, heralding a new era of renal cancer treatment.
As there is no study showing PTP4A3 expression and its potential roles in renal cancer, further experiments were conducted to evaluate whether PTP4A3 play as an oncogene in RCC. In our validation cohort, PTP4A3 expression is significantly elevated in the tumor tissues by the immunohistochemical staining of PTP4A3 in 49 kidney cancers and 6 normal tissues. As for the functional study, ectopic expression of PTP4A3 promotes proliferation, migration and invasion in RCC cell lines. Taken collectively, the above findings firstly suggested that PTP4A3 serve as an oncogene in renal cancer. The metastasis phenotype could be partially inhibited by IL-1alpha inhibitors, and MMP11 and FOXN3 were elevated in the PTP4A3 overexpressing colon cells (22). These immune related cytokines may connect the oncogene PTP4A3 with the immune infiltration. We also tried to verify the correlation of PTP4A3 and IL-1 alpha in RCC cells, but the mRNA level of IL-1 alpha showed different changes in 786-O and Caki-2 cell. Considering the possible inconsistent of mRNA level and secretion level in supernatant, the cytokines in supernatant should be detected by the ELISA to further confirm the correlations. Luckily, we found that PTP4A3 upregulated the mRNA of TGFB1 in both 786-O and Caki-2 cells. TGFB1 expressed on the Treg cell surface and maintain the immune hemostasis. The correlation of PTP4A3 and TGFB1 were also been confirmed in Table 3 by the TIMER database (R=0.370, p<0.001). Therefore, it indicated that PTP4A3 may also regulate the immune microenvironment through cytokines TGFB1 in RCC.
In this study, we focused on the association between PTP4A3 and cytotoxic T cells, which have not been jointly reported. CD3+, CD8+ T cells represent T cells, cytotoxic T cells. Our explorations examined the immune infiltrates' density and the location of immune cell populations in intratumor and invasion margins. In KIRC, we observed negative association between PTP4A3 expression and the density of CD3+ TILs intratumor and CD8+ TILs intratumor in the tumor microenvironment. There was no significant relationship between PTP4A3 expression and the number of CD3+ TILs invasion margins and CD8+ TILs invasion margins. However, we did not assess whether the combination of PTP4A3 and CD3+ and CD8+   (51,52). These findings suggested that immune infiltrates of T cells play a protective role. The role of CD8+ T cells in renal cancer is controversial, maybe attribute to technical factors such as antibodies used to examine CD8+ T cells and different tumor location detected (intratumor vs invasion margins TILs) (52,53). In general, the tumor-infiltrating T cells play vital roles in renal cancer, including cytotoxic T cells, memory T cells, regulatory T cells, and so on (54,55). In the future work, the relationship between PTP4A3 expression and B cells and other T cells including CD4+ T cells, CD45RO+ memory T cells, regulatory T cell need to be explored in renal cancer.
In this study, the role of PTP4A3 in renal cell carcinoma immune microenvironments was explored using online databases, including the Oncomine, TCGA, GEPIA, TIMER and Progscan databases. Consistent with previous findings, PTP4A3 expression was elevated in many tumor types compared with normal tissues, including renal cancer. Kaplan-Meier survival analysis suggested that PTP4A3 expression was correlated with poor DFS but not with OS in KIRC and KIRP. This may be related to the cutoff percentage setting (50% cut-off) and case numbers in the GEPIA database. PTP4A3 expression is higher in gastric cancer than in colon cancer, and the cutoff settings in various cancer types should be adjusted base on the immunohistochemistry of each cancer type. The association of PTP4A3 expression with B cells, CD8+ T cells, CD4+T cells and neutrophil infiltration was confirmed in renal cancer ( Figure 5 and Supplementary Table 3) using the GEPIA and TIMER databases. These results were consistent with PRL-3-zumab work from Zeng et al. (50). PTP4A3 is an oncogene related to immune cell infiltration in KIRC and KIRP.
Tumorigenesis and tumor development are complicated processes involving thousands of genes and pathways. Because of this, signal gene signatures are rarely used to predict tumor prognosis. Due to the association of PTP4A3 expression with immune cell infiltration, the relationship of immune markers and PTP4A3 expression was analyzed using the TIMER database to identify gene profile signatures that may be able to predict prognosis and immunotherapy response in KIRP. Correlations were observed between PTP4A3 expression and M1 macrophage immune markers such as INOS and COX2 in both KIRP and KIRC (Supplementary Table 3 and Table 3), suggesting that PTP4A3 may regulate the polarization of M1 macrophages in renal cancer. A positive correlation was observed between PTP4A3 expression and B cell immune marker (CD79A) in KIRP only (Supplementary Table 3 and Table 3), and this result was confirmed with the TIMER database. Dendritic cells can play a vital role in tumor growth and metastasis by inhibiting CD8+ T cell cytotoxicity and promoting Treg cells (56). PTP4A3 expression was correlated with HLA-DRA, BDCA-1, and BDCA-4 in dendritic cells, and with CCR8 and TGFB1 in Treg cells, indicating that PTP4A3 promoted tumor growth and metastasis through the classical PI3K-AKT-PTEN signal pathway and also by regulating the tumor immune microenvironment (18)(19)(20). Interestingly, exhausted T cell markers such as PD1, LAG3 and GZMB were also significantly correlated with PTP4A3 expression in KIRP (Supplementary Table 3). PD1 is the immune checkpoint receptor expressed on the activated T cells; blocking the PD1/PD-L1 interaction would release immune suppression and trigger immune clearance in tumors. This study is the first reported association of PTP4A3 and PD1 in renal cancer. The combination of PTP4A3 with immune cells is proposed as a predictor for immunotherapy response and prognosis in renal cancers, with potential applications in other cancer types such as NSCLC or gastric cancer as well (57,58). These findings further confirm the relationship of PTP4A3 expression to immune infiltration in KIRP, suggesting that PTP4A3 influences immune escape in the renal cancer microenvironment. We are confidence in our results base on the various online datasets, therefore we are working on the collection of renal cancer samples and focusing on the immune-related role of PTP4A3 in renal cancer.
In conclusion, PTP4A3 may be a regulator in the tumor immune microenvironment and a useful biomarker for predicting prognosis and immunotherapy response in renal cancer. These results will be verified in our further workings.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the Ethics Committee of Peking University Cancer Hospital & Institute (Beijing, China). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

AUTHOR CONTRIBUTIONS
QS, ZL, and SL conceived the project and wrote the manuscript. JW and SW managed data acquisition. QF, JW, and SW participated in the data analysis. QS and SL participated in the discussion and language editing. ZL reviewed the manuscript. YZ did the IHC, cell proliferation migration, and invasion experiments during the revised of the manuscript. LZ and QY counted the CD3+ and CD8+ TILs. All authors contributed to the article and approved the submitted version.