- Research article
- Open Access
Analysis of tumor environmental response and oncogenic pathway activation identifies distinct basal and luminal features in HER2-related breast tumor subtypes
Breast Cancer Research volume 13, Article number: R62 (2011)
Breast cancer heterogeneity occurs as a consequence of the dysregulation of numerous oncogenic pathways as well as many non-genetic factors, including tumor microenvironmental stresses such as hypoxia, lactic acidosis, and glucose deprivation. Although the importance of these non-genetic factors is well recognized, it is not clear how to integrate these factors within the genetic framework of cancer as the next logical step in understanding tumor heterogeneity.
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 genes.
This study demonstrates differences in microenvironmental conditions in HER2-related subgroups defined by distinct oncogenic pathway activities, and provides a mechanistic explanation for differences in the observed hypoxia response between these subgroups. Collectively, these data demonstrate the potential of a pathway-based classification strategy as a framework to integrate genetic and non-genetic factors to investigate the basis of tumor heterogeneity.
Breast cancer is a collection of distinct diseases characterized by differences in oncogenic mechanisms and clinical characteristics including prognosis and response to therapeutic regimens. Clinically, tumors are classified on the basis of tumor size, visual characteristics, and a limited number of histochemical markers including estrogen receptor, progesterone receptor and HER2 receptor status. While the strongest and most consistent division in breast cancer is between the basal and luminal types and systematic differences strongly suggest a different cell of origin for these two dominant categories [1–5] additional molecular or mechanistic heterogeneity can be identified within classes.
Amplification of HER2 (ERBB2) occurs in approximately 15 to 25% of human breast tumors and defines a clinically unique subgroup of breast tumors . Despite the introduction of Herceptin (trastuzumab), a monoclonal antibody that therapeutically targets HER2, the prognosis for these patients remains poor since a substantial number of patients either fail to respond to this therapy or develop resistance over time [7–10]. Therefore, it is clear that additional heterogeneity exists within this class of tumors. A number of recent studies have demonstrated the genomic complexity of HER2+ tumor heterogeneity, reporting that these tumors have varied and complex patterns of copy number alterations, global gene expression, and DNA methylation [11, 12]. Moreover, molecular subtypes based on analysis of gene expression patterns  as well as more complex analyses of patterns of oncogenic and tumor suppressor pathway activity  have identified further subtypes that provide additional insight into the molecular and clinical heterogeneity of HER2+ tumors beyond the original ERBB2+ classification [1–5]. These studies also suggest that further investigation of the basis of HER2+ tumor heterogeneity, including the interaction of genetic and non-genetic factors, is necessary to fully elucidate the molecular mechanisms driving HER2-mediated oncogenesis.
Inherited germline mutations as well as acquired somatic mutations, unique to each individual patient, provide a source of variation at the genetic level. However, additional variation in the environmental conditions further influence tumor phenotypes and select for tumors adapted for growth within these microenvironments. The most visible of these environmental factors relate to common changes in the physical and chemical alterations of cancers, including reduced oxygen tension (hypoxia), high lactate, extracellular acidosis (lactic acidosis), and glucose depletion [15–18]. These changes may result from poor tissue perfusion, abnormal tumor vasculature, and/or genetic dysregulation of metabolism in the cancer cell.
Hypoxia (or low pO2) is recognized as a risk factor for poor clinical outcomes and increased metastasis, partly due to the gene expression response to hypoxia . The cellular response to hypoxia is triggered by the stabilization of hypoxia regulators, HIF-1α or EPAS1 (HIF-2α) proteins, under limited oxygen. However, the HIF protein stabilization and hypoxia response are not limited to low pO2 and may also result from a wide range of genetic alterations and signaling malfunctions. These conditions represent "pseudo-hypoxia" and include the loss of VHL , TP53 , or PTEN  or activation of PI3-kinase/Akt  and HER2  pathways. In addition, the HIF genes can also be expressed at different levels among different cells and contribute to the varying degrees of hypoxia response [24, 25]. Hypoxia also favors glycolysis and leads to lactic acidosis with the accumulation of lactate and acidity in solid tumors. Avid glucose uptake and glycolysis can lead to glucose deprivation, another stress found in solid tumors [26, 27]. Glucose deprivation triggers activation of AMPK and LKB-1, which in turn activates TSC1/TSC2 and inhibits the central energy sensor mTOR . In addition, long term exposure to glucose deprivation has been shown to contribute to KRAS mutation in human cancers .
Given the importance of these microenvironmental stresses, many methods have been developed to measure their occurrence in solid tumors. These methods include direct measurements of oxygen tension (EF5, Eppendorf polarographic probe), acidity (pH probes), and lactate/glucose levels (bioluminescence technique) [30–33]. These approaches allow us to obtain a calibrated map of the respective metabolites' distributions within the tumor tissues that indicates relative abundance and spatial distribution. The level of each of these stresses in viable tumor areas shows large variations between tumors and can be related to tumor behavior and different clinical risks [30, 34–38]. Therefore, an ability to assess these physiological parameters prior to treatment has the potential to aid in the development of novel therapeutic strategies for individual patients.
Drawbacks to these in vivo tumor measurements are that they are frequently invasive, technically challenging, or require tumors to be snap-frozen in a sophisticated laboratory setting; they are not amenable to routine and optimal clinical care. We have been working towards developing a conceptual framework to evaluate and incorporate the influences of these stresses to better understand the basis of tumor heterogeneity. The influence of tumor environment can be manifested by the expression of "endogenous" response genes; for example, tissue hypoxia can be inferred by elevated levels of CA9 . When the cellular response genes are captured by microarrays in the form of "gene signatures", these in vitro derived gene expression signatures are capable of quantitatively assessing and classifying tumors based on their environmental profile analogous to oncogenic pathway activation and other phenotypes [24, 40–45].
We recently reported the development of a tumor classification strategy based on patterns of oncogene and tumor suppressor activity as measured by gene expression signatures of each oncogenic pathway . This study identified 17 breast tumor subgroups that are associated with the intrinsic subtypes of breast cancers and have shown an ability to further delineate potential mechanisms of disease by identifying novel patterns of oncogenic pathway activity and copy number variation between subgroups [1, 46]. In the current study, we apply a similar framework towards understanding the nature of environmental and metabolic stresses in the context of defined lineage and oncogene specific tumor characteristics. Our findings indicate that a strong hypoxic response is the most characteristic of basal type tumors, which extends to a subset of HER2 related cancers that express characteristics consistent with derivation from the basal lineage. These differences in the hypoxia responses are likely due to the cell-type specific expression levels of HIF mRNAs which are associated with mammary differentiation programs.
Materials and methods
Gene expression signatures and tumor samples
The methods and a description of the training data used to develop the oncogenic pathway gene expression signatures have been published  and a detailed description of the methods, including the training data, for the development of the microenvironment signatures as well as for signatures to measure basal-luminal and subgroup 7/10 characteristics are provided in the Additional file 1, Supplemental methods section. Signature parameters, genes, and regression weights for the microenvironment signatures are reported in Additional file 2 Table S1. Briefly, a signature represents a group of genes that collectively demonstrate a consistent pattern of expression in relation to a given phenotype. Each signature is derived from the first principal component or the factor corresponding to the largest singular value as determined by singular value decomposition. A binary probability regression model is then estimated using Bayesian methods based on the vectors representing the two phenotypic states of the training data. The gene selection, identification, and regression model is based solely on the training data and maintains statistical independence from the validation dataset. This enables evaluation of predictive probabilities of each of the two phenotypic states in the training data for each sample in the validation datasets. Two previously described meta-datasets containing expression data from 1,143 and 547 human breast tumor samples respectively, were collected from 10 and 2 independent datasets for which Affymetrix HG-U133 CEL files were publicly available . A summary of publicly available samples utilized in the study is provided in Additional file 3, Table S2. Samples were first RMA or MAS5 normalized using Affymetrix Expression Console ver1.0. MAS5 data were log2 transformed and both RMA and log2MAS5 datasets were filtered to include only those probes on the Affymetrix (Santa Clara, CA, USA) HG-U133A array. Samples were then normalized by Bayesian Factor Regression Analysis (BFRM) to remove technical variation as previously reported . All studies using human data were performed in compliance with the Helsinki Declaration.
AU565 and HCC202 cells were maintained in RPMI with 15% FBS and 1% penicillin/streptomycin while MCF7, BT474, CAMA1, BT20, MDAMB231, and MDAMB157 were cultured in DMEM with 10% FBS and 1% penicillin/streptomycin. Primary luminal and basal cells are obtained from human breast organoids and isolated as described . For hypoxia, cells were cultured under 1% O2 for 4 or 24 hours.
RNA interference and overexpression
The cells were plated in six-well plates at a density of 3 × 105 cells per well. For RNAi, cells were transfected with control non-targeting or siRNA against HIF-1α or HIF-2α with lipofectamine (Invitrogen, Carlsbad, CA, USA). RNA was extracted 48 hours after transfection and the levels of indicated transcripts were examined by real-time RT-PCR.
Total RNAs from cells, including cancer cell lines and primary breast cells, under normal culture condition or 48 h of hypoxia (1% O2) were isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and then reverse-transcribed to cDNAs with SuperScript II reverse transcription kit (Invitrogen, Carlsbad, CA, USA) and used for real-time PCR with Power SYBRGreen PCR Mix (Applied Biosystem, Foster City, CA, USA) and primers for HIF1A (TGCTCATCAGTTGCCACTTC, CAGAAGTTTCCTCACACGCA), EPAS1 (CCATGTCTCCACCTTCAGA, GCTTCAGCTTCAGCTTGTTG), ACTB (β-actin) (CTCTTCCAGCCTTCCTTCCT, AGCACTGTGTTGGCGTACAG), CA9 (TCCTCAAGAACCCCAGAATAA, CCTCCATAGCGCCAATGACT), GLUT1 (AACTCTTCAGCCAGGGTCCAC, CACAGTGAAGATGATGAAGAC), VEGFA (TGCTCTACCTCCACCATGCCAA, TGATGATTCTGCCCTCCTCCTTC). The relative levels of each gene from real-time PCR were normalized with the level of β-actin.
Analysis of gene expression levels
To compare relative levels of expression for individual genes or the average of several related genes, RMA normalized Affymetrix U133A expression data from the primary or the validation datasets were compared. The Affymetrix probes used for each gene or condition are reported in the Supplemental methods. For the analysis of the average basal-like or hypoxia-related gene expression, probes associated with each condition, reported in the Supplemental methods, were averaged across each sample. An unpaired t-test or Mann-Whitney U-test, as reported in the text, was used to compare between groups and a linear regression was used to compare expression values on a continuous variable scale.
Analysis of co-regulation of oncogenic signaling and microenvironment pathways in breast cancer
To investigate the influence of genetic and non-genetic factors in human breast cancer, we first constructed gene expression signatures to define the influences of a series of microenvironmental stresses and then analyzed them in the context of previously defined oncogenic pathways. We developed a series of novel gene expression signatures capable of quantitatively assessing the influence of hypoxia, lactic acidosis, acidosis, and glucose depletion (Glu(-)) using previously published gene expression studies of breast cells [42, 44, 48] (Additional file 4 Figure S1 and Additional file 3, Table S1). The pathway activities of these stress signatures were integrated with a series of 18 previously reported oncogenic pathway signatures using a Bayesian binary regression strategy to quantitatively assess the stress pathway activities in a collection of 1,143 human breast tumors that were collected from 10 independent datasets and normalized by Bayesian Factor Regression Modeling (BFRM) to remove technical variation . The predicted hypoxia and lactic acidosis pathway activities were further validated by their respective significant correlations with the average expression levels of previously defined hypoxia-inducible genes  and the expression of TXNIP, which is induced by lactic acidosis and is the putative marker of the lactic acidosis response  (Additional file 4, Figure S1). Concurrently, the predicted probabilities of 18 oncogenic signaling pathways were determined in this dataset  (Figure 1) to enable the integrative analysis.
Two-way hierarchical clustering (Figure 1A) and a Pearson correlation (Figure 1B) was calculated to define the statistical correlation between co-activated pathways These analyses identify three separate clusters of pathways with a clear clustering among the hypoxia, Glu(-) and acidosis pathways together with the EGFR, TGFβ, HER2, TNFα and STAT3 pathways (Figure 1A). A significant relationship was also identified between the ER, PR, and p53 pathways which were found to be grouped with lactic acidosis pathway activity in a second cluster (Figure 1). Our previous study on the patterns of oncogenic pathway activity in human breast tumors  demonstrated clear correlations between the interferon pathways as well as between MYC and RAS, which are consistent with the findings of many other studies [49–51]. As illustrated in Figure 1A, these correlations persist with the inclusion of environmental conditions. Moreover, we are able to confirm the previously identified positive and negative relationships between lactic acidosis with p53 and PI3K/AKT pathways, respectively . Finally, the positive correlation and clustering of the hypoxia pathway with EGFR, TGFβ and STAT3 was consistent with their reported reciprocal positive regulation [52, 53]. The re-discovery of these connections among different pathways in human tumors provides a measure of validation for the expression signatures and the predictive approach to integrate genetic and non-genetic sources of tumor phenotypes.
Analysis of microenvironment conditions in pathway-defined subgroups of breast cancer identifies differences between HER2 related subgroups
We recently reported the identification of 17 breast cancer subgroups based on patterns of 18 different oncogene and tumor suppressor pathway activities . These subgroups, while not simply a refinement of the molecular subtypes defined by the intrinsic gene list (that is, basal, luminal A, luminal B, ERBB2+, normal-like), have been shown to be significantly associated with these intrinsic subtypes to further delineate biological and clinical heterogeneity within these groups. Therefore, we next investigated microenvironmental stress pathway activities in the context of these 17 subgroups (Additional file 4, Figure S2). Interestingly, in subgroups 7 and 10, which are characterized by high levels of HER2 pathway activity (Additional file 4, Figure S2) and have been shown to be enriched for ERBB2+ tumors , the pathway activities of hypoxia/EGFR/TNFα/TGFβ/STAT3 showed consistently higher levels in subgroup 7 than subgroup 10 tumors (Figures 2A, 2B). As illustrated in Figure 2Aand quantitatively assessed in Figure 2B, clear differences in patterns of predicted microenvironmental pathway activities, specifically hypoxia response, glucose depletion, and acidosis, as well as a number of hypoxia-associated oncogenic signaling pathways, including EGFR/TNFα/TGFβ/STAT3, are characteristically different between subgroups 7 and 10. To quantitatively assess these differences, an unpaired t-test was used to compare patterns of predicted pathway activity in each subgroup (Figure 2B and Additional file 4, Figure S3). We found that subgroup 7 has significantly (P <0.0001, unpaired t-test) higher levels of predicted hypoxia response and many other pathways that were shown to be clustered with hypoxia, including EGFR, TGFβ, TNFα, STAT3, Glu(-) and acidosis (Figure 2A and Additional file 4 Figure S3). In contrast, tumors assigned to subgroup 10 had significantly higher levels of ER, PR, IFNα, IFNγand SRC pathway activity (Additional file 4, Figure S3). In order to validate the differences in patterns of microenvironmental pathway activity, a validation dataset comprised of 547 independent breast tumor samples was assessed for patterns of oncogenic pathway activity and assigned to subgroups . Patterns of microenvironmental pathway activity were then assessed (Figure 2C) and examined quantitatively (Figure 2D) demonstrating patterns of microenvironmental pathway activity in subgroups 7 and 10 consistent with the discovery dataset of 1,143 samples. Specifically, this independent, validation dataset also demonstrates a significant difference (unpaired t-test) in patterns of hypoxia response, glucose depletion, and acidosis (Figure 2D).
Subgroup 7 tumors exhibit a strong hypoxia response
Since tumors in subgroup 7 were characterized by a higher level of hypoxia response (Figure 2 and Additional file 4, Figure S3) when compared to subgroup 10 tumors in both a primary and validation dataset, we next examined whether these tumors showed differences in expression of known hypoxia-induced genes. As shown in Figure 3A and 3B (and in higher resolution in Additional file 4, Figures S4A and S4B), subgroup 7 tumors demonstrated higher expression levels of a panel of previously-defined common hypoxia-induced genes  (list of 89 Affymetrix probes provided in Supplemental methods) in both the primary (Figure 3A) and validation (Figure 3B) datasets. Consistent with these observations, significantly higher mRNA levels of many hypoxia-induced genes including VEGFA (P = 0.0004, Mann-Whitney U-test), EGLN3 (P <0.0001), GLUT1 (P = 0.0121), and DEC1 (P <0.0001) were observed in subgro up 7 tumors (Additional file 4, Figure S5). Collectively these data suggest that tumors in subgroups 7 and 10, while both sharing high HER2 pathway activity, differ not only in patterns of oncogenic pathway signaling but are also characterized by differences in the degrees of hypoxia response.
In order to understand the basis for the exaggerated hypoxia response observed in subgroup 7 tumors, we next analyzed gene expression levels of hypoxia regulators in these tumors. It was determined that subgroup 7 tumors demonstrated a greater than two-fold higher level of mRNA expression (P <0.0001, Mann-Whitney U-test) for both HIF1A (HIF-1α) (Figures 3C, D) and EPAS1 (HIF-2α) (Figures 3F, G) as compared to subgroup 10 tumors in both the primary (Figures 3C, F) and validation datasets (Figures 3D, G). While HIF-1α and HIF-2α are usually regulated at the protein level via oxygen-dependent degradation, several studies have reported a similar mRNA induction associated with a stronger hypoxia response [24, 25]. To further validate the hypoxia response in subgroup 7 tumors, a number of additional known hypoxia regulators were identified from BioCarta, and it was determined that subgroup 7 tumors exhibited significantly higher levels of expression of ARNT (P = 0.0078), EP300 (P = 0.009), JUN (P <0.0001), and LDHA (P = 0.0092) (Additional file 4, Figure S6). Because several studies have also linked the level of HER2 amplification to the hypoxia pathway [23, 54], we compared the mRNA levels of HER2 among these two HER2 related subgroups and found no significant difference in either the primary or validation datasets (Additional file 4, Figure S6).
As we noted previously , a clear advantage of utilizing the pathway-based classification scheme is an ability to assign new samples including cancer cell lines to corresponding subgroups based on the similarity of pathway composition. Therefore, we took advantage of this ability to identify in vitro cell culture model systems to test differences between subgroups 7 and 10. Using a previously analyzed dataset of 50 breast cancer cell lines , we identified two HER2+ breast cancer cell lines representative of subgroup 7 (AU565) and subgroup 10 (HCC202), each with a probability of subgroup membership greater than 0.90. Consistent with tumor data, we determined that the subgroup 7 cell line (AU565) has significantly higher mRNA levels of HIF1A (HIF-1α) (Figure 3E) and EPAS1 (HIF-2α) (Figure 3H) when compared to the subgroup 10 cell line HCC202. Moreover, when both cell lines were grown under hypoxic conditions, AU565 cells exhibited a much stronger hypoxia response as measured by higher mRNA levels of several hypoxia-inducible genes (Additional file 4, Figure S7). This is consistent with an enhanced hypoxia response and higher hypoxia-inducible factor levels in the subgroup 7 breast tumors.
HER2 related subtypes exhibit basal and luminal type cell features
It has been reported that basal-like breast tumors have a significantly more robust hypoxia response compared to luminal-type breast cancers , similar to the differential hypoxia response between the two identified HER2 related subgroups. In addition, basal breast tumors are characterized by low ER and PR expression as well as high levels of EGFR [1, 46], elevated expression of AKT3 , CD44 and MET  and decreased GATA3 expression ; features that are also consistent with each of the HER2 related subgroups investigated in the current study. Recent studies have also suggested that HER2 related tumors may be derived from late luminal progenitor cells, an intermediate differentiation stage between luminal progenitor cells from which basal-like tumors arise and differentiated luminal cells from which luminal A and B tumors arise [58, 59]. Therefore, we next investigated whether the two HER2 related subgroups might have additional characteristics consistent with basal and luminal tumors in addition to being HER2+.
We first measured differences in expression patterns of pathways and genes known to be differentially expressed in basal and luminal breast tumors. As shown in Figure 2 (and Additional file 4, Figure S3), tumors in subgroup 7 exhibit significantly higher hypoxia and EGFR pathway activity and lower levels of ER and PR pathway activity [55, 60] when compared to tumors in subgroup 10, suggesting that tumors in subgroup 7 and 10 tumors, in addition to being HER2+, have basal-like and luminal-like patterns of pathway activities, respectively. Next, we compared the expression of several basal and luminal specific genes in the two HER2 related subgroups. As shown in Figure 4 (and validated in Additional file 4, Figure S8), subtype 7 tumors have significantly higher expression of basal related genes (Supplemental methods) including AKT3, MET, EGFR, CD44 and basal type cytokeratins (KRT5 and KRT17) (Figure 4A, and Additional file 4, Figures S8, S9, S10). In contrast, subgroup 10 tumors have higher expression of the luminal-specific master regulator GATA3 (Figure 4B) [61, 62]. Consistent with this observation, HCC202 cells (subgroup 10) had significantly higher mRNA levels of GATA3 (Figure 4C) than AU565 cells (subgroup 7). Although these results demonstrate that subgroup 10 tumors do not express the majority of basal-type markers and do exhibit certain luminal-type features, these tumors did not have higher levels of luminal-type cytokeratins (KRT18 and KRT19) or several other luminal-specific markers (Additional file 4, Figure S11).
The differential expression of a number of luminal and basal-specific genes and pathways in subgroups 7 and 10 prompted us to further investigate this possibility at the genome-wide level. We used the expression of a well-annotated dataset of 24 luminal and 24 basal breast cancer cell lines to develop a gene expression signature with the ability to differentiate between basal and luminal tumor characteristics (Additional file 4, Figure S12A). The basal-luminal signature was first validated using a leave-one-out cross-validation among the cancer cell lines (Additional file 4, Figure S12B). The accuracy of the signature was then investigated by predicting the basal or luminal characteristics of tumors and then comparing the probability of basal or luminal status against the intrinsic subtype based on previously reported classifications of tumors  in both the primary and validation datasets. Our analysis determined that most basal tumors were strongly predicted to be basal-like (mean predicted probability: 0.933 and 0.992, primary and validation datasets, respectively). The majority of luminal tumors were predicted to have luminal-like characteristics (mean predicted probability: 0.321 and 0.438, respectively), albeit with a less uniform distribution of predicted probabilities as compared to basal tumors, likely owing to the high degree of heterogeneity in luminal breast tumors. However, the difference in the predicted probability between the intrinsic subtypes was still highly significant for both the primary and validation datasets (P <0.0001) (Figure 4D, 4E). In contrast to basal and luminal tumors, ERBB2+ tumors showed a wide-range of predicted basal-luminal characteristics (Figure 4D, 4E). However, when we segregated these tumors into the two pathway-based HER2 related subgroups (subgroups 7 and 10), we found that subgroup 10 tumors were predicted to be predominantly luminal in both the primary (congruence rate: 83%) and validation (congruence rate: 60%) datasets, while most tumors in subgroup 7 were predicted to be basal in the primary (congruence rate: 72%) and validation (congruence rate: 81.1%) datasets when a predicted probability of 0.5 is used as the cut-off for class assignment (Figure 4F, 4G) and to determine the level of agreement between 7/10 subgroup identity and basal/luminal characteristics (congruence rate).
To further examine the similarity between characteristics of each of the two HER2 related subgroups and basal/luminal cancer cell lines, we also performed the reciprocal analysis. We developed a gene expression signature, using the samples from the primary dataset as the training data, to differentiate between subgroup 7 and 10 tumor characteristics (Additional file 4, Figure S13). Among the genes identified in this signature, we noted, consistent with our earlier analyses, GATA3 and MET enrichment in subgroup 10 and 7 tumors, respectively, both of which are known to be preferentially expressed in luminal and basal tumor, respectively (Additional file 4, Figure S13). Although the intent of this signature is not to assign samples to subgroups 7 and 10, but instead to investigate differences in group characteristics, we first validated the ability of this signature to identify subgroup 7 and 10 characteristics using samples in the validation dataset. As illustrated in Figure 4H, even with the inclusion of samples with a lower probability (<0.70) of subgroup assignment, we demonstrate that the majority of tumors assigned to subgroups 7 were predicted to have subgroup 7-like characteristics (mean predicted probability: 0.7950) while samples assigned to subgroup 10 based on patterns of pathway activity were determined to have subgroup 10-like characteristics (mean predicted probability: 0.2022) with an overall congruence rate of 84.5% (P <0.0001). This signature and training model were then used to investigate the subgroup 7 and subgroup 10 characteristics of the breast cancer cell lines with established basal and luminal status. We found that the majority of the luminal and basal cancer cells were characteristic of subgroups 10 and 7, respectively (congruence rate 81.3%, P <0.0001), when a cutoff of 0.5 was used for class assignment and to calculate congruence rate (Figure 4I). Taken together, these data strongly suggested that subgroup 10 and subgroup 7 tumors have molecular characteristics that are consistent with luminal and basal type breast cancers, respectively. These results, consistent with several previous studies [3, 11, 63], suggest that HER2 related tumor heterogeneity may arise as a consequence of basal and luminal cell features and provide evidence supporting the hypothesis that the cell of origin for subgroup 7 tumors may have basal-like features, whereas subgroup10 tumors may arise from progenitor cells that have further differentiated and have features more consistent with luminal-like tumors.
HER2 related tumors demonstrate significant positive correlations between hypoxia response, basal gene expression, and HIF gene expression
Given that our results demonstrate a clear correlation among basal-like features, hypoxia response and HIF gene expression when comparing tumors in subgroups 7 and 10, we next investigated whether these observations can be further extended to the heterogeneity among HER2 related tumors irrespective of the classification strategy. Since different classification strategies lead to consistent but imperfect agreements sensitive to varying thresholds, we sought to further determine the relationship between the continuous variables of these properties. To investigate these correlations, we first compared the average expression of basal associated genes (detailed in Figure 4A) with the average expression of hypoxia-related genes (detailed in Figure 3A). As expected from our previous results, a significant, positive relationship was found to exist between average hypoxia-and basal-related gene expression in tumors assigned to subgroups 7 and 10 in both the primary (Figure 5A, P <0.0001) and validation (Figure 5F, P <0.0001) datasets. In order to investigate this relationship among HER2 related tumors, we next examined tumors classified as HER2 related (or ERBB2+) by the intrinsic gene list [61, 62] (Figure 5B, G) as well as in those tumors that were reported to be HER2+ by IHC (Figure 5C, H). As illustrated in Figures 5B (P <0.0001) and 5C (P <0.0001) for the primary dataset, and validated in Figures 5G (P = 0.0012) and 5H (P = 0.0069), this relationship is evident regardless of which of the three independent methods is used to classify tumors as HER2+ or HER2 related.
Furthermore, our data suggest that differential tumor hypoxia response is associated in subgroups 7 and 10 tumors with expression of HIF mRNA (Figure 3). Among the tumors classified as HER2+ by IHC, there is a significant correlation between both HIF1A (HIF-1α) (Figure 5D, P = 0.0256) and EPAS1 (HIF-2α) (Figure 5E, P = 0.0038) expression and hypoxia-related gene expression. These results are further confirmed in the validation dataset as illustrated in Figures 5I (P = 0.0032) and 5J (P = 0.0094). Collectively, these data suggest that patterns of basal-related gene expression, as well as HIF levels, directly correlate with the expression of hypoxia-related genes in HER2 related tumors regardless of whether tumors are defined as HER2 related by a HER2 gene expression signature (subgroups 7/10), by the intrinsic gene list (ERBB2+), or by HER2 IHC status.
Basal-like breast tumors exhibit an exaggerated hypoxia response and increased HIF-1α expression but not lower pO2
Our data suggest that tumors assigned to subgroup 7, which express high levels of HIF-1α and hypoxia-related genes, and subgroup 10, which express lower levels of HIF-1α and hypoxia-related genes, have characteristics in common with basal and luminal breast tumors, respectively. Since our results demonstrate that the correlation between average basal gene expression and average hypoxia-related gene expression is consistent across HER2 related tumors irrespective of classification strategy, we next investigated whether these relationships were evident across all breast tumors given that basal-like tumors have been previously reported to be characterized by an enhanced hypoxia response [55, 64]. As illustrated in Figure 6, a direct correlation was observed in both the primary (Figure 6A, P <0.0001) and validation (Figure 6D, P <0.0001) datasets between basal-like and hypoxia-related gene expression. We further examined the correlation between HIF-1α expression and both hypoxia-related gene expression (Figure 6B, E) and basal-like gene expression (Figure 6C, F) across both the primary and validation dataset and demonstrated a significant (P <0.0001) relationship between both patterns of gene expression. These results, consistent with previous studies [55, 64] collectively suggest that as tumors express more basal-like features, irrespective of HER2 status, they also demonstrate higher levels of hypoxia-related gene expression as well as higher levels of HIF-1α. Consistent with the latter point, we demonstrate that tumors assigned to subgroups 7 or that are classified as basal using the intrinsic gene list express higher levels of HIF-1α as compared to tumors assigned to subgroup 10 or classified as luminal, respectively (Additional file 4, Figure S14, P <0.0001).
Given that significantly higher mRNA levels of HIF-1α expression were noted in both HER2 related tumors and basal tumors with high levels of predicted hypoxia response, we next investigated whether the observed levels of predicted hypoxia response correlated with physical levels of oxygenation (pO2). We took advantage of data from our previously published study that measured pO2 levels in situ in human breast cancer using an Eppendorf electrode . As shown in Figure S15, basal tumors in this cohort also showed significantly higher levels of hypoxia-induced gene expression (P = 0.0065) (Additional file 4, Figure. S15A) and HIF-1α mRNA (P = 0.0261) (Additional file 4 Figure S15B) when compared to luminal tumors. Importantly, no significant difference (P = 0.28) was observed for pO2 levels (Additional file 4, Figure S15C). Therefore, these data demonstrate that the exaggerated hypoxia response observed in basal-like tumors cannot be solely attributed to decreased oxygen levels. Instead, we postulated that intrinsically higher mRNA expression of hypoxia regulators, including HIF-1α, was responsible for the exaggerated hypoxia response in basal-type tumors.
Since basal-like breast tumors and HER2 related tumors with basal-like characteristics, including subgroup 7 tumors, demonstrated increased levels of HIF-1α, we investigated whether such differences are a cell-lineage phenomenon which can also delineate normal basal and luminal mammary epithelial cells. Primary breast epithelial cells were separated based on surface expression of EPCAM (TACSTD1) from reduction mammoplasty specimens and gene expression levels were determined by microarray . Both arrays and real-time PCR showed that basal epithelial cells had a greater than three-fold higher level of both HIF-1α and HIF-2α expression compared to isogenic luminal epithelial cells (Figure 6G). When placed under hypoxic conditions, primary basal epithelial cells also demonstrated a significantly greater hypoxia response (Figure 6H). From these data, we hypothesize that the stronger hypoxia response is associated with higher expression of hypoxia regulators shared by both basal-like and HER2 related tumors with basal-like characteristics, including subgroup 7 tumors. Given that normal basal epithelial cells also demonstrate enhanced HIF-1α expression, the association may be a cell-lineage distinction that can be traced to systematic differences in the composition of hypoxia regulators between normal, non-transformed basal and luminal breast epithelial cells.
HIF-1α mRNA is responsible for the strong hypoxia response in basal cells
Since our results demonstrated a correlation between predicted hypoxia response and HIF-1α mRNA expression in both basal breast tumors and HER2 related tumors with basal-like characteristics, including subgroup 7 tumors, we next investigated whether HIF-1α mRNA expression directly regulated the exaggerated hypoxia-induced gene expression in these cells. We analyzed HIF1A (HIF-1α expression in six common breast cancer cell lines to determine whether basal breast cell lines, which have been previously reported to have a higher hypoxia response , also express higher levels of HIF1A (HIF-1α) and EPAS1 (HIF-2α) mRNA as compared to luminal breast cancer cell lines. We determined that the basal breast cancer cell lines (BT20, MDA-MB231, MDA-MB157) had consistently higher mRNA levels of both HIF1A (HIF-1α) (Figure 7A) and EPAS1 (HIF-2α) (Additional file 4, Figure S16A) when compared to the luminal cell lines (MCF7, BT474, CAMA1). Furthermore, we found that the basal breast cancer cell lines also exhibited elevated expression levels of three hypoxia-inducible genes (VEGF, GLUT1, and CA9) even under normoxia (Figure 7B and Additional file 4, Figure S16B, C). These results are consistent with previous studies reporting a strong constitutive hypoxia response and glycolytic phenotype in basal breast cancer cell lines [64, 65].
To investigate the roles of HIF1A (HIF-1α) and EPAS1 (HIF-2α) in the hypoxia response, we reduced expression of both genes by siRNA-mediated gene silencing in MCF7 (luminal) and MDA-MB231 (basal) cells. This gene silencing significantly reduced both transcripts in MDA-MB231 cells to a level comparable to those in MCF7 cells (Additional file 4, Figure S17A, B). Silencing of HIF-1α, but not HIF-2α, mRNA significantly reduced the exaggerated induction of VEGF (Figure 7C), GLUT1, and CA9 expression (Additional file 4, Figure S17C, D) in MDA-MB231 cells under hypoxia (1% pO2). Importantly, the silencing of HIF1A (HIF-1α), but not EPAS1 (HIF-2α), in MDA-MB231 reduced the induction of these hypoxia-inducible genes to the levels comparable to MCF7 cells (Figure 7C and Additional file 4, Figure S17). Finally, we examined whether a similar mechanism may be also important for the stronger hypoxia response in the breast cancer cell lines in subgroup 7. Consistent with our findings in the basal-breast cancer cell line, AU565 cells (subgroup 7) demonstrated that silencing of HIF1A (HIF-1α), but not EPAS1 (HIF-2α), led to the abolishment of the exaggerated hypoxia response seen in these cells (Figure 7D and Additional file 4, Figure S18). Taken together, these data are consistent with the role of high levels of HIF-1α regulating the exaggerated hypoxia response in both basal-like cancers and HER2 related tumors with basal-like characteristics and further suggest that there is an apparent strong lineage-specific component to this phenotype.
Breast tumors are currently characterized in the clinic based on tumor size, visual characteristics, and a limit number of histochemical markers including estrogen receptor, progesterone receptor, and HER2 receptor status. Although these measures provide information about tumor properties, they offer minimal insight into the underlying biological mechanisms and provide only limited guidance in the development of therapeutic strategies. This can be best illustrated by the differential response of patients to treatment with the HER2 inhibitor Herceptin; while HER2 negative patients do not response to Herceptin, only fraction of HER2+ patients respond to this treatment suggesting that additional unrecognized heterogeneity must exist within these patients and that more sophisticated strategies must be employed to investigate tumor heterogeneity and to develop therapeutic regimens.
Consistent with this hypothesis, a number of recent studies have illustrated the heterogeneity of tumors driven by HER2 amplification. These studies have demonstrated variability in patterns of copy number variation, DNA methylation, and global gene expression [11, 12]. Moreover, molecular subtypes based on the analysis of gene expression patterns  as well as more complex analyses of patterns of oncogenic and tumor suppressor pathway activities  have identified classes that can subdivide HER2 related tumors and provide further insight into the molecular and clinical heterogeneity of HER2 related tumors beyond the ERBB2+ molecular classification. As previously discussed, we recently reported a quantitative tumor classification scheme based on patterns of pathway-associated signatures that can better match the complexity of the tumor, and in doing so, identified two subgroups characterized by high levels of HER2 pathway activity . In the current study, we have utilized this framework to analyze non-genetic factors in the context of genetic factors to define the relationship between these two distinct sources of tumor heterogeneity. We found a high degree of specific co-regulation between non-genetic stresses and oncogenic signaling pathways. Many of the co-regulated and correlated pathways observed in the tumors are consistent with previous findings to validate our approach. For example, the strong co-regulation of hypoxia/EGFR/TNFα/TGFβ/STAT3 pathways clearly distinguishes between two groups of HER2 related tumors with significant differences in the pathway activities of these pathways. By analyzing patterns of non-genetic stress factors, we demonstrate significant variation in patterns of predicted tumor microenvironment activity in two HER2 related subgroups classified by patterns of oncogenic pathway activity. The HER2 related tumors that are characterized by high levels of predicted hypoxia response display enhanced levels of hypoxia-related genes including HIF-1α, exhibit similar patterns of pathway activity and gene expression as compared to basal-like breast tumors. In contrast, the HER2 related tumors that express low levels of hypoxia have oncogenic features in common with luminal-like tumors. As illustrated in Figure 5, the correlation between hypoxia-related gene expression, HIF gene expression, and basal-like features is common across HER2 related tumors irrespective of the classification strategy used to identify these tumors; tumors classified by patterns of oncogenic pathway activity (subgroups 7/10), by the intrinsic gene list (ERBB2+) or by IHC (HER2+) all demonstrate this characteristic. Moreover, we demonstrate that breast tumors in general demonstrate this pattern of coincidental gene expression where the increased expression of basal-like genes correlates with the expression of HIF and hypoxia-related genes. These results are supported by, and are consistent with, several studies that suggest that the observed heterogeneity in ER status among HER2 related tumors is consistent with differences in the basal and luminal origins of HER2-driven breast cancers [3, 11, 63] as well as additional studies which report that basal-like tumors have higher levels of the hypoxia gene expression program compared to luminal breast tumors [55, 64]. Finally we demonstrate that HIF-1α expression levels are higher in basal-like tumors, HER2 related tumors with high hypoxia, and normal basal breast epithelial cells and that HIF-1α regulates the expression of hypoxia-induced genes. Our data support the idea that the differences in the hypoxic response and mRNA level of hypoxia regulators may be an intrinsic cell lineage property associated with the origin of the tumor progenitor cell. Consistent with this idea, it was recently suggested that HER2 related tumors arise from late luminal progenitor cells, an intermediate differentiation stage between luminal progenitor cells from which basal-like tumors arise and differentiated luminal cells from which luminal A and B tumors arise . Our results, in conjunction with others  suggest that some HER2 related tumors may arise from early luminal progenitor cells that are more basal-like in nature whereas other HER2 related tumors may arise from more differentiated, late luminal progenitor cells resulting in tumors with more luminal-like characteristics.
The positive correlation between both HIF-1α and EPAS1 (HIF-2α) mRNA levels and overall breast tumor hypoxia response suggests a role for transcriptional regulation of these hypoxia regulators. Our findings which demonstrate that elevated expression of HIF-1α mRNA in the basal and HER2 related breast tumors with basal-like features, including subgroup 7 tumors, correlates with hypoxia response are consistent with recent studies reporting higher mRNA levels of HIF-1α and EPAS1 (HIF-2α) in glioma stem cells and renal epithelial cells [24, 25]. The consistent pattern of higher mRNA levels of hypoxia-inducible factors in the less differentiated glioma stem cells , hematopoietic stem cells [67, 68] and basal mammary epithelial cells (this study) suggests that the differentiation status of these cells can be coupled with the fine-tuning of the hypoxia response and that these changes may be responsible for their distinct sensitivity to oxygen and metabolic requirements. It has also been reported that pO2 and hypoxia-inducible factors play important roles in the differentiation processes of several different cell types [69, 70] and has been shown to induce de-differentiation of breast cancer cells promoting invasion .
Our findings also have translation potential and important therapeutic implications. Gaining a more complete understanding of the biological mechanisms regulating oncogenesis will enable the development of more rational therapeutic strategies to enhance therapeutic efficacy. The cellular origin of HER2 related tumors, the patterns of oncogenic pathway activity, and the contributions of microenvironmental stresses may be clinically relevant since differences in levels of PI3K, EGFR, hypoxia, or MET activities in tumors have been reported to contribute to various therapeutic responses to HER2-targeting and other therapies, such as Trastuzumab (Herceptin). The influence of these factors on response can be direct as in the case of PI3K [59, 72], or indirect, as in the case of development of drug resistance under hypoxia . Aberrant MET expression has also been shown to confer resistance to the EGFR inhibitor, Gefitinib . Therefore, co-activation of these pathways in subgroup 7 tumors provides further justification for the development and investigation of therapeutic regimens that incorporate compounds targeting more than one oncogenic pathway in order to increase efficacy and reduce resistance.
We have incorporated non-genetic (microenvironmental stresses) with genetic (oncogenic and tumor suppressor genes) factors together in a pathway-based classification framework to further dissect and identify biological mechanisms underlying tumor heterogeneity. We validate this approach by investigating the basis for the distinct patterns of oncogenic pathway activity and microenvironmental conditions of two HER2 related tumor subgroups and found significant similarities with many features distinguishing between the basal and luminal intrinsic subtypes of breast cancers. By examining tumors identified as HER2 related by patterns of oncogenic pathway activity (subgroups7/10), tumors defined as HER2 related by the intrinsic gene list (ERBB2+), or tumors that are HER2+ by IHC, we demonstrate that these features are common to HER2 related tumors irrespective of the strategy used to characterize them as HER2+. By revealing the role of HIF-1α in regulating the hypoxia response and showing that this mechanism, as well as the expression of additional features, link HER2 related tumors in subgroup 7 with basal-like breast tumors and subgroup 10 tumors with luminal like breast tumors, respectively, we identify additional oncogenic mechanisms and HER2 related tumor heterogeneity that was not previously evident. While our current studies focus only on a subset of chemical stresses, a similar approach is likely to be applicable for incorporating the mechanical factors, physical forces, interstitial pressure, paracrine milieu and other non-genetic conditions that may impact tumor phenotypes. Such strategies are scalable to incorporate additional sources of information to investigate and discover cross-talk between different pathways.
aryl hydrocarbonreceptor nuclear translocator
Baysian Factor Regression Modeling
carbonic anhydrase IX
deleted in esophageal cancer
egl nin homolog 3
epidermal growth factor receptor
e1A binding protein p300
endothelial PAS domain containing protein 1
v-erb-B2 erythroblastic leukemia viral oncogene homolog 2
GATA binding protein 3
glucose transporter type 1
Human epidermal growth factor receptor 2
hypoxia-inducible factor 1 alpha
hypoxia-inducible factor 2 alpha
insulin growth factor 1 receptor
Jun proto oncogene
17, 18, 19, cytokeratin 5, 17, 18, 19
signal transducer and activator of transduction 3
tumor growth factor beta
tumor necrosis factor alpha
thioredoxin interacting protein
vascular endothelial growth factor A.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lonning PE, Borresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752. 10.1038/35021093.
Allinen M, Beroukhim R, Cai L, Brennan C, Lahti-Domenici J, Huang H, Porter D, Hu M, Chin L, Richardson A, Schnitt S, Sellers WR, Polyak K: Molecular characterization of the tumor microenvironment in breast cancer. Cancer Cell. 2004, 6: 17-32. 10.1016/j.ccr.2004.06.010.
Chin K, DeVries S, Fridlyand J, Spellman PT, Roydasgupta R, Kuo WL, Lapuk A, Neve RM, Qian Z, Ryder T, Chen F, Feiler H, Tokuyasu T, Kingsley C, Dairkee S, Meng Z, Chew K, Pinkel D, Jain A, Ljung BM, Esserman L, Albertson DG, Waldman FM, Gray JW: Genomic and transcriptional aberrations linked to breast cancer pathophysiologies. Cancer Cell. 2006, 10: 529-541. 10.1016/j.ccr.2006.10.009.
Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, Clark L, Bayani N, Coppe JP, Tong F, Speed T, Spellman PT, DeVries S, Lapuk A, Wang NJ, Kuo WL, Stilwell JL, Pinkel D, Albertson DG, Waldman FM, McCormick F, Dickson RB, Johnson MD, Lippman M, Ethier S, Gazdar A, Gray JW: A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006, 10: 515-527. 10.1016/j.ccr.2006.10.008.
Visvader JE: Keeping abreast of the mammary epithelial hierarchy and breast tumorigenesis. Genes Dev. 2009, 23: 2563-2577. 10.1101/gad.1849509.
Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL: Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science. 1987, 235: 177-182. 10.1126/science.3798106.
Baselga J, Tripathy D, Mendelsohn J, Baughman S, Benz CC, Dantis L, Sklarin NT, Seidman AD, Hudis CA, Moore J, Rosen PP, Twaddell T, Henderson IC, Norton L: Phase II study of weekly intravenous trastuzumab (Herceptin) in patients with HER2/neu-overexpressing metastatic breast cancer. Semin Oncol. 1999, 26: (4 Suppl 12):78-83.
Piccart-Gebhart MJ, Procter M, Leyland-Jones B, Goldhirsch A, Untch M, Smith I, Gianni L, Baselga J, Bell R, Jackisch C, Cameron D, Dowsett M, Barrios CH, Steger G, Huang CS, Andersson M, Inbar M, Lichinitser M, Láng I, Nitz U, Iwata H, Thomssen C, Lohrisch C, Suter TM, Rüschoff J, Suto T, Greatorex V, Ward C, Straehle C, McFadden E, et al: Trastuzumab after adjuvant chemotherapy in HER2-positive breast cancer. N Engl J Med. 2005, 353: 1659-1672. 10.1056/NEJMoa052306.
Romond EH, Perez EA, Bryant J, Suman VJ, Geyer CE, Davidson NE, Tan-Chiu E, Martino S, Paik S, Kaufman PA, Swain SM, Pisansky TM, Fehrenbacher L, Kutteh LA, Vogel VG, Visscher DW, Yothers G, Jenkins RB, Brown AM, Dakhil SR, Mamounas EP, Lingle WL, Klein PM, Ingle JN, Wolmark N: Trastuzumab plus adjuvant chemotherapy for operable HER2-positive breast cancer. N Engl J Med. 2005, 353: 1673-1684. 10.1056/NEJMoa052122.
Smith IE: New drugs for breast cancer. Lancet. 2002, 360: 790-792. 10.1016/S0140-6736(02)09903-8.
Staaf J, Jonsson G, Ringner M, Vallon-Christersson J, Grabau D, Arason A, Gunnarsson H, Agnarsson BA, Malmstrom PO, Johannsson OT, Loman N, Barkardottir RB, Borg A: High-resolution genomic and expression analyses of copy number alterations in HER2-amplified breast cancer. Breast Cancer Res. 2010, 12: R25-10.1186/bcr2568.
Holm K, Hegardt C, Staaf J, Vallon-Christersson J, Jonsson G, Olsson H, Borg A, Ringner M: Molecular subtypes of breast cancer are associated with characteristic DNA methylation patterns. Breast Cancer Res. 2010, 12: R36-10.1186/bcr2590.
Staaf J, Ringner M, Vallon-Christersson J, Jonsson G, Bendahl PO, Holm K, Arason A, Gunnarsson H, Hegardt C, Agnarsson BA, Luts L, Grabau D, Ferno M, Malmstrom PO, Johannsson OT, Loman N, Barkardottir RB, Borg A: Identification of subtypes in human epidermal growth factor receptor 2--positive breast cancer reveals a gene signature prognostic of outcome. J Clin Oncol. 2010, 28: 1813-1820. 10.1200/JCO.2009.22.8775.
Gatza ML, Lucas JE, Barry WT, Kim JW, Wang Q, Crawford MD, Datto MB, Kelley M, Mathey-Prevot B, Potti A, Nevins JR: A pathway-based classification of human breast cancer. Proc Natl Acad Sci USA. 2010, 107: 6994-6999. 10.1073/pnas.0912708107.
Vaupel P: Tumor microenvironmental physiology and its implications for radiation oncology. Semin Radiat Oncol. 2004, 14: 198-206. 10.1016/j.semradonc.2004.04.008.
Vaupel P, Hockel M: Blood supply, oxygenation status and metabolic micromilieu of breast cancers: characterization and therapeutic relevance. Int J Oncol. 2000, 17: 869-879.
Gatenby RA, Gillies RJ: A microenvironmental model of carcinogenesis. Nat Rev Cancer. 2008, 8: 56-61. 10.1038/nrc2255.
Harris AL: Hypoxia--a key regulatory factor in tumour growth. Nat Rev Cancer. 2002, 2: 38-47. 10.1038/nrc704.
Motzer RJ, Bander NH, Nanus DM: Renal-cell carcinoma. N Engl J Med. 1996, 335: 865-875. 10.1056/NEJM199609193351207.
Blagosklonny MV, An WG, Romanova LY, Trepel J, Fojo T, Neckers L: p53 inhibits hypoxia-inducible factor-stimulated transcription. J Biol Chem. 1998, 273: 11995-11998. 10.1074/jbc.273.20.11995.
Zundel W, Schindler C, Haas-Kogan D, Koong A, Kaper F, Chen E, Gottschalk AR, Ryan HE, Johnson RS, Jefferson AB, Stokoe D, Giaccia AJ: Loss of PTEN facilitates HIF-1-mediated gene expression. Genes Dev. 2000, 14: 391-396.
Zhong H, Chiles K, Feldser D, Laughner E, Hanrahan C, Georgescu MM, Simons JW, Semenza GL: Modulation of hypoxia-inducible factor 1alpha expression by the epidermal growth factor/phosphatidylinositol 3-kinase/PTEN/AKT/FRAP pathway in human prostate cancer cells: implications for tumor angiogenesis and therapeutics. Cancer Res. 2000, 60: 1541-1545.
Laughner E, Taghavi P, Chiles K, Mahon PC, Semenza GL: HER2 (neu) signaling increases the rate of hypoxia-inducible factor 1alpha (HIF-1alpha) synthesis: novel mechanism for HIF-1-mediated vascular endothelial growth factor expression. Mol Cell Biol. 2001, 21: 3995-4004. 10.1128/MCB.21.12.3995-4004.2001.
Chi JT, Wang Z, Nuyten DS, Rodriguez EH, Schaner ME, Salim A, Wang Y, Kristensen GB, Helland A, Borresen-Dale AL, Giaccia A, Longaker MT, Hastie T, Yang GP, van de Vijver MJ, Brown PO: Gene expression programs in response to hypoxia: cell type specificity and prognostic significance in human cancers. PLoS Med. 2006, 3: e47-10.1371/journal.pmed.0030047.
Li Z, Bao S, Wu Q, Wang H, Eyler C, Sathornsumetee S, Shi Q, Cao Y, Lathia J, McLendon RE, Hjelmeland AB, Rich JN: Hypoxia-inducible factors regulate tumorigenic capacity of glioma stem cells. Cancer Cell. 2009, 15: 501-513. 10.1016/j.ccr.2009.03.018.
Schroeder T, Yuan H, Viglianti BL, Peltz C, Asopa S, Vujaskovic Z, Dewhirst MW: Spatial heterogeneity and oxygen dependence of glucose consumption in R3230Ac and fibrosarcomas of the Fischer 344 rat. Cancer Res. 2005, 65: 5163-5171. 10.1158/0008-5472.CAN-04-3900.
Sonveaux P, Vegran F, Schroeder T, Wergin MC, Verrax J, Rabbani ZN, De Saedeleer CJ, Kennedy KM, Diepart C, Jordan BF, Kelley MJ, Gallez B, Wahl ML, Feron O, Dewhirst MW: Targeting lactate-fueled respiration selectively kills hypoxic tumor cells in mice. J Clin Invest. 2008, 118: 3930-3942.
Shaw RJ, Kosmatka M, Bardeesy N, Hurley RL, Witters LA, DePinho RA, Cantley LC: The tumor suppressor LKB1 kinase directly activates AMP-activated kinase and regulates apoptosis in response to energy stress. Proc Natl Acad Sci USA. 2004, 101: 3329-3335. 10.1073/pnas.0308061100.
Yun J, Rago C, Cheong I, Pagliarini R, Angenendt P, Rajagopalan H, Schmidt K, Willson JK, Markowitz S, Zhou S, Diaz LA, Velculescu VE, Lengauer C, Kinzler KW, Vogelstein B, Papadopoulos N: Glucose deprivation contributes to the development of KRAS pathway mutations in tumor cells. Science. 2009, 325: 1555-1559. 10.1126/science.1174229.
Mueller-Klieser W, Walenta S: Geographical mapping of metabolites in biological tissue with quantitative bioluminescence and single photon imaging. Histochem J. 1993, 25: 407-420. 10.1007/BF00157805.
Kallinowski F, Schlenger KH, Runkel S, Kloes M, Stohrer M, Okunieff P, Vaupel P: Blood flow, metabolism, cellular microenvironment, and growth rate of human tumor xenografts. Cancer Res. 1989, 49: 3759-3764.
Dewhirst MW, Klitzman B, Braun RD, Brizel DM, Haroon ZA, Secomb TW: Review of methods used to study oxygen transport at the microcirculatory level. Int J Cancer. 2000, 90: 237-255. 10.1002/1097-0215(20001020)90:5<237::AID-IJC1>3.0.CO;2-T.
Vujaskovic Z, Rosen EL, Blackwell KL, Jones EL, Brizel DM, Prosnitz LR, Samulski TV, Dewhirst MW: Ultrasound guided pO2 measurement of breast cancer reoxygenation after neoadjuvant chemotherapy and hyperthermia treatment. Int J Hyperthermia. 2003, 19: 498-506. 10.1080/0265673031000121517.
Brizel DM, Schroeder T, Scher RL, Walenta S, Clough RW, Dewhirst MW, Mueller-Klieser W: Elevated tumor lactate concentrations predict for an increased risk of metastases in head-and-neck cancer. Int J Radiat Oncol Biol Phys. 2001, 51: 349-353.
Schwickert G, Walenta S, Sundfor K, Rofstad EK, Mueller-Klieser W: Correlation of high lactate levels in human cervical cancer with incidence of metastasis. Cancer Res. 1995, 55: 4757-4759.
Walenta S, Salameh A, Lyng H, Evensen JF, Mitze M, Rofstad EK, Mueller-Klieser W: Correlation of high lactate levels in head and neck tumors with incidence of metastasis. Am J Pathol. 1997, 150: 409-415.
Walenta S, Wetterling M, Lehrke M, Schwickert G, Sundfor K, Rofstad EK, Mueller-Klieser W: High lactate levels predict likelihood of metastases, tumor recurrence, and restricted patient survival in human cervical cancers. Cancer Res. 2000, 60: 916-921.
Walenta S, Mueller-Klieser WF: Lactate: mirror and motor of tumor malignancy. Semin Radiat Oncol. 2004, 14: 267-274. 10.1016/j.semradonc.2004.04.004.
Koukourakis MI, Giatromanolaki A, Sivridis E, Simopoulos K, Pastorek J, Wykoff CC, Gatter KC, Harris AL: Hypoxia-regulated carbonic anhydrase-9 (CA9) relates to poor vascularization and resistance of squamous cell head and neck cancer to chemoradiotherapy. Clin Cancer Res. 2001, 7: 3399-3403.
Winter SC, Buffa FM, Silva P, Miller C, Valentine HR, Turley H, Shah KA, Cox GJ, Corbridge RJ, Homer JJ, Musgrove B, Slevin N, Sloan P, Price P, West CM, Harris AL: Relation of a hypoxia metagene derived from head and neck cancer to prognosis of multiple cancers. Cancer Res. 2007, 67: 3441-3449. 10.1158/0008-5472.CAN-06-3322.
Bild AH, Yao G, Chang JT, Wang Q, Potti A, Chasse D, Joshi MB, Harpole D, Lancaster JM, Berchuck A, Olson JA, Marks JR, Dressman HK, West M, Nevins JR: Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature. 2006, 439: 353-357. 10.1038/nature04296.
Chen JL, Lucas JE, Schroeder T, Mori S, Wu J, Nevins J, Dewhirst M, West M, Chi JT: The genomic analysis of lactic acidosis and acidosis response in human cancers. PLoS Genet. 2008, 4: e1000293-10.1371/journal.pgen.1000293.
Chi JT, Rodriguez EH, Wang Z, Nuyten DS, Mukherjee S, van de Rijn M, van de Vijver MJ, Hastie T, Brown PO: Gene expression programs of human smooth muscle cells: tissue-specific differentiation and prognostic significance in breast cancers. PLoS Genet. 2007, 3: 1770-1784.
Chen JL, Merl D, Peterson CW, Wu J, Liu PY, Yin H, Muoio DM, Ayer DE, West M, Chi JT: Lactic acidosis triggers starvation response with paradoxical induction of TXNIP through MondoA. PLoS Genet. 2010, 6: e1001093-10.1371/journal.pgen.1001093.
Lucas JE, Kung HN, Chi JT: Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers. PLoS Comput Biol. 2010, 6: e1000920-10.1371/journal.pcbi.1000920.
Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Eystein Lonning P, Borresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001, 98: 10869-10874. 10.1073/pnas.191367098.
Huper G, Marks JR: Isogenic normal basal and luminal mammary epithelial isolated by a novel method show a differential response to ionizing radiation. Cancer Res. 2007, 67: 2990-3001. 10.1158/0008-5472.CAN-06-4065.
Saito S, Furuno A, Sakurai J, Sakamoto A, Park HR, Shin-Ya K, Tsuruo T, Tomida A: Chemical genomics identifies the unfolded protein response as a target for selective cancer cell killing during glucose deprivation. Cancer Res. 2009, 69: 4225-4234. 10.1158/0008-5472.CAN-08-2689.
Sinn E, Muller W, Pattengale P, Tepler I, Wallace R, Leder P: Coexpression of MMTV/v-Ha-ras and MMTV/c-myc genes in transgenic mice: synergistic action of oncogenes in vivo. Cell. 1987, 49: 465-475. 10.1016/0092-8674(87)90449-1.
Nielsen TO, Hsu FD, Jensen K, Cheang M, Karaca G, Hu Z, Hernandez-Boussard T, Livasy C, Cowan D, Dressler L, Akslen LA, Ragaz J, Gown AM, Gilks CB, van de Rijn M, Perou CM: Immunohistochemical and clinical characterization of the basal-like subtype of invasive breast carcinoma. Clin Cancer Res. 2004, 10: 5367-5374. 10.1158/1078-0432.CCR-04-0220.
D'Cruz CM, Gunther EJ, Boxer RB, Hartman JL, Sintasath L, Moody SE, Cox JD, Ha SI, Belka GK, Golant A, Cardiff RD, Chodosh LA: c-MYC induces mammary tumorigenesis by means of a preferred pathway involving spontaneous Kras2 mutations. Nat Med. 2001, 7: 235-239. 10.1038/84691.
Franovic A, Gunaratnam L, Smith K, Robert I, Patten D, Lee S: Translational up-regulation of the EGFR by tumor hypoxia provides a nonmutational explanation for its overexpression in human cancer. Proc Natl Acad Sci USA. 2007, 104: 13092-13097. 10.1073/pnas.0702387104.
Niu G, Briggs J, Deng J, Ma Y, Lee H, Kortylewski M, Kujawski M, Kay H, Cress WD, Jove R, Yu H: Signal transducer and activator of transcription 3 is required for hypoxia-inducible factor-1alpha RNA expression in both tumor cells and tumor-associated myeloid cells. Mol Cancer Res. 2008, 6: 1099-1105. 10.1158/1541-7786.MCR-07-2177.
Yamamoto Y, Ibusuki M, Okumura Y, Kawasoe T, Kai K, Iyama K, Iwase H: Hypoxia-inducible factor 1alpha is closely linked to an aggressive phenotype in breast cancer. Breast Cancer Res Treat. 2008, 110: 465-475. 10.1007/s10549-007-9742-1.
Tan EY, Yan M, Campo L, Han C, Takano E, Turley H, Candiloro I, Pezzella F, Gatter KC, Millar EK, O'Toole SA, McNeil CM, Crea P, Segara D, Sutherland RL, Harris AL, Fox SB: The key hypoxia regulated gene CAIX is upregulated in basal-like breast tumours and is associated with resistance to chemotherapy. Br J Cancer. 2009, 100: 405-411. 10.1038/sj.bjc.6604844.
Nakatani K, Thompson DA, Barthel A, Sakaue H, Liu W, Weigel RJ, Roth RA: Up-regulation of Akt3 in estrogen receptor-deficient breast cancers and androgen-independent prostate cancer lines. J Biol Chem. 1999, 274: 21528-21532. 10.1074/jbc.274.31.21528.
Charafe-Jauffret E, Ginestier C, Monville F, Finetti P, Adelaide J, Cervera N, Fekairi S, Xerri L, Jacquemier J, Birnbaum D, Bertucci F: Gene expression profiling of breast cell lines identifies potential new basal markers. Oncogene. 2006, 25: 2273-2284. 10.1038/sj.onc.1209254.
Dairkee SH, Puett L, Hackett AJ: Expression of basal and luminal epithelium-specific keratins in normal, benign, and malignant breast tissue. J Natl Cancer Inst. 1988, 80: 691-695. 10.1093/jnci/80.9.691.
Berns K, Horlings HM, Hennessy BT, Madiredjo M, Hijmans EM, Beelen K, Linn SC, Gonzalez-Angulo AM, Stemke-Hale K, Hauptmann M, Beijersbergen RL, Mills GB, van de Vijver MJ, Bernards R: A functional genetic approach identifies the PI3K pathway as a major determinant of trastuzumab resistance in breast cancer. Cancer Cell. 2007, 12: 395-402. 10.1016/j.ccr.2007.08.030.
Livasy CA, Karaca G, Nanda R, Tretiakova MS, Olopade OI, Moore DT, Perou CM: Phenotypic evaluation of the basal-like subtype of invasive breast carcinoma. Mod Pathol. 2006, 19: 264-271. 10.1038/modpathol.3800528.
Asselin-Labat ML, Sutherland KD, Barker H, Thomas R, Shackleton M, Forrest NC, Hartley L, Robb L, Grosveld FG, van der Wees J, Lindeman GJ, Visvader JE: Gata-3 is an essential regulator of mammary-gland morphogenesis and luminal-cell differentiation. Nat Cell Biol. 2007, 9: 201-209. 10.1038/ncb1530.
Kouros-Mehr H, Slorach EM, Sternlicht MD, Werb Z: GATA-3 maintains the differentiation of the luminal cell fate in the mammary gland. Cell. 2006, 127: 1041-1055. 10.1016/j.cell.2006.09.048.
Sircoulomb F, Bekhouche I, Finetti P, Adelaide J, Ben Hamida A, Bonansea J, Raynaud S, Innocenti C, Charafe-Jauffret E, Tarpin C, Ben Ayed F, Viens P, Jacquemier J, Bertucci F, Birnbaum D, Chaffanet M: Genome profiling of ERBB2-amplified breast cancers. BMC Cancer. 2010, 10: 539-10.1186/1471-2407-10-539.
Kao J, Salari K, Bocanegra M, Choi YL, Girard L, Gandhi J, Kwei KA, Hernandez-Boussard T, Wang P, Gazdar AF, Minna JD, Pollack JR: Molecular profiling of breast cancer cell lines defines relevant tumor models and provides a resource for cancer gene discovery. PLoS One. 2009, 4: e6146-10.1371/journal.pone.0006146.
Robey IF, Lien AD, Welsh SJ, Baggett BK, Gillies RJ: Hypoxia-inducible factor-1alpha and the glycolytic phenotype in tumors. Neoplasia. 2005, 7: 324-330. 10.1593/neo.04430.
Lim E, Vaillant F, Wu D, Forrest NC, Pal B, Hart AH, Asselin-Labat ML, Gyorki DE, Ward T, Partanen A, Feleppa F, Huschtscha LI, Thorne HJ, Fox SB, Yan M, French JD, Brown MA, Smyth GK, Visvader JE, Lindeman GJ: Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers. Nat Med. 2009, 15: 907-913. 10.1038/nm.2000.
Takubo K, Goda N, Yamada W, Iriuchishima H, Ikeda E, Kubota Y, Shima H, Johnson RS, Hirao A, Suematsu M, Suda T: Regulation of the HIF-1alpha level is essential for hematopoietic stem cells. Cell Stem Cell. 2010, 7: 391-402. 10.1016/j.stem.2010.06.020.
Simsek T, Kocabas F, Zheng J, Deberardinis RJ, Mahmoud AI, Olson EN, Schneider JW, Zhang CC, Sadek HA: The distinct metabolic profile of hematopoietic stem cells reflects their location in a hypoxic niche. Cell Stem Cell. 2010, 7: 380-390. 10.1016/j.stem.2010.07.011.
Caniggia I, Mostachfi H, Winter J, Gassmann M, Lye SJ, Kuliszewski M, Post M: Hypoxia-inducible factor-1 mediates the biological effects of oxygen on human trophoblast differentiation through TGFbeta(3). J Clin Invest. 2000, 105: 577-587. 10.1172/JCI8316.
Ezashi T, Das P, Roberts RM: Low O2 tensions and the prevention of differentiation of hES cells. Proc Natl Acad Sci USA. 2005, 102: 4783-4788. 10.1073/pnas.0501283102.
Helczynska K, Kronblad A, Jogi A, Nilsson E, Beckman S, Landberg G, Pahlman S: Hypoxia promotes a dedifferentiated phenotype in ductal breast carcinoma in situ. Cancer Res. 2003, 63: 1441-1444.
Nagata Y, Lan KH, Zhou X, Tan M, Esteva FJ, Sahin AA, Klos KS, Li P, Monia BP, Nguyen NT, Hortobagyi GN, Hung MC, Yu D: PTEN activation contributes to tumor inhibition by trastuzumab, and loss of PTEN predicts trastuzumab resistance in patients. Cancer Cell. 2004, 6: 117-127. 10.1016/j.ccr.2004.06.022.
Engelman JA, Zejnullahu K, Mitsudomi T, Song Y, Hyland C, Park JO, Lindeman N, Gale CM, Zhao X, Christensen J, Kosaka T, Holmes AJ, Rogers AM, Cappuzzo F, Mok T, Lee C, Johnson BE, Cantley LC, Janne PA: MET amplification leads to gefitinib resistance in lung cancer by activating ERBB3 signaling. Science. 2007, 316: 1039-1043. 10.1126/science.1141478.
We recognize research support from the NIH (NCI R01CA125618 to J.T.C., CA084955 to J.R.M.) and Komen Foundation grant KG090869 to JTC and JRM. MLG is supported by post-doctoral fellowship NIH-CA139890-01. We also appreciate the technical assistance of Ms. Gudrun Huper and helpful discussions with members of Chi and Nevins labs. The funding sources have no role in the study design; in the collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.
The authors declare that they have no competing interests.
MLG, HNK, JRM and JTC conceived and designed the study. MLG and HNK generated and analyzed data. KLB, MWD, and JRM provided unpublished reagents and data. MLG, HNK and JTC wrote the manuscript. All authors read and approved the final manuscript.
Michael L Gatza, Hsiu-Ni Kung contributed equally to this work.
Electronic supplementary material
Additional File 1: Supplementary methods. This document provides additional information about the methods used in the manuscript. (PDF 1 MB)
Additional File 2: Supplemental Table S1. Summary of microenvironment signatures. This document provides a summary of the microenvironment signatures and parameters. (PDF 192 KB)
Additional File 3: Supplemental Table S2. Summary of tumor samples. This document provides a summary of the tumor samples used in the study. (XLSX 34 KB)
Additional File 4: Supplementary Figures S1 to S18. This document contains Supplementary Figures S1 to S18. (XLSX 61 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Gatza, M.L., Kung, HN., Blackwell, K.L. et al. Analysis of tumor environmental response and oncogenic pathway activation identifies distinct basal and luminal features in HER2-related breast tumor subtypes. Breast Cancer Res 13, R62 (2011). https://doi.org/10.1186/bcr2899
- Validation Dataset
- Hypoxia Response
- Luminal Progenitor Cell
- Microenvironmental Stress
- Congruence Rate