Statistics and network-based approaches to identify molecular mechanisms that drive the progression of breast cancer

Breast cancer (BC) is one of the most malignant tumors and the leading cause of cancer-related death in women worldwide. So, an in-depth investigation on the molecular mechanisms of BC progression is required for diagnosis, prognosis and therapies. In this study, we identified 127 common differentially expressed genes (cDEGs) between BC and control samples by analyzing five gene expression profiles with NCBI accession numbers GSE139038, GSE62931, GSE45827, GSE42568 and GSE54002, based-on two statistical methods LIMMA and SAM. Then we constructed protein-protein interaction (PPI) network of cDEGs through the STRING database and selected top-ranked 7 cDEGs (BUB1, ASPM, TTK, CCNA2, CENPF, RFC4, and CCNB1) as a set of key genes (KGs) by cytoHubba plugin in Cytoscape. Several BC-causing crucial biological processes, molecular functions, cellular components, and pathways were significantly enriched by the estimated cDEGs including at-least one KGs. The multivariate survival analysis showed that the proposed KGs have a strong prognosis power of BC. Moreover, we detected some transcriptional and post-transcriptional regulators of KGs by their regulatory network analysis. Finally, we suggested KGs-guided three repurposable candidate-drugs (Trametinib, selumetinib, and RDEA119) for BC treatment by using the GSCALite online web tool and validated them through molecular docking analysis, and found their strong binding affinities. Therefore, the findings of this study might be useful resources for BC diagnosis, prognosis and therapies.