Abstract Purpose: Breast cancer most often recurs and metastasizes to the distal organs that had their primary tumors surgically excised. Due to heterogeneous nature of the breast cancer, metastasis organotropism has poorly understood. In this study, we assessed the specific cancer-related gene expression changes occurring with metastatic breast cancer recurrence to distant organs comparing with non-metastatic breast cancer. Experimental design: We develop a clinically relevant metastatic syngeneic mouse model using two cell lines 4T1-luc2 that can only metastasize to lung and 4T1.2-luc2 metastasized to bone and lung. After the primary mammary tumor surgically excised from mice, we monitor mice using IVIS imaging for metastatic tumor lesion development in the distant organs followed by MRI scan to detect spinal bone metastases. Animals were sacrificed to collect spines, bones, and lung lesions for molecular analysis. We also validate the gene expression data using the RNA-seq cohort data from the Roswell Park Cancer Institute patient samples with advanced breast cancer bone metastases and matched non-metastatic breast cancer. Results: Our cell lines, RNA-seq data analysis suggested that there are at least 50 genes which are statistically highly expressed in 4T1.2-luc2 compared to 4T1-luc2 cells. Out of which only a few genes are well established for cancer metastasis biology, include ANGPTL77, SERPINE2, TSPAN11, ESM1, LCN2 which are expressed more than 8 fold in 4T1.2-luc2 compared to 4T1-luc2 cells. To our surprise, we found that most of these genes are differentially expressed when metastasized to distant organ. Our data revealed that LCN2 (7 fold) and CD133 (above 20 fold) are overexpressed in the spine and bone compared to the primary or lung met lesions formed by 4T1.2-luc2. Conversely, LCN2 and CD133 genes are downregulated in breast cancer lung metastasis tissues, whereas CD151, EPHA2, and TWIST1 genes are highly overexpressed in lung metastatic lesion compared to primary, bone or spine. Further, the RNA-seq data of patient samples explained that LCN2, CD133 expressions are significantly higher (8 out of 10 patients) in advanced breast cancer bone metastases compared with matching non-metastatic breast cancer. Conclusion: Our data suggested that LCN2, CD133 would be a prognostic marker for breast cancer spinal bone metastasis, whereas CD151, TWIST1, EPHA2 would be a prognostic marker for lung metastasis and organ-specific relapse. Citation Format: Aparna Maiti, Nitai C. Hait. LCN2, CD133 expressions in breast cancer spinal bone metastasis and CD151, TWIST1, EPHA2 in lung metastasis: Prognostic value and organ-specific relapse [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 39.