Identification of a Nomogram from Ferroptosis-Related Long Noncoding RNAs Signature to Analyze Overall Survival in Patients with Bladder Cancer

Cui, Yuanshan; Zhou, Zhongbao; Chai, Yumeng; Che, Xuanyan; Zhang, Yong

doi:https://doi.org/10.1155/2021/8533464

Journal of Oncology

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Prognostic Models Based on Machine Learning for Clinical Cancer Research

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 8533464 | https://doi.org/10.1155/2021/8533464

Identification of a Nomogram from Ferroptosis-Related Long Noncoding RNAs Signature to Analyze Overall Survival in Patients with Bladder Cancer

Yuanshan Cui,^1,2Zhongbao Zhou,¹Yumeng Chai,¹Xuanyan Che,¹and Yong Zhang¹

Academic Editor: Zhixiong Liu

Received02 Jun 2021

Revised09 Jul 2021

Accepted13 Aug 2021

Published26 Aug 2021

Abstract

Purpose. This study aimed to establish a nomogram to predict the overall survival (OS) of patients with bladder cancer (BC) by ferroptosis-related long noncoding RNAs (FRlncRNAs) signature. Methods. We obtained FRlncRNAs expression profiles and clinical data of patients with BC from the Cancer Genome Atlas database. The patients were divided into the training set, testing set, and overall set. Lasso regression and multivariate Cox regression were used to establish the FRlncRNAs signature, the prognosis of each group was compared by Kaplan–Meier (K-M) analysis, and the receiver operating characteristic (ROC) curve evaluated the accuracy of the model. The Gene Set Enrichment Analysis (GSEA) was used for the visualization of the functional enrichment for FRlncRNAs. The databases of GEPIA and K-M Plotter were used for subsequent functional analysis of major FRlncRNAs. Results. Thirteen prognostic FRlncRNAs (LINC00942, MAFG-DT, AL049840.3, AL136084.3, OCIAD1-AS1, AC062017.1, AC008074.2, AC018653.3, AL031775.1, USP30-AS1, LINC01767, AC132807.2, and AL354919.2) were identified to be significantly different, constituting an FRlncRNAs signature. Patients with BC were divided into low-risk group and high-risk group by this signature in the training, testing, and overall sets. K-M analysis showed that the prognosis of patients in the high-risk group was poor and the difference in the subgroup analyses was statistically significant. ROC analysis revealed that the predictive ability of the model was more accurate than traditional assessment methods. A risk score based on FRlncRNAs signature was an independent prognostic factor for the patients with BC (HR = 1.388, 95%CI = 1.228–1.568, ). Combining the FRlncRNAs signature and clinicopathological factors, a predictive nomogram was constructed. The nomogram can accurately predict the overall survival of patients and had high clinical practicability. The GSEA analysis showed that the primary pathways were WNT, MAPK, and cell-matrix adhesion signaling pathways. The major FRlncRNAs (MAFG-DT) were associated with poor prognosis in the GEPIA and K-M Plotter database. Conclusion. Thirteen prognostic FRlncRNAs and their nomogram were accurate tools for predicting the OS of BC, which might be molecular biomarkers and therapeutic targets.

1. Introduction

Bladder cancer (BC) is one of the most common cancers worldwide, with nearly 430,000 newly diagnosed patients each year [1]. Approximately 90% of BC belong to urothelial carcinoma according to the histological origin [2]. However, unlike other tumors with a single direction of progression, BC is divided into two types, including nonmuscle invasive bladder cancer (NMIBC) and muscle invasive bladder cancer (MIBC) [3]. The two are not different just in the depth infiltration but have distinctly different characteristics and outcomes [4]. The obvious characteristic of NMIBC is high local recurrence rate [5]. Although NMIBC can be treated by transurethral resection, about 20% of patients with NMIBC progress to MIBC within 5 years, followed by lymphatic metastasis and distant organ metastasis [6, 7]. Although the patients undergo radical cystectomy, the 5-year survival rate of patients with MIBC is still less than 50% [8]. Therefore, it is very meaningful to identify molecular markers of tumor progression and early metastasis.

Ferroptosis, as iron-dependent programmed cell death, is mainly characterized by lipid peroxidation [9]. With the development of molecular mechanism research, ferroptosis has been proved to be involved in a variety of important pathophysiological processes, mainly including the occurrence of cancer and the formation of ischemia-reperfusion injury and neurodegenerative diseases [10]. Thus, modulating ferroptosis in tumor cells may be a novel therapeutic modality [11]. An increasing number of studies support the involvement of ferroptosis in the pathophysiology of BC development and progression [12, 13].

Long noncoding RNAs (lncRNAs) were found to play a wide range of roles in a variety of important biological processes, including cell proliferation and differentiation, regulation of gene expression, RNA translation, and regulation of microRNAs [14, 15]. In BC, ferroptosis had dual and contradictory roles in the process of tumorigenesis, but the exact mechanisms that lead to ferroptosis during cancer have not been clearly studied [16]. Ferroptosis-related long noncoding RNAs (FRlncRNAs) can be involved in the invasion, metastasis, prognosis, and chemoresistance of BC via modulating ferroptosis. The FRlncRNAs expression profiles of TCGA database were not performed to explore novel biomarkers for forecasting the prognosis of BC. Therefore, we aimed to utilize TCGA database to establish FRlncRNAs signature and seek new biomarkers to predict the prognosis of patients with BC.

2. Materials and Methods

2.1. Data Collection

We followed the methods of Dr. Zhou et al. [15]. We obtained RNA sequencing (RNA-seq) of patients with BC from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/) database. The inclusion criteria were as follows: (a) patients were diagnosed as BC; (b) patients had complete lncRNA expression profiles and clinical data. According to the inclusion criteria, 414 patients were included in our study. In addition, the clinical data (mainly including age, gender, pathological stage, TNM stage, survival status, survival time, and survival prognosis) were also downloaded from TCGA. We excluded the patients who had a follow-up time of less than 30 days and incomplete RNA-seq and clinical data. The study did not need the approval of the ethics committee because the data involved in this study were all from TCGA database and strictly complied with TCGA guidelines (http://cancergenome.nih.gov/abouttcga/policies/publicationguidelines).

2.2. Identification of FRlncRNAs and Differentially Expressed Genes

We obtained the profiles of lncRNAs from all the RNA-seq data and standardized the RNA-seq data by log2 transformation. A total of 14,142 lncRNAs were isolated from the TCGA-BLCA dataset. The Gene Set Enrichment Analysis (GSEA, http://www.gsea-msigdb.org/gsea/index.jsp) was used to download the list of 60 ferroptosis-related genes (FRGs), of which 59 FRGs were expressed in BC. Pearson analysis was used to estimate the correlation degree between lncRNAs and FRGs. The square of correlation coefficient ∣R2∣> 0.3, and was considered as FRlncRNAs. Finally, 1,810 FRlncRNAs were selected.

2.3. Identification of the FRlncRNAs Signature

We divided 395 patients into training set (197 cases) and testing set (198 cases) according to a ratio of 1 : 1. The basic characteristics of each group are shown in Table 1. In the training set, univariate Cox regression analysis was used to evaluate the prognostic value of FRlncRNAs. If FRlncRNAs showed , they would be included in the least absolute shrinkage and selection operator (lasso) regression. We visualized the coexpression of FRlncRNAs and FRGs using Cytoscape 3.6.0 and the correlation between FRlncRNAs using “igraph package” and “reshape2 package” in R software. Then, the results of lasso regression were included in the multivariate Cox regression model to establish the risk score of each patient. A risk score (∗(expression of lncRNA_i)) was established by the expression level of FRlncRNA multiplied by regression coefficient (β). Based on the median of risk score, the patients were divided into high-risk group and low-risk group, and the survival rates between the two groups were compared by log-rank test.

In addition, the same formula was used to calculate the risk score of each patient in the testing set and overall set to verify the stability of the established model. The Kaplan-Meier (K-M) curve was conducted to analyze the survival outcomes of each set. The receiver operating characteristic (ROC) curve and its area under the curve (AUC) value were applied to evaluate the specificity and sensitivity of the established model by “ROC package” in R software.

2.4. Construction of the Nomogram

An independent prognostic model was established by Cox regression analysis. The 1 -, 3 -, and 5-year survival rates were predicted by establishing nomogram. The stability of the model was evaluated by the index of concordance (C-index), calibration curves, and ROC curves. The basic characteristics of patients were included in the multivariate Cox regression to determine whether the risk score was an independent predictor of prognosis.

2.5. Functional Analysis of FRlncRNAs

The functional enrichment of FRlncRNAs was explained by the Gene Set Enrichment Analysis (GSEA) (http://www.broadinstitute.org/gsea/index.jsp). The study accessed the functional enrichment of FRlncRNAs and visualized the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to ferroptosis. In the enrichment analysis, a two-sided value of <0.05 was considered significant.

2.6. Verification of Major FRlncRNAs in External Databases

The Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/) database contained RNA-seq and clinical data compiled by TCGA and GTEx databases after standardized transformation. The K-M Plotter database (http://kmplot.com/analysis/) included RNA-seq and prognostic data of 405 BLCA patients. The above two databases were used to verify the relationship between major FRlncRNAs and patient’s prognosis.

2.7. Statistical Analysis

The K-M analysis was utilized to generate the survival curve in this study, and the log-rank test was used to compare whether there was a difference between the two groups of patients. The Cox regression and lasso regression were applied to access the prognostic effect of FRlncRNAs signature and clinicopathological data. R software (version 3.6) was used for statistical analysis. Statistical tests were two-sided, and studies were considered statistically significant when .

3. Results

3.1. Identification of the FRlncRNAs Signature

The flowchart of this research is shown in Figure 1. By univariate Cox regression, we screened a total of 177 FRlncRNAs with prognostic value for patients with BC (). Then, we identified 29 FRlncRNAs by lasso regression (Figures 2(a) and 2(b) and Table S1). By multivariate Cox regression, thirteen FRlncRNAs were independent prognostic factors and their coefficients are also shown in Table 2. An FRlncRNAs signature was established based on these 13 FRlncRNAs. The formula of risk score was as follows: risk score = (0.0153LINC00942) + (0.0550MAFG-DT) + (0.1824AL049840.3) + (0.1928 AL136084.3) – (0.3582OCIAD1-AS1) – (0.3362AC062017.1) – (0.2606 AC008074.2) – (0.2419AC018653.3) – (0.2375 AL031775.1) – (0.1942USP30-AS1) – (0.1094 LINC01767) – (0.0730 AC132807.2) – (0.0662AL354919.2).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 2

Construction and evaluation of the FRlncRNAs signature in the training set. Lasso coefficient values and vertical dashed lines were calculated at the best log (lambda) value (a), and coefficients (b) of prognostic-related lncRNAs are displayed. (c) The K-M curve showed that the high-risk group had a worse survival rate than the low-risk group (). (d) The ROC curve is given for this signature and its AUC value. (e) Heatmap of the 13 FRlncRNAs profiles showed the expression of FRlncRNAs in the high-risk and the low-risk group. (f) Scatter plot showed the correlation between the survival status and risk score of patients. (g) Risk score distribution plot showed the distribution of high-risk and low-risk patients. (h) ROC curves and their AUC value represented 1-, 3-, and 5-year OS.

3.2. Prognostic Influence of the FRlncRNAs Signature

To access the sensitivity and stability of this signature, the patients in the training set were divided into low-risk group (99 cases) and high-risk group (98 cases) based on the median of risk scores. The K-M curve showed a shorter overall survival (OS) in the high-risk group compared with the low-risk group () (Figure 2(c)). Then, the study used ROC curve to assess the accuracy of the signature and found an AUC value of 0.776 in the training set (Figure 2(d)). The heatmap showed a significant difference in the expression of 13 FRlncRNAs between the high-risk group and the low-risk group (Figure 2(e)). The scatter plot identified that patients with high-risk score had a lower OS than patients with low-risk score (Figure 2(f)). In addition, the distribution map of risk score demonstrated that the high-risk group had a higher risk score than the low-risk group (Figure 2(g)). The AUC value corresponding to 1, 3, and 5 years of OS was 0.787, 0.779, and 0.807 (Figure 2(h)). Besides, we used K-M curves to evaluate the prognostic effect of 13 FRlncRNAs. Among them, four FRlncRNAs (LINC00942, MAFG-DT, AL049840.3, and AL136084.3) were detrimental to the prognosis of patients and nine FRlncRNAs (OCIAD1-AS1, AC062017.1, AC008074.2, AC018653.3, AL031775.1, USP30-AS1, LINC01767, AC132807.2, and AL354919.2) were favorable to the prognosis of patients (Figure 3). In short, this prognostic-related signature had good stability and sensitivity in predicting the OS of patients.

3.3. Validation of the FRlncRNAs Signature

To evaluate the predictive power of the FRlncRNAs signature, we calculated the risk score of each patient in the testing set and overall set using the same method and divided patients into the low-risk group and the high-risk group. We used the K-M curve to analyze the OS of patients in the testing set () (Figure 4(a)) and overall set () (Figure 4(c)), and the final results were consistent with the training set. Subsequently, the ROC curve showed that this signature has a strong predictive ability in the testing set (AUC = 0.866) (Figure 4(b)) and overall set (AUC = 0.811) (Figure 4(d)) for the OS of patients. Furthermore, ROC time curves and its AUC value also confirmed that this signature had a better prognostic ability for patients in the testing set (1-year AUC = 0.753, 3-year AUC = 0.803, and 5-year AUC = 0.812) (Figure 4(e)) and overall set (1-year AUC = 0.770, 3-year AUC = 0.787, and 5-year AUC = 0.807) (Figure 4(g)). Besides, heatmaps showed that the expression profiles of 13 FRlncRNAs were also consistent with those in the training set (Figures 4(f) and 4(h)). Also, the scatter plot identified that patients with high-risk score had a lower survival rate than patients with low-risk score, and the distribution map of risk score demonstrated that the high-risk group had a higher risk score than the low-risk group (Figures 4(i)–4(l)). These results indicated that the FRlncRNAs signature had a stable prognostic-predictive ability.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 4

Validation of the FRlncRNAs signature for BC patients in the testing set and overall set. K-M curves showed that the high-risk group had the worse OS than the low-risk group in the testing set (a) and overall set (c). ROC curves and its AUC value in the testing set (b) and overall set (d). ROC curves and their AUC value represented 1-, 3-, and 5-year OS in the testing set (e) and overall set (g). Heatmap of 13 FRlncRNAs profiles showed the expression of FRlncRNAs in the high-risk group and the low-risk group in the testing set (f) and overall set (h). Scatter plot showed the outcomes between the survival status and risk score in the testing set (i) and overall set (j). Risk score distribution plot showed the distribution of high-risk and low-risk BC patients in the testing set (k) and overall set (l).

3.4. Construction and Evaluation of the Prognostic Nomogram

Based on the result of univariate Cox regression in the overall set, the risk score and stage were independent prognostic indicators in patients with BC, in which the HR of risk score was 1.414 (95% CI: 1.272–1.573, , Figure 5(a) and Table 3). The risk score remained an independent prognostic indicator in multivariate Cox regression after controlling for clinical characteristics (HR = 1.388, 95% CI = 1.228–1.568, , Figure 5(b)). Subsequently, we included the risk score, patient age, and tumor stage in the nomogram. The nomogram showed that risk score and tumor stage had the greatest contribution to 1 -, 3 -, and 5-year OS in patients with BC (Figure 5(c)). The higher the risk score of samples, the worse the prognosis of patients. The multivariate ROC curve of risk score based on the FRlncRNAs signature and clinicopathologic characteristics indicated that AUC value was 0.811, which was the higher than AUC value of age (0.544), gender (0.429), stage (0.672), T stage (0.644), M stage (0.524), and N stage (0.654), indicating that the predictive power of nomogram was more accurate than that of stage and TNM stage (Figure 5(d)). These results showed that the FRlncRNAs signature can be a good indicator in predicting the prognosis of patients compared with the predictive power of existing lncRNA-related signatures reported in recent studies (Table 4) [17–22]. The C-index of the nomogram was 0.755 (se = 0.032). By AUC of 1-, 3-, and 5-year survival rate, the nomogram had a better prediction for the prognosis of patients with BC (Figure 5(e)).

(a)

(b)

(c)

(d)

(e)

Figure 5

Clinical value of the FRlncRNAs signature in BC patients. The univariate Cox regression showed that risk score and clinicopathological features including age, stage, T stage, and N stage were prognostic-related variables (a). The multivariate Cox regression analysis showed that the risk score was independent prognostic factors (b). Construction of a prognostic nomogram based on risk score and clinicopathological indexes to predict 1-, 3-, and 5-year OS of BC patients (c). The multivariate ROC curve showed that predictive accuracy of risk score was higher than other clinicopathological indexes (d). Calibration curves displayed the concordance between predicted and observed 1-, 3-, and 5-year OS (e).

3.5. Stratification Analysis of the Clinicopathological Index

To further evaluate the FRlncRNAs signature and verify its stability in predicting the OS of patients in the high-risk group and the low-risk group, we conducted a stratified analysis based on clinicopathological indexes, including gender (female and male), age (≤65 years and >65 years), stage (I-II and III-IV stages), T stage (T1-2 stage and T3-4 stage), and N stage (N0 stage and N1-3 stages).

The result of the K-M curve showed that the OS of the high-risk group was worse than that of the low-risk group in different clinical stratification () (Figures 6(a)–6(j)). Clinical influences of risk score for patients in the training, testing, and overall sets are shown in Table 5.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

3.6. Gene Functional Analysis

We totally extracted 936 GO pathways (Table S2) and 51 KEGG pathways (Table S3) from GSEA 4.1.0. KEGG analysis revealed that the signaling pathway such as WNT signaling pathway, mitogen-activated protein kinase (MAPK), and ERBB signaling pathway was significantly enriched in the high-risk group (Figures 7(a)–7(d)). Go analysis found that the functions of FRlncRNAs mainly focused on physiological processes, such as cell-matrix adhesion and positive regulation of cell division (Figures 7(e) and 7(f)).

(a)

(b)

(c)

(d)

(e)

(f)

3.7. Construction of the Coexpression Network and Verification of Major FRlncRNAs

The Sankey diagram showed the association between FRlncRNAs, FRGs, and risk types (Figure 8(a)). The correlation of 13 FRlncRNAs is shown in Figure 8(b). The coexpression network between prognostic FRlncRNAs and FRGs is shown in Figure 8(c). After literature search, we selected MAFG-DT (MAFG-AS1) for further study. We verified the expression levels and survival outcomes of MAFG-DT in two external databases (GEPIA and K-M Plotter). In the cohort from the GEPIA database, the expression of MAFG-DT in tumor tissues was higher than that in normal tissues (Figure 9(a)), and the expression level of MAFG-DT differed across the three pathological stages, which may indicate that MAFG-DT was closely related to the prognosis of patients (Figure 9(b)). As shown in the GEPIA database, high MAFG-DT expression levels were associated with poor prognosis () (Figures 9(c) and 9(d)). In the 405 patients from the K-M Plotter database, 5-year OS in the group with high expression of MAFG-DT was lower than that in the group with low expression () (Figures 9(e) and 9(f)). In order to intuitively show the expression correlation of MAFG-DT and FRGs, we made a linear correlation graph including GPX4, NCOA4, and SLC1A5 (Figure 9(g)). In conclusion, these results showed that the expression of MAFG-DT was higher in BC tissues, but with the increase of the stage, the expression of MAFG-DT also increased, and there was a significant correlation with survival outcomes.

(a)

(b)

(c)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

4. Discussion

Despite advances in surgery and chemotherapy, the prognosis of patients with advanced and metastatic BC remained unsatisfactory [23]. Furthermore, because of the different causative molecules, although some patients have the same TNM stage or similar risk factors, they may have different clinical outcomes. Therefore, molecular biomarkers to predict tumor prognosis are of great importance [24, 25]. Ferroptosis may be associated with biological behaviors such as proliferation, invasion, and metastasis in BC [6, 7]. And lncRNAs may play key regulatory roles in ferroptosis-related biological processes of malignant tumor cells [26, 27]. However, a prognostic tool based on FRlncRNAs for patients with BC is still lacking.

In our study, FRlncRNAs were collected by establishing a coexpression network of lncRNAs and FRGs. Initially, we identified 29 FRlncRNAs related to prognosis and constructed a prognostic model comprising 13 FRlncRNAs (LINC00942, MAFG-DT, AL049840.3, AL136084.3, OCIAD1-AS1, AC062017.1, AC008074.2, AC018653.3, AL031775.1, USP30-AS1, LINC01767, AC132807.2, and AL354919.2) via multivariate Cox regression and lasso regression. K-M curve showed that the OS of patients with high-risk score were shorter than those of patients with low-risk score. In addition, the ROC curve proved that the 13 FRlncRNAs signature was highly sensitive and specific prognostic markers in BC. Besides, this result was further verified in the testing set and overall set. The 13 FRlncRNAs’ signature was also associated with poor OS of BC patients in different subgroups, especially in age, gender, AJCC stage, T stage, and N stage. Next, a nomogram was used to count the risk score and can predict the OS of patients. The calibration curve showed that the signature had a higher sensitivity and clinical applicability than the traditional standard. Then, we established an FRlncRNA-mRNAs coexpression network to analyze the effect of the 13 FRlncRNAs’ signature. And we also analyzed the correlation between 13 FRlncRNAs. The GSEA analysis revealed that signaling pathways such as WNT, MAPK, ERBB, cell-matrix adhesion, and positive regulation of cell division were significantly enriched in the high-risk group. Finally, we verified the major FRlncRNAs in two external databases. The results showed that the expression of MAFG-DT was higher in BC tissues, but with the increase of the stage, the expression of MAFG-DT also increased, and there was a significant correlation with survival outcomes. Correlation analysis found that MAFG-DT and protein-coding genes (GPX4, NCOA4, and SLC1A5) were obviously related.

Currently, multiple FRlncRNAs have been reported to be associated with poor prognosis in a variety of tumors. Among them, LINC00942 exerts its functions as an oncogene in promoting METTL14-mediated m6A methylation and regulating the expression and stability of its target genes CXCR4 and CYP1B1 in breast cancer (BRCA) initiation and progression, which provides new targets and crosstalk m6A epigenetic modification mechanism for BRCA prevention and treatment [28]. Six FRlncRNAs were included in the relevant clinical prediction models, which showed a good prediction effect for patients with BC. Tong et al. found that AL049840.3, AL031775.1, and USP30-AS1 may predict the prognosis and progression of patients with BC as an epithelial-mesenchymal transition- (EMT-) related lncRNAs [29]. Wang et al. found that OCIAD1-AS1 and AL354919.2 may predict the prognosis and progression of patients with BC as immune-related lncRNAs [30]. USP30-AS1 has also been found to be associated with autophagy, which may be involved in the diagnosis and prognosis of BC [31].

As a major FRlncRNA, MAFG-DT is a tumorigenic lncRNA in a variety of cancers. Firstly, MAFG-DT as a part of lncRNAs signature of tumor-infiltrating B lymphocytes had some predictive value for the survival outcomes and immunotherapy response of patients with antiprogrammed death-1 (PD-1) therapy and added significant predictive power to current immune checkpoint gene markers [32]. Next, Li et al. reported that silencing of MAFG-AS1 inhibited BC cell proliferation, metastasis, and invasion, while overexpression of MAFG-AS1 in BC cell had opposite biological effects, and it was further confirmed that MAFG-AS1/miR-143-3p/COX-2 axis was involved in the progression of BC [33]. Besides, Xiao et al. showed that MAFG-AS1 can promote BC proliferation, invasion, metastasis, and EMT in vitro and in vivo. Mechanistically, MAFG-AS1 directly binding to Hu antigen R (HuR) could recruit ubiquitin-specific proteinase 5 (USP5) to prevent HuR from degrading by ubiquitination [34]. For six remaining FRlncRNAs (AL136084.3, AC062017.1, AC008074.2, AC018653.3, LINC01767, and AC132807.2), there were no studies to report their prognostic roles in cancer. Thus, more researches were necessary to explore how these lncRNAs affect the prognosis of patients with BC through ferroptosis exactly.

Currently, multiple studies have highlighted the role of lncRNAs in the pathogenesis of BC. These lncRNAs were downregulated or upregulated during the origin, proliferation, and migration of BC. Therefore, it is very necessary to clarify the physiological roles of lncRNAs and their contributions to the development and progression of BC. For example, lncRNAs can contribute to mRNA degradation and failure of protein translation by mediating posttranscriptional gene silencing [35]. In addition, lncRNA can also remodel chromatin structure by directing the formation of heterochromatin, thereby affecting the apparent phenotype of cells [36]. LncRNAs can also be involved in both cis- and transgene regulation, thereby enhancing or repressing gene expression [35]. These effects are crucial for the proliferation, differentiation, regeneration, and apoptosis of human cells, thereby helping to maintain the balance of the normal physiological function of the body [37, 38]. However, relative to the wide application of lncRNAs in BC diagnosis, researches of lncRNA-targeting therapies are relatively scarce. Since lncRNAs diagnostic is noninvasive and can detect the occurrence of BC more rapidly compared with cystoscopy, a commonly used diagnostic tool [39].

Ferroptosis can synergistically enhance antitumor activity in combination with immune checkpoint inhibitors (ICIs), even in ICI-resistant types [40]. Currently, there are fewer studies analyzing the relationship between ICIs and ferroptosis. Accumulating evidence has found that microRNAs (miRNAs) and lncRNA were crucial in mediating the regulation of ferroptosis. Nrf2 reduces ROS generation by inhibiting iron absorption [41]. As such, miRNAs can disrupt the process of ferroptosis by regulating the expression of Nrf2 [42, 43]. Meanwhile, miRNAs were involved in regulating iron absorption, transport, storage, and metabolism. In recent years, multiple regulators of ferroptosis have been discovered, including ATF3/4, SLC7A11, and ACSl4. Importantly, lncRNAs play a key role in regulating the expression of these factors [44].

Ferroptosis, as a new form of cell death, may provide a new direction for tumor therapy. However, many key questions remain unresolved, such as how ferroptosis intersects with other cell deaths and host immunogenicity is less well studied. Therefore, this study explored ferroptosis-related biomarkers that can be used to predict BC prognosis, which may provide references for therapeutic modalities of cancer. However, the current study has some weaknesses. First, the data source of this analysis was single, and the amount of data included was small, so the results may have some bias. Second, the study was retrospective, and prospective studies may be needed to demonstrate the prognostic function of FRlncRNAs. Third, to further validate the stability and accuracy of this prognostic model, our established prognostic model needs further analysis in other independent cohorts. Fourth, relevant functional experiments should be performed to further analyze the underlying molecular mechanisms by which FRlncRNAs affect the development and progression of BC.

5. Conclusions

The coexpression network of FRlncRNA-mRNA provided a valuable source for revealing the function of FRlncRNAs in BC. Thirteen FRlncRNAs were considered to be significantly associated with OS of patients with BC. An FRlncRNAs signature that was composed of thirteen FRlncRNAs was used to differentiate patients at different risks, and it was a significantly independent factor for patients with BC. Therefore, the thirteen FRlncRNAs and their signature might be molecular biomarkers and therapeutic targets for patients with BC.

Data Availability

The datasets used to support the findings of this study are available from the TCGA (https://portal.gdc.cancer.gov/repository) database.

Disclosure

Yuanshan Cui and Zhongbao Zhou are considered the co-first authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors contributed to conception and design, manuscript writing, and final approval of manuscript. Yong Zhang provided administrative support. Yuanshan Cui and Zhongbao Zhouc contributed to provision of study materials or patients. Yuanshan Cui and Zhongbao Zhou carried out collection and assembly of data. Yuanshan Cui, Zhongbao Zhou, Yumeng Chai, and Xuanyan Che contributed to data analysis and interpretation. Yuanshan Cui and Zhongbao Zhou contributed equally to this work.

Acknowledgments

This work was supported by the Beijing Municipal Administration of Hospitals’ Ascent Plan, under Code DFL20190502; Beijing Municipal Administration of Hospitals Clinical Medicine Development of Special Funding Support, under Code ZYLX201820; and National Nature Science Foundation of China, under Code 81801429.

Supplementary Materials

Supplementary Table S1: twenty-nine FRlncRNAs identified after univariate Cox regression and lasso regression. Supplementary Table S2: gene set enrichment GO analysis results according to the signature of thirteen FRlncRNAs (top 20 values). Supplementary Table S3: gene set enrichment KEGG analysis results based on the signature of thirteen FRlncRNAs (top 20 values). . (Supplementary Materials)

References

A. M. Kamat, N. M. Hahn, J. A. Efstathiou et al., “Bladder cancer,” The Lancet, vol. 388, no. 10061, pp. 2796–2810, 2016.
View at: Publisher Site | Google Scholar
A. T. Lenis, P. M. Lec, K. Chamie, and M. Mshs, “Bladder cancer,” JAMA, vol. 324, no. 19, pp. 1980–1991, 2020.
View at: Publisher Site | Google Scholar
M. G. K. Cumberbatch, I. Jubber, P. C. Black et al., “Epidemiology of bladder cancer: a systematic review and contemporary update of risk factors in 2018,” European Urology, vol. 74, no. 6, pp. 784–795, 2018.
View at: Publisher Site | Google Scholar
P. G. Yousef and M. Y. Gabril, “An update on the molecular pathology of urinary bladder tumors,” Pathology, Research & Practice, vol. 214, no. 1, pp. 1–6, 2018.
View at: Publisher Site | Google Scholar
V. Soukup, O. Čapoun, D. Cohen et al., “Prognostic performance and reproducibility of the 1973 and 2004/2016 world health organization grading classification systems in non-muscle-invasive bladder cancer: a European association of Urology non-muscle invasive bladder cancer guidelines panel systematic review,” European Urology, vol. 72, no. 5, pp. 801–813, 2017.
View at: Publisher Site | Google Scholar
J. Taylor, E. Becher, and G. D. Steinberg, “Update on the guideline of guidelines: non-muscle-invasive bladder cancer,” BJU International, vol. 125, no. 2, pp. 197–205, 2020.
View at: Publisher Site | Google Scholar
S. A. Funt and J. E. Rosenberg, “Systemic, perioperative management of muscle-invasive bladder cancer and future horizons,” Nature Reviews Clinical Oncology, vol. 14, no. 4, pp. 221–234, 2017.
View at: Publisher Site | Google Scholar
R. Li, M. Metcalfe, J. Kukreja, and N. Navai, “Role of radical cystectomy in non-organ confined bladder cancer: a systematic review,” Bladder Cancer, vol. 4, no. 1, pp. 31–40, 2018.
View at: Publisher Site | Google Scholar
B. Hassannia, P. Vandenabeele, and T. Vanden Berghe, “Targeting ferroptosis to iron out cancer,” Cancer Cell, vol. 35, no. 6, pp. 830–849, 2019.
View at: Publisher Site | Google Scholar
J. P. Friedmann Angeli, D. V. Krysko, and M. Conrad, “Ferroptosis at the crossroads of cancer-acquired drug resistance and immune evasion,” Nature Reviews Cancer, vol. 19, no. 7, pp. 405–414, 2019.
View at: Publisher Site | Google Scholar
S. V. Torti and F. M. Torti, “Iron and cancer: 2020 vision,” Cancer Research, vol. 80, no. 24, pp. 5435–5448, 2020.
View at: Publisher Site | Google Scholar
P. Sharma, T. Shimura, J. K. Banwait, and A. Goel, “Andrographis-mediated chemosensitization through activation of ferroptosis and suppression of β-catenin/Wnt-signaling pathways in colorectal cancer,” Carcinogenesis, vol. 41, no. 10, pp. 1385–1394, 2020.
View at: Publisher Site | Google Scholar
X. Xu, X. Zhang, C. Wei et al., “Targeting SLC7A11 specifically suppresses the progression of colorectal cancer stem cells via inducing ferroptosis,” European Journal of Pharmaceutical Sciences, vol. 152, Article ID 105450, 2020.
View at: Publisher Site | Google Scholar
S. Panni, R. C. Lovering, P. Porras, and S. Orchard, “Non-coding RNA regulatory networks,” Biochimica et Biophysica Acta (BBA)—Gene Regulatory Mechanisms, vol. 1863, no. 6, Article ID 194417, 2020.
View at: Publisher Site | Google Scholar
W. Zhou, S. Zhang, H. B. Li et al., “Development of prognostic indicator based on autophagy-related lncRNA analysis in colon adenocarcinoma,” BioMed Research International, vol. 2020, Article ID 9807918, 14 pages, 2020.
View at: Publisher Site | Google Scholar
X. Chen, R. Kang, G. Kroemer, and D. Tang, “Broadening horizons: the role of ferroptosis in cancer,” Nature Reviews Clinical Oncology, vol. 18, no. 5, pp. 280–296, 2021.
View at: Publisher Site | Google Scholar
H. Wu, Z.-Y. Zhang, Z. Zhang et al., “Prediction of bladder cancer outcome by identifying and validating a mutation-derived genomic instability-associated long noncoding RNA (lncRNA) signature,” Bioengineered, vol. 12, no. 1, pp. 1725–1738, 2021.
View at: Publisher Site | Google Scholar
K. Zhao, Q. Zhang, T. Zeng, J. Zhang, N. Song, and Z. Wang, “Identification and validation of a prognostic immune-related lncRNA signature in bladder cancer,” Translational Andrology and Urology, vol. 10, no. 3, pp. 1229–1240, 2021.
View at: Publisher Site | Google Scholar
W. J. Luo, X. Tian, W. H. Xu et al., “Construction of an immune‐related LncRNA signature with prognostic significance for bladder cancer,” Journal of Cellular and Molecular Medicine, vol. 25, no. 9, pp. 4326–4339, 2021.
View at: Publisher Site | Google Scholar
Y. Du, B. Wang, X. Jiang et al., “Identification and validation of a stromal EMT related LncRNA signature as a potential marker to predict bladder cancer outcome,” Frontiers in oncology, vol. 11, Article ID 620674, 2021.
View at: Publisher Site | Google Scholar
L. Qing, P. Gu, M. Liu et al., “Extracellular matrix-related six-lncRNA signature as a novel prognostic biomarker for bladder cancer,” OncoTargets and Therapy, vol. 13, pp. 12521–12538, 2020.
View at: Publisher Site | Google Scholar
Z. Sun, C. Jing, C. Xiao, and T. Li, “An autophagy-related long non-coding RNA prognostic signature accurately predicts survival outcomes in bladder urothelial carcinoma patients,” Aging, vol. 12, no. 15, pp. 15624–15637, 2020.
View at: Publisher Site | Google Scholar
R. Nadal and J. Bellmunt, “Management of metastatic bladder cancer,” Cancer Treatment Reviews, vol. 76, pp. 10–21, 2019.
View at: Publisher Site | Google Scholar
C. Zhang, J. Hu, H. Li et al., “Emerging biomarkers for predicting bladder cancer lymph node metastasis,” Frontiers in oncology, vol. 11, Article ID 648968, 2021.
View at: Publisher Site | Google Scholar
P. K. Parizi, F. Yarahmadi, H. M. Tabar et al., “MicroRNAs and target molecules in bladder cancer,” Medical Oncology, vol. 37, no. 12, p. 118, 2020.
View at: Publisher Site | Google Scholar
Q. Ma, X. Dai, W. Lu, X. Qu, N. Liu, and C. Zhu, “Silencing long non-coding RNA MEG8 inhibits the proliferation and induces the ferroptosis of hemangioma endothelial cells by regulating miR-497-5p/NOTCH2 axis,” Biochemical and biophysical research communications, vol. 556, pp. 72–78, 2021.
View at: Publisher Site | Google Scholar
J. Lu, F. Xu, and H. Lu, “LncRNA PVT1 regulates ferroptosis through miR-214-mediated TFR1 and p53,” Life Sciences, vol. 260, Article ID 118305, 2020.
View at: Publisher Site | Google Scholar
T. Sun, Z. Wu, X. Wang et al., “LNC942 promoting METTL14-mediated m6A methylation in breast cancer cell proliferation and progression,” Oncogene, vol. 39, no. 31, pp. 5358–5372, 2020.
View at: Publisher Site | Google Scholar
H. Tong, T. Li, S. Gao, H. Yin, H. Cao, and W. He, “An epithelial-mesenchymal transition-related long noncoding RNA signature correlates with the prognosis and progression in patients with bladder cancer,” Bioscience Reports, vol. 41, no. 1, 2021.
View at: Publisher Site | Google Scholar
J. Wang, C. Shen, D. Dong, X. Zhong, Y. Wang, and X. Yang, “Identification and verification of an immune-related lncRNA signature for predicting the prognosis of patients with bladder cancer,” International Immunopharmacology, vol. 90, Article ID 107146, 2021.
View at: Publisher Site | Google Scholar
J. Wan, C. Guo, H. Fang, Z. Xu, Y. Hu, and Y. Luo, “Autophagy-related long non-coding RNA is a prognostic indicator for bladder cancer,” Frontiers in oncology, vol. 11, Article ID 647236, 2021.
View at: Publisher Site | Google Scholar
M. Zhou, Z. Zhang, S. Bao et al., “Computational recognition of lncRNA signature of tumor-infiltrating B lymphocytes with potential implications in prognosis and immunotherapy of bladder cancer,” Briefings in Bioinformatics, vol. 22, no. 3, 2021.
View at: Publisher Site | Google Scholar
D. Li, S. Zhong, Z. Zhu et al., “LncRNA MAFG-AS1 promotes the progression of bladder cancer by targeting the miR-143-3p/COX-2 Axis,” Pathobiology, vol. 87, no. 6, pp. 345–355, 2020.
View at: Publisher Site | Google Scholar
M. Xiao, J. Liu, L. Xiang et al., “MAFG-AS1 promotes tumor progression via regulation of the HuR/PTBP1 axis in bladder urothelial carcinoma,” Clinical and Translational Medicine, vol. 10, no. 8, p. e241, 2020.
View at: Publisher Site | Google Scholar
S. Toden, T. J. Zumwalt, and A. Goel, “Non-coding RNAs and potential therapeutic targeting in cancer,” Biochimica et Biophysica Acta (BBA)—Reviews on Cancer, vol. 1875, no. 1, Article ID 188491, 2021.
View at: Publisher Site | Google Scholar
K. Hirota, T. Miyoshi, K. Kugou, C. S. Hoffman, T. Shibata, and K. Ohta, “Stepwise chromatin remodelling by a cascade of transcription initiation of non-coding RNAs,” Nature, vol. 456, no. 7218, pp. 130–134, 2008.
View at: Publisher Site | Google Scholar
C. Gulìa, S. Baldassarra, F. Signore et al., “Role of non-coding RNAs in the etiology of bladder cancer,” Genes, vol. 8, no. 11, 2017.
View at: Publisher Site | Google Scholar
C. Pop-Bica, D. Gulei, R. Cojocneanu-Petric, C. Braicu, B. Petrut, and I. Berindan-Neagoe, “Understanding the role of non-coding RNAs in bladder cancer: from dark matter to valuable therapeutic targets,” International Journal of Molecular Sciences, vol. 18, no. 7, 2017.
View at: Publisher Site | Google Scholar
M. G. K. Cumberbatch, A. P. Noon, and A. P. Noon, “Epidemiology, aetiology and screening of bladder cancer,” Translational Andrology and Urology, vol. 8, no. 1, pp. 5–11, 2019.
View at: Publisher Site | Google Scholar
R. Tang, J. Xu, B. Zhang et al., “Ferroptosis, necroptosis, and pyroptosis in anticancer immunity,” Journal of Hematology & Oncology, vol. 13, no. 1, p. 110, 2020.
View at: Publisher Site | Google Scholar
Y. Tang, C. Li, Y.-J. Zhang, and Z.-H. Wu, “Ferroptosis-Related Long Non-Coding RNA signature predicts the prognosis of Head and neck squamous cell carcinoma,” International Journal of Biological Sciences, vol. 17, no. 3, pp. 702–711, 2021.
View at: Publisher Site | Google Scholar
D. Xu, H. Zhu, C. Wang et al., “microRNA-455 targets cullin 3 to activate Nrf2 signaling and protect human osteoblasts from hydrogen peroxide,” Oncotarget, vol. 8, no. 35, pp. 59225–59234, 2017.
View at: Publisher Site | Google Scholar
S. Kabaria, D. C. Choi, A. D. Chaudhuri, M. R. Jain, H. Li, and E. Junn, “MicroRNA-7 activates Nrf2 pathway by targeting Keap1 expression,” Free Radical Biology and Medicine, vol. 89, pp. 548–556, 2015.
View at: Publisher Site | Google Scholar
X. Song and D. Long, “Nrf2 and ferroptosis: a new research direction for neurodegenerative diseases,” Frontiers in Neuroscience, vol. 14, p. 267, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Yuanshan Cui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

431

Downloads

645

Citations