Analysis of tumor environmental response and oncogenic pathway activation identifies distinct basal and luminal features in HER2-related breast tumor subtypes

We report here the development of a series of gene expression signatures to measure the influences of microenvironmental stresses. The pathway activities of hypoxia, lactic acidosis, acidosis and glucose deprivation were investigated in a collection of 1,143 breast tumors, which have been separated into 17 breast tumor subgroups defined by their distinct patterns of oncogenic pathways. A validation dataset comprised of 547 breast tumors was also used to confirm the major findings, and representative breast cancer cell lines were utilized to validate in silico results and mechanistic studies.Through the integrative pathway analysis of microenvironmental stresses and oncogenic events in breast tumors, we identified many known and novel correlations between these two sources of tumor heterogeneity. Focusing on differences between two human epidermal growth factor receptor 2 (HER2)-related subgroups, previously identified based on patterns of oncogenic pathway activity, we determined that these subgroups differ with regards to tumor microenvironmental signatures, including hypoxia. We further demonstrate that each of these subgroups have features consistent with basal and luminal breast tumors including patterns of oncogenic signaling pathways, expression of subtype specific genes, and cellular mechanisms that regulate the hypoxia response. Importantly, we also demonstrate that the correlated pattern of hypoxia-related gene expression and basal-associated gene expression are consistent across HER2-related tumors whether we analyze the tumors as a function of our pathway-based classification scheme, using the intrinsic gene list (ERBB2+), or based on HER2 IHC status. Our results demonstrate a cell lineage-specific phenomenon in which basal-like tumors, HER2-related tumors with high hypoxia, as well as normal basal epithelial cells express increased mRNA levels of HIF-1α compared to luminal types and silencing of HIF-1α results in decreased expression of hypoxia-induced g