Major regulators of the multi-step metastatic process are potential therapeutic targets for breast cancer management.

Metastasis is a multi-step process that leads to the dissemination of tumor cells to new sites and, consequently, to multi-organ neoplasia. Although most lethal breast cancer cases are related to metastasis occurrence, little is known about the dysregulation of each step, and clinicians still lack reliable therapeutic targets for metastasis impairment. To fill these gaps, we constructed and analyzed gene regulatory networks for each metastasis step (cell adhesion loss, epithelial-to-mesenchymal transition, and angiogenesis). Through topological analysis, we identified E2F1, EGR1, EZH2, JUN, TP63, and miR-200c-3p as general hub-regulators, FLI1 for cell-adhesion loss specifically, and TRIM28, TCF3, and miR-429 for angiogenesis. Applying the FANMOD algorithm, we identified 60 coherent feed-forward loops regulating metastasis-related genes associated with distant metastasis-free survival prediction. miR-139-5p, miR-200c-3p, miR-454-3p, and miR-1301-3p, among others, were the FFL's mediators. The expression of the regulators and mediators was observed to impact overall survival and to go along with metastasis occurrence. Lastly, we selected 12 key regulators and observed that they are potential therapeutic targets for canonical and candidate antineoplastics and immunomodulatory drugs, like trastuzumab, goserelin, and calcitriol. Our results highlight the relevance of miRNAs in mediating feed-forward loops and regulating the expression of metastasis-related genes. Altogether, our results contribute to understanding the multi-step metastasis complexity and identifying novel therapeutic targets and drugs for breast cancer management.