Tumor metastasis causes over 90% of cancer related death and no currently available therapies target it. However, there is limited understanding regarding the epigenetic regulation of genes during this complex process. Here by integrating single-cell ATAC-seq (scATAC-seq), single-cell RNA-seq (scRNA-seq), microarray, bulk RNA-seq, immunohistochemistry (IHC) staining, as well as proteomics datasets from paired primary and liver metastatic colorectal cancer (CRC) patient-derived xenograft (PDX) model and patients, we discovered that liver metastatic CRC cells lose their colon-specific chromatin accessible sites yet gain liver-specific ones. Importantly, we observed elevated accessibility of HNF4A, a liver-specific transcription factor, in liver metastatic CRC cells. Subsequently, we performed clustering analysis of liver metastatic CRC cells together with cells involved in liver development, revealing significant heterogeneity among the liver metastatic CRC cells. Over 50% of the liver metastatic CRC cells exhibited characteristics similar to those of erythroid progenitors and hepatocytes, showing increased expression of genes involved in oxidative phosphorylation and glycolysis. Moreover, our discovery further revealed that the MHC and IFN response genes in these cells exhibit moderate epigenetic activity, which is significantly associated with the low objective response rates in checkpoint blockade immunotherapy. Our findings uncovered the critical roles of HNF4A and the cell populations within liver metastatic CRC cells might serve as crucial therapeutic targets for addressing liver metastasis and improving the immunotherapy response in patients with CRC.Copyright © 2023. Published by Elsevier Ltd.