Abstract
Interferon regulatory factor-1 (IRF1) is a tumor suppressor that regulates cell fate in several cell types. Here, we report an inverse correlation in expression of nuclear IRF1 and the autophagy regulator ATG7 in human breast cancer cells that directly affects their cell fate. In mice harboring mutant Atg7, nuclear IRF1 was increased in mammary tumors, spleen, and kidney. Mechanistic investigations identified ATG7 and the cell death modulator beclin-1 (BECN1) as negative regulators of IRF1. Silencing ATG7 or BECN1 caused estrogen receptor-α to exit the nucleus at the time when IRF1 nuclear localization occurred. Conversely, silencing IRF1 promoted autophagy by increasing BECN1 and blunting IGF1 receptor and mTOR survival signaling. Loss of IRF1 promoted resistance to antiestrogens, whereas combined silencing of ATG7 and IRF1 restored sensitivity to these agents. Using a mathematical model to prompt signaling hypotheses, we developed evidence that ATG7 silencing could resensitize IRF1-attenuated cells to apoptosis through mechanisms that involve other estrogen-regulated genes. Overall, our work shows how inhibiting the autophagy proteins ATG7 and BECN1 can regulate IRF1-dependent and -independent signaling pathways in ways that engender a new therapeutic strategy to attack breast cancer. Cancer Res; 75(6); 1046–55. ©2015 AACR.
Introduction
Breast cancer remains one of the most prevalent cancers among American women (1). With more than 70% of these breast cancers expressing the estrogen receptor-α (ER+), many tumors depend on estrogen to drive cell proliferation and promote cell survival. The routine use of antiestrogens, such as tamoxifen and fulvestrant (Faslodex; ICI 182,780; ICI), and aromatase inhibitors, including letrozole and anastrozole, are now established as effective therapies for the treatment of breast cancer (2, 3). The widespread use of endocrine therapy is a major contributor leading to the decline in breast cancer–related mortality in women with ER+ disease (4). However, clinical efficacy of these therapies is often limited by aberrations in prosurvival and prodeath signaling leading to drug resistance. Precisely, how signaling within breast cancer cells regulates the cell-fate decision in response to endocrine therapies remains largely unknown.
Interferon regulatory factor-1 (IRF1) is a nuclear transcription factor that is activated by many immune effector molecules, including type I and type II interferons (IFNγ), TNFα, retinoic acid, IL1, IL6, and in response to viral infection (5). Upon activation, IRF1 alters the transcriptional activity of genes involved in immunomodulation, antiviral responses, and tumor suppression (5, 6). Previous studies have established a role for IRF1 in inhibiting breast cancer cell growth, including by opposing the actions of IRF2 and inducing apoptosis (7, 8). Furthermore, gene-expression microarrays comparing antiestrogen-responsive (MCF7/LCC1) and -resistant (MCF7/LCC9) human breast cancer cell lines report that IRF1 is downregulated in resistant cells, suggesting that loss of IRF1 may contribute to endocrine resistance through its reduced ability to induce cell death (9–11). Indeed, IRF1 mRNA levels correlate with a pathologic complete response and reduced risk of recurrence and death in some breast cancers (12, 13). Recent evidence suggests that autophagy proteins have a major effect on the regulation of inflammatory transcriptional responses. Embryonic mice lacking the autophagy-related 5 (Atg5) gene have increased inflammation in tissue with impaired clearance of apoptotic cells, raising the possibility that autophagy has a role in inflammation and autoimmunity (14).
A critical decision point in the determination of cell fate in response to endocrine therapies appears to be affected by the balance between autophagy (prosurvival) and apoptosis (prodeath; ref. 15). The activities of IRF1 in multiple cell types suggest that it may be a key component in regulating cell-fate decisions and perhaps a critical regulator of the switch between apoptosis and autophagy. Thus, we investigated the biologic function of IRF1 in autophagy in the context of endocrine responsiveness in breast cancer. We established an important (inverse) link between ATG7 and IRF1 expression/function, which we confirmed using immunohistochemical studies in human breast carcinomas and tissue from Atg7+/− mice. We also showed that inhibition of the autophagy proteins ATG7 and beclin-1 (BECN1) promoted IRF1 signaling and induced apoptotic cell death. Conversely, inhibition of IRF1 prolonged cell survival and promoted autophagy. Thus, our data suggest that IRF1 inhibits breast cancer cell growth through its ability to regulate both autophagy and apoptosis. To better understand the molecular interactions, we constructed an influence diagram of the novel signaling and then used this as a guide to build mathematical models using our experimental data. Model simulations suggested additional experimentation to distinguish whether ATG7 knockout promotes endocrine responsiveness through either IRF1-dependent or -independent signaling. Collectively, these data indicate a major role for IRF1 in regulating cell-fate decisions in breast cancer and also establish an IRF1-independent path as a further signaling component in ATG7-mediated cell death.
Materials and Methods
Cell culture, reagents, and siRNA treatments
MCF7, T47D, BT-474, and MDA-MB-231 cells were maintained in improved minimal essential media (IMEM) with phenol red and supplemented with 5% FBS (Life Technologies). MCF7/LCC1 (LCC1) and MCF7/LCC9 (LCC9) cells were grown in phenol red-free IMEM supplemented with 5% charcoal-stripped calf serum. All cells were maintained in a humidified atmosphere at 37°C and 95% air/5% CO2. ATG7 (SignalSilence; Cell Signaling Technology), BECN1 and STAT1 (three unique siRNAs for each target; OriGene), IRF1 (Silencer Select; consisting of three different siRNA for same target; Life Technologies), or control (Ctrl) siRNA were transfected in cells using Lipofectamine RNAiMAX (Life Technologies) according to the manufacturer's instructions. ERα and ATG7 cDNA were from Origene; ICI 182,780 (ICI; Faslodex; Fulvestrant) was from Tocris Bioscience; hydroxychloroquine (HCQ), 3-methyladenine (3-MA), and N-acetylcysteine (NAC) were from Sigma-Aldrich and the JAK inhibitor 1 (C19H16FN3O) was obtained from EMD Millipore.
Cell proliferation
Cells were transfected with IRF1, ATG7, BECN1, or Ctrl siRNA and seeded at a density of 0.2 × 106 per well in 24-well plates. One day after plating, cells were treated with the indicated concentration of ICI or vehicle control. Cells were incubated with ICI for 6 days and then cell density was measured as previously described (10, 16).
Western blot analysis
Lysates were harvested from transfected cell monolayers and protein analysis was measured as previously described (17). Membranes were probed for proteins of interest at 4°C overnight. To confirm equal loading of the gels, membranes were reprobed for β-actin (1:1,000; Santa Cruz Biotechnology).
MMP, autophagosome formation, and ROS assays
Cells were reverse transfected with ATG7, BECN1, IRF1, or Ctrl siRNA and plated in 6-well tissue culture plates. The following day, cells were treated with 100 nmol/L ICI or vehicle control. Cells were harvested 48 hours later and stained as described in the Mitochondrial Permeability Detection Kit for flow cytometry (Enzo). Accumulation of autophagic vesicles was measured using a modified monodansylcadaverine according to the manufacturer's instructions (Enzo Cyto-ID Autophagy detection kit). Total reactive oxygen species (ROS) were stained according to the Enzo's Total ROS Detection Kit instructions. Stained cells were detected and appropriate signals measured by FACS (LCCC FACS Shared Resource).
Autophagosome maturation and localization studies
Cells (1 × 105) were reverse transfected with IRF1, ATG7, BECN1, or Ctrl siRNA and seeded onto 18 × 18-mm glass coverslips. The following day, IRF1 siRNA cells were transfected with LC3 tagged with a GFP. Twenty-four hours later, cells were treated with 500 nmol/L ICI or vehicle, and then fixed and stained for IRF1 as previously described (11). IRF1 and ERα were also measured in LCC1 cells following knockdown of ATG7, BECN1, or IRF1 to determine subcellular localization.
Mathematical modeling and data fitting
Matlab (version 7.9.0) was used to build the mathematical model and perform simulations. An ordinary differential equation formalism (18–20) was used to model the outcomes observed in experiments. Experimental data on cell proliferation under different conditions were used to fit the model by a least-squares approach. Further details of the model are described in Supplementary Fig. S1.
In vivo experiments
Tumors from female Atg7+/− and wild-type (WT) mice (21) were treated with medroxyprogesterone acetate (DepoProvera; Pfizer) and 7,12-dimethylbenz(a)anthracene (DMBA; Sigma) to induce mammary tumors, as described earlier (22). Mice were euthanized once tumors reached 10% of the mouse body size, or when a single tumor volume reached 1,000 mm3; mammary tumors and organs were removed at necroscopy, fixed in neutral buffered formalin, and processed using routine histologic methods.
Human tissue microarrays
Human tissue microarrays (TMA) were obtained from the Familial Cancer Registry (FCR) at Georgetown Lombardi Comprehensive Cancer Center. TMAs contained 107 cores from breast tumors of FCR participants (2 cores/subject) with a known BRCA1/2 mutation status (known BRCA1/2 mutation carriers and noncarriers). The histologic diagnosis was confirmed by a pathologist at the time of the IHC scoring. Patient/tumor characteristics are shown in Supplementary Table S1.
Immunohistochemistry
Sections (5 μm) from paraffin-embedded human TMAs and mouse mammary tissue were stained with mouse anti-IRF1 (1:100; BD Biosciences) or rabbit anti-APG7 (1:130; Santa Cruz Biotechnology) antibodies as previously described (17). Assessments of immunohistochemical staining was performed by a pathologist blinded both to the clinical and molecular data. Cytoplasmic and nuclear staining of ATG7 and IRF1 were described in terms of intensity (0, negative; 1+, weak; 2+, intermediate; 3+, strong) and distribution (0, negative; +1, <10%; 2+, 11%–50%; 3+, >50%). A simplified total IHC score (the sum of the intensity and distribution scores) was calculated for each core. For in vivo staining, a computer-assisted counting technique with a grid filter to select cells was used to quantify the immunohistochemical staining of IRF1 (23). Cells stained positive were expressed as a percentage of the total cell number examined (100 cells sampled from each tissue site within each breast tumor section).
Statistical analysis
Prism 5.0 was used for statistical analysis. Differences between two groups were analyzed by Student t tests, and multiple group comparisons were assessed by one-way ANOVA followed by Bonferroni post hoc tests. The Pearson correlation coefficient was calculated to measure the correlation between ATG7 and IRF1 in human tumor tissue. Differences were considered statistically significant at P < 0.05.
Results
Nuclear IRF1 and ATG7 are inversely correlated in vivo and in human breast cancer tissue
We used IHC to study IRF1 expression in mammary tumors from WT mice, and mice with a mutated Atg7 allele (+/−); Atg7 null mice (−/−) are not viable (21). Ninety-eight percent of breast adenocarcinoma tumors generated by this model were ER+. Mammary tumors, spleen, and kidney from Atg7+/− mice exhibited increased nuclear IRF1 staining (arrows in Fig. 1A) compared with their WT controls (Fig. 1A; P < 0.001). We then examined whether nuclear IRF1 protein expression was associated with ATG7 protein expression in human disease by staining TMAs containing n = 107 cores from breast tumors. Because only nuclear IRF1 is likely to be transcriptionally active, we scored separately nuclear and cytoplasmic IRF1 staining (24). We found the primary form of IRF1 to be cytoplasmic (expressed in 80% of cases), whereas nuclear IRF1 was only detected in 13% (14/107) of tumors. Furthermore, we found that >96% of samples (90/94) with positive ATG7 expression had low or undetectable nuclear IRF1 expression (Fig. 1B and C). Conversely, >70% of samples (10/14) with nuclear IRF1 expression had low ATG7 expression. An inverse correlation was found between ATG7 and nuclear IRF1 expression in these human breast tumors (Fig. 1C; r = −0.28, P < 0.01). We then examined the expression of IRF1 and ATG7 in ER+ and ER− human breast tumors using IHC. ER status was identified in 52 of the 107 human breast cancer cores. ER+ tumors exhibited significantly greater ATG7 (P < 0.001) and decreased nuclear IRF1 (P < 0.01) expression when compared with ER− tumors (Fig. 1D and Supplementary Table S2), suggesting that the reduced nuclear IRF1 expression in ER+ breast tumors may be functionally related to increased ATG7 expression.
Nuclear IRF1 and ATG7 are inversely correlated in vivo and in human breast cancers. A, immunohistochemical staining of IRF1 in WT and ATG7+/− mouse mammary tumors (top), kidney (middle), and spleen (bottom). Positive cells were quantified using ImageJ software; n = 3. B, immunohistochemical staining of ATG7 (left) and IRF1 (right) in matched human breast cancer tissues. Arrows, nuclear IRF1 staining. C, tissue sections of 107 breast cancer cores were immunostained for the expression of ATG7 and IRF1, and their correlation was analyzed by Pearson rank correlation (P < 0.01). D, average immunohistochemical score; n = 52; **, P < 0.01; ***, P < 0.001 between indicated groups.
Nuclear IRF1 and ATG7 are inversely correlated in vivo and in human breast cancers. A, immunohistochemical staining of IRF1 in WT and ATG7+/− mouse mammary tumors (top), kidney (middle), and spleen (bottom). Positive cells were quantified using ImageJ software; n = 3. B, immunohistochemical staining of ATG7 (left) and IRF1 (right) in matched human breast cancer tissues. Arrows, nuclear IRF1 staining. C, tissue sections of 107 breast cancer cores were immunostained for the expression of ATG7 and IRF1, and their correlation was analyzed by Pearson rank correlation (P < 0.01). D, average immunohistochemical score; n = 52; **, P < 0.01; ***, P < 0.001 between indicated groups.
ATG7 and BECN1 knockdown stimulates IRF1 expression
Because a putative link was established between IRF1 and ATG7, we used various inhibitors and siRNA targeting of ATG7 and other key autophagy proteins in breast cancer cell lines and measured IRF1 protein expression. Antiestrogen-sensitive LCC1 cells were treated with vehicle (control), 100 nmol/L ICI, 5 mmol/L 3-MA, 10 μmol/L HCQ, or reverse transfected with ATG5, ATG7, or BECN1 siRNA (Fig. 2A) for 48 hours before total protein was harvested. ATG5 and ATG7 are both involved in autophagosome membrane formation and elongation (14, 15). Levels of IRF1 protein were measured in breast cancer cell lines with deficient autophagy and were found to be significantly elevated in LCC1 cells with reduced ATG7 and BECN1 (Fig. 2A; P < 0.001). When further studied in other ER+ breast cancer cell lines, including MCF7 (antiestrogen-sensitive), LCC9 (antiestrogen-resistant), T47D (antiestrogen-sensitive and p53 mutant), and BT-474 (antiestrogen-sensitive, HER2 amplified), we observed that ATG7- and BECN1-knockout cells (except ATG7 knockout in BT-474 cells) had elevated expression of IRF1 compared with control cells (Fig. 2B and C). We further determined that ATG7 and BECN1 siRNA induced apoptosis in LCC1 and LCC9 cells as measured by an increase in mitochondrial membrane permeability (Fig. 2D). However, the ER− cell line MDA-MB-231 did not exhibit enhanced IRF1 expression following transfection with either ATG7 or BECN1 siRNA (Fig. 2E). We then overexpressed ATG7 in LCC1 cells and found that ectopic expression of ATG7 significantly reduced IRF1 expression (Fig. 2F; P < 0.05). These results suggest that ER+ breast cancer cells may use autophagy machinery (ATG7 and BECN1) to regulate IRF1 expression and block apoptosis.
ATG7 and BECN1 inhibit IRF1 and apoptosis in human ER+ breast cancer cells. A, Western blot images of IRF1, ATG5, ATG7, and BECN1 in LCC1 cells treated with vehicle, 100 nmol/L ICI, 5 mmol/L 3-MA, 10 μmol/L HCQ, or transfected with ATG5, ATG7, or BECN1 siRNA for 48 hours. B and C, MCF7, LCC9, T47D, and BT-474 were transfected with Ctrl, ATG7 siRNA (B), or BECN1 siRNA (C) and the following day, treated with 100 nmol/L ICI for 48 hours. Western blot hybridization was used to measure IRF1, ATG7, and BECN1 protein expression. D, mitochondrial permeability assay measured by flow cytometry in LCC1 and LCC9 cells transfected with ATG7, BECN1, or control (Ctrl) siRNA. E, IRF1 expression was measured in MDA-MB-231 cells transfected with ATG7, BECN1, or Ctrl siRNA. F, LCC1 cells transfected with ATG7 cDNA or empty vector were treated with 100 nmol/L ICI for 48 hours and expressions of IRF1 and ATG7 were analyzed by Western blot analysis. β-Actin served as the loading control; n = 3 independent experiments; *, P < 0.05; **, P < 0.01; ***, P < 0.001 versus control/vehicle experiment.
ATG7 and BECN1 inhibit IRF1 and apoptosis in human ER+ breast cancer cells. A, Western blot images of IRF1, ATG5, ATG7, and BECN1 in LCC1 cells treated with vehicle, 100 nmol/L ICI, 5 mmol/L 3-MA, 10 μmol/L HCQ, or transfected with ATG5, ATG7, or BECN1 siRNA for 48 hours. B and C, MCF7, LCC9, T47D, and BT-474 were transfected with Ctrl, ATG7 siRNA (B), or BECN1 siRNA (C) and the following day, treated with 100 nmol/L ICI for 48 hours. Western blot hybridization was used to measure IRF1, ATG7, and BECN1 protein expression. D, mitochondrial permeability assay measured by flow cytometry in LCC1 and LCC9 cells transfected with ATG7, BECN1, or control (Ctrl) siRNA. E, IRF1 expression was measured in MDA-MB-231 cells transfected with ATG7, BECN1, or Ctrl siRNA. F, LCC1 cells transfected with ATG7 cDNA or empty vector were treated with 100 nmol/L ICI for 48 hours and expressions of IRF1 and ATG7 were analyzed by Western blot analysis. β-Actin served as the loading control; n = 3 independent experiments; *, P < 0.05; **, P < 0.01; ***, P < 0.001 versus control/vehicle experiment.
IRF1 localizes to the nucleus following ATG7 and BECN1 knockdown
To examine further the relationships among IRF1, ATG7/BECN1, and ER status, we used confocal microscopy to establish the subcellular localization of IRF1 and ERα following ATG7 and BECN1 knockdown. ER+ LCC1 cells were transfected with IRF1, ATG7, BECN1, or control (Ctrl) siRNA for 48 hours and stained for IRF1 (green) and ERα (red); DAPI (blue) staining indicates nuclear localization (Fig. 3A). In LCC1 control cells, IRF1 and ERα were predominantly localized in the nucleus and a significant amount of each colocalized (Fig. 3A). Transfection with ATG7 or BECN1 siRNA enhanced the nuclear expression of IRF1; however, ERα was no longer detected in nuclei that had IRF1 induced by siATG7 and siBECN1 (Fig. 3A). These results were further confirmed by detecting increased IRF1 in the nuclear fraction of LCC1 and LCC9 cells transfected with ATG7 or BECN1 siRNA (Supplementary Fig. S2).
ATG7 and BECN1 knockdown causes IRF1 nuclear localization. A, LCC1 cells transfected with IRF1, ATG7, BECN1, or Ctrl siRNA were used to determine IRF1 and ERα localization by immunofluorescent confocal microscopy. B and C, LCC1 and LCC9 cells were transfected with Ctrl, ATG7 (B), or BECN1 siRNA (C) and treated with 100 nmol/L ICI for 48 hours. Western blot hybridization of protein homogenates was used to measure IRF1, ATG7, BECN1, STAT1, P-STAT1, and BCL2. D, MDA-MB-231 cells were transfected with ERα or control cDNA and treated with vehicle or 100 nmol/L ICI for 48 hours. Western blot hybridization of protein homogenates was used to measure ERα and IRF1 expression. E, Western blot images of LCC1 cells transfected with Ctrl, ATG7, STAT1, and ATG7+STAT1 siRNA and treated with vehicle, 100 nmol/L ICI, or 15 nmol/L JAK inhibitor (JI1) for 48 hours. F, Western blot images of LCC1 cells transfected with Ctrl, BECN1, STAT1, and BECN1+STAT1 siRNA and treated with vehicle, 100 nmol/L ICI, or 15 nmol/L JAK inhibitor (JI1) for 48 hours. β-Actin served as the loading control; n = 3 independent experiments.
ATG7 and BECN1 knockdown causes IRF1 nuclear localization. A, LCC1 cells transfected with IRF1, ATG7, BECN1, or Ctrl siRNA were used to determine IRF1 and ERα localization by immunofluorescent confocal microscopy. B and C, LCC1 and LCC9 cells were transfected with Ctrl, ATG7 (B), or BECN1 siRNA (C) and treated with 100 nmol/L ICI for 48 hours. Western blot hybridization of protein homogenates was used to measure IRF1, ATG7, BECN1, STAT1, P-STAT1, and BCL2. D, MDA-MB-231 cells were transfected with ERα or control cDNA and treated with vehicle or 100 nmol/L ICI for 48 hours. Western blot hybridization of protein homogenates was used to measure ERα and IRF1 expression. E, Western blot images of LCC1 cells transfected with Ctrl, ATG7, STAT1, and ATG7+STAT1 siRNA and treated with vehicle, 100 nmol/L ICI, or 15 nmol/L JAK inhibitor (JI1) for 48 hours. F, Western blot images of LCC1 cells transfected with Ctrl, BECN1, STAT1, and BECN1+STAT1 siRNA and treated with vehicle, 100 nmol/L ICI, or 15 nmol/L JAK inhibitor (JI1) for 48 hours. β-Actin served as the loading control; n = 3 independent experiments.
Protein expression of antiapoptotic BCL2 is induced by 17β-estradiol, whereas BCL2 promoter activity is reduced by IRF1 (11). We measured BCL2 protein expression in LCC1 and LCC9 cells transfected with ATG7, BECN1, or Ctrl siRNA. Knockdown of BECN1 significantly inhibited BCL2 expression in LCC1 and LCC9 cells, whereas knockdown of ATG7 only reduced BCL2 in LCC1 cells (Fig. 3B and C). The association between ERα and IRF1 was further investigated by measuring basal ERα and IRF1 protein levels in MDA-MB-231 (231-wt), and MDA-MB-231 cells transfected with a full-length ERα cDNA (231-ER). 231-ER cells had one third less IRF1 expression compared with their 231-wt counterparts (P < 0.05; Fig. 3D). These data strongly support a functional link between ER-mediated IRF1 regulation and BCL2 repression in autophagy-deficient breast cancer cells.
We also investigated whether inhibition of the JAK–STAT pathway would block IRF1 expression in LCC1 cells deficient in ATG7 and BECN1. Despite inhibition of STAT1 and STAT2, IRF1 remained elevated in ATG7 (Fig. 3E) and BECN1 knockdown cells (Fig. 3F). IRF1 gene expression is also regulated by ROS formation (25). LCC1 cells transfected with ATG7 siRNA had increased levels of ROS compared with control-transfected cells (Supplementary Fig. S3). Treatment with 5 mmol/L NAC, an antioxidant that suppresses ROS, had no effect on IRF1 suppression in ATG7-deficient cells (Supplementary Fig. S3). Collectively, these data demonstrate that ATG7 and BECN1 regulate IRF1 through STAT1/2- and ROS-independent mechanisms.
Knockdown of IRF1 stimulates autophagy
LCC1 and LCC9 cells were transfected with IRF1 or Ctrl siRNA and treated with vehicle or 100 nmol/L ICI for 48 hours. Consistent with previous reports (11), we found increased levels of BCL2 when IRF1 was repressed (Fig. 4A). Western blotting also demonstrated that the total expression of LC3B-II (the lipidated form of LC3 found in autophagosome membranes) was increased in LCC1 and LCC9 cells expressing little to no detectable IRF1 (Fig. 4A). Furthermore, expression of p62 (an autophagic cargo marker degraded in the final stages of autophagy) was reduced in LCC1 cells transfected with IRF1 siRNA (Fig. 4A).
IRF1 knockdown induces autophagy. A, LCC1 and LCC9 cells were transfected with IRF1 siRNA and treated with 100 nmol/L ICI for 48 hours. Western blot analysis of protein homogenates was used to measure expression of IRF1, BCL2, LC3B, and p62. B, autophagosome formation assay in LCC1 and LCC9 cells transfected with IRF1 or Ctrl siRNA and treated with 100 nmol/L ICI. Data are presented as the percentage of total cells positive for green fluorescence; n = 3; *, P < 0.05; ***, P < 0.001 versus control/vehicle experiment. C, LCC1 cells were transfected with IRF1/Ctrl siRNA and LC3–GFP and then treated with ICI. D, LCC9 cells were transfected with IRF1 or control cDNA and treated with 100 nmol/L ICI. Representative Western blot images of indicated proteins. Equivalent protein loading was assessed by measuring β-actin; n = 3 independent experiments.
IRF1 knockdown induces autophagy. A, LCC1 and LCC9 cells were transfected with IRF1 siRNA and treated with 100 nmol/L ICI for 48 hours. Western blot analysis of protein homogenates was used to measure expression of IRF1, BCL2, LC3B, and p62. B, autophagosome formation assay in LCC1 and LCC9 cells transfected with IRF1 or Ctrl siRNA and treated with 100 nmol/L ICI. Data are presented as the percentage of total cells positive for green fluorescence; n = 3; *, P < 0.05; ***, P < 0.001 versus control/vehicle experiment. C, LCC1 cells were transfected with IRF1/Ctrl siRNA and LC3–GFP and then treated with ICI. D, LCC9 cells were transfected with IRF1 or control cDNA and treated with 100 nmol/L ICI. Representative Western blot images of indicated proteins. Equivalent protein loading was assessed by measuring β-actin; n = 3 independent experiments.
We used confocal microscopy and flow cytometry to confirm that breast cancer cells with reduced IRF1 have enhanced autophagic vacuole formation. Although ICI slightly enhanced autophagosome formation (P < 0.05), cells transfected with IRF1 siRNA had significantly higher numbers of autophagic vacuoles (P < 0.001) compared with Ctrl siRNA–treated cells (Fig. 4B). By immunofluorescence, we measured LC3 punctae formation in LCC1 cells transfected with IRF1 or Ctrl siRNA (Fig. 4C). IRF1 knockdown cells showed elevated LC3 punctae compared with control cells. In contrast, overexpression of IRF1 in LCC9 cells resulted in an accumulation of p62 (Fig. 4D). Taken together, these results suggest that IRF1 inhibits autophagy through modulation of autophagy.
To identify autophagy effectors downstream of IRF1 action, we measured the expression of selected autophagy and growth factor proteins in LCC1 and LCC9 cells transfected with IRF1 or Ctrl siRNA. BECN1 contains an interferon-stimulated response element (ISRE) in its promoter and is negatively regulated by IRF1 in LCC1 cells (Supplementary Fig. S3A–S3C). Stimulation of mTOR inhibits the initiation of autophagy and acts downstream of insulin-like growth factor 1 receptor (IGF1R; ref. 15). Knockdown of IRF1 with siRNA inhibited IGF1R and mTOR protein expression (Supplementary Fig. S4A–S4C). These findings suggest that IRF1 may inhibit autophagy through alteration of the IGF1R–mTOR and BECN1 pathways.
Mathematical modeling of IRF1-related signaling data
We next determined how silencing ATG7 and/or IRF1 affects ICI responsiveness in antiestrogen-sensitive breast cancer cells. Loss of antiestrogen sensitivity with dominant-negative IRF1 (dnIRF1) was confirmed by crystal violet assay (Supplementary Fig. S5). To establish that IRF1 siRNA had the same effect on ICI responsiveness, we measured the effect of IRF1 knockdown on ICI-mediated inhibition of cell density (Fig. 5A). Western blot analysis confirmed that IRF1 and ATG7 successfully inhibited expression of the appropriate target (Fig. 5A and B). LCC1 cells transfected with IRF1 siRNA were markedly less sensitive to growth inhibition by 10 to 1,000 nmol/L ICI treatment compared with Ctrl siRNA–transfected cells (Fig. 5A; P < 0.05). Cells were then transfected with ATG7 siRNA alone and in combination with IRF1 siRNA to determine whether ATG7 knockdown could resensitize IRF1-deficient cells to ICI. Although IRF1 siRNA blocked ICI responsiveness, LCC1 cells were growth inhibited by ATG7 siRNA alone, which was further decreased by 100 nmol/L ICI treatment (Fig. 5B). LCC1 cells with silenced IRF1 and ATG7 had restored sensitivity to 100 nmol/L ICI compared with cells transfected with IRF1 siRNA (Fig. 5B; P < 0.001).
Mathematical modeling of proliferation data and experimental validation. A, LCC1 cells transfected with IRF1 or Ctrl siRNA were treated with vehicle control or ICI for 6 days and cell density measured by crystal violet. B, LCC1 cells transfected with Ctrl, IRF1, ATG7, or IRF1+ATG7 siRNA were treated with vehicle control or 100 nmol/L ICI for 6 days before analysis by crystal violet assay. Relative cell densities of each case are plotted. Western blot analysis shows each knockdown target of interest; n = 3 independent experiments; *, P < 0.05; ***, P < 0.001 versus indicated groups. C, wiring diagram of model. D, model predictions for increasing the concentration of siATG7 by a factor of 2. E, predictions for increasing the concentration of siIRF1 by a factor of 3. F, experimental results for increasing the concentration of siIRF1 by a factor of 3; n = 3 independent experiments; ***, P < 0.001 versus indicated groups.
Mathematical modeling of proliferation data and experimental validation. A, LCC1 cells transfected with IRF1 or Ctrl siRNA were treated with vehicle control or ICI for 6 days and cell density measured by crystal violet. B, LCC1 cells transfected with Ctrl, IRF1, ATG7, or IRF1+ATG7 siRNA were treated with vehicle control or 100 nmol/L ICI for 6 days before analysis by crystal violet assay. Relative cell densities of each case are plotted. Western blot analysis shows each knockdown target of interest; n = 3 independent experiments; *, P < 0.05; ***, P < 0.001 versus indicated groups. C, wiring diagram of model. D, model predictions for increasing the concentration of siATG7 by a factor of 2. E, predictions for increasing the concentration of siIRF1 by a factor of 3. F, experimental results for increasing the concentration of siIRF1 by a factor of 3; n = 3 independent experiments; ***, P < 0.001 versus indicated groups.
To both capture and quantitatively understand the experimental data detailed above, we created a phenomenologic mathematical model of LCC1 cells based on the key molecular components (ERα, IRF1, and ATG7) and their interactions (Fig. 5C). The model was used to predict the effect of further knocking down ATG7 and IRF1, respectively. The model predictions in Fig. 5D and E show that an increase in siATG7 concentration would not significantly affect the results, whereas increasing the siIRF1 concentration should increase the proliferation in response to ICI for both the siIRF1 and siIRF1+siATG7 cases (Fig. 5E). We then experimentally validated these model predictions by repeating the study with an increased concentration of siIRF1. Supplementary Figure S6 shows that the higher concentration of siIRF1 further reduced IRF1 protein expression without effecting cell density. Figure 5F shows that siIRF1+siATG7 still results in ICI sensitivity, suggesting that ATG7 knockdown activates cell death through both IRF1-dependent and -independent mechanisms.
High IRF1 expression is associated with prolonged survival
Because IRF1 promotes antiestrogen sensitivity, we investigated the relationship between high and low IRF1 expression and breast cancer–specific survival among several patient datasets. Using the dataset by Wang and colleagues (26), we found that high IRF1 expression is associated with prolonged survival in the total population of patients with breast cancer analyzed (Fig. 6A) and in the ER+ subgroup only (Fig. 6B). We also analyzed a combination of clinical datasets and found that high IRF1 expression is associated with prolonged survival in more than 1,600 patients (Fig. 6C) and in a smaller subset of ER− patients (Fig. 6D). These findings suggest that activation of IRF1 may be important for improved patient survival.
High IRF1 expression is associated with prolonged survival in human datasets. A and B, Kaplan–Meier plots show survival analysis for IRF1 expression using the Wang and colleagues (26). dataset. Of note, 286 samples were split into quartiles according to IRF1 expression level (A); 209 samples with ER+ status were split and analyzed the same way (B). C and D, several datasets were combined using the Kaplan–Meier Plotter, and IRF1 expression was analyzed in the total population (C) and ER− subgroup (D).
High IRF1 expression is associated with prolonged survival in human datasets. A and B, Kaplan–Meier plots show survival analysis for IRF1 expression using the Wang and colleagues (26). dataset. Of note, 286 samples were split into quartiles according to IRF1 expression level (A); 209 samples with ER+ status were split and analyzed the same way (B). C and D, several datasets were combined using the Kaplan–Meier Plotter, and IRF1 expression was analyzed in the total population (C) and ER− subgroup (D).
Discussion
The results presented here reveal that IRF1 is a critical signaling protein that contributes to the switch between apoptosis and autophagy to determine breast cancer cell fate. Upregulation of autophagy has been shown to protect breast cancer cells from antiestrogens (15). Previous studies also report that IRF1 mediates the antiproliferative and proapoptotic effects in cancer cells (9, 10, 27–29). Thus, we investigated the relationship between nuclear IRF1 and the key autophagy protein, ATG7, in vivo and in human breast cancer. Atg7-deficient mice exhibited increased nuclear IRF1 staining in DMBA-induced mammary tumors, kidney, and spleen, suggesting a broad role for IRF1 in autophagy signaling. Immunohistochemical staining of IRF1 in human breast tumors was broadly consistent with a study by Zhu and colleagues (24), further supporting a tumor-suppressive role for IRF1. ATG7 expression was detected in nearly all tumors. Because autophagy has a prosurvival role in breast cancer, ATG7 expression could enhance the aggressiveness of breast cancer cells and their ability to block proapoptotic stimuli. Furthermore, ATG7 could negatively regulate the activation state of IRF1 to alter breast cancer aggressiveness and response to therapy. Interestingly, we found that ER− tumor cells expressed detectable (1+ and higher) nuclear IRF1 staining compared with ER+ breast tumor cells. These findings suggested that some breast tumors may differentially regulate the activation state/subcellular localization of IRF1 through modulation of ATG7 and ERα.
We are the first to show that silencing the autophagy proteins, ATG7 and BECN1, induced IRF1 and restored drug sensitivity to ICI. Previous studies have reported that nude mice with atg5-deficient thymi showed massive inflammatory infiltrates, suggesting that autophagy negatively regulated T-cell activation (30, 31). Although we did not observe IRF1 induction with the autophagy inhibitors 3-MA, HCQ, or ATG5 siRNA, ATG7 and BECN1 may target different immunologic signaling molecules, specifically IRF1. We also found little to no IRF1 induction in BT-474 or MDA-MB-231 cells when either ATG7 or BECN1 was silenced (Fig. 2). BT-474 cells have amplified HER2, and in 3T3 cells transduced with Her2/neu, IRF1 became constitutively active (5). Thus, BT-474 cells may have constitutively active IRF1, and knockdown of ATG7 with siRNA had no additional effect on IRF1 activation. WT MDA-MB-231 cells do not express ERα; thus, IRF1 induction by ATG7 and BECN1 may be partly due to ERα regulation.
We found that knockdown of ATG7 sensitized antiestrogen-resistant breast cancer cells to ICI through mechanisms involving ERα and IRF1. Silencing of ATG7 and BECN1 caused ERα to exit and IRF1 to enter the nucleus (Fig. 3A). It is possible that ERα was inactivated or degraded whereas IRF1 was activated in response to ATG7 and BECN1 knockdown. However, the role of altered receptor coactivator or corepressor expression in this situation is still unclear. For example, the transcription factor p300 can regulate IRF1 promoter activity and acts synergistically with ERα (32). Examining the role of other coactivators/corepressors is ongoing. IRF1 also inhibited BCL2 protein expression when either ATG7 or BECN1 was silenced (Fig. 3B and C). This observation is clinically important because tamoxifen-resistant breast tumors have sustained BCL2 expression (33, 34). We previously showed that inhibiting BCL2 with the BH3 mimetic, GX15-070 (obatoclax), in combination with an antiestrogen successfully inhibited breast cancer cell growth through modulation of apoptosis and autophagy (16).
Although IRF1 promotes apoptosis (10), we established that loss of IRF1 enhanced autophagy in human breast cancer cells. IRF1 inhibited autophagic vacuole formation and reduced BECN1 protein expression (Fig. 3B–D). BECN1 contains an interferon-regulated response element sequence and is regulated by BCL2 (35), suggesting a possible interaction between IRF1 and BECN1. These data suggested that IRF1 might play a role in enhancing cell death by preventing autophagic cargo degradation. To further delineate how IRF1 inhibited the initiation of autophagy, we examined the expression of selected proteins to find downstream targets modified by loss of IRF1. IGF1R is implicated in autophagy signaling and can also promote endocrine resistance (36). IRF1 was necessary for IGF1R protein stability (Supplementary Fig. S4B–S4D), suggesting that IRF1 can disrupt IGF1R/mTOR–mediated autophagy (Supplementary Fig. S4B–S4D). IRF1 was previously shown to inhibit the autophagic response in lipopolysaccharide-stimulated macrophages by modulating the activation of mTOR (37).
Using our experimental data, we constructed two mathematical models fitting the observation that silencing ATG7 or BECN1 can block the effect of IRF1 siRNA on ICI responsiveness. To select the most correct model, an experiment for which the models predicted different results was identified and run, illustrating the usefulness of mathematical modeling to help evaluate complex biologic behaviors and prioritize experiments. The new experimental results supported the view that ATG7 knockdown could also resensitize IRF1-deficient cells to ICI through an IRF1-independent mechanism. This could be due to a switch from autophagy to apoptosis and be mediated by other proteins that are shared between pathways, such as BCL2 (38). It is also possible that silencing ATG7 activated other cell death pathways such as necroptosis or an ER-mediated stress pathway (39, 40). Our colocalization studies also supported the modeling by showing that ERα left the nucleus and stimulated nuclear IRF1 expression following ATG7 knockdown.
Previous reports suggest that high mRNA levels of IRF1 are a marker for good prognosis for patients with breast cancer given neoadjuvant treatment (12, 41). We now show that high IRF1 expression contributes to prolonged survival in more than 1,600 cases, highlighting the potential for IRF1 activation in the clinical setting. Furthermore, restoring IRF1 expression or controlling its modulation may be useful for the treatment of both ER+ and ER− breast tumors. Because IFNγ has produced mixed results as an antitumor agent in the clinic (5), using newer compounds that activate IRF1, such as baicalein, might restore IRF1 activity and lead to better outcomes for some patients (42). We now show that inhibiting key autophagy proteins, ATG7 and BECN1 also induces IRF1, which could provide a possible therapeutic option. This novel ATG7–IRF1 axis provides new insight into the mechanisms underlying antiestrogen resistance, and implicates restoration of IRF1 as a potential therapeutic strategy for the treatment of breast cancer.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: J.L. Schwartz-Roberts, A.N. Shajahan-Haq, R.B. Riggins, R. Clarke
Development of methodology: J.L. Schwartz-Roberts, A.N. Shajahan-Haq, B.R. Haddad
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J.L. Schwartz-Roberts, K.L. Cook, M. Axelrod, A. Wärri, B.V. Kallakury
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J.L. Schwartz-Roberts, C. Chen, L. Jin, B.R. Haddad, B.V. Kallakury, W.T. Baumann, R. Clarke
Writing, review, and/or revision of the manuscript: J.L. Schwartz-Roberts, K.L. Cook, A.N. Shajahan-Haq, A. Wärri, R.B. Riggins, L. Jin, B.V. Kallakury, W.T. Baumann, R. Clarke
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J.L. Schwartz-Roberts
Study supervision: J.L. Schwartz-Roberts, A.N. Shajahan-Haq, R. Clarke
Acknowledgments
The authors acknowledge the assistance provided by the Nontherapeutic Subject Registry (Familial Cancer Registry), Microscopy and Imaging, Flow Cytometry, Histopathology and Tissue, and Tissue Culture Shared Resources at Georgetown University Medical Center. These shared resources are partially supported by the NIH/NCI Cancer Center Support Grant (CCSG) P30-CA051008. The authors thank Dr. J.H. Yim (City of Hope) for critically reviewing this article. The authors also thank Drs. M. Liu, T. Sherman, and P. Furth for their guidance and support.
Grant Support
This work was supported by Public Health Service awards U54-CA149147, R01-CA131465, and Susan G. Komen grant KG090245 (R. Clarke). J.L. Schwartz-Roberts is the recipient of NIH training grant F31CA165514-01A1 that also supported this research study.
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.