Abstract
Androgen receptor (AR) is expressed in 90% of estrogen receptor alpha–positive (ER+) breast tumors, but its role in tumor growth and progression remains controversial. Use of two anti-androgens that inhibit AR nuclear localization, enzalutamide and MJC13, revealed that AR is required for maximum ER genomic binding. Here, a novel global examination of AR chromatin binding found that estradiol induced AR binding at unique sites compared with dihydrotestosterone (DHT). Estradiol-induced AR-binding sites were enriched for estrogen response elements and had significant overlap with ER-binding sites. Furthermore, AR inhibition reduced baseline and estradiol-mediated proliferation in multiple ER+/AR+ breast cancer cell lines, and synergized with tamoxifen and fulvestrant. In vivo, enzalutamide significantly reduced viability of tamoxifen-resistant MCF7 xenograft tumors and an ER+/AR+ patient-derived model. Enzalutamide also reduced metastatic burden following cardiac injection. Finally, in a comparison of ER+/AR+ primary tumors versus patient-matched local recurrences or distant metastases, AR expression was often maintained even when ER was reduced or absent. These data provide preclinical evidence that anti-androgens that inhibit AR nuclear localization affect both AR and ER, and are effective in combination with current breast cancer therapies. In addition, single-agent efficacy may be possible in tumors resistant to traditional endocrine therapy, as clinical specimens of recurrent disease demonstrate AR expression in tumors with absent or refractory ER.
Implications: This study suggests that AR plays a previously unrecognized role in supporting E2-mediated ER activity in ER+/AR+ breast cancer cells, and that enzalutamide may be an effective therapeutic in ER+/AR+ breast cancers. Mol Cancer Res; 14(11); 1054–67. ©2016 AACR.
This article is featured in Highlights of This Issue, p. 1031
Introduction
AR is more frequently expressed in breast cancer than estrogen receptor alpha (ER) or progesterone receptor (PR) (1); however, the role of AR is complex, dependent on the hormonal milieu, and remains controversial. AR positivity is associated with better prognosis in ER+ breast cancer (2–4), possibly due to the fact that like ER, AR positivity is indicative of a more well-differentiated state. In the presence of estradiol (E2), the androgen dihydrotestosterone (DHT) decreased E2-induced proliferation (2) and ER transcriptional activity (5), leading to the conclusion that AR is protective in breast cancer. However, there is accumulating evidence that androgen signaling and AR are involved in resistance to ER-directed endocrine therapies. De novo or acquired resistance to anti-estrogen therapies is a frequent occurrence, and ultimately all metastatic ER+ breast cancers are resistant (6, 7). In ER+ tumors responsive to neoadjuvant aromatase inhibitor (AI) therapy, AR mRNA and nuclear AR protein decreased, whereas in nonresponsive tumors it remained elevated (8, 9). AR overexpression in breast cancer cell lines resulted in resistance to tamoxifen and AIs in vitro and in vivo (10, 11). One mechanism of resistance to anti-estrogen therapies may therefore be tumor adaptation from estrogen to androgen dependence.
AIs block the conversion of androgens to estrogens, and free testosterone and dehydroepiandrosterone sulfate (DHEA-S) increased in patients on AIs (12). Furthermore, high levels of the adrenal androgen DHEA-S are predictive of failure on AIs, and circulating DHEA-S increased during treatment in patients with tumors that progressed during AI treatment (13). Patients with tumors exhibiting a high ratio of percent cells positive for AR versus ER protein are more likely to have recurrent disease while on tamoxifen and also have a worse overall prognosis compared with those with a more equal ratio of these two receptors, as is found in normal breast epithelium (14). So although AR, like ER, is associated with a better prognosis, anti-androgen therapies may benefit patients with AR+ breast cancers if the tumors are dependent on activated AR.
We previously reported that the new generation AR antagonist enzalutamide, which inhibits AR nuclear localization, decreased estrogen-induced tumor growth, while the first-generation AR antagonist bicalutamide did not (14). However, the mechanism by which enzalutamide affected ER activity was not known. Herein, we demonstrate for the first time that in response to E2, nuclear localization of AR supports maximum ER genomic binding, and that AR inhibition with the pure antagonist enzalutamide significantly decreases E2-induced growth of ER+/AR+ cell lines and patient-derived xenografts, as well as tamoxifen-resistant tumors in vivo, and also decreases metastatic burden. Importantly, these data suggest that patients with ER+/AR+ breast cancer may benefit from combining anti-androgen therapy with anti-estrogen therapy, and that tumors resistant to traditional ER-directed therapies may be responsive to AR-directed drugs.
Materials and Methods
Cell lines
All cell lines were authenticated by short tandem repeat analysis using AmpFLSTR Identifiler PCR Amplification Kit (Life Technologies) and tested negative for mycoplasma in January 2015. MCF7 cells were obtained from Dr. Kate Horwitz at the University of Colorado Anschutz Medical Campus. MCF7-TamR cells obtained from Dr. Doug Yee at the University of Minnesota were generated by chronic treatment of MCF7 cells with 100 nmol/L tamoxifen. All other cell lines were obtained from the ATCC. Additional cell culture details are included in Supplementary Material. The BCK4 cell line is an ER+/AR+ breast cancer line recently derived from a pleural effusion (15), and the PT12 breast cancer cell line is ER+/AR+ and created from a patient-derived xenograft (PDX; ref. 16). Originally the PT12 PDX was described as AR-negative (16), but upon staining of the original passage with a more sensitive AR antibody (SP107 from Cell Marque), the PDX was found to be AR+ (Supplementary Fig. S6A).
Cellular assays and reagents
Cells were treated with 10 nmol/L estradiol (E2, Sigma Aldrich) and 10 nmol/L dihydrotestosterone (DHT, Sigma Aldrich). Androgen concentrations have been previously examined in breast cancer (17) and intratumoral DHT concentrations (249 pg/g) were significantly higher than in blood. The DHT concentration of the current study is consistent with other in vitro studies of DHT in breast cancer (18, 19), and approximates levels of circulating testosterone in obese, postmenopausal women (12) as well as DHT levels in FBS used during routine tissue culture. 10 μmol/L enzalutamide (Medivation) approximates the IC50 of the three cell lines studied and is a clinically achievable, well-tolerated treatment concentration (NCT01889238).
Proliferation assays
Proliferation assays were performed using the IncuCyte ZOOM live cell imaging system (Essen BioSciences) or crystal violet as described previously (20). For synergy experiments, percent inhibition was calculated compared with vehicle control, and the combination index was calculated for each dose combination by CalcuSyn (21) (Biosoft). Soft agar assays were performed in 6-well plates using 0.5% bottom and 0.25% top layer agar (Difco Agar Noble, BD Biosciences). Wells were photographed and colony number and size was determined by ImageJ software (NIH, Bethesda, MD).
Tumor studies
Xenograft experiments were approved by the University of Colorado Institutional Animal Care and Use Committee [IACUC protocol 83614(01)1E] and were conducted in accordance with the NIH Guidelines of Care and Use of Laboratory Animals. A total of 1 × 106 MCF7-GFP-Luc cells were mixed with growth factor–reduced Matrigel (BD Biosciences) and injected bilaterally into the mammary fat pad of female ovariectomized athymic nu/nu mice (Taconic). E2 pellets (60-day release, 1.5 mg/pellet, Innovative Research of America) were implanted subcutaneously (SQ) at the back of the neck. Once tumors were established, mice were randomized into groups based on total tumor burden as measured by in vivo imaging. Mice received enzalutamide in their chow (∼50 mg/kg daily dose). Enzalutamide was mixed with ground mouse chow (Research Diets Inc.) at 0.43 mg/g chow. Control mice received the same chow without enzalutamide. All mice were given free access to enzalutamide-formulated chow or control chow (CTRL chow) during the study. Mice were euthanized by carbon dioxide asphyxiation followed by cervical dislocation, and tumors were harvested. The MCF7-TamR xenograft experiment was performed as described above without estrogen pellets. For the PT-12 xenograft study, 1 × 106 cells were injected bilaterally into the mammary fat pad of NOD-SCID-IL2Rgc−/− female mice. Mice were implanted with a DHT (8 mg) or E2 (2 mg) pellet. For the metastasis experiment, 2.5 × 105 GFP-Luciferase labeled PT12 cells were injected intracardially in NOD-SCID-IL2Rgc−/− mice implanted with E2 pellets (2 mg). PT12 experiments with DHT were performed in ovariectomized females, while PT12 experiments with E2 were performed in nonovariectomized females as the E2 pellet overrides the estrus cycle.
Immunoblotting
Whole-cell protein extracts (50 μg) were denatured, separated on SDS-PAGE gels and transferred to PVDF membranes. After blocking in 3% BSA in TBS-T, membranes were probed overnight at 4°C. Primary antibodies used were: ERalpha (Neomarkers Ab-16, 1:500 dilution), AR (EMD Millipore PG-21, 1:500 dilution), Topo 1 (Santa Cruz Biotechnology C-21, 1:100 dilution) and α-tubulin (clone B-5-1-2 from Sigma, 1:30,000 dilution). After incubation with appropriate secondary antibody, results were detected using Western Lightning Chemiluminescence Reagent Plus (Perkin Elmer).
Nuclear–cytoplasmic fractionation
A total of 1 × 106 cells were seeded in 10-cm dishes in medium supplemented with 5% charcoal-stripped serum (CSS). After 3 days, the cells were pretreated with vehicle or 10 μmol/L enzalutamide for 3 hours and then cotreated with either DHT for 3 hours plus or minus enzalutamide, or E2 for 1 hour plus or minus enzalutamide. Cells were washed with PBS and cellular fractionation was performed using the NE-PER Nuclear and Cytoplasmic Extraction Kit (Life Technologies) as per manufacturer's instructions.
Proximity ligation assay
PLA was performed using the Duolink kit according to manufacturer's instructions (Olink Bioscience). Briefly, 1.5 × 104 cells were plated in 8-well chamber slides and hormone starved in phenol red–free media with 5% CSS for 72 hours. Cells were then pretreated with vehicle or enzalutamide for 3 hours, then treated with hormones ± enzalutamide as described for 1 hour. After fixation with 4% paraformaldehyde, cells were permeabilized then blocked. Samples were then incubated with primary antibodies AR D6F11 (Cell Signaling Technology), and ERα clone 6F11 (Vector Laboratories) overnight at 4°C. Samples were then incubated with secondary antibodies linked to PLA probes and ligase was added. Detection reagent red was added and DAPI mounting media was added to visualize nuclei. Images were captured using a 20× objective. DAPI-labeled nuclei and red ERα/AR complexes were quantified using CellProfiler (22).
ChIP-seq
ChIP-seq was performed by Active Motif. Briefly, 1 × 106 MCF7 cells were seeded in 15-cm dishes in phenol red–free medium supplemented with 5% CSS for 72 hours. Cells were pretreated with vehicle, 10 μmol/L enzalutamide, or 30 μmol/L MJC13 for 3 hours. E2 was then added for 1 hour in continued presence of vehicle, enzalutamide, or MJC13. For AR ChIP-seq, an additional sample was treated with DHT for 4 hours. The cells were washed with PBS then fixed as per the manufacturer's instructions (Active Motif). AR antibody H-280 (Santa Cruz Biotechnology) or ER antibody HC-20 (Santa Cruz Biotechnology) were utilized. Peak calls were made by MACS2 (23) with default parameters using the sequence alignments obtained from Active Motif. Motif discovery was performed on 100 base pairs surrounding the peak summit using BioProspector (24). Patser (25) was used to determine significant matches to AREs and EREs.
RNA-seq
RNA Libraries were constructed using Illumina TruSEQ stranded mRNA Sample Prep Kit (cat# RS-122-2101). Total RNA was combined with RNA purification beads to bind PolyA RNA to oligodT magnetic beads. mRNA was eluted and converted to double stranded DNA. A Tailing, adapter ligation, and PCR amplification using 15 cycles was used to complete the library construction. Libraries were quantitated via Qubit, analyzed on a Bioanalyzer Tape Station and diluted to appropriate concentration to run on an Illumina HiSEQ 2500 High Throughput Flow Cell. Reads were mapped to the human genome (hg19) by gSNAP, expression (FPKM) derived by Cufflinks, and differential expression analyzed with ANOVA in R (26, 27).
Statistical analyses
For most analyses, statistical significance was evaluated using a two-tailed Student t test or ANOVA with Bonferroni or Dunnett multiple comparisons test or nonparametric equivalents in GraphPad Prism (Ver 6, GraphPad Software) or SAS (ver 9.4, SAS Institute). Test assumptions were checked for all analyses. If data distributions were skewed, data transformations were attempted to allow the use of parametric tests. If data transformations failed, than a nonparametric test was used. For the PT12 xenograft experiments, due to unequal time measurements, the repeated measures mixed model approach was used rather than a standard repeated-measures ANOVA. The data met the assumption of normality (Shapiro–Wilk test P > 0.05 and frequency distribution graphs were symmetrical without evidence of outliers). For the PT12 cardiac injection experiment, we were not able to use a repeated measures approach as the data did not meet the assumptions of a normal distribution despite different data transformations. Therefore, a single Wilcoxon rank-sum test was used to determine difference between E2 and E2 plus enzalutamide at week 12. P ≤ 0.05 was considered statistically significant, with P values indicated in figures as*, P ≤ 0.05; **, P ≤ 0.01; ***, P ≤ 0.001; ****, P ≤ 0.0001. Error bars represent SEM unless otherwise noted.
Results
AR inhibition impairs ER+/AR+ breast cancer cell proliferation
The role of AR in ER+/AR+ breast cancer remains controversial, with conflicting data suggesting either proliferative or protective effects on breast cancer cells in vitro (2, 5, 28–31). Most studies of AR function in breast cancer have focused on the effect of androgen stimulation in the presence of E2 in hormone-depleted media. However, we believe it is more relevant to study the effects of activating or inhibiting AR either (i) in the absence of E2, to model postmenopausal women with breast cancer treated with AIs, or (ii) in full serum, which contains androgens as well as sufficient estrogens to induce ER activity and genomic binding (32).
Enzalutamide, which inhibits AR nuclear translocation and DNA binding (14, 33), significantly decreased growth of MCF7 cells grown in full serum (Fig. 1A) as well as two additional ER+/AR+ cell lines, T47D and ZR-75-1, (Supplementary Fig. S1A) in a concentration-dependent manner. This shows that AR activity is necessary for ER+/AR+ cell growth under typical culture conditions. Enzalutamide also decreased colony size of MCF7 cells (Fig. 1B) and T47D cells (Supplementary Fig. S1B) grown in soft agar using complete culture media, similar in magnitude to the effect of the anti-estrogen tamoxifen. We next decreased AR expression in MCF7 cells using two different shRNA constructs, and AR protein was confirmed to be decreased by Western blot analysis (Fig. 1C). AR knockdown led to a significant decrease in MCF7 cell growth over the course of 7 days (Fig. 1C), further demonstrating that AR is required for baseline proliferation of ER+/AR+ breast cancer cells in hormone-replete conditions.
AR inhibition decreases ER+/AR+ breast cancer growth. A, proliferation of MCF7 cells treated with increasing concentrations of enzalutamide (Enza) was monitored by IncuCyte. B, MCF7 cells were grown in soft agar with enzalutamide or tamoxifen and colony size was measured by ImageJ. C, immunoblotting for AR, ER, and Tubulin in MCF7 cells expressing a nontargeting (shNeg) or AR-targeting (shAR15 and shAR17) shRNA constructs (top). Proliferation was monitored by IncuCyte (bottom). D, MCF7 cells were grown in media with CSS for 72 hours then treated with vehicle (Veh), E2, or E2 + enzalutamide or MJC13 and cell number was measured by crystal violet. E, MCF7 cells were grown in media with CSS for 72 hours then treated with Veh, E2, or E2 + enzalutamide for 24 hours followed by cell-cycle analysis. F, MCF7 cells expressing shNeg, shAR15, or shAR17 were cultured in media with CSS for 72 hours then treated with vehicle or E2 and growth was measured by crystal violet. Error bars, SEM. *, P < 0.05; ****, P < 0.0001 by ANOVA with Dunnett multiple comparison test.
New-generation AR inhibitors decrease E2-induced proliferation
We previously showed that enzalutamide, which does not bind to ER by ligand binding assay, inhibits E2-induced growth of ER+/AR+ breast cancer cells in vitro and in vivo (14). To demonstrate that this is AR-dependent and not specific to enzalutamide, we also utilized MJC13, which inhibits AR nuclear localization by targeting ligand-induced dissociation of AR from FKBP52 in the cytosol (34). Both enzalutamide and MJC13 inhibited E2-induced proliferation in MCF7 cells (Fig. 1D). Enzalutamide decreased E2-induced growth in a concentration-dependent manner in additional luminal cell lines T47D and ZR-75-1, as well as in ER+/AR+ PT12 cells recently created from a patient-derived xenograft (ref. 16; Supplementary Fig. S1C). EC50 values for enzalutamide-mediated inhibition of E2-induced growth in MCF7 and T47D cells were determined to be 19.0 μmol/L and 17.1 μmol/L, respectively (Supplementary Fig. S1D), which are concentrations readily achieved in patients.
To specifically determine whether enzalutamide affected E2-induced proliferation, cell-cycle analysis of E2-treated MCF7 and T47D cells was performed. Enzalutamide significantly decreased the percent cells in S and G2–M phases compared with E2 treatment alone (Fig. 1E and Supplementary Fig. S1E). Silencing AR using shRNA also significantly decreased E2-induced proliferation of MCF7 cells compared with cells transduced with nontargeting shRNA (Fig. 1F). Together, these data confirm that pharmacologic AR inhibition or AR knockdown similarly diminish E2-driven proliferation of ER+/AR+ breast cancer cells.
AR inhibitors diminish ER genome binding
AR is capable of interacting with ER and estrogen response elements (EREs; refs. 2, 5, 35), thus, we postulated that inhibitors of AR nuclear localization might diminish baseline and E2-induced growth by altering ER genomic binding. To test this hypothesis, MCF7 cells were pretreated for 3 hours with vehicle, enzalutamide, or MJC13 then treated with E2 1 hour and global ER ChIP-seq was performed. Surprisingly, the anti-androgens enzalutamide or MJC13 dramatically decreased E2-induced ER genomic binding (Fig. 2A). The majority of sites displayed an approximate 50% decrease in ER binding (Fig. 2A–C), with no appreciable shift in the location of ER-binding sites upon enzalutamide or MJC13 treatment. The decrease in ER-binding intensity by enzalutamide or MJC13 was confirmed by qPCR after ChIP at previously characterized ER-binding sites including GREB1, GATA3, and PGR (Fig. 2D and E). Together, this suggests that the interaction of AR and ER is necessary for efficient ER genomic binding in response to E2, and that inhibition of nuclear AR localization decreases E2-induced ER activity by diminishing ER genome binding.
AR inhibitors diminish ER genomic binding. ChIP-seq for ER in MCF7 cells grown in CSS for 3 days then treated with E2 ± enzalutamide or MJC-13. A, heatmap of ER binding. The heatmap is shown with a horizontal window of ± 2 kb. B, the number of binding sites identified by MACS2, using vehicle treatment as the control. C, the ER ChIP-seq signal at individual sites with E2 alone (x-axis) versus E2 + enzalutamide (blue) or MJC13 (red; y-axis). D and E, ChIP-qPCR (D) and ChIP-seq read depth (E) at well-characterized ER-binding sites. Error bars, SEM. **, P < 0.01; ***, P < 0.001; ****, P < 0.0001 by ANOVA with Dunnett multiple comparison test.
Enzalutamide decreases nuclear localization of both AR and ER
As E2-induced ER genome binding was globally decreased by anti-androgens, we speculated that ER nuclear localization in response to E2 might be affected. Immunofluorescent staining of MCF7 cells grown in CSS revealed nuclear localization of both ER and AR following E2 treatment (Fig. 3A). Notably, treatment with enzalutamide decreased nuclear localization of both receptors (Fig. 3A), while bicalutamide did not, suggesting that the mechanism by which enzalutamide globally inhibits ER genomic binding may be by decreasing ER nuclear localization. We previously showed that enzalutamide decreased E2-driven growth of MCF7 xenograft tumors equally as well as tamoxifen (14). IHC for ER performed on tumors from mice on enzalutamide -containing chow had significantly decreased nuclear localization of ER compared with tumors from mice on CTRL chow (14).
AR and ER colocalize in the nucleus in response to E2. A, MCF7 cells were grown in media with CSS for 72 hours then pretreated with vehicle (veh), 10 μmol/L enzalutamide, or 1 μmol/L bicalutamide (bic). Following pretreatment, cells were treated with vehicle or 10 nmol/L E2 ± enzalutamide or bicalutamide as shown for an additional 3 hours. Cells were then fixed and ICC was performed for AR (green) and ER (red). B and C, MCF7 cells were grown in media with CSS for 72 hours then treated with the indicated treatment for 3 hours, and nuclear extracts were immunoblotted for AR and TOPO1. D, ER+/AR+ ZR-75-1 or (E) ER−/AR+ MDA-453 cells were grown in media with CSS for 72 hours, then pretreated for 3 hours with enzalutamide or vehicle control. Following pretreatment, cells were treated with vehicle, 10 nmol/L DHT, or 10 nmol/L E2 ± enzalutamide as shown for 3 additional hours. Nuclear extracts were then obtained and subjected to Western blotting for AR and TopoI. F, MCF7 cells were grown in media with CSS for 72 hours then treated with E2 ± enzalutamide for 1 hour followed by fixation and PLA staining for AR and ER (red). Nuclei were stained with DAPI (blue). G, fluorescent intensity per nuclei was measured by CellProfiler. Error bars, SEM. ****, P < 0.0001 by ANOVA with Dunnett multiple comparison test.
To further examine AR nuclear localization in response to E2, MCF7 cells were treated with E2 or DHT plus or minus enzalutamide for 3 hours, and nuclear and cytoplasmic protein fractions were isolated. DHT induced a strong increase in AR nuclear localization as expected, which was largely blocked by cotreatment with enzalutamide, but not bicalutamide (Fig. 3B). E2 treatment also increased AR nuclear localization, and this effect was blocked by enzalutamide, but not bicalutamide (Fig. 3C). E2-induced nuclear localization of AR was also observed in ZR-75-1 (Fig. 3D). However, E2 did not induce AR nuclear localization in ER−/AR+ MDA-MB-453 cells (Fig. 3E) or MDA-MB-231 cells (Supplementary Fig. S2A), suggesting that the observed AR nuclear localization is not due to promiscuous binding of E2 to AR, but rather that AR becomes localized to the nucleus in an ER-dependent manner upon E2 stimulation in ER+/AR+ breast cancer cells.
As both ER and AR were localized to the nucleus following E2 treatment, we next tested whether E2 induced AR and ER to colocalize using the proximity ligation assay (PLA). MCF7 cells treated with 10 nmol/L E2 for 1 hour demonstrated a strong increase in PLA signal when probed for ER and AR compared with vehicle control or enzalutamide treatment alone. This E2-induced increase in PLA signal was dramatically inhibited by pretreatment with enzalutamide (Fig. 3F and G). Similar results were observed in T47D cells (Supplementary Fig. S2B–S2D), suggesting that AR colocalizes with ER in the nucleus in response to E2.
E2 induces AR DNA binding distinct from DHT
To examine whether the observed nuclear localization of AR in response to E2 was associated with AR genome binding, hormone-deprived MCF7 cells were treated with DHT or E2 followed by global AR ChIP-seq analysis. As expected, DHT treatment induced a significant increase in AR genome binding compared with vehicle treatment (Fig. 4A and B). Among the 1,813 DHT-induced AR-binding sites identified in MCF7 cells, 49% were previously identified as bound by AR in LNCaP, a prostate cancer cell line, while 73.6% were bound by AR in MDA-MB-453, an ER−/AR+ breast cancer cell line (ref. 36; Supplementary Fig. S3A). This indicated that DHT-induced AR binding may be more similar between luminal breast cancer cell lines than between breast and prostate cancer cell lines, and is similar to previously reported findings in ZR-75-1 cells (35).
E2 induces AR genome binding that overlaps with ER binding. ChIP-seq for AR in MCF7 cells grown in CSS for 3 days then treated with E2 for 1 hour or DHT for 4 hours. A, heatmap of binding showing a horizontal window of ± 2 kb and enriched motifs from each category. B, the number of binding sites identified by MACS2, using vehicle treatment as the control. C, the number of AR-binding sites that are unique to DHT (red), unique to E2 (blue), or shared (overlap) are shown. D and E, ChIP-qPCR (D) and ChIP-seq read depth (E) results show AR binding at well-characterized ER-binding sites following E2 treatment. F, the percentage of AR-binding sites in response to DHT (left) or E2 (right) that were also identified as ER-binding sites (blue) is shown.
Surprisingly, E2 also induced AR genome binding, with 1,380 AR binding events identified in E2-treated MCF7 cells (Fig. 4A and B). Enzalutamide abolished E2-induced AR genomic binding, consistent with inhibition of AR nuclear localization and previously published reports in prostate cancer (33). Only 25% of all AR-bound sites overlapped between the two hormone treatments, indicating a large shift in AR genomic binding between DHT and E2 (Fig. 4C). For example, qPCR after ChIP demonstrated that DHT, but not E2, induced a robust increase in AR binding at previously characterized AR targets FKBP5 and ZBTB16 (Fig. 4D). Both E2 and DHT treatments resulted in AR binding to previously characterized ER targets GREB1 and GATA3, but only E2 treatment resulted in AR binding at a different ER target, PGR (Fig. 4D and E).
The most highly enriched motif among AR-binding sites in response to DHT was a FOXA1 motif (Fig. 4A), consistent with previous studies demonstrating strong overlap between AR and FOXA1-binding sites in breast cancer cells (36). However, the most highly enriched motif among AR binding sites unique to E2 treatment was a slightly degenerate estrogen response element (ERE; Fig. 4A), suggesting that AR was bound within 200 bp of ER-binding sites in the presence of E2. Indeed, full palindromic EREs were highly enriched among these sites, compared with sites bound by AR in response to DHT (Supplementary Fig. S3B and S3C). Validated nuclear ERα network was the most highly enriched pathway among genes near AR-binding sites unique to E2 treatment, whereas this network was not enriched among genes near AR-binding sites in response to DHT (Supplementary Table S1). Thus, in response to E2, AR binds to many sites correlated with ER regulation.
Finally we compared AR and ER binding following E2 treatment and found that 75% of E2-induced AR-binding sites overlapped with ER-binding sites (Fig. 4F). Notably, ER genome binding was more strongly inhibited by enzalutamide or MJC13 at these overlapping sites compared with nonoverlapping sites (Supplementary Fig. S3D and S3E), suggesting that AR might be facilitating ER binding at these loci. Taken together, these data demonstrate that in response to E2, AR and ER bind a significant number of overlapping loci and suggest that new generation antiandrogens which inhibit AR nuclear localization decrease ER activity and E2-mediated tumor growth by diminishing ER genome binding.
Enzalutamide synergizes with anti-estrogens
Because enzalutamide inhibited baseline and E2-induced growth by a different mechanism than currently used anti-estrogens, we hypothesized that it might act synergistically with anti-estrogens such as tamoxifen or fulvestrant in ER+/AR+ breast cancer cells. T47D cells were treated with varying concentrations of enzalutamide and/or tamoxifen, and all combinations showed synergistic inhibition of E2-induced growth as determined by CalcuSyn (Fig. 5A). Enzalutamide and tamoxifen also showed synergy or additive inhibition of E2-induced proliferation in MCF7 cells (Supplementary Fig. S4A), and the combination of enzalutamide plus tamoxifen reduced MCF7 growth in soft agar more significantly than either drug alone (Supplementary Fig. S4B). We also tested for synergy between enzalutamide and fulvestrant. In BCK4 cells, these drugs showed synergy in 10 of 15 dose combinations (Fig. 5B), with similar results also observed in PT12 cells (Fig. 5C) and ZR-75-1 cells (Supplementary Fig. S4C). Together, this shows that enzalutamide effectively synergizes with anti-estrogens to inhibit both baseline and E2-induced growth of ER+/AR+ cells, likely due to the ability of enzalutamide to inhibit AR as well as to indirectly inhibit ER.
Enzalutamide (Enza) synergizes with tamoxifen (tam) and fulvestrant in vitro. A, T47D cells were grown in media with complete serum and enzalutamide and/or tamoxifen, and cell number was monitored by IncuCyte. Percent inhibition was compared to vehicle after 5 days, and synergy was calculated using CalcuSyn software. A combination index (CI) value < 1 is indicative of synergy. B and C, BCK4 or PT12 cells were grown in phenol red–free media with CSS for 1 day then treated with E2 with enzalutamide and/or fulvestrant and cell number was monitored by IncuCyte.
Enzalutamide inhibits tamoxifen-resistant tumor growth
Resistance to currently-used endocrine therapies is a common occurrence facing ER+ breast cancer patients. Therefore, we also tested whether enzalutamide could inhibit growth of tamoxifen-resistant MCF7 (MCF7-TamR) cells (37). In vitro, both enzalutamide and MJC13 significantly decreased growth of MCF7-TamR cells (Fig. 6A). Enzalutamide also decreased growth of MCF7-TamR cells in soft agar, and the combination of enzalutamide + tamoxifen was more effective than enzalutamide alone (Fig. 6B).
Enzalutamide (Enza) inhibits tamoxifen-resistant tumor growth in vitro and in vivo, and AR is expressed in recurrent breast cancers. A, growth of MCF7-TamR cells treated with vehicle, tamoxifen (Tam), enzalutamide, or MJC13 for 7 days. B, MCF7-TamR cells were plated in soft agar and the number of colonies was counted after 14 days. C and D, MCF7-TamR cells were implanted into the mammary glands of nude mice with estrogen pellets and were matched into groups to receive either control chow (CTRL), tamoxifen pellets, enzalutamide-containing chow, or both (tam+enza). C, tumor growth was measured over time by luminescence. D, final tumor weights of mice from each group. E, IHC for AR and ER in clinical samples of patient-matched primary tumor and recurrence 110 months later. F, IHC for AR and ER in clinical samples of patient-matched primary tumor and metastasis 167 months later (400×). *, P < 0.05; ***, P < 0.001; ****, P < 0.0001 by ANOVA with Dunnett multiple comparison test.
We next tested whether enzalutamide could inhibit growth of tamoxifen-resistant tumor xenografts in vivo using GFP-luciferase–labeled MCF7-TamR cells. Once tumors were established, mice were matched into groups to receive CTRL chow, tamoxifen pellets, enzalutamide -containing chow, or enzalutamide + tamoxifen. Twenty days after beginning treatment, the enzalutamide -treated mice demonstrated a significant decrease in tumor viability by IVIS compared with those in the CTRL group (Fig. 6C). Each treatment resulted in a significant decrease in tumor weight compared with control-treated tumors, with enzalutamide + tamoxifen resulting in the smallest tumors by weight at the end of the experiment (Fig. 6D). TUNEL staining revealed increased apoptosis in each of the treatment groups compared with CTRL (Supplementary Fig. S5A). Interestingly, the combination resulted in a significant decrease in ER expression compared with CTRL or either drug alone (Supplementary Fig. S5B and S5C).
AR is expressed in recurrent ER+ breast cancers
To validate the potential clinical utility of anti-androgens as a therapy for advanced ER+ tumors refractory to traditional anti-estrogen–directed therapy, we examined AR expression in primary tumors compared with the same patient's local recurrence or metastatic disease. Sections of formalin-fixed, paraffin-embedded breast tumors from a cohort of 192 female patients (median age of 68 years) diagnosed with breast cancer at the Massachusetts General Hospital (Partners) between 1977 and 1993, treated with adjuvant tamoxifen and followed through 1998 were stained for AR (14). Of 49 patients with ER+/AR+ primary tumors that developed local recurrence, 96% retained AR positivity (>1% cells positive) in the recurrence. Furthermore, in more than half of these cases, the ratio of AR to ER expression (percent cells positive) was higher in the recurrence compared with the primary tumor.
Of 55 patients that developed distant metastasis, 67% retained AR positivity in the metastatic lesion. Notably, one patient with an ER+/AR− primary tumor developed an ER−/AR+ metastasis. Nearly half of these metastases showed an increased ratio of AR to ER expression compared with the primary tumor. Two examples of cases in which the recurrence or metastasis displayed increased percent cells positive for AR, but decreased percent cells positive for ER compared with the primary tumor are shown in Fig. 6E and 6F. Our findings are consistent with other studies demonstrating that AR status is highly conserved in recurrences and metastases (38), and that AR is more highly expressed in metastases than ER and PR (39). Collectively, this suggests that anti-androgens may be a useful therapeutic strategy for patients with anti-estrogen–refractory disease, as AR is frequently expressed in recurrences and metastases, often at even higher levels than in the primary tumor.
Enzalutamide inhibits primary and metastatic tumor growth in vivo
To assess the effect of enzalutamide on E2- and DHT-induced growth in vivo, we utilized AR+/ER+ PT12 cells, recently cultured from a patient-derived xenograft and expressing GFP-luciferase (16). Cells were injected orthotopically in mice implanted with either E2 or DHT pellets, and once tumors were established, mice were matched on the basis of tumor burden to receive either enzalutamide-containing or CTRL chow. Although E2 induced more rapid tumor growth, DHT also stimulated tumor growth (Fig. 7A and B). This demonstrates that DHT, in the absence of E2, promotes ER+/AR+ tumor growth in vivo, similar to our previous finding with MCF7 xenografts (14). As shown in Fig. 7A, enzalutamide significantly reduced the growth rate of E2-driven tumors when compared with E2 alone [difference between treatment groups, F(1,23) = 37.41, P < 0.0001; and group*time interaction, F(3,23) = 13.75, P < 0.0001]. Enzalutamide also significantly reduced growth rate of DHT-driven tumors compared with DHT alone [difference between treatment groups, F(1,8) = 27.80, P = 0.001; and group*time interaction, F(3,24) = 11.34, P < 0.0001]. In E2-driven tumors, bromodeoxyuridine (BrdUrd) staining demonstrated that enzalutamide significantly decreased proliferation (Fig. 7C), but had no effect on apoptosis (not shown). Conversely, in DHT-driven tumors enzalutamide significantly increased apoptosis as measured by cleaved caspase-3 staining (Fig. 7D), but had no effect on proliferation.
Enzalutamide (Enza) decreases hormone-driven growth of PT12 primary tumors and metastases. A–D, 1 × 106 GFP-luciferase expressing PT12 cells were injected orthotopically into the mammary fat pad of NOD-SCID-IL2Rgc−/− mice followed by implantation of either an E2 or DHT pellet. When tumors reached an average of 39 mm3, mice were matched into the following groups: E2 with control chow (CTRL; n = 10) or enzalutamide chow (n = 10), or DHT with control chow (n = 5) or enzalutamide chow (n = 5). Tumor viability was measured by IVIS for mice with E2 pellets (A) or DHT pellets (B). C, BrdUrd staining of tumors from mice with E2 pellets. D, Cleaved caspase staining of tumors from mice with DHT pellets. E, PT12 GFP-luciferase cells were injected intracardially in NOD-SCID-IL2Rgc−/− mice. Metastatic burden of mice treated with E2 or E2 + enzalutamide was monitored using IVIS (photons/second) over 12 weeks (total signal supine+prone). F, IVIS signal from mice in supine position at 2 weeks versus 12 weeks in E2 versus E2 + enzalutamide mice. G, IVIS image of mice in the supine position in the E2 (n = 11) or E2 + enzalutamide (n = 12) groups after 12 weeks, red denotes high IVIS signal. Error bars, SEM. *, P < 0.05; ***, P < 0.001; ****, P < 0.0001 by repeated-measures mixed model approach for (A and B), ANOVA with Bonferroni's multiple comparison test for (C and D), and Wilcoxon rank sum test for E.
To identify the molecular mechanisms by which enzalutamide decreased E2-induced tumor growth, we performed RNA-seq on PT12 tumors from E2-treated mice. Enzalutamide significantly altered 484 genes (P < 0.05, fold change > 1.2); 144 upregulated and 340 downregulated compared with E2 alone (Supplementary Table S2). Of these, 107 (22.1%) of the genes affected by enzalutamide were previously identified as regulated by estradiol in the original PT12 xenograft model (Supplementary Fig. S6B; ref. 16). Metacore analysis of the 484 genes altered by enzalutamide treatment identified AR and ER as among the transcription factors most highly implicated as upstream regulators (Supplementary Table S3A). Gene set over-representation analysis also identified AR regulation as a highly enriched pathway, as well as the HIF-1α and HIF-2α networks (Supplementary Table S3B), which have previously been associated with AR in prostate cancer (40). Finally, of the 340 genes downregulated by enzalutamide in PT12 tumors, 56 were also identified in our AR ChIP-seq experiment as being the nearest gene to sites bound by AR in response to E2 treatment in MCF7 cells. Notably, several genes decreased by enzalutamide are reported to be both ER targets and critical for ER activity including GREB1, an E2-responsive ER coactivator (41), and the histone demethylases KDM3A and KDM4B, which mediate ER binding to target gene promoters (42, 43). Together, these data confirm our in vitro observations and show that enzalutamide alters expression of ER target genes in vivo.
Finally, as we found that AR is frequently expressed in metastases of ER+ breast cancers, we tested whether enzalutamide could inhibit metastatic growth in vivo. PT12 cells were injected intracardially into mice implanted with E2 pellets, and mice were randomized onto CTRL chow or chow containing enzalutamide. Mice were monitored weekly by IVIS imaging of luciferase over 12 weeks in both the supine and prone positions. Tumors in the enzalutamide-treated group were significantly smaller at week 12 (z = −3.82, P = 0.0001, two-sided test; Fig. 7E). Next we analyzed IVIS signal of mice at week 2 (first detectable luciferase signal) versus week 12 (Fig. 7F). While control mice showed a significant increase in tumor burden over time, there was no significant increase in tumor burden in the enzalutamide-treated mice (Fig. 7F and G). This held true whether IVIS signal was measured in the supine (shown) or prone position (not shown) or both added together. These data demonstrate clearly, in a model of ER+/AR+ breast cancer recently derived from a patient, that enzalutamide is effective in reducing the growth of metastatic disease.
Discussion
AR was previously thought to antagonize ER activity because androgens such as DHT diminished the transcriptional and proliferative response of breast cancer cells to E2 (2), likely because AR and ER compete for some of the same binding sites on chromatin. However, our previous and current studies indicate that DHT is proliferative in the context of no E2, as would be the case in a postmenopausal woman treated with AIs (14). The AR antagonist bicalutamide has also been shown to increase E2-induced ER activity (2). But unlike bicalutamide, enzalutamide and MJC13 are newer generation anti-androgens that inhibit AR nuclear localization, and this is the first study to test the effects of new-generation AR inhibitors on ER chromatin binding. By inhibiting AR nuclear localization or decreasing AR expression by shRNA, we have discovered that AR supports ER genome binding and activity in breast cancer.
In response to E2, AR translocated to the nucleus in ER+ cell lines and bound chromatin at sites that overlap with ER-binding sites and are enriched for EREs. Inhibiting nuclear localization of AR with enzalutamide or MJC13 dramatically decreased E2-induced ER chromatin binding, with the greatest effects observed at sites also bound by AR. Both the AR antagonist enzalutamide and AR knockdown decreased baseline and E2-induced breast cancer cell proliferation in vitro, and enzalutamide decreased both DHT- and E2-stimulated growth of ER+/AR+ xenografts as well as metastatic burden in vivo.
These results further underscore the crosstalk between AR and ER in dual-positive breast cancer cells. There is evidence that DHT metabolites can have estrogenic effects and stimulate breast cancer growth through ER activation (44). However, enzalutamide inhibits growth differently in E2-driven tumors, where it decreases proliferation, compared with DHT-driven tumors, where it increases apoptosis. This suggests that DHT-driven growth in ER+/AR+ xenografts is not mediated through ER, but rather directly through AR. Conversely, the expression of the AR cofactor ARA70 can result in E2 having a weak agonist effect on AR (45, 46), which could explain how AR is translocated to the nucleus in response to E2 in our studies. Another possibility is that AR is activated by growth factor pathways subsequent to ER activation. However, E2 only drove nuclear localization of AR in ER+ cell lines, suggesting an ER-dependent mechanism. Further studies are ongoing to determine the mechanism of AR nuclear translocation in response to E2, and whether ARA70 is necessary for AR DNA binding in response to E2.
Although our findings may seem contradictory to prior studies, in actuality they are not mutually exclusive. Ligand-bound AR interfered with E2-mediated ER activity (2) and diminished E2-induced upregulation of a subset of ER target genes (35), likely due to competition between AR and ER for some of the same genome-binding sites; this is consistent with our observation that AR can bind to ER-binding sites in the genome. In addition, a recent study using AI-resistant MCF7 cells found that ER and AR cooperate on known androgen- and estrogen-responsive gene promoters (10).
We propose that in ER+/AR+ breast cancer cells, AR supports ER nuclear localization and genome binding, possibly by increasing chromatin availability of ER-bound loci (47), by stabilizing ER binding to chromatin, and/or interacting directly with ER or as part of an ER-containing complex, as the proximity ligation assay suggests. While this challenges the current view of AR as antagonizing ER activity, a similar effect has been observed with retinoic acid receptor-α (RARα), which interacts with ER-binding sites in an ER-dependent manner in ER+ breast cancer cells (48). This interaction is required for E2-induced proliferation and ER transcriptional activity (48). In addition, glucocorticoid receptor (GR), which is highly similar to AR (49, 50), increases chromatin availability and subsequent ER binding at response elements bound by both receptors, a mechanism termed “assisted loading” (47). Similarly, our data show that ER chromatin binding is most inhibited by anti-androgens at sites where E2 induces binding of AR and ER. Thus, even though androgen-bound AR can diminish E2-stimulated ER activity, we show that anti-androgens that prevent AR nuclear translocation have the same effect, suppression of ER activity, via an entirely different molecular mechanism.
Importantly, the assays herein were performed using endogenous AR and ER in cells that naturally express both receptors. In light of recent data from our laboratory and others suggesting that the ratio of AR:ER protein expression is a predictor of response to traditional ER-directed endocrine therapy (14) and DCIS progression (51), it is likely that the interplay of these receptors may depend on their relative expression level, the levels of their respective ligands in circulation and within tumors, and levels of shared cofactors such as FOXA1. Our data show that across multiple cell lines and preclinical models of ER+/AR+ breast cancer, AR antagonists such as enzalutamide and MJC13 that inhibit AR nuclear translocation also inhibit ER activity indirectly. This combined effect on AR and ER may account for the synergy demonstrated between enzalutamide and the anti-estrogens tamoxifen and fulvestrant in vitro. Further analysis of these data is ongoing to determine the mechanisms of this synergy.
Collectively, these data strongly suggest that in the most common form of breast cancer, ER+/AR+ disease, primary tumors and recurrent disease may become reliant on AR and that AR may serve as an effective therapeutic target either in combination with traditional ER-directed therapies (particularly in tumors that have a high AR:ER protein ratio) or upon resistance to ER-directed therapies. Our current and prior (14) studies on the role of AR in ER+ breast cancer contribute to a deeper understanding of the complex molecular interplay between the two most widely expressed hormone receptors in breast cancer (AR and ER), and have already led to clinical trials testing the efficacy of enzalutamide in combination with the AI exemestane in patients with advanced ER+ disease (NCT02007512) and in combination with fulvestrant (NCT01597193).
Finally, our data demonstrate that enzalutamide effectively inhibits growth of ER+/AR+ metastases in vivo. A recent study of ER− PDX models that metastasize from the orthotopic site found that AR mRNA was increased in circulating tumor cells and micrometastases compared with the primary tumors (52), indicating AR may be an important target for inhibition of metastasis. Likewise, we show in clinical specimens of patient-matched ER+/AR+ primary tumors compared with local or distant recurrences occurring during tamoxifen treatment, AR expression is often maintained, and sometimes increased, in breast cancers refractory to anti-estrogen therapy. We are actively investigating the specific role that AR plays in facilitating the process of metastasis in both ER+ and ER− breast cancer.
Disclosure of Potential Conflicts of Interest
V.T. Phan is a senior scientist and has ownership interest (including patents) in Medivation, Inc. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: N.C. D'Amato, A. Elias, B.M. Jacobsen, J.K. Richer
Development of methodology: N.C. D'Amato, N.S. Spoelstra, V.T. Phan, B.M. Jacobsen, J.K. Richer
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): N.C. D'Amato, B. Babbs, N.S. Spoelstra, K.T. Carson Butterfield, T.J. Rogers, C.A. Sartorius, B.M. Jacobsen, J.K. Richer
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): N.C. D'Amato, M.A. Gordon, N.S. Spoelstra, K.C. Torkko, V.T. Phan, T.J. Rogers, A. Elias, J. Gertz, B.M. Jacobsen
Writing, review, and/or revision of the manuscript: N.C. D'Amato, M.A. Gordon, N.S. Spoelstra, K.C. Torkko, V.T. Phan, T.J. Rogers, A. Elias, J. Gertz, B.M. Jacobsen, J.K. Richer
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): N.C. D'Amato, M.A. Gordon, K.T. Carson Butterfield
Study supervision: B.M. Jacobsen, J.K. Richer
Grant Support
This work was supported by grants BC120183 W81XWH-13-1-0090 DOD BCRP Clinical Translational Award (to J.K. Richer and A. Elias); R01 CA187733-01A1 (to J.K. Richer); American Cancer Society Award124475-PF-13-314-01-CDD (to N.C. D'Amato); and NIH/NCI Cancer Center Support Grant P30CA046934.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Acknowledgments
The authors thank A. Protter at Medivation, Inc. for comments on the manuscript. We thank Dr. Marc Cox for graciously providing MJC13 for these studies. We thank Active Motif, Inc. for their assistance with AR and ER ChIPseq experiments. The authors also acknowledge the Genomics and Microarray Core and other shared resources of Colorado's NIH/NCI Cancer Center Support Grant P30CA046934.
Footnotes
Note: Supplementary data for this article are available at Molecular Cancer Research Online (http://mcr.aacrjournals.org/).
B.M. Jacobsen and J.K. Richer are co-senior authors of this article.
- Received May 11, 2016.
- Revision received July 29, 2016.
- Accepted August 2, 2016.
- ©2016 American Association for Cancer Research.