M2-like tumor-associated macrophages (M2-like TAMs) have important roles in the progression and therapeutics of cancers. We aimed to detect novel M2-like TAM-related biomarkers in hepatocellular carcinoma (HCC) via integrative analysis of single-cell RNA-seq (scRNA-seq) and bulk RNA-seq data to construct a novel prognostic signature, reveal the “immune landscape”, and screen drugs in HCC.
M2-like TAM-related genes were obtained by overlapping the marker genes of TAM identified from scRNA-seq data and M2 macrophage modular genes identified by weighted gene co-expression network analysis (WGCNA) using bulk RNA-seq data. Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses were carried out to screen prognostic genes from M2-like TAM-related genes, followed by a construction of a prognostic signature, delineation of risk groups, and external validation of the prognostic signature. Analyses of immune cells, immune function, immune evasion scores, and immune-checkpoint genes between high- and low-risk groups were done to further reveal the immune landscape of HCC patients. To screen potential HCC therapeutic agents, analyses of gene–drug correlation and sensitivity to anti-cancer drugs were conducted.
A total of 127 M2-like TAM-related genes were identified by integrative analysis of scRNA-seq and bulk-seq data. PDLIM3, PAM, PDLIM7, FSCN1, DPYSL2, ARID5B, LGALS3, and KLF2 were screened as prognostic genes in HCC by univariate Cox regression and LASSO regression analyses. Then, a prognostic signature was constructed and validated based on those genes for predicting the survival of HCC patients. In terms of drug screening, expression of PAM and LGALS3 was correlated positively with sensitivity to simvastatin and ARRY-162, respectively. Based on risk grouping, we predicted 10 anticancer drugs with high sensitivity in the high-risk group, with epothilone B having the lowest half-maximal inhibitory concentration among all drugs tested.