IntroductionHypoxia plays an important role in the heterogeneity, relapse, metastasis, and drug resistance of breast cancer. In this study, we explored the hypoxia-related biological signatures in different subtypes of breast cancer and identified the key prognostic factors by bioinformatics methods.MethodsBased on The Cancer Genome Atlas (TCGA) Breast Cancer datasets, we divided the samples into immune-activated/suppressed populations by single-sample gene set enrichment analysis (ssGSEA) and then used hierarchical clustering to further identify hypoxic/non-hypoxic populations from the immune-suppressed samples. A hypoxia related risk model of breast cancer was constructed.ResultsNuclear factor interleukin-3 regulated (NFIL3), serpin family E member 1 (SERPINE1), FOS, biglycan (BGN), epidermal growth factor receptor (EGFR), and sushi-repeat-containing protein, X-linked (SRPX) were identified as key hypoxia-related genes. Margin status, American Joint Committee on Cancer (AJCC) stage, hypoxia status, estrogen receptor/progesterone receptor (ER/PR) status, NFIL3, SERPINE1, EGFR, and risk score were identified as independent prognostic indicators for breast cancer patients. The 3- and 5-year survival curves of the model and immunohistochemical staining on the breast cancer microarray verified the statistical significance and feasibility of our model. Among the different molecular types of breast cancer, ER/PR+ and HER2+ patients might have higher hypoxia-related risk scores. ER/PR-negative samples demonstrated more activated immune-related pathways and better response to most anticancer agents.DiscussionOur study revealed a novel risk model and potential feasible prognostic factors for breast cancer and might provide new perspectives for individual breast cancer treatment.