Construction and validation of a metabolism-associated gene signature for predicting the prognosis, immune landscape, and drug sensitivity in bladder cancer.

Tumor Metabolism is strongly correlated with prognosis. Nevertheless, the prognostic and therapeutic value of metabolic-associated genes in BCa patients has not been fully elucidated. First, in this study, metabolism-related differential expressed genes DEGs with prognostic value in BCa were determined. Through the consensus clustering algorithm, we identified two molecular clusters with significantly different clinicopathological features and survival prognosis. Next, a novel metabolism-related prognostic model was established. Its reliable predictive performance in BCa was verified by multiple external datasets. Multivariate Cox analysis exhibited that risk score were independent prognostic factors. Interestingly, GSEA enrichment analysis of GO, KEGG, and Hallmark gene sets showed that the biological processes and pathways associated with ECM and collagen binding in the high-risk group were significantly enriched. Notely, the model was also significantly correlated with drug sensitivity, immune cell infiltration, and immunotherapy efficacy prediction by the wilcox rank test and chi-square test. Based on the 7 immune infiltration algorithm, we found that Neutrophils, Myeloid dendritic cells, M2 macrophages, Cancer-associated fibroblasts, etc., were more concentrated in the high-risk group. Additionally, in the IMvigor210, GSE111636, GSE176307, or our Truce01 (registration number NCT04730219) cohorts, the expression levels of multiple model genes were significantly correlated with objective responses to anti-PD-1/anti-PD-L1 immunotherapy. Finally, the expression of interested model genes were verified in 10 pairs of BCa tissues and para-carcinoma tissues by the HPA and real-time fluorescent quantitative PCR. Altogether, the signature established and validated by us has high predictive power for the prognosis, immunotherapy responsiveness, and chemotherapy sensitivity of BCa.

Shen C ，Bi Y ，Chai W ，Zhang Z ，Yang S ，Liu Y ，Wu Z ，Peng F ，Fan Z ，Hu H ... -  《BMC Medical Genomics》

Identification and validation of a dysregulated TME-related gene signature for predicting prognosis, and immunological properties in bladder cancer.

During tumor growth, tumor cells interact with their tumor microenvironment (TME) resulting in the development of heterogeneous tumors that promote tumor occurrence and progression. Recently, there has been extensive attention on TME as a possible therapeutic target for cancers. However, an accurate TME-related prediction model is urgently needed to aid in the assessment of patients' prognoses and therapeutic value, and to assist in clinical decision-making. As such, this study aimed to develop and validate a new prognostic model based on TME-associated genes for BC patients. Transcriptome data and clinical information for BC patients were extracted from The Cancer Genome Atlas (TCGA) database. Gene Expression Omnibus (GEO) and IMvigor210 databases, along with the MSigDB, were utilized to identify genes associated with TMEs (TMRGs). A consensus clustering approach was used to identify molecular clusters associated with TMEs. LASSO Cox regression analysis was conducted to establish a prognostic TMRG-related signature, with verifications being successfully conducted internally and externally. Gene ontology (GO), KEGG, and single-sample gene set enrichment analyses (ssGSEA) were performed to investigate the underlying mechanisms. The potential response to ICB therapy was estimated using the Tumor Immune Dysfunction and Exclusion (TIDE) algorithm and Immunophenoscore (IPS). Additionally, it was found that the expression level of certain genes in the model was significantly correlated with objective responses to anti-PD-1 or anti-PD-L1 treatment in the IMvigor210, GSE111636, GSE176307, or Truce01 (registration number NCT04730219) cohorts. Finally, real-time PCR validation was performed on 10 paired tissue samples, and in vitro cytological experiments were also conducted on BC cell lines. In BC patients, 133 genes differentially expressed that were associated with prognosis in TME. Consensus clustering analysis revealed three distinct clinicopathological characteristics and survival outcomes. A novel prognostic model based on nine TMRGs (including C3orf62, DPYSL2, GZMA, SERPINB3, RHCG, PTPRR, STMN3, TMPRSS4, COMP) was identified, and a TMEscore for OS prediction was constructed, with its reliable predictive performance in BC patients being validated. MultiCox analysis showed that the risk score was an independent prognostic factor. A nomogram was developed to facilitate the clinical viability of TMEscore. Based on GO and KEGG enrichment analyses, biological processes related to ECM and collagen binding were significantly enriched among high-risk individuals. In addition, the low-risk group, characterized by a higher number of infiltrating CD8+ T cells and a lower burden of tumor mutations, demonstrated a longer survival time. Our study also found that TMEscore correlated with drug susceptibility, immune cell infiltration, and the prediction of immunotherapy efficacy. Lastly, we identified SERPINB3 as significantly promoting BC cells migration and invasion through differential expression validation and in vitro phenotypic experiments. Our study developed a prognostic model based on nine TMRGs that accurately and stably predicted survival, guiding individual treatment for patients with BC, and providing new therapeutic strategies for the disease.

Shen C ，Chai W ，Han J ，Zhang Z ，Liu X ，Yang S ，Wang Y ，Wang D ，Wan F ，Fan Z ，Hu H ... -  《Frontiers in Immunology》

Identification of a dysregulated CircRNA-associated gene signature for predicting prognosis, immune landscape, and drug candidates in bladder cancer.

Increasing evidences have demonstrated that circular RNA (circRNAs) plays a an essential regulatory role in initiation, progression and immunotherapy resistance of various cancers. However, circRNAs have rarely been studied in bladder cancer (BCa). The purpose of this research is to explore new circRNAs and their potential mechanisms in BCa. A novel ceRNA-regulated network, including 87 differentially expressed circRNAs (DE-circRNAs), 126 DE-miRNAs, and 217 DE-mRNAs was constructed to better understanding the biological processes using Cytoscape 3.7.1 based on our previously high-throughput circRNA sequencing and five GEO datasets. Subsequently, five randomly selected circRNAs (upregulated circ_0001681; downregulated circ_0000643, circ_0001798, circ_0006117 and circ_0067900) in 20 pairs of BCa and paracancerous tissues were confirmed using qRT-PCR. Functional analysis results determined that 772 GO functions and 32 KEGG pathways were enriched in the ceRNA network. Ten genes (PFKFB4, EDNRA, GSN, GAS1, PAPPA, DTL, TGFBI, PRSS8, RGS1 and TCF4) were selected for signature construction among the ceRNA network. The Human Protein Atlas (HPA) expression of these genes were consistent with the above sequencing data. Notably, the model was validated in multiple external datasets (GSE13507, GSE31684, GSE48075, IMvigor210 and GSE32894). The immune-infiltration was evaluated by 7 published algorithms (i.e., TIMER, CIBERSORT, CIBERSORT-ABS, QUANTISEQ, MCPCOUNTER, XCELL and EPIC). Next, Correlations between riskscore or risk groups and clinicopathological data, overall survival, recognized immunoregulatory cells or common chemotherapeutic agents of BCa patients were performed using wilcox rank test, chi-square test, cox regression and spearman's correlation analysis; and, these results are significant. According to R package "GSVA" and "clusterProfiler", the most significantly enriched HALLMARK and KEGG pathway was separately the 'Epithelial Mesenchymal Transition' and 'Ecm Receptor Interaction' in the high- vs. low-risk group. Additionally, the functional experiments in vitro also revealed that the overexpression of has_circ_0067900 significantly impaired the proliferation, migration, and invasion capacities of BCa cells. Collectively, the results of the current study provide a novel landscape of circRNA-associated ceRNA-regulated network in BCa. The ceRNA-associated gene model which was constructed presented a high predictive performance for the prognosis, immunotherapeutic responsiveness, and chemotherapeutic sensitivity of BCa. And, has_circ_0067900 was originally proposed as tumor suppressor for patients with BCa.

Shen C ，Li Z ，Zhang Y ，Zhang Z ，Wu Z ，Da L ，Yang S ，Wang Z ，Zhang Y ，Qie Y ，Zhao G ，Lin Y ，Huang S ，Zhou M ，Hu H ... -  《Frontiers in Oncology》

Therapeutic Benefits and Prognostic Value of a Model Based on 7 Immune-associated Genes in Bladder Cancer.

Cao M ，Cao Y ，Xue S ，Zhang Q ，Zhang H ，Xue W ... -  《ALTERNATIVE THERAPIES IN HEALTH AND MEDICINE》

Construction and validation of a bladder cancer risk model based on autophagy-related genes.

Autophagy has an important association with tumorigenesis, progression, and prognosis. However, the mechanism of autophagy-regulated genes on the risk prognosis of bladder cancer (BC) patients has not been fully elucidated yet. In this study, we created a prognostic model of BC risk based on autophagy-related genes, which further illustrates the value of genes associated with autophagy in the treatment of BC. We first downloaded human autophagy-associated genes and BC datasets from Human Autophagy Database and The Cancer Genome Atlas (TCGA) database, and finally obtained differential prognosis-associated genes for autophagy by univariate regression analysis and differential analysis of cancer versus normal tissues. Subsequently, we downloaded two datasets from Gene Expression Omnibus (GEO), GSE31684 and GSE15307, to expand the total number of samples. Based on these genes, we distinguished the molecular subtypes (C1, C2) and gene classes (A, B) of BC by consistent clustering analysis. Using the genes merged from TCGA and the two GEO datasets, we conducted least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis to obtain risk genes and construct autophagy-related risk prediction models. The accuracy of this risk prediction model was assessed by receiver operating characteristic (ROC) and calibration curves, and then nomograms were constructed to predict the survival of bladder cancer patients at 1, 3, and 5 years, respectively. According to the median value of the risk score, we divided BC samples into the high- and low-risk groups. Kaplan-Meier (K-M) survival analysis was performed to compare survival differences between subgroups. Then, we used single sample gene set enrichment analysis (ssGSEA) for immune cell infiltration abundance, immune checkpoint genes, immunotherapy response, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, and tumor mutation burden (TMB) analysis for different subgroups. We also applied quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry (IHC) techniques to verify the expression of these six genes in the model. Finally, we chose the IMvigor210 dataset for external validation. Six risk genes associated with autophagy (SPOCD1, FKBP10, NAT8B, LDLR, STMN3, and ANXA2) were finally screened by LASSO regression algorithm and multivariate Cox regression analysis. ROC and calibration curves showed that the model established was accurate and reliable. Univariate and multivariate regression analyses were used to verify that the risk model was an independent predictor. K-M survival analysis indicated that patients in the high-risk group had significantly worse overall survival than those in the low-risk group. Analysis by algorithms such as correlation analysis, gene set variation analysis (GSVA), and ssGSEA showed that differences in immune microenvironment, enrichment of multiple biologically active pathways, TMB, immune checkpoint genes, and human leukocyte antigens (HLAs) were observed in the different risk groups. Then, we constructed nomograms that predicted the 1-, 3-, and 5-year survival rates of different BC patients. In addition, we screened nine sensitive chemotherapeutic drugs using the correlation between the obtained expression status of risk genes and drug sensitivity results. Finally, the external dataset IMvigor210 verified that the model is reliable and efficient. We established an autophagy-related risk prognostic model that is accurate and reliable, which lays the foundation for future personalized treatment of bladder cancer.

Shen C ，Yan Y ，Yang S ，Wang Z ，Wu Z ，Li Z ，Zhang Z ，Lin Y ，Li P ，Hu H ... -  《-》