Comprehensive analysis of macrophage-related genes in prostate cancer by integrated analysis of single-cell and bulk RNA sequencing

Macrophages, as essential components of the tumor immune microenvironment (TIME), could promote growth and invasion in many cancers. However, the role of macrophages in tumor microenvironment (TME) and immunotherapy in PCa is largely unexplored at present. Here, we investigated the roles of macrophage-related genes in molecular stratification, prognosis, TME, and immunotherapeutic response in PCa. Public databases provided single-cell RNA sequencing (scRNA-seq) and bulk RNAseq data. Using the Seurat R package, scRNA-seq data was processed and macrophage clusters were identified automatically and manually. Using the CellChat R package, intercellular communication analysis revealed that tumor-associated macrophages (TAMs) interact with other cells in the PCa TME primarily through MIF - (CD74+CXCR4) and MIF - (CD74+CD44) ligand-receptor pairs. We constructed coexpression networks of macrophages using the WGCNA to identify macrophage-related genes. Using the R package ConsensusClusterPlus, unsupervised hierarchical clustering analysis identified two distinct macrophage-associated subtypes, which have significantly different pathway activation status, TIME, and immunotherapeutic efficacy. Next, an 8-gene macrophage-related risk signature (MRS) was established through the LASSO Cox regression analysis with 10-fold cross-validation, and the performance of the MRS was validated in eight external PCa cohorts. The high-risk group had more active immune-related functions, more infiltrating immune cells, higher HLA and immune checkpoint gene expression, higher immune scores, and lower TIDE scores. Finally, the NCF4 gene has been identified as the hub gene in MRS using the “mgeneSim” function.