Skip to main content
Download PDF

Ferroptosis-related long non-coding RNA signature predicts the prognosis of bladder cancer

BMC Cancer volume 22, Article number: 719 (2022) Cite this article

Abstract

Background

Ferroptosis is an iron-dependent programmed cell death modality that may have a tumor-suppressive function. Therefore, regulating ferroptosis in tumor cells could serve as a novel therapeutic approach. This article focuses on ferroptosis-associated long non-coding RNAs (lncRNAs) and their potential application as a prognostic predictor for bladder cancer (BCa).

Methods

We retrieved BCa-related transcriptome information and clinical information from the TCGA database and ferroptosis-related gene sets from the FerrDb database. Least absolute shrinkage and selection operator regression (LASSO) and Cox regression models were used to identify and develop predictive models and validate the model accuracy. Finally, we explored the inter-regulatory relationships between ferroptosis-related genes and immune cell infiltration, immune checkpoints, and m6A methylation genes.

Results

Kaplan–Meier analyses screened 11 differentially expressed lncRNAs associated with poor BCa prognosis. The signature (AUC = 0.720) could be utilized to predict BCa prognosis. Additionally, GSEA revealed immune and tumor-related pathways in the low-risk group. TCGA showed that the p53 signaling pathway, ferroptosis, Kaposi sarcoma − associated herpesvirus infection, IL − 17 signaling pathway, MicroRNAs in cancer, TNF signaling pathway, PI3K − Akt signaling pathway and HIF − 1 signaling pathway were significantly different from those in the high-risk group. Immune checkpoints, such as PDCD-1 (PD-1), CTLA4, and LAG3, were differentially expressed between the two risk groups. m6A methylation-related genes were significantly differentially expressed between the two risk groups.

Conclusion

A new ferroptosis-associated lncRNAs signature developed for predicting the prognosis of BCa patients will improve the treatment and management of BCa patients.

Peer Review reports

Introduction

Bladder cancer (BCa) is one of the most common malignancies in the urogenital system and one of the top ten predominant malignancies worldwide [1]. Bladder cancer is divided into muscle-invasive bladder cancer (MIBC) and non-muscle-invasive bladder cancer (NMIBC), according to whether the tumor invades the muscle layer of the bladder [2]. Although surgical treatment and postoperative Bacillus Calmette-Guerin (BCG) perfusion and some other immunotherapy are applied to the clinical treatment of BCa [3], about 20% of BCa cases still show an invasion into the bladder muscle. Despite the available treatments, MIBC recurrence, progression, and mortality are high [4]. The five-year overall survival (OS) rate for all stages of urothelial cancer patients is approximately 66–68% [5]. Effective clinical management of BCa is greatly limited by the preclinical models and a lack of accurate biomarkers for early diagnosis. Therefore, it is crucial to explore other forms of cell death to overcome the resistance of tumor cells and discover new and effective prognostic biomarkers for early BCa.

Ferroptosis is a newly discovered form of programmed cell death, operating differently from apoptosis and autophagy [6]. Ferroptosis is iron-dependent because it is triggered by the accumulation of intracellular iron, lipid peroxides, and reactive oxygen species (ROS). The primary mechanism of ferroptosis is the induction of cell death through lipid peroxidation in the presence of divalent iron or ester oxygenase, which catalyzes the high expression of unsaturated fatty acids in cell membrane [7]. Ferroptosis is involved in various critical biological processes, including cancer and neurodegenerative diseases [8, 9]. Furthermore, recent studies have increasingly confirmed that the regulation of ferroptosis could serve as a new therapeutic tool [10], which requires a closer investigation on the link between ferroptosis and cancer.

Long non-coding RNAs (lncRNAs) are a class of non-coding RNA longer than 200 nucleotides in length. LncRNAs can regulate different physiological and biochemical cellular processes via mediating chromosomal modifications, transcriptional activation, and interference [11]. In addition to gene regulation, lncRNAs are involved in various bioregulatory processes, including those related to tumorigenesis, progression, and metastasis [12]. However, current knowledge on the association between ferroptosis, lncRNAs, and cancer is far from comprehensive. Wang et al. showed that in lung cancer, lncRNA LINC00336 regulates tumor progression by inhibiting ferroptosis mechanisms through interaction with ELAVL1 [13]. LncRNA LINK-A, an oncogene, plays a vital role in endogenous tumor suppression and presentation of cancer cell antigens [14]. In addition, the lncRNAMT1DP showed increased sensitivity of non-small cell lung cancer to ferroptosis via regulating the miR-365a-3p/NRF2 signaling pathway [15]. Therefore, lncRNAs can act as an independent prognostic factor for tumors, providing new directions for individualized tumor treatment.

Currently, there are no studies reporting the association of ferroptosis-related lncRNAs with BCa overall survival. This study developed the first prognostic model of differentially expressed ferroptosis-related lncRNAs with prognostic lncRNAs for BCa.

Method

Data collection

BCa transcriptome expression data were retrieved from the TCGA portal (https://cancergenome.nih.gov/). The data included 414 tumor samples and 19 healthy samples. Clinical characteristics of the BCa patients obtained included age, stage, TNM stage, survival time, and survival status. Patients with incomplete information were excluded from our analysis. Samples with OS ≤ 30 days were excluded for non-neoplastic death (Table 1). Corresponding ferroptosis-related genes were downloaded from the FerrDb database [16].FerrDb is a comprehensive, manually curated, and up-to-date database for studying ferroptosis markers and regulators in health and disease. In this study, we identified 247 genes related to triggering effects (Table S1). The relationship between the ferroptosis-related lncRNAs and BCa was assessed using Pearson correlation. The association was considered significant if the correlation coefficient |R2|> 0.3 at P < 0.001. Statistical significance of differential expression of ferroptosis-related lncRNAs was set at a fold-change (FC) value of > 1.0 and FDR-corrected value of P < 0.01.

Table 1 The clinical characteristics of patients in the TCGA dataset
Full size table

Enrichment Analysis of ferroptosis-related DEGs

Functional enrichment analysis of differentially expressed genes (DEGs) was performed using Metascape (http://metascape.org) [17] and the Database for Annotation, Visualization and Integrated Discovery (DAVID) [18]. In addition, functional analysis of biological processes (BP), molecular functions (MF), and cellular components (CC) regulated by the differentially expressed ferroptosis-related lncRNAs were analyzed based on Kyoto Encyclopedia of Genes and Genomes (KEGG) data [19] using R software and Metascape database. The P-value was set at P < 0.05 as a critical value.

Development of the ferroptosis-related lncRNAs prognostic signature

Least absolute shrinkage and selection operator (LASSO, Tibshirani,1996) method is a compression estimation method, which shapes a more refined model by constructing a penalty function that compresses some coefficients and sets some coefficients to zero. LASSO method retains the advantage of subset shrinkage, and it is also a biased estimation for multicollinear data. Thus, it can realize the selection of variables while estimating parameters, and better solves the multicollinearity problem in regression analysis. We employed LASSO-penalized Cox regression analysis and Univariate Cox regression analysis using the "glmnet" package in R to develop the ferroptosis-related lncRNAs signature. The risk score was calculated with the below formula [20], and each BCa patient's risk score was evaluated. With the median value as a bound, the RNAs were divided into low-risk (< median value) and high-risk (≥ median value) groups.

(Coefficient lncRNA1 × expression of lncRNA1) + (Coefficient lncRNA2 × expression of lncRNA2) + ⋯ + (Coefficient lncRNAn × expression lncRNAn).

The predictive nomogram

We performed Gene set enrichment analyses (GSEA [21] to define the lncRNAs signatures in the KEGG pathways and searched in the TCGA-BLCA database. Statistical significance was set at P < 0.05 and false discovery rate (FDR) of q < 0.25. In order to enable clinicians to easily use the prognostic model to evaluate the 1-year, 3-year and 5-year OS of patients with BCa, We combined univariate and multivariate clinical features(gender、grader、age and stage) with significant prognosis and prognostic models, and established nomogram with R software package regplot (https://github.com/cran/regplot).

Immunity analysis and gene expression

The CIBERSORT [22], ESTIMATE [23], MCPcounter [24], single-sample gene set enrichment analysis (ssGSEA) [25] and TIMER [26] algorithms were compared to assess cellular components or cell immune responses between high-risk and low-risk groups based on the ferroptosis-related lncRNAs signature. Differences in immune response under different algorithms were revealed using a Heatmap. In addition, ssGSEA was used to quantify tumor-infiltrating immune cell subgroups between the two groups and assess their immune function. The potential immune checkpoint was also acquired from previous literature.

Cell culture

Bladder cancer cell lines T24 and EJ and normal bladder epithelial cells line SV-HUC were obtained from the Shanghai Branch, Chinese Academy of Sciences. The cells were cultured in RPMI-1640 (ThermoFischer Scientific, Waltham, MA, USA) with 10% fetal calf serum (Sigma-Aldrich, St. Louis, MO, USA) and passaged by 0.25% trypsinization with EDTA (Invitrogen, Grand Island, NY). All the cells were cultured at 37° C in 5% CO2.Tumor cells in logarithmic phase were selected for experiment.

Validation of the diff-lncRNAs

The identified diff-LncRNAs were further validated by real-time-quantitative PCR (RT-qPCR) analysis using the following human cell lines (Bladder cancer cells (T24 and EJ) and normal bladder epithelial cells line SV-HUC). The cell lines were purchased from the Shanghai Institute of Cell Science, Chinese Academy of Sciences. The RT-qPCR was conducted according to the procedures. Briefly, TRIzol (Invitrogen; Thermo Fisher Scientific, Inc., Carlsbad, CA, USA) was used to extract total RNA from the cells. The lncRNAs primers and real-time fluorescent quantitative PCR testing kit were obtained from FulenGen Co., Ltd. (Guangzhou, China). 18S served as an endogenous control. The relative quantification of LncRNA levels was determined by the △△C t method. The synthesis of our first cDNA was synthesized using EntiLink™ 1st Strand cDNA Synthesis Kit (ELK Biotechnology, EQ003), and real-time fluorescent quantitative PCR was conducted on the StepOne™ Real-Time PCR instrument (Life Technologies) using EnTurbo™ SYBR Green PCR SuperMix kit (ELK Biotechnology, EQ001), and 3 double holes were set up for each sample. The specific primer sequence list was shown in Table 2.

Table 2 The specific primer sequence list
Full size table

Drug sensitivity analysis

We use the pRRophetic algorithm [27] to predict the IC50 value of drugs by constructing a relevant ridge regression model. The model takes the expression profile of GDSC cell line(https://www.cancerrxgene.org/) as the training set and TCGA queue as the validation set. we predicted the IC50 values of axitinib, bortezomib, cisplatin, gefitinib, sorafenib, sunitinib, temsirolimus and vinblastine drugs in each sample of TCGA data set. Spearman correlation test was used to analyze the correlation between the expression of lncRNA and the IC50 values of these drugs and cisplatin.

Statistical analysis

Data were analyzed using Bioconductor packages in R software (version 4.0.2). Normal and non-normal distributed variables were analyzed by the unpaired student's t-test and the Wilcoxon test, respectively. Benjamini–Hochberg method was used to identify the differentially expressed lncRNAs based on FDR. The ssGSEA-normalized BCa DEGs were compared with a human genome using "GSVA" (R-package). The sensitivity and specificity of the derived prognostic signatures for BCa in comparison to other clinicopathological was assessed by the receiver operating characteristic (ROC) curve and decision curve analysis (DCA) [28]. The relationship between ferroptosis-related lncRNAs and clinicopathological manifestations was evaluated using logistic regression analyses and a heatmap graph. Finally, the survival analysis of BCa patients based on the ferroptosis-related lncRNAs signature was analzyed using the Kaplan–Meier survival analysis. For each analysis, statistical significance was set at P < 0.05.

Results

Extraction and functional enrichment analysis of differential genes associated with ferroptosis

We began with analyzing the BCa transcriptome expression data obtained from TCGA. By differential enrichment analysis using DAVID database, 61 differentially expressed genes (DEGs) were found to be associated with ferroptosis. Of these genes, 25 were downregulated, and 36 were upregulated (Table S2). BP functional enrichment analysis revealed that these genes were involved in intrinsic apoptotic signaling pathway, multicellular organismal homeostasis, and response to toxic substances. Among the MF specific terms, we found “iron ion binding”, “cargo recrptor activity”, “oxidoreductase activity acting on single donors with incorporation of molecular oxygen”, and “protein kinase inhibitor activity”. Among the CC terms, we mainly found “lipid droplet”, “caveola”, and “chromatin”. KEGG-based analysis showed that overexpressed genes mainly involved the p53 signaling pathway, ferroptosis, Kaposi sarcoma-associated herpesvirus infection, IL-17 signaling pathway, MicroRNAs in cancer, the TNF signaling pathway, PI3K-Akt signaling pathway, and HIF-1 signaling pathway (Fig. 1A-D).

Fig. 1
figure 1

GO and KEGG analyses for ferroptosis-related differentially expressed genes. (A) KEGG. (B) GO. (C-D) GO and KEGG analyses for ferroptosis-related differentially expressed genes by metascape

Full size image

Prognostic features and survival analysis based on ferroptosis-associated lncRNAs

In this study, 518 ferroptosis-related lncRNAs were identified(Table S3), and 34 lncRNAs associated with the prognosis of bladder cancer were screened out by univariate Cox analysis. The result showed that except for AL583785.1 and LINC02762 which were high-risk prognostic LncRNAs, others were low-risk LncRNAs( P < 0.05, Fig. 3C). Multivariate Cox analysis of these obtained lncRNAs showed that 7 of the 11 differentially expressed lncRNAs (Table S4) were independent prognostic indicators of BCa.

We next calculated risk scores and constructed a BCa prognostic model using the lncRNAs. Kaplan–Meier analysis showed that patients with high-risk lncRNAs expression had poorer survival than lncRNAs with low-risk group (P < 0.001, Fig. 2A). The risk model (AUC = 0.720) showed a stronger performance and predictive power than traditional clinicopathological features (Fig. 2B). Using the patients' risk survival status maps, we found that patients' risk scores were inversely correlated with the survival of BCa patients. Interestingly, from the heat map, we found that most of the novel lncRNAs identified in this study were negatively correlated with the risk model (Fig. 2C). The AUC predictive values of these new IncRNA models for 1-year, 3-year, and 5-year survival were 0.720, 0.697, and 0.706, respectively (Fig. 2D).The DCA plot showed the optimal predictive performance of our risk model (Fig. 2E). The heat map showed the visualization of the expression of the 11 ferroptosis-related lncRNAs included in the risk model (Fig. 2C). Univariate and multifactorial Cox analyses showed that lncRNAs model (HR: 1.05, 95CI: 1.03–1.07) and tumor stage (HR: 1.66, 95CI: 1.37–2.03) were independent prognostic factors for OS of BCa patients (Fig. 3A-B). Figure 3D shows the interactions of prognosis-related LncRNAs regulating genes. A heatmap of the association between prognostic model and clinicopathological manifestations of lncRNAs associated with ferroptosis was also analyzed (Fig. 4A). A hybrid column line plot combining clinicopathological characteristics and the novel ferroptosis-related lncRNAs prognostic signature (Fig. 4B) was stable and accurate, therefore it could be applied in clinical management of BCa patients.

Fig. 2
figure 2

Ferroptosis-related lncRNAs signature based on TCGA. A Kaplan–Meier curves result. B The AUC values of the risk factors. C Risk survival status plot. D The AUC of the for the prediction of 1, 3, 5-year survival rate of BCa. E The DCA of the risk factors

Full size image
Fig. 3
figure 3

Univariate and multivariate COX analysis for the expression of ferroptosis-related lncRNAs. (A) Univariate. (B) Multivariate. (C) Univariate COX analysis for the expression of ferroptosis-related lncRNAs. (D) The relationship between the novel lncRNA and mRNA expression

Full size image
Fig. 4
figure 4

A Heatmap for ferroptosis-related lncRNAs prognostic signature and clinicopathological manifestations. B A nomogram for both clinic-pathological factors and prognostic ferroptosis-related lncRNAs

Full size image

Gene set enrichment analysis

The gene set enrichment analysis (GSEA) revealed various immune and tumor-related pathways, which are prognostic signature regulators of most novel lncRNAs associated with ferroptosis, such as Adhesion junction, ECM receptor interaction, Chemokine signaling pathway, B cell receptor signaling pathway, TGF-β signaling pathway, MAPK receptor signaling pathway, Notch signaling pathway, and Bladder cancer. (Fig. 5A). Using the R software package pRRophetic, we predicted the IC50 values of axitinib, bortezomib, cisplatin, gefitinib, sorafenib, sunitinib, temsirolimus and vinblastine drugs in each sample of TCGA data set, and further analyzed the Spearman rank correlation coefficient between the expression of these lncRNAs and these drugs. It can be observed that multiple lncRNAs have significant correlation with the IC50 of various drugs, Bortezomib showed significant negative correlation with eight lncRNAs (Fig. 5B).

Fig. 5
figure 5

A Gene enrichment analysis for ferroptosis-related lncRNAs based on TCGA. B Correlation between the expression of 11 lncRNAs and drugs

Full size image

Expression of immune-related genes

Next, we generated a heat map of the immune response based on CIBERSORT, ESTIMATE, MCPcounter, single sample gene set enrichment analysis (ssGSEA), and TIMER algorithm (Fig. 6A). From ssGSEA on the TCGA-BLCA data and correlation analysis between immune cell subsets and related functions, it was found that T cell functions included checkpoint (suppression), lysis, HLA, inflammatory regulation, co-stimulation, co-inhibition, and type II INF response. Significant differences between low-risk and high-risk patients were detected (Fig. 6B). As checkpoint inhibitor-based immunotherapy strategies are emerging as one of the most promising cancer treatment tools for some drug-resistant tumors, we further explored the differences in immune checkpoint expression between the two groups, and observed significant differences in PDCD-1 (PD-1), CTLA4, LAG3, and BTLA expression (Fig. 6C). In addition, the comparison of m6A-related mRNA expression and the expression of METTL3, RBM15, ZC3H13, YTHDC1, YTHDF1, YTHDF2, HNRNP and FTO in the high-risk and low-risk groups were significant (Fig. 6D).

Fig. 6
figure 6

A Heatmap for immune responses based on CIBERSORT, ESTIMATE, MCPcounter, ssGSEA, and TIMER algorithms among high- and low-risk group. B ssGSEA for the association between immune cell subpopulations and related functions C. Expression of immune checkpoints among high and low BCa risk groups. D The expression of m6A-related genes between high and low BCa risk group

Full size image

Validation of the identified Diff-lncRNAs

Figure 7 showed the results of the qRT-PCR. The expression of lncRNAs (AL031775.1, AC018653.3, AC011468.1, AL583785.1, AC021321.1, AP003352.1, `ETV7-AS1`, U47924.1, AC010326.3) were significantly downregulated in the SV-HUC cell line as compared with the T24 cell line, while LINC02762 expression did not show significant differences in the two cells line. Similarly, the expression of lncRNAs (AL031775.1, AC018653.3, AC011468.1, AL583785.1, AC021321.1, `ETV7-AS1`, U47924.1, AC010326.3) was significantly upregulated in EJ cell lines, different from SV-HUC cells. However, the expressions of AP003352.1 and LINC02762 did not show significant differences in the two cell lines. The result indicated that the eight lncRNAs (AL031775.1, AC018653.3, AC011468.1, AL583785.1, AC021321.1, `ETV7-AS1`, U47924.1, AC010326.3) could be used as ferroptosis-related biomarkers to predict the prognosis of BCa patients. Those results provided new targets for the treatment and management of BCa patients. (The specific primer sequences are listed in Table 2).

Fig. 7
figure 7

The results of qRT-PCR shows the differential expression of iron death-associated LncRNAs associated with bladder cancer prognosis between bladder cancer cells (T24 and EJ cell lines) and normal bladder epithelial cells line SV-HUC. (* Means p < 0.05, ** means p < 0.001, *** means p < 0.0001)

Full size image

Discussion

Ferroptosis or iron-dependent cell death can regulate tumor proliferation, invasion, and progression [29]. Therefore, ferroptosis induction is emerging as a potential anti-cancer therapeutic strategy via triggering tumor cell death, especially for patients with drug-resistant tumors [29]. In this study, we identified 47 DEGs associated with ferroptosis. GO analysis showed that these DEGs were involved in the BP intrinsic apoptotic signaling pathway, multicellular organismal homeostasis, and response to toxic substances. MF terms were mainly related to “regulate iron ion binding”, “oxidoreductase activity”, “acting on single donors with incorporation”, and CC terms were related to chromatin, receptor complex, and endoplasmic reticulum lumen. KEGG signaling pathway prediction analysis further demonstrated that these overexpressed genes were mainly involved in the p53 signaling pathway, Ferroptosis, Kaposi sarcoma-associated herpesvirus infection, IL -17 signaling pathway, MicroRNAs in cancer, TNF signaling pathway, PI3K-Akt signaling pathway, and HIF-1. A recent study showed that p53 could enhance SLC7A11 (solute carrier family 7 member 11), SAT1 (spermidine/spermine N1-acetyltransferase 1), and GLS2 (glutaminase 2) expression by suppressing ferroptosis. On the other hand, p53 inhibits ferroptosis by directly inhibiting DPP4 (dipeptidyl peptidase 4) activity or inducing CDKN1A/p21 (cell cycle protein-dependent kinase inhibitor 1A) expression [30]. Another independent study also confirmed that ferroptosis could be regulated by p53 signaling and tumor-associated mutant p53 (mutp53). The primary manifestation is that the regulation of ferroptosis via p53 contributes to the tumor-suppressive function of p53 and in addition, the accumulation of mutp53 protein in cancer cells increases the sensitivity of cancer cells to ferroptosis [31]. Pretreatment with FG-4592, an inhibitor of prolyl hydroxylase of HIF, reduces renal injury via AKT/GSK-3β-mediated activation of Nrf2 in advanced stages of ferroptosis [32]. In addition, Yadong Sun et al. reported that cancer spheroids proliferate using mammalian targets of rapamycin (mTOR) and utilize the lipid peroxidase GPX4 against ferroptosis [33].

To date, several studies have shown that lncRNAs can act as anti-cancer targets by regulating ferroptosis [34,35,36]. Ferroptosis-related lncRNAs could predict the prognosis of colon cancer patients [37]. LINC00618 expedites ferroptosis through adding lipid reactive oxygen and iron in leukemia and reduces the level of SLC7A11, which accelerates ferroptosis by inducing apoptosis [38]. A recent independent study suggested that iron-dependent cell death-associated lncRNAs can be a prognostic factor for colorectal cancer and HNSCC patients [39, 40]. Therefore, it is of great significance to developed a ferroptosis-related lncRNA prediction signature for BCa patients. In our present study, we constructed a model characterized by 11 lncRNAs (AL031775.1, AC024060.2, AC018653.3, AC011468.1, AL583785.1, AC021321.1, AP003352.1, `ETV7-AS1`, U47924.1, AC010326.3, and LINC02762) associated with iron-dependent cell death to predict the prognosis of BCa patients. Currently, multiple ferroptosis-related lncRNAs have been reported to be associated with a poor prognosis in a variety of tumors. Mei Chen et al. found that AL031775.1, AP003352.1 could estimate the prognosis and development of bladder cancer patients [41]. AC018653.3 was a gene in fifteen kinds of lncRNAs to predict prognosis in colorectal cancer [42]. Six immune-related lncRNAs, including AC011468.1, showed an underlying significance in the prognosis of bladder cancer patients [43]. For seven remaining ferroptosis-related lncRNAs (AC024060.2, AL583785.1, AC021321.1, `ETV7-AS1`, U47924.1, AC010326.3 and LINC02762), there were no studies reporting their prognostic roles in cancers. Therefore, more studies are needed to explore how these lncRNAs affect the prognosis of BC patients through iron failure.

Subsequently, using GSEA, we revealed potential signaling pathways, for example, Adhesion junction, ECM receptor interaction, Chemokine signaling pathway, B cell receptor signaling pathway, TGF-β signaling pathway, MAPK receptor signaling pathway, Notch signaling pathway, and Bladder cancer, for the 11 ferroptosis-associated lncRNAs. It was reported that ferroptosis promotes neutrophil adhesion to coronary vascular endothelial cells via the TLR4/Trif/I type IFN signaling pathway, thereby coordinating neutrophil recruitment to damaged myocardium and promoting programmed cell death [44]. Several recent studies have confirmed the interaction among ferroptosis and immune checkpoint inhibitors and immune cell infiltration [45]. Similarly, lncRNAs play critical roles in ferroptosis. Previous study found that LINC00618 promotes VCR-induced ferroptosis and apoptosis, while LINC00618 accelerates ferroptosis in an apoptosis-dependent manner. LINC00618 attenuated lymphatic-specific decapping enzyme (LSH) expression and LSH enhanced SLC7A11 promoter region after recruitment to SLC7A11 transcription, further inhibiting ferroptosis [38]. LncRNA OIP5-AS1 promotes PCa progression and ferroptosis resistance through miR-128-3p/SLC7A11 signaling [46]. More interestingly, we also discovered a relationship between ferroptosis-associated lncRNAs and m6A methylation genes.

Despite the encouraging data, certain key questions such as the interconnection of ferroptosis with other types of cell death and host immunogenicity remain unclear. This study discovered novel ferroptosis biomarkers for BCa prognosis. These biomarkers could be used to predict and potentially treat BCa. The signature profile should be further investigated using a different cohort as our findings were not validated using clinical samples. Therefore, an in-deep validation with more clinical data from BCa patients is needed before translating our results into clinical practice.

Conclusion

In summary, we developed a BCa prognostic model base on eleven ferroptosis-related genes (AL031775.1, AC024060.2, AC018653.3, AC011468.1, AL583785.1, AC021321.1, AP003352.1`ETV7-AS1`, U47924.1, AC010326.3 and LINC02762). The novel model lays a foundation for further developing new research strategies to explore the mechanisms of ferroptosis and predicting the prognosis of BCa patients.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the [https://portal.gdc.cancer.gov/].

Abbreviations

BCa:

Bladder cancer

TCGA:

The cancer genome atlas

diff-LncRNAs:

Differentially expressed LncRNAs

qRT-PCR:

Real-time quantitative PCR

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

CC:

Cellular component

MF:

Molecular function

BP:

Biological process

FC:

Fold-change

DEGs:

Differentially expressed genes

GSEA:

Gene Set Enrichment Analysis

DAVID:

The Database for Annotation, Visualization and Integrated Discovery

OS:

Overall Survival

ssGSEA:

Single-sample gene set enrichment analysis

EDTA:

Ethylene Diamine Tetraacetic Acid

DCA:

Decision curve analysis

References

  1. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, Bray F. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359-386.

    Article  CAS  PubMed  Google Scholar 

  2. Burger M, Catto JW, Dalbagni G, Grossman HB, Herr H, Karakiewicz P, Kassouf W, Kiemeney LA, La Vecchia C, Shariat S, et al. Epidemiology and risk factors of urothelial bladder cancer. Eur Urol. 2013;63(2):234–41.

    Article  PubMed  Google Scholar 

  3. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  4. Crispen PL, Kusmartsev S. Mechanisms of immune evasion in bladder cancer. Cancer immunology, immunotherapy : CII. 2020;69(1):3–14.

    Article  PubMed  Google Scholar 

  5. Pettenati C, Ingersoll MA. Mechanisms of BCG immunotherapy and its outlook for bladder cancer. Nat Rev Urol. 2018;15(10):615–25.

    Article  CAS  PubMed  Google Scholar 

  6. Huang L, McClatchy DB, Maher P, Liang Z, Diedrich JK, Soriano-Castell D, Goldberg J, Shokhirev M, Yates JR 3rd, Schubert D, et al. Intracellular amyloid toxicity induces oxytosis/ferroptosis regulated cell death. Cell Death Dis. 2020;11(10):828.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, Patel DN, Bauer AJ, Cantley AM, Yang WS, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149(5):1060–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Ren JX, Sun X, Yan XL, Guo ZN, Yang Y. Ferroptosis in Neurological Diseases. Front Cell Neurosci. 2020;14:218.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Jiang M, Qiao M, Zhao C, Deng J, Li X, Zhou C. Targeting ferroptosis for cancer therapy: exploring novel strategies from its mechanisms and role in cancers. Transl Lung Cancer Res. 2020;9(4):1569–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Liang C, Zhang X, Yang M, Dong X. Recent Progress in Ferroptosis Inducers for Cancer Therapy. Adv Mater (Deerfield Beach, Fla). 2019;31(51): e1904197.

    Article  CAS  Google Scholar 

  11. Zhu Y, Wang S, Xi X, Zhang M, Liu X, Tang W, Cai P, Xing S, Bao P, Jin Y, et al. Integrative analysis of long extracellular RNAs reveals a detection panel of noncoding RNAs for liver cancer. Theranostics. 2021;11(1):181–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, Tsai MC, Hung T, Argani P, Rinn JL, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010;464(7291):1071–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wang M, Mao C, Ouyang L, Liu Y, Lai W, Liu N, Shi Y, Chen L, Xiao D, Yu F, 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. Lin A, Li C, Xing Z, Hu Q, Liang K, Han L, Wang C, Hawke DH, Wang S, Zhang Y, et al. The LINK-A lncRNA activates normoxic HIF1α signalling in triple-negative breast cancer. Nat Cell Biol. 2016;18(2):213–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Gai C, Liu C, Wu X, Yu M, Zheng J, Zhang W, Lv S, Li W. MT1DP loaded by folate-modified liposomes sensitizes erastin-induced ferroptosis via regulating miR-365a-3p/NRF2 axis in non-small cell lung cancer cells. Cell Death Dis. 2020;11(9):751.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zhou N, Bao J: FerrDb: a manually curated resource for regulators and markers of ferroptosis and ferroptosis-disease associations. Database : the journal of biological databases and curation 2020, 2020.

  17. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):P3.

    Article  PubMed  Google Scholar 

  19. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545-d551.

    Article  CAS  PubMed  Google Scholar 

  20. Tibshirani R. The lasso method for variable selection in the Cox model. Stat Med. 1997;16(4):385–95.

    Article  CAS  PubMed  Google Scholar 

  21. Powers RK, Goodspeed A, Pielke-Lombardo H, Tan AC, Costello JC. GSEA-InContext: identifying novel and common patterns in expression experiments. Bioinformatics (Oxford, England). 2018;34(13):i555–64.

    Article  CAS  Google Scholar 

  22. Kawada JI, Takeuchi S, Imai H, Okumura T, Horiba K, Suzuki T, Torii Y, Yasuda K, Imanaka-Yoshida K, Ito Y. Immune cell infiltration landscapes in pediatric acute myocarditis analyzed by CIBERSORT. J Cardiol. 2021;77(2):174–8.

    Article  PubMed  Google Scholar 

  23. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, Selves J, Laurent-Puig P, Sautès-Fridman C, Fridman WH, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17(1):218.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160(1–2):48–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Xiao B, Liu L, Li A, Xiang C, Wang P, Li H, Xiao T. Identification and Verification of Immune-Related Gene Prognostic Signature Based on ssGSEA for Osteosarcoma. Front Oncol. 2020;10: 607622.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, Li B, Liu XS. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509-w514.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Geeleher P, Cox N, Huang RS. pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PLoS One. 2014;9(9): e107468.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Van Calster B, Wynants L, Verbeek JFM, Verbakel JY, Christodoulou E, Vickers AJ, Roobol MJ, Steyerberg EW. Reporting and Interpreting Decision Curve Analysis: A Guide for Investigators. Eur Urol. 2018;74(6):796–804.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Manz DH, Blanchette NL, Paul BT, Torti FM, Torti SV. Iron and cancer: recent insights. Ann N Y Acad Sci. 2016;1368(1):149–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kang R, Kroemer G, Tang D. The tumor suppressor protein p53 and the ferroptosis network. Free Radical Biol Med. 2019;133:162–8.

    Article  CAS  Google Scholar 

  31. Zhang C, Liu J, Liang Y, Wu R, Zhao Y, Hong X, Lin M, Yu H, Liu L, Levine AJ, et al. Tumour-associated mutant p53 drives the Warburg effect. Nat Commun. 2013;4:2935.

    Article  PubMed  CAS  Google Scholar 

  32. Li X, Zou Y, Xing J, Fu YY, Wang KY, Wan PZ, Zhai XY. Pretreatment with Roxadustat (FG-4592) Attenuates Folic Acid-Induced Kidney Injury through Antiferroptosis via Akt/GSK-3β/Nrf2 Pathway. Oxid Med Cell Longev. 2020;2020:6286984.

    PubMed  PubMed Central  Google Scholar 

  33. Sun Y, Berleth N, Wu W, Schlütermann D, Deitersen J, Stuhldreier F, Berning L, Friedrich A, Akgün S, Mendiburo MJ, et al. Fin56-induced ferroptosis is supported by autophagy-mediated GPX4 degradation and functions synergistically with mTOR inhibition to kill bladder cancer cells. Cell Death Dis. 2021;12(11):1028.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Mao C, Wang X, Liu Y, Wang M, Yan B, Jiang Y, Shi Y, Shen Y, Liu X, Lai W, et al. A G3BP1-Interacting lncRNA Promotes Ferroptosis and Apoptosis in Cancer via Nuclear Sequestration of p53. Can Res. 2018;78(13):3484–96.

    Article  CAS  Google Scholar 

  35. Qi W, Li Z, Xia L, Dai J, Zhang Q, Wu C, Xu S. LncRNA GABPB1-AS1 and GABPB1 regulate oxidative stress during erastin-induced ferroptosis in HepG2 hepatocellular carcinoma cells. Sci Rep. 2019;9(1):16185.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Yang Y, Tai W, Lu N, Li T, Liu Y, Wu W, Li Z, Pu L, Zhao X, Zhang T, et al. lncRNA ZFAS1 promotes lung fibroblast-to-myofibroblast transition and ferroptosis via functioning as a ceRNA through miR-150-5p/SLC38A1 axis. Aging. 2020;12(10):9085–102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Cai HJ, Zhuang ZC, Wu Y, Zhang YY, Liu X, Zhuang JF, Yang YF, Gao Y, Chen B, Guan GX. Development and validation of a ferroptosis-related lncRNAs prognosis signature in colon cancer. Bosn J Basic Med Sci. 2021;21(5):569–76.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Wang Z, Chen X, Liu N, Shi Y, Liu Y, Ouyang L, Tam S, Xiao D, Liu S, Wen F, et al. A Nuclear Long Non-Coding RNA LINC00618 Accelerates Ferroptosis in a Manner Dependent upon Apoptosis. Molecular therapy : the journal of the American Society of Gene Therapy. 2021;29(1):263–74.

    Article  CAS  Google Scholar 

  39. Tang Y, Li C, Zhang YJ, Wu ZH. Ferroptosis-Related Long Non-Coding RNA signature predicts the prognosis of Head and neck squamous cell carcinoma. Int J Biol Sci. 2021;17(3):702–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Zhang W, Fang D, Li S, Bao X, Jiang L, Sun X. Construction and Validation of a Novel Ferroptosis-Related lncRNA Signature to Predict Prognosis in Colorectal Cancer Patients. Front Genet. 2021;12: 709329.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Tong H, Li T, Gao S, Yin H, Cao H, He W. An epithelial-mesenchymal transition-related long noncoding RNA signature correlates with the prognosis and progression in patients with bladder cancer. Biosci Rep. 2021;41(1):1–13.

  42. Li N, Shen J, Qiao X, Gao Y, Su HB, Zhang S. Long Non-Coding RNA Signatures Associated with Ferroptosis Predict Prognosis in Colorectal Cancer. Int J Gen Med. 2022;15:33–43.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Zhao K, Zhang Q, Zeng T, Zhang J, Song N, Wang Z. Identification and validation of a prognostic immune-related lncRNA signature in bladder cancer. Transl Androl Urol. 2021;10(3):1229–40.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Li W, Feng G, Gauthier JM, Lokshina I, Higashikubo R, Evans S, Liu X, Hassan A, Tanaka S, Cicka M, et al. Ferroptotic cell death and TLR4/Trif signaling initiate neutrophil recruitment after heart transplantation. J Clin Investig. 2019;129(6):2293–304.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Tang R, Liu X, Liang C, Hua J, Xu J, Wang W, Meng Q, Liu J, Zhang B, Yu X, et al. Deciphering the Prognostic Implications of the Components and Signatures in the Immune Microenvironment of Pancreatic Ductal Adenocarcinoma. Front Immunol. 2021;12: 648917.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Zhang Y, Guo S, Wang S, Li X, Hou D, Li H, Wang L, Xu Y, Ma B, Wang H, et al. LncRNA OIP5-AS1 inhibits ferroptosis in prostate cancer with long-term cadmium exposure through miR-128-3p/SLC7A11 signaling. Ecotoxicol Environ Saf. 2021;220: 112376.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Special thanks to the researchers that collected, curated, and maintain the TCGA data, whose high-quality work and effort have made studies like this possible.

Funding

This work was supported in part by grants from Natural Science Foundation of Hunan Province (#2021JJ50069).

Author information

Author notes

  1. Jian Hou, Zhenquan Lu, Xiaobao Cheng and Runan Dong contributed equally to this work.

Authors and Affiliations

  1. Department of Urology, Zhuzhou Central Hospital, Zhuzhou, 412007, China

    Genyi Qu & Yong Xu

  2. Department of Surgery, Division of Urology, The University of Hongkong-ShenZhen Hospital, Shenzhen, 518000, China

    Jian Hou, Zhenquan Lu, Xiaobao Cheng, Runan Dong, Yi Jiang & Guoqing Wu

Authors
  1. Jian Hou
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zhenquan Lu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xiaobao Cheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Runan Dong
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yi Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Guoqing Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Genyi Qu
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Yong Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Jian Hou, Zhenquan Lu, Xiaobao Cheng and Runan Dong wrote the main manuscript text, GenYi Qu and Yong Xu performed experiments, Yi Jiang and Guoqing Wu collected data, All the authors reviewed the manuscript and discussed the results and edited the manuscript.

Corresponding authors

Correspondence to Genyi Qu or Yong Xu.

Ethics declarations

Ethics approval and consent to participate

Not applicable, data was collected from public data repositories.

Consent for publication

Not applicable.

Competing interests

All authors declare that they have no conflict of interests.

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

Hou, J., Lu, Z., Cheng, X. et al. Ferroptosis-related long non-coding RNA signature predicts the prognosis of bladder cancer. BMC Cancer 22, 719 (2022). https://doi.org/10.1186/s12885-022-09805-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09805-9

Keywords

BMC Cancer

ISSN: 1471-2407

Contact us