Abstract
Immune checkpoint blockade (ICB) therapy revolutionized cancer treatment, but many patients with impaired MHC-I expression remain refractory. Here, we combined FACS-based genome-wide CRISPR screens with a data-mining approach to identify drugs that can upregulate MHC-I without inducing PD-L1. CRISPR screening identified TRAF3, a suppressor of the NFκB pathway, as a negative regulator of MHC-I but not PD-L1. The Traf3-knockout gene expression signature is associated with better survival in ICB-naïve patients with cancer and better ICB response. We then screened for drugs with similar transcriptional effects as this signature and identified Second Mitochondria-derived Activator of Caspase (SMAC) mimetics. We experimentally validated that the SMAC mimetic birinapant upregulates MHC-I, sensitizes cancer cells to T cell–dependent killing, and adds to ICB efficacy. Our findings provide preclinical rationale for treating tumors expressing low MHC-I expression with SMAC mimetics to enhance sensitivity to immunotherapy. The approach used in this study can be generalized to identify other drugs that enhance immunotherapy efficacy.
MHC-I loss or downregulation in cancer cells is a major mechanism of resistance to T cell–based immunotherapies. Our study reveals that birinapant may be used for patients with low baseline MHC-I to enhance ICB response. This represents promising immunotherapy opportunities given the biosafety profile of birinapant from multiple clinical trials.
This article is highlighted in the In This Issue feature, p. 1307
Introduction
Immune checkpoint blockade (ICB), primarily anti–PD-1, anti–PD-L1, and anti-CTLA4, has been shown to induce remarkable response across different cancer types in many clinical studies (1–8). However, only a minority of patients respond to ICB, and resistance can develop in patients who initially respond (9). The mechanisms of resistance have been under extensive investigation. Multiple cancer cell–intrinsic mechanisms have been identified, including low neoantigen load (2, 3, 10, 11), loss of antigen presentation (12, 13), loss of HLA heterozygosity (14, 15), impaired response or prolonged exposure to IFNγ (16–19), and activation of certain oncogenic pathways/signatures (20–22), as well as other mechanisms (22, 23). Among these ICB resistance mechanisms, insufficient antigen presentation to activate T cells comprises a large proportion of patients (24–26). MHC class I (MHC-I), the key component in antigen presentation, presents intracellular peptide antigens to the cell surface for recognition by the specific CD8+ T cells (27). Loss of MHC-I by genetic mutation or epigenetic silencing impairs T-cell recognition and activation, and compromises antitumor activity. Furthermore, ICB treatment of CD8+ T cells in the absence of antigen stimulation can induce, rather than alleviate, T-cell dysfunction (28). Therefore, enhancing MHC-I levels in cancer cells, especially in cancers with low baseline MHC-I levels, is a promising strategy to improve ICB efficacy.
Genes encoding the MHC-I components (for example, HLA-A/B/C and B2M) are expressed widely yet are under tight regulation (29). In tumors, IFNγ secreted by tumor-infiltrating T cells can induce the expression of the MHC-I component genes through the JAK–STAT pathway (30) and boost the CD8+ T cell–dependent cancer cell elimination (31). However, recent reports have shown that IFNγ can also exert immunosuppressive functions, partially through its induction of immune checkpoint genes such as CD274 (encoding PD-L1) to inhibit T-cell activation and promote T-cell dysfunction (19, 32, 33). In addition to IFNγ, multiple other cytokines were reported to upregulate both MHC-I and PD-L1 with opposite effects on tumor immunity (31). Therefore, we sought to identify gene perturbations and drugs that specifically induce MHC-I expression without increasing PD-L1 expression.
CRISPR screening is a powerful technology with increasing popularity for target identification in cancer and immunology, and FACS-based CRISPR screens have previously identified regulators of MHC-I or PD-L1 (34–36). However, it is difficult to identify factors that regulate MHC-I but not PD-L1 based on existing data, due to variations in sorting strategies and limitations in statistical power. Likewise, several exciting studies have reported chemotherapy or targeted therapy with immunomodulatory effects by regulating MHC-I and/or PD-L1 (37–40), but a systematic approach to identifying drugs that preferentially modulate MHC-I has proved elusive.
In this study, we used a CRISPR screening approach with dual-marker FACS sorting to identify factors that decouple the regulation of MHC-I and PD-L1. The experimentally validated target was used to generate a knockout (KO) differential expression signature. Using this signature, we analyzed transcriptome data from drug perturbation studies to identify drugs that regulate MHC-I but not PD-L1. Finally, we validated the effect of the identified drug to enhance ICB response in a T cell–dependent manner in vivo.
Results
An Integrated Workflow to Identify Drugs with MHC-I–Modulatory Effects
Analysis of bulk tumor gene expression data from The Cancer Genome Atlas (TCGA) revealed a strong positive correlation between genes encoding MHC-I and PD-L1 in multiple cancer types (Supplementary Fig. S1A). Furthermore, analysis of Cancer Cell Line Encyclopedia (CCLE) cell line gene expression data in vitro revealed a similar trend (Supplementary Fig. S1B). Despite the overall positive correlation, the ratio of MHC-I over PD-L1 expression is positively correlated with patient overall survival in the majority of cancer types, although TCGA tumors were not treated with immunotherapy (Supplementary Fig. S1C).
To search for drugs that can specifically upregulate MHC-I without inducing PD-L1, we designed an integrated experimental and computational workflow, which consists of five steps (Supplementary Fig. S1D). First, we performed genome-wide CRISPR screens with dual-marker sorting to identify candidate genes whose KO can upregulate MHC-I without upregulating PD-L1. Second, we performed experiments on the identified gene candidates to validate the MHC-I–specific effect and test the sensitivity to ICB therapy. Third, we characterized the target KO gene expression signature in vitro. Next, we searched the Gene Expression Omnibus (GEO) to identify drugs with differential gene expression signature similar to that from the candidate gene KO. Finally, we validated the drug for its immunomodulatory effects in vitro and in vivo.
CRISPR Screen Identifies TRAF3 as a Negative Regulator of MHC-I
To systematically identify regulators of MHC-I and PD-L1, we performed FACS-based genome-wide CRISPR screens using the mouse melanoma cell line B16F10 (Fig. 1A), which can be induced by IFNγ to express MHC-I and PD-L1. Low-dose (0.1 ng/mL) IFNγ showed minimal induction of MHC-I and PD-L1 and enabled us to identify negative regulators whose deficiency can increase MHC-I and/or PD-L1 protein level (Fig. 1B). In comparison, high-dose (10 ng/mL) IFNγ strongly induced MHC-I and PD-L1 and enabled us to identify the positive regulators of MHC-I and/or PD-L1 protein (Fig. 1B). We transduced the B16F10 cells with our mouse genome-wide CRISPR KO library, expanded the transduced cells, and treated the cells with either a low or high dose of IFNγ for 2 days. We then performed FACS to isolate the MHC-IhiPD-L1hi, MHC-IhiPD-L1lo, MHC-IloPD-L1hi, and MHC-IloPD-L1lo subpopulations and compared the single guide RNA (sgRNA) frequencies in each subpopulation with the unsorted population (Fig. 1A). There was a large decrease in the number of detected sgRNAs in the regulator-enriched quadrants compared with the unsorted population, indicating strong selections in these subpopulations (Supplementary Fig. S2A).
CRISPR screens identify novel regulators of MHC-I. A, Workflow of using CRISPR screens to identify the positive or negative regulators of MHC-I and/or PD-L1. We transduced B16F10 mouse melanoma cells with an in-house–designed genome-wide sgRNA library, expanded the transduced cells, and stimulated the cells with 0.1 ng/mL (low dose) or 10 ng/mL IFNγ (high dose) for different levels of MHC-I/PD-L1 induction. We then performed FACS to isolate the MHC-IhiPD-L1hi, MHC-IhiPD-L1lo, MHC-IloPD-L1hi, and MHC-IloPD-L1lo subpopulations. We amplified and sequenced the gRNAs in these subpopulations as well as the bulk presorting population to identify genes that were enriched in each sorted subpopulation. MHC-I and/or PD-L1 regulators are expected to be enriched in the subpopulations as indicated. B, Titration of IFNγ concentration to test its effect on MHC-I and PD-L1 expression. Histogram of H2-Kb and PD-L1 levels assessed by flow cytometry following 2-day treatment with different concentrations of IFNγ. C, CRISPR screen reveals known and novel candidate regulators of MHC-I and/or PD-L1. Ranked dot plots of gene enrichment in each sorted subpopulation compared with the unsorted population are shown. The X axis shows the rank of each gene, and the Y axis shows the log2 enrichment of sgRNAs for each gene in the indicated subpopulation compared with the unsorted population. D, Immunoblot of B16F10 cells with indicated genotypes shows good efficiency of Traf3-KO or overexpression. E and F, Validation of TRAF3 as a negative regulator of MHC-I with IFNγ induction. B16F10 cells transduced with control sgRNA or sgTraf3 were cultured for 48 hours, with IFNγ concentrations as indicated, and then assessed on their MHC-I and PD-L1 levels. E, Typical histogram of H2-Kb and PD-L1 FACS plot of control or Traf3-deficient B16F10 cells in each treatment condition. F, Quantification of median fluorescence intensity (MFI) of H2-Kb or PD-L1 from E. Values are normalized to the sgControl group with vehicle treatment (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). G and H, Validation of TRAF3 as a negative regulator of MHC-I with type-I IFN induction. B16F10 cells transduced with control sgRNA or sgTraf3 were cultured for 48 hours with treatment as indicated, and then assessed on their MHC-I and PD-L1 levels. 500 U/mL IFNα or IFNβ was used. G, Typical histogram of H2-Kb and PD-L1 FACS plot of control or Traf3-deficient B16F10 cells in each treatment condition. H, Quantification of MFI of H2-Kb or PD-L1 from G. Values are normalized to the sgControl group with vehicle treatment (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition).
CRISPR screens identify novel regulators of MHC-I. A, Workflow of using CRISPR screens to identify the positive or negative regulators of MHC-I and/or PD-L1. We transduced B16F10 mouse melanoma cells with an in-house–designed genome-wide sgRNA library, expanded the transduced cells, and stimulated the cells with 0.1 ng/mL (low dose) or 10 ng/mL IFNγ (high dose) for different levels of MHC-I/PD-L1 induction. We then performed FACS to isolate the MHC-IhiPD-L1hi, MHC-IhiPD-L1lo, MHC-IloPD-L1hi, and MHC-IloPD-L1lo subpopulations. We amplified and sequenced the gRNAs in these subpopulations as well as the bulk presorting population to identify genes that were enriched in each sorted subpopulation. MHC-I and/or PD-L1 regulators are expected to be enriched in the subpopulations as indicated. B, Titration of IFNγ concentration to test its effect on MHC-I and PD-L1 expression. Histogram of H2-Kb and PD-L1 levels assessed by flow cytometry following 2-day treatment with different concentrations of IFNγ. C, CRISPR screen reveals known and novel candidate regulators of MHC-I and/or PD-L1. Ranked dot plots of gene enrichment in each sorted subpopulation compared with the unsorted population are shown. The X axis shows the rank of each gene, and the Y axis shows the log2 enrichment of sgRNAs for each gene in the indicated subpopulation compared with the unsorted population. D, Immunoblot of B16F10 cells with indicated genotypes shows good efficiency of Traf3-KO or overexpression. E and F, Validation of TRAF3 as a negative regulator of MHC-I with IFNγ induction. B16F10 cells transduced with control sgRNA or sgTraf3 were cultured for 48 hours, with IFNγ concentrations as indicated, and then assessed on their MHC-I and PD-L1 levels. E, Typical histogram of H2-Kb and PD-L1 FACS plot of control or Traf3-deficient B16F10 cells in each treatment condition. F, Quantification of median fluorescence intensity (MFI) of H2-Kb or PD-L1 from E. Values are normalized to the sgControl group with vehicle treatment (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). G and H, Validation of TRAF3 as a negative regulator of MHC-I with type-I IFN induction. B16F10 cells transduced with control sgRNA or sgTraf3 were cultured for 48 hours with treatment as indicated, and then assessed on their MHC-I and PD-L1 levels. 500 U/mL IFNα or IFNβ was used. G, Typical histogram of H2-Kb and PD-L1 FACS plot of control or Traf3-deficient B16F10 cells in each treatment condition. H, Quantification of MFI of H2-Kb or PD-L1 from G. Values are normalized to the sgControl group with vehicle treatment (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition).
We used the MAGeCK (refs. 41, 42; Fig. 1C; Supplementary Fig. S2B; Supplementary Table S1) and RIGER (ref. 43; Supplementary Fig. S2C) computational algorithms to analyze the CRISPR screen data. Comparison of the MHC-IloPD-L1lo subpopulation with the unsorted population in the 10 ng/mL IFNγ cohort identified known positive regulators of both MHC-I and PD-L1, including Jak/Stat and Ifngr genes (Fig. 1C). Comparison of the MHC-IhiPD-L1lo subpopulation with the unsorted population in the 10 ng/mL IFNγ cohort identified Cd274 (encoding PD-L1), as expected, as the top PD-L1–specific positive regulator (Supplementary Fig. S2B). Comparison of the MHC-IhiPD-L1hi subpopulation with the unsorted population in the 0.1 ng/mL IFNγ cohort identified novel candidate common negative regulators, including Stub1 and Ube2n (Fig. 1C). Stub1 was reported to downregulate PD-L1 (35), so its role in regulating MHC-I merits further study.
The comparisons most central to the goals of this study are those differentially regulating MHC-I and PD-L1 (MHC-IloPD-L1hi and MHC-IhiPD-L1lo). Comparison of the MHC-IloPD-L1hi subpopulation with the unsorted population in the 10 ng/mL IFNγ condition revealed genes known to be involved in MHC-I transcription and assembly, including B2m, Nlrc5, Tap1, Tap2, Tapbp, and H2-K1, which were ranked among the top MHC-I–specific positive regulators (Fig. 1C). Comparison of the MHC-IhiPD-L1lo subpopulation with the unsorted population in the 0.1 ng/mL IFNγ condition identified previously reported (Ezh2; refs. 36, 44) and novel candidate negative regulators of MHC-I, including Traf3, Tada3, and Med13 (Fig. 1C). TRAF3, the top-ranking putative MHC-I–specific negative regulator in our screens, is an adaptor protein involved in multiple signaling pathways leading to NFκB activation (45). NFκB is known to regulate MHC-I component genes by directly binding to their enhancers (46). A recent study in neuroblastoma also implicated TRAF3 as mediating the MHC-I–repressive effects of N4BP1 on MHC-I (47). These previous studies comport with our screen results that TRAF3 is a negative regulator of MHC-I.
We then validated the effect of TRAF3 perturbation on MHC-I and PD-L1 levels using CRISPR-mediated single gene KO of Traf3 with two different sgRNAs, and overexpression of Traf3 (Fig. 1D). As expected, Traf3-KO increased the expression of MHC-I, but not PD-L1, in response to various doses of IFNγ (Fig. 1E and F). Conversely, overexpression of Traf3 specifically decreased the expression of MHC-I, but not PD-L1, in response to IFNγ (Supplementary Fig. S2D). To test whether MHC-I is upregulated by cell-extrinsic autocrines, we compared the sgControl or sgTraf3 cells cultured in conditioned media from Traf3–wild-type (WT) or Traf3-KO cells. We did not detect a significant difference in MHC-I expression between different media conditions (Supplementary Fig. S2E). Furthermore, we compared the coculture of GFP+ WT B16F10 cells with GFP− Traf3-WT or Traf3-KO B16F10 cells, and observed no significant difference in MHC-I expression level in GFP+ WT cells between different coculture conditions either (Supplementary Fig. S2F). These data suggest that Traf3-KO specifically upregulates MHC-I in a cell-autonomous manner.
To test the generality of this regulation, we measured the effect of Traf3 deficiency in the context of induction by other cytokines besides IFNγ. We found that Traf3 deficiency also led to specific upregulation of MHC-I in response to TNFα and type I IFNs (Fig. 1G and H; Supplementary Fig. S2G and S2H). We also assessed the role of TRAF3 in MHC-I regulation in an independent mouse cell line (Supplementary Fig. S2I–S2M), the colorectal cancer line CT26. Despite higher baseline MHC-I expression in CT26, Traf3-KO also led to specific upregulation of MHC-I with or without cytokine treatments (Supplementary Fig. S2I–S2M). We further tested the function of TRAF3 in human cancer cells and found that its KO can preferentially upregulate MHC-I in multiple human lines (Supplementary Fig. S2N–S2Q), including SKMEL5 (melanoma), K028 (melanoma), HT29 (colorectal cancer), and MCF7 (breast cancer). These results support the general role of TRAF3 in suppressing MHC-I in human and mouse cancer cell lines.
TRAF3 Depletion Activates MHC-I through NFκB
To study the mechanism by which TRAF3 negatively regulates MHC-I levels, we performed RNA sequencing (RNA-seq) and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) on Traf3-WT or Traf3-KO B16F10 cells that were treated with vehicle or 1 ng/mL IFNγ. As expected, RNA-seq revealed multiple MHC-I components, such as H2-K1 and H2-D1, among the top upregulated genes in the Traf3-KO cells under baseline condition or after IFNγ induction (Fig. 2A; Supplementary Fig. S3A; Supplementary Table S2). Other immune-related genes were also upregulated, including the proinflammatory chemokines Cxcl10 and Ccl9 and components of the NFκB pathway (Fig. 2A; Supplementary Fig. S3A; Supplementary Table S2). By contrast, Traf3-KO did not alter the transcription of Cd274 (PD-L1; Supplementary Table S2). Gene Ontology (GO) analysis shows that genes upregulated upon Traf3-KO were enriched in antigen presentation, the NFκB pathway, and TNF-mediated signaling (Fig. 2B; Supplementary Fig. S3B). This is consistent with previous reports that TRAF3 is a negative regulator of NFκB in response to various extracellular signals (48–50).
TRAF3 deficiency upregulates MHC-I through NFκB. A, RNA-seq of B16F10 cells transduced with sgControl or sgTraf3 shows upregulation of MHC-I–related genes in the absence of Traf3. Heat map of differential expression of genes induced by TRAF3 deficiency with IFNγ treatment. B, Gene set enrichment analysis of upregulated pathways (GO biological pathway) in sgTraf3 cells compared with sgControl cells with IFNγ treatment. Multiple pathways, such as antigen presentation and NFκB signaling, were upregulated by the deletion of Traf3. C, ATAC-seq of Traf3-normal or -deficient B16F10 cells revealed that TRAF3 deficiency leads to higher chromatin accessibility near genes encoding components of the MHC-I complex. **, P < 0.01; ***, P < 0.001. D, Cistrome-GO analysis of the more accessible regions in sgTraf3 compared with sgControl cells with IFNγ treatment. E, Cistrome toolkit analysis of ATAC-seq data revealed that DNA-binding sites of RELA were more open in the Traf3-deficient cells. F, Enrichment of motifs in the accessible chromatin regions specific to Traf3-deficient cells. The top enriched motifs (RELA) are shown. G, Typical immunoblot of NFκB signaling components in Traf3-normal or Traf3-deficient B16F10 cells in response to IFNγ induction. B16F10 cells transduced with sgControl or sgTraf3 were induced by 1 ng/mL IFNγ for 0, 0.5, or 24 hours, and then harvested for immunoblot. Data from one typical experiment out of three biological replicates is shown. H, Quantification of immunoblot signals from G based on three biological replicates. TRAF3 deletion leads to upregulated NFκB signaling. I and J, B16F10 cells transduced with control sgRNA or sgTraf3 were treated with vehicle control, IFNγ (1 ng/mL), and/or TPCA-1 (1 μmol/L) for 48 hours, and then assessed on their MHC-I and PD-L1 levels. I, Typical histogram of H2-Kb and PD-L1 FACS plot of control- or Traf3-deficient B16F10 cells in each treatment condition. J, Quantification of MFI of H2-Kb or PD-L1 from I. Values are normalized to sgControl group with vehicle treatment (***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing IFNγ and IFNγ+TPCA groups).
TRAF3 deficiency upregulates MHC-I through NFκB. A, RNA-seq of B16F10 cells transduced with sgControl or sgTraf3 shows upregulation of MHC-I–related genes in the absence of Traf3. Heat map of differential expression of genes induced by TRAF3 deficiency with IFNγ treatment. B, Gene set enrichment analysis of upregulated pathways (GO biological pathway) in sgTraf3 cells compared with sgControl cells with IFNγ treatment. Multiple pathways, such as antigen presentation and NFκB signaling, were upregulated by the deletion of Traf3. C, ATAC-seq of Traf3-normal or -deficient B16F10 cells revealed that TRAF3 deficiency leads to higher chromatin accessibility near genes encoding components of the MHC-I complex. **, P < 0.01; ***, P < 0.001. D, Cistrome-GO analysis of the more accessible regions in sgTraf3 compared with sgControl cells with IFNγ treatment. E, Cistrome toolkit analysis of ATAC-seq data revealed that DNA-binding sites of RELA were more open in the Traf3-deficient cells. F, Enrichment of motifs in the accessible chromatin regions specific to Traf3-deficient cells. The top enriched motifs (RELA) are shown. G, Typical immunoblot of NFκB signaling components in Traf3-normal or Traf3-deficient B16F10 cells in response to IFNγ induction. B16F10 cells transduced with sgControl or sgTraf3 were induced by 1 ng/mL IFNγ for 0, 0.5, or 24 hours, and then harvested for immunoblot. Data from one typical experiment out of three biological replicates is shown. H, Quantification of immunoblot signals from G based on three biological replicates. TRAF3 deletion leads to upregulated NFκB signaling. I and J, B16F10 cells transduced with control sgRNA or sgTraf3 were treated with vehicle control, IFNγ (1 ng/mL), and/or TPCA-1 (1 μmol/L) for 48 hours, and then assessed on their MHC-I and PD-L1 levels. I, Typical histogram of H2-Kb and PD-L1 FACS plot of control- or Traf3-deficient B16F10 cells in each treatment condition. J, Quantification of MFI of H2-Kb or PD-L1 from I. Values are normalized to sgControl group with vehicle treatment (***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing IFNγ and IFNγ+TPCA groups).
Consistent with the RNA-seq data, ATAC-seq showed that DNA accessibility is higher in Traf3-KO cells in the promoter regions of MHC-I components (Fig. 2C; Supplementary Fig. S3C; Supplementary Table S2). We also performed Cistrome-GO analysis (51) of the ATAC-seq data, which ranks each gene's regulatory potential with nearby peaks' strength and distance, and then examines the enrichment of GO pathways on the top-ranking genes. Genes with increased regulatory potential in Traf3-KO cells showed significant enrichment of antigen processing and presentation, TNF signaling pathway, and NFκB signaling pathway (Fig. 2D; Supplementary Fig. S3D). In addition, HOMER motif enrichment analysis (52) found NFκB family motif to be the most enriched in the sites with increased chromatin accessibility upon Traf3 KO (Fig. 2E). Furthermore, we used the CistromeDB Toolkit function (53, 54) to identify transcription factors (TF) whose published binding sites have the highest overlap with our differential ATAC-seq peaks, quantified by GIGGLE score (55). This analysis found binding sites of RELA, a component of NFκB, to be the most enriched in the increased accessible sites in Traf3-KO cells (Fig. 2F; Supplementary Fig. S3E). Other TFs with enriched binding sites, such as AIRE, STAT3, IRF8, and STAT1 (Fig. 2F), are also known to be involved in immune response.
The NFκB pathway is involved in multiple biological processes, such as immune response, developmental processes, cell proliferation, and cell death (56–58). The NFκB family of TFs can be activated by “canonical” or “noncanonical” signaling pathways, which differ from each other in upstream regulators, TF components, and signaling kinetics (59). The canonical pathway predominantly involves RELA, cREL, and NFKB1 (p105/p50), whereas the noncanonical pathway mainly involves RELB and NFKB2 (p100/p52). Immunoblot of Traf3-WT or Traf3-KO cells with IFNγ treatment revealed that TRAF3 deficiency led to activation of both canonical and noncanonical NFκB pathways, with increased phospho-RELA, RELB, cREL, and p52 (the active form of NFKB2; Fig. 2G and H). Analysis of existing chromatin immunoprecipitation sequencing (ChIP-seq) datasets in Cistrome DB (53, 54) revealed that both canonical and noncanonical NFκB TFs can bind to the promoters/enhancers of MHC-I in human and mouse cells (refs. 60–62; Supplementary Fig. S3F). In addition, consistent with the reported role of TRAF3 as an adaptor for cIAP-mediated degradation of NIK (48, 49), an activator of the noncanonical NFκB pathway, TRAF3-deficient cells showed higher NIK levels (Supplementary Fig. S3G). To test whether the NFκB pathway underlies the upregulation of MHC-I in the absence of TRAF3, we treated Traf3-KO cells with NFκB inhibitors, TPCA-1 and IKK-16, and tested whether the upregulation of MHC-I in Traf3-KO cells could be reversed. Indeed, the ability of IFNγ to upregulate MHC-I was drastically attenuated when Traf3-KO cells were treated with TPCA-1 (Fig. 2I and J) or IKK-16 (Supplementary Fig. S3H and S3I). This finding is consistent with recent studies that activation of NFκB-related pathways can induce MHC-I expression in melanoma cells (63, 64). Taken together, these data suggest that TRAF3 deficiency leads to upregulation of MHC-I through elevated NFκB activity.
TRAF3 Deficiency Facilitates T Cell–Driven Killing of Cancer Cells
Because TRAF3 loss upregulated MHC-I but not PD-L1 in response to IFNγ, we hypothesized that it may also render cancer cells more susceptible to T cell–driven cytotoxicity. We tested this hypothesis by coculturing B16F10 cells with either Pmel-1 T cells, which have low-affinity T-cell receptors (TCR) for an endogenous melanoma antigen, or OT-I T cells, which have high affinity for an exogenously introduced ovalbumin antigen (Fig. 3A). Traf3-KO resulted in fewer viable cancer cells and higher MHC-I expression than Traf3-WT after coculture in both models (Fig. 3B–D; Supplementary Fig. S4A–S4E). In the coculture with Pmel-1 T cells, the killing efficiency in WT B16F10 cells is low at baseline due to the low expression of MHC-I. Therefore, it is common to pretreat cancer cells with IFNγ, which boosts MHC-I expression, prior to coculture with T cells in this model. We used this method and observed significantly more efficient killing of Traf3-KO cells (relative cell number in Fig. 3C; raw killing% in Supplementary Fig. S4C). Furthermore, we found that the Traf3-KO B16F10 cells could be efficiently killed by Pmel-1 T cells even without IFNγ pretreatment (relative cell number in Supplementary Fig. S4A; raw killing% in Supplementary Fig. S4D), and exhibited a substantial induction of MHC-I expression (Supplementary Fig. S4E). Conversely, overexpression of TRAF3 resulted in reduced sensitivity of cancer cells to T cell–driven cytotoxicity (relative cell number in Fig. 3E, raw killing% in Supplementary Fig. S4F).
TRAF3 deletion sensitizes cancer cells to T cell–driven cytotoxicity. A, Workflow of testing the role of TRAF3 in regulating the response of cancer cells to T cell–driven cytotoxicity through in vitro coculture. Traf3-normal or Traf3-deficient B16F10 cells were cultured with antigen-specific CD8+ T cells (Pmel-1 or OT-I) for 1 to 3 days. Cell number and cell surface MHC-I/PD-L1 expression were quantified by FACS. B and C, Relative B16F10 cell number after coculture with (B) OT-I T cells or (C) Pmel-1 T cells at different Effector:Target (E:T) ratios revealed a higher sensitivity of Traf3-deficient B16F10 cells to T-cell–mediated cytotoxicity. For coculture with Pmel-1 T cells, B16F10 cells were pretreated with 1 ng/mL IFNγ for 12 hours prior to the coculture. The bar plots present the relative cell number in each group, normalized to the cell number in the sgControl group in each E:T condition. Mean ± SD and individual replicate values are shown for each group (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). D, Relative MFI of H2-Kb or PD-L1 of B16F10 cells cocultured with OT-I T cells. Values are normalized to the sgControl group in E:T = 0.25 condition. Mean ± SD and individual replicate values are shown for each group (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). E, Relative B16F10 cell number after coculture with Pmel-1 T cells at different E:T ratios revealed a lower sensitivity of Traf3-overexpressing B16F10 cells to T cell–mediated cytotoxicity. The bar plots present the relative cell number in each group, normalized to the cell number in the pEF1α group in each E:T condition. Mean ± SD and individual replicate values are shown for each group (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing pEF1α-Traf3 and pEF1α in each condition). F, Workflow of testing the role of TRAF3 in modulating the ICB response in vivo. 4 × 105 B16F10 cells (Traf3-normal or Traf3-deficient) were transplanted subcutaneously into syngeneic recipient mice. Starting on day 6 after transplantation, we treated the recipients with control IgG or combination ICB every third day for a total of four doses. We monitored tumor size and recipient survival. G, Longitudinal tumor size of sgControl or sgTraf3 tumors treated by control IgG or ICB. Mean ± SEM is shown for each group at each time point (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). H, Kaplan–Meier curves of recipients of sgControl or sgTraf3 tumors treated by control IgG or ICB. The sgTraf3 cohort with ICB treatment survived significantly longer than the other groups (**, P < 0.01; log-rank test with Benjamini–Hochberg adjustment of multiple comparisons).
TRAF3 deletion sensitizes cancer cells to T cell–driven cytotoxicity. A, Workflow of testing the role of TRAF3 in regulating the response of cancer cells to T cell–driven cytotoxicity through in vitro coculture. Traf3-normal or Traf3-deficient B16F10 cells were cultured with antigen-specific CD8+ T cells (Pmel-1 or OT-I) for 1 to 3 days. Cell number and cell surface MHC-I/PD-L1 expression were quantified by FACS. B and C, Relative B16F10 cell number after coculture with (B) OT-I T cells or (C) Pmel-1 T cells at different Effector:Target (E:T) ratios revealed a higher sensitivity of Traf3-deficient B16F10 cells to T-cell–mediated cytotoxicity. For coculture with Pmel-1 T cells, B16F10 cells were pretreated with 1 ng/mL IFNγ for 12 hours prior to the coculture. The bar plots present the relative cell number in each group, normalized to the cell number in the sgControl group in each E:T condition. Mean ± SD and individual replicate values are shown for each group (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). D, Relative MFI of H2-Kb or PD-L1 of B16F10 cells cocultured with OT-I T cells. Values are normalized to the sgControl group in E:T = 0.25 condition. Mean ± SD and individual replicate values are shown for each group (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). E, Relative B16F10 cell number after coculture with Pmel-1 T cells at different E:T ratios revealed a lower sensitivity of Traf3-overexpressing B16F10 cells to T cell–mediated cytotoxicity. The bar plots present the relative cell number in each group, normalized to the cell number in the pEF1α group in each E:T condition. Mean ± SD and individual replicate values are shown for each group (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing pEF1α-Traf3 and pEF1α in each condition). F, Workflow of testing the role of TRAF3 in modulating the ICB response in vivo. 4 × 105 B16F10 cells (Traf3-normal or Traf3-deficient) were transplanted subcutaneously into syngeneic recipient mice. Starting on day 6 after transplantation, we treated the recipients with control IgG or combination ICB every third day for a total of four doses. We monitored tumor size and recipient survival. G, Longitudinal tumor size of sgControl or sgTraf3 tumors treated by control IgG or ICB. Mean ± SEM is shown for each group at each time point (**, P < 0.01; ***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing sgTraf3 and sgControl in each condition). H, Kaplan–Meier curves of recipients of sgControl or sgTraf3 tumors treated by control IgG or ICB. The sgTraf3 cohort with ICB treatment survived significantly longer than the other groups (**, P < 0.01; log-rank test with Benjamini–Hochberg adjustment of multiple comparisons).
We further tested the effect of Traf3-KO on tumor response to ICB treatment in vivo. Specifically, we transplanted Traf3-WT or Traf3-KO B16F10 cancer cells in vivo, treated the recipient mice with control IgG or combined anti–PD-1/anti-CTLA4 therapy, and monitored disease progression over time (Fig. 3F). Consistent with the in vitro coculture results, Traf3-KO cancer cells showed a better response to ICB treatment than Traf3-WT cells in vivo, and recipients of Traf3-KO cancer cells survived significantly longer upon ICB treatment (Fig. 3G and H). Interestingly, Traf3-KO tumors also grew more slowly than TRAF3-normal tumors even without ICB (Fig. 3G). In addition, depletion of CD8 T cells by anti-CD8 treatment ablated the slower tumor growth caused by Traf3-KO, with control IgG or ICB treatment (Supplementary Fig. S4G). This indicated that the effect of Traf3-KO on tumor growth is dependent on the presence of CD8 T cells. To provide further support of our experimental observations, we analyzed multiple published coculture or in vivo CRISPR screens. In general, Traf3-KO tended to sensitize cancer cells to T cell–driven cytotoxicity in those studies (Supplementary Fig. S4H), in agreement with our finding. In contrast, in natural killer (NK)–cell coculture CRISPR screens, Traf3-KO seemed to desensitize cancer cells to NK-cell killing (Supplementary Fig. S4I), consistent with the role of MHC-I in inhibiting NK-cell cytotoxicity (65).
Traf3-KO Signature Is Correlated with Higher MHC-I in Primary Patient Samples
Next, we evaluated whether TRAF3 function was correlated with MHC-I expression or T-cell infiltration in clinical cohorts where bulk tumor RNA-seq has been performed. Because TRAF3 regulates downstream TF activity, rather than directly using TRAF3 expression itself to measure TRAF3 functional deficiency, we derived a Traf3-KO gene expression signature from RNA-seq (Supplementary Fig. S5A; Supplementary Table S3). Specifically, we used the top differential genes from TRAF3 deficiency as signature genes, and evaluated the Traf3-KO signature score in each clinical tumor RNA-seq sample by a weighted sum of signature gene expression. We further used it to assess the association of TRAF3 deficiency with MHC-I antigen presentation, cytotoxic T-cell infiltration, overall survival, and immunotherapy response in clinical cohorts as appropriate.
We first examined TCGA tumors, where baseline antitumor immunity can influence prognosis, although patients were not treated with ICB (66, 67). We found the Traf3-KO signature level to be positively correlated with MHC-I expression, CD8+ T-cell infiltration, and patient survival in melanoma (Fig. 4A–C). We also applied an orthogonal algorithm TIDE (66), which we developed to predict ICB response based on expression signatures of T-cell dysfunction and exclusion, to TCGA tumors as if they were to be treated with ICB. The predicted responders also showed higher Traf3-KO signature score than the predicted nonresponders (Fig. 4D). Furthermore, we observed positive correlation of Traf3-KO signature and MHC-I expression, CD8 T-cell infiltration, and patient survival in multiple other cancer types in TCGA (Supplementary Fig. S6A–S6D). Finally, to assess the MHC-I–specific effect of Traf3-KO, we created a restricted Traf3-KO signature to include just the differential genes in the MHC-I component genes. This MHC-I–restricted Traf3-KO signature still shows positive correlation with patient survival in many cancer types in TCGA (Supplementary Fig. S6E). These results suggest that the MHC-I–induction effect by TRAF3 inhibition is associated with better clinical outcome.
Traf3-KO signature is correlated with higher MHC-I in primary patient samples. A–C, Traf3-KO signature is positively correlated with (A) MHC-I expression, (B) CD8+ T-cell infiltration, and (C) patient survival in the TCGA-SKCM dataset. PCC, Pearson correlation coefficient; TPM, Transcripts Per Kilobase Million. D, TIDE-predicted ICB responders showed higher Traf3-KO signature values in the TCGA SKCM (skin cutaneous melanoma) dataset. E, Expression of MHC-I component or related genes in MHC-I–high or MHC-I–low samples. F, Gene set enrichment analysis of differentially expressed genes in MHC-I–high versus MHC-I–low samples. G, MHC-I–high samples show higher expression of genes involved in antigen presentation, NFκB signaling, TNF signaling, and Toll-like receptor pathways. H, Traf3-KO signature score in MHC-I–high or MHC-I–low RNA-seq samples. I, H3K27ac ChIP-seq results for the HLA-A/B/C loci in MHC-I–high or MHC-I–low samples. J, Top 200 genes with the higher regulatory potential value in MHC-I–high versus MHC-I–low H3K27ac ChIP-seq samples.K, Traf3-KO signature score of H3-K27ac ChIP-seq regulatory potential values in MHC-I–high or MHC-I–low samples. L, Cistrome-GO enrichment analysis of H3K27ac ChIP-seq peaks with stronger signal in MHC-I–high compared with MHC-I–low samples. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
Traf3-KO signature is correlated with higher MHC-I in primary patient samples. A–C, Traf3-KO signature is positively correlated with (A) MHC-I expression, (B) CD8+ T-cell infiltration, and (C) patient survival in the TCGA-SKCM dataset. PCC, Pearson correlation coefficient; TPM, Transcripts Per Kilobase Million. D, TIDE-predicted ICB responders showed higher Traf3-KO signature values in the TCGA SKCM (skin cutaneous melanoma) dataset. E, Expression of MHC-I component or related genes in MHC-I–high or MHC-I–low samples. F, Gene set enrichment analysis of differentially expressed genes in MHC-I–high versus MHC-I–low samples. G, MHC-I–high samples show higher expression of genes involved in antigen presentation, NFκB signaling, TNF signaling, and Toll-like receptor pathways. H, Traf3-KO signature score in MHC-I–high or MHC-I–low RNA-seq samples. I, H3K27ac ChIP-seq results for the HLA-A/B/C loci in MHC-I–high or MHC-I–low samples. J, Top 200 genes with the higher regulatory potential value in MHC-I–high versus MHC-I–low H3K27ac ChIP-seq samples.K, Traf3-KO signature score of H3-K27ac ChIP-seq regulatory potential values in MHC-I–high or MHC-I–low samples. L, Cistrome-GO enrichment analysis of H3K27ac ChIP-seq peaks with stronger signal in MHC-I–high compared with MHC-I–low samples. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
In addition to publicly available databases, we also collected primary samples from treatment-naïve patients with melanoma treated at the Dana-Farber Cancer Institute. Assessment of MHC-I levels by immunohistochemistry (IHC) of the formalin-fixed paraffin-embedded (FFPE) samples identified MHC-I–high and MHC-I–low tumors, from which we macrodissected cancer cell sections for RNA-seq and H3K27 acetylation fixed tissue ChIP-seq (FitAc-seq; ref. 68) to examine the molecular features associated with MHC-I expression. RNA-seq results on MHC-I components (HLA-A, HLA-B, HLA-C, and B2M) and other genes in the same pathway were consistent with the IHC results (Fig. 4E). Examination of the differentially expressed genes revealed 673 upregulated and 425 downregulated genes in the MHC-I–high compared with the MHC-I–low samples, with the upregulated genes being enriched in NFκB and TNF signaling pathways (Fig. 4F and G; Supplementary Table S4). Indeed, MHC-I–high samples showed a significantly higher Traf3-KO signature score (Fig. 4H). Consistent with the RNA-seq data, ChIP-seq revealed higher levels of H3K27ac near MHC-I and higher regulatory potential of MHC-I–related genes (Fig. 4I and J; Supplementary Table S4). The regulatory potential of Traf3-KO signature genes was also higher in the MHC-I–high samples (Fig. 4K). Cistrome GO analysis of the differentially active chromosomal regions from ChIP-seq revealed that NFκB binding sites are more active in the MHC-I–high samples (Fig. 4L). These results suggest that the negative regulation of MHC-I by TRAF3 is also reflected in the primary melanoma tumors from patients.
Traf3-KO Signature Is Correlated with Better ICB Response
In patients treated with immunotherapy, antigen expression and/or presentation levels are known to be a critical factor of treatment response and an issue addressed in multiple response biomarker studies (12, 69). Higher MHC-I expression is associated with better ICB response in most published ICB cohorts (Supplementary Fig. S7A). Based on the role of TRAF3 in MHC-I regulation, we tested the correlation between TRAF3 functional deficiency and immunotherapy response in the published ICB treatment cohorts. We examined the RNA-seq data of pretreatment biopsies in a recent melanoma cohort with 63 biopsies treated with anti–PD-1 monotherapy and 57 treated with combined anti–PD-1 and anti-CTLA4 (70). We confirmed that the Traf3-KO signature was positively correlated with tumor MHC-I, cytotoxic T-cell infiltration, ICB response, and progression-free survival (Fig. 5A–D). Because the MHC-I–TCR interaction is a conserved mechanism for T-cell recognition of cancer cells, we expect the Traf3-KO signature to be associated with ICB response in various cancer types.
Traf3-KO signature is correlated with better response to ICB. A–D, Traf3-KO signature is positively correlated with (A) MHC-I expression, (B) intratumoral cytotoxic T lymphocyte (CTL) infiltration, (C) ICB response, and (D) overall survival and progression-free survival in patients treated by anti–PD-1 or combined anti–PD-1 and anti-CTLA4 in the Gide and colleagues (70) study in melanoma. E–H, Traf3-KO signature is positively correlated with (E) MHC-I expression, (F) intratumoral CTL infiltration, (G) ICB response, and (H) overall survival and progression-free survival in patients treated by anti–PD-1 or combined anti–PD-1 and anti-CTLA4 in the Mariathasan and colleagues (71) study in urothelial carcinoma. I, Traf3-KO signature is positively correlated with MHC-I expression and cytotoxic T-cell infiltration and negatively correlated with survival hazard in most ICB treatment clinical trials. *, P < 0.05; **, P < 0.01; ***, P < 0.001. PD, progressive disease; SD, stable disease; PR, partial response; CR, complete response.
Traf3-KO signature is correlated with better response to ICB. A–D, Traf3-KO signature is positively correlated with (A) MHC-I expression, (B) intratumoral cytotoxic T lymphocyte (CTL) infiltration, (C) ICB response, and (D) overall survival and progression-free survival in patients treated by anti–PD-1 or combined anti–PD-1 and anti-CTLA4 in the Gide and colleagues (70) study in melanoma. E–H, Traf3-KO signature is positively correlated with (E) MHC-I expression, (F) intratumoral CTL infiltration, (G) ICB response, and (H) overall survival and progression-free survival in patients treated by anti–PD-1 or combined anti–PD-1 and anti-CTLA4 in the Mariathasan and colleagues (71) study in urothelial carcinoma. I, Traf3-KO signature is positively correlated with MHC-I expression and cytotoxic T-cell infiltration and negatively correlated with survival hazard in most ICB treatment clinical trials. *, P < 0.05; **, P < 0.01; ***, P < 0.001. PD, progressive disease; SD, stable disease; PR, partial response; CR, complete response.
We therefore applied our findings to another study with approximately 300 urothelial carcinoma biopsies (71), the largest published cohort with available expression data. Again, Traf3-KO signature was positively correlated with tumor MHC-I, cytotoxic T-cell infiltration, ICB response, and progression-free survival (Fig. 5E–H). Furthermore, the complete or MHC-I–restricted Traf3-KO signature showed a similar trend (Supplementary Table S5; Fig. 5I; Supplementary Fig. S7B–S7E) in most published ICB clinical studies across different cancer types. Together, our data support the general relevance of Traf3-KO signature in predicting ICB response.
Identification of a SMAC Mimetic, Birinapant, as a Regulator of MHC-I to Enhance ICB Efficacy
Based on the effect of Traf3-KO on MHC-I levels and sensitivity to T cell–driven killing, we hypothesized that drug treatment that mimics the transcriptional effect of Traf3-KO may achieve similar outcomes. We searched existing transcriptomic studies on drug-treatment effects in the GEO database, with differential expression data involving more than 200 targeted therapy/chemotherapy drugs. Different drug treatments showed distinct transcriptional regulation of MHC-I components and PD-L1/PD-L2 (Fig. 6A; Supplementary Fig. S8A; Supplementary Table S6). We computed the Traf3-KO signature score of each drug treatment, and, reassuringly, it is significantly positively correlated with MHC-I expression (Supplementary Fig. S8B). Focusing on drugs that increase Traf3-KO signature and MHC-I expression, we found several drugs previously reported to upregulate MHC-I, including the CDK4/6 inhibitor palbociclib (38) and the EZH2 inhibitor EPZ6438 (ref. 44; Supplementary Table S6). Furthermore, we found multiple Second Mitochondria-derived Activator of Caspase (SMAC) mimetics as novel candidates for MHC-I–specific regulation (Fig. 6B). SMAC mimetics are a family of proapoptotic compounds that can inhibit the cIAP/XIAP E3 ligases, leading to activated NFκB activity. We therefore examined the immunomodulatory effects of a SMAC mimetic, birinapant, using in vitro and in vivo assays.
Birinapant can specifically upregulate MHC-I and add to the efficacy of ICB treatment. A, Heat map of the regulation of transcription of MHC-I components and PD-L1/PD-L2 in each drug treatment condition. Data are ordered by the magnitude of MHC-I regulation from the most upregulated (left) to the most downregulated (right). B, Characterization of each drug treatment for effect on MHC-I expression and Traf3-KO signature value. C, Flow cytometry of B16F10 cells treated by IFNγ and/or different concentrations of the SMAC mimetic birinapant (Bir). Typical histograms of H2-Kb and PD-L1 levels in each condition are shown. D, The relative MFI of samples in C. Mean ± SD is shown for each group. Values are normalized to the mean MFI at 0 ng/mL IFNγ and 0 μmol/L birinapant. E, Relative B16F10 cell number after coculture with OT-I T cells at different E:T ratios revealed a higher sensitivity of B16F10 cells to T cell–mediated cytotoxicity under birinapant treatment. The bar plots present the relative cell number in each group, normalized to the cell number in the vehicle treatment group in each E:T condition. Mean ± SD and individual replicate values are shown for each group (***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test). F, Relative MFI of H2-Kb or PD-L1 of B16F10 cells cocultured with OT-I T cells. Values are normalized to the vehicle treatment group in E:T = 0.25 condition. Mean ± SD and individual replicate values are shown for each group (***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing vehicle and SMAC treatment groups in each condition). G, Longitudinal tumor size of tumors under different treatments. Mean ± SEM is shown for each group at each time point (*, P < 0.05; **, P < 0.01; two-way ANOVA with Benjamini–Hochberg post test). H, Gene set enrichment analysis of the bulk tumor RNA-seq data evaluating the gene sets that were upregulated in response to SMAC mimetic. I, Quantification of Traf3-KO signature in each treatment group revealed a trend toward upregulation in response to birinapant treatment.
Birinapant can specifically upregulate MHC-I and add to the efficacy of ICB treatment. A, Heat map of the regulation of transcription of MHC-I components and PD-L1/PD-L2 in each drug treatment condition. Data are ordered by the magnitude of MHC-I regulation from the most upregulated (left) to the most downregulated (right). B, Characterization of each drug treatment for effect on MHC-I expression and Traf3-KO signature value. C, Flow cytometry of B16F10 cells treated by IFNγ and/or different concentrations of the SMAC mimetic birinapant (Bir). Typical histograms of H2-Kb and PD-L1 levels in each condition are shown. D, The relative MFI of samples in C. Mean ± SD is shown for each group. Values are normalized to the mean MFI at 0 ng/mL IFNγ and 0 μmol/L birinapant. E, Relative B16F10 cell number after coculture with OT-I T cells at different E:T ratios revealed a higher sensitivity of B16F10 cells to T cell–mediated cytotoxicity under birinapant treatment. The bar plots present the relative cell number in each group, normalized to the cell number in the vehicle treatment group in each E:T condition. Mean ± SD and individual replicate values are shown for each group (***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test). F, Relative MFI of H2-Kb or PD-L1 of B16F10 cells cocultured with OT-I T cells. Values are normalized to the vehicle treatment group in E:T = 0.25 condition. Mean ± SD and individual replicate values are shown for each group (***, P < 0.001; two-way ANOVA with Benjamini–Hochberg post test comparing vehicle and SMAC treatment groups in each condition). G, Longitudinal tumor size of tumors under different treatments. Mean ± SEM is shown for each group at each time point (*, P < 0.05; **, P < 0.01; two-way ANOVA with Benjamini–Hochberg post test). H, Gene set enrichment analysis of the bulk tumor RNA-seq data evaluating the gene sets that were upregulated in response to SMAC mimetic. I, Quantification of Traf3-KO signature in each treatment group revealed a trend toward upregulation in response to birinapant treatment.
We first tested the effect of birinapant treatment on MHC-I levels in mouse B16F10 and CT26 cells. Indeed, birinapant treatment specifically upregulated MHC-I in response to IFNγ treatment, with minimal effect on PD-L1 (Fig. 6C and D; Supplementary Fig. S8C). Furthermore, birinapant can also specifically upregulate MHC-I in multiple human cancer cell lines (Supplementary Fig. S8D–S8G). Similar to the effect of Traf3-KO, birinapant treatment led to reduced inactive forms of NFκB (p100 and p105) and increased active form of NFκB (p52; Supplementary Fig. S8H). Consistent with the reported role of SMAC mimetics to induce cIAP degradation (72, 73), we found lower cIAP and higher NIK levels in birinapant-treated cells (Supplementary Fig. S8I). When added to the cancer cells cocultured with OT-I T cells, birinapant resulted in specific upregulation of MHC-I in cancer cells and facilitated better cancer cell killing by T cells (Fig. 6E and F). We further tested whether this immunomodulatory effect of birinapant functions through TRAF3 in cancer cells by using Traf3-WT or Traf3-KO cancer cells cocultured with antigen-specific T cells, treated by vehicle or birinapant. TRAF3-deficient cancer cells showed comparable high levels of MHC-I expression and sensitivity to T cell–mediated cytotoxicity as birinapant-treated cells, and were not further sensitized by birinapant (Supplementary Fig. S8J and S8K). These data suggest that birinapant's effect on MHC-I expression is dependent on the presence of TRAF3.
We further tested the effect of birinapant in vivo by treating established B16F10 melanoma tumors with birinapant and/or combination anti–PD-1/CTLA4. We observed moderately reduced tumor growth in vivo with birinapant treatment alone, and further reduced tumor progression with combination therapy of birinapant and ICB (Fig. 6G, left). In addition, CD8 T-cell depletion by anti-CD8 treatment ablated the slower tumor growth caused by birinapant treatment (Fig. 6G, right), indicating the effect of birinapant treatment on tumor growth to be CD8 T-cell–dependent. We performed RNA-seq of the bulk tumor in each treatment group to identify the pathways regulated by birinapant. Indeed, birinapant-upregulated genes, whether with ICB treatment or not, are involved in inflammatory response, innate/adaptive immune response, chemokine/cytokine-mediated signaling, NFκB signaling, and MHC-I antigen presentation (Fig. 6H; Supplementary Table S7). Inference of tumor immune infiltration by TIMER (74, 75) predicted a trend for higher infiltration of dendritic cells and CD8+ T cells, though the difference is not statistically significant (Supplementary Fig. S8L and S8M). In addition, birinapant treatment increased the Traf3-KO signature in the tumor (Fig. 6I, left), consistent with our in silico prediction and in vitro epistatic experiment showing that SMAC mimetic treatment phenocopies Traf3-KO. In contrast, we did not observe a significant difference in Traf3-KO signature between birinapant + ICB and ICB groups (Fig. 6I, right). Because we evaluated tumor RNA-seq after ICB treatment, the cancer cells that show higher Traf3-KO signature were likely more susceptible to immune killing, making it difficult to detect such signature through bulk tumor RNA-seq. This ICB-induced selection for immune-resistant clones was characterized in our recent work applying clonal tracing of cancer cells following ICB treatment (76). Furthermore, the presence of other cell types in the tumor microenvironment might obscure the enrichment of Traf3-KO signature.
Discussion
Immunotherapy, especially ICB, has brought paradigm shifts to cancer treatment and improved patient survival in many cancer types. Its therapeutic effect relies on the recognition of cancer-specific antigens by cytotoxic T cells, and defect in cancer cell antigen presentation is a major ICB-resistance mechanism. Multiple studies have sought to identify regulators of MHC-I (36, 44). However, due to the common regulatory pathways shared by MHC-I and PD-L1, upregulating MHC-I often leads to upregulation of PD-L1, which suppresses antitumor immunity. In this study, we established a workflow to identify drugs with MHC-I–specific enhancing effects. We first performed dual-marker FACS-based CRISPR screens to identify regulators of MHC-I and/or PD-L1, and found TRAF3 as an MHC-I–specific negative regulator through its modulation of the NFκB pathway. TRAF3 depletion sensitized cancer cells to antigen-specific T cell–driven cytotoxicity. Traf3-KO signature is significantly correlated with higher expression of genes in the antigen presentation pathway, higher cytotoxic T-cell infiltration, and better overall survival in patients with cancer from many TCGA cohorts and ICB clinical trials. Furthermore, data mining of public transcriptomic profiles of drug treatment followed by experimental validation identified that birinapant, a SMAC mimetic, phenocopies the effect of TRAF3 deficiency on MHC-I expression and adds to the efficacy of ICB treatment in suppressing tumor growth. Our study identifies a novel treatment strategy that can potentially be used in MHC-I–low cancers to enhance ICB therapeutic efficacy.
TRAF3 has been reported to negatively regulate the NFκB pathway activity, and multiple ChIP-seq studies revealed NFκB binding sites in the promoter/enhancer regions of MHC-I components. Consistent with these previous studies, we found that Traf3-KO upregulates MHC-I in an NFκB-dependent manner. There are no published TRAF3 small-molecule inhibitors, so we searched published studies to identify drugs that phenocopy the effect of TRAF3 inhibition. The SMAC mimetic family ranked among the top of the list. SMAC mimetics can trigger the autoubiquitination and degradation of Cellular Inhibitor of Apoptosis Proteins (cIAP), thereby promoting NFκB activation via canonical or noncanonical pathways (49, 77). We then experimentally validated that birinapant, a SMAC mimetic, can promote the immunosensitivity of cancer cells through specifically inducing MHC-I expression in a TRAF3/NFκB-dependent manner. A recent study also showed that SMAC mimetics can synergize with ICB in a CD8 T cell–dependent and TNFα-dependent manner (78), although SMAC mimetics did not upregulate MHC-I in their glioblastoma model (CT-2A). Our preliminary analyses of published single cell RNA-seq (scRNA-seq) data from various cancer types (79) found that glioblastoma cells tend to have high expression of MHC-I and NFκB component genes. We speculate that cells already high in MHC-I and NFκB might not increase MHC-I when treated with SMAC mimetics. Indeed, in A375 and 786O cell lines with higher baseline MHC-I and NFκB expression, birinapant also failed to increase MHC-I. These results suggest context-dependent MHC-I regulation and antitumor response induced by SMAC mimetics, which merit further mechanistic studies.
In addition to SMAC mimetics, mining the GEO database identified multiple other targeted therapy or chemotherapy agents with the potential to cause a transcriptional effect similar to TRAF3 inhibition. This comports with multiple other studies (37–39, 80, 81) in revealing the immunomodulatory effects of chemotherapy and targeted therapies. For example, palbociclib, a CDK4/6 inhibitor reported to induce MHC-I via type-III IFN (38), and EPZ6438, an EZH2 inhibitor reported to induce MHC-I via chromatin modification (36, 44), ranked among the top candidates. Regimens containing azacitidine, another top-ranking candidate, were also reported to induce antigen presentation in non–small cell lung cancer cells (82). Our approach using differential gene expression signatures can be extended to identify drugs with other immunomodulatory effects.
Our study has limitations that future studies could address. First, our finding can potentially assist the treatment of tumors whose MHC-I is suppressed through gene regulatory mechanisms, but not those with genetic deletion of MHC-I components or loss of heterozygosity of HLA genes. Nevertheless, the clinical observation in melanoma treated by ICB (13) suggests that a large number of patients can still benefit from this approach. Second, in addition to upregulation of MHC-I, other mechanisms may also contribute to the immune sensitivity of Traf3-KO tumors. For example, Traf3-KO cells had higher expression of chemokines which may facilitate T-cell infiltration. Further studies are needed to elucidate such alternative mechanisms and their impact on the immune sensitivity of Traf3-KO tumors. Third, when mining drug treatment data to predict drugs that may upregulate MHC-I through TRAF3/NFκB, our analysis may be limited by the amount of drug treatment data available in GEO. More publicly available transcriptional profiling data on drug screening would improve the power of our approach.
In summary, we present a workflow integrating CRISPR screens with functional genomics characterization and computational mining of publicly accessible data to identify novel regulators of tumor immunity. We showed that the SMAC mimetic birinapant can phenocopy TRAF3 deficiency and preferentially induce MHC-I expression, leading to enhanced sensitivity of cancer cells to T cell–driven cytotoxicity. The combination of SMAC mimetic and ICB is likely to benefit patients with a low baseline level of MHC-I expression. This is potentially exciting given the biosafety of profiles of birinapant from multiple early-phase clinical studies. Furthermore, our general approach can be applied to identify other drugs with immunomodulatory effects.
Methods
TCGA Cohorts
We downloaded all the TCGA tumor data from the Broad GDAC Firehose, including transcriptomic profiles and clinical information of tumor samples, from 33 cancer types. All of the gene expression levels were measured by logarithmic Transcripts Per Kilobase Million (TPM). We downloaded the IHC-derived tumor purity data from the study by Aran and colleagues (83), and we also estimated the tumor purity using the CHAT algorithm (84). For all the tumors, we estimated the immune cell infiltration in the tumor microenvironment using TIMER algorithm, with default parameters and TPM as input. Because of the distinct clinical behavior and molecular profile of cancer subtypes, we divided breast cancer datasets into the PAM50 subtypes of basal, lumA, lumB, and HER2 positive. We focused on cancer types where higher MHC-I is associated with better overall survival (CoxPH regression of patient survival proportional hazard and MHC-I expression with a z-score < –0.5) for downstream analyses.
CCLE Cohorts
We downloaded the RNA-seq data of more than 1,000 cancer cell lines from the DepMap Portal, which contains processed gene expression data quantified according to the Genotype-Tissue Expression pipelines.
Flow Cytometry
We performed FACS of B16F10 or CT26 cells treated under different conditions as indicated to quantify the level of H2-Kb/H2-Kd and PD-L1. B16F10 or CT26 cells were dissociated by TrypLE treatment, washed with PBS-2% FBS, and incubated with DAPI (1:10,000 dilution), anti–H2-Kb (clone AF6-88.5, BioLegend, 1:400 dilution), and anti–PD-L1 (clone MIH5, BD Biosciences, 1:400 dilution) for 1 hour on ice. Cells were then washed and resuspended in PBS-2% FBS and analyzed on BD LSR-Fortessa instrument. FACS data were then analyzed by Flowjo software.
Cell Lines and Cell Culture
Cell lines B16F10, CT26, SKMEL5, and MCF7 were purchased from the ATCC and authenticated using standard short tandem repeat analysis in 2019. Cell line K028 was derived from a patient with melanoma under Institutional Review Board–approved protocols. Cell line HT29 was obtained from the CCLE core facility, which obtained them directly from commercial sources and authenticated the lines using standard short tandem repeat analysis. 293FT, B16F10, SKMEL5, K029, and MCF7 were cultured in DMEM supplemented with 10% FBS, 1% l-glutamine, and 1% penicillin and streptomycin. CT26 and HT29 were cultured in RPMI supplemented with 10% FBS, 1% l-glutamine, and 1% penicillin and streptomycin. All cells were used at low passage numbers and were tested for Mycoplasma using the MycoAlert Mycoplasma Detection Kit (Lonza).
CRISPR Screens for Regulators of MHC-I
We performed CRISPR screens to identify regulators of MHC-I induced by IFNγ. We first transduced 1 × 108 B16F10 cells with our M1/M2 mouse genome-wide CRISPR library at multiplicity of infection of 0.3 to ensure that most transduced cells received only 1 virion. We cultured the transduced cells for 2 days in DMEM full media, selected with 2 μg/mL puromycin for 2 days, and cultured for an additional 2 days to allow depletion of residual protein of the targeted genes. We then treated the cells with high- or low-dosage IFNγ to identify MHC-I regulators.
For identification of positive regulators, we treated the transduced cells with 10 ng/mL IFNγ for 2 days and flow-sorted 2 × 107 cells for the following populations: MHC-Ihi PD-L1hi, MHC-Ihi PD-L1lo, MHC-Ilo PD-L1hi, MHC-Ilo PD-L1lo. The sgRNAs that are present in the MHC-Ilo PD-L1hi population are expected to enrich for MHC-I–specific positive regulators. For identification of negative regulators, we treated the transduced cells with 0.1 ng/mL IFNγ for 2 days, and flow-sorted 2 × 107 cells for the following populations: MHC-Ihi PD-L1hi, MHC-Ihi PD-L1lo, MHC-Ilo PD-L1hi, MHC-Ilo PD-L1lo. The gRNAs that are present in the MHC-Ihi PD-L1lo population are expected to enrich for MHC-I–specific negative regulators.
Genomic DNA extraction and gRNA library construction were performed as described before (85). Briefly, we extracted genomic DNA from each sorted population or from the unsorted population by phenol–chloroform extraction, and then performed two rounds of PCR to construct the barcode-indexed sequencing library for each sample. We sequenced each library at an approximately 300× average coverage over the CRISPR library.
Analysis of CRISPR Screen
MAGeCK software was used to quantify and test for sgRNA and gene enrichment (41). The sequence reads were trimmed to remove a constant portion of the sgRNA sequences and mapped to the M1/M2 mouse genome-wide CRISPR library with MAGeCK “count” module, which computes read counts for each sgRNA. We used MAGeCK “test” module to identify MHC-I/PD-L1 regulators, for which targeting sgRNA was significantly enriched in the sorted populations compared with the control unsorted populations grown in parallel. MAGeCK returned sgRNA and gene log fold changes, which represent the enrichment level of the sgRNAs and genes in each cell population. Because negative selection is meaningless in FACS-based CRISPR screens, we assigned zero log fold change for genes whose log fold changes were negative. Finally, we calculated the mean log fold change of genes in M1 and M2 CRISPR screens and visualized the results using functions in R package MAGeCKFlute (42).
We analyzed the data from the cohort treated by 10 ng/mL IFNγ to identify the positive regulators of MHC-I and/or PD-L1. Specifically, we compared the gRNA distribution between the MHC-Ilo PD-L1lo population and unsorted population to identify positive regulators for both; between the MHC-Ilo PD-L1hi population and the unsorted population to identify MHC-I–specific positive regulators; and between the MHC-Ihi PD-L1lo population and the unsorted population to identify PD-L1–specific positive regulators.
We analyzed the data from the cohort treated by 0.1 ng/mL IFNγ to identify the negative regulators of MHC-I and/or PD-L1. Specifically, we compared the gRNA distribution between the MHC-Ihi PD-L1hi population and unsorted population to identify negative regulators for both; between the MHC-Ihi PD-L1lo population and the unsorted population to identify MHC-I–specific negative regulators; and between the MHC-Ilo PD-L1hi population and the unsorted population to identify PD-L1–specific negative regulators.
Cloning of Screen Candidates for Validation
We cloned four gRNAs of Traf3 into the LentiCRISPRv2-Puro construct to validate the role of TRAF3 in the regulation of MHC-I. We used the cloning protocol as previously established (86, 87). Out of the four gRNAs, two successfully depleted TRAF3 protein level >80%. We used these two gRNAs for our follow-up validation studies.
We also cloned the open reading frame of Traf3 into pER1α-Puro vector for overexpression of Traf3 and examination of its role in MHC-I/PD-L1 regulation.
Validation of the Role of TRAF3 in MHC-I/PD-L1 Regulation
We performed CRISPR/Cas9–mediated deletion of Traf3 in multiple cell lines, including B16F10 and CT26. For each line, we cultured Traf3-normal or Traf3-deficient cells with different dosages of IFNγ, TNFα, IFNα, or IFNβ for 2 days to test its regulation of MHC-I and/or PD-L1 with these induction methods.
Virus Packaging
Viral packaging was performed as previously described (85). Briefly, we cultured 293FT cells at 40% to 50% confluence the day before transfection. Transfection was performed using X-tremeGENE HD (Sigma-Aldrich). For each 10-cm dish, 8 μg of lentivectors, 2.4 μg of pMD2.G, and 6 μg of psPAX2 were added into 2 mL Opti-MEM (Life Technologies). Forty microliters of X-tremeGENE HD was added to the plasmid mixture. The complete mixture was incubated for 15 minutes and then added to 1.1 × 107 293FT cells per 10-cm dish. After overnight culture, the media were changed to 8 mL DMEM + 15% FBS for virus collection. We harvested virus at 48, 60, and 72 hours after the start of transfection.
Coculture of Cancer Cells and T Cells for T-cell Cytotoxicity Assay
Coculture of cancer cells and T cells was performed as previously described (66). Briefly, B16F10 cells were maintained in complete DMEM (10% FBS and 50 U/mL of penicillin–streptomycin). CD8 T cells isolated from spleen and lymph nodes from Pmel-1 or OT-I mice were stimulated with anti-CD3/CD28 beads (Thermo Fisher, 11452-D) and then cultured in complete RPMI 1640 media (10% FBS, 20 mmol/L HEPES, 1 mmol/L sodium pyruvate, 0.05 mmol/L 2-mercaptoethanol, 2 mmol/L l-glutamine, 50 U/mL streptomycin and penicillin, and 20 ng/mL recombinant mouse IL2).
To test the sensitivity of cancer cells to T cell–driven cytotoxicity (OT-I or Pmel-1 model), we plated B16F10 cells (sgRosa26, sgTraf3, pEF1a-Empty, or pEF1a-Traf3) at equal density in all wells and added T cells at ratios to cancer cells. With the OT-I model, we first incubated the B16F10 cells with 1 nmol/L SIINFEKL peptide for 2 hours prior to coculture with T cells. With the Pmel-1 model, B16F10 cells were either pretreated with 1 ng/mL IFNγ overnight or untreated prior to the coculture. There are two to four cell-culture replicates for each condition. After a 1-day or 3-day coculture with T cells, we counted the remaining cancer cells by FACS using the precision count beads (BioLegend, 424902). T cells present in these cultures were gated out based on antibodies specific for CD45 (BioLegend, clone 30-F11) or CD8 (BioLegend, clone 53-6.7).
Immunoblot
We treated B16F10 cells (sgRosa26 or sgTraf3) by 1 ng/mL IFNγ for 0, 0.5, or 24 hours and lysed cells in RIPA buffer. Lysates were resolved by SDS-PAGE, and immunoblotting was performed as described (88). To visualize NIK levels, cells were pretreated by 0 or 5 μmol/L MG132 for 4 hours before lysis. Antibodies specific for TRAF3, STAT1, phospho-STAT1(Y705), RELA, RELB, phospho-RELA(S536), cREL, p100/p52, p105/p50, IKK, and NIK were from Cell Signaling Technology. Antibody specific for VCL was from Santa Cruz Biotechnology. Antibody specific for cIAP was from R&D Systems. Antibodies were used at each manufacturer's recommended concentrations.
RNA-seq
For comparison of transcription profile between control and Traf3−/− B16F10 cells, we treated sgControl or sgTraf3-transduced B16F10 cells with vehicle control or 1 ng/mL IFNγ for 2 days, and extracted total RNA using the miRNeasy Plus Mini Kit (QIAGEN, #17004) following the manufacturer's protocol. For comparison of tumors treated by birinapant and/or ICB, we harvested bulk tumors from recipients in each group on day 15 after transplantation and extracted total RNA. Total RNA samples in each experiment were submitted to Novogene Inc. for sequencing. Standard mRNA library preparation kit was used for library preparation. Paired-end 150 bp sequencing was performed on Illumina HiSeq 2500. Sequencing reads were mapped to the mm10 genome by RSEM. Statistics for differentially expressed genes were calculated by DESeq2.
Traf3-KO Signature
To study the clinical relevance of TRAF3, we established a Traf3-KO signature by extracting the top 200 upregulated and 200 downregulated genes and taking the normalized DESeq2 Wald statistics as weights. The weights of genes are normalized to –1 to 1 range by the equation ki = wi/max(w), where ki indicates the weight of the i th gene, and wi indicates the DESeq2 Wald statistics of the i th gene. For each input expression profile, we computed a Traf3-KO signature score to estimate the Traf3-KO level by calculating the weighted sum expression of the signature genes following the equation ) where S denotes the signature score, and Xi denotes the expression level of the i th gene. Finally, we evaluated the association of TRAF3 deficiency with MHC-I expression, tumor-infiltrating T cells, patient outcome, and response to ICB.
ATAC-seq
For comparison of chromatin accessibility between control and Traf3−/− B16F10 cells, we treated sgControl or sgTraf3-transduced B16F10 cells with vehicle control or 1 ng/mL IFNγ for 2 days. ATAC-seq was performed in duplicates (200,000 cells) by the Center for Functional Cancer Epigenetics at Dana-Farber Cancer Institute as previously described (76, 89). For data analysis, we used Burrows–Wheeler Aligner (90) to map sequencing reads to the reference genome and MACS2 (91) for peak calling. DESeq2 (92) was applied to identify the differentially accessible regions between ICB-resistant lines and the parental control line.
Mouse Experiments
All mice were housed in standard cage in Dana-Farber Cancer Institute Animal Resources Facility (ARF). All animal procedures were carried out under the ARF Institutional Animal Care and Use Committee (IACUC) protocol and were in accordance with the IACUC standards for the welfare of animals. WT C57BL/6 recipient mice were purchased from Charles River laboratory.
To test the function of TRAF3 in tumor's response to combination anti–PD-1/anti-CTLA4 treatment, we transplanted 4 × 105 sgRosa26 or sgTraf3 B16F10 cells subcutaneously into the right flank of 6- to 8-week-old female C57BL/6 mice (from Charles River Laboratory). Starting day 6 after transplantation, we began treatment with control IgG [clone 2A3, 200 μg per mouse in 200 μL Hank's Balanced Salt Solution (HBSS) buffer] or combination anti–PD-1 (clone 1A12, 100 μg per mouse) and anti-CTLA4 (clone 9D9, 100 μg per mouse) by i.p. injection every third day for a total of 4 times. We monitored tumor growth along the treatment.
To test the efficacy of ICB and SMAC mimetic combination treatment, we transplanted 4 × 105 parental B16F10 cells subcutaneously into the left and right flanks of 6- to 8-week-old female C57BL/6 mice (from Charles River Laboratory). Starting day 5 after transplantation, we began treatment with vehicle or birinapant (dissolved in 15% captisol, 600 μg per mouse) by i.p. injection every third day for a total of 3 times. Starting day 6 after transplantation, we began treatment with control IgG (clone 2A3, 200 μg per mouse in 200 μL HBSS buffer) or combination anti–PD-1 (clone 1A12, 100 μg per mouse) and anti-CTLA4 (clone 9D9, 100 μg per mouse) by i.p. injection every third day for a total of 3 times. We monitored tumor growth along the treatment and the survival of recipients.
Mice were euthanized using CO2 inhalation. All mouse experiments were performed in compliance with institutional guidelines as approved by the IACUC of Dana-Farber Cancer Institute. The maximum tumor diameter permitted under the relevant animal protocols is 20 mm, which was not exceeded in any experiment.
Profiling of Primary MHC-I–High or MHC-I–Low Melanoma Samples
For comparison of transcription and epigenetic profile between MHC-I–high and MHC-I–low primary melanoma tumors, we performed RNA-seq and FiTAc-seq of human tumor biopsies. Samples were grouped into MHC-I–high or MHC-I–low based on IHC of tumor biopsy from the same sampling.
FiTAc-seq: H3K27Ac ChIP-seq on FFPE Samples
FiTAc-seq (68) was used to profile H3K27ac ChIP-seq on FFPE samples. Sample selection was performed with dual-IHC staining assessing tumor enrichment with SOX10 antibody and MHC-I expression status. Accordingly, areas of interest were marked, and samples were classified into low or high MHC-I categories. The selected areas of interest were manually macrodissected from 10 unstained sections (10 μmol/L each) per case. Briefly, the macrodissected tissues were washed with xylenes to remove paraffin, rehydrated in an ethanol/water series, and incubated overnight with 200 μL of lysis buffer [1% SDS, 50 mmol/L Tris-HCl (pH 8), and 10 mmol/L EDTA] on a thermomixer set at 65°C and shaking at 1,000 rpm. Next, samples were put on ice, supplemented with protease and deacetylase inhibitors, and sonicated with Covaris E220 instrument for 5 minutes (PIP105, 5% DF, 200 cbp) in 1 mL AFA Fiber milliTUBEs. Solubilized chromatin (ranging from 450 ng to 3 μg) was precleared, diluted, and immunoprecipitated overnight with 3 μg of H3K27ac antibody (Diagenode, cat. #A1723-0041-D, lot #C15410196) following a low-input protocol; and finally washed, eluted, and DNA purified. Libraries were prepared with the Accel-NGS 2S Plus DNA Library Kit (Swift Biosciences, cat. no. 21024) on an automated next-generation sequencing (NGS) workstation, amplified for 14 cycles, and sequenced on a Next-seq instrument (PE35).
Immunotherapy Trials
We collected 12 cancer patient cohorts with available RNA-seq data, patient survival, and immunotherapy response from published studies, including more than 900 ICB-treated tumors from 61 gastric cancer, 420 melanoma, 35 renal cell carcinoma, and 348 urothelial cancer. For each cohort, we standardized the transcriptome data across patients by variance stabilizing transformation.
Drug Treatment Transcriptomic Data
By searching for transcriptomic studies on drug perturbation from GEO, we manually collected 719 expression profiles, including both microarray and RNA-seq datasets. Microarray profiles were normalized by the normalizeQuantiles function in limma. For each dataset, we then performed differential expression analysis with limma (for microarray) or DESeq2 (for RNA-seq). Finally, for each drug treatment condition, we got the pattern of gene expression change, which served as the input gene expression signature into signature-based analysis for mining immune-modulating drugs.
Testing the Effect of SMAC Mimetic Birinapant on MHC-I Expression and Sensitivity to T Cell–Driven Killing
To test the effect of birinapant on MHC-I and PD-L1 expression, we treated B16F10 or CT26 cells with IFNγ (0 or 1 ng/mL) and/or birinapant (0, 1, 5, 10, or 20 μmol/L) for 48 hours and tested the effects on MHC-I or PD-L1 expression by flow cytometry.
To test the effect of birinapant on sensitivity of cancer cells to T cell–driven cytotoxicity, we started treatment of B16F10 cells (parental, sgRosa26, or sgTraf3) with vehicle control or 1 μmol/L birinapant 24 hours prior to the coculture and maintained the same condition throughout the coculture.
Testing the combination efficacy of SMAC mimetic and ICB in vivo is detailed in the section “Mouse Experiments” above.
Correlation Analysis
We took the average log2TPM of HLA-A, HLA-B, HLA-C, and B2M as MHC-I expression, and then tested the correlation between MHC-I and PD-L1 expression by using Pearson's correlation adjusted by tumor purity estimated by IHC or the CHAT algorithm (83) in 25 cancer types in TCGA. MHC-I and PD-L1 show strong positive correlation in the majority of cancer types, where they also show strong correlation with CD8 T-cell infiltration. We also computed the correlation between MHC-I and PD-L1 expression in cancer cell lines, which shows the similar results.
We calculated the Pearson correlation of TRAF3-KO signature score with MHC-I expression and CTL level (average expression of CD8A, CD8B, GZMA, GZMB, and PRF1) in multiple tumor cohorts. The TRAF3-KO signature score shows notable positive correlation with MHC-I expression and CTL level in most of the TCGA cancer types and immunotherapy trials.
Survival Analysis
We used the Cox proportional hazards regression analysis to test the association between MHC-I expression and patient survival and selected 16 cancer types, in which the MHC-I expression is positively correlated with patient survival (coxph z-statistics < –0.5), for downstream analysis. To determine the decoupled effect of MHC-I and PD-L1 on patient survival, we tested the association between MHC-I/PD-L1 differential expression and patient survival, and found patients with higher MHC-I/PD-L1 differential expression have longer overall survival in all the 16 cancer types.
The clinical relevance of TRAF3 regulating MHC-I expression is confirmed by testing the association between TRAF3-KO signature and overall survival or progressive survival of patients in TCGA and immunotherapy trials with cox regression.
Gene Set Enrichment/Overrepresentation Analysis
All pathway/GO enrichment analyses were performed using the EnrichAnalyze function in MAGeCKFlute R package, and the enriched pathways are visualized using functions from MAGeCKFlute.
Cistrome-GO and Toolkit for CistromeDB Analysis
We selected top 1,000 enriched peaks from ATAC-seq or H3K27ac ChIP-seq and entered the BED file into Cistrome-GO website to perform pathway enrichment analysis and Toolkit to identify the TFs that have a significant binding overlap with the input peak sets.
Data Availability
All sequencing data were uploaded to GEO with accession number GSE149826.
Code Availability
All codes are available upon request.
Authors' Disclosures
S.S. Gu reports grants from Sara Elizabeth O'Brien Trust during the conduct of the study. A. Lako reports personal fees from Bristol Myers Squibb outside the submitted work. A.K. Tewari reports grants from NIH during the conduct of the study. A. Das Sahu reports personal fees from C-reveal Therapeutics and Checkmate Pharma, and nonfinancial support from Lead Pharma outside the submitted work. G.J. Freeman reports grants from NCI during the conduct of the study; and personal fees from Roche, Bristol-Myers Squibb, Xios, Origimed, Triursus, iTeos, NextPoint, IgM, Jubilant, GV20, and Trillium outside the submitted work; in addition, G.J. Freeman has a patent for PD-L1/PD-1 pathway issued, licensed, and with royalties paid from Roche, Merck MSD, Bristol-Myers Squibb, Merck KGA, Boehringer-Ingelheim, AstraZeneca, Dako, Mayo Clinic, and Novartis; and G.J. Freeman has equity in Nextpoint, Triursus, Xios, iTeos, IgM, GV20, and Geode. S. Rodig has research grants from Merck, Bristol-Myers Squibb, Affimed, and KITE/Gilead for work outside this investigation. F.S. Hodi reports grants and personal fees from Bristol Myers Squibb and Novartis, and personal fees from Merck, EMD Serono, Sanofi, Aduro, Apricity, Surface, Compass Therapeutics, Pionyr, Torque, Rheos, Bicara, Psioxus, Pieris, Eisai, Checkpoint Therapeutics, Idera, Takeda, Genentech, BioEntre, Gossamer, and Iovance outside the submitted work; in addition, F.S. Hodi has a patent for Methods for Treating MICA-Related Disorders pending and with royalties paid, a patent for Tumor antigens and uses thereof issued, a patent for Angiopoietin-2 Biomarker pending, a patent for Compositions and Methods for Identification, Assessment, prevention, and treatment of Melanoma using PD-1 Isoforms pending, a patent for Therapeutic Peptides pending, a patent for Methods of using pembrolizumab and trebananib pending, a patent for Vaccine compositions and methods for restoring NKG2D pathway function against cancers pending and with royalties paid, a patent for Antibodies that bind to MHC class I polypeptide-related sequence A with royalties paid, and a patent for Anti-Galectin antibody biomarker predictive of antiimmune checkpoint and antiangiogenesis responses pending. M. Brown reports personal fees from GV20 Oncotherapy, Kronos Bio, and H3 Biomedicine, and grants and personal fees from Novartis outside the submitted work. X.S. Liu reports grants from NIH during the conduct of the study; and X.S. Liu is a cofounder, board member, and consultant of GV20 Oncotherapy and its subsidiaries, SAB of 3DMedCare, consultant for Genentech, and stockholder of AMGN, JNJ, MRK, PFE, and receives sponsored research funding from Takeda and Sanofi outside the submitted work. No disclosures were reported by the other authors.
Authors' Contributions
S.S. Gu: Conceptualization, formal analysis, validation, investigation, methodology, writing–original draft, writing–review and editing. W. Zhang: Conceptualization, data curation, formal analysis, investigation, methodology, writing–original draft, writing–review and editing. X. Wang: Formal analysis, validation, investigation, methodology, writing–review and editing. P. Jiang: Conceptualization, formal analysis, investigation, methodology, writing–review and editing. N. Traugh: Formal analysis, validation, investigation. Z. Li: Formal analysis, validation, investigation. C. Meyer: Conceptualization, methodology. B. Stewig: Validation, investigation. Y. Xie: Formal analysis. X. Bu: Investigation. M.P. Manos: Investigation.A. Font-Tello: Investigation. E. Gjini: Investigation. A. Lako: Investigation. K. Lim: Investigation. J. Conway: Investigation. A.K. Tewari:Investigation. Z. Zeng: Investigation. A. Das Sahu: Investigation.C. Tokheim: Investigation. J.L. Weirather: Investigation. J. Fu: Investigation. Y. Zhang: Investigation. B. Kroger: Validation, investigation. J.H. Liang: Investigation. P. Cejas: Investigation. G.J. Freeman: Supervision, writing–review and editing. S. Rodig: Resources, investigation. H.W. Long: Resources, investigation. B.E. Gewurz: Supervision, investigation, writing–review and editing. F.S. Hodi: Resources, supervision, investigation. M. Brown: Conceptualization, resources, supervision, writing–review and editing. X.S. Liu: Conceptualization, resources, supervision, writing–review and editing.
Acknowledgments
This study was supported by grants from the NIH (R01CA234018 to X.S. Liu, R01AI137337 to B.E. Gewurz, and P50CA101942-12 and P50CA206963 to G.J. Freeman), Breast Cancer Research Foundation (BCRF-19-100 to X.S. Liu), Burroughs Wellcome Career Award in Medical Sciences (to B.E. Gewurz), and Sara Elizabeth O'Brien Trust Fellowship (to S.S. Gu). The authors thank Drs. Kai Wucherpfennig and Deng Pan for their insightful suggestions on this study.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.