Identification of a prognostic gene signature based on an immunogenomic landscape analysis of bladder cancer

Development of a prognostic signature for immune-associated genes in bladder cancer and exploring potential drug findings

BACKGROUND

Bladder cancer, predominantly affecting men, is a prevalent malignancy of the urinary system. Although platinum-based chemotherapy has demonstrated certain enhancements in overall survival when compared to surgery alone, the efficacy of treatments is impeded by the unfavorable side effects of conventional chemotherapy medications. Nonetheless, immunotherapy exhibits potential in the treatment of bladder cancer.

METHODS

To create an immune-associated prognostic signature for bladder cancer, bioinformatics analyses were performed utilizing The Cancer Genome Atlas (TCGA) database in this study. By identifying differential gene expressions between the high-risk and low-risk groups, a potential therapeutic drug was predicted using the Connectivity Map database. Subsequently, the impact of this drug on the growth of T24 cells was validated through MTT assay and 3D cell culture techniques.

RESULTS

The signature included 1 immune-associated LncRNA (NR2F1-AS1) and 16 immune-associated mRNAs (DEFB133, RBP7, PDGFRA, CGB3, PDGFD, SCG2, ADCYAP1R1, OPRL1, PGR, PSMD1, TANK, PRDX1, ADIPOR2, S100A8, AHNAK, EGFR). Based on the assessment of risk scores, the patients were classified into cohorts of low-risk and high-risk individuals. The cohort with low risk demonstrated a considerably higher likelihood of survival in comparison to the group with high risk. Furthermore, variations in immune infiltration were noted among the two categories. Cephaeline, a possible medication, was discovered by analyzing variations in gene expression. It exhibited promise in suppressing the viability and growth of T24 bladder cancer cells.

CONCLUSION

The novel predictive pattern allows for efficient categorization of patients with bladder cancer, enabling focused and rigorous treatment for those expected to have a worse prognosis. The discovery of a possible curative medication establishes a basis for forthcoming immunotherapy trials in bladder cancer.

A transcriptome study of p53-pathway related prognostic gene signature set in bladder cancer

p53 pathway is important in tumorigenesis. However, no study has been performed to specifically investigate the role of p53 pathway genes in bladder cancer (BLCA). In this study, transcriptomics data of muscle invasive bladder cancer patients (n = 411) from The Cancer Genome Atlas (TCGA) were investigated. Using the hallmark p53 pathway gene set, the Non-Negative Matrix factorization (NMF) analysis identified two subtypes (C1 and C2). Clinical, survival, and immunological analysis were done to validate distinct characteristics of the subtypes. Pathway enrichment analysis showed the subtype C1 with poor prognosis having enrichment in genes of the immunity related pathways, where C2 subtype with better prognosis being enriched in genes of the steroid synthesis and drug metabolism pathways. A signature gene set consisting of MDGA2, GNLY, GGT2, UGT2B4, DLX1, and DSC1 was created followed by a risk model. Their expressions were analyzed in RNA extracted from the blood and matched tumor tissues of BLCA patients (n = 10). DSC1 had significant difference of expression (p = 0.005) between the blood and tumor tissues in our BLCA samples. Contrary to the usual normal bladder tissue to blood ratio, DLX1 expression was lower (p = 0.02734) in tumor tissues than in blood. Being the first research of p53 pathway related signature gene set in bladder cancer, this study potentially has a substantial impact on the development of biomarkers for BLCA.

Expression pattern, tumor immune landscape, and prognostic value of N7‑methylguanosine regulators in bladder urothelial carcinoma

N7-Methylguanosine (m7G) modification is important in post-transcriptional regulation. dysregulation of m7G RNA modification has been reported to be markedly associated with cancer. However, its importance in bladder urothelial carcinoma (BLCA) remains poorly characterized. The present study systematically analyzed mRNA gene expression data and clinical information from The Cancer Genome Atlas and further constructed robust risk signatures for the four regulators of m7G RNA modification (nudix hydrolase 11, gem nuclear organelle-associated protein 5, eukaryotic translation initiation factor 3 subunit D and cytoplasmic FMR1 interacting protein 1). The differential expression and cell function of m7G-related genes in bladder cancer cells were verified by reverse transcription-quantitative PCR, Cell Counting Kit-8 and colony formation assays. The four-gene-based model could accurately predict the prognosis of BLCA. Nomogram-based clinical decisions had a higher net benefit compared with that of individual predictors. Through immune infiltration analysis, it was found that immune cell infiltration affected the prognosis of patients with BLCA. Finally, the present study identified potential therapeutics that differ between high and low-risk groups based on four genes. In summary, the current findings revealed an essential role for m7G RNA modification regulators in BLCA, and developed risk signatures as promising prognostic markers in patients with BLCA.

Clinical Significance and Tumor Microenvironment Characterization of a Novel Immune-Related Gene Signature in Bladder Cancer

Cancer immunotherapy plays a crucial role in bladder cancer (BC) progression. Increasing evidence has elucidated the clinicopathologic significance of the tumor microenvironment (TME) in predicting outcomes and therapeutic efficacy. This study sought to establish a comprehensive analysis of the immune-gene signature combined with TME to assist in BC prognosis. We selected sixteen immune-related genes (IRGs) after a weighted gene co-expression network and survival analysis. Enrichment analysis revealed that these IRGs were actively involved in Mitophagy and Renin secretion pathways. After multivariable COX analysis, the IRGPI comprising NCAM1, CNTN1, PTGIS, ADRB3, and ANLN was established to predict the overall survival of BC, which was validated in both TCGA and GSE13507 cohorts. In addition, a TME gene signature was developed for molecular and prognosis subtyping with unsupervised clustering, followed by a panoramic landscape characterization of BC. In summary, the IRGPI model developed in our study provided a valuable tool with an improved prognosis for BC.

An epithelial-mesenchymal transition-related mRNA signature associated with the prognosis, immune infiltration and therapeutic response of colon adenocarcinoma

Background: Epithelial-mesenchymal transition (EMT) is closely associated with cancer cell metastasis. Colon adenocarcinoma (COAD) is one of the most common malignancies in the world, and its metastasis leading to poor prognosis remains a challenge for clinicians. The purpose of this study was to explore the prognostic value of EMT-related genes (EMTRGs) by bioinformatics analysis and to develop a new EMTRGs prognostic signature for COAD. Methods: The TCGA-COAD dataset was downloaded from the TCGA portal as the training cohort, and the GSE17538 and GSE29621 datasets were obtained from the GEO database as the validation cohort. The best EMTRGs prognostic signature was constructed by differential expression analysis, Cox, and LASSO regression analysis. Gene set enrichment analysis (GSEA) is used to reveal pathways that are enriched in high-risk and low-risk groups. Differences in tumor immune cell levels were analyzed using microenvironmental cell population counter and single sample gene set enrichment analysis. Subclass mapping analysis and Genomics of Drug Sensitivity in Cancer were applied for prediction of immunotherapy response and chemotherapy response, respectively. Results: A total of 77 differentially expressed EMTRGs were identified in the TCGA-COAD cohort, and they were significantly associated with functions and pathways related to cancer cell metastasis, proliferation, and apoptosis. We constructed EMTRGs prognostic signature with COMP, MYL9, PCOLCE2, SCG2, and TIMP1 as new COAD prognostic biomarkers. The high-risk group had a poorer prognosis with enhanced immune cell infiltration. The GSEA demonstrated that the high-risk group was involved in "ECM Receptor Interaction," "WNT Signaling Pathway" and "Colorectal Cancer." Furthermore, patients with high risk scores may respond to anti-CTLA4 therapy and may be more resistant to targeted therapy agents BI 2536 and ABT-888. Conclusion: Together, we developed a new EMTRGs prognostic signature that can be an independent prognostic factor for COAD. This study has guiding implications for individualized counseling and treatment of COAD patients.

Distinct immune and inflammatory response patterns contribute to the identification of poor prognosis and advanced clinical characters in bladder cancer patients

Due to the molecular heterogeneity, most bladder cancer (BLCA) patients show no pathological responses to immunotherapy and chemotherapy yet suffer from their toxicity. This study identified and validated three distinct and stable molecular clusters of BLCA in cross-platform databases based on personalized immune and inflammatory characteristics. H&E-stained histopathology images confirmed the distinct infiltration of immune and inflammatory cells among clusters. Cluster-A was characterized by a favorable prognosis and low immune and inflammatory infiltration but showed the highest abundance of prognosis-related favorable immune cell and inflammatory activity. Cluster-B featured the worst prognosis and high immune infiltration, but numerous unfavorable immune cells exist. Cluster-C had a favorable prognosis and the highest immune and inflammatory infiltration. Based on machine learning, a highly precise predictive model (immune and inflammatory responses signature, IIRS), including FN1, IL10, MYC, CD247, and TLR2, was developed and validated to identify the high IIRS-score group that had a poor prognosis and advanced clinical characteristics. Compared to other published models, IIRS showed the highest AUC in 5 years of overall survival (OS) and a favorable predictive value in predicting 1- and 3- year OS. Moreover, IIRS showed an excellent performance in predicting immunotherapy and chemotherapy's response. According to immunohistochemistry and qRT-PCR, IIRS genes were differentially expressed between tumor tissues with corresponding normal or adjacent tissues. Finally, immunohistochemical and H&E-stained analyses were performed on the bladder tissues of 13 BLCA patients to further demonstrate that the IIRS score is a valid substitute for IIR patterns and can contribute to identifying patients with poor clinical and histopathology characteristics. In conclusion, we established a novel IIRS depicting an IIR pattern that could independently predict OS and acts as a highly precise predictive biomarker for advanced clinical characters and the responses to immunotherapy and chemotherapy.

Immune-related gene risk score predicting the effect of immunotherapy and prognosis in bladder cancer patients

Background: Immune checkpoint inhibitor therapy has changed the treatment model of metastatic bladder cancer. However, only approximately 20% of patients benefit from this therapy, and robust biomarkers to predict the effect of immunotherapy are still lacking. In this study, we aimed to investigate whether immune-related genes could be indicators for the prognosis of bladder cancer patients and the effect of immunotherapy.

Methods: Based on bladder cancer dataset from the Cancer Genome Atlas (TCGA) and GSE48075, 22 immune microenvironment-related cells were identified by CIBERSORT. After performing a series of bioinformatic and machine learning approaches, we identified distinct tumor microenvironment clusters and three bladder cancer specific immune-related genes (EGFR, OAS1 and MST1R). Then, we constructed immune-related gene risk score (IRGRS) by using the Cox regression method and validated it with the IMvigor210 dataset.

Results: IRGRS-high patients had a worse overall survival than IRGRS-low patients, which was consistent with the result in the IMvigor210 dataset. Comprehensive analysis shows that patients with high IRGRS scores are mainly enriched in basal/squamous type (Ba/Sq), and tumor metabolism-related pathways are more Active, with higher TP53 and RB1 gene mutation rates, lower CD4+/CD8+ T cell infiltration, higher M0 macrophage infiltration, and lower immunotherapy efficacy. In contrast, Patients with low IRGRS scores are mainly enriched in the luminal papillary type (LumP), which is associated with the activation of IL-17 and TNF signaling pathways, higher mutation rates of FGFR3 and CDKN1A genes, higher CD4+/CD8+ T cell infiltration content, and The level of M0 macrophage infiltration was relatively low, and the immunotherapy was more probably effective.

Conclusion: Our study constructed an IRGRS for bladder cancer and clarified the immune and molecular characteristics of IRGRS-defined subgroups of bladder cancer to investigate the association between IRGRS and its potential implications for prognosis and immunotherapy.

Values of OAS gene family in the expression signature, immune cell infiltration and prognosis of human bladder cancer

Background

Bladder cancer (BLCA) is one of the most common genitourinary malignancies in the world, but its pathogenic genes have not been fully identified and the treatment outcomes are still unsatisfactory. Although the members of 2', 5'-oligoadenylate synthetase (OAS) gene family are known involved in some tumorous biological processes, the roles of the OAS gene family in BLCA are still undetermined.

Methods

By combining vast bioinformatic datasets analyses of BLCA and the experimental verification on clinical BLCA specimen, we identified the expressions and biological functions of OAS gene family members in BLCA with comparison to normal bladder tissues.

Results

The expression levels of OAS gene family members were higher in BLCA than in normal bladder tissues. The expression levels of most OAS genes had correlations with genomic mutation and methylation, and with the infiltration levels of CD4 + T cells, CD8 + T cells, neutrophils, and dendritic cells in the microenvironment of BLCA. In addition, high expressions of OAS1, OAS2, OAS3, and OASL predicted better overall survival in BLCA patients.

Conclusions

The highly expressed OAS genes in BLCA can reflect immune cells infiltration in the tumor microenvironment and predict the better overall survival of BLCA, and thus may be considered as a signature of BLCA. The study provides new insights into the diagnosis, treatment, and prognosis of BLCA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10102-8.

Identification of immune subtypes to guide immunotherapy and targeted therapy in clear cell renal cell carcinoma

Accumulating pieces of evidence suggested that immunotypes may indicate the overall immune landscape in the tumor microenvironment, which were closely related to therapeutic response. The purpose of this study was to classify and define the immune subtypes of clear cell renal cell carcinoma (ccRCC), so as to authenticate the potential immune subtypes that respond to immunotherapy. Transcriptome expression profile and mutation profile data of ccRCC, as well as clinical characteristics used in this study were obtained from TCGA database. There were significant differences in the infiltration of immune cells, immune checkpoints, and antigens between ccRCC and para-cancerous tissues. According to immune components, patients with ccRCC were divided into three immune subtypes, with different clinical and molecular characteristics. Compared with other subtypes, IS2 showed cold immune phenotype, and was associated with better survival. IS1 represented complex immune populations and was associated with poor overall survival (OS) and progression free survival (PFS). Further analysis indicated that expression of immune checkpoints also differed among the three subtypes, and was abnormally up-regulated in IS3. Pathway enrichment analysis indicated that the mTOR signaling pathway was abnormally enriched in IS3, while the TGF_BETA, ANGIOGENESIS and receptor tyrosine kinase signaling pathways were abnormally enriched in IS2. Furthermore, there was an abnormal enrichment of the epithelial-to-mesenchymal transition (EMT) signaling pathway in IS1, which may be associated with a higher rate of metastasis. Finally, SCG2 was screened as a specific antigen of ccRCC, which was not only related to poor prognosis, but also significantly associated with immune cells and immune checkpoints. In conclusion, the immune subtypes of ccRCC may provide new insights into the tumor biology and the precise clinical management of this disease.

Integrated Analysis and Identification of Critical RNA-Binding Proteins in Bladder Cancer

Simple Summary

The role of RNA-binding proteins (RBPs) in bladder cancer (BC) remains unclear. Therefore, we analyzed the clinical information and RNA sequencing data from patients with BC and identified RBPs that may be promising predictors of BC.

Abstract

RBPs in the development and progression of BC remains unclear. Here, we elucidated the role of RBPs in predicting the survival of patients with BC. Clinical information and RNA sequencing data of the training and validation cohorts were downloaded from the Cancer Genome Atlas and Gene Expression Omnibus databases, respectively. Survival-related differentially expressed RBPs were identified using Cox regression analyses. A total of 113 upregulated and 54 downregulated RBPs were observed, with six showing prognostic values (AHNAK, MAP1B, LAMA2, P4HB, FASN, and GSDMB). In both the GSE32548 and GSE31684 datasets, patients with low-risk scores in survival-related six RBPs-based prognostic model showed longer overall survival than those with high-risk scores. AHNAK, MAP1B, P4HB, and FASN expression were significantly upregulated in both BC tissues and cell lines. BC tissues from high-risk group showed higher proportions of naive CD4+ T cells, M0 and M2 macrophages, and neutrophils and lower proportions of plasma cells, CD8+ T cells, and T-cell follicular helper compared to low-risk group. AHNAK knockdown significantly inhibited the proliferation, invasion, and migration of BC cells in vitro and inhibited the growth of subcutaneous tumors in vivo. We thus developed and functionally validated a novel six RBPs-based prognostic model for BC.

Immunogenomic Characteristics of Cell-Death-Associated Genes with Prognostic Implications in Bladder Cancer

Bladder cancer is one of the most common genitourinary malignant cancers worldwide. Cell death processes, including apoptosis, ferroptosis, and necrosis, provide novel clinical and immunological insights promoting the management of precision medicine. Therefore, this study aimed to evaluate the transcriptomic profile of signatures in cell death pathways with significant prognostic implications in patients with bladder cancer from multiple independent cohorts (n = 1999). First, genes involved in apoptosis (n = 19), ferroptosis (n = 31), and necrosis (n = 6) were analyzed to evaluate the prognostic implications in bladder cancer. Significant genes were included to establish the cell-death index (CDI) of 36 genes that distinguished patients according to high and low risks. Survival analysis using the Kaplan-Meier curves clustered patients based on overall survival (18.8 vs. 96.7 months; hazard model [HR] = 3.12, P<00001). Cox proportional hazard model was significantly associated with a higher risk of mortality using 10 external independent cohorts in patients with CDIhigh (HR = 1.31, 95% CI: 1.04-1.62). To explore immune parameters associated with CDI, microenvironment cell-population-counter algorithms indicated increased intratumoral heterogeneity and macrophage/monocyte infiltration and CD8+ T cells in patients with CDIhigh group. Besides, the CDIhigh group showed an increased expression of the following immune checkpoints: CD276, PD-L1, CTLA-4, and T-cell exhaustion signatures. Cytokine expression analysis revealed the highest association of IL-9R, IL-17A, IL-17F, GDF7, and IFNW1 with the high-risk group. In addition, 42 patients with BCa receiving immunotherapies were enrolled from a real-world cohort, and expression patterns of three CDI hub genes (DRD5, SCL2A14, and IGF1) were detected using immunohistochemical staining. Patients with triple-negative staining of tumor tissues had significantly higher tumor-associated macrophage abundance, PD-L1 expression, predicted immunocompromised microenvironment, and prominently progressive progression (HR = 4.316, P = 0.0028). In conclusion, this study highlights the immunoevasive tumor microenvironment characterized by the higher tumor-associated macrophage infiltration with the presence of immune checkpoint and T-cell exhaustion genes in patients with BCa at CDIhigh risk who might suffer progression and be more suitable to benefit from immune checkpoint inhibitors or other immunotherapies.

Overexpression of OAS1 Is Correlated With Poor Prognosis in Pancreatic Cancer

Background

OAS1 expression in pancreatic cancer has been confirmed by many studies. However, the prognostic value and mechanism of OAS1 in pancreatic cancer have not been analyzed.

Methods

The RNA-seq in pancreatic cancer were obtained by UCSC XENA and GEO database. In addition, immunohistochemical validation and analysis were performed using samples from the 900th hospital. The prognosis of OAS1 was evaluated by timeROC package, Cox regression analysis, and Kaplan-Meier survival curves. Then, the main functional and biological signaling pathways enrichment and its relationship with the abundance of immune cells were analyzed by bioinformatics.

Results

OAS1 was highly expressed in pancreatic cancer compared with normal pancreatic tissue. High OAS1 expression was associated with poor overall survival (p<0.05). The OAS1 was significantly correlated to TNM staging (p=0.014). The timeROC analysis showed that the AUC of OAS1 was 0.734 for 3-year OS. In addition, the expression of OAS1 was significantly correlated with the abundance of a variety of immune markers. GSEA showed that enhanced signaling pathways associated with OAS1 include Apoptosis, Notch signaling pathway, and P53 signaling pathway.

Conclusions

OAS1 is a valuable prognostic factor in pancreatic cancer. Moreover, it may be a potential immunotherapeutic target.

Identification of an immune gene-associated prognostic signature in patients with bladder cancer

A deeper understanding of the interaction between tumor cell and the immune microenvironment in bladder cancer may help select predictive and prognostic biomarkers. The current study aims to construct a prognostic signature for bladder cancer by analysis of molecular characteristics, as well as tumor-immune interactions. RNA-sequencing and clinical information from bladder cancer patients were downloaded from the TCGA database. The single sample Gene Sets Enrichment Analysis (ssGSEA) and Cell type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) were employed to separate the samples into two clusters. Lasso Cox regression was performed to construct an immune gene signature for bladder cancer. The correlation between key target genes of immune checkpoint blockade and the prognostic signature was also analyzed. Dataset from Gene Expression Omnibus (GEO) was retrieved for validation. Two immunophenotypes and immunological characteristics were identified, and a 17-immune gene signature was constructed to provide an independent prognostic signature for bladder cancer. The signature was verified through external validation and correlated with genomic characteristics and clinicopathologic features. Finally, a nomogram was generated from the clinical characteristics and immune signature. Our study reveals a tumor-immune microenvironment signature useful for prognosis in bladder cancer. The results provide information on the potential development of treatment strategies for bladder cancer patients. Prospective studies are warranted to validate the prognostic capability of this model, but these data highlight the role of the microenvironment in the clinical outcome of patients.

Prognosis Risk Model Based on Pyroptosis-Related lncRNAs for Bladder Cancer

Objective. Bladder cancer (BC) is the most common malignancy in the urinary system and is prone to recurrence and metastasis. Pyroptosis is a kind of cell necrosis that is triggered by the gasdermin protein family. lncRNAs are noncoding RNAs that are more than 200 nucleotides long. Both pyroptosis and lncRNAs are associated with tumor development and progression. This study is aimed at exploring and establishing a prognostic signature of BC based on pyroptosis-related lncRNAs. Methods. In this study, The Cancer Genome Atlas (TCGA) database provided us with the RNA sequencing transcriptome data of bladder cancer patients, and we identified differentially expressed pyroptosis-related lncRNAs in bladder cancer. Then, the prognostic significance of these lncRNAs was assessed using univariate Cox regression analysis and LASSO regression analysis. Subsequently, 4 pyroptosis-related lncRNAs, namely, AL121652.1, AL161729.4, AC007128.1, and AC124312.3, were identified by multivariate Cox regression analysis, thus constructing the prognostic risk model. Then, we compared the levels of immune infiltration, differences in cell function, immune checkpoints, and m6A-related gene expression levels between the high- and low-risk groups. Result. Patients were divided into low-risk or high-risk groups based on the median risk score. Kaplan–Meier survival analysis indicated that the overall survival of bladder cancer patients in the low-risk group was substantially superior to that in the high-risk group ( $p<0.001$ ). The receiver operating characteristic (ROC) curve further confirmed the credibility of our model. Moreover, gene set enrichment analysis (GSEA) indicated that these were different signal pathways significantly enriched between the two groups. Immune infiltration, immune checkpoint, and N6-methyladenosine-related gene analysis also reflected that there were notable differences between the two groups. Conclusion. Therefore, this prognostic risk model is based on the level of pyroptotic lncRNAs, which is conducive to individualized assessment of the risk of patients and provides a reference for clinical treatment. This will also help provide insights into the prognosis and treatment of bladder cancer.

SCG2 is a Prognostic Biomarker Associated With Immune Infiltration and Macrophage Polarization in Colorectal Cancer

Colorectal cancer (CRC) is the second most lethal malignancy around the world. Limited efficacy of immunotherapy creates an urgent need for development of novel treatment targets. Secretogranin II (SCG2) is a member of the chromogranin family of acidic secretory proteins, has a role in tumor microenvironment (TME) of lung adenocarcinoma and bladder cancer. Besides, SCG2 is a stroma-related gene in CRC, its potential function in regulating tumor immune infiltration of CRC needs to be fully elucidated. In this study, we used western blot, real-time PCR, immunofluorescence and public databases to evaluate SCG2 expression levels and distribution. Survival analysis and functional enrichment analysis were performed. We examined TME and tumor infiltrating immune cells using ESTIMATE and CIBERSORT algorithm. The results showed that SCG2 expression was significantly decreased in CRC tumor tissues, and differentially distributed between tumor and adjacent normal tissues. SCG2 was an independent prognostic predictor in CRC. High expression of SCG2 correlated with poor survival and advanced clinical stage in CRC patients. SCG2 might regulate multiple tumor- and immune-related pathways in CRC, influence tumor immunity by regulating infiltration of immune cells and macrophage polarization in CRC.

Construction and verification of a prognostic risk model based on immunogenomic landscape analysis of bladder caner

This study was designed to construct a prognostic risk model to predict prognosis and immunotherapy response of bladder cancer (BCa) patinets. 350 differential expressed immune-related genes (DEIRGs) were obtained according to the transcriptome profiling and immune-related genes from the Cancer Genome Atlas (TCGA) database and ImmPort database, respectively. A prognostic risk model was constructed based on 15 hub genes through univariate, multivariate, and LASSO Cox regression analyses. The area under the receiver operating characteristic (ROC) curve was 0.743, indicating the superiority of the model. The scatter plot showed that as the risk score increased, the overall survival decreased significantly. In addition, all results were internally verified by the TCGA cohort. The model showed that the higher the grade, clinical stage, and TNM stage of BCa, the higher the risk score of patients. The tumor mutation burden of the low-risk group was generally higher than that of the high-risk group. Immune cell infiltration analysis showed that CD8 T cells, naive CD4 T cells, follicular helper T cells and M0 Macrophage were significantly different between the two groups. Several key immune checkpoint genes were found to be significantly different between the two groups, such as CTLA4, PD-L1, CD47, CD276, CXCL8, and HAVCR2/TIM3. Finally, the analysis of immunotherapy revealed that the efficacy of CTLA4 or PD1 blockers alone was better in the low-risk group than in the high-risk group. Taken together, we developed and validated a prognostic risk model based on 15 hub genes, which performed well in predicting prognosis and immunotherapy response of BCa patients.

Scoring System Based on RNA Modification Writer-Related Genes to Predict Overall Survival and Therapeutic Response in Bladder Cancer

Introduction

It’s widely reported the “writer” enzymes mediated RNA adenosine modifications which is known as a crucial mechanism of epigenetic regulation in development of tumor and the immunologic response in many kinds of cancers. However, the potential roles of these writer genes in the progression of bladder cancer (BLCA) remain unclear.

Materials and Methods

We comprehensively described the alterations of 26 RNA modification writer genes in BLCA from the genetic and transcriptional fields and identified writer-related genes from four independent datasets. Utilizing least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression, we constructed a ten writer-related gene signature. After that, we confirmed the predictive and prognostic value of this signature on another six independent datasets and established a nomogram to forecast the overall survival (OS) and mortality odds of BLCA patients clinically.

Results

The writer-related genes signature showed good performance in predicting the OS for BLCA patients. Moreover, the writer-related gene signature was related to EMT-related pathways and immune characteristics. Furthermore, the immune cell infiltration levels of CD8 T cells, cytotoxic cells, M1/2 macrophage cells and tumor mutation burden might be able to predict which patients will benefit from immunotherapy. This could also be reflected by the writer-related gene signature.

Conclusions

This signature might play an important role in precision individualized immunotherapy. The present work highlights the crucial clinical implications of RNA modifications and may help developing individualized therapeutic strategies for patients with BLCA.

MINDY1 promotes bladder cancer progression by stabilizing YAP

Background

Bladder cancer is one of the most commonly diagnosed urological malignant tumor. The Hippo tumor suppressor pathway is highly conserved in mammals and plays an important role in carcinogenesis. YAP is one of major key effectors of the Hippo pathway. However, the mechanism supporting abnormal YAP expression in bladder cancer remains to be characterized.

Methods

Western blot was used to measure the expression of MINDY1 and YAP, while the YAP target genes were measured by real-time PCR. CCK8 assay was used to detect the cell viability. The xeno-graft tumor model was used for in vivo study. Protein stability assay was used to detect YAP protein degradation. Immuno-precipitation assay was used to detect the interaction domain between MINDY1 and YAP. The ubiquitin-based Immuno-precipitation assays were used to detect the specific ubiquitination manner happened on YAP.

Results

In the present study, we identified MINDY1, a DUB enzyme in the motif interacting with ubiquitin-containing novel DUB family, as a bona fide deubiquitylase of YAP in bladder cancer. MINDY1 was shown to interact with, deubiquitylate, and stabilize YAP in a deubiquitylation activity-dependent manner. MINDY1 depletion significantly decreased bladder cancer cell proliferation. The effects induced by MINDY1 depletion could be rescued by further YAP overexpression. Depletion of MINDY1 decreased the YAP protein level and the expression of YAP/TEAD target genes in bladder cancer, including CTGF, ANKRD1 and CYR61.

Conclusion

In general, our findings establish a previously undocumented catalytic role for MINDY1 as a deubiquitinating enzyme of YAP and provides a possible target for the therapy of bladder cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02095-4.

Identification of a prognostic signature based on immune-related genes in bladder cancer

Bladder cancer (BLCA) has a high incidence and recurrence rate, and the effect of immunotherapy varies from person to person. Immune-related genes (IRGs) have been shown to be associated with immunotherapy and prognosis in many other cancers, but their role in immunogenic BLCA is less well defined. In this study, we constructed an eight-IRG risk model, which demonstrated strong prognostic and immunotherapeutic predictive power. The signature was significantly related to tumor clinicopathological characteristics, tumor class, immune cell infiltration and mutation status. Additionally, a nomogram containing the risk score and other potential risk factors could effectively predict the long-term overall survival probability of BLCA patients. The enriched mechanisms identified by gene set enrichment analysis suggested that the reason why this signature can accurately distinguish high- and low-risk populations may be closely related to the different degrees of innate immune response and T cell activation in different patients.

STAT1 is a key gene in a gene regulatory network related to immune phenotypes in bladder cancer: An integrative analysis of multi-omics data

The immunophenotype of bladder cancer plays a pivotal role in the prognosis of cancer, but the effect of different epigenetic factors on different immunophenotypes in bladder tumours remains unclear. This study used multi-omics data analysis to provide molecular basis support for different immune phenotypes. Unsupervised cluster analysis revealed distinct subclusters with higher (subcluster B2) or lower cytotoxic immune phenotypes (subcluster A1) related to PD-L1 and IFNG expression. Mutational landscape analyses showed that the mutation level of TP53 in subcluster B1 was highest than other subclusters, and subcluster B1 had a lower frequency of concurrent mutation than subcluster A2. A total of 2364 differentially expressed genes were identified between subclusters A2 and B1, and the main functions of the up-regulated genes in subcluster B1 were enriched in the activation of T cells and other related pathways. We found that STAT1 was a key gene in a gene regulatory network related to immune phenotypes in bladder cancer. Finally, we constructed a prognostic prediction model by LASSO Cox regression which could distinguish high-risk and low-risk cases significantly. In conclusion, the present study addressed a field synopsis between genetic and epigenetic events in immune phenotypes of bladder cancer.