Skip to main content
Download PDF
  • Primary research
  • Open access
  • Published:

Construction and validation of a ferroptosis-related long noncoding RNA signature in clear cell renal cell carcinoma

Cancer Cell International volume 22, Article number: 283 (2022) Cite this article

Abstract

Background

Clear cell renal cell carcinoma (ccRCC) is characterized by the accumulation of lipid-reactive oxygen species. Ferroptosis, due to the lipid peroxidation, has been reported to be strongly correlated with tumorigenesis and progression. However, the functions of the ferroptosis process in ccRCC remain unclear.

Methods

After sample cleaning, data integration, and batch effect removal, we used the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases to screen out the expression and prognostic value of ferroptosis-related lncRNAs and then performed the molecular subtyping using the K-means method. Then, the functional pathway enrichment and immune microenvironment infiltration between the different clusters were carried out. The results showed a significant difference in immune cell infiltration between the two clusters and the associated marker responded to individualized differences in treatment. Then, least absolute shrinkage and selection operator (LASSO) Cox regression was used to establish a prognostic signature based on 5 lncRNAs. This signature could accurately predicted patient prognosis and served as an independent clinical risk factor. We then combined significant clinical parameters in multivariate Cox regression and the prognostic signature to construct a clinical predictive nomogram, which provides appropriate guidance for predicting the overall survival of ccRCC patients.

Results

The prognostic differentially expressed ferroptosis-related LncRNAs (DEFRlncRNAs) were found, and 5 lncRNAs were finally used to establish the prognostic signature in the TCGA cohort, with subsequently validation in the internal and external cohorts. Moreover, we conducted the molecular subtyping and divided the patients in the TCGA cohort into two clusters showing differences in Hallmark pathways, immune infiltration, immune target expression, and drug therapies. Differences between clusters contributed to individualizing treatment. Furthermore, a nomogram was established to better predict the clinical outcomes of the ccRCC patients.

Conclusions

Our study conducted molecular subtyping and established a novel predictive signature based on the ferroptosis-related lncRNAs, which contributed to the prognostic prediction and individualizing treatment of ccRCC patients.

Introduction

Renal cell carcinoma (RCC) is one of the most common malignant tumors of the urinary system [1]. As the major pathologic subtype of RCC, clear cell renal cell carcinoma (ccRCC), especially metastatic ccRCC, often has high morbidity and mortality [2]. Approximately 25–30% of ccRCC patients have metastasis at initial diagnosis, which usually indicates poor prognosis. Moreover, the tumor-node-metastasis (TNM) staging system, which is currently applied in clinical practice, is considered less accurate in evaluating the prognosis and progression of ccRCC patients [3]. Meanwhile, 73–75% of identified ccRCC driver aberrations were subclonal, which might contribute to different clinical outcomes [4, 5]. Hence, performing suitable molecular subtyping and exploring new prognostic signatures to diagnose and evaluate the prognosis of ccRCC patients remains significant.

During the past decades, ferroptosis has been gradually identified as an iron-dependent, nonapoptotic cell death mode characterized by the accumulation of lipid reactive oxygen species [6, 7]. An increasing number of studies have shown that dysregulation of ferroptosis-related genes (FRGs) plays an important role in the occurrence and development of many diseases, especially cancer [8,9,10]. Studies related to ferroptosis have become the focus in treating and detecting related diseases. Moreover, as the noncoding RNAs play an emerging role in cancer, long noncoding RNAs (lncRNAs) have been fully studied [11]. Interestingly, an increasing number of studies have demonstrated that lncRNAs can participate in the ferroptosis process and then influence tumor development [12,13,14]. Therefore, exploring ferroptosis-associated lncRNAs might provide new ideas and insights for ccRCC treatment and prediction.

In this study, we first explored the potential biological functions and correlation of FRGs. Afterward, we screened the prognostic ferroptosis-related lncRNAs (FRlncRNAs) and conducted the molecular subtyping of ccRCC patients. The correlation between the molecular clusters and immune cell infiltration was explored. Afterward, we conducted the LASSO Cox regression to establish a 5-ferroptosis-related LncRNAs signature in ccRCC patients and validated it in the ICGC database. Furthermore, a nomogram was constructed integrating the prognostic signature and significant clinical parameters. The results showed that a good predictive performance for the overall survival (OS) of ccRCC patients.

Materials and methods

Data selection and processing

The ccRCC sequencing data (HTSeq-FPKM) and the latest corresponding clinical information (Additional file 1: Table S1) were downloaded from the TCGA database (https://cancergenome.nih.gov/), including 539 ccRCC samples and 72 normal controls. Moreover, we downloaded the transcriptome profile and corresponding survival data from the ICGC database (http://daco.icgc.org/) as the validation cohort, including 91 RCC samples. Then, we distinguished between lncRNAs and mRNAs using the human GTF annotation file. The Matrix processing and batch effect removal were conducted via the limma and sva packages in R (v 4.0.3).

Cell lines and clinical specimens

Seven kidney and RCC cell lines (kidney cell lines: HK-2; RCC cell lines: 786-O, 769-P, OS-RC-2, A498, ACHN, and Caki-1) were obtained from the American Type Culture Collection (ATCC). All cells were maintained in RPMI-1640 (Corning, USA) or high-glucose DMEM medium (Gibco, USA) with 10% fetal bovine serum (BI, Israel) and 1% penicillin and streptomycin (Gibco, USA) at 37 °C and 5% CO2. Ten paired ccRCC and adjacent normal tissues were obtained from ccRCC patients undergoing surgical resection at Peking University First Hospital. Detailed information on the 10 paired tissue specimens is shown in Additional file 2: Table S2. The Ethics Committee approved this study of PUFH, and all patients signed informed consent forms. All procedures were performed according to the World Medical Association Declaration of Helsinki.

Real-time quantitative PCR (qPCR)

Total RNA of 10 paired clinical samples and 7 cell lines were extracted by Takara kit according to the manufacturer's protocol. Then, the RNA was reverse transcribed to cDNA in a 20 µl reaction system. All gene transcripts were quantified by qPCR using SYBR Premix ExTaq kit, and TUBA was used as a normalization control. The primer sequences are listed in the Additional file 3: Table S3. Each reaction was performed four times, and the 2^−△△CT method was used to calculate the relative mRNA expression level.

Identification of prognostic ferroptosis-related differentially expressed lncRNAs

According to previous studies [15,16,17], we obtained 259 FRGs, and the list is shown in Additional file 4: Table S4. After that, we screened ferroptosis-related lncRNAs with a filter (correlation  > 0.5, p < 0.01), including 2854 lncRNAs. The limma package was used to determine the differentially expressed lncRNAs (DElncRNAs) between the ccRCC patients and normal controls, including 1333 lncRNAs [18]. We intersected the DELncRNAs and ferroptosis-related LncRNAs to obtain the ferroptosis-related DElncRNAs (FRDElncRNAs). Univariate Cox regression of OS was performed on the FRDElncRNAs. Those with p < 0.01 were considered prognostic FRDElncRNAs.

Molecular subtyping in ccRCC based on prognostic FRDElncRNAs

After obtaining prognostic FRDELncRNAs, we performed consensus clustering to identify the molecular subtypes of ccRCC by using the ConsensusClusterplus R package [19]. We selected 80% of the prognostic FRDElncRNAs resampling 100 times and determined clusterings of specified cluster counts (k). After this, the pairwise consensus values were calculated and stored in a symmetrical consensus matrix for each k. The k, at which there was no appreciable increase, was determined by the cumulative distribution function (CDF) plot and delta area plot. The alteration in immune infiltration between different clusters was estimated using the CIBERSORT method (Additional file 5: Figure S1).

Potential biological functional enrichment

To gain insights into the cellular functions directly regulated by FRG transcriptional control, we compared the list of FRGs to the biological pathways annotated by the Kyoto Encyclopedia of Genes and Genomes (KEGG) [20]. Afterward, according to the two clusters, Gene Set Variation Analysis (GSVA) was conducted using the GSVA package in R software v.4.0.3 to investigate the enrichment of HALLMARK pathways with the h.all.v7.4.symbols.gmt gene set from the Molecular Signature Database [21].

Analysis of potential therapeutic targets in different clusters

According to the two divided clusters, we determined whether there were differences in treatment effects between groups based on relevant databases. Since targeted drugs are commonly used to treat advanced kidney cancer, we used the R pRRophetic package to estimate drug response as determined by the half-maximal inhibitory concentration (IC50) for each kidney cancer patient on the Genomics of Drug Sensitivity in Cancer (GDSC) website [22]. The drug predictive model was applied to the processed, standardized, and filtered clinical tumor expression data, and then it estimated the drug sensitivity for each patients. Furthermore, based on the Cancer Immunome Altas (TCIA) database (http://tcia.at/), the Immunophenoscore (IPS) was obtained [23, 24]. To predict sensitivity to immunotherapy between different clusters, we compared the IPS of the two clusters in different immunotherapy decisions.

Construction of the prognostic predictive risk signature

First, the TCGA cohort patients were randomly divided a the training set and an internal validation set. Meanwhile, the patients in the ICGC cohort were used as the external validation cohort. Based on the prognostic FRDELncRNAs, we performed LASSO Cox regression using the glmnet R package. We calculated each patient's risk score using the regression coefficient score of the individual lncRNAs and their expression values. We defined the formula for calculating the prognostic risk score as follows: Risk score = coef(Lnc1)*Exp(Lnc1) + coef(Lnc2)* Exp (Lnc2) + … + coef(Lncn)* Exp (Lncn), where “coef” represents the coefficient score estimated by LASSO Cox regression, and “Exp” represents the expression value of the individual LncRNAs. The detailed information of the signature is shown (Table 1). Then, we classified the ccRCC patients into the high- and low-risk groups, according to the median risk score of the training group as the cutoff [25].

Table 1 Detailed information of the LncRNAs in the prognostic signature
Full size table

Validation of the prognostic risk signature

We conducted the Kaplan–Meier and receiver operating characteristic (ROC) curve analyses to assess the prognostic risk signature’s validity. According to the calculated median risk score, all samples in each group were divided into high- and low-risk groups, and the survminer and timeROC packages were used to validate the predictive accuracy in the training and validation sets. The area under the curve (AUC) values corresponding to 1-, 3-, and 5-years were calculated. The time-dependent ROC curve was used to validate the predictive performance of the signature. An AUC value of 0.75 or higher was considered the significant predictive value, and the value of 0.60 or higher was regarded as acceptable for prediction. Furthermore, univariate and multivariate Cox regression was conducted to explore whether the ferroptosis-signature (FerroSig) could serve as an independent factor.

Construction and validation of the nomogram

To better predict the prognosis of patients with ccRCC, we established a predictive nomogram based on clinical parameters and prognostic signature [26]. In brief, we first performed univariate and multivariate Cox regression analyses to identify clinical parameters and risk scores that could be used as independent risk factors. Subsequently, the significant factors were used to construct the predictive nomogram. We then evaluated the nomogram effect using calibration curves and time-dependent ROC curves. An AUC value of 0.75 or higher was considered a significant predictive value, and a value of 0.60 or higher was regarded as acceptable for prediction.

Results

Identification of the prognostic FRDElncRNAs

The data processing was performed as described in the methods above. The flow chart of the whole process is shown (Fig. 1). To explore the potential functions of the FRGs, we first conducted KEGG pathway enrichment analysis. The results showed that FRGs were mainly enriched in ferroptosis and autophagy pathways (Fig. 2A). Afterward, we explored the FRlncRNAs with a correlation  > 0.5 and p < 0.01, using the limma package. Next, 1333 DElncRNAs between the ccRCC and normal samples from the TCGA set were screened, and the 5 most obvious upregulated- and downregulated lncRNAs were identified (Fig. 2B). Theoverlapping lncRNAs in both DElncRNAs and FRlncRNAs were identified as the FRDElncRNAs. There were 723 lncRNAs for OS selected as prognostic FRDELncRNAs for subsequent analyses (Fig. 2C).

Fig. 1
figure 1

Flow chart of the whole analysis processes of this study

Full size image
Fig. 2
figure 2

Function annotations of ferroptosis-related genes and Identification of the ferroptosis-related LncRNAs. A KEGG enrichment analysis of the ferroptosis-related genes, the larger the shape of the dot, the more corresponding genes are represented. B Volcano plot of the differentially expressed ferroptosis-related LncRNAs between the ccRCC and normal samples. The 5 most significant up-and down-regulated LncRNAs were labeled separately. C The Venn plot of the overlapped LncRNAs between two cohorts

Full size image

Molecular subtyping showed differences in therapeutic choice and immune microenvironments

The prognostic FRDElncRNAs above were used to screen the molecular subtypes of ccRCC by using the ConsensusClusterPlus packages. The K-means method was performed for clustering, and 80% of the lncRNAs were sampled 100 times using the resampling method. The consensus CDF and delta area were calculated to determine the clustering outcomes, as shown in Fig. 3B. When the cluster number was 2, there was no significant increase in the area under the CDF curve. Hence, we finally divided the samples in the TCGA cohort into cluster 1 and cluster 2. The representative consensus matrix of two clusters (Fig. 3A) displayed a well-defined 2-block structure. Principal component analysis (PCA) showed that the samples from two clusters could be well separated (Fig. 3C). Then, a heatmap integrating the expression of the prognostic FRDElncRNAs and clinical parameters in each subtype was generated, as is shown in Fig. 3D. The results showed that pM stage was higher in cluster1 than in cluster 2, which also clarified that the patients in cluster 1 had a worse prognosis.

Fig. 3
figure 3

Molecular subtyping of the ccRCC patients based on the prognostic ferroptosis-related LncRNAs. A The heatmap corresponding to the consensus matrix for k = 2 was obtained by applying consensus clustering. Color gradients represent consensus values from 0–1; white corresponds to 0 and dark blue to 1. B Consensus among clusters for each category number k. C Principal Component Analysis and D K-M survival analysis of the two clusters. E The composite heatmap corresponding to the cluster and mRNA expression, TNM stage, AJCC stage, and ISUP grade, and Age as the annotations. **P < 0.01

Full size image

Comprehensive analyses of the molecular clusters

Subsequently, we explored HALLMARK pathway enrichments alterations between the two clusters using the GSVA method. The results showed that several oncogenic pathways were significantly altered between the two clusters, such as hypoxia and apoptosis (Fig. 4A). Furthermore, using the pRRophetic packages, we calculated and analyzed the IC50 between the two clusters for clinical drugs used for advanced ccRCC patients. According to the median IC50 value, the outcome indicated that the patients in cluster 2 were more sensitive to Daunorubicin and Tipifarnib. In contrast, using the Dasatinib, Paclitaxel, Sorafenib, and Pazopanib would be more effective for patients in cluster 1 (Fig. 4B).

Fig. 4
figure 4

Biological functions and immune infiltration between two clusters. A Gene set variation analysis (GSVA) was performed to compute HALLMARK pathways between two clusters. B The IC50 data on drugs for ccRCC differentially expressed in the two clusters obtained by applying pRRophetic were shown

Full size image

Since the tumor microenvironment (TME) plays a vital role in the development of tumors, we further explored the correlation between the TME and molecular clusters. As shown in Fig. 5A, we first calculated the stromal score, immune score, and ESTIMATE score. Only the immune score showed a significant difference between the two clusters. Hence, we then explored the immune cell infiltration between two clusters using the CIBERSORT method. Cluster 1 exhibited a higher infiltration of CD8 T cells, regulatory CD4 T cells, and Neutrophils. In comparison, cluster 2 showed a higher infiltration of macrophages (Fig. 5B). Meanwhile, the expression levels of some potential immunotherapy targets changed significantly in both classifications. The mRNA expression levels of potential immune therapeutic targets, such as BRAF and PD-1, were markedly higher in cluster 1 (Fig. 5C). Next, we explored the TCIA database for differences in the presence of immunotherapeutic targets. Previous studies reported the role of IPS in predicting the immunotherapy response of ccRCC patients. According to the IPS, the outcome showed that patients in Cluster 1 obtained a higher IPS of CTLA4 and CTLA4 + PD-1 than Cluster 2 (Fig. 5D). These results suggest that the cluster 1 group may be more sensitive to immunotherapy.

Fig. 5
figure 5

Tumor environment especially immune infiltration and potential immune therapeutic target between two clusters. A Tumor environment scores between two clusters. B The differentially expressed immune infiltrated cells in the two clusters obtained by applying Cibersoft were shown. C Common potential immune therapeutic targets between two clusters were shown. *P < 0.05, **P < 0.01, *** P < 0.001 D The different expressed of four immune status targets between two clusters were explored from the TCIA database

Full size image

Construction and validation of the prognostic signatures based on the prognostic FRDELncRNAs

Based on the prognostic FRDElncRNAs obtained from the univariate Cox regression analysis, we constructed a 5 FRDElLncRNAs based prognostic signature using the LASSO regression (Fig. 6A and B). Detailed information on the lncRNAs from the signature is shown in Table 1. The risk score was constructed according to the following formula: risk score = 0.058 (LINC00460) + 0.088 (LINC00894) + 0.118 (VPS9D1-AS1) + 0.013 (CYTOR) + 0.075 (FOXD2-AS1).  Then, we explored the prognosis of TCGA-KIRC patients affected by the expression levels of the five LncRNAs, including overall survival and disease-free survival (Additional file 5: Fig. S1). Then, based on the calculated median risk score cutoff, patients were divided into the high- and low-risk groups. The risk score distribution, survival status, and expression of ncRNAs from the signatures are exhibited in the TCGA training cohort, TCGA validation set, and ICGC validation set are exhibited in Fig. 6C–E. The Kaplan–Meier log-rank test and the time-dependent ROC curve were used to evaluate the predictive ability and accuracy of the prognostic signature. The results of the Kaplan–Meier log-rank test showed that the high-risk group had a significantly worse OS than the low-risk group in the TCGA training set (Fig. 7A), TCGA validation set (Fig. 7B), and ICGC validation set (Fig. 7C). Moreover, the time-dependent ROC curve proved the 1-year, 3-year, and 5-year predictive accuracy of the signature for OS (Fig. 7A–C).

Fig. 6
figure 6

Establishment of the 5-LncRNAs based prognostic signature. A LASSO coefficient profiles of the prognostic DEFRGs. B Partial likelihood deviance was plotted versus log (Lambda). The vertical dotted line indicates the lambda value with the minimum error and the largest lambda value. (C-E) LncRNA expression patterns and the distribution of survival status increased risk score in the TCGA training set, TCGA internal validation set, and ICGC external validation set

Full size image
Fig. 7
figure 7

Validation of the prognostic signature. K-M plot analyses between the high-and low-risk patients in the TCGA training cohort A, TCGA validation cohort B, and ICGC validation cohort C. The 1-, 3-, 5-year time-dependent ROC curves in the TCGA training cohort (A), TCGA validation cohort (B), and ICGC validation cohort (C)

Full size image

Construction and validation of the prognostic nomogram

After establishing and validating the signatures based on the prognostic FRDElncRNAs, univariate and multivariate Cox regression analyses were used to explore independent risk factors in the TCGA dataset. As shown in Table 2, univariate Cox regression analysis showed that AJCC stage, ISUP grade, age, and the risk signature were significantly correlated with OS (Table 2). Moreover, the multivariate Cox regression analyses of the clinical parameters above showed that AJCC stage, ISUP grade, age, and the risk signature, were significantly correlated with OS (Table 2).

Table 2 Univariate and multivariate Cox analyses of clinical parameters and risk signature
Full size table

To better assess patient prognosis and guide clinical decision-making, we established a nomogram integrating the risk signature and significant clinical parameters in the multivariate Cox regression analyses (Fig. 8A). The C-index of 0.773 showed a good agreement, and the established nomogram was shown (Fig. 8B). The calibration curves showed that the predictive nomogram could well predict the survival status of patients at 1, 3, and 5 years (Fig. 8C). Additionally, the nomogram showed better predictive value than clinical indicators and the FRlncRNAs signatures (Fig. 8D).

Fig. 8
figure 8

Construction and validation of the prognostic nomogram. A Forest plot of the significant clinical parameters in the multivariate Cox regression. B The nomogram based on the significant clinical parameters and risk signature C Calibration curves of the nomogram for 1-, 3-, and 5-year survival prediction. D The predictive value of the nomogram, risk signature, and clinical parameters

Full size image

Validation of the lncRNAs in cell lines and clinical specimens

We finally validated the 5 lncRNAs in signature in ccRCC clinical samples and RCC cell lines. Compared to the normal cell line HK-2, the expression levels of LncRNAs in RCC cell lines were inconsistent, with higher or lower levels present (Fig. 9A). The results may be due to the inability of a single cell line to mimic the overall situation of RCC and adjacent normal tissues. We then examined the expression levels of the corresponding mRNAs in 10 paired clinical samples, in which LNC000894, VPS9D-AS1, and CYTOR were highly expressed in cancer, while LNC000460 and FOXD2-AS1 were not evident (Fig. 9B).

Fig. 9
figure 9

Validation of expression levels of LncRNAs in signaturein cell lines and samples. A Levels of mRNA expression of LncRNAs in cell lines. B The mRNA expression levels of LncRNAs in 10 pairs of paired clinical samples. *P < 0.05

Full size image

Discussion

The ccRCC is a molecularly heterogeneous tumor characterized as radiotherapy and chemotherapy-resistant [2, 5]. With the development of diagnostic and therapeutic techniques, the 5-year survival rate of ccRCC patients has been significantly improved. However, 25–30% of ccRCC patients have metastases at initial diagnosis [4]. The 5-year survival rate is merely 10%. Most ccRCC patients have no apparent symptoms such as pain and hematuria, usually resulting in diagnostic difficulty in the early stage [27]. Moreover, the current TNM staging system used in clinical practice lacks accuracy for prognostic evaluation [28, 29]. The early diagnosis and accurate assessment of ccRCC patients remain challenging for the reasons mentioned above. Therefore, it is still vital to screen out new clinical and molecular biomarkers for diagnosis and treatment.

Ferroptosis plays essential roles in the progression and tumorigenesis of RCC [30]. Notably, the expression levels of various ferroptosis-related genes were significantly correlated with the prognosis of ccRCC patient, suggesting that targeting ferroptosis-related pathway might be an effective option for ccRCC treatment [31, 32]. Moreover, lncRNAs plays important roles in regulating the expression of FRGs and the process of ferroptosis [33, 34]. Hence, in this study, we comprehensively analyzed the expression and prognosis of ferroptosis-related lncRNAs. And then, we performed molecularly typing and developed a prognostic model based on 5 lncRNAs in patients with ccRCC, using the method above.

Although many molecular subtypes of ccRCC based on gene expression have been proposed in recent years, the cluster of lncRNAs associated with ferroptosis has not been fully explored [35]. Therefore, we divided ccRCC patients into two clusters based on prognostic DEFRlncRNAs using the NMF algorithm. PCA analysis showed significant differences between the two clusters. The Kaplan–Meier plot showed that cluster 1 had a worse prognosis than cluster2. GSVA analysis showed that many tumor-related pathways were significantly altered between the two clusters (such as hypoxia and coagulation), suggesting that patients in the two clusters may have different sensitivities to some clinical drugs. We then tested the sensitivity of different cluster of patients to some drugs commonly used in advanced kidney cancer. The results showed promising differences, which contributed to the implementation of individualized treatment.

The immune components of the tumor microenvironment and immune cells significantly regulate tumor development [36]. We further compared the differences in immune cell infiltration between the two clusters. We found high levels of the CD8 + T cells, T cell regulatory, follicular helper T cells, and memory B cells were presented high expression levels in cluster 1. In contrast, the expression levels of neutrophils and macrophages were significantly increased in cluster 2. Unlike most cancer types, previous studies have shown that the density of CD8 T cells is correlated with poor prognosis in patients with RCC [37, 38]. This also explains why patients in cluster 1 have a worse prognosis. Moreover, the infiltration of mesenchymal cells and neutrophils may serve as a protective factor for RCC. We then examined at the mRNA expression levels of common therapeutic targets, and patients in cluster 1 had significantly upregulated expression levels of PD1, CTLA4, and other genes. We also examined IPS under different immunotherapy modalities through the TCIA database. We found higher CTLA4 and CTLA4 + PD1 IPS in cluster1, suggesting that patients in cluster1 may be more sensitive to immunotherapy. However, cluster 1 patients showed a worse prognosis in previous results. This suggests that the sensitivity of immunotherapy may not have a significant impact on prognosis.

Afterward, we established a predictive signature using LASSO regression. The Kaplan–Meier plot and time-dependent ROC curves showed that the predictive signature exhibited good predictive performance. Furthermore, most of the lncRNAs in our signature have been reported in various cancer types. LINC00460 has been extensively studied in cancer and has been shown in several studies to be a prognostic target in renal clear cell carcinoma [39,40,41]. LINC00894 was reported to promote breast cancer metastasis by regulating ZEB1 [42]. CYTOR and VPS9D1-AS1 are associated with the prognosis of multiple cancers and can regulate the progression of multiple cancers by sponging miRNAs [43,44,45]. Similarly, FOXD2-AS1 was strongly associated with the prognosis and progression of cancer patients in various cancers [46, 47]. These validation results in multiple datasets and literature mining results indicate that the prognostic signature predicts the prognosis of ccRCC patients and may function as the regulator of ccRCC progression.

Nevertheless, this study has several certain limitations. First, our data are based on the TCGA and ICGC databases, and more independent datasets are needed for testing and validation. Second, some lncRNAs in our signature play an essential role in cancer and need to be validated in future experiments. In conclusion, we first systematically analyzed the expression and prognostic value of ferroptosis-related lncRNAs and assessed immune infiltration and potential prognostic targets by molecular subtyping in ccRCC patients.

Conclusion

In conclusion, our study identified the FRDELncRNAs and successfully constructed an individualized ccRCC signature (riskScore), which proved to be significantly correlated with OS in both the training and validation cohorts. We also estimated the potential relationships among immune cell infiltration, immunotherapy-related targets, and potential therapeutic drugs between the two molecular subtyping clusters. Our research is anticipated to provide new insights into ferroptosis-related lncRNAs for future work.

Availability of data and materials

Public databases analyzed in this study could be found here: TCGA (https://portal.gdc.cancer.gov/) and ICGC database (http://daco.icgc.org/). The authors will provide the original data supporting the conclusions of this paper without undue reservation.

Abbreviations

RCC:

Renal cell carcinoma

ccRCC:

Clear cell renal cell carcinoma

TNM:

Tumor node metastasis

LncRNAs:

Long non-coding RNAs

FRGs:

Ferroptosis-related genes

TCGA:

The Cancer Genome Atlas

GEO:

Gene expression omnibus

DELncRNAs:

Differentially expressed LncRNAs

FRDELncRNAs:

Ferroptosis-related differentially expressed LncRNAs

KEGG:

Kyoto encyclopedia of genes and genomes

GSVA:

Gene set variation analysis

GDSC:

Genomics of drug sensitivity in cancer

LASSO:

Least absolute shrinkage and selection operator

ROC:

Receiver operating characteristic

AUC:

Area under the curve

OS:

Overall survival

References

  1. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71(1):7–33.

    Article  PubMed  Google Scholar 

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

  3. Warren AY, Harrison D. WHO/ISUP classification, grading and pathological staging of renal cell carcinoma: standards and controversies. World J Urol. 2018;36(12):1913–26.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Turajlic S, Swanton C. Metastasis as an evolutionary process. Science. 2016;352(6282):169–75.

    Article  CAS  PubMed  Google Scholar 

  5. Sanfrancesco JM, Cheng L. Complexity of the genomic landscape of renal cell carcinoma: Implications for targeted therapy and precision immuno-oncology. Crit Rev Oncol Hematol. 2017;119:23–8.

    Article  PubMed  Google Scholar 

  6. Li J, Cao F, Yin H-L, Huang Z-J, Lin Z-T, Mao N, et al. Ferroptosis: past, present and future. Cell Death Dis. 2020;11(2):88.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Hirschhorn T, Stockwell BR. The development of the concept of ferroptosis. Free Radic Biol Med. 2019;133:130–43.

    Article  CAS  PubMed  Google Scholar 

  8. Mou Y, Wang J, Wu J, He D, Zhang C, Duan C, et al. Ferroptosis, a new form of cell death: opportunities and challenges in cancer. J Hematol Oncol. 2019;12(1):34.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Sun Y, Chen P, Zhai B, Zhang M, Xiang Y, Fang J, et al. The emerging role of ferroptosis in inflammation. Biomed Pharmacother. 2020;127:110108.

    Article  CAS  PubMed  Google Scholar 

  10. Liang C, Zhang X, Yang M, Dong X. Recent progress in ferroptosis inducers for cancer therapy. Adv Mater. 2019;31(51):e1904197.

    Article  PubMed  CAS  Google Scholar 

  11. Ransohoff JD, Wei Y, Khavari PA. The functions and unique features of long intergenic non-coding RNA. Nat Rev Mol Cell Biol. 2018;19(3):143–57.

    Article  CAS  PubMed  Google Scholar 

  12. Mao C, Wang X, Liu Y, Wang M, Yan B, Jiang Y, et al. A G3BP1-interacting lncRNA promotes ferroptosis and apoptosis in cancer via nuclear sequestration of p53. Cancer Res. 2018;78(13):3484–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wang M, Mao C, Ouyang L, Liu Y, Lai W, Liu N, et al. Long noncoding RNA LINC00336 inhibits ferroptosis in lung cancer by functioning as a competing endogenous RNA. Cell Death Differ. 2019;26(11):2329–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wang Z, Chen X, Liu N, Shi Y, Liu Y, Ouyang L, et al. A nuclear long non-coding RNA LINC00618 accelerates ferroptosis in a manner dependent upon apoptosis. Mol Ther. 2021;29(1):263–74.

    Article  CAS  PubMed  Google Scholar 

  15. Zhu L, Yang F, Wang L, Dong L, Huang Z, Wang G, et al. Identification the ferroptosis-related gene signature in patients with esophageal adenocarcinoma. Cancer Cell Int. 2021;21(1):124.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Hassannia B, Vandenabeele P, Vanden BT. Targeting ferroptosis to iron out cancer. Cancer Cell. 2019;35(6):830–49.

    Article  CAS  PubMed  Google Scholar 

  17. Doll S, Freitas FP, Shah R, Aldrovandi M, da Silva MC, Ingold I, et al. FSP1 is a glutathione-independent ferroptosis suppressor. Nature. 2019;575(7784):693–8.

    Article  CAS  PubMed  Google Scholar 

  18. Greenberg SA, Sanoudou D, Haslett JN, Kohane IS, Kunkel LM, Beggs AH, et al. Molecular profiles of inflammatory myopathies. Neurology. 2002;59(8):1170–82.

    Article  CAS  PubMed  Google Scholar 

  19. Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26(12):1572–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  22. Geeleher P, Cox NJ, Huang RS. Clinical drug response can be predicted using baseline gene expression levels and in vitro drug sensitivity in cell lines. Genome Biol. 2014;15(3):R47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.

    Article  CAS  PubMed  Google Scholar 

  24. Hackl H, Charoentong P, Finotello F, Trajanoski Z. Computational genomics tools for dissecting tumour-immune cell interactions. Nat Rev Genet. 2016;17(8):441–58.

    Article  CAS  PubMed  Google Scholar 

  25. Hong W, Liang L, Gu Y, Qi Z, Qiu H, Yang X, et al. Immune-related lncRNA to construct novel signature and predict the immune landscape of human hepatocellular carcinoma. Mol Ther Nucleic Acids. 2020;22:937–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol. 2008;26(8):1364–70.

    Article  PubMed  Google Scholar 

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

  28. Nguyen CT, Campbell SC. Staging of renal cell carcinoma: past, present, and future. Clin Genitourin Cancer. 2006;5(3):190–7.

    Article  PubMed  Google Scholar 

  29. Elkassem AA, Allen BC, Sharbidre KG, Rais-Bahrami S, Smith AD. Update on the role of imaging in clinical staging and restaging of renal cell carcinoma based on the AJCC 8th edition from the AJR special series on cancer staging. AJR Am J Roentgenol. 2021;217(3):541–55.

    Article  PubMed  Google Scholar 

  30. Yu H, Guo P, Xie X, Wang Y, Chen G. Ferroptosis, a new form of cell death, and its relationships with tumourous diseases. J Cell Mol Med. 2017;21(4):648–57.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  33. Lu J, Xu F, Lu H. LncRNA PVT1 regulates ferroptosis through miR-214-mediated TFR1 and p53. Life Sci. 2020;260:118305.

    Article  CAS  PubMed  Google Scholar 

  34. Song X, Long D. Nrf2 and ferroptosis: a new research direction for neurodegenerative diseases. Front Neurosci. 2020;14:267.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Zhong W, Zhang F, Huang C, Lin Y, Huang J. Identification of epithelial-mesenchymal transition-related lncRNA with prognosis and molecular subtypes in clear cell renal cell carcinoma. Front Oncol. 2020;10:591254.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Xiong Y, Wang Z, Zhou Q, Zeng H, Zhang H, Liu Z, et al. Identification and validation of dichotomous immune subtypes based on intratumoral immune cells infiltration in clear cell renal cell carcinoma patients. J Immunother Cancer. 2020;8(1):e000447.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Nakano O, Sato M, Naito Y, Suzuki K, Orikasa S, Aizawa M, et al. Proliferative activity of intratumoral CD8(+) T-lymphocytes as a prognostic factor in human renal cell carcinoma: clinicopathologic demonstration of antitumor immunity. Cancer Res. 2001;61(13):5132–6.

    CAS  PubMed  Google Scholar 

  38. Giraldo NA, Becht E, Pagès F, Skliris G, Verkarre V, Vano Y, et al. Orchestration and prognostic significance of immune checkpoints in the microenvironment of primary and metastatic renal cell cancer. Clin Cancer Res. 2015;21(13):3031–40.

    Article  CAS  PubMed  Google Scholar 

  39. Yang H, Xiong X, Li H. Development and interpretation of a genomic instability derived lncRNAs based risk signature as a predictor of prognosis for clear cell renal cell carcinoma patients. Front Oncol. 2021;11:678253.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Wang Y, Yan K, Wang L, Bi J. Genome instability-related long non-coding RNA in clear renal cell carcinoma determined using computational biology. BMC Cancer. 2021;21(1):727.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhang S, Zhang F, Niu Y, Yu S. Aberration of lncRNA LINC00460 is a promising prognosis factor and associated with progression of clear cell renal cell carcinoma. Cancer Manag Res. 2021;13:6489–97.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Meng D-F, Shao H, Feng C-B. LINC00894 enhances the progression of breast cancer by sponging miR-429 to regulate ZEB1 expression. Onco Targets Ther. 2021;14:3395–407.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Zhu H, Shan Y, Ge K, Lu J, Kong W, Jia C. LncRNA CYTOR promotes pancreatic cancer cell proliferation and migration by sponging miR-205-5p. Pancreatology. 2020;20(6):1139–48.

    Article  CAS  PubMed  Google Scholar 

  44. Hu B, Yang X-B, Yang X, Sang X-T. LncRNA CYTOR affects the proliferation, cell cycle and apoptosis of hepatocellular carcinoma cells by regulating the miR-125b-5p/KIAA1522 axis. Aging (Albany NY). 2020;13(2):2626–39.

    Article  Google Scholar 

  45. Liu H, Zhang X, Jin X, Yang Y, Liang G, Ma Y, et al. Long noncoding RNA VPS9D1-AS1 sequesters microRNA-525-5p to promote the oncogenicity of colorectal cancer cells by upregulating HMGA1. Cancer Manag Res. 2020;12:9915–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Liang X, Chen Z, Wu G. FOXD2-AS1 predicts dismal prognosis for oral squamous cell carcinoma and regulates cell proliferation. Cell Transplant. 2020;29:963689720964411.

    PubMed  Google Scholar 

  47. Xue W, Shen Z, Li L, Zheng Y, Yan D, Kan Q, et al. Long non-coding RNAs MACC1-AS1 and FOXD2-AS1 mediate NSD2-induced cisplatin resistance in esophageal squamous cell carcinoma. Mol Ther Nucleic Acids. 2021;23:592–602.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank the researchers and study participants at the Peking University Institute of Urology for their contributions.

Funding

This study was supported by grants from the National Natural Science Foundation of China [81670617, 81672546, 81602253, 81772703, 81872083], Natural Science Foundation of Beijing [7152146, 7172219], and Wuxi “Taihu Talents Program” Medical and Health High-level Talents Project. All authors would like to thank all contributors to the TCGA and GEO projects. Moreover, all authors thank all researchers at the Institution of Urology, Peking University.

Author information

Author notes

  1. Zhenpeng Zhu, Cuijian Zhang have contributed equally to this work

Authors and Affiliations

  1. Department of Urology, Peking University First Hospital, Beijing, 100034, China

    Zhenpeng Zhu, Cuijian Zhang, Jinqin Qian, Weijie Zhu, Yanqing Gong, Xuesong Li, Jian Lin & Liqun Zhou

  2. Institution of Urology, Peking University, Beijing, 100034, China

    Zhenpeng Zhu, Cuijian Zhang, Jinqin Qian, Weijie Zhu, Yanqing Gong, Xuesong Li, Jian Lin & Liqun Zhou

  3. National Urological Cancer Center, Beijing, 100034, China

    Zhenpeng Zhu, Cuijian Zhang, Jinqin Qian, Weijie Zhu, Yanqing Gong, Xuesong Li, Jian Lin & Liqun Zhou

  4. Beijing Key Laboratory of Urogenital Diseases (Male) Molecular Diagnosis and Treatment Center, Beijing, 100034, China

    Zhenpeng Zhu, Cuijian Zhang, Jinqin Qian, Weijie Zhu, Yanqing Gong, Xuesong Li, Jian Lin & Liqun Zhou

  5. The Affiliated Wuxi Second People’s Hospital of Nanjing Medical University, Wuxi, 214000, China

    Ninghan Feng & Yang Wang

Authors
  1. Zhenpeng Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Cuijian Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jinqin Qian
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ninghan Feng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Weijie Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yang Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yanqing Gong
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xuesong Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jian Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Liqun Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LZ, JL, XL, SH, YG, and ZZ designed and supervised the study, drafted the manuscript, and interpreted the data. LZ, JL, and XL provided the funding. SH and YG provided advice and supervision. ZZ, XJ, WZ, CX, TC, CH performed the experiments, obtained, and analyzed the data. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xuesong Li, Jian Lin or Liqun Zhou.

Ethics declarations

Ethics approval and consent to participate

Not Applicable.

Consent for publication

Not Applicable.

Competing interests

The authors declare that there are no competing of interest.

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. Clinical information of 530 ccRCC patients.

Additional file 2:

Table S2. Detailed information of the 10 paired clinical specimens.

Additional file 3:

Table S3. Paired primers of the lncRNAs in the FerroSig.

Additional file 4:

Table S4. The gene list of the ferroptosis-related genes.

Additional file 5: Figure S1.

KM-Plot for overall survival and disease-free survival, based on 5 LncRNAs in signature.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Zhu, Z., Zhang, C., Qian, J. et al. Construction and validation of a ferroptosis-related long noncoding RNA signature in clear cell renal cell carcinoma. Cancer Cell Int 22, 283 (2022). https://doi.org/10.1186/s12935-022-02700-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-022-02700-0

Keywords

Cancer Cell International

ISSN: 1475-2867

Contact us