AUTHOR=Wang Qiaona , Lin Yunshou , Yu Wenguan , Chen Xiaogang , He Qingqing , Ye Zhiyu 

TITLE=The core role of macrophages in hepatocellular carcinoma: the definition of molecular subtypes and the prognostic risk system

JOURNAL=Frontiers in Pharmacology

VOLUME=14

YEAR=2023

URL=https://www.frontiersin.org/journals/pharmacology/articles/10.3389/fphar.2023.1228052

DOI=10.3389/fphar.2023.1228052

ISSN=1663-9812

ABSTRACT=<p><bold>Background:</bold> In patients with hepatocellular carcinoma (HCC), the tumor microenvironment (TME) is resistant to immunotherapy because of its specificity. It is meaningful to explore the role of macrophage, which is one of the most abundant immune cells in the TME, in cellular communication and its effect on the prognosis and immunotherapy of HCC.</p><p><bold>Methods:</bold> Dimensionality reduction and clustering of the single-cell RNA-seq data from the GSE149614 dataset were carried out to identify the cellular composition of HCC. CellChat was used to analyze the communication between different cells. The specifically highly expressed genes of macrophages were extracted for univariate Cox regression analysis to obtain prognostic genes for HCC cluster analysis, and the risk system of macrophage-specifically highly expressed genes was developed by random forest analysis and multivariate Cox regression analysis. Prognosis, TME infiltration, potential responses to immunotherapy, and antineoplastic drugs were compared among molecular subtypes and between risk groups.</p><p><bold>Results:</bold> We found that HCC included nine identifiable cell types, of which macrophages had the highest communication intensity with each of the other eight cell types. Of the 179 specifically highly expressed genes of macrophage, 56 were significantly correlated with the prognosis of HCC, which classified HCC into three subtypes, which were reproducible and produced different survival outcomes, TME infiltration, and immunotherapy responses among the subtypes. In the integration of four macrophage-specifically highly expressed genes for the development of a risk system, the risk score was significantly involved in higher immune cell infiltration, poor prognosis, immunotherapy response rate, and sensitivity of six drugs.</p><p><bold>Conclusion:</bold> In this study, through single-cell RNA-seq data, we identified nine cell types, among which macrophage had the highest communication intensity with the rest of the cell types. Based on specifically highly expressed genes of macrophage, we successfully divided HCC patients into three clusters with distinct prognosis, TME, and therapeutic response. Additionally, a risk system was constructed, which provided a potential reference index for the prognostic target and preclinical individualized treatment of HCC.</p>