BackgroundThe past decade has witnessed unprecedented scientific breakthroughs, including immunotherapy, which has great potential in clinical applications for liver cancer.MethodsPublic data were obtained from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases and analyzed with R software.ResultsThe LASSO and SVM-RFE machine learning algorithms identified 16 differentially expressed genes (DEGs) related to immunotherapy, namely, GNG8, MYH1, CHRNA3, DPEP1, PRSS35, CKMT1B, CNKSR1, C14orf180, POU3F1, SAG, POU2AF1, IGFBPL1, CDCA7, ZNF492, ZDHHC22, and SFRP2. Moreover, a logistic model (CombinedScore) was established based on these DEGs, showing an excellent prediction performance for liver cancer immunotherapy. Patients with a low CombinedScore might respond better to immunotherapy. Gene Set Enrichment Analysis showed that many metabolism pathways were activated in patients with a high CombinedScore, including butanoate metabolism, bile acid metabolism, fatty acid metabolism, glycine serine and threonine metabolism, and propanoate metabolism. Our comprehensive analysis showed that the CombinedScore was negatively correlated with the levels of most tumor-infiltrating immune cells and the activities of key steps of cancer immunity cycles. Continually, the CombinedScore was negatively associated with the expression of most immune checkpoints and immunotherapy response-related pathways. Moreover, patients with a high and a low CombinedScore exhibited diverse genomic features. Furthermore, we found that CDCA7 was significantly correlated with patient survival. Further analysis showed that CDCA7 was positively associated with M0 macrophages and negatively associated with M2 macrophages, suggesting that CDCA7 could influence the progression of liver cancer cells by affecting macrophage polarization. Next, single-cell analysis showed that CDCA7 was mainly expressed in prolif T cells. Immunohistochemical results confirmed that the staining intensity of CDCA7 was prominently increased in the nucleus in primary liver cancer tissues compared to adjacent non-tumor tissues.ConclusionsOur results provide novel insights into the DEGs and factors affecting liver cancer immunotherapy. Meanwhile, CDCA7 was identified as a potential therapeutic target in this patient population.