- Open Access
A novel ferroptosis-related gene signature associated with cell cycle for prognosis prediction in patients with clear cell renal cell carcinoma
BMC Cancer volume 22, Article number: 1 (2022)
It is of great urgency to explore useful prognostic markers for patients with clear cell renal cell carcinoma (ccRCC). Prognostic models based on ferroptosis-related gene (FRG) in ccRCC is poorly reported for now.
Comprehensive analysis of 22 FRGs were performed in 629 ccRCC samples from two independent patient cohorts. We carried out least absolute shrinkage and selection operator analysis to screen out prognosis-related FRGs and constructed prognosis model for patients with ccRCC. Weighted gene co-expression network analysis was also carried out for potential functional enrichment analysis.
Based on the TCGA cohort, a total of 11 prognosis-associated FRGs were selected for the construction of the prognosis model. Significantly differential overall survival (hazard ratio = 3.61, 95% CI: 2.68–4.87, p < 0.0001) was observed between patients with high and low FRG score in the TCGA cohort, which was further verified in the CPTAC cohort with hazard ratio value of 5.13 (95% CI: 1.65–15.90, p = 0.019). Subgroup survival analysis revealed that our FRG score could significantly distinguish patients with high survival risk among different tumor stages and different tumor grades. Functional enrichment analysis illustrated that the process of cell cycle, including cell cycle-mitotic pathway, cytokinesis pathway and nuclear division pathway, might be involved in the regulation of ccRCC through ferroptosis.
We developed and verified a FRG signature for the prognosis prediction of patients with ccRCC, which could act as a risk factor and help to update the tumor staging system when integrated with clinicopathological characteristics. Cell cycle-related pathways might be involved in the regulation of ccRCC through ferroptosis.
It is reported that in the United States, there will be 76,080 new renal malignancy cases and about 13,780 cases of death related to renal carcinoma in 2021 . While in China, 74,000 new cases and a crude mortality of 1.97/105 for renal cancer were estimated in 2015 . Significant heterogeneity is found in different types of renal malignancy, among which, clear cell renal cell carcinoma (ccRCC) accounts for about 80% of malignant cases in renal . As one of the most aggressive malignancies, ccRCC is responsible for most of the death cases caused by renal tumor . Even for localized cases, about 25% of patients with ccRCC could also be troubled by tumor recurrence after receiving operative treatment . Tumor staging system is currently the most fashionable method for survival prediction of patients with ccRCC. However, different survival outcomes could also be found in patients with similar tumor staging. Therefore, it is of great urgency to explore useful prognostic markers and develop novel prognostic models for patients with ccRCC.
Ferroptosis is a newfound process of programmed cell death, which differs from the traditional cell death processes since it is caused by the lethal accumulation of iron-dependent lipid hydroperoxides . Current studies have also reported the important role of ferroptosis-related gene (FRG) in ccRCC. Through facilitating ferroptosis, SUV39H1 deficiency could restrain cell growth of ccRCC in vitro and in vivo . Reduced expression of NCOA4, which is one of the FRG, was reported to be associated with tumor progression and poor prognosis of ccRCC . In addition, cell density-regulated ferroptosis was found to be regulated via TAZ in cell death of renal cancer . However, the prognostic model based on FRG in ccRCC is poorly reported for now.
Here, we preformed comprehensive analysis of FRG from two independent patient cohorts to develop and verify a prognostic model based on FRG and explored the potential mechanism underlying the FRG signature.
Patient cohorts and data sources
Two patient cohorts from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/) and Clinical Proteomic Tumor Analysis Consortium (CPTAC)  were recruited in this study. Only patients with complete gene expression data and clinical prognosis information were selected for analysis. The TCGA cohort included 531 ccRCC patients with pre-processed RNA-sequencing data and corresponding clinical data, which were retrieved from TCGA database. Ninety-eight patients in the CPTAC cohort with complete clinical data and processed RNA-seq sequencing data were also download from CPTAC database, including RNA-seq data of matched normal renal tissues. Basic clinicopathologic features were shown in Table 1. In addition, twenty-two ferroptosis-related genes were retrieved from the previous study .
Development and verification of the prognosis model based on FRGs
Based on 22 FRGs in the TCGA cohort, we carried out least absolute shrinkage and selection operator (LASSO)-cox regression analysis by using glmnet package in R environment, which had been widely used in prognostic studies related to survival [12,13,14,15], to screen out survival-related FRGs. Their respective coefficients were also calculated corresponding to their weights in the prognosis model. We set the numbers of lambda as 1000 to ensure the robustness in the LASSO analysis. The FRG score was calculated as follow:
Coefi represents the coefficient of each prognosis-associated FRG, while Feri refers to the corresponding mRNA expression. We further verified the prognosis model in the independent CAPAC cohort, with cut-off value of the median value for each cohort.
Constructing and evaluating a predictive nomogram combining FRG score and clinicopathologic factors
In order to construct a predictive nomogram for patients with ccRCC, we combined FRG score and clinicopathologic factors via nomogramEx and rms packages. Calibration with bootstrapping, receiver operating characteristic (ROC) curves with area under curve (AUC) value and decision curve were further carried out to evaluate the performance of the nomogram.
Functional enrichment analysis
To explore the potential mechanism of the FRG score, we carried out weighted gene co-expression network analysis (WGCNA) to develop co-expression gene networks based on valid differentially expressed genes (DEGs, |fold change| ≥ 2, p < 0.01) in ccRCC, which were summarized in GEPIA2 . Correlations between gene modules and predicted clinicopathologic features were calculated to find out the optimum module that was significantly associated with the FRG score. Finally, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Genetic Ontology (GO) analysis were carried out to explore the potential biological mechanisms in which the FRG signature might be involved via Metascape .
In this study, R (3.6.2) were used for statistical analysis. Comparation of continuous variable between two groups was carried out through the Mann-Whitney U test. Comparation of continuous variable among more than two groups was carried out through the analysis of variance (ANOVA). Kaplan–Meier curve analysis was performed by using log-rank test to compare overall survival (OS). Pearson correlation analysis was also performed to calculate the correlation coefficient between two variables.
Developed and verified the prognosis model based on FRGs
We firstly identify the important roles of FRGs in ccRCC from the CAPAC cohort. As shown in Fig. 1A, 7 FRGs (HSPB1, FANCD2, TFRC, RPL8, CARS, CDKN1A and SLC7A11) were significantly up-regulated in tumor samples, while 8 FRGs, including MT1G, CISD1, FDFT1, SLC1A5, GLS2, ATP5G3, ACSL4 and HSPA5, were found to be significantly down-regulated in ccRCC (Fig. 1B).
Based on the TCGA cohort, LASSO-cox regression analysis screened out 11 prognosis-associated FRGs for the construction of the prognosis model, including CARS, FANCD2, ACSL4, CISD1, SLC1A5, SLC7A11, MT1G, CDKN1A, FDFT1, GLS2 and NCOA4 (Fig. 2a, b). The selected genes and their respective coefficients were shown in Table 2. The FRG score was calculated as mentioned in the method part. Significantly differential OS (hazard ratio = 3.61, 95% CI: 2.68–4.87, p < 0.0001) was observed between patients with high and low FRG score in the TCGA cohort (Fig. 2c), which was further verified in the CPTAC cohort with hazard ratio value of 5.13 (95% CI: 1.65–15.90, p = 0.019, Fig. 2d). As shown in Fig. 2e, the heatmap and correlation analysis indicated that higher FRG score was significantly associated with higher rate of death (r = 0.42, p < 0.0001), higher tumor stage (r = 0.41, p < 0.0001), and higher tumor stage (r = 0.45, p < 0.0001).
The FRG score could act as a prognostic factor for patients with ccRCC
We next to perform univariate cox regression analysis to further evaluate the prediction performance of the FRG score. As shown in Fig. 3a, the FRG score could act as a risk factor for survival prediction of ccRCC patient in the TCGA cohort, which was further verified in the CPTAC cohort (Fig. 3b). Further comparative analyses revealed that the patients with high tumor grade (P < 0.0001, Fig. 3c) or high tumor stage (P < 0.0001, Fig. 3d) tended to have higher level of FRG score. In addition, higher level of FRG score was also found in patients with tumor lymph node metastasis status (P < 0.0001, Fig. 3e) or tumor distant metastasis status (P < 0.0001, Fig. 3f).
Subgroup survival analysis of the prognosis model in the TCGA cohort revealed that our FRG score could significantly distinguish patients with high survival risk among different tumor stages (Fig. 4a) and different tumor grades (Fig. 4b). New tumor staging system based on current staging system and the FRG score performed well in distinguishing ccRCC patients with different clinical prognoses (Fig. 4c). Patients with stage i-ii/high risk score tumors had similar survival outcomes compared to patients with stage iii-iv/low risk score tumors (p = 0.0824). In addition, patients with grade 1–2/high risk score tumors also had similar survival outcomes when compared to patients with grade 3–4/low risk score tumors (p = 0.9163, Fig. 4d).
Improved prognostic accuracy of the FRG score integrated with clinicopathologic features
To explore whether the accuracy of the prognosis model could be improved through combining our FRG score and clinicopathologic features, we developed an integrated nomogram based on the FRG score, patient age, tumor grade and tumor stage (Fig. 5a). The calibration analysis indicated that the survival rate predicted by the nomogram had excellent agreement with actual observations at 1-, 3- and 5-year follow up (Fig. 5b). Further decision curve analysis verified the improved prognostic accuracy via the FRG score (Fig. 5c). ROC curve analyses illustrated that that AUC of the nomogram for survival prediction in 3- and 5-year reached to 84.5 and 83.2%, respectively (Fig. 5d, e).
Cell cycle-related pathways were associated with the FRG score in ccRCC
A total of 2353 DEGs were analyzed through WGCNA in this study. According to the recommendation of pickSoftThreshold, the soft-thresholding power of β value was set as 18 (Fig. 6a). All the DEGs associated with ccRCC were then hierarchically clustered into 5 gene modules (Fig. 6b). As shown in Fig. 6c, correlation analysis indicated that the green model (MEgreen) seemed to have the highest correlation with FRG score. The heatmap also illustrated the relationships of the 43 classified genes of the green model and the FRG score (Fig. 6d). Further functional enrichment analysis revealed that our FRG score might involve in pathways associated with the process of cell cycle, including cell cycle-mitotic pathway, cytokinesis pathway and nuclear division pathway (Fig. 6e).
Due to the significant heterogeneity and aggressiveness of ccRCC, different clinical outcomes could still be found in patients with similar tumor stage or grade. Therefore, it is of great urgency to find out novel prognostic markers for clinical practices. Fortunately, high-throughput genetic techniques for oncology have revolutionized the development of prognostic biomarkers for malignant tumor [10, 18].
Ferroptosis is closely related to tumor invasion and metastasis. It was reported that tumor-infiltrating CD8+ T cells with CD36 deficiency had low expression of FRGs, and CD36 deficiency had been confirmed to be associated with reduced ferroptosis in tumor-infiltrating CD8+ T cells . A recent study has also revealed that oleic acid could protect melanoma cells from ferroptosis through acsl3-dependent manner. In addition, melanoma cells from lymph nodes were more resistant to ferroptosis .
Here, we developed and verified a prognostic model based on FRGs from two independent patient cohorts. Cox regression analysis revealed that the FRG score could act as a survival risk factor for patients with ccRCC. Improved prognostic accuracy was also found in the nomogram integrated with FRG score and clinicopathologic features, which had excellent agreements with the actual survival rates at 1-, 3- and 5-year follow up.
The new tumor staging system based on current staging system and the ferroptosis-related prognosis model was found to act better in distinguishing ccRCC patients with worse prognosis. Patients with low tumor stage (stage i-ii) might be faced with similar survival risk to patients with high tumor stage (stage iii-iv) if they were accompanied with high FRG score, which might account for the clinical observations that different survival outcomes could sometime be found in patients with similar tumor staging, indicating the potential application value of our ferroptosis-based prognosis model in clinical practices.
A total of 43 core genes clustered in the green model were found to be significantly associated with our FRG score. Further functional enrichment analysis revealed that our FRG score might involve in pathways associated with the process of cell cycle in ccRCC. It was reported that ferroptosis could be regulated by p53, which was an indispensable regulator of the cell cycle and could enhance ferroptosis by inhibiting SLC7A11. In addition, p53 could also inhibit ferroptosis by directly inhibiting the activity of DPP4 or inducing the expression of CDKN1A/p21 .
Several limitations could still be found in this study. Firstly, cross-validations among two independent patient cohorts were carried out in this study, however, potential bias might still exist since retrospective public cohorts were used for analyses. Secondly, prospective single- or multi-center studies are still wanted for further verifying the ferroptosis-related prognostic model. Finally, even though our study revealed that cell cycle-related pathways were associated with the FRG score in ccRCC, experimental studies for potential mechanism exploring and function verification are still needed for subsequent analyses.
We developed and verified a FRG signature for the prognosis prediction of patients with ccRCC, which could act as a risk factor and help to update the tumor staging system when integrated with clinicopathological characteristics. Cell cycle-related pathways might be involved in the regulation of ccRCC through ferroptosis, which still need further experimental studies for function verifications of the study.
Clear-cell renal cell carcinoma
The Cancer Genome Atlas
Clinical Proteomic Tumor Analysis Consortium
Least absolute shrinkage and selection operator
Receiver operating characteristic
Area under curve
Weighted gene co-expression network analysis
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71(1):7–33.
Zhang S, Sun K, Zheng R, Zeng H, He J. Cancer incidence and mortality in China, 2015. J Natl Cancer Center. 2020;1(1):2–11.
Nabi S, Kessler ER, Bernard B, Flaig TW, Lam ET. Renal cell carcinoma: a review of biology and pathophysiology. F1000Res. 2018;7:307.
Reuter VE. The pathology of renal epithelial neoplasms. Semin Oncol. 2006;33(5):534–43.
De P, Otterstatter MC, Semenciw R, Ellison LF, Marrett LD, et al. Trends in incidence, mortality, and survival for kidney cancer in Canada, 1986-2007. Cancer Causes Control. 2014;25(10):1271–81.
Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149(5):1060–72.
Wang J, Yin X, He W, Xue W, Zhang J, Huang Y. SUV39H1 deficiency suppresses clear cell renal cell carcinoma growth by inducing ferroptosis. Acta Pharm Sin B. 2021;11(2):406–19.
Mou Y, Wu J, Zhang Y, Abdihamid O, Duan C, Li B. Low expression of ferritinophagy-related NCOA4 gene in relation to unfavorable outcome and defective immune cells infiltration in clear cell renal carcinoma. BMC Cancer. 2021;21(1):18.
Yang WH, Ding CC, Sun T, Rupprecht G, Lin CC, Hsu D, et al. The hippo pathway effector TAZ regulates Ferroptosis in renal cell carcinoma. Cell Rep. 2019;28(10):2501–8.
Clark DJ, Dhanasekaran SM, Petralia F, Pan J, Song X, Hu Y, et al. Integrated Proteogenomic characterization of clear cell renal cell carcinoma. Cell. 2019;179(4):964–83.
Liu Z, Zhao Q, Zuo ZX, Yuan SQ, Yu K, Zhang Q, et al. Systematic analysis of the aberrances and functional implications of Ferroptosis in Cancer. iScience. 2020;23(7):101302.
Hirata T, Arai Y, Yuasa S, Abe Y, Takayama M, Sasaki T, et al. Associations of cardiovascular biomarkers and plasma albumin with exceptional survival to the highest ages. Nat Commun. 2020;11(1):3820.
Cao R, Yuan L, Ma B, Wang G, Qiu W, Tian Y. An EMT-related gene signature for the prognosis of human bladder cancer. J Cell Mol Med. 2020;24(1):605–17.
Liu Y, Zhang X, Zhang J, Tan J, Li J, Song Z. Development and validation of a combined Ferroptosis and immune prognostic classifier for hepatocellular carcinoma. Front Cell Dev Biol. 2020;8:596679.
Liang JY, Wang DS, Lin HC, Chen XX, Yang H, Zheng Y, et al. A novel Ferroptosis-related gene signature for overall survival prediction in patients with hepatocellular carcinoma. Int J Biol Sci. 2020;16(13):2430–41.
Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–60.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.
Wang LB, Karpova A, Gritsenko MA, Kyle JE, Cao S, Li Y, et al. Proteogenomic and metabolomic characterization of human glioblastoma. Cancer Cell. 2021;39(4):509–28.
Ma X, Xiao L, Liu L, Ye L, Su P, Bi E, et al. CD36-mediated ferroptosis dampens intratumoral CD8 + T cell effector function and impairs their antitumor ability. Cell Metab. 2021;33(5):1001-12.
Ubellacker JM, Tasdogan A, Ramesh V, Shen B, Mitchell EC, Martin-Sandoval MS, et al. Lymph protects metastasizing melanoma cells from ferroptosis. Nature. 2020;585(7823):113–8.
Kang R, Kroemer G, Tang D. The tumor suppressor protein p53 and the ferroptosis network. Free Radic Biol Med. 2019;133:162–8.
We are grateful to the TCGA and the CPTAC project teams.
This work was supported by the National Natural Science Foundation of China (81972393 and 82002665).
Ethics approval and consent to participate
Consent for publication
Authors confirmed that this work can be published. The content of this manuscript is original that has not yet been accepted or published elsewhere.
The author declares that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, S., Zhang, E., Guo, T. et al. A novel ferroptosis-related gene signature associated with cell cycle for prognosis prediction in patients with clear cell renal cell carcinoma. BMC Cancer 22, 1 (2022). https://doi.org/10.1186/s12885-021-09033-7
- Clear cell renal cell carcinoma
- Cell cycle