Dataset Information

Tumor Microenvironment Analysis Identified Subtypes Associated With the Prognosis and the Tumor Response to Immunotherapy in Bladder Cancer

ABSTRACT: Background: The efficiency of immune checkpoint inhibitors (ICIs) in bladder cancer (BLCA) treatment has been widely validated; however, the tumor response to ICIs was generally low. It is critical and urgent to find biomarkers that can predict tumor response to ICIs. The tumor microenvironment (TME), which may play important roles to either dampen or enhance immune responses, has been widely concerned. Methods: The cancer genome atlas BLCA (TCGA-BLCA) cohort (n = 400) was used in this study. Based on the proportions of 22 types of immune cells calculated by CIBERSORT, TME was classified by K-means Clustering and differentially expressed genes (DEGs) were determined. Based on DEGs, patients were classified into three groups, and cluster signature genes were identified after reducing redundant genes. Then TMEscore was calculated based on cluster signature genes, and the samples were classified to two subtypes. We performed somatic mutation and copy number variation analysis to identify the genetic characteristics of the two subtypes. Correlation analysis was performed to explore the correlation between TMEscore and the tumor response to ICIs as well as the prognosis of BLCA. Results: According to the proportions of immune cells, two TME clusters were determined, and 1,144 DEGs and 138 cluster signature genes were identified. Based on cluster signature genes, samples were classified into TMEscore-high (n = 199) and TMEscore-low (n = 201) subtypes. Survival analysis showed patients with TMEscore-high phenotype had better prognosis. Among the 45 differentially expressed micro-RNAs (miRNAs) and 1,033 differentially expressed messenger RNAs (mRNAs) between the two subtypes, 16 miRNAs and 287 mRNAs had statistically significant impact on the prognosis of BLCA. Furthermore, there were 94 genes with significant differences between the two subtypes, and they were enriched in RTK-RAS, NOTCH, WNT, Hippo, and PI3K pathways. The Tumor Immune Dysfunction and Exclusion (TIDE) score of TMEscore-high BLCA was statistically lower than that of TMEscore-low BLCA. Receiver operating characteristic (ROC) curve analysis showed that the area under the curve (AUC) of TMEscore and tumor mutation burden (TMB) is 0.6918 and 0.5374, respectively. Conclusion: We developed a method to classify BLCA patients to two TME subtypes, TMEscore-high and TMEscore-low, and we found TMEscore-high subtype of BLCA had a good prognosis and a good response to ICIs.

SUBMITTER: Zhang H

PROVIDER: S-EPMC7957069 | biostudies-literature | 2021 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundBladder cancer (BLCA) is prone to metastasis and has poor prognosis with unsatisfactory treatment responsiveness. Disulfidptosis is a recently discovered, novel mode of cell death that is closely associated with human cancers. However, a comprehensive analysis of the relationship between disulfidptosis and BLCA is lacking. Therefore, this study aimed to explore the potential effect of disulfidptosis on BLCA and identify a biomarker for evaluating the prognosis and immunotherapy of patients with BLCA.Material and methodsWe acquired BLCA RNA sequencing data from The Cancer Genome Atlas Urothelial Bladder Carcinoma (TCGA-BLCA) cohort (containing 19 normal samples and 409 tumor samples) and the GES39281 cohort (containing 94 tumor samples) which were used for external validation of the signature. Initially, we performed unsupervised consensus clustering to explore disulfidptosis-related subgroups. We then conducted functional enrichment analysis on these subgroups to gain insights into their biological significance and evaluate their immunotherapy response and chemotherapy sensitivity. Next, we conducted Least Absolute Shrinkage and Selection Operator (LASSO) regression and multivariate Cox regression to construct a prognostic signature in the TCGA training set for prognosis-related differentially expressed genes (DEGs) in the disulfidptosis-related subgroups. Subsequently, we used a receiver operating characteristic (ROC) curve and independent prognostic analysis to validate the predictive performance of the signature in the TCGA testing and the GES39281 cohorts. Finally, we explored the therapeutic value of this signature in patients with BLCA, in terms of immunotherapy and chemotherapy.ResultIn this study, we obtained two subgroups: DRG-high (238 samples) and DRG-low (160 samples). The DRG-high group exhibited a poor survival rate compared to the DRG-low group and had a significant association with tumor grade, stage, and metastasis. Additionally, several pathways related to cancer and the immune system were enriched in the high-DRG group. Moreover, the DRG-high group exhibited higher expression of PD1 and CTLA4 and had a better response to immunotherapy in patients with both PD1 and CTLA4 positivity. Conversely, the DRG-high group was more sensitive to common chemotherapeutic agents. A prognostic signature was created, consisting of COL5A1, DIRAS3, NKG7, and POLR3G and validated as having a robust predictive capability. Patients in the low-risk-score group had more immune cells associated with tumor suppression and better immunotherapy outcomes.ConclusionThis study contributes to our understanding of the characteristics of disulfidptosis-related subgroups in BLCA. Disulfidptosis-related signatures can be used to assess the prognosis and immunotherapy of patients with BLCA.

Project description:BackgroundBiomarkers based on immune context may guide prognosis prediction. T-cell inactivation, exclusion, or dysfunction could cause unfavorable tumor microenvironments, which affect immunotherapy and prognosis. However, none of the immuno-biomarkers reported to date can differentiate colorectal-cancer (CRC) patients. Thus, we aimed to classify CRC patients according to the levels of T-cell activation, exclusion, and dysfunction in the tumor microenvironment.MethodsRNAseq data of 618 CRC patients from The Cancer Genome Atlas and microarray data of 316 CRC patients from Gene Expression Omnibus were analysed using the Tumor Immune Dysfunction and Exclusion algorithm. Unsupervised clustering was used to classify patients.ResultsBased on the expression signatures of myeloid-derived suppressor cells, cancer-associated fibroblasts, M2-like tumor-associated macrophages, cytotoxic T-lymphocytes, and PD-L1, all patients were clustered into four subtypes: cluster 1 had a high level of immune dysfunction, cluster 2 had a low level of immune activation, cluster 3 had intense immune exclusion, and cluster 4 had a high level of immune activation and a moderate level of both dysfunction and exclusion signatures. Compared with cluster 1, the hazard ratios and 95% confidential intervals for overall survival were 0.63 (0.35-1.13) for cluster 2, 0.55 (0.29-1.03) for cluster 3, and 0.30 (0.14-0.64) for cluster 4 in multivariate Cox regression. Similar immune clustering and prognosis patterns were obtained upon validation in the GSE39582 cohort. In subgroup analysis, immune clustering was significantly associated with overall survival among stage I/II patients, microsatellite stable/instability-low patients, and patients not treated with adjuvant therapy.ConclusionsOur findings demonstrated that classifying CRC patients into different immune subtypes serves as a reliable prognosis predictor and may help to refine patient selection for personalized cancer immunotherapy.

Project description:Background: Tumor necrosis factor (TNF) family members play vital roles in cancer development and antitumor immune responses. However, the expression patterns, prognostic values, and immunological characteristics of TNF members in bladder carcinoma (BLCA) remain unclear. Methods: The training cohort, TCGA-BLCA, was downloaded from The Cancer Genome Atlas; another two Gene Expression Omnibus datasets (GSE13507 and GSE32894) and the Xiangya cohort (RNA-sequencing cohort collected from our hospital) were used as the external validation cohort. The least absolute shrinkage and selection operator (LASSO) algorithm and cross-validation were used to screen variables. Cox regression model and random survival forest (RSF) were used to develop the risk score, respectively. Then, we systematically correlated the TNF risk score with the tumor microenvironment (TME) cell infiltration, molecular subtypes of BLCA, and the potential value for predicting the efficacy of immunotherapy. Results: We developed two TNF-based patterns, named TNF cluster 1 and TNF cluster 2. TNF cluster 1 exhibited poorer survival outcome and an inflamed TME characteristic compared with TNF cluster 2. We then filtered out 196 differentially expressed genes between the two TNF clusters and applied the LASSO algorithm and cross-validation to screen out 22 genes to build the risk score. For risk score, we found that RSF exhibited higher efficacy than the Cox regression model, and we chose the risk score developed by RSF for the following analysis. BLCA patients in the higher risk score group showed significantly poorer survival outcomes. Moreover, these results could be validated in the external validation cohorts, including the GSE13507, GSE32894, and Xiangya cohorts. Then, we systematically correlated the risk score with TME cell infiltration and found that it was positively correlated with the infiltration of a majority of immune cells. Also, a higher risk score indicated a basal subtype of BLCA. Notably, the relationship between risk score, TME cell infiltration, and molecular subtypes could be validated in the Xiangya cohort. Conclusion: We developed and validated a robust TNF-based risk score, which could predict prognostic outcomes, TME, and molecular subtypes of BLCA. However, the value of risk score predicting the efficacy of immunotherapy needs further research.

Project description:Objective: Because of the modest immunotherapeutic response among ovarian carcinoma (OC) patients, it is significant to evaluate antitumor immune response and develop more effective precision immunotherapeutic regimens. Here, this study aimed to determine diverse immune subtypes of OC. Methods: This study curated the expression profiles of prognostic immunologically relevant genes and conducted consensus clustering analyses for determining immune subtypes among OC patients in TCGA cohort. With Boruta algorithm, characteristic genes were screened for conducting an immune scoring system through principal component analysis algorithm. The single-sample gene set enrichment analysis and ESTIAMTE methods were adopted for quantifying the immune infiltrates and responses to chemotherapeutic agents were estimated with pRRophetic algorithm. Two immunotherapeutic cohorts were used for investigating the efficacy of immune score in predicting therapeutic benefits. Results: Two immune subtypes were conducted among 377 OC patients. Immune subtype 2 was characterized by worse clinical prognosis, more frequent genetic variations and mutations, enhanced immune infiltrates, and increased expression of MHC molecules and programmed cell death protein 1/programmed death ligand 1 (PD-1/PD-L1). In total, 30 prognosis-relevant characteristic immune subtype-derived genes were identified for constructing the immune score of OC patients. High immune score was linked with more dismal prognosis, decreased immune infiltrations, and expression of MHC molecules. High immune score presented favorable sensitivity to doxorubicin and vinorelbine and reduced sensitivity to cisplatin. In addition, immune score possessed the potential in predicting benefits from anti-PD-1/PD-L1 therapy. Conclusion: Collectively, our findings propose two complex and diverse immune subtypes of OC. Quantitative assessment of immune subtypes in individual patients strengthens the understanding of tumor microenvironment features and promotes more effective immunotherapeutic regimens.

Dataset Information

Tumor Microenvironment Analysis Identified Subtypes Associated With the Prognosis and the Tumor Response to Immunotherapy in Bladder Cancer

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets