Identification of immune-related genes that predict prognosis and risk of bladder cancer: bioinformatics analysis of TCGA database.
Ontology highlight
ABSTRACT: Background
Bladder cancer (BLCA) is the major tumor of the urinary system, and immune-related genes (IRGs) contribute significantly to its initiation and prognosis.Results
A total of 51 prognostic IRGs significantly associated with overall survival were identified. Functional enrichment analysis revealed that these genes were actively involved in tumor-related functions and pathways. Using multivariate Cox regression analysis, we detected 11 optimal IRGs (ADIPOQ, PPY, NAMPT, TAP1, AHNAK, OLR1, PDGFRA, IL34, MMP9, RAC3, and SH3BP2). We validated the prognostic value of this signature in two validation cohorts: GSE13507 (n = 165) and GSE32894 (n = 224). Furthermore, we performed a western blot and found that the expression of these IRGs matched their mRNA expression in TCGA. Moreover, correlations between risk score and immune-cell infiltration indicated that the prognostic signature reflected infiltration by several types of immune cells.Conclusion
We identified and validated an 11-IRG-based risk signature that may be a reliable tool to evaluate the prognosis of BLCA patients and help to devise individualized immunotherapies.Methods
Bioinformatics analysis was performed using TCGA and ImmPort databases. Cox regression was used to identify prognostic signatures. Two external GEO cohorts and western blotting of samples were performed to validate the IRG signature.
SUBMITTER: Guo L
PROVIDER: S-EPMC8386543 | biostudies-literature |
REPOSITORIES: biostudies-literature
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.
| S-EPMC7658688 | biostudies-literature
Project description:BackgroundMuscle-invasive bladder cancer (MIBC) with high tumor stages accounts for most bladder cancer patient mortality. Platinum-based chemotherapy provides insufficient survival benefits; however, immunotherapy is a promising option for MIBC.ResultsThere were 31 differentially expressed IRGs that significantly correlated with the clinical outcomes of MIBC patients. A prognostic signature based on 12 IRGs (MMP9, RBP7, ADIPOQ, AHNAK, OAS1, RAC3, SLIT2, EDNRA, IL34, PDGFD, PPY, IL17RD) performed moderately in prognostic predictions with area under the curve (AUC) equal to 0.76. The high-risk patient group presented worse survival outcomes (hazard ratio 1.197, 95% confidence interval 1.103-1.299, p < 0.001). Furthermore, immune cell infiltration analysis showed increased tumor infiltration of macrophages in the high-risk group.ConclusionThis novel prognostic signature can effectively divide MIBC patients into different risk groups, allowing for intensive treatment of high-risk individuals who have worse predicted survival outcomes.MethodsBioinformatics analyses were conducted using the Cancer Genome Atlas (TCGA) database. Differentially expressed genes and survival-associated immune-related genes (IRGs) were analyzed through a computational algorithm and Cox regression. The potential mechanisms of IRG expression were explored with transcription factors, and a prognosis classification based on IRG expression was developed to stratify patients into distinct risk groups.
| S-EPMC7880322 | biostudies-literature
Project description:Bladder cancer is the most common cancer of the urinary system and its treatment has scarcely progressed for nearly 30 years. Advances in checkpoint inhibitor research have seemingly provided a new approach for treatment. However, there have been issues predicting immunotherapeutic biomarkers and identifying new therapeutic targets. We downloaded the gene expression profile and clinical data of 408 cases bladder urinary cancer from the Cancer Genome Atlas (TCGA) portal, and the abundance ratio of immune cells for each sample was obtained via the "Cell Type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT)" algorithm. Then, four survival-related immune cells were obtained via Kaplan-Meier survival analysis, and 933 immune-related genes were obtained via a variance analysis. Enrichment, protein-protein interaction, and co-expression analyses were performed for these genes. Lastly, 4 survival-related immune cells and 24 hub genes were identified, four of which were related to overall survival. More importantly, these immune cells and genes were closely related to the clinical features. These cells and genes may have research value and clinical application in bladder cancer immunotherapy. Our study not only provides cell and gene targets for bladder cancer immunotherapy, but also provides new ideas for researchers to explore the immunotherapy of various tumors.
| S-EPMC6974676 | biostudies-literature
Project description:BackgroundRecurrence 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.MethodsThe 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.ResultsResults 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).ConclusionsWe 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.
| S-EPMC8749062 | biostudies-literature
Project description:PurposeIdentify immune-related lncRNA (IRL) signature related to the prognosis and immunotherapeutic efficiency for bladder cancer (BLCA) patients.MethodsA total of 397 samples, which contained RNA-seq and clinical information from The Cancer Genome Atlas (TCGA) database, were used for the following study. Then the Lasso penalized Cox proportional hazards regression model was used to construct prognostic signature. According to the optimal cut-off value determined by time-dependent ROC curve, low and high-risk groups were set up. One immunotherapy microarray dataset as validation set was used to verify the ability of predicting immunotherapy efficacy. Furthermore, more evaluation between two risk groups related clinical factors were conducted. Finally, external validation of IRL-signature was conducted in Zhengzhou cohort.ResultFour IRLs (HCP5, IPO5P1, LINC00942, and LINC01356) with significant prognostic value (P<0.05) were distinguished. This signature can accurately predict the overall survival of BLCA patients and was verified in the immunotherapy validation set. IRL-signatures can be used as independent prognostic risk factor in various clinical subgroups. According to the results of GSVA and MCP algorithm, we found that IRL-signature risk score is strikingly negative correlated with tumor microenvironment (TME) CD8+T cells and Cytotoxic lymphocytes infiltration, indicating that the better prognosis and immunotherapy might be caused partly by these. Then, the results from the TIDE analysis revealed that IRL could efficiently predict the response of immunotherapy in BLCA. External validation had similar results with TCGA-BLCA cohort.ConclusionsThe novel IRL-signature has a significant prognostic value for BLCA patients might facilitate predicting the efficacy of immunotherapy.
| S-EPMC7855860 | biostudies-literature
Project description:Background Bladder cancer is the tenth most common cancer worldwide. Valuable biomarkers in the field of diagnostic bladder cancer are urgently required. Method Here, the gene expression matrix and clinical data were obtained from The Cancer Genome Atlas (TCGA), GSE13507, GSE32894, and Mariathasan et al. Five prognostic genes were identified by the univariate, robust, and multivariate Cox's regression and were used to develop a prognosis-related model. The Kaplan–Meier survival curves and receiver operating characteristics were used to evaluate the model's effectiveness. The potential biological functions of the selected genes were analyzed using CIBERSORT and ESTIMATE algorithms. Cancer Therapeutics Response Portal (CTRP) and PRISM datasets were used to identify drugs with high sensitivity. Subsequently, using the bladder cancer (BLCA) cell lines, the role of TNFRSF14 was determined by Western blotting, cell proliferation assay, and 5-ethynyl-20-deoxyuridine assay. Results GSDMB, CLEC2D, APOL2, TNFRSF14, and GBP2 were selected as prognostic genes in bladder cancer patients. The model's irreplaceable reliability was validated by the training and validation cohorts. CD8+ T cells were highly infiltrated in the high-TNFRSF14-expression group, and M2 macrophages were the opposite. Higher expression of TNFRSF14 was associated with higher expression levels of LCK, interferon, MHC-I, and MHC-II, while risk score was the opposite. Many compounds with higher sensitivity for treating bladder cancer patients in the low-TNFRSF14-expression group were identified, with obatoclax being a potential drug most likely to treat patients in the low-TNFRSF14-expression group. Finally, the proliferation of BLCA cell lines was increased in the TNFRSF14-reduced group, and the differential expression was identified. TNFRSF14 plays a role in bladder cancer progression through the Wnt/β-catenin-dependent pathway. TNFRSF14 is a potential protective biomarker involved in cell proliferation in BLCA. Conclusion We conducted a study to establish a 5-gene score model, providing reliable prediction for the outcome of bladder cancer patients and therapeutic drugs to individualize therapy. Our findings provide a signature that might help determine the optimal treatment for individual patients with bladder cancer.
| S-EPMC8545590 | biostudies-literature
Project description:IntroductionBladder cancer (BC) is one of the most prevalent urinary cancers, and its management is still a problem causing recurrence and progression, elevating mortality.Materials and methodsWe aimed at the nuclear mitochondria-related genes (MTRGs), collected from the MITOMAP: A Human Mitochondrial Genome Database. Meanwhile, the expression profiles and clinical information of BC were downloaded from the Cancer Genome Atlas (TCGA) as a training group. The univariate, multivariate, and the least absolute shrinkage and selection operator (LASSO) Cox regression analyses were used to construct a nuclear mitochondrial-related multi-genes signature and the prognostic nomogram.ResultsA total of 17 nuclear MTRGs were identified to be correlated with the overall survival (OS) of BC patients, and a nuclear MTRGs signature based on 16 genes expression was further determined by the LASSO Cox regression analysis. Based on a nuclear MTRGs scoring system, BC patients from the TCGA cohort were divided into high- and low- nuclear MTRGs score groups. Patients with a high nuclear MTRGs score exhibited a significantly poorer outcome (median OS: 92.90 vs 20.20 months, p<0.0001). The nuclear MTRGs signature was further verified in three independent datasets, namely, GSE13507, GSE31684, and GSE32548, from the Gene Expression Omnibus (GEO). The BC patients with a high nuclear MTRGs score had significantly worse survival (median OS in GSE13507: 31.52 vs 98.00 months, p<0.05; GSE31684: 32.85 months vs unreached, p<0.05; GSE32548: unreached vs unreached, p<0.05). Furthermore, muscle-invasive bladder cancer (MIBC) patients had a significantly higher nuclear MTRGs score (p<0.05) than non-muscle-invasive bladder cancer (NMIBC) patients. The integrated signature outperformed each involved MTRG. In addition, a nuclear MTRGs-based nomogram was constructed as a novel prediction prognosis model, whose AUC values for OS at 1, 3, 5 years were 0.76, 0.75, and 0.75, respectively, showing the prognostic nomogram had good and stable predicting ability. Enrichment analyses of the hallmark gene set and KEGG pathway revealed that the E2F targets, G2M checkpoint pathways, and cell cycle had influences on the survival of BC patients. Furthermore, the analysis of tumor microenvironment indicated more CD8+ T cells and higher immune score in patients with high nuclear MTRGs score, which might confer sensitivity to immune checkpoint inhibitors.ConclusionsNot only could the signature and prognostic nomogram predict the prognosis of BC, but it also had potential therapeutic guidance.
| S-EPMC8528313 | biostudies-literature
Project description:AimsIn the cancer-related research field, there is currently a major need for a greater number of valuable biomarkers to predict the prognosis of hepatocellular carcinoma (HCC). In this study, we aimed to screen hub genes related to immune cell infiltration and explore their prognostic value for HCC.MethodsWe analyzed five datasets (GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630) from the Gene Expression Omnibus database to screen the differentially expressed genes (DEGs). A protein-protein interaction network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes; then, the hub genes were identified. Functional enrichment of the genes was performed on the Metascape website. Next, the expression of these hub genes was validated in several databases, including Oncomine, Gene Expression Profiling Interactive Analysis 2 (GEPIA2), and Human Protein Atlas. We explored the correlations between the hub genes and infiltrated immune cells in the TIMER2.0 database. The survival curves were generated in GEPIA2, and the univariate and multivariate Cox regression analyses were performed using TIMER2.0.ResultsThe top ten hub genes [DNA topoisomerase II alpha (TOP2A), cyclin B2 (CCNB2), protein regulator of cytokinesis 1 (PRC1), Rac GTPase-activating protein 1 (RACGAP1), aurora kinase A (AURKA), cyclin-dependent kinase inhibitor 3 (CDKN3), nucleolar and spindle-associated protein 1 (NUSAP1), cell division cycle-associated 5 (CDCA5), abnormal spindle microtubule assembly (ASPM), and non-SMC condensin I complex subunit G (NCAPG)] were identified in subsequent analysis. These genes are most markedly enriched in cell division, suggesting their close association with tumorigenesis. Multi-database analyses validated that the hub genes were upregulated in HCC tissues. All hub genes positively correlated with several types of immune infiltration, including B cells, CD4+ T cells, macrophages, and dendritic cells. Furthermore, these hub genes served as independent prognostic factors, and the expression of these hub genes combing with the macrophage levels could help predict an unfavorable prognosis of HCC.ConclusionIn sum, these hub genes (TOP2A, CCNB2, PRC1, RACGAP1, AURKA, CDKN3, NUSAP1, CDCA5, ASPM, and NCAPG) may be pivotal markers for prognostic prediction as well as potentially work as targets for immune-based intervention strategies in HCC.
| S-EPMC7831279 | biostudies-literature
Project description:Bladder cancer is one of the most common malignant tumors of the urinary system that seriously threatens the health of a population. In recent years, the application of immunotherapy has significantly changed the treatment of bladder cancer, but only some patients can benefit from the treatment with immune-checkpoint inhibitors. Many problems are unsolved in the field of bladder cancer immunotherapy, especially in the search for genes that are critical to the level of immune cell infiltration and new effective therapeutic targets. We attempted to use bioinformatics analysis to identify immune gene markers related to the prognosis of bladder cancer and to establish a prognostic signature for bladder cancer patients based on their immune gene expression profiles. We used univariate Cox proportional hazards regression analysis, the least absolute shrinkage and selection operator (LASSO) Cox regression, and multivariate Cox proportional hazards regression analysis from The Cancer Genome Atlas bladder cancer cohort (TCGA-BLCA). Fifteen genes related to prognosis were screened using the survival analysis, correlation analysis, cancer and adjacent cancer differential expression analysis, and mutation analysis. The potential biological role of these genes was determined using survival analysis and principal component analysis (PCA). The receiver operating characteristic (ROC) curve assesses the prognostic value of the predictive signature. The gene ontology (GO), Kyoto Encyclopedia of Gene and Genome (KEGG), Gene set enrichment analysis (GSEA), and other methods were used to reveal the differential gene enrichment in the signaling pathways and cellular processes of high- and low-risk groups. The single-sample GSEA (ssGSEA) method was used to quantify the infiltration levels of 24 immune cells in the tumor immune microenvironment and these immune genes were found to be closely related to the tumor immune microenvironment. In summary, we screened 15 immune genes that were closely related to bladder cancer overall survival (OS) and may be potential prognostic indicators of bladder cancer. They may have research and clinical application value in bladder cancer immunotherapy. We used 15 immune genes to construct a new immune-related gene signature that was verified and could be helpful in improving individualized prognosis of patients with bladder cancer.
| S-EPMC7324668 | biostudies-literature
Project description:Breast cancer is one of the most common malignancies in women worldwide. Many studies have shown that tumor microenvironment cells, immune cells, and stromal cell infiltration have an important impact on prognosis, so it is important to identify biomarkers for achieving better treatment and prognosis.To better understand the relationship between immune and stromal cell-related genes and prognosis, we screened patients with breast cancer in The Cancer Genome Atlas (TCGA) database and divided them into high and low groups based on immune/stromal scores. We next identified differentially expressed immune-related genes that are significantly associated with the prognosis of patients with breast cancer for functional enrichment analysis and protein-protein interaction networks, respectively. Finally, we selected a separate breast cancer cohort in gene expression synthesis (GEO) for validation.Both immune scores and stromal scores are meaningful in the correlation of subtype classification. Disease-free survival of cases with the high score group of immune scores is statistically longer than the cases in the low score group. Differentially expressed immune-related genes extracted from the comparison can effectively evaluate the prognosis of patients with breast cancer and these genes are primarily involved in immune responses, extracellular matrix, and chemokine activity. At last, we obtained a series of verified tumor immune-related genes that predict the prognosis of patients with breast cancer.Combining the Estimation of Stromal and Immune Cells in Malignant Tumor Tissues using Expression database and the TCGA database to extract the list of tumor microenvironment related genes which may help to outline the prognosis of patients with breast cancer. Some previously overlooked genes have the potential to become additional biomarkers for breast cancer. Further research on these genes can reveal a new understanding of the potential relationship between tumor microenvironment and breast cancer prognosis.
| S-EPMC7220520 | biostudies-literature