Abstract
Breast cancer (BC) is the most common tumor in women, and the molecular mechanism underlying its pathogenesis remains unclear. In this study, we aimed to investigate gene modules related to the phenotypes of BC, and identify representative candidate biomarkers for clinical prognosis of BC patients. Using weighted gene co-expression network analysis, we here identified NPY5R as a hub gene in BC. We further found that NPY5R was frequently downregulated in BC tissues compared with adjacent tumor-matched control tissues, due to its aberrant promoter CpG methylation which was confirmed by methylation analysis and treatment with demethylation agent. Higher expression of NPY5R was closely associated with better prognosis for BC patients. Gene set enrichment analysis showed that transcriptome signatures concerning apoptosis and cell cycle were critically enriched in specimens with elevated NPY5R. Ectopic expression of NPY5R significantly curbed breast tumor cell growth, induced cell apoptosis and G2/M arrest. Moreover, NPY5R also promoted the sensitivity of BC cells to doxorubicin. Mechanistically, we found that NPY5R restricted STAT3 signaling pathway activation through interacting with IL6, which may be responsible for the antitumor activity of NPY5R. Collectively, our findings indicate that NPY5R functions as a tumor suppressor but was frequently downregulated in BC.