Abstract
The purpose of the present study was to screen the prognostic targets for breast cancer based on a co-expression modules analysis. The microarray dataset GSE73383 was downloaded from the Gene Expression Omnibus (GEO) database, including 15 breast cancer samples with good prognosis and 9 breast cancer samples with poor prognosis. The differentially expressed genes (DEGs) were identified with the limma package. The Database for Annotation, Visualization and Integrated Discovery was used to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Furthermore, the co‑expression analysis of DEGs was conducted with weighted correlation analysis. The interaction associations were analyzed with the Human Protein Reference Database and BioGRID. The protein‑protein interactions (PPI) network was constructed and visualized by Cytoscape software. A total of 491 DEGs were identified in breast cancer samples with poor prognosis compared with those with good prognosis, and they were enriched in 85 GO terms and 4 KEGG pathways. 368 DEGs were co‑expressed with others, and they were clustered into 10 modules. Module 6 was the most relevant to the clinical features, and 21 genes and 273 interaction pairs were selected out. Abnormal expression levels of required for meiotic nuclear division 5 homolog A (RMND5A) and angiopoietin‑like protein 1 (ANGPTL1) were associated with a poor prognosis. It was indicated that SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily D, member 1, SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily D, member 1, dihydropyrimidinase‑like 2, RMND5A and ANGPTL1 were potential prognostic markers in breast cancer, and the cell cycle may be involved in the regulation of breast cancer.