The human oncoprotein, mucin 1 (MUC1), drives tumorigenesis in breast carcinomas by promoting epithelial-to-mesenchymal transition (EMT), epigenetic reprogramming, and evasion of immune response. MUC1 interacts with STAT1, through JAK/STAT signaling, and stimulates transcription of IFN-stimulated genes, specifically IFN-induced transmembrane protein 1 (IFITM1). Our laboratory has previously shown that IFITM1 overexpression in aromatase inhibitor (AI)-resistant breast cancer cells promotes aggressiveness. Here, we demonstrate that differential regulation of MUC1 in AI-sensitive (MCF-7 and T-47D) compared with AI-resistant (MCF-7:5C) cells is critical in mediating IFITM1 expression. A tumor microarray of 94 estrogen receptor–positive human breast tumors correlated coexpression of MUC1 and IFITM1 with poor recurrence-free survival, poor overall survival, and AI-resistance. In this study, we investigated the effects of MUC1/IFITM1 on cell survival and proliferation. We knocked down MUC1 levels with siRNA and pharmacologic inhibitors, which abrogated IFITM1 mRNA and protein expression and induced cell death in AI-resistant cells. In vivo, estrogen and ruxolitinib significantly reduced tumor size and decreased expression of MUC1, P-STAT1, and IFITM1.

Implications:

MUC1 and IFITM1 overexpression drives AI resistance and can be targeted with currently available therapies.

Visual Overview:http://mcr.aacrjournals.org/content/molcanres/17/5/1180/F1.large.jpg.

An estimated 70% of breast cancers express the estrogen receptor alpha (ERα), the progesterone receptor (PgR), or both, and estrogen stimulation of these receptors is a significant factor in the development and growth of breast cancer (1). In the context of ER-positive breast cancer, disrupting estrogen activity using aromatase inhibitors (AI) represents the primary and most effective approach to therapy. AIs inhibit the growth of ER-positive breast cancer by blocking the aromatase enzyme, which regulates the conversion of androgen to estrogen in peripheral tissues, thus resulting in an almost complete loss of circulating estrogen (2). Despite the effectiveness of these agents, however, 30% to 40% of patients eventually develop resistance to AI treatments. Clinical recurrence of breast cancer after AI treatment usually occurs within 3 years of the conclusion of therapy (3, 4). A better understanding of the mechanisms of AI resistance may contribute to the development of new therapeutic strategies and aid in the search for new therapeutic targets and agents.

Mucin 1 (MUC1) is a transmembrane heterodimeric glycoprotein that is localized to the apical border of epithelial cells and is aberrantly overexpressed in approximately 90% of breast cancers (5–8). The MUC1 N-terminal (MUC1-N) subunit is the heavily glycosylated mucin component of the heterodimer and it supports the mucous layer, functions in cell adhesion and immune suppression, and aids in the trapping of microbes (9). The MUC1 C-terminal (MUC1-C) subunit functions as an oncoprotein by promoting signaling pathways that control cell proliferation and survival (10, 11). Of relevance to breast cancer, MUC1-C is known to associate with ERα upon stimulation by 17β-estradiol (E2; ref. 12). MUC1-C binds directly to the ERα DNA-binding domain and stabilizes ERα by blocking its ubiquitination and degradation. MUC1-C also enhances ERα promoter occupancy, increases recruitment of coactivators, and stimulates ERα-mediated transcription (12). In contrast, loss of MUC1 confers inhibition of ERα translational activity and cell proliferation (13). Clinically, MUC1 expression correlates with ERα levels in breast tumors and its overexpression is associated with resistance to tamoxifen therapy (6, 12, 14). Its role in AI resistance, however, is not known.

MUC1 expression has previously been shown to be upregulated by IFNγ in cell lines derived from breast carcinomas (15). Most recently, we published that the IFNα signaling pathway is hyperactivated in AI-resistant breast cancer cells and that these resistant cells constitutively overexpress several IFNα-stimulated genes (ISG), which promote an aggressive phenotype (16). Of note, we found that IFN-induced transmembrane protein 1 (IFITM1), a critical downstream target of the IFNα signaling pathway, was overexpressed in AI-resistant cells and AI-resistant breast tumors. Targeting IFITM1 expression induced cell death in the resistant cells in vitro and in vivo (16, 17). Of particular interest was the observation that estrogen treatment suppressed IFITM1 expression in AI-resistant MCF-7:5C cells, which was associated with cell death (16). IFITM1 expression is regulated by IFNs through activation of the JAK/STAT signaling pathway. Previous studies have reported that MUC1 can interact with STAT1 protein in breast cancer cells to maintain activating phosphorylation and that this interaction is associated with poor overall survival (18). We hypothesize that MUC1 augments IFITM1 expression in AI-resistant MCF-7:5C cells and that crosstalk between MUC1/IFITM1 promotes an aggressive phenotype in the resistant cells.

In this study, we investigated the interaction between MUC1 and IFITM1 in AI-resistant MCF-7:5C cells and AI-sensitive T-47D and MCF-7 cells. We examined whether estrogen treatment or inhibition of JAK/STAT signaling can disrupt this crosstalk and thus contribute to death in the resistant cells. We found that siRNA knockdown of MUC1 markedly reduced IFITM1 expression in AI-resistant MCF-7:5C cells, which was associated with cell death, and that treatment with E2 further enhanced this phenotype. In addition, we found that MUC1, STAT1, and STAT2 directly bound to the IFITM1 promoter induce its expression through modulation of JAK/STAT1 signaling in the resistant cells. Finally, using tissue microarray of 94 ER+ breast tumors and in silico analysis of 4,151 patient samples in publicly available databases, we demonstrated that high MUC1 and IFITM1 expression correlate with poor patient outcome and overall survival. Together, these findings demonstrate a critical role for MUC1/IFITM1 crosstalk in promoting the aggressiveness of AI-resistant breast cancer and they suggest that inhibiting MUC1 or IFITM1 expression either with low-dose estrogen therapy or JAK/STAT inhibition has potential as a viable treatment option for AI-resistant breast cancer.

Cell lines

The MCF-7 cell line was obtained from Dr. V. Craig Jordan (University of Texas MD Anderson Cancer Center, Houston, TX) and maintained in RPMI1640 medium supplemented as described previously (19). The long-term estrogen-deprived human breast cancer cell line; MCF-7:5C was cloned from parental MCF-7 cells following long-term (>12 months) culture in estrogen-free medium composed of phenol red–free RPMI1640 supplemented with 10% FBS treated three times with dextran-coated charcoal (Thermo Fisher Scientific, catalog no. 7440-44-0; ref. 19). The T-47DA:18 (20) cell line (hereafter referred to as T-47D) was derived from T-47D (21) cells originally obtained from ATCC and cultured in the same media as the MCF-7 cells. All cell lines were cultured at 37°C under 5% CO2.

Western blotting

Cells were seeded in 6-well plates and allowed to acclimatize overnight. Following 24-hour treatment as indicated with 1 nmol E2 (Sigma, catalog no. E8875), 72-hour treatment with siMUC1 (Santa Cruz Biotechnology, catalog no. SC-37266), or 48-hour treatment with ruxolitinib/Jakafi (Rux; Cayman Chemical, catalog no. 11609) cells were harvested. Total cell lysate was separated by gel electrophoresis and transferred to polyvinylidene difluoride membrane as described previously (17, 22). Target proteins were detected using either anti-MUC1 (Santa Cruz Biotechnology, catalog no. SC-7313), anti-ERα (Santa Cruz Biotechnology, catalog no. SC-544), anti-P-STAT1 (Santa Cruz Biotechnology, catalog no. SC-8394), anti-P-STAT2 (Cell Signaling Technology, catalog no. 88410S), anti-STAT1 (Santa Cruz Biotechnology, catalog no. SC-464), anti-STAT2 (Santa Cruz Biotechnology, catalog no. SC-514193), anti-PARP-1 (Santa Cruz Biotechnology, catalog no. SC-8007), anti-IFITM1 (Santa Cruz Biotechnology, catalog no. SC-374026), or anti-β-actin (Cell Signaling Technology, catalog no. 3700S) antibodies. The appropriate horseradish peroxidase (HRP)-conjugated secondary antibody (Cell Signaling Technology, catalog no. 7076S and catalog no. 7074S) was applied and the positive bands were detected on autoradiography film as described previously (17, 22).

RNA isolation and real-time PCR

Cells were seeded in 6-well plates and allowed to acclimatize overnight. Following 24-hour treatment with 1 nmol E2 or siMUC1, RNA was isolated using the RNeasy Mini kit (Qiagen, catalog no. 74104). First-strand cDNA synthesis was performed from 3 μg total RNA using M-MLV Reverse Transcriptase (Invitrogen, catalog no. 28025-013) on a Bio-Rad MyCycler. RT-PCR was conducted using the ViiA 7 Real-Time PCR system (Applied Biosystems) and SYBR Green Reagent (Applied Biosystems, catalog no. 4367659) with 25 pmol primers specific for human MUC1 (sense: 5′-ACCTACCATCCTATGAGCGAG-3′; antisense: 5′-GGTTTGTGTAAGAGAGGCTGC-3′), IFITM1 (sense: 5′-GGATTTCGGCTTGTCCCGAG-3′; antisense: 5′- CCATGTGGAAGGGAGGGCTC-3′), ERα (sense: 5′-AAGAGGGTGCCAGGCTTTGT-3′; antisense: 5′-CAGGATCTCTAGCCAGGCACAT-3′), STAT1 (sense: 5′- CCGCCATGTTTACAGCAGAT-3′; antisense: 5′-GTCCCCTAGGACCTCCTCAT -3′), and STAT2 (sense: 5′-GCAGCACCATTTGCGGAA -3′; antisense: 5′-ACAGGTGTTTCGAGAACTGGC-3′). PUM1 was used as the internal control (sense: 5′-TCACCGAGGCCCCTCTGAACCCTA-3′; antisense: 5′-GGCAGTAATCTCCTTCTGCATCC T-3′). Relative mRNA expression level was determined as the ratio of the signal intensity to that of PUM1 using the formula: 2−ΔCt. When cells were treated, fold change in gene expression was normalized to PUM1 and then compared with the untreated value for that cell line using the formula: 2−ΔΔCt.

Immunofluorescent staining

IF was performed as described previously (17). Briefly, cells were seeded in 2-well slides and allowed to acclimatize overnight. Following 24-hour treatment with 1 nmol E2, cells were fixed with methanol. Because of the use of mouse antibodies on mouse tissue, blocking and antibody dilution were performed using the Mouse on Mouse (MOM) Kit (Vector Labs, catalog no. FMK-2201) following manufacturer's instructions. Sections were stained using antibodies against anti-MUC1 (Santa Cruz Biotechnology, catalog no. SC-7313), anti-ERα (Santa Cruz Biotechnology, catalog no. SC-544). Secondary antibodies were FITC (Santa Cruz Biotechnology, catalog no. SC-2359) or Texas Red (Santa Cruz Biotechnology, catalog no. SC-2781) conjugated. Slides were visualized on a Leica TCS SPE confocal microscope in the Confocal Imaging Core at The University of Kansas Medical Center (Kansas City, KS). Images were collected and analyzed using the Leica LAS AF Lite software (Leica Biosystems).

Small interfering RNA (siRNA) transfections

Cells were transiently transfected with siRNA for MUC1 (Santa Cruz Biotechnology, catalog no. sc-37266) or a scrambled negative control (Santa Cruz Biotechnology, catalog no. sc-37007). The MUC1 and control siRNAs were pools of three target-specific 20–25 nt siRNAs. Cells were seeded the night before transfection and allowed to reach 60% confluence by the time of transfection. Twenty nanomoles of each siRNA was introduced using Lipofectamine 2000 (Invitrogen, catalog no. 11668-027) in OptiMEM Reduced-Serum Medium (Gibco, catalog no. 11058-021) according to the manufacturer's instructions. After overnight incubation, the transfection mixture was replaced with normal culture medium, containing E2 only where indicated.

Cell counting for proliferation

Cells were assayed for viability and proliferation in 24-well plates in triplicate in estrogen-free medium, supplemented with E2 only where indicated for 72 hours. After treatment, cells were counted by Trypan blue (Sigma, catalog no. T8154) exclusion direct cell counts.

TUNEL staining

Cells were seeded in 2-well slides and allowed to acclimatize overnight. After the indicated treatment conditions, TUNEL staining was conducted using the Click-iT Plus TUNEL Assay Kit (Invitrogen, catalog no. C10618) following manufacturer's instructions. Tumor samples were deparaffinized by clearance in xylene, rehydrated through graded ethanol series, and then subjected to TUNEL staining as above. The average TUNEL intensity was quantified using the red color channel on ImageJ software for a minimum of three images.

Dual luciferase reporter assay

For IFITM1 promoter assays, 0.8 μg of plasmid DNA and the pRL CMV Renilla vector were used as described previously (22). For analysis of IFITM1 promoter activity, the pGL3 plasmid with the first 750 nucleotides of the IFITM1 promoter inserted (pGL3-IFITM1 [−750/−1]), was used (22). The pGL3-Basic-IRES was a kind gift from Joshua Mendell (Addgene, catalog no. 64784; ref. 23). After 24 hours, transfection reagent was replaced with normal cell culture media containing Rux where indicated. Luciferase and Renilla activities were measured 24 hours later using the Dual-Luciferase Reporter Assay Kit (Promega, catalog no. E1910) according to the manufacturer's instructions on a BioTek Synergy 4 microplate reader using the Gen 5 data analysis software (BioTek Instruments).

Coimmunoprecipitation

For coimmunoprecipitation (coIP) experiments, cell lysates were collected in RIPA buffer (Thermo Scientific, catalog no. 89901) with protease inhibitor cocktail (Roche Diagnostics, catalog no. 11836-153-001) and phosphatase inhibitor (Sigma, catalog no. P0044), and sonicated on ice. Cell lysates containing equivalent protein concentrations (5 μg) were precleared with 50:50 Protein A/G–coated magnetic beads (Invitrogen, catalog no. 10001D/catalog no. 10003D) then incubated overnight at 4°C with 2 μg appropriate antibody or control IgG. 50:50 Protein A/G–coated magnetic beads were then added for the final 1 hour of incubation time. Immune complexes were washed three times with PBS, resuspended in Laemmli sample buffer containing dithiothreitol and β-mercaptoethanol (Invitrogen, catalog no. NP0007), boiled for 5 minutes, and subjected to Western blotting analysis.

Chromatin immunoprecipitation assay

Chromatin immunoprecipitation (ChIP) was performed using the ChIP-IT Express Kit (Active Motif, catalog no. 53008) according to the manufacturer's instructions using sonication as the method for chromatin shearing. Lysates were immunoprecipitated (IP) overnight (18 hours) with the following antibodies anti-STAT1 (Santa Cruz Biotechnology, catalog no. SC-464), anti-STAT2 (Santa Cruz Biotechnology, catalog no. SC-514193), anti-MUC1 (Santa Cruz Biotechnology, catalog no. SC-7313) or an equal amount of mouse (Santa Cruz Biotechnology, catalog no. SC-2762) or rabbit IgG (Santa Cruz Biotechnology, catalog no. SC-2027). Resulting DNA was analyzed using qPCR as described previously (24), and data are represented as a percentage of input DNA.

Orthotopic cell line transplantation

Eight- to 10-week-old virgin female athymic nude (nu/nu) or NOD-SCID IL2Rgammanull (NSG) mice were purchased from Jackson Laboratories and used for in vivo experiments. Animal experiments were conducted following protocols approved by the University of Kansas School of Medicine Animal Care and Use (ACUP#: 2016-2341). Cells were suspended in 50:50 PBS/Matrigel (Corning, catalog no. 354262) and bilaterally injected into fourth mammary fat pads as described previously (25). A total of 3 × 106 cells were delivered per injection in a volume of 100 μL. Tumor volume was calculated weekly as described previously (25). When tumors reached a mean volume of 0.20 cm3, groups of 5–15 mice were randomly assigned to treatment groups. Where indicated, mice were administered 50 μg/kg body weight of ruxolitinib (LC Laboratories, catalog no. R6688) suspended in methylcellulose (Acros Organics, catalog no. 9004-67-5) by oral gavage every other day. Experiments were terminated after 43 days due to gastrointestinal side-effects in the treatment group. When using MCF-7 and MCF-7:5C cells for estrogen experiments, capsules containing 1:4 β-estradiol/cholesterol, which produce a mean serum estradiol level of approximately 80 pg/mL, were implanted subcutaneously as described previously (26).

IHC staining

IHC staining was performed after tissue deparaffinization by clearance in xylene and hydration through graded ethanol series as described previously (17). Sections were stained using primary human antibodies targeted against anti-IFITM1 (Santa Cruz Biotechnology, catalog no. SC-374026), anti-MUC1 (Santa Cruz Biotechnology, catalog no. SC-7313), anti-P-STAT1 (Santa Cruz Biotechnology, catalog no. SC-8394), and horseradish peroxidase–conjugated biotinylated secondary antibodies (Cell Signaling Technology, catalog nos. 7076S and 7074S). Immunoperoxidase signal was produced using 3,3′-Diaminobenzidine (DAB; Vector Laboratories, catalog no. SK-4100) and amplified using the Vectastain Elite ABC Kit (Vector Laboratories, catalog no. PK-4000). Tissue sections were counter stained using hematoxylin and mounted in xylene. Slides were imaged on a Nikon Eclipse 80i Upright Microscope in the Imaging Core of The University of Kansas Medical Center (Kansas City, KS).

Tumor microarray analysis

A tumor microarray of 104 patient samples as described previously (17) was analyzed for MUC1 and IFITM1 expression and scored by staining intensity by three independent analysts blinded to the patient details. The MUC1 and IFITM1 antibodies were validated using samples of human colorectal cancer as described previously. The MUC1 antibody was also validated using normal human breast tissues samples taken from reduction mammoplasties at KUMC (Supplementary Fig. S7A). Clinicopathologic data including age, race, clinical stage, and Her2 staining were published previously (17). Two cores from each tumor were analyzed and averaged in an attempt to account for the impact of tumor heterogeneity on protein expression. MUC1 staining intensity was quantified manually on a scale of 0–3 where 0 means no staining, 1+ is faint staining, 2+ is moderate staining and 3+ is strong staining (Supplementary Fig. S7B). Cores were scored separately for IFITM1 and MUC1 expression by three independent individuals before accessing patient medical records. Any discrepancies were resolved by group consensus. Demonstration of IFITM1 staining intensity and final distribution of IFITM1 expression was published previously on the same scale (17). For Kaplan–Meier plotting, a score of 0 to 1 was labeled as negative (−) and 2 to 3 as positive (+). Histopathologic parameters were extracted from patient pathology records. Clinical parameters including disease recurrence and survival data were obtained from electronic medical records for 94 of 104 patients.

E2 treatment reduces MUC1 levels and induces apoptosis in AI-resistant breast cancer cells

We previously reported the development of an AI-resistant breast cancer cell line, MCF-7:5C, which was derived from parental MCF-7 cells following long-term estrogen deprivation (19, 27). A unique phenotype of MCF-7:5C cells is that they express high levels of functional ERα but lack PgR and undergo apoptosis in the presence of 17β-estradiol (E2; refs. 27, 28). Our study begins by examining the effect of E2 treatment on MUC1 expression in AI-sensitive T-47D and MCF-7 cells and AI-resistant MCF-7:5C cells. First, we measured baseline expression of MUC1 and ERα protein and mRNA in all three cell lines and we found that MUC1 protein and mRNA levels were highest in T-47D cells compared with MCF-7 and MCF-7:5C cells; however, ERα protein and mRNA levels were highest in AI-resistant MCF-7:5C cells (Fig. 1A). Next, we confirmed that estrogen differentially regulates MUC1 expression in AI-sensitive T-47D and MCF-7 cells compared with AI-resistant MCF-7:5C cells. Estrogen starvation causes upregulation of ERα in the AI-sensitive cell lines and subsequent estrogen treatment (1 nmol/L E2) significantly increased MUC1 protein and mRNA level in MCF-7 and T-47D cells within 24 hours; however, in AI-resistant MCF-7:5C cells, estrogen treatment significantly reduced MUC1 protein and mRNA level within 12 hours (Fig. 1B; Supplementary Fig. S1A). Notably, in MCF-7:5C cells, MUC1 levels remained reduced up to 72 hours, which was followed by increased PARP cleavage, and reduced P-STAT1 and IFITM1 expression (Supplementary Fig. S1C). MUC1 downregulation in AI-resistant MCF-7:5C cells was confirmed with immunofluorescence (Fig. 1C). For comparison, similar experiments were performed with the pure antiestrogen ICI 182,780 and we found that ICI 182,780 decreased MUC1 levels in T-47D and MCF-7 cells, but increased MUC1 expression in AI-resistant MCF-7:5C cells (Fig. 1B), further suggesting differential hormonal regulation of MUC1 in AI-sensitive versus AI-resistant breast cancer cells.

Figure 1.

E2 treatment reduces MUC1 levels in AI-resistant breast cancer cells. A, Whole-cell lysates from 80% confluent T-47D, MCF-7, and MCF-7:5C cells were immunoblotted for MUC1 and ERα protein expression. MUC1 and ERα mRNA expression were determined by RT-PCR and are displayed relative to MCF-7 cells. Data represents three independent experiments run in triplicate. B, T-47D, MCF-7, and MCF-7:5C cells were treated for 0, 12, and 24 hours with 1 nmol β-estradiol (E2) or 1mmol ICI 182,780 (ICI) and immunoblotted for MUC1 and ERα expression. C, After 72 hours of E2 treatment, MUC1 and ERα protein expression in T-47D, MCF-7, and MCF-7:5C cells was determined using immunofluorescent staining (*, P < 0.05; **, P < 0.01).

Figure 1.

E2 treatment reduces MUC1 levels in AI-resistant breast cancer cells. A, Whole-cell lysates from 80% confluent T-47D, MCF-7, and MCF-7:5C cells were immunoblotted for MUC1 and ERα protein expression. MUC1 and ERα mRNA expression were determined by RT-PCR and are displayed relative to MCF-7 cells. Data represents three independent experiments run in triplicate. B, T-47D, MCF-7, and MCF-7:5C cells were treated for 0, 12, and 24 hours with 1 nmol β-estradiol (E2) or 1mmol ICI 182,780 (ICI) and immunoblotted for MUC1 and ERα expression. C, After 72 hours of E2 treatment, MUC1 and ERα protein expression in T-47D, MCF-7, and MCF-7:5C cells was determined using immunofluorescent staining (*, P < 0.05; **, P < 0.01).

Close modal

Loss of MUC1 alone and in combination with E2 treatment inhibits cell proliferation and induces apoptosis

To test whether MUC1 expression held any functional significance on the aggression of AI-resistant cells, we first knocked down MUC1 with three unique MUC1-specific siRNAs (Supplementary Fig. S1B) and a combined pool of all three siRNAs (Fig. 2A). Trypan blue direct cell counting confirmed that E2 treatment significantly enhanced the growth of MCF-7 and T-47D cells, but reduced the growth of AI-resistant MCF-7:5C cells (Fig. 2B), as previously reported (27). Notably, loss of MUC1 resulted in lower cell proliferation in AI-resistant MCF-7:5C cells; however, it did not alter the growth of AI-sensitive MCF-7 and T-47D cells (Fig. 2C). To further investigate the effects of MUC1 loss on AI-resistant MCF-7:5C cells we combined MUC1 knockdown with estrogen treatment. Phase contrast images demonstrated an increase in dead, floating cells (Fig. 2D) and TUNEL staining verified induction of apoptosis most significantly in the MUC1 knockdown plus estrogen-treated cells compared with siMUC1 or estrogen treatment alone (Fig. 2E).

Figure 2.

Loss of MUC1 alone and in combination with E2 treatment inhibits cell proliferation and induces apoptosis. A, Cell lysates from T-47D, MCF-7, and MCF-7:5C cells transfected with scrambled control siRNA (siCon) or a pool of three different MUC1 siRNA sequences (siMUC1) were immunoblotted for MUC1 expression. B, Cell proliferation was measured by Trypan blue exclusion cells counts after 72 hours of E2 treatment in T-47D, MCF-7, and MCF-7:5C cells. Data represent means ± SD from three experiments conducted in triplicate. C, Cell proliferation was measured by Trypan blue exclusion cells counts after 72 hours of siMUC1 transfection in T-47D, MCF-7, and MCF-7:5C cells. Data represent means ± SD from three experiments conducted in triplicate. D, MCF-7:5C cells were transiently transfected with siCon or siMUC1 alone or in combination with E2 and photos taken 72 hours after transfection. Phase contrast images were taken to visualize the number of floating cells indicating apoptosis. E, MCF-7:5C cells were transiently transfected with siCon or siMUC1 alone or in combination with E2 and measured for apoptosis by TUNEL staining. The intensity of TUNEL (red) staining was used to assess the effect of IFITM1 knockdown on cell survival and was quantified using ten separate images with Image J Software (right; *, P < 0.05; **, P < 0.01; ***, P < 0.001).

Figure 2.

Loss of MUC1 alone and in combination with E2 treatment inhibits cell proliferation and induces apoptosis. A, Cell lysates from T-47D, MCF-7, and MCF-7:5C cells transfected with scrambled control siRNA (siCon) or a pool of three different MUC1 siRNA sequences (siMUC1) were immunoblotted for MUC1 expression. B, Cell proliferation was measured by Trypan blue exclusion cells counts after 72 hours of E2 treatment in T-47D, MCF-7, and MCF-7:5C cells. Data represent means ± SD from three experiments conducted in triplicate. C, Cell proliferation was measured by Trypan blue exclusion cells counts after 72 hours of siMUC1 transfection in T-47D, MCF-7, and MCF-7:5C cells. Data represent means ± SD from three experiments conducted in triplicate. D, MCF-7:5C cells were transiently transfected with siCon or siMUC1 alone or in combination with E2 and photos taken 72 hours after transfection. Phase contrast images were taken to visualize the number of floating cells indicating apoptosis. E, MCF-7:5C cells were transiently transfected with siCon or siMUC1 alone or in combination with E2 and measured for apoptosis by TUNEL staining. The intensity of TUNEL (red) staining was used to assess the effect of IFITM1 knockdown on cell survival and was quantified using ten separate images with Image J Software (right; *, P < 0.05; **, P < 0.01; ***, P < 0.001).

Close modal

MUC1 stabilizes JAK/STAT signaling which is necessary for IFITM1 expression

We previously reported that IFITM1, an IFN-stimulated gene, is overexpressed in AI-resistant MCF-7:5C cells and that its overexpression promotes an aggressive phenotype. (16, 17) We therefore examined whether MUC1 crosstalk with JAK/STAT signaling plays a role in regulating IFITM1 expression in MCF-7:5C cells. In Fig. 3A and D, we demonstrate that knockdown of MUC1 with siRNA significantly reduced IFITM1 and P-STAT1 expression in AI-resistant MCF-7:5C cells with a less pronounced effect seen in AI-sensitive T-47D cells. Notably, the AI-sensitive MCF-7 cell line does not express P-STAT1, P-STAT2, or IFITM1. Next, we targeted JAK/STAT signaling using the JAK1/2 inhibitor, ruxolitinib, and found reduced P-STAT1 and IFITM1 expression at the protein and mRNA level in MCF-7:5C and T-47D cells (Fig. 3B and D). Notably, E2 treatment also reduced P-STAT2 and IFITM1 expression in both cell lines (Fig. 3C). Physical binding of MUC1 with P-STAT1 and P-STAT2 was confirmed by coimmunoprecipitation using ERα as a positive control. Under normal conditions (no treatment), MUC1 is bound to P-STAT1 and P-STAT2 in AI-resistant MCF-7:5C cells, but only P-STAT2 in T-47D cells. Inhibition of STAT phosphorylation with ruxolitinib markedly reduced this association in both cell lines (Fig. 3E). In addition, we confirmed the role of MUC1 in IFITM1 expression using the MUC1 inhibitor, GO-201, in T-47D, MCF-7, and MCF-7:5C cells (Supplementary Fig. S1D). These data support a role for MUC1 upstream of IFITM1 expression. Notably, inhibition of JAK/STAT signaling with ruxolitinib did not significantly enhance the antiproliferative effect of either ICI 182,780 or estrogen in the AI-sensitive or AI-resistant cells, respectively (Supplementary Fig. S2).

Figure 3.

MUC1 stabilizes JAK/STAT signaling which is necessary for IFITM1 expression. A, T-47D, MCF-7, and MCF-7:5C cells were transiently transfected with siCon or siMUC1 for 72 hours and immunoblotted for MUC1, P-STAT2, STAT2, P-STAT1, STAT1, and IFITM1 expression. B, T-47D, MCF-7, and MCF-7:5C cells were treated for 48 hours with Ruxolitinib (Rux) and immunoblotted for MUC1, P-STAT2, STAT2, P-STAT1, STAT1, and IFITM1 expression. C, T-47D, MCF-7, and MCF-7:5C cells were treated for 24 hours with E2 and immunoblotted for MUC1, P-STAT2, STAT2, P-STAT1, STAT1, and IFITM1 expression. D, T-47D and MCF-7:5C cells were transiently transfected with control or MUC1 siRNA or treated with ruxolitinib (Rux). After 24 hours, mRNA expression of IFITM1 was determined by RT-PCR. E, Whole-cell lysates from T-47D and MCF-7:5C cells were treated with vehicle (Con) or 10 μmol Rux for 8 hours were immunoprecipitated with anti-MUC1 antibody or rabbit IgG and samples immunoblotted for MUC1, P-STAT2, P-STAT1, and ERα protein interaction. F, Fixed T-47D and MCF-7:5C cell lysates were subjected to chromatin immunoprecipitation (ChIP) with antibodies against STAT1, STAT2, MUC1 or species-specific IgG control. As indicated, fixed whole-cell lysates from T-47D and MCF-7:5C cells treated with 10 μmol/L ruxolitinib for 24 hours. qPCR was performed on the isolated DNA using primers designed to amplify the ISRE regulatory regions. Recruitment of the indicated proteins to the ISRE site was compared to input and IgG DNA and displayed as mean ± SD of technical triplicates in two independent experiments. Fold change between control and treatment is displayed above the indicated protein (**, P < 0.01; ***, P < 0.001).

Figure 3.

MUC1 stabilizes JAK/STAT signaling which is necessary for IFITM1 expression. A, T-47D, MCF-7, and MCF-7:5C cells were transiently transfected with siCon or siMUC1 for 72 hours and immunoblotted for MUC1, P-STAT2, STAT2, P-STAT1, STAT1, and IFITM1 expression. B, T-47D, MCF-7, and MCF-7:5C cells were treated for 48 hours with Ruxolitinib (Rux) and immunoblotted for MUC1, P-STAT2, STAT2, P-STAT1, STAT1, and IFITM1 expression. C, T-47D, MCF-7, and MCF-7:5C cells were treated for 24 hours with E2 and immunoblotted for MUC1, P-STAT2, STAT2, P-STAT1, STAT1, and IFITM1 expression. D, T-47D and MCF-7:5C cells were transiently transfected with control or MUC1 siRNA or treated with ruxolitinib (Rux). After 24 hours, mRNA expression of IFITM1 was determined by RT-PCR. E, Whole-cell lysates from T-47D and MCF-7:5C cells were treated with vehicle (Con) or 10 μmol Rux for 8 hours were immunoprecipitated with anti-MUC1 antibody or rabbit IgG and samples immunoblotted for MUC1, P-STAT2, P-STAT1, and ERα protein interaction. F, Fixed T-47D and MCF-7:5C cell lysates were subjected to chromatin immunoprecipitation (ChIP) with antibodies against STAT1, STAT2, MUC1 or species-specific IgG control. As indicated, fixed whole-cell lysates from T-47D and MCF-7:5C cells treated with 10 μmol/L ruxolitinib for 24 hours. qPCR was performed on the isolated DNA using primers designed to amplify the ISRE regulatory regions. Recruitment of the indicated proteins to the ISRE site was compared to input and IgG DNA and displayed as mean ± SD of technical triplicates in two independent experiments. Fold change between control and treatment is displayed above the indicated protein (**, P < 0.01; ***, P < 0.001).

Close modal

Because inhibition of JAK/STAT signaling with ruxolitinib reduced IFITM1 expression in MCF-7:5C and T-47D cells at the mRNA level, we next assessed the effect of ruxolitinib on the IFITM1 promoter. We performed ChIP assay to assess the binding of MUC1, STAT1, and STAT2, to the IFN-stimulated response element (ISRE) located at the -1 position in the IFITM1 promoter. Baseline levels showed greater binding of MUC1, STAT1, and STAT2 to the IFITM1 promoter in AI-resistant MCF-7:5C and AI-sensitive T-47D cells compared with MCF-7 cells (Supplementary Fig. S3B). Upon treatment with ruxolitinib, binding of MUC1, STAT1, and STAT2 was significantly decreased compared with controls in both MCF-7:5C and T-47D cells indicating a role for these proteins in promoting IFITM1 expression (Fig. 3F). Dual luciferase reporter assay supported these findings (Supplementary Fig. S3A). Inhibition of IFITM1 expression had no effect on upstream MUC1 expression (Supplementary Fig. S3C). These studies found a more prominent role for P-STAT2 in driving IFITM1 expression in AI-sensitive T-47D cells, while indicating that IFITM1 expression in AI-resistant MCF-7:5C cells is more dependent on P-STAT1. Overall, these data validate the importance of MUC1 in transcriptional regulation of IFITM1 expression.

E2 treatment inhibits the growth of AI-resistant MCF-7:5C tumors and suppresses P-STAT1 and IFITM1 levels in vivo

Because our in vitro data indicated that estrogen (E2) treatment significantly reduced MUC1 and IFITM1 expression in AI-resistant MCF-7:5C cells (Fig. 3C) and that this reduction was associated with cell death (Fig. 2), we next determined the effect of E2 treatment on AI-resistant tumor growth and MUC1 and IFITM1 expression in vivo. For this experiment, AI-sensitive MCF-7 cells and AI-resistant MCF-7:5C cells (3 × 106) were bilaterally injected into the fourth mammary fat pad of NSG mice. As expected, treatment with E2 significantly enhanced the growth of MCF-7 tumors; however, it markedly reduced the growth of AI-resistant MCF-7:5C tumors (Fig. 4A and B). TUNEL staining of MCF-7:5C tumors revealed that the reduction in tumor growth was due to an increase in apoptosis (Fig. 4C). Similar experiments were performed using athymic nude mice and the results were comparable with those seen with NSG mice (Supplementary Fig. S4A–S4C). Of importance, IHC staining of MCF-7:5C tumors revealed high expression of MUC1, P-STAT1, and IFITM1 under control conditions (-E2); however, treatment with E2 completely reduced the expression of all three proteins (Fig. 4D). A similar finding was observed for MCF-7:5C tumors generated in Nude mice (Supplementary Fig. S5A and S5B).

Figure 4.

E2 treatment reduces tumor size and reduces MUC1, P-STAT1, and IFITM1 levels in AI-resistant cells in vivo. A, 3 million MCF-7 and MCF-7:5C cells were bilaterally injected into the fourth mammary fat pad of ovariectomized female NSG mice with and without E2 capsule implants respectively. 35 days post-injection, mice were randomized into two groups with half of the E2 capsule implants removed from MCF-7 injected mice and half of the MCF-7:5C injected mice received estrogen capsule implants. Tumors were measured using digital calipers and measurements were used to calculate tumor volume (mm3) over time. B, At the end of the experiment tumors were excised and final tumor volume was measured. In the MCF-7 injected mice, 12 tumors were collected from the -E2 group and 11 tumors were collected from the +E2 group. In the MCF-7:5C injected mice, 14 tumors were collected from the -E2 group and 16 tumors were collected from the +E2 group. C, Cell death was analyzed by TUNEL staining and quantified by ImageJ software. D, Hematoxylin and Eosin (H+E), IgG (negative staining control) and MUC1, P-STAT1, and IFITM1 staining in NSG mice from untreated and E2-treated mice (*, P < 0.05; **, P < 0.01; and ***, P < 0.001).

Figure 4.

E2 treatment reduces tumor size and reduces MUC1, P-STAT1, and IFITM1 levels in AI-resistant cells in vivo. A, 3 million MCF-7 and MCF-7:5C cells were bilaterally injected into the fourth mammary fat pad of ovariectomized female NSG mice with and without E2 capsule implants respectively. 35 days post-injection, mice were randomized into two groups with half of the E2 capsule implants removed from MCF-7 injected mice and half of the MCF-7:5C injected mice received estrogen capsule implants. Tumors were measured using digital calipers and measurements were used to calculate tumor volume (mm3) over time. B, At the end of the experiment tumors were excised and final tumor volume was measured. In the MCF-7 injected mice, 12 tumors were collected from the -E2 group and 11 tumors were collected from the +E2 group. In the MCF-7:5C injected mice, 14 tumors were collected from the -E2 group and 16 tumors were collected from the +E2 group. C, Cell death was analyzed by TUNEL staining and quantified by ImageJ software. D, Hematoxylin and Eosin (H+E), IgG (negative staining control) and MUC1, P-STAT1, and IFITM1 staining in NSG mice from untreated and E2-treated mice (*, P < 0.05; **, P < 0.01; and ***, P < 0.001).

Close modal

Blockade of JAK/STAT signaling using ruxolitinib reduces the growth of AI-resistant MCF-7:5C tumors in vivo

Because our in vitro data revealed that P-STAT1 and IFITM1 expression were markedly elevated in AI-resistant MCF-7:5C cells compared with parental MCF-7 cells and that treatment with the JAK inhibitor ruxolitinib markedly reduced P-STAT1 and IFITM1 expression in these cells (Fig. 5A), we next determined whether ruxolitinib has therapeutic benefits in vivo. MCF-7 and MCF-7:5C cells were bilaterally injected into the fourth mammary fat pad of NSG mice and once the tumors reached approximately 200 mm3 in size the mice were randomized to either the control group (-Rux) or the ruxolitinib-treated group (50 mg/kg). Our data showed that ruxolitinib treatment significantly reduced the growth of AI-resistant MCF-7:5C tumors within the first 2 weeks (Fig. 5B); however, it did not alter the growth of MCF-7 tumors (Fig. 5B) or T-47D tumors (data not shown) throughout the duration of the study. TUNEL staining showed increased apoptosis after ruxolitinib treatment in MCF-7:5C tumors (Fig. 5C), but not in MCF-7 tumors (Supplementary Fig. S6) and this increase in cell death was associated with a reduction in P-STAT1 and IFITM1 expression, as confirmed by Western blot analysis (Fig. 5D) and IHC staining (Fig. 5E).

Figure 5.

Ruxolitinib (Rux) treatment reduces tumor size and reduces MUC1, P-STAT1, and IFITM1 levels in AI-resistant cells in vivo. A, Whole-cell lysates from MCF-7 and MCF-7:5C cells treated with 10 μmol ruxolitinib Rux for 24 hours were immunoblotted for P-STAT1, STAT1, and IFITM1 protein expression. B, 3 million MCF-7 or MCF-7:5C cells were injected into the fourth mammary fat pad of female NSG mice. After 32 days of tumor growth, mice were randomized to treatment groups and half were given 50 mg/kg ruxolitinib by oral gavage every other day. At the end of the experiment tumors were excised and tumor volumes were determined by digital caliper measurement. In the MCF-7 injected mice, 8 tumors were collected from the −Rux group and 7 tumors were collected from the +Rux group. In the MCF-7:5C injected mice, 5 tumors were collected from the −Rux group and 6 tumors were collected from the +Rux group. C, Tumors were subjected to TUNEL staining to quantify apoptotic cells. The intensity of TUNEL (red) staining was quantified using ten separate images with ImageJ Software (bottom). D, Tumor expression of P-STAT1, STAT1 and IFITM1 after 21 days of treatment was determined by immunoblot. E, Tumors were fixed, embedded in paraffin, and sectioned onto glass slides. H+E staining revealed tumor architecture and MUC1, IFITM1, and P-STAT1 expression was determined by IHC (*, P < 0.05; **, P < 0.01; and ***, P < 0.001).

Figure 5.

Ruxolitinib (Rux) treatment reduces tumor size and reduces MUC1, P-STAT1, and IFITM1 levels in AI-resistant cells in vivo. A, Whole-cell lysates from MCF-7 and MCF-7:5C cells treated with 10 μmol ruxolitinib Rux for 24 hours were immunoblotted for P-STAT1, STAT1, and IFITM1 protein expression. B, 3 million MCF-7 or MCF-7:5C cells were injected into the fourth mammary fat pad of female NSG mice. After 32 days of tumor growth, mice were randomized to treatment groups and half were given 50 mg/kg ruxolitinib by oral gavage every other day. At the end of the experiment tumors were excised and tumor volumes were determined by digital caliper measurement. In the MCF-7 injected mice, 8 tumors were collected from the −Rux group and 7 tumors were collected from the +Rux group. In the MCF-7:5C injected mice, 5 tumors were collected from the −Rux group and 6 tumors were collected from the +Rux group. C, Tumors were subjected to TUNEL staining to quantify apoptotic cells. The intensity of TUNEL (red) staining was quantified using ten separate images with ImageJ Software (bottom). D, Tumor expression of P-STAT1, STAT1 and IFITM1 after 21 days of treatment was determined by immunoblot. E, Tumors were fixed, embedded in paraffin, and sectioned onto glass slides. H+E staining revealed tumor architecture and MUC1, IFITM1, and P-STAT1 expression was determined by IHC (*, P < 0.05; **, P < 0.01; and ***, P < 0.001).

Close modal

Clinical significance of high MUC1 and IFITM1 coexpression in ER-positive breast cancer

Normal breast tissue from 6 reduction mammoplasties and tumor samples from a tumor microarray of 94 patients with ER+ breast cancer were analyzed by IHC and scored according to staining intensity (Supplementary Fig. S7). χ2 analysis confirmed a strong correlation between MUC1 and IFITM1 expression (Fig. 6A). We have previously reported an association between recurrence-free survival and overall survival based on IFITM1 expression alone in this patient cohort (17). Here, MUC1 intensity alone did not predict patient outcomes (Fig. 6B and C); however, positive coexpression of both MUC1 and IFITM1 significantly correlated with recurrence of disease and overall survival (Fig. 6D and E). We sought to verify this correlation in a larger cohort of breast cancers and mined data from BCCRC xenograft, METABRIC, TCGA, and TCGA 2015 datasets. We found that MUC1 and IFITM1 were coexpressed in 13%–27% of breast tumors (Fig. 7A) and that their coexpression was significantly associated with decreased recurrence-free and overall survival (Fig. 7B). These findings suggest that the presence of MUC1 and IFITM1 in breast tumors confers a poor prognosis for patients.

Figure 6.

Clinical significance of MUC1 and IFITM1 co-expression in ER-positive breast cancer. A-E, A tumor microarray of 104 patient samples was analyzed for MUC1 and IFITM1 expression and scored for staining intensity by three independent analysts. Low expression was given a rating of 0–1, intermediate expression was given a rating of 2, and high expression was given a rating of 3 and score distributions are displayed (A). MUC1 intensity score is corelated with breast cancer recurrence (B) and overall survival (C). An intensity score of (0-1) and (2–3) indicates positive (+) or negative (−) expression of MUC1 and IFITM1. Recurrence-free survival (D) and overall survival (E) of patients based on IFITM1 and MUC1 expression is displayed using Kaplan–Meier plots. Statistical analysis was performed using a log-rank Mantel–Cox test indicating significant difference in the curves.

Figure 6.

Clinical significance of MUC1 and IFITM1 co-expression in ER-positive breast cancer. A-E, A tumor microarray of 104 patient samples was analyzed for MUC1 and IFITM1 expression and scored for staining intensity by three independent analysts. Low expression was given a rating of 0–1, intermediate expression was given a rating of 2, and high expression was given a rating of 3 and score distributions are displayed (A). MUC1 intensity score is corelated with breast cancer recurrence (B) and overall survival (C). An intensity score of (0-1) and (2–3) indicates positive (+) or negative (−) expression of MUC1 and IFITM1. Recurrence-free survival (D) and overall survival (E) of patients based on IFITM1 and MUC1 expression is displayed using Kaplan–Meier plots. Statistical analysis was performed using a log-rank Mantel–Cox test indicating significant difference in the curves.

Close modal
Figure 7.

MUC1 and IFITM1 co-occur in breast cancer. A and B,In silico analysis of tumor gene expression using the cbioportal.org database for deposited breast tumor expression for both MUC1 and IFITM1 expression resulted in four datasets with significant overexpression of both proteins: BCCRC Xenograft, METABRIC, TCGA, and TCGA 2015. Fisher exact test on each of the four datasets was conducted to determine likelihood of co-occurrence of MUC1 and IFITM1 in breast cancer. All datasets indicated association towards co-occurrence of both MUC1 and IFITM1 in tumors and two demonstrated a statistically significant relationship. Log OR > 0: association toward cooccurrence, log OR ≤0: association towards mutual exclusivity, P < 0.05: significant association.

Figure 7.

MUC1 and IFITM1 co-occur in breast cancer. A and B,In silico analysis of tumor gene expression using the cbioportal.org database for deposited breast tumor expression for both MUC1 and IFITM1 expression resulted in four datasets with significant overexpression of both proteins: BCCRC Xenograft, METABRIC, TCGA, and TCGA 2015. Fisher exact test on each of the four datasets was conducted to determine likelihood of co-occurrence of MUC1 and IFITM1 in breast cancer. All datasets indicated association towards co-occurrence of both MUC1 and IFITM1 in tumors and two demonstrated a statistically significant relationship. Log OR > 0: association toward cooccurrence, log OR ≤0: association towards mutual exclusivity, P < 0.05: significant association.

Close modal

While AIs remain the standard of care for patients with ER-positive breast cancer, approximately 30% of patients will develop resistance to these agents. Mechanisms of resistance to AIs have been extensively studied, yet few novel, targetable regulators of AI resistance have been discovered and used clinically. Resistance can occur through ER-independent signaling, upregulated growth factor signaling, and changes to cell-cycle regulators and apoptotic mechanisms (29). We have previously reported that JAK/STAT signaling drives IFITM1 expression in our AI-resistant model, which directly promotes the aggressive phenotype of these cells. We found that targeting IFITM1 causes the death of AI-resistant cells in vivo and in vitro (16, 17). In this study, we demonstrate for the first time that MUC1 augments the underlying JAK/STAT signaling, enhances IFITM1 expression, and promotes an aggressive phenotype in AI-resistant breast cancer cells. Clinically, we found that coexpression of MUC1 and IFITM1 positively correlated with aggressive disease and resistance to therapy. Notably, in vivo studies using NSG mice demonstrated that treatment with ruxolitinib, a JAK/STAT inhibitor, and estrogen (E2) effectively reduced the growth of AI-resistant breast tumors with no significant effect seen for AI-sensitive tumors. Overall, these findings suggest that targeting MUC1 and the JAK/STAT signaling pathway simultaneously might be a novel strategy to treat a subset of patients with AI-resistant ER+ breast cancer.

We first investigated the correlation between MUC1 and IFITM1 in promoting breast cancer aggressiveness. MUC1, a known oncoprotein, is overexpressed in many types of cancer, including breast (30). Its C-terminus is a coactivator in many tumor-promoting signaling pathways in both ER-positive and -negative breast cancers (31, 32). Here, we elucidate a role for MUC1 in augmenting IFNα-mediated JAK/STAT signaling in a subset of AI-resistant tumors. Following constitutive activation of the type 1 IFNα receptor (IFNAR) by IFNα, MUC1 binds to STAT1 to maintain phosphorylation and regulate the expression of IFN-stimulated genes (ISG), including IFITM1 (18). Previously, our laboratory has shown that IFITM1 drives the aggressive phenotype of AI-resistant breast cancer (17). The IFITM1 promoter contains both IFN-stimulated response elements (ISRE) and STAT response elements. We show for the first time through in vitro (Fig. 3) and in vivo (Fig. 5) studies that MUC1 interacts with STAT1/2 to drive expression of IFITM1 and promote an aggressive phenotype in AI-resistant cells. Specifically, we found that MUC1 enhances STAT binding to the IFITM1 promoter through maintenance of STAT activation, which is consistent with reports of MUC1 cooperativity with STAT1 signaling (ref. 18; see visual overview). In addition, MUC1 may interact with chromatin remodelers and recruit coactivators to cause upregulation of IFITM1 expression, as it has a known role in enhancing recruitment of CBP, EzH2, c-Src, and DNMT (33–36). Future research should investigate the role of MUC1 in regulating SWI/SNF complexes, which are known to regulate ISG expression, including IFITM1 (37). Downregulation of MUC1 abrogates IFITM1 protein and mRNA expression (Fig. 3). MUC1 also associates with the IFITM1 ISRE promoter element in conjunction with STAT1/2 (Fig. 3), playing a role in constitutive activation of IFITM1 to maintain aggressiveness and AI resistance. IFITM1 expression causes tumor progression by downregulating the CDK inhibitor p21 (17). p21 increases the expression of proteins involved in cell progression and proliferation, like PCNA, Ki67, cyclin D1 and E, pRb, and E2F1 (38). The key role of MUC1 in regulating IFITM1 is of primary importance in breast cancer as high coexpression of MUC1 and IFITM1 together may predict patient outcomes, demonstrated through tumor microarray analysis of patients with ER-positive breast cancer (Fig. 6). Not only does IFITM1 expression predict overall outcome (17), as previously reported, but high expression of MUC1 and IFITM1 was associated with disease recurrence and poor overall survival in large patient datasets, emphasizing the opportunity to use MUC1 and IFITM1 as biomarkers of aggressive disease. Overexpression of MUC1 and IFITM1 was also examined through data mining and MUC1 and IFITM1 were found to co-occur in 13% to 27% of breast cancers from four distinct databases (Fig. 7). This resistance-driving crosstalk led us to further probe the interaction between MUC1 and IFITM1 and investigate strategies for targeting these tumors in the clinic.

Mechanistically, we uncovered a role for MUC1 in inducing IFITM1 expression and found that disrupting this interaction leads to cell death (Visual Overview). Notably, we observed differential regulation of MUC1 by estrogen in AI-resistant cells compared to AI-sensitive cells. In AI-sensitive cells, ERα is only active upon E2 treatment, which then upregulates MUC1 expression. MUC1 is influenced by estrogen treatment and ERα expression due to estrogen response elements (ERE) within its promoter. In addition, MUC1 can directly bind ERα which prevents its ubiquitination and degradation (30). This interaction is enhanced in the presence of estrogen and promotes the binding of ERα with its coactivators (12). Loss of MUC1 confers inhibition of ERα translational activity and cell proliferation (13). In contrast, treatment of AI-resistant cells and tumors with E2 results in MUC1 suppression and induction of death (Figs. 1 and 4). This dichotomy may tie back into initial mechanisms of AI resistance.

It is well-established that estrogen promotes breast cancer development and progression; however, preclinical data from our laboratory (16, 25, 27) and other investigators have previously shown that estrogen is also capable of inducing cell death in ER-positive breast cancers that have acquired resistance to estrogen deprivation (i.e., AI-resistant breast cancer). Clinically, low-dose estrogen therapy (6 mg/daily E2) has also been shown to be effective in treating a subset of postmenopausal women (30%–40%) with ER-positive AI-resistant breast cancer (39, 40). The ability of estrogen to induce death in AI-resistant cells is partly attributed to activation of the mitochondrial death pathway. However, our current study highlights a novel pathway in which estrogen also downregulates MUC1 expression, which reduces IFITM1 expression in AI-resistant cells. The loss of MUC1 and IFITM1 leads to cell death in AI-resistant MCF-7:5C cells but not AI-sensitive MCF-7 and T47D cells. AI-resistant cells can rely on ligand-independent signaling where ERα is constitutively active and promotes upregulation of MUC1. Estrogen may interfere with ERα stability, which in turn may abrogate the transcription and/or stabilization of MUC1 and its rate of turnover. Alternatively, it is possible that E2 treatment may disrupt ERα-coactivator binding at the MUC1 promoter thereby reducing MUC1 and consequently IFITM1 expression in the resistant cells (41). MUC1 can also be regulated by other nuclear hormone receptors, including PgR and the androgen receptor (AR), which may be inhibited by activation of ERα (42–44). Uncovering this mechanism allowed us to isolate MUC1 and consequently, IFITM1, as targets of E2-induced cell death in AI-resistant cells. MUC1 can also be constitutively overexpressed in ER-negative breast cancer through amplification of the MUC1 gene and is associated with patient prognosis (31, 45). Recently, high MUC1 expression has been associated with targetable CD274/PDL1 gene activation, PD-L1 expression, immune evasion and epithelial to mesenchymal transition in ER-negative breast cancer (32).

In addition to our in vitro data, we also investigated the in vivo effects of estrogen treatment on AI resistance and IFITM1 expression. Estrogen acts paradoxically in AI-resistant cells and instead of promoting their growth, induces death. We observed this effect in AI-resistant MCF-7:5C cells upon treatment with E2 (Fig. 2D and E). Perhaps this occurs because of the downregulation of MUC1 and IFITM1 levels after E2 treatment (Fig. 3A), which leads to an increase in p21, a decrease in phospho-p21, and cell death (17). Significantly, E2 treatment caused tumor regression and downregulation of these proteins in vivo, exclusively in the AI-resistant cell line (Fig. 4).

As an alternative therapy to estrogen treatment or potentially as a combination therapy, we investigated the effects of ruxolitinib on MUC1 and IFITM1 expression and AI resistance. Ruxolitinib is a small-molecule inhibitor of JAK1 and JAK2 and so widely dampens the JAK/STAT signaling pathway. Ruxolitinib is currently under clinical investigation in combination treatment with trastuzumab in metastatic HER2+ breast cancer (32), with chemotherapy in triple-negative inflammatory breast cancer (46), and with pembrolizumab in metastatic triple-negative breast cancer (TNBC; ref. 47). However, a phase II clinical trial of ruxolitinib in patients with metastatic TNBC demonstrated that patients treatment was ineffective in this subpopulation of patients as a single agent (48). These patients were selected on the basis of high P-STAT3 expression; alternatively, more successful results may have been achieved by selecting patients with high MUC1 and IFITM1 expression. The dichotomy of treatment efficacy between HER2+ and triple-negative breast cancers and our research in AI-resistant breast cancer cells may be due to the effects of ruxolitinib on downregulating IFITM1, a linchpin in maintaining AI resistance, specific to ER-positive breast cancer. Our data indicate that ruxolitinib, as a single agent or in combination with other therapies (i.e., estrogen), could be effective in treating patients with ER-positive AI-resistant or IFITM1-expressing breast cancer. This could be due to the downregulation of MUC1 and IFITM1 expression upon ruxolitinib treatment that occurs both in vitro and in vivo (Figs. 3C–F and 5). We should note that the pure antiestrogen, fulvestrant (ICI 182,780), which is known to degrade ER, is currently a second-line therapy for ER-positive breast cancers that fail tamoxifen and/or AI treatment. ICI 182,780 profoundly inhibited the proliferation of AI-sensitive MCF-7 and T47D cells; however, its inhibitory effect was not augmented when combined with ruxolitinib. Similarly, ruxolitinib alone blocks IFITM1 expression and significantly inhibited the proliferation of AI-resistant MCF-7:5C cells; however, the inhibitory effect of ruxolitinib was not further enhanced when combined with fulvestrant (ICI 182,780). Although ruxolitinib blocks IFITM1 expression in the AI-sensitive T-47D cells, the JAK/STAT inhibitor did not significantly inhibit cell proliferation in this model, suggesting that AI sensitivity and ER expression is a main determinant of ruxolitinib sensitivity. These data suggest that the clinical benefit of ruxolitinib may be limited to high IFITM1-expressing AI-resistant breast tumors.

Finally, we demonstrated that MUC1 knockdown in combination with estrogen treatment causes the most effective increase in cell death in AI-resistant cells. We also examined the effects of a MUC1 inhibitor, GO-201 (45, 49), on IFITM1 expression and observed decreases in IFITM1 expression in vitro (Supplementary Fig. S1D). Currently, a MUC1 third-generation inhibitor, GO-203, is undergoing phase I/II clinical trials for patients with acute myeloid leukemia (45, 49). However, in vivo experiments within our laboratory have shown limited response to GO-201 and GO-201 demonstrates a cytotoxic effect on normal human fibroblasts (50). In addition to small-molecule inhibitors, MUC1 vaccines are also being investigated for their efficacy in multiple tumor types (51, 52). Our research elucidated two novel mechanisms to target IFITM1 expression both in vitro and in vivo through treatment with estrogen or ruxolitinib. The combination of these therapies will be explored in future research.

In conclusion, this study demonstrated that MUC1 is a key regulator of JAK/STAT signaling and IFITM1 transcription, which promotes AI-resistant breast cancer cell aggression and survival. MUC1 and IFITM1 coexpression could be used as a biomarker to isolate patients who may successfully be treated with estrogen, Rux, or combination therapy. Future studies are needed to determine the exact role that IFITM1 has in tumorigenesis and how combination therapies can effectively treat AI-resistant patients. Our data suggests that MUC1 and IFITM1 crosstalk may be a marker of aggressive breast cancer and that current therapies can be adapted to target resistant breast cancer.

No potential conflicts of interest were disclosed.

Conception and design: T.E. Escher, A.J. Lui, E.S. Geanes, O. Tawfik, C.R. Hagan, J. Lewis-Wambi

Development of methodology: T.E. Escher, A.J. Lui, E.S. Geanes, O. Tawfik, C.R. Hagan, J. Lewis-Wambi

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): T.E. Escher, A.J. Lui, E.S. Geanes, K.R. Walter, O. Tawfik, J. Lewis-Wambi

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): T.E. Escher, A.J. Lui, E.S. Geanes, C.R. Hagan, J. Lewis-Wambi

Writing, review, and/or revision of the manuscript: T.E. Escher, A.J. Lui, E.S. Geanes, K.R. Walter, O. Tawfik, C.R. Hagan, J. Lewis-Wambi

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): T.E. Escher, E.S. Geanes

Study supervision: J. Lewis-Wambi

The authors would like to thank the KUMC Confocal Imaging Facility and Imaging Core. The confocal core is supported, in part, by NIH/NIGMS COBRE grant P20GM104936. We acknowledge support from the University of Kansas (KU) Cancer Center's Biospecimen Repository Core Facility staff for helping obtain human specimens and performing histological work. The authors also acknowledge support from the KU Cancer Center's Support Grant (P30 CA168524). This study was supported, in part, by grants from the Department of Defense (W81XWH-12-1-0139, to J. Lewis-Wambi), the National Cancer Institute (K01CA120051, to J. Lewis-Wambi), start-up funds from the University of Kansas Medical Center (KUMC; to J. Lewis-Wambi), the KUMC Biomedical Research Training Program (BRTP, to A.J. Lui), the National Cancer Institute (1F30CA203160-01, to A.J. Lui), the American Medical Association (AMA) Foundation (to A.J. Lui), the University of Kansas Cancer Center (CCSG grants P30 CA168524-0, to J. Lewis-Wembi, C.R. Hagan, and K.R. Walter), and the University of New Mexico Comprehensive Cancer Center (GMaP grant 3P30CA118100-12S2, to J. Lewis-Wembi). We acknowledge the Flow Cytometry Core Laboratory at the University of Kansas Medical Center, which is sponsored, in part, by the NIH/NIGMS COBRE grant P30 GM10332.

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.

1.
Masood
S
. 
Estrogen and progesterone receptors in cytology: A comprehensive review
.
Diagn Cytopathol
1992
;
8
:
475
91
.
2.
Miller
WR
. 
Biological rationale for endocrine therapy in breast cancer
.
Best Pract Res Clin Endocrinol Metab
2004
;
18
:
1
32
.
3.
Blackwell
KL
. 
Are all aromatase inhibitors alike?
Breast Cancer Res Treat
2008
;
112
:35–43
.
4.
Chen
S
,
Masri
S
,
Wang
X
,
Phung
S
,
Yuan
YC
,
Wu
X
. 
What do we know about the mechanisms of aromatase inhibitor resistance?
J Steroid Biochem Mol Biol
2006
;
102
:
232
40
.
5.
Khodarev
NN
,
Pitroda
SP
,
Beckett
MA
,
MacDermed
DM
,
Huang
L
,
Kufe
DW
, et al
MUC1-induced transcriptional programs associated with tumorigenesis predict outcome in breast and lung cancer
.
Cancer Res
2009
;
69
:
2833
7
.
6.
Pitroda
SP
,
Khodarev
NN
,
Beckett
MA
,
Kufe
DW
,
Weichselbaum
RR
. 
MUC1-induced alterations in a lipid metabolic gene network predict response of human breast cancers to tamoxifen treatment
.
Proc Natl Acad Sci U S A
2009
;
106
:
5837
41
.
7.
MacLean
GD
,
Reddish
MA
,
Longenecker
BM
. 
Prognostic significance of preimmunotherapy serum CA27.29 (MUC-1) mucin level after active specific immunotherapy of metastatic adenocarcinoma patients
.
J Immunother
1997
;
20
:
70
8
.
8.
Kufe
DW
. 
Mucins in cancer: function, prognosis and therapy
.
Nat Rev Cancer
2009
;
9
:
874
85
.
9.
Geng
Y
,
Yeh
K
,
Takatani
T
,
King
MR
. 
Three to tango: MUC1 as a ligand for both e-selectin and ICAM-1 in the breast cancer metastatic cascade
.
Front Oncol
2012
;
2
:
76
.
10.
Kufe
D
,
Inghirami
G
,
Abe
M
,
Hayes
D
,
Justi-Wheeler
H
,
Schlom
J
. 
Differential reactivity of a novel monoclonal antibody (DF3) with human malignant versus benign breast tumors
.
Hybridoma
1984
;
3
:
223
32
.
11.
Ahmad
R
,
Alam
M
,
Rajabi
H
,
Kufe
D
. 
The MUC1-C oncoprotein binds to the BH3 domain of the pro-apoptotic BAX protein and blocks BAX function
.
J Biol Chem
2012
;
287
:
20866
75
.
12.
Wei
X
,
Xu
H
,
Kufe
D
. 
MUC1 oncoprotein stabilizes and activates estrogen receptor alpha
.
Mol Cell
2006
;
21
:
295
305
.
13.
Klinge
CM
,
Radde
BN
,
Imbert-Fernandez
Y
,
Teng
Y
,
Ivanova
MM
,
Abner
SM
, et al
Targeting the intracellular MUC1 C-terminal domain inhibits proliferation and estrogen receptor transcriptional activity in lung adenocarcinoma cells
.
Mol Cancer Ther
2011
;
10
:
2062
71
.
14.
Ma
XJ
,
Wang
Z
,
Ryan
PD
,
Isakoff
SJ
,
Barmettler
A
,
Fuller
A
, et al
A two-gene expression ratio predicts clinical outcome in breast cancer patients treated with tamoxifen
.
Cancer Cell
2004
;
5
:
607
16
.
15.
Clark
S
,
McGuckin
MA
,
Hurst
T
,
Ward
BG
. 
Effect of interferon-gamma and TNF-alpha on MUC1 mucin expression in ovarian carcinoma cell lines
.
Dis Markers
1994
;
12
:
43
50
.
16.
Choi
HJ
,
Lui
A
,
Ogony
J
,
Jan
R
,
Sims
PJ
,
Lewis-Wambi
J
. 
Targeting interferon response genes sensitizes aromatase inhibitor resistant breast cancer cells to estrogen-induced cell death
.
Breast Cancer Res
2015
;
17
:
6
.
17.
Lui
AJ
,
Geanes
ES
,
Ogony
J
,
Behbod
F
,
Marquess
J
,
Valdez
K
, et al
IFITM1 suppression blocks proliferation and invasion of aromatase inhibitor-resistant breast cancer in vivo by JAK/STAT-mediated induction of p21
.
Cancer Lett
2017
;
399
:
29
43
.
18.
Khodarev
N
,
Ahmad
R
,
Rajabi
H
,
Pitroda
S
,
Kufe
T
,
McClary
C
, et al
Cooperativity of the MUC1 oncoprotein and STAT1 pathway in poor prognosis human breast cancer
.
Oncogene
2010
;
29
:
920
9
.
19.
Jiang
SY
,
Wolf
DM
,
Yingling
JM
,
Chang
C
,
Jordan
VC
. 
An estrogen receptor positive MCF-7 clone that is resistant to antiestrogens and estradiol
.
Mol Cell Endocrinol
1992
;
90
:
77
86
.
20.
Murphy
CS
,
Pink
JJ
,
Jordan
VC
. 
Characterization of a receptor-negative, hormone-nonresponsive clone derived from a T47D human breast cancer cell line kept under estrogen-free conditions
.
Cancer Res
1990
;
50
:
7285
92
.
21.
Keydar
I
,
Chen
L
,
Karby
S
,
Weiss
FR
,
Delarea
J
,
Radu
M
, et al
Establishment and characterization of a cell line of human breast carcinoma origin
.
Eur J Cancer
1979
;
15
:
659
70
.
22.
Ogony
J
,
Choi
HJ
,
Lui
A
,
Cristofanilli
M
,
Lewis-Wambi
J
. 
Interferon-induced transmembrane protein 1 (IFITM1) overexpression enhances the aggressive phenotype of SUM149 inflammatory breast cancer cells in a signal transducer and activator of transcription 2 (STAT2)-dependent manner
.
Breast Cancer Res
2016
;
18
:
25
.
23.
Chang
TC
,
Wentzel
EA
,
Kent
OA
,
Ramachandran
K
,
Mullendore
M
,
Lee
KH
, et al
Transactivation of miR-34a by p53 broadly influences gene expression and promotes apoptosis
.
Mol Cell
2007
;
26
:
745
52
.
24.
Hagan
CR
,
Regan
TM
,
Dressing
GE
,
Lange
CA
. 
CK2-dependent phosphorylation of progesterone receptors (PR) on Ser81 regulates PR-B isoform-specific target gene expression in breast cancer cells
.
Mol Cell Biol
2011
;
31
:
2439
52
.
25.
Lewis
JS
,
Meeke
K
,
Osipo
C
,
Ross
EA
,
Kidawi
N
,
Li
T
, et al
Intrinsic mechanism of estradiol-induced apoptosis in breast cancer cells resistant to estrogen deprivation
.
J Natl Cancer Inst
2005
;
97
:
1746
59
.
26.
O'Regan
RM
,
Cisneros
A
,
England
GM
,
MacGregor
JI
,
Muenzner
HD
,
Assikis
VJ
, et al
Effects of the antiestrogens tamoxifen, toremifene, and ICI 182,780 on endometrial cancer growth
.
J Natl Cancer Inst
1998
;
90
:
1552
8
.
27.
Lewis
JS
,
Osipo
C
,
Meeke
K
,
Jordan
VC
. 
Estrogen-induced apoptosis in a breast cancer model resistant to long-term estrogen withdrawal
.
J Steroid Biochem Mol Biol
2005
;
94
:
131
41
.
28.
Lewis-Wambi
JS
,
Cunliffe
HE
,
Kim
HR
,
Willis
AL
,
Jordan
VC
. 
Overexpression of CEACAM6 promotes migration and invasion of oestrogen-deprived breast cancer cells
.
Eur J Cancer
2008
;
44
:
1770
9
.
29.
Chumsri
S
,
Howes
T
,
Bao
T
,
Sabnis
G
,
Brodie
A
. 
Aromatase, aromatase inhibitors, and breast cancer
.
J Steroid Biochem Mol Biol
2011
;
125
:
13
22
.
30.
Zaretsky
JZ
,
Barnea
I
,
Aylon
Y
,
Gorivodsky
M
,
Wreschner
DH
,
Keydar
I
. 
MUC1 gene overexpressed in breast cancer: structure and transcriptional activity of the MUC1 promoter and role of estrogen receptor alpha (ERalpha) in regulation of the MUC1 gene expression
.
Mol Cancer
2006
;
5
:
57
.
31.
Siroy
A
,
Abdul-Karim
FW
,
Miedler
J
,
Fong
N
,
Fu
P
,
Gilmore
H
, et al
MUC1 is expressed at high frequency in early-stage basal-like triple-negative breast cancer
.
Hum Pathol
2013
;
44
:
2159
66
.
32.
Maeda
T
,
Hiraki
M
,
Jin
C
,
Rajabi
H
,
Tagde
A
,
Alam
M
, et al
MUC1-C induces PD-L1 and immune evasion in triple-negative breast cancer
.
Cancer Res
2018
;
78
:
205
15
.
33.
Rajabi
H
,
Tagde
A
,
Alam
M
,
Bouillez
A
,
Pitroda
S
,
Suzuki
Y
, et al
DNA methylation by DNMT1 and DNMT3b methyltransferases is driven by the MUC1-C oncoprotein in human carcinoma cells
.
Oncogene
2016
;
35
:
6439
45
.
34.
Tagde
A
,
Rajabi
H
,
Stroopinsky
D
,
Gali
R
,
Alam
M
,
Bouillez
A
, et al
MUC1-C induces DNA methyltransferase 1 and represses tumor suppressor genes in acute myeloid leukemia
.
Oncotarget
2016
;
7
:
38974
87
.
35.
Wei
X
,
Xu
H
,
Kufe
D
. 
Human mucin 1 oncoprotein represses transcription of the p53 tumor suppressor gene
.
Cancer Res
2007
;
67
:
1853
8
.
36.
Cascio
S
,
Faylo
JL
,
Sciurba
JC
,
Xue
J
,
Ranganathan
S
,
Lohmueller
JJ
, et al
Abnormally glycosylated MUC1 establishes a positive feedback circuit of inflammatory cytokines, mediated by NF-kappaB p65 and EzH2, in colitis-associated cancer
.
Oncotarget
2017
;
8
:
105284
98
.
37.
Cui
K
,
Tailor
P
,
Liu
H
,
Chen
X
,
Ozato
K
,
Zhao
K
. 
The chromatin-remodeling BAF complex mediates cellular antiviral activities by promoter priming
.
Mol Cell Biol
2004
;
24
:
4476
86
.
38.
Abbas
T
,
Dutta
A
. 
P21 in cancer: intricate networks and multiple activities
.
Nat Rev Cancer
2009
;
9
:
400
14
.
39.
Ellis
MJ
,
Gao
F
,
Dehdashti
F
,
Jeffe
DB
,
Marcom
PK
,
Carey
LA
, et al
Lower-dose vs high-dose oral estradiol therapy of hormone receptor-positive, aromatase inhibitor-resistant advanced breast cancer: a phase 2 randomized study
.
JAMA
2009
;
302
:
774
80
.
40.
Iwase
H
,
Yamamoto
Y
,
Yamamoto-Ibusuki
M
,
Murakami
KI
,
Okumura
Y
,
Tomita
S
, et al
Ethinylestradiol is beneficial for postmenopausal patients with heavily pre-treated metastatic breast cancer after prior aromatase inhibitor treatment: a prospective study
.
Br J Cancer
2013
;
109
:
1537
42
.
41.
Klinge
CM
. 
Estrogen receptor interaction with co-activators and co-repressors
.
Steroids
2000
;
65
:
227
51
.
42.
Mitchell
S
,
Abel
P
,
Madaan
S
,
Jeffs
J
,
Chaudhary
K
,
Stamp
G
, et al
Androgen-dependent regulation of human MUC1 mucin expression
.
Neoplasia
2002
;
4
:
9
18
.
43.
Dharmaraj
N
,
Wang
P
,
Carson
DD
. 
Cytokine and progesterone receptor interplay in the regulation of MUC1 gene expression
.
Mol Endocrinol
2010
;
24
:
2253
66
.
44.
Smirnova
OV
. 
Competitive agonists and antagonists of steroid nuclear receptors: evolution of the concept or its reversal
.
Biochemistry Biokhimiia
2015
;
80
:
1227
34
.
45.
Kufe
DW
. 
MUC1-C oncoprotein as a target in breast cancer: activation of signaling pathways and therapeutic approaches
.
Oncogene
2013
;
32
:
1073
81
.
46.
ClinicalTrials.gov.
Study of ruxolitinib (INCB018424) with preoperative chemotherapy for triple negative inflammatory breast cancer
.
Available from:
https://ClinicalTrials.gov/show/NCT02876302.
47.
ClinicalTrials.gov.
Pembrolizumab and ruxolitinib phosphate in treating patients with metastatic stage IV triple negative breast cancer
.
Available from:
https://ClinicalTrials.gov/show/NCT03012230.
48.
Stover
DG
,
Gil Del Alcazar
CR
,
Brock
J
,
Guo
H
,
Overmoyer
B
,
Balko
J
, et al
Phase II study of ruxolitinib, a selective JAK1/2 inhibitor, in patients with metastatic triple-negative breast cancer
.
NPJ Breast Cancer
2018
;
4
:
10
.
49.
Joshi
S
,
Kumar
S
,
Choudhury
A
,
Ponnusamy
MP
,
Batra
SK
. 
Altered mucins (MUC) trafficking in benign and malignant conditions
.
Oncotarget
2014
;
5
:
7272
84
.
50.
Tyuryaeva
II
,
Lyublinskaya
OG
,
Podkorytov
IS
,
Skrynnikov
NR
. 
Origin of anti-tumor activity of the cysteine-containing GO peptides and further optimization of their cytotoxic properties
.
Sci Rep
2017
;
7
:
40217
.
51.
Lakshminarayanan
V
,
Supekar
NT
,
Wei
J
,
McCurry
DB
,
Dueck
AC
,
Kosiorek
HE
, et al
MUC1 vaccines, comprised of glycosylated or non-glycosylated peptides or tumor-derived MUC1, can circumvent immunoediting to control tumor growth in MUC1 transgenic mice
.
PLoS One
2016
;
11
:
e0145920
.
52.
Vasir
B
,
Zarwan
C
,
Ahmad
R
,
Crawford
KD
,
Rajabi
H
,
Matsuoka
K
, et al
Induction of antitumor immunity ex vivo using dendritic cells transduced with fowl pox vector expressing MUC1, CEA, and a triad of costimulatory molecules (rF-PANVAC)
.
J Immunother
2012
;
35
:
555
69
.