- Research article
- Open Access
Identification of an independent immune-genes prognostic index for renal cell carcinoma
BMC Cancer volume 21, Article number: 746 (2021)
Considerable evidence has indicated an association between the immune microenvironment and clinical outcome in ccRCC. The purpose of this study is to extensively figure out the influence of immune-related genes of tumors on the prognosis of patients with ccRCC.
Files containing 2498 immune-related genes were obtained from the Immunology Database and Analysis Portal (ImmPort), and the transcriptome data and clinical information relevant to patients with ccRCC were identified and downloaded from the TCGA data-base. Univariate and multivariate Cox regression analyses were used to screen out prognostic immune genes. The immune risk score model was established in light of the regression coefficient between survival and hub immune-related genes. We eventually set up a nomogram for the prediction of the overall survival for ccRCC. Kaplan-Meier (K-M) and ROC curve was used in evaluating the value of the predictive risk model. A P value of < 0.05 indicated statistically significant differences throughout data analysis.
Via differential analysis, we found that 556 immune-related genes were expressed differentially between tumor and normal tissues (p < 0. 05). The analysis of univariate Cox regression exhibited that there was a statistical correlation between 43 immune genes and survival risk in patients with ccRCC (p < 0.05). Through Lasso-Cox regression analysis, we established an immune genetic risk scoring model based on 18 immune-related genes. The high-risk group showed a bad prognosis in K-M analysis. (p < 0.001). ROC curve showed that it was reliable of the immune risk score model to predict survival risk (5 year over survival, AUC = 0.802). The model indicated satisfactory AUC and survival correlation in the validation data set (5 year OS, Area Under Curve = 0.705, p < 0.05). From Multivariate regression analysis, the immune-risk score model plays an isolated role in the prediction of the prognosis of ccRCC. Under multivariate-Cox regression analysis, we set up a nomogram for comprehensive prediction of ccRCC patients’ survival rate. At last, it was identified that 18 immune-related genes and risk scores were not only tremendously related to clinical prognosis but also contained in a variety of carcinogenic pathways.
In general, tumor immune-related genes play essential roles in ccRCC development and progression. Our research established an unequal 18-immune gene risk index to predict the prognosis of ccRCC visually. This index was found to be an independent predictive factor for ccRCC.
Renal cell carcinoma (RCC) is among one of the most prevalent malignancies affecting humanity worldwide. Its incidence rate has increased during the past 10 years, consisting of 2–3% of the whole newly diagnosed carcinoma cases . In histology, clear cell RCC (ccRCC) is the predominant RCC subtype, responsible for nearly 75% of total renal carcinoma cases . As considerable progress has been achieved in screening, diagnosing, and treating a variety of types of tumors through surgery and drug therapy [3,4,5], the clinical prognosis of ccRCC remains unsatisfactory [2, 6]. Thus, identifying several prognostic factors and targets is crucial to making the therapy and clinical outcomes of ccRCC patients better.
Immune evasion has recently attracted great interest as one of the fundamental characteristics of carcinoma . Immunotherapies, including immune checkpoint blockade, have produced astonishing results in the management of malignancies. Accumulating evidence shows immune-related components, including immune genes, antigens, and immune cells, contribute greatly to the occurrence and malignant progression of cancer and are valuable markers for cancer diagnosis and prognosis . Additionally, immune genes in tumor TME have great potential as prognostic biomarkers . However, IRG predictive models still need extensive study when it comes to ccRCC biology.
This study aims to reveal the distribution and pedigree of immune-related genes in patients with ccRCC and explore the influence of immune-related genes on the prognosis of ccRCC cases. In addition, we established an immune genetic risk score model and set the nomogram, which was used for predicting the prognosis of ccRCC.
First, through the ImmPort data-base, we gained the table of 2498 immune genes, then downloaded the transcriptome records. A total of 72 paracancerous tissues and 507 cancerous tissues were included in ccRCC cases from the TCGA data-base. Furthermore, the clinicopathological data of 507 ccRCC patients were gained as well, which include age, sex, pathological grading, tumor staging and TNM staging, vital status, and survival time (Table 1). At last, the correction of transcriptome records depends on the “LIMMA” software package in R software.
Function analysis of related genes
For exploring the principal biological processes of the selected hub genes, we conducted the Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) analysis. The enriched KEGG and GO terms were identified by DAVID (https://david.ncifcrf.gov/).
Survival analysis of hub genes and comparison of their expression levels
RCC cases’ clinical records from the TCGA data-base include survival time, vital status, and TNM staging (remove missing information cases). The survival analysis of the hub genes was carried out by the survival R software package. A log-rank test was used to detect the difference in overall survival. Survival curves were demonstrated by using the Kaplan-Meier method. P < 0.05 was considered statistically significant.
Gene functional-enrichment analysis
For studying the biological characteristics of renal carcinoma, we performed gene enrichment analysis (GSEA, version 3.0, the broad institute of MIT and Harvard, http://software.broadinstitute.org/gsea/downloads.jsp) between the cancerous tissues and paracancerous tissues. The number of permutations is 1000; collapse dataset to gene symbols is “false”, and permutation type is “phenotype”. Additional options selected included weighted enrichment statistics, and Signal2Noise metric was applied to ranking genes. The experimental group was composed of a high-expression set, and the control group was composed of a low-expression set. Gene set databases c2.cp.kegg.v7.0.symbols.gmt was applied to enrichment analysis. Cut-off criteria included gene set size > 500 and < 15, while nominal P value of < 0.05 and an FDR of < 0.25 were considered significant.
The analysis was carried out entirely through R statistical language version 3.6.1 (https://www.R-project.org). All of the tests had two sides, and a level of P < 0.05 was accepted as statistically significant. The continuous variables following normal distribution were compared by independent t-test, while those in skewness were compared by Mann-Whitney U test. In the light of the Pearson correlation coefficient, correlation matrices were schemed using R-software. We study the connection between OS and immune cell infiltration on the basis of the Kaplan-Meier curve, assessed with a log-rank test. The relevance between OS and immune cell infiltration was visualized by the K-M curve and further evaluated by log-rank test. Sensitivity and specificity in the predictive model of recurrence were analyzed by time-dependent ROC curves. The univariate regression model was applied in analyzing the influence of single-variable on survival. The LASSO-Cox regression models were used to identify the independent factors for survival. According to the Cox analysis, we used regression coefficients to build a nomogram.
Differentially expressed genes of ccRCC
Files containing 2498 immune-related genes were obtained from the ImmPort data-base. For analyzing differential expression, transcriptome records of 72 para cancer and 507 tumor tissues were obtained from the TCGA data-base. By differential expression analysis between cancer tissues and normal tissues, 556 differentially expressed immune genes were recognized, among which 402 immune genes were upregulated and the remaining 154 immune genes were downregulated (P < 0.05, Fig. 1A). Heat map of the topmost ten up-and down-regulated differentially expressed genes are displayed in Fig. 1B.
Functional annotation of differentially expressed genes in renal carcinoma
We learned about the biological properties of 556 DEIGS by KEGG and GO analysis. David’s results showed that the topmost 3 enrichment GO items of the upregulated genes were cAMP-mediated signals, humoral immune response and negative regulation of ERBB signaling pathway, while the topmost 3 enrichment GO components of down-regulated genes are lymphocyte activation, humoral immune response, and regulation of leukocyte mediated immunity (Fig. 2A). Through pathway enrichment analysis, it was found that the top 3 biological pathways enriched with upregulated genes were the JAK − STAT signaling pathway, PI3K − Akt signaling pathway and Rap1 signaling pathway. In comparison, the top 3 biological pathways enriched with down-regulated genes were the cytokine-cytokine receptor interaction, Th1 and Th2 cell differentiation, and JAK-STAT signaling pathway (Fig. 2B).
Construction of the immune-related prognostic model
We established a PPI network that was based on the differentially expressed genes and recognized 496 genes with over 50 para cancer nodes (Figure S1). To reveal relationships of these 496 DEIGs with the prognoses of patients with ccRCC, we identified 43 prognostic DEIGs through single variable Cox regression analysis (Table 2). KIRC data, downloaded from TCGA, were separated into 2 groups at random (training group: validation group, 1:1). After that, a lasso regression study was performed for the purpose of increasing robustness and selecting isolated indexes for survival in all according to the training group. At last, 18 DEIGs were obtained for the establishment of a prognostic indicator (Fig. 3A, B, Table 3). In the light of the risk index established, we divide patients into high-risk groups or low-risk ones (Fig. 3C). The differential expression of the model genes between high-and low-risk sufferers in the training set of ccRCC is shown in this heatmap (Fig. 3D). The Kaplan-Meier analysis indicated shorter overall survival among those high-risk sufferers in the validation and training set (p < 0.05, Fig. 4A, B, C). The ROC curve indicated better sensitivity and specificity in the risk model when used to predict survival risk (the AUC values of 5-year overall survival in the validation group and training validation group are 0.705 and 0.802, separately, Fig. 4D, E, F). In light of multivariate and univariate Cox regression analysis of age, sex, tumor staging, pathological grade, TNM stage and risk score, we determined if the immune risk score model was isolated of age, sex, tumor stage and other clinicopathological parameters. Among the single variable Cox models, age, pathological grading, tumor staging, T, M stage and high-risk scores are related to low survival rates (Fig. 5A). Of the multivariate Cox model, age, pathological grade and risk score were the only isolated predictors (Fig. 5B). For forecasting the prognosis of patients with ccRCC at 3 and 5 years after operations, we established another nomogram (Fig. 5C, D, E) based on the variables related to the overall survival rate ((OS)), namely age, sex, tumor staging pathological grade, TNM staging and risk score (Fig. 5C, D, E).
Correlation between 18 immune genes together with a risk score of immune gene and clinical prognosis
The scale of every gene model was determined at different pathological stagings. EDN1, GNAL1, and ICAM1 were remarkably correlated with the progression of ccRCC (Fig. 6). IFNG’s and XCL’s expressions were connected to the infiltration of CD4+, CD8+, and myeloid dendritic cells (Figure S2). The immune gene risk score was closely related to the grade, pathological stage, and clinical TNM stage (Fig. 7).
Gene set enrichment analysis results with hallmark genes of risk scores
To investigate the biological connection of risk scores in ccRCC development, we carried out a GSEA on the risk scores from a TCGA renal carcinoma group. It showed that high-risk scores are connected to IL6 JAK STAT3 SIGNALING, EPITHELIAL MESENCHYMAL TRANSITION, and WNT BETA CATENIN SIGNALING (Fig. 8).
ccRCC is a heterogeneous disease with different ethnic characteristics resulting from renal epithelial cells  and accounts for most RCC-related deaths . Although radical nephrectomy has been proven to be an effective treatment for local renal cancer. Many patients may experience development and metastasis after surgical resection. Given that targeted treatment for advanced and metastatic ccRCC has been fully developed, response to treatment is diverse . Identification of molecular mechanisms and relevant prognostic factors may be critical to the treatment of ccRCC . The prognosis of the tumor is closely connected with TME, particularly when considering the tumor immune microenvironment [14, 15]. Different types of cancers have diverse immune gene subpopulations. Therefore, investigating the immune gene subsets is vital for evaluating the risk and prognosis of ccRCC.
In the research, a large amount of specimen data was conducted to assess the immune genes of ccRCC in comprehensiveness and detail. We analyzed the expression of 2498 immune genes that were gained from the ImmPort database in ccRCC and normal tissues. Furthermore, we recognized and set up a risk score model for ccRCC via a single variable and LASSO-Cox regression analyses. The model was composed of 18 DEIGs (AVP, CCL22, CRP, CSF2, CXCL5, EDN1, GAL, GNAI1, GNRH1, HMOX1, ICAM1, IFNG, IL4, PDGFRA, PPY, TGFB1, TNFSF11 and XCL1). The K-M analysis showed that high-risk scores not only were associated with lower overall survival but also predicted advanced stage and higher pathological grade.
Several genes in the model have been studied in renal cell carcinoma. Arginine vasopressin (AVP) and its type 2 receptor (V2R) play an important role in regulating salt and water homeostasis. Activation of V2R can stimulate the proliferation of renal cell carcinoma (RCC) cell line in vitro . It has been shown that the increase of CXCL5 cytokines is associated with sunitinib resistance in renal cell carcinoma . CCL22, CRP, ICAM1, IFNG, PDGFRA, TGFB1 and TNFSF11 have been confirmed to be involved in the malignant progression and metastasis of renal cell carcinoma through different biological mechanisms [18,19,20,21,22,23,24]. EDN1 may serve as a promising prognostic and diagnostic biomarker for ccRCC . In addition, CXCL5, IL4 may be involved in the regulation of the immune microenvironment in renal cell carcinoma .
ccRCC immune models are typically established by screening immune-related lncRNAs. A new prognostic gene marker based on immune lncRNA in patients with KIRC patients was found . Zhao et al. integrated multiple levels of data to construct immune, inflammatory, or KIRC-oriented neighbor networks and KIRC-related gene directed networks. Their analysis showed that genes related to immune and inflammation have unique topological characteristics and related KIRC expression patterns in the networks. Furthermore, they identified five core clusters for constructing specific prognostic biomarkers for KIRC . Another study evaluated the prognostic value of individual gene expression by using TCGA data and ccRCC patient data . In the study, a predictive nomogram was generated. Independent prognostic factors were identified not only for examining the functional involvement of individual genes in vitro and in vivo RCC models but also for the assessment of OS and progression-free survival of patients with ccRCC in the first, fifth, and eighth year . Our study shed light on the role of immune-related genes in tumorigenesis and malignant development of ccRCC. Based on immune genes, we established a novel risk-score model consisting of immune genes and verified it to predict ccRCC prognosis. Our risk model showed excellent predictive performance in terms of prediction and may thus make contribution to developing novel prognostic indexes of ccRCC. Furthermore, we analyzed the expression profiles of the model genes in the pathological grade and stage of RCC. The model-associated immune genes strongly showed an association with immuno-infiltrating cells, which may be used for targeting clinical immunotherapy.
However, our research still has limitations. First, our study only included the expression profiles of a part of the Western population. Extensive sample sequencing data from other countries and races are needed to enhance our conclusions. Second, our results were based on the RNA sequencing results of entire tumor tissue, and the diversity of different cell compositions in the TME was not considered. Third, we only focused on transcriptional expression profile data. Gene methylation level, mutation level, and other equally essential data in tumor progression were not considered. These data are pivotal to exploring tumor progression.
The research provides a basis for the application of immune genes in the prognosis of ccRCC. It is credible for the immune gene risk score model to forecast the prognosis of ccRCC, which serves as an independent prognostic element for ccRCC patients. Our results may be helpful to in personalized treatment for patients with ccRCC and in exploring novel biomarkers for the targeted therapy of ccRCC.
Renal cell carcinoma
Clear cell renal cell carcinoma
The Cancer Genome Atlas
Kyoto Encyclopedia of Genes and Genomes
Database for Annotation, Visualization, and Integrated Discovery
Kidney renal clear cell carcinoma
Gene enrichment analysis
False Discovery Rate
Receiver Operating Characteristic
Least absolute shrinkage and selection operator
Differentially expressed immune genes
Area Under Curve
Ljungberg B, Bensalah K, Canfield S, Dabestani S, Hofmann F, Hora M, et al. EAU guidelines on renal cell carcinoma: 2014 update. Eur Urol. 2015;67(5):913–24.
Barata PC, Rini BI. Treatment of renal cell carcinoma: current status and future directions. CA Cancer J Clin. 2017;67(6):507–24. https://doi.org/10.3322/caac.21411.
Dagher J, Kammerer-Jacquet SF, Dugay F, Beaumont M, Lespagnol A, Cornevin L, et al. Clear cell renal cell carcinoma: a comparative study of histological and chromosomal characteristics between primary tumors and their corresponding metastases. Virchows Arch. 2017;471(1):107–15. https://doi.org/10.1007/s00428-017-2124-0.
Fernandez-Pello S, Hofmann F, Tahbaz R, Marconi L, Lam TB, Albiges L, et al. A systematic review and meta-analysis comparing the effectiveness and adverse effects of different systemic treatments for non-clear cell renal cell carcinoma. Eur Urol. 2017;71(3):426–36. https://doi.org/10.1016/j.eururo.2016.11.020.
Lin YW, Lee LM, Lee WJ, Chu CY, Tan P, Yang YC, et al. Melatonin inhibits MMP-9 transactivation and renal cell carcinoma metastasis by suppressing Akt-MAPKs pathway and NF-kappaB DNA-binding activity. J Pineal Res. 2016;60(3):277–90. https://doi.org/10.1111/jpi.12308.
Garcia JA, Rini BI. Recent progress in the management of advanced renal cell carcinoma. CA Cancer J Clin. 2007;57(2):112–25. https://doi.org/10.3322/canjclin.57.2.112.
Aparicio LMA, Fernandez IP, Cassinello J. Tyrosine kinase inhibitors reprogramming immunity in renal cell carcinoma: rethinking cancer immunotherapy. Clin Transl Oncol. 2017;19(10):1175–82. https://doi.org/10.1007/s12094-017-1657-7.
Yang J, Li H, Hu S, Zhou Y. ACE2 correlated with immune infiltration serves as a prognostic biomarker in endometrial carcinoma and renal papillary cell carcinoma: implication for COVID-19. Aging (Albany NY). 2020;12(8):6518–35. https://doi.org/10.18632/aging.103100.
Murphy KA, James BR, Guan Y, Torry DS, Wilber A, Griffith TS. Exploiting natural anti-tumor immunity for metastatic renal cell carcinoma. Hum Vaccin Immunother. 2015;11(7):1612–20. https://doi.org/10.1080/21645515.2015.1035849.
Wang XM, Lu Y, Song YM, Dong J, Li RY, Wang GL, et al. Integrative genomic study of Chinese clear cell renal cell carcinoma reveals features associated with thrombus. Nat Commun. 2020;11(1):739. https://doi.org/10.1038/s41467-020-14601-9.
Scelo G, Larose TL. Epidemiology and risk factors for kidney cancer. J Clin Oncol. 2018;7:JCO2018791905.
Li JK, Chen C, Liu JY, Shi JZ, Liu SP, Liu B, et al. Long noncoding RNA MRCCAT1 promotes metastasis of clear cell renal cell carcinoma via inhibiting NPR3 and activating p38-MAPK signaling. Mol Cancer. 2017;16(1):111. https://doi.org/10.1186/s12943-017-0681-0.
Wang C, Wang Y, Hong T, Cheng B, Gan S, Chen L, et al. Blocking the autocrine regulatory loop of Gankyrin/STAT3/CCL24/CCR3 impairs the progression and pazopanib resistance of clear cell renal cell carcinoma. Cell Death Dis. 2020;11(2):117. https://doi.org/10.1038/s41419-020-2306-6.
Bremnes RM, Al-Shibli K, Donnem T, Sirera R, Al-Saad S, Andersen S, et al. The role of tumor-infiltrating immune cells and chronic inflammation at the tumor site on cancer development, progression, and prognosis: emphasis on non-small cell lung cancer. J Thorac Oncol. 2011;6(4):824–33. https://doi.org/10.1097/JTO.0b013e3182037b76.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10. https://doi.org/10.1158/0008-5472.CAN-17-0307.
Sinha S, Dwivedi N, Tao S, Jamadar A, Kakade VR, Neil MO, et al. Targeting the vasopressin type-2 receptor for renal cell carcinoma therapy. Oncogene. 2020;39(6):1231–45. https://doi.org/10.1038/s41388-019-1059-0.
Giuliano S, Dufies M, Ndiaye PD, Viotti J, Borchiellini D, Parola J, et al. Resistance to lysosomotropic drugs used to treat kidney and breast cancers involves autophagy and inflammation and converges in inducing CXCL5. Theranostics. 2019;9(4):1181–99. https://doi.org/10.7150/thno.29093.
Berlato C, Khan MN, Schioppa T, Thompson R, Maniati E, Montfort A, et al. A CCR4 antagonist reverses the tumor-promoting microenvironment of renal cancer. J Clin Invest. 2017;127(3):801–13. https://doi.org/10.1172/JCI82976.
Jin C, Shi L, Li Z, Liu W, Zhao B, Qiu Y, et al. Circ_0039569 promotes renal cell carcinoma growth and metastasis by regulating miR-34a-5p/CCL22. Am J Transl Res. 2019;11(8):4935–45.
Liu Y, Shang D. Transforming growth factor-beta1 enhances proliferative and metastatic potential by up-regulating lymphoid enhancer-binding factor 1/integrin alphaMbeta2 in human renal cell carcinoma. Mol Cell Biochem. 2020;465(1–2):165–74. https://doi.org/10.1007/s11010-019-03676-8.
Robinson RL, Sharma A, Bai S, Heneidi S, Lee TJ, Kodeboyina SK, et al. Comparative STAT3-regulated gene expression profile in renal cell carcinoma subtypes. Front Oncol. 2019;9:72. https://doi.org/10.3389/fonc.2019.00072.
Shah A, Lal P, Toorens E, Palmer MB, Schwartz L, Vergara N, et al. Acquired cystic kidney disease-associated renal cell carcinoma (ACKD-RCC) harbor recurrent mutations in KMT2C and TSC2 genes. Am J Surg Pathol. 2020;44(11):1479–86. https://doi.org/10.1097/PAS.0000000000001530.
Steven A, Leisz S, Fussek S, Nowroozizadeh B, Huang J, Branstetter D, et al. Receptor activator of NF-kappaB (RANK)-mediated induction of metastatic spread and association with poor prognosis in renal cell carcinoma. Urol Oncol. 2018;36(11):502.e515–24.
Wen L, Guo L, Zhang W, Li Y, Jiang W, Di X, et al. Cooperation between the inflammation and coagulation systems promotes the survival of circulating tumor cells in renal cell carcinoma patients. Front Oncol. 2019;9:504. https://doi.org/10.3389/fonc.2019.00504.
Wang S, Yu ZH, Chai KQ. Identification of EGFR as a novel key gene in clear cell renal cell carcinoma (ccRCC) through bioinformatics analysis and meta-analysis. Biomed Res Int. 2019;2019:6480865.
Foguer K, Braga Mde S, Peron JP, Bortoluci KR, Bellini MH. Endostatin gene therapy inhibits intratumoral macrophage M2 polarization. Biomed Pharmacother. 2016;79:102–11. https://doi.org/10.1016/j.biopha.2016.01.035.
Khadirnaikar S, Kumar P, Pandi SN, Malik R, Dhanasekaran SM, Shukla SK. Immune associated LncRNAs identify novel prognostic subtypes of renal clear cell carcinoma. Mol Carcinog. 2019;58(4):544–53. https://doi.org/10.1002/mc.22949.
Zhao E, Li L, Zhang W, Wang W, Chan Y, You B, et al. Comprehensive characterization of immune- and inflammation-associated biomarkers based on multi-omics integration in kidney renal clear cell carcinoma. J Transl Med. 2019;17(1):177. https://doi.org/10.1186/s12967-019-1927-y.
Cancer Genome Atlas Research N. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013;499(7456):43–9. https://doi.org/10.1038/nature12222.
This work was supported by the Jiangsu Province “Six Talent Peaks Project” (WSN-011) by the National Natural Science Foundation of China (grant number 81672531, 81972386). The funders had no role in the design of this study, data collection and analysis, the decision to publish, or preparation of this manuscript.
Ethics approval and consent to participate
Consent for publication
The authors have no conflicts of interest to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, G., Wei, X., Su, S. et al. Identification of an independent immune-genes prognostic index for renal cell carcinoma. BMC Cancer 21, 746 (2021). https://doi.org/10.1186/s12885-021-08367-6
- Immune-related genes
- Clear cell renal cell carcinoma
- Immune risk score model