Breast cancer metastasis is the major reason for deaths from breast cancer. Identification of breast cancer metastasis is of great importance for the management and prediction of cancer progression. However, the key genes and signaling pathways remain unclear in metastatic breast cancers. Our objective is to find the key molecules and signaling pathways by analyzing the RNA-sequence data. The GSE189411 was constructed by the Illumina NovaSeq 6000 (Mus musculus). The KEGG and GO analyses showed the cytokine−cytokine receptor interaction and human papillomavirus infection are the two major processes during the liver metastasis of breast cancer cells. Moreover, we discovered ten key relevant molecules including ITGB2, FCGR3A, CD86, CD80, FOXP3, SYK, CCR5, VCAM1, RAC2, ICAM1. Our study may provide novel insights for the early diagnosis of breast cancer metastasis.