Metastasis is the main cause of mortality in colorectal cancer (CRC) patients. Exploring the mechanisms of metastasis is of great importance in both clinical and fundamental CRC research. CRC is a highly heterogeneous disease with variable therapeutic outcomes of treatment. In this study, we applied spatial transcriptomics (ST) to generate a tissue-wide transcriptome from two primary colorectal cancer tissues and their matched liver metastatic tissues. Spatial RNA information showed intratumoral heterogeneity (ITH) of both primary and metastatic tissues. The comparison of gene expressions across tissues revealed an apparent enrichment of cancer stem cells (CSCs) in metastatic tissues and identified FOXD1 as a novel metastatic CSC marker. Trajectory and pseudo-time analyses revealed distinct evolutionary trajectories and a dedifferentiation-differentiation process during metastasis. CellphoneDB analysis suggested a dominant interaction of CD74-MIF with tumor cells in metastatic tissues. Further analysis confirmed FOXD1 as a maker of CSCs and the predictor of patient survival, especially in metastatic diseases. Our study found ITH of primary and metastatic tissues and provides novel insights into the cellular mechanisms underlying liver metastasis of CRC and foundations for therapeutic strategies for CRC metastasis.