Perturbations of pathway co-expression network identify a core network in metastatic breast cancer

Metastases are the main cause of death in advanced breast cancer (BC) patients. Although chemotherapy and hormone therapy are current treatment strategies, drug resistance is frequent and still not completely understood. In this study, a bioinformatics analysis was performed on BC patients to explore the molecular mechanisms associated with BC metastasis. Microarray gene expression profiles of metastatic and non metastatic BC patients were downloaded from Gene Expression Omnibus (GEO) dataset. Raw data were normalized and merged using the Combat tool. Pathways enriched with differently expressed genes were identified and a pathway co-expression network was generated using Pearson's correlation. We identified from this network, which includes 17 pathways and 128 interactions, the pathways that most influence the network efficiency. Moreover, protein interaction network was investigated to identify hub genes of the pathway network. The prognostic role of the network was evaluated with a survival analysis using an independent dataset. In conclusion, the pathway co-expression network could contribute to understanding the mechanism and development of BC metastases.