Research article | Open | Published:
Induction of cell cycle arrest and inflammatory genes by combined treatment with epigenetic, differentiating, and chemotherapeutic agents in triple-negative breast cancer
Breast Cancer Researchvolume 20, Article number: 145 (2018)
A combination of entinostat, all-trans retinoic acid, and doxorubicin (EAD) induces cell death and differentiation and causes significant regression of xenografts of triple-negative breast cancer (TNBC).
We investigated the mechanisms underlying the antitumor effects of each component of the EAD combination therapy by high-throughput gene expression profiling of drug-treated cells.
Microarray analysis showed that entinostat and doxorubicin (ED) altered expression of genes related to growth arrest, inflammation, and differentiation. ED downregulated MYC, E2F, and G2M cell cycle genes. Accordingly, entinostat sensitized the cells to doxorubicin-induced growth arrest at G2. ED induced interferon genes, which correlated with breast tumors containing a higher proportion of tumor-infiltrating lymphocytes. ED also increased the expression of immune checkpoint agonists and cancer testis antigens. Analysis of TNBC xenografts showed that EAD enhanced the inflammation score in nude mice. Among the genes differentially regulated between the EAD and ED groups, an all-trans retinoic acid (ATRA)-regulated gene, DHRS3, was induced in EAD-treated xenografts. DHRS3 was expressed at lower levels in human TNBC metastases compared to normal breast or primary tumors. High expression of ED-induced growth arrest and inflammatory genes was associated with better prognosis in TNBC patients.
Entinostat potentiated doxorubicin-mediated cell death and the combination induced inflammatory signatures. The ED-induced immunomodulation may improve immunotherapy. Addition of ATRA to ED may potentiate inflammation and contribute to TNBC regression.
All-trans retinoic acid (ATRA) induces breast cancer differentiation [1, 2], cell death , and inflammation . However, its limited treatment success in solid tumors  may be attributed to the frequent epigenetic silencing of the retinoic acid receptor beta (RAR-β) [6, 7]. We and others have shown that RAR-β is underexpressed and epigenetically silenced in breast cancer, and that histone deacetylase inhibitors (HDACi) cause reexpression of RAR-β and sensitize the cells to treatment [8, 9].
HDACi are currently employed in the clinic for the treatment of a wide variety of solid and hematological malignancies , including breast cancer . The proposed anticancer activities of HDACi include: the induction of apoptosis and cell cycle arrest [12, 13], the inhibition of angiogenesis [14, 15], and the stimulation of cancer cell differentiation . Importantly, HDACi have been shown to enhance the immunogenicity of cancer cells .
Although HDACi showed limited effect as single agents to treat breast cancer, their use in combination with other anticancer agents is currently being evaluated . Studies in advanced solid tumors in which HDACi were combined either with doxorubicin  or retinoic acid  suggested enhanced antitumor activity. Previously, we showed that a combination of entinostat, ATRA, and low-dose doxorubicin (EAD) effectively induced cell death and differentiation and decreased tumor size in xenografts of TNBC cell lines . In this study, we provide insights into additional mechanisms underlying the combined drug effect on the decrease of tumor volume. Entinostat and doxorubicin combination (ED) induced cell cycle growth arrest, interferon response, and other inflammatory signatures. Addition of ATRA to the entinostat and doxorubicin combination (EAD) most effectively induced inflammatory features in nude mice.
Patient samples, cell lines, constructs, and reagents
RNA and cDNA were generated from normal and breast cancer tissue, as described previously . MDA-MB-231 and HCC1937 cells were obtained from the American Type Culture Collection and SUM-149 and SUM-159 cells were obtained from Dr. S. Ethier (Wayne State University, Detroit, MI, USA). ATRA and doxorubicin were purchased from Sigma Chemicals. Entinostat was provided by Syndax Pharmaceuticals, LLC.
Transcriptome array and bioinformatics analysis
MDA-MB-231 cells were treated for 48 h with entinostat (2.5 μM), ATRA (1 μM), and doxorubicin (0.2 μM) singly or in combination. RNA was extracted using RNeasy Mini Kit (Qiagen) and Illumina HumanHT12v4 gene expression array was performed by the Microarray Core at Johns Hopkins. Illumina microarray data were preprocessed by background subtraction followed by quantile normalization using GenomeStudio. After preprocessing using GenomeStudio, data were imported and analyzed using R (using base and Bioconductor packages) and Ingenuity® Pathway Analysis (IPA). Detailed methods are provided in Additional file 1.
Fluorescence-activated cell quantification
For cell cycle determinations, cells were permeabilized with cold 70% ethanol and stained with propidium iodide (Sigma). Samples were run on the BD FACSCalibur system (Becton Dickinson).
Xenografts of MDA-MB-231 cells were established in six to eight athymic nude mice by injecting 2 × 106 tumor cells per flank, two flanks per mouse. The mice were treated for 4 weeks, receiving entinostat (2.5 mg/kg) for 5 days/weeks oral, ATRA (5 mg/kg) 5 days/weeks i.p., and doxorubicin (2 mg/kg) once a week i.v.
Correlation of gene expression and patient outcome
Kaplan–Meier analysis was performed by utilizing an extended version of the previously established breast cancer transcriptomic database . The Kaplan–Meier survival plot, the hazard ratio with 95% confidence intervals, and the logrank p value were calculated and plotted in R using the package “survival” .
The cell line results were expressed as mean ± standard errors of mean (SEM). Two-tailed Student’s t tests (95% confidence interval) were performed on pairwise combinations of data to determine statistical significance defined as p < 0.05, p < 0.01, and p < 0.001. Results of qRT-PCR using tumor xenograft cDNAs were analyzed using the median and two-tailed Mann–Whitney test. Statistical analyses were performed using GraphPad Prism version 5.0 (GraphPad Software, Inc.).
Entinostat and doxorubicin combination induces gene reprogramming
Previously, we showed that a combination of entinostat, ATRA, and doxorubicin (EAD) resulted in increased cell death and differentiation, and as a consequence, a decrease of tumor volume . To explore additional mechanisms of tumor regression caused by EAD, and to identify drug targets other than RAR-β, we performed high-throughput gene expression profiling analysis of MDA-MB-231 cells treated with entinostat, ATRA, and doxorubicin as single, double, and triple combinations.
An unsupervised cluster analysis using the top 10% most differentially expressed genes across all eight different treatment conditions revealed clustering of vehicle with ATRA and ED with EAD treatment groups (Fig. 1a). Principal component analysis (PCA) revealed similar grouping of samples on the first component and a minor batch effect on the second component (Additional file 2: Figure S1A) which was corrected in further analysis.
Differential expression analysis identified gene expression changes (≥ 2-fold change and false discovery rate ≤ 0.05) across all seven different combinations compared to vehicle treatment (Fig. 1b). EAD and ED combination treatments elicited the greatest number of differentially expressed genes, at 2098 and 2152 respectively. Consistent with the behavior of HDAC inhibitors, most of the genes altered by entinostat treatment alone are upregulated (228/237, or 96.2%).
Gene set analysis reveals independent and synergistic biological roles for each treatment modality
A comprehensive analysis of the pathways regulated by the treatment groups in comparison to vehicle was performed using gene set variant analysis (GSVA)  on Hallmark gene sets defined by the Molecular Signatures Database (MSigDB). ANOVA with multiple test correction was performed for each gene set across all treatment groups to identify differentially regulated gene sets (Fig. 1c).
ATRA treatment alone had a gene set profile similar to vehicle while doxorubicin downregulated genes related to mitotic spindle and UV response. As described previously for other HDACi , entinostat increased expression of genes associated with angiogenesis, as well as genes related to TNF-α, inflammatory response, apoptosis, early estrogen receptor (ER) response, and others (Fig. 1c). Doxorubicin and entinostat single treatments downregulated MYC and E2F targets and progression through the G2M checkpoint, in comparison to vehicle (Fig. 1c). However, the ED/EAD combination most effectively decreased expression of these genes related to proliferation and cell-cycle progression. ED/EAD also downregulated genes related to TGF-β, protein secretion, and unfolded protein response. Interestingly, ED/EAD combinations most effectively upregulated IFN-α, IFN-γ, complement, TP53, and UV-response genes (Fig. 1c). Although some of these latter pathways were also upregulated to a lesser extent by doxorubicin or entinostat, IFN-α and IFN-γ showed the strongest synergistic enrichment following treatment with ED (Fig. 1c).
Combined entinostat and doxorubicin regulates growth arrest, inflammation, and differentiation pathways
Due to the large gene expression effect of ED treatment, we focused our analysis on identifying differences between ED versus single treatments with entinostat or doxorubicin. In this analysis, 868 genes were differentially expressed (Additional file 2: Figure S1B, Additional file 3: Table S1, and Additional file 4: Table S2), including upregulation of cell cycle, interferon, and CTA genes and downregulation of cell cycle genes (Fig. 1d), which suggests a synergistic effect of ED combination treatment on inducing cell death and inflammation (Fig. 1c).
To determine the activation of specific biological functions potentially related to a decrease of tumor volume due to EAD treatment, we performed pathway enrichment analysis using Ingenuity® Pathway Analysis (IPA) on the 868 synergistically differentially expressed genes following ED treatment (Fig. 1e). Cancer-related pathways showed negative scores, while inflammatory response (Additional file 1: Figure S1C), cellular movement (Additional file 2: Figure S1D), and immune cell trafficking pathways had positive scores (Fig. 1e and Additional file 5: Table S3). Other pathways related to cell death and survival, growth and proliferation, and development were also identified, with mixed enrichment scores (Fig. 1e and Additional file 5: Table S3).
Consistent with IPA analysis, gene set enrichment analysis (GSEA) identified increased inflammation-related (interferon alpha and gamma), cell death (upregulation of TP53, apoptosis), and growth arrest and proliferation (decreased G2M progression, E2F and myc targets) gene sets in the ED gene signature (Additional file 6: Table S4).
ED and EAD therapies induce cell growth arrest
Among the genes related to cell growth arrest and death identified by array analysis, we validated the following ED upregulated genes by qRT-PCR: BTG2, GADD45G, CEBPA, CCNA1, and CDKN1C. Their mRNA expression levels following ED treatment were 2–4.5-fold higher than entinostat-treated MDA-MB-231 cells (Fig. 2a). In contrast to the decrease of cyclin D1 by entinostat, observed mainly in ER-positive cell lines , treatment of TNBC cells with entinostat induced cyclin D1 and its target, CDK6, mRNA in comparison to vehicle-treated cells (Additional file 7: Figure S2A). However, ED and EAD treatment significantly decreased cyclin D1 and CDK6 mRNA, in comparison to entinostat (Additional file 7: Figure S2A). Similarly, cyclin-D1 protein levels were higher in entinostat-treated cells, and significantly decreased by doxorubicin, present in ED and EAD (Fig. 2b and Additional file 7: Figure S2B, left). On the other hand, the level of the tumor suppressor cyclin-A (CCNA1) protein was decreased by entinostat treatment, and restored by doxorubicin in ED and EAD (Fig. 2b and Additional file 7: Figure S2B, right panel).
Previously, we reported that ED therapy induces cell death by an increase in apoptosis in cell culture and in tumor xenografts . EAD most effectively induced apoptosis and necrosis and cell death . Since we observed an increase in cell growth arrest genes by ED and EAD treatment (Fig. 2a, b), we further analyzed whether cell growth arrest is involved in the ED and EAD-induced cell death. As expected, high doses of doxorubicin (200 nM) induced arrest at the G2 phase of the cell cycle (Fig. 2c). Interestingly, ED increased arrest at the G2 phase in comparison to doxorubicin alone (6.25–50 nM) in MDA-MB-231 cells (Fig. 2c), SUM-149 cells (Fig. 2d), and HCC1937 cells (Additional file 8: Table S5), suggesting that entinostat increases doxorubicin sensitivity and growth arrest. Entinostat also significantly decreased the S phase of the cell cycle in the three cell lines (Fig. 2c, d and Additional file 8: Table S5), consistent with its effect on decreasing cell proliferation .
ED and EAD induce interferon genes associated with higher immune cell infiltration
Since we identified regulation of interferon (IFN) response genes as one of the most significant pathways associated with ED treatment (Fig. 1c and Additional file 6: Table S4), we analyzed the expression of type-1 IFN genes in the array, and validated expression by qRT-PCR. Hierarchical clustering analysis of these gene sets revealed upregulation of IFN-α (Fig. 3A and Additional file 9: Table S6) and IFN-γ (Additional file 10: Figure S3A and Additional file 11: Table S7) genes by ED and EAD in comparison to doxorubicin and the other groups in MDA-MB-231 cells. Within the IFN genes, GBP1, CXCL10, IRF1, and STAT1 (Fig. 3A) were described previously to be part of a tumor relapse-free signature in breast cancer patients . By qRT-PCR, we showed that CXCL10, IRF1, and STAT1 were induced in MDA-MB-231 by EAD (at low doses of doxorubicin), in comparison to ED and single treatments (Fig. 3Ba). Interferon-responsive genes, the tripartite motif (TRIM) proteins TRIM48 and TRIM51 (Fig. 3Bb), and interferon gamma (IFNG) (Additional file 10: Figure S3Ba) were induced at higher levels in MDA-MB-231 cells by ED and EAD treatments in comparison to single treatments. We observed a significant increase of CXCL10 and a trend to increase TRIM48 expression in mice xenografts of MDA-MB-231 cells treated with ED and EAD (Fig. 3Bc). The induction of proinflammatory genes such as CXCL10 and TRIM48 by epigenetic therapy was also observed in another TNBC cell line, SUM-159 (Additional file 10: Figure S3Bb).
Lastly, we sought to identify the relationship of immune infiltration in breast cancer with respect to the expression of IFN-α (Fig. 3A) and IFN-γ (Additional file 10: Figure S3A) genes induced by ED. Immune infiltration groups were identified in The Cancer Genome Atlas (TCGA) breast cancer dataset using a core set of immune genes  into clusters of low, medium, and high infiltration markers (Additional file 10: Figure S3C). Hierarchical clustering was performed in the same TCGA samples using these ED-induced IFN genes. There was a significant enrichment of samples with high immune infiltration in a cluster with high expression of ED-induced IFN genes, with stepwise decrease of immune infiltration groups with lower ED-induced IFN gene expression (Fig. 3Ca and Additional file 10: Figure S3D, Additional file 9: Table S6, and Additional file 11: Table S7). There is a strong correlation between ED-induced IFN-α (Fig. 3Cb, Cc) and IFN-γ (Additional file 10: Figure S3E, F) gene scores and immune infiltration scores (Spearman correlation = 0.72 and 0.85, respectively). A 10,000-fold gene-wise permutation analysis to identify significant association between immune infiltration genes and IFN genes revealed statistically significant association (p < 0.001) (Fig. 3Cb, Cc and Additional file 10: Figure S3E, F). These analyses suggest that the increase of IFN genes by ED and EAD treatment may increase the infiltration of immune cells into the breast tumor site.
ED induces proinflammatory genes
In addition, we identified by array analysis that ED regulated different inflammatory genes, including the CCL2, CXCL16 (Fig. 4a), and TNFSF9 (Fig. 4c) cytokines, the EGR2 and DLX3 transcription factors (Fig. 4a), the PADI4 enzyme (Fig. 4a), the CCR3 receptor (Fig 4a), and the SPANXN1 and SPANXN5 cancer testis antigens (CTAs) (Fig. 4b), which are related to immune cells recruitment and activation (Additional file 2: Figure S1C, D). qRT-PCR validation showed that entinostat induced the expression of these inflammatory genes in MDA-MB-231 cells and this effect is further potentiated by ED (Fig. 4). Since we observed the ED induction of the CTAs SPANXN5 in MDA-MB-231 cells (Fig. 4b) and SPANXN1 in both MDA-MB-231 (Fig. 4b) and SUM-159 cells (Additional file 12: Figure S4A), and the immune checkpoint agonist TNFSF9/CD137L in MDA-MB-231 cells (Fig. 4c) in the array and qRT-PCR, we studied the expression of other members of these gene families across our treatment groups.
The CTAs CTCFL/BORIS, IL13RA2, and PRAME are significantly induced by entinostat and further regulated by the combination of entinostat and doxorubicin (ED) (Fig. 4b). ED also induced the expression of different checkpoint agonists, including the ligands ICOSL and GITRL (Fig. 4c) in MDA-MB-231 cells. The immune checkpoint CD80 is also induced by entinostat, and further regulated by doxorubicin. The ligand PDL1 is induced by entinostat and significantly decreased by doxorubicin (Fig. 4c). The checkpoint agonist GITRL and the cancer testis antigen IL13RA2 proteins were also induced in MDA-MB-231 cells by ED and EAD treatment (Fig. 4d and Additional file 12: Figure S4B). These patterns suggest activation of the immune system by ED as possible contributors to tumor cell death in an in-vivo system.
EAD combination induces inflammatory features in xenografts in nude mice
Although immune deficient, athymic nude mice retain increased natural killer (NK) cell activity and tumoricidal macrophages . Therefore, we investigated whether the increase in inflammatory genes in the tumor cells by ED and EAD treatment would also be followed by an increase in inflammation. The study pathologist’s quantification of inflammation in the tumor xenografts following the eight different treatments revealed that EAD treatment, followed by ED, resulted in significantly higher inflammatory scores (Fig. 5A and Additional file 13: Table S8). The inflammatory score consisted of primarily neutrophils and macrophages, some edema, and lesser numbers of lymphocytes.
These results in tumor xenografts showed that although EAD and ED equally regulated many proinflammatory genes (Fig. 4), EAD treatment in vivo showed a higher effect on the induction of inflammatory features in nude mice.
EAD-specific inflammatory genes and regulation of DHRS3
As shown in Fig. 1a (and Additional file 2: Figure S1A) there is a high degree of overlap of the ED and EAD signatures and very few genes are unique to EAD’s triple drug effect in MDA-MB-231 cells. In order to identify genes that potentially contribute to the EAD’s effect on cell death and differentiation  and inflammation, we investigated differentially expressed genes between EAD and ED (Fig. 5B).
EAD treatment induced ELF3 and IL-1β in comparison to ED (Fig. 5B). Previously, we have described ELF3 as an EAD-induced gene with a role in differentiation ( and Fig. 5C). In addition, EAD upregulates DHRS3, CCL26, TNFα, CD14, IL-1α, and others (Fig. 5B and Additional file 14: Figure S5A). qRT-PCR validation showed that mRNA levels for the cytokine CCL26 were induced in MDA-MB-231 following EAD treatment in comparison to ED and other treatments (Fig. 5Ca). Interestingly, the dehydrogenase/reductase member 3 (DHRS3, also known as retSDR1) which is involved in maintaining the cellular supply of retinol metabolites  was identified in the array as one the few genes induced by ATRA in comparison to vehicle. DHRS3 levels were higher in EAD-treated MDA-MB-231 (Fig. 5Cb) and SUM-159 (Additional file 14: Figure S5B) cells in comparison to ED and were also induced by all combinations of ATRA treatment, including the individual treatment. Similarly, ATRA and to a lesser extent entinostat increased DHRS3 protein levels in MDA-MB-231 cells (Fig. 5Cc and Additional file 14: Figure S5C). Treatment of MDA-MB-231 cells with EAD, in the presence of a low dose of doxorubicin (12.5 nM), was more effective than single agents and ED in induction of DHRS3 expression (Fig. 5Cc and Additional file 14: Figure S5C). Treatment of cells with a high dose of doxorubicin (200 nM) also increased DHRS3 expression and the combination with ATRA (AD) most effectively caused increased expression of this protein (Fig. 5Cc and Additional file 14: Figure S5C). DHRS3 mRNA was also significantly higher in EAD-treated tumors in comparison to ED or vehicle (Fig. 5D). DHRS3 levels also tended to increase in the EA-treated tumors in comparison to vehicle (Fig. 5D). We also investigated DHRS3 expression in normal breast, primary, and metastatic TNBC to correlate their levels with tumor progression. Metastatic TNBCs showed lower levels of DHRS3 in comparison to primary tumors and normal breast (Fig. 5E). The role of ELF3, CCL26, and DHRS3 in the induction of inflammation by EAD remains to be determined.
ED-regulated genes are associated with better disease outcome in TNBC patients
We next investigated whether there is a correlation of the epigenetic therapy-induced genes with survival in patients with TNBC. Higher expression of ED-induced genes such as BTG2, CEBPA, CCNA1, and GADD45G related to growth arrest, and IFNG, TNFSF9, IL1B, and TRIM48 related to inflammatory response, significantly correlated with a better prognosis of patients with basal/TNBC. Longer overall survival (Fig. 6a, n = 241 patients), relapse-free survival (Additional file 15: Figure S6A, n = 614 patients), and metastases-free survival (Additional file 15: Figure S6B, n = 232 patients) were observed in TNBC patients with higher expression of ED-induced genes.
Previously, we have demonstrated that a combination of entinostat, ATRA, and doxorubicin (EAD) resulted in the greatest induction of cell death and cancer stem cell differentiation, and in consequence inhibition of xenografts of TNBC cells . In this article, we show additional mechanisms likely to be involved in EAD-mediated decrease of tumor volume. Global gene expression analysis of TNBC cells following EAD and single treatments revealed that ED is a potent gene reprogramming therapy, and differentially regulated the expression of many genes more efficiently than single agents (Fig. 1).
We speculate that this ED effect on gene expression could be attributed to the inhibition of topoisomerase II-beta by both agents, as we have described previously in the case of RAR-β induction . Pathway analysis showed that ED treatment regulated genes related to inflammation, growth arrest, and differentiation (Fig. 1). ED downregulated MYC and E2F targets and genes related to progression through the G2M checkpoint (Fig. 1). Recently, Topper et al.  showed that the combination of HDAC and DNA methyltransferase inhibitors decreased MYC-driven cell proliferation in lung cancer. In line with this finding related to the regulation of cell cycle arrest and death genes, we observed that entinostat treatment sensitized TNBC cells to doxorubicin-induced G2 cell cycle arrest (Fig. 2). Retinoic acid , doxorubicin , and HDACi  as single agents were also shown to induce cell arrest. Previously, we showed that ED and EAD increased apoptosis and EAD is the most effective to induce necrosis . Collectively, these data showed that the increase in G2 cell arrest by ED and EAD (Fig. 2) led to cancer cell apoptosis and necrosis  and contributed to the decrease of tumor volume .
Among the inflammatory genes, ED increased the expression of interferon (IFN) genes, which correlated to the higher levels of infiltrated immune cells in patient breast tumors (Fig. 3). These findings suggest that the induction of IFN genes by the epigenetic treatment may favor immune cell recruitment to the tumor site. Type I IFN possesses the potent ability to activate several immune cell types [32, 33]. The presence of lymphocytic infiltration in early-stage breast cancer was associated with good prognosis and high response rates to neoadjuvant chemotherapy , especially in ER−/HER2− tumors . Interestingly, CXCL10, IRF1, and STAT1, shown to be expressed in breast cancers of patients who did not relapse , were also induced by EAD (Fig. 3). In addition to the regulation of interferon responsive genes, such as members of the tripartite motif (TRIM) family (Fig. 3), ED treatment also significantly induced the expression of members of the cancer/testis antigens (CTA) family in TNBC cells (Fig. 4). Epigenetic modifications including promoter hypomethylation and histone deacetylation have important roles in CTA gene activation . CTAs are protein antigens normally expressed in a wide variety of malignant tumors but not in normal adult tissues, except for testis, and therefore have been viewed as attractive targets for cancer immunotherapy . Similarly, 5-azacitidine regulates interferon signaling, antigen processing and presentation, cytokines/chemokines, and CTA genes [30, 38, 39]. In addition, we showed that several inflammatory genes related to immune cell activation and migration, including immune checkpoints, are also epigenetically regulated and induced by entinostat and further potentiated by the ED combination (Fig. 4). In fact, HDACi have been shown to enhance immunogenicity of cancer cells. Several groups have reported the upregulation of natural killer (NK) cell-activating ligands, MHC class I and II molecules, components of the machinery for antigen presentation, and costimulatory molecules on the surface of cancer cells exposed to HDACi [40,41,42]. Entinostat potentiates the effect of immune checkpoint-blocking antibodies, and the combination decreased regulatory T cells (Tregs) and myeloid derived suppressor cells (MDSC) in mice , and in patients . Treatment with epigenetic therapy prevented exhaustion of CD8+ T cells and increased their expansion after immune-checkpoint blockade . Chemotherapy was also shown to augment tumor immunity . Previously we observed that EAD induced almost two times higher levels of necrosis compared to ED in vitro, and in tumor xenografts . The increase in the number of dying cells in tumors of EAD-treated mice  may release damage-associated molecular patterns, release antigens, and initiate an immune response . We found that entinostat in combination with doxorubicin (ED) most effectively upregulated interferon response, tumor antigens, cytokines, costimulatory molecules, and other inflammatory genes (Figs. 1, 3, and 4), and therefore may improve immunotherapies. Although we observed a similar effect of ED and EAD on the regulation of inflammatory genes, in-vivo EAD treatment was the most effective to increase the recruitment of immune cells and edema at the tumor site of nude mice (Fig. 5), the least immune-deficient mouse model .
ATRA is critical in maintaining immune homeostasis [4, 48] to differentiate myeloid-derived suppressor cells (MDSCs) into dendritic cells (DCs) and to improve their immunostimulatory capacity [49, 50]. Treatment of renal cell carcinoma patients with ATRA substantially decreased the presence of MDSCs in peripheral blood . A recent study demonstrated that in lung cancer patients, p53 vaccine-generated immune responses were improved if patients received a short course of ATRA . Some of the genes induced by EAD in comparison to ED, such as ELF3 [53, 54], DHRS3, IL-1β, CCL26 , TNF-α, CD14, and IL-1α (Fig. 5 and Additional file 12: Figure S4), play a role in inflammation. DHRS3 is involved in maintaining the cellular supply of retinol metabolites  and was described to be induced by the retinoid X receptor (RXR) rexinoid ligand bexarotene in MMTV-erbB2 mice . We also observed a correlation of the epigenetic therapy-induced genes, related to growth arrest and inflammation, with survival in patients with TNBC. Collectively, these data suggested that ED and EAD treatment likely potentiates tumor immune surveillance.
Some of the mechanisms involved in the EAD-induced decrease in tumor volume (Fig. 6b) include an increase in differentiation and cell death , due to an increase in cell cycle arrest. Entinostat potentiated doxorubicin-mediated cell death and the combination induced inflammatory signatures. The ED-induced immunomodulation may improve immunotherapy. Addition of ATRA to ED may potentiate inflammation and contribute to TNBC regression.
All-trans retinoic acid
Cancer testes antigens
Dehydrogenase/reductase member 3
Gene set enrichment analysis
Histone deacetylase inhibitors
The Cancer Genome Atlas
Triple-negative breast cancer
Ginestier C, Wicinski J, Cervera N, Monville F, Finetti P, Bertucci F, Wicha MS, Birnbaum D, Charafe-Jauffret E. Retinoid signaling regulates breast cancer stem cell differentiation. Cell Cycle. 2009;8(20):3297–302.
Bhat-Nakshatri P, Goswami CP, Badve S, Sledge GW Jr, Nakshatri H. Identification of FDA-approved drugs targeting breast cancer stem cells along with biomarkers of sensitivity. Sci Rep. 2013;3:2530.
Seewaldt VL, Johnson BS, Parker MB, Collins SJ, Swisshelm K. Expression of retinoic acid receptor beta mediates retinoic acid-induced growth arrest and apoptosis in breast cancer cells. Cell Growth Differ. 1995;6(9):1077–88.
Hall JA, Grainger JR, Spencer SP, Belkaid Y. The role of retinoic acid in tolerance and immunity. Immunity. 2011;35(1):13–22.
Connolly RM, Nguyen NK, Sukumar S. Molecular pathways: current role and future directions of the retinoic acid pathway in cancer prevention and treatment. Clin Cancer Res. 2013;19(7):1651–9.
Anso E, Mullen AR, Felsher DW, Mates JM, Deberardinis RJ, Chandel NS. Metabolic changes in cancer cells upon suppression of MYC. Cancer Metab. 2013;1(1):7.
Tang XH, Gudas LJ. Retinoids, retinoic acid receptors, and cancer. Annu Rev Pathol. 2011;6:345–64.
Sirchia SM, Ferguson AT, Sironi E, Subramanyan S, Orlandi R, Sukumar S, Sacchi N. Evidence of epigenetic changes affecting the chromatin state of the retinoic acid receptor beta2 promoter in breast cancer cells. Oncogene. 2000;19(12):1556–63.
Nacht M, Ferguson AT, Zhang W, Petroziello JM, Cook BP, Gao YH, Maguire S, Riley D, Coppola G, Landes GM, et al. Combining serial analysis of gene expression and array technologies to identify genes differentially expressed in breast cancer. Cancer Res. 1999;59(21):5464–70.
Li Y, Seto E. HDACs and HDAC inhibitors in cancer development and therapy. Cold Spring Harbor Perspect Med. 2016;6(10). https://doi.org/10.1101/cshperspect.a026831.
Connolly RM, Rudek MA, Piekarz R. Entinostat: a promising treatment option for patients with advanced breast cancer. Future Oncol. 2017;13(13): 1137–148. https://doi.org/10.2217/fon-2016-0526.
Vigushin DM, Ali S, Pace PE, Mirsaidi N, Ito K, Adcock I, Coombes RC. Trichostatin A is a histone deacetylase inhibitor with potent antitumor activity against breast cancer in vivo. Clin Cancer Res. 2001;7(4):971–6.
Huang X, Gao L, Wang S, Lee CK, Ordentlich P, Liu B. HDAC inhibitor SNDX-275 induces apoptosis in erbB2-overexpressing breast cancer cells via down-regulation of erbB3 expression. Cancer Res. 2009;69(21):8403–11.
Ray A, Alalem M, Ray BK. Loss of epigenetic Kruppel-like factor 4 histone deacetylase (KLF-4-HDAC)-mediated transcriptional suppression is crucial in increasing vascular endothelial growth factor (VEGF) expression in breast cancer. J Biol Chem. 2013;288(38):27232–42.
Kaluza D, Kroll J, Gesierich S, Yao TP, Boon RA, Hergenreider E, Tjwa M, Rossig L, Seto E, Augustin HG, et al. Class IIb HDAC6 regulates endothelial cell migration and angiogenesis by deacetylation of cortactin. EMBO J. 2011;30(20):4142–56.
Salvador MA, Wicinski J, Cabaud O, Toiron Y, Finetti P, Josselin E, Lelievre H, Kraus-Berthier L, Depil S, Bertucci F, et al. The histone deacetylase inhibitor abexinostat induces cancer stem cells differentiation in breast cancer with low Xist expression. Clin Cancer Res. 2013;19(23):6520–31.
West AC, Smyth MJ, Johnstone RW. The anticancer effects of HDAC inhibitors require the immune system. Oncoimmunology. 2014;3(1):e27414.
Munster PN, Marchion D, Thomas S, Egorin M, Minton S, Springett G, Lee JH, Simon G, Chiappori A, Sullivan D, et al. Phase I trial of vorinostat and doxorubicin in solid tumours: histone deacetylase 2 expression as a predictive marker. Br J Cancer. 2009;101(7):1044–50.
Pili R, Salumbides B, Zhao M, Altiok S, Qian D, Zwiebel J, Carducci MA, Rudek MA. Phase I study of the histone deacetylase inhibitor entinostat in combination with 13-cis retinoic acid in patients with solid tumours. Br J Cancer. 2012;106(1):77–84.
Merino VF, Nguyen N, Jin K, Sadik H, Cho S, Korangath P, Han L, Foster YM, Zhou XC, Zhang Z, et al. Combined treatment with epigenetic, differentiating, and chemotherapeutic agents cooperatively targets tumor-initiating cells in triple-negative breast cancer. Cancer Res. 2016;76(7):2013–24.
Gyorffy B, Lanczky A, Eklund AC, Denkert C, Budczies J, Li Q, Szallasi Z. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010;123(3):725–31.
Gyorffy B, Surowiak P, Budczies J, Lanczky A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer. PLoS One. 2013;8(12):e82241.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Jin G, Bausch D, Knightly T, Liu Z, Li Y, Liu B, Lu J, Chong W, Velmahos GC, Alam HB. Histone deacetylase inhibitors enhance endothelial cell sprouting angiogenesis in vitro. Surgery. 2011;150(3):429–35.
Alao JP, Lam EW, Ali S, Buluwela L, Bordogna W, Lockey P, Varshochi R, Stavropoulou AV, Coombes RC, Vigushin DM. Histone deacetylase inhibitor trichostatin A represses estrogen receptor alpha-dependent transcription and promotes proteasomal degradation of cyclin D1 in human breast carcinoma cell lines. Clin Cancer Res. 2004;10(23):8094–104.
Ascierto ML, Kmieciak M, Idowu MO, Manjili R, Zhao Y, Grimes M, Dumur C, Wang E, Ramakrishnan V, Wang XY, et al. A signature of immune function genes associated with recurrence-free survival in breast cancer patients. Breast Cancer Res Treat. 2012;131(3):871–80.
Yao J, Caballero OL, Yung WK, Weinstein JN, Riggins GJ, Strausberg RL, Zhao Q. Tumor subtype-specific cancer-testis antigens as potential biomarkers and immunotherapeutic targets for cancers. Cancer Immunol Res. 2014;2(4):371–9.
Zhou Q, Facciponte J, Jin M, Shen Q, Lin Q. Humanized NOD-SCID IL2rg−/− mice as a preclinical model for cancer research and its potential use for individualized cancer therapies. Cancer Lett. 2014;344(1):13–9.
Haeseleer F, Huang J, Lebioda L, Saari JC, Palczewski K. Molecular characterization of a novel short-chain dehydrogenase/reductase that reduces all-trans-retinal. J Biol Chem. 1998;273(34):21790–9.
Topper MJ, Vaz M, Chiappinelli KB, DeStefano Shields CE, Niknafs N, Yen RC, Wenzel A, Hicks J, Ballew M, Stone M, et al. Epigenetic therapy ties MYC depletion to reversing immune evasion and treating lung cancer. Cell. 2017;171(6):1284–300 e1221.
Bar-On O, Shapira M, Hershko DD. Differential effects of doxorubicin treatment on cell cycle arrest and Skp2 expression in breast cancer cells. Anti-Cancer Drugs. 2007;18(10):1113–21.
Zitvogel L, Galluzzi L, Kepp O, Smyth MJ, Kroemer G. Type I interferons in anticancer immunity. Nat Rev Immunol. 2015;15(7):405–14.
Minn AJ, Wherry EJ. Combination cancer therapies with immune checkpoint blockade: convergence on interferon signaling. Cell. 2016;165(2):272–5.
Denkert C, Loibl S, Noske A, Roller M, Muller BM, Komor M, Budczies J, Darb-Esfahani S, Kronenwett R, Hanusch C, et al. Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol. 2010;28(1):105–13.
Loi S, Sirtaine N, Piette F, Salgado R, Viale G, Van Eenoo F, Rouas G, Francis P, Crown JP, Hitre E, et al. Prognostic and predictive value of tumor-infiltrating lymphocytes in a phase III randomized adjuvant breast cancer trial in node-positive breast cancer comparing the addition of docetaxel to doxorubicin with doxorubicin-based chemotherapy: BIG 02-98. J Clin Oncol. 2013;31(7):860–7.
De Smet C, Lurquin C, Lethe B, Martelange V, Boon T. DNA methylation is the primary silencing mechanism for a set of germ line- and tumor-specific genes with a CpG-rich promoter. Mol Cell Biol. 1999;19(11):7327–35.
Whitehurst AW. Cause and consequence of cancer/testis antigen activation in cancer. Annu Rev Pharmacol Toxicol. 2014;54:251–72.
Roulois D, Loo Yau H, Singhania R, Wang Y, Danesh A, Shen SY, Han H, Liang G, Jones PA, Pugh TJ, et al. DNA-demethylating agents target colorectal cancer cells by inducing viral mimicry by endogenous transcripts. Cell. 2015;162(5):961–73.
Li H, Chiappinelli KB, Guzzetta AA, Easwaran H, Yen RW, Vatapalli R, Topper MJ, Luo J, Connolly RM, Azad NS, et al. Immune regulation by low doses of the DNA methyltransferase inhibitor 5-azacitidine in common human epithelial cancers. Oncotarget. 2014;5(3):587–98.
Khan AN, Magner WJ, Tomasi TB. An epigenetic vaccine model active in the prevention and treatment of melanoma. J Transl Med. 2007;5:64.
Magner WJ, Kazim AL, Stewart C, Romano MA, Catalano G, Grande C, Keiser N, Santaniello F, Tomasi TB. Activation of MHC class I, II, and CD40 gene expression by histone deacetylase inhibitors. J Immunol. 2000;165(12):7017–24.
Skov S, Pedersen MT, Andresen L, Straten PT, Woetmann A, Odum N. Cancer cells become susceptible to natural killer cell killing after exposure to histone deacetylase inhibitors due to glycogen synthase kinase-3-dependent expression of MHC class I-related chain A and B. Cancer Res. 2005;65(23):11136–45.
Kim K, Skora AD, Li Z, Liu Q, Tam AJ, Blosser RL, Diaz LA Jr, Papadopoulos N, Kinzler KW, Vogelstein B, et al. Eradication of metastatic mouse cancers resistant to immune checkpoint blockade by suppression of myeloid-derived cells. Proc Natl Acad Sci U S A. 2014;111(32):11774–9.
Tomita Y, Lee MJ, Lee S, Tomita S, Chumsri S, Cruickshank S, Ordentlich P, Trepel JB. The interplay of epigenetic therapy and immunity in locally recurrent or metastatic estrogen receptor-positive breast cancer: correlative analysis of ENCORE 301, a randomized, placebo-controlled phase II trial of exemestane with or without entinostat. Oncoimmunology. 2016;5(11):e1219008.
Ghoneim HE, Fan Y, Moustaki A, Abdelsamed HA, Dash P, Dogra P, Carter R, Awad W, Neale G, Thomas PG, et al. De novo epigenetic programs inhibit PD-1 blockade-mediated T cell rejuvenation. Cell. 2017;170(1):142–157 e119.
Emens LA. Chemotherapy and tumor immunity: an unexpected collaboration. Front Biosci. 2008;13:249–57.
Kono K, Mimura K, Kiessling R. Immunogenic tumor cell death induced by chemoradiotherapy: molecular mechanisms and a clinical translation. Cell Death Dis. 2013;4:e688.
Raverdeau M, Mills KH. Modulation of T cell and innate immune responses by retinoic Acid. J Immunol. 2014;192(7):2953–8.
Kusmartsev S, Cheng F, Yu B, Nefedova Y, Sotomayor E, Lush R, Gabrilovich D. All-trans-retinoic acid eliminates immature myeloid cells from tumor-bearing mice and improves the effect of vaccination. Cancer Res. 2003;63(15):4441–9.
Kusmartsev S, Su Z, Heiser A, Dannull J, Eruslanov E, Kubler H, Yancey D, Dahm P, Vieweg J. Reversal of myeloid cell-mediated immunosuppression in patients with metastatic renal cell carcinoma. Clin Cancer Res. 2008;14(24):8270–8.
Mirza N, Fishman M, Fricke I, Dunn M, Neuger AM, Frost TJ, Lush RM, Antonia S, Gabrilovich DI. All-trans-retinoic acid improves differentiation of myeloid cells and immune response in cancer patients. Cancer Res. 2006;66(18):9299–307.
Iclozan C, Antonia S, Chiappori A, Chen DT, Gabrilovich D. Therapeutic regulation of myeloid-derived suppressor cells and immune response to cancer vaccine in patients with extensive stage small cell lung cancer. Cancer Immunol Immunother. 2013;62(5):909–18.
Oliver JR, Kushwah R, Hu J. Multiple roles of the epithelium-specific ETS transcription factor, ESE-1, in development and disease. Lab Invest. 2012;92(3):320–30.
Kushwah R, Oliver JR, Wu J, Chang Z, Hu J. Elf3 regulates allergic airway inflammation by controlling dendritic cell-driven T cell differentiation. J Immunol. 2011;187(9):4639–53.
Kitaura M, Suzuki N, Imai T, Takagi S, Suzuki R, Nakajima T, Hirai K, Nomiyama H, Yoshie O. Molecular cloning of a novel human CC chemokine (Eotaxin-3) that is a functional ligand of CC chemokine receptor 3. J Biol Chem. 1999;274(39):27975–80.
Li Y, Zhang Y, Hill J, Kim HT, Shen Q, Bissonnette RP, Lamph WW, Brown PH. The rexinoid, bexarotene, prevents the development of premalignant lesions in MMTV-erbB2 mice. Br J Cancer. 2008;98(8):1380–8.
Karn T, Pusztai L, Ruckhaberle E, Liedtke C, Muller V, Schmidt M, Metzler D, Wang J, Coombes KR, Gatje R, et al. Melanoma antigen family A identified by the bimodality index defines a subset of triple negative breast cancers as candidates for immune response augmentation. Eur J Cancer. 2012;48(1):12–23.
Cancer Genome Atlas Network, Comprehensive molecular portraits of human breast tumours. Nature. 2012;490(7418):61–70.
The authors thank Dr Q. Zhao for the IPA data processing and interpretation, and Dr David Huso for the analysis of tumor histology.
This work was supported by funding from the DOD BCRP Center of Excellence (Grant W81XWH-04-1-0595 to SS and Grant W81XWH-09-1-0499 to VFM). BG was supported by the NVKP_16-1-2016-0037 grant of the National Research, Development and Innovation Office, Hungary.
Availability of data and materials
The GEO accession number for the data reported in this article is GSE63351.
Ethics approval and consent to participate
Freshly resected breast tissue of women undergoing reduction mammoplasty, primary tumors, and pleural effusion from women undergoing treatment and from the Rapid Autopsy Program, were provided by the Johns Hopkins Surgical Pathology Department, under protocols approved by the institutional review board. All animal studies were performed according to the guidelines and approval of the Animal Care Committee of the Johns Hopkins School of Medicine, Baltimore.
Consent for publication
All authors approved the final version of this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional methods. (DOCX 45 kb)
Figure S1. Gene expression array profile of TNBC cells identified entinostat and doxorubicin (ED) as a gene reprogramming combination. (A) Principal component analysis (PCA) 3D projection of gene expression data obtained from microarray analysis onto first three principal components. Each ball represents a different sample; different treatments indicated. 1 and 2 are duplicates from different batches. (B) Venn diagram showing number of genes common (intersection) and unique to each indicated treatment combination. Each treatment gene signature derived following normalization with genes present in MDA-MB-231 cells after vehicle (DMSO) treatment. Ingenuity® Pathway Analysis (IPA) generated network of (C) inflammatory and (D) cellular movement signaling upregulated in MDA-MB-231 cells treated with ED combination compared to single treatments. Color indicates genes upregulated (red) and downregulated (green) in MDA-MB-231 cells by ED treatment. Open and closed edges indicate genes with direct and indirect relationships respectively. E entinostat, A all-trans retinoic acid, D doxorubicin (different combinations). (PPTX 579 kb)
Table S1. ED genes in comparison to entinostat and doxorubicin treatments. (DOCX 41 kb)
Table S2. Genes differentially expressed by ED treatment and validated by qRT-PCR. (DOCX 14 kb)
Table S3. Ingenuity® Pathway Analysis of ED genes (DOCX 13 kb)
Table S4. Gene set analysis on ED genes. (DOCX 12 kb)
Figure S2. ED and EAD induce cell growth arrest. (A) qRT-PCR of cyclin D1 (CCND1) and cyclin-dependent kinase 6 (CDK6) (identified by array analysis) in MDA-MB-231 cells treated with entinostat (2.5 μM), ATRA (1 μM), and doxorubicin (0.2 μM) singly, and combinations, for 48 h. (B) ImageJ quantification of cyclin D1 (left) and cyclin A (right) protein expression in MDA-MB-231 cells treated with entinostat (2.5 μM), ATRA (1 μM), and doxorubicin (0.2 μM) singly, and combinations, for 48 h. *Compared to entinostat in qRT-PCR: *p < 0.05, **p < 0.01. t test used to compare mean level of mRNA expression (± SEM), after RPL39 normalization. (PPTX 75 kb)
Table S5. ED and EAD induce growth arrest in HCC1937 TNBC cells. (DOCX 15 kb)
Table S6. IFN-α genes induced by ED in MDA-MB-231 cells and correlated with immune infiltration. (DOCX 14 kb)
Figure S3. ED induces interferon gamma genes associated with an increase in tumor lymphocytes. (A) Hierarchical supervised clustering of expression of interferon-gamma (IFN-G) genes against signatures of MDA-MB-231 cells following treatments. (B) qRT-PCR of (a) IFN-G in MDA-MB-231 and (b) CXCL10 and TRIM48 in SUM-159 cells treated with EAD singly and in combinations (doxorubicin 12.5 and 200 nM). (C) Unsupervised hierarchical cluster analysis of tumor-infiltrating lymphocyte genes , used in Fig. 3C to classify immune infiltration (low, medium, and high) in TCGA RNA-seq breast cancer patient dataset . (D) Hierarchical supervised clustering of expression of IFN-γ genes against TCGA RNA-seq breast cancer patient dataset. Bars above identify different tumor subtypes (PAM50) and inflammatory cell content (immune, low–high) identified in (C). (E) One-way ANOVA showed significant difference across one or more groups (#1 low, #2 medium, #3 high immune cells) and post-hoc pairwise Student t test revealed statistically significant differences across all groups (p < 0.05). (F) IFN-γ score correlation with immune infiltration. *p < 0.05, **p < 0.01, ***p < 0.001. (PPTX 538 kb)
Table S7. IFN-γ genes induced by ED in MDA-MB-231 cells and correlated with immune infiltration. (DOCX 15 kb)
Figure S4. ED regulates expression of inflammatory genes. (A) qRT-PCR of SPANXN1 in SUM-159 cells treated as described in text (doxorubicin 200 nM). t test used to compare mean level of mRNA expression (± SEM) after RPL39 normalization. **p < 0.01, ***p < 0.001. (B) ImageJ quantification of GITRL, IL13RA2, and housekeeping β-actin protein in MDA-MB-231 cells treated as described (doxorubicin 12.5 and 200 nM). (PPTX 74 kb)
Table S8. Inflammation score significance in mouse xenografts. (DOCX 14 kb)
Figure S5. genesEAD regulates inflammatory genes. qRT-PCR of TNF-α (a) and CD14 and IL1a (b) in MDA-MB-231 cells (A) and DHRS3 in SUM-159 cells (B) treated as described in text (doxorubicin 12.5 and 200 nM). t test used to compare mean level of mRNA expression (± SEM) after RPL39 normalization. *p < 0.05, **p < 0.01. (C) ImageJ quantification of DHRS3 and housekeeping β-actin proteins in MDA-MB-231 cells treated as described. (PPTX 105 kb)
Figure S6. ED-induced genes correlate with a better prognosis in TNBC patients. Kaplan–Meier curves of relapsefree survival (RFS) (A) and metastases-free survival (DMFS) (B) showing correlation of ED-induced gene expression and prognosis in basal/TNBC patients, over a period of 12 years. (PPTX 378 kb)