Abstract
Current evidence indicate that cancer-associated fibroblasts (CAFs) play an important role in prostate cancer (PCa) development and progression. In this study, we identified CAF-related molecular subtypes and prognostic index for PCa patients undergoing radical prostatectomy through integrating single-cell and bulk RNA sequencing data. We completed analyses using software R 3.6.3 and its suitable packages. Through single-cell and bulk RNA sequencing analysis, NDRG2, TSPAN1, PTN, APOE, OR51E2, P4HB, STEAP1 and ABCC4 were used to construct molecular subtypes and CAF-related gene prognostic index (CRGPI). These genes could clearly divide the PCa patients into two subtypes in TCGA database and the BCR risk of subtype 1 was 13.27 times higher than that of subtype 2 with statistical significance. Similar results were observed in MSKCC2010 and GSE46602 cohorts. In addtion, the molucular subtypes were the independent risk factor of PCa patients. We orchestrated CRGPI based on the above genes and divided 430 PCa patients in TCGA database into high- and low- risk groups according to the median value of this score. We found that high-risk group had significant higher risk of BCR than low-risk group (HR: 5.45). For functional analysis, protein secretion was highly enriched in subtype 2 while snare interactions in vesicular transport was highly enriched in subtype 1. In terms of tumor heterogeneity and stemness, subtype 1 showd higher levels of TMB than subtype 2. In addition, subtype 1 had significant higher activated dendritic cell score than subtype 2. Based on eight CAF-related genes, we developed two prognostic subtypes and constructed a gene prognostic index, which could predict the prognosis of PCa patients very well.