Background: Hepatocellular carcinoma (HCC) is among the deadliest cancers worldwide, and advanced HCC is difficult to treat. Identifying specific cell subpopulations in the tumor microenvironment and exploring interactions between the cells and their environment are crucial for understanding the development, prognosis, and treatment of tumors.Methods: In this study, we constructed a tumor ecological landscape of 14 patients with HCC from 43 tumor tissue samples and 14 adjacent control samples. We used bioinformatics analysis to reveal cell subpopulations with potentially specific functions in the tumor microenvironment and to explore the interactions between tumor cells and the tumor microenvironment.Results: Immune cell infiltration was evident in the tumor tissues, and BTG1⁺RGS1⁺ central memory T cells (Tcms) interact with tumor cells through CCL5-SDC4/1 axis. HSPA1B may be associated with remodeling of the tumor ecological niche in HCC. Cancer-associated fibroblasts (CAFs) and macrophages (TAMs) were closely associated with tumor cells. APOC1⁺SPP1⁺ TAM secretes SPP1, which binds to ITGF1 secreted by CAFs to remodel the tumor microenvironment. More interestingly, FAP⁺ CAF interacts with naïve T cells via the CXCL12–CXCR4 axis, which may lead to resistance to immune checkpoint inhibitor therapy.Conclusion: Our study suggests the presence of tumor cells with drug-resistant potential in the HCC microenvironment. Among non-tumor cells, high NDUFA4L2 expression in fibroblasts may promote tumor progression, while high HSPA1B expression in central memory T cells may exert anti-tumor effects. In addition, the CCL5–SDC4/1 interaction between BTG1⁺RGS1⁺ Tcms and tumor cells may promote tumor progression. Focusing on the roles of CAFs and TAMs, which are closely related to tumor cells, in tumors would be beneficial to the progress of systemic therapy research.

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. 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. 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. 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.

Background: The tumor suppressor gene TP53 is frequently mutated or inactivated in bladder cancer (BLCA), which is implicated in the pathogenesis of tumor. However, the cellular mechanisms of TP53 mutations are complicated, yet well-defined, but their clinical prognostic value in the management of BLCA remains controversial. This study aimed to explore the role of TP53 mutation in regulating the tumor microenvironment (TME), elucidate the effects of TP53 activity on BLCA prognosis and immunotherapy response.Methods: A TP53-related signature based on TP53-induced and TP53-repressed genes was used to construct a TP53 activity-related score and classifier. The abundance of different immune cell types was determined using CIBERSORT to estimate immune cell infiltration. Moreover, the heterogeneity of the tumor immune microenvironment between the high and low TP53 score groups was further evaluated using single-cell mass cytometry (CyTOF) and imaging mass cytometry (IMC). Moreover, pathway enrichment analysis was performed to explore the differential biological functions between tumor epithelial cells with high and low TP53 activity scores. Finally, the receptor–ligand interactions between immune cells and tumor epithelial cells harboring distinct TP53 activity were analyzed by single-cell RNA-sequencing.Results: The TP53 activity-related gene signature differentiated well between TP53 functional retention and inactivation in BLCA. BLCA patients with low TP53 scores had worse survival prognosis, more TP53 mutations, higher grade, and stronger lymph node metastasis than those with high TP53 scores. Additionally, CyTOF and IMC analyses revealed that BLCA patients with low TP53 scores exhibited a potent immunosuppressive TME. Consistently, single-cell sequencing results showed that tumor epithelial cells with low TP53 scores were significantly associated with high cell proliferation and stemness abilities and strongly interacted with immunosuppressive receptor–ligand pairs.Conclusion: BLCA patients with low TP53 scores have a worse prognosis and a more immunosuppressive TME. This TP53 activity-related signature can serve as a potential prognostic signature for predicting the immune response, which may facilitate the development of new strategies for immunotherapy in BLCA.

Background: The alternative usage of promoters provides a way to regulate gene expression, has a significant influence on the transcriptome, and contributes to the cellular transformation of cancer. However, the function of alternative promoters (APs) in hepatocellular carcinoma (HCC) has not been systematically studied yet. In addition, the potential mechanism of regulation to the usage of APs remains unclear. DNA methylation, one of the most aberrant epigenetic modifications in cancers, is known to regulate transcriptional activity. Whether DNA methylation regulates the usage of APs needs to be explored. Here, we aim to investigate the effects of DNA methylation on usage of APs in HCC.Methods: Promoter activities were calculated based on RNA-seq data. Functional enrichment analysis was implemented to conduct GO terms. Correlation tests were used to detect the correlation between promoter activity and methylation status. The LASSO regression model was used to generate a diagnostic model. Kaplan–Meier analysis was used to compare the overall survival between high and low methylation groups. RNA-seq and whole-genome bisulfite sequencing (WGBS) in HCC samples were performed to validate the correlation of promoter activity and methylation.Results: We identified 855 APs in total, which could be well used to distinguish cancer from normal samples. The correlation of promoter activity and DNA methylation in APs was observed, and the APs with negative correlation were defined as methylation-regulated APs (mrAPs). Six mrAPs were identified to generate a diagnostic model with good performance (AUC = 0.97). Notably, the majority of mrAPs had CpG sites that could be used to predict clinical outcomes by methylation status. Finally, we verified 85.6% of promoter activity variation and 92.3% of methylation changes in our paired RNA-seq and WGBS samples, respectively. The negative correlation between promoter activity and methylation status was further confirmed in our HCC samples.Conclusion: The aberrant methylation status plays a critical role in the precision usage of APs in HCC, which sheds light on the mechanism of cancer development and provides a new insight into cancer screening and treatment.

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.

Background: Muscle-invasive bladder cancer (MIBC) is a prevalent cancer characterized by molecular and clinical heterogeneity. Assessing the spatial heterogeneity of the MIBC microenvironment is crucial to understand its clinical significance.Methods: In this study, we used imaging mass cytometry (IMC) to assess the spatial heterogeneity of MIBC microenvironment across 185 regions of interest in 40 tissue samples. We focused on three primary parameters: tumor (T), leading-edge (L), and nontumor (N). Cell gating was performed using the Cytobank platform. We calculated the Euclidean distances between cells to determine cellular interactions and performed single-cell RNA sequencing (scRNA-seq) to explore the molecular characteristics and mechanisms underlying specific fibroblast (FB) clusters. scRNA-seq combined with spatial transcriptomics (ST) facilitated the identification of ligand–receptor (L–R) pairs that mediate interactions between specific FB clusters and endothelial cells. Machine learning algorithms were used to construct a prognostic gene signature.Results: The microenvironments in the N, L, and T regions of MIBC exhibited spatial heterogeneity and regional diversity in their components. A distinct FB cluster located in the L region—identified as S3—is strongly associated with poor prognosis. IMC analyses demonstrated a close spatial association between S3 and endothelial cells, with S3-positive tumors exhibiting increased blood vessel density and altered vascular morphology. The expression of vascular endothelial growth factor receptor and active vascular sprouting were significant in S3-positive tumors. scRNA-seq and ST analyses indicated that the genes upregulated in S3 were associated with angiogenesis. NOTCH1–JAG2 signaling pathway was identified as a significant L–R pair specific to S3 and endothelial cell interactions. Further analysis indicated that YAP1 was a potential regulator of S3. Machine learning algorithms and Gene Set Variation Analysis were used to establish an S3-related gene signature that was associated with the poor prognosis of tumors including MIBC, mesothelioma, glioblastoma multiforme, lower-grade glioma, stomach adenocarcinoma, uveal melanoma, kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, and lung squamous cell carcinoma.Conclusions: We assessed the spatial landscape of the MIBC microenvironment and revealed a specific FB cluster with prognostic potential. These findings offer novel insights into the spatial heterogeneity of the MIBC microenvironment and highlight its clinical significance.