Low ATM protein expression in malignant tumor as well as cancer-associated stroma are independent prognostic factors in a retrospective study of early-stage hormone-negative breast cancer

Introduction The serine/threonine protein kinase ataxia telangiectasia mutated (ATM) is critical in maintaining genomic integrity. Upon DNA double-strand breaks, ATM phosphorylates key downstream proteins including p53 and BRCA1/2, thereby orchestrating complex signaling pathways involved in cell cycle arrest, DNA repair, senescence and apoptosis. Although sporadic mutation of ATM occurs rarely in breast cancer, the status of its protein expression and its clinical significance in breast cancer remain not well established. Our study was designed to investigate the influence of ATM protein in both tumor and cancer-associated stroma on clinical outcome in hormone-positive (HPBC) and hormone-negative (HNBC) early-stage breast cancer (EBC). Methods Tissue microarrays (TMAs), containing formalin-fixed, paraffin-embedded resected tumors from two cohorts of patients (HPBC cohort: n = 130; HNBC cohort: n = 168) diagnosed at the Tom Baker Cancer Centre, Calgary, Canada, were analyzed for ATM protein expression using fluorescence immunohistochemistry (IHC) and automated quantitative analysis (AQUA). ATM expression levels were measured within the tumor as a whole (tATM) as indicated by pan-cytokeratin expression, tumor nuclear compartment (nATM) as indicated by both DAPI and pan-cytokeratin-positive results, and cancer-associated stroma (csATM) as indicated by vimentin-positive and pan-cytokeratin-negative results. ATM expression levels within these compartments were correlated with clinical outcome. Results While tATM and nATM were significantly lower in tumors compared to normal breast epithelial tissues, csATM was significantly higher than the corresponding normal tissue compartment. In addition, the median expression level of both tATM and nATM were two- to threefold lower (P <0.001) in HNBC than in HPBC. In both HNBC and HPBC cohorts, patients with low tATM, nATM and csATM tumors had significantly poorer survival outcomes than those with a high tATM, nATM and csATM, but this effect was more pronounced in HNBC. A multivariate analysis demonstrates that these biomarkers predict survival independent of tumor size and lymph node status, but only in the HNBC cohort (P <0.001). Conclusions Low ATM protein expression in both malignant tumor and stromal compartments likely contributes to the aggressive nature of breast cancer and is an independent prognostic factor associated with worse survival in HNBC patients. Electronic supplementary material The online version of this article (doi:10.1186/s13058-015-0575-2) contains supplementary material, which is available to authorized users.


Introduction
Breast cancer continues to be the most frequent cancer type and the second leading cause of cancer death in females in Western Europe and North America. Although we have significantly improved diagnostic tools and therapeutic interventions in breast cancer over the last decade, breast oncologists still face paramount challenges in hormone-negative breast cancer, especially triple-negative breast cancer, the most aggressive breast cancer subtype, for which we do not yet have any molecular-targeted treatments apart from cytotoxic chemotherapy. Several large studies have identified RNA expression signatures that are prognostic for disease recurrence and metastasis early-stage hormone-positive breast cancer (HPBC) [1][2][3]. No such biomarkers exist currently for the more aggressive hormone-negative breast cancer (HNBC). In addition, mounting evidence indicates that levels of proteins that are immediately relevant to cell growth and metabolism are often not particularly well correlated to mRNA levels [4], suggesting that it may be more informative to directly measure protein levels as a indicator of biological behavior. These biomarkers likely reflect genetic signatures of early cancer formation and evolution of cancer aggressiveness.
Ataxia telangiectasia mutated (ATM) belongs to a family of so-called phosphatidyl inositol-3 kinase (PI3K)-like serine/threonine protein kinases (PIKKs), most of which are involved in the cellular response to various stresses [5]. In response to DNA double-strand breaks, ATM phosphorylates multiple key downstream proteins including p53 and BRCA1/2, the most common sporadic and inheritable mutated genes in breast cancer [6], thereby orchestrating complex signaling pathways involved in cell cycle arrest, DNA repair, senescence and apoptosis [7]. Germline mutation of ATM genes causes ataxia telangiectasia, a rare pleiotropic autosomal recessive disorder characterized by heightened radiation sensitivity and a predisposition to cancers including breast cancer [6,8]. The incidence of sporadic ATM mutations in breast cancer is low (2 to 3%), as compared to p53 and PIK3A mutations whose incidence rates are 30 to 40% [9]. Although we have some evidence supporting ATM loss in malignant breast tumors [10][11][12][13][14][15][16][17], and its prognostic significance in breast cancer was recently uncovered [16,17], the protein level of ATM in breast cancer-associated stromal tissues remains unknown.
In this study, using fluorescent immunohistochemistry (IHC) and automated quantitative analysis (AQUA), we measured ATM protein expression in both malignant tumor and stromal compartments of HPBC and HNBC. We hypothesized that reduced ATM expression correlates with increased tumor aggressiveness and poorer clinical outcome. Remarkably, we also identified a strong association between stromal expression of ATM and clinical outcome, suggesting that the tumor microenvironment contributes significantly to the overall course of disease.

Case selection and clinical data collection
This study was approved by the University of Calgary Conjoint Faculties Research Ethics Board (Ethic ID #24525), in accordance with the Tri-Council Policy Statement on Research with Human Subjects. Consents directly from patients were not required in this study as per Ethics Board guidelines. Two cohorts of clinical data were collected retrospectively through chart review. One cohort consists of 130 hormone-positive resected earlystage breast cancer (EBC) patients diagnosed at the Tom Baker Cancer Centre (TBCC) from 1989 to 2000. The second cohort consists of 168 hormone-negative resected EBC patients diagnosed at the TBCC from 1999 to 2004. All patients had stage I to III breast cancer with no evidence of metastatic disease at diagnosis and had undergone either lumpectomy or mastectomy. Hormone status including estrogen receptor (ER) and progesterone receptor (PR) were initially assessed by local breast pathologists according to the American Society of Clinical Oncology and the College of American Pathologists guidelines at the time and further evaluated on tissue microarray (TMA) using an independent staining method. Samples with discordant results between initial pathological diagnosis and further evaluation on TMA were excluded.

Tissue microarray generation
All archived formalin-fixed paraffin-embedded (FFPE) resected EBC tumor samples were retrieved from Calgary Laboratory Services and reviewed by a pathologist (AMM). TMAs were constructed as previously described [18], consisting of two cohorts of patients (cohort 1: n = 130 HPBC; cohort 2: n = 168 HNBC) as described in Case Selection above. Normal breast epithelial tissue specimens (n = 8), normal tonsil tissue, an ATMdeficient human lymphoblastoid cell line (L3) and an age-matched ATM-expressing human lymphoblastoid cell line (BT/C3ABR) were also included as controls [19].

Four-color fluorescence immunohistochemistry
TMA sections (4 μm) were deparaffinized in xylene, rinsed in ethanol, and rehydrated as previously described [18,20]. Heat-induced epitope retrieval was performed by heating slides to 121°C for 6 minutes in a Tris/ EDTA-based buffer (pH 9.0) (Target Retrieval Solution, DAKO, Mississauga, ON, Canada). IHC was performed using a DAKO Autostainer Link 48. Primary antibodies used were a rabbit monoclonal antibody to ATM (clone Y170, 1:1000, Abcam, Cambridge, MA, USA), a mouse monoclonal antibody to vimentin, which marks the stromal compartment (clone V9, 1:5000, DAKO), and a guinea pig monoclonal pan-cytokeratin antibody identifying the tumor compartment (1:100, Acris, San Diego, CA, USA). Briefly, endogenous peroxidase activity was quenched with a 10-minute incubation of peroxidase-blocking reagent (a component of the DAKO EnVision™ + System) followed by a 15-minute protein block (Signal Stain; Cell Signaling, Danvers, MA, USA) to eliminate nonspecific antibody binding. All incubations were performed at room temperature. Slides were washed with a wash buffer containing Tris-buffered saline and 0.05% Tween-20 (TBST, DAKO) and then incubated with primary antibodies for 30 minutes. After three washes, a goat anti-rabbit horseradish peroxidase (HRP)-conjugated secondary antibody from the DAKO EnVision™ + System was applied for 30 minutes. Slides were again washed in TBST and treated for 5 minutes with TSA-Plus Cy3 tyramide signal amplification reagent (1:100, PerkinElmer, Woodbridge, ON, Canada). To stop enzyme activity in preparation for vimentin staining, slides were treated with peroxidase block for 30 minutes, followed by four washes of TBST, a 10-minute protein block, and three more washes. To identify stromal cells, slides were incubated for 30 minutes with the vimentin antibody, washed in TBST, and then treated for 30 minutes with a goat anti-mouse HRP-conjugated secondary antibody from the DAKO EnVision™ + System. After three TBST washes, slides were incubated for 5 minutes with TSA-Plus Cy5 tyramide signal amplification reagent (PerkinElmer). To identify tumor cells, slides were further incubated for 30 minutes with the pan-cytokeratin antibody, followed by three washes and a 30-minute incubation of a 1:200 dilution of Alexa-488 conjugated goat anti-guinea pig antibody (Life Technologies, Burlington, ON, Canada). After three washes in TBST, the TMA slides were mounted with ProLong™ Gold anti-fade mounting medium containing diamidino-2-phenylindole (DAPI) (Life Technologies) and stored at 4°C overnight to set before scanning.

Automated image acquisition and analysis
Automated image acquisition was performed using an Aperio Scanscope FL (Aperio Inc., Vista, CA, USA). Briefly, seamless high-resolution slide images were acquired using the Scanscope FL 10-bit monochrome TDI line-image capture camera using filters specific for DAPI to define the nuclear compartment, fluoroscein isothiocyanate (FITC) to define cytokeratin for the tumor cytosolic compartment, Cy3 to define the target marker ATM, and Cy5 to define vimentin-positive nonmalignant stromal cells. Images were then analyzed using the HistoRX AQUAnalysis™ program, version 2.4.4.1 as described [13]. Briefly, a tumor-specific binary mask was generated to distinguish the cancer cells from surrounding stromal tissue by thresholding the pan-cytokeratin images. Thresholding created a binary mask that identified the presence or absence of tumor cells by the presence of a pixel that was 'on' or 'off ' , respectively. A similar thresholding technique was applied to generate a nonmalignant stroma-specific mask (from the pancytokeratin-negative regions of the vimentin images, a proven superior method for identifying stromal regions [12]. To isolate tumor nuclei, mutually exclusive tumor nuclear and cyoplasmic masks were created by subtracting cytokeratin-positive areas from DAPI-positive areas within the tumor mask. All images were processed using optimized threshold values and all subsequent image manipulations involved only image information from the masked areas. Images were manually validated according to the following: 1) >10% of the tissue area is pancytokeratin positive, 2) >50% of the image was usable (that is not compromised due to overlapping or out-offocus tissue). Unusable areas within each image were manually cropped so that they were excluded from the final analysis.
ATM compartment-specific quantification ATM expression was calculated as the mean Cy3 pixel intensity within each defined compartment. The tumorspecific, tumor nuclear-specific and stroma-specific AQUA scores represent ATM protein expression in the total malignant area, in nuclei within the malignant area, and in the tumor-associated stromal area, respectively. For each patient sample, the average compartmentspecific AQUA score over triplicate cores was used to define the ATM expression tumor score (tATM), tumor nuclear score (nATM), and the cancer-associated stromal score (csATM).

Statistical analysis
Summary and descriptive statistics were used to describe demographic and clinicopathological characteristics of patients in the two cohorts. Student's t test was used to compare the continuous variable age between low and high ATM groups. Fisher's exact test was used to compare the categorical variables between low and high ATM groups. Histograms were used to evaluate the distribution of ATM expression (tATM, nATM, and csATM) of the two cohorts. Pairwise group comparisons of ATM expression (tATM, nATM, and csATM) between groups (HPBC, HNBC, human epidermal growth factor receptor 2 (HER2)-positive and triple-negative breast cancer (TNBC)) were performed using parametric (Student's t test) and nonparametric (Wilcoxon) tests. Cut-points to define low and high ATM expression were modeled using X-tile software [21] in the two cohorts independently. Disease-specific survival (DSS) and disease-free survival (DFS) estimates were calculated using the Kaplan-Meier (KM) survival analysis and logrank tests were used to compare between low and high ATM groups. Cox proportional hazards regression methods were used to assess the prognostic effect of ATM expression (tATM, nATM, and csATM) for DSS multivariably. A multivariate model from selected clinical known prognostic factors such as tumor size, grade, lymph node (LN) status, lymphovascular invasion (LVI) and age was created. The prognostic ability of ATM expression (tATM, nATM, and csATM) was tested using a backward selection approach to obtain the optimal model. All tests were two-sided, and a P value of 0.05 was considered statistically significant. All statistical analyses were conducted using R statistics software (Version 3.1.2).

Patient characteristics
Our study adheres to the REMARK criteria for the study of biomarkers [22]. The clinical and tumor characteristics of patients in the two cohorts (HPBC and HNBC) are shown in Table 1. Briefly, in the HNBC cohort, mean age at diagnosis was 53 years. The majority of cases (72%) were grade 3, consistent with the aggressive biology of hormone-negative tumors. More than half (58%) the cases were LN negative, and a majority of cases (81%) were of stages I or II, indicative of the early staging of this cohort. Two-thirds of the cases were TNBC, and the remaining third were hormone-negative, HER2positive cases. Eleven cases in the HNBC cohort were deemed hormone positive based on analysis of the TMA spots, and were excluded from final statistical analysis. Sixty percent of patients had adjuvant radiation treatment (RT). Seventy-five percent of patients had adjuvant chemotherapy; the majority of these patients (83%) had anthracycline-based chemotherapy, while the remainder (17%) received taxane-based chemotherapy. In the HPBC cohort, the mean age at diagnosis was 67 years, older than our HNBC cohort that is consistent with clinical observation. Both cohorts shared similar patterns of tumor size, LVI, stage, LN status and adjuvant RT. Almost all patients (99%) in the HPBC cohort received adjuvant hormone treatment instead of chemotherapy.

ATM expression analysis
ATM protein expression levels were measured by automated quantitative fluorescence IHC using AQUAnalysis software [23] (Figure 1). Antibody specificity was validated using the L3 ATM-deficient human lymphoblastoid cell line and the BT/C3ABR age-matched ATMexpressing human lymphoblastoid cell line as controls ( Figure 1A, left two panels) [19]. Positive and negative (by omission of the primary antibody) tonsil tissue controls were included in each run ( Figure 1A, right two panels). Representative examples of the ATM-staining patterns in normal breast epithelium and invasive ductal carcinomas representing preserved or aberrantly reduced ATM expression are shown in Figure 1B. Normal breast epithelial tissue (pan-cytokeratin-positive) and hormonepositive tumors with high AQUA scores demonstrated a relatively uniform and intense nuclear expression pattern of ATM, consistent with the reported cellular distribution of ATM [7] ( Figure 1B, top panel). In contrast, hormone-negative tumors with low AQUA scores exhibited only faint and/or sporadic nuclear staining pattern in malignant cells ( Figure 1B, bottom panel).
To determine whether ATM expression was altered in EBC compared to normal breast tissue, we compared the median ATM AQUA scores of our two EBC cohorts to the median scores from normal breast tissue ( Figures 1C and D). tATM was calculated within each patient tissue core as the average Cy3 pixel intensity within the pan-cytokeratin-positive malignant cell area ( Figure 1C). nATM was calculated within each patient tissue core as the average Cy3 pixel intensity within both the pan-cytokeratin and DAPI-positive malignant cell area ( Figure 1D). Median tATM and nATM scores in normal breast epithelium were 148 (95% confidence interval (CI): 84 to 19) and 294 (95% CI: 208 to 548) respectively ( Figure 1C Figure 1D, hashed red line) respectively. The median tATM and nATM scores were two-to threefold lower in the HNBC cohort than in the HPBC cohort (P <0.001) (Figure 2A and B). We further classified our HNBC cohort into TNBC and hormone-negative, HER2-positive groups based on HER2 status. We found a persistently significant difference in tATM and nATM scores between the HPBC and the TNBC cohorts, and between the HPBC and the hormone-negative HER2positive groups (the median tATM scores in TNBC and hormone-negative HER2-positive groups were 52 (95% CI: 37.4 to 69.2) and 85 (95% CI: 44.2 to 117.7) respectively; the median nATM scores in TNBC and hormone-negative HER2-positive groups were 94 (95% CI: 66.9 to 142.2) and 187 (95% CI: 80.7 to 299.0) respectively; Figure 2C and D). It seemed that the TNBC group had the lowest tATM and nATM expression numerically; statistically, the different between the TNBC and the hormonenegative HER2-positive group was not significant, likely due to the small sample sizes of both groups ( Figure 2C and D). We further asked whether there were differences in ATM protein between the cancer-associated stromal compartment adjacent to the malignant tumor compared to the stromal compartment around normal breast epithelium. The csATM score was calculated within each patient tissue core as the average Cy3 pixel intensity of the vimentin-positive but pan-cytokeratin-negative stromal cell area. While median stromal ATM score adjacent to normal breast epithelium was 45 (95% CI: 32 to 72) ( Figure 3A, hashed black line), median csATM score in the HPBC cohort was 85 (95% CI: 79 to 95) ( Figure 3A, hashed green line) and median csATM was in the HNBC cohort was 62 (95% CI: 51 to 81) ( Figure 3A, hashed red line). Interestingly, we found that csATM expression was higher in significant portions of patients in both of our two EBC cohorts, more pronounced in the HPBC than in the HNBC cohort, compared to normal stromal ATM expression surrounding normal breast epithelium ( Figure 3B). Our result is consistent with a recent study demonstrating that an abnormally activated ATM-associated pathway had been identified in breast cancer-associated fibroblasts (CAFs) [24], although its biological significance remains elusive. We were not able to identify the biological significance of the increase of csATM compared to normal stroma ATM expression in our study as high csATM expression did not associate with tumor aggressiveness. Rather, low csATM scores appear to be associated with tumor-aggressive features (see below).

Association between ATM protein expression and clinical outcome
Patients were dichotomized within each cohort into those with high and low tATM, nATM and csATM scores, using cut-points identified by X-Tile software [21]. In the HPBC cohort, low tATM and csATM scores correlated with larger tumor size (P = 0.002 and only with a trend P = 0.08 respectively), LN positivity (P = 0.03 and P = 0.004 respectively), higher grade (P = 0.029 and P = 0.033 respectively), and higher stage (P = 0.002 and P = 0.014 respectively), suggesting that low tATM and csATM scores are associated with more aggressive breast tumors in HPBC (Table 1 and Table S1 in Additional file 1). Interestingly, these correlations were not apparent in the HNBC cohort.
Patients were dichotomized into low and high tATM, nATM and csATM scores as described above for further survival analysis using the KM method. As a whole (the HPBC plus the HNBC cohort), patients with low tATM, nATM and csATM scores, compared to patients with high tATM, nATM and csATM scores, had a worse DSS ( Figure S2A-C in Additional file 2). Then, patients with HNBC and HPBC were examined separately. In both cohorts, compared to patients with a high tATM, nATM and csATM expression, patients with low tATM, nATM and csATM scores had worse DSS outcome with median follow-up of 7.5 and 15 years in the HNBC and HPBC cohorts respectively (Figures 4A-F). This effect seems greater in the HNBC cohort compared with the HPBC cohort ( Figures 4A-F).
Cox proportional hazards regression models were used to determine if reduced ATM expression identified by tATM, nATM and csATM scores was an independent prognostic factor for DSS. A multivariate model from selected known clinical prognostic factors such as tumor size, grade, LN status, LVI and age was created. As expected, tumor size and LN status were significantly associated with survival. In the HNBC cohort, the significant association between these biomarkers and DSS persisted in the multivariate model after adjusting for all these known prognostic variables in the optimal model (P <0.001), but not in the HPBC cohort (Table 2A-C and  data not shown). Furthermore, within the HNBC cohort, the prognostic value of ATM expression appeared to be similar in the LN-negative and LN-positive subsets ( Figure S3A and B in Additional file 3). However, it seemed to have a greater effect in the T1/T2 subset compared to the T3/T4 subset ( Figure S3C and D in Additional file 3), although the latter result may be confounded by a very small sample size in the T3/T4 subgroup. One hundred and sixty-eight HNBC and 130 HPBC patients diagnosed with EBC were included in the study. Clinical and histopathological features were summarized and described in numeric numbers and proportions in both cohorts. Sixty-seven cases in the HNBC cohort and 28 cases in the HPBC cohort were excluded for failing to obtain ATM staining due to technical issues. Therefore, a total of 101 cases in the HNBC cohort and 102 cases in the HPBC cohort were dichotomized into high and low ATM groups based on tATM AQUA scores. The cut-point was set at 22% in the HNBC cohort and 38% in the HPBC cohort respectively based on X-tile prediction. Student's t test was used to compare the continuous variable age between low and high ATM groups. Fisher's exact test was used to compare the clinicopathological features of breast tumors in high and low ATM groups. HPBC, hormone-positive breast cancer; HNBC, hormone-negative breast cancer; tATM, tumor ATM; LVI, lymphvascular invasion; LN, lymph node status; DCIS, ductal carcinoma in situ; ER, estrogen receptor; PR, progesterone receptor; HER2, human epidermal growth factor receptor 2; RT, radiation treatment; CMT, combined modality treatment; 5FU, 5-fluorouracil; EBC, early-stage breast cancer; ATM, ataxia telangiectasia mutated; AQUA, automated quantitative immunofluorescence analysis.

Discussion
In this study, we show that ATM expression was reduced in HNBC, as compared to either HPBC or normal breast epithelium. In both HPBC and HNBC cohorts, reduced ATM expression was associated with poorer disease-specific outcomes. Interestingly, stromal ATM expression was higher in both cohorts as compared to normal epithelium, though lower expression in this compartment was also associated with poorer clinical outcome. These expression studies were enhanced by fluorescent IHC and AQUAnalysis, which allowed us to describe ATM expression as a continuous variable across a broader range than is possible with brightfield IHC. The dynamic range of protein quantification that is possible in the fluorescent realm allows the analysis, in a research context, of samples with widely varying degrees of protein expression. It should be noted that fluorescent IHC and the associated automated image analysis workflows are not currently used in clinical practice; therefore some data presented in this study may not be directly applicable. It is however possible that they could be translated into clinically relevant assays upon identification of the levels of ATM expression that offer prognostic value. To the best of our knowledge, this is the first study to examine ATM protein expression in malignant tumor and cancer-associated stromal tissues simultaneously. About 60 to 80% of the EBC patients in both cohorts had ATM scores (tATM or nATM) lower than the mean values for the same compartment in normal breast epithelium. This suggests that loss of ATM expression is an important early event signaling malignant transformation of breast epithelial cells as the majority of EBC has some degree of ATM deficiency. In addition, ATM loss continues when breast cancer evolves and becomes more aggressive, to a certain threshold that significantly impacts on cancer-specific survival. Reduced ATM expression within malignant breast epithelial tissues was predominant in HNBC, especially TNBC, the most aggressive breast cancer type. These results, along with similar outcomes reported by others [10][11][12][13][14][15][16][17], suggest that ATM loss not only occurs early during breast tumorigenesis, but also may promote cancer progression by permitting the accumulation of genetic mutations due to genome instability as a result of ATM loss, thereby selecting cancer cells bearing aggressive phenotypes.
In HPBC, low ATM within malignant breast epithelial cells and stroma is associated with aggressive characteristics of breast cancer including large size, high grade and LN involvement. It is an interesting observation that this association is not seen in HNBC, which is a more aggressive tumor type. It is likely that ATM loss predominates the aggressive nature regardless of other worse features in this aggressive tumor type, which is consistent with our later finding of independent prognostic ability of ATM in only the HNBC cohort. Previous studies [16,17] demonstrated the significant association between low ATM expression and worse clinicopathological features in an undifferentiated breast cancer cohort containing both HPBC and HNBC cases. We believe that our data, consistent with others [10][11][12][13][14][15][16][17], support the notion that ATM loss contributes to the aggressive nature of breast tumors as they evolve. It is interesting that ATM is higher in cancerassociated stromal tissues compared with stromal tissues around normal breast epithelium. Although this is consistent with another study [24], validation is needed using a larger distinct cohort. To the best of our knowledge, the mechanisms of regulation and functional roles of ATM in malignant stroma are presently unknown. Although it may seem contradictory that the ATM level is increased in cancer-associated stroma compared to normal stroma, albeit the ATM level is decreased in malignant epithelium compared to normal epithelium, we believe that this could be a reflection of the complexity of ATM regulation in diverse cell types. We speculate that the heightened level of DNA damage and oxidative stress [25] that are predominant in malignancy may initially lead to significant induction of ATM in the malignant epithelial tumor as well as in the stroma, Subsequently, a range of other factors may lead to differentially active downregulation of ATM, more so within the epithelial tumor and less so within the stroma. For example, one of such factors may include caspases that are activated particularly in the necrotic cores of aggressive tumors, which were recently shown to proteolytically cleave ATM [26]. Other regulatory factors are discussed later. Furthermore, we speculate that ATM may have context-dependent roles in stroma that are distinct from its well-known conventional function in maintaining genomic stability in response to DNA damage within the epithelial component. For example, it may regulate inflammatory response such as chemokine release (for example interleukin (IL)-6, IL-8) to promote/inhibit tumor growth in a paracrine-mediated fashion [27]. In addition, it may serve as a contributor to extracellular matrix stiffness (in addition to many other cytokines) and regulates tumor invasion and metastasis through bidirectional signaling between tumor cells and stroma microenvironment [28]. Therefore, it is not surprising that low level of ATM in cancer-associated stroma is also associated with a poor survival ( Figure 4E and F). It should be noted that cut-points used to define positive and negative ATM tumors varied widely in previous studies despite all having used conventional IHC as there is no well-defined method for choosing cut-points for biomarker analysis. For example, Bueno's group [16] defined ATM IHC scores as negative or positive according to the absence or presence of nuclear staining in epithelial cells evaluated by two independent pathologists who had been blinded to the outcome. Abdel-Fatah's group [17] used a cutoff of <25% cells being classed as low and ≥25% as high for nuclear ATM protein level. Because ATM expression in our study is a continuous variable that is generated automatically by the analysis software, we chose to use the X-tile program [21], which is a validated methodology to define our cut-point in two independent cohorts. The other advantage of using such digital image analysis as compared to conventional IHC is that it can potentially avoid human bias (inter-and intra-observer variability) as the intensity scores are objective representations of protein expression. In addition, this platform allows us to define the specific localization of ATM expression within the tumor. Therefore, we can evaluate tumor nuclear and cytoplasmic expression of ATM concurrently, in a similar fashion as the above tumor and stromal ATM analysis. We examined nuclear versus cytoplasmic ratio of ATM expression (ncATM) in both EBC cohorts as we speculated that aberrant localization of ATM may also contribute to loss of ATM function thereby resulting in tumor aggressiveness and poor prognosis. Although we did observe a significantly decreased ncATM in more aggressive breast tumor, HNBC compared with HPBC ( Figure S1 in Additional file 4), the prognostic value of ncATM for clinical outcome (DFS and DSS) was not as significant as tATM, nATM and csATM in both cohorts (data not shown). We speculate that this may be due to our small sample size.
Lastly, our study illustrates the independent prognostic ability of ATM protein expression levels for DSS in both malignant tumor and stromal tissues in HNBC. Bueno et al. found that ATM protein expression associated significantly with both DFS and DSS [16]. We also observed a strong trend that low ATM group (tATM, nATM and csATM) had worse DFS than high ATM (tATM, nATM and csATM) group with a P value of 0.05 to 0.09 based on log-rank tests of KM analysis (data not shown). The lack of statistical significance may be an outcome of relative small sample size as well as the different methodologies of evaluating ATM protein expression and defining cut-points to divide ATM low and high groups. We, and others [16,17], clearly demonstrated that ATM protein has potential use as a novel biomarker in EBC, in addition to the known clinicopathological prognostic factors such as tumor size, grade, LVI and LN status. This could be a much simpler tool than several biomarker tools deployed in HPBC such as Oncotype DX, PAM50 and MammaPrint, which are all RNA-based assays [1][2][3].
Our study did not seek to investigate the mechanisms underlying the regulation of ATM expression in EBC. Independent prognostic value of ATM expression levels (tATM, nATM and csATM) in both malignant tumor and stromal tissues was evaluated in a multivariate model. A multivariate model was created using selected known prognostic factors such as tumor size, lymph node (LN) status, lymphovascular invasion (LVI), grade and age. Cox proportional hazards regression methods were used to calculate P value, hazard ratio (HR) and confidence interval (CI) of each variable in the model. ATM, ataxia telangiectasia mutated; HNBC, hormone-negative breast cancer; tATM, tumor ATM; nATM, tumor nuclear ATM; csATM, cancer-associated stromal ATM.
ATM gene mutations have been described in the literature, but since the incidence of ATM mutations in breast cancer is low [9], we would not expect this to explain the full range of our results. Other mechanisms such as gene copy number loss [16,29], posttranscriptional microRNA regulation [16,30,31], epigenetic modifications [32,33] and direct caspase cleavage [26] for altering ATM levels have been described. It should also be noted that a simple measure of ATM protein expression by IHC cannot establish the functional capability of ATM within the tumor. Measuring of phosphorylated ATM protein and its targets within the ATM signaling pathways may provide more insight with regard to ATM functions in breast cancer. Our results taken in the context of what is known about the biological role of ATM in protecting the genome suggest that it is very likely that a substantial proportion of breast cancers, which have lost ATM expression at the protein level, are in fact deficient also in ATM function [34]. We attempted to investigate the predictive value of ATM protein for the efficacy of adjuvant RT and chemotherapy in our cohorts. Due to small sample size, we cannot draw any conclusions. We noted that other studies have showed that ATM loss seems to associate with resistance to anthracycline chemotherapy in breast cancer [16,17,29], which seems counterintuitive to many preclinical studies revealing that ATM-deficient cells are sensitive to cytotoxic chemotherapy [7]. More interestingly, we and others have recently shown that ATM loss can predict sensitivity to poly (ADP-ribose) polymerase (PARP) inhibitors, a novel drug that may exploit this deficiency, in breast cancer [35] as well as other cancers [19,36,37]. Prospective clinical trials are warranted to validate its predicative value.

Conclusions
In summary, we report that low ATM protein expression in malignant breast epithelial tissues can be detected in a significant proportion of EBC, and is more pronounced in HNBC compared to HPBC. Our analysis also demonstrates that low ATM protein expression can independently predict a poor DSS in early-stage HNBC. Furthermore, we show that ATM protein expression is increased in cancer-associated stroma compared to normal stroma, although its biological significance remains to be elucidated.