Breast cancer (BC) constitutes a significant peril to global women's health. Contemporary research progressively suggests that mitochondrial dysfunction plays a pivotal role in both the inception and advancement of BC. However, investigations delving into the correlation between mitochondrial-related genes (MRGs) and the prognosis and metastasis of BC are still infrequent. Utilizing data from the TCGA database, we employed the "limma" R package for differential expression analysis. Subsequently, both univariate and multivariate Cox regression analyses were executed, alongside LASSO Cox regression analysis, to pinpoint prognostic MRGs and to further develop the prognostic model. External validation (GSE88770 merged GSE425680) and internal validation were further conducted. Our investigation delved into a broad spectrum of analyses that included functional enrichment, metabolic and immune characteristics, immunotherapy response prediction, intratumor heterogeneity (ITH), mutation, tumor mutational burden (TMB), microsatellite instability (MSI), cellular stemness, single-cell, and drug sensitivity analysis. We validated the protein and mRNA expressions of prognostic MRGs in tissues and cell lines through immunohistochemistry and qRT-PCR. Moreover, leveraging the GSE102484 dataset, we conducted differential gene expression analysis to identify MRGs related to metastasis, subsequently developing metastasis models via 10 distinct machine-learning algorithms and then selecting the best-performing model. The division between training and validation cohorts was set at 70% and 30%, respectively. A prognostic model was constructed by 9 prognostic MRGs, which were DCTPP1, FEZ1, KMO, NME3, CCR7, ISOC2, STAR, COMTD1, and ESR2. Patients within the high-risk group experienced more adverse outcomes than their counterparts in the low-risk group. The ROC curves and constructed nomogram showed that the model exhibited an excellent ability to predict overall survival (OS) for patients and the risk score was identified as an independent prognostic factor. The functional enrichment analysis showed a strong correlation between metabolic progression and MRGs. Additional research revealed that the discrepancies in outcomes between the two risk categories may be attributed to a variety of metabolic and immune characteristics, as well as differences in intratumor heterogeneity (ITH), tumor mutational burden (TMB), and cancer stemness indices. ITH, TIDE, and IPS analyses suggested that patients possessing a low-risk score may exhibit enhanced responsiveness to immunotherapy. Additionally, distant metastasis models were established by PDK4, NRF1, DCAF8, CHPT1, MARS2 and NAMPT. Among these, the XGBoost model showed the best predicting ability. In conclusion, MRGs significantly influence the prognosis and metastasis of BC. The development of dual clinical prediction models offers crucial insights for tailored and precise therapeutic strategies, and paves the way for exploring new avenues in understanding the pathogenesis of BC.