Abstract
Prostate cancer (PCa) is a multifocal disease characterized by genomic and phenotypic heterogeneity within a single gland. In this study, Visium spatial transcriptomics (ST) analysis was applied to PCa tissues with different histological structures to infer the molecular events involved in Gleason score (GS) progression. The spots in tissue sections were classified into various groups using Principal Component Analysis (PCA) and Louvain clustering analysis based on transcriptome data. Anotation of the spots according to GS revealed notable similarities between transcriptomic profiles and histologically identifiable structures. The accuracy of macroscopic GS determination was bioinformatically verified through malignancy-related feature analysis, specifically inferred copy number variation (inferCNV), as well as developmental trajectory analyses, such as diffusion pseudotime (DPT) and partition-based graph abstraction (PAGA). Genes related to GS progression were identified from the differentially expressed genes (DEGs) through pairwise comparisons of groups along a GS gradient. The proteins encoded by the representative oncogenes UQCRB and LBH were found to be highly expressed in advanced-stage PCa tissues. Knockdown of their mRNAs significantly suppressed PCa cell proliferation and invasion. These findings were validated using The Cancer Genome Atlas Prostate Adenocarcinoma (TCGA-PRAD) dataset, as well as through histological and cytological experiments. The results presented here establish a foundation for ST-based evaluation of GS progression and provide valuable insights into the GS progression-related genes UQCRB and LBH.