Abstract
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.
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
Introduction
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.
Results
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).
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. 1E–G; 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. 2A–C). 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.
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).
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. 4B–D; 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.
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. 5A–C; 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).
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.
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.
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.
Discussion
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).
Methods
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+CD11b−F4/80−CD11c+ CD103+; T cells: live CD45+CD11b− F4/80−NKp46−CD3+; CD4+ T cells: live CD45+CD11b−F4/80−NKp46−CD3+CD4+; CD8+ T cells: live CD45+CD11b−F4/80−NKp46−CD3+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.
Authors' Disclosures
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.
Authors' Contributions
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.
Acknowledgments
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.