Effective therapies for non–small cell lung cancer (NSCLC) remain challenging despite an increasingly comprehensive understanding of somatically altered oncogenic pathways. It is now clear that therapeutic agents with potential to impact the tumor immune microenvironment potentiate immune-orchestrated therapeutic benefit. Herein, we evaluated the immunoregulatory properties of histone deacetylase (HDAC) and bromodomain inhibitors, two classes of drugs that modulate the epigenome, with a focus on key cell subsets that are engaged in an immune response. By evaluating human peripheral blood and NSCLC tumors, we show that the selective HDAC6 inhibitor ricolinostat promotes phenotypic changes that support enhanced T-cell activation and improved function of antigen-presenting cells. The bromodomain inhibitor JQ1 attenuated CD4+FOXP3+ T regulatory cell suppressive function and synergized with ricolinostat to facilitate immune-mediated tumor growth arrest, leading to prolonged survival of mice with lung adenocarcinomas. Collectively, our findings highlight the immunomodulatory effects of two epigenetic modifiers that, together, promote T cell–mediated antitumor immunity and demonstrate their therapeutic potential for treatment of NSCLC.

Significance: Selective inhibition of HDACs and bromodomain proteins modulates tumor-associated immune cells in a manner that favors improved T-cell function and reduced inhibitory cellular mechanisms. These effects facilitated robust antitumor responses in tumor-bearing mice, demonstrating the therapeutic potential of combining these epigenetic modulators for the treatment of NSCLC. Cancer Discov; 7(8); 852–67. ©2017 AACR.

This article is highlighted in the In This Issue feature, p. 783

The tumor microenvironment contains a variety of stromal/supporting cells, including those of hematopoietic origin. Specific cell types such as T lymphocytes, natural killer (NK) cells, and various myeloid cells play diverse roles in shaping the course of the immune response against cancer (1–3). Of particular importance are antigen-presenting cells (APC), which process tumor-associated antigens (TAA) and present them to potential tumor-reactive T cells in a process that leads to T-cell priming. Despite their potential to generate an antitumor response, T-cell antigen–specific responses are often inhibited by immunosuppressive cells such as myeloid-derived suppressor cells (MDSC) and CD4+FOXP3+ regulatory T cells (Treg; refs. 4–6). Thus, agents that promote stimulatory immunophenotypic and functional changes in APCs and effector T cells or disrupt inhibitory cellular mechanisms may generate a more permissive tumor microenvironment favoring enhanced antitumor responses.

Histone deacetylases (HDAC) are a family of enzymes that modulate the expression of genes or their products by removing acetyl groups from lysine residues on histone and nonhistone proteins. As a consequence, they regulate numerous cellular processes, some of which have been implicated in aspects of tumor development and behavior such as differentiation and cell-cycle progression in malignant cells, as well as apoptosis (7–9). Class I and II HDACs are among the most active regulators of histone and nonhistone proteins in cancer (10), and their pharmacologic inhibition is seeing increasing interest for treatment of solid cancers. The contribution of the various classes of the HDAC family to tumor pathogenesis is well documented. Class I HDACs have been reported to contribute to the transformation of epithelial cells and/or metastasis (11, 12), whereas members of class II HDACs play critical roles in regulating autophagy, proliferation, and differentiation in tumor cells (13–16). The roles of class III (Sirtuins) and class IV HDAC family members are still emerging, although existing reports implicate their involvement in DNA damage response pathways, chromatin regulation, genome integrity, cell survival under stress, tumor metabolism, and cell death (17–19). Given their clear importance in cancer pathogenesis, pharmacologic inhibition of HDACs has become an area of significant interest for cancer therapy.

Several studies have evaluated the relative utility of both pan- and isozyme-selective HDAC inhibitors for cancer therapy (20, 21). Although promising therapeutic outcomes have been demonstrated using pan-HDAC inhibitors, especially when combined with chemotherapeutic drugs, their use is associated with adverse effects and dose-limiting toxicity in patients (22–25). Thus, isozyme-specific inhibitors may offer improved benefits arising from their selective inhibition of some but not all HDACs. Although many translational studies exploring the therapeutic efficacy of HDAC inhibitors have focused on their effects on malignant cells (26–28), multiple studies have highlighted the immunomodulatory activity of HDAC inhibitors on various immune cellular subsets when used alone or in combination with other immunotherapeutic agents (8, 9, 20, 21, 29–33). Furthermore, HDAC inhibitors have been shown to promote tumor cell recognition by immune cells via induction of MHC constituents and antigen-presentation machinery on tumor cells (34). Despite these data, the effects of HDAC inhibitors on the cellular constituents of the tumor immune microenvironment remain poorly understood.

Bromodomain inhibitors are a second class of agents that alter cellular epigenetic and transcriptional programs through regulating the recognition of acetylated lysine residues by transcriptional machinery. JQ1 is a well-characterized inhibitor of the bromodomain and extraterminal (BET) family of bromodomain-containing proteins and may indirectly regulate epigenetic footprints by modulating the interactions of histone acetyl transferases (HAT) and HDACs with transcription factors and proteins involved in gene expression (35). Although JQ1, like other BET inhibitors, has demonstrated efficacy in hematologic malignancies, its efficacy in solid cancers has only recently been explored. For example, JQ1 was shown to synergize with the pan-HDAC inhibitor SAHA to suppress tumor growth in models of pancreatic ductal and lung adenocarcinoma (36). Although the mechanisms governing therapeutic synergy were not fully explored in that report, it was suggested that these agents function to repress MYC expression and inflammatory cascades that regulate IL6 levels (36). Thus, the effect of JQ1 on tumor-associated immune cells from this and other studies remains largely undefined.

Recent preclinical and clinical data generated with antibodies that target immune checkpoints, namely PD-1/PD-L1 and CTLA4, demonstrate that antitumor immunity can be driven in multiple tumor types by targeting inhibitory interactions among T cells, APCs, and tumor cell populations (37, 38). Although the use of these agents represents a major advance in the treatment of cancer, it remains the case that in most tumor types, only a minority of patients respond and the barriers to response are diverse and poorly characterized. As the durability of antitumor immunity is increasingly appreciated to hinge upon immune cell contributions, we sought to identify the immunomodulatory effects of small-molecule epigenetic modulators including selective HDAC inhibitors and JQ1. We describe significant dynamic changes with respect to phenotype and function of T cells and monocytes/macrophages in peripheral tissues as well as in the tumor microenvironment by ricolinostat (ACY-1215), an HDAC inhibitor selective for the class II HDAC6. We also describe the Treg-specific effects of JQ1, which serves as the basis for its use as a partner agent with ricolinostat, a combination that facilitated sustained antitumor responses in immunocompetent, genetically engineered mouse models of non–small cell lung cancer (NSCLC).

Ricolinostat Reduced the Frequency of CD4+FOXP3+ Tregs and Promoted Transient Activation of Peripheral and Tumor-Associated T Lymphocytes

HDAC inhibition results in tumor growth arrest in immunocompetent animal models (39–41). As HDACs are ubiquitously expressed in many tissues, including immune cells, an outstanding issue is whether some of the reported therapeutic benefits of HDAC inhibitors stem from their effects on immune cells. In this regard, we sought to determine the effects of selective inhibition of members of class I and class II HDACs on leukocytes. First, we cultured peripheral blood mononuclear cells (PBMC) from healthy donors with ricolinostat and compared it with a panel of HDAC inhibitors with varying selectivity for class I HDACs 1/2/3 and class II HDAC6, including ACY-241, ACY-738, ACY-775, and tubastatin A. This panel also included entinostat (MS-275), a class I HDAC-selective agent currently being evaluated in clinical trials, in order to assess the relative impact of selective inhibition of class I HDAC enzymes (Supplementary Table S1). Notably, viability of cells cultured for 24 hours with these HDAC inhibitors was largely similar (Supplementary Fig. S1A). Furthermore, the proportions (Supplementary Fig. S1B) and phenotype (not shown) of leukocyte subsets, including monocytes, B cells, NK cells, CD4+FOXP3 and CD8+ conventional T cells, were not significantly affected by these drugs. Although the percentage of CD4+FOXP3+ Tregs within PBMCs was reduced when cultured with ACY-738 or entinostat, this effect was most striking with ricolinostat or ACY-241 treatment (Supplementary Fig. S1C and S1D). Interestingly, the highly selective HDAC6 inhibitors ACY-775 and tubastatin A did not significantly affect Treg frequency. This finding suggests that inhibition of class I HDAC and class II HDAC6 by compounds such as ricolinostat and ACY-241 may affect Treg frequency more so than single class-selective agents such as entinostat.

Next, PBMCs from patients with NSCLC (Supplementary Table S2) or bulk single-cell suspensions from dissociated freshly resected patient lung tumors were cultured with ricolinostat or entinostat for 24 or 72 hours. Flow cytometry analysis of leukocyte subsets within the cell cultures revealed that ricolinostat significantly reduced Treg proportions relative to DMSO control or entinostat [Fig. 1A (P = 0.04) and D (P = 0.05); Supplementary Fig. S1E], similar to what was observed for healthy donor PBMCs (Supplementary Fig. S1C and S1D). In parallel, expression of the CD69 activation marker was substantially upregulated on conventional CD8+ and CD4+ T cells in both patient [Fig. 1B (P = 0.003) and C (P = 0.01)] and healthy donor (Supplementary Fig. S2A and S2B) PBMC cultures in the presence of ricolinostat but not entinostat, suggesting that an increased activation profile is promoted by ricolinostat. In addition, we observed a substantial increase in phosphorylated STAT3 levels that accompanied increased CD69 expression in T cells within healthy donor PBMCs exposed to ricolinostat (Supplementary Fig. S2C and S2D).

Figure 1.

Reduction of CD4+FOXP3+ Tregs and upregulation of CD69 on T cells in PBMCs from patients with NSCLC and dissociated tumor cultures in the presence of ricolinostat. PBMCs from patients with NSCLC (A–C, G) or dissociated tumors (D–F) were cultured for 24 or 72 hours, respectively, with ricolinostat or entinostat, after which the frequency and phenotype of T-cell subsets were assessed by FACS. DMSO was used as a control. A and D, Percentage of CD4+FOXP3+ Tregs in (A) PBMCs cultured for 24 hours or (D) tumors of patients with NSCLC cultured for 72 hours with 2.5 μm of indicated HDAC inhibitors. B, C, E, F, Summary (left), and representative histograms (right) of expression levels of CD69 on gated (B and E) CD8+ and (C and F) CD4+ T cells within the PBMC and tumor cultures. G, Patient PBMCs cultured and treated as described above were washed and restimulated with PMA + ionomycin in the presence of golgi plug for 6 hours. Percentage of IFNγ-positive CD3+CD8+ T cells as assessed by intracellular cytokine staining. Data represent the mean ±SEM of samples analyzed from 5 patients. *, P < 0.05; ***, < 0.001.

Figure 1.

Reduction of CD4+FOXP3+ Tregs and upregulation of CD69 on T cells in PBMCs from patients with NSCLC and dissociated tumor cultures in the presence of ricolinostat. PBMCs from patients with NSCLC (A–C, G) or dissociated tumors (D–F) were cultured for 24 or 72 hours, respectively, with ricolinostat or entinostat, after which the frequency and phenotype of T-cell subsets were assessed by FACS. DMSO was used as a control. A and D, Percentage of CD4+FOXP3+ Tregs in (A) PBMCs cultured for 24 hours or (D) tumors of patients with NSCLC cultured for 72 hours with 2.5 μm of indicated HDAC inhibitors. B, C, E, F, Summary (left), and representative histograms (right) of expression levels of CD69 on gated (B and E) CD8+ and (C and F) CD4+ T cells within the PBMC and tumor cultures. G, Patient PBMCs cultured and treated as described above were washed and restimulated with PMA + ionomycin in the presence of golgi plug for 6 hours. Percentage of IFNγ-positive CD3+CD8+ T cells as assessed by intracellular cytokine staining. Data represent the mean ±SEM of samples analyzed from 5 patients. *, P < 0.05; ***, < 0.001.

Close modal

In contrast to patient PBMCs, we found that tumor-infiltrating immune cells in the two-dimensional (2-D) cultures of NSCLC patient tumors were significantly less viable in the presence of entinostat when compared with ricolinostat (Supplementary Fig. S3), precluding further analysis of entinostat-treated tumor cultures. Although many tumor-infiltrating T cells already exhibited an activated profile based on CD69 positivity ex vivo (data not shown), exposure to ricolinostat led to a modest increase in the frequency of CD69+CD8+ T cells (Fig. 1E, P = 0.39) and a significant increase in CD69+CD4+ cells (Fig. 1F, P = 0.04). This observation suggests that ricolinostat enhanced activation of T cells in a setting where tumor-relevant antigens were likely present. Furthermore, after removal of drugs from primary cultures, CD8+ T cells from patient PBMCs that were re-stimulated in secondary cultures with phorbol myristate acetate (PMA) and ionomycin showed a higher frequency of IFNγ positivity (Fig. 1G, P = 0.05; Supplementary Fig. S4A–S4C), indicating enhanced effector cytokine secretion and cytotoxic capability following exposure to ricolinostat.

Ricolinostat Promotes Increased Expression of MHC and Costimulatory Molecules on Human Monocytes and Tumor-Associated Macrophages

Existing reports demonstrate that HDAC6 inhibition can effect functional changes in APCs through a mechanism that involves regulation of inflammatory cytokine production (42). We therefore explored whether ricolinostat alters APC phenotype and function. Specifically, we evaluated CD14+CD11b+ monocytes in the peripheral blood of patients with NSCLC and CD14CD45+CD68+CD11b+ macrophages within disaggregated NSCLC tumors. Upon 24-hour culture with ricolinostat, the CD14+ monocyte fraction significantly upregulated surface expression of MHC class II molecules (P = 0.04), as well as CD86 (P = 0.01) but not CD80 (data not shown), phenotypic changes that are associated with increased priming ability of APCs (43). In contrast, entinostat only promoted upregulation of CD86 while not affecting MHC class II expression (Fig. 2A and B). Ricolinostat elicited a similar pattern of increased expression of MHC class II and CD86 on the CD14CD45+CD68+CD11b+ macrophages within 2-D tumor cultures (Fig. 2C and D) and in PBMCs obtained from healthy donors (Supplementary Fig. S5A and S5B). These findings demonstrate a unique modulatory effect of ricolinostat on human monocytic cells and tumor-infiltrating macrophages, and suggest that ricolinostat promotes phenotypic changes that support enhanced antigen presentation and costimulatory capabilities. Consistent with this notion, ricolinostat-exposed CD14+ monocytes were superior at inducing allogeneic T-cell proliferative responses in mixed lymphocyte reactions (Fig. 2E).

Figure 2.

Increased expression of MHC class II and CD86 on monocytes/macrophages in PBMCs from patients with NSCLC and dissociated tumor cultures in the presence of ricolinostat. PBMCs from patients with NSCLC (A and B) or dissociated tumors (C and D) were cultured for 24 or 72 hours, respectively, with ricolinostat or entinostat (2.5 μm), after which the phenotype of CD14+CD11b+ monocytes or CD45+CD68+CD11b+ tumor macrophages were assessed by FACS. DMSO was used as a control. Summary (left) or representative histograms (right) for the expression levels of (A and C) HLA-DR and (B and D) costimulatory molecule CD86 on gated CD3CD14+ monocytes from (A and B) PBMCs or (C and D) CD14CD45+CD68+CD11b+ macrophages in dissociated tumors after culture with indicated HDAC inhibitors. MFI, mean fluorescent intensity. E, Purified CD14+ cells from patient PBMCs that had been cultured with ricolinostat or entinostat for 24 hours were incubated with cell trace violet (CTV)–labeled purified T cells from allogeneic donor PBMCs for 6 days in the presence of 20 IU/mL of recombinant human IL2. The percentage of proliferation by the responder T cells was determined by CTV dilution in response to stimulation by CD14+ cells. Data represent the mean ±SEM of 5 patients (A–D) or 2 independent experiments (E). *, P < 0.05; **, P < 0.01. ns, not significant.

Figure 2.

Increased expression of MHC class II and CD86 on monocytes/macrophages in PBMCs from patients with NSCLC and dissociated tumor cultures in the presence of ricolinostat. PBMCs from patients with NSCLC (A and B) or dissociated tumors (C and D) were cultured for 24 or 72 hours, respectively, with ricolinostat or entinostat (2.5 μm), after which the phenotype of CD14+CD11b+ monocytes or CD45+CD68+CD11b+ tumor macrophages were assessed by FACS. DMSO was used as a control. Summary (left) or representative histograms (right) for the expression levels of (A and C) HLA-DR and (B and D) costimulatory molecule CD86 on gated CD3CD14+ monocytes from (A and B) PBMCs or (C and D) CD14CD45+CD68+CD11b+ macrophages in dissociated tumors after culture with indicated HDAC inhibitors. MFI, mean fluorescent intensity. E, Purified CD14+ cells from patient PBMCs that had been cultured with ricolinostat or entinostat for 24 hours were incubated with cell trace violet (CTV)–labeled purified T cells from allogeneic donor PBMCs for 6 days in the presence of 20 IU/mL of recombinant human IL2. The percentage of proliferation by the responder T cells was determined by CTV dilution in response to stimulation by CD14+ cells. Data represent the mean ±SEM of 5 patients (A–D) or 2 independent experiments (E). *, P < 0.05; **, P < 0.01. ns, not significant.

Close modal

Ricolinostat Promotes Quantitative and Phenotypic Changes That Facilitate Tumor-Infiltrating T-cell Activation and Function in an NSCLC Model

Given the importance of tumor-associated immune cells in shaping the course of tumor progression and antitumor immunity, we next investigated whether the phenotypic and functional changes seen in vitro on immune cells upon ricolinostat treatment are recapitulated in vivo, especially in the context of tumor antigens. To address this issue, we examined the effects of ricolinostat on immune cells infiltrating lung adenocarcinomas that spontaneously develop in two fully immunocompetent mouse models: (i) genetically engineered mice (GEM) harboring an activating KrasG12D mutation concurrent with Trp53 loss (designated KP), or (ii) GEM with T790M/L858R EGFR mutations (designated TL), as previously described (44).

Mice with established adenocarcinomas (as confirmed by MRI) were treated once daily with ricolinostat for 7 days prior to evaluation of dissociated tumors by comprehensive multi-parametric flow cytometry. Under this short-term treatment, we observed modest increases in the frequency of tumor-infiltrating CD8+ T cells and significant decreases in CD4+FOXP3+ Treg cells among lymphoid cells (Fig. 3A and B), resulting in a significant elevation of CD8:Treg ratios in treated tumors (Fig. 3C). Such increases in CD8:Treg ratio are often associated with enhanced antitumor responses (45). Although there were no substantial phenotypic changes with respect to inhibitory receptor molecule expression (Supplementary Fig. S6A–S6D), tumor-infiltrating CD8+ T cells in ricolinostat-treated mice exhibited increased expression of the activation marker CD69 (Fig. 3D) and enhanced secretion of the effector cytokine IFNγ following ex vivo stimulation compared with equivalent cells in the vehicle control group (Fig. 3E).

Figure 3.

Ricolinostat promotes phenotypic and functional changes in tumor-infiltrating T-cell subsets and macrophages that are consistent with immune activation. Single-cell suspensions of lung tumor nodules of KP and TL mice treated with ricolinostat or vehicle for 7 days were stained and then subjected to FACS analysis to assess proportions, phenotype, and function of CD45+ immune cell subsets. The single-cell suspensions were also stimulated ex vivo for intracellular cytokine production by T cells. Proportion of indicated lymphoid cell subsets in (A) KP or (B) TL tumors, and (C) ratio of CD8 to CD4+FOXP3+ Tregs in tumors of KP and TL mice. Proportion of (D) tumor-infiltrating CD8+ T cells expressing CD69 activation marker and (E) producing IFNγ after ex-vivo stimulation. F, The expression levels of MHC class II and G, CD86 costimulatory molecules on tumor-associated CD11c+CD11blo macrophages in vehicle- and ricolinostat-treated KP and TL mice. MFI, mean fluorescent intensity. Data are mean ±SEM of 5 to 9 mice per group. *, P < 0.05; ***, P < 0.001.

Figure 3.

Ricolinostat promotes phenotypic and functional changes in tumor-infiltrating T-cell subsets and macrophages that are consistent with immune activation. Single-cell suspensions of lung tumor nodules of KP and TL mice treated with ricolinostat or vehicle for 7 days were stained and then subjected to FACS analysis to assess proportions, phenotype, and function of CD45+ immune cell subsets. The single-cell suspensions were also stimulated ex vivo for intracellular cytokine production by T cells. Proportion of indicated lymphoid cell subsets in (A) KP or (B) TL tumors, and (C) ratio of CD8 to CD4+FOXP3+ Tregs in tumors of KP and TL mice. Proportion of (D) tumor-infiltrating CD8+ T cells expressing CD69 activation marker and (E) producing IFNγ after ex-vivo stimulation. F, The expression levels of MHC class II and G, CD86 costimulatory molecules on tumor-associated CD11c+CD11blo macrophages in vehicle- and ricolinostat-treated KP and TL mice. MFI, mean fluorescent intensity. Data are mean ±SEM of 5 to 9 mice per group. *, P < 0.05; ***, P < 0.001.

Close modal

Consistent with this finding, single-cell RNA sequencing (RNA-seq) revealed that a number of genes associated with T-cell activation and T-cell receptor (TCR) signaling including Cd69, Cd44, Cd247, and Zap70 were upregulated in tumor-infiltrating T cells from KP mice treated with ricolinostat when compared with equivalent cells derived from the tumors of the vehicle group (Supplementary Fig. S7). Among myeloid cell populations, the proportions of CD11c+CD11blo tumor-associated macrophages (TAM; Supplementary Fig. S8A) were comparable between the treatment group and vehicle controls (data not shown). However, under this short-term ricolinostat treatment, TAMs exhibited significantly elevated expression of MHC class II as well as costimulatory molecule CD86 (Fig. 3F and G; Supplementary Fig. S8B). Consistent with this phenotype, ricolinostat-exposed TAMs analyzed by single-cell RNA-seq showed higher expression of key genes associated with MHC class II expression including Cd74 and H2-Aa relative to their counterparts in the control mice (Supplementary Fig. S8C).

As an increase in acetylation of tubulin is considered to be one measure of pharmacologic inhibition of HDAC6 (42, 46, 47), we observed that ricolinostat treatment led to increased levels of acetyl α-tubulin (lysine 40) in the tumor (Supplementary Fig. S9), providing supporting evidence that changes occurring consequent to ricolinostat treatment stem from drug activity in tumor-associated immune cells.

Collectively, these findings from in vivo treatment of tumor-bearing mice are consistent with our in vitro observations in PBMCs and indicate that, upon ricolinostat treatment, tumor-infiltrating macrophages underwent changes analogous to a more mature APC phenotype that supports a more activated, effector profile of surrounding T cells.

BET Bromodomain Inhibitor JQ1 Disrupts Treg Signature Protein Expression and Attenuates Suppressive Function

Based on the positive immunodynamic and phenotypic changes mediated by ricolinostat in immune cell subsets, we reasoned that these effects should translate to augmented antitumor immune response. In initial efficacy studies, we found that ricolinostat treatment in immunocompetent KP mice provided only marginal therapeutic benefit; that is, minimal inhibition of tumor growth was observed relative to vehicle-treated controls (Supplementary Fig. S10). In an attempt to identify a partner agent to augment antitumor effect, we focused on JQ1, an inhibitor of the BET family of bromodomain proteins. Our choice of JQ1 was based on our hypothesis that, due to its indirect regulation of epigenetic “reader” proteins, which are key to immune cell biology, JQ1 should also potentiate immunomodulatory effects in the tumor microenvironment (48–50). To this end, we treated KP mice with JQ1 in order to determine whether it promotes changes in tumor-associated immune cells that are beneficial for tumor immunity. Analysis of lung tumor–infiltrated T cells revealed that CD4+FOXP3+ Tregs showed significant downregulation of FOXP3, CTLA4, and PD-1 upon JQ1 treatment, whereas little phenotypic changes were observed in the conventional T-cell subsets (Fig. 4A and B; Supplementary Fig. S11A and S11B). In addition, phospho-STAT5, a bromodomain-dependent target gene (51–53), was downregulated in tumors of JQ1-treated mice, specifically in Tregs, supporting the notion that JQ1 activity is dominant in Tregs relative to conventional T cells (Supplementary Fig. S12A–S12D).

Figure 4.

BET bromodomain inhibitor JQ1 disrupts signature protein expression and attenuates suppressive function of Tregs within lung tumors of treated mice. Single-cell suspensions generated from lung tumor nodules excised from KP mice treated with JQ1 or ricolinostat for 1 week were subjected to FACS analysis. Mice that received vehicle served as controls. A, Representative histograms and B, summary of expression levels of FOXP3, CTLA4, PD-1, and CD69 on tumor-infiltrating CD4+FOXP3+ Tregs. CD4+CD25hi Tregs sorted from the tumors of vehicle or JQ1-treated KP mice were cocultured with CFSE-labeled CD4+CD25 T cells isolated from the spleen of the same mouse. Cells were stimulated with α-CD3 in the presence of T-depleted splenocytes as APCs for 3 days. C, Representative CFSE profiles of proliferating T cells and D, summary of percentage of cells proliferating in the presence of tumor (left) or splenic (right) Tregs at indicated Treg:T-cell ratios. Data in B and D are mean ±SEM of 5 to 6 mice per group. *, P < 0.05.

Figure 4.

BET bromodomain inhibitor JQ1 disrupts signature protein expression and attenuates suppressive function of Tregs within lung tumors of treated mice. Single-cell suspensions generated from lung tumor nodules excised from KP mice treated with JQ1 or ricolinostat for 1 week were subjected to FACS analysis. Mice that received vehicle served as controls. A, Representative histograms and B, summary of expression levels of FOXP3, CTLA4, PD-1, and CD69 on tumor-infiltrating CD4+FOXP3+ Tregs. CD4+CD25hi Tregs sorted from the tumors of vehicle or JQ1-treated KP mice were cocultured with CFSE-labeled CD4+CD25 T cells isolated from the spleen of the same mouse. Cells were stimulated with α-CD3 in the presence of T-depleted splenocytes as APCs for 3 days. C, Representative CFSE profiles of proliferating T cells and D, summary of percentage of cells proliferating in the presence of tumor (left) or splenic (right) Tregs at indicated Treg:T-cell ratios. Data in B and D are mean ±SEM of 5 to 6 mice per group. *, P < 0.05.

Close modal

Furthermore, the phenotypic changes in Tregs noted above were accompanied by decreased suppressive function as Tregs isolated from the tumors of JQ1-treated mice were less efficient at inhibiting the proliferation of responding effector T cells when compared with their counterparts in vehicle-treated mice (Fig. 4C and D). In contrast, Tregs isolated from the spleens of these tumor-bearing mice were comparable in their inhibitory function (Supplementary Fig. S13), demonstrating that tumor–Tregs were more susceptible to JQ1 effect, possibly due to their unique and heightened expression of PD-1, CTLA4, and CD69 (Supplementary Fig. S14A and S14B). Consistent with these findings in the GEM, FOXP3+ Tregs present within dissociated tumors from patients with NSCLC cultured with JQ1 showed diminished expression levels of FOXP3, CTLA4, and PD-1 relative to those cultured in the presence of ricolinostat or control DMSO (Supplementary Fig. S15). In summary, these results indicate that JQ1 preferentially modulates Treg phenotype and function in tumor-infiltrating immune cells.

BET Bromodomain Inhibitor JQ1 Synergizes with Ricolinostat to Enhance Immune-Dependent Suppression of Tumor Growth

A recent report by Mazur and colleagues showed that JQ1 has an antitumor effect in KP mice in part due to downregulation of MYC and induction of apoptosis in tumor cells (36). The preceding data also demonstrate that JQ1 exhibits a unique effect on tumor-associated Tregs. We therefore hypothesized that these tumor and immune cell–specific JQ1 effects, when combined with the immunomodulatory effects of ricolinostat, should foster a tumor immune microenvironment that favors stronger stimulation of tumor-infiltrating T cells and an augmented antitumor response. Thus, mice with lung adenocarcinomas were treated with ricolinostat and/or JQ1 in long-term efficacy studies (Supplementary Fig. S16). Immunocompetent GEMs with established tumors (150–250 mm3) were administered a daily dose of ricolinostat and/or JQ1, and tumor growth kinetics were evaluated. In KP mice, ricolinostat or JQ1 monotherapy led to minimal or moderate delay in tumor growth, respectively, but combination treatment with both agents resulted in significant arrest of tumor growth (Fig. 5A–C). In alignment with these outcomes, vehicle-treated mice harbored large tumors with high-grade adenocarcinomas, whereas single and combination therapy tumors were mostly low-grade adenocarcinomas often exhibiting the lowest histopathology grade and increased cellular infiltrates (Supplementary Fig. S17A). In addition, tumors in the treatment groups, particularly the combination, displayed reduced Ki67 and elevated cleaved caspase-3 expression, signaling reduced proliferation and increased cell death, respectively (Supplementary Fig. S17B and S17C).

Figure 5.

JQ1 synergizes with ricolinostat to promote antitumor immunity. KP mice with tumor burdens of approximately 200 to 400 mm3 were injected i.p. once daily with ricolinostat alone, JQ1 alone, or the combination of the two drugs for 5 to 6 weeks with or without depleting antibodies against CD4 or CD8. Tumor growth was monitored weekly by MRI. A, Tumor growth kinetics. B, Change in tumor size after 2 weeks of treatment with indicated drugs and antibodies. C, Representative tumor MRI of mice treated with single agents or combination as indicated. D, Progression-free survival of KP mice in each treatment condition. For immune profiling, cohorts of KP mice were euthanized after 4 to 6 weeks of treatment, and tumor cell suspensions were subjected to FACS analysis. E, Percentage of CD8+ T cells that expressed CD69 within tumor-infiltrating CD45+ leukocytes. F, Immune cells isolated from tumors of mice treated with indicated drugs were stimulated ex vivo for 6 hours in the presence of golgi plug. Percentage of CD8+ T cells within tumor-infiltrating leukocytes that secreted IFNγ or G, expressed CD107α based on mean fluorescent intensity (MFI). H, CD25CD3+ T cells, CD45EPCAM+ tumor cells, and CD11bloCD11c+ TAMs were sorted from tumors of treated mice, and equivalent numbers of each population were cocultured for 2 days in the presence of tumor cell lysates. Graph represents percentage of dead EPCAM+ tumor cells as determined by populations that stained positive for viability dye. I, Summary of expression levels of MHC class II and J, CD86 on TAMs in the tumors of treated mice. Control mice received vehicle and Rat IgG isotype in the depletion studies. Data are mean ±SEM of 6 to 8 mice per group. *, P < 0.05; **, P <0.01; ***, P < 0.001.

Figure 5.

JQ1 synergizes with ricolinostat to promote antitumor immunity. KP mice with tumor burdens of approximately 200 to 400 mm3 were injected i.p. once daily with ricolinostat alone, JQ1 alone, or the combination of the two drugs for 5 to 6 weeks with or without depleting antibodies against CD4 or CD8. Tumor growth was monitored weekly by MRI. A, Tumor growth kinetics. B, Change in tumor size after 2 weeks of treatment with indicated drugs and antibodies. C, Representative tumor MRI of mice treated with single agents or combination as indicated. D, Progression-free survival of KP mice in each treatment condition. For immune profiling, cohorts of KP mice were euthanized after 4 to 6 weeks of treatment, and tumor cell suspensions were subjected to FACS analysis. E, Percentage of CD8+ T cells that expressed CD69 within tumor-infiltrating CD45+ leukocytes. F, Immune cells isolated from tumors of mice treated with indicated drugs were stimulated ex vivo for 6 hours in the presence of golgi plug. Percentage of CD8+ T cells within tumor-infiltrating leukocytes that secreted IFNγ or G, expressed CD107α based on mean fluorescent intensity (MFI). H, CD25CD3+ T cells, CD45EPCAM+ tumor cells, and CD11bloCD11c+ TAMs were sorted from tumors of treated mice, and equivalent numbers of each population were cocultured for 2 days in the presence of tumor cell lysates. Graph represents percentage of dead EPCAM+ tumor cells as determined by populations that stained positive for viability dye. I, Summary of expression levels of MHC class II and J, CD86 on TAMs in the tumors of treated mice. Control mice received vehicle and Rat IgG isotype in the depletion studies. Data are mean ±SEM of 6 to 8 mice per group. *, P < 0.05; **, P <0.01; ***, P < 0.001.

Close modal

These antitumor effects were largely dependent on CD8+ and CD4+ T cells, because tumor growth arrest was abrogated upon incorporation of CD8- or CD4-depleting antibodies in the combination treatment regimen (Fig. 5B). Furthermore, combination therapy led to longer progression-free survival compared with either agent alone (Fig. 5D), suggesting that JQ1 synergizes with ricolinostat to drive T cell–mediated antitumor response. Similar to KP mice, the combination was also efficacious in mice bearing mutant EGFR (TL) lung tumors, as well as Lewis lung carcinomas, as evidenced by marked inhibition of tumor growth in each of these models (Supplementary Fig. S18A and S18B). These data suggest that the therapeutic benefits of this combinatorial approach extend beyond lung adenocarcinomas driven by ectopic expression of driver mutations to those arising by spontaneous development.

Immune profiling of tumor nodules further revealed that, relative to the single-agent treatments, tumor-infiltrating CD8+ T cells under the combination therapy were more activated, as evidenced by increased expression of CD69 (Fig. 5E). In addition, these tumor-infiltrating CD8+ T cells demonstrated significantly enhanced secretion of effector cytokine IFNγ (Fig. 5F) as well as increased degranulation as indicated by CD107α staining (Fig. 5G).

Furthermore, when equivalent numbers of T cells, TAMs, and EPCAM+ tumor cells sorted from tumor nodules of ricolinostat and/or JQ1-treated mice were cocultured together with tumor cell lysate, we found a higher proportion of tumor cell death in the presence of T cells isolated from the tumors of dual agent–treated mice (Fig. 5H). Thus, combination treatment enhanced the tumor-killing potential of tumor-infiltrating T cells. With respect to the TAMs, ricolinostat treatment, either alone or in combination with JQ1, markedly promoted upregulation of MHC class II molecules and CD86 (Fig. 5I and J).

Ricolinostat and JQ1 Treatment Promotes Quantitative Changes in Tumor-Infiltrating Immune Cell Subsets That Signal a Reshaping of Tumor-Immune Cell Dynamics

Although we detected very little change in the proportions and phenotypes of T cells under acute (7 days) drug treatments with the exception of Tregs (Fig. 3A–C; Supplementary Fig. S6), we wondered whether the immune cell dynamics change over time, especially upon long-term treatment where we observed change in tumor size after prolonged drug exposure. Although the proportion of conventional T cells did not change dramatically within the tumors across treatment groups (Fig. 6A), tumors of KP mice treated with JQ1 alone or in combination with ricolinostat exhibited substantially reduced Treg numbers, yielding significantly increased CD8:Treg ratios, an effect that was not evident in peripheral Tregs (Fig. 6A and B; Supplementary Fig. S19A and S19C). In contrast, CD11c+ TAMs increased with more pronounced clusters of these cells within the tumors of mice treated with JQ1 and/or ricolinostat (Fig. 6C; Supplementary Fig. S19B), demonstrating peripheral reshaping of immune cell dynamics within the tumor in the presence of these drugs. Although their frequencies did not differ dramatically, tumor-infiltrating T cells, particularly CD8+ cells, in JQ1- and JQ1/ricolinostat-treated mice, exhibited significantly reduced expression levels of the immune checkpoint receptors PD-1, CTLA4, TIM3, and BTLA relative to the vehicle- and ricolinostat-treated groups (Supplementary Fig. S20A–S20D), indicating a correlation between the expression levels of inhibitory proteins on these T cells and the outcome of tumor growth kinetics.

Figure 6.

Dynamics of immune cells infiltrating the tumors of ricolinostat and/or JQ1-treated KP mice. Cell suspensions generated from tumor nodules of KP mice that were treated with vehicle, ricolinostat, and/or JQ1 for 5 to 6 weeks were subjected to FACS analysis to assess proportions of CD45+ T lymphocyte subsets. A, Proportion of CD4+, CD8+ conventional T cells, or CD4+FOXP3+ Tregs. B, Ratio of CD8 to CD4+FOXP3+ Tregs within the tumors. Frozen sections of fresh tumor nodules from these treated KP mice were stained for TAMs (CD11c+; red) and T cells (CD3+; green), and counterstained with DAPI (nuclei; blue). C, Representative immunofluorescent staining of tumor sections from mice treated as indicated. Images were captured on a Nikon Eclipse 80i fluorescence microscope equipped with CoolSNAP CCD camera and merged images created with NIS elements imaging software. Scale bar, 25 μm. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 6.

Dynamics of immune cells infiltrating the tumors of ricolinostat and/or JQ1-treated KP mice. Cell suspensions generated from tumor nodules of KP mice that were treated with vehicle, ricolinostat, and/or JQ1 for 5 to 6 weeks were subjected to FACS analysis to assess proportions of CD45+ T lymphocyte subsets. A, Proportion of CD4+, CD8+ conventional T cells, or CD4+FOXP3+ Tregs. B, Ratio of CD8 to CD4+FOXP3+ Tregs within the tumors. Frozen sections of fresh tumor nodules from these treated KP mice were stained for TAMs (CD11c+; red) and T cells (CD3+; green), and counterstained with DAPI (nuclei; blue). C, Representative immunofluorescent staining of tumor sections from mice treated as indicated. Images were captured on a Nikon Eclipse 80i fluorescence microscope equipped with CoolSNAP CCD camera and merged images created with NIS elements imaging software. Scale bar, 25 μm. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Close modal

Next, we wondered whether, despite the lack of substantial changes in immune cell dynamics under acute drug exposure, other changes were already being imprinted in tumor–immune cell populations. In this regard, we performed gene expression profiling to address this issue and as an approach to elucidate the mechanistic underpinnings for each drug as a sole agent. Our results revealed that TAMs sorted from the tumors of ricolinostat-treated mice, but not JQ1-treated mice, show increased gene transcripts for the MHC class II transactivator Ciita and MHC class II H2-DMb relative to vehicle controls (Supplementary Fig. S21A and S21B). On the other hand, JQ1 but not ricolinostat treatment remarkably led to downregulation of Treg-associated genes and/or genes encoding inhibitory coreceptors including Foxp3, Ctla4, Cd25, Ccr4, Pdcd1, and Btla in sorted tumor T cells (Supplementary Fig. S22A and S22B). This effect is not due to a reduction in Treg numbers within the analyzed T-cell samples, as these were comparable between ricolinostat and JQ1 treatments, and only slightly lower than the vehicle group (Supplementary Fig. S23). Interestingly, this gene expression pattern was largely absent in ricolinostat-exposed tumor T cells, which displayed a distinct profile with minimal changes in the gene transcripts for the above-mentioned inhibitory receptors (Supplementary Fig. S22A and S22B).

Taken together, these findings demonstrate that selective inhibition of HDAC6 by ricolinostat and bromodomain inhibition by JQ1 evokes a distinct gene transcriptional landscape that yields striking immunomodulatory features which, when combined, synergize to potentiate robust immune-mediated antitumor response in our described NSCLC models.

The utility of HDAC inhibitors as anticancer agents is gaining traction in oncology in part due to reported cytostatic effects on tumor cells (54). However, their effects on immune cells recruited to tumors are still not well understood. Such knowledge is critical, as cross-talk between immune cells and cancer cells shapes the course of tumor progression. Although pan-HDAC inhibitors are currently being evaluated as therapeutic options for a number of cancers, isozyme-selective HDAC inhibitors may possess unique properties. Indeed, we observed that ricolinostat, which is highly selective for HDAC6 but retains some inhibitory activities on class I HDACs as well, has striking immunomodulatory effects. In the present study, we characterized the effects of ricolinostat on immune cells and evaluated the therapeutic potential of combining this HDAC inhibitor with an inhibitor of BET bromodomain proteins in models of NSCLC.

The increased expression of MHC class II and costimulatory molecule CD86 that we observed on TAMs upon ricolinostat treatment in lung tumor–bearing mice indicates that ricolinostat possesses immunomodulatory properties that promote phenotypic changes favoring improved delivery of costimulatory signals and antigen presentation. Furthermore, the observed increase in CD69 expression on ricolinostat-exposed T cells also suggests that it may facilitate cellular events that lead to T-cell activation. Such an activating effect may be beneficial for amplification of TCR-induced antigen-driven T-cell priming, especially in the context of tumor-associated antigens. Although they did not focus on T cells, Sotomayor and colleagues also reported that HDAC6 inhibition improved APC function in a manner involving the regulation of STAT3/IL10 tolerogenic axis in APCs (42), raising the possibility that selective HDAC6 inhibition may differentially regulate activation and function in multiple cell types.

Our results suggest that ricolinostat has the potential to facilitate T-cell priming in patients with cancer. This notion is supported by the increased activation and enhanced effector function of tumor-infiltrating CD8+ T cells observed in tumor-bearing mice treated with ricolinostat. Although ricolinostat treatment alone did not robustly dampen tumor growth in our NSCLC models, it substantially increased T cell–mediated antitumor efficacy when coupled with JQ1. Thus, the enhancement of the therapeutic activity of ricolinostat by JQ1 cotreatment suggests that promoting antigen presentation and costimulation alone is not sufficient to achieve sustained antitumor responses. Rather, a significant reduction in inhibitory cellular mechanisms operative within the tumor bed as achieved by JQ1 therapy is fundamental to augmenting the therapeutic benefit of ricolinostat treatment. Nevertheless, the enhanced APC function conferred by ricolinostat warrants further exploration, as does its application in NSCLC treatment in tandem with vaccination strategies.

Given that CD4+FOXP3+ Tregs accumulate in many solid cancers (55), agents that target these cells may offer therapeutic benefits. In this regard, we found that JQ1, more than ricolinostat treatment, led to substantial reduction in tumor-infiltrating Treg fractions evidenced by reduced FOXP3 expression level, as well as reduced Treg function. These findings suggest that Tregs may be more susceptible to epigenetic modifications by bromodomain proteins than their conventional T-cell counterparts. Given that FOXP3, CTLA4, and PD-1 are critical for Treg stability and/or function (56–58), their diminished expression or transcriptional repression may compromise tumor–Treg survival and suppressive function as we observed in lung tumor–bearing, JQ1-treated mice. Although existing reports demonstrate that certain HDAC inhibitors affect Treg biology (35, 59, 60), and ricolinostat clearly evoked deleterious effects on patient tumor–Tregs in vitro, the more profound anti-Treg effect in our in vivo lung tumor model is attributed to JQ1.

Our gene expression profiling provides a snapshot of potential mechanisms of action by these drugs in our NSCLC investigational model: Ricolinostat promoted increased transcription of genes that are key components of MHC class II expression in tumor macrophages, molecular changes that support improved antigen presentation and T-cell priming (61). On the other hand, bromodomain inhibition by JQ1 remarkably led to a decrease in transcription of gene signatures associated with Treg maintenance/function which significantly overlap with those associated with immune checkpoint coreceptors (62–64). We postulate that ricolinostat plays a pleiotropic role on immune cells and functions by increasing activation of T cells and promoting APC function, whereas JQ1 contributes to reduction of Treg proportions in addition to disrupting the suppressive function of residual Tregs present in the tumor bed. The net effect is promotion of leukocyte-orchestrated antitumor response by way of optimizing antigen delivery, increasing costimulatory signals, and lowering inhibitory cellular mechanisms (Fig. 7).

Figure 7.

Proposed mechanism of action of ricolinostat and JQ1 in GEM models of NSCLC. In treatment-naïve lung tumors that develop in GEM models of NSCLC, inhibitory cellular mechanisms such as Tregs outweigh (co)-stimulatory signals resulting in impaired T-cell function. Ricolinostat promoted increased MHC and costimulatory molecules, which favors enhanced antigen presentation, whereas JQ1 led to a reduction in Treg numbers and function, effects that facilitate T-cell activation and function. The net result is a cooperative immunotherapeutic effect orchestrated by both agents to favor enhanced T-cell function, promoting robust antitumor immunity.

Figure 7.

Proposed mechanism of action of ricolinostat and JQ1 in GEM models of NSCLC. In treatment-naïve lung tumors that develop in GEM models of NSCLC, inhibitory cellular mechanisms such as Tregs outweigh (co)-stimulatory signals resulting in impaired T-cell function. Ricolinostat promoted increased MHC and costimulatory molecules, which favors enhanced antigen presentation, whereas JQ1 led to a reduction in Treg numbers and function, effects that facilitate T-cell activation and function. The net result is a cooperative immunotherapeutic effect orchestrated by both agents to favor enhanced T-cell function, promoting robust antitumor immunity.

Close modal

As the race continues toward finding effective therapeutic modalities for solid cancers, such as NSCLC, our findings show that the combination of ricolinostat and JQ1 leads to robust immune-mediated antitumor effects in lung tumors driven by different mutation genotypes, which is timely and important. This crucial contribution by the immune cells, particularly T cells, is underscored by the abrogation of tumor growth arrest seen under ricolinostat and JQ1 therapy upon T-cell depletion. In a nutshell, our findings highlight previously unknown immunomodulatory effects of these two classes of drugs in combination and support their therapeutic testing in lung adenocarcinoma. Although not the focus of this study, both agents have direct antitumor effects (36, 41, 48, 65). This raises the exciting possibility that ricolinostat and JQ1 could be utilized in combination as “hybrid” therapies that harness not just direct antitumor effects, but also indirect effects via peripheral reshaping of immune cell dynamics and function, converging to foster strong anticancer benefits. With JQ1 currently being evaluated in clinical trials for solid cancers, and ricolinostat, which has a good safety profile, presently in clinical development, both drugs are promising therapeutic agents that could be combined as epigenetics-based immunotherapeutics for the treatment of lung adenocarcinomas.

Mice

Genetically engineered mice harboring Kras+/LSL-G12DTrp53L/L or EGFRLSL-T790M/L858R have been previously described (44). For lung tumor induction, mice received 1 × 106 colony-forming units of Cre-encoding adenovirus intranasally at 5 to 6 weeks of age. Tumor formation (typically 6–8 weeks after inoculation) and volume were confirmed by MRI and quantified using 3-D Slicer software. C57BL/6 mice were purchased from the Jackson Laboratory. All mice were maintained under pathogen-free conditions in accordance with institutional guidelines for animal welfare, and experiments involving human tissues were conducted under approved, Dana-Farber/Harvard Cancer Center (DF/HCC) Institutional Review Board (IRB) protocol 02-180 with informed consent by patients.

Drugs and Reagents

Ricolinostat (ACY-1215), ACY-241, ACY-738, and ACY-775 were provided by Acetylon Pharmaceuticals. Tubastatin was obtained from ApexBio Technology, and Entinostat was obtained from Selleckchem. Ricolinostat was formulated in 10% DMSO in 5% dextrose and intraperitoneally administered to mice at 50 mg/kg daily. JQ1, kindly provided by James Bradner (Dana-Farber Cancer Institute, Boston, MA), was formulated in 10% DMSO in 10% 2-hydroxylpropyl β-cyclodextrin (Sigma-Aldrich) and intraperitoneally administered to mice at 25 mg/kg daily. Each drug's vehicle control represents DMSO diluted into the vehicle agent. Recombinant human cytokines IL2, IL4, and GM-CSF were purchased from Peprotech.

Antibodies

In vivo antibodies anti-CD4 (GK1.5), anti-CD8 (53-6.72), and rat IgG2b (LTF-2) were purchased from Bio X Cell. Mice were first administered two consecutive doses (400 μg/mouse) of antibodies at day −2 and day −1 and twice per week thereafter along with ricolinostat and JQ1 combination treatment. Antibodies utilized for flow cytometric analyses are listed in Supplementary Information.

Tumor Profiling and Flow Cytometry

To generate cell suspensions, tumor nodules were excised from lungs of KP and TL mice, cut into small pieces, and further dissociated in RPMI-1640 buffer containing 5% FBS, 100 IU/mL collagenase type IV (Invitrogen), and 50 μg/mL DNAse I (Roche) for 45 minutes at 37°C. After incubation, cells were treated with red blood cell lysis buffer and filtered through a 70-μm cell strainer. After centrifugation, cell pellets were resuspended in 1X PBS/2% FBS. About 0.5 to 1 × 106 cells were stained for surface markers in 1X PBS/2% FBS for 15 minutes at 4°C. Intracellular staining was performed for FOXP3 using the FOXP3 staining kit (eBioscience) according to the manufacturer's instructions. The Cytofix/Cytoperm Kit (BD Biosciences) was used for intracellular cytokine staining. Briefly, cells were stained for surface markers, including CD4, CD8, and CD3, and fixed and permeabilized, followed by intracellular staining with PE-conjugated anti-IFNγ, or isotype-matched mAbs. In all stained samples, dead cells were excluded using Live/Dead Fixable Dead Cell Staining Kit (Invitrogen). Cells were acquired on the LSR Fortessa (BD Biosciences) and analyzed with FlowJo software (Treestar). Gating strategy for FACS analyses is shown in Supplementary Fig. S24.

Cell Purifications and In Vitro Studies

PBMCs were isolated by Ficoll gradient from fresh or frozen blood samples obtained from healthy donors and patients. All subjects provided consent, and all procedures were followed under IRB-approved institutional protocols. When frozen, PBMCs were thawed and rested in 10% complete medium (RPMI-1640 supplemented with 10% pooled human serum, 100 μg/mL penicillin, 100 μg/mL streptomycin, and 2 mmol/L l-glutamine) at 37°C for 6 hours, after which cells (0.5 × 106) were cultured with DMSO as a control, or 2.5 μmol/L of individual HDAC inhibitors (ricolinostat, ACY-241, ACY-738, ACY-775, tubastatin, and entinostat) as indicated for 24 hours. For tumor cell cultures, freshly resected tumors from consented patients with NSCLC were dissociated into single-cell suspensions and cultured similarly for 72 hours in culture media supplemented with 20 IU/mL of IL2 and 20 ng/mL of GM-CSF. Cells were then washed and subjected to phenotypic analysis by FACS. In some experiments, cultures were maintained in 100 IU/mL of IL2 for 48 hours prior to analysis of Treg populations.

For human proliferation assays, T cells were isolated from PBMCs by positive selection using a T-cell enrichment kit (Miltenyi Biotec) per the manufacturer's instructions. Purified T cells (typically ∼90% pure) were then labeled with 2.5 μmol/L carboxyfluorescein succinimidyl ester (CFSE) or cell trace violet (Life Technologies) at 37°C for 7 minutes. Cells were washed thrice with RPMI + 10% FBS, counted, and used as responder cells. Purified CD14+ cells isolated from patient PBMCs by positive selection using CD14 microbeads (Miltenyi Biotec) were cultured with ricolinostat or entinostat for 24 hours, after which cells were washed thrice before coculture with allogeneic T cells in mixed lymphocyte reactions.

For mouse Treg suppression assays, CD4+CD25hi Tregs were sorted from spleen or tumor cell suspensions and cultured at graded numbers with CFSE-labeled CD4+CD25 cells (0.5 × 105) isolated from the spleen of the same mouse. T cell–depleted, mitomycin C–treated splenocytes (0.25 × 105) were added as APCs, and cultures were stimulated with soluble α-CD3 (clone 2C11; eBioscience) at 0.5 μg/mL concentration for 3 days.

For intracellular cytokine detection assays, immune cells from tumors of treated or control mice were obtained after Ficoll gradient separation as previously described (66). Cells (1 × 106) were cultured with PMA (50 ng) and ionomycin (500 ng) for 6 hours at 37°C. GolgiPlug (BD Pharmingen) and FITC-conjugated CD107α (BioLegend) were added for the last 5 hours of culture. For equivalent assays utilizing human PBMCs, 0.5 × 106 cells were similarly stimulated before analysis.

Immunohistochemistry

Tumors were extracted from lungs of tumor-bearing mice and fixed in 10% formalin overnight, after which they were embedded in paraffin and sectioned at the Specialized Histopathology core at Harvard Medical School. Prior to immunohistochemical staining, tumor sections were washed twice with Histo-Clear II (National Diagnostics), followed by two washes with 100% ethanol, and subsequent hydration with washes of 90%, 80%, 70%, and 50% ethanol. Tumor sections were heated in 10 mmol/L sodium citrate buffer (pH 6.0) for antigen unmasking. After cooling, sections were washed in de-ionized water, then incubated in 3% hydrogen peroxide for 10 minutes at room temperature followed by washes in deionized water and 1X PBS. For antigen blocking, sections were incubated in PBS buffer containing 0.5% Tween, 1% BSA plus 5% serum for 1 hour at room temperature. Sections were then stained in block buffer containing primary antibodies: anti-alpha tubulin (acetyl K40), clone EPR16772 at 1:100 (Abcam), or phospho-STAT5 alpha (Tyr694) at 1:200 (Invitrogen) overnight in a wet chamber at 4°C in the dark. The following day, sections were washed 3 times in 1X PBS and then stained with secondary biotinylated antibody in PBS blocking buffer for 1 hour at room temperature. Sections were washed 3 times with 1X PBS, and Elite Vectastain ABC Kit (Vector Laboratories) was applied for 30 minutes following the manufacturer's instructions at room temperature in the dark. Sections were washed with 1X PBS, developed with DAB reagent (Peroxidase Substrate Kit, Vector Laboratories), and counterstained with hematoxylin. Sections were then washed twice with deionized water followed by one wash with 1X PBS, and additional washes of increasing ethanol concentration for dehydration followed by incubation in Histo-Clear II. Slides were mounted with VectaMount Permanent Mounting Medium (Vector Laboratories) and covered with glass coverslips (Denville Scientific). After 24 hours, sections were viewed with an Olympus BX43 Trinocular Microscope. For all IHC quantitations, 10 randomly selected fields from at least 4 different tumors in each treatment group were used to quantitate the percentage of tissue positive for each marker using ImageJ software. Images were converted to a grayscale red-green-blue stack. Positive stain in the grayscale image was quantified at the appropriate threshold as %total image area positive for stain. Quantitation as a percentage of total tissue is shown where indicated.

Immunofluorescence

Tumors were excised from lungs of tumor-bearing mice, and embedded in and immediately snap-frozen in optimal cutting temperature (O.C.T.) compound (Fisher Healthcare). Frozen fresh tissues were cryosectioned at the Specialized Histopathology core at Harvard Medical School. Sections were fixed for 10 minutes in precooled (–20°C) acetone, washed 3 times in ice-cold 1X PBS, and blocked in 1X PBS with 10% BSA for 1 hour at room temperature. Sections were then stained with fluorochrome-conjugated antibodies: anti-mouse CD3 Alexa Fluor 488, clone 17A2 at 1:200, anti-mouse CD11c Alexa Fluor 594, clone N418 at 1:100, (BioLegend) overnight in a wet chamber at 4°C in the dark. Sections were subsequently washed 3 times in ice-cold 1X PBS and counterstained with DAPI (FluorPure Grade, Life Technologies) at 0.5 μg/mL for 7 minutes. Following this, sections were further washed thrice with 1X PBS and mounted with ProlongGold mounting media and coverslips (Corning) and allowed to cure for 24 hours at room temperature in the dark. Sections were imaged using a Nikon Eclipse 80i fluorescence microscope equipped with CoolSNAP CCD camera (Roper Scientific). NIS Elements Imaging software was used to create merged images.

Single-Cell RNA Sequencing of Tumor-Infiltrating TAMs and T Cells

TAMs (CD45+CD11c+CD11blo) and CD45+CD3+CD25lo/− T cells were sorted into 96-well plates containing 50 μL of 1X PBS with recombinant ribonuclease inhibitor (Life Technologies). Briefly, tumor cell suspensions were incubated with live/dead fixable dye after FcγR blocking and subsequently stained with antibodies against CD45, CD3, CD25, CD11b, and CD11c. The SmartSeq2 libraries were prepared according to the SmartSeq2 protocol (67, 68) with some modifications (69), the details of which are provided in Supplementary Information. Correlation between expression levels for genes of interest highlighted in this study is shown in Supplementary Fig. S25.

Statistical Analysis

Statistical significance was evaluated by the Student t test for comparisons between two groups (DMSO/vehicle vs. depicted drug) or one-way ANOVA for multigroup comparisons using GraphPad Prism software. P < 0.05 was considered statistically significant (*); P < 0.01 was considered very significant (**); and P < 0.001 was considered highly significant (***).

S.S. Jones has ownership interest (including patents) in Acetylon Pharmaceuticals. J. Bradner is President of NIBR at Novartis. S.N. Quayle has ownership interest (including patents) in Acetylon Pharmaceuticals. K.-K. Wong reports receiving commercial research support from Acetylon Pharmaceuticals. No potential conflicts of interest were disclosed by the other authors.

Conception and design: D.O. Adeegbe, Y. Liu, J.L. Guerriero, J. Bradner, S.N. Quayle, K.-K. Wong

Development of methodology: D.O. Adeegbe, Y. Liu, A.R. Aref, Y. Li, J. Castrillon, J.L. Guerriero, S. Han, K.-K. Wong

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): D.O. Adeegbe, Y. Liu, P.H. Lizotte, Y. Kamihara, Y. Li, S. Liu, X. Wang, T. Warner-Hatten, J. Castrillon, N. Poudel-Neupane, H. Zhang, S. Han, M.M. Awad, J. Ritz, K.-K. Wong

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): D.O. Adeegbe, Y. Liu, P.H. Lizotte, Y. Kamihara, R. Dries, J. Castrillon, G.-C. Yuan, H. Zhang, J.L. Guerriero, D.A. Barbie, J. Ritz, P.S. Hammerman, J. Bradner, K.-K. Wong

Writing, review, and/or revision of the manuscript: D.O. Adeegbe, Y. Liu, P.H. Lizotte, Y. Kamihara, C. Almonte, X. Wang, H. Zhang, J.L. Guerriero, M.M. Awad, D.A. Barbie, J. Ritz, S.S. Jones, P.S. Hammerman, S.N. Quayle, K.-K. Wong

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): Y. Liu, C. Almonte, N. Poudel-Neupane, P.S. Hammerman, K.-K. Wong

Study supervision: D.O. Adeegbe, Y. Liu, G.-C. Yuan, P.S. Hammerman, J. Bradner, K.-K. Wong

We thank Glenn Dranoff for critical reading of the manuscript and Elena Ivanova for technical assistance with immunofluorescence. The data presented in this article are tabulated in the main paper and in Supplementary Materials.

This work was supported by the National Cancer Institute R01 CA195740, CA163896, CA166480, CA122794, and CA140594 (to K.-K. Wong), CA205150 (to P.S. Hammerman), and 5RO1CA183560 (to J. Ritz); the Thoracic Foundation (to K.-K. Wong); the Gross-Loh Family Fund for Lung Cancer Research and Susan Spooner Family Lung Cancer Research Fund at Dana-Farber Cancer Institute (to K.-K. Wong); and research funding from Acetylon Pharmaceuticals.

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.
Adam
JK
,
Odhav
B
,
Bhoola
KD
. 
Immune responses in cancer
.
Pharmacol Ther
2003
;
99
:
113
32
.
2.
Domagala-Kulawik
J
. 
The role of the immune system in non-small cell lung carcinoma and potential for therapeutic intervention
.
Transl Lung Cancer Res
2015
;
4
:
177
90
.
3.
Quail
DF
,
Joyce
JA
. 
Microenvironmental regulation of tumor progression and metastasis
.
Nat Med
2013
;
19
:
1423
37
.
4.
Facciabene
A
,
Motz
GT
,
Coukos
G
. 
T-regulatory cells: key players in tumor immune escape and angiogenesis
.
Cancer Res
2012
;
72
:
2162
71
.
5.
Nishikawa
H
,
Sakaguchi
S
. 
Regulatory T cells in cancer immunotherapy
.
Curr Opin Immunol
2014
;
27
:
1
7
.
6.
Marvel
D
,
Gabrilovich
DI
. 
Myeloid-derived suppressor cells in the tumor microenvironment: Expect the unexpected
.
J Clin Invest
2015
;
125
:
3356
64
.
7.
Haberland
M
,
Montgomery
RL
,
Olson
EN
. 
The many roles of histone deacetylases in development and physiology: Implications for disease and therapy
.
Nat Rev Genet
2009
;
10
:
32
42
.
8.
Licciardi
PV
,
Karagiannis
TC
. 
Regulation of immune responses by histone deacetylase inhibitors
.
ISRN Hematol
2012
;
2012
:
690901
.
9.
Ropero
S
,
Esteller
M
. 
The role of histone deacetylases (HDACs) in human cancer
.
Mol Oncol
2007
;
1
:
19
25
.
10.
Li
Y
,
Seto
E
. 
HDACs and HDAC inhibitors in cancer development and therapy
.
Cold Spring Harb Perspect Med
2016
;
6
.
11.
Peinado
H
,
Ballestar
E
,
Esteller
M
,
Cano
A
. 
Snail mediates E-cadherin repression by the recruitment of the Sin3A/histone deacetylase 1 (HDAC1)/HDAC2 complex
.
Mol Cell Biol
2004
;
24
:
306
19
.
12.
von Burstin
J
,
Eser
S
,
Paul
MC
,
Seidler
B
,
Brandl
M
,
Messer
M
,
von Werder
A
, et al
E-cadherin regulates metastasis of pancreatic cancer in vivo and is suppressed by a SNAIL/HDAC1/HDAC2 repressor complex
.
Gastroenterology
2009
;
137
:
361
71
,
371 e1–5
.
13.
Clocchiatti
A
,
Florean
C
,
Brancolini
C
. 
Class IIa HDACs: from important roles in differentiation to possible implications in tumourigenesis
.
J Cell Mol Med
2011
;
15
:
1833
46
.
14.
Duong
V
,
Bret
C
,
Altucci
L
,
Mai
A
,
Duraffourd
C
,
Loubersac
J
, et al
Specific activity of class II histone deacetylases in human breast cancer cells
.
Mol Cancer Res
2008
;
6
:
1908
19
.
15.
Oehme
I
,
Linke
J-P
,
Böck
BC
,
Milde
T
,
Lodrini
M
,
Hartenstein
B
, et al
Histone deacetylase 10 promotes autophagy-mediated cell survival
.
Proc Natl Acad Sci U S A
2013
;
110
:
E2592
601
.
16.
Pandey
UB
,
Nie
Z
,
Batlevi
Y
,
McCray
BA
,
Ritson
GP
,
Nedelsky
NB
, et al
HDAC6 rescues neurodegeneration and provides an essential link between autophagy and the UPS
.
Nature
2007
;
447
:
859
63
.
17.
Gorospe
M
,
de Cabo
R
. 
AsSIRTing the DNA damage response
.
Trends Cell Biol
2008
;
18
:
77
83
.
18.
Mostoslavsky
R
,
Chua
KF
,
Lombard
DB
,
Pang
WW
,
Fischer
MR
,
Gellon
L
, et al
Genomic instability and aging-like phenotype in the absence of mammalian SIRT6
.
Cell
2006
;
124
:
315
29
.
19.
Wang
RH
,
Sengupta
K
,
Li
C
,
Kim
HS
,
Cao
L
,
Xiao
C
, et al
Impaired DNA damage response, genome instability, and tumorigenesis in SIRT1 mutant mice
.
Cancer Cell
2008
;
14
:
312
23
.
20.
Kroesen
M
,
Gielen
P
,
Brok
IC
,
Armandari
I
,
Hoogerbrugge
PM
,
Adema
GJ
. 
HDAC inhibitors and immunotherapy; a double edged sword?
Oncotarget
2014
;
5
:
6558
72
.
21.
West
AC
,
Smyth
MJ
,
Johnstone
RW
. 
The anticancer effects of HDAC inhibitors require the immune system
.
Oncoimmunology
2014
;
3
:
e27414
.
22.
San-Miguel
JF
,
Hungria
VT
,
Yoon
SS
,
Beksac
M
,
Dimopoulos
MA
,
Elghandour
A
, et al
Panobinostat plus bortezomib and dexamethasone versus placebo plus bortezomib and dexamethasone in patients with relapsed or relapsed and refractory multiple myeloma: A multicentre, randomised, double-blind phase 3 trial
.
Lancet Oncol
2014
;
15
:
1195
206
.
23.
Balasubramanian
S
,
Verner
E
,
Buggy
JJ
. 
Isoform-specific histone deacetylase inhibitors: The next step?
Cancer Lett
2009
;
280
:
211
21
.
24.
Gryder
BE
,
Sodji
QH
,
Oyelere
AK
. 
Targeted cancer therapy: Giving histone deacetylase inhibitors all they need to succeed
.
Future Med Chem
2012
;
4
:
505
24
.
25.
Ganai
SA
. Novel approaches towards designing of isoform-selective inhibitors against class II histone deacetylases: The acute requirement for targeted anticancer therapy
.
Curr Top Med Chem
2016
;
16
:
2441
52
.
26.
Bhadury
J
,
Nilsson
LM
,
Muralidharan
SV
,
Green
LC
,
Li
Z
,
Gesner
EM
, et al
BET and HDAC inhibitors induce similar genes and biological effects and synergize to kill in Myc-induced murine lymphoma
.
Proc Natl Acad Sci U S A
2014
;
111
:
E2721
30
.
27.
Ellis
L
,
Bots
M
,
Lindemann
RK
,
Bolden
JE
,
Newbold
A
,
Cluse
LA
, et al
The histone deacetylase inhibitors LAQ824 and LBH589 do not require death receptor signaling or a functional apoptosome to mediate tumor cell death or therapeutic efficacy
.
Blood
2009
;
114
:
380
93
.
28.
Lindemann
RK
,
Newbold
A
,
Whitecross
KF
,
Cluse
LA
,
Frew
AJ
,
Ellis
L
, et al
Analysis of the apoptotic and therapeutic activities of histone deacetylase inhibitors by using a mouse model of B cell lymphoma
.
Proc Natl Acad Sci U S A
2007
;
104
:
8071
6
.
29.
Chinnaiyan
P
,
Vallabhaneni
G
,
Armstrong
E
,
Huang
SM
,
Harari
PM
. 
Modulation of radiation response by histone deacetylase inhibition
.
Int J Radiat Oncol Biol Phys
2005
;
62
:
223
9
.
30.
Di Fazio
P
,
Schneider-Stock
R
,
Neureiter
D
,
Okamoto
K
,
Wissniowski
T
,
Gahr
S
, et al
The pan-deacetylase inhibitor panobinostat inhibits growth of hepatocellular carcinoma models by alternative pathways of apoptosis
.
Cell Oncol
2010
;
32
:
285
300
.
31.
Kim
K
,
Skora
AD
,
Li
Z
,
Liu
Q
,
Tam
AJ
,
Blosser
RL
, 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
:
11774
9
.
32.
Armeanu
S
,
Bitzer
M
,
Lauer
UM
,
Venturelli
S
,
Pathil
A
,
Krusch
M
, et al
Natural killer cell-mediated lysis of hepatoma cells via specific induction of NKG2D ligands by the histone deacetylase inhibitor sodium valproate
.
Cancer Res
2005
;
65
:
6321
9
.
33.
Berghuis
D
,
Schilham
MW
,
Vos
HI
,
Santos
SJ
,
Kloess
S
,
Buddingh
EP
, et al
Histone deacetylase inhibitors enhance expression of NKG2D ligands in Ewing sarcoma and sensitize for natural killer cell-mediated cytolysis
.
Clin Sarcoma Res
2012
;
2
:
8
.
34.
Magner
WJ
,
Kazim
AL
,
Stewart
C
,
Romano
MA
,
Catalano
G
,
Grande
C
, et al
Activation of MHC class I, II, and CD40 gene expression by histone deacetylase inhibitors
.
J Immunol
2000
;
165
:
7017
24
.
35.
van Loosdregt
J
,
Vercoulen
Y
,
Guichelaar
T
,
Gent
YY
,
Beekman
JM
,
van Beekum
O
, et al
Regulation of Treg functionality by acetylation-mediated Foxp3 protein stabilization
.
Blood
2010
;
115
:
965
74
.
36.
Mazur
PK
,
Herner
A
,
Mello
SS
,
Wirth
M
,
Hausmann
S
,
Sánchez-Rivera
FJ
, et al
Combined inhibition of BET family proteins and histone deacetylases as a potential epigenetics-based therapy for pancreatic ductal adenocarcinoma
.
Nat Med
2015
;
21
:
1163
71
.
37.
Das
R
,
Verma
R1
,
Sznol
M1
,
Boddupalli
CS1
,
Gettinger
SN1
,
Kluger
H
, et al
Combination therapy with anti-CTLA-4 and anti-PD-1 leads to distinct immunologic changes in vivo
.
J Immunol
2015
;
194
:
950
9
.
38.
Pardoll
DM
. The blockade of immune checkpoints in cancer immunotherapy
.
Nat Rev Cancer
2012
;
12
:
252
64
.
39.
Crisanti
MC
,
Wallace
AF
,
Kapoor
V
,
Vandermeers
F
,
Dowling
ML
,
Pereira
LP
, et al
The HDAC inhibitor panobinostat (LBH589) inhibits mesothelioma and lung cancer cells in vitro and in vivo with particular efficacy for small cell lung cancer
.
Mol Cancer Ther
2009
;
8
:
2221
31
.
40.
Pratap
J
,
Akech
J
,
Wixted
JJ
,
Szabo
G
,
Hussain
S
,
McGee-Lawrence
ME
, et al
The histone deacetylase inhibitor, vorinostat, reduces tumor growth at the metastatic bone site and associated osteolysis, but promotes normal bone loss
.
Mol Cancer Ther
2010
;
9
:
3210
20
.
41.
Woan
KV
,
Lienlaf
M
,
Perez-Villaroel
P
,
Lee
C
,
Cheng
F
,
Knox
T
, et al
Targeting histone deacetylase 6 mediates a dual anti-melanoma effect: Enhanced antitumor immunity and impaired cell proliferation
.
Mol Oncol
2015
;
9
:
1447
57
.
42.
Cheng
F
,
Lienlaf
M
,
Wang
HW
,
Perez-Villarroel
P
,
Lee
C
,
Woan
K
, et al
A novel role for histone deacetylase 6 in the regulation of the tolerogenic STAT3/IL-10 pathway in APCs
.
J Immunol
2014
;
193
:
2850
62
.
43.
Whiteside
TL
. 
Methods to monitor immune response and quality control
.
Dev Biol (Basel)
2004
;
116
:
219
28
;
discussion 229–36
.
44.
DuPage
M
,
Dooley
AL
,
Jacks
T
. 
Conditional mouse lung cancer models using adenoviral or lentiviral delivery of Cre recombinase
.
Nat Protoc
2009
;
4
:
1064
72
.
45.
Sato
E
,
Olson
SH
,
Ahn
J
,
Bundy
B
,
Nishikawa
H
,
Qian
F
, et al
Intraepithelial CD8+ tumor-infiltrating lymphocytes and a high CD8+/regulatory T cell ratio are associated with favorable prognosis in ovarian cancer
.
Proc Natl Acad Sci U S A
2005
;
102
:
18538
43
.
46.
Asthana
J
,
Kapoor
S
,
Mohan
R
,
Panda
D
. 
Inhibition of HDAC6 deacetylase activity increases its binding with microtubules and suppresses microtubule dynamic instability in MCF-7 cells
.
J Biol Chem
2013
;
288
:
22516
26
.
47.
Zhang
Y
,
Li
N
,
Caron
C
,
Matthias
G
,
Hess
D
,
Khochbin
S
, et al
HDAC-6 interacts with and deacetylates tubulin and microtubules in vivo
.
EMBO J
2003
;
22
:
1168
79
.
48.
Filippakopoulos
P
,
Knapp
S
. 
Targeting bromodomains: epigenetic readers of lysine acetylation
.
Nat Rev Drug Discov
2014
;
13
:
337
56
.
49.
Belkina
AC
,
Nikolajczyk
BS
,
Denis
GV
. 
BET protein function is required for inflammation: Brd2 genetic disruption and BET inhibitor JQ1 impair mouse macrophage inflammatory responses
.
J Immunol
2013
;
190
:
3670
8
.
50.
Prinjha
R
,
Tarakhovsky
A
. 
Chromatin targeting drugs in cancer and immunity
.
Genes Dev
2013
;
27
:
1731
8
.
51.
Liu
S
,
Walker
SR
,
Nelson
EA
,
Cerulli
R
,
Xiang
M
,
Toniolo
PA
, et al
Targeting STAT5 in hematologic malignancies through inhibition of the bromodomain and extra-terminal (BET) bromodomain protein BRD2
.
Mol Cancer Ther
2014
;
13
:
1194
205
.
52.
Ott
CJ
,
Kopp
N
,
Bird
L
,
Paranal
RM
,
Qi
J
,
Bowman
T
, et al
BET bromodomain inhibition targets both c-Myc and IL7R in high-risk acute lymphoblastic leukemia
.
Blood
2012
;
120
:
2843
52
.
53.
Toniolo
PA
,
Liu
S
,
Yeh
JE
,
Moraes-Vieira
PM
,
Walker
SR
,
Vafaizadeh
V
, et al
Inhibiting STAT5 by the BET bromodomain inhibitor JQ1 disrupts human dendritic cell maturation
.
J Immunol
2015
;
194
:
3180
90
.
54.
Falkenberg
KJ
,
Johnstone
RW
. 
Histone deacetylases and their inhibitors in cancer, neurological diseases and immune disorders
.
Nat Rev Drug Discov
2014
;
13
:
673
91
.
55.
Nishikawa
H
,
Sakaguchi
S
. 
Regulatory T cells in tumor immunity
.
Int J Cancer
2010
;
127
:
759
67
.
56.
Wing
K
,
Onishi
Y
,
Prieto-Martin
P
,
Yamaguchi
T
,
Miyara
M
,
Fehervari
Z
, et al
CTLA-4 control over Foxp3+ regulatory T cell function
.
Science
2008
;
322
:
271
5
.
57.
Williams
LM
,
Rudensky
AY
. 
Maintenance of the Foxp3-dependent developmental program in mature regulatory T cells requires continued expression of Foxp3
.
Nat Immunol
2007
;
8
:
277
84
.
58.
Zhang
B
,
Chikumaa
S
,
Hori
S
,
Fagarasan
S
,
Honjo
T
. 
Nonoverlapping roles of PD-1 and FoxP3 in maintaining immune tolerance in a novel autoimmune pancreatitis mouse model
.
Proc Natl Acad Sci U S A
2016
;
113
:
8490
5
.
59.
Shen
L
,
Ciesielski
M
,
Ramakrishnan
S
,
Miles
KM
,
Ellis
L
,
Sotomayor
P
, et al
Class I histone deacetylase inhibitor entinostat suppresses regulatory T cells and enhances immunotherapies in renal and prostate cancer models
.
PLoS One
2012
;
7
:
e30815
.
60.
Beier
UH
,
Wang
L
,
Han
R
,
Akimova
T
,
Liu
Y
,
Hancock
WW
. 
Histone deacetylases 6 and 9 and sirtuin-1 control Foxp3+ regulatory T cell function through shared and isoform-specific mechanisms
.
Sci Signal
2012
;
5
:
ra45
.
61.
den Haan
JM
,
Arens
R
,
van Zelm
MC
. 
The activation of the adaptive immune system: Cross-talk between antigen-presenting cells, T cells and B cells
.
Immunol Lett
2014
;
162
:
103
12
.
62.
Fuertes Marraco
SA
,
Neubert
NJ
,
Verdeil
G
,
Speiser
DE
. 
Inhibitory receptors beyond T cell exhaustion
.
Front Immunol
2015
;
6
:
310
.
63.
Liu
C
,
Workman
CJ
,
Vignali
DA
. 
Targeting regulatory T cells in tumors
.
FEBS J
2016
;
283
:
2731
48
.
64.
Montler
R
,
Bell
RB
,
Thalhofer
C
,
Leidner
R
,
Feng
Z
,
Fox
BA
, et al
OX40, PD-1 and CTLA-4 are selectively expressed on tumor-infiltrating T cells in head and neck cancer
.
Clin Transl Immunology
2016
;
5
:
e70
.
65.
Zhao
Y
,
Yu
D
,
Wu
H
,
Liu
H
,
Zhou
H
,
Gu
R
, et al
Anticancer activity of SAHA, a potent histone deacetylase inhibitor, in NCI-H460 human large-cell lung carcinoma cells in vitro and in vivo
.
Int J Oncol
2014
;
44
:
451
8
.
66.
Akbay
EA
,
Koyama
S
,
Carretero
J
,
Altabef
A
,
Tchaicha
JH
,
Christensen
CL
, et al
Activation of the PD-1 pathway contributes to immune escape in EGFR-driven lung tumors
.
Cancer Discov
2013
;
3
:
1355
63
.
67.
Picelli
S
,
Björklund
ÅK
,
Faridani
OR
,
Sagasser
S
,
Winberg
G
,
Sandberg
R
. 
Smart-seq2 for sensitive full-length transcriptome profiling in single cells
.
Nat Methods
2013
;
10
:
1096
8
.
68.
Picelli
S
,
Faridani
OR
,
Björklund
ÅK
,
Winberg
G
,
Sagasser
S
,
Sandberg
R
. 
Full-length RNA-seq from single cells using Smart-seq2
.
Nat Protoc
2014
;
9
:
171
81
.
69.
Trombetta
JJ
,
Gennert
D
,
Lu
D
,
Satija
R
,
Shalek
AK
,
Regev
A
. 
Preparation of single-cell RNA-Seq libraries for next generation sequencing
.
Curr Protoc Mol Biol
2014
;
107
:
4 22
1
17
.

Supplementary data