X. Liu, A. Jamroze, E. Cortes, Y. Ji, K. Zhao, F. Zhao, J. Ho, Y. S. Liu, E. Davicioni, S. Wu, W. Li, F. Y. Feng, J. J. Alumkal, D. E. Spratt, C. J. Sweeney, H. Yu, Q. Hu, C. Zou, D. Zhang, K. Lin, Y. Lu, G. Chatta, K. L. Nastiuk, D. W. Goodrich, K. Rycaj, J. S. Kirk, I. Puzanov, Z. Xia, S. Liu, J. Wang, D. G. Tang
Background: Cancer progression is often accompanied by dedifferentiation and acquisition of stem cell-like properties (stemness). In prostate cancer (PCa), lineage plasticity and therapy resistance remain major clinical challenges, yet a unified quantitative transcriptomic framework connecting stemness, androgen receptor (AR) signaling, castration resistance, and disease progression across the PCa continuum is lacking. Methods: We integrated 87,339 transcriptomic profiles from 33 preclinical and clinical datasets spanning the PCa continuum from normal prostate and treatment-naive primary PCa (Pri-PCa) to PCa treated with neoadjuvant ADT (nADT) and metastatic castration-resistant PCa (mCRPC), with single-cell RNA-seq analyses encompassing 115,197 cells. Cancer stemness was quantified using a transcriptome-derived mRNA-based Stemness Index (mRNAsi; hereafter Stemness), and a 12-gene PCa-Stem signature was developed to capture PCa-specific stemness. Stemness, PCa-Stem, canonical AR activity (c_AR-A), castration-reprogrammed AR activity (cr_AR-A), RB1-loss, PTEN-loss, and MYC activity signatures were quantified across cohorts. Functional validation included MYC inhibition and representative PCa-Stem signature gene depletion in PCa models. Clinical prognostic significance was evaluated in independent patient cohorts. Results: The Stemness score and c_AR-A increased concordantly during early prostate tumorigenesis but diverged with PCa progression: as Gleason grade increased, c_AR-A declined while Stemness continually increased. mCRPC exhibited the highest Stemness and lowest c_AR-A, a pattern recapitulated in Pten/Rb1/Trp53-deficient mouse models. Global Stemness increased progressively across the PCa continuum, was enriched in aggressive PAM50-LumB and PCS1 subtypes, associated with proliferative and lineage plasticity programs, and predicted poor patient survival. The newly derived 12-gene PCa-Stem signature provided a PCa-specific molecular representation of Stemness and tracked disease progression and poor patient survival. Network analyses identified a coordinated mitotic regulatory program linking MYC activity, RB1-loss, cr_AR-A, and the PCa-Stem signature. Spatial and single-cell transcriptomic analyses localized the PCa-Stem program to lineage plasticity-related epithelial cells and demonstrated progressive expansion of PCa-Stem+ epithelial cells during PCa progression. Functional perturbation of representative PCa-Stem signature genes, as well as genetic and pharmacological MYC inhibition, consistently suppressed Stemness-associated phenotypes in diverse PCa models. Conclusions: Cancer Stemness quantitatively captures PCa aggressiveness, lineage plasticity, treatment resistance, disease progression, and poor patient survival. cr_AR-A, RB1 loss, and MYC activation cooperate to reinforce the high-Stemness state and therapy resistance in mCRPC. Collectively, our work establishes a trajectory-integrated transcriptomic framework defining cancer Stemness as a quantifiable molecular and clinical determinant of PCa aggressiveness, lineage plasticity, disease progression, therapy resistance and patient survival.