Abstract
Breast cancer (BC) has high morbidity, with significant relapse and mortality rates in women worldwide. Therefore, further exploration of its pathogenesis is of great significance. This study selected therapy genes and possible biomarkers to predict BC using bioinformatic methods. To this end, the study examined 21 healthy breasts along with 457 BC tissues in two Gene Expression Omnibus (GEO) datasets and then identified differentially expressed genes (DEGs). Survival-associated DEGs were screened using the Kaplan-Meier curve. Based on Gene Ontology (GO) annotation, survival-associated DEGs were mostly associated with cell division and cellular response to hormone stimulus. The enriched Kyoto Encyclopedia of Gene and Genome (KEGG) pathway was mostly correlated with cell cycle and tyrosine metabolism. Using overlapped survival-associated DEGs, a survival-associated PPI network was constructed. PPI analysis revealed three hub genes (EZH2, CCNB1, and PPARG) by their degree of connection. These hub genes were confirmed using The Cancer Genome Atlas (TCGA)-BRCA dataset and BC tissue samples. Through Gene Set Enrichment Analysis (GSEA), the molecular mechanism of the potential therapy and prognostic genes were evaluated. Thus, hub genes were shown to be associated with KEGG_CELL_CYCLE and VANTVEER_BREAST_CANCER_POOR_PROGNOSIS gene sets. Finally, based on integrated bioinformatics analysis, this study identified three hub genes as possible prognostic biomarkers and therapeutic targets for BC. The results obtained further understanding of the underground molecular mechanisms related to BC occurrence and prognostic outcomes.