Dataset Information

Identifying a hypoxia related score to predict the prognosis of bladder cancer: a study with The Cancer Genome Atlas (TCGA) database.

ABSTRACT:

Background

Recurrence is common in bladder cancer, with a hypoxic tumor microenvironment (TME) playing a role in genetic instability and prognosis of bladder cancer. However, we still lack practical hypoxia related model for predicting the prognosis of bladder cancer. In this study, we identified new prognosis-related hypoxia genes and established a new hypoxia score related signature.

Methods

The Gene Set Variation Analysis (GSVA) algorithm was utilized to calculate the hypoxia score of bladder cancer cases found on the The Cancer Genome Atlas (TCGA) database on the gene expression profiles. The cases were first divided into low- and high-hypoxia score groups and then differentially expressed genes (DEGs) expression analysis was conducted. Hypoxia-related genes were identified using weighted gene co-expression network analysis (WGCNA). We then conducted a protein-protein interaction (PPI) network and carried out functional enrichment analysis of the genes that overlapped between DEGs and hypoxia-related genes. LASSO Cox regression analysis was used to establish a hypoxia-related prognostic signature, which was validated using the GSE69795 dataset downloaded from GEO database.

Results

Results from Kaplan-Meier analysis showed that patients with a high hypoxia score had significantly poor overall survival compared to patients with low hypoxia score. We selected 270 DEGs between low- and high-hypoxia score groups, while WGCNA analysis identified 1,313 genes as hypoxia-related genes. A total of 170 genes overlapped between DEGs and hypoxia-related genes. LASSO algorithms identified 29 genes associated with bladder cancer prognosis, which were used to construct a novel 29-gene signature model. The prognostic risk model performed well, since the receiver operating characteristic (ROC) curve showed an accuracy of 0.802 (95% CI: 0.759-0.844), and Cox proportional hazards regression analysis proved the model an independent predictor with hazard ratio (HR) =1.789 (95% CI: 1.585-2.019) (P<0.001). The low-risk score patients had remarkably longer overall survival than patients with a higher score (survival rate 71.06% vs. 23.66%) in the The Cancer Genome Atlas (TCGA) cohort (P<0.0001) and in the dataset GSE69795 (P=0.0079).

Conclusions

We established a novel 29-gene hypoxia-related signature model to predict the prognosis of bladder cancer cases. This model and identified hypoxia-related genes may further been used as biomarkers, assisting the evaluation of prognosis of bladder cancer cases and decision making in clinical practice.

SUBMITTER: Zhang Z

PROVIDER: S-EPMC8749062 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundBladder cancer (BCa) is a common urothelial malignancy. The Cancer Genome Atlas (TCGA) database allows for an opportunity to analyze the relationship between gene expression and clinical outcomes in bladder cancer patients. This study is aimed at identifying prognosis-related genes in the bladder cancer microenvironment.MethodsImmune scores and stromal scores were calculated by applying the ESTIMATE algorithm. We divided bladder cancer patients into high and low groups based on their immune/stromal scores. Then, differentially expressed genes (DEGs) were identified in bladder cancer patients based on the TCGA database. We evaluated the correlation between immune/stromal scores and clinical characteristics as well as prognosis. Finally, we validated identified genes associated with bladder cancer prognosis through a cohort study in the Gene Expression Omnibus (GEO) database.ResultsA higher stromal score was associated with female (vs. malep = 0.037), age > 65 (vs.age ≤ 65 p = 0.015), T3/4 (vs. T1/2,p < 0.001), N status(p = 0.016), and pathological high grade (vs. low gradeP < 0.001). By analyzing DEGs, there were 1125 genes commonly upregulated, and 209 genes were commonly downregulated. Protein-protein interaction networks further showed the important protein that may be involved in the biological behavior and prognosis of BCa, such as FN1, CXCL12, CD3E, LCK, and ZAP70. A total of 14 DEGs were found to be associated with overall survival of bladder cancer. After validation by a cohort of 165 BCa cases with detailed follow-up information from GSE13507, 10 immune-associated DEGs were demonstrated to be predictive of prognosis in BCa. Among them, 5 genes have not been reported previously associated with the prognosis of BCa, including BTBD16, OLFML2B, PRRX1, SPINK4, and SPON2.ConclusionsOur study elucidated tight associations between stromal score and clinical characteristics as well as prognosis in BCa. Moreover, we obtained a group of genes closely related to the prognosis of BCa in the tumor microenvironment.

Project description:BackgroundThe incidence rate of prostate cancer is increasing rapidly. This study aims to explore the gene-associated mechanism of prostate cancer biochemical recurrence (BCR) after radical prostatectomy and to construct a biochemical recurrence of prostate cancer prognostic model.MethodsThe DEseq2 R package was used for the differential expression of mRNA. The ClusterProfiler R package was used to analyze the functional enrichment of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to explore related mechanisms. The Survival, Survminer, and My.stepwise R packages were used to construct the prognostic model to predict the biochemical recurrence-free probability. The RMS R package was used to draw the nomogram. For evaluating the prognostic model, the timeROC R package was used to draw the time-dependent ROC curve (receiver operating characteristic curve).ResultTo investigate the association between mRNA and prostate cancer, we performed differential expression analysis on the TCGA (The Cancer Genome Atlas) database. Seven protein-coding genes (VWA5B2, ARC, SOX11, MGAM, FOXN4, PRAME, and MMP26) were picked as independent prognostic genes by regression analysis. Based on their Cox coefficient, a risk score formula was proposed. According to the risk scores, patients were divided into high- and low-risk groups based on the median score. Kaplan-Meier plot curves showed that the low-risk group had a better biochemical recurrence-free probability compared to the high-risk group. The 1-year, 3-year, and 5-year AUCs (areas under the ROC curve) of the model were 77%, 81%, and 86%, respectively. In addition, we built a nomogram based on the result of multivariate Cox regression analysis. Furthermore, we select the GSE46602 dataset as our external validation. The 1-year, 3-year, and 5-year AUCs of BCR-free probability were 83%, 82%, and 80%, respectively. Finally, the levels of seven genes showed a difference between PRAD tissues and adjacent non-tumorous tissues.ConclusionsThis study shows that establishing a biochemical recurrence prediction prognostic model comprising seven protein-coding genes is an effective and precise method for predicting the progression of prostate cancer.

Project description:Background and aimHepatocellular carcinoma is a common malignant tumor of the digestive system with a poor prognosis. The high recurrence rate and metastasis after surgery reduce the survival time of patients. Therefore, assessing the overall survival of patients with hepatocellular carcinoma after hepatectomy is critical to clinicians' clinical decision-making. Conventional hepatocellular carcinoma assessment systems (such as tumor lymph node metastasis and Barcelona clinical hepatocellular carcinoma) are obviously insufficient in assessing the overall survival rate of patients. This research is devoted to the development of nomogram assessment tools to assess the overall survival probability of patients undergoing liver resection.MethodsWe collected the clinical and pathological information of 438 hepatocellular carcinoma patients undergoing surgery from The Cancer Genome Atlas (TCGA) database, then excluded 87 patients who did not meet inclusion criteria. Univariate and multivariate analyses were performed on patient characteristics and related pathological factors. Finally, we developed a nomogram model to predict patient's prognosis.ResultsA retrospective analysis of 438 consecutive samples from the TCGA database of patients with hepatocellular carcinoma who underwent potentially curative liver resections. Six risk factors were included in the final model. In the training set, the discriminative ability of the nomogram was very good (concordance index = 0.944), and the external verification method (concordance index = 0.962) was used for verification. At the same time, the internal and external calibration of the model was verified, showing that the model was well calibrated. The calibration between the evaluation of the nomogram and the actual observations was good. According to the patient's risk factors, we determined the patient's Kaplan-Meyer survival analysis curve. Finally, the clinical decision curve was used to compare the benefits of two different models in evaluating patients' clinical outcomes.ConclusionsThe nomogram can be used to evaluate the post-hepatectomy 1-, 3-, and 5-year survival rates of patients with hepatocellular carcinoma. The Kaplan-Meyer curve can intuitively display the survival differences among patients with various risk factors. The clinical decision curve is a good reference guide for clinical application.

Project description:BackgroundBladder cancer (BLCA) is the ninth most common cancer globally, as well as the fourth most common cancer in men, with an incidence of 7%. However, few effective prognostic biomarkers or models of BLCA are available at present.MethodsThe prognostic genes of BLCA were screened from one cohort of The Cancer Genome Atlas (TCGA) database through univariate Cox regression analysis and functionally annotated by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The intersecting genes of the BLCA gene set and focal adhesion-related gene were obtained and subjected to the least absolute shrinkage and selection operator regression (LASSO) to construct a prognostic model. Gene set enrichment analysis (GSEA) of high- and low-risk patients was performed to explore further the biological process related to focal adhesion genes. Univariate and multivariate Cox analysis, receiver operating characteristic (ROC) curve analysis, and Kaplan-Meier survival analysis (KM) were used to evaluate the prognostic model. DNA methylation analysis was presented to explore the relationship between prognosis and gene methylation. Furthermore, immune cell infiltration was assessed by CIBERSORT, ESTIMATE, and TIMER. The model was verified in an external GSE32894 cohort of the Gene Expression Omnibus (GEO) database, and the Prognoscan database presented further validation of genes. The HPA database validated the related protein level, and functional experiments verified significant risk factors in the model.ResultsVCL, COL6A1, RAC3, PDGFD, JUN, LAMA2, and ITGB6 were used to construct a prognostic model in the TCGA-BLCA cohort and validated in the GSE32894 cohort. The 7-gene model successfully stratified the patients into both cohorts' high- and low-risk groups. The higher risk score was associated with a worse prognosis.ConclusionsThe 7-gene prognostic model can classify BLCA patients into high- and low-risk groups based on the risk score and predict the overall survival, which may aid clinical decision-making.