Int J Med Sci 2020; 17(18):2917-2925. doi:10.7150/ijms.51126 This issue

Research Paper

The effect of CCL5 on the immune cells infiltration and the prognosis of patients with kidney renal clear cell carcinoma

Shuheng Bai1,3, YinYing Wu2, Yanli Yan1, Haojing Kang1, Jiangzhou Zhang3, Wen Ma3, Ying Gao1, Beina Hui1, Rong Li1, Xiaozhi Zhang1, Juan Ren1 Corresponding address

1. Department of Radiotherapy, Oncology Department, First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, Shaanxi Province, China, 710061;
2. Department of Chemotherapy, Oncology Department, First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, Shaanxi Province, China, 710061;
3. Medical School, Xi'an Jiaotong University Xi'an, Shaanxi Province, China, 710061

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Bai S, Wu Y, Yan Y, Kang H, Zhang J, Ma W, Gao Y, Hui B, Li R, Zhang X, Ren J. The effect of CCL5 on the immune cells infiltration and the prognosis of patients with kidney renal clear cell carcinoma. Int J Med Sci 2020; 17(18):2917-2925. doi:10.7150/ijms.51126. Available from https://www.medsci.org/v17p2917.htm

File import instruction

Abstract

Background: Kidney renal clear cell carcinoma (KIRC) is the most representative subtype of renal cancer. Immune infiltration was associated with the survival time of patients with tumors. C-C chemokine ligand 5 (CCL5) can promote the malignant process of tumor and be related to infiltration immune cells in some cancers, but not reported in KIRC.

Methods: The expression profile and clinical data were obtained from The Cancer Genome Atlas (TCGA) database. The correlation between the expression level of CCL5 and clinical features in KIRC was analyzed. Gene Set Enrichment Analysis (GSEA) was utilized to explore the functions and pathways of CCL5 in KIRC. Then, the analysis between the survival and immune infiltration cells was carried out, as well as the non-parametric tests between the CCL5 expression and the ratios of immune infiltration cells.

Results: The correlations between the expression levels of CCL5 in KIRC and clinical features including survival time, pathological stage, grade, and status of the patient, have been identified. Meanwhile, GSEA analysis has shown relationships between the expression of CCL5 and immune pathways. The immune infiltrated cells were correlated with the prognosis of KIRC, especially regulatory T cells (Tregs), mast cells, and dendritic cells. And Tregs was associated with the CCL5 expression.

Conclusion: The increased expression of CCL5 is related to poor prognosis and clinical features. Meanwhile, CCL5 is related to Tregs ratios and CCL5 may act as a typical chemokine to recruit Tregs in KIRC. CCL5 could be used as a biomarker for the prognosis prediction and a potential therapeutic target for patients with KIRC.

Introduction

Renal cancer is a common cancer, and its incidence rates in males and females are 5% and 3%, respectively[1]. Kidney renal clear cell carcinoma (KIRC) accounts for 70%-80% of renal cancer, which is the most representative subtype, with the increased incidence rate year by year. Early diagnosis of KIRC is more difficult compared with other cancers, because kidney cancer-related clinical symptoms are atypical. KIRC has poor responses to conventional chemotherapy and radiation therapy[2]. All of these lead to a low five-year survival rate of advanced patients, which is only 10% to 20% [3, 4]. Targeted therapy had a positive effect on prolonging the survival time of patients, including VEGF-TKIs, mTOR inhibitors, and monoclonal antibodies as Bevacizumab [5, 6]. But the drug-resistance due to long-term use has not been resolved. Novel treatment target is urgently needed. Immune therapy, especially immune checkpoint inhibitors, is a very promising treatment measure in patients with KIRC [7, 8]. But not all of the patients can benefit from it. Objective response rate to anti PD-L1 therapy is only around 20%. And patients who had a favorable response to immune checkpoint inhibitors did not show long time remission [7, 9, 10]. So, some new immune checkpoints are needed to be discovered urgently.

Chemokines are a type of small secreted protein with a molecular weight of 8-11KDa, composed of 70-125 amino acids. They are classified into CXC, CC, C and CX3C subfamily, according to the position and remaining structure of the first two cysteine residues in the N segment. And The N-terminal of CC subfamily chemokines contains two adjacent Cys residues. CCL5 belongs to the CC family of chemokine and is mainly expressed in T cells, macrophages and some tumor cells [11]. Some researche had demonstrated that CCL5 promotes the development of malignant tumors including lung cancer, colorectal cancer, esophageal cancer, breast cancer and so on [12-15]. So CCL5 may be utilized as a biomarker to predict prognosis and act as a new target to treat KIRC.

Immune infiltration in the tumor microenvironment has been demonstrated that it is not only associated with the survival time of cancer patients, but also has great potential for cancer treatment[16, 17]. Studies revealed some immune infiltration cells, including regulatory T cells, mast cells, CD8+T cells, and dendritic cells, had great effects on KIRC progression[18-21]. Moreover, some studies found that CCL5 is related to some immune diseases and regulates immune infiltrating cells which can yield immune escape [22-25]. Yet up to now, there are no studies which explore its prognostic value and its relationship with immune infiltration cells in KIRC.

In this study, Relationships between the mRNA expression levels of CCL5 and the clinical-pathological features, infiltration immune cell landscape of KIRC patients were analyzed based on TCGA database and bioinformatics analysis including CIBERSORT and GSEA. The purpose of this study was to evaluate the prognostic value of CCL5 and to explore the potential relationship between CCLR and infiltration immune cells in KIRC patients.

Methods

Data acquiring and related bioinformatics analysis

The transcriptome profiling data about Kidney Renal Clear Cell Carcinoma and corresponding clinical information were obtained from TCGA official website (https://cancergenome.nih.gov/), which is a public database demonstrated major cancer related genomic alterations [26]. Gene expression data (level 3) based on HTSeq-FPKM workflow, which including 539 tumor samples and 72 control samples (normal tissue adjacent to tumor in the same patient), and corresponding clinical data were obtained. Then boxplots were used to show the differences of genes expression between normal groups and tumor groups, and some clinical variables (stage, age, gender, pathologic grade, TNM). Finally, the survival curve (K-M) was applied to explore the relationships between CCL5 expression level and the overall survival time of patients with KIRC.

Gene set enrichment analysis

Gene set enrichment analysis (GSEA) is a method to determine that if predefined sets of genes are differentially expressed in different phenotypes. Predefined gene sets may be genes in a known metabolic pathway, located in the same cytogenetic band, sharing the same Gene Ontology category, or any user-defined set[27]. In this study, GSEA was carried out to explore the significant signal pathways related to CCL5 groups. Gene set permutations were performed 1000 times for each analysis. The c2.cp.kegg.v7.0.symbols.gmt data set was obtained from the GSEA website. The expression level of CCL5 was used as a phenotype label. The p value represents the credibility of the enrichment result, q value is the p value after multiple hypothesis testing correction. The nominal p value and normalized enrichment score (NES) were used to sort the pathways enriched in each phenotype.

Assessment of immune infiltration

CIBERSORT (http://cibersort.stanford.edu/) is an analytical tool to provide an estimation of the abundances of immune cells in mixed cell populations, using gene expression data [28]. It was used to evaluate the relative ratios of 22 types of infiltration immune cells in each cancer example, which including M1 macrophages, M2 macrophages, M0 macrophages, T follicular helper cells, resting memory CD4+ T cells, activated memory CD4+ T cells, γδT cells, CD8+ T cells, Tregs CD4+ T cells, naive CD4+ T cells, resting natural killer (NK) cells, activated NK cells, resting mast cells, activated mast cells, memory B cells, resting dendritic cells (DC), activated DC, naïve B cells, monocytes, plasma cells, neutrophils and eosinophils [29]. p‑value of CIBERSORT, which can evaluate the statistical significance of the deconvolution consequences, was used as set p<0.05.

Identification of the relationship between the profile of CCL5 and immune infiltration cells

Mann-Whitney test was used to evaluate the relationships between the expression of CCL5 and some types of immune infiltration cells which had effects on the overall survival time. The significant fitting accuracy was set as p<0.01.

Gene Expression Profiling Interactive Analysis (GEPIA) is a web-based tool to deliver fast and customized functionalities based on TCGA and GTEx data. It provides key interactive and customizable functions including differential expression analysis, profiling plotting, and the correlation analysis[30]. The “correlation” module was utilized to confirm the relationships between CCL5 expression and the possible gene markers of tumor-infiltrating immune cells which had been selected through the above process [16].

Statistical analysis

All statistical analyses were conducted by using R (v.3.6.1), GraphPad_Prism (V 8.0) and SPSS (V 25). Non-normally distributed data used Wilcoxon test and Chi-square test. Kaplan-Meier method was utilized to analyze the overall survival. Mann-Whitney test (U-test) was used to analyze the correlation between two groups.

Results

The relationship between the mRNA expression levels of CCL5 and the clinical features of patients with KIRC

The CCL5 mRNA levels of 539 KIRC tissues were significantly higher than that of 72 normal tissues (P<0.01) as shown in Figure 1A, and consistent results were showed in the 72 matched tissue samples from the KIRC patients (p<0.001) (Figure 1B). In the group of KIRC, CCL5 expression ranged from 0.87 to 269.103, and the median expression was 30.5, which was a criterion to define low and high expression groups. Survival analysis showed that the patients with high expression levels of CCL5 exhibited significantly shorter overall survival time (OS) (p=0.011) than the patients with low CCL5 expression levels, shown in Figure 1C.

Then, boxplots were performed to show the relationships between the expression of CCL5 and clinical characters including the age (p=0.373), gender (p=0.289), stage (p<0.01), pathologic T (p<0.01), pathologic N (p=0.039), pathologic M (p<0.01), grade (p<0.01), status (p<0.01), shown in Fig 1D-K.

 Fig 1 

The relationship between CCL5 and the clinical features. (A) Different expression profiles of CCL5 between 539 KIRC tissues and 72 normal tissues. (B) Different expression profiles of CCL5 between 72 KIRC tissues and matched normal tissue from the same patients. (C) Effect of CCL5 expression on OS in all patients, p=0.011. (D-K) Relationships between CCL5 expression and clinical characters of patients with KIRC including the age, gender, stage, grade, pathologic T, pathologic N, pathologic M, and status.

Int J Med Sci Image

(View in new window)

 Fig 2 

Enrichment plots from GSEA. (A) Chemokine signaling pathway. (B) Cytokine-cytokine receptor interaction pathway. (C) T cell receptor signaling pathway separately.

Int J Med Sci Image

(View in new window)

High expression levels of CCL5 of patients with KIRC were related to some clinical characters as shown in Table 1. Stage, pathological T stage, pathological M stage, the status of patients, tumor pathological grade were all significantly related with CCL5 high expression in patients with KIRC (all meet p<0.01); However, some clinical characters including age, gender, and pathologic N stage were not associated with the expression levels of CCL5 in patients with KIRC.

 Table 1 

The clinical characteristics of KIRC patients between high expression of CCL5 group and low expression group.

Clinical characteristicsVariableNo. patientsCCL5 expressionχ2P
low(<30.5)high(≥30.5)
Age<602461241220.030.862
≥60284141143
Genderfemale186101852.120.154
male344164180
Statusalive35719616110.51<0.01
dead17369104
GradeG11411335.46<0.01
G222713988
G320687119
G4752154
pathologic TT127116210921.28<0.01
T2692841
T317971108
T41147
pathologic NN02391191204.290.117
N116412
NX275142133
pathologic MM042122120010.25<0.01
M1782652
MX301713
pathologic Stagestage I26516010527.77<0.01
stage II572532
stage III1235271
stage IV822755

KEGG pathways related to CCL5

To identify the CCL5 related signaling transduction pathways in KIRC, GSEA analysis was performed between the data of high CCL5 expression levels of and that of low CCL5 expression. GSEA revealed significant differences in the enrichment of the KEGG Collection between CCL5 high expression group and low CCL5 expression groups.

Relied on the false discovery rate (FDR) and the normalized enrichment score (NES), top 20 signaling transduction pathways which were satisfied with p<0.01 were ranked by NES, as shown in Table 2. It contained chemokine signaling pathways. The top 5 pathways were cytokine-cytokine receptor interaction, autoimmune thyroid disease, antigen processing and antigen presentation, natural killer cell mediated cytotoxicity, and immune network for IgA production. Among these pathways, some pathways were related to immune response. So, a speculation that the expression levels of CCL5 may play an important role in immune response in patients with KIRC was raised. Figure 2 showed chemokine signaling pathway, cytokine-cytokine receptor interaction pathway and T cell receptor signaling pathway analyzed by GSEA.

Overview of immune infiltration in RCC and normal samples

CIBERSORT, which provides an estimation of the abundance of immune cells in a mixed cell populations, was utilized to assess the proportion of 22 immune cell types in different samples. The proportions difference of all 22 types of immune infiltration cells between tumor tissue and normal tissue was analyzed. The results illustrated that 5 types of immune cells were significantly different between the two groups (p<0.01) as shown in Fig 3. Among these cells, naive B cell, CD4+ memory resting T cells, resting Dendritic cells were significantly lower in tumor group compared with the normal group, While CD8+ T cells, Tregs were significantly higher in tumor group (p<0.01).

 Fig 3 

The proportions difference of all 22 types of immune infiltration cells between tumor tissue and normal tissue. Comparisons of 22 immune cells between tumor (red parts) and normal (blue parts) tissues showed that naive B cell, CD4+ memory resting T cells, resting Dendritic cells were lower in the tumor group, and CD8+ T cells, Tregs were higher in tumor group, compared with the normal group (p<0.01).

Int J Med Sci Image

(View in new window)

 Table 2 

Gene sets enriched in the samples of KIRC in the high expression of CCL5 group.

Gene set nameSIZEESNESNOM p-valFDR q-val
KEGG_CYTOKINE_CYTOKINE_RECEPTOR_INTERACTION2640.640152.722257900
KEGG_AUTOIMMUNE_THYROID_DISEASE500.907812.713472100
KEGG_ANTIGEN_PROCESSING_AND_PRESENTATION810.821372.643913500
KEGG_NATURAL_KILLER_CELL_MEDIATED_CYTOTOXICITY1320.697612.630446200
KEGG_INTESTINAL_IMMUNE_NETWORK_FOR_IGA_PRODUCTION460.886172.620240200
KEGG_CELL_ADHESION_MOLECULES_CAMS1310.699552.574541600
KEGG_VIRAL_MYOCARDITIS680.774042.480061506.15E-05
KEGG_HEMATOPOIETIC_CELL_LINEAGE850.690172.442231201.08E-04
KEGG_TYPE_I_DIABETES_MELLITUS410.92422.416135502.39E-04
KEGG_PRIMARY_IMMUNODEFICIENCY350.876742.368313803.42E-04
KEGG_CYTOSOLIC_DNA_SENSING_PATHWAY550.643242.367143603.11E-04
KEGG_ASTHMA280.868542.340905203.86E-04
KEGG_CHEMOKINE_SIGNALING_PATHWAY1880.598632.2787430.0019342366.73E-04
KEGG_ALLOGRAFT_REJECTION350.9452.259700506.87E-04
KEGG_LEISHMANIA_INFECTION700.756272.24858120.0019920326.64E-04
KEGG_GRAFT_VERSUS_HOST_DISEASE370.929162.239085706.98E-04
KEGG_T_CELL_RECEPTOR_SIGNALING_PATHWAY1080.62432.1902990.0019417480.00127215
KEGG_SYSTEMIC_LUPUS_ERYTHEMATOSUS1350.720482.184673300.001251814
KEGG_JAK_STAT_SIGNALING_PATHWAY1550.486352.08528850.0019120460.003184194

Overall survival time of immune infiltration cells in KIRC

The relationships between the 22 types of immune cells and the overall survival time of patients were analyzed by survival analysis (K-M). It was found that three types of immune cells had a significant impact on the overall survival time. As showed in Fig 4, T cells regulatory, mast resting cells and dendritic resting cells had significant effects on the overall survival time (p<0.01). Regulatory T cells had a negative regulation on the survival time of patients with KIRC which is consistent with their function, but resting mast cells and resting dendritic cells had positive effects on it.

Relationship between CCL5 expression and immune infiltration cells

Regulatory T cells, resting Mast cells and resting dendritic cells had shown significant effects on the OS. Then, non-parametric tests including Kruskal-Walls test, Mann-Whitney test (U-test) and independent sample median test were utilized to explore the relationships between the 3 types of immune infiltration cells and the expression levels of CCL5. The group of high expression levels of CCL5 had positive effects on the proportions of regulatory T cells (median 0.183 vs 0.064, p<0.01). The significant correlation between the expression levels of CCL5 and the proportions of resting mast cells were found (median 0.308 vs 0.049, p<0.01) shown in Fig 5. It was found that the proportion of resting dendritic cells was not influenced by the expression levels of CCL5.

Then, “correlation” module (Spearman test) of GEPIA was used to confirm the result. The gene biomarkers of mast cells, including CPA3, TPSB2, TPSAB1, MS4A2, HDC, were not correlated with CCL5. But the gene biomarkers of Tregs, including CD2, CD25, Foxp3, CTLA-4, were correlated with CCL5 significantly (p<0.01), shown in Table 3. So, CCL5 may play an important role in regulating the abundance of Tregs in KIRC.

 Table 3 

The results of the correlation between CCL5 and gene markers of mast cells and Tregs by the “correlation” module of GEPIA.

Gene markersPR
Mast cellsCPA30.720.016
TPSB20.110.07
TPSAB10.0850.076
MS4A20.580.024
HDC0.570.025
TregsCD4<0.010.88
CD25<0.010.27
Foxp3<0.010.62
CTLA-4<0.010.64

Discussion

In this study, the correlation between the expression levels of CCL5 in KIRC and some clinical features including the survival time, pathological stage, and so on, was analyzed to explore the value of CCL5 as a tumor biomarker in KIRC. Meanwhile, KEGG pathways analysis for CCL5 showed an interesting relationship between the expression of CCL5 and some immune signaling transduction pathways, such as T cells receptor pathway. CIBERSORT was utilized to explore the relationship between the difference of ratios of 22 types of immune infiltration cells and the different expression levels of CCL5 in KIRC. The relationships were verified by the “correlation” module of GEPIA. The results showed that the expression of CCL5 is related to regulatory T cells, which led to poor survival time in patients with KIRC.

 Fig 4 

Relationships between the Overall survival time and three types of immune infiltration cells. A, regulatory T cells. B, resting mast cells. C, resting dendritic cell

Int J Med Sci Image

(View in new window)

 Fig 5 

Relationship between CCL5 expression and immune infiltration cells. Mann-Whitney test was utilized, and the low and high expression groups of CCL5 were grouped by CCL5 median expression (30.5). A, regulatory T cells. B, resting mast cells. C, resting dendritic cell.

Int J Med Sci Image

(View in new window)

Kidney Renal Clear Cell Carcinoma (KIRC) is the most representative subtype of kidney cancers. However, the clinical treatment effect is poor, and the 5-year survival of patients with the advanced stage is as low as 12%[6]. So, identifying more biomarkers for diagnosis and prognosis is important and urgent. CCL5 is widely secreted by natural killer cells, T cells, fibroblasts, epithelial cells, platelets, and some certain types of tumor cells [31, 32]. It has been demonstrated as a biomarker for the prognosis in esophageal cancer, colorectal cancer, breast cancer and so on [33-35]. But in kidney renal clear cell carcinoma, there was only one study which implied that CCL5 may be a crucial potential gene associated with carcinogenesis of KIRC based on the literature review [36]. This present study found that CCL5 mRNA was significantly overexpressed in KIRC tissues compared with normal tissues adjacent to cancer. This study also revealed significant connection between the expression levels of CCL5 and disease stage, pathological grade, and pathological TNM stage. In addition, survival analysis also revealed that the high expression level of CCL5 induced a poor survival time compared with the low expression group. These results indicate that CCL5 gene plays an important role in the development of KIRC.

Multiple studies have confirmed that chemokines, which contain 4 subgroups CXC, CC, CX3C and C, play important roles in regulating tumor microenvironment. These chemokines including CCL5 can bind to seven-transmembrane G protein receptors on the target cell membrane, thereby activating downstream effector molecules and promoting the malignant process of tumors. CCL5 promotes tumor cell proliferation through the mTOR pathway or increases glucose uptake. It can also indirectly promote tumor cell proliferation by recruiting macrophages, fibroblasts, T cells, etc [37, 38]. At the same time, CCL5 can promote angiogenesis through promoting endothelial cell migration, neovascularization, and the secretion of vascular endothelial growth factor (VEGF) [39-41]. In addition, it plays an important role in tumor cell migration through upregulating MMP2 and MMP9 expression or activating integrin [42-44].

In this study, KEGG pathways selected by NES scores from GSEA results contained chemokine signaling pathway and cytokine-cytokine receptor interaction, which may also verify the functions of CCL5. Meanwhile, many signaling transduction pathways related to immune functions, such as T cell receptor pathways, were contained in these KEGG pathways. Studies had revealed that CCL5 promoted tumor migration and invasion through recruiting M2 type macrophages (TAMs) in breast cancer [32]. Meanwhile, Tregs had been found that it was recruited by CCL5-CCR5 axis in ovarian cancer [45]. CCL5 may have effects on tumor infiltration cells. So, CIBERSORT was utilized to explore tumor infiltration immune cells in KIRC. Tumor infiltration immune cells are a major component of the tumor microenvironment and are closely related to tumor cell proliferation, treatment response and prognosis [46, 47]. Different proportions of Regulatory T cells (Tregs) including naive B cell, resting CD4+ T cells memory, resting dendritic cells, CD8+ T cells were identified between tumor tissues and normal tissue. And among these five types of tumor infiltration immune cells, Tregs was the only one, which not only significantly had effects on survival time but also was significantly related to the expression levels of CCL5 in KIRC.

Tregs are a group of lymphocytes which negatively regulate the body's immune response and play an important role in maintaining self-tolerance and avoiding the immune response from self-damage. It is simply divided into two types by different sources, one is natural Treg (nTreg) which is derived from thymus and mostly expressed FOXP3, another is inducible Treg (iTreg) which is induced peripherally[48]. Tregs have been studied because of its function of suppressing the bodies' anti-tumor immune response [49]. At present, it is thought that the immunosuppressive effect is exerted mainly through the following ways. Firstly, Tregs are activated by T cell receptor (TCR) to directly inhibit the formation of effector T cells. Secondly, Treg inhibit the activity of immune cells through cross-linking between cell surface cytotoxic T lymphocyte-associated antigen-4 (cTLA-4) and antigen-presenting cells (APC). Thirdly, Tregs secrete cytokines such as IL-10 and TGF-B which can inhibit the function of effector T cells [50, 51]. Some research reported that an increased number of Tregs no matter in peripheral blood or in tumors have a negative correlation with survival time in KIRC [52-54]. Some studies utilized various methods to deplete Treg cells or control Treg cell function and these methods were effective in treating tumors by increasing immune ability of anti-tumor[55, 56]. But the lack of special targets for depletion and functional impairment of Treg cells is a vast obstacle to overcome[57].

CCL5 acts as an important factor in recruiting Tregs and assisting Tregs to implement its function suggested by some studies. A research had demonstrated that tumor-infiltrating Tregs were increased by doing intratumoral injection of CCL4 or CCL5 in mice models. It also shown an interesting result that tumor-infiltrating MO-MDSCs (myeloid-derived suppressor cells) produce high levels of the CCL3, CCL4, and CCL5 which is CCR5 ligands, and directly attract Tregs in a CCR5-dependent manner [58]. Meanwhile, a study identified that zoledronic acid can significantly affect the interaction between Tregs and cancer cells by reducing the expression of CCL5 and CCL2[59]. Furthermore, Tregs, which were recruited by CCL5, promoted the invasion of ovarian cancer cells through matrix metalloproteinase-9 (MMP9), which enhanced the degradation of the extracellular matrix and enabled the invasion of tumor cells, also been identified [45, 60]. So, this present study demonstrated that a higher expression level of CCL5 is associated with the more abundance of Tregs in KIRC. CCL5 may act as a typical chemokine to recruit Tregs to mediated immunosuppression and lead to immune escape in KIRC. Furthermore, down regulation of CCL5 in KIRC may improve the immune effect by decreasing recruitment of Tregs. CCL5 may act as a target of immune therapy in the future.

In summary, the increased expression of CCL5 is related to poor prognosis and some clinical features including TNM stage, pathological grade. So CCL5 could act as a biomarker for the prognosis and potential therapy target for patients with KIRC. Meanwhile, the poor prognosis may be caused by the function of CCL5 to recruit more Tregs which induced the immune suppression in patients with KIRC.

All these results were based on a series of bioinformatics algorithms and databases. Experiments, as immunohistochemistry, are needed to confirm these results.

Availability of data and materials

The datasets used and analyzed during the current study are available from the TCGA database.

Funding

This manuscript is supported by the National Natural Science Foundations of China (Juan Ren, 81772793/H1621, Juan Ren, 31201060/C0709; Juan Ren, 30973175/H1621; and Juan Ren, 81772793/H1621); Supported by Scientific and Technological Research Foundation of Shaan'xi Province (Juan Ren, 2020JM-368); Supported by Program for New Century Excellent Talents in University (Juan Ren, NCET-12-0440); The Fundamental Research Funds for the Central Universities (Juan Ren, 2012).

Author Contributions

JR conceived and supervised the study; JR, SHB, YYW, YYL, HJK, JZZ, WM, YG, RL, and BNH analyzed data; JR and SHB wrote the manuscript; JR, SHB and XZZ made manuscript revisions. All authors have read and approved the final version of this submission.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7-34

2. Siska PJ, Beckermann KE, Rathmell WK, Haake SM. Strategies to overcome therapeutic resistance in renal cell carcinoma. Urol Oncol. 2017;35:102-10

3. Hsieh JJ, Le VH, Oyama T, Ricketts CJ, Ho TH, Cheng EH. Chromosome 3p Loss-Orchestrated VHL, HIF, and Epigenetic Deregulation in Clear Cell Renal Cell Carcinoma. J Clin Oncol. 2018: JCO2018792549.

4. Li X, Meng X, Wei C, Zhou Y, Chen H, Huang H. et al. Dissecting LncRNA Roles in Renal Cell Carcinoma Metastasis and Characterizing Genomic Heterogeneity by Single-Cell RNA-seq. Mol Cancer Res. 2018;16:1879-88

5. Bedke J, Gauler T, Grunwald V, Hegele A, Herrmann E, Hinz S. et al. Systemic therapy in metastatic renal cell carcinoma. World J Urol. 2017;35:179-88

6. Atkins MB, Tannir NM. Current and emerging therapies for first-line treatment of metastatic clear cell renal cell carcinoma. Cancer Treat Rev. 2018;70:127-37

7. Choueiri TK, Fishman MN, Escudier B, McDermott DF, Drake CG, Kluger H. et al. Immunomodulatory Activity of Nivolumab in Metastatic Renal Cell Carcinoma. Clin Cancer Res. 2016;22:5461-71

8. Flippot R, Escudier B, Albiges L. Immune Checkpoint Inhibitors: Toward New Paradigms in Renal Cell Carcinoma. Drugs. 2018;78:1443-57

9. Li Y, Yang S, Yue H, Yuan D, Li L, Zhao J. et al. Unraveling LGALS1 as a Potential Immune Checkpoint and a Predictor of the Response to Anti-PD1 Therapy in Clear Cell Renal Carcinoma. Pathol Oncol Res. 2019

10. Sharpe AH, Pauken KE. The diverse functions of the PD1 inhibitory pathway. Nat Rev Immunol. 2018;18:153-67

11. Leighton SP, Nerurkar L, Krishnadas R, Johnman C, Graham GJ, Cavanagh J. Chemokines in depression in health and in inflammatory illness: a systematic review and meta-analysis. Mol Psychiatr. 2018;23:48-58

12. Singh R, Singh S, Mishra M, Lillard J. Biological significance of CCR5/CCL5 axis in liver cancer. Eur J Immunol. 2019;49:148 -

13. Walens A, DiMarco AV, Lupo R, Kroger BR, Damrauer JS, Alvarez JV. CCL5 promotes breast cancer recurrence through macrophage recruitment in residual tumors. Elife. 2019 8

14. Suenaga M, Cao S, Zhang W, Ning Y, Okazaki S, Berger MD. et al. Genetic variants in CCL5 and CCR5 genes and serum VEGF-A levels to predict efficacy of bevacizumab in metastatic colorectal cancer patients receiving first-line chemotherapy. J Clin Oncol. 2017 35

15. Wu YC, Shen YC, Chang JWC, Hsieh JJ, Chu Y, Wang CH. Autocrine CCL5 promotes tumor progression in esophageal squamous cell carcinoma in vitro. Cytokine. 2018;110:94-103

16. Wang W, Ren S, Wang Z, Zhang C, Huang J. Increased expression of TTC21A in lung adenocarcinoma infers favorable prognosis and high immune infiltrating level. Int Immunopharmacol. 2019;78:106077

17. Gooden MJM, de Bock GH, Leffers N, Daemen T, Nijman HW. The prognostic influence of tumour-infiltrating lymphocytes in cancer: a systematic review with meta-analysis. Brit J Cancer. 2011;105:93-103

18. Xiong Y, Liu L, Xia Y, Qi Y, Chen Y, Chen L. et al. Tumor infiltrating mast cells determine oncogenic HIF-2alpha-conferred immune evasion in clear cell renal cell carcinoma. Cancer Immunol Immunother. 2019;68:731-41

19. Siska PJ, Beckermann KE, Mason FM, Andrejeva G, Greenplate AR, Sendor AB. et al. Mitochondrial dysregulation and glycolytic insufficiency functionally impair CD8 T cells infiltrating human renal cell carcinoma. JCI Insight. 2017 2

20. Fu H, Zhu Y, Wang Y, Liu Z, Zhang J, Wang Z. et al. Tumor Infiltrating Mast Cells (TIMs) Confers a Marked Survival Advantage in Nonmetastatic Clear-Cell Renal Cell Carcinoma. Ann Surg Oncol. 2017;24:1435-42

21. Chen YH, Chen SH, Hou J, Ke ZB, Wu YP, Lin TT. et al. Identifying hub genes of clear cell renal cell carcinoma associated with the proportion of regulatory T cells by weighted gene co-expression network analysis. Aging (Albany NY). 2019;11:9478-91

22. Lapteva N, Huang XF. CCL5 as an adjuvant for cancer immunotherapy. Expert Opin Biol Th. 2010;10:725-33

23. Czaja AJ. Review article: chemokines as orchestrators of autoimmune hepatitis and potential therapeutic targets. Aliment Pharm Ther. 2014;40:261-79

24. Araujo JM, Gomez AC, Aguilar A, Salgado R, Balko JM, Bravo L. et al. Effect of CCL5 expression in the recruitment of immune cells in triple negative breast cancer. Sci Rep-Uk. 2018 8

25. Chang LY, Lin YC, Mahalingam J, Huang CT, Chen TW, Kang CW. et al. Tumor-derived chemokine CCL5 enhances TGF-beta-mediated killing of CD8(+) T cells in colon cancer by T-regulatory cells. Cancer Res. 2012;72:1092-102

26. Liu J, Lichtenberg T, Hoadley KA, Poisson LM, Lazar AJ, Cherniack AD. et al. An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell. 2018;173:400-16 e11

27. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545-50

28. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y. et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12:453-7

29. Zhang S, Zhang E, Long J, Hu Z, Peng J, Liu L. et al. Immune infiltration in renal cell carcinoma. Cancer Sci. 2019;110:1564-72

30. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:W98-W102

31. Donlon TA, Krensky AM, Wallace MR, Collins FS, Lovett M, Clayberger C. Localization of a Human T-Cell-Specific Gene, Rantes (D17s136e), to Chromosome 17q11.2-Q12. Genomics. 1990;6:548-53

32. An G, Wu F, Huang S, Feng L, Bai J, Gu S. et al. Effects of CCL5 on the biological behavior of breast cancer and the mechanisms of its interaction with tumorassociated macrophages. Oncol Rep. 2019;42:2499-511

33. Goto M, Liu M. Chemokines and their receptors as biomarkers in esophageal cancer. Esophagus. 2019

34. Araujo JM, Gomez AC, Aguilar A, Salgado R, Balko JM, Bravo L. et al. Effect of CCL5 expression in the recruitment of immune cells in triple negative breast cancer. Sci Rep. 2018;8:4899

35. Ucuncu M, Serilmez M, Sari M, Bademler S, Karabulut S. The Diagnostic Significance of PDGF, EphA7, CCR5, and CCL5 Levels in Colorectal Cancer. Biomolecules. 2019 9

36. Song E, Song W, Ren M, Xing L, Ni W, Li Y. et al. Identification of potential crucial genes associated with carcinogenesis of clear cell renal cell carcinoma. J Cell Biochem. 2018;119:5163-74

37. Aldinucci D, Colombatti A. The Inflammatory Chemokine CCL5 and Cancer Progression. Mediat Inflamm. 2014

38. Aldinucci D, Casagrande N. Inhibition of the CCL5/CCR5 Axis against the Progression of Gastric Cancer. International Journal of Molecular Sciences. 2018 19

39. Sax MJ, Gasch C, Athota VR, Freeman R, Rasighaemi P, Westcott DE. et al. Cancer cell CCL5 mediates bone marrow independent angiogenesis in breast cancer. Oncotarget. 2016;7:85437-49

40. Wang SW, Liu SC, Sun HL, Huang TY, Chan CH, Yang CY. et al. CCL5/CCR5 axis induces vascular endothelial growth factor-mediated tumor angiogenesis in human osteosarcoma microenvironment. Carcinogenesis. 2015;36:104-14

41. Wang SW, Liu SC, Tang CH. CCL5 promotes vascular endothelial growth factor-dependent angiogenesis in human osteosarcoma. Ann Oncol. 2015;26:129 -

42. Huang CY, Fong YC, Lee CY, Chen MY, Tsai HC, Hsu HC. et al. CCL5 increases lung cancer migration via PI3K Akt and NF-kappa B pathways. Biochem Pharmacol. 2009;77:794-803

43. Long HX, Xie RK, Xiang T, Zhao ZQ, Lin S, Liang ZQ. et al. Autocrine CCL5 Signaling Promotes Invasion and Migration of CD133+Ovarian Cancer Stem-Like Cells via NF-?B-Mediated MMP-9 Upregulation. Stem Cells. 2012;30:2309-19

44. Kato T, Fujita Y, Nakane K, Mizutani K, Terazawa R, Ehara H. et al. CCR1/CCL5 interaction promotes invasion of taxane-resistant PC3 prostate cancer cells by increasing secretion of MMPs 2/9 and by activating ERK and Rac signaling. Cytokine. 2013;64:251-7

45. You Y, Li Y, Li M, Lei M, Wu M, Qu Y. et al. Ovarian cancer stem cells promote tumour immune privilege and invasion via CCL5 and regulatory T cells. Clin Exp Immunol. 2018;191:60-73

46. Fridman WH, Pages F, Sautes-Fridman C, Galon J. The immune contexture in human tumours: impact on clinical outcome. Nat Rev Cancer. 2012;12:298-306

47. Finotello F, Trajanoski Z. Quantifying tumor-infiltrating immune cells from transcriptomics data. Cancer Immunol Immun. 2018;67:1031-40

48. Kohno M, Murakami J, Wu L, Chan M, Yun Z, De Perrot M. Selective Depletion of Foxp3(+) Regulatory T Cells After Local Radiotherapy Induces Abscopal Effects in Murine Malignant Mesothelioma. J Thorac Oncol. 2019;14:S477-S8

49. Wing JB, Tanaka A, Sakaguchi S. Human FOXP3(+) Regulatory T Cell Heterogeneity and Function in Autoimmunity and Cancer. Immunity. 2019;50:302-16

50. Rothstein DM, Camirand G. New insights into the mechanisms of Treg function. Curr Opin Organ Tran. 2015;20:376-84

51. Tanaka A, Sakaguchi S. Regulatory T cells in cancer immunotherapy. Cell Research. 2017;27:109-18

52. Griffiths RW, Elkord E, Gilham DE, Ramani V, Clarke N, Stern PL. et al. Frequency of regulatory T cells in renal cell carcinoma patients and investigation of correlation with survival. Cancer Immunol Immun. 2007;56:1743-53

53. Finke JH, Rini B, Ireland J, Rayman P, Richmond A, Golshayan A. et al. Sunitinib Reverses Type-1 Immune Suppression and Decreases T-Regulatory Cells in Renal Cell Carcinoma Patients. Clin Cancer Res. 2008;14:6674-82

54. Zhu G, Pei L, Yin H, Lin F, Li X, Zhu X. et al. Profiles of tumor-infiltrating immune cells in renal cell carcinoma and their clinical implications. Oncol Lett. 2019;18:5235-42

55. Tanaka A, Sakaguchi S. Targeting Treg cells in cancer immunotherapy. Eur J Immunol. 2019;49:1140-6

56. Tanaka A, Sakaguchi S. Regulatory T cells in cancer immunotherapy. Cell Res. 2017;27:109-18

57. Ohue Y, Nishikawa H. Regulatory T (Treg) cells in cancer: Can Treg cells be a new therapeutic target?. Cancer Sci. 2019;110:2080-9

58. Schlecker E, Stojanovic A, Eisen C, Quack C, Falk CS, Umansky V. et al. Tumor-Infiltrating Monocytic Myeloid-Derived Suppressor Cells Mediate CCR5-Dependent Recruitment of Regulatory T Cells Favoring Tumor Growth. J Immunol. 2012;189:5602-11

59. Liu H, Wang SH, Chen SC, Chen CY, Lin TM. Zoledronic acid blocks the interaction between breast cancer cells and regulatory T-cells. BMC Cancer. 2019;19:176

60. Qiu H, Zmina PM, Huang AY, Askew D, Bedogni B. Inhibiting Notch1 enhances immunotherapy efficacy in melanoma by preventing Notch1 dependent immune suppressive properties. Cancer Lett. 2018;434:144-51

Author contact

Corresponding address Corresponding author: Juan Ren, Ph.D., P.R. China, Tel: +86 29 85324029, E-mail address: 869491533com (J. Ren)


Received 2020-7-26
Accepted 2020-9-27
Published 2020-10-18