Integrative analysis of 1q23.3 copy-number gain in metastatic urothelial carcinoma

Comprehensive analysis of tumor immune-related gene signature for predicting prognosis, immunotherapy, and drug sensitivity in bladder urothelial carcinoma

Background

Bladder urothelial carcinoma (BLCA) is globally recognized as a prevalent malignancy. Its treatment remains challenging due to the extensive morbidity, high mortality rates, and compromised quality of life from postoperative complications and the lack of specific molecular targets. Our aim was to establish a prognostic model to evaluate the prognostic significance, assess immunotherapy responses, and determine drug susceptibility in patients with BLCA.

Methods

From The Cancer Genome Atlas (TCGA) datasets, we obtained BLCA clinical details and expression data of immune-related genes. These data were analyzed using R and related packages. Differential expression analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, weighted gene co-expression network analysis (WGCNA), univariate and multivariate Cox regression analysis, prognostic analysis, model establishment and evaluation, gene set variation analysis (GSVA), immune function and checkpoint analysis, immunotherapy response prediction, and prediction of drug sensitivity were conducted.

Results

A total of 11 differentially expressed immune genes (DEIGs) were selected to establish the bladder carcinoma immune-related gene signature for BLCA prognosis prediction. In both the training and testing groups, the high-risk cohort showed a lower overall survival (OS) than the low-risk cohort. The area under the receiver operating characteristic curve (AUC) was 0.712 in the training group and 0.631 in the testing group, highlighting its predictive capacity. In the external validation datasets GSE39281 and IMvigor210, the OS of the high-risk group was significantly lower than that of the low-risk group, with AUC values of 0.609 and 0.563, respectively. Patients in the training group were categorized into low- and high-risk groups based on the bladder carcinoma immune gene signature (BCIGS) median risk score. GSVA showed 21 KEGG pathways positively correlated with model risk scores. The high-risk group presented with elevated stromal score, immune score, ESTIMATE score, and T cell exclusion score level. Conversely, the low-risk group displayed heightened cytotoxic T-lymphocyte antigen 4 (CTLA4) expression, indicative of a better response to immune checkpoint inhibitors (ICIs). Notably, significant disparities were found in immune subtypes, immune-related function, and immune-related survival between the two risk groups. The AUC values of our model are 0.765 and 0.660, respectively, surpassing those of other models, such as the tumor inflammation signature (TIS), tumor immune dysfunction and exclusion (TIDE), and various clinical factors. We also presented a nomogram, with the AUCs for predicting 1-, 2-, and 3-year OS at 0.727, 0.772, and 0.765 respectively, suggesting the signature's robust predictive power. Finally, 20 small molecular compounds were identified, with the TW.37 drug's half maximum inhibitory concentration (IC50) value difference being the most pronounced between the high- and low-risk patient groups, indicating its potential as a treatment option.

Conclusions

Our constructed immune-related gene signature model forecasts BLCA patient prognosis and potentially guides individualized immunotherapy and chemotherapeutic drug choices.

Establish TIIC signature score based the machine learning fusion in bladder cancer

BACKGROUND

Bladder cancer is a prevalent malignant tumor with high heterogeneity. Current treatments, such as transurethral resection of bladder tumor (TURBT) and intravesical Bacillus Calmette-Guérin (BCG) therapy, still have limitations, with approximately 30% of non-muscle-invasive bladder cancer (NMIBC) progressing to muscle-invasive bladder cancer (MIBC), and a substantial number of MIBC patients experiencing recurrence after surgery. Immunotherapy has shown potential benefits, but accurate prediction of its prognostic effects remains challenging.

METHODS

We analyzed bladder cancer RNA-seq data and clinical information from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and used various machine learning algorithms to screen for feature RNAs related to tumor-infiltrating immune cells (TIICs) from single-cell data. Based on these RNAs, we established a TIIC signature score and evaluated its relationship with overall survival (OS) and immunotherapy response in bladder cancer patients.

RESULTS

The study identified 171 TIIC-RNAs and selected 11 TIIC-RNAs with prognostic value through survival analysis. The TIIC signature score established using a machine learning fusion method was significantly associated with OS and showed good predictive performance in different datasets. Additionally, the signature score was negatively correlated with immunotherapy response, with patients with low TIIC feature scores showing better survival outcomes after immunotherapy. Further biological functional analysis revealed a close association between the TIIC signature score and immune regulation processes, cellular metabolism, and genetic variations.

CONCLUSION

This study successfully constructed and validated an RNA signature scoring system based on tumor-infiltrating immune cell (TIIC) features, which can effectively predict OS and the effectiveness of immunotherapy in bladder cancer patients.

Advancing bladder cancer management: development of a prognostic model and personalized therapy

Background

Bladder cancer (BLCA) was recognized as a significant public health challenge due to its high incidence and mortality rates. The influence of molecular subtypes on treatment outcomes was well-acknowledged, necessitating further exploration of their characterization and application. This study was aimed at enhancing the understanding of BLCA by mapping its molecular heterogeneity and developing a robust prognostic model using single-cell and bulk RNA sequencing data. Additionally, immunological characteristics and personalized treatment strategies were investigated through the risk score.

Methods

Single-cell RNA sequencing (scRNA-seq) data from GSE135337 and bulk RNA-seq data from several sources, including GSE13507, GSE31684, GSE32894, GSE69795, and TCGA-BLCA, were utilized. Molecular subtypes, particularly the basal-squamous (Ba/Sq) subtype associated with poor prognosis, were identified. A prognostic model was constructed using LASSO and Cox regression analyses focused on genes linked with the Ba/Sq subtype. this model was validated across internal and external datasets to ensure predictive accuracy. High- and low-risk groups based on the risk score derived from TCGA-BLCA data were analyzed to examine their immune-related molecular profiles and treatment responses.

Results

Six molecular subtypes were identified, with the Ba/Sq subtype being consistently associated with poor prognosis. The prognostic model, based on basal-squamous subtype-related genes (BSSRGs), was shown to have strong predictive performance across diverse clinical settings with AUC values at 1, 3, and 5 years indicating robust predictability in training, testing, and entire datasets. Analysis of the different risk groups revealed distinct immune infiltration and microenvironments. Generally higher tumor mutation burden (TMB) scores and lower tumor immune dysfunction and exclusion (TIDE) scores were exhibited by the low-risk group, suggesting varied potentials for systemic drug response between the groups. Finally, significant differences in potential systemic drug response rates were also observed between risk groups.

Conclusions

The study introduced and validated a new prognostic model for BLCA based on BSSRGs, which was proven effective in prognosis prediction. The potential for personalized therapy, optimized by patient stratification and immune profiling, was highlighted by our risk score, aiming to improve treatment efficacy. This approach was promised to offer significant advancements in managing BLCA, tailoring treatments based on detailed molecular and immunological insights.

NECTIN4 Amplification Is Frequent in Solid Tumors and Predicts Enfortumab Vedotin Response in Metastatic Urothelial Cancer

PURPOSE

The anti-NECTIN4 antibody-drug conjugate enfortumab vedotin (EV) is approved for patients with metastatic urothelial cancer (mUC). However, durable benefit is only achieved in a small, yet uncharacterized patient subset. NECTIN4 is located on chromosome 1q23.3, and 1q23.3 gains represent frequent copy number variations (CNVs) in urothelial cancer. Here, we aimed to evaluate NECTIN4 amplifications as a genomic biomarker to predict EV response in patients with mUC.

MATERIALS AND METHODS

We established a NECTIN4-specific fluorescence in situ hybridization (FISH) assay to assess the predictive value of NECTIN4 CNVs in a multicenter EV-treated mUC patient cohort (mUC-EV, n = 108). CNVs were correlated with membranous NECTIN4 protein expression, EV treatment responses, and outcomes. We also assessed the prognostic value of NECTIN4 CNVs measured in metastatic biopsies of non-EV-treated mUC (mUC-non-EV, n = 103). Furthermore, we queried The Cancer Genome Atlas (TCGA) data sets (10,712 patients across 32 cancer types) for NECTIN4 CNVs.

RESULTS

NECTIN4 amplifications are frequent genomic events in muscle-invasive bladder cancer (TCGA bladder cancer data set: approximately 17%) and mUC (approximately 26% in our mUC cohorts). In mUC-EV, NECTIN4 amplification represents a stable genomic alteration during metastatic progression and associates with enhanced membranous NECTIN4 protein expression. Ninety-six percent (27 of 28) of patients with NECTIN4 amplifications demonstrated objective responses to EV compared with 32% (24 of 74) in the nonamplified subgroup (P < .001). In multivariable Cox analysis adjusted for age, sex, and Bellmunt risk factors, NECTIN4 amplifications led to a 92% risk reduction for death (hazard ratio, 0.08 [95% CI, 0.02 to 0.34]; P < .001). In the mUC-non-EV, NECTIN4 amplifications were not associated with outcomes. TCGA Pan-Cancer analysis demonstrated that NECTIN4 amplifications occur frequently in other cancers, for example, in 5%-10% of breast and lung cancers.

CONCLUSION

NECTIN4 amplifications are genomic predictors of EV responses and long-term survival in patients with mUC.

Tumor immune dysfunction and exclusion subtypes in bladder cancer and pan-cancer: a novel molecular subtyping strategy and immunotherapeutic prediction model

BACKGROUND

Molecular subtyping is expected to enable precise treatment. However, reliable subtyping strategies for clinical application remains defective and controversial. Given the significance of tumor immune dysfunction and exclusion (TIDE), we aimed to develop a novel TIDE-based subtyping strategy to guide personalized immunotherapy in the bladder cancer (BC).

METHODS

Transcriptome data of BC was used to evaluate the heterogeneity and the status of TIDE patterns. Subsequently, consensus clustering was applied to classify BC patients based on TIDE marker-genes. Patients' clinicopathological, molecular features and signaling pathways of the different TIDE subtypes were well characterized. We also utilize the deconvolution algorithms to analyze the tumor microenvironment, and further explore the sensitivity and mechanisms of each subtype to immunotherapy. Furthermore, BC patient clinical information, real-world BC samples and urine samples were collected for the validation of our findings, which were used for RNA-seq analysis, H&E staining, immunohistochemistry and immunofluorescence staining, and enzyme-linked immunosorbent assay. Finally, we also explored the conservation of our novel TIDE subtypes in pan-cancers.

RESULTS

We identified 69 TIDE biomarker genes and classified BC samples into three subtypes using consensus clustering. Subtype I showed the lowest TIDE status and malignancy with the best prognosis and highest sensitivity to immune checkpoint blockade (ICB) treatment, which was enriched of metabolic related signaling pathways. Subtype III represented the highest TIDE status and malignancy with the poorest prognosis and resistance to ICB treatment, resulting from its inhibitory immune microenvironment and T cell terminal exhaustion. Subtype II was in a transitional state with intermediate TIDE level, malignancy, and prognosis. We further confirmed the existence and characteristics of our novel TIDE subtypes using real-world BC samples and collected patient clinical data. This subtyping method was proved to be more efficient than previous known methods in identifying non-responders to immunotherapy. We also propose that combining our TIDE subtypes with known biomarkers can potentially improve the sensitivity and specificity of these biomarkers. Moreover, besides guiding ICB treatment, this classification approach can assist in selecting the frontline or recommended drugs. Finally, we confirmed that the TIDE subtypes are conserved across the pan-tumors.

CONCLUSIONS

Our novel TIDE-based subtyping method can serve as a powerful clinical tool for BC and pan-cancer patients, and potentially guiding personalized therapy decisions for selecting potential beneficiaries and excluding resistant patients of ICB therapy.

Investigating the impact of regulatory B cells and regulatory B cell-related genes on bladder cancer progression and immunotherapeutic sensitivity

BACKGROUND

Regulatory B cells (Bregs), a specialized subset of B cells that modulate immune responses and maintain immune tolerance in malignant tumors, have not been extensively investigated in the context of bladder cancer (BLCA). This study aims to elucidate the roles of Bregs and Breg-related genes in BLCA.

METHODS

We assessed Breg infiltration levels in 34 pairs of BLCA and corresponding paracancerous tissues using immunohistochemical staining. We conducted transwell and wound healing assays to evaluate the impact of Bregs on the malignant phenotype of SW780 and T24 cells. Breg-related genes were identified through gene sets and transcriptional analysis. The TCGA-BLCA cohort served as the training set, while the IMvigor210 and 5 GEO cohorts were used as external validation sets. We employed LASSO regression and random forest for feature selection and developed a risk signature using Cox regression. Primary validation of the risk signature was performed through immunohistochemical staining and RT-qPCR experiments using the 34 local BLCA samples. Additionally, we employed transfection assays and flow cytometry to investigate Breg expansion ability and immunosuppressive functions.

RESULTS

Breg levels in BLCA tissues were significantly elevated compared to paracancerous tissues (P < 0.05) and positively correlated with tumor malignancy (P < 0.05). Co-incubation of SW780 and T24 cells with Bregs resulted in enhanced invasion and migration abilities (all P < 0.05). We identified 27 Breg-related genes, including CD96, OAS1, and CSH1, which were integrated into the risk signature. This signature demonstrated robust prognostic classification across the 6 cohorts (pooled HR = 2.25, 95% CI = 1.52-3.33). Moreover, the signature exhibited positive associations with advanced tumor stage (P < 0.001) and Breg infiltration ratios (P < 0.05) in the local samples. Furthermore, the signature successfully predicted immunotherapeutic sensitivity in three cohorts (all P < 0.05). Knockdown of CSH1 in B cells increased Breg phenotype and enhanced suppressive ability against CD8 + T cells (all P < 0.05).

CONCLUSIONS

Bregs play a pro-tumor role in the development of BLCA. The Breg-related gene signature established in this study holds great potential as a valuable tool for evaluating prognosis and predicting immunotherapeutic response in BLCA patients.

Collaborating single-cell and bulk RNA sequencing for comprehensive characterization of the intratumor heterogeneity and prognostic model development for bladder cancer

INTRODUCTION

Gaining a deeper insight into the single-cell RNA sequencing (scRNA-seq) results of bladder cancer (BLCA) provides a transcriptomic profiling of individual cancer cells, which may disclose the molecular mechanisms involved in BLCA carcinogenesis.

METHODS

scRNA data were obtained from GSE169379 dataset. We used the InferCNV software to determine the copy number variant (CNV) with normal epithelial cells serving as the reference, and performed the pseudo-timing analysis on subsets of epithelial cell using Monocle3 software. Transcription factor analysis was conducted using the Dorothea software. Intercellular communication analysis was performed using the Liana software. Cox analysis and LASSO regression were applied to establish a prognostic model.

RESULTS

We investigated the heterogeneity of tumors in four distinct cell types of BLCA cancer, namely immune cells, endothelial cells, epithelial cells, and fibroblasts. We evaluated the transcription factor activity of different immune cells in BLCA and identified significant enrichment of TCF7 and TBX21 in CD8+ T cells. Additionally, we identified two distinct subtypes of cancer-associated fibroblasts (CAFs), namely iCAFs and myoCAFs, which exhibited distinct communication patterns. Using sub-cluster and cell trajectory analyses, we identified different states of normal-to-malignant cell transformation in epithelial cells. TF analysis further revealed high activation of MYC and SOX2 in tumor cells. Finally, we identified five model genes (SLCO3A1, ANXA1, TENM3, EHBP1, LSAMP) for the development of a prognostic model, which demonstrated high effectiveness in stratifying patients across seven different cohorts.

CONCLUSIONS

We have developed a prognostic model that has demonstrated significant efficacy in stratifying patients with BLCA.

Single-cell RNA sequencing reveals the cellular and molecular characteristics of high-grade and metastatic bladder cancer

PURPOSE

Metastatic bladder cancer (BC) has the highest somatic mutation frequency and recurrence rate of all tumors. However, the cellular and molecular characteristics of BC remain unclear.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) on the samples of paracancerous normal tissue (PNT), primary tumor (PT) and lymph node metastasis (LNM). The proportions and gene expression profiles of different cell types in the tumor microenvironment (TME) were investigated.

RESULTS

In total, 50,158 cells were classified into six populations. Malignant cells of PT and LNM exhibited large mutant DNA fragments, while the cell phenotypes and gene expression profiles differed during differentiation. Metastasis was associated with a poorer prognosis than PT. Tumor-associated stromal cells and inhibitory immune cells were the main cell populations in PT and LNM. Cell-cell communication analysis revealed the roles of signaling pathways of inflammatory cancer-associated fibroblast (iCAF) and tumor-associated macrophage (TAM) in exhaustion of T cells. In addition, iCAF may recruit TAM to promote formation of the TME earlier than the differentiation of tumor cells.

CONCLUSION

This study through scRNA-seq enhanced our understanding of new features about the cellular and molecular similarities and differences of high-grade and metastatic bladder cancer, which might provide potential therapeutic targets in future treatment.

Elucidating the Associated Biological Function and Clinical Significance of RHOJ Expression in Urothelial Carcinoma

Urothelial cancer, a common urinary system malignancy, often presents treatment challenges due to metastasis and chemotherapy side effects. Angiogenesis, crucial for tumor growth, has become a target for drug development. This study explores the expression, prognostic value, and clinical correlation of RHOJ in the TCGA BLCA, GSE31684, and GSE32894 datasets. We identify common differentially expressed genes across these databases and utilize g:Profiler and Cytoscape ClueGO for functional assessment. Further, we perform a gene set enrichment analysis (GSEA) using Hallmark gene sets and use the imsig package for immune cell infiltration analysis. Our analysis indicates that RHOJ expression levels significantly impact survival rates, tumor progression, and immune response in urothelial tumors. High RHOJ expression correlated with poor prognosis, advanced disease stages, and an increase in monocyte population within the tumor microenvironment. This aligns with current literature indicating a key role of immune infiltration in bladder cancer progression and treatment response. Moreover, the GSEA and imsig results further suggest a potential mechanistic link between RHOJ expression and immune-related pathways. Considering the increasing emphasis on immunotherapeutic strategies in bladder cancer management, our findings on RHOJ's potential as a diagnostic biomarker and its association with immune response open new avenues for therapeutic interventions.

Development and experimental validation of a folate metabolism-related gene signature to predict the prognosis and immunotherapeutic sensitivity in bladder cancer

Folate metabolism is critical for the maintenance of genomic stability due to its regulatory ability to methylation, nucleotide metabolism, and reduction capabilities in cancer cells. However, the prognostic value of folate metabolism-related genes has not been clarified, especially in bladder cancer (BLCA). 91 folate metabolism-related genes were retrieved from the public database. TCGA-BLCA cohort, obtained from the Cancer Genome Atlas, was selected for training, while GSE13507, GSE31684, and GSE32894, downloaded from the Gene Expression Omnibus, and 35 BLCA samples collected from the local hospital were used for external validation. Through genomic difference detection, protein-protein interaction network analysis, LASSO regression, and Cox regression, a three-gene signature, including ATIC, INS, and MTHFD1L, was constructed. The signature was a reliable prognosis predictor across multiple independent cohorts (pooled hazard ratio = 2.79, 95% confidence interval = 1.79-4.33). The signature was associated with the BLCA malignant degree, which was validated in the local clinical samples (P < 0.01) and multiple cell lines (all P < 0.05). Additionally, the TIDE algorithm, GSE111636 cohort, and IMvigor210 cohort indicated that the signature was a promising tool to evaluate the immunotherapeutic response. Collectively, a folate metabolism-related gene signature was constructed to predict the prognosis and immunotherapeutic sensitivity in BLCA, which was verified in multiple large-scale cohorts, clinical samples, and cellular experiments, providing novel insights into the biological mechanisms.

A robust gene prognostic index composed of GZMB, IRF1, and TP63 can stratify the risk of two metastatic urothelial carcinoma cohorts based on immune checkpoint blockade therapy

BACKGROUND

Immune checkpoint blockade (ICB) therapy has become a first-line treatment option for metastatic urothelial carcinoma (mUC) patients who do not meet the criteria of cisplatin. Still, only a few people can benefit from it, so useful predictive markers are needed.

METHODS

Download the ICB-based mUC and chemotherapy-based bladder cancer cohorts, and extract the expression data of pyroptosis-related genes (PRG). The LASSO algorithm was used to construct the PRG prognostic index (PRGPI) in the mUC cohort, and we verified the prognostic ability of PRGPI in two mUC and two bladder cancer cohorts.

RESULTS

Most of the PRG in the mUC cohort were immune-activated genes, and a few were immunosuppressive genes. The PRGPI composed of GZMB, IRF1, and TP63 can stratify the risk of mUC. In IMvigor210 and GSE176307 cohorts, the P-values of Kaplan Meier analysis was < 0.01 and 0.002, respectively. PRGPI could also predict ICB response, and the chi-square test of the two cohorts had P-values of 0.002 and 0.046, respectively. In addition, PRGPI can also predict the prognosis of two bladder cancer cohorts without ICB therapy. The PRGPI and the expression of PDCD1/CD274 had a high degree of synergistic correlation. The Low PRGPI group showed prominent characteristics of immune infiltration and was enriched in the immune signal activation pathway.

CONCLUSION

The PRGPI we constructed can effectively predict the treatment response and overall survival rate of mUC patients treated with ICB. The PRGPI can help mUC patients achieve individualized and accurate treatment in the future.

Homologous recombination mRNAs (RAD21, RAD50 and BARD1) have a potentially poor prognostic role in ERBB2-low bladder cancer patients

Human epidermal growth factor receptor 2 (HER2/ERBB2) factor is known to be implicated in many malignancies and the potential of it as a prognostic biomarker was reported years ago. Molecular subtypes of HER2/ERBB2 negative and positive with distinct clinical outcomes have been identified in recent years; however, it is still under investigation for bladder cancer. This study evaluates the biological and prognostic significance of RAD21, RAD50 and BARD1 (homologous recombination biomarkers) mRNA levels with ERBB2 low and high expression to explore their impact on bladder cancer patient survival and cancer aggressiveness. The expression of ERBB2, RAD21, RAD50 and BARD1 mRNA levels was assessed in The Cancer Genome Atlas (TCGA) bladder cancer dataset along with four validation cohorts. Outcome analysis was evaluated using disease-free survival (DFS) and overall survival (OS). Univariate and multivariate analysis were used to evaluate the relationship between RAD21, RAD50, BARD1 and ERBB2 expression and clinicopathological variables. A significant increase in mRNA expression levels of RAD21, RAD50 and BARD1 was noticed in ERBB2-low patients compared to ERBB2-high patients. This overexpression of the homologous recombination repair transcripts was associated with poor outcome in ERBB2-low tumors, not in ERBB2-high tumors. Furthermore, the combined expression of high RAD21/RAD50, high RAD21/BARD1 or high RAD50/BARD1 were significantly associated with worse DFS and a better outcome for those with low co-expression in the ERBB2-low cohort. High expression of either RAD21/RAD50 or RAD21/BARD1 in ERBB2-low cohort associated with higher chance of metastasis. In addition, gene expression of BARD1 alone or in combination with RAD50 acted as an independent prognostic factor for worst survival. The data presented in this study reveal a connection between RAD21, RAD50, BARD1 and ERBB2 and patient survival. Importantly, it provided novel findings and potential prognostic markers, particularly in ERBB2-low bladder cancer.

Construction and experimental validation of a B cell-related gene signature to predict the prognosis and immunotherapeutic sensitivity in bladder cancer

BACKGROUND

B cells are essential components of tumor microenvironment and exert important functions in anti-tumor immune response. However, the prognosis value of B cell-related genes in bladder cancer (BLCA) remains obscure.

MATERIALS AND METHODS

The infiltrating levels of B cells were measured via the CD20 staining in the local samples and the computational biology analyses in the TCGA-BLCA cohort. The single-cell RNA sequencing analysis, gene-pair strategy, LASSO regression, random forest, and Cox regression were used for B cell-related signature construction. TCGA-BLCA cohort was chosen as the training cohort, and three independent cohorts from GEO and the local cohort were used for external validation. 326 B cells were adopted to explore the association between the model and B cells' biological processes. TIDE algorithm and two BLCA cohorts receiving anti-PD1/PDL1 treatment were utilized to detect its predictive ability to the immunotherapeutic response.

RESULTS

High infiltration levels of B cells heralded favorable prognosis, both in the TCGA-BLCA cohort and the local cohort (all P < 0.05). A 5-gene-pair model was established and served as a significant prognosis predictor across multiple cohorts (pooled hazard ratio = 2.79, 95% confidence interval = 2.22-3.49). The model could evaluate the prognosis effectively in 21 of 33 cancer types (P < 0.05). The signature was negatively associated with B cells' activation, proliferation, and infiltrating levels, and could serve as a potential predictor of immunotherapeutic outcomes.

CONCLUSIONS

A B cell-related gene signature was constructed to predict the prognosis and immunotherapeutic sensitivity in BLCA, helping to guide the personalized treatment.

The role of E3 ubiquitin ligases and deubiquitinases in bladder cancer development and immunotherapy

Bladder cancer is one of the common malignant urothelial tumors. Post-translational modification (PTMs), including ubiquitination, acetylation, methylation, and phosphorylation, have been revealed to participate in bladder cancer initiation and progression. Ubiquitination is the common PTM, which is conducted by E1 ubiquitin-activating enzyme, E2 ubiquitin-conjugating enzyme and E3 ubiquitin-protein ligase. E3 ubiquitin ligases play a key role in bladder oncogenesis and progression and drug resistance in bladder cancer. Therefore, in this review, we summarize current knowledge regarding the functions of E3 ubiquitin ligases in bladder cancer development. Moreover, we provide the evidence of E3 ubiquitin ligases in regulation of immunotherapy in bladder cancer. Furthermore, we mention the multiple compounds that target E3 ubiquitin ligases to improve the therapy efficacy of bladder cancer. We hope our review can stimulate researchers and clinicians to investigate whether and how targeting E3 ubiquitin ligases acts a novel strategy for bladder cancer therapy.

Apoptosis related genes mediated molecular subtypes depict the hallmarks of the tumor microenvironment and guide immunotherapy in bladder cancer

Apoptosis has been discovered as a mechanism of cell death. The purpose of this study is to identify the diagnostic signature factors related to bladder cancer (BLCA) through apoptosis related genes (ARGs). Clinicopathological parameters and transcriptomics data of 1,440 BLCA patients were obtained from 7 datasets (GSE13507, GSE31684, GSE32548, GSE32894, GSE48075, TCGA-BLCA, and IMvigor210). We first identified prognosis-related ARGs in BLCA and used them to construct two ARGs molecular subtypes by using consensus clustering algorithm. By using principal component analysis algorithms, a ARGscore was constructed to quantify the index of individualized patient. High ARGscore correlated with progressive malignancy and poor outcomes in BLCA patients. High ARGscore was associated with higher immune cell, higher estimate scores, higher stromal scores, higher immune scores, higher immune checkpoint, and lower tumor purity, which was consistent with the "immunity tidal model theory". Preclinically, BLCA immunotherapy cohorts confirmed patients with low ARGscore demonstrated significant therapeutic advantages and clinical benefits. These findings contribute to our understanding of ARGs and immunotherapy in BLCA. The ARGscore is a potentially useful tool to predict the prognosis and immunotherapy in BLCA.

LINC01589 serves as a potential tumor-suppressor and immune-related biomarker in endometrial cancer: A review

Currently, increasing attention is being paid to biomarkers in endometrial cancer. Immune infiltration of the tumor microenvironment has been shown to significantly affect the overall survival (OS) of uterine corpus endometrial carcinoma (UCEC) patients. LINC01589 is a long non-coding RNA (lncRNA) that is rarely reported in cancer and is assumed to play a role in immune regulation. We therefore evaluated the role of LINC01589 in UCEC using the Cancer Genome Atlas (TCGA) database. We analyzed the expression of LINC01589 using the gene expression profiles of LINC01589 in the UCEC projects in TCGA. Comparisons between the differentially expressed genes (DEGs) of the cancer and adjacent normal tissues of the UCEC projects revealed that LINC01589 expression was decreased in UCEC tissues. A multivariate cox regression analysis indicated that LINC01589 upregulation could serve as an independent prognostic factor for survival. Furthermore, there was a positive correlation between LINC01589 expression and B cell, T cell, NK cell, monocytic lineage, and myeloid dendritic cell infiltration in UCEC patients. In addition, 5 clusters of hub genes were detected by comparison of different expression levels of LINC01589 in the UCEC groups. The analysis of the reactome pathway using gene set enrichment analysis (GSEA) revealed immune-related pathways, including CD22-mediated B cell receptor (BCR) regulation and antigen-activated BCRs, leading to the generation of second messengers and complement cascade pathways that were significantly enriched in the high LINC01589 expression group. Thus, LINC01589 may serve as a prognostic biomarker, as it is associated with immune infiltration in UCEC.

Prognostic value of COL10A1 and its correlation with tumor-infiltrating immune cells in urothelial bladder cancer: A comprehensive study based on bioinformatics and clinical analysis validation

Introduction

Bladder cancer (BLCA) is one of the most lethal diseases. COL10A1 is secreted small-chain collagen in the extracellular matrix associated with various tumors, including gastric, colon, breast, and lung cancer. However, the role of COL10A1 in BLCA remains unclear. This is the first research focusing on the prognostic value of COL10A1 in BLCA. In this research, we aimed to uncover the association between COL10A1 and the prognosis, as well as other clinicopathological parameters in BLCA.

Methods

We obtained gene expression profiles of BLCA and normal tissues from the TCGA, GEO, and ArrayExpress databases. Immunohistochemistry staining was performed to investigate the protein expression and prognostic value of COL10A1 in BLCA patients. GO and KEGG enrichment along with GSEA analyses were performed to reveal the biological functions and potential regulatory mechanisms of COL10A1 based on the gene co-expression network. We used the “maftools” R package to display the mutation profiles between the high and low COL10A1 groups. GIPIA2, TIMER, and CIBERSORT algorithms were utilized to explore the effect of COL10A1 on the tumor immune microenvironment.

Results

We found that COL10A1 was upregulated in the BLCA samples, and increased COL10A1 expression was related to poor overall survival. Functional annotation of 200 co-expressed genes positively correlated with COL10A1 expression, including GO, KEGG, and GSEA enrichment analyses, indicated that COL10A1 was basically involved in the extracellular matrix, protein modification, molecular binding, ECM-receptor interaction, protein digestion and absorption, focal adhesion, and PI3K-Akt signaling pathway. The most commonly mutated genes of BLCA were different between high and low COL10A1 groups. Tumor immune infiltrating analyses showed that COL10A1 might have an essential role in recruiting infiltrating immune cells and regulating immunity in BLCA, thus affecting prognosis. Finally, external datasets and biospecimens were used, and the results further validated the aberrant expression of COL10A1 in BLCA samples.

Conclusions

In conclusion, our study demonstrates that COL10A1 is an underlying prognostic and predictive biomarker in BLCA.

Comprehensive analysis of a novel RNA modifications-related model in the prognostic characterization, immune landscape and drug therapy of bladder cancer

Background: Bladder cancer (BCa) is the leading reason for death among genitourinary malignancies. RNA modifications in tumors closely link to the immune microenvironment. Our study aimed to propose a promising model associated with the "writer" enzymes of five primary RNA adenosine modifications (including m⁶A, m⁶A_m, m¹A, APA, and A-to-I editing), thus characterizing the clinical outcome, immune landscape and therapeutic efficacy of BCa. Methods: Unsupervised clustering was employed to categorize BCa into different RNA modification patterns based on gene expression profiles of 34 RNA modification "writers". The RNA modification "writers" score (RMS) signature composed of RNA phenotype-associated differentially expressed genes (DEGs) was established using the least absolute shrinkage and selection operator (LASSO), which was evaluated in meta-GEO (including eight independent GEO datasets) training cohort and the TCGA-BLCA validation cohort. The hub genes in the RMS model were determined via weighted gene co-expression network analysis (WGCNA) and were further validated using human specimen. The potential applicability of the RMS model in predicting the therapeutic responsiveness was assessed through the Genomics of Drug Sensitivity in Cancer database and multiple immunotherapy datasets. Results: Two distinct RNA modification patterns were determined among 1,410 BCa samples from a meta-GEO cohort, showing radically varying clinical outcomes and biological characteristics. The RMS model comprising 14 RNA modification phenotype-associated prognostic DEGs positively correlated with the unsatisfactory outcome of BCa patients in meta-GEO training cohort (HR = 3.00, 95% CI = 2.19-4.12) and TCGA-BLCA validation cohort (HR = 1.53, 95% CI = 1.13-2.09). The infiltration of immunosuppressive cells and the activation of EMT, angiogenesis, IL-6/JAK/STAT3 signaling were markedly enriched in RMS-high group. A nomogram exhibited high prognostic prediction accuracy, with a concordance index of 0.785. The therapeutic effect of chemotherapeutic agents and antibody-drug conjugates was significantly different between RMS-low and -high groups. The combination of the RMS model and conventional characteristics (TMB, TNB and PD-L1) achieved an optimal AUC value of 0.828 in differentiating responders from non-responders to immunotherapy. Conclusion: We conferred the first landscape of five forms of RNA modifications in BCa and emphasized the excellent power of an RNA modifications-related model in evaluating BCa prognosis and immune landscape.

Identification and validation of a novel prognostic model based on platinum Resistance-related genes in bladder cancer

BACKGROUND

The depth of response to platinum in urothelial neoplasm tissues varies greatly. Biomarkers that have practical value in prognosis stratification are increasingly needed. Our study aimed to select a set of BC (bladder cancer)-related genes involved in both platinum resistance and survival, then use these genes to establish the prognostic model.

MATERIALS AND METHODS

Platinum resistance-related DEGs (differentially expressed genes) and tumorigenesis-related DEGs were identified. Ten most predictive co-DEGs were acquired followed by building a risk score model. Survival analysis and ROC (receiver operating characteristic) plot were used to evaluate the predictive accuracy. Combined with age and tumor stages, a nomogram was generated to create a graphical representation of survival rates at 1-, 3-, 5-, and 8-year in BC patients. The prognostic performance was validated in three independent BC datasets with platinum-based chemotherapy. The potential mechanism was explored by enrichment analysis.

RESULTS

PPP2R2B, TSPAN7, ATAD3C, SYT15, SAPCD1, AKR1B1, TCHH, AKAP12, AGLN3, and IGF2 were selected for our prognostic model. Patients in high- and low-risk groups exhibited a significant survival difference with HR (hazard ratio) = 2.7 (p < 0.0001). The prognostic nomogram of predicting 3-year OS (overall survival) for BC patients could yield an AUC (area under the curve) of 0.819. In the external validation dataset, the risk score also has a robust predictive ability.

CONCLUSION

A prognostic model derived from platinum resistance-related genes was constructed, we confirmed its value in predicting platinum-based chemotherapy benefits and overall survival for BC patients. The model might assist in therapeutic decisions for bladder malignancy.

Comprehensive analysis of the prognostic value and functions of prefoldins in hepatocellular carcinoma

Prefoldins (PFDNs), a group of proteins known to be associated with cytoskeletal rearrangement, are involved in tumor progression in various cancer types. However, little is known about the roles of PFDNs in hepatocellular carcinoma (HCC). Herein, we investigated the transcriptional and survival data of PFDNs from The Cancer Genome Atlas (TCGA) database. Gene Ontology (GO), Gene Set Enrichment Analysis (GSEA), and single-sample gene set enrichment analysis (ssGSEA) were used to evaluate the potential functions of PFDN1/2/3/4. We also detected the expression of PFDN1/2/3/4 via immunohistochemistry (IHC), Western blotting, and real-time PCR in our clinical samples. We found that the PFDN family showed elevated expression in HCC tissues, while only PFDN1/2/3/4 were found to be significantly correlated with poor prognosis of patients with HCC in the TCGA database. Further investigation was associated with PFDN1-4. We found that the expression of PFDN1/2/3/4 was significantly associated with advanced clinicopathologic features. Apart from the TCGA database, IHC, real-time PCR, and immunoblotting identified the overexpression of PFDN1/2/3/4 in HCC tissues and HCC cell lines. Taken together, these results indicated that PFDN1/2/3/4 might be novel prognostic biomarkers and treatment targets for patients with HCC.

Insights Into the Properties, Biological Functions, and Regulation of USP21

Deubiquitylating enzymes (DUBs) antagonize ubiquitination by removing ubiquitin from their substrates. The role of DUBs in controlling various physiological and pathological processes has been extensively studied, and some members of DUBs have been identified as potential therapeutic targets in diseases ranging from tumors to neurodegeneration. Ubiquitin-specific protease 21 (USP21) is a member of the ubiquitin-specific protease family, the largest subfamily of DUBs. Although USP21 was discovered late and early research progress was slow, numerous studies in the last decade have gradually revealed the importance of USP21 in a wide variety of biological processes. In particular, the pro-carcinogenic effect of USP21 has been well elucidated in the last 2 years. In the present review, we provide a comprehensive overview of the current knowledge on USP21, including its properties, biological functions, pathophysiological roles, and cellular regulation. Limited pharmacological interventions for USP21 have also been introduced, highlighting the importance of developing novel and specific inhibitors targeting USP21.

CD93 orchestrates the tumor microenvironment and predicts the molecular subtype and therapy response of bladder cancer

BACKGROUND

CD93 is newly reported to normalize vasculature and attenuate pancreatic cancer therapy response, but its role in bladder cancer (BLCA) is unknown.

METHOD

The immunologic role of CD93 is analyzed across TCGA pan-cancers. The correlation between CD93 and BLCA clinical and tumor microenvironment features, predicted immunotherapy pathways, molecular subtypes, therapeutic signatures and mutation status was evaluated in TCGA-BLCA and other two BLCA cohorts. The impact of CD93 on immunotherapy response was validated by five real-world cohorts, and chemotherapy response was assessed with IC50. CD93-based risk model was constructed with LASSO regression and validated by seven independent cohorts.

RESULT

CD93 is positively correlated with immunomodulators, tumor-infiltrating lymphocytes (TILs) and immune checkpoints across pan-cancers. In BLCA, CD93 leads to higher T cell inflamed score and expression of immune checkpoints. However, CD93 is indicative of more aggressive clinical features, worse survival, more tumor-associated macrophages and regulatory T cells recruitment, less recognition and killing of cancer cells by T cells, lower predicted chemotherapy and immunotherapy response, which is further validated by immunotherapy cohorts (IMvigor210: 16.11% vs 29.53%; GSE176307: 15.56% vs 20.93%). Notably, CD93 correlates with enriched neuroendocrine subtype and epithelial-mesenchymal transition differentiation, while CD93-low group has enriched luminal subtype. Pathways including hypoxia and Wnt-β-catenin are enriched along with CD93 expression, and more frequent FGFR3 mutation is also observed. Lastly, the CD93-based risk model, validated by seven independent cohorts, is powerful in distinguishing the survival probability of BLCA (3-year AUC 0.808).

CONCLUSION

CD93 plays a critical role in tumor immune regulation. CD93 expression indicates more aggressive clinicopathological status and molecular subtypes of BLCA and worse therapy response, which implies that combing anti-CD93 therapy with immunotherapy (or chemotherapy) may be potentially beneficial for BLCA in clinical practice.

Classification of Muscle Invasive Bladder Cancer to Predict Prognosis of Patients Treated with Immunotherapy

Background

Recently, immunotherapies have been approved for advanced muscle invasive bladder cancer (MIBC) treatment, but only a small fraction of MIBC patients could achieve a durable drug response. Our study is aimed at identifying tumor microenvironment (TME) subtypes that have different immunotherapy response rates.

Methods

The mRNA expression profiles of MIBC samples from seven discovery datasets (GSE13507, GSE31684, GSE32548, GSE32894, GSE48075, GSE48276, and GSE69795) were analyzed to identify TME subtypes. The identified TME subtypes were then validated by an independent dataset (TCGA-MIBC). The subtype-related biomarkers were discovered using computational analyses and then utilized to establish a random forest predictive model. The associations of TME subtypes with immunotherapy therapeutic responses were investigated in a group of patients who had been treated with immunotherapy. A prognostic index model was constructed using the subtype-related biomarkers. Two nomograms were built by the subtype-related biomarkers or the clinical parameters.

Results

Two TME subtypes, including ECM-enriched class (EC) and immune-enriched class (IC), were found. EC was associated with greater extracellular matrix (ECM) pathways, and IC was correlated with immune pathways, respectively. Overall survival was significantly greater for tumors classified as IC, whereas the EC subtype had a worse prognosis. A total of nine genes (AKAP12, APOL3, CXCL13, CXCL9, GBP4, LRIG1, PEG3, PODN, and PTPRD) were selected by computational analyses to construct the random forest model. The area under the curve (AUC) values for this model were 0.827 and 0.767 in the testing and external validation datasets, respectively. Therapeutic response rates were greater in IC patients than in EC patients (28 percent vs. 18 percent). Patients with a high prognostic index had a poorer prognosis than those with a low prognostic index. The nomogram constructed from nine genes and stage achieved a C-index of 0.71.

Conclusion

The present investigation defined two distinct TME subtypes and developed models to assess immunotherapeutic treatment outcomes.

Pyroptosis-Related Signature Predicts Prognosis and Immunotherapy Efficacy in Muscle-Invasive Bladder Cancer

Pyroptosis has profound impacts on tumor cell proliferation, invasion, and metastasis and is of great clinical significance for different cancers. However, the role of pyroptosis in the progression and prognosis of muscle invasive bladder cancer (MIBC) remains poorly characterized. Here, we collected multicenter MIBC data and performed integrated analysis to dissect the role of pyroptosis in MIBC and provide an optimized treatment for this disease. Based on transcriptomic data, we developed a novel prognostic model named the pyroptosis-related gene score (PRGScore), which summarizes immunological features, genomic alterations, and clinical characteristics associated with the pyroptosis phenotype. Samples with high PRGScore showed enhancement in CD8⁺ T cell effector function, antigen processing machinery and immune checkpoint and better response to immunotherapy by programmed cell death 1 (PD-1) and programmed cell death ligand 1 (PD-L1) inhibitors, which indicates that PRGScore is a valuable signature in the identification of populations sensitive to immune checkpoint inhibitors. Collectively, our study provides insights into further research targeting pyroptosis and its tumor immune microenvironment (TME) and offers an opportunity to optimize the treatment of MIBC.

Prefoldin Function in Cellular Protein Homeostasis and Human Diseases

Cellular functions are largely performed by proteins. Defects in the production, folding, or removal of proteins from the cell lead to perturbations in cellular functions that can result in pathological conditions for the organism. In cells, molecular chaperones are part of a network of surveillance mechanisms that maintains a functional proteome. Chaperones are involved in the folding of newly synthesized polypeptides and assist in refolding misfolded proteins and guiding proteins for degradation. The present review focuses on the molecular co-chaperone prefoldin. Its canonical function in eukaryotes involves the transfer of newly synthesized polypeptides of cytoskeletal proteins to the tailless complex polypeptide 1 ring complex (TRiC/CCT) chaperonin which assists folding of the polypeptide chain in an energy-dependent manner. The canonical function of prefoldin is well established, but recent research suggests its broader function in the maintenance of protein homeostasis under physiological and pathological conditions. Interestingly, non-canonical functions were identified for the prefoldin complex and also for its individual subunits. We discuss the latest findings on the prefoldin complex and its subunits in the regulation of transcription and proteasome-dependent protein degradation and its role in neurological diseases, cancer, viral infections and rare anomalies.

Identification and Validation of Key Genes of Differential Correlations in Gastric Cancer

Background: Gastric cancer (GC) is aggressive cancer with a poor prognosis. Previously bulk transcriptome analysis was utilized to identify key genes correlated with the development, progression and prognosis of GC. However, due to the complexity of the genetic mutations, there is still an urgent need to recognize core genes in the regulatory network of GC.

Methods: Gene expression profiles (GSE66229) were retrieved from the GEO database. Weighted correlation network analysis (WGCNA) was employed to identify gene modules mostly correlated with GC carcinogenesis. R package ‘DiffCorr’ was applied to identify differentially correlated gene pairs in tumor and normal tissues. Cytoscape was adopted to construct and visualize the gene regulatory network.

Results: A total of 15 modules were detected in WGCNA analysis, among which three modules were significantly correlated with GC. Then genes in these modules were analyzed separately by “DiffCorr”. Multiple differentially correlated gene pairs were recognized and the network was visualized by the software Cytoscape. Moreover, GEMIN5 and PFDN2, which were rarely discussed in GC, were identified as key genes in the regulatory network and the differential expression was validated by real-time qPCR, WB and IHC in cell lines and GC patient tissues.

Conclusions: Our research has shed light on the carcinogenesis mechanism by revealing differentially correlated gene pairs during transition from normal to tumor. We believe the application of this network-based algorithm holds great potential in inferring relationships and detecting candidate biomarkers.

Combination Biomarker of Immune Checkpoints Predict Prognosis of Urothelial Carcinoma

In contrast to Western counties, the incidence of urothelial carcinoma (UC) remains mar-edly elevated in Taiwan. Regulatory T cells (Tregs) play a crucial role in limiting immune responses within the tumor microenvironment. To elucidate the relationship between immune checkpoints in the tumor immune microenvironment and UC progression, we utilize the Gene Expression Omnibus (GEO) to analyze a microarray obtained from 308 patients with UC. We observed that the expression level of CD276 or TIM-3 was positively correlated with late-stage UC and poor prognosis. Patients with simultaneously high CD276 and TIM-3 expression in tumors have significantly reduced both univariate and multivariate survival, indicating that mRNA levels of these immune checkpoints could be independent prognostic biomarkers for UC overall survival and recurrence. Our cohort study showed rare CD8+ cytotoxic T-cells and Tregs infiltration during early-stage UC-known as cold tumors. Approximately 30% of late-stage tumors exhibited highly infiltrated cytotoxic T cells with high PD-1 and FOXP3 expression, which implied that cytotoxic T cells were inhibited in the advanced UC microenvironment. Collectively, our findings provide a better prognosis prediction by combined immune checkpoint biomarkers and a basis for early-stage UC standard treatment to convert cold tumors into hot tumors, followed by immune checkpoint therapy.

Alterations of Chromatin Regulators in the Pathogenesis of Urinary Bladder Urothelial Carcinoma

Simple Summary

Urinary bladder cancer is one of the ten major cancers worldwide, with higher incidences in males, in smokers, and in highly industrialized countries. New therapies beyond cytotoxic chemotherapy are urgently needed to improve treatment of these tumors. A better understanding of the mechanisms underlying their development may help in this regard. Recently, it was discovered that a group of proteins regulating the state of chromatin and thus gene expression is exceptionally and frequently affected by gene mutations in bladder cancers. Altered function of these mutated chromatin regulators must therefore be fundamental in their development, but how and why is poorly understood. Here we review the current knowledge on changes in chromatin regulators and discuss their possible consequences for bladder cancer development and options for new therapies.

Abstract

Urothelial carcinoma (UC) is the most frequent histological type of cancer in the urinary bladder. Genomic changes in UC activate MAPK and PI3K/AKT signal transduction pathways, which increase cell proliferation and survival, interfere with cell cycle and checkpoint control, and prevent senescence. A more recently discovered additional category of genetic changes in UC affects chromatin regulators, including histone-modifying enzymes (KMT2C, KMT2D, KDM6A, EZH2), transcription cofactors (CREBBP, EP300), and components of the chromatin remodeling complex SWI/SNF (ARID1A, SMARCA4). It is not yet well understood how these changes contribute to the development and progression of UC. Therefore, we review here the emerging knowledge on genomic and gene expression alterations of chromatin regulators and their consequences for cell differentiation, cellular plasticity, and clonal expansion during UC pathogenesis. Our analysis identifies additional relevant chromatin regulators and suggests a model for urothelial carcinogenesis as a basis for further mechanistic studies and targeted therapy development.

Construction and Validation of an Autophagy-Related Prognostic Signature and a Nomogram for Bladder Cancer

Objective

Bladder cancer (BC) is one of the top ten cancers endangering human health but we still lack accurate tools for BC patients’ risk stratification. This study aimed to develop an autophagy-related signature that could predict the prognosis of BC. In order to provide clinical doctors with a visual tool that could precisely predict the survival probability of BC patients, we also attempted to establish a nomogram based on the risk signature.

Methods

We screened out autophagy-related genes (ARGs) combining weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) in BC. Based on the screened ARGs, we performed survival analysis and Cox regression analysis to identify potential prognostic biomarkers. A risk signature based on the prognostic ARGs by multivariate Cox regression analysis was established, which was validated by using seven datasets. To provide clinical doctors with a useful tool for survival possibility prediction, a nomogram assessed by the ARG-based signature and clinicopathological features was constructed, verified using four independent datasets.

Results

Three prognostic biomarkers including BOC (P = 0.008, HR = 1.104), FGF7(P = 0.030, HR = 1.066), and MAP1A (P = 0.001, HR = 1.173) were identified and validated. An autophagy-related risk signature was established and validated. This signature could act as an independent prognostic feature in patients with BC (P = 0.047, HR = 1.419). We then constructed two nomograms with and without ARG-based signature and subsequent analysis indicated that the nomogram with ARG signature showed high accuracy for overall survival probability prediction of patients with BC (C-index = 0.732, AUC = 0.816). These results proved that the ARG signature improved the clinical net benefit of the standard model based on clinicopathological features (age, pathologic stage).

Conclusions

Three ARGs were identified as prognosis biomarkers in BC. An ARG-based signature was established for the first time, showing strong potential for prognosis prediction in BC. This signature was proven to improve the clinical net benefit of the standard model. A nomogram was established using this signature, which could lead to more effective prognosis prediction for BC patients.

A comprehensive molecular characterization of the 8q22.2 region reveals the prognostic relevance of OSR2 mRNA in muscle invasive bladder cancer

Technological advances in molecular profiling have enabled the comprehensive identification of common regions of gene amplification on chromosomes (amplicons) in muscle invasive bladder cancer (MIBC). One such region is 8q22.2, which is largely unexplored in MIBC and could harbor genes with potential for outcome prediction or targeted therapy. To investigate the prognostic role of 8q22.2 and to compare different amplicon definitions, an in-silico analysis of 357 patients from The Cancer Genome Atlas, who underwent radical cystectomy for MIBC, was performed. Amplicons were generated using the GISTIC2.0 algorithm for copy number alterations (DNA_Amplicon) and z-score normalization for mRNA gene overexpression (RNA_Amplicon). Kaplan-Meier survival analysis, univariable, and multivariable Cox proportional hazard ratios were used to relate amplicons, genes, and clinical parameters to overall (OS) and disease-free survival (DFS). Analyses of the biological functions of 8q22.2 genes and genomic events in MIBC were performed to identify potential targets. Genes with prognostic significance from the in silico analysis were validated using RT-qPCR of MIBC tumor samples (n = 46). High 8q22.2 mRNA expression (RNA-AMP) was associated with lymph node metastases. Furthermore, 8q22.2 DNA and RNA amplified patients were more likely to show a luminal subtype (DNA_Amplicon_core: p = 0.029; RNA_Amplicon_core: p = 0.01). Overexpression of the 8q22.2 gene OSR2 predicted shortened DFS in univariable (HR [CI] 1.97 [1.2; 3.22]; p = 0.01) and multivariable in silico analysis (HR [CI] 1.91 [1.15; 3.16]; p = 0.01) and decreased OS (HR [CI] 6.25 [1.37; 28.38]; p = 0.0177) in RT-qPCR data analysis. Alterations in different levels of the 8q22.2 region are associated with manifestation of different clinical characteristics in MIBC. An in-depth comprehensive molecular characterization of genomic regions involved in cancer should include multiple genetic levels, such as DNA copy number alterations and mRNA gene expression, and could lead to a better molecular understanding. In this study, OSR2 is identified as a potential biomarker for survival prognosis.

The biology and rationale of targeting nectin-4 in urothelial carcinoma

Bladder cancer is the tenth most common cancer type worldwide. Urothelial carcinoma is the most common type of bladder cancer and accounts for 90% of bladder cancer cases in the USA and Europe. Novel approaches are needed to improve patient outcomes. Nectin-4 is a tumour-associated antigen found on the surface of most urothelial carcinoma cells. In the antibody-drug conjugate enfortumab vedotin, human anti-nectin-4 antibody is linked to the cytotoxic microtubule-disrupting agent monomethyl auristatin E. In ongoing phase I, II and III clinical trials, enfortumab vedotin has been evaluated as a monotherapy and in combination with a checkpoint inhibitor and/or chemotherapy in locally advanced and metastatic urothelial carcinoma. Encouraging data from the phase II study resulted in the FDA granting accelerated approval for enfortumab vedotin in December 2019 for patients with locally advanced or metastatic urothelial carcinoma who were previously treated with platinum and a checkpoint inhibitor therapy. Moreover, data from a phase I study led to the FDA granting breakthrough therapy designation to enfortumab vedotin combined with pembrolizumab as a first-line treatment in February 2020 for cisplatin-ineligible patients with locally advanced or metastatic urothelial carcinoma. Results of ongoing and future combination studies of enfortumab vedotin with immunotherapy and other novel agents are eagerly awaited.

The Role of Prefoldin and Its Subunits in Tumors and Their Application Prospects in Nanomedicine

Prefoldin (PFDN) is a hexameric chaperone complex that is widely found in eukaryotes and archaea and consists of six different subunits (PFDN1-6). Its main function is to transfer actin and tubulin monomers to the eukaryotic cell cytoplasmic chaperone protein (c-CPN) specific binding during the assembly of the cytoskeleton, to stabilize the newly synthesized peptides so that they can be folded correctly. The current study found that each subunit of PFDN has different functions, which are closely related to the occurrence, development and prognosis of tumors. However, the best characteristics of each subunit have not been fully affirmed. The connection between research and tumors can change the understanding of PFDN and further extend its potential prognostic role and structural function to cancer research and clinical practice. This article mainly reviews the role of canonical PFDN and its subunits in tumors and other diseases, and discusses the potential prospects of the unique structure and function of PFDN in nanomedicine.

The functions and mechanisms of prefoldin complex and prefoldin-subunits

The correct folding is a key process for a protein to acquire its functional structure and conformation. Prefoldin is a well-known chaperone protein that regulates the correct folding of proteins. Prefoldin plays a crucial role in the pathogenesis of common neurodegenerative diseases (Alzheimer's disease, Parkinson's disease, and Huntington's disease). The important role of prefoldin in emerging fields (such as nanoparticles, biomaterials) and tumors has attracted widespread attention. Also, each of the prefoldin subunits has different and independent functions from the prefoldin complex. It has abnormal expression in different tumors and plays an important role in tumorigenesis and development, especially c-Myc binding protein MM-1. MM-1 can inhibit the activity of c-Myc through various mechanisms to regulate tumor growth. Therefore, an in-depth analysis of the complex functions of prefoldin and their subunits is helpful to understand the mechanisms of protein misfolding and the pathogenesis of diseases caused by misfolded aggregation.

SWI/SNF chromatin remodeling complex and glucose metabolism are deregulated in advanced bladder cancer

Bladder cancer (BC) is a frequently diagnosed malignancy affecting predominantly adult and elderly populations. It is expected that due to the longer life time, BC will become even more frequent in the future; thus in consequence, it will represent serious health problem of older society part. The treatment of advanced BC is mostly ineffective due to its very aggressive behavior. So far, no effective targeted therapy is used for BC treatment. Here, we found that BC is characterized by lower protein levels of BRM, INI1, and BAF155 main subunits of SWI/SNF chromatin remodeling complex (CRC) which is involved in global control of gene expression and influences various important cellular processes like: cell cycle control, apoptosis, DNA repair, etc. Moreover, the expression of SMARCA2, a BRM encoding gene, strongly correlated with BC metastasis and expression of such metabolic genes as PKM2 and PRKAA1. Furthermore, the analysis of T24 and 5637 commonly used BC cell lines revealed different expression levels of metabolic genes including FBP1 gene encoding Frutose-1,6-Bisphosphatase, an enzyme controlling glycolysis flux and gluconeogenesis. The tested BC cell lines exhibited various molecular and metabolic alterations as well as differential glucose uptake, growth rate, and migration potential. We have shown that BRM subunit is involved in the transcriptional control of genes encoding metabolic enzymes. Moreover, we found that the FBP1 expression level and the SWI/SNF CRCs may serve as markers of molecular subtypes of BC. Collectively, this study may provide a new knowledge about the molecular and metabolic BC subtypes which likely will be of high importance for the clinic in the future.

Downregulation of ST3GAL5 is associated with muscle invasion, high grade and a poor prognosis in patients with bladder cancer

In patients with bladder cancer (BC), the association between ST3 β-galactoside α-2,3-sialyltransferase 5 (ST3GAL5) expression and clinical outcomes, particularly regarding muscle-invasive disease, high tumor grade and prognosis, remain unknown. In the present study, the expression of ST3GAL5 and its association with clinical outcomes in patients with BC was analyzed using various public bioinformatics databases. The difference in ST3GAL5 expression between BC and healthy bladder tissues was also evaluated using data from the Oncomine database, The Cancer Genome Atlas and Gene Expression Omnibus database. The differences in ST3GAL5 expression between muscle invasive BC (MIBC) and non-muscle invasive BC (NMIBC), and high- and low-grade BC were also analyzed. Furthermore, genes that were positively co-expressed with ST3GAL5 in patients with BC were identified from the intersection between the Oncomine, Gene Expression Profiling Interactive Analysis 2 and UALCAN databases. Enrichment analysis by Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Reactome pathway enrichment analyses and a gene-concept network was performed using R package. Gene set enrichment analysis was also performed to assess the signaling pathways influenced by the high and low expression of ST3GAL5 in BC. The results indicated that ST3GAL5 expression was significantly lower in BC tissues compared with normal bladder tissues (P<0.05). Furthermore, ST3GAL5 expression in MIBC and high-grade BC was significantly lower compared with NMIBC and low-grade BC (P<0.05), respectively. The results from Kaplan-Meier survival analysis result demonstrated that ST3GAL5 downregulation was associated with poor survival in patients with BC (P<0.05). Taken together, these findings suggested that ST3GAL5 may be considered as an anti-oncogene in BC, could represent a potential predictive and prognostic biomarker for BC and may be a molecular target for tumor therapy.

RNAseq expression patterns of canine invasive urothelial carcinoma reveal two distinct tumor clusters and shared regions of dysregulation with human bladder tumors

BACKGROUND

Invasive urothelial carcinoma (iUC) is highly similar between dogs and humans in terms of pathologic presentation, molecular subtypes, response to treatment and age at onset. Thus, the dog is an established and relevant model for testing and development of targeted drugs benefiting both canine and human patients. We sought to identify gene expression patterns associated with two primary types of canine iUC tumors: those that express a common somatic mutation in the BRAF gene, and those that do not.

METHODS

We performed RNAseq on tumor and normal tissues from pet dogs. Analysis of differential expression and clustering, and positional and individual expression was used to develop gene set enrichment profiles distinguishing iUC tumors with and without BRAFV595E mutations, as well as genomic regions harboring excessive numbers of dysregulated genes.

RESULTS

We identified two expression clusters that are defined by the presence/absence of a BRAFV595E (BRAFV600E in humans) somatic mutation. BRAFV595E tumors shared significantly more dysregulated genes than BRAF wild-type tumors, and vice versa, with 398 genes differentiating the two clusters. Key genes fall into clades of limited function: tissue development, cell cycle regulation, immune response, and membrane transport. The genomic site with highest number of dysregulated genes overall lies in a locus corresponding to human chromosome 8q24, a region frequently amplified in human urothelial cancers.

CONCLUSIONS

These data identify critical sets of genes that are differently regulated in association with an activating mutation in the MAPK/ERK pathway in canine iUC tumors. The experiments also highlight the value of the canine system in identifying expression patterns associated with a common, shared cancer.

Systematic screening of protein-coding gene expression identified HMMR as a potential independent indicator of unfavorable survival in patients with papillary muscle-invasive bladder cancer

Papillary and non-papillary are two histological patterns of bladder carcinogenesis and are considered as dual-track oncogenic pathways, which have different genetic alterations. The TCGA-bladder cancer (BLCA) database contains clinicopathological, genomic and survival data from over 400 muscle-invasive bladder cancer patients. In this study, using data from this database, we performed a systematic screening of gene expression to identify the protein-coding gene that might have prognostic value in papillary and non-papillary muscle-invasive bladder cancer (MIBC). The data of patients with primary MIBC in TCGA-BLCA was acquired from the UCSC Xena project (http://xena.ucsc.edu) for re-analysis. By setting |log2 fold change|≥2 and adjusted p value <0.01 as the screening criteria, we found 751 significantly dysregulated genes, including 183 overexpressed and 568 downregulated genes. HMMR was identified as a potential prognostic marker with unique expression. Multivariate analysis showed that its expression was an independent prognostic indicator of shorter progression-free survival (PFS) (HR: 1.400, 95%CI: 1.021-1.920, p = 0.037) in the papillary subtype. ENST00000393915.8 and ENST00000358715.3, two transcripts that contain all 18 exons and encode the full length of HMMR, were significantly upregulated in cancer tissues compared with normal bladder tissues. None of the 17 CpG sites in its DNA locus was relevant to HMMR expression. 26/403 (6.5%) MIBC cases had HMMR gene-level amplification, which was associated with upregulated HMMR expression compared with the copy-neutral and deletion groups. Gene set enrichment analysis (GSEA) in papillary MIBC found that the high HMMR expression group was associated with upregulated genes enriched in multiple gene sets with well-established role in BC development, including G2M checkpoint, E2 F Targets, Myc Targets V1, Myc Targets V2 and Glycolysis. Based on these findings, we infer that HMMR expression might be a specific prognostic marker in terms of PFS in papillary MIBC. DNA amplification might be an important mechanism of its elevation.

YWHAZ amplification/overexpression defines aggressive bladder cancer and contributes to chemo-/radio-resistance by suppressing caspase-mediated apoptosis

The objective of this study was to characterize the oncogenic actions of a recently identified cancer-associated gene YWHAZ (also named as 14-3-3 ζ/δ) in urothelial carcinomas of the urinary bladder (UCUB). A genome-wide study revealed YWHAZ to be involved in the amplicon at 8q22.3, and its genetic amplification was detected predominantly in muscle-invasive bladder cancer (MIBC). Immunohistochemical staining confirmed the association of YWHAZ overexpression with higher tumor stages, lymph node/vascular invasion, and mitotic activity. Univariate and multivariate analyses further indicated the prognostic potential of YWHAZ for more aggressive cancer types. Both gene set enrichment analysis and STRING network studies suggested involvement of YWHAZ in regulating caspase-mediated apoptosis. Ectopic expression of YWHAZ in bladder cells with low endogenous YWHAZ levels boosted cell resistance to doxorubicin and cisplatin, as well as to ionizing radiation. Conversely, YWHAZ-knockdown using specific shRNA in cells with high endogenous YWHAZ levels diminished survival activity, suppressing cell growth and increasing cell death. Our findings confirm the essential role played by YWHAZ in sustaining cell proliferation during chemo/radiotherapy. Treatments based on anti-YWHAZ strategies may thus be beneficial for UCUB patients overexpressing YWHAZ. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Collagen-rich airway smooth muscle cells are a metastatic niche for tumor colonization in the lung

Metastases account for the majority of cancer deaths. While certain steps of the metastatic cascade are well characterized, identification of targets to block this process remains a challenge. Host factors determining metastatic colonization to secondary organs are particularly important for exploration, as those might be shared among different cancer types. Here, we showed that bladder tumor cells expressing the collagen receptor, CD167a, responded to collagen I stimulation at the primary tumor to promote local invasion and utilized the same receptor to preferentially colonize at airway smooth muscle cells (ASMCs)—a rich source of collagen III in lung. Morphologically, COL3-CD167a-driven metastatic foci are uniquely distinct from typical lung alveolar metastatic lesions and exhibited activation of the CD167a-HSP90-Stat3 axis. Importantly, metastatic lung colonization could be abrogated using an investigational drug that attenuates Stat3 activity, implicating this seed-and-soil interaction as a therapeutic target for eliminating lung metastasis.

Collagen is a dynamic component of both the tumor and metastatic niche. Here, the authors show that airway smooth muscle cells are a collagen III rich niche bladder cancer cells expressing CD167a, and Stat3 is a downstream target for abrogating  these collagen III/CD167a-driven metastatic foci.

Four novel biomarkers for bladder cancer identified by weighted gene coexpression network analysis

Bladder cancer (BC) is one of the most malignancies in terms of incidence and recurrence worldwide. The aim of this study is to find out novel and prognostic biomarkers for patients with BC. First, we identified 258 differentially expressed genes by using GSE19915 from Gene Expression Omnibus database. Second, a total of 33 modules were identified by constructing a coexpression network by using weighted gene coexpression network analysis and yellow module was regarded as the key module. Furthermore, by constructing protein-protein interaction networks, we preliminarily picked out 13 genes. Among them, four hub genes (CCNB1, KIF4A, TPX2, and TRIP13) were eventually identified by using five different methods (survival analysis, one-way analysis of variance, the Spearman correlation analysis, receiver operating characteristic curve, and expression value comparison), which were significantly correlated with the prognosis of BC. The validation of transcriptional and translational levels made sense (based on Oncomine and The Human Protein Atlas database). Moreover, functional enrichment analysis suggested that all the hub genes played crucial roles in chromosome segregation, sister chromatid segregation, nuclear chromosome segregation, mitotic nuclear division, nuclear division, and organelle fission during cell mitosis. In addition, three of the hub genes (KIF4A, TPX2, and TRIP13) might be potential targets of cancer drugs according to the results of the genetical alteration. In conclusion, this study indicates that four hub genes have great predictive value for the prognosis of BC, and may contribute to the exploration of the further and more in-depth research of BC.

MCL1 and DEDD Promote Urothelial Carcinoma Progression

Focal amplification of chromosome 1q23.3 in patients with advanced primary or relapsed urothelial carcinomas is associated with poor survival. We interrogated chromosome 1q23.3 and the nearby focal amplicon 1q21.3, as both are associated with increased lymph node disease in patients with urothelial carcinoma. Specifically, we assessed whether the oncogene MCL1 that resides in 1q21.3 and the genes that reside in the 1q23.3 amplicon were required for the proliferation or survival of urothelial carcinoma. We observed that suppressing MCL1 or the death effector domain-containing protein (DEDD) in the cells that harbor amplifications of 1q21.3 or 1q23.3, respectively, inhibited cell proliferation. We also found that overexpression of MCL1 or DEDD increased anchorage independence growth in vitro and increased experimental metastasis in vivo in the nonamplified urothelial carcinoma cell line, RT112. The expression of MCL1 confers resistance to a range of apoptosis inducers, while the expression of DEDD led to resistance to TNFα-induced apoptosis. These observations identify MCL1 and DEDD as genes that contribute to aggressive urothelial carcinoma. IMPLICATIONS: These studies identify MCL1 and DEDD as genes that contribute to aggressive urothelial carcinomas.

The Role of COL5A2 in Patients With Muscle-Invasive Bladder Cancer: A Bioinformatics Analysis of Public Datasets Involving 787 Subjects and 29 Cell Lines

Bladder cancer (BC) is one of the most common malignancies. Two previous studies identified collagen type V alpha 2 (COL5A2) as a potential biomarker in BC, both are simple reanalysis of a single transcriptomic dataset without subgroup analysis for muscle-invasive BC (MIBC). We focused in MIBC patients and explored the role of COL5A2 from an integration perspective, using refined methodology covering individual participant data meta-analysis and bioinformatics analysis. Eight transcriptomic datasets of 787 MIBC patients (including one dataset containing genomic mutation information) and two drug sensitivity datasets of 29 cell lines in which more than 250 compounds were analyzed. We found subjects with increased COL5A2 gene expression exhibited poorer prognosis, and the power analysis confirmed adequate sample size. FGFR3 was the only gene differential mutated between the COL5A2 high and low expression groups. Differential expression and co-expression network analysis suggested potential association between COL5A2 expression and essential pathways involved in cancer invasion and dissemination, including cell adhesion, extracellular matrix organization, and epithelial-mesenchymal transition. Coordinately, analysis of drug screening datasets and gene-drug interaction also revealed COL5A2 expression linked to cell morphogenesis, angiogenesis, blood vessel development, and urogenital development. The utility and feasibility of COL5A2 for clinically applicable prognosis prediction and risk classification and the exact underlying molecular mechanism should be further investigated in subsequent studies.

A prospective study to examine the accuracies and efficacies of prediction systems for response to neoadjuvant chemotherapy for muscle invasive bladder cancer

The present study established systems to predict the chemo-sensitivity of muscle invasive bladder cancer (MIBC) for neoadjuvant chemotherapy (NAC) with methotrexate, vinblastine, doxorubicin plus cisplatin (M-VAC) and carboplatin plus gemcitabine (CaG) by analyzing microarray data. The primary aim of the study was to investigate whether the clinical response would increase by combining these prediction systems. Treatment of each MIBC case was allocated into M-VAC NAC, CaG NAC, surgery, or radiation therapy groups by their prediction score (PS), which was calculated using the designed chemo-sensitivity prediction system. The therapeutic effect of the present study was compared with the results of historical controls (n=76 patients) whose treatments were not allocated using the chemo-sensitivity prediction system. In addition, the overall survival between the predicted to be responder (positive PS) group and predicted to be non-responder (negative PS) group was investigated in the present study. Of the 33 patients with MIBC, 25 cases were positive PS and 8 were negative PS. Among the 25 positive PS cases, 7 were allocated to receive M-VAC NAC and 18 were allocated to receive CaG NAC according to the results of the prediction systems. Of the 8 negative PS cases, 3 received CaG NAC, 1 received surgery without NAC and 4 received radiation therapy. The total clinical response to NAC was 88.0% (22/25), which was significantly increased compared with the historical controls [56.6% (43/76) P=0.0041]. Overall survival of the positive PS group in the study was significantly increased compared with the negative PS group (P=0.027). In conclusion, the combination of the two prediction systems may increase the treatment efficacy for patients with MIBC by proposing the optimal NAC regimen. In addition, the positive PS group would have a better prognosis compared with the negative PS group. These results suggest that the two prediction systems may lead to the achievement of ‘precision medicine’.

The deubiquitinase USP21 stabilizes MEK2 to promote tumor growth

Deubiquitinases (DUBs) play essential roles in normal cell proliferation and tumor growth. However, the molecular mechanisms of DUBs on hepatocellular carcinoma (HCC) remains largely unknown. In this study, based on analysis of several HCC datasets, we found that the USP21 gene, which encodes a member of the ubiquitin-specific protease family, is highly amplified and overexpressed in HCCs, with the extent of this up-regulation significantly correlating with poor clinical outcomes. Inhibition of USP21 in HCC cell lines decreased cell proliferation, anchorage-independent growth, cell cycle progression, and in vivo tumor growth. Conversely, ectopic expression of USP21 transformed the normal human hepatocyte line HL-7702 and increased the tumorigenicity of the HCC cell line MHCC97L. Mechanistically, USP21 stabilized MEK2 by decreasing its polyubiquitination at Lys48, thereby activating the ERK signaling pathway. Importantly, MEK2 partially mediated the optimal expression of USP21-mediated oncogenic phenotypes. These findings indicate that USP21-mediated deubiquitination and stabilization of MEK2 play a critical role in HCC development.

The clinical significance of COL5A2 in patients with bladder cancer: A retrospective analysis of bladder cancer gene expression data

The present study was aimed to investigate the relationship between the expression of collagen type V alpha 2 chain (COL5A2) and clinical outcomes of patients with bladder cancer.Chi-square test and log-rank-based survival analysis were performed to assess the correlation of COL5A2 expression with clinical characteristics and survivals of patients with bladder cancer using GSE13507. Gene set enrichment analysis was conducted to study the relevant mechanisms.Bladder cancer patients in COL5A2 low expression group were associated with better invasiveness (P < .0001), tumor grade (P = .001), T staging (P < .0001), N staging (P = .002), cancer specific survival (P < .0001), overall survival (P < .0001), and a trend of better M staging (P = .053) than those in COL5A2 high expression group.COL5A2 might affect the progression of bladder cancer through "Coagulation," "Hypoxia," "Apical junction," "Ultraviolet response," "Epithelial mesenchymal transition," "Angiogenesis," "KRAS (KRAS proto-oncogene, GTPase) signaling,""Complement,""IL2-STAT5-signaling," "Inflammatory response," "IL6-JAK-STAT3-signaling," "Myogenesis," "TNF α signaling," "Apoptosis," and "Hedgehog-signaling."Our results demonstrated that COL5A2 was correlated with poor clinical outcomes and survivals of patients with bladder cancer, suggesting that it could be regarded as a biomarker of bladder cancer.

Ubiquitin specific peptidase 21 regulates interleukin-8 expression, stem-cell like property of human renal cell carcinoma

USP family proteins play essential roles in cancer cell proliferation and apoptosis and represent as candidate targets for cancer therapeutics. However, the effects and underlying mechanism of USP21 on renal cell carcinomas (RCC) remain unclear. In the present study, we investigate the effects of USP21 on proliferation, invasion and cancer stem cells (CSCs) property of RCC cell lines. As a result, siRNA-mediated depletion of USP21 inhibits cell proliferation, invasion ability and decreases the CSCs percentage of RCC cell lines. Complementarily, forced expression of USP21 leads to increase of tumorigenic properties. In addition, CSCs properties assessed by sphere formation assays demonstrated that depletion of USP21 impair the self-renewal capability of CSCs. Furthermore, decrease USP21 levels is associated with repression of interleukin 8 (IL-8), a chemokine that regulates CSCs characteristics in RCC. Mechanistically, USP21 binds to the promoter region of IL-8 and mediates transcriptional initiation. These data suggest that USP21/IL-8 could be a pair of the critical molecular targets for the development of therapeutic strategies for RCC.

Identification of Biomarkers Correlated with the TNM Staging and Overall Survival of Patients with Bladder Cancer

Objective: To identify candidate biomarkers correlated with clinical prognosis of patients with bladder cancer (BC). Methods: Weighted gene co-expression network analysis was applied to build a co-expression network to identify hub genes correlated with tumor node metastasis (TNM) staging of BC patients. Functional enrichment analysis was conducted to functionally annotate the hub genes. Protein-protein interaction network analysis of hub genes was performed to identify the interactions among the hub genes. Survival analyses were conducted to characterize the role of hub genes on the survival of BC patients. Gene set enrichment analyses were conducted to find the potential mechanisms involved in the tumor proliferation promoted by hub genes. Results: Based on the results of topological overlap measure based clustering and the inclusion criteria, top 50 hub genes were identified. Hub genes were enriched in cell proliferation associated gene ontology terms (mitotic sister chromatid segregation, mitotic cell cycle and, cell cycle, etc.) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (cell cycle, Oocyte meiosis, etc.). 17 hub genes were found to interact with ≥5 of the hub genes. Survival analysis of hub genes suggested that lower expression of MMP11, COL5A2, CDC25B, TOP2A, CENPF, CDCA3, TK1, TPX2, CDCA8, AEBP1, and FOXM1were associated with better overall survival of BC patients. BC samples with higher expression of hub genes were enriched in gene sets associated with P53 pathway, apical junction, mitotic spindle, G2M checkpoint, and myogenesis, etc. Conclusions: We identified several candidate biomarkers correlated with the TNM staging and overall survival of BC patients. Accordingly, they might be used as potential diagnostic biomarkers and therapeutic targets with clinical utility.

Tobacco-Specific Carcinogens Induce Hypermethylation, DNA Adducts, and DNA Damage in Bladder Cancer

Smoking is a major risk factor for the development of bladder cancer; however, the functional consequences of the carcinogens in tobacco smoke and bladder cancer-associated metabolic alterations remain poorly defined. We assessed the metabolic profiles in bladder cancer smokers and non-smokers and identified the key alterations in their metabolism. LC/MS and bioinformatic analysis were performed to determine the metabolome associated with bladder cancer smokers and were further validated in cell line models. Smokers with bladder cancer were found to have elevated levels of methylated metabolites, polycyclic aromatic hydrocarbons, DNA adducts, and DNA damage. DNA methyltransferase 1 (DNMT1) expression was significantly higher in smokers than non-smokers with bladder cancer. An integromics approach, using multiple patient cohorts, revealed strong associations between smokers and high-grade bladder cancer. In vitro exposure to the tobacco smoke carcinogens, 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone and benzo[a]pyrene (BaP) led to increase in levels of methylated metabolites, DNA adducts, and extensive DNA damage in bladder cancer cells. Cotreatment of bladder cancer cells with these carcinogens and the methylation inhibitor 5-aza-2'-deoxycytidine rewired the methylated metabolites, DNA adducts, and DNA damage. These findings were confirmed through the isotopic-labeled metabolic flux analysis. Screens using smoke-associated metabolites and DNA adducts could provide robust biomarkers and improve individual risk prediction in bladder cancer smokers. Noninvasive predictive biomarkers that can stratify the risk of developing bladder cancer in smokers could aid in early detection and treatment. Cancer Prev Res; 10(10); 588-97. ©2017 AACR.

RASAL2, a RAS GTPase-activating protein, inhibits stemness and epithelial-mesenchymal transition via MAPK/SOX2 pathway in bladder cancer

Muscle-invasive or metastatic bladder cancer (BCa) is associated with a very poor prognosis, and the underlying mechanism remains poorly understood. In this study, we demonstrate RASAL2, a RAS GTPase-activating protein (RAS GAP), acts as a tumor suppressor in BCa. First, RASAL2 was downregulated in BCa specimens and inversely correlated with pathological grades and clinical stages. Furthermore, we observed that RASAL2 could inhibit BCa stemness and epithelial–mesenchymal transition (EMT) based on our gain-of-function and loss-of-function experiments. Mechanistically, we found that mitogen-activated protein kinase/SOX2 signaling had a critical role for maintaining the stemness and mesenchymal properties of RASAL2-deficient BCa cells because inhibition of ERK activity or knockdown of SOX2 could reverse these phenotypes. Also, RASAL2 could inhibit BCa tumorigenesis and distant metastasis in vivo. Moreover, there was an inverse correlation between RASAL2 expression and the stemness/EMT status in subcutaneous xenograft and human BCa specimens. Taken together, our data indicate that RASAL2 is a tumor suppressor in BCa, and modulates cancer stemness and EMT for BCa recurrence and metastasis.