- 1Department of Urology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
- 2Guangdong Provincial Clinical Research Center for Urological Diseases, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
Background: Some studies have revealed a close relationship between metabolism-related genes and the prognosis of bladder cancer. However, the relationship between metabolism-related long non-coding RNAs (lncRNA) regulating the expression of genetic material and bladder cancer is still blank. From this, we developed and validated a prognostic model based on metabolism-associated lncRNA to analyze the prognosis of bladder cancer.
Methods: Gene expression, lncRNA sequencing data, and related clinical information were extracted from The Cancer Genome Atlas (TCGA). And we downloaded metabolism-related gene sets from the human metabolism database. Differential expression analysis is used to screen differentially expressed metabolism-related genes and lncRNAs between tumors and paracancer tissues. We then obtained metabolism-related lncRNAs associated with prognosis by correlational analyses, univariate Cox analysis, and logistic least absolute shrinkage and selection operator (LASSO) regression. A risk scoring model is constructed based on the regression coefficient corresponding to lncRNA calculated by multivariate Cox analysis. According to the median risk score, patients were divided into a high-risk group and a low-risk group. Then, we developed and evaluated a nomogram including risk scores and Clinical baseline data to predict the prognosis. Furthermore, we performed gene-set enrichment analysis (GSEA) to explore the role of these metabolism-related lncRNAs in the prognosis of bladder cancer.
Results: By analyzing the extracted data, our research screened out 12 metabolism-related lncRNAs. There are significant differences in survival between high and low-risk groups divided by the median risk scoring model, and the low-risk group has a more favorable prognosis than the high-risk group. Univariate and multivariate Cox regression analysis showed that the risk score was closely related to the prognosis of bladder cancer. Then we established a nomogram based on multivariate analysis. After evaluation, the modified model has good predictive efficiency and clinical application value. Furthermore, the GSEA showed that these lncRNAs affected bladder cancer prognosis through multiple links.
Conclusions: A predictive model was established and validated based on 12 metabolism-related lncRNAs and clinical information, and we found these lncRNA affected bladder cancer prognosis through multiple links.
Introduction
Bladder cancer, as a common tumor of the urinary system, has always been the focus of research. Recent studies have shown that the metabolism of glycogen, lipid, amino acid, and other substances is closely related to the diagnosis and prognosis of tumors (1, 2). For example, in cancer cells, blocking the supply of lipids can have serious effects on bioenergetics, membrane biosynthesis, and intracellular signaling processes (3). And the specificity of glycogen metabolism in tumor cells, the Warburg effect, has also been widely studied and discussed (4–7). These studies have opened up potential therapeutic targets for tumors. And the tumor metabolism-based therapy treatment has led to improved cancer outcomes (2, 8). Despite the success of previous studies in improving the prognosis of bladder cancer, there were 213,000 deaths from bladder cancer worldwide in 2020 (9). And little reformation in treatment has occurred in recent years. Therefore, we still need to explore new biomarkers to aid in diagnosis and treatment.
The lncRNAs play a crucial role in transcriptional regulation, epigenetic gene regulation, and disease, a kind of non-coding RNAs greater than 200 nucleotides, with no protein-coding function (10, 11). With the deepening of research, the relationship between lncRNA and progression, metastasis, and prognosis of various tumors has been found, such as the digestive, nervous, bone, and urinary system (12–17).
For bladder cancer, immune-associated and autophagy-associated lncRNAs were found as a marker of early diagnostic and prognostic treatment (18–20). Recent studies show that lncRNAs play important roles in glucose and lipid (21, 22).
Although we know the nonnegligible fact that metabolism-related genes and lncRNAs are involved in tumor prognosis. At present, the research on metabolism-related lncRNA is not abundant, which is the starting point of this dissertation. It is of great underlying value to study metabolism-related lncRNA to accurately predict the prognosis of bladder cancer. Hence, to address the clinical values, we tried to establish a nomogram containing metabolism-related lncRNAs and clinical data and explore the possible functions of metabolism-related lncRNAs beings through the GSEA.
Methods
Data Collection and Processing
The RNA-Seq, lncRNA sequencing data, and related clinical information of bladder cancer were extracted from TCGA (https://portal.gdc.cancer.gov). And metabolism-related gene sets were downloaded from the GSEA database (https://www.gsea-msigdb.org/gsea/index.jsp). Taking | log2FC | > 0.5 and P < 0.05 as standards, the “edgeR” package was used to construct the volcano map and distinguish the differentially expressed metabolism-related gene and lncRNAs in tumors and paracancer tissues. The correlation between metabolism-related genes and lncRNAs is measured by the Pearson correlation coefficient. Metabolism-related gene and lncRNAs whose correlation coefficients satisfy | R2 | > 0.5 and P < 0.05 are considered related and used for further analysis.
Identification of Prognostic Related LncRNAs and Construction and Validation Risk Score Model
Patients with a follow-up time no less than 30 days were divided into the training group and the validation group in a ratio of 2 to 1. Univariate Cox regression analysis was used to analyze clinical survival data and lncRNAs. These lncRNAs associated with prognosis were further screened by LASSO regression analysis. The n selected lncRNAs were analyzed by multivariate Cox regression analysis, and the risk scoring model was based on coefficients. The AUC was used to evaluate the scoring model.
We can calculate the risk score for each patient using the scoring model. The patients were divided into high and low-risk groups based on the median risk score. The Kaplan-Meier survival curve was plotted to show the difference in prognosis between the two groups. The risk score and survival status of each sample were presented using risk curves and scatter plots. Heat maps were used to show the expression state of selected lncRNAs.
Construction and Aessement Nomogram
Risk scores, age, gender, and TNM were included in univariate and multivariate Cox analysis to obtain variables associated with prognosis. Then, a nomogram was established based on the results of multivariate Cox analysis. We evaluated the nomogram by calculating the C-index and plotting the area under the receiver operating characteristic curve (AUC), calibration plot, and decision curve analysis (DCA).
Gene-Set Enrichment Analysis
The differentially expressed metabolism-related genes in the high-risk group and low-risk groups were analyzed for gene enrichment. To explore what biological functions or pathways these expressed different genes might be involved in.
Statistical Analysis
All statistical computations were conducted using the R software, version 4.0.2 (The R Foundation for Statistical Computing, Vienna, Austria http://www.R-project.org). And P < 0.05 was considered statistically significant.
Results
Differentially Expressed Genes and LncRNAs
We grabbed RNA-seqs, lncRNA sequences, and clinical data of bladder cancer tissues and paracancer tissues data from the TCGA database. And 937 metabolism-related gene sets were downloaded from the GSEA database. 169 differentially expressed metabolism-related genes (Additional File: Table 1) and 2301 differentially expressed lncRNAs (Additional File: Table 2) were obtained by differential expression analysis (Figure 1).
Figure1 Identification of differentially expressed Metabolism-related genes [volcano plots (A)] and lncRNAs [volcano plots (B)].
Identification of Metabolism-Related LncRNAs Associated With Prognosis
We identified 105 metabolism-related lncRNAs by correlation Pearson analysis, as shown in Additional File Table 3. And univariate Cox regression was performed, we screened out 20 metabolism-related lncRNAs associated with prognosis. To control overfitting, LASSO regression analysis was used, and finally, 12 metabolism-related lncRNAs associated with prognosis were obtained. These selected lncRNAs were CIRBP.AS1, AC018653.3, AL357033.4, LINC02004, DUXAP8, AC010331.1, PWAR6, AC025575.2, AL355353.1, AL731567.1, AC074117.1 and AC073335.2. Among them, LINC02004 is unfavorable to the prognosis of bladder cancer (Figure 2).
Figure 2 (A) Univariate Cox regression analysis was used to screen prognostic-related metabolism-related lncRNAs. (B, C) The screened 20 prognostic-related metabolism-related lncRNAs were incorporated into the Lasso regression model and penalties were used for controlling the overfitting effects of the model.
Construction and Validation Risk Score Model
393 patients were divided into training group (262) and validation group (131). In the training group, the risk score model was determined by multivariate cox regression analysis to calculate the coefficients for selected lncRNAs.
The median risk score of patients calculated based on this risk score model was used to be divided into high and low-risk groups. In the training group and validation group, scores and gene expression, and survival prognosis for high and low risk are shown (Figure 3). The high-risk group had an unfavorable prognosis than the low-risk group (P < 0.05). And in both the training group and validation group, the AUC was greater than 0.71 (Figure 3). Furthermore, in Additional File 1, the good predictive effectiveness of the risk scoring model in each subgroup (age, sex, TNM stage) was presented.
Figure 3 The performance of risk score model based on 12 metabolism-related lncRNAs in training (A–E) and validation groups (F–J). Kaplan–Meier curves showed that patients with high-risk scores had a worse overall survival than patients in the low-risk group. (A, F) Patients with low-risk scores generally had longer survival years, and lower mortality than the other patients. (B, C, G, H) Heat maps of the risk scores were based on metabolism-related lncRNAs in the training group (D) and validation group (I). ROC curves and AUC values of 3-, 5- year of risk score model in the training group (E) and validation group (J).
Construction and Aessement Nomogram
Clinical data and scores were included in the univariate and multivariate analysis (Figure 4). We constructed a nomogram based on the results of multivariate Cox analysis including sick score, age, and T stage. The c-index of the training group and validation group are 0.741 and 0.689 respectively. The ROC and DCA Curves showed satisfactory results (Figure 5). These results indicate that our nomogram is reliable and has clinical application value.
Figure 4 The result of univariate (A) and multivariate (B) Cox regression performed in the risk score and clinical data. Both in the two forests above suggested that the risk score was an independent prognostic factor.
Figure 5 The nomogram (A) of the risk score model and its predictive power and reliability in the training group (B, D, F) and validation group (C, E, G).
Gene-Set Enrichment Analysis
We performed Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis in differentially expressed metabolism-related genes between high and low-risk groups. Then we presented 9 KEGG molecular pathways in which different genes were enriched (Figure 6). Some metabolic pathways enriched into the high-risk group included galactose metabolism and sugar and nucleotide sugar metabolism. And glycine serine and threonine metabolism, glycerophospholipid metabolism, oxidative phosphorylation, and linoleic acid metabolism enriched into the low-risk group.
Figure 6 Gene-set enrichment analysis. The analysis results showed significant enrichment in the high-risk group (A–H), and in the low-risk group (I–L). And the metabolism-related genes are enriched in the high-risk group (G, H), and the low-risk group (I–L).
Discussion
Worldwide, bladder cancer is characterized by high incidence, poor prognosis, and high economic burden. How to accurately predict prognosis and explore new therapeutic targets are urgently needed in the clinic. Previous studies have found metabolism-related genes and lncRNAs affect the prognosis of bladder cancer in different aspects (1, 23). However, there is little research on the relationship between metabolism-related lncRNAs and bladder cancer. This thesis fills the research gap by identifying metabolism-related lncRNAs associated with bladder cancer, exploring their function channels, and developing a nomogram to predict prognosis individually.
Although there are few studies on the prognostic relationship between metabolism-related related lncRNAs and bladder cancer, the prediction model based on other types of lncRNAs has a good performance in predicting the prognosis of bladder cancer. In tumor immunity, an 8 immune-related LncRNA classifier for prognostic prediction, could predict prognosis and immunotherapeutic response (24). And patients with bladder cancer grouped according to identified immune-related lncRNA signature showed different immune states (25, 26). Ferroptosis, an iron-dependent form of nonapoptotic cell death has been proved to be closely related to the development of cancer (27). The ferroptosis-related lncRNA signatures can predict the prognosis of bladder cancer patients, which still needs further experimental verification in the future (28).
From the Kaplan-Meier survival curve we drew according to subgroups, we could see that this model still has good predictive performance in most subgroups age > 65, female, male, T1-2, T3-4, N0, and M0. In the remaining three subgroups, no significant differences were observed, and the insufficient sample size was considered as the main reason.
Our study developed a risk score model with good reliability. This reconfirmed that these selected metabolism-related lncRNA (CIRBP.AS1, AC018653.3, AL357033.4, LINC02004, DUXAP8, AC010331.1, PWAR6, AC025575.2, AL355353.1, AL731567.1, AC074117.1 and AC073335.2.) are dysregulated in bladder cancer.
LncRNA DUXAP8 was located on chromosome 20q11 with 2307 bp in length8 (29). DUXAP8 can regulate PTEN to alter the prognosis of bladder cancer (30). Similar to Liu’s research, DUXAP8 might be a useful lncRNAs resource for prognostic or diagnostic markers for bladder cancer (31). Recent studies have found that silencing DUXAP8 expression can inhibit cell proliferation and promote apoptosis by targeting miR-26b-5P (32). DUXAP8 is involved in the initiation of many tumors and has the potential as a therapeutic target (33). LncRNA PWAR6 regulates proliferation and migration by epigenetically silencing YAP1 in tumorigenesis of pancreatic ductal adenocarcinoma (34). One study has adopted a machine-learning Algorithm and Laplacian score feature selection, and an advanced over-sampling, they found AC074117.1 is the target of multiple cancer-related miRNAs and interacts with adjacent protein-coding genes (35). RNA is associated with prognosis in metastatic lung Adenocarcinomaand rectal adenocarcinoma (36, 37). Also according to existing research findings, AC010331.1, a favorable prognostic factor, can improve prognosis by upregulating autophagy-related gene expression (38). AC018653.3 and CIRBP.AS1 may also work through regulating Ferroptosis-related genes and N6-methyladenosine related genes respectively (39, 40). These suggest that abnormal regulation of some lncRNA may affect multiple aspects of bladder tumors cell. The rest of the lncRNAs (AL357033.4, LINC02004, AC025575.2, AL355353.1, AL731567.1, and AC073335.2) have not been studied in bladder cancer yet. We believe these lncRNAs could be a new focus.
We found through GSEA that different metabolism-related pathways were significantly enriched both in high and low-risk groups. The galactose metabolism and amino sugar and nucleotide sugar metabolism enrich in patients with high-risk scores, which might mean a poor prognosis in these patients. Consistent with our research, Overexpression of the galactose transporter and galactose binding lectin, may contribute to tumor progression (41, 42). On the contrary, serine and threonine metabolism, glycerophospholipid metabolism, oxidative phosphorylation, and linoleic acid metabolism are enriched in patients who would have a favorable prognosis. The potential value of serine and threonine metabolism and linoleic acid metabolism as a therapeutic target has been demonstrated in some studies (43, 44).
We were the first to build a nomogram based on metabolism-related lncRNA. After verification and evaluation, the model can accurately predict the prognosis of bladder cancer. Using GSEA, we initial research the potential function of 12 metabolism-related lncRNA. Admittedly, there are also some limitations to our study. First, further experiments on cells and animals are needed to verify the function of this lncRNA. Second, contained clinical data lacks some important data, such as comorbidities, therapies, smoking, and the cause of death.
Conclusion
We identified 12 metabolism-related lncRNAs associated with the prognosis of bladder cancer. And we found these lncRNA affecting bladder cancer prognosis through multiple links could be the focus of future research. From this, the first nomogram was established and validated based on metabolism-related lncRNAs and clinical information.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author Contributions
JTH, JLH, and JK conceived and designed the study, participated in the collection of data and data analysis, and drafted the manuscript. CL and ZS assisted in the design of this research and project development. HY, JL, WX, and HS analyzed the data and reviewed the article. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Science and Technology Planning Project of Guangdong Province (Grant No. 2020A1515111119) and Guangdong Provincial Clinical Research Center for Urological Diseases (Grant No.2020B1111170006). The funders played no role in design of the study, collection, analysis, and interpretation of data or in writing the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.833763/full#supplementary-material
Supplementary Table 1 | 169 metabolism-related genes differentially expressed in tumors and paracancer tissues.
Supplementary Table 2 | 2301 lncRNAs differentially expressed in tumors and paracancer tissues.
Supplementary Table 3 | 105 metabolism-related lncRNAs were identified.
Supplementary Figure 1 | The survival outcomes of bladder cancer patients with different risk scores in subgroups. (A) Age < 65, (B) Age ≥ 65, (C) Female, (D) Male, (E) T1-2, (F) T3-4, (G) N0, (H) N1-3, (I) M0, (J) M1.
References
1. Putluri N, Shojaie A, Vasu VT, Vareed SK, Nalluri S, Putluri V, et al. Metabolomic Profiling Reveals Potential Markers and Bioprocesses Altered in Bladder Cancer Progression. Cancer Res (2011) 71(24):7376–86. doi: 10.1158/0008-5472.Can-11-1154
2. Bellmunt J, de Wit R, Vaughn DJ, Fradet Y, Lee JL, Fong L, et al. Pembrolizumab as Second-Line Therapy for Advanced Urothelial Carcinoma. N Engl J Med (2017) 376(11):1015–26. doi: 10.1056/NEJMoa1613683
3. Röhrig F, Schulze A. The Multifaceted Roles of Fatty Acid Synthesis in Cancer. Nat Rev Cancer (2016) 16(11):732–49. doi: 10.1038/nrc.2016.89
4. Wang G, Wang JJ, Guan R, Du L, Gao J, Fu XL. Strategies to Target Glucose Metabolism in Tumor Microenvironment on Cancer by Flavonoids. Nutr Cancer (2017) 69(4):534–54. doi: 10.1080/01635581.2017.1295090
5. Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg Effect: The Metabolic Requirements of Cell Proliferation. Sci (New York NY) (2009) 324(5930):1029–33. doi: 10.1126/science.1160809
6. Massari F, Ciccarese C, Santoni M, Iacovelli R, Mazzucchelli R, Piva F, et al. Metabolic Phenotype of Bladder Cancer. Cancer Treat Rev (2016) 45:46–57. doi: 10.1016/j.ctrv.2016.03.005
7. Cheng S, Wang G, Wang Y, Cai L, Qian K, Ju L, et al. Fatty Acid Oxidation Inhibitor Etomoxir Suppresses Tumor Progression and Induces Cell Cycle Arrest via Pparγ-Mediated Pathway in Bladder Cancer. Clin Sci (Lond Engl: 1979) (2019) 133(15):1745–58. doi: 10.1042/cs20190587
8. Zhu J, Thompson CB. Metabolic Regulation of Cell Growth and Proliferation. Nat Rev Mol Cell Biol (2019) 20(7):436–50. doi: 10.1038/s41580-019-0123-5
9. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
10. Kumar MM, Goyal R. Lncrna as a Therapeutic Target for Angiogenesis. Curr Topics Med Chem (2017) 17(15):1750–7. doi: 10.2174/1568026617666161116144744
11. Chandra Gupta S, Nandan Tripathi Y. Potential of Long Non-Coding Rnas in Cancer Patients: From Biomarkers to Therapeutic Targets. Int J Cancer (2017) 140(9):1955–67. doi: 10.1002/ijc.30546
12. Yang Z, Li X, Yang Y, He Z, Qu X, Zhang Y. Long Noncoding Rnas in the Progression, Metastasis, and Prognosis of Osteosarcoma. Cell Death Dis (2016) 7(9):e2389. doi: 10.1038/cddis.2016.272
13. Zou W, Wang Z, Wang F, Li L, Liu R. Hu M. A Metabolism-Related 4-Lncrna Prognostic Signature and Corresponding Mechanisms in Intrahepatic Cholangiocarcinoma. BMC Cancer (2021) 21(1):608. doi: 10.1186/s12885-021-08322-5
14. Lu Y, Wang W, Liu Z, Ma J, Zhou X, Fu W. Long Non-Coding RNA Profile Study Identifies a Metabolism-Related Signature for Colorectal Cancer. Mol Med (Cambridge Mass) (2021) 27(1):83. doi: 10.1186/s10020-021-00343-x
15. Zhao J, Wang L, Wei B. Identification and Validation of an Energy Metabolism-Related Lncrna-Mrna Signature for Lower-Grade Glioma. BioMed Res Int (2020) 2020:3708231. doi: 10.1155/2020/3708231
16. Aurilia C, Donati S, Palmini G, Miglietta F, Iantomasi T, Brandi ML. The Involvement of Long Non-Coding Rnas in Bone. Int J Mol Sci (2021) 22(8):3909. doi: 10.3390/ijms22083909
17. Quan J, Pan X, Zhao L, Li Z, Dai K, Yan F, et al. Lncrna as a Diagnostic and Prognostic Biomarker in Bladder Cancer: A Systematic Review and Meta-Analysis. OncoTargets Ther (2018) 11:6415–24. doi: 10.2147/ott.S167853
18. Du Y, Wang B, Jiang X, Cao J, Yu J, Wang Y, et al. Identification and Validation of a Stromal EMT Related Lncrna Signature as a Potential Marker to Predict Bladder Cancer Outcome. Front Oncol (2021) 11:620674. doi: 10.3389/fonc.2021.620674
19. Wu J, Li W, Ning J, Yu W, Rao T, Cheng F. Long Noncoding RNA UCA1 Targets Mir-582-5p and Contributes to the Progression and Drug Resistance of Bladder Cancer Cells Through ATG7-Mediated Autophagy Inhibition. OncoTargets Ther (2019) 12:495–508. doi: 10.2147/ott.S183940
20. Eissa S, Safwat M, Matboli M, Zaghloul A, El-Sawalhi M, Shaheen A. Measurement of Urinary Level of a Specific Competing Endogenous RNA Network (FOS and RCAN Mrna/Mir-324-5p, Mir-4738-3p,/Lncrna Mir-497-HG) Enables Diagnosis of Bladder Cancer. Urol Oncol (2019) 37(4):292. doi: 10.1016/j.urolonc.2018.12.024
21. Tang J, Yan T, Bao Y, Shen C, Yu C, Zhu X, et al. Lncrna GLCC1 Promotes Colorectal Carcinogenesis and Glucose Metabolism by Stabilizing C-Myc. Nat Commun (2019) 10(1):3499. doi: 10.1038/s41467-019-11447-8
22. Zeng Y, Ren K, Zhu X, Zheng Z, Yi G. Long Noncoding Rnas: Advances in Lipid Metabolism. Adv Clin Chem (2018) 87:1–36. doi: 10.1016/bs.acc.2018.07.001
23. Zhou Q, Zhan H, Lin F, Liu Y, Yang K, Gao Q, et al. Lincrna-P21 Suppresses Glutamine Catabolism and Bladder Cancer Cell Growth Through Inhibiting Glutaminase Expression. Biosci Rep (2019) 39(4). doi: 10.1042/bsr20182372
24. Wu Y, Zhang L, He S, Guan B, He A, Yang K, et al. Identification of Immune-Related Lncrna for Predicting Prognosis and Immunotherapeutic Response in Bladder Cancer. Aging (2020) 12(22):23306–25. doi: 10.18632/aging.104115
25. Wang J, Shen C, Dong D, Zhong X, Wang Y, Yang X. Identification and Verification of an Immune-Related Lncrna Signature for Predicting the Prognosis of Patients With Bladder Cancer. Int Immunopharmacol (2021) 90:107146. doi: 10.1016/j.intimp.2020.107146
26. Hua S, Xie Z, Wang W, Wan Z, Chen M, Zhao S, et al. Identification and Validation of a Novel Immune-Related Lncrna Signature for Bladder Cancer. Front Oncol (2021) 11:704946. doi: 10.3389/fonc.2021.704946
27. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: An Iron-Dependent Form of Nonapoptotic Cell Death. Cell (2012) 149(5):1060–72. doi: 10.1016/j.cell.2012.03.042
28. Chen M, Nie Z, Li Y, Gao Y, Wen X, Cao H, et al. A New Ferroptosis-Related Lncrna Signature Predicts the Prognosis of Bladder Cancer Patients. Front Cell Dev Biol (2021) 9:699804. doi: 10.3389/fcell.2021.699804
29. Sun M, Nie FQ, Zang C, Wang Y, Hou J, Wei C, et al. The Pseudogene DUXAP8 Promotes Non-Small-Cell Lung Cancer Cell Proliferation and Invasion by Epigenetically Silencing EGR1 and RHOB. Mol Ther: J Am Soc Gene Ther (2017) 25(3):739–51. doi: 10.1016/j.ymthe.2016.12.018
30. Lin MG, Hong YK, Zhang Y, Lin BB, He XJ. Mechanism of Lncrna DUXAP8 in Promoting Proliferation of Bladder Cancer Cells by Regulating PTEN. Eur Rev Med Pharmacol Sci (2018) 22(11):3370–7. doi: 10.26355/eurrev_201806_15158
31. Jiang B, Hailong S, Yuan J, Zhao H, Xia W, Zha Z, et al. Identification of Oncogenic Long Noncoding RNA SNHG12 and DUXAP8 in Human Bladder Cancer Through a Comprehensive Profiling Analysis. Biomed Pharmacother = Biomed Pharmacother (2018) 108:500–7. doi: 10.1016/j.biopha.2018.09.025
32. Liu Y, Zhang G, Chen H, Wang H. Silencing Lncrna DUXAP8 Inhibits Lung Adenocarcinoma Progression by Targeting Mir-26b-5p. Biosci Rep (2021) 41(1). doi: 10.1042/bsr20200884
33. Xue C, Cai X, Jia J. Long Non-Coding RNA Double Homeobox a Pseudogene 8: A Novel Oncogenic Propellant in Human Cancer. Front Cell Dev Biol (2021) 9:709069. doi: 10.3389/fcell.2021.709069
34. Huang S, Li Y, Hu J, Li L, Liu Z, Guo H, et al. Lncrna PWAR6 Regulates Proliferation and Migration by Epigenetically Silencing YAP1 in Tumorigenesis of Pancreatic Ductal Adenocarcinoma. J Cell Mol Med (2021) 25(9):4275–86. doi: 10.1111/jcmm.16480
35. Zhang X, Li T, Wang J, Li J, Chen L, Liu C. Identification of Cancer-Related Long Non-Coding Rnas Using Xgboost With High Accuracy. Front Genet (2019) 10:735. doi: 10.3389/fgene.2019.00735
36. Li Z, Wang D, Yin H. A Seven Immune-Related Lncrna Signature Predicts the Survival of Patients With Colon Adenocarcinoma. Am J Trans Res (2020) 12(11):7060–78.
37. Li M, Yang B, Li X, Ren H, Zhang L, Li L, et al. Identification of Prognostic Factors Related to Super Enhancer-Regulated Cerna Network in Metastatic Lung Adenocarcinoma. Int J Gen Med (2021) 14:6261–75. doi: 10.2147/ijgm.S332317
38. Wan J, Guo C, Fang H, Xu Z, Hu Y, Luo Y. Autophagy-Related Long Non-Coding RNA is a Prognostic Indicator for Bladder Cancer. Front Oncol (2021) 11:647236. doi: 10.3389/fonc.2021.647236
39. Cui Y, Zhou Z, Chai Y, Che X, Zhang Y. Identification of a Nomogram From Ferroptosis-Related Long Noncoding Rnas Signature to Analyze Overall Survival in Patients With Bladder Cancer. J Oncol (2021) 2021:8533464. doi: 10.1155/2021/8533464
40. Wu J, Cai Y, Zhao G, Li M. A Ten N6-Methyladenosine-Related Long Non-Coding Rnas Signature Predicts Prognosis of Triple-Negative Breast Cancer. J Clin Lab Anal (2021) 35(6):e23779. doi: 10.1002/jcla.23779
41. Schoeppner HL, Raz A, Ho SB, Bresalier RS. Expression of an Endogenous Galactose-Binding Lectin Correlates With Neoplastic Progression in the Colon. Cancer (1995) 75(12):2818–26. doi: 10.1002/1097-0142(19950615)75:12<2818::aid-cncr2820751206>3.0.co;2-#
42. Kouznetsova VL, Kim E, Romm EL, Zhu A, Tsigelny IF. Recognition of Early and Late Stages of Bladder Cancer Using Metabolites and Machine Learning. Metabol: Off J Metabol Soc (2019) 15(7):94. doi: 10.1007/s11306-019-1555-9
43. Cai B, Ma L, Nong S, Wu Y, Guo X, Pu J. β-Elemene Induced Anticancer Effect in Bladder Cancer Through Upregulation of PTEN and Suppression of AKT Phosphorylation. Oncol Lett (2018) 16(5):6019–25. doi: 10.3892/ol.2018.9401
Keywords: bladder cancer, TCGA, metabolism, lncRNA, nomogram, prognosis
Citation: Hu J, Lai C, Shen Z, Yu H, Lin J, Xie W, Su H, Kong J and Han J (2022) A Prognostic Model of Bladder Cancer Based on Metabolism-Related Long Non-Coding RNAs. Front. Oncol. 12:833763. doi: 10.3389/fonc.2022.833763
Received: 12 December 2021; Accepted: 03 February 2022;
Published: 25 February 2022.
Edited by:
Christos K. Kontos, National and Kapodistrian University of Athens, GreeceCopyright © 2022 Hu, Lai, Shen, Yu, Lin, Xie, Su, Kong and Han. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jinli Han, aGFuamlubGlAbWFpbC5zeXN1LmVkdS5jbg==; Jianqiu Kong, a29uZ2pxNUBtYWlsLnN5c3UuZWR1LmNu
†These authors have contributed equally to this work and share first authorship