Skip to main content

A novel nine-microRNA-based model to improve prognosis prediction of renal cell carcinoma



With the improved knowledge of disease biology and the introduction of immune checkpoints, there has been significant progress in treating renal cell carcinoma (RCC) patients. Individual treatment will differ according to risk stratification. As the clinical course varies in RCC, it has developed different predictive models for assessing patient’s individual risk. However, among other prognostic scores, no transparent preference model was given. MicroRNA as a putative marker shown to have prognostic relevance in RCC, molecular analysis may provide an innovative benefit in the prophetic prediction and individual risk assessment. Therefore, this study aimed to establish a prognostic-related microRNA risk score model of RCC and further explore the relationship between the model and the immune microenvironment, immune infiltration, and immune checkpoints. This practical model has the potential to guide individualized surveillance protocols, patient counseling, and individualized treatment decision for RCC patients and facilitate to find more immunotherapy targets.


Downloaded data of RCC from the TCGA database for difference analysis and divided it into a training set and validation set. Then the prognostic genes were screened out by Cox and Lasso regression analysis. Multivariate Cox regression analysis was used to establish a predictive model that divided patients into high-risk and low-risk groups. The ENCORI online website and the results of the RCC difference analysis were used to search for hub genes of miRNA. Estimate package and TIMER database were used to evaluate the relationship between risk score and tumor immune microenvironment (TME) and immune infiltration. Based on Kaplan-Meier survival analysis, search for immune checkpoints related to the prognosis of RCC.


There were nine miRNAs in the established model, with a concordance index of 0.702 and an area under the ROC curve of 0.701. Nine miRNAs were strongly correlated with the prognosis (P < 0.01), and those with high expression levels had a poor prognosis. We found a common target gene PDGFRA of hsa-miR-6718, hsa-miR-1269b and hsa-miR-374c, and five genes related to ICGs (KIR2DL3, TNFRSF4, LAG3, CD70 and TNFRSF9). The immune/stromal score, immune infiltration, and immune checkpoint genes of RCC were closely related to its prognosis and were positively associated with a risk score.


The established nine-miRNAs prognostic model has the potential to facilitate prognostic prediction. Moreover, this model was closely related to the immune microenvironment, immune infiltration, and immune checkpoint genes of RCC.

Peer Review reports


In the latest update of GLOBOCAN worldwide cancer statistics, it was estimated that approximately 431,288 new cases of kidney cancer were diagnosed in 2020, with about 179,368 patients succumbing to the disease ( Renal cell carcinoma (RCC) accounts for 90% of kidney cancers. RCC can be broken down into three major histological subtypes: clear cell RCC (CCRCC), papillary RCC (PRCC), and chromophobe RCC (ChRCC) [1,2,3] . These three types of RCC account for 65–70%, 15–20%, and 5–7%, respectively. Patients with CCRCC have the worst disease-specific survival because it is discovered at a more advanced stage [4]. Considerable progress has been made in treating patients with RCC, with innovative immunotherapy making dramatic advances in the management of this disease [5]. At present, the main treatment methods for RCC are: first-line combined treatment of metastatic RCC combined with immune checkpoint blocking, combining immune checkpoint blockade with tyrosine kinase inhibitors, and combination therapy with anti–vascular endothelial growth factor antibodies [6]. Immune checkpoints (IC) are inhibitory regulators of the immune system that are crucial for maintaining self-tolerance. The IC pathways are mechanisms adopted by cancer cells to disguise themselves as regular components of the human body to undergo immune evasion [7]. At the same time, immune infiltration of tumors is closely associated with clinical outcomes in RCC [8, 9].

MicroRNAs (miRNAs) can regulate various target genes and are thus involved in the regulation of various of physiological and pathological processes, including cell proliferation, fundamental biological processes, signal pathways, tumor formation and development and so on [10,11,12,13]. According to their effects on tumor development, miRNAs are classified as carcinogenic or tumor-suppressive. As an essential determinant of post-transcriptional regulation, miRNA participates in the occurrence and progression of cancer and regulates anti-cancer immune response [14].

In this study, we are committed to finding miRNAs related to the prognosis of RCC and constructing a prognostic model through lasso and multivariate cox regression, further looking for the relationship between the prognostic model and immune infiltration/ checkpoints. This research results should contribute to the immunotherapy and prognosis of RCC.

Materials and methods

Data preparation

Gene expression profiles and miRNA expression profiles of RCC samples were downloaded from the TCGA data portal (, April 2021). Clinical data of RCC samples were downloaded from the website ( After excluding samples with incomplete information, 865 samples were included in the study, including 512 KIRC, 288 KIRP and 65 KICH. A list of 79 immune checkpoint genes (ICGs) was obtained from the literature (PMID: 32814346) [15]. The downloaded mRNA and miRNA raw data were integrated and standardized using the Edger software package, and differences were analyzed to obtain differentially expressed genes and their expression level. The data provided by TCGA are public and open-ended, and therefore does not require the approval of a local ethics committee.

Construction of prognostic model for RCC

Based on the results of the difference analysis, univariate cox analysis was carried out using the R (vision 4.0.4) package ‘survival’. L1-penalized (lasso) characterized by simultaneous variable selection and shrinkage is a valuable method for determining interpretable prediction rules in high-dimensional data [16]. Based on the ‘glmnet’ package in R, lasso cox regression analysis was applied to build an optimal prognostic signature for RCC samples. The optimal values of the penalty parameter λ were determined through 26 cross-validations. Finally, multivariate cox analyses were conducted to compare their predictive power of prognostic. An equation for calculating risk score was generated based on the expression levels of these prognostic miRNAs and their regression coefficients from the multivariate cox analyses as follows:

Risk score = coefmiRNA1 x exprmiRNA1 + coefmiRNA2 x exprmiRNA2+ · ···· + coefmiRNAn x exprmiRNAn.

According to the median, the RS values for each sample were calculated, and the samples were divided into high-risk and low-risk groups. The forest plot was drawn by ‘survminer’ in R language.

Survival analysis of prognostic model

Overall survival (OS) time and recurrence-free survival (RFS) time were used as variables to calculate the scores of 3 years and 5 years. Nomogram and calibration charts were drawn using ‘rms’ package in R language. The risk groups’ discrepancies in OS time and RFS time were analyzed using Kaplan-Meier (KM) survival analysis and the log-rank test. Survival analysis uses GraphPad Prism 8.4.2 analysis, and p-value < 0.05 was considered significant. ROC curve was used to estimate the prediction performance of the prognosis index and carried out using the R package ‘timeROC’.

Subgroup analysis of clinicopathological factors

The TNM stage, age, gender and other indicators of the patients were included in the study to analyze the influence of these factors on the prognosis of the patients by log-rank-test. And univariate and multivariate cox analyses were conducted to compare their predictive power of prognostic. These analyses used the R language and GraphPad Prism, and p-value < 0.05 was considered significant.

Construction of miRNA-mRNA regulatory network

To further study the interactions between miRNAs and hub genes, the Encyclopedia of RNA Interactomes (ENCORI)(; version 3.0) was utilized to predict the targeted miRNAs of hub genes. In our study, the targeted miRNAs of hub genes were defined according to the positive results of ≥3 miRNA-target predicting databases, including PicTar, TargetScan, PITA, and miRanda. After taking the intersection with the hub genes selected by the difference analysis in the training set, the Venn diagram was drawn on the Venny online website (Venny 2.1.0 ( Finally, the interaction network of miRNAs and hub genes was constructed using Cytoscape.

Correlation between risk score and immune infiltration

TIMER includes 10,897 samples across 32 cancer types from TCGA to estimate the abundance of six TIIC subsets (B cells, CD4 T cells, CD8 T cells, macrophages, neutrophils, and dendritic cells) [17]. The information of the six immune cells of RCC was downloaded from the TIMER website ( At the same time, the stromal/ immune score of RCC was downloaded from the ESTIMATE website ( The Wilcoxon test was used to compare the immune infiltration differences between high and low-risk groups. The log-rank test was used to analyze the relationship between stromal/immune score and OS of RCC. Finally, we analyzed the correlation between the risk score and the six immune cells.

Correlation between risk score and immune checkpoint genes

Immune checkpoint genes (ICGs) play critical roles in circumventing self-reactivity and represent a novel target to develop treatments for cancers [18]. Using the ‘vlookup’ function in Excel, 76 immune checkpoint genes obtained from the literature were matched with the gene expression profiles of RCC. Differential analysis was done using the Edger package with p-value < 0.05 and HR > 1. The log-rank test was used to analyze the relationship between ICGs and OS of RCC, and the Wilcoxon test was used to compare the ICGs between high and low-risk groups.


Difference analysis and patients grouping

Based on the TCGA database, 63 differentially expressed miRNAs (DEMs) (false discovery rate (FDR) < 0.05 and |log2fold change (log2FC) | ≥ 2) were identified using 1021 miRNAs expression profiles between 893 RCC samples and 128 adjacent normal kidney tissues (Fig. 1). These 1021 samples were matched with the clinical data downloaded from cBioPortal, and finally, 865 patients with complete clinical data were included in the study. The patients were randomly divided into a training set containing 451 samples and a validation set containing 414 samples using the caret package.

Fig. 1
figure 1

Difference analysis. A. The heatmap of top 20 miRNAs by difference analysis. B. Volcano plot of miRNAs in RCC patients. Red color indicates up-regulated expression, and green color represents down-regulated expression

Prognostic model construction and risk score formula

First, 26 miRNAs related to prognosis were selected through univariate cox analysis (Table S1). Based on the expression of these 26 overlapped miRNAs in the training set, the lasso analysis identified a 9-miRNA signature that was significantly associated with survival rate based on the optimal lambda value (Fig. 2). The Concordance index of the model was 0.702, and the 9 miRNAs were hsa-miR-105-2, hsa-miR-122, hsa-miR-1269a, hsa-miR-1269b, hsa-miR-1293, hsa-miR-155, hsa-miR-224, hsa-miR-374c and hsa-miR-6718. These nine miRNAs were all up-regulated expressed in RCC. Hsa-miR-105-2 was the most statistically significant with p < 0.001(Table 1). To facilitate the utility of the identified prognostic miRNAs in routine clinical practice, the following formula was developed to generate risk score for each patient: Risk score = 1.42033 × exprhsa-miR-105-2 + 1.01781 × exprhsa-miR-122 + 1.03488 × exprhsa-miR-1269a + 1.09167 × exprhsa-miR-1269b + 1.09789 × exprhsa-miR-1293 + 1.07460 × exprhsa-miR-155 + 0.96770 × exprhsa-miR-224+ 1.11368 × exprhsa-miR-374c + 1.17308 × exprhsa-miR-6718. Therefore, patients were divided into low-risk and high-risk groups by the same median risk score as the cut-off points in the two independent cohorts.

Fig. 2
figure 2

Prognostic model construction A. 26 cross-validations in Lambda plot. B. min plot. C. Forest plot of 9 miRNAs. ‘***’ p < 0.001, ‘**’ p < 0.01, ‘*’ p < 0.05

Table 1 Nine miRNAs in the prognostic model

Survival analysis

Based on KM survival analysis, the high-risk group has a significantly worse prognosis (Fig. 3F). To further analyze the effect of the prognostic model, a 3-year and 5-year nomogram was drawn (Fig. 3A). In addition, calibration plots of the prognostic model fitted well in the training set, which indicated good calibration ability (Fig. 3B-C). The area under the ROC curve for 3-year survival in the prognostic model was 0.686. The area under the 5-year ROC curve was 0.705, which indicated that the prognostic model with low to moderate diagnostic ability (Fig. 3D). At the same time, we also conducted a separate survival analysis for these nine miRNAs. The results were that the up-regulated expressed miRNAs have a poor prognosis for RCC (Fig. 3E).

Fig. 3
figure 3

Survival analysis related to prognostic model. A. Nine-miRNA-based prognostic model to predict 3 and 5-year overall survival (OS) in RCC patients. B-C. Calibration plots of the nine-miRNA-based prognostic model of 3 and 5-year OS. D. ROC curve of 3 and 5-year OS. E. Impact of nine miRNAs on OS in RCC based on KM analysis. F. Impact of a nine-miRNA-based prognostic model in training set on OS in RCC based on KM analysis

Correlation between prognostic index and clinicopathological factors

The clinical data of 451 patients were presented in Table 2, and based on KM survival analysis, the relationship between clinicopathological factors and OS was drawn (Fig. 4A). Clinicopathological factors had an impact on the prognosis, except for gender. Univariate and multivariate analysis further proved that our index could be used as an independent prognostic index and had a higher predictive value than the traditional tumor pathological stage (p < 0.05, Fig. 4B-C). (To ensure the accuracy of the results, Mstage was not included in the Cox regression analysis as lacking too much data.)

Table 2 Clinicopathological characteristics of RCC patients in training set
Fig. 4
figure 4

Correlation with clinicopathological factors. A. Impact of clinicopathological factors on OS in RCC based on KM analysis. B. Univariate analysis of RS and clinicopathological factors. C. Multivariate analysis of RS and clinicopathological factors

miRNA-mRNA regulatory network

Four miRNAs were selected after multivariate analysis. First, we found 1488 target genes for hsa-miR-374c, 659 target genes for hsa-miR-1269b, 2536 target genes for hsa-miR-6718, and no target genes for hsa-mirR-105-2 from ENCORI. Then, venny2.1 was used to obtain the intersection between the target gene of miRNA and the up-regulated and down-regulated genes in RCC (Fig. 5A). Finally, Cytoscape was used to make the miRNA-mRNA regulatory network, and the top 50 hub genes were drawn on the network (Fig. 5B-C). PDGFRA was the common hub gene between the three miRNAs and down-regulated genes.

Fig. 5
figure 5

miRNA-mRNA network. A. The intersection of the hub gene predicted by the ENCORI website and the difference analysis result of RCC patients. B-C. The network between miRNA and the up-regulated / down-regulated genes in the intersection

Relationship between prognostic index and tumor microenvironment (TME)

We can find that the immune score was related to survival. The median survival time of the high immune/stromal score group was shorter than that of the low immune score group (Fig. 6B). And the immune/stromal score in the high-risk group was higher than that of the low-risk group (Fig. 6A). The box plot shows that the content of B cells, T cells, CD8 +, neutrophils, and myeloid dendritic cells are significantly positively correlated with RS (p < 0.05, Fig. 6C). Neutrophils and myeloid dendritic cells have the strongest correlation (Fig. 6D).

Fig. 6
figure 6

Correlation between immune infiltration and risk score. A. The relationship between RS and immune/stromal score. B. Impact of immune/ stromal score on OS in RCC based on KM analysis. C. The relationship between RS and immune infiltration cells. D. Correlation between RS and Neutrophil/ Myeloid dendritic cell

Association between prognostic index and immune checkpoint gene

After taking the intersection between 76 ICGs and the mRNA of RCC patients, The Edger package was used for differential analysis. A total of six ICGs were screened out, KIR2DL3, TNFRSF4, LAG3, CD70 and TNFRSF9 were up-regulated genes, while ICOSLG was a down-regulated gene. Based on KM survival analysis, five up-regulated genes were found to be related to the OS, and up-regulation of gene expression indicated poor prognosis (Fig. 7A). In addition, the high-risk group of ICGs expression was higher than that of the low-risk group, which was consistent with the results of the prognosis model (Fig. 7B).

Fig. 7
figure 7

Correlation between Immune checkpoint genes and risk score. A. Impact of five immune checkpoint genes on OS in RCC based on KM analysis. B. The relationship between RS and Immune checkpoint genes

Validation of the miRNA signature by the validation set

To further test the ability of the prognostic model, the data in the validation set was calculated using the risk score formula and divided into high and low score groups based on the median. Based on Kaplan-Meier survival analysis, the high-risk group has a significantly worse prognosis (Fig. 8). The results were consistent with the training set, confirming the accuracy of the prognostic model.

Fig. 8
figure 8

Impact of a nine-miRNA-based prognostic model in the validation set on OS in RCC based on KM analysis


There is an urgent need for a molecular marker-based method to accurately predict RCC patients’ survival rates. The major established risk factors for RCC based on clinical factors include hypertension, excess body weight, and smoking [19]. Increasing evidence shows that miRNAs play a vital role in RCC prognosis. In this study, we established a prognostic model of RCC containing nine miRNAs, and all these miRNAs have been shown to have a harmful effect on survival time. Their high expression leads to a shorter survival time.

In one literature report, we only found one target gene of hsa-miR-105-2 as CSAG1 [20]. Chondrosarcoma-associated gene one protein (CSAG1) depletion results in delayed mitotic progression and multipolar mitotic exit. Depletion of CSAG1 disrupts centrosomes and leads to multipolar spindles, particularly in cells with compromised p53 function [21]. It has been reported that high hsa-miR-122 expression promotes the malignant phenotype of clear cell RCC (ccRCC) by targeting occludin [22]. hsa-miR-1269b can promote the proliferation and migration of hepatocellular carcinoma, drives Cisplatin Resistance of Human Non-Small Cell Lung Cancer (NSCLC) [23,24,25]. Few studies related to hsa-miR-1269a and cancer were found. It was reported that hsa-miR-155 and hsa-miR-224 worked as cancer suppressors in RCC [26,27,28,29]. hsa-miR-374c-5p plays a crucial role in the invasion and migration of cervical cancer by acting on the Foxc1/snail pathway [30]. hsa-miR-6718 may be related to DNMT3A-mutant acute myeloid leukemia patients [31]. Due to the inhibitory and regulatory effects of miRNAs on downstream mRNAs, in recent decades, scientists have been committed to developing robust and safe targeting methods to restore these inhibitory miRNAs in cancer cells [32]. miRNAs play the role of a sponge to silence mRNA. Perhaps these miRNAs can be used as therapeutic targets to inhibit tumor growth.

Based on KM survival analysis, it can be found that the model can predict the prognosis of RCC. The survival time of the high-risk group was much shorter than that of the low-risk group. Moreover, the high expression of the nine miRNAs in the model all predicts a poor prognosis. The area under the ROC curve suggests that the model with moderate diagnostic ability. Clinicopathological factors had an impact on the prognosis, except for gender. Univariate and multivariate analysis further proved that our index could be an independent prognostic index and had a higher predictive value than the traditional tumor pathological stage. Moreover, the three-year and five-year calibration curves had good accuracy and consistency, which all prove that the prediction performance of this model was better. Based on the results of multivariate analysis, we expected the target genes of 4 miRNAs with p < 0.05 and found that PDGFRA was the common down-regulated gene of hsa-miR-374c, hsa-miR-6718, and hsa-miR-1269b. This gene may be used as a breakthrough point in the mechanism of miRNA-mRNA for RCC, and improve the prognosis of patients by inhibiting upstream miRNAs. Platelet-derived growth factor receptor alpha (PDGFRA) encodes a cell surface tyrosine kinase receptor for the platelet-derived growth factor family members. Studies suggest that this gene plays a role in organ development, wound healing, and tumor progression, whereas mutant-PDGFRA is potently oncogenic [33, 34].

RCC can be treated with ablation (destruction of the malignant tissue with heat or cold), partial or radical nephrectomy (removal of the kidney), or active surveillance (monitoring of tumor growth with periodic radiographic studies) [35]. In recent years, the approval of immunotherapeutic drugs and immunotherapy-based combination strategies has revolutionized the treatment of patients with advanced renal cell carcinoma (aRCC) [36]. In order to study the impact of TME and immune infiltration on RCC, we conducted KM analysis and correlation analysis. It was found that the group with a high immune/stromal score indicated a poor prognosis, which was positively correlated with RS. Furthermore, Infiltration levels of B cell, T cell CD8 +, Neutrophil, and Myeloid dendritic cell were higher in a high-risk group and were significantly positively correlated with RS. Neutrophils are recruited into the tumor microenvironment by the chemokine receptor system and play a role in tumor progression and metastasis [37]. The cytokines and chemokines secreted by the TME and tumor-associated immune cells promote a sustained and non-resolving tumor-associated inflammation. In addition, regulatory B cells increase tumor activity by regulating immune cells [38]. These may be the reason why the increase of these cells was negatively correlated with the prognosis. It was consistent with previous studies on immune infiltration.

ICGs play a key role in evading self-reactions and represent a new goal for t developing of cancer treatments. For example, Nivolumab, a programmed death 1 (PD-1) immune checkpoint inhibitor monoclonal antibody, was approved as a monotherapy for aRCC after treatment with a VEGF-targeted agent [39]. In our study, a total of 5 immune checkpoints related to RCC were found, and all of them were up-regulated genes. Among them, LAG3 and TNFRSF9 have the strongest correlation with prognosis. Lymphocyte activating 3 (LAG3) belongs to the Ig superfamily and contains four extracellular Ig-like domains and represents an inhibitory receptor, mainly found on activated immune cells and involved in the exhaustion of T cells in malignant diseases. In cancer, T cells are continuously exposed to antigens, leading to the gradual loss of cytokine production and the ability of CD8+ T cells to specifically kill tumor cells [40]. LAG3 is expressed on regulatory T cells in tumors and induces the production of IL-10 and TGF-β1, which helps tumor immune escape. Therefore, blocking LAG3 on regulatory T cells may reduce their inhibitory function and lead to a revival of CD8 + tumor-infiltrating lymphocytes activity [41]. Therefore, LAG3 is suitable as a clinical target to block cancer progression [42]. In our study, the high expression of LAG3 indicates a poor prognosis, and it is positively correlated with RS, which is consistent with these studies. TNF receptor superfamily member 9(TNFRSF9) plays a vital role in regulating the function of CD8+ T cells. It has been reported that TNFRSF9+ CD8+ T cells have exhaustion and effector phenotypes, identified as poor prognostic factors for ccRCC [43]. In our study, TNFRSF9 and T CD8+ cell scores were higher in the high-risk group than in the low-risk group. This may be the reason for the poor prognosis of the high-risk group. The influence of CD70 molecule (CD70), TNF receptor superfamily member 4 (TNFRSF4) and killer cell immunoglobulin like receptor, two Ig domains and long cytoplasmic tail 3 (KIR2DL3) on the prognosis of RCC is not so obvious as the first two genes. In most studies, these three genes need to work together with other genes, and they need to be verified as immunotherapy targets.

In the study, we constructed a nine-miRNAs prognostic model, which has been successfully verified in the validation data set and shows reliable predictive performance and significant clinical value in a comprehensive analysis with clinicopathological factors. We found PDGFRA, a common target gene of hsa-miR-6718, hsa-miR-1269b and hsa-miR-374c, and five genes related to ICGs.

In conclusion, we successfully constructed a nine-miRNAs prognostic model with good predictive power. Moreover, this model was closely related to the TME, immune infiltration and ICGs of RCC, and it is hoped that it can provide new ideas for RCC immunotherapy.

Availability of data and materials

The data that support the findings of this study are available in The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO).



Renal cell carcinoma


The cancer genome atlas


Tumor microenvironment




Risk score


Overall survival


Hazard ratio


Immune checkpoint genes


  1. Motzer RJ, Jonasch E, Agarwal N, Bhayani S, Bro WP, Chang SS, et al. Kidney cancer, version 2.2017, NCCN clinical practice guidelines in oncology. J Natl Compr Canc Netw. 2017;15(6):804–34.

  2. Capitanio U, Bensalah K, Bex A, Boorjian SA, Bray F, Coleman J, et al. Epidemiology of renal cell carcinoma. Eur Urol. 2019;75(1):74–84.

    PubMed  Google Scholar 

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

    Article  PubMed Central  Google Scholar 

  4. Jonasch E, Walker CL, Rathmell WK. Clear cell renal cell carcinoma ontogeny and mechanisms of lethality. Nat Rev Nephrol. 2021;17(4):245–61.

    Article  CAS  PubMed  Google Scholar 

  5. Barata PC, Rini BI. Treatment of renal cell carcinoma: current status and future directions. CA Cancer J Clin. 2017;67(6):507–24.

    Article  PubMed  Google Scholar 

  6. Chowdhury N, Drake CG. Kidney cancer: an overview of current therapeutic approaches. Urol Clin North Am. 2020;47(4):419–31.

    Article  PubMed  Google Scholar 

  7. Zhang JC, Chen WD, Alvarez JB, Jia K, Shi L, Wang Q, et al. Cancer immune checkpoint blockade therapy and its associated autoimmune cardiotoxicity. Acta Pharmacol Sin. 2018;39(11):1693–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Bussard KM, Mutkus L, Stumpf K, Gomez-Manzano C, Marini FC. Tumor-associated stromal cells as key contributors to the tumor microenvironment. Breast Cancer Res. 2016;18(1):84.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Rupaimoole R, Slack FJ. MicroRNA therapeutics: towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. 2017;16(3):203–22.

    Article  CAS  PubMed  Google Scholar 

  11. Ferragut Cardoso AP, Udoh KT, States JC. Arsenic-induced changes in miRNA expression in cancer and other diseases. Toxicol Appl Pharmacol. 2020;409:115306.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zhao XQ, Ao Y, Chen HY, Wang H. The role of miRNA in kidney development. Yi chuan = Hereditas. 2020;42(11):1062–72.

    PubMed  Google Scholar 

  13. Ghafouri-Fard S, Shirvani-Farsani Z, Branicki W, Taheri M. MicroRNA signature in renal cell carcinoma. Front Oncol. 2020;10:596359.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Yi M, Xu L, Jiao Y, Luo S, Li A, Wu K. The role of cancer-derived microRNAs in cancer immune escape. J Hematol Oncol. 2020;13(1):25.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Hu FF, Liu CJ, Liu LL, Zhang Q, Guo AY. Expression profile of immune checkpoint genes and their roles in predicting immunotherapy response. Brief Bioinform. 2021;22(3):bbaa176.

  16. Zhang Y, Li H, Zhang W, Che Y, Bai W, Huang G. LASSObased CoxPH model identifies an 11lncRNA signature for prognosis prediction in gastric cancer. Mol Med Rep. 2018;18(6):5579–93.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Marin-Acevedo JA, Dholaria B, Soyano AE, Knutson KL, Chumsri S, Lou Y. Next generation of immune checkpoint therapy in cancer: new developments and challenges. J Hematol Oncol. 2018;11(1):39.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Chow WH, Dong LM, Devesa SS. Epidemiology and risk factors for kidney cancer. Nat Rev Urol. 2010;7(5):245–57.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Liu GW, Qin ZM, Shen QH. An ensemble method integrated with miRNA expression data for predicting miRNA targets in stomach adenocarcinoma. Cancer Biomark. 2017;20(4):617–25.

    Article  CAS  PubMed  Google Scholar 

  21. Sapkota H, Wren JD, Gorbsky GJ. CSAG1 maintains the integrity of the mitotic centrosome in cells with defective p53. J Cell Sci. 2020;133(10):jcs239723.

  22. Jingushi K, Kashiwagi Y, Ueda Y, Kitae K, Hase H, Nakata W, et al. High miR-122 expression promotes malignant phenotypes in ccRCC by targeting occludin. Int J Oncol. 2017;51(1):289–97.

    Article  CAS  PubMed  Google Scholar 

  23. Yang W, Xiao W, Cai Z, Jin S, Li T. miR-1269b drives cisplatin resistance of human non-small cell lung Cancer via modulating the PTEN/PI3K/AKT signaling pathway. OncoTargets Ther. 2020;13:109–18.

    Article  CAS  Google Scholar 

  24. Kong XX, Lv YR, Shao LP, Nong XY, Zhang GL, Zhang Y, et al. HBx-induced MiR-1269b in NF-κB dependent manner upregulates cell division cycle 40 homolog (CDC40) to promote proliferation and migration in hepatoma cells. J Transl Med. 2016;14(1):189.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Wang X, Gao J, Zhou B, Xie J, Zhou G, Chen Y. Identification of prognostic markers for hepatocellular carcinoma based on miRNA expression profiles. Life Sci. 2019;232:116596.

    Article  CAS  PubMed  Google Scholar 

  26. Meng L, Xing Z, Guo Z, Qiu Y, Liu Z. Hypoxia-induced microRNA-155 overexpression in extracellular vesicles promotes renal cell carcinoma progression by targeting FOXO3. Aging. 2021;13(7):9613–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. García-Donas J, Beuselinck B, Inglada-Pérez L, Graña O, Schöffski P, Wozniak A, et al. Deep sequencing reveals microRNAs predictive of antiangiogenic drug response. JCI insight. 2016;1(10):e86051.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zhang J, Ye Y, Chang DW, Lin SH, Huang M, Tannir NM, et al. Global and targeted miRNA expression profiling in clear cell renal cell carcinoma tissues potentially links miR-155-5p and miR-210-3p to both tumorigenesis and recurrence. Am J Pathol. 2018;188(11):2487–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Pan Y, Hu J, Ma J, Qi X, Zhou H, Miao X, et al. MiR-193a-3p and miR-224 mediate renal cell carcinoma progression by targeting alpha-2,3-sialyltransferase IV and the phosphatidylinositol 3 kinase/Akt pathway. Mol Carcinog. 2018;57(8):1067–77.

    Article  CAS  PubMed  Google Scholar 

  30. Huang Y, Huang H, Li M, Zhang X, Liu Y, Wang Y. MicroRNA-374c-5p regulates the invasion and migration of cervical cancer by acting on the Foxc1/snail pathway. Biomed Pharmacother = Biomedecine & pharmacotherapie. 2017;94:1038–47.

    Article  CAS  Google Scholar 

  31. Lauber C, Correia N, Trumpp A, Rieger MA, Dolnik A, Bullinger L, et al. Survival differences and associated molecular signatures of DNMT3A-mutant acute myeloid leukemia patients. Sci Rep. 2020;10(1):12761.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Mollaei H, Safaralizadeh R, Rostami Z. MicroRNA replacement therapy in cancer. J Cell Physiol. 2019;234(8):12369–84.

    Article  CAS  PubMed  Google Scholar 

  33. CCL C, Deshmukh S, Jessa S, Hadjadj D, Lisi V, Andrade AF, et al. Histone H3.3G34-Mutant Interneuron Progenitors Co-opt PDGFRA for Gliomagenesis. Cell. 2020;183(6):1617–1633.e1622.

    Article  Google Scholar 

  34. Evans EK, Gardino AK, Kim JL, Hodous BL, Shutes A, Davis A, et al. A precision therapy against cancers driven by KIT/PDGFRA mutations. Sci Transl Med. 2017;9(414):eaao1690.

  35. Hsieh JJ, Purdue MP, Signoretti S, Swanton C, Albiges L, Schmidinger M, et al. Renal cell carcinoma. Nat Rev Dis Primers. 2017;3:17009.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Rini BI, Battle D, Figlin RA, George DJ, Hammers H, Hutson T, et al. The society for immunotherapy of cancer consensus statement on immunotherapy for the treatment of advanced renal cell carcinoma (RCC). J Immunother Cancer. 2019;7(1):354.

    Article  PubMed  PubMed Central  Google Scholar 

  37. SenGupta S, Subramanian BC, Parent CA. Getting TANned: how the tumor microenvironment drives neutrophil recruitment. J Leukoc Biol. 2019;105(3):449–62.

    Article  CAS  PubMed  Google Scholar 

  38. Tokunaga R, Naseem M, Lo JH, Battaglin F, Soni S, Puccini A, et al. B cell and B cell-related pathways for novel cancer treatments. Cancer Treat Rev. 2019;73:10–9.

    Article  CAS  PubMed  Google Scholar 

  39. Hayashi H, Nakagawa K. Combination therapy with PD-1 or PD-L1 inhibitors for cancer. Int J Clin Oncol. 2020;25(5):818–30.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Solinas C, Migliori E, De Silva P, Willard-Gallo K. LAG3: the biological processes that motivate targeting this immune checkpoint molecule in human cancer. Cancers. 2019;11(8):1213.

  42. Puhr HC, Ilhan-Mutlu A. New emerging targets in cancer immunotherapy: the role of LAG3. ESMO Open. 2019;4(2):e000482.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Li Y, Wang Z, Jiang W, Zeng H, Liu Z, Lin Z, et al. Tumor-infiltrating TNFRSF9(+) CD8(+) T cells define different subsets of clear cell renal cell carcinoma with prognosis and immunotherapeutic response. Oncoimmunology. 2020;9(1):1838141.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The results published here are in whole or part based upon data generated by The Cancer Genome Atlas managed by the NCI and NHGRI. Information about TCGA can be found at


This work was supported by grants from National Natural Science Foundation of China (81272372) and by Zhongnan Hospital of Wuhan University Science, Technology and Innovation Seed Fund (znpy2016033).

Author information

Authors and Affiliations



XC, FW and XL took part in research design. XC, HZ and XL drafted the manuscript. XC, JF, WH, XY and SL participated in the data analysis. XL took part in research funding. All the authors reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Fubing Wang or Xinghua Long.

Ethics declarations

Ethics approval and consent to participate

Not applicable. All data in this study are publicly available. All methods were carried out in accordance with the Declaration of Helsinki. All utilized public data sets were generated by others who obtained ethical approval.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest in relation to the submission.

Additional information

Publisher’s Note

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

Supplementary Information

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, C., Zeng, H., Fan, J. et al. A novel nine-microRNA-based model to improve prognosis prediction of renal cell carcinoma. BMC Cancer 22, 264 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: