Gene expression profiling of primary male breast cancers reveals two unique subgroups and identifies N-acetyltransferase-1 (NAT1) as a novel prognostic biomarker
© Johansson et al.; licensee BioMed Central Ltd. 2012
Received: 4 August 2011
Accepted: 14 February 2012
Published: 14 February 2012
Male breast cancer (MBC) is a rare and inadequately characterized disease. The aim of the present study was to characterize MBC tumors transcriptionally, to classify them into comprehensive subgroups, and to compare them with female breast cancer (FBC).
A total of 66 clinicopathologically well-annotated fresh frozen MBC tumors were analyzed using Illumina Human HT-12 bead arrays, and a tissue microarray with 220 MBC tumors was constructed for validation using immunohistochemistry. Two external gene expression datasets were used for comparison purposes: 37 MBCs and 359 FBCs.
Using an unsupervised approach, we classified the MBC tumors into two subgroups, luminal M1 and luminal M2, respectively, with differences in tumor biological features and outcome, and which differed from the intrinsic subgroups described in FBC. The two subgroups were recapitulated in the external MBC dataset. Luminal M2 tumors were characterized by high expression of immune response genes and genes associated with estrogen receptor (ER) signaling. Luminal M1 tumors, on the other hand, despite being ER positive by immunohistochemistry showed a lower correlation to genes associated with ER signaling and displayed a more aggressive phenotype and worse prognosis. Validation of two of the most differentially expressed genes, class 1 human leukocyte antigen (HLA) and the metabolizing gene N-acetyltransferase-1 (NAT1), respectively, revealed significantly better survival associated with high expression of both markers (HLA, hazard ratio (HR) 3.6, P = 0.002; NAT1, HR 2.5, P = 0.033). Importantly, NAT1 remained significant in a multivariate analysis (HR 2.8, P = 0.040) and may thus be a novel prognostic marker in MBC.
We have detected two unique and stable subgroups of MBC with differences in tumor biological features and outcome. They differ from the widely acknowledged intrinsic subgroups of FBC. As such, they may constitute two novel subgroups of breast cancer, occurring exclusively in men, and which may consequently require novel treatment approaches. Finally, we identified NAT1 as a possible prognostic biomarker for MBC, as suggested by NAT1 positivity corresponding to better outcome.
Male breast cancer (MBC) is a rare cancer form accounting for only 0.6% of all breast cancer cases in the Nordic countries . MBC is similar to female breast cancer (FBC) in many ways and is often likened to post-menopausal breast cancer in women due to the high prevalence of estrogen receptor (ER) positivity and relatively high age at onset. There are nevertheless also distinct differences; there is an ongoing debate regarding the level of similarity between FBC and MBC, and whether MBC may be a unique tumor type with biological features and clinicopathological parameters distinct from FBC [2–4]. MBC tumors are more frequently hormone receptor positive than FBC tumors (ER positivity 91% vs. 76% and progesterone receptor (PR) positivity 81% vs. 67%, respectively). Human epidermal growth factor receptor 2 (HER2) over-expression and/or amplification appear less frequent in MBC and the mean age at diagnosis is approximately five years older than for women [2, 3, 5, 6]. Risk factors include hormonal imbalances (for example, caused by liver disease, Klinefelter's syndrome or obesity), genetic predisposition (mainly due to BRCA2 mutations) and environmental factors (for example, exposure to chronic heat or radiation) [7, 8]. Survival rates have been debated, with some studies finding that men diagnosed with breast cancer have a worse prognosis than women [9, 10], whereas other studies have reported similar prognoses [11, 12]. The rarity of the disease has, however, precluded randomized trials for optimizing patient management; thus, recommendations for treating MBC are extrapolated from small retrospective trials and prior knowledge of FBC .
Importantly, no major progress has been made in the treatment of MBC since the introduction of hormonal therapy; survival rates have not improved over the last decades, unlike the female counterpart. Refined, comprehensive classification and identification of novel biomarkers will greatly increase our understanding of the pathobiology of the disease, and enable personalized clinical management as well as rationales for targeted therapy. We have previously described two genomic subgroups of MBC by array-based comparative genomic hybridization (aCGH) , and a few smaller studies have described specific differences between MBC and FBC based on gene expression (GEX) , microRNA [16, 17] and genomic profiles [14, 18], respectively.
In the present study, we aimed to understand MBC on the transcriptional level and to subclassify tumors into comprehensive subgroups. We also wanted to further validate the previously identified genomic subgroups that were based on the same cases , and to compare MBC with FBC. To this end, molecular profiling has been extensively applied to FBC by numerous independent researchers, resulting in the subdivision into gene expression-based 'intrinsic' subgroups associated with differences in survival as well as biological phenotypes [19–22]. Herein, we describe two stable subgroups of MBC, luminal M1 and luminal M2, respectively, highly correlated to the recently described MBC genomic subgroups . Remarkably, these subgroups were distinct from the well-established intrinsic subgroups of FBC, and may, as such, represent unique subtypes of breast cancer arising exclusively in males. The largest subgroup (luminal M1), comprising two-thirds of all cases, displayed a more aggressive phenotype and worse prognosis compared to the other cases, while high expression of immune response and ER-related genes was seen in the smaller subgroup (luminal M2). Finally, we identified N-acetyltransferase-1 (NAT1) as a potential prognostic biomarker in MBC.
Materials and methods
Clinicopathological data for the fresh frozen and paraffin-embedded MBC tumors, respectively
Fresh frozen tumors N (%)
Paraffin-embedded tumors N (%)
Age at diagnosis
42 to 93
23 to 98
BRCA2 mutation status
Invasive cancer in combination with DCIS
Follow-up time (years)
0.20 to 15
0.04 to 15
Adjuvant endocrine therapy
Gene expression (GEX) analysis
Tumor cellularity was determined on H&E stained sections and only tumors with high (> 70%) tumor cell content were included. Total RNA was extracted from fresh frozen tumors using the RNeasy Lipid Tissue Mini Kit (QIAGEN, Valencia, CA, USA) and RNA integrity was assessed on an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). RNA quantification was performed using a NanoDrop ND-1000 (NanoDrop Products, Wilmington, DE, USA). Sixty-six samples with RIN values ≥ 7 were hybridized to Human HT-12 v3.0 Expression BeadChips (Illumina Inc, San Diego, CA, USA) in three batches at the SCIBLU Genomics Center at Lund University. Data normalization and management were performed using BioArray Software Environment (BASE)  and R . Data were normalized using quantile normalization in BASE and were thereafter log2 transformed. To handle potential platform related biases, four samples each from hybridization batches one and two were re-hybridized in the third batch, resulting in a total of 74 experiments. A principal component analysis (PCA) was run and associations between principal components and technical and biological annotations were evaluated, whereupon a batch effect was detected as the main principal component (Additional file 2A). To correct for technical biases a supervised empirical Bayes method (ComBat) was thus applied . PCA was then carried out on ComBat corrected data, whereupon no technical variation was found among the main principal components (Additional file 2B). The gene expression data have been published in NCBI's Gene Expression Omnibus (GEO) database (GSE31259) .
Unsupervised discovery of MBC GEX subgroups
Probes with low signals (mean < 5.8 across all experiments) were filtered away and probes that varied the most across experiments were selected for use in unsupervised hierarchical clustering (HCL). Probes were mean centered across experiments. Pearson correlation distance and complete linkage were used for HCL. To assess the robustness of the initial HCL analysis, a multiscale bootstrap resampling was performed on the probes using the R-package Pvclust . To further validate the stability of the clusters, resampling was performed on the samples to obtain 10,000 bootstrapped datasets. The co-clustering frequencies of sample pairs across the datasets were calculated. HCL was then carried out on the co-clustering frequencies and the dendrogram clusters were compared with the clusters from the initial HCL analysis as described . When stable clusters were detected the procedure was repeated for each subcluster until no more stable clusters could be detected.
The Illumina probes were re-annotated using Re-annotation and Mapping for Oligonucleotide Array Technologies (ReMOAT) , and for the ontology studies only probes with good or perfect quality were used. A two-class unpaired significance analysis of microarray (SAM) was performed for MBC subgroups to identify differentially expressed genes, and the false discovery rate (FDR) 0 was used as a cut-off for significance. Up- and down-regulated genes were run separately in the database for annotation, visualization and integrated discovery (DAVID) v6.7 to identify possible enrichment of genes with specific biological themes separating the subgroups [30, 31].
Module signatures from FBC
where x i is the expression of gene i in the module and w i is either +1 or -1 depending on the up- or down-regulation of each gene in the original FBC study . The module scores were also calculated for a reference dataset representing all subgroups of FBC  as well as for an external dataset of MBC .
Independent validation and comparison with FBC using external datasets
An external GEX dataset on custom made cDNA microarray slides containing 16,457 sequence-verified I.M.A.G.E. clones (Research Genetics, Invitrogen, Carlsbad, CA, USA) with 37 MBCs was downloaded from ArrayExpress (ID: E-TABM-810) . Data quality assessment and normalization were performed in R . Limma packages were used for background correction; a within-array method for data centering followed by a between-array method using quantile normalization. A normalized FBC dataset containing 359 FBC samples representing all intrinsic subgroups was downloaded from GEO (GSE22133) [22, 26]. Samples in the validation sets were classified into MBC subgroups using nearest centroid classification. Centroids were calculated in our MBC data using the top 124 genes from the SAM analysis. Samples were classified based on to which centroid they showed the highest correlation and were unclassified if the correlation was < 0.2. We also performed nearest centroid classification for the intrinsic subgroups of FBC using the genes from Hu et al. on the MBC samples in our study as well as the MBC validation samples . The Hu classifier relies on expression levels relative to the spectrum of FBC, that is, including both ER negative (ER-) and positive (ER+) samples. Since the vast majority of MBCs are ER+, we constructed an ER+ specific FBC subtype classifier. Briefly, ER+ samples were extracted from the FBC reference dataset and genes were mean centered across these samples. A SAM analysis was performed between the luminal A and B samples among these ER+ FBC samples, whereupon 300 significant genes with FDR = 0 were selected. Centroids were calculated for luminal A and B tumors separately using these 300 genes. Both MBC datasets were classified using these ER+ FBC luminal centroids. Finally, for the MBC validation dataset, the detection of stable subgroups was performed using the same unsupervised approach used for our dataset, and these subgroups were then compared with the subgroups from the centroid classification.
Validation immunohistochemistry (IHC)
A tissue microarray (TMA) with two 1 mm cores from each of 220 MBC tumors was constructed as described . Sections of 3 to 4 μm were cut, transferred to glass slides, dried at room temperature and then baked in a heat chamber for two hours at 60°C. The DAKO Envision horseradish peroxidase rabbit/mouse kit (DAKO, Glostrup, Denmark) and a Dakocytomation Autostainer (DAKO) were used for the staining procedure. A monoclonal mouse antibody to the polymorphic heavy chain of human MHC Class I (HC10, diluted in 1:1,000 in high pH), with preferential binding to HLA-B and HLA-C alleles and some HLA-A was generously provided by Prof. Dr. J. Neefjes [33, 34], and the primary NAT1 antibody (diluted 1:1,000 in low pH) has been previously described [35, 36]. The evaluation of NAT1 and HLA was performed by one reader (IJ) in a blinded manner. The intensity of the staining in the tumor cells was scored on a scale as: 0 (absent), 1 (weak), 2 (moderate) or 3 (strong). The percentage of positively stained tumor cells was scored as: 0 (< 5%), 1 (5 to 25%), 2 (26 to 50%), 3 (51 to 75%) or 4 (> 75%).
All figures and statistical calculations were generated in R . For the survival analyses the survival and survcomp packages were used with distant metastasis free survival (DMFS) as end-point. All P-values are two-sided.
Discovery of two stable subgroups of MBC
Differences in key biological processes between MBC subgroups
Gene ontology indicates that luminal M1 tumors are more aggressive than luminal M2 tumors
A SAM analysis was performed, resulting in 544 up-regulated genes and 370 down-regulated genes in the luminal M2 compared to the luminal M1 subgroup (FDR = 0). The up- and down-regulated genes were uploaded separately into DAVID and some of the most relevant gene ontology (GO) terms associated with genes up-regulated in luminal M1 tumors included: cell migration, cell adhesion, angiogenesis, cell cycle, cell division and HOX genes (Additional file 3). Luminal M2 tumors, on the other hand, displayed up-regulation of genes associated with the GO term class I histocompatibility antigen, which is involved in the immune system (Additional file 4). Taken together, this suggests that luminal M1 tumors may be more aggressive than luminal M2 tumors, and may thus be associated with inferior outcome.
Validation of the MBC subgroups with an independent external dataset
The MBC subgroups differ from the intrinsic subgroups of FBC
In an effort to identify the degree of similarity between MBCs and the commonly used intrinsic subtypes of FBC [19, 21], we applied a centroid-based approach based on the Hu genes to our MBC dataset. This classification left 55% of the samples unclassified. Taking into account that MBCs are generally ER+, we also classified the MBCs using ER+ FBC luminal subgroup centroids. Even using this approach, 36% of the samples remained unclassified. Conversely, when ER+ FBC samples were classified using the MBC subgroup centroids, 151 samples (63%) were unclassified. Interestingly, luminal M1 tumors showed significantly higher correlations to the HER2 and basal centroids (Wilcoxon test, P = 1.1 × 10-3 and P = 6.8 × 10-5, respectively), while luminal M2 tumors were significantly higher correlated to both luminal A and B centroids (Wilcoxon test, P = 1.9 × 10-5 and P = 0.018, respectively). Furthermore, pronounced differences in the GEX patterns were observed when the Hu genes and the ER+ FBC luminal subgroup centroid genes were compared across male and female breast cancers (Additional files 5 and 6). When the MBC tumors were clustered with the ER+ FBC subgroup centroid genes, the two MBC subgroups were mixed between the main clusters found (Additional file 7). These findings indicate that the MBC subgroups differ substantially from the intrinsic subgroups of FBC.
NAT1 protein expression is prognostic for MBC
Uni- and multi-variate analysis of five-year breast cancer specific survival
1.0 to 5.9
1.0 to 7.2
No vs . Yes
0.33 to 5.8
< 45 vs . ≥ 45
1.3 to 6.6
0.82 to 6.8
T2 vs . T1
1.4 to 7.9
1.1 to 8.3
Pos vs . Neg
1.1 to 5.0
0.40 to 2.3
3 vs . 1 to 2
0.22 to 1.1
Yes vs . No
The correlation between immune response and prognosis
Significant up-regulation of immune-related genes was observed in luminal M2 compared to luminal M1 tumors (Additional file 4), and a higher correlation to the immune module was also seen (Figure 2A, C). More specifically, class 1 HLA genes were strongly up-regulated in luminal M2 tumors. To explore the prognostic impact of immune-related genes we therefore assessed the protein levels of HLA in the extended cohort. Tumors were defined as positive, moderate or negative when > 50%, 5 to 50% or < 5% of the cancer cells were positive, respectively (Additional file 8D-F). Ninety-three (42%) tumors were positive, 51 (23%) displayed moderate staining and 58 (26%) were negative (data were missing for 21 (9%) tumors). A significant difference in the protein expression was observed, with the luminal M2 subgroup containing more HLA positive tumors than the luminal M1 subgroup (Fisher's exact test, P = 0.039). Most interestingly, HLA positivity was associated with significantly better DMFS than moderate HLA expression (HR 3.6 (95% CI 1.6 to 7.9) P = 0.002). The DMFS was similar in the HLA negative and the HLA positive groups during the beginning of the follow-up, but late events (> 8 years) among HLA negative cases nevertheless resulted in poor DMFS in the latter group, comparable to the moderate HLA group (Figure 4B).
While it is generally accepted that FBC is a heterogeneous disease both in terms of transcriptional profiles, genomic aberrations and survival [19–22], whether MBC can be classified into comprehensive subgroups associated with differences in clinicopathological variables has not yet been elucidated. The inferior outcome reported clearly indicates the requirement to better understand the pathobiology of MBC , and the need to stratify patients based on tumor characteristics, potentially in need of alternative treatment strategies.
In this study, we have performed GEX profiling on 66 MBC tumors to study the disease transcriptionally and to subclassify tumors based on mRNA profiles. We were thus able to stratify tumors into two stable subgroups, luminal M1 and luminal M2, respectively, associated with different biological features and clinicopathological characteristics. While additional subgroups may exist, the sample size was too small for further subdivisions. As evidenced from the SAM analysis, the two GEX subgroups displayed a large number of differentially expressed genes. Patients tended to be diagnosed at a younger age in the luminal M1 subgroup, but no differences in histological grade or tumor size were observed. Interestingly, the two subgroups correlated with the genomic subgroups (male-simple and male-complex, respectively) that we previously defined based on genomic aberrations within the same patient cohort . Specifically, the luminal M2 subgroup described in the present study contained the majority of the male-simple tumors, while the luminal M1 subgroup harbored most of the male-complex tumors. There was a tendency towards better DMFS among patients with male-simple tumors with few genomic aberrations compared to patients with male-complex tumors with highly rearranged genomes and survival comparable to luminal B FBCs. Male-simple tumors were less frequently aneuploid, displayed a lower fraction of genome altered as well as lower S-phase fractions , further supporting better outcome among patients with luminal M2 (male-simple) vs. luminal M1 (male-complex) tumors.
To better understand the biological differences between the two subgroups, we performed an ontology analysis, whereupon up-regulation of genes involved in cell migration, cell adhesion, angiogenesis, cell cycle, cell division and HOX genes was identified among luminal M1 tumors. Two of the well-known hallmarks of cancer development are activation of invasion and metastasis and induction of angiogenesis, respectively . Up-regulation of genes involved in these processes in the luminal M1 subgroup indicates that these tumors may have a more aggressive phenotype than luminal M2 tumors. HOX genes are DNA binding factors involved in the transcriptional regulation of many key development factors, and de-regulated expression of HOX genes has been found to be involved in carcinogenesis and metastasis in many different cancer forms, including breast cancer . In the present study, HOXB7, which has been shown to be involved in epithelial-mesenchymal transition, migration and invasion in FBC , was significantly up-regulated in luminal M1 tumors compared to luminal M2 tumors. This finding also corresponds to the higher frequency of distant metastases in the former group.
An additional hallmark of cancer pertains to activation of the immune system, specifically by T and B lymphocytes, macrophages and natural killer cells . To this end, luminal M2 tumors displayed high expression of genes associated with class I histocompatibility antigens, which are involved in regulating the immune response, further supporting a more favorable outcome among luminal M2 tumors. A Kaplan-Meier survival analysis indeed suggested better DMFS among luminal M2 tumors, lending additional support to this notion.
To further understand the biological differences between the MBC subgroups, we investigated seven key biological processes associated with differences in survival among the intrinsic subtypes of FBC by calculating scores for each module in the two MBC subgroups. The luminal M2 subgroup demonstrated higher scores for the immune response and ER modules, while luminal M1 tumors displayed higher scores for the tumor invasion and metastasis, proliferation and HER2 modules, again indicating differences in tumor aggressiveness between the subgroups. A surprising finding was the low score for the ER module among luminal M1 tumors, despite almost all MBC tumors in the present study being ER+. In FBC, only ER- tumors display a low ER module score, suggesting that luminal M1 MBC tumors, although positive by IHC, in fact differ from ER+ FBCs. Luminal M2 MBC tumors also appear to differ from the FBC intrinsic subgroups, as the correlations between module scores and ER+ FBC intrinsic subgroups were not recapitulated. While some similarities to both luminal A and B subgroups were observed among luminal M2 tumors, a high module score for immune response was also seen, a feature only associated with the HER2 and basal intrinsic subtypes of FBC. These findings underscore the difficulty in capturing the complexity of molecular alterations associated with MBC subtypes using single protein markers. To date, only two other studies have attempted to subclassify MBC into the major FBC intrinsic subtypes. Applying IHC and the commonly used protein markers for FBC subtyping, approximately 80% of the tumors were classified as luminal A and approximately 20% as luminal B [41, 42]. The discrepancy with our findings is probably due to the inability of a small number of protein markers to fully capture the differences in transcriptional profiles between subtypes. This is illustrated by, for example, luminal M1 tumors being ER+ by IHC, while displaying less active ER signaling, more similar to ER- FBC.
Several studies of FBC have underlined the importance of regulation of the immune system in ER- and HER2 positive (HER2+) tumors [43–45]. Teschendorf et al. identified an immune response related seven-gene signature in ER- tumors correlated to risk of distant metastases [43, 44]. Further, Staaf et al. defined a predictor prognostic of outcome for HER2+ FBC tumors that included genes associated with immune response, tumor invasion and metastasis. The better prognosis group displayed up-regulation of the immune response and low invasive ability. Of interest, this predictor also performed well in ER- FBC, but not in ER+/HER2- FBC . Given the correlation to the immune module in the luminal M2 subgroup of ER+ MBCs, and the lack of association with the ER signaling module in the luminal M1 subgroup, these findings indicate that the two subgroups of MBC described herein may constitute two new subgroups of breast cancer, with unique biological and clinical features, occurring only in males. Hypothetically, these patients may require novel treatment strategies. Specifically, despite the majority of tumors in the luminal M1 subgroup being ER+ they had a low ER signaling module score, suggesting that the ER pathway may not be active. A recent comprehensive study of steroid hormone receptors in breast cancer revealed gender specific differences, suggesting differential hormonal dependency . Hence, whether these MBC patients respond to endocrine therapy like tamoxifen may be questioned, and needs to be further investigated in prospective randomized studies. Intriguingly, a recent study implicated HOXB7, one of the genes up-regulated in luminal M1 MBC tumors, in rendering FBC cells resistant to tamoxifen. In addition, high expression of HOXB7 in tamoxifen treated FBC patients correlated with poor disease free survival .
Importantly, we were able to validate the MBC subgroups in an independent dataset ; the centroid classification in this dataset resulted in the same distribution into the two subgroups as in our dataset, and unsupervised clustering revealed two stable subgroups with similar characteristics as those found in our dataset. Specifically, the GEX patterns of the centroid genes were similar and the module scores showed the same trends in the validation dataset, indicating similar associations with biological processes. Unfortunately, however, no information on outcome was available from that study.
A comparison between the MBC subgroups and a dataset representing all intrinsic subgroups of FBC revealed pronounced differences; 55% and 36% of the MBC samples were unclassified when applying the Hu gene centroid and FBC ER+ luminal gene centroid classifications, respectively. The fraction of unclassified tumors among FBCs has previously been reported to be 0 to 20% . The GEX patterns based on these genes also differed between male and female breast tumors, further indicating that the subgroups of MBC identified herein are not represented by the intrinsic subgroups of FBC. When our MBC samples were clustered based on the ER+ FBC luminal centroid genes, the two MBC subgroups were intermixed between the clusters. In support of this finding, our previous study of genomic profiles in MBC also revealed that the two MBC subgroups largely differed from genomic subgroups described in FBC . Although some of the centroid genes in ER+ luminal FBC may differ between the MBC subgroups, they are clearly not the most dominant feature.
To further characterize the subgroups and investigate the clinical relevance of our findings, we investigated NAT1 and HLA protein levels in a series of 220 MBCs, as differences in mRNA levels of the corresponding genes were found between the subgroups. A significantly worse DMFS was observed for the moderate HLA group compared with the positive HLA group. Curiously, however, the HLA negative group initially had a prognosis similar to the HLA positive group, but late recurrences resulted in worse long-time survival. Due to the fairly small sample size and the fact that not many patients remained alive after eight years, this finding needs to be interpreted cautiously and more studies are needed to validate these findings. Interestingly, women with node-negative breast cancer displaying a mixed HLA class I expression pattern had a worse prognosis according to a study by Gudmundsdóttir et al. . Tumors with high expression of HLA class I can be recognized by T lymphocytes of the specific immune system. On the other hand, tumors lacking expression of HLA class I may be targeted by NK cells of the non-specific immune system [50, 51]. Tumors with a mixed expression of HLA class I may thus hypothetically be able to avoid the specific immune system by displaying too few antigens, while evading the non-specific immune system by inhibiting NK cell activity [49, 51].
Several studies have reported higher expression of NAT1 on both protein and mRNA levels in ER+ FBC compared to ER- FBC [20, 52–54], and high expression of NAT1 has been shown to correlate with better outcome among ER+ FBCs [55, 56]. The NAT1 antibody used in the present study has been used over the past 20 years with consistent results on cellular localization [35, 36, 55], and has been demonstrated to be uniquely specific following Western blot analysis . However, there has also been an indication of a nuclear location , which has been linked to turnover of the NAT1 protein. The precise role of the small proportion of nuclear staining of NAT1 has not been established, but may well be related to protein turnover and gene regulation. Patients whose tumors were positive for NAT1 displayed a significantly better prognosis than those with NAT1 negative tumors in the present study, a finding that remained significant in a multivariate analysis. ER status was not included in the multivariate analysis, because only seven tumors were ER-. It has, however, been shown that ER status provides independent prognostic information in MBC . Luminal M2 tumors displayed higher NAT1 levels compared to luminal M1 tumors, thus further supporting an association between the luminal M2 subgroup and better outcome. Bièche et al. found high NAT1 to be predictive of response to tamoxifen in women with ER+ breast cancer. In general, altered tamoxifen metabolism and bioavailability may contribute to tamoxifen resistance, and the xenobiotic-metabolizing enzyme NAT1 may be part of this explanation . Unfortunately, we had too few cases to be able to detect any association between NAT1 and survival among only tamoxifen treated patients, but this finding warrants further investigation and suggests that the low expression of NAT1 among luminal M1 tumors may lead to tamoxifen resistance, and that these patients may hence require alternative treatment approaches.
We have detected two unique and stable subgroups of MBC with differences in tumor biological features and outcome, luminal M1 tumors being more aggressive and associated with worse prognosis, while luminal M2 tumors, on the other hand, displayed up-regulated immune response and activated ER signaling, generally favorable features. Importantly, both MBC subgroups differed from the established intrinsic subgroups of FBC, indicating that they constitute two novel subgroups of breast cancer occurring only in males. Consequently, men diagnosed with breast cancer may require other management and treatment strategies than women. Finally, we identified NAT1 as a possible prognostic biomarker for MBC, as suggested by NAT1 positivity corresponding to better outcome.
array-based comparative genomic hybridization
BioArray Software Environment
comparative genomic hybridization
The Database for Annotation, Visualization and Integrated Discovery
distant metastasis free survival
female breast cancer
false discovery rate
Gene Expression Omnibus
human epidermal growth factor receptor 2
class 1 human leukocyte antigen
male breast cancer
Nottingham Histological Grade
principal component analysis
Re-annotation and Mapping for Oligonucleotide Array Technologies
significance analysis of microarray
The authors wish to thank all involved Pathology Departments for providing tissue for the study, and Kristina Lövgren for technical assistance. This study was supported by grants from the Swedish Cancer Society, the G Nilsson Cancer Foundation, the Mrs. B Kamprad Foundation, and the Lund University Hospital Research Foundation. IH was supported by the Swedish Cancer Society, and the SCIBLU Genomics center is supported by governmental funding of clinical research within the national health services (ALF) and by Lund University.
- Engholm G, Ferlay J, Christensen N, Bray F, Gjerstorff ML, Klint Å, Køtlum JE, Ólafsdóttir E, Pukkala E, Storm HH: NORDCAN: Cancer Incidence, Mortality, Prevalence and Prediction in the Nordic Countries. Acta Oncologica. 2009, 49: 532-544.Google Scholar
- Giordano SH: A review of the diagnosis and management of male breast cancer. Oncologist. 2005, 10: 471-479. 10.1634/theoncologist.10-7-471.View ArticleGoogle Scholar
- Agrawal A, Ayantunde AA, Rampaul R, Robertson JF: Male breast cancer: a review of clinical management. Breast Cancer Res Treat. 2007, 103: 11-21. 10.1007/s10549-006-9356-z.View ArticleGoogle Scholar
- Korde LA, Zujewski JA, Kamin L, Giordano S, Domchek S, Anderson WF, Bartlett JM, Gelmon K, Nahleh Z, Bergh J, Cutuli B, Pruneri G, McCaskill-Stevens W, Gralow J, Hortobagyi G, Cardoso F: Multidisciplinary meeting on male breast cancer: summary and research recommendations. J Clin Oncol. 2010, 28: 2114-2122. 10.1200/JCO.2009.25.5729.View ArticleGoogle Scholar
- Rudlowski C, Friedrichs N, Faridi A, Fuzesi L, Moll R, Bastert G, Rath W, Buttner R: Her-2/neu gene amplification and protein expression in primary male breast cancer. Breast Cancer Res Treat. 2004, 84: 215-223. 10.1023/B:BREA.0000019953.92921.7e.View ArticleGoogle Scholar
- Nilsson C, Johansson I, Ahlin C, Thorstensson s, Amini RM, Holmqvist M, Hedenfalk I, Fjällskog M-L: Molecular subtyping of male breast cancer using alternative definitions and its prognostic impact. Submitted,
- Ottini L, Rizzolo P, Zanna I, Falchetti M, Masala G, Ceccarelli K, Vezzosi V, Gulino A, Giannini G, Bianchi S, Sera F, Palli D: BRCA1/BRCA2 mutation status and clinical-pathologic features of 108 male breast cancer cases from Tuscany: a population-based study in central Italy. Breast Cancer Res Treat. 2009, 116: 577-586. 10.1007/s10549-008-0194-z.View ArticleGoogle Scholar
- Brinton LA, Richesson DA, Gierach GL, Lacey JV, Park Y, Hollenbeck AR, Schatzkin A: Prospective evaluation of risk factors for male breast cancer. J Natl Cancer Inst. 2008, 100: 1477-1481. 10.1093/jnci/djn329.View ArticleGoogle Scholar
- Nilsson C, Holmqvist M, Bergkvist L, Hedenfalk I, Lambe M, Fjällskog M-L: Similarities and differences in the characteristics and treatment of breast cancer in men and women-a population based study (Sweden). Acta Oncol. 2011, 50: 1083-1088. 10.3109/0284186X.2011.602114.View ArticleGoogle Scholar
- Anderson WF, Devesa SS: In situ male breast carcinoma in the Surveillance, Epidemiology, and End Results database of the National Cancer Institute. Cancer. 2005, 104: 1733-1741. 10.1002/cncr.21353.View ArticleGoogle Scholar
- Thalib L, Hall P: Survival of male breast cancer patients: Population-based cohort study. Cancer Sci. 2008, 100 (2): 292-295.View ArticleGoogle Scholar
- Giordano SH, Cohen DS, Buzdar AU, Perkins G, Hortobagyi GN: Breast carcinoma in men: a population-based study. Cancer. 2004, 101: 51-57. 10.1002/cncr.20312.View ArticleGoogle Scholar
- Giordano SH: Male breast cancer: it's time for evidence instead of extrapolation. Onkologie. 2008, 31: 505-506.View ArticleGoogle Scholar
- Johansson I, Nilsson C, Berglund P, Strand C, Jonsson G, Staaf J, Ringner M, Nevanlinna H, Barkardottir RB, Borg A, Olsson H, Luts L, Fjallskog ML, Hedenfalk I: High-resolution genomic profiling of male breast cancer reveals differences hidden behind the similarities with female breast cancer. Breast Cancer Res Treat. 2011, 129: 747-760. 10.1007/s10549-010-1262-8.View ArticleGoogle Scholar
- Callari M, Cappelletti V, De Cecco L, Musella V, Miodini P, Veneroni S, Gariboldi M, Pierotti MA, Daidone MG: Gene expression analysis reveals a different transcriptomic landscape in female and male breast cancer. Breast Cancer Res Treat. 2011, 127: 601-610. 10.1007/s10549-010-1015-8.View ArticleGoogle Scholar
- Lehmann U, Streichert T, Otto B, Albat C, Hasemeier B, Christgen H, Schipper E, Hille U, Kreipe HH, Langer F: Identification of differentially expressed microRNAs in human male breast cancer. BMC Cancer. 2010, 10: 109-10.1186/1471-2407-10-109.View ArticleGoogle Scholar
- Fassan M, Baffa R, Palazzo JP, Lloyd J, Crosariol M, Liu CG, Volinia S, Alder H, Rugge M, Croce CM, Rosenberg A: MicroRNA expression profiling of male breast cancer. Breast Cancer Res. 2009, 11: R58-10.1186/bcr2348.View ArticleGoogle Scholar
- Tommasi S, Mangia A, Iannelli G, Chiarappa P, Rossi E, Ottini L, Mottolese M, Zoli W, Zuffardi O, Paradiso A: Gene copy number variation in male breast cancer by aCGH. Anal Cell Pathol (Amst). 2010, 33 (3): 113-119.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Hu Z, Fan C, Oh DS, Marron JS, He X, Qaqish BF, Livasy C, Carey LA, Reynolds E, Dressler L, Nobel A, Parker J, Ewend MG, Sawyer LR, Wu J, Liu Y, Nanda R, Tretiakova M, Ruiz Orrico A, Dreher D, Palazzo JP, Perreard L, Nelson E, Mone M, Hansen H, Mullins M, Quackenbush JF, Ellis MJ, Olopade OI, Bernard PS, Perou CM: The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics. 2006, 7: 96-10.1186/1471-2164-7-96.View ArticleGoogle Scholar
- Jonsson G, Staaf J, Vallon-Christersson J, Ringner M, Holm K, Hegardt C, Gunnarsson H, Fagerholm R, Strand C, Agnarsson BA, Kilpivaara O, Luts L, Heikkila P, Aittomaki K, Blomqvist C, Loman N, Malmstrom P, Olsson H, Johannsson OT, Arason A, Nevanlinna H, Barkardottir RB, Borg A: Genomic subtypes of breast cancer identified by array comparative genomic hybridization display distinct molecular and clinical characteristics. Breast Cancer Res. 2010, 12: R42-10.1186/bcr2596.View ArticleGoogle Scholar
- Vallon-Christersson J, Nordborg N, Svensson M, Hakkinen J: BASE--2nd generation software for microarray data management and analysis. BMC Bioinformatics. 2009, 10: 330-10.1186/1471-2105-10-330.View ArticleGoogle Scholar
- The R Project for Statistical Computing. [http://www.r-project.org]
- Johnson WE, Li C, Rabinovic A: Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007, 8: 118-127.View ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.View ArticleGoogle Scholar
- Shimodaira H: Approximately unbiased tests of regions using multistep-multiscale bootstrap resampling. Ann Stat. 2004, 32 (6): 26-View ArticleGoogle Scholar
- Lindgren D, Frigyesi A, Gudjonsson S, Sjodahl G, Hallden C, Chebil G, Veerla S, Ryden T, Mansson W, Liedberg F, Hoglund M: Combined gene expression and genomic profiling define two intrinsic molecular subtypes of urothelial carcinoma and gene signatures for molecular grading and outcome. Cancer Res. 2010, 70: 3463-3472. 10.1158/0008-5472.CAN-09-4213.View ArticleGoogle Scholar
- Barbosa-Morais NL, Dunning MJ, Samarajiwa SA, Darot JF, Ritchie ME, Lynch AG, Tavare S: A re-annotation pipeline for Illumina BeadArrays: improving the interpretation of gene expression data. Nucleic Acids Res. 2010, 38: e17-10.1093/nar/gkp942.View ArticleGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.View ArticleGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticleGoogle Scholar
- Desmedt C, Haibe-Kains B, Wirapati P, Buyse M, Larsimont D, Bontempi G, Delorenzi M, Piccart M, Sotiriou C: Biological processes associated with breast cancer clinical outcome depend on the molecular subtypes. Clin Cancer Res. 2008, 14: 5158-5165. 10.1158/1078-0432.CCR-07-4756.View ArticleGoogle Scholar
- Stam NJ, Vroom TM, Peters PJ, Pastoors EB, Ploegh HL: HLA-A- and HLA-B-specific monoclonal antibodies reactive with free heavy chains in western blots, in formalin-fixed, paraffin-embedded tissue sections and in cryo-immuno-electron microscopy. Int Immunol. 1990, 2: 113-125. 10.1093/intimm/2.2.113.View ArticleGoogle Scholar
- Stam NJ, Spits H, Ploegh HL: Monoclonal antibodies raised against denatured HLA-B locus heavy chains permit biochemical characterization of certain HLA-C locus products. J Immunol. 1986, 137: 2299-2306.Google Scholar
- Stanley LA, Copp AJ, Pope J, Rolls S, Smelt V, Perry VH, Sim E: Immunochemical detection of arylamine N-acetyltransferase during mouse embryonic development and in adult mouse brain. Teratology. 1998, 58: 174-182. 10.1002/(SICI)1096-9926(199811)58:5<174::AID-TERA3>3.0.CO;2-Q.View ArticleGoogle Scholar
- Stanley LA, Coroneos E, Cuff R, Hickman D, Ward A, Sim E: Immunochemical detection of arylamine N-acetyltransferase in normal and neoplastic bladder. J Histochem Cytochem. 1996, 44: 1059-1067. 10.1177/44.9.8773572.View ArticleGoogle Scholar
- Suzuki R, Shimodaira H: Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics. 2006, 22: 1540-1542. 10.1093/bioinformatics/btl117.View ArticleGoogle Scholar
- Hanahan D, Weinberg RA: Hallmarks of cancer: the next generation. Cell. 2011, 144: 646-674. 10.1016/j.cell.2011.02.013.View ArticleGoogle Scholar
- Shah N, Sukumar S: The Hox genes and their roles in oncogenesis. Nat Rev Cancer. 2010, 10: 361-371. 10.1038/nrc2826.View ArticleGoogle Scholar
- Wu X, Chen H, Parker B, Rubin E, Zhu T, Lee JS, Argani P, Sukumar S: HOXB7, a homeodomain protein, is overexpressed in breast cancer and confers epithelial-mesenchymal transition. Cancer Res. 2006, 66: 9527-9534. 10.1158/0008-5472.CAN-05-4470.View ArticleGoogle Scholar
- Kornegoor R, Verschuur-Maes AH, Buerger H, Hogenes MC, de Bruin PC, Oudejans JJ, van der Groep P, Hinrichs B, van Diest PJ: Molecular subtyping of male breast cancer by immunohistochemistry. Mod Pathol. 2012, 25: 398-404. 10.1038/modpathol.2011.174.View ArticleGoogle Scholar
- Ge Y, Sneige N, Eltorky MA, Wang Z, Lin E, Gong Y, Guo M: Immunohistochemical characterization of subtypes of male breast carcinoma. Breast Cancer Res. 2009, 11: R28-10.1186/bcr2258.View ArticleGoogle Scholar
- Teschendorff AE, Miremadi A, Pinder SE, Ellis IO, Caldas C: An immune response gene expression module identifies a good prognosis subtype in estrogen receptor negative breast cancer. Genome Biol. 2007, 8: R157-10.1186/gb-2007-8-8-r157.View ArticleGoogle Scholar
- Teschendorff AE, Caldas C: A robust classifier of high predictive value to identify good prognosis patients in ER-negative breast cancer. Breast Cancer Res. 2008, 10: R73-10.1186/bcr2138.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Shaaban AM, Ball GR, Brannan RA, Cserni G, Benedetto AD, Dent J, Fulford L, Honarpisheh H, Jordan L, Jones JL, Kanthan R, Maraqa L, Litwiniuk M, Mottolese M, Pollock S, Provenzano E, Quinlan PR, Reall G, Shousha S, Stephens M, Verghese ET, Walker RA, Hanby AM, Speirs V: A comparative biomarker study of 514 matched cases of male and female breast cancer reveals gender-specific biological differences. Breast Cancer Res Treat. 2011,Google Scholar
- Jin K, Kong X, Shah T, Penet MF, Wildes F, Sgroi DC, Ma XJ, Huang Y, Kallioniemi A, Landberg G, Bieche I, Wu X, Lobie PE, Davidson NE, Bhujwalla ZM, Zhu T, Sukumar S: The HOXB7 protein renders breast cancer cells resistant to tamoxifen through activation of the EGFR pathway. Proc Natl Acad Sci USA. 2011,Google Scholar
- Millikan RC, Newman B, Tse CK, Moorman PG, Conway K, Dressler LG, Smith LV, Labbok MH, Geradts J, Bensen JT, Jackson S, Nyante S, Livasy C, Carey L, Earp HS, Perou CM: Epidemiology of basal-like breast cancer. Breast Cancer Res Treat. 2008, 109: 123-139. 10.1007/s10549-007-9632-6.View ArticleGoogle Scholar
- Gudmundsdottir I, Gunnlaugur Jonasson J, Sigurdsson H, Olafsdottir K, Tryggvadottir L, Ogmundsdottir HM: Altered expression of HLA class I antigens in breast cancer: association with prognosis. Int J Cancer. 2000, 89: 500-505. 10.1002/1097-0215(20001120)89:6<500::AID-IJC6>3.0.CO;2-#.View ArticleGoogle Scholar
- Karre K: Express yourself or die: peptides, MHC molecules, and NK cells. Science. 1995, 267: 978-979. 10.1126/science.7863341.View ArticleGoogle Scholar
- Garrido F, Ruiz-Cabello F, Cabrera T, Perez-Villar JJ, Lopez-Botet M, Duggan-Keen M, Stern PL: Implications for immunosurveillance of altered HLA class I phenotypes in human tumours. Immunol Today. 1997, 18: 89-95. 10.1016/S0167-5699(96)10075-X.View ArticleGoogle Scholar
- Wakefield L, Robinson J, Long H, Ibbitt JC, Cooke S, Hurst HC, Sim E: Arylamine N-acetyltransferase 1 expression in breast cancer cell lines: a potential marker in estrogen receptor-positive tumors. Genes Chromosomes Cancer. 2008, 47: 118-126. 10.1002/gcc.20512.View ArticleGoogle Scholar
- Tozlu S, Girault I, Vacher S, Vendrell J, Andrieu C, Spyratos F, Cohen P, Lidereau R, Bieche I: Identification of novel genes that co-cluster with estrogen receptor alpha in breast tumor biopsy specimens, using a large-scale real-time reverse transcription-PCR approach. Endocr Relat Cancer. 2006, 13: 1109-1120. 10.1677/erc.1.01120.View ArticleGoogle Scholar
- Adam PJ, Berry J, Loader JA, Tyson KL, Craggs G, Smith P, De Belin J, Steers G, Pezzella F, Sachsenmeir KF, Stamps AC, Herath A, Sim E, O'Hare MJ, Harris AL, Terrett JA: Arylamine N-acetyltransferase-1 is highly expressed in breast cancers and conveys enhanced growth and resistance to etoposide in vitro. Mol Cancer Res. 2003, 1: 826-835.Google Scholar
- Dolled-Filhart M, Ryden L, Cregger M, Jirstrom K, Harigopal M, Camp RL, Rimm DL: Classification of breast cancer using genetic algorithms and tissue microarrays. Clin Cancer Res. 2006, 12: 6459-6468. 10.1158/1078-0432.CCR-06-1383.View ArticleGoogle Scholar
- Bieche I, Girault I, Urbain E, Tozlu S, Lidereau R: Relationship between intratumoral expression of genes coding for xenobiotic-metabolizing enzymes and benefit from adjuvant tamoxifen in estrogen receptor alpha-positive postmenopausal breast carcinoma. Breast Cancer Res. 2004, 6: R252-263. 10.1186/bcr784.View ArticleGoogle Scholar
- Butcher NJ, Arulpragasam A, Minchin RF: Proteasomal degradation of N-acetyltransferase 1 is prevented by acetylation of the active site cysteine: a mechanism for the slow acetylator phenotype and substrate-dependent down-regulation. J Biol Chem. 2004, 279: 22131-22137. 10.1074/jbc.M312858200.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.