Proteomic profiling identifies key coactivators utilized by mutant ERα proteins as potential new therapeutic targets


Identification of co-regulators recruited to ERα mutants

Since Y537S, D538G, and ESR1-YAP1 ERα (Fig. 1a) promote E2-independent transcriptional activation of an ERE-dependent reporter (Fig. 1b), we tested whether this activation is through recruiting or repelling distinct co-regulators (coactivators and co-repressors). Using our 4xERE DNA pulldown assay to identify co-regulators recruited to ERα [21], we first utilized recombinant WT and Y537S ERα, along with HeLa S3 nuclear extract (HNE) to form complexes. Washed complexes were then subjected to label-free, quantitative MS and bound proteins were normalized by the amount of ESR1 N-terminal peptides bound (Fig. 1c, top; Supplementary Table 2).

Fig. 1
Fig. 1

Proteomics of co-regulators recruited to the Y537S ERα LBD point mutant. a Schematic diagram of WT, Y537S, and D538G ERα point mutants (mutations indicated by arrows), and the ESR1-YAP1 fusion protein. Numbers refer to amino-acid residues in ERα and YAP1. AF1 activation function 1, DBD DNA-binding domain. Hinge, region between DBD and ligand-binding domain (LBD); AF2 Activation function 2, WW WW domain, TAD transcription activation domain. In ESR1-YAP1, blue represents ERα residues (1–365); yellow indicates YAP1 residues (230–504). b Mutant ERα proteins display E2-independent transcriptional activity. Vectors expressing YFP or YFP-tagged WT, Y537S, D538G, or ESR1-YAP1 proteins were co-expressed with an ERE-dependent luciferase reporter (pERE-E1b-luc) in HeLa cells grown in charcoal-stripped fetal bovine serum (FBS). Cells were then treated with/without 10 nM E2 for overnight. Luciferase activity (RLU relative light units) was assayed from whole-cell extracts of the cells transfected in triplicate. Data are represented as mean ± SEM (n = 3); ***p < 0.001. c MS data depicted as a heatmap for WT or Y537S ERα-dependent coactivators recruited to EREs. (Top) Schematic of ERE DNA pulldown assay using HeLa S3 NE as the source of co-regulators (modified from ref. [21]) and purified ERα proteins. (Bottom) MS data were analyzed from duplicate reactions using a label-free method and depicted as in ref. [21]. Peptides number of peptides detected; amount (vsESR1) amount normalized by sum of area under the curve for six N-terminal ESR1 peptides (see Supplementary Table 2). Fold change represents the ratio of amount detected normalized to unliganded WT ERα. Fold change cutoff used was ≥1.5. SRC-1 to -3 (gene symbols: NCOA1, NCOA2, NCOA3), p300 (gene symbol: EP300), CBP (gene symbol: CREBBP). d Immunoblotting validation of KMT2D and SRC-3 enrichment with purified Y537S ERα bound to EREs using independent DNA pulldown samples in the absence of E2. ERα protein binding was assayed by an N-terminal (N-term) antibody. TBP serves as a loading control. 2% input represents 2% of the starting HNE employed in the ERE DNA pulldown. e Y537S ERα and the KMT2D complex interact directly in an enhanced manner relative to unliganded WT ERα. ERE DNA pulldown assays were performed with purified ERα proteins and a purified KMT2D “fusion” complex [22], and then analyzed by immunoblotting. (i) Representative immunoblot probing for select KMT2D complex members and ERα. (ii) Quantification of KMT2D signal across three independent pulldown assays. Binding was quantified using Image J and normalized to ERα signal

Compared with unliganded WT ERα, a subset of coactivators were recruited in an enhanced manner to Y537S ERα (Fig. 1c, bottom). Namely, the histone H3 lysine 4 (H3K4) methyltransferase KMT2D complex displayed the greatest enrichment with Y537S ERα, along with SRC-1, -2, and -3, p300, CBP, and KMT2D’s paralog, KMT2C. Immunoblotting validated the enhanced KMT2D and SRC-3 recruitment (Fig. 1d).

As SRC-3 directly binds Y537S ERα [14, 19], we tested whether a KMT2D complex [22] would directly interact with Y537S ERα. Indeed, the binding of the KMT2D complex was enhanced with Y537S, along with E2-bound WT, compared to unliganded WT ERα (Fig. 1e). We found other potential coactivators with enhanced binding to Y537S ERα (Supplementary Figure 1a). We further validated PELP1 recruitment by immunoblotting, as an inhibitor disrupting this interaction is described [23] (Supplementary Figure 1b). We found very few co-repressors had reduced recruitment to Y537S ERα (Supplementary Figure 1c).

As purified D538G ERα failed to recruit SRC-3 despite binding EREs (data not shown), we resorted to using extracts from transfected 293T cells as sources of WT, Y537S, and D538G ERα proteins. We again observed enhanced recruitment of KMT2D and SRCs to Y537S compared to unliganded WT ERα (Supplementary Figure 2). The D538G mutant also recruited these coactivators, but three to four times less than that of Y537S. Additional potential co-regulators displayed enhanced binding to both ERα mutants. However, for subsequent functional characterization, we chose to focus on SRCs and the KMT2D complex.

SRCs are critical for ERα LBD mutant activity and cell growth

We next tested the functional consequence of enhanced SRC-ERα mutant interactions on transcription. Knockdown of SRC-3 using published siRNAs [24] significantly reduced both Y537S and D538G ERα-mediated transcriptional activity on the ERE-Luc reporter (Fig. 2a). Treatment with a small molecule inhibitor (SMI), SI-1, which inhibits the activities of all three SRCs (IC50 = 0.2 μM) [25] and reduces SRC protein levels (Fig. 2b), severely reduced the transcriptional activities of WT and mutant ERα (Fig. 2c). Cell viability was minimally affected (Supplementary Figure 3a). We also tested whether the combination of an oral SERD and SI-1 would further reduce LBD mutant ERα transcriptional activity. We focused on AZD9496 [26] (AZD) as it: (1) reduced endogenous ERα protein, (2) was significantly more potent than ICI182,780 (ICI, fulvestrant) in reducing mutant ERα transcriptional activities (unlike another SERD GDC-0810 [27] (GDC; Supplementary Figure 3b-d), and (3) AZD was reported as more effective than ICI at inhibiting tumor growth promoted by Y537S ERα [14]. The combination of SI-1 and AZD synergistically reduced both Y537S and D538G activities on the ERE-luciferase (Luc) reporter (Fig. 2d, e, Supplementary Table 4).

Fig. 2
Fig. 2

SRCs are key co-regulators of WT, Y537S, and D538G ERα transcriptional activity and are essential for breast cancer cell growth. a Knockdown of SRC-3 by co-transfection of HeLa cells with an siRNA targeting pool reduced transcriptional activity of LBD point mutant ERα proteins, as assayed by an ERE-luciferase reporter, in hormone-depleted media. Data are represented as mean ± SEM (n = 3); ***p < 0.001. RLU relative light units. b A “pan-SRC” inhibitor, SI-1, reduced all three SRC protein levels in MCF-7 breast cancer cells treated overnight, as assayed by immunoblotting. β-actin serves as a loading control. c SRC inhibitor SI-1 reduced WT, Y537S, and D538G ERα transcriptional activity at concentrations higher than the established IC50. Vectors expressing YFP-tagged WT, Y537S, and D538G ERα proteins were co-expressed with pERE-E1b-luc in HeLa cells, and then cells were treated with dimethyl sulphoxide (DMSO, vehicle control) or SI-1 overnight. Luciferase activity was measured as in Fig. 1b. Data are represented as mean ± SEM (n = 3). d,e Combination of an SRC inhibitor (SI-1) and an oral SERD (AZD9496) synergistically reduce Y537S (d) and D538G (e) mutant ERα transcriptional activity. Experiment was done as in (c), except that AZD9496 was added with or without SI-1 to co-transfected HeLa cells. Data are represented as mean (n = 3). f The oral SERD AZD9496 is more effective than the SRC inhibitor SI-1 in reducing cell viability of MCF-7 lines expressing WT or Y537S ERα. The lentiviral transduced MCF-7 stably expressing cell lines were treated with vehicle (DMSO) or different concentrations of AZD9496 or SI-1 as indicated. After 6 days of treatment, viability was assayed by a MTT assay. Data are represented as mean ± SEM (n = 3). Synergism was observed with combination (Combo) treatments 4 and 5 (12.5/200 nM AZD9496/SI-1; 25/400 nM AZD9496/SI-1) in Y537S ERα-expressing cells (shown as red arrows; CI values are shown in Supplementary Table 5)

We tested whether SRC inhibition would affect growth of stably expressing WT or Y537S ERα MCF-7 cell lines [28]. SI-1 at 400 nM reduced viability in both WT and Y537S ERα-expressing cells by 91% (Fig. 2f). When combinations of SI-1 with AZD or ICI were tested, a synergistic reduction in viability of Y537S ERα-expressing cells was found at the two highest combined doses or highest combined dose, respectively (Fig. 2f and Supplementary Figure 3e; Supplementary Tables 5 and 6). Thus, SI-1 combined with AZD was most effective in reducing both transcription and cell growth mediated by mutant ERα proteins.

Inhibiting SRCs and mutant ERα most effectively reduces patient-derived xenograft tumor growth

We next tested the efficacy of an improved pan-SRC inhibitor (SI-2) [25], AZD, or the combination in a patient-derived xenograft (PDX) expressing Y537S ERα (WHIM 20 [11]) for tumor reduction (Fig. 3a). SI-2 was chosen instead of SI-1, given its reduced IC50 (3.4 nM) and ability to reduce ER tumor growth [25]. After tumors grew to 350 mm3, mice were randomized, E2 was withdrawn to mimic AI treatment, and mice were then treated with control vehicle, SI-2, AZD, or the combination. After 4.5 weeks, SI-2 alone significantly reduced tumor volume, AZD gave a larger reduction, but the combination gave the most significant reduction in tumor growth.

Fig. 3
Fig. 3

Inhibition of SRCs and ERα reduces tumor burden in a PDX model of Y537S ERα breast cancer. Significance was determined using one-way ANOVA and Tukey’s test to correct for multiple comparisons. *p < 0.05; **p < 0.01; ***p < 0.001. a Tumor volume (n = 8) is reduced with treatment of either SI-2 or AZD9496, and further reduced upon treatment with a combination of both inhibitors. Data are represented as mean ± SEM. b–e Quantification of tumor immunoblotting (n = 6/treatment) performed with Image J analysis relative to GAPDH expression or imaging data (n = 8) unless otherwise indicated. Data are represented as mean ± SD. ud undetected. Color legend is the same as in (a). b Quantification of ERα, SRC-1, and SRC-3 protein expression across tumors. By Dixon’s Q test, the tumor samples in the last lane probed for SRC-1 and SRC-3 (Supplementary Figure 4a) were outliers at 99% confidence and thus excluded from the Image J analysis. c The ERα target gene PR displays reduced expression with AZD9496 and combination treatment. d BrdU is decreased in the AZD9496 and combination therapy groups (n = 8 mice per group, with mean ± SD of three slides for each mouse counted). e Cleaved PARP is increased with SI-1 treatment and decreased with AZD9496 and combination therapies

We next confirmed that the drugs indeed affected their intended targets and tested for effects on apoptosis and proliferation. First, as expected, AZD treatment significantly reduced ERα expression in tumor lysates. Unexpectedly, AZD treatment upregulated SRC expression (Fig. 3b, Supplementary Figure 4a), which may have relevance for patients receiving AZD monotherapy in clinical trials (NCT02248090/NCT03236974). Second, we tested the expression of a classical ER target gene, PR (Fig. 3c, Supplementary Figure 4a). While SI-2 did not affect PR expression, AZD clearly did. Third, we found that SI-2 increased an apoptosis marker (cleaved PARP protein), while AZD decreased proliferation as measured by BrdU incorporation (Fig. 3d, e, Supplementary Figure 4b). Finally, we did not observe any significant toxicity with any drug treatment after examining recipient mouse livers by histochemistry and measuring body weights (Supplementary Figures 4c-d). Thus, our PDX data support a potential new treatment regime for breast cancers bearing ESR1 LBD mutations, which is to combine a SRC inhibitor with an oral SERD.

KMT2C/2D are novel coactivators for Y537S ERα

From above, we found that the KMT2C/2D complexes were preferentially enriched with Y537S ERα (Fig. 1c, Supplementary Figure 2a). To determine the functional role of KMT2C/2D, we tested whether their depletion would affect Y537S ERα-mediated reporter expression. Knocking down KMT2C, KMT2D, or both together reduced Y537S ERα transcriptional activity (Fig. 4a). Upon double KMT2C/2D knockdown, WT ERα transcriptional activity was also reduced (Fig. 4b). However, ESR1-YAP1 transcriptional activity was not affected, ruling out a general transcriptional effect.

Fig. 4
Fig. 4

Knockdown of KMT2C/2D reduces WT and Y537S ERα transcriptional activity and breast cancer cell growth. a HeLa cells grown in phenol red-free, charcoal-stripped media were co-transfected with pERE-E1b-luc, YFP-tagged Y537S ERα, and different siRNAs (25 nM each for 50 nM total). Cell lysates were assayed for luciferase activity (RLU). Data are represented as mean ± SEM (n = 3); ***p < 0.001. NC#1 siRNA served as the negative control for KMT2D targeting siRNA, while Control-A siRNA pool served as the negative control for KMT2C targeting siRNA pool. b Transfection of siRNAs targeting KMT2C and KMT2D reduced expression of the pERE-E1b-luc reporter, as compared to non-targeting siRNAs (siControl), in HeLa cells co-transfected with YFP-tagged WT or Y537S ERα vectors, but not an YFP-tagged ESR1-YAP1 fusion. Luciferase activity was measured as in Fig. 1b. Data are represented as mean ± SEM (n = 3); **p < 0.01; ***p < 0.001. c Knockdown of KMT2C and KMT2D in lentiviral transduced stably expressing WT or Y537S MCF-7 cells results in reduced anchorage-independent growth in soft agar and confers sensitivity to anti-estrogens. siRNAs (same as above) were transfected into the two cell lines at a final concentration of 100 nM (50 nM each), and then re-plated in 24 well plates 24 h later. After 1 week with either vehicle or anti-estrogen (100 or 1000 nM 4-hydroxytamoxifen (Tam) or ICI) treatment, colonies formed in soft agar were counted and quantified. Data are represented as mean ± SEM (n = 4); *p < 0.05; **p < 0.01; ***p < 0.001

We additionally found that KMT2C/2D knockdown reduced anchorage-independent growth of WT, Y537S, and D538G ERα-expressing cells (Supplementary Figure 5), with significant differences observed between WT and mutant ERα data. More importantly, knockdown of KMT2C/2D significantly sensitized the partially resistant Y537S ERα cells to anti-estrogens currently given in the clinic (Fig. 4c). Overall, our data reveal an important functional role for KMT2C/2D with Y537S ERα in transcription and cell growth.

SRCs and KMT2D are crucial for growth of inducible LBD mutant ERα cells

As our above cell lines stably overexpressed LBD ERα mutants [28], we created MCF-7 cell lines supporting conditional (doxycycline (Dox)-inducible) expression of FLAG-tagged WT, Y537S, or D538G ERα. The FLAG tag did not impair the transcriptional activities of ERα proteins (Supplementary Figure 6). Dox addition to cells grown in charcoal-stripped media induced expression of these ERα proteins, with Y537S and D538G ERα accumulating to a similar extent but less than WT (Fig. 5a). We next performed co-immunoprecipitations to test whether SRC-3 and KMT2D displayed enhanced association with inducible mutant ERα proteins under hormone-depleted conditions. We found greatest KMT2D and SRC-3 association with the Y537S mutant, with less recruited to D538G ERα (Fig. 5b). Importantly, the Dox-inducible mutant ERα proteins activated endogenous ERα target genes (GREB1 and TFF1) in a hormone-independent manner and D538G ERα had a weaker effect than Y537S, in accordance with other models [7, 9, 29,30,31] (Fig. 5c and Supplementary Figure 7a). We next asked whether ablation of these key co-regulators would inhibit the viability of these cells. Consistently, SI-1 significantly reduced the viability of all ERα-expressing cells (Fig. 5d), while knockdown of KMT2C/2D selectively affected WT and Y537S ERα-expressing cells (Fig. 5e).

Fig. 5
Fig. 5

SRCs and KMT2C/2D promote cell growth and transcription of ERE-containing target genes in conditionally expressed (Dox-inducible) WT and mutant ERα MCF-7 cell lines. Prior to all assays, cell lines were grown in charcoal-stripped, phenol red-free media with/without Dox for at least 2 days. a Confirmation of FLAG-tagged WT, Y537S, and D538G ERα Dox-inducible expression in MCF-7 cells transduced with specific lentiviruses. Immunoblotting was performed on whole-cell extracts with FLAG, ERα (HC-20 antibody), or β-actin (loading control) antibodies. *, Nonspecific band. →, Position of endogenous ERα. b KMT2D and SRC-3 co-immunoprecipitate with mutant ERα. Whole-cell extracts from Dox-inducible MCF-7 cells were subject to immunoprecipitation of FLAG, followed by immunoblotting. For analysis, intensities were normalized for different precipitated ERα levels by Image J. c Dox-inducible expression of Y537S and D538G ERα in MCF-7 cells results in activation of two canonical ERα target genes independent of E2. Cell lines were treated with Dox and the WT ERα line was treated (±) 10 nM E2 overnight, followed by RNA isolation. Relative levels of GREB1 or TFF1 mRNAs were determined by RT-qPCR with ESR1 mRNA as the normalizer. Data are represented as mean ± SEM (n = 3); **p < 0.01; ***p < 0.001. d, e Cell viability (measured by Cell Titer Glo as RLU) of Dox-inducible WT, Y537S, and D538G ERα expressing MCF-7 lines (induced as above). *p < 0.05; **p < 0.01; ***p < 0.001. d Cell viability is reduced in all three cell lines upon treatment with SRC inhibitor SI-1 (1 μM) after 3 days’ exposure. DMSO served as the vehicle control. Data are represented as mean ± SEM (n = 3). e Cell viability of Dox-inducible WT and Y537S ERα-expressing MCF-7 lines, but not the D538G-expressing ERα line, is reduced upon knockdown of KMT2C and KMT2D after 3 days exposure to 100 nM total siRNA (50 nM each). Data are represented as mean ± SEM (n = 3). f Knockdown of KMT2C and KMT2D reduces Y537S ERα-enhanced expression of GREB1 and TFF1. Cells were transfected with siRNAs (100 nM total, 50 nM each). Relative levels of GREB1 or TFF1 mRNAs were determined by RT-qPCR using ACTB mRNA as the normalizer. Data are represented as mean ± SEM (n = 3); *p < 0.05; **p < 0.01; ***, p < 0.001. g Dox-induced FLAG-tagged WT, Y537S, and D538G ERα proteins occupy EREs of GREB1 and TFF1 genes in MCF-7 cells, implying direct transcriptional regulation. In contrast, ERα proteins minimally occupy a negative control region from intron 4 of the CCND1 gene, which lacks endogenous ERα binding [36]. Where indicated, WT ERα cells were treated with 100 nM E2 for 45 min. ChIP assays employed an antibody against FLAG to IP the FLAG-tagged ERα proteins and associated DNA. Isolated DNA was assayed by ChIP-qPCR. Representative data were plotted relative to percentage of starting input chromatin and are represented as mean of triplicate qPCR reactions ± SEM. Supplementary Figure 7d shows a repeated ChIP assay. h KMT2D occupies EREs of GREB1 and TFF1 genes in a Dox-dependent manner correlating with increased Y537S ERα occupancy. ChIP-qPCR was performed using an antibody to KMT2D. Representative data were plotted as above and the CCND1 gene intron 4 served as a negative control region. Supplementary Figure 7e shows a repeated ChIP assay

Knockdown of KMT2C/2D modulates Y537S ERα direct target gene expression

We next tested the effect of KMT2C/2D depletion on GREB1 and TFF1 gene expression in our Dox-inducible MCF-7 cells. KMT2C/2D knockdown reduced Y537S ERα-mediated transcription of both genes (Fig. 5f and Supplementary Figure 7b). In WT ERα cells, the loss of KMT2C/2D also reduced TFF1 mRNA, but GREB1 or ESR1 were unaffected. We extended our analysis to 10 total ERα target genes by depleting only KMT2D using a validated siRNA [32], as KMT2D had greater recruitment than KMT2C to Y537S ERα (Fig. 1c). We observed that Y537S mutant expression regulates select ERα targets, a subset of which is reduced by KMT2D depletion (Supplementary Figure 7c).

Chromatin immunoprecipitation (ChIP) was used to examine whether the ERα mutants bound EREs [33,34,35] located upstream of the GREB1 and TFF1 gene transcription start sites (TSS) in a Dox-dependent manner (Fig. 5g and Supplementary Figure 7d). All ERα proteins displayed Dox-dependent enrichment on EREs upstream of the GREB1 and TFF1 TSS, but not on a negative control region [36]. Thus, the binding of the LBD mutant ERα proteins to multiple EREs suggests direct regulatory roles in regulation of these genes.

Chromatin occupancy of Y537S ERα and KMT2D are positively correlated

As SRC-3 and p300 are co-localized with Y537S ERα on chromatin [29], we next tested whether KMT2D chromatin occupancy correlated with Y537S ERα. We performed ChIP using a validated KMT2D antibody [37, 38] (Fig. 5h and Supplementary Figure 7e). KMT2D occupancy increased in Y537S ERα cells in a Dox-dependent manner for all EREs assayed, but not for the control region, without increased KMT2D protein expression (Supplementary Figure 7f). Thus, Dox-induced occupancy of EREs by Y537S ERα is correlated with increased recruitment of KMT2D and enhanced transcription of GREB1 and TFF1.

Co-regulators recruited to the ESR1-YAP1 fusion

We next investigated co-regulator recruitment to ESR1-YAP1, which contains the N terminus and DBD, but lacks the LBD, of ERα fused in-frame to the C terminus of the Yes-associated protein 111 (Fig. 1a). ESR1-YAP1 promotes high levels of expression of an ERE-dependent reporter without E2 (Fig. 1b and Supplementary Figure 6a). We compared co-regulator recruitment to either purified WT ERα or ESR1-YAP1 bound to ERE DNA in the absence of E2 (Fig. 6a). After normalization for ERα-binding differences, the ESR1-YAP1 protein did not recruit more SRC-3 and KMT2D as compared to WT ERα, unlike the LBD ERα mutants (validated in Fig. 1d). Instead, the ESR1-YAP1 protein recruited many subunits of the 26S proteasome (Fig. 6a). Immunoblotting validated enhanced proteasome recruitment (SUG1/PSMC5 and 20S C2/PSMA1) of ESR1-YAP1 vs. WT ERα, even with E2 (Fig. 6b). We further validated that the 26S proteasome was recruited to ESR1-YAP1 in E2-deprived T47D stable cell line [11] extracts (Supplementary Figure 8a).

Fig. 6
Fig. 6

The ESR1-YAP1 fusion protein recruits the 26S proteasome for activated transcription of an ERE-driven luciferase reporter. a MS data depicted as a heatmap (displayed as in Fig. 1c) for ESR1-YAP1-dependent coactivators recruited to EREs from HNE. Recombinant purified WT or ESR1-YAP1 proteins were added to duplicate ERE DNA pulldown reactions. Fold change cutoff was ≥1.5 for enrichment over unliganded WT ERα. *ESR1 normalized ESR1 (see Supplementary Table 2). **YAP1 YAP1 corrected for ESR1. b Immunoblotting validation of the 26S proteasome being recruited to ESR1-YAP1. Independent ERE DNA pulldown samples were used to detect proteins recruited to ESR1-YAP1 compared to WT with/without E2 (i) or compared to WT or Y537S ERα (ii). Levels of ERα bound to the EREs were determined with an ERα antibody recognizing an N-terminal (N-term) epitope. TBP served as a loading control. 3% HNE, 3% of the starting HeLa S3 NE employed in the ERE DNA pulldown. c A 26S proteasome inhibitor, MG132, reduces ESR1-YAP1 transcriptional activity on an ERE-driven luciferase reporter. HeLa cells grown in phenol red-free, charcoal-stripped media were co-transfected with a vector expressing YFP-tagged ESR1-YAP1 protein and pERE-E1b-luc. Cells were then treated with vehicle (0.1% DMSO) or 1 μM MG132 for overnight. Luciferase activity (RLU relative light units) was assayed from whole-cell lysates. Data are represented as mean ± SEM (n = 3); ***p < 0.001

The 26S proteasome plays a role in WT ERα degradation that is linked with the receptor’s ability to activate ERE-driven reporter genes [39,40,41]. We thus tested whether a proteasome inhibitor, MG132, would similarly inhibit ESR1-YAP1 transcriptional activity. Indeed, MG132 treatment of cells transfected with an ERE-driven reporter and an ESR1-YAP1 expression plasmid reduced Luc activity, as compared to the vehicle control (Fig. 6c). Interestingly, MG132 or the FDA-approved proteasome inhibitor, bortezomib, could inhibit the transcriptional activity of a GAL4 DBD-YAP1 fusion (Supplementary Figure 8b).

Proteasome activity is important for ESR1-YAP1-mediated cell growth and gene expression

We next tested the functional significance of the ESR1-YAP1: 26S proteasome interaction in T47D stable lines expressing HA-tagged YFP, WT ERα, or ESR1-YAP1 proteins (Fig. 7aii and Supplementary Figure 8c). Increasing concentrations of bortezomib treatment for 3 days efficiently reduced the growth of all three T47D lines compared to the vehicle control (Fig. 7a,i).

Fig. 7
Fig. 7

Proteasome inhibition reduces viability of breast cancer cells expressing the ESR1-YAP1 fusion and modulates ESR1-YAP1 target gene expression. T47D cell lines expressing these constructs were grown in phenol red-free, charcoal-stripped media for at least 2 days. a Increasing concentrations of bortezomib reduce cell viability of HA epitope-tagged YFP, WT ERα, and ESR1-YAP1 stably expressing T47D cell lines. (i) Cell viability was measured (via Cell Titer Glo as RLU) after 3 days of bortezomib or vehicle (DMSO) treatments. Data are represented as mean ± SEM (n = 3). (ii) Expression levels of ERα proteins were assayed in whole-cell extracts made from HA-tagged YFP, WT ERα, and ESR1-YAP1 expressing T47D cell lines. The N-terminal ERα antibody was used to detect endogenous ERα (denoted by →) as well as the ESR1-YAP1 fusion. β-actin served as a loading control. b Expression of HA-tagged ESR1-YAP1 activates expression of two classical ERE-containing target genes as compared to E2-deprived HA-tagged WT ERα. WT ERα-expressing cells were treated with/without 10 nM E2 for 24 h. Relative levels of TFF1 or PGR mRNAs were determined by RT-qPCR using ACTB mRNA as the normalizer. Data are represented as mean ± SEM (n = 3); **p < 0.01; ***p < 0.001. c ESR1-YAP1 directly occupies certain EREs of the TFF1 and PGR genes. Where indicated, HA-tagged WT ERα cells were treated with 100 nM E2 for 45 min before ChIP assays. ChIP-qPCR assays employed an antibody against HA to IP the HA-tagged ERα proteins and associated DNA. Representative data were plotted relative to percentage of starting input chromatin, which was represented as the mean of triplicate qPCR reactions ± SEM. CCND1 gene intron 4 served as a negative control region. Supplementary Figure 8d shows a repeated ChIP assay. d Proteasome inhibitor treatment of T47D cells expressing HA-tagged ESR1-YAP1 modulates ERE-containing target gene expression. Cells were treated with vehicle (0.1% DMSO) or 4 or 16 nM bortezomib for 17 h. Relative levels of PGR, TFF1, ESR1, or ESR1-YAP1 mRNAs were determined by RT-qPCR using GAPDH mRNA as the normalizer. Data are represented as mean ± SEM (n = 3); *p < 0.05; ***p < 0.001

We next wanted to define the effect of ESR1-YAP1 expression in T47D cells grown in an E2-deprived state on ERα target gene expression. The ESR1-YAP1 fusion promoted the expression of two target genes, TFF1 and PGR, significantly above the level of unliganded WT ERα-expressing cells (Fig. 7b), even though much less ESR1-YAP1 fusion was expressed (Fig. 7a, ii and Supplementary Figure 8c). The regulation of these two genes is likely direct, as ChIP assays revealed ESR1-YAP1 occupancy at two defined ERE enhancer-like sequences [35, 42] (Fig. 7c and Supplementary Figure 8d).

The effect of proteasome inhibition on TFF1 and PGR gene expression in E2-deprived ESR1-YAP1-expressing cells was tested by treating cells with bortezomib (Fig. 7d). Proteasome inhibition had both a dose-dependent and gene-specific effect on the ESR1-YAP1 targets, decreasing PGR and increasing TFF1, which resembles the effect of proteasome inhibitors on these genes after E2 treatment of MCF-7 cells [43,44,45,46]. Bortezomib decreased endogenous ESR1 mRNA expression, consistent with prior reports [43, 47]. Finally, we observed that bortezomib stimulated ESR1-YAP1 mRNA expression (driven by the CMV promoter), which may explain why bortezomib did not severely reduce ESR1-YAP1 cell viability as compared to overexpressed WT ERα (Fig. 7a). Thus, the proteasome modulates ESR1-YAP1 target genes and cell growth, suggesting a new approach for treating tumors bearing this class of resistance mutation.

Potential clinical relevance of mutant ERα-binding coactivators

To investigate whether the expression levels of the coactivators identified in this study correlate with patient outcomes, we queried two existing expression data sets—the Symmans Breast 2 ER-positive tamoxifen-treated patients [48] or ER-positive breast cancer patients treated with endocrine therapy (in KM plotter [49]; Supplementary Figure 9). In the Symmans data set, we found significantly higher KMT2D mRNA levels in patients who had a metastatic occurrence after 3 years of tamoxifen treatment vs. those that had not (panel a). We also found that higher SRC-3 and proteasomal subunit mRNA levels correlated with reduced survival from distant metastasis (panels b, c). In the KM plotter analysis, we found that higher KMT2D, but not KMT2C, mRNA significantly correlated with reduced recurrence-free survival (panel d). Finally, we assayed a metastatic breast cancer mutation/amplification database [50] and found that the ESR1 LBD mutations were not present in patients with mutations in either KMT2C or KMT2D genes (panel e).

Source link


Please enter your comment!
Please enter your name here