Metastasis is the primary cause of cancer mortality. The primary tumors of colorectal cancer (CRC) often metastasize to the liver. In this study, we have collected 122 samples from 45 CRC patients. Among them, 32 patients have primary tumors, adjacent normal tissues, and matched liver metastases. Thirteen patients have primary tumors without distant metastasis and matched normal tissues. Characterization of these samples was conducted by whole-exome and RNA sequencing and SNP6.0 analysis. Our results revealed no significant difference in genetic alterations including common oncogenic mutations, whole genome mutations and copy number variations between primary and metastatic tumors. We then assembled gene co-expression networks and identified metastasis-correlated gene networks of immune-suppression, epithelial-mesenchymal transition (EMT) and angiogenesis as the key events and potentially synergistic drivers associated with CRC metastasis. Further independent cohort validation using published datasets has verified that these specific gene networks are up regulated throughout the tumor progression. The gene networks of EMT, angiogenesis, immune-suppression and T cell exhaustion are closely correlated with the poor patient outcome and intrinsic anti-PD-1 resistance. These results offer insights of combinational strategy for the treatment of metastatic CRC.