Skip to main content

CTLA4 has a profound impact on the landscape of tumor-infiltrating lymphocytes with a high prognosis value in clear cell renal cell carcinoma (ccRCC)



Cytotoxic T-lymphocyte associated protein 4 (CTLA4) inhibitors have been shown to significantly prolong the overall survival (OS) in a wide range of cancers. However, its application in clear cell renal cell carcinoma (ccRCC) is limited due to the therapy response, and the prognostic value of CTLA4 in ccRCC has not been investigated in detail.


By using immunohistochemistry, Kaplan–Meier (K–M) analysis, uni- and multi-variate Cox analysis, we comprehensively and systematically studied the prognostic value of CTLA4 in ccRCC. Then, we applied Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) and CIBERSORT, ESTIMATE algorithm, ssGSEA and somatic mutation analyses to reveal the impact of CTLA4 on the landscape of tumor-infiltrating lymphocytes (TILs) infiltration and genetic mutation. Besides, given current concerns caused by combined immunotherapy, we also investigated the relationship between CTLA4 and other immune checkpoints.


In vitro experiment and data mining showed that, CTLA4 was up-regulated in ccRCC tissues and closely related to the disease progression as well as a poor prognosis. Deeper researches demonstrated that CTLA4 regulates T cell activation and was significantly linked to TIL-abundant tumor microenvironment (TME), but was accompanied by an immunosuppressed phenotype. Mutation analysis showed that CTLA4 was associated with more frequent BRCA-associated protein 1 (BAP1) mutation. Moreover, we found that CTLA4 was markedly correlated with multiple immune checkpoints, which suggested that ccRCC patients with high expressed CTLA4 may benefit more from immune checkpoint blockades (ICBs) combined therapy.


CTLA4 has a profound impact on the landscape of TILs and genetic mutation, and can be used as the biomarker with high prognosis value in ccRCC.


Clear cell renal cell carcinoma (ccRCC) is the most common and fatal histological subtype of renal cell carcinoma in adults, accounting for about 65–70% in RCC, characterized by abundant tumor-infiltrating lymphocytes (TILs) infiltration within the tumor microenvironment (TME) [1, 2]. Since the symptoms were not obvious, 20% of patients were initially diagnosed with metastases and nearly 30% relapsed with metastasis after surgical excision [3]. Therefore, the development of novel therapeutic strategies for ccRCC is necessary.

Along with advances in cancer immunology, the role of TME in ccRCC has attracted increased attention in recent years, consisting of cancer cells, fibroblasts, myofibroblasts, endothelial cells, TILs, and extracellular matrix [4]. The composition of TILs within TME determines whether its phenotype is anti-tumor immunity or immune evasion. Anti-tumor immunity is characterized by CD8+ T cells, M1 macrophages, while immune evasion is characterized by mast cells, T cells regulatory (Tregs), and M2-macrophages. Macrophages in tumor tissues are biased towards M2 subtype [5]. Furthermore, immune checkpoint within the TILs is a crucial factor in maintaining immune evasion of TME [6]. E.g. Tregs can inhibit the activation of CD8+ T cells through Cytotoxic T-lymphocyte associated protein 4 (CTLA4), triggering tumor immunosuppression [7]. Numerous clinical studies have reported that the immune checkpoint would be an ideal target for driving T cell mediated anti-tumor immunity [8, 9]. With tremendous progress, immune checkpoint blockades (ICBs) such as programmed death-ligand 1 (PD-L1) inhibitor Durvalumab and CTLA4 inhibitor Tremelimumab et al. have been shown to significantly prolong overall survival (OS) in a wide range of cancers [10,11,12].

CTLA4, immune checkpoint protein, has received extensive attention in immunotherapy. CTLA4 was highly expressed in CD8+ T cells and inhibited the T cells activation through competitively blocking the binding of CD28 with B7, leading to immune evasion [13]. Since the CTLA4 inhibitor ipilimumab (the first ICB approved by the FDA) has significantly prolonged the OS of patients with metastatic melanoma, CTLA4 inhibitors have proven to be effective agents for many cancers [14]. However, the clinical application of CTLA4 inhibitors in ccRCC is strictly limited, for that the therapeutic response and prognostic value of CTLA4 in ccRCC have not been investigated in detail. In the current work, we attempted to comprehensively analyze the prognostic value of CTLA4 in ccRCC and its impact on TILs and genetic landscape through bioinformatics and in vitro experiment, which may be beneficial to the therapeutic response of CTLA4 inhibitors in ccRCC patients.

Materials and methods


The mRNA-seq data from 533 ccRCC and 72 normal tissues and corresponding clinical information, as well as 336 ccRCC somatic mutation data, were retrieved from The Cancer Genome Atlas (TCGA,, patients with OS more than 30 days were retained. Ultimately, 494 ccRCC and 68 normal tissues were included for our analyses. RNA-array data GSE40435 (GPL10558) and GSE46699 (GPL570) were obtained from the Gene Expression Omnibus (GEO,, including 101 ccRCC and paired normal kidney tissues, 67 ccRCC and 63 normal tissues, respectively.


Formalin Fixed Paraffin Embedded (FFPE) ccRCC (n = 5) and corresponding normal specimens (n = 5) were collected from the department of pathology, Renmin Hospital of Wuhan University (China). Our study was approved by the Ethics Committee of Renmin Hospital of Wuhan University.

Immunohistochemistry with the rabbit monoclonal to CTLA4 (abcam, Shanghai, China) and the secondary antibody (Aspen, Wuhan, China) was performed following the manufacturer’s instructions [15]. Images were obtained using a BX63 microscope (Olympus, Japan), and the density of CTLA4 expression was evaluated by the percentage of CTLA4 positive cells in the total cells.

Analysis of the functions of CTLA4

First, the CTLA4 related genes were screened based on Spearman correlation method with the absolute value of the correlation coefficient > 0.6 and p value < 0.05 as the threshold. Then, Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out through R package clusterProfiler [16] to reveal the potential function of CTLA4 in ccRCC.

Genetic mutation analysis

The somatic mutation data of ccRCC patients from TCGA were analyzed by R package Maftool [17] to identify the impact of genetic mutation by CTLA4.

Immune landscape of TME

CIBERSORT was a gene-based deconvolution algorithm developed by Newman et al. [18] and was applied to predict the abundance of immune cells using complex gene expression data in this investigation. Furthermore, the immune score was calculated by ESTIMATE algorithm, and the immunosuppression score was obtained from ssGSEA [19, 20].


Data analysis was conducted by using R 3.6.2 together with Graphpad prism 6. Difference of continuous variables between the two groups was calculated by the Wilcoxon test or student’s t test according to Shapiro-test and Bartlett-test. Chi-square (χ2) test was applied to assess the difference between categorical variables. Kaplan–Meier (K–M) analysis (Log-rank test based), uni- and multi-variate Cox analyses were used to evaluate the survival prognosis. Spearman method was utilized in correlation analysis. Data are presented as mean ± SD and P value < 0.05 was considered statistically significant (p < 0.05*, p < 0.01**, p < 0.001***).


High prognosis value of CTLA4 in ccRCC

In TCGA dataset, CTLA4 was highly expressed in tumor tissues (N = 494) compared to normal tissues (N = 68, p < 0.001), and the paired test also confirmed that CTLA4 was highly expressed in ccRCC (N = 68, p < 0.001) (Fig. 1a, b). Moreover, CTLA4 overexpression in ccRCC had also been reproduced in the GEO datasets GSE40435 and GSE46699 (Fig. 1c).

Fig. 1

CTLA4 has a high prognosis value in ccRCC. The expression of CTLA4 between ccRCC and normal renal tissues in a TCGA dataset and c GSE46699. The expression of CTLA4 between paired ccRCC and normal tissues in b TCGA dataset and c GSE40435. d Verifying the expression of CTLA4 in ccRCC tissues through immunohistochemistry. Scale bars, 100 mm. Survival analysis in ccRCC patients from e TCGA dataset and f GEPIA database. ccRCC

The immunohistochemistry was utilized to further validate the above results, and the result showed that the density of CTLA4 expression was higher in ccRCC tissues compared with normal tissues (p < 0.001, Fig. 1d). It seems that CTLA4 was steadily up-regulated in ccRCC both in data mining and in vitro experiment.

In TCGA dataset, patients were initially classified into high and low groups based on CTLA4 level. The K-M curves showed that CTLA4 shortened the OS in ccRCC, with the median survival time being 63.73 and 91.73 months in high and low group, respectively (p < 0.001) (Fig. 1e). In addition, 512 cases from GEPIA also confirmed that CTLA4 was a risk gene with Hazard ratio (HR) > 1 (HR = 1.5, p = 0.013, Fig. 1f).

We next sought to investigate the role of CTLA4 in cancer progression, and the result suggested that the overexpressed CTLA4 was related to high grade (χ2 = 12.465, p < 0.001), advanced stage (χ2 = 22.510, p < 0.001), patient with tumor state (χ2 = 7.874, p = 0.005) and death (χ2 = 9.965, p = 0.002, Table 1). We also found that CTLA4 had a positive correlation with high grade in GSE40435 cohorts (χ2 = 3.971, p = 0.046), indicating that CTLA4 may function as an oncogene in the progress of ccRCC. Subsequently, uni- and multi-variate Cox analyses found that some variables, including CTLA4, age, tumor pathological grade and stage, were independent risk factors in ccRCC (p < 0.01) (Table 2). These results proved that CTLA4 contributed to the progression of ccRCC with a high prognosis value.

Table 1 Correlation between CTLA4 and clinicopathological features in ccRCC
Table 2 Univariate and multivariate Cox analyses

CTLA4 indicated a higher density of TILs in ccRCC tumor microenvironment, but an immunosuppressed phenotype.

To outline the corresponding function of the CTLA4 in ccRCC TME, we performed KEGG and GO analysis based on 200 CTLA4 related protein coding genes (Additional file 1: Table S1). 35 KEGG items were identified, including T cell mediated immune related pathway, natural killer cell mediated cytotoxicity, T helper cells differentiation, PD-L1 expression and programmed death‐1 receptor (PD-1) checkpoint pathway in cancer (Table 3), and CTLA4 was positively correlated with T cell receptor signaling pathway (cor = 0.80), Natural killer cell mediated cytotoxicity (cor = 0.75), Th1 and Th2 cell differentiation (cor = 0.75) and Th17 cell differentiation (cor = 0.78) (Fig. 2a). The biological processes of CTLA4 yielded from the GO analysis were associated with the activation and differentiation of T cells (Table 3) and positively correlated with T cell activation (cor = 0.79), regulation of lymphocyte activation (cor = 0.79) and T cell differentiation(cor = 0.80) (Fig. 2b).

Table 3 Top 10 items in functional enrichment analysis
Fig. 2

CTLA4 indicated a higher density of TILs in ccRCC tumor microenvironment, but an immunosuppressed phenotype. a CTLA4 is positively correlated with T cell receptor signaling pathway, Natural killer cell mediated cytotoxicity, Th1 and Th2 cell differentiation and Th17 cell differentiation. b CTLA4 is positively correlated with T cell activation and T cell differentiation. c Highly expressed CTLA4 has a higher Immune score, indicating that there are more TILs within the TME. d CTLA4 is positively correlated with infiltrating CD8+ T cells and Tregs. e, f Highly expressed CTLA4 is accompanied by a higher Immunosuppression score and higher expression of T cell exhaustion markers, indicating an immunosuppressed phenotype. TILs, tumor-infiltrating lymphocytes; TME, tumor microenvironment; Tregs, T regulatory cells

Previous studies have demonstrated that TILs within the TME can be regarded as a prognostic indicator in ccRCC [21]. Besides, the results of GO and KEGG promoted us to continue to investigate the role of CTLA4 in TILs infiltration, which affected the ICBs’ response. Our results showed that CTLA4 was associated with a higher immune score, which was calculated by the ESTIMATE algorithm and represented the level of TILs, indicating that CTLA4 promoted the recruitment of immune cells into the TME (Fig. 2c). Furthermore, there was a great difference in the composition of TILs between high and low CTLA4 groups. CTLA4 increased the infiltration of T Cells CD8+ , Tregs, Macrophage M1, whereas Plasma cells, NK cells activated, Monocytes, Macrophage M2, Dendritic cells activated were less infiltrated in CTLA4 high group (Additional file 1: Table S2). The correlation analysis result was presented in Fig. 2d, showing that CTLA4 was positively correlated with CD8+ T cells (cor = 0.50, p < 0.001), Tregs (cor = 0.28, p < 0.001) (Fig. 2d). However, the immunosuppression score as well as the expression of CD8 + T cell exhaustion markers Hepatitis A virus cellular receptor 2 (HAVCR2), lymphocyte activation gene-3 (LAG3), and T cell immunoglobulin and ITIM domain (TIGIT) was higher in CTLA4 high group (Fig. 2e, f). All in all, CTLA4 changed the landscape of TILs in ccRCC TME, and indicated a higher density of TILs, especially the CD8+ T cells and Tregs, but faced an immunosuppressed phenotype.

Genetic altered by CTLA4 in ccRCC

Genetic changes include non-synonymous mutations, which are mainly composed of missense mutation, synonymous mutation, insertion or deletion, and copy number gain or loss [22,23,24,25]. Tumor mutation burden (TMB) can be used as a biomarker to predict the efficacy of ICBs [26]. Some studies have shown that the RCC was sensitive to ICBs, although the TMB in RCC was moderate [27, 28]. To identify the somatic mutations that were altered by CTLA4 in ccRCC, we performed the mutation analysis and the result showed that overexpressed CTLA4 was correlated with BRCA-associated protein 1 (BAP1) mutation (p < 0.05, Fig. 3). The TMB in the high CTLA4 expression group tended to be higher than the low expression group, although it was not statistically significant. Moreover, Nonsense Mutation and In Frame Ins in the high CTLA4 expression group were higher than those in the low group (Table 4). BAP1 is a deubiquitinating enzyme and considered to be a tumor suppressor, and the loss of BAP1 contributes to the metastasis and poor prognosis in various cancers [29].

Fig. 3

Genetic changes by CTLA4 in ccRCC. Genetic mutation analysis of CTLA4 between high and low CTLA4 groups, the result shows that the BAP1 mutation is more frequent in the CTLA4 high group

Table 4 Classification of variants in different groups

CTLA4 was highly related to other immune checkpoint molecules

Recently, the combined inhibition of PD-L1 and CTLA4 has attracted much attention [30]. Planchard et al. [31] reported that combination immunotherapy of PD-L1 and CTLA4 considerably prolonged the OS in advanced refractory colorectal cancer. Combination immunotherapy tends to replace monotherapy, for that the combinational usage of ICBs can produce higher synergistic anti-tumor efficiency and reduce side effects [32]. Therefore, we continued to explore the correlation between CTLA4 and other immune checkpoint molecules, including PDCD1 (PD-1), CD274 (PD-L1), LAG3, indoleamine-2,3-dioxygenase-1 (IDO1), and TIGIT [33, 34]. The results showed that CTLA4 was highly and positively related to PD-1, PD-L1, LAG3, IDO1, and TIGIT (Fig. 4).

Fig. 4

CTLA4 was highly related to other immune checkpoint molecules. CTLA4 is highly and positively related to other immune checkpoint molecules PDCD1(PD-1), CD274(PD-L1), LAG3, IDO1, TIGIT


CTLA4, as a transmembrane protein expressed in activated CD4+ T and CD8+ T cells, has received a lot of attention for its interaction with cancer. CTLA4 negatively regulates T cell activation by blocking the function of costimulatory signal and differentiation cluster CD28:B7 binding [35]. CTLA4 inhibitors reverse inhibitory immune signal and restore the anti-cancer response by blocking the interaction between CTLA4 and the ligand expressed by antigen presenting cells [32]. With the approval of CTLA4 inhibitor Ipilimumab for clinical applications, it has been used for metastatic melanoma after the first-line treatment [14]. Furthermore, Romano et al. [36] proved that Ipilimumab can exert a therapeutic effect by targeting Tregs in tumors.

Several studies have shown that CD8+ T cells are in a state of abnormal activation when they turn to exhausted phenotypes due to the elimination of tumor cells, which can not only up-regulate the expression of immunosuppressive cytokines, but also directly lead to immunosuppression [37]. Exhausted CD8+ T cells continue to activate the expression of CTLA4 and other immune checkpoint receptors under the chronic stimulation of tumor antigens, which further promoting tumor invasion [38]. Here, we confirmed that CTLA4 was up-regulated in ccRCC tissues through data mining and in vitro experiment, and revealed that CTLA4 was correlated with poor prognosis. Corresponding to other studies, this study confirmed that CTLA4 played an important role in the regulation of T cells and represented more TILs infiltration to the TME, especially the CD8+ T cells and Tregs. However, the phenotype of the TME trended to immunosuppression, and the infiltrating CD8+ T cells biased to exhaustion in the CTLA4 high group, synergizing with Tregs, ultimately leading to tumor metastasis and progression.

We noticed that tumors with high TMB were sensitive to ICBs, contributing to a better outcome [26]. Therefore, we tried to outline the relationship between somatic mutation and CTLA4. We found that the Nonsense Mutation and In Frame Ins were markedly higher, and the BAP1 mutation was more frequent in the CTLA4 high group. As a tumor suppressor gene, the loss of BAP1 tended to cause poor prognosis and higher TMB, suggesting that patients with high CTLA4 expression might be more sensitive to ICBs. Furthermore, the somatic mutation of BAP1 is more abundant in highly metastatic tumors, such as uveal melanoma [39]. BAP1 is located at chromosome 3p21, adjacent to 3p25 where VHL is located [40]. Considering that 3p deletions in ccRCC are common, we think that 3p deletions might cause the inactivation of these tumor suppressor genes. Nevertheless, further explorations are needed to reveal the underlying mechanisms between CTLA4 and BAP1.

Finally, we revealed that CTLA4 was highly related to other immune checkpoints: PD-1, PD-L1, LAG3, IDO1, and TIGIT. Activation of the PD-1/PD-L1 signaling pathway contributes to TME with immune evasion, and its inhibitors are representative, which have been used in lots of solid tumors [41]. LAG3 can negatively regulate the activation and function of T cells, and its antagonists have been applied clinically [42]. At present, some studies have been devoted to the synergism between LAG3 and PD-1 in enhancing the efficacy of immunotherapy [43]. IDO1 is overexpressed in cancer cells, inhibiting the function of effector T cells and promoting the infiltration of Tregs. Studies have demonstrated that IDO1 is a promising target for improving patient outcomes in the field of immune- oncology [44]. The above results suggested that the CTLA4 inhibitor combined with other ICBs like PD-1 inhibitor nivolumab or LAG3 inhibitor may obtain a better therapeutic response in ccRCC, since that preclinical and clinical studies have provided evidence that combination inhibitor of CTLA4 and other ICBs can enhance the anti-tumor efficiency of CD8+ T cells [45, 46].


We comprehensively and systematically studied the prognostic value of CTLA4 in ccRCC and its impact on the landscape of TME TILs infiltration and genetic mutation, finding that CTLA4, acted like an oncogene, can accelerate the progression of ccRCC with a high prognostic value, and that CTLA4 was associated with more TILs infiltrated TME but had an immunosuppressed phenotype. Besides, patients with high CTLA4 levels may benefit more from the combined ICBs therapy. However, the potential role of CTLA4 in the progression of ccRCC needs further verification in vitro and in vivo.

Availability of data and materials

All analyzed data are included in this published article and its additional file. The original data used during the current study are available from the corresponding author on reasonable request.



Cytotoxic T-lymphocyte associated protein 4


Clear cell renal cell carcinoma


Tumor microenvironment


Tumor-infiltrating lymphocytes


Overall survival


BRCA-associated protein 1


Immune checkpoint blockades


T cells regulatory


Programmed death-1 receptor


Programmed death-ligand 1


The Cancer Genome Atlas


Gene Expression Omnibus


Formalin fixed paraffin embedded


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes




Hazard ratio


Hepatitis A virus cellular receptor 2


Lymphocyte activation gene-3


T cell immunoglobulin and ITIM domain


Tumor mutation burden




  1. 1.

    Inamura K. Renal cell tumors: understanding their molecular pathological epidemiology and the 2016 WHO classification. Int J MolSci. 2017;18(10):2195.

    Article  CAS  Google Scholar 

  2. 2.

    Owens B. Kidney cancer. Nature. 2016;537(7620):S97.

    CAS  Article  Google Scholar 

  3. 3.

    Posadas EM, Limvorasak S, Figlin RA. Targeted therapies for renal cell carcinoma. Nat Rev Nephrol. 2017;13(8):496–511.

    Article  Google Scholar 

  4. 4.

    Hui L, Chen Y. Tumor microenvironment: sanctuary of the devil. Cancer Lett. 2015;368(1):7–13.

    CAS  Article  Google Scholar 

  5. 5.

    Cervantes-Villagrana RD, Albores-Garcia D, Cervantes-Villagrana AR, Garcia-Acevez SJ. Tumor-induced neurogenesis and immune evasion as targets of innovative anti-cancer therapies. Signal Transd Target Ther. 2020;5(1):99.

    CAS  Article  Google Scholar 

  6. 6.

    Eun Y, Kim IY, Sun JM, Lee J, Cha HS, Koh EM, et al. Risk factors for immune-related adverse events associated with anti-PD-1 pembrolizumab. Sci Rep. 2019;9(1):14039.

    Article  CAS  Google Scholar 

  7. 7.

    Topalian SL, Taube JM, Anders RA, Pardoll DM. Mechanism-driven biomarkers to guide immune checkpoint blockade in cancer therapy. Nat Rev Cancer. 2016;16(5):275–87.

    CAS  Article  Google Scholar 

  8. 8.

    Pardoll DM. Immunology beats cancer: a blueprint for successful translation. Nat Immunol. 2012a;13(12):1129–32.

    CAS  Article  Google Scholar 

  9. 9.

    Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012b;12(4):252–64.

    CAS  Article  Google Scholar 

  10. 10.

    Antonia SJ, Villegas A, Daniel D, Vicente D, Murakami S, Hui R, et al. Overall survival with durvalumab after chemoradiotherapy in stage III NSCLC. N Engl J Med. 2018;379(24):2342–50.

    CAS  Article  Google Scholar 

  11. 11.

    Nishijima TF, Muss HB, Shachar SS, Moschos SJ. Comparison of efficacy of immune checkpoint inhibitors (ICIs) between younger and older patients: a systematic review and meta-analysis. Cancer Treat Rev. 2016;45:30–7.

    CAS  Article  Google Scholar 

  12. 12.

    van de Ven K, Borst J. Targeting the T-cell co-stimulatory CD27/CD70 pathway in cancer immunotherapy: rationale and potential. Immunotherapy. 2015;7(6):655–67.

    Article  CAS  Google Scholar 

  13. 13.

    Gough SC, Walker LS, Sansom DM. CTLA4 gene polymorphism and autoimmunity. Immunol Rev. 2005;204:102–15.

    CAS  Article  Google Scholar 

  14. 14.

    Peggs KS, Quezada SA, Chambers CA, Korman AJ, Allison JP. Blockade of CTLA-4 on both effector and regulatory T cell compartments contributes to the antitumor activity of anti-CTLA-4 antibodies. J Exp Med. 2009;206(8):1717–25.

    CAS  Article  Google Scholar 

  15. 15.

    Lee SJ, Jun SY, Lee IH, Kang BW, Park SY, Kim HJ, et al. CD274, LAG3, and IDO1 expressions in tumor-infiltrating immune cells as prognostic biomarker for patients with MSI-high colon cancer. J Cancer Res ClinOncol. 2018;144(6):1005–14.

    CAS  Article  Google Scholar 

  16. 16.

    Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    CAS  Article  Google Scholar 

  17. 17.

    Mayakonda A, Lin DC, Assenov Y, Plass C, Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 2018;28(11):1747–56.

    CAS  Article  Google Scholar 

  18. 18.

    Hao Y, Yan M, Heath BR, Lei YL, Xie Y. Fast and robust deconvolution of tumor infiltrating lymphocyte from expression profiles using least trimmed squares. PLoSComputBiol. 2019;15(5):e1006976.

    CAS  Google Scholar 

  19. 19.

    Kardos J, Chai S, Mose LE, Selitsky SR, Krishnan B, Saito R, et al. Claudin-low bladder tumors are immune infiltrated and actively immune suppressed. JCI Insight. 2016;1(3):e85902.

    Article  Google Scholar 

  20. 20.

    Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.

    Article  Google Scholar 

  21. 21.

    Zhang C, Li Z, Qi F, Hu X, Luo J. Exploration of the relationships between tumor mutation burden with immune infiltrates in clear cell renal cell carcinoma. Ann Transl Med. 2019;7(22):648.

    CAS  Article  Google Scholar 

  22. 22.

    Lawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499(7457):214–8.

    CAS  Article  Google Scholar 

  23. 23.

    Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9(1):34.

    Article  CAS  Google Scholar 

  24. 24.

    Zehir A, Benayed R, Shah RH, Syed A, Middha S, Kim HR, et al. Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nat Med. 2017;23(6):703–13.

    CAS  Article  Google Scholar 

  25. 25.

    Cheng Q, Li J, Fan F, Cao H, Dai ZY, Wang ZY, et al. Identification and analysis of glioblastoma biomarkers based on single cell sequencing. Front BioengBiotechnol. 2020;8:167.

    Google Scholar 

  26. 26.

    Chan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, et al. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol. 2019;30(1):44–56.

    CAS  Article  Google Scholar 

  27. 27.

    Gunjur A. Nivolumab plus ipilimumab in advanced renal-cell carcinoma. Lancet Oncol. 2018;19(5):e232.

    Article  Google Scholar 

  28. 28.

    Hammers HJ, Plimack ER, Infante JR, Rini BI, McDermott DF, Lewis LD, et al. Safety and efficacy of nivolumab in combination with ipilimumab in metastatic renal cell carcinoma: the CheckMate 016 Study. J ClinOncol. 2017;35(34):3851–8.

    CAS  Article  Google Scholar 

  29. 29.

    Sahtoe DD, van Dijk WJ, Ekkebus R, Ovaa H, Sixma TK. BAP1/ASXL1 recruitment and activation for H2A deubiquitination. Nat Commun. 2016;7:10292.

    CAS  Article  Google Scholar 

  30. 30.

    Planchard D, Reinmuth N, Orlov S, Fischer JR, Sugawara S, Mandziuk S, et al. ARCTIC: durvalumab with or without tremelimumab as third-line or later treatment of metastatic non-small-cell lung cancer. Ann Oncol. 2020;31(5):609–18.

    CAS  Article  Google Scholar 

  31. 31.

    Chen EX, Jonker DJ, Loree JM, Kennecke HF, Berry SR, Couture F, et al. Effect of combined immune checkpoint inhibition vs best supportive care alone in patients with advanced colorectal cancer: The Canadian Cancer Trials Group CO.26 Study. JAMA Oncol. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Liu F, Huang J, Liu X, Cheng Q, Luo C, Liu Z. CTLA-4 correlates with immune and clinical characteristics of glioma. Cancer Cell Int. 2020;20:7.

    CAS  Article  Google Scholar 

  33. 33.

    Blackburn SD, Shin H, Haining WN, Zou T, Workman CJ, Polley A, et al. Coregulation of CD8+ T cell exhaustion by multiple inhibitory receptors during chronic viral infection. Nat Immunol. 2009;10(1):29–37.

    CAS  Article  Google Scholar 

  34. 34.

    Munn DH. Indoleamine 2,3-dioxygenase, Tregs and cancer. Curr Med Chem. 2011;18(15):2240–6.

    CAS  Article  Google Scholar 

  35. 35.

    Huang J, Liu F, Liu Z, Tang H, Wu H, Gong Q, et al. Immune checkpoint in glioblastoma: promising and challenging. Front Pharmacol. 2017;8:242.

    Article  CAS  Google Scholar 

  36. 36.

    Romano E, Kusio-Kobialka M, Foukas PG, Baumgaertner P, Meyer C, Ballabeni P, et al. Ipilimumab-dependent cell-mediated cytotoxicity of regulatory T cells ex vivo by nonclassical monocytes in melanoma patients. ProcNatlAcadSci USA. 2015;112(19):6140–5.

    CAS  Article  Google Scholar 

  37. 37.

    Ma J, Zheng B, Goswami S, Meng L, Zhang D, Cao C, et al. PD1(Hi) CD8(+) T cells correlate with exhausted signature and poor clinical outcome in hepatocellular carcinoma. J Immunother Cancer. 2019;7(1):331.

    Article  Google Scholar 

  38. 38.

    Speiser DE, Ho PC, Verdeil G. Regulatory circuits of T cell function in cancer. Nat Rev Immunol. 2016;16(10):599–611.

    CAS  Article  Google Scholar 

  39. 39.

    Harbour JW, Onken MD, Roberson ED, Duan S, Cao L, Worley LA, et al. Frequent mutation of BAP1 in metastasizing uveal melanomas. Science. 2010;330(6009):1410–3.

    CAS  Article  Google Scholar 

  40. 40.

    Hakimi AA, Chen YB, Wren J, Gonen M, Abdel-Wahab O, Heguy A, et al. Clinical and pathologic impact of select chromatin-modulating tumor suppressors in clear cell renal cell carcinoma. EurUrol. 2013;63(5):848–54.

    Google Scholar 

  41. 41.

    Farhood B, Najafi M, Mortezaee K. CD8(+) cytotoxic T lymphocytes in cancer immunotherapy: a review. J Cell Physiol. 2019;234(6):8509–21.

    CAS  Article  Google Scholar 

  42. 42.

    Ruffo E, Wu RC, Bruno TC, Workman CJ, Vignali DAA. Lymphocyte-activation gene 3 (LAG3): the next immune checkpoint receptor. SeminImmunol. 2019;42:101305.

    CAS  Google Scholar 

  43. 43.

    Andrews LP, Marciscano AE, Drake CG, Vignali DA. LAG3 (CD223) as a cancer immunotherapy target. Immunol Rev. 2017;276(1):80–96.

    CAS  Article  Google Scholar 

  44. 44.

    Komiya T, Huang CH. Updates in the clinical development of epacadostat and other indoleamine 2,3-dioxygenase 1 inhibitors (IDO1) for human cancers. Front Oncol. 2018;8:423.

    Article  Google Scholar 

  45. 45.

    Mahoney KM, Rennert PD, Freeman GJ. Combination cancer immunotherapy and new immunomodulatory targets. Nat Rev Drug Discov. 2015;14(8):561–84.

    CAS  Article  Google Scholar 

  46. 46.

    Postow MA, Callahan MK, Wolchok JD. Immune checkpoint blockade in cancer therapy. J ClinOncol. 2015;33(17):1974–82.

    CAS  Article  Google Scholar 

Download references


We would like to thank all the researchers who participated in the study for their contributions.


This project was supported by the National Natural Science Foundation of China (Grant No. 81860276) and the Scientific Research Projects Fund Project of Hubei Health Commission (Grant No. WJ2019M179).

Author information




Conception: SL, FW; Acquisition of data: LZ, FD; Analysis, Interpretation, Validation of data: SL, FW, YW, YF, MY, DY; Statistical analysis: YZ, ZD; Methodology: SL, MY, DY; Writing—original draft: FW, WT; Writing—review & editing: SL, FW; Revising: YL, SL, FW; Supervision and Management: YL, YC. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yeqiang Liu or Yanxiang Cheng.

Ethics declarations

Ethics approval and consent to participate

Our study was approved by the Ethics Committee of Renmin Hospital of Wuhan University (NO.2018017).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests in this work.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

CTLA4-related protein coding genes by correlation analysis. Table S2. Comparison of CIBERSORT immune cell fractions.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liu, S., Wang, F., Tan, W. et al. CTLA4 has a profound impact on the landscape of tumor-infiltrating lymphocytes with a high prognosis value in clear cell renal cell carcinoma (ccRCC). Cancer Cell Int 20, 519 (2020).

Download citation


  • CTLA4
  • Clear cell renal cell carcinoma
  • Tumor microenvironment
  • Tumor-infiltrating lymphocytes
  • Prognosis
  • CD8+ T cells
  • Immune checkpoints