Characterization of the Lipid Metabolism in Bladder Cancer to Guide Clinical Therapy

Yang, Yuan-Yuan; Hong, Sen-Yuan; Xun, Yang; Liu, Chen-Qian; Sun, Jian-Xuan; Xu, Jin-Zhou; Xu, Meng-Yao; An, Ye; He, Deng; Xia, Qi-Dong; Wang, Shao-Gang

doi:https://doi.org/10.1155/2022/7679652

Journal of Oncology

On this page

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

Special Issue

The Role of Lipid Metabolism in Tumor Progression

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 7679652 | https://doi.org/10.1155/2022/7679652

Characterization of the Lipid Metabolism in Bladder Cancer to Guide Clinical Therapy

Yuan-Yuan Yang,¹Sen-Yuan Hong,¹Yang Xun,¹Chen-Qian Liu,¹Jian-Xuan Sun,¹Jin-Zhou Xu,¹Meng-Yao Xu,¹Ye An,¹Deng He,²Qi-Dong Xia,¹and Shao-Gang Wang¹

Academic Editor: Yingkun Xu

Received22 Jul 2022

Revised02 Aug 2022

Accepted09 Aug 2022

Published12 Sept 2022

Abstract

Background. Bladder cancer is one of the most common malignancies of the urinary system with an unfavorable prognosis. More and more studies have suggested that lipid metabolism could influence the progression and treatment of tumors. However, there are few studies exploring the relationship between lipid metabolism and bladder cancer. This study aimed to explore the roles that lipid metabolism-related genes play in patients with bladder cancer. Methods. TCGA_BLCA cohort and GSE13507 cohort were included in this study, and transcriptional and somatic mutation profiles of 309 lipid metabolism-related genes were analyzed to discover the critical lipid metabolism-related genes in the incurrence and progression of bladder cancer. Furthermore, the TCGA_BLCA cohort was randomly divided into training set and validation set, and the GSE13507 cohort was served as an external independent validation set. We performed the LASSO regression and multivariate Cox regression in training set to develop a prognostic signature and further verified this signature in TCGA_BLCA validation set and GSE13507 external validation set. Finally, we systematically investigated the association between this signature and tumor microenvironment, drug response, and potential functions and then verified the differential expression status of signature genes in the protein level by immunohistochemistry. Results. A novel 6-lipidmetabolism-related gene signature was identified and validated, and this risk score model could predict the prognosis of patients with bladder cancer. In addition, the prognostic model was tightly related to immune cell infiltration and tumor mutation burden. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) showed that mTOR signaling pathway, G2M checkpoint, fatty acid metabolism, and hypoxia were enriched in patients in the high-risk score groups. Furthermore, 3 therapies specific for bladder cancer patients in different risk scores were identified. Conclusions. In conclusion, we investigated the lipid metabolism-related genes in bladder cancer through comprehensive bioinformatic analysis. A novel 6-gene signature associated with lipid metabolism for predicting the outcomes of patients with bladder cancer was conducted and validated. Furthermore, the risk score model could be utilized to indicate the choice of therapy in bladder cancer.

1. Introduction

Bladder cancer is one of the most common malignancies of the urinary system; it has the 13^th highest mortality among all cancers, and the mortality is still rising despite tremendous efforts that have been made for the treatment [1]. The progression of bladder cancer is a multistage process including environmental and genetic factors [2]. Previous studies indicate that tobacco smoking and occupational exposure to carcinogens are major factors [3]. The primary treatment for bladder cancer is transurethral resection of the bladder (TURB) accompanied by intravesical chemotherapy or immunotherapy [3]. However, the prognosis of patients still remains unfavorable. Therefore, it is essential to identify prognostic biomarkers to guide the selection of treatment for improving the curative effects.

Accumulating evidence demonstrates that clinical outcome, epigenetic status, and treatment resistance are associated with tumor metabolism [4, 5]. Metabolic reprogramming has been considered to be a new hallmark of malignant tumors [6]. Cancer cells usually require more energy to meet their biological needs than normal cells [7], while fatty acid oxidation is an important energy source for cancer cells, so lipid metabolism in cancer cells has been recognized as playing an important role in cancer progression [8].

Recent studies have shown that fatty acid metabolism has a close connection with bladder cancer [9]. Epidemiologic studies have shown that free fatty acid (FFA) level was increased in bladder cancer compared with paracancerous tissues. Cheng et al. demonstrated that inhibition of fatty metabolism is important in sustaining tumor growth through PPAR-γ-mediated pathway [10]. Besides, the molecular mechanism of drug resistance of cancer therapy might include lipid metabolism reprogramming [11]. Rysman E et al. indicated that altered lipid composition of cellular membranes could disrupt the uptake of chemotherapeutic agents and lead to chemotherapeutic resistance [12]. What is more, the tumor microenvironment (TME) plays an important role in the progression of bladder cancer, while metabolic disorders including FA metabolism changes have a crucial impact on cancer [13]. Fatty acid secreted in the microenvironment could affect the function and phenotype of infiltrating immune cells [14].

In this study, we explored the lipid metabolism-related genes in bladder cancer to conduct a model to predict patient prognosis. The prognostic risk score model independently predicted the survival outcome of bladder cancer patients. What is more, the relationship between risk score model and TME cell-infiltrating characteristics was investigated and suitable therapy treatment could be selected through the risk score model. These findings can provide a new sight into exploring the metabolic mechanism and treatment for bladder cancer.

2. Methods

2.1. Data Acquisition

We included two datasets in this study, and they are TCGA_BLCA cohort and GSE13507 cohort. The transcriptome profiles, somatic mutation profiles, and clinical information about TCGA_BLCA were downloaded from the GDC_Portal (https://portal.gdc.cancer.gov/), and the transcriptome profiles and corresponding clinical information of GSE13507 cohort were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/gds/?term=GSE13507) and served as independent external validation sets. Notably, the batch effects between TCGA_BLCA and GSE13507 cohorts were normalized by the R package “sva.”

2.2. Landscape of Lipid Metabolism-Related Genes in Bladder Cancer

We firstly investigated the role of lipid metabolism-related genes in TCGA_BLCA cohort, performed differential expression analysis between tumor and normal samples, and screened differentially expressed lipid metabolism-related genes with a threshold of FC (fold change, FC) > 1.5 and adjusted value < 0.05. Following this, we performed univariate cox regression of these differentially expressed lipid metabolism-related genes in TCGA_BLCA cohort to further identify the critical differentially expressed lipid metabolism-related genes with significant prognostic value. Subsequently, we explored the somatic mutation of these critical genes and draw the mutation atlas.

2.3. Identification of the Prognostic Signature

Having systematically summarized the role of lipid metabolism-related genes in TCGA_BLCA cohort, we would like to establish a prognostic signature by these critical genes. Thus, we first randomly split all the patients in TCGA_BLCA cohort with a ratio of 1 : 1 that one for the construction and the other for the verification. Moreover, the GSE13507 cohort was served as an external independent validation cohort. Then, we separately conducted LASSO regression to screen appropriate variables and multivariate cox regression to establish the prognostic signature in the training set, and a formula of risk score based on these lipid metabolism-related genes was established:

i represents each gene in the prognostic signature, coef (i) represents the coefficient of this gene, and exp (i) is the expression value of it. Thus, each sample in training set, validation set, and GSE13507 acquired a risk score according to this formula. Moreover, we set the medium value of the risk score in training set as the threshold, and each patient received a risk level that the higher is high risk and the lower is low risk.

2.4. Further Verification of the Prognostic Signature

We first conducted survival analysis in all three sets to confirm the prognostic value of this signature, then performed univariate Cox regression to calculate the hazard ratio of the risk score, and carried out a meta-analysis to summarize the HR of risk score in three different sets. Besides, the differences in clinical information between high-risk patients and low-risk patients were explored by the chi-square test. Also, the ROC curves of the risk score in three sets were plotted and the area under the curves was calculated. Furthermore, we combined the indicator of tumor mutation burden and our risk score to predict the survival of patients together and explored the difference in TMB between high-/low-risk patients and the correlation between TMB and risk score. Finally, we further investigated the mutation differences between high-risk patients and low-risk patients in all TCGA_BLCA patients.

2.5. Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA)

Having constructed and verified the prognostic signature, we wonder about the further potential mechanisms behind the risk score. Thus, we separately conducted the GSEA and GSVA in all TCGA_BLCA patients by the R package “clusterProfiler” and “GSVA.” The gene ontology (GO) gene sets, KEGG gene sets, Hallmarks gene sets, metabolism-related gene sets, and cell death-related gene sets were used to do the corresponding analysis.

2.6. Tumor Microenvironment, Drug Response, and Immunohistochemistry

Notably, we carried out eight different algorithms to quantify the immune cells or immune-related function of each patient according to its transcriptome profile, and they were XCELL, TIMER, QUANTISEQ, MCPCOUNTER, EPIC, CIBERSORT-ABS, CIBERSORT, and ssGSEA. Then, both differential infiltration analysis and correlation tests were conducted. Following this, we separately estimate the drug response of each patient to the commonly used drugs including cisplatin, gemcitabine, and others small molecule drugs by R package “proPhetics.” Also, the TIDE score was compared between high-risk patients and low-risk patients to predict the drug response to immunotherapy. Finally, we searched the immunohistochemistry profiles of these signature genes in the Human Protein Atlas database (HPA, https://www.proteinatlas.org/) and compared and verified the differential expression status in the protein expressed level.

2.7. Statistical Analysis

All the data were processed or analyzed by the R program version 4.1.1 and Microsoft Office Excel. was regarded with significant statistical differences.

3. Results

3.1. Identification of 18 Vital Differentially Expressed Lipid Metabolism-Related Genes in Bladder Cancer

Expression data of 309 lipid metabolism-related genes (LMRGs) were collected from the GEO and TCGA cohorts. The flow chart of this research is presented in Figure 1. 89 differently expressed LMRGs were found between normal and bladder cancer tissues, with 57 genes upregulated and 32 genes downregulated in cancer samples when the cutoff was set to FC (fold change, FC) > 1.5 and FDR<0.05 (Figures 2(a) and 2(b)). Univariate Cox regression analysis was conducted on 89 differently expressed LMRGs in TCGA_BLCA cohort. A total of 18 genes with prognostic value were identified with a value < 0.05 (Figure 2(c)). The somatic mutation profile of 18 LMRGs associated with prognosis was first summarized. A total of 54 of 412 bladder cancer samples experienced mutations of LMRGs, with a frequency of 13.11%. ACOX2, SLC27A2, and ACLY had the highest mutation frequency (Figure 2(d)). Further analyses demonstrated a mutation co-occurrence relationship between ACSF2 and PTGIS, ACOX2 and GPX1, SLC27A2 and CYP1B1, and ACLY and DHCR24 (Figure 2(e)).

(a)

(b)

(c)

(d)

(e)

Figure 2

Identification of 18 vital differentially expressed lipid metabolism-related genes (LMRGs) in bladder cancer. (a) Volcano map of the expression patterns of 89 LMRGs in bladder cancer, with 57 genes upregulated and 32 genes downregulated. (b) Heat map showing the differently expressed LMRGs between normal and tumor tissues. False discovery rate (FDR) < 0.05. (c) Forrest plot of 18 LMRGs related to prognosis based on univariate analysis. (d) The mutation frequency of 18 LMRGs in 412 patients with bladder cancer from the TCGA cohort. (e) The mutation co-occurrence and exclusion analyses for 18 LMRGs.

3.2. Construction and Verification of the Prognostic Index

To construct a prognostic index of lipid metabolism-related genes, we obtained 404 samples of bladder cancer from TCGA database and divided them into two groups as training cohort (N = 204) and testing cohort (N = 200). The basic characteristics of the patients included are shown in Table 1. Then, the least absolute shrinkage and selection operator (LASSO) Cox regression analysis was conducted to narrow the number of genes. Finally, six genes (CYP1B1, ACOT13, NUDT19, SCD, IL4I1, and DECR1) were used for the construction of prognostic risk score model (Figures 3(a) and 3(b)). Then, six genes were shifted from them through multivariate Cox regression (Figure 3(d)). The coefficient of these genes is displayed in Table 2. The chord diagram showed the correlation between the six genes filtered by LASSO Cox regression analysis (Figure 3(c)). To estimate the effectiveness of this model, we also collected 165 samples from GEO database (GSE13507) as validation. We then calculated the risk score of samples and divided the patients into high- and low-risk score groups according to median risk score acquired from training cohort. We also found that patients in the high-risk score group had the worse survival status than those in the low-risk group (Figures 3(f)–3(k)). We found that clinical characteristics such as stage and status were positively related to high-risk group (Figure 3(e)). The Kaplan–Meier survival analyses indicated that patients in the high-risk score group had worse survival outcome (Figures 3(l)–3(n)). The area under curve (AUC) showed the effectiveness of the prognostic index (Figures 3(o)–3(q)). Furthermore, we conducted meta-analysis based on the three cohorts to demonstrate that the risk score was in good validity (HR = 1.49, 95% CI = 1.15–1.94, p = 0.01) (Figure 3(r)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

(p)

(q)

(r)

Figure 3

Construction and verification of the prognostic index. (a) LASSO coefficients of the 18 lipid metabolism-related genes (LMRGs). (b) Identification of genes for development of prognostic risk score model. (c) Chord diagram showed the correlation between the 6 LMRGs selected through LASSO Cox regression analysis. (d) Hazard ratio of each gene after multivariate Cox regression analysis. (e) The distribution of clinical characteristics of patients in high- and low-risk score groups. (f–k) Risk scores of patients in the training cohort (TCGA), testing cohort (TCGA), and verification cohort (GEO). (l-n) Survival status of each patient in the training cohort, testing cohort, and verification cohort. (k-m) Kaplan–Meier survival analyses of patients in the training cohort (, log-rank test), testing cohort (, log-rank test), and verification cohort (, log-rank test). (o-q) AUC curve of patients in the training cohort, testing cohort, and verification cohort. (r) Meta-analysis of the training cohort, testing cohort, and verification cohort, and the results showed that the risk score was in good validity (HR = 1.49, 95% CI = 1.15–1.94, ).

3.3. The Relationship between Tumor Mutation Burden (TMB) and Risk Score

We compared the tumor mutation burden (TMB) in high- and low-risk score patients and found the TMB of patients in the low-risk score group was higher (Figures 4(a) and 4(b)). The mutation spectrum of patients with high-risk score (Figure 4(c)) and low-risk score (Figure 4(d)) was mapped, and 14 significantly mutant genes including (ATAD5, PIK3CA, FGFR3, HUWE1, SPTAN1, GON4L, ALMS1, RELN, STAG2, AHNAK, MED13, UTRN, C2orf16, and ATR) were obtained (Table 3). Survival curves demonstrated that the prognosis of patients with high TMB was better (Figure 4(e)); furthermore, patients with high TMB and low-risk score had the best prognosis (Figure 4(f)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4

Relationship between tumor mutation burden (TMB) and risk score. (a, b) Comparison of TMB between patients in the high- and low-risk score groups, and high tumor mutation burden was associated with low-risk score. (c, d) Mutation spectrum of high-risk and low-risk patients. (e) Survival analyses for low- and high-TMB patient groups using the Kaplan–Meier curves (, log-rank test). (f) Survival analyses for four groups of TMB and risk score using the Kaplan–Meier curves. The high-TMB and low-risk score groups indicated better overall survival than the other three groups (, log-rank test).

3.4. The Relationship between Risk Score and Tumor Microenvironment (TME) in Bladder Cancer

We estimated the relationship between risk score and the immune checkpoints (ICBs), and the results indicated that NRP1, CD44, CD276, and TNFSF9 were significantly higher expressed in high-risk score patients (Figure 5(a)). Most ICBs were highly expressed in low-risk score groups (Figure 5(b)). Immune-related function analysis demonstrated that the expression of HLA, inflammation-promoting factors, and cytolytic activity was lower in high-risk score patients compared with low-risk score patients (Figure 5(c)). Then, we explored the composition of immune infiltration cells through seven methods. We found that CD8+T, Treg, and NK-activated cells were higher in the low-risk score group compared with the high-risk group (Figures 5(d) and 5(e)). In addition, CD8+T cell was negatively associated with the risk score (Figure 5(f)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5

Relationship between risk score and tumor microenvironment (TME) in bladder cancer. (a) Immune checkpoint expression in high- and low-risk score patients (; ; ). (b) Correlation between immune checkpoint and six genes (; ; ). (c) Immune-related function of high- and low-risk score patients. (d) Correlation of risk score and immune cell infiltration detected by seven different methods. (e) Immune cell infiltration of patients in high- and low-risk score groups; CD8+T, Treg, and NK-activated cells were higher in low-risk score group compared with high-risk group. (f) Correlation between CD8+T-cell infiltration and risk scores.

3.5. Gene Set Variation Analysis (GSVA) and Gene Set Enrichment Analysis (GSEA)

We performed GSEA and GSVA to verify the correlation between risk score and pathways involved in the formation of bladder cancer. We conducted GSVA using Hallmark gene sets, and the results showed that risk score was positively associated with 31 hallmark pathways including mTOR signaling pathway, G2M checkpoint, fatty acid metabolism, and hypoxia (Figure 6(a)). We performed GSEA in different risk score groups using the Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways (Figures 6(b) and 6(c)). Besides, we confirmed that risk score was associated with pathways in cell death including autophagy and ferroptosis (Figure 6(d)). Also, we performed metabolism analysis and discovered risk score was positively related to fat-soluble vitamins (Figure 6(e)).

(a)

(b)

(c)

(d)

(e)

3.6. Response to Chemotherapy and Immunotherapy for Bladder Cancer of High- and Low-Risk Score Patients

We estimated the response to chemotherapy drugs in different risk score groups and found that patients in the low-risk score were more sensitive to methotrexate sensitivity; furthermore, metformin could be used in bladder cancer patients of low-risk score (Figure 7(a) and 7(b)). We also conducted immunotherapy analysis between the high- and low-risk score groups and found that anti-CTLA4 and anti-PD-1 therapy could make a difference in therapeutic effect between two groups no matter whether used alone or in combination, and patients in the high-risk score groups showed better sensitivity compared with those in the low-risk score groups (Figure 7(c)).

(a)

(b)

(c)

3.7. Immunohistochemistry (IHC) Verification of ACOT13, CYP1, DECR1, IL4I1, and SCD

We obtained the results of immunohistochemical staining of ACOT13, CYP1, DECR1, IL4I1, and SCD in both normal tissues and bladder cancer tissues and found that the expression of CYP1 was higher in normal tissues, while other genes were higher in tumor tissues than in normal tissues (Figure 8).

(a)

(b)

(c)

(d)

(e)

4. Discussion

Bladder cancer is a common malignant tumor with high rates of recurrence [15]. The treatments of bladder cancer have advanced a lot during the past decade; however, the high morbidity and mortality remain unchanged [16]. Effective therapies and prognostic markers still need to be identified.

Accumulating studies have indicated that lipid metabolism is a crucial step in metabolic reprogramming, while metabolic reprogramming is a new hallmark of malignant tumors [17]. Despite the importance of lipid metabolism in bladder cancer, few studies were conducted to explore the association between lipid metabolism and bladder cancer.

In this study, we identified the potential mechanism and prognostic value of lipid metabolism-related genes in bladder cancer via bioinformatic analysis. A 6-gene prognostic risk model was constructed by LASSO Cox analysis. In addition, these lipid metabolism-relatedgene-based signatures were tightly associated with TNM stage, T stage, N stage, and status. The patients with low-risk scores were found better outcomes than those with high-risk scores. Our findings showed that this risk model was an independent prognostic prediction for survival and was tightly associated with tumor mutation burden (TMB) and tumor microenvironment (TME).

Recently, more and more studies have suggested that the progression and prognosis of bladder cancer were tightly related to immune cell infiltration [18]. Furthermore, the functions of immune cell could be influenced by metabolic reprogramming [19]. Therefore, there must exist a relationship between risk scores and immune cell infiltration. Our results showed that patients in the high-risk score group had higher expression of M0 macrophages, while Tregs and CD8+T cells were upregulated in the low-risk group, demonstrating a differential infiltration pattern between the subgroups. In addition, patients in high-risk groups had higher expression of NRP1, CD44, CD276, and TNFSF9 than those in low-risk groups. These results indicated that the unfavorable prognosis of patients in the high-risk score groups might be due to the immunosuppressive environment and elevated expression of immune checkpoint genes [20]. We also found that tumor mutation burden was negatively associated with the risk score, and the prognosis of patient with high tumor mutation burden and low score was the best. It might be due to that tumor mutations could generate immunogenic neoantigens, thus leading to immune checkpoint blockade [21].

GSEA and GSVA of hallmarks suggested that G2M checkpoint, EMT pathway, and hypoxia were highly expressed in the high-risk score groups. MT Dillon et al. suggested that inhibition of G2M checkpoint could slow down the progression of tumors [22]. Epithelial-mesenchymal transition (EMT) is a process of epithelial cells acquiring mesenchymal features. It is associated with tumor initiation, invasion, metastasis, and resistance to therapy [23]. Hypoxia has emerged as a crucial factor in tumor pathophysiology, and microenvironment promotes altered cellular metabolism including lipid metabolism. Furthermore, reports suggested that hypoxia could trigger EMT in bladder cancer [24]. The pathways of autophagy and ferroptosis were also involved in patients of high-risk score groups, and Enyong Dai et al. suggested that autophagy-dependent ferroptosis could drive tumor-associated macrophage polarization [25]. The poorer prognosis of patients in the high-risk score groups might be tightly associated with mechanisms above.

The risk signature was constructed with ACOT13, CYP1, DECR1, IL4I1, NUDT19, and SCD ACOT13, and a member of acy1-CoA thioesterase (ACOT) enzymes can catalyze hydrolysis of fatty acyl-CoA into free fatty acids. It is usually enriched in oxidative tissues and tightly related to mitochondria [26]. Previous studies found that it was associated with many diseases, including lung cancer, pheochromocytomas, and paragangliomas [27]. CYP1 enzymes could catalyze the metabolic activation of procarcinogens and deactivation of certain anticancer drugs. Inhibition of CYP1 is an effective approach for chemoprevention, and many studies have suggested that inhibitors and prodrug target CYP1 are promising anticancer strategies [28]. DECR1 is an auxiliary enzyme of beta-oxidation, and it participates in redox homeostasis by controlling the balance between saturated and unsaturated phospholipids. Deletion of DECR1 can impair lipid metabolism and reduce tumor growth; therefore, DECR1 is important in the progression of tumor growth and treatment resistance [29]. NUDT19 has been identified to promote the proliferation, migration, and EMT process of tumor via mTORC1/P70S6K signaling pathway [30]. IL4I1 frequently associates with AHR (aryl hydrocarbon receptor) activity and activates the AHR through the generation of indole metabolites and kynurenic acid. In summary, it associates with reduced survival in patients with tumor and enhances the progression of tumor [31]. The previous study indicated that upregulation of SCD could proliferate cancer cells in a lipid-depleted environment for it could synthesize monounsaturated fatty acids. Decreased tumor SCD activity could slow tumor growth [32].

Furthermore, we identified the molecule drugs highly related to lipid metabolism genes for the treatment of bladder cancer. Methotrexate, a well-established antimetabolite, has been used separately or in combination for antitumoral activity for a while [33], and we found it suitable to treat patients with low-risk score. Excepting chemotherapy, immunotherapy is another important treatment for bladder cancer; in our study, we found that anti-CTLA4 and anti-PD1 were sensitive to patients with high-risk score of bladder cancer no matter whether used separately or in combination. These two drugs are immune checkpoint inhibition and have been licensed for the treatment of bladder cancer [34]. Metformin, often used for diabetes, is known to induce apoptosis in many types of cancers and has the feasibility as a drug repositioning used for the treatment of bladder cancer [35]. Our study indicated that patients in the low-risk score groups were suitable for the treatment of metformin compared with patients in the high-risk score groups.

Our study constructed and validated a prognostic signature model based on lipid metabolism genes, which could predict the prognosis of patients with bladder cancer well and guide the treatment for patients with bladder cancer. Our study also has limitations, and the main limitation of the study is that we do not have experimental studies in vivo and in vitro. Further studies would be conducted to validate what roles the lipid metabolism-related genes play in bladder cancer.

5. Conclusions

In conclusion, we investigated the lipid metabolism-related genes in bladder cancer through comprehensive bioinformatic analysis. A novel 6-gene signature associated with lipid metabolism for predicting the outcomes of patients with bladder cancer was conducted and validated. Furthermore, the risk score model could be utilized to indicate the choice of therapy in bladder cancer.

Data Availability

The source data of this study were derived from the public repositories, as indicated in the section of “Methods” of the manuscript. All data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

YYY, HD, XQD, and WSG designed the study, analyzed and interpreted the data, drafted the manuscript, and critically revised the manuscript; XQD and SJX involved in statistical analysis; YYY, XQD, XY, LCQ, SJX, HD, and WSG designed methodology; HD, XQD, and WSG administered project; all authors wrote the original draft; XQD and WSG wrote, reviewed, and edited the manuscript; and HD, XQD, and WSG are co-correspondence authors and contributed equally to this work.

Acknowledgments

The authors thank all the R programming package developer. This study was supported by the National Natural Science Foundation of China (81974092).

References

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
J. A. Witjes, H. M. Bruins, R. Cathomas et al., “European association of urology guidelines on muscle-invasive and metastatic bladder cancer: summary of the 2020 guidelines,” European Urology, vol. 79, no. 1, pp. 82–104, 2021.
View at: Publisher Site | Google Scholar
M. Babjuk, M. Burger, E. M. Compérat et al., “European association of urology guidelines on non-muscle-invasive bladder cancer (tat1 and carcinoma in situ) - 2019 update,” European Urology, vol. 76, no. 5, pp. 639–657, 2019.
View at: Publisher Site | Google Scholar
J. M. Karasinska, J. T. Topham, S. E. Kalloger et al., “Altered gene expression along the glycolysis-cholesterol synthesis axis is associated with outcome in pancreatic cancer,” Clinical Cancer Research, vol. 26, no. 1, pp. 135–146, 2020.
View at: Publisher Site | Google Scholar
S. Cristea, G. L. Coles, D. Hornburg et al., “The MEK5-ERK5 Kinase axis controls lipid metabolism in small-cell lung cancer,” Cancer Research, vol. 80, no. 6, pp. 1293–1303, 2020.
View at: Publisher Site | Google Scholar
C. Thakur and F. Chen, “Connections between metabolism and epigenetics in cancers,” Seminars in Cancer Biology, vol. 57, pp. 52–58, 2019.
View at: Publisher Site | Google Scholar
S. Beloribi-Djefaflia, S. Vasseur, and F. Guillaumond, “Lipid metabolic reprogramming in cancer cells,” Oncogenesis, vol. 5, no. 1, Article ID e189, 2016.
View at: Publisher Site | Google Scholar
H. Iwamoto, M. Abe, Y. Yang et al., “Cancer lipid metabolism confers antiangiogenic drug resistance,” Cell Metabolism, vol. 28, no. 1, pp. 104–117, 2018.
View at: Publisher Site | Google Scholar
F. Massari, C. Ciccarese, M. Santoni et al., “Metabolic phenotype of bladder cancer,” Cancer Treatment Reviews, vol. 45, pp. 46–57, 2016.
View at: Publisher Site | Google Scholar
S. Cheng, G. Wang, Y. Wang et al., “Fatty acid oxidation inhibitor etomoxir suppresses tumor progression and induces cell cycle arrest via PPARγ-mediated pathway in bladder cancer,” Clinical Science, vol. 133, no. 15, pp. 1745–1758, 2019.
View at: Publisher Site | Google Scholar
Y. Cao, “Adipocyte and lipid metabolism in cancer drug resistance,” Journal of Clinic Investigation, vol. 129, no. 8, pp. 3006–3017, 2019.
View at: Publisher Site | Google Scholar
A. J. Hoy, S. R. Nagarajan, and L. M. Butler, “Tumour fatty acid metabolism in the context of therapy resistance and obesity,” Nature Reviews Cancer, vol. 21, no. 12, pp. 753–766, 2021.
View at: Publisher Site | Google Scholar
I. Elia and M. C. Haigis, “Metabolites and the tumour microenvironment: from cellular mechanisms to systemic metabolism,” Nature Metabolism, vol. 3, no. 1, pp. 21–32, 2021.
View at: Publisher Site | Google Scholar
K. Man, A. Kallies, and A. Vasanthakumar, “Resident and migratory adipose immune cells control systemic metabolism and thermogenesis,” Cellular Molecular Immunol, vol. 19, no. 3, pp. 421–431, 2022.
View at: Publisher Site | Google Scholar
V. G. Patel, W. K. Oh, and M. D. Galsky, “Treatment of muscle-invasive and advanced bladder cancer in 2020,” CA: A Cancer Journal for Clinicians, vol. 70, no. 5, pp. 404–423, 2020.
View at: Publisher Site | Google Scholar
I. H. Kim and H. J. Lee, “Perioperative systemic treatment for muscle-invasive bladder cancer: current evidence and future perspectives,” International Journal of Molecular Sciences, vol. 22, no. 13, Article ID 7201, 2021.
View at: Publisher Site | Google Scholar
J. Liu, Y. Peng, L. Shi et al., “Skp2 dictates cell cycle-dependent metabolic oscillation between glycolysis and TCA cycle,” Cell Research, vol. 31, no. 1, pp. 80–93, 2021.
View at: Publisher Site | Google Scholar
S. Shi, T. Ma, and Y. Xi, “Characterization of the immune cell infiltration landscape in bladder cancer to aid immunotherapy,” Archives of Biochemistry and Biophysics, Article ID 108950, p. 708, 2021.
View at: Publisher Site | Google Scholar
S. M. Morrissey, F. Zhang, C. Ding et al., “Tumor-derived exosomes drive immunosuppressive macrophages in a pre-metastatic niche through glycolytic dominant metabolic reprogramming,” Cell Metabolism, vol. 33, no. 10, pp. 2040–2058.e10, 2021.
View at: Publisher Site | Google Scholar
X. Chen, R. Xu, D. He et al., “CD8⁺ T effector and immune checkpoint signatures predict prognosis and responsiveness to immunotherapy in bladder cancer,” Oncogene, vol. 40, no. 43, pp. 6223–6234, 2021.
View at: Publisher Site | Google Scholar
J. Pei, Y. Li, T. Su et al., “Identification and validation of an immunological expression-based prognostic signature in breast cancer,” Frontiers in Genetics, vol. 11, p. 912, 2020.
View at: Publisher Site | Google Scholar
M. T. Dillon, J. S. Good, and K. J. Harrington, “Selective targeting of the G2/M cell cycle checkpoint to improve the therapeutic index of radiotherapy,” Clinical Oncology, vol. 26, no. 5, pp. 257–65, 2014.
View at: Publisher Site | Google Scholar
I. Pastushenko and C. Blanpain, “EMT transition states during tumor progression and metastasis,” Trends in Cell Biology, vol. 29, no. 3, pp. 212–226, 2019.
View at: Publisher Site | Google Scholar
J. P. Joseph, M. K. Harishankar, A. A. Pillai, and A. Devi, “Hypoxia induced EMT: a review on the mechanism of tumor progression and metastasis in OSCC,” Oral Oncology, vol. 80, pp. 23–32, 2018.
View at: Publisher Site | Google Scholar
E. Dai, L. Han, J. Liu et al., “Autophagy-dependent ferroptosis drives tumor-associated macrophage polarization via release and uptake of oncogenic KRAS protein,” Autophagy, vol. 16, no. 11, pp. 2069–2083, 2020.
View at: Publisher Site | Google Scholar
N. Du, D. Dong, L. Sun et al., “Identification of ACOT13 and PTGER2 as novel candidate genes of autosomal dominant polycystic kidney disease through whole exome sequencing,” European Journal of Medical Research, vol. 26, no. 1, p. 142, 2021.
View at: Publisher Site | Google Scholar
P. L. M. Dahia, “Novel hereditary forms of pheochromocytomas and paragangliomas,” Front Horm Res, vol. 41, pp. 79–91, 2013.
View at: Publisher Site | Google Scholar
J. Cui and S. Li, “Inhibitors and prodrugs targeting CYP1: a novel approach in cancer prevention and therapy,” Current Medicinal Chemistry, vol. 21, no. 5, pp. 519–52, 2014.
View at: Publisher Site | Google Scholar
A. Blomme, C. A. Ford, E. Mui et al., “2, 4-dienoyl-CoA reductase regulates lipid homeostasis in treatment-resistant prostate cancer,” Nature Communication, vol. 11, no. 1, p. 2508, 2020.
View at: Publisher Site | Google Scholar
C. Lan, Y. Wang, X. Su, J. Lu, and S. Ma, “LncRNA LINC00958 Activates mTORC1/P70S6K signalling pathway to promote epithelial-mesenchymal transition process in the hepatocellular carcinoma,” Cancer Investigation, vol. 39, no. 6-7, pp. 539–549, 2021.
View at: Publisher Site | Google Scholar
A. Sadik, L. F. Somarribas Patterson, S. Öztürk et al., “IL4I1 is a metabolic immune checkpoint that activates the AHR and promotes tumor progression,” Cell, vol. 182, no. 5, pp. 1252–1270, 2020.
View at: Publisher Site | Google Scholar
E. C. Lien, A. M. Westermark, Y. Zhang et al., “Low glycaemic diets alter lipid metabolism to influence tumour growth,” Nature, vol. 599, no. 7884, pp. 302–307, 2021.
View at: Publisher Site | Google Scholar
S. Kern, I. Truebenbach, M. Höhn et al., “Combined antitumoral effects of pretubulysin and methotrexate,” Pharmacol Res Perspect, vol. 7, no. 1, Article ID e00460, 2019.
View at: Publisher Site | Google Scholar
L. Dyck and K. H. G. Mills, “Immune checkpoints and their inhibition in cancer and infectious diseases,” European Journal of Immunology, vol. 47, no. 5, pp. 765–779, 2017.
View at: Publisher Site | Google Scholar
J. H. Jang, E. G. Sung, I. H. Song, T. J. Lee, and J. Y. Kim, “Metformin induces caspase-dependent and caspase-independent apoptosis in human bladder cancer T24 cells,” Anticancer Drugs, vol. 31, no. 7, pp. 655–662, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Yuan-Yuan Yang 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