Background: Cholangiocarcinoma was a highly malignant liver cancer with poor prognosis, and immune infiltration status was considered an important factor in response to immunotherapy. In this investigation, we tried to locate immune infiltration related genes of cholangiocarcinoma through combination of bulk-sequencing and single-cell sequencing technology.Methods: Single sample gene set enrichment analysis was used to annotate immune infiltration status in datasets of TCGA CHOL, GSE32225, and GSE26566. Differentially expressed genes between high- and low-infiltrated groups in TCGA dataset were yielded and further compressed in other two datasets through backward stepwise regression in R environment. Single-cell sequencing data of GSE138709 was loaded by Seurat software and was used to examined the expression of infiltration-related gene set. Pathway changes in malignant cell populations were analyzed through scTPA web tool.Results: There were 43 genes differentially expressed between high- and low-immune infiltrated patients, and after further compression, PNOC and LAIR2 were significantly correlated with high immune infiltration status in cholangiocarcinoma. Through analysis of single-cell sequencing data, PNOC was mainly expressed by infiltrated B cells in tumor microenvironment, while LAIR2 was expressed by Treg cells and partial GZMB+ CD8 T cells, which were survival related and increased in tumor tissues. High B cell infiltration levels were related to better overall survival. Also, malignant cell populations demonstrated functionally different roles in tumor progression.Conclusion: PNOC and LAIR2 were biomarkers for immune infiltration evaluation in cholangiocarcinoma. PNOC, expressed by B cells, could predict better survival of patients, while LAIR2 was a potential marker for exhaustive T cell populations, correlating with worse survival of patients.

Anti–PD-L1–based combination immunotherapy has become the first-line treatment for unresectable hepatocellular carcinoma (HCC). However, the objective response rate is lower than 40%, highlighting the need to identify mechanisms of tolerance to immune checkpoint inhibitors and accurate biomarkers of response. In this study, we used next-generation sequencing to analyze HCC samples from 10 patients receiving anti–PD-L1 therapy. Activation of the renin–angiotensin system was elevated in nonresponders compared with responders, and angiotensin-converting enzyme 2 (ACE2) expression was significantly downregulated in nonresponders. ACE2 deficiency promoted HCC development and anti–PD-L1 resistance, whereas ACE2 overexpression inhibited HCC progression in immune-competent mice. Mass cytometry by time of flight revealed that ACE2-deficient murine orthotopic tumor tissues featured elevated M2-like tumor-associated macrophages, displayed a CCR5+PD-L1+ immunosuppressive phenotype, and exhibited high VEGFα expression. ACE2 downregulated tumor-intrinsic chemokine (C–C motif) ligand 5 expression by suppressing NF-κB signaling through the ACE2/angiotensin-(1–7)/Mas receptor axis. The lower chemokine (C–C motif) ligand 5 levels led to reduced activation of the JAK–STAT3 pathway and suppressed PD-L1 and VEGFα expression in macrophages, blocking macrophage infiltration and M2-like polarization. Pharmacologic targeting of CCR5 using maraviroc enhanced the tumor-suppressive effect of anti–PD-L1 therapy. Together, these findings suggest that activation of the ACE2 axis overcomes the immunosuppressive microenvironment of HCC and may serve as an immunotherapeutic target and predictive biomarker of response to PD-L1 blockade. Significance: ACE2 regulates the immune landscape of hepatocellular carcinoma by abrogating M2-like macrophage polarization and sensitizes tumors to anti–PD-L1, suggesting that harnessing the ACE2 axis could be a promising strategy to improve immunotherapy efficacy.