TP53-related signature for predicting prognosis and tumor microenvironment characteristics in bladder cancer: A multi-omics study

Leveraging Single-Cell Multi-Omics to Decode Tumor Microenvironment Diversity and Therapeutic Resistance

Recent developments in single-cell multi-omics technologies have provided the ability to identify diverse cell types and decipher key components of the tumor microenvironment (TME), leading to important advancements toward a much deeper understanding of how tumor microenvironment heterogeneity contributes to cancer progression and therapeutic resistance. These technologies are able to integrate data from molecular genomic, transcriptomic, proteomics, and metabolomics studies of cells at a single-cell resolution scale that give rise to the full cellular and molecular complexity in the TME. Understanding the complex and sometimes reciprocal relationships among cancer cells, CAFs, immune cells, and ECs has led to novel insights into their immense heterogeneity in functions, which can have important consequences on tumor behavior. In-depth studies have uncovered immune evasion mechanisms, including the exhaustion of T cells and metabolic reprogramming in response to hypoxia from cancer cells. Single-cell multi-omics also revealed resistance mechanisms, such as stromal cell-secreted factors and physical barriers in the extracellular matrix. Future studies examining specific metabolic pathways and targeting approaches to reduce the heterogeneity in the TME will likely lead to better outcomes with immunotherapies, drug delivery, etc., for cancer treatments. Future studies will incorporate multi-omics data, spatial relationships in tumor micro-environments, and their translation into personalized cancer therapies. This review emphasizes how single-cell multi-omics can provide insights into the cellular and molecular heterogeneity of the TME, revealing immune evasion mechanisms, metabolic reprogramming, and stromal cell influences. These insights aim to guide the development of personalized and targeted cancer therapies, highlighting the role of TME diversity in shaping tumor behavior and treatment outcomes.

Loss of NUMB drives aggressive bladder cancer via a RHOA/ROCK/YAP signaling axis

Advances in bladder cancer (BCa) treatment have been hampered by the lack of predictive biomarkers and targeted therapies. Here, we demonstrate that loss of the tumor suppressor NUMB promotes aggressive bladder tumorigenesis and worsens disease outcomes. Retrospective cohort studies show that NUMB-loss correlates with poor prognosis in post-cystectomy muscle-invasive BCa patients and increased risk of muscle invasion progression in non-muscle invasive BCa patients. In mouse models, targeted Numb ablation induces spontaneous tumorigenesis and sensitizes the urothelium to carcinogenic insults, accelerating tumor onset and progression. Integrative transcriptomic and functional analyses in mouse and human BCa models reveal that upregulation of YAP transcriptional activity via a RHOA/ROCK-dependent pathway is a hallmark of NUMB-deficient BCa. Pharmacological or genetic inhibition of this molecular pathway selectively inhibits proliferation and invasion of NUMB-deficient BCa cells in 3D-Matrigel organoids. Thus, NUMB-loss could serve as a biomarker for identifying high-risk patients who may benefit from targeted anti-RHOA/ROCK/YAP therapies.

Multi-omics analysis unveils the predictive value of IGF2BP3/SPHK1 signaling in cancer stem cells for prognosis and immunotherapeutic response in muscle-invasive bladder cancer

BACKGROUND

Muscle invasive bladder cancer (MIBC) is a life-threatening malignant tumor characterized by high metastasis rates, poor prognosis, and limited treatment options. Immune checkpoint inhibitors (ICIs) targeting PD-1 and PD-L1 represent an emerging treatment for MIBC immunotherapy. However, the characteristics of patients likely to benefit from immunotherapy remain unclear.

METHODS

We performed single-cell mass cytometry (CyTOF) analysis of 179,483 single cells to characterize potential immunotherapy-related cancer stem cells (CSCs)-like populations in the tumor microenvironment of 38 MIBC tissues. The upregulated expression of IGF2BP3 in CD274 + ALDH + CSC-like cells, which was associated with poor clinical prognosis, was analyzed by bulk RNA-sequencing data from an in-house cohort. The functional role of IGF2BP3 was determined through cell proliferation, colony formation, cell apoptosis and sphere formation assays. The regulation of SPHK1 expression by IGF2BP3 was  investigated using methylated RNA immunoprecipitation sequencing (MeRIP-seq) and bulk RNA-sequencing (bulk RNA-seq). We further utilized single-nucleus RNA sequencing (snRNA-seq) data from 67,988 cells of 25 MIBC tissues and single-cell RNA sequencing (scRNA-seq) data from MIBC patient-derived organoids to characterize the molecular features of bladder cancer cells co-expressing IGF2BP3 and SPHK1. Spatial transcriptomics (ST) and co-detection by indexing (CODEX) analysis were used to describe the spatial distribution and interactions of IGF2BP3 + SPHK1 + bladder cancer cells and immune cells.

RESULTS

A subset of CD274 + ALDH + CSC-like cells was identified, associating with immunosuppression and low survival rates in MIBC patients. IGF2BP3, an m6A reader gene, was found to be upregulated in the CD274 + ALDH + CSC-like cell population and linked to poor clinical prognosis in MIBC. Knockout of IGF2BP3 dramatically promoted cell apoptosis and reduced cell proliferation in T24 cells. By integrating MeRIP-seq and bulk RNA-seq analyses, we identified SPHK1 served as a substrate for IGF2BP3 in an m6A-dependent manner. Further snRNA-seq, scRNA-seq, ST, and CODEX analysis revealed a closer topographical distance between IGF2BP3 + SPHK1 + bladder cancer cells and exhausted CD8 + T cells, providing one explanation for the superior response to immunotherapy in IGF2BP3 + SPHK1 + bladder cancer cells-enriched patients. Finally, an ICI-associated signature was developed based on the enriched genes of IGF2BP3 + SPHK1 + bladder cancer cells, and its potential ability to predict the response to immunotherapy was validated in two independent immunotherapy cohort.

CONCLUSIONS

Our study highlighted the critical involvement of the IGF2BP3/SPHK1 signaling in maintaining the stemness of CSCs and promoting MIBC progression. Additionally, these findings suggested that the IGF2BP3/SPHK1 signaling might serve as a biomarker for prognosis and immunotherapy response in MIBC.

Distinct genomic and immunologic tumor evolution in germline TP53-driven breast cancers

Pathogenic germline TP53 alterations cause Li-Fraumeni Syndrome (LFS), and breast cancer is the most common cancer in LFS females. We performed first of its kind multimodal analysis of LFS breast cancer (LFS-BC) compared to sporadic premenopausal BC. Nearly all LFS-BC underwent biallelic loss of TP53 with no recurrent oncogenic variants except ERBB2 (HER2) amplification. Compared to sporadic BC, in situ and invasive LFS-BC exhibited a high burden of short amplified aneuploid segments (SAAS). Pro-apoptotic p53 target genes BAX and TP53I3 failed to be up-regulated in LFS-BC as was seen in sporadic BC compared to normal breast tissue. LFS-BC had lower CD8+ T-cell infiltration compared to sporadic BC yet higher levels of proliferating cytotoxic T-cells. Within LFS-BC, progression from in situ to invasive BC was marked by an increase in chromosomal instability with a decrease in proliferating cytotoxic T-cells. Our study uncovers critical events in mutant p53-driven tumorigenesis in breast tissue.

Clinical implications of single cell sequencing for bladder cancer

Bladder cancer (BC) is the 10th most common cancer worldwide, with about 0.5 million reported new cases and about 0.2 million deaths per year. In this scoping review, we summarize the current evidence regarding the clinical implications of single-cell sequencing for bladder cancer based on PRISMA guidelines. We searched PubMed, CENTRAL, Embase, and supplemented with manual searches through the Scopus, and Web of Science for published studies until February 2023. We included original studies that used at least one single-cell technology to study bladder cancer. Forty-one publications were included in the review. Twenty-nine studies showed that this technology can identify cell subtypes in the tumor microenvironment that may predict prognosis or response to immune checkpoint inhibition therapy. Two studies were able to diagnose BC by identifying neoplastic cells through single-cell sequencing urine samples. The remaining studies were mainly a preclinical exploration of tumor microenvironment at single cell level. Single-cell sequencing technology can discriminate heterogeneity in bladder tumor cells and determine the key molecular properties that can lead to the discovery of novel perspectives on cancer management. This nascent tool can advance the early diagnosis, prognosis judgment, and targeted therapy of bladder cancer.

A gene signature of cancer-associated fibroblasts predicts prognosis and treatment response in bladder cancer

OBJECTIVE

Due to the pivotal role cancer-associated fibroblasts (CAFs) play in tumor progression, our study aimed to develop a signature of CAFs-related gene (CRG) to predict the survival outcomes and treatment response of bladder cancer (BLCA).

METHODS

The transcriptome data and relevant clinical information about BLCA were collected from publicly available databases, including The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. Weighted gene co-expression network analysis was utilized to uncover CAFs-associated hub genes, and subsequently, a risk model for survival prognosis was constructed using LASSO-Cox regression. The immune microenvironment, immune infiltration, immunotherapy response, and drug sensitivity were explored using ESTIMATE, CIBERSORT, TIDE, and oncoPredict algorithms. To verify the expression of the CRGs, additional analyses were performed using online databases (HPA, CCLE, TIMER, cBioPortal, and TISCH).

RESULTS

Our study developed a CRG signature and constructed a prognostic model. Significant differences in overall survival were observed between the two risk stratifications. The risk score increased with the infiltration of CAFs and tumor staging progression, while closely correlating with immune checkpoint expression and infiltration of CD8 T cells, follicular helper T cells, regulatory T cells, activated dendritic cells, M0 macrophages, M2 macrophages, and resting mast cells. Furthermore, a higher proportion of patients in the low-risk stratification exhibited responsiveness to immunotherapy, and significant variances in sensitivity to multiple chemotherapy medications were observed between the two risk stratifications.

CONCLUSION

The construction of the risk model based on the CRG signature offers new avenues for the prognosis evaluation and development of personalized treatment strategies for BLCA.

Small Cajal Body-Specific RNA12 Promotes Carcinogenesis through Modulating Extracellular Matrix Signaling in Bladder Cancer

Background: Small Cajal body-specific RNAs (scaRNAs) are a specific subset of small nucleolar RNAs (snoRNAs) that have recently emerged as pivotal contributors in diverse physiological and pathological processes. However, their defined roles in carcinogenesis remain largely elusive. This study aims to explore the potential function and mechanism of SCARNA12 in bladder cancer (BLCA) and to provide a theoretical basis for further investigations into the biological functionalities of scaRNAs. Materials and Methods: TCGA, GEO and GTEx data sets were used to analyze the expression of SCARNA12 and its clinicopathological significance in BLCA. Quantitative real-time PCR (qPCR) and in situ hybridization were applied to validate the expression of SCARNA12 in both BLCA cell lines and tissues. RNA sequencing (RNA-seq) combined with bioinformatics analyses were conducted to reveal the changes in gene expression patterns and functional pathways in BLCA patients with different expressions of SCARNA12 and T24 cell lines upon SCARNA12 knockdown. Single-cell mass cytometry (CyTOF) was then used to evaluate the tumor-related cell cluster affected by SCARNA12. Moreover, SCARNA12 was stably knocked down in T24 and UMUC3 cell lines by lentivirus-mediated CRISPR/Cas9 approach. The biological effects of SCARNA12 on the proliferation, clonogenic, migration, invasion, cell apoptosis, cell cycle, and tumor growth were assessed by in vitro MTT, colony formation, wound healing, transwell, flow cytometry assays, and in vivo nude mice xenograft models, respectively. Finally, a chromatin isolation by RNA purification (ChIRP) experiment was further conducted to delineate the potential mechanisms of SCARNA12 in BLCA. Results: The expression of SCARNA12 was significantly up-regulated in both BLCA tissues and cell lines. RNA-seq data elucidated that SCARAN12 may play a potential role in cell adhesion and extracellular matrix (ECM) related signaling pathways. CyTOF results further showed that an ECM-related cell cluster with vimentin+, CD13+, CD44+, and CD47+ was enriched in BLCA patients with high SCARNA12 expression. Additionally, SCARNA12 knockdown significantly inhibited the proliferation, colony formation, migration, and invasion abilities in T24 and UMUC3 cell lines. SCARNA12 knockdown prompted cell arrest in the G0/G1 and G2/M phase and promoted apoptosis in T24 and UMUC3 cell lines. Furthermore, SCARNA12 knockdown could suppress the in vivo tumor growth in nude mice. A ChIRP experiment further suggested that SCARNA12 may combine transcription factors H2AFZ to modulate the transcription program and then affect BLCA progression. Conclusions: Our study is the first to propose aberrant alteration of SCARNA12 and elucidate its potential oncogenic roles in BLCA via the modulation of ECM signaling. The interaction of SCARNA12 with the transcriptional factor H2AFZ emerges as a key contributor to the carcinogenesis and progression of BLCA. These findings suggest SCARNA12 may serve as a diagnostic biomarker and potential therapeutic target for the treatment of BLCA.

Comprehensive identification of onco-exaptation events in bladder cancer cell lines revealed L1PA2-SYT1 as a prognosis-relevant event

Transposable elements (TEs) can provide ectopic promoters to drive the expression of oncogenes in cancer, a mechanism known as onco-exaptation. Onco-exaptation events have been extensively identified in various cancers, with bladder cancer showing a high frequency of onco-exaptation events (77%). However, the effect of most of these events in bladder cancer remains unclear. This study identified 44 onco-exaptation events in 44 bladder cancer cell lines in 137 RNA-seq datasets from six publicly available cohorts, with L1PA2 contributing the most events. L1PA2-SYT1, L1PA2-MET, and L1PA2-XCL1 had the highest frequency not only in cell lines but also in TCGA-BLCA samples. L1PA2-SYT1 showed significant tumor specificity and was found to be activated by CpG island demethylation in its promoter. The upregulation of L1PA2-SYT1 enhances the in vitro invasion of bladder cancer and is an independent risk factor for patient's overall survival, suggesting L1PA2-SYT1 being an important event that promotes the development of bladder cancer.

Identification of tumor immunophenotypes associated with immunotherapy response in bladder cancer

OBJECTIVES

This study aims to reveal immunophenotypes associated with immunotherapy response in bladder cancer, identify the signature genes of immune subtypes, and provide new molecular targets for improving immunotherapy response.

METHODS

Bladder cancer immunophenotypes were characterized in the bulk RNA sequencing dataset GSE32894 and Imvigor210, and gene expression signatures were established to identify the immunophenotypes. Expression of gene signatures were validated in single-cell RNA sequencing dataset GSE145140 and human proteins expression data source. Investigation of Immunotherapy Response was performed in IMvigor210 dataset. Prognosis of tumor immunophenotypes was further analyzed.

RESULTS

Inflamed and immune-excluded immunophenotypes were characterized based on the tumor immune cell scores. Risk score models that were established rely on RNA sequencing profiles and overall survival of bladder cancer cohorts. The inflamed tumors had lower risk scores, and the low-risk tumors were more likely to respond to atezolizumab, receiving complete response/partial response (CR/PR). Patients who responded to atezolizumab had higher SRRM4 and lower NPHS1 and TMEM72 expression than the non-responders. SRRM4 expression was a protective factor for bladder cancer prognosis, while the NPHS1 and TMEM72 showed the opposite pattern.

CONCLUSION

This study provided a novel classification method for tumor immunophenotypes. Bladder cancer immunophenotypes can predict the response to immune checkpoint blockade. The immunophenotypes can be identified by the expression of signature genes.