Characterizing ligand-receptor interactions and unveiling the pro-tumorigenic role of CCL16-CCR1 axis in the microenvironment of hepatocellular carcinoma

Background: The heterogeneity of the tumor microenvironment significantly influences the prognosis of hepatocellular carcinoma (HCC) patients, with cell communication through ligand-receptor complexes playing a central role.

Methods: We conducted single-cell transcriptomic analysis on ten HCC tissues to identify ligand-receptor genes involved in malignant HCC cell communication using CellChat. Leveraging RNA-Seq data from the TCGA Liver Cancer (TCGA-LIHC) and Liver Cancer - RIKEN, JP (LIRI-JP) cohorts, we employed Cox regression analysis to screen for prognosis-related genes. Prognostic risk models were constructed through unsupervised clustering and differential gene expression analysis. Subsequently, a co-culture system involving tumor cells and macrophages was established. A series of experiments, including Transwell assays, immunofluorescence staining, immunoprecipitation, flow cytometry, and immunohistochemistry, were conducted to elucidate the mechanism through which HCC cells recruit macrophages via the CCL16-CCR1 axis.

Results: Single-cell analysis unveiled significant interactions between malignant HCC cells and macrophages, identifying 76 related ligand-receptor genes. Patients were classified into three subtypes based on the expression patterns of eight prognosis-related ligand-receptor genes. The subtype with the worst prognosis exhibited reduced infiltration of T cell-related immune cells, downregulation of immune checkpoint genes, and increased M2-like tumor-associated macrophage scores. In vitro experiments confirmed the pivotal role of the CCL16-CCR1 axis in the recruitment and M2 polarization of tumor-associated macrophages. Clinical samples demonstrated a significant association between CCL16 protein expression levels and advanced stage, lymph node metastasis, and distant metastasis. Immunohistochemistry and immunofluorescence staining further confirmed the correlation between CCL16 and CCR1, CD68, and CD206, as well as CD68+CCR1+ macrophage infiltration.

Conclusions: Our study identified molecular subtypes, a prognostic model, and immune microenvironment features based on ligand-receptor interactions in malignant HCC cell communication. Moreover, we revealed the pro-tumorigenic role of HCC cells in recruiting M2-like tumor-associated macrophages through the CCL16-CCR1 axis.