AbstractCurrent prognostic tools cannot clearly distinguish indolent and aggressive prostate cancer (PC). We hypothesized that analyzing individual contributions of epithelial and stromal components in localized PC (LPC) could improve risk stratification, as stromal subtypes may have been overlooked due to the emphasis on malignant epithelial cells. Hence, we derived molecular subtypes of PC using gene expression analysis of LPC samples from prostatectomy patients (cohort 1, n = 127) and validated these subtypes in two independent prostatectomy cohorts (cohort 2, n = 406, cohort 3, n = 126). Stroma and epithelium-specific signatures were established from laser-capture microdissection data and non-negative matrix factorization was used to identify subtypes based on these signatures. Subtypes were functionally characterized by gene set and cell type enrichment analyses, and survival analysis was conducted. Three epithelial (E1-E3) and three stromal (S1-S3) PC subtypes were identified. While subtyping based on epithelial signatures showed inconsistent associations to biochemical recurrence (BCR), subtyping by stromal signatures was significantly associated with BCR in all three cohorts, with subtype S3 indicating high BCR risk. Subtype S3 exhibited distinct features, including significantly decreased cell-polarity and myogenesis, significantly increased infiltration of M2-polarized macrophages and CD8 + T-cells compared to subtype S1. For patients clinically classified as CAPRA-S intermediate risk, S3 improved prediction of BCR. This study demonstrates the potential of stromal signatures in identification of clinically relevant PC subtypes, and further indicated that stromal characterization may enhance risk stratification in LPC and may be particularly promising in cases with high prognostic ambiguity based on clinical parameters.