[A risk scoring model based on M2 macrophage-related genes for predicting prognosis of HBV-related hepatocellular carcinoma]

Abstract

Objective: To investigate the prognostic value of M2 macrophage-related genes (MRG) in hepatitis B virus (HBV)- related hepatocellular carcinoma (HCC).

Methods: The transcriptome data of 73 patients with HBV-related HCC were obtained from TCGA database, and the MRG modules were identified by WGCNA. The MRG-based risk scoring model was constructed by LASSO regression analysis and validated using an external dataset. The correlation of the risk score with immune cell infiltration and drug sensitivity of HCC were analyzed with CIBERSORT and R. pRRophetic. The signaling pathways of the differential genes between the high- and low-risk groups were investigated using GSVA and GSEA enrichment analyses, and MRG expressions at the single cell level were validated using R.Seurat. The cell interaction intensity was analyzed by R.Cellchat to identify important cell types related to HCC progression. MRG expression levels were detected by RT-qPCR in THP-1 cells with HCC-conditioned medium-induced M2 polarization and in HBV-positive HCC cells.

Results: A high M2 macrophage infiltration level was significantly correlated with a poor prognosis of HCC, and 5 hub MRG (VTN, GCLC, PARVB, TRIM27, and GMPR) were identified. The overall survival of HCC patients was significantly lower in the high-risk than in the low-risk group. The high- and the low-risk groups showed significant enrichment of M2 macrophages and na?ve B cells, respectively, and were sensitive to BI. 2536 and to AG. 014699, AKT. inhibitor. Ⅷ, AZD. 0530, AZD7762, and BMS. 708163, respectively. The proliferation-related and metabolism-related pathways were enriched in the high-risk group, where monocytes showed the most active cell interactions during HCC progression. VTN was significantly upregulated in HCC cell lines, while GCLC, PARVB, TRIM27, and GMPR were upregulated in M2 THP-1 cells.

Conclusion: The MRG-based risk scoring model can accurately predict the prognosis of HBV-related HCC and reveal the differences in tumor microenvironment to guide precision treatment of the patients.

Keywords: M2 macrophages; hepatitis B virus; hepatocellular carcinoma; prognostic model.

图1

M2巨噬细胞高浸润的HCC患者预后不良 Fig.1 HCC patients with high M2 macrophage infiltration level have poor prognosis. A: Proportion of immune cells in HBV-related HCC patients. B-D: Prognosis of patients stratified by the number of M0 macrophage (B), M1 macrophage (C) and M2 macrophage (D).

图2

风险评分的构建 Fig.2 Construction of prognostic risk scoring model. A-C: Identification of M2 macrophage-related gene (MRG) modules by WGCNA. D-F: Five hub MRGs (VTN, GCLC, PARVB, TRIM27 and GMPR) identified by LASSO regression analysis for constructing the risk scoring model.

图3

风险评分的验证 Fig.3 Validation of the risk scoring model. A-D: Different patterns of survival status and survival time between the high-risk group (A, B) and low-risk group (C, D). E, F: Validation of the MRG prognostic model in the training dataset (E) and testing dataset (F). G, H: ROC curves of the risk score for predicting overall survival (OS) in the training dataset (G) and testing dataset (H). I: Validation of the MRG prognostic model using the external dataset.

图4

列线图的构建 Fig.4 Construction of the nomogram. A: Nomogram of the risk scores and clinical characteristics. B: Calibration curves for evaluating OS predictions at 3 and 5 years. C, D: ROC and DCA curves for determining the accuracy of the nomogram for OS at 1, 3 and 5 years, respectively. E-G: Correlation of the risk scores with clinical stages. H, I: Univariate analysis and multivariate analysis for validating the independent prognostic value of the risk scores.

图5

风险评分的临床预测价值 Fig.5 Clinical predictive value of the risk scoring model. A, B: Infiltrating level of immune cells in the high- and low-risk groups. C: Sensitivity of anti-tumor immunotherapy in the high- and low-risk groups. D: IC₅₀ values of common chemotherapy drugs.

图6

风险评分相关的功能富集分析 Fig.6 Functional enrichment analysis of signaling pathways related to the risk scoring model. A: GSVA in the high- and low-risk groups. B: GSEA in the high- and low-risk groups. C: Molecular interaction networks between the pathways.

图7

MRG与HCC致病基因的关系 Fig.7 Relationship between the 5 hub MRGs and HCC pathogenic genes. A: Relationship between the hub genes and the top 20 HCC pathogenic genes. B: A bubble plot illustrating the Pearson correlation between 5 hub MRGs (VTN, GCLC, PARVB, TRIM27 and GMPR) and the top 20 HCC pathogenic genes.

图8

表达MRG的细胞类型 Fig.8 Expression of MRGs in different cell types. A-C: 5 hub MRGs are mainly expressed in hepatocytes, monocyte, T cells and NK cell. D, E: Cell chat intensity in high- and low-risk groups. F: UMAP projection showing the immune landscape of HCC, colored by cluster. G: Expressions of GCLC, PARVB, GMPR and TRIM27 in HCC immune microenvironment.

图9

MRG在肝癌细胞和巨噬细胞中的表达 Fig.9 Expressions of MRGs in HCC cell lines and macrophages. A: Expression of CD163 and CD206 on THP-1 cells detected by flow cytometry. B-F: Expression of MRGs in PLC/PRF/5 and HCM-THP-1 cells. **P<0.01, ***P<0.001, ****P<0.0001.