Although immunotherapy has revolutionized cancer care, patients with pancreatic ductal adenocarcinoma (PDA) rarely respond to these treatments, a failure that is attributed to poor infiltration and activation of T cells in the tumor microenvironment (TME). We performed an in vivo CRISPR screen and identified lysine demethylase 3A (KDM3A) as a potent epigenetic regulator of immunotherapy response in PDA. Mechanistically, KDM3A acts through Krueppel-like factor 5 (KLF5) and SMAD family member 4 (SMAD4) to regulate the expression of the epidermal growth factor receptor (EGFR). Ablation of KDM3A, KLF5, SMAD4, or EGFR in tumor cells altered the immune TME and sensitized tumors to combination immunotherapy, whereas treatment of established tumors with an EGFR inhibitor, erlotinib, prompted a dose-dependent increase in intratumoral T cells. This study defines an epigenetic–transcriptional mechanism by which tumor cells modulate their immune microenvironment and highlights the potential of EGFR inhibitors as immunotherapy sensitizers in PDA.

Significance:

PDA remains refractory to immunotherapies. Here, we performed an in vivo CRISPR screen and identified an epigenetic–transcriptional network that regulates antitumor immunity by converging on EGFR. Pharmacologic inhibition of EGFR is sufficient to rewire the immune microenvironment. These results offer a readily accessible immunotherapy-sensitizing strategy for PDA.

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

Immune-checkpoint blockade (ICB) has affected cancer care for multiple malignancies, including melanoma, non–small cell lung cancer, and renal cell carcinoma (1, 2). Nevertheless, most patients with these malignancies fail to achieve potent and durable responses to ICB-based immunotherapy. Moreover, many cancer types, including pancreatic ductal adenocarcinoma (PDA), exhibit negligible responses to ICB in the absence of microsatellite instability (1–3). This poor response has been attributed to immunosuppressive conditions in the tumor microenvironment (TME), including the scarcity and dysfunction of infiltrating CD8 T cells, as increased abundance and activity of tumor-infiltrating CD8 T cells enhance immunotherapy across different types of cancers (1, 4–6). It has become increasingly clear that tumor cells play a central role in establishing a permissive immune TME by releasing secreted factors that can actively drive the exclusion of T cells and resistance to immunotherapy (7–9). Although several cancer cell–intrinsic pathways have been reported to influence the composition of the immune TME (10–19), the key genetic and epigenetic drivers of this process remain largely unknown.

In the present study, we performed a targeted in vivo CRISPR screen and identified lysine demethylase 3A (KDM3A) as a potent epigenetic regulator of the immune TME. Cancer cells lacking KDM3A gave rise to tumors with marked alterations in lymphoid and myeloid infiltration and profound improvements in their response to combination immunotherapy. Unbiased molecular profiling suggested that KDM3A acts through Krueppel-like factor 5 (KLF5) and SMAD family member 4 (SMAD4) to regulate the expression of the epidermal growth factor receptor (EGFR), and functional studies revealed that ablation of tumor cell–intrinsic KLF5, SMAD4, or EGFR resulted in increased T-cell infiltration and improved immunotherapy responsiveness in previously refractory tumors. Moreover, treatment of established tumors with the EGFR inhibitor erlotinib was sufficient to promote a T cell–inflamed TME. Given the known requirement for EGFR signaling during pancreatic tumorigenesis (20–23) and erlotinib's status as an FDA-approved therapy for PDA (24), our results highlight the potential use of this drug as an immunotherapy-sensitizing agent for this treatment-refractory disease.

An In Vivo CRISPR Screen Identifies Epigenetic Regulators of Tumor Immunity in PDA

We previously established a group of clonal congenic primary tumor cell lines from the KPCY mouse model of spontaneous PDA (12, 25). Upon implantation into immunocompetent C57BL/6 mice, these clones recapitulate the heterogeneity of the immune TME seen in human PDA, with some clones giving rise to T cell–inflamed tumors and others giving rise to non–T cell–inflamed tumors (12). Importantly, genomic and epigenetic differences, rather than tumor mutational burden, account for these observed differences in immune cell infiltration (12). To better understand the epigenetic differences distinguishing T cell–inflamed and non–T cell–inflamed samples, we used ATAC-seq (assay for transposase-accessible chromatin coupled with sequencing) to compare the two groups. This analysis revealed different chromatin accessibility patterns between T cell–inflamed and non–T cell–inflamed samples, and that differences in the epigenetic landscape of cancer cells correlate with differences in the corresponding TME (Supplementary Fig. S1A and S1B). This finding is consistent with recent studies describing functions of epigenetic regulators in modulating antitumor immunity in other types of cancers (26–30).

T cell–inflamed and non–T cell–inflamed tumors derived from our panel of PDA tumor cell clones exhibit heterogeneous responses to a combination immunotherapy regimen consisting of gemcitabine (G), abraxane (A), CD40 agonistic antibody (F), CTLA4-blocking antibody (C), and PD-1–blocking antibody (P), hereafter referred to as GAFCP (10–12). GAFCP therapy improves clinical outcomes in preclinical mouse models of PDA (31, 32), and a similar regimen has shown promising results and is currently being tested in a phase II clinical trial of patients with metastatic PDA (NCT03214250; ref. 33). When challenged with GAFCP, PDA tumor cell clone 6419c5—a non–T cell-inflamed tumor with abundant myeloid cells and a paucity of tumor-infiltrating T cells (12)—responded modestly to GAFCP therapy, with no tumor regressions (Supplementary Fig. S1C). Nevertheless, RNA sequencing (RNA-seq) of tumor cells sorted from subcutaneously implanted tumors that were treated with either control or GAFCP therapy revealed significant transcriptional differences (Supplementary Fig. S1D and S1E). Specifically, tumor cells that persisted following GAFCP treatment were enriched for gene signatures similar to those enriched in untreated T cell–low clones compared with T cell–high clones (Supplementary Fig. S1F). These results indicate that the transcriptional networks underlying T-cell infiltration status mirror those that define sensitivity to immunotherapy. Thus, we concluded that this model would be a suitable platform for investigating other epigenetic and transcriptional regulators that shape the immune TME and influence response to immunotherapy in PDA.

To systematically interrogate tumor cell–intrinsic epigenetic factors modulating antitumor immunity, we devised an in vivo–focused CRISPR screen strategy (Fig. 1A). First, we engineered the 6419c5 tumor cell clone to stably express Cas9 and then transduced it with a single guide RNA (sgRNA) library (34) consisting of 5,100 sgRNAs targeting 850 epigenetic factors and RNA-binding factors, as well as 100 nontargeting control sgRNAs. We passaged the transduced cells in vitro for 12 days to allow for genome editing and then transplanted the cells into immunocompetent C57BL/6 mice. After 11 days, mice were treated with either vehicle or GAFCP therapy for an additional 12 to 14 days. Samples were collected at two time points prior to implantation (“T0” and “preinjection”) and from control- and GAFCP-treated tumors at the end of the experiment. Following sequencing, the representation of each sgRNA was measured for each sample and compared pairwise with all other samples (Supplementary Table S1).

Figure 1.

A genetic screen identifies KDM3A as a tumor cell–intrinsic epigenetic regulator suppressing antitumor immunity in pancreatic cancer. A, Diagram showing the strategy of an in vivo CRISPR-based genetic screen. B, Volcano plot showing the comparison between GAFCP and Ctrl groups. Each dot represents a gene whose sgRNAs were enriched in Ctrl (blue) or GAFCP (red) groups, with P < 0.05. Top candidate epigenetic regulators whose sgRNAs were enriched in the Ctrl group (P < 0.025 and >3 sgRNAs meeting a statistical cutoff) are highlighted. C, Top, representative immunofluorescent staining images of CD3+ T cells and Gr1+ myeloid cells in Kdm3a-WT or Kdm3a-KO tumors stained for CD3 or Gr1 (red) and YFP (green). Bottom: quantification of CD3+ T cells and Gr1+ myeloid cells comparing Kdm3a-WT and Kdm3a-KO tumors (counts taken from n = 5–8 images/group; data combined from two independent 6419c5 knockout clones). D, Flow-cytometric analysis of tumor-infiltrating immune cells in subcutaneously implanted Kdm3a-WT and Kdm3a-KO 6419c5 and 6694c2 tumors (n = 9–13/group). E, Flow-cytometric analysis of tumor-infiltrating immune cells in orthotopically implanted Kdm3a-WT and Kdm3a-KO 6419c5 and 6694c2 tumors (n = 5–10/group). F and G, Representative flow plots showing the abundance of myeloid cells (F) and T cells (G) within total CD45+ leukocytes in Kdm3a-WT (blue) or Kdm3a-KO (red) tumors. In C, statistical differences were calculated using unpaired Student t test. In D–E, statistical differences were calculated using multiple t test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01;***, P< 0.001. gMDSC, granulocytic myeloid-derived suppressor cell.

Figure 1.

A genetic screen identifies KDM3A as a tumor cell–intrinsic epigenetic regulator suppressing antitumor immunity in pancreatic cancer. A, Diagram showing the strategy of an in vivo CRISPR-based genetic screen. B, Volcano plot showing the comparison between GAFCP and Ctrl groups. Each dot represents a gene whose sgRNAs were enriched in Ctrl (blue) or GAFCP (red) groups, with P < 0.05. Top candidate epigenetic regulators whose sgRNAs were enriched in the Ctrl group (P < 0.025 and >3 sgRNAs meeting a statistical cutoff) are highlighted. C, Top, representative immunofluorescent staining images of CD3+ T cells and Gr1+ myeloid cells in Kdm3a-WT or Kdm3a-KO tumors stained for CD3 or Gr1 (red) and YFP (green). Bottom: quantification of CD3+ T cells and Gr1+ myeloid cells comparing Kdm3a-WT and Kdm3a-KO tumors (counts taken from n = 5–8 images/group; data combined from two independent 6419c5 knockout clones). D, Flow-cytometric analysis of tumor-infiltrating immune cells in subcutaneously implanted Kdm3a-WT and Kdm3a-KO 6419c5 and 6694c2 tumors (n = 9–13/group). E, Flow-cytometric analysis of tumor-infiltrating immune cells in orthotopically implanted Kdm3a-WT and Kdm3a-KO 6419c5 and 6694c2 tumors (n = 5–10/group). F and G, Representative flow plots showing the abundance of myeloid cells (F) and T cells (G) within total CD45+ leukocytes in Kdm3a-WT (blue) or Kdm3a-KO (red) tumors. In C, statistical differences were calculated using unpaired Student t test. In D–E, statistical differences were calculated using multiple t test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01;***, P< 0.001. gMDSC, granulocytic myeloid-derived suppressor cell.

Close modal

To confirm the robustness of the screening platform and strategy, we first identified genes whose sgRNAs were depleted from tumor cells in both the control (Ctrl) and GAFCP groups compared with the preinjection sample, as such genes would likely be important for tumor growth and progression in vivo. As sgRNAs targeting DNA methyltransferase 1 (Dnmt1) were the most highly depleted (Supplementary Fig. S2A and S2B; Supplementary Tables S2 and S3), we then knocked out Dnmt1 (Dnmt1-KO) in two PDA tumor cell clones, 6419c5 and 6694c2 (Supplementary Fig. S2C). This resulted in a decreased ability to form tumors (Supplementary Fig. S2D), confirming the importance of Dnmt1 in tumor growth. Moreover, flow cytometry revealed increased infiltration of T cells and decreased infiltration of myeloid cells in Dnmt1-KO tumors (Supplementary Fig. S2E and S2F). T-cell depletion partially restored the Dnmt1-KO cells' tumorigenic and growth capacities (Supplementary Fig. S2G and S2H), consistent with recent findings in other cancers showing that DNMT1 promotes tumor growth in part by suppressing antitumor immunity (35). Thus, our screening strategy can detect genes that control antitumor immunity to modulate tumor growth in vivo.

Tumor Cell–Intrinsic KDM3A Suppresses Antitumor Immunity in PDA

We next identified genes targeted by sgRNAs that are depleted from tumor cells in the GAFCP group compared with the Ctrl group. After filtering by P value (< 0.025) and the presence of >3 sgRNAs meeting a statistical cutoff, we nominated seven epigenetic regulators as potential suppressors of immunotherapy response (Fig. 1B; Supplementary Fig. S3A; Supplementary Table S4). To determine whether the hits in our screen regulate immune cell composition in the TME, we individually knocked out all seven genes in 6419c5 (Supplementary Fig. S3B). Given that the makeup of the immune TME, particularly T cells and myeloid cells, determines the responsiveness of PDA tumors to GAFCP treatment, we screened the tumors derived from the mutant clones for changes in immune infiltrates by flow cytometry. In five of the seven mutants (Prdm8, Atf2, Kdm3a, Ep400, and Prdm5), we observed a significant increase in intratumoral T cells and a significant decrease in intratumoral myeloid cells (Supplementary Fig. S3C), suggestive of a reprogramming of the immune landscape. We then deleted the top three hits from the screen (Prdm8, Atf2, and Kdm3a) in a second PDA tumor cell clone (6694c2) and observed a similar change in immune infiltration in Kdm3a mutants but not in the other two genes (Supplementary Fig. S3D–S3F). KDM3A has recently been shown to promote PDA progression in animal models (36). Moreover, an analysis of data from The Cancer Genome Atlas (TCGA) revealed KDM3A as the only gene among the seven candidates whose elevated expression in PDA correlates with poor clinical outcome (Supplementary Fig. S3G and S3H). Thus, we focused our attention on the immune TME regulatory activities of KDM3A.

Next, we performed a more detailed characterization of KDM3A's impact on the tumor immune microenvironment. In subcutaneous tumors, Kdm3a loss was associated with an increase in tumor-infiltrating T cells and dendritic cells (DC) and a concomitant decrease in myeloid cells compared with wild-type (WT) tumors (Fig. 1C and D; Supplementary Fig. S3I). This trend was confirmed in an orthotopic transplantation model (Fig. 1EG; Supplementary Fig. S3J). Recent studies identified subtypes of fibroblasts in PDA tumors that may have important implications for tumor immunity (37, 38). Utilizing surface markers identified by one of these reports, we analyzed the abundance of different fibroblast populations in Kdm3a-WT and Kdm3a-KO tumors and observed a decrease of myofibroblastic cancer-associated fibroblasts (myCAF) following deletion of tumor cell–intrinsic KDM3A (Supplementary Fig. S3K). Importantly, subcutaneously implanted Kdm3a-KO tumors acquired a marked sensitivity to GAFCP therapy compared with WT tumors, with most exhibiting a complete regression that led to a marked extension of survival (Fig. 2AC). These results demonstrate that KDM3A is an epigenetic suppressor of antitumor immunity and that loss of KDM3A reprograms the immune TME to a T cell–inflamed state, thereby conferring sensitivity to combination immunotherapy.

Figure 2.

Genetic ablation of tumor cell–intrinsic KDM3A promotes a T cell–inflamed TME in PDA. A, Growth curves of Kdm3a-WT or Kdm3a-KO tumors from two non–T cell–inflamed lines (6419c5 and 6694c2) with or without GAFCP treatment. B, Overall survival curves of animals with Kdm3a-WT or Kdm3a-KO tumors with or without GAFCP treatment. C, Waterfall plots showing change in tumor size of Kdm3a-WT or Kdm3a-KO tumors relative to the start of treatment measured with or without GAFCP treatment. In A–C, tumor cells were subcutaneously implanted into C57BL/6 mice (two knockout tumor cell clones for each parental PDA tumor cell clone were used, n = 5–16/group). In A, statistical differences were calculated using two-way ANOVA with multiple comparisons. In B, statistical differences were calculated using log-rank (Mantel–Cox) test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Figure 2.

Genetic ablation of tumor cell–intrinsic KDM3A promotes a T cell–inflamed TME in PDA. A, Growth curves of Kdm3a-WT or Kdm3a-KO tumors from two non–T cell–inflamed lines (6419c5 and 6694c2) with or without GAFCP treatment. B, Overall survival curves of animals with Kdm3a-WT or Kdm3a-KO tumors with or without GAFCP treatment. C, Waterfall plots showing change in tumor size of Kdm3a-WT or Kdm3a-KO tumors relative to the start of treatment measured with or without GAFCP treatment. In A–C, tumor cells were subcutaneously implanted into C57BL/6 mice (two knockout tumor cell clones for each parental PDA tumor cell clone were used, n = 5–16/group). In A, statistical differences were calculated using two-way ANOVA with multiple comparisons. In B, statistical differences were calculated using log-rank (Mantel–Cox) test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Close modal

Transcriptional and Epigenetic Profiling Identifies KLF5 and SMAD4 Transcription Factors as Potential Mediators of KDM3A Activity

In previous work, transcriptional profiling of sorted tumor cells identified several gene signatures that distinguish T cell–high clones from T cell–low clones (12). For example, T cell–low clones are enriched for TGFβ signaling and MYC signatures, whereas T cell–high clones are enriched for interferon response signatures (12). When we performed RNA-seq and gene set enrichment analysis (GSEA; ref. 39) on sorted tumor cells from WT versus Kdm3a-KO tumors (Supplementary Fig. S4A), we found that they enriched for molecular features similar to those that defined T cell–low and T cell–high tumor cells, including TGFβ signaling and MYC signatures in WT tumor cells and interferon response signatures in Kdm3a-KO tumor cells (Supplementary Fig. S4B and S4C). Kdm3a-KO tumor cells also expressed higher levels of PD-L1 and MHC class I in vivo (Supplementary Fig. S4D and S4E). When treated with interferon-gamma in culture, however, induction of MHC I and PD-L1 occurred to a similar degree in Kdm3a-KO cells and WT cells, suggesting that increased MHC class I expression by Kdm3a-KO cells is a consequence of the T cell–inflamed TME (Supplementary Fig. S4F). Collectively, these results suggest that KDM3A regulates gene networks like those that arise spontaneously during the formation of T cell–inflamed and non–T cell–inflamed tumors.

To identify transcriptional regulators mediating KDM3A's effects, we probed the epigenetic changes associated with KDM3A loss. Previous studies have identified KDM3A as a histone demethylase that colocalizes with and removes methyl groups from dimethylated histone H3 lysine residue 9 (H3K9me2; refs. 40–45), a histone modification associated with transcriptional repression (46–50). To unbiasedly examine global changes in histone modifications, we performed histone mass spectrometry analysis. No significant differences were observed between WT and Kdm3a-KO cells (Supplementary Fig. S4G), suggesting that KDM3A loss affects H3K9me2 levels at a local rather than a global level.

To characterize such locus-specific effects, we performed chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) for both H3K9me2 and KDM3A. We observed that both signals were largely distributed in distal intergenic regions (Supplementary Fig. S4H), as previously reported (51–54). Motif analysis revealed a strong signal for the consensus DNA binding sequence recognized by members of the KLF family of transcriptional regulators. KDM3A peaks that were lost in Kdm3a-KO cells enriched for the KLF motif (Fig. 3A), which is highly conserved across multiple KLF family members (55). Accordingly, genomic regions enriching KDM3A ChIP-seq signals also included the KLF motif, consistent with a previous study in prostate cancer cells (ref. 51; Fig. 3B). Furthermore, we found that peaks associated with H3K9ac, a classic marker of active promoters (56, 57), that were decreased in Kdm3a-KO cells also enriched for the KLF motif (Fig. 3C).

Figure 3.

Transcriptional and epigenetic profiling identifies KLF5 and SMAD4 as potential transcriptional regulators downstream KDM3A. A–C, Graphs showing the top enriched motifs of transcriptional regulators in genomic regions with lost H3K9me2, gained H3K9ac, and gained KDM3A ChIP-seq signals in Kdm3a-WT tumor cells compared with Kdm3a-KO tumor cells (q-value < 0.01 following the DREME- and Tomtom-based motif analysis, using ChIP-seq analysis of both 6419c5 and 6694c2 tumor cell clones). Two knockout clones from each parental tumor cell clone were pooled together for ChIP-seq analysis (n = 2/group, 8 samples in total). D, Volcano plot illustrating the differential gene expression between sorted YFP+ tumor cells from Kdm3a-WT or Kdm3a-KO tumors (6419c5 and 6694c2 tumor cell clones, two knockout clones of each parental tumor cell clone were used, n = 3/group and 18 samples in total). Each dot represents a gene [Kdm3a-WT (blue) or Kdm3a-KO (red), with Padj < 0.05]. Several factors known to drive a non–T cell–inflamed TME in PDA are highlighted. E, Dot plot showing the expression of KDM3A and KLF5 in the TCGA PAAD data set. F, Leading-edge plot of GSEA result highlighting Hallmark_TGF_beta_signaling. G, Venn diagram showing the overlapping gene targets of KDM3A, KLF5, and SMAD2/3. H, Leading-edge plot of GSEA result highlighting overlap genes of KDM3A, KLF5, and SMAD2/3 target genes.

Figure 3.

Transcriptional and epigenetic profiling identifies KLF5 and SMAD4 as potential transcriptional regulators downstream KDM3A. A–C, Graphs showing the top enriched motifs of transcriptional regulators in genomic regions with lost H3K9me2, gained H3K9ac, and gained KDM3A ChIP-seq signals in Kdm3a-WT tumor cells compared with Kdm3a-KO tumor cells (q-value < 0.01 following the DREME- and Tomtom-based motif analysis, using ChIP-seq analysis of both 6419c5 and 6694c2 tumor cell clones). Two knockout clones from each parental tumor cell clone were pooled together for ChIP-seq analysis (n = 2/group, 8 samples in total). D, Volcano plot illustrating the differential gene expression between sorted YFP+ tumor cells from Kdm3a-WT or Kdm3a-KO tumors (6419c5 and 6694c2 tumor cell clones, two knockout clones of each parental tumor cell clone were used, n = 3/group and 18 samples in total). Each dot represents a gene [Kdm3a-WT (blue) or Kdm3a-KO (red), with Padj < 0.05]. Several factors known to drive a non–T cell–inflamed TME in PDA are highlighted. E, Dot plot showing the expression of KDM3A and KLF5 in the TCGA PAAD data set. F, Leading-edge plot of GSEA result highlighting Hallmark_TGF_beta_signaling. G, Venn diagram showing the overlapping gene targets of KDM3A, KLF5, and SMAD2/3. H, Leading-edge plot of GSEA result highlighting overlap genes of KDM3A, KLF5, and SMAD2/3 target genes.

Close modal

We focused our attention on Klf5 as it was the most downregulated member of the KLF family in Kdm3a-KO cells (Fig. 3D; Supplementary Fig. S4I) and its expression positively correlated with the expression of KDM3A in human PDA samples (Fig. 3E). An analysis of TCGA data also showed that KLF5 expression in PDA correlates with poor clinical outcome (Supplementary Fig. S4J). Further examination of the RNA-seq data also identified SMAD4, which has recently been implicated in suppressing antitumor immunity in PDA and other types of cancers (11, 58, 59), as a possible mediator of KDM3A's immune-regulatory activity. Specifically, genes whose expression decreased upon KDM3A loss were enriched for both KLF and SMAD4 DNA binding motifs (Supplementary Fig. S4K), as well as the SMAD4-activating TGFβ Hallmark gene set (Fig. 3F). Interestingly, Kdm3a-KO cells expressed lower levels of factors that have previously been shown to promote a T cell–low microenvironment in PDA (Fig. 3D), including Epha2, Ptgs2, and Csf2 (11, 60, 61).

To determine whether KDM3A functionally interacts with KLF5 or SMAD4, we performed several experiments. First, we performed an immunoprecipitation (IP)–mass spectrometry experiment to identify KDM3A-interacting proteins. However, we did not observe any direct interaction between KDM3A and KLF5 or SMAD4 (Supplementary Table S5 and Supplementary Fig. S4L). To investigate whether KLF5, SMAD4, and KDM3A modulate a shared group of genes, we analyzed recently published ChIP-seq data sets of KLF5 and SMAD2/3 (closely related to SMAD4) in mouse pancreatic tumor cells (62). We observed that KDM3A, KLF5, and SMAD2/3 share significant overlap of gene targets (Fig. 3G). Importantly, GSEA indicated that the shared targets by these proteins are downregulated following ablation of KDM3A (Fig. 3H). Taken together, our transcriptomic and epigenomic analyses nominated KLF5 and SMAD4 as transcriptional regulators that may coordinate the output of KDM3A's epigenetic program.

Loss of Tumor Cell–Intrinsic KLF5 or SMAD4 Promotes a T Cell–Inflamed TME

To functionally examine the role of KLF5 and SMAD4 on antitumor immunity, we used CRISPR to knock out Klf5 and Smad4—along with several other transcription factors downregulated with Kdm3a loss—in the T cell–low clone 6694c2 (Supplementary Fig. S5A). These cells were implanted into immunocompetent mice, and immune infiltrates in the resultant tumors were examined by flow cytometry. Here, similar to our findings with KDM3A, we observed a significant increase of T cells and a concomitant decrease of myeloid cells with loss of KLF5 or SMAD4 (Fig. 4A; Supplementary Fig. S5B). Importantly, none of the other transcription factors representing potential KDM3A mediators yielded consistent changes in immune infiltration other than SMAD3, a binding factor of SMAD4 (Supplementary Fig. S5C). The effect of KLF5 on immune infiltration was also confirmed using an orthotopic model and in a second T cell–low clone (Fig. 4BD; Supplementary Fig. S5D–S5F). These results demonstrate that KLF5 and SMAD4 act within tumor cells to promote intratumoral infiltration of myeloid cells and exclusion of T cells.

Figure 4.

Genetic ablation of tumor cell–intrinsic KLF5 and SMAD proteins promotes a T cell–inflamed TME in PDA. A, Flow-cytometric analysis of immune cells in subcutaneously implanted WT, Klf5-KO, or Smad4-KO 6694c2 tumors (n = 10/group, two knockout clones/gene). B, Flow-cytometric analysis of immune cells in orthotopically implanted WT and Klf5-KO 6694c2 tumors (n = 5–10/group, two knockout clones/gene). C, Flow-cytometric analysis of immune cells in subcutaneously implanted WT and Klf5-KO 6419c5 tumors (n = 4–6/group, two knockout clones/gene). D, Representative flow plots showing the abundance of myeloid cells and T cells within total CD45+ leukocytes in Klf5-WT (blue) or Klf5-KO 6694c2 (red) tumors. E, Leading-edge plots from the GSEA highlighting the enrichment of KLF5 and SMAD4 gene signatures based on the RNA-seq analysis of Kdm3a-WT or Kdm3a-KO tumor cells (gene signatures were generated based on Padj < 0.01 and fold change > 2). F, Leading-edge plots from the GSEA highlighting the enrichment of the SMAD4 gene signature based on the RNA-seq analysis of Klf5-WT or Klf5-KO tumor cells (gene signatures were generated based on Padj < 0.01 and fold change > 2). G, Leading-edge plots from the GSEA highlighting the enrichment of the Hallmark_TGF_beta_signaling gene set. In A, statistical differences were calculated using one-way ANOVA with multiple comparisons. In B and C, statistical differences were calculated using an unpaired Student t test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Figure 4.

Genetic ablation of tumor cell–intrinsic KLF5 and SMAD proteins promotes a T cell–inflamed TME in PDA. A, Flow-cytometric analysis of immune cells in subcutaneously implanted WT, Klf5-KO, or Smad4-KO 6694c2 tumors (n = 10/group, two knockout clones/gene). B, Flow-cytometric analysis of immune cells in orthotopically implanted WT and Klf5-KO 6694c2 tumors (n = 5–10/group, two knockout clones/gene). C, Flow-cytometric analysis of immune cells in subcutaneously implanted WT and Klf5-KO 6419c5 tumors (n = 4–6/group, two knockout clones/gene). D, Representative flow plots showing the abundance of myeloid cells and T cells within total CD45+ leukocytes in Klf5-WT (blue) or Klf5-KO 6694c2 (red) tumors. E, Leading-edge plots from the GSEA highlighting the enrichment of KLF5 and SMAD4 gene signatures based on the RNA-seq analysis of Kdm3a-WT or Kdm3a-KO tumor cells (gene signatures were generated based on Padj < 0.01 and fold change > 2). F, Leading-edge plots from the GSEA highlighting the enrichment of the SMAD4 gene signature based on the RNA-seq analysis of Klf5-WT or Klf5-KO tumor cells (gene signatures were generated based on Padj < 0.01 and fold change > 2). G, Leading-edge plots from the GSEA highlighting the enrichment of the Hallmark_TGF_beta_signaling gene set. In A, statistical differences were calculated using one-way ANOVA with multiple comparisons. In B and C, statistical differences were calculated using an unpaired Student t test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Close modal

Next, we compared the transcriptional profiles of tumor cells sorted from subcutaneously implanted Kdm3a-KO, Smad4-KO, and Klf5-KO tumors (Supplementary Fig. S5G). GSEA revealed that loss of Klf5 or Smad4 enriched for the gene signature defining Kdm3a-KO cells (Fig. 4E). Furthermore, ablation of Klf5 resulted in transcriptional alterations that recapitulated those seen with loss of Smad4 (Fig. 4F). Moreover, deletion of either KLF5 or SMAD4 resulted in a negative enrichment of the TGFβ signature, as was observed for KDM3A (Fig. 4G). Collectively, these results suggest that KLF5, SMAD4, and KDM3A coordinately regulate transcription in tumor cells to suppress antitumor immunity.

The KDM3A–KLF5–SMAD4 Transcriptional Axis Converges on the EGFR

Because KDM3A, KLF5, and SMAD4 broadly regulate a shared transcriptional network, we sought to identify downstream pathways and gene targets on which all three converge. The 75 genes that were consistently decreased following ablation of Kdm3a, Klf5, or Smad4 enriched most significantly for the KRAS signaling gene set (Fig. 5AC; Supplementary Fig. S5H). Interestingly, genes encoding members of the EGF family of growth factors, as well as the EGFR itself, were among those showing the most consistent downregulation following ablation of these genes (Fig. 5B; Supplementary Fig. S6A). Several additional findings led us to focus on EGFR as a potential mediator of the KDM3A–KLF5–SMAD4 axis: (i) TGFβ treatment, which activates SMAD4 transcriptional activity, induced the expression of Egfr (Supplementary Fig. S6B); (ii) analysis of the Cancer Dependency Map Portal (DepMap; ref. 63) revealed EGFR to exhibit the greatest codependency with KLF5, but not SMAD4 or KDM3A, for survival across tumor lines as compared with any other gene (Fig. 5D); and (iii) in PDA patient samples, EGFR expression correlated with KDM3A and KLF5 expression (Fig. 5E) and worse clinical outcome (Supplementary Fig. S6C).

Figure 5.

Molecular profiling identifies EGFR as a tumor cell–intrinsic factor regulating the immune TME and sensitivity to combination immunotherapy in PDA. A, Venn diagram showing the overlap of differentially expressed genes in Klf5-KO, Kdm3a-KO, and Smad4-KO tumor cells compared with WT. B, Heat map showing expression of the 75 overlapped genes in A. C, Graph showing the top enriched hallmark gene sets based on the GSEA of overlapped genes identified in A. D, Dot plot showing the dependency scores for EGFR and KLF5 across all tumor cell lines in the Project Achilles/Cancer Dependency Map Portal (DepMap). E, Dot plot showing the expression of EGFR, KDM3A, and KLF5 in the TCGA PAAD data set. F, Flow-cytometric analysis of immune cells in subcutaneously implanted WT and Egfr-KO 6419c5 and 6694c2 tumors (n = 7–16/group, two knockout clones from each parental tumor cell clone). G, Flow-cytometric analysis of immune cells in orthotopically implanted WT and Egfr-KO 6419c5 and 6694c2 tumors (n = 5–18/group, two knockout clones from each parental tumor cell clone). H, Leading-edge plots from the GSEA highlighting the enrichment of EGFR gene signature based on RNA-seq analysis of Kdm3a-KO, Klf5-KO, and Smad4-KO tumor cells (gene signatures were generated based on Padj < 0.01 and fold change > 2). I, Graph showing top enriched functional pathways from GSEA using Hallmark gene sets. In F and G, statistical differences were calculated using multiple t test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Figure 5.

Molecular profiling identifies EGFR as a tumor cell–intrinsic factor regulating the immune TME and sensitivity to combination immunotherapy in PDA. A, Venn diagram showing the overlap of differentially expressed genes in Klf5-KO, Kdm3a-KO, and Smad4-KO tumor cells compared with WT. B, Heat map showing expression of the 75 overlapped genes in A. C, Graph showing the top enriched hallmark gene sets based on the GSEA of overlapped genes identified in A. D, Dot plot showing the dependency scores for EGFR and KLF5 across all tumor cell lines in the Project Achilles/Cancer Dependency Map Portal (DepMap). E, Dot plot showing the expression of EGFR, KDM3A, and KLF5 in the TCGA PAAD data set. F, Flow-cytometric analysis of immune cells in subcutaneously implanted WT and Egfr-KO 6419c5 and 6694c2 tumors (n = 7–16/group, two knockout clones from each parental tumor cell clone). G, Flow-cytometric analysis of immune cells in orthotopically implanted WT and Egfr-KO 6419c5 and 6694c2 tumors (n = 5–18/group, two knockout clones from each parental tumor cell clone). H, Leading-edge plots from the GSEA highlighting the enrichment of EGFR gene signature based on RNA-seq analysis of Kdm3a-KO, Klf5-KO, and Smad4-KO tumor cells (gene signatures were generated based on Padj < 0.01 and fold change > 2). I, Graph showing top enriched functional pathways from GSEA using Hallmark gene sets. In F and G, statistical differences were calculated using multiple t test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Close modal

To determine whether EGFR has a role in antitumor immunity, we knocked out EGFR in two PDA tumor cell clones and assessed the immune infiltration of the resulting tumors from both subcutaneous and orthotopic tumor models. As observed previously with KDM3A, SMAD4, and KLF5, loss of EGFR rewired the immune TME from a non–T cell–inflamed state to a T cell–inflamed state (Fig. 5F and G; Supplementary Fig. S6D–S6F). Remarkably, transcriptomic analysis of Egfr-KO cells revealed dysregulation of the same genes and pathways that were perturbed in Kdm3a, Klf5, and Smad4-KO cells (Fig. 5H and I), including downregulation of the TGFβ gene set and induction of interferon response gene sets. Importantly, among the 75 genes that were decreased by loss of Kdm3a, Klf5, and Smad4, 70 genes were also decreased following ablation of EGFR (Supplementary Fig. S6G). In addition, the tumor cell–intrinsic expression of two immune-modulatory secreted factors (CSF2 and PTGS2) that we and others have previously shown to promote a non–T cell–inflamed TME (11, 60, 61) in PDA was also decreased across all four mutants compared with WT tumor cells (Fig. 3A and 5B; Supplementary Fig. S6G). Taken together, these results suggest that KDM3A, KLF5, SMAD4, and EGFR converge on a conserved transcriptional program in cancer cells, which in turn influences the composition of the tumor immune microenvironment.

EGFR Inhibition Sensitizes Tumors to Combination Immunotherapy

Clones lacking Klf5, Smad4, or Egfr gave rise to tumors with a reduction in myeloid cells and an increase in T cells. Hence, we sought to determine whether these tumors responded to GAFCP immunotherapy, as we had earlier shown for Kdm3a mutants (Fig. 2). We found that tumors derived from Klf5-, Smad4-, or Egfr-KO tumor cells all became sensitized to GAFCP therapy, leading to tumor regression and improved animal survival, whereas parental cells showed minimal response (Fig. 6A). Therefore, the changes in immune infiltrates caused by loss of KLF5, SMAD4, or EGFR are sufficient to sensitize tumors to the effects of combination immunotherapy.

Figure 6.

Ablation of tumor cell–intrinsic KLF5, SMAD4, and EGFR increases the sensitivity of PDA tumors to a combination immunotherapy. Growth curves, waterfall plots, and overall survival curves of WT, Klf5-KO, Smad4-KO, and Egfr-KO tumors with or without GAFCP treatment are shown. Tumor cells were subcutaneously implanted into C57BL/6 mice (two knockout tumor cell clones were pooled together for this experiment, n = 5/group). Statistical differences were calculated using two-way ANOVA with multiple analysis (tumor growth curves) and log-rank (Mantel–Cox) test (survival curves).

Figure 6.

Ablation of tumor cell–intrinsic KLF5, SMAD4, and EGFR increases the sensitivity of PDA tumors to a combination immunotherapy. Growth curves, waterfall plots, and overall survival curves of WT, Klf5-KO, Smad4-KO, and Egfr-KO tumors with or without GAFCP treatment are shown. Tumor cells were subcutaneously implanted into C57BL/6 mice (two knockout tumor cell clones were pooled together for this experiment, n = 5/group). Statistical differences were calculated using two-way ANOVA with multiple analysis (tumor growth curves) and log-rank (Mantel–Cox) test (survival curves).

Close modal

EGFR signaling has been shown to play a central role in PDA biology (20–24, 64, 65). Erlotinib, a tyrosine kinase inhibitor (TKI) of EGFR, is approved for the treatment of patients with PDA based on a small but significant clinical benefit (24). Importantly, these modest clinical effects were observed in the absence of any adjunctive immunotherapies. Given that deletion of the Egfr gene results in a T cell–inflamed TME, we hypothesized that pharmacologic inhibition of EGFR might similarly shift the tumor immune landscape. Consistent with this hypothesis, we found that erlotinib treatment of subcutaneous or orthotopic 6694c2 tumors prompted a dramatic and dose-dependent increase in intratumoral T cells and decrease in intratumoral myeloid cells (Fig. 7A and B; Supplementary Fig. S7A and S7B). A more modest trend was observed in tumors derived from clone 6419c5 (Supplementary Fig. S7C and S7D), a difference that may be explained by the high abundance of granulocytic myeloid-derived suppressor cells (gMDSC) in these tumors, a feature that has been associated with resistance to EGFR inhibitors (66). Furthermore, we examined the immune infiltration of 6694c2 tumors following a short-term treatment of erlotinib (five days) and observed increased infiltration of T cells (Supplementary Fig. S7E). Therefore, pharmacologic inhibition of EGFR modulates the immune TME and may create a permissive environment for effective immunotherapy in non–T cell–inflamed (“cold”) PDA tumors.

Figure 7.

Erlotinib treatment increases T-cell infiltration in PDA tumors and enhances the sensitivity of PDA tumors to a combination immunotherapy. A, Flow-cytometric analysis of immune cells in subcutaneously implanted 6694c2 tumors treated with either vehicle or erlotinib (n = 5/group). B, Flow-cytometric analysis of immune cells in orthotopically implanted 6694c2 tumors treated with either vehicle or erlotinib (n = 5/group). For subcutaneous tumors, treatment started on day 10 after implantation, and tumors were collected for flow analysis 10 days later. For orthotopic tumors, treatment started on day 7 after implantation, and tumors were harvested for flow analysis 10 days later. C, Growth curve, overall survival curve, and waterfall plot showing change in tumor size relative to the start of treatment of tumors from a non–T cell–inflamed tumor cell clone, 6694c2, treated with vehicle, erlotinib (E), FCP, or erlotinib combined with FCP treatment (EFCP). In A, statistical differences were calculated using one-way ANOVA with multiple comparisons. In B, statistical differences were calculated using an unpaired Student t test. In C, statistical differences were calculated using two-way ANOVA with multiple comparisons and log-rank (Mantel–Cox) test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Figure 7.

Erlotinib treatment increases T-cell infiltration in PDA tumors and enhances the sensitivity of PDA tumors to a combination immunotherapy. A, Flow-cytometric analysis of immune cells in subcutaneously implanted 6694c2 tumors treated with either vehicle or erlotinib (n = 5/group). B, Flow-cytometric analysis of immune cells in orthotopically implanted 6694c2 tumors treated with either vehicle or erlotinib (n = 5/group). For subcutaneous tumors, treatment started on day 10 after implantation, and tumors were collected for flow analysis 10 days later. For orthotopic tumors, treatment started on day 7 after implantation, and tumors were harvested for flow analysis 10 days later. C, Growth curve, overall survival curve, and waterfall plot showing change in tumor size relative to the start of treatment of tumors from a non–T cell–inflamed tumor cell clone, 6694c2, treated with vehicle, erlotinib (E), FCP, or erlotinib combined with FCP treatment (EFCP). In A, statistical differences were calculated using one-way ANOVA with multiple comparisons. In B, statistical differences were calculated using an unpaired Student t test. In C, statistical differences were calculated using two-way ANOVA with multiple comparisons and log-rank (Mantel–Cox) test. P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P< 0.001.

Close modal

To directly examine whether erlotinib treatment could enhance the efficacy of combination immunotherapy, we combined erlotinib (E) with CD40 agonist (F), PD-1 blockade (P), and CTLA4 blockade (C). We observed a significant delay in tumor growth and improved clinical outcomes in mice treated with EFCP compared with the mice treated with only erlotinib (E) or FCP (Fig. 7C). This result directly demonstrates that erlotinib can enhance the responsiveness of a non–T cell–inflamed PDA tumors to combination immunotherapy. This delay in tumor growth translated into a significant survival benefit for animals treated with EFCP compared with vehicle control and a trend toward improved survival from the addition of erlotinib to FCP immunotherapy (P = 0.07; Fig. 7C). The effects on tumor growth persisted at later time points, as an analysis of animals surviving 20 days after initiation of therapy revealed persistently smaller tumors and one regression (Supplementary Fig. S7F). Together, these data suggest a potential new therapeutic opportunity for an otherwise treatment-refractory tumor type, providing a rationale for further preclinical, and even clinical, testing of this combination regimen.

Tumor heterogeneity describes distinct features that include both tumor cell–intrinsic and tumor cell–extrinsic properties (67). Functionally, heterogeneous immune TMEs lead to distinct responses to immunotherapy (4). Recent studies have highlighted the role of tumor cell–intrinsic factors in modulating both the abundance and activation of T cells and the consequences for productive antitumor immunity (7–9). A full understanding of these molecular regulators could inform novel therapeutic strategies to enhance the efficacy of existing immunotherapies. This principle is particularly relevant to PDA, which is refractory to most existing treatment options, including ICB-based immunotherapy (68). CRISPR loss-of-function genetic screening has been used to discover novel molecular regulators of antitumor immunity. Specifically, two recent studies utilized in vivo CRISPR screen strategies to unbiasedly identify potential therapeutic targets in melanoma and lung cancer (30, 69). The current study is, to our knowledge, the first in vivo screen to identify tumor cell–intrinsic factors that can mediate antitumor immunity in pancreatic cancer.

We and others have shown that genomic mutations are unlikely to explain immune TME heterogeneity (10–12, 19, 60, 61, 70, 71). Based on our observation that T cell–inflamed and non–T cell–inflamed tumor cells possess different epigenetic landscapes, we designed a screen to identify epigenetic factors important for immunotherapy responsiveness in vivo. Using this approach, together with transcriptional and epigenetic profiling and functional validation studies, we conclude that KDM3A, KLF5, and SMAD4 transcriptionally converge on EGFR to maintain a T cell–low TME and drive resistance to immunotherapies. Loss of SMAD4 is frequently observed in human PDA samples, suggesting that it functions as a tumor suppressor (70), yet our study suggests that ablation of tumor cell–intrinsic SMAD4 promotes antitumor immunity. These seemingly discrepant observations may be due to the different functions of SMAD4 during tumor initiation versus progression. Nevertheless, the fact that SMAD4 loss is common in PDA, whereas responses to ICB are virtually nonexistent in patients with microsatellite-stable PDA, indicates that SMAD4 loss is insufficient to sensitize PDA tumors to either PD-1 or CTLA4 inhibition. In the future, it will be important to evaluate whether immunotherapy-responsive PDA tumors are enriched for SMAD4 alterations. In addition, we observed that KDM3A loss resulted in a depletion of myCAFs from the TME, which could have an indirect effect on immune populations. Future work is needed to define the precise mechanism by which KDM3A, SMAD4, and KLF5 regulate antitumor immunity and whether this epigenetic–transcriptomic axis plays a role in other cancer types.

Our findings are in line with recent reports that small-molecule inhibitors of EGFR enhance the efficacy of immunotherapy. For example, one recent study demonstrated that EGFR inhibition increased the sensitivity of EGFR-mutated non–small cell lung cancer to ICB-based immunotherapy (72). Mechanistically, we found that EGFR inhibition led to a downregulation of known immunosuppressive molecules, including PTGS2 and CSF2 (Supplementary Fig S6G; refs. 11, 30, 60, 61, 73), as well as an induction of type I interferon signaling in tumor cells (Fig. 5I), matching the observation from another recent study (74). Thus, there are at least two possible mechanisms by which EGFR may alter T-cell infiltration: (i) inhibition of immunosuppressive molecules and (ii) activation of IFN signaling.

Importantly, the EGFR inhibitor erlotinib is approved by the FDA for treating pancreatic cancer in combination with gemcitabine (24). Although EGFR signaling has been shown to play an important role in pancreatic cancer progression (20–24), erlotinib's clinical benefit in PDA has been limited (22, 24). The current study provides a rationale for further examining whether the combination of EGFR inhibition with immunotherapies will elicit a more potent response, not just in pancreatic cancer but also in other KRAS-driven cancers such as lung cancer. In this respect, it is noteworthy that the combination therapy used in this study (GAFCP) is similar to one being evaluated in a phase II clinical trial for PDA (NCT03214250) and for which promising phase I data have been reported (33). As resistance to EGFR-targeting TKIs is frequently observed (75–77), a deeper understanding of the molecular mechanisms underlying EGFR-mediated control of antitumor immunity is needed. In particular, the finding that different PDA tumor clones exhibit variable immune responses following EGFR inhibition suggests that intertumoral heterogeneity could lead to varied responses in patients. Consequently, elucidating the mechanisms underlying EGFR-dependent immune infiltration may help identify those patients most likely to respond to the immune sensitizing effects of EGFR inhibitors as well as facilitate the development of alternative therapies for those nonresponders. Furthermore, it is also important to better understand the molecular changes in immune cells, especially CD8 T cells, following the treatment of EGFR inhibitors.

Accumulating evidence shows that different tumors exploit distinct cellular and molecular mechanisms to establish non–T cell–inflamed TMEs (12, 78). Consequently, further subtyping of non–T cell–inflamed tumors and identification of subtype-specific immunosuppressive modules will be important for precision medicine and improved efficacy of immunotherapy. This might be particularly relevant to our observation that different T cell–low PDA tumors vary in their sensitivity to erlotinib treatment. Because it has been previously shown that an increased abundance of gMDSCs is associated with therapy resistance, non–T cell–inflamed tumors having an overabundance of this cell type might derive particular benefit from treatment with anti-gMDSC agents (79).

Animals

All mouse experiment procedures used in this study were performed following the NIH guidelines. All mouse procedure protocols utilized in this study were in accordance with, and with the approval of, the Institutional Animal Care and Use Committee of the University of Pennsylvania. All WT C57BL/6 mice were purchased from The Jackson Laboratory and/or bred and maintained at the University of Pennsylvania.

PDA Tumor Cell Clones and Cell Lines

Mouse pancreatic tumor cell clones were generated in our lab recently (12) and were examined by the Research Animal Diagnostic Laboratory at the University of Missouri, using the Infectious Microbe PCR Amplification Test (IMPACT) II. PDA tumor cells were cultured in standard cell culture medium including DMEM (high glucose without sodium pyruvate, Corning, 10-017-CV) with 10% heat-inactivated FBS (Genesee Scientific; 25-514) and glutamax (Thermo 35050061). These tumor cell clones were used for fewer than 20 passages. YFP labeling of these tumor cell clones was derived from their genetic background (these tumor cell clones were derived from the KPCY mice as described in our previous report; ref. 25). We used 293T cells (Clontech, 632180) for lentivirus packaging.

Implantation of Tumor Cells

Pancreatic tumor cells were dissociated into single-cell suspension with 0.25% trypsin for three minutes at room temperature (Gibco), washed with serum-free DMEM twice, and counted in preparation for subcutaneous or orthotopic implantations. PDA tumor cells (2.0 × 105) were implanted subcutaneously and 5.0 × 104 PDA tumor cells were implanted orthotopically into the pancreas of six- to eight-week-old female C57BL/6 mice as previously described. Tumors were harvested 17 to 24 days following implantation for flow analysis, fixation, and other described experiments. Endpoint criteria included tumor volume exceeding 500 mm3, severe cachexia, or weakness and inactivity.

Subcutaneous Tumor Growth and Regression Assessments

For tumor growth assessment, tumors were measured every three to four days. Tumor length and width were examined with calipers, and tumor volumes were then calculated as length*width2/2. Tumor volumes of 500 mm3 were used as an endpoint for overall survival analysis following our mouse protocols. Tumor progression and/or regressions were calculated using the initial tumor size at the start of treatment to tumor size at the endpoint of analysis for calculation.

Treatments and T-cell Depletions in Mice

For therapy treatment, Gemcitabine (Hospira) and nab-paclitaxel (Abraxane, Celgene) were both purchased from the Hospital of the University of Pennsylvania Pharmacy. Gemcitabine (G) was procured as pharmaceutical-grade suspension at 38 mg/mL, and then diluted to 12 mg/mL in PBS and administered at 120 mg/kg via intraperitoneal (i.p.) injection. Abraxane (A) was purchased as a pharmaceutical-grade powder resuspended at 12 mg/mL in PBS for i.p. injection at dose of 120 mg/kg. Vehicle control mice received the equivalent to nab-paclitaxel dose of human albumin (huAlb; Sigma). For anti-CD40 agonist treatment, mice were injected i.p. with 100 μg of either agonistic CD40 rat anti-mouse IgG2a (clone FGK45, endotoxin-free) or the isotype control IgG2a (clone 2A3, Bio X Cell) 48 hours after chemotherapy. For checkpoint blockade treatment, mice were injected i.p. with 200 μg of anti–PD-1 (clone RMP 1-14; Bio X Cell) and 200 μg anti-CTLA4 (clone 9H10; Bio X Cell), beginning from described therapy start time point, with six and three doses, respectively. Control mice received the isotype control IgG2a (clone 2A3; Bio X Cell) on treatment days. Therapy was started when tumor was 3 to 5 mm size. CD4+ and CD8+ T cells were depleted using i.p. injections of 200 mg anti-CD4 (clone GK1.5; Bio X Cell) and anti-CD8 (clone 2.43; Bio X Cell), three days prior to tumor implantation, and every three days for the duration of the experiment. Control groups received IgG2b isotype control (Bio X Cell). For erlotinib treatment, we prepared 500 mg erlotinib (HY-50896) in DMSO and then diluted the drug in corn oil (5 mL DMSO + 45 mL corn oil) for a final concentration of 10 mg/mL. We delivered 100 mL and 200 mL of stock by gavage for 50 mg/kg and 100 doses, respectively, daily.

CRISPR/Cas9 Screen Analysis

Cas9 protein was introduced to the 6419c5 PDA tumor cell clone by transduction of the lentiCas9-Blast plasmid (Addgene plasmid #52962) and selected with 10 μg/mL blasticidin. Tumor cells were transduced with sgRNA lentivirus library at multiplicity of infection of 0.2. 1,000× coverage of the sgRNA library in tumor cells was maintained during in vitro cell culture. Enough tumor cells were implanted into mice to maintain 500× coverage of the sgRNA library. Tumor cells were collected for library preparation to maintain at least 500× coverage of the sgRNA library in each sample. Genomic DNA from collected tumor cells was harvested using the QIAamp DNA mini kit (Qiagen), and libraries were prepared following previously described protocols (80). Briefly, DNA containing sgRNAs were amplified over 25 cycles with Herculase II fusion DNA polymerase (Agilent) following the manufacturer-provided protocol with PCR.1 forward and reverse primer. Three micrograms of genomic DNA was used for each PCR reaction, and several PCR reactions were run in parallel for each collected sample for maintaining library coverage. Product of PCR.1 reaction was pooled for each sample, and 5 μL of each pooled PCR.1 sample was used as a template for PCR.2. For PCR.2, template was amplified over seven cycles with PCR.2 forward and reverse primers, and the resulting product was purified with the QIAquick PCR Purification Kit (Qiagen) followed by gel extraction with the QIAquick gel extraction kit (Qiagen). The barcoded libraries were pooled at an equal molar ratio and sequenced on a NextSeq500/550 (Illumina, 75 cycles High Output kit v2.0) to generate 75-bp single end reads. MAGeCK software was used for screen analysis (81). Briefly, the sequencing data were de-barcoded and merged, and the 20-bp sgRNA sequence was aligned to the reference sgRNA library without allowing for any mismatches. The read counts were calculated for each sgRNA using the method normalizing to the nontargeting sgRNAs. Differential analysis of sgRNA and targeted genes was also done following the MAGeCK instructions with standard parameters.

RNA-seq Analysis

RNA samples were extracted from sorted YFP+ tumor cells from subcutaneous tumors, using the Qiagen RNeasy Micro Kit following the manufacturer's instructions. RNA was sent out to a commercial company, Novogene, for library preparation and high-throughput sequencing using Illumina sequencers to generated paired-end results. Raw counts of gene transcripts were obtained using an alignment-independent tool, Salmon (https://combine-lab.github.io/salmon/), using standard settings. The raw count matrix was subsequently imported into R-studio (R version 3.3.3) and used as input for DESeq2 following the vignette of the package for normalization and differential gene-expression analysis. Salmon was used to normalize and quantitate gene expression in transcripts-per-million (tpm) through quasi-alignment. Differentially expressed genes were used as input for GSEA and other functional analysis. Differentially expressed genes were also used as input for online EnrichR analysis following the provided instructions (82, 83). Detailed scripts and parameters used for each step of the analysis provided by reasonable request to the authors.

ATAC-seq Analysis

For library construction, 50,000 tumor cells sorted from subcutaneously implanted tumors were used as input as previously described (12). The libraries were prepared as previously described with minor modifications (84). Briefly, nuclei were isolated from sorted cells using a solution of 10 mmol/L Tris-HC pH 7.4, 10 mmol/L NaCl, 3 mmol/L MgCl2, and 0.1% IGEPAL CA-630. Immediately following nuclei isolation, the transposition reaction was conducted using Tn5 transposase and TD buffer (Illumina) for 45 minutes at 37°C. Transposed DNA fragments were purified using a Qiagen MinElute Kit, barcoded with primers based on Illumina TruSeq indices, and PCR amplified for 12 cycles using NEBNext High Fidelity 2× PCR master mix (New England Biolabs). Libraries were purified by extraction from a 6% TBE gel, followed by column purification with the Qiagen PCR Cleanup Kit. Sequencing was performed using a HiSeq2500 (Illumina) on rapid run mode to generate 50-bp paired-end reads. For analysis, after adapter trimming with cutadapt, reads were aligned to the mouse reference genome (mm10) by Bowtie2. MACS2 was applied to identify accessible chromatin (85). Peaks were subsequently merged using BEDTools, and ATAC-seq read counts were calculated in the merged peaks for every sample (86). The resulting count table was used to identify differentially accessible loci with edgeR, which were then used for functional analysis by GREAT analysis (version 3.0.0; refs. 87–89). Detailed scripts and parameters used for each step of the analysis provided by reasonable request to the authors.

ChIP-seq Analysis

For ChIP, about 1.5 × 107 cells were cultured on 15-cm plates with normal culture medium. Medium was changed to DMEM without serum first, and then cells were cross-linked with 1.1% formaldehyde for 10 minutes at room temperature. Then samples were quenched with 125 mmol/L glycine for 5 minutes at room temperature on plate. Cells were then washed with cold PBS three times and then lysed with 1 mL nuclei isolation buffer for 15 minutes on ice (20 mmol/L PIPES, 85 mmol/L KCl, and 0.5% NP-40). Nuclei were spun down at 300 × g for 5 minutes at 4°C and resuspended in 500 μL ChIP lysis buffer with proteinase inhibitor (1% SDS, 10 mmol/L EDTA pH 8, 50 mmol/L Tris-HCl pH 8.0). Samples were then sonicated with a Diagenode Bioruptor at high-intensity setting for 15 × 5 minute cycles of 30 seconds on/off, to yield about 300-bp fragments. After centrifugation at 4°C with 13,500 rpm for 10 minutes, soluble chromatin was then diluted using ChIP dilution buffer (0.01% SDS, 1.2 mmol/L EDTA pH 8, 16.7 mmol/L Tris-HCl pH 8, 167 mmol/L NaCl, and 1.1% Triton X 100), and precleared with magnetic Protein G Dynabeads (Invitrogen, 10004D). IP with anti-H3K9me2 (Abcam, ab32521), anti-H3K9ac (Abcam, ab4441), anti-KDM3A (Novus Biologicals, NB100-77282) was performed overnight at 4°C with rotation, with 10% sample kept as input DNA. Products were recovered with magnetic Protein G Dynabeads and washed sequentially using ChIP Low Salt Buffer (0.1% SDS, 2.0 mmol/L EDTA pH 8.0, 20 mmol/L Tris-HCl pH 8.0, 150 mmol/L NaCl, 1.0% Triton X 100), ChIP High Salt Buffer (0.1% SDS, 2.0 mmol/L EDTA pH 8.0, 20 mmol/L Tris-HCl pH 8.0, 500 mmol/L NaCl, 1.0% Triton X 100), Lithium Chloride Buffer (1 mmol/L EDTA pH 8.0, 10 mmol/L Tris-HCl pH 8.0, 250 mmol/L LiCl, 1% NP-40, 1% sodium deoxycholate), and TE buffer (1 mmol/L EDTA pH 8.0, 10 mmol/L Tris-HCl pH 8.0) 10 minutes each buffer at room temperature. Samples were eluted twice using ChIP Elution Buffer (2% SDS, 100 mmol/L NaHCO3) at room temperature for 15 minutes each elution with rotation. Eluted sample DNA was incubated with RNase A (0.2 mg/mL) overnight at 65°C, then incubated with proteinase K (0.2 mg/mL) for 2 hours at 45°C. DNA samples were purified using the QIAquick PCR Purification Kit (Qiagen) and used for preparation of ChIP-seq libraries. Libraries were prepared following the NEBNext Ultra II DNA Library Prep Kit (New England Biolabs) protocol, using NEBNext Multiplex Oligos Index Primers Sets. Libraries were pooled and sequenced on the NextSeq500/550 (Illumina, 75 cycles High Output kit v2.0) to generate 75-bp single end reads. Sequencing reads were aligned to the mouse (mm10) using Bowtie2 with default settings (90). Resulting sam files were processed and marked with duplicates using SAMtools (91). Peak calling was performed using MACS2 (85) with default settings. Peaks were subsequently merged using BEDTools, and read counts were calculated in the merged peaks for every sample (86). The resulting count table was then used to find differentially genomic regions with edgeR (87, 88). Differential genomic regions were then transformed to bed files and used for genomic region annotation using the ChIP-seeker analysis tool (92) as well as the motif analysis using the MEME suite (93). The DREME analysis and the Tomtom analysis were used in the MEME suite for motif discovery and motif comparison (94, 95). The HOCOMOCO core collection of motifs (version 11) was used for the motif enrichment analysis. Detailed scripts and parameters used for each step of the analysis provided by reasonable request to the authors.

Cancer Dependency Map Portal Data Analysis

The depmap portal (https://depmap.org/portal/) and the CRISPR (Avana) Public 19Q3 data set were used for this analysis. No samples were excluded from the data set in this analysis. Following the depmap instruction, the dependency scores of genes were downloaded and then plotted as dot plots. Pearson correlation was calculated for all the plots using Prism.

Lentiviral Transduction of Tumor Cells for CRISPR-Mediated Ablation

The CRISPR vector, lentiCRISPR v2, was a gift from Feng Zhang (Addgene plasmid # 52961). The vector and pVSVg and psPAX2 lentiviral packaging plasmids (Addgene) were cotransfected into 293T cells (Clontech, 632180) using PEI reagent (Polysciences, 23966-2). Lentiviral particles were collected 48 hours after transfection and filtered for usage. Tumor cells transduced either with Cas9-Puro (control, ctrl) or Cas9-guide-Puro (knockout, KO; Ctrl is from Addgene, #52961, and KOs were cloned following the instruction from Addgene) were selected with 8 μg/mL puromycin (Invitrogen; A1113803). Single-cell clones were picked from bulk knockout cell line using single-cell sorting using BD Jazz FACS sorting machine in to 96-well plates. Knockout efficiencies were assessed by gene-specific qPCR analysis of target gene. CRISPR sgRNA sequences used were:

  • Dnmt1-sg1-TGA-AAC-TTC-ACC-TAG-TTC-CG; Dnmt1-sg2-CAA-ATA-GAT-CCC-CAA-GAT-CC;

  • Atf2-sg1-TGG-ACG-AAC-GAT-AGC-TGA-TG; Atf2-sg2-CGT-TCG-TCC-AGC-ATC-ATT-AC;

  • Brd1-sg1-TAA-TTG-AGC-TGC-TGC-GCA-AG; Brd1-sg2-GAG-CGT-AGC-AGC-ACA-GTT-AG;

  • Cbx7-sg1-GGA-AGA-GAG-GTC-CGA-AAC-CC; Cbx7-sg2-CCA-GAC-GTT-GTG-CAG-ACC-CC;

  • Ep400-sg1-TCG-GAA-CAT-GTA-GGG-CCG-GC; Ep400-sg2-CAG-AAG-CCG-ACC-CCT-TTA-AG;

  • Prdm5-sg1-GCT-CGA-TTC-ACA-CTG-TGG-AC; Prdm5-sg2-GAG-TTC-AAA-TGC-GAG-AAC-TG;

  • Prdm8-sg1-GCT-CTT-CCG-CTC-GTC-CGA-TG; Prdm8-sg2-TGT-TCT-ACC-GCT-CTC-TCC-GC;

  • Kdm3a-sg1-TAC-AGG-ATG-TTA-ACA-GTC-TT; Kdm3a-sg2-AAC-TCT-TCA-AGT-CAA-CTG-TG;

  • Tead2-sg1-GAA-GAC-GAG-AAC-GCG-AAA-GC; Tead2-sg2-TTC-GAG-CCA-AAA-CCT-GAA-TA;

  • Mycn-sg1-CGA-GTA-CGT-GCA-CGC-CCT-AC; Mycn-sg2-GCC-GAC-TCT-CGC-TTG-TTC-AC;

  • Prdm16-sg1-CTA-CGA-GAG-TCC-TCC-ATA-CC; Prdm16-sg2-CCC-GAT-TTC-CAT-CTT-CCG-CT;

  • Trim28-sg1-GGA-CCT-GCT-AAG-ACT-CGA-GA; Trim28-sg2-ATA-ATT-CTC-CAC-GAT-GTC-TT;

  • Klf5-sg1-CGC-GTG-TTT-CAG-ATC-GTC-TC; Klf5-sg2-TGC-GAA-CCC-GGC-CCG-CGA-CG;

  • Smad3-sg1-GCT-CCA-TGG-CCC-GTA-ATT-CA; Smad3-sg2-ACC-TAC-CTG-GAA-TAT-TGC-TC;

  • Smad4-sg1-GCC-AAG-TAA-TCG-CGC-ATC-AA; Smad4-sg2-TCC-GTT-GAT-GCG-CGA-TTA-CT;

  • Egfr-sg1-CGG-TCA-GAG-ATG-CGA-CCC-TC; Egfr-sg2-ACT-GCC-CAT-GCG-GAA-CTT-AC.

Quantitative PCR (qPCR) Analysis for Gene Expression

RNA was prepared from cultured PDA tumor cells using the RNeasy Mini Kit or RNeasy Micro Kit (Qiagen). cDNA was generated using the High-Capacity cDNA Reverse Transcription Kit from 1 μg RNA in 20 μL reaction volume and diluted 1:10 for qPCR analysis (Life Technologies). qPCR analysis was performed using 2 μL diluted cDNA with biological and technical replicates using SsoAdvanced SYBR reagent (Bio-Rad) and Bio-Rad qPCR platform, and results were normalized to the expression of Tbp using the Bio-Rad software. Primer sequences utilized for qPCR were:

  • Dnmt1-F-AAG-AAT-GGT-GTT-GTC-TAC-CGA-C; Dnmt1-R-CAT-CCA-GGT-TGC-TCC-CCT-TG;

  • Atf2-F-CCG-TTG-CTA-TTC-CTG-CAT-CAA; Atf2-R-TTG-CTT-CTG-ACT-GGA-CTG-GTT;

  • Brd1-F-AAC-ACT-GAC-CTA-CGC-ACA-AGC; Brd1-R-GCC-TCT-CGC-TGT-TCT-CCT-TAT-T;

  • Cbx7-F-TGC-GGA-AGG-GCA-AAG-TTG-AAT; Cbx7-R-ACA-AGG-CGA-GGG-TCC-AAG-A;

  • Ep400-F-CGG-TTC-TCA-GGA-TAA-ACT-GGC; Ep400-R-CAC-CTC-CGC-TCT-TGA-GCA-A;

  • Prdm5-F-GAG-AAG-CGA-ATG-CCT-GAA-GAC; Prdm5-R-CTC-CCA-CGT-ACC-TCC-CAC-A;

  • Prdm8-F-ATG-GAG-GAT-TCA-GGC-ATC-CAG; Prdm8-R-GGA-CCG-AAT-ATG-GCG-TTC-TCT;

  • Kdm3a-F-GTG-ACA-CAA-CCA-TTT-TCA-ACC-TG; Kdm3a-R-CAC-CCT-GTT-GGC-AGT-TCT-TC;

  • Tead2-F-GAA-GAC-GAG-AAC-GCG-AAA-GC; Tead2-R-GAT-GAG-CTG-TGC-CGA-AGA-CA;

  • Mycn-F-ACC-ATG-CCG-GGG-ATG-ATC-T; Mycn-R-AGC-ATC-TCC-GTA-GCC-CAA-TTC;

  • Prdm16-F-CCC-CAC-ATT-CCG-CTG-TGA-T; Prdm16-R-CTC-GCA-ATC-CTT-GCA-CTC-A;

  • Trim28-F-CGG-CGC-TAT-GGT-GGA-TTG-T; Trim28-R-GGT-TAG-CAT-CCT-GGG-AAT-CAG-AA;

  • Klf5-F-CCG-GAG-ACG-ATC-TGA-AAC-ACG; Klf5-R-GTT-GAT-GCT-GTA-AGG-TAT-GCC-T;

  • Smad3-F-CAC-GCA-GAA-CGT-GAA-CAC-C; Smad3-R-GGC-AGT-AGA-TAA-CGT-GAG-GGA;

  • Smad4-F-ACA-CCA-ACA-AGT-AAC-GAT-GCC; Smad4-R-GCA-AAG-GTT-TCA-CTT-TCC-CCA;

  • Egfr-F-GCC-ATC-TGG-GCC-AAA-GAT-ACC; Egfr-R-GTC-TTC-GCA-TGA-ATA-GGC-CAA-T;

  • Tbp-F-AGA-ACA-ATC-CAG-ACT-AGC-AGC-A; Tbp-R-GGG-AAC-TTC-ACA-TCA-CAG-CTC.

Flow Cytometry of Implanted Tumors

For flow-cytometric analyses, subcutaneous or orthotopic tumors were chopped into small pieces and digested in collagenase (1 mg/mL in DMEM; Sigma-Aldrich) at 37°C for 45 minutes and filtered through a 70-μm cell strainer to generate single-cell suspensions. Single-cell suspensions were then stained with fluorescence labeled antibodies on ice for 40 minutes and washed twice with cold PBS with 5% FBS for flow-cytometric analysis. Cells were then analyzed by flow cytometry using BD FACS (BD Biosciences) and FlowJo software (Treestar). CD44 was used as a marker for T-cell activation given our prior studies demonstrating that CD44 is consistently differentially expressed in CD8 T cells from naturally occurring T cell–low and T cell–high tumors (12).

Gating Strategies for Immune Cells.

Myeloid cells: live CD45+CD11b+; granulocytic (g)MDSCs/neutrophils: live CD45+CD11b+Gr1+; macrophages: live CD45+F4/80+ CD11b+; CD11c+ DCs: live CD45+F4/80 CD11c+; CD103+ DCs: live CD45+CD11bF4/80CD11c+ CD103+; T cells: live CD45+CD11b F4/80NKp46CD3+; CD4+ T cells: live CD45+CD11bF4/80NKp46CD3+CD4+; CD8+ T cells: live CD45+CD11bF4/80NKp46CD3+CD8+.

Antibodies Used for Flow Analysis.

CD279 (PD-1) FITC (29F.1A12) BioLegend 135214, CD335 (NKp46) PE (29A1.4) BioLegend 137604, CD103 PE/Dazzle 594 (2E7) BioLegend 121430, CD3e PE/Cy5 (145-2C11) BioLegend 100310, CD45 AF700 (30-F11) BioLegend 103128, CD8a PE/Cy7 (53-6.7) BioLegend 100722, Ly-6G V450 (1A8) BD 560603, H-2Kb/H-2Db (MHCI) AF647 (28–8-6) BioLegend 114612, F4/80 APC/Cy7 (BM8) BioLegend 123118, CD11b PerCP-Cy5.5 (M1/70) BD 550993, CD11c BV605 (N418) BioLegend 117334, Ly-6C BV570 (HK1.4) BioLegend 128030, CD4 BV650 (RM4-5) BioLegend 100546, CD44 APC BioLegend 103012.

TGFβ Treatment

PDA tumor cells were plated into 6-well plates, and 24 hours later were treated with either citric acid (control) or 10 ng/mL TGFβ (CST) over three days for RNA collection for RNA-seq analysis.

In Vitro Treatment with IFNγ for MHC-I Expression Analysis

Kdm3a-WT and Kdm3a-KO PDA tumor cells were treated with 100 ng/mL of IFNγ (PeproTech) in cell culture medium for 24 hours. Tumor cells were then trypsinized from culture plates and resuspended in PBS with 5% FBS for staining of antibodies. Tumor cells were then stained with MHC-I (BioLegend 114612) antibody on ice for 40 minutes and washed twice with cold PBS with 5% FBS for flow-cytometric analysis using BD FACS LSR machine (BD Biosciences) and FlowJo software (Treestar).

Immunofluorescent and IHC Staining

For CD3 and Gr1 staining, collected implanted tumor tissues were fixed in Zn-formalin for 24 hours and embedded in paraffin. Sections were deparaffinized, rehydrated, and prepared by antigen retrieval for six minutes each, and then blocked in 5% donkey serum (Sigma, D9663) for one hour at room temperature, incubated with primary antibodies overnight at 4°C, washed with 0.1% PBST (PBS with Tween-20), incubated with secondary antibodies for one hour at room temperature, and then washed and mounted. Slides were visualized and imaged using an Olympus IX71 inverted multicolor fluorescent microscope and a DP71 camera. For CD3 and Gr1 staining quantification, stained cells were counted for CD3+ T cells manually in five to eight fields per sample.

Primary antibodies included CD3 (Abcam ab5690, 1:100 dilution), Gr1 (eBioscience 14-5931-85, 1:50 dilution), and YFP (Abcam, ab6673). Secondary antibodies were purchased from Invitrogen (A-11055, A-21207, and A-21209) and were used as 1:250 dilution for all staining.

Immunoblotting

For whole-cell lysates, cells were washed with PBS and lysed in RIPA lysis buffer. Extracted proteins were separated by SDS-PAGE, transferred to PVDF membrane, blocked in 5% nonfat milk in PBS with 0.1% Tween-20, probed with primary antibodies, and detected with horseradish peroxidase–conjugated secondary antibodies (Jackson ImmunoResearch). Primary antibodies used include KDM3A (Novus Biologicals, NB100-77282) and alpha-TUBULIN (CST, 3873S).

Histone Extraction, Derivatization, Mass Spectrometry, and PTM Mark Quantification

For extraction, histones were acid-extracted according to standard protocols. Briefly, nuclei were isolated from cells lysed with a hypotonic lysis buffer. Histones were extracted from the nuclei with 0.2N H2SO4, precipitated in 33% trichloroacetic acid (TCA), washed with acetone, and resuspended in diH2O. For derivatization, histones were derivatized according to standard protocol (96). Briefly, 5 to 10 μg of acid-extracted histones were resuspended in 100 mmol/L ammonium bicarbonate (pH 8.0) and mixed with freshly prepared propionic anhydride with acetonitrile for 15 minutes at 37°C. Histones were then digested with trypsin (enzyme: sample ratio 1:20) overnight at 37°C. The peptides were then desalted and stored dried. They were resuspended in 0.1% formic acid just before mass spectrometry. For direct injection–MS, samples were placed in a TriVersa NanoMate (Advion) and acquired either manually or by using a sequence coordinated with MS acquisition by a contact closure. The NanoMate was set up with a spray voltage of 1.7 kV and a gas pressure of 0.5 psi. Samples were acquired in the Orbitrap Fusion Tribrid (Thermo Scientific). All scans were acquired in the orbitrap, at 240,000 resolution for the full MS and at 120,000 resolution for MS-MS. The AGC target for the tSIM-MSX scans was set to 10E6. The full description of the DI-MS acquisition method has been previously described (63). For histone peptide quantification, raw files were searched with a modified version of the software EpiProfile 2.064. Histone peptides are collected in MS scans, and isobaric peptides are collected in targeted preset MS-MS scans. The software reads the intensities from MS scans to calculate the percentage of all peptides with the same amino acid sequence. The unique fragment ions in the MS-MS scans are extracted to discriminate isobaric peptide intensities from the MS scans. The software EpiProfileLite is available at https://github.com/zfyuan/EpiProfileLite, with user guide.

Immunoprecipitation

For IP experiments, Dynabeads Protein G (Invitrogen) were first cross-linked with either rabbit IgG or KDM3A (Novus Biologicals, NB100-77282) antibodies at 6 μg antibodies/50 μl beads according to standard protocols. Briefly, antibodies and beads were mixed for 1 hour at room temperature, and then incubated with 20 mmol/L DMP in 0.2M triethanolamine for 40 minutes at room temperature. The reaction was quenched with 100 mmol/L Tris-HCl pH 8.0 for 1 hour at room temperature. The cross-linked beads were then washed 3× with 0.1M glycine pH 2.5, then washed 2× with PBST, and stored in PBST at 4°C for future use. For the IP, cultured cells were washed in cold PBS, and nuclei were isolated as previously described. Nuclei were then lysed in lysis buffer (20 mmol/L Tris pH 7.4, 0.15 M NaCl, 1.5 mmol/L MgCl2, 1 mmol/L EDTA, 1% glycerol, 1%NP-40) with an EDTA-free protease inhibitor mixture (Roche Applied Science) and precleared by incubating with Dynabeads Protein G (Invitrogen) for 1 hour at room temperature. Precleared nuclear lysates were then incubated with cross-linked beads overnight at 4°C. After washing with lysis buffer, protein was eluted 2× from the beads with 100 mmol/L glycine pH 2.5, which was then neutralized with 4 μL 1M Tris-HCl pH 8.5. Protein was then denatured by boiling with SDS protein loading buffer. Prepared samples were then subjected to either immunoblotting, as previously described, or LC/MS-MS analysis.

LC/MS-MS Analyses and Data Processing

LC/MS-MS analysis was performed as previously described (97) by the Proteomics and Metabolomics Facility at the Wistar Institute using a Q Exactive HF mass spectrometer (ThermoFisher Scientific) coupled with a Nano-ACQUITY UPLC system (Waters). Samples were reduced with TCEP, alkylated with iodoacetamide, digested in-gel with trypsin, and injected onto a UPLC Symmetry trap column (180 μm i.d. × 2 cm packed with 5 μm C18 resin; Waters). Tryptic peptides were separated by reverse-phase HPLC on a BEH C18 nanocapillary analytical column (75 μm i.d. × 25 cm, 1.7 μm particle size; Waters) using a 95-minute gradient formed by solvent A (0.1% formic acid in water) and solvent B (0.1% formic acid in acetonitrile). A 30-minute blank gradient was run between sample injections to minimize carryover. Eluted peptides were analyzed by the mass spectrometer set to repetitively scan m/z from 400 to 2,000 in positive ion mode. The full MS scan was collected at 60,000 resolution followed by data-dependent MS-MS scans at 15,000 resolution on the 20 most abundant ions exceeding a minimum threshold of 10,000. Peptide match was set as preferred; exclude isotopes option and charge-state screening were enabled to reject singly and unassigned charged ions.

Peptide sequences were identified using MaxQuant v1.6.8.0 (98). MS-MS spectra were searched against a UniProt mouse protein database (October 1, 2019) using full tryptic specificity with up to two missed cleavages, static carboxamidomethylation of Cys, and variable oxidation of MET and protein N-terminal acetylation. Consensus identification lists were generated with false discovery rates of 1% at protein and peptide levels.

Data Resources

The RNA-seq and ChIP-seq data, fastq files, and processed count tables have been deposited at the Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) with following access numbers: GSE146832, GSE146833, GSE146834, and GSE156889.

Software and Statistical Analysis

PRISM software and R were used for data processing, statistical analysis, and result visualization (http://www.graphpad.com). The R language and environment for graphics (https://www.r-project.org) was used in this study for the bioinformatics analysis of CRISPR screen, RNA-seq, and ChIP-seq data. The R packages used for all analysis described in this manuscript were from the Bioconductor and CRAN. Statistical comparisons between two groups were performed using the Student unpaired t test. For comparisons between multiple groups, one-way ANOVA with Tukey HSD post-test was used. For survival comparison between two groups, log-rank P values of Kaplan–Meier curves were determined in GraphPad Prism 8 (GraphPad). On graphs, bars represent either range or standard error of mean, as indicated in legends. For all figures, P < 0.05 was considered statistically significant; *, P < 0.05; **, P < 0.01; ***, P < 0.001.

R.H. Vonderheide reports grants and personal fees from Apexigen, Janssen, and Lilly, grants from Fibrogen and Inovio, and personal fees from AstraZeneca, Celgene, Genentech, Medimmune, Merck, and Verastem during the conduct of the study; in addition, R.H. Vonderheide is an inventor on a patent for Universal Cancer Antigens, an inventor on a patent for Cellular Therapy issued to formerly Novartis, and an inventor on a research antibody licensed to multiple merchants. B.Z. Stanger reports grants from National Institutes of Health during the conduct of the study; personal fees from iTeos Therapeutics and grants from Boehringer-Ingelheim outside the submitted work. No disclosures were reported by the other authors.

J. Li: Conceptualization, data curation, formal analysis, investigation, methodology, writing–review and editing. S. Yuan: Conceptualization, data curation, formal analysis, investigation, methodology, writing–review and editing. R.J. Norgard: Investigation. F. Yan: Conceptualization and investigation. Y.H. Sun: Resources and investigation. I.-K. Kim: Investigation and methodology. A.J. Merrell: Investigation. Y. Sela: Investigation. Y. Jiang: Investigation. N.V. Bhanu: Resources, formal analysis, and investigation. B.A. Garcia: Conceptualization and resources. R.H. Vonderheide: Conceptualization and resources. A. Blanco: Conceptualization, supervision, methodology, writing–review and editing. B.Z. Stanger: Conceptualization, supervision, funding acquisition, writing–review and editing.

This work was supported by grants from NIH (CA229803, DK083355, and CA224970), the Abramson Family Cancer Research Institute, the Abramson Cancer Center, and the NIH/Penn Center for Molecular Studies in Digestive and Liver Diseases. B.A. Garcia acknowledges NIH grants CA196539 and U54CA232568, as well as St. Jude Collaborative A11576. We thank members of the Stanger laboratories for technical help and scientific discussions. We also thank members of the flow cytometry and cell sorting resource laboratory for technical help with flow analysis and sorting.

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.
Ribas
A
,
Wolchok
JD
. 
Cancer immunotherapy using checkpoint blockade
.
Science
2018
;
359
:
1350
5
.
2.
Sanmamed
MF
,
Chen
L
. 
A paradigm shift in cancer immunotherapy: from enhancement to normalization
.
Cell
2018
;
175
:
313
26
.
3.
Le
DT
,
Durham
JN
,
Smith
KN
,
Wang
H
,
Bartlett
BR
,
Aulakh
LK
, et al
Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade
.
Science
2017
;
357
:
409
13
.
4.
Binnewies
M
,
Roberts
EW
,
Kersten
K
,
Chan
V
,
Fearon
DF
,
Merad
M
, et al
Understanding the tumor immune microenvironment (TIME) for effective therapy
.
Nat Med
2018
;
24
:
541
50
.
5.
Tumeh
PC
,
Harview
CL
,
Yearley
JH
,
Shintaku
IP
,
Taylor
EJM
,
Robert
L
, et al
PD-1 blockade induces responses by inhibiting adaptive immune resistance
.
Nature
2014
;
515
:
568
71
.
6.
Spranger
S
. 
Mechanisms of tumor escape in the context of the T-cell–inflamed and the non-T-cell-inflamed tumor microenvironment
.
Int Immunol
2016
;
28
:
383
91
.
7.
Li
J
,
Stanger
BZ
. 
The tumor as organizer model
.
Science
2019
;
363
:
1038
9
.
8.
Spranger
S
,
Gajewski
TF
. 
Impact of oncogenic pathways on evasion of antitumour immune responses
.
Nat Rev Cancer
2018
;
18
:
139
47
.
9.
Wellenstein
MD
,
de Visser
KE
. 
Cancer-cell-intrinsic mechanisms shaping the tumor immune landscape
.
Immunity
2018
;
48
:
399
416
.
10.
Li
J
,
Yuan
S
,
Norgard
RJ
,
Yan
F
,
Yamazoe
T
,
Blanco
A
, et al
Tumor cell–intrinsic USP22 suppresses antitumor immunity in pancreatic cancer
.
Cancer Immunol Res
2019
;
8
:
282
91
.
11.
Markosyan
N
,
Li
J
,
Sun
YH
,
Richman
LP
,
Lin
JH
,
Yan
F
, et al
Tumor cell-intrinsic EPHA2 suppresses anti-tumor immunity by regulating PTGS2 (COX-2)
.
J Clin Invest
2019
;
129
:
3594
609
.
12.
Li
J
,
Byrne
KT
,
Yan
F
,
Yamazoe
T
,
Chen
Z
,
Baslan
T
, et al
Tumor cell-intrinsic factors underlie heterogeneity of immune cell infiltration and response to immunotherapy
.
Immunity
2018
;
49
:
178
93
.
13.
Peng
W
,
Chen
JQ
,
Liu
C
,
Malu
S
,
Creasy
C
,
Tetzlaff
MT
, et al
Loss of PTEN promotes resistance to T cell–mediated immunotherapy
.
Cancer Discov
2016
;
6
:
202
16
.
14.
Welte
T
,
Kim
IS
,
Tian
L
,
Gao
X
,
Wang
H
,
Li
J
, et al
Oncogenic mTOR signalling recruits myeloid-derived suppressor cells to promote tumour initiation
.
Nat Cell Biol
2016
;
18
:
632
44
.
15.
Spranger
S
,
Bao
R
,
Gajewski
TF
. 
Melanoma-intrinsic β-catenin signalling prevents anti-tumour immunity
.
Nature
2015
;
523
:
231
5
.
16.
Liao
W
,
Overman
MJ
,
Boutin
AT
,
Shang
X
,
Zhao
D
,
Dey
P
, et al
KRAS-IRF2 axis drives immune suppression and immune therapy resistance in colorectal cancer
.
Cancer Cell
2019
;
35
:
559
72
.
17.
Wang
G
,
Lu
X
,
Dey
P
,
Deng
P
,
Wu
CC
,
Jiang
S
, et al
Targeting YAP-dependent MDSC infiltration impairs tumor progression
.
Cancer Discov
2016
;
6
:
80
95
.
18.
Kortlever
RM
,
Sodir
NM
,
Wilson
CH
,
Burkhart
DL
,
Pellegrinet
L
,
Brown Swigart
L
, et al
Myc cooperates with ras by programming inflammation and immune suppression
.
Cell
2017
;
6
:
1301
15
.
19.
Sodir
NM
,
Kortlever
RM
,
Barthet
VJA
,
Campos
T
,
Pellegrinet
L
,
Kupczak
S
, et al
Myc instructs and maintains pancreatic adenocarcinoma phenotype
.
Cancer Discov
2020
;
10
:
588
607
.
20.
Ardito
CM
,
Grüner
BM
,
Takeuchi
KK
,
Lubeseder-Martellato
C
,
Teichmann
N
,
Mazur
PK
, et al
EGF receptor is required for KRAS-induced pancreatic tumorigenesis
.
Cancer Cell
2012
;
22
:
304
17
.
21.
Navas
C
,
Hernández-Porras
I
,
Schuhmacher
AJ
,
Sibilia
M
,
Guerra
C
,
Barbacid
M
. 
EGF receptor signaling is essential for K-ras oncogene-driven pancreatic ductal adenocarcinoma
.
Cancer Cell
2012
;
22
:
318
30
.
22.
Nedaeinia
R
,
Avan
A
,
Manian
M
,
Salehi
R
,
Ghayour-Mobarhan
M
. 
EGFR as a potential target for the treatment of pancreatic cancer: dilemma and controversies
.
Curr Drug Targets
2014
;
15
:
1293
301
.
23.
Ng
SSW
,
Tsao
MS
,
Nicklee
T
,
Hedley
DW
. 
Effects of the epidermal growth factor receptor inhibitor OSI-774, Tarceva, on downstream signaling pathways and apoptosis in human pancreatic adenocarcinoma
.
Mol Cancer Ther
2002
;
1
:
777
83
.
24.
Moore
MJ
,
Goldstein
D
,
Hamm
J
,
Figer
A
,
Hecht
JR
,
Gallinger
S
, et al
Erlotinib plus gemcitabine compared with gemcitabine alone in patients with advanced pancreatic cancer: a phase III trial of the National Cancer Institute of Canada Clinical Trials Group
.
J Clin Oncol
2007
;
25
:
1960
6
.
25.
Rhim
AD
,
Mirek
ET
,
Aiello
NM
,
Maitra
A
,
Bailey
JM
,
McAllister
F
, et al
EMT and dissemination precede pancreatic tumor formation
.
Cell
2012
;
148
:
349
61
.
26.
Burr
ML
,
Sparbier
CE
,
Chan
KL
,
Chan
YC
,
Kersbergen
A
,
Lam
EYN
, et al
An evolutionarily conserved function of polycomb silences the MHC class I antigen presentation pathway and enables immune evasion in cancer
.
Cancer Cell
2019
;
36
:
385
401
.
27.
Peng
D
,
Kryczek
I
,
Nagarsheth
N
,
Zhao
L
,
Wei
S
,
Wang
W
, et al
Epigenetic silencing of TH1-type chemokines shapes tumour immunity and immunotherapy
.
Nature
2015
;
527
:
249
53
.
28.
Sheng
W
,
LaFleur
MW
,
Nguyen
TH
,
Chen
S
,
Chakravarthy
A
,
Conway
JR
, et al
LSD1 ablation stimulates anti-tumor immunity and enables checkpoint blockade
.
Cell
2018
;
174
:
549
63
.
29.
Pan
D
,
Kobayashi
A
,
Jiang
P
,
De Andrade
LF
,
Tay
RE
,
Luoma
AM
, et al
A major chromatin regulator determines resistance of tumor cells to T cell-mediated killing
.
Science
2018
;
359
:
770
5
.
30.
Li
F
,
Huang
Q
,
Luster
TA
,
Hu
H
,
Zhang
H
,
Ng
WL
, et al
In vivo epigenetic CRISPR screen identifies asf1a as an immunotherapeutic target in kras-mutant lung adenocarcinoma
.
Cancer Discov
2020
;
10
:
270
87
.
31.
Winograd
R
,
Byrne
KT
,
Evans
RA
,
Odorizzi
PM
,
Meyer
ARL
,
Bajor
DL
, et al
Induction of T-cell immunity overcomes complete resistance to PD-1 and CTLA-4 blockade and improves survival in pancreatic carcinoma
.
Cancer Immunol Res
2015
;
3
:
399
411
.
32.
Byrne
KT
,
Vonderheide
RH
. 
CD40 stimulation obviates innate sensors and drives T Cell Immunity in Cancer
.
Cell Rep
2016
;
15
:
2719
32
.
33.
O'Hara
MH
,
O'Reilly
EM
,
Rosemarie
M
,
Varadhachary
G
,
Wainberg
ZA
,
Ko
A
, et al
Abstract CT004: a phase Ib study of CD40 agonistic monoclonal antibody APX005M together with gemcitabine (Gem) and nab-paclitaxel (NP) with or without nivolumab (Nivo) in untreated metastatic ductal pancreatic adenocarcinoma (PDAC) patients
.
Cancer Res
2019
;
79
:
CT004
.
34.
Shi
Z
,
Xu
S
,
Xing
S
,
Yao
K
,
Zhang
L
,
Xue
L
, et al
Mettl17, a regulator of mitochondrial ribosomal RNA modifications, is required for the translation of mitochondrial coding genes
.
FASEB J
2019
;
33
:
13040
50
.
35.
Goel
S
,
Decristo
MJ
,
Watt
AC
,
Brinjones
H
,
Sceneay
J
,
Li
BB
, et al
CDK4/6 inhibition triggers anti-tumour immunity
.
Nature
2017
;
548
:
471
5
.
36.
Dandawate
P
,
Ghosh
C
,
Palaniyandi
K
,
Paul
S
,
Rawal
S
,
Pradhan
R
, et al
The histone demethylase KDM3A, increased in human pancreatic tumors, regulates expression of DCLK1 and promotes tumorigenesis in mice
.
Gastroenterology
2019
;
157
:
1646
59
.
37.
Elyada
E
,
Bolisetty
M
,
Laise
P
,
Flynn
WF
,
Courtois
ET
,
Burkhart
RA
, et al
Cross-species single-cell analysis of pancreatic ductal adenocarcinoma reveals antigen-presenting cancer-associated fibroblasts
.
Cancer Discov
2019
;
9
:
1102
23
.
38.
Chakravarthy
A
,
Khan
L
,
Bensler
NP
,
Bose
P
,
De Carvalho
DD
. 
TGF-β-associated extracellular matrix genes link cancer-associated fibroblasts to immune evasion and immunotherapy failure
.
Nat Commun
2018
;
9
:
4692
.
39.
Subramanian
A
,
Tamayo
P
,
Mootha
VK
,
Mukherjee
S
,
Ebert
BL
,
Gillette
MA
, et al
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
.
Proc Natl Acad Sci U S A
2005
;
102
:
15545
50
.
40.
Yamane
K
,
Toumazou
C
,
Tsukada
Y ichi
,
Erdjument-Bromage
H
,
Tempst
P
,
Wong
J
, et al
JHDM2A, a JmjC-containing H3K9 demethylase, facilitates transcription activation by androgen receptor
.
Cell
2006
;
125
:
483
95
.
41.
Goda
S
,
Isagawa
T
,
Chikaoka
Y
,
Kawamura
T
,
Aburatani
H
. 
Control of histone H3 lysine 9 (H3K9) methylation state via cooperative two-step demethylation by jumonji domain containing 1A (JMJD1A) homodimer
.
J Biol Chem
2013
;
288
:
36948
56
.
42.
Tateishi
K
,
Okada
Y
,
Kallin
EM
,
Zhang
Y
. 
Role of Jhdm2a in regulating metabolic gene expression and obesity resistance
.
Nature
2009
;
458
:
757
61
.
43.
Kuroki
S
,
Matoba
S
,
Akiyoshi
M
,
Matsumura
Y
,
Miyachi
H
,
Mise
N
, et al
Epigenetic regulation of mouse sex determination by the histone demethylase Jmjd1a
.
Science
2013
;
341
:
1106
9
.
44.
Li
J
,
Yu
B
,
Deng
P
,
Cheng
Y
,
Yu
Y
,
Kevork
K
, et al
KDM3 epigenetically controls tumorigenic potentials of human colorectal cancer stem cells through Wnt/β-catenin signalling
.
Nat Commun
2017
;
8
:
15146
.
45.
Cloos
PAC
,
Christensen
J
,
Agger
K
,
Helin
K
. 
Erasing the methyl mark: histone demethylases at the center of cellular differentiation and disease
.
Genes Dev
2008
;
22
:
1115
40
.
46.
Hublitz
P
,
Albert
M
,
Peters
AHFM
. 
Mechanisms of transcriptional repression by histone lysine methylation
.
Int J Dev Biol
2009
;
53
:
335
54
.
47.
Peters
AHFM
,
Kubicek
S
,
Mechtler
K
,
O'Sullivan
RJ
,
Derijck
AAHA
,
Perez-Burgos
L
, et al
Partitioning and plasticity of repressive histone methylation states in mammalian chromatin
.
Mol Cell
2003
;
12
:
1577
89
.
48.
Rice
JC
,
Briggs
SD
,
Ueberheide
B
,
Barber
CM
,
Shabanowitz
J
,
Hunt
DF
, et al
Histone methyltransferases direct different degrees of methylation to define distinct chromatin domains
.
Mol Cell
2003
;
12
:
1591
8
.
49.
Tachibana
M
,
Sugimoto
K
,
Fukushima
T
,
Shinkai
Y
. 
SET domain-containing protein, G9a, is a novel lysine-preferring mammalian histone methyltransferase with hyperactivity and specific selectivity to lysines 9 and 27 of histone H3
.
J Biol Chem
2001
;
276
:
25309
17
.
50.
Tachibana
M
,
Sugimoto
K
,
Nozaki
M
,
Ueda
J
,
Ohta
T
,
Ohki
M
, et al
G9a histone methyltransferase plays a dominant role in euchromatic histone H3 lysine 9 methylation and is essential for early embryogenesis
.
Genes Dev
2002
;
16
:
1779
91
.
51.
Wilson
S
,
Fan
L
,
Sahgal
N
,
Qi
J
,
Filipp F
V
. 
The histone demethylase KDM3A regulates the transcriptional program of the androgen receptor in prostate cancer cells
.
Oncotarget
2017
;
8
:
30328
43
.
52.
Barski
A
,
Cuddapah
S
,
Cui
K
,
Roh
TY
,
Schones
DE
,
Wang
Z
, et al
High-resolution profiling of histone methylations in the human genome
.
Cell
2007
;
129
:
823
37
.
53.
Filion
GJ
,
Van Steensel
B
. 
Reassessing the abundance of H3K9me2 chromatin domains in embryonic stem cells
.
Nat Genet
2010
;
42
:
4
.
54.
Wen
B
,
Wu
H
,
Shinkai
Y
,
Irizarry
RA
,
Feinberg
AP
. 
Large histone H3 lysine 9 dimethylated chromatin blocks distinguish differentiated from embryonic stem cells
.
Nat Genet
2009
;
41
:
246
50
.
55.
McConnell
BB
,
Yang
VW
. 
Mammalian Krüppel-Like factors in health and diseases
.
Physiol Rev
2010
;
90
:
1337
81
.
56.
Liang
G
,
Lin
JCY
,
Wei
V
,
Yoo
C
,
Cheng
JC
,
Nguyen
CT
, et al
Distinct localization of histone H3 acetylation and H3-K4 methylation to the transcription start sites in the human genome
.
Proc Natl Acad Sci U S A
2004
;
101
:
7357
62
.
57.
Bernstein
BE
,
Kamal
M
,
Lindblad-Toh
K
,
Bekiranov
S
,
Bailey
DK
,
Huebert
DJ
, et al
Genomic maps and comparative analysis of histone modifications in human and mouse
.
Cell
2005
;
120
:
169
81
.
58.
Tauriello
DVF
,
Palomo-Ponce
S
,
Stork
D
,
Berenguer-Llergo
A
,
Badia-Ramentol
J
,
Iglesias
M
, et al
TGFβ drives immune evasion in genetically reconstituted colon cancer metastasis
.
Nature
2018
;
554
:
538
43
.
59.
Mariathasan
S
,
Turley
SJ
,
Nickles
D
,
Castiglioni
A
,
Yuen
K
,
Wang
Y
, et al
TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells
.
Nature
2018
;
554
:
544
8
.
60.
Bayne
LJ
,
Beatty
GL
,
Jhala
N
,
Clark
CE
,
Rhim
AD
,
Stanger
BZ
, et al
Tumor-derived granulocyte-macrophage colony-stimulating factor regulates myeloid inflammation and T cell immunity in pancreatic cancer
.
Cancer Cell
2012
;
21
:
822
35
.
61.
Pylayeva-Gupta
Y
,
Lee
KE
,
Hajdu
CH
,
Miller
G
,
Bar-Sagi
D
. 
Oncogenic kras-induced GM-CSF production promotes the development of pancreatic neoplasia
.
Cancer Cell
2012
;
21
:
836
47
.
62.
David
CJ
,
Huang
Y-H
,
Chen
M
,
Su
J
,
Zou
Y
,
Bardeesy
N
, et al
TGF-β tumor suppression through a lethal EMT
.
Cell
2016
;
164
:
1015
30
.
63.
Meyers
RM
,
Bryan
JG
,
McFarland
JM
,
Weir
BA
,
Sizemore
AE
,
Xu
H
, et al
Computational correction of copy number effect improves specificity of CRISPR-Cas9 essentiality screens in cancer cells
.
Nat Genet
2017
;
49
:
1779
84
.
64.
Zhang
Y
,
Velez-Delgado
A
,
Mathew
E
,
Li
D
,
Mendez
FM
,
Flannagan
K
, et al
Myeloid cells are required for PD-1/PD-L1 checkpoint activation and the establishment of an immunosuppressive environment in pancreatic cancer
.
Gut
2017
;
66
:
124
36
.
65.
Bruns
CJ
,
Solorzano
CC
,
Harbison
MT
,
Ozawa
S
,
Tsan
R
,
Fan
D
, et al
Blockade of the epidermal growth factor receptor signaling by a novel tyrosine kinase inhibitor leads to apoptosis of endothelial cells and therapy of human pancreatic carcinoma
.
Cancer Res
2000
;
60
:
2926
35
.
66.
Feng
PH
,
Yu
CT
,
Chen
KY
,
Luo
CS
,
Wu
SM
,
Liu
CY
, et al
S100A9+ MDSC and TAM-mediated EGFR-TKI resistance in lung adenocarcinoma: the role of RELB
.
Oncotarget
2018
;
9
:
7631
43
.
67.
Runa
F
,
Hamalian
S
,
Meade
K
,
Shisgal
P
,
Gray
PC
,
Kelber
JA
. 
Tumor microenvironment heterogeneity: challenges and opportunities
.
Curr Mol Biol Reports
2017
;
3
:
218
29
.
68.
Morrison
AH
,
Byrne
KT
,
Vonderheide
RH
. 
Immunotherapy and prevention of pancreatic cancer
.
Trends Cancer
2018
;
4
:
418
28
.
69.
Manguso
RT
,
Pope
HW
,
Zimmer
MD
,
Brown
FD
,
Yates
KB
,
Miller
BC
, et al
In vivo CRISPR screening identifies Ptpn2 as a cancer immunotherapy target
.
Nature
2017
;
547
:
413
8
.
70.
Balli
D
,
Rech
AJ
,
Stanger
BZ
,
Vonderheide
RH
. 
Immune cytolytic activity stratifies molecular subsets of human pancreatic cancer
.
Clin Cancer Res
2017
;
23
:
3129
38
.
71.
Spranger
S
,
Luke
JJ
,
Bao
R
,
Zha
Y
,
Hernandez
KM
,
Li
Y
, et al
Density of immunogenic antigens does not explain the presence or absence of the T-cell-inflamed tumor microenvironment in melanoma
.
Proc Natl Acad Sci U S A
2016
;
113
:
E7759
68
.
72.
Sugiyama
E
,
Togashi
Y
,
Takeuchi
Y
,
Shinya
S
,
Tada
Y
,
Kataoka
K
, et al
Blockade of EGFR improves responsiveness to PD-1 blockade in EGFR-mutated non-small cell lung cancer
.
Sci Immunol
2020
;
5
:
eaav3937
.
73.
Zelenay
S
,
van der Veen
AG
,
Böttcher
JP
,
Snelgrove
KJ
,
Rogers
N
,
Acton
SE
, et al
Cyclooxygenase-dependent tumor growth through evasion of immunity
.
Cell
2015
;
162
:
1257
70
.
74.
Gong
K
,
Guo
G
,
Panchani
N
,
Bender
ME
,
Gerber
DE
,
Minna
JD
, et al
EGFR inhibition triggers an adaptive response by co-opting antiviral signaling pathways in lung cancer
.
Nat Cancer
2020
;
1
:
394
409
.
75.
Jänne
PA
,
Chih-Hsin Yang
J
,
Kim
DW
,
Planchard
D
,
Ohe
Y
,
Ramalingam
SS
, et al
AZD9291 in EGFR inhibitor-resistant non-small-cell lung cancer
.
N Engl J Med
2015
;
372
:
1689
99
.
76.
Rosell
R
,
Carcereny
E
,
Gervais
R
,
Vergnenegre
A
,
Massuti
B
,
Felip
E
, et al
Erlotinib versus standard chemotherapy as first-line treatment for European patients with advanced EGFR mutation-positive non-small-cell lung cancer (EURTAC): a multicentre, open-label, randomised phase 3 trial
.
Lancet Oncol
2012
;
12
:
735
42
.
77.
Sequist
LV
,
Soria
JC
,
Goldman
JW
,
Wakelee
HA
,
Gadgeel
SM
,
Varga
A
, et al
Rociletinib in EGFR-mutated non-small-cell lung cancer
.
N Engl J Med
2015
;
372
:
1700
9
.
78.
Kim
IS
,
Gao
Y
,
Welte
T
,
Wang
H
,
Liu
J
,
Janghorban
M
, et al
Immuno-subtyping of breast cancer reveals distinct myeloid cell profiles and immunotherapy resistance mechanisms
.
Nat Cell Biol
2019
;
21
:
1113
26
.
79.
Tavazoie
MF
,
Pollack
I
,
Tanqueco
R
,
Ostendorf
BN
,
Reis
BS
,
Gonsalves
FC
, et al
LXR/ApoE activation restricts innate immune suppression in cancer
.
Cell
2018
;
172
:
825
40
.
80.
Shalem
O
,
Sanjana
NE
,
Hartenian
E
,
Shi
X
,
Scott
DA
,
Mikkelsen
TS
, et al
Genome-scale CRISPR-Cas9 knockout screening in human cells
.
Science
2014
;
343
:
84
7
.
81.
Li
W
,
Xu
H
,
Xiao
T
,
Cong
L
,
Love
MI
,
Zhang
F
, et al
MAGeCK enables robust identification of essential genes from genome-scale CRISPR/Cas9 knockout screens
.
Genome Biol
2014
;
15
:
554
.
82.
Chen
EY
,
Tan
CM
,
Kou
Y
,
Duan
Q
,
Wang
Z
,
Meirelles
G V.
, et al
Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool
.
BMC Bioinformatics
2013
;
14
:
128
.
83.
Kuleshov
MV
,
Jones
MR
,
Rouillard
AD
,
Fernandez
NF
,
Duan
Q
,
Wang
Z
, et al
Enrichr: a comprehensive gene set enrichment analysis web server 2016 update
.
Nucleic Acids Res
2016
;
44
:
W90
W97
.
84.
Buenrostro
JD
,
Giresi
PG
,
Zaba
LC
,
Chang
HY
,
Greenleaf
WJ
. 
Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position
.
Nat Methods
2013
;
10
:
1213
8
.
85.
Zhang
Y
,
Liu
T
,
Meyer
CA
,
Eeckhoute
J
,
Johnson
DS
,
Bernstein
BE
, et al
Model-based analysis of ChIP-Seq (MACS)
.
Genome Biol
2008
;
9
:
R137
.
86.
Quinlan
AR
,
Hall
IM
. 
BEDTools: a flexible suite of utilities for comparing genomic features
.
Bioinformatics
2010
;
26
:
841
2
.
87.
Robinson
MD
,
McCarthy
DJ
,
Smyth
GK
. 
edgeR: a bioconductor package for differential expression analysis of digital gene expression data
.
Bioinformatics
2009
;
26
:
139
40
.
88.
McCarthy
DJ
,
Chen
Y
,
Smyth
GK
. 
Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation
.
Nucleic Acids Res
2012
;
40
:
4288
97
.
89.
McLean
CY
,
Bristor
D
,
Hiller
M
,
Clarke
SL
,
Schaar
BT
,
Lowe
CB
, et al
GREAT improves functional interpretation of cis-regulatory regions
.
Nat Biotechnol
2010
;
28
:
495
501
.
90.
Langmead
B
,
Salzberg
SL
. 
Fast gapped-read alignment with bowtie 2
.
Nat Methods
2012
;
9
:
357
9
.
91.
Li
H
,
Handsaker
B
,
Wysoker
A
,
Fennell
T
,
Ruan
J
,
Homer
N
, et al
The sequence alignment/map format and SAMtools
.
Bioinformatics
2009
;
25
:
2078
9
.
92.
Yu
G
,
Wang
LG
,
He
QY
. 
ChIP seeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization
.
Bioinformatics
2015
;
31
:
2382
3
.
93.
Bailey
TL
,
Boden
M
,
Buske
FA
,
Frith
M
,
Grant
CE
,
Clementi
L
, et al
MEME Suite: tools for motif discovery and searching
.
Nucleic Acids Res
2009
;
37
:
W202
W208
.
94.
Bailey
TL
. 
DREME: motif discovery in transcription factor ChIP-seq data
.
Bioinformatics
2011
;
27
:
1653
9
.
95.
Gupta
S
,
Stamatoyannopoulos
JA
,
Bailey
TL
,
Noble
WS
. 
Quantifying similarity between motifs
.
Genome Biol
2007
;
8
:
R24
.
96.
Sidoli
S
,
Bhanu
NV
,
Karch
KR
,
Wang
X
,
Garcia
BA
. 
Complete workflow for analysis of histone post-translational modifications using bottom-up mass spectrometry: from histone extraction to data analysis
.
J Vis Exp
2016
;
111
:54112
.
97.
Seo
JH
,
Agarwal
E
,
Bryant
KG
,
Caino
MC
,
Kim
ET
,
Kossenkov
A V
, et al
Syntaphilin ubiquitination regulates mitochondrial dynamics and tumor cell movements
.
Cancer Res
2018
;
78
:
4215
28
.
98.
Cox
J
,
Mann
M
. 
MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification
.
Nat Biotechnol
. 
2008
;
26
:
1367
72
.