Dataset Information

Metabolic classification of bladder cancer based on multi-omics integrated analysis to predict patient prognosis and treatment response.

ABSTRACT:

Background

Currently, no molecular classification is established for bladder cancer based on metabolic characteristics. Therefore, we conducted a comprehensive analysis of bladder cancer metabolism-related genes using multiple publicly available datasets and aimed to identify subtypes according to distinctive metabolic characteristics.

Methods

RNA-sequencing data of The Cancer Genome Atlas were subjected to non-negative matrix fractionation to classify bladder cancer according to metabolism-related gene expression; Gene Expression Omnibus and ArrayExpress datasets were used as validation cohorts. The sensitivity of metabolic types to predicted immunotherapy and chemotherapy was assessed. Kaplan-Meier curves were plotted to assess patient survival. Differentially expressed genes between subtypes were identified using edgeR. The differences among identified subtypes were compared using the Kruskal-Wallis non-parametric test. To better clarify the subtypes of bladder cancer, their relationship with clinical characteristics was examined using the Fisher's test. We also constructed a risk prediction model using the random survival forest method to analyze right-censored survival data based on key metabolic genes. To identify genes of prognostic significance, univariate Cox regression, lasso analysis, and multivariate regression were performed sequentially.

Results

Three bladder cancer subtypes were identified according to the expression of metabolism-related genes. The M1 subtype was characterized by high metabolic activity, low immunogenicity, and better prognosis. M2 exhibited moderate metabolic activity, high immunogenicity, and the worst prognosis. M3 was associated with low metabolic activity, low immunogenicity, and poor prognosis. M1 showed the best predicted response to immunotherapy, whereas patients with M1 were predicted to be the least sensitive to cisplatin. By contrast, M2 showed the worst predicted response to immunotherapy but was predicted to be more sensitive to cisplatin, doxorubicin, and other first-line anticancer drugs. M3 was the most sensitive to gemcitabine. The risk model based on metabolic genes effectively predicted the prognosis of bladder cancer patients.

Conclusions

Metabolic classification of bladder cancer has potential clinical value and therapeutic feasibility by inhibiting the associated pathways. This classification can provide valuable insights for developing precise bladder cancer treatment.

SUBMITTER: Tang C

PROVIDER: S-EPMC8117567 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundRecently, immunotherapies have been approved for advanced muscle invasive bladder cancer (MIBC) treatment, but only a small fraction of MIBC patients could achieve a durable drug response. Our study is aimed at identifying tumor microenvironment (TME) subtypes that have different immunotherapy response rates.MethodsThe mRNA expression profiles of MIBC samples from seven discovery datasets (GSE13507, GSE31684, GSE32548, GSE32894, GSE48075, GSE48276, and GSE69795) were analyzed to identify TME subtypes. The identified TME subtypes were then validated by an independent dataset (TCGA-MIBC). The subtype-related biomarkers were discovered using computational analyses and then utilized to establish a random forest predictive model. The associations of TME subtypes with immunotherapy therapeutic responses were investigated in a group of patients who had been treated with immunotherapy. A prognostic index model was constructed using the subtype-related biomarkers. Two nomograms were built by the subtype-related biomarkers or the clinical parameters.ResultsTwo TME subtypes, including ECM-enriched class (EC) and immune-enriched class (IC), were found. EC was associated with greater extracellular matrix (ECM) pathways, and IC was correlated with immune pathways, respectively. Overall survival was significantly greater for tumors classified as IC, whereas the EC subtype had a worse prognosis. A total of nine genes (AKAP12, APOL3, CXCL13, CXCL9, GBP4, LRIG1, PEG3, PODN, and PTPRD) were selected by computational analyses to construct the random forest model. The area under the curve (AUC) values for this model were 0.827 and 0.767 in the testing and external validation datasets, respectively. Therapeutic response rates were greater in IC patients than in EC patients (28 percent vs. 18 percent). Patients with a high prognostic index had a poorer prognosis than those with a low prognostic index. The nomogram constructed from nine genes and stage achieved a C-index of 0.71.ConclusionThe present investigation defined two distinct TME subtypes and developed models to assess immunotherapeutic treatment outcomes.