In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70. © 2021 The Authors