Dataset Information

Development and Validation of a Hypoxia-Related Signature for Predicting Survival Outcomes in Patients With Bladder Cancer.

ABSTRACT:

Objectives

This study aimed to develop and validate a hypoxia signature for predicting survival outcomes in patients with bladder cancer.

Methods

We downloaded the RNA sequence and the clinicopathologic data of the patients with bladder cancer from The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov/repository?facetTab=files) and the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) databases. Hypoxia genes were retrieved from the Molecular Signatures Database (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp). Differentially expressed hypoxia-related genes were screened by univariate Cox regression analysis and Lasso regression analysis. Then, the selected genes constituted the hypoxia signature and were included in multivariate Cox regression to generate the risk scores. After that, we evaluate the predictive performance of this signature by multiple receiver operating characteristic (ROC) curves. The CIBERSORT tool was applied to investigate the relationship between the hypoxia signature and the immune cell infiltration, and the maftool was used to summarize and analyze the mutational data. Gene-set enrichment analysis (GSEA) was used to investigate the related signaling pathways of differentially expressed genes in both risk groups. Furthermore, we developed a model and presented it with a nomogram to predict survival outcomes in patients with bladder cancer.

Results

Eight genes (AKAP12, ALDOB, CASP6, DTNA, HS3ST1, JUN, KDELR3, and STC1) were included in the hypoxia signature. The patients with higher risk scores showed worse overall survival time than the ones with lower risk scores in the training set (TCGA) and two external validation sets (GSE13507 and GSE32548). Immune infiltration analysis showed that two types of immune cells (M0 and M1 macrophages) had a significant infiltration in the high-risk group. Tumor mutation burden (TMB) analysis showed that the risk scores between the wild types and the mutation types of TP53, MUC16, RB1, and FGFR3 were significantly different. Gene-Set Enrichment Analysis (GSEA) showed that immune or cancer-associated pathways belonged to the high-risk groups and metabolism-related signal pathways were enriched into the low-risk group. Finally, we constructed a predictive model with risk score, age, and stage and validated its performance in GEO datasets.

Conclusion

We successfully constructed and validated a novel hypoxia signature in bladder cancer, which could accurately predict patients' prognosis.

SUBMITTER: Zhang F

PROVIDER: S-EPMC8188560 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Recently, growing evidence has revealed the significant effect of reprogrammed metabolism on pancreatic cancer in relation to carcinogenesis, progression, and treatment. However, the prognostic value of metabolism-related genes in pancreatic cancer has not been fully revealed. We identified 379 differentially expressed metabolic-related genes (DEMRGs) by comparing 178 pancreatic cancer tissues with 171 normal pancreatic tissues in The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression project (GTEx) databases. Then, we used univariate Cox regression analysis together with Lasso regression for constructing a prognostic model consisting of 15 metabolic genes. The unified risk score formula and cutoff value were taken into account to divide patients into two groups: high risk and low risk, with the former exhibiting a worse prognosis compared with the latter. The external validation results of the International Cancer Genome Consortium (IGCC) cohort and the Gene Expression Omnibus (GEO) cohort further confirm the effectiveness of this prognostic model. As shown in the receiver operating characteristic (ROC) curve, the area under curve (AUC) values of the risk score for overall survival (OS), disease-specific survival (DSS), and progression-free survival (PFS) were 0.871, 0.885, and 0.886, respectively. Based on the Gene Set Enrichment Analysis (GSEA), the 15-gene signature can affect some important biological processes and pathways of pancreatic cancer. In addition, the prognostic model was significantly correlated with the tumor immune microenvironment (immune cell infiltration, and immune checkpoint expression, etc.) and clinicopathological features (pathological stage, lymph node, and metastasis, etc.). We also built a nomogram based on three independent prognostic predictors (including individual neoplasm status, lymph node metastasis, and risk score) for the prediction of 1-, 3-, and 5-year OS of pancreatic cancer, which may help to further improve the treatment strategy of pancreatic cancer.

Project description:Background: Pancreatic ductal adenocarcinoma (PDAC) is a highly fatal, aggressive cancer characterized by invasiveness and metastasis. In this study, we aimed to propose a gene prediction model based on metastasis-related genes (MTGs) to more accurately predict PDAC prognosis. Methods: Differentially expressed MTGs (DE-MTGs) were identified via integrated analysis of gene expression omnibus (GEO) datasets and Human Cancer Metastasis Database (HCMDB). Overall survival (OS) related DE-MTGs were then identified and a prognostic gene signature was established using Lasso-Cox regression with TCGA-PAAD datasets. Tumor immunity was analyzed using ESTIMATE and CIBERSORT algorithms. Finally, a nomogram predicting 1-year, 2-year, and 3-year OS of PDAC patients was established based on the prognostic gene signature and relevant clinical parameters using a stepwise Cox regression model. Results: A total of 36 DE-MTGs related to OS were identified in PDAC. Consequently, an MTG-based gene signature comprising of RACGAP1, RARRES3, TPX2, MMP28, GPR87, KIF14, and TSPAN7 was established to predict the OS of PDAC. The MTG-based gene signature was able to distinguish high-risk patients with significantly poorer prognosis and accurately predict OS of PDAC in both the training and external validation datasets. Cox regression analysis indicated that the MTG-based gene signature was an independent prognostic factor in PDAC. The gene set enrichment analysis (GSEA) showed that molecular alterations in the high-risk group were associated with multiple oncological pathways. Moreover, analysis of tumor immunity revealed significantly higher levels of follicular helper T cells and M0 macrophage infiltration, and lower levels of infiltrating naïve B cells, CD8 T cells, monocytes, and resting dendritic cells in the high-risk group. Immune cell infiltration levels were significantly associated with the expression of the seven DE-MTGs. Finally, a nomogram was established by incorporating the prognostic gene signature and clinical parameters, which was superior to the AJCC staging system in predicting the OS of PDAC patients. Conclusions: The DE-MTGs we identified were closely associated with the progress and prognosis of PDAC and are potential therapeutic targets. The MTG-based gene signature and nomogram may serve to improve the individualized prediction of survival, assisting in clinical decision-making.

Project description:BackgroundBoth lncRNAs and glycolysis are considered to be key influencing factors in the progression of bladder cancer (BCa). Studies have shown that glycolysis-related lncRNAs are an important factor affecting the overall survival and prognosis of patients with bladder cancer. In this study, a prognostic model of BCa patients was constructed based on glycolysis-related lncRNAs to provide a point of reference for clinical diagnosis and treatment decisions.MethodsThe transcriptome, clinical data, and glycolysis-related pathway gene sets of BCa patients were obtained from The Cancer Genome Atlas (TCGA) database and the Gene Set Enrichment Analysis (GSEA) official website. Next, differentially expressed glycolysis-related lncRNAs were screened out, glycolysis-related lncRNAs with prognostic significance were identified through LASSO regression analysis, and a risk scoring model was constructed through multivariate Cox regression analysis. Then, based on the median of the risk scores, all BCa patients were divided into either a high-risk or low-risk group. Kaplan-Meier (KM) survival analysis and the receiver operating characteristic (ROC) curve were used to evaluate the predictive power of the model. A nomogram prognostic model was then constructed based on clinical indicators and risk scores. A calibration chart, clinical decision curve, and ROC curve analysis were used to evaluate the predictive performance of the model, and the risk score of the prognostic model was verified using the TCGA data set. Finally, Gene Set Enrichment Analysis (GSEA) was performed on glycolysis-related lncRNAs.ResultsA total of 59 differentially expressed glycolysis-related lncRNAs were obtained from 411 bladder tumor tissues and 19 pericarcinomatous tissues, and 9 of those glycolysis-related lncRNAs (AC099850.3, AL589843.1, MAFG-DT, AC011503.2, NR2F1-AS1, AC078778.1, ZNF667-AS1, MNX1-AS1, and AC105942.1) were found to have prognostic significance. A signature was then constructed for predicting survival in BCa based on those 9 glycolysis-related lncRNAs. ROC curve analysis and a nomogram verified the accuracy of the signature.ConclusionThrough this study, a novel prognostic prediction model for BCa was established based on 9 glycolysis-related lncRNAs that could effectively distinguish high-risk and low-risk BCa patients, and also provide a new point of reference for clinicians to make individualized treatment and review plans for patients with different levels of risk.

Project description:BackgroundAutophagy, a highly conserved self-digesting process, has been deeply involved in the development and progression of oral squamous cell carcinoma (OSCC). However, the prognostic value of autophagy-related genes (ARGs) for OSCC still remains unclear. Our study set out to develop a multigene expression signature based on ARGs for individualized prognosis assessment in OSCC patients.MethodsBased on The Cancer Genome Atlas (TCGA) database, we identified prognosis-related ARGs through univariate COX regression analysis. Then we performed the least absolute shrinkage and selection operator (LASSO) regression analysis to identify an optimal autophagy-related multigene signature with the subsequent validation in testing set, GSE41613 and GSE42743 datasets.ResultsWe identified 36 prognosis-related ARGs for OSCC. Subsequently, the multigene signature based on 13 prognostic ARGs was constructed and successfully divided OSCC patients into low and high-risk groups with significantly different overall survival in TCGA training set (p < 0.0001). The autophagy signature remained as an independent prognostic factor for OSCC in univariate and multivariate Cox regression analyses. The area under the curve (AUC) values of the receiver operating characteristic (ROC) curves for 1, 3, and 5-year survival were 0.758, 0.810, 0.798, respectively. Then the gene signature was validated in TCGA testing set, GSE41613 and GSE42743 datasets. Moreover, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and single-sample gene set enrichment analysis (ssGSEA) revealed the underlying biological characteristics and signaling pathways associated with this signature in OSCC. Finally, we constructed a nomogram by combining the gene signature with multiple clinical parameters (age, gender, TNM-stage, tobacco, and alcohol history). The concordance index (C-index) and calibration plots demonstrated favorable predictive performance of our nomogram.ConclusionIn summary, we identified and verified a 13-ARGs prognostic signature and nomogram, which provide individualized prognosis evaluation and show insight for potential therapeutic targets for OSCC.