Identification of two immune-related risk score signatures through integrated analysis of multi-omics data in hepatocellular carcinoma

•Two risk scores (ICS and ICRGS) were constructed through integrated analysis of bulk and scRNA sequencing data in HCC.•The functional roles of the DEGs between the prognostic immune cells were clarified by enrichment analyses.•The interaction among the prognostic immune cells was preliminarily investigated by cell–cell communication analysis.•Potential therapeutic drugs targeting prognostic DEGs were screened by CellMiner. Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related death worldwide. Immunotherapy has become a major treatment for advanced HCC, but the therapeutic effects remain unsatisfactory. In this study, we constructed an immune cell risk score (ICS) and an immune cell-related gene risk score (ICRGS) for the prognosis prediction of HCC through integrated analysis of bulk and single-cell RNA (scRNA) sequencing data. These two risk score signatures both showed good predictive values in the training and validation cohorts. The potential interactions among these prognostic immune cell types were elucidated by cell–cell communication analysis. The results of enrichment analysis and gene set enrichment analysis (GSEA) of the prognostic genes showed that metabolic-related processes were involved in the immune response of HCC. Furthermore, the results of correlation analyses further confirmed the hub genes that were strongly correlated with immune cells. Finally, potential therapeutic drugs targeting these hub genes were screened by CellMiner based on NCI-60 cell line set. Taken together, two useful models for the prognosis prediction of HCC patients were constructed in this study. The functional differences between the two groups of HCC patients separated by ICS or ICRGS provide fundamental knowledge for finding synergistic therapeutic targets for HCC immunotherapy.