Next Article in Journal
Innovative MRI Techniques in Neuroimaging Approaches for Cerebrovascular Diseases and Vascular Cognitive Impairment
Next Article in Special Issue
Cadmium Exposure and Risk of Breast Cancer by Histological and Tumor Receptor Subtype in White Caucasian Women: A Hospital-Based Case-Control Study
Previous Article in Journal
Impact and Relevance of the Unfolded Protein Response in HNSCC
Previous Article in Special Issue
Lipid Metabolism and Endocrine Resistance in Prostate Cancer, and New Opportunities for Therapy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 11
10.3390/ijms20112655
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estrogen Receptor Positive Breast Cancer with High Expression of Androgen Receptor has Less Cytolytic Activity and Worse Response to Neoadjuvant Chemotherapy but Better Survival

by
Maiko Okano
1,
Masanori Oshi
1,
Ali Linsk Butash
1,
Mariko Asaoka
1,
Eriko Katsuta
1,
Xuan Peng
2,
Qianya Qi
2,
Li Yan
2 and
Kazuaki Takabe
1,3,4,5,6,7,*
1
Breast Surgery, Department of Surgical Oncology, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
2
Department of Biostatistics & Bioinformatics, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
3
Department of Surgery, University at Buffalo Jacobs School of Medicine and Biomedical Sciences, The State University of New York, Buffalo, NY 14203-1121, USA
4
Department of Breast Surgery and Oncology, Tokyo Medical University, Tokyo 160-8402, Japan
5
Department of Surgery, Yokohama City University, Yokohama 236-0004, Japan
6
Department of Surgery, Niigata University Graduate School of Medical and Dental Sciences, Niigata 951-8510, Japan
7
Department of Breast Surgery, Fukushima Medical University School of Medicine, Fukushima 960-1295, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(11), 2655; https://doi.org/10.3390/ijms20112655
Submission received: 22 April 2019 / Revised: 24 May 2019 / Accepted: 27 May 2019 / Published: 30 May 2019
(This article belongs to the Special Issue Hormones-Dependent Cancers: New Aspects on Biochemistry and Molecular Pathology (II))
Browse Figures
Review Reports Versions Notes

Abstract

:
Estrogen receptor (ER) positive breast cancer (BC), the most abundant BC subtype, is notorious for poor response to neoadjuvant chemotherapy (NAC). The androgen receptor (AR) was reported to support estradiol-mediated ER activity in an in vitro system. Recently, ER-positive BC with fewer tumor infiltrating lymphocytes (TILs) was shown to have a better prognosis, opposite to the trend seen with ER-negative BC. We hypothesized that ER-positive BC with high expression of AR will have fewer TILs and an inferior response to NAC, but with a better prognosis. In both TCGA and METABRIC cohorts, AR expression was significantly higher in ER-positive BCs compared to ER-negatives (p < 0.001, p < 0.001, respectively) and it correlated with ER expression (R = 0.630, R = 0.509, respectively). In ER-positive tumors, AR high tumors enriched UV response down (NES = 2.01, p < 0.001), and AR low tumors enriched DNA repair (NES = −2.02, p < 0.001). AR high tumors were significantly associated with procancer regulatory T-cells, and AR low tumors were associated with anticancer immune cells, such as CD4, CD8, and Gamma-Delta T-cells and memory B-cells in ER-positive BC (p < 0.01). Further, cytolytic activity was significantly lower in AR high BC in both cohorts. Finally, AR high tumors had a significantly lower rate of attaining pathological complete response to NAC (GSE22358), but better survival. In conclusion, our results demonstrated that high AR has fewer tumor infiltrating lymphocytes as well as cytolytic activity and an inferior response to NAC, but better survival in ER-positive BC.

1. Introduction

Breast cancer (BC) is the most common cancer among women in the United States excluding skin cancers. About one in eight U.S. women (about 12%) will develop invasive BC in their lifetime [1]. BC is categorized into clinical subtypes by receptor expression statuses with distinctive phenotypes, including varying response to certain treatments [2]. Estrogen receptor (ER) positive BC is the most abundant subtype that accounts for 70% of BC. Although ER-positive BC has a relatively better prognosis compared with the other subtypes [3,4], it is well known to have less tumor infiltrating lymphocytes (TIL) [5,6], and less responsiveness to neoadjuvant chemotherapy (NAC) [7,8]. Interestingly, there have been multiple reports that demonstrated better survival with less TILs in ER-positive BC, which is quite opposite to what is known in ER-negative tumors [9,10].
The androgen receptor (AR) is one of the members of the steroid nuclear receptor family, which includes ER and progesterone receptor (PR). AR is expressed in 50–90% of breast cancers [11,12] and some claim that it is more common than ER or PR [13]. Unlike ER or PR, however, the role of AR in BC progression is mechanistically complex and remains highly controversial.
AR would rather be famous as the potential target in the therapy for triple negative breast cancer (TNBC), which does not express ER, PR, and HER2/neu and tends to be more aggressive than other subtypes of breast cancers. TNBC is considered not to be a single disease, and is classified into several molecular subtypes, among which there are luminal AR subtypes (LAR) [14]. In vitro experiments have shown that AR is associated with the growth of TNBC cell lines via Wnt/β-catenin signaling and G-protein coupled estrogen receptor activation [15,16]. AR inhibitors are expected to be effective for TNBC patients otherwise lacking molecular targets and several clinical trials have been performed using the AR antagonist for TNBC patients.
On the other hand, several studies have reported the role of AR in ER-positive BC. A study had identified a crosstalk between AR and epidermal growth factor receptor (EGFR), which is involved in cancer cell growth, in ER-positive BC cells [17]. Another report showed that specific domains of AR and ERα physically interact and demonstrated the functional consequences of such interaction using yeast and mammalian two-hybrid systems [18]. AR was also demonstrated to support estradiol-mediated ER activity in BC cells that express both ER and AR [13], which led to clinical trials currently underway using AR-targeting agents (NCT00468715, NCT02007512) [19]. Further, ER-positive BC with high AR expression showed resistance to tamoxifen and aromatase inhibitors in both in vitro and in vivo systems [20,21]. However, several papers that analyzed many patient data concluded that AR is a favorable marker for survival [22,23].
We hypothesized that ER-positive BC with high levels of AR expression attracts TILs and has a worse response to NAC, but better patient survival. In this study, we investigated the associations of AR mRNA expression with ER expression, the number of TILs, the response to NAC, and survival using gene expression data from the publicly available large cohorts.

2. Results

2.1. Expression of AR mRNA was Significantly Higher in ER-Positive BC

AR expression has been reported to be higher in ER-positive compared with ER-negative BC; however, previous studies used immunohistochemistry (IHC) for AR measurement, which is semiquantitative and the antibodies used were not standardized [12,24,25]. We used transcriptomic data that is rigorously quantitative. AR mRNA expression was significantly lower in basal subtypes among PAM50 classification determined by the gene expression profile (Figure 1A). As expected, ER-positive tumors classified by IHC in TCGA, expressed significantly higher levels of AR mRNA expression when compared to ER-negative patients (Figure 1B). Strikingly, these results were mirrored in the completely independent validation cohort, METABRIC (p < 0.001, Figure 1C,D).

2.2. AR mRNA Expression Correlates with ER mRNA Expression

Since several papers have shown the interaction of AR and ER, it was of interest whether high expression of AR mRNA correlates with that of ER mRNA (Figure 2). We found it not only in TCGA (r = 0.630, p < 0.001), but also in METABRIC (r = 0.509, p < 0.001) as well.

2.3. UV Response Down and DNA Repair Gene Sets Were Enriched in AR High and AR Low Tumors, Respectively, in ER-Positive BC

We conducted GSEA to investigate which gene sets are associated with AR high tumors or AR low tumors in ER-positive BC. High AR expression tumors enriched UV response down (NES = 2.01, p < 0.001) hallmark gene sets in GSEA (Figure 3 and Table 1). On the other hand, AR low expression tumors significantly enriched DNA repair (NES = −2.03, p < 0.001) hallmark gene sets in GSEA (Figure 3B and Table 2).

2.4. Activated Memory CD4 T Cells and Gamma Delta T cells Were Significantly Lower in AR High Tumors Compared with AR Low Tumors in ER-Positive BC

In order to investigate the relationship between TIL and AR expression in ER-positive BC, immune cell composition was estimated using the CIBERSORT algorithm in TCGA and METABRIC cohorts (Figure 4). We found that AR high tumors associated with significantly low infiltration of CD4 T cells and Gamma Delta T cells (p < 0.001) consistently in both TCGA and METABRIC cohorts (Figure 4). The relationship between AR expression and CD8 T cells or B cells were inconsistent in those cohorts. The infiltration of regulatory T cells, which is known to suppress anticancer immunity and negatively impact patient survival, was also inconsistent between the two cohorts.

2.5. Tumor Heterogeneity and Cytolytic Activity Score (CYT) Were Significantly Lower in AR High Compared with AR Low Tumors in ER-Positive BC

The current dogma is that heterogenous tumors generate neoantigens that attract TILs, which are expected to attack cancer with cytolytic activity. Although ER-positive BC is known to have less TILs when compared to triple negative BC, it was of interest to investigate whether this mechanistic model applies in AR high tumors where we found lower infiltration of TILs among ER-positive BC. We found that AR high tumors were significantly associated with low mutant allele tumor heterogeneity (MATH) scores that assess tumor heterogeneity in TCGA whole cohort, but not in any other subtypes, including ER-positive BC (Figure 5 top row). This result implies that infiltration of TILs may not be due to tumor heterogeneity and neoantigen production in AR high tumors among ER-positive BC. On the other hand, CYT was significantly lower in AR high tumors not only in the whole cohort, but also in ER-positive BC and not in the other subtypes, and consistently in both TCGA and METABRIC cohorts (Figure 5 middle and lower rows). This result demonstrates that despite some discrepancy in the infiltration of CD8 T cells and regulatory T cells, AR high tumors possess less overall anticancer immunity.

2.6. AR High Tumors Demonstrated Lower Rates of Attaining Pathological Complete Response (pCR) to Neoadjuvant Chemotherapy (NAC), but Better Survival

In order to explore the relationship between AR expression and response to NAC and subsequent survival, ten NAC cohorts from the GEO database that reported gene expression data and pCR or distant relapse-free survival were analyzed (Figure 6). Although there was only one cohort that reached statistical significance (GSE22358), the pCR rate after NAC was unequivocally lower in AR high tumors in all ten ER-positive BC NAC cohorts (GSE20194, GSE21974, GSE22093, GSE22226, GSE23988, GSE25055, GSE25066, GSE42822, GSE50948). The two ER-positive BC NAC cohorts that have associated survival data, both demonstrated better survival in AR high tumors in ER-positive BC although statistical significance was achieved in only one (GSE25066), which is in agreement with previous reports [26,27].

3. Discussion

In the present study, we found that AR mRNA expression was higher in ER-positive tumors and the level of AR expression correlated with that of ER expression. AR high tumors enriched UV response down gene sets, whereas AR low tumors enriched DNA repair gene sets. AR high tumors consistently had less infiltration of CD4 and Gamma Delta T cells, but there was a discrepancy with infiltration of CD8 and regulatory T cells between the cohorts. AR high tumors had low tumor heterogeneity only in the whole cohort, whereas CYT was consistently low in the whole cohort and the ER-positive BC cohort in both TCGA and METABRIC. Finally, while AR high tumors demonstrated low pCR rates, they were associated with better survival after NAC.
Several studies have reported that AR is expressed higher in ER-positive primary BC than ER-negative BC [24,28]. D’Amato et al. demonstrated that estradiol-induced AR-binding sites were enriched for estrogen response elements with significant overlap with ER-binding sites, and that AR is required for maximum ER genomic binding [13]. Our finding that AR expression correlates with ER expression is in agreement with this notion that AR not only interacts with ER but also supports estradiol-mediated ER activity.
AR was demonstrated to protect prostate cancer cells from DNA damage in vitro [29,30]. AR was also shown to maintain production of DNA damage response molecules in some breast cancer cells and it prevents the accumulation of damaged DNA in vitro [31]. Recent investigations have demonstrated that high AR expression was significantly associated with worse local recurrence-free survival after radiation therapy, which is known to induce DNA damage, and AR can be as a mediator of radioresistance in triple negative breast cancer [32]. Our results are similar to these notions that AR is strongly associated with DNA damage, and AR high tumors associated with chromatin-remodeling genes which were downregulated by UV stress. Our result that AR low tumor significantly enriched DNA repair gene sets aligns with this notion that tumors with AR low expression have more DNA damage because of the lack of the protection by AR and it enriched DNA repair gene sets as a response. This notion is also consistent with our result that MATH score was high in AR low tumor.
Although there are not many articles investigating the relation between TIL and AR, it has been proposed that there is negative correlation between AR expression and immune cell infiltration in patients with HER2 and ER-positive tumors [33]. In the present study, among the infiltrating TILs, CD8 and regulatory T cells, which are known to be the major players, were not consistently low in AR high tumors. Together with the fact that tumor heterogeneity was not significantly lower in AR high tumors in the ER-positive BC cohort, the relationship between AR expression and TIL does not follow the dogma that heterogeneous tumor produces neoantigens that attract TILs [34]. On the other hand, CYT was consistently low in AR high tumors, which was associated with poor response to NAC. To this end, it is speculated that interaction between AR and ER may have some effect on tumor immune microenvironment separate from neoantigen-TIL dogma that suppresses cytolytic activity.
High infiltration of TILs is a known indicator of a good response to NAC and a better prognosis in the vast majority of cancers, including ER-negative BC [5,6]. This is partially explained by lymphocytes helping chemotherapy to fight cancer cells, and enhancement of anticancer immunity long after chemotherapy is completed. In contrast, recent studies have demonstrated that the opposite trend is the case for the ER-positive BC, where patients’ prognosis after NAC is better with fewer TILs [9,10]. Some argue that the prognosis of breast cancer with high AR is better because these cancers are well differentiated [27]. Our result is strikingly in agreement with these previous reports that AR high patients have better survival despite lower TIL infiltration and worse response to NAC. Based upon our results, we cannot help but speculate that high expression of AR that interacts with ER promotes a mitotic phenotype, less TIL, and less cytolytic activity that leads to NAC resistance. However, enhanced ER activity may have resulted in improvement of response to adjuvant hormonal therapy that consequently prolonged survival.
Our study has some limitations. The first and biggest limitation is that our research has demonstrated only the association between AR expression, TIL, and NAC, but there is no data that indicates causality. The second limitation is that TCGA, METABRIC, and GEO databases have limited numbers of patients and have substantial missing values. In the future, further prospective studies with larger sample volume are needed in order to overcome these issues.

4. Materials and Methods

4.1. The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) Patient Cohorts

The gene expression level quantification data (mRNA Z-score of RNA-sequence) of 1093 breast cancer patients, and level three Z-score normalized gene expression data of 1904 patients from METABRIC data were downloaded through cBioPortal as described previously [35,36,37,38]. Gene expression microarray datasets were obtained from ten Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/). Patients were classified as either AR expression high or low group using the median of AR mRNA expression of the individual group. Distant relapse-free survival was determined from time of diagnosis to clinical confirmation of tumor recurrence.

4.2. Gene Set Enrichment Analysis (GSEA)

GSEA was calculated by the Broad Institute software (http://software.broadinstitute.org/gsea/index.jsp), as previously described [35,36,37,39]. We classified the patients into two groups according to AR expression using the median of each group.

4.3. CIBERSORT, the Immune Cytolytic Activity (CYT) Score, and the Mutant-Allele Tumor Heterogeneity (MATH)

We utilized CIBERSORT, a bioinformatic algorithm, in order to calculate immune cell composition from gene expression profiles, in order to estimate tumor infiltrating immune cells in tumors [40]. Immune cell fraction data was downloaded through The Cancer Imaging Archive (TCIA) (https://tcia.at/home) [41]. The immune cytolytic activity (CYT) score was defined as the geometric mean of granzyme A (GZMA) and perforin (PRF1), which are pivotal cytolytic effectors, with expression values in transcripts per million [42]. The mutant-allele tumor heterogeneity (MATH) level was calculated from the median absolute deviation and the median of its mutant-allele fractions at tumor-specific mutated loci as described previously [43,44].

4.4. Statistical Analysis

The Kaplan–Meier method with log-rank test was used to determine distant relapse-free survival as previously described [35,36,37,45]. Fisher’s exact test was used to determine if there were nonrandom associations between two categorical variables. Student’s t-test was used to analyze the differences between continuous values. In all analyses, a two-sided p < 0.05 was considered statistically significant. All statistical analyses were performed using Microsoft Excel 2010, R software (http:///www.rproject.org/), and Bioconductor (http://bioconductor.org/).

5. Conclusions

Our results demonstrated that high AR has fewer tumor infiltrating lymphocytes and a worse response to NAC, but better survival in ER-positive BC. The fact that AR expression plays some cancer-immunological role in ER-positive breast cancer is enhanced by this report, which can add a novel aspect to future research and treatment of ER-positive breast cancer.

Author Contributions

M.O. (Maiko Okano) conceptualized study design, performed analysis and prepared the article. M.O. (Masanori Oshi), M.A., E.K., X.P., Q.Q., L.Y., performed analysis. A.L.B. prepared the article. K.T. provided supervision of analysis and prepared the article.

Funding

This work was supported by NIH grant R01CA160688 to K.T., and National Cancer Institute (NCI) grant P30CA016056 involving the use of Roswell Park Cancer Institute’s Resources.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BCBreast Cancer
EREstrogen Receptor
ARAndrogen Receptor
NACNeo-adjuvant Chemotherapy
TILTumor Infiltrating Lymphocytes
TCGAThe Cancer Genome Atlas
METABRICMolecular Taxonomy of Breast Cancer International Consortium
GSEAGene Set Enrichment Analysis
CYTCytolytic Activity
MATHMutant-Allele Tumor Heterogeneity

References

  1. SEER Stat Fact Sheets: Female Breast Cancer. Available online: https://seer.cancer.gov/statfacts/html/breast.html (accessed on 1 March 2019).
  2. Goldhirsch, A.; Wood, W.C.; Coates, A.S.; Gelber, R.D.; Thurlimann, B.; Senn, H.J. Strategies for subtypes—Dealing with the diversity of breast cancer: Highlights of the St. Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. Ann. Oncol. Off. J. Eur. Soc. Med. Oncol. 2011, 22, 1736–1747. [Google Scholar] [CrossRef] [PubMed]
  3. Fisher, B.; Redmond, C.; Fisher, E.R.; Caplan, R. Relative worth of estrogen or progesterone receptor and pathologic characteristics of differentiation as indicators of prognosis in node negative breast cancer patients: Findings from National Surgical Adjuvant Breast and Bowel Project Protocol B-06. Off. J. Am. Soc. Clin. Oncol. 1988, 6, 1076–1087. [Google Scholar] [CrossRef]
  4. Crowe, J.J.; Gordon, N.H.; Hubay, C.A.; Shenk, R.R.; Zollinger, R.M.; Brumberg, D.J.; McGuire, W.L.; Shuck, J.M. Estrogen receptor determination and long term survival of patients with carcinoma of the breast. Surg. Gynecol. Obstet. 1991, 173, 273–278. [Google Scholar] [PubMed]
  5. Denkert, C.; Loibl, S.; Noske, A.; Roller, M.; Muller, B.M.; Komor, M.; Budczies, J.; Darb-Esfahani, S.; Kronenwett, R.; Hanusch, C.; et al. Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer. Off. J. Am. Soc. Clin. Oncol. 2010, 28, 105–113. [Google Scholar] [CrossRef]
  6. Savas, P.; Salgado, R.; Denkert, C.; Sotiriou, C.; Darcy, P.K.; Smyth, M.J.; Loi, S. Clinical relevance of host immunity in breast cancer: From TILs to the clinic. Nat. Rev. Clin. Oncol. 2016, 13, 228–241. [Google Scholar] [CrossRef]
  7. Fasching, P.A.; Heusinger, K.; Haeberle, L.; Niklos, M.; Hein, A.; Bayer, C.M.; Rauh, C.; Schulz-Wendtland, R.; Bani, M.R.; Schrauder, M.; et al. Ki67, chemotherapy response, and prognosis in breast cancer patients receiving neoadjuvant treatment. BMC Cancer 2011, 11, 486. [Google Scholar] [CrossRef]
  8. Colleoni, M.; Viale, G.; Zahrieh, D.; Pruneri, G.; Gentilini, O.; Veronesi, P.; Gelber, R.D.; Curigliano, G.; Torrisi, R.; Luini, A.; et al. Chemotherapy is more effective in patients with breast cancer not expressing steroid hormone receptors: A study of preoperative treatment. Off. J. Am. Assoc. Cancer Res. 2004, 10, 6622–6628. [Google Scholar] [CrossRef]
  9. Denkert, C.; von Minckwitz, G.; Darb-Esfahani, S.; Lederer, B.; Heppner, B.I.; Weber, K.E.; Budczies, J.; Huober, J.; Klauschen, F.; Furlanetto, J.; et al. Tumour-infiltrating lymphocytes and prognosis in different subtypes of breast cancer: A pooled analysis of 3771 patients treated with neoadjuvant therapy. Lancet Oncol. 2018, 19, 40–50. [Google Scholar] [CrossRef]
  10. Miyoshi, Y.; Shien, T.; Ogiya, A.; Ishida, N.; Yamazaki, K.; Horii, R.; Horimoto, Y.; Masuda, N.; Yasojima, H.; Inao, T.; et al. Associations in tumor infiltrating lymphocytes between clinicopathological factors and clinical outcomes in estrogen receptor-positive/human epidermal growth factor receptor type 2 negative breast cancer. Oncol. Lett. 2019, 17, 2177–2186. [Google Scholar] [CrossRef] [PubMed]
  11. Moinfar, F.; Okcu, M.; Tsybrovskyy, O.; Regitnig, P.; Lax, S.F.; Weybora, W.; Ratschek, M.; Tavassoli, F.A.; Denk, H. Androgen receptors frequently are expressed in breast carcinomas: Potential relevance to new therapeutic strategies. Cancer 2003, 98, 703–711. [Google Scholar] [CrossRef]
  12. Peters, A.A.; Buchanan, G.; Ricciardelli, C.; Bianco-Miotto, T.; Centenera, M.M.; Harris, J.M.; Jindal, S.; Segara, D.; Jia, L.; Moore, N.L.; et al. Androgen receptor inhibits estrogen receptor-alpha activity and is prognostic in breast cancer. Cancer Res. 2009, 69, 6131–6140. [Google Scholar] [CrossRef]
  13. D’Amato, N.C.; Gordon, M.A.; Babbs, B.; Spoelstra, N.S.; Butterfield, K.T.C.; Torkko, K.C.; Phan, V.T.; Barton, V.N.; Rogers, T.J.; Sartorius, C.A.; et al. Cooperative Dynamics of AR and ER Activity in Breast Cancer. Mol. Cancer Res. MCR 2016, 14, 1054–1067. [Google Scholar] [CrossRef]
  14. Lehmann, B.D.; Bauer, J.A.; Chen, X.; Sanders, M.E.; Chakravarthy, A.B.; Shyr, Y.; Pietenpol, J.A. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J. Clin. Investig. 2011, 121, 2750–2767. [Google Scholar] [CrossRef] [Green Version]
  15. Huang, R.; Han, J.; Liang, X.; Sun, S.; Jiang, Y.; Xia, B.; Niu, M.; Li, D.; Zhang, J.; Wang, S.; et al. Androgen Receptor Expression and Bicalutamide Antagonize Androgen Receptor Inhibit beta-Catenin Transcription Complex in Estrogen Receptor-Negative Breast Cancer. Cell Physiol. Biochem. 2017, 43, 2212–2225. [Google Scholar] [CrossRef]
  16. Shen, Y.; Yang, F.; Zhang, W.; Song, W.; Liu, Y.; Guan, X. The Androgen Receptor Promotes Cellular Proliferation by Suppression of G-Protein Coupled Estrogen Receptor Signaling in Triple-Negative Breast Cancer. Cell Physiol. Biochem. 2017, 43, 2047–2061. [Google Scholar] [CrossRef]
  17. Migliaccio, A.; Di Domenico, M.; Castoria, G.; Nanayakkara, M.; Lombardi, M.; de Falco, A.; Bilancio, A.; Varricchio, L.; Ciociola, A.; Auricchio, F. Steroid receptor regulation of epidermal growth factor signaling through Src in breast and prostate cancer cells: Steroid antagonist action. Cancer Res. 2005, 65, 10585–10593. [Google Scholar] [CrossRef]
  18. Panet-Raymond, V.; Gottlieb, B.; Beitel, L.K.; Pinsky, L.; Trifiro, M.A. Interactions between androgen and estrogen receptors and the effects on their transactivational properties. Mol. Cell. Endocrinol. 2000, 167, 139–150. [Google Scholar] [CrossRef]
  19. Schwartzberg, L.S.; Yardley, D.A.; Elias, A.D.; Patel, M.; LoRusso, P.; Burris, H.A.; Gucalp, A.; Peterson, A.C.; Blaney, M.E.; Steinberg, J.L.; et al. A Phase I/Ib Study of Enzalutamide Alone and in Combination with Endocrine Therapies in Women with Advanced Breast Cancer. Off. J. Am. Assoc. Cancer Res. 2017, 23, 4046–4054. [Google Scholar] [CrossRef]
  20. Rechoum, Y.; Rovito, D.; Iacopetta, D.; Barone, I.; Andò, S.; Weigel, N.L.; O’Malley, B.W.; Brown, P.H.; Fuqua, S.A. AR collaborates with ERalpha in aromatase inhibitor-resistant breast cancer. Breast Cancer Res. Treat. 2014, 147, 473–485. [Google Scholar] [CrossRef]
  21. De Amicis, F.; Thirugnansampanthan, J.; Cui, Y.; Selever, J.; Beyer, A.; Parra, I.; Weigel, N.L.; Herynk, M.H.; Tsimelzon, A.; Lewis, M.T.; et al. Androgen receptor overexpression induces tamoxifen resistance in human breast cancer cells. Breast Cancer Res. Treat. 2010, 121, 1–11. [Google Scholar] [CrossRef]
  22. Castellano, I.; Allia, E.; Accortanzo, V.; Vandone, A.M.; Chiusa, L.; Arisio, R.; Durando, A.; Donadio, M.; Bussolati, G.; Coates, A.S.; et al. Androgen receptor expression is a significant prognostic factor in estrogen receptor positive breast cancers. Breast Cancer Res. Treat. 2010, 124, 607–617. [Google Scholar] [CrossRef] [Green Version]
  23. Bozovic-Spasojevic, I.; Zardavas, D.; Brohée, S.; Ameye, L.; Fumagalli, D.; Ades, F.; De Azambuja, E.; Bareche, Y.; Piccart, M.; Paesmans, M.; et al. The Prognostic Role of Androgen Receptor in Patients with Early-Stage Breast Cancer: A Meta-analysis of Clinical and Gene Expression Data. Clin. Cancer Res. 2017, 23, 2702–2712. [Google Scholar] [CrossRef]
  24. Gonzalez-Angulo, A.M.; Stemke-Hale, K.; Palla, S.L.; Carey, M.; Agarwal, R.; Meric-Berstam, F.; Traina, T.A.; Hudis, C.; Hortobagyi, G.N.; Gerald, W.L.; et al. Androgen receptor levels and association with PIK3CA mutations and prognosis in breast cancer. Clin. Cancer Res. 2009, 15, 2472–2478. [Google Scholar] [CrossRef]
  25. Witzel, I.; Graeser, M.; Karn, T.; Schmidt, M.; Wirtz, R.; Schütze, D.; Rausch, A.; Jänicke, F.; Milde-Langosch, K.; Müller, V. Androgen receptor expression is a predictive marker in chemotherapy-treated patients with endocrine receptor-positive primary breast cancers. J. Cancer Res. Clin. Oncol. 2013, 139, 809–816. [Google Scholar] [CrossRef]
  26. Vera-Badillo, F.E.; Templeton, A.J.; de Gouveia, P.; Diaz-Padilla, I.; Bedard, P.L.; Al-Mubarak, M.; Seruga, B.; Tannock, I.F.; Ocana, A.; Amir, E. Androgen receptor expression and outcomes in early breast cancer: A systematic review and meta-analysis. J. Natl. Cancer Inst. 2014, 106, 319. [Google Scholar] [CrossRef]
  27. Tsang, J.Y.; Ni, Y.B.; Chan, S.K.; Shao, M.M.; Law, B.K.; Tan, P.H.; Gary, M.T. Androgen receptor expression shows distinctive significance in ER-positive and negative breast cancers. Ann. Surg. Oncol. 2014, 21, 2218–2228. [Google Scholar] [CrossRef]
  28. Park, S.; Koo, J.S.; Kim, M.S.; Park, H.S.; Lee, J.S.; Lee, J.S.; Kim, S.I.; Park, B.W.; Lee, K.S. Androgen receptor expression is significantly associated with better outcomes in estrogen receptor-positive breast cancers. Ann. Oncol. 2011, 22, 1755–1762. [Google Scholar] [CrossRef]
  29. Polkinghorn, W.R.; Parker, J.S.; Lee, M.X.; Kass, E.M.; Spratt, D.E.; Iaquinta, P.J.; Arora, V.K.; Yen, W.F.; Cai, L.; Zheng, D.; et al. Androgen receptor signaling regulates DNA repair in prostate cancers. Cancer Discov. 2013, 3, 1245–1253. [Google Scholar] [CrossRef]
  30. Goodwin, J.F.; Schiewer, M.J.; Dean, J.L.; Schrecengost, R.S.; de Leeuw, R.; Han, S.; Ma, T.; Den, R.B.; Dicker, A.P.; Feng, F.Y.; et al. A hormone-DNA repair circuit governs the response to genotoxic insult. Cancer Discov. 2013, 3, 1254–1271. [Google Scholar] [CrossRef]
  31. Min, A.; Jang, H.; Kim, S.; Lee, K.H.; Kim, D.K.; Suh, K.J.; Yang, Y.; Elvin, P.; O’Connor, M.J.; Im, S.A. Androgen Receptor Inhibitor Enhances the Antitumor Effect of PARP Inhibitor in Breast Cancer Cells by Modulating DNA Damage Response. Mol. Cancer Ther. 2018, 17, 2507–2518. [Google Scholar] [CrossRef] [Green Version]
  32. Speers, C.; Zhao, S.G.; Chandler, B.; Liu, M.; Wilder-Romans, K.; Olsen, E.; Nyati, S.; Ritter, C.; Alluri, P.G.; Kothari, V.; et al. Androgen receptor as a mediator and biomarker of radioresistance in triple-negative breast cancer. NPJ Breast Cancer 2017, 3, 29. [Google Scholar] [CrossRef]
  33. van Rooijen, J.M.; Qiu, S.Q.; Timmer-Bosscha, H.; van der Vegt, B.; Boers, J.E.; Schröder, C.P.; de Vries, E.G. Androgen receptor expression inversely correlates with immune cell infiltration in human epidermal growth factor receptor 2-positive breast cancer. Eur. J. Cancer 2018, 103, 52–60. [Google Scholar] [CrossRef]
  34. Coulie, P.G.; Van den Eynde, B.J.; van der Bruggen, P.; Boon, T. Tumour antigens recognized by T lymphocytes: At the core of cancer immunotherapy. Nat. Rev. Cancer 2014, 14, 135–146. [Google Scholar] [CrossRef]
  35. Terakawa, T.; Katsuta, E.; Yan, L.; Turaga, N.; McDonald, K.A.; Fujisawa, M.; Guru, K.A.; Takabe, K. High expression of SLCO2B1 is associated with prostate cancer recurrence after radical prostatectomy. Oncotarget 2018, 9, 14207–14218. [Google Scholar] [CrossRef] [Green Version]
  36. Young, J.; Kawaguchi, T.; Yan, L.; Qi, Q.; Liu, S.; Takabe, K. Tamoxifen sensitivity-related microRNA-342 is a useful biomarker for breast cancer survival. Oncotarget 2017, 8, 99978–99989. [Google Scholar] [CrossRef] [Green Version]
  37. Kawaguchi, T.; Yan, L.; Qi, Q.; Peng, X.; Gabriel, E.M.; Young, J.; Liu, S.; Takabe, K. Overexpression of suppressive microRNAs, miR-30a and miR-200c are associated with improved survival of breast cancer patients. Sci. Rep. 2017, 7, 15945. [Google Scholar] [CrossRef]
  38. Kim, S.Y.; Kawaguchi, T.; Yan, L.; Young, J.; Qi, Q.; Takabe, K. Clinical Relevance of microRNA Expressions in Breast Cancer Validated Using the Cancer Genome Atlas (TCGA). Ann. Surg. Oncol. 2017, 24, 2943–2949. [Google Scholar] [CrossRef] [Green Version]
  39. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  40. Newman, A.M.; Liu, C.L.; Green, M.R.; Gentles, A.J.; Feng, W.; Xu, Y.; Hoang, C.D.; Diehn, M.; Alizadeh, A.A. Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 2015, 12, 453–457. [Google Scholar] [CrossRef] [Green Version]
  41. Charoentong, P.; Finotello, F.; Angelova, M.; Mayer, C.; Efremova, M.; Rieder, D.; Hackl, H.; Trajanoski, Z. Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cell Rep. 2017, 18, 248–262. [Google Scholar] [CrossRef] [Green Version]
  42. Narayanan, S.; Kawaguchi, T.; Yan, L.; Peng, X.; Qi, Q.; Takabe, K. Cytolytic Activity Score to Assess Anticancer Immunity in Colorectal Cancer. Ann. Surg. Oncol. 2018, 25, 2323–2331. [Google Scholar] [CrossRef]
  43. Mroz, E.A.; Rocco, J.W. MATH, a novel measure of intratumor genetic heterogeneity, is high in poor-outcome classes of head and neck squamous cell carcinoma. Oral Oncol. 2013, 49, 211–215. [Google Scholar] [CrossRef] [Green Version]
  44. Rocco, J.W. Mutant allele tumor heterogeneity (MATH) and head and neck squamous cell carcinoma. Head Neck Pathol. 2015, 9, 1–5. [Google Scholar] [CrossRef]
  45. Ramanathan, R.; Olex, A.L.; Dozmorov, M.; Bear, H.D.; Fernandez, L.J.; Takabe, K. Angiopoietin pathway gene expression associated with poor breast cancer survival. Breast Cancer Res. Treat. 2017, 162, 191–198. [Google Scholar] [CrossRef]
Ijms 20 02655 g001 550
Figure 1. Expression levels of AR mRNA in TCGA and METABRIC cohorts. AR mRNA levels were quantified by mRNA Z-scores of RNA-sequences. Bold lines demonstrate the median and box plot is 95% interval. * p < 0.001. (A) PAM50 subtype in TCGA; (B) ER positivity defined by IHC in TCGA; (C) PAM50 subtype in METABRIC; (D) ER positivity defined by IHC in METABRIC.
Figure 1. Expression levels of AR mRNA in TCGA and METABRIC cohorts. AR mRNA levels were quantified by mRNA Z-scores of RNA-sequences. Bold lines demonstrate the median and box plot is 95% interval. * p < 0.001. (A) PAM50 subtype in TCGA; (B) ER positivity defined by IHC in TCGA; (C) PAM50 subtype in METABRIC; (D) ER positivity defined by IHC in METABRIC.
Ijms 20 02655 g001
Ijms 20 02655 g002 550
Figure 2. Pearson correlation coefficient (r) and p-value (p) between AR mRNA expression and ER mRNA expression. (A) A scatter plot of TCGA cohort; (B) A scatter plot of METABRIC cohort.
Figure 2. Pearson correlation coefficient (r) and p-value (p) between AR mRNA expression and ER mRNA expression. (A) A scatter plot of TCGA cohort; (B) A scatter plot of METABRIC cohort.
Ijms 20 02655 g002
Ijms 20 02655 g003 550
Figure 3. Gene set enrichment analysis (GSEA) demonstrated significant enrichment of the following hallmark gene sets to AR mRNA expression in ER-positive BC in TCGA. (A) UV response down gene set enriched in high expression of AR mRNA; (B) DNA repair gene set enriched in low expression of AR mRNA.
Figure 3. Gene set enrichment analysis (GSEA) demonstrated significant enrichment of the following hallmark gene sets to AR mRNA expression in ER-positive BC in TCGA. (A) UV response down gene set enriched in high expression of AR mRNA; (B) DNA repair gene set enriched in low expression of AR mRNA.
Ijms 20 02655 g003
Ijms 20 02655 g004 550
Figure 4. In CIBERSORT analysis, activated memory CD4 T cells and Gamma Delta T cells in AR high tumors are lower than those in AR low tumors in ER-positive BC.
Figure 4. In CIBERSORT analysis, activated memory CD4 T cells and Gamma Delta T cells in AR high tumors are lower than those in AR low tumors in ER-positive BC.
Ijms 20 02655 g004
Ijms 20 02655 g005 550
Figure 5. Tumor heterogeneity and cytolytic activity score in AR high BC is lower than that in AR low ER-positive BC.
Figure 5. Tumor heterogeneity and cytolytic activity score in AR high BC is lower than that in AR low ER-positive BC.
Ijms 20 02655 g005
Ijms 20 02655 g006 550
Figure 6. pCR rate and Kaplan–Meier survival curve by AR mRNA levels of GEO cohorts. High and low expressions are represented by the red and blue lines, respectively.
Figure 6. pCR rate and Kaplan–Meier survival curve by AR mRNA levels of GEO cohorts. High and low expressions are represented by the red and blue lines, respectively.
Ijms 20 02655 g006
Table
Table 1. Significant gene sets that associated with high expression of AR using Gene Set Enrichment Analysis (GSEA)(TCGA).
Table 1. Significant gene sets that associated with high expression of AR using Gene Set Enrichment Analysis (GSEA)(TCGA).
NameSizeESNESNOM p-valFDR q-valFWER p-val
HALLMARK_UV_RESPONSE_DN1400.55212.01130.00000.01260.0210
1 ES, enrichment score; 2 NES, normalized enrichment score; 3 NOM-p-val, normalized p-value; 4 FDR q-val, false discovery rate q-value; 5 FWER p-val, family-wise error rate p-value.
Table
Table 2. Significant gene sets that associated with low expression of AR using Gene Set Enrichment Analysis (GSEA)(TCGA).
Table 2. Significant gene sets that associated with low expression of AR using Gene Set Enrichment Analysis (GSEA)(TCGA).
NameSizeESNESNOM p-valFDR q-valFWER p-val
HALLMARK_DNA_REPAIR142−0.5265−2.02980.00000.02810.0150
1 ES, enrichment score; 2 NES, normalized enrichment score; 3 NOM-p-val, normalized p-value; 4 FDR q-val, false discovery rate q-value; 5 FWER p-val, family-wise error rate p-value.

Share and Cite

MDPI and ACS Style

Okano, M.; Oshi, M.; Butash, A.L.; Asaoka, M.; Katsuta, E.; Peng, X.; Qi, Q.; Yan, L.; Takabe, K. Estrogen Receptor Positive Breast Cancer with High Expression of Androgen Receptor has Less Cytolytic Activity and Worse Response to Neoadjuvant Chemotherapy but Better Survival. Int. J. Mol. Sci. 2019, 20, 2655. https://doi.org/10.3390/ijms20112655

AMA Style

Okano M, Oshi M, Butash AL, Asaoka M, Katsuta E, Peng X, Qi Q, Yan L, Takabe K. Estrogen Receptor Positive Breast Cancer with High Expression of Androgen Receptor has Less Cytolytic Activity and Worse Response to Neoadjuvant Chemotherapy but Better Survival. International Journal of Molecular Sciences. 2019; 20(11):2655. https://doi.org/10.3390/ijms20112655

Chicago/Turabian Style

Okano, Maiko, Masanori Oshi, Ali Linsk Butash, Mariko Asaoka, Eriko Katsuta, Xuan Peng, Qianya Qi, Li Yan, and Kazuaki Takabe. 2019. "Estrogen Receptor Positive Breast Cancer with High Expression of Androgen Receptor has Less Cytolytic Activity and Worse Response to Neoadjuvant Chemotherapy but Better Survival" International Journal of Molecular Sciences 20, no. 11: 2655. https://doi.org/10.3390/ijms20112655

APA Style

Okano, M., Oshi, M., Butash, A. L., Asaoka, M., Katsuta, E., Peng, X., Qi, Q., Yan, L., & Takabe, K. (2019). Estrogen Receptor Positive Breast Cancer with High Expression of Androgen Receptor has Less Cytolytic Activity and Worse Response to Neoadjuvant Chemotherapy but Better Survival. International Journal of Molecular Sciences, 20(11), 2655. https://doi.org/10.3390/ijms20112655

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop