Identification and validation of cancer-associated fibroblast-related subtypes and the prognosis model of biochemical recurrence in prostate cancer based on single-cell and bulk RNA sequencing.

BACKGROUND

Cancer-associated fibroblasts (CAFs) are an essential component of the tumor immune microenvironment that are involved in extracellular matrix (ECM) remodeling. We aim to investigate the characteristics of CAFs in prostate cancer and develop a biochemical recurrence (BCR)-related CAF signature for predicting the prognosis of PCa patients.

METHODS

The bulk RNA-seq and relevant clinical information were obtained from the TCGA and GEO databases, respectively. The infiltration scores of CAFs in prostate cancer patients were calculated using the MCP counter and EPIC algorithms. The single-cell RNA sequencing (scRNA-seq) was downloaded from the GEO database. Subsequently, univariate Cox regression analysis was employed to identify prognostic genes associated with CAFs. We identified two subtypes (C1 and C2) of prostate cancer that were associated with CAFs via non-negative matrix factorization (NMF) clustering. In addition, the BCR-related CAF signatures were constructed using Lasso regression analysis. Finally, a nomogram model was established based on the risk score and clinical characteristics of the patients.

RESULTS

Initially, we found that patients with high CAF infiltration scores had shorter biochemical recurrence-free survival (BCRFS) times. Subsequently, CAFs in four pairs of tumors and paracancerous tissues were identified. We discovered 253 significantly differentially expressed genes, of which 13 had prognostic significance. Using NMF clustering, we divided PCa patients into C1 and C2 subgroups, with the C1 subgroup having a worse prognosis and substantially enriched cell cycle, homologous recombination, and mismatch repair pathways. Furthermore, a BCR-related CAFs signature was established. Multivariate COX regression analysis confirmed that the BCR-related CAFs signature was an independent prognostic factor for BCR in PCa. In addition, the nomogram was based on the clinical characteristics and risk scores of the patient and demonstrated high accuracy and reliability for predicting BCR. Lastly, our findings indicate that the risk score may be a useful tool for predicting PCa patients' sensitivity to immunotherapy and drug treatment.

CONCLUSION

NMF clustering based on CAF-related genes revealed distinct TME immune characteristics between groups. The BCR-related CAF signature accurately predicted prognosis and immunotherapy response in prostate cancer patients, offering a promising new approach to cancer treatment.