ALKBH5 regulates anti-PD-1 therapy response by modulating lactate and suppressive immune cell accumulation in tumor microenvironment.

Although immune checkpoint blockade (ICB) therapy has revolutionized cancer treatment, many patients do not respond or develop resistance to ICB. N6 -methylation of adenosine (m6A) in RNA regulates many pathophysiological processes. Here, we show that deletion of the m6A demethylase Alkbh5 sensitized tumors to cancer immunotherapy. Alkbh5 has effects on m6A density and splicing events in tumors during ICB. Alkbh5 modulates Mct4/Slc16a3 expression and lactate content of the tumor microenvironment and the composition of tumor-infiltrating Treg and myeloid-derived suppressor cells. Importantly, a small-molecule Alkbh5 inhibitor enhanced the efficacy of cancer immunotherapy. Notably, the ALKBH5 gene mutation and expression status of melanoma patients correlate with their response to immunotherapy. Our results suggest that m6A demethylases in tumor cells contribute to the efficacy of immunotherapy and identify ALKBH5 as a potential therapeutic target to enhance immunotherapy outcome in melanoma, colorectal, and potentially other cancers.

The adaptive immune response is tightly regulated through immune checkpoint pathways that serve to inhibit T cell activation, thereby maintaining self-tolerance and preventing autoimmunity. The two major checkpoints involve interactions between cytotoxic T lymphocyte antigen 4 (CTLA-4) and programmed cell death protein 1 (PD-1) on T cells and their ligands CD80/CD86 and PD-L1, respectively, which are expressed on various immune cells under physiological conditions. However, expression of these proteins on tumor cells inhibits the T cell activation and enables immune evasion and tumor cell survival. The development of antibodies (Abs) and fusion proteins against PD-1, PD-L1, and CTLA-4, which block negative signaling and enhance the T cell response to tumor antigens, has proven to be a breakthrough in the treatment of solid tumors. Nevertheless, such immune checkpoint blockade (ICB) is ineffective against some tumor types, and many patients who initially respond develop resistance and relapse after ICB. Consequently, understanding the mechanisms of tumor sensitivity, evasion, and resistance to ICB is under intense investigation (1). One of the proposed mechanisms for the failure of ICB is ineffective T cell infiltration and activation due to immunosuppressive conditions within the tumor microenvironment (TME). There is thus an urgent need to develop approaches to increase the sensitivity of tumors to ICBs through combination treatment with molecules that convert an immune-suppressive to an immune-active TME.

Epitranscriptomics is an emerging field that seeks to identify and understand chemical modifications in RNA; the enzymes that deposit, remove, and interpret the modifications (writers, erasers, and readers, respectively); and their effects on gene expression via regulation of RNA metabolism, function, and localization (2, 3). N-methyladenosine (m6A) is the most prevalent internal RNA modification in many species, including mammals. In eukaryotic mRNAs, m6A is abundant in 5′UTRs, 3′UTRs, and stop codons (4–6). The m6A modification is catalyzed by a large RNA methyltransferase complex composed of a catalytic subunit METTL3 and its interacting proteins METTL14, a splicing factor (WTAP), a novel protein (KIAA1429), and other as yet unidentified proteins (2, 3). Conversely, removal of m6A is catalyzed by the RNA demethylases FTO and ALKBH5 (7, 8). In addition, FTO demethylates N,2′-O-dimethyladenosine (m6Am) to reduce the stability of target mRNAs and small nuclear RNA (snRNA) biogenesis (9, 10). The m6A RNA reader proteins, YTH domain-containing proteins (e.g., YTHDF1, YTHDF2, and YTHDF3), specifically bind modified RNA and mediate its effects on RNA stability and translation (11, 12).

In addition to the physiological roles of m6A in regulating RNA metabolism in such crucial processes as stem cell differentiation, circadian rhythms, spermatogenesis, and the stress response (2, 13), increasing evidence supports a pathological role for perturbed m6A metabolism in several disease states. For example, recent studies have shown that the m6A status of mRNA is involved in the regulation of T cell homeostasis (14), viral infection (15), and cancer (16–21).

Here, we employed well-established ICB mouse models of melanoma and colorectal carcinoma to investigate the roles of tumor cell intrinsic Alkbh5 and Fto functions in modulating the response to immunotherapy. We found that CRISPR-mediated deletion of Alkbh5 or Fto in the B16 mouse melanoma (22) or CT26 colorectal carcinoma (23–25) cell line had no effect on tumor growth in untreated mice, but Alkbh5 knockout (KO) significantly reduced tumor growth and prolonged mouse survival during immunotherapy. Alkbh5 deficiency altered immune cell infiltration and metabolite composition in the TME. In addition, the efficacy of cancer immunotherapy was enhanced by pharmacological inhibition of Alkbh5. Finally, we show that gene mutation or down-regulation of the ALKBH5 in melanoma patients correlates with a positive response to PD-1 blockade with pembrolizumab or nivolumab. Thus, our results identify a major role for tumor m6A demethylase in controlling the efficacy of immunotherapy and suggest that combination treatment with ALKBH5 inhibitors may be an approach to sensitize immunotherapy or to overcome tumor resistance to ICB.

Results

Deletion of the m6A RNA Demethylase Alkbh5 Enhances the Efficacy of Anti–PD-1 Treatment.

To determine the role of m6A demethylation enzymes in tumor cells in the response to anti–PD-1 therapy, we employed a mouse model using the poorly immunogenic murine melanoma cell line B16 or modestly immunogenic colorectal cancer cell line CT26. In the standard protocol (Fig. 1), B16 cells were deleted of Alkbh5 or Fto by CRISPR/Cas9 editing and subcutaneously injected into wild-type syngeneic C57BL/6 mice, which were then vaccinated on days 1 and 4 with GVAX (26), composed of irradiated B16 cells secreting granulocyte-macrophage colony-stimulating factor (GM-CSF) to induce an antitumor T cell response. The mice were then treated with anti–PD-1 Ab on days 6, 9, and 12 (or as indicated for individual experiments). In the CT26 model, control or KO cells were subcutaneously injected into BALB/c mice, and mice were then treated with anti–PD-1 Ab on days 11, 14, 17, 20, and 23 (Fig. 1). Gene editing was performed with up to four distinct Alkbh5- or Fto-targeting single-guide RNAs (sgRNAs) per gene (or nontargeting control sgRNAs, NTC), and B16 lines with complete deletion were selected for further experiments (). Compared with NTC-B16 tumors, growth of Alkbh5-KO and Fto-KO tumors was significantly reduced by GVAX/anti–PD-1 treatment (Fig. 1 and ) and the survival of Alkbh5- but not Fto-deficient tumor-bearing mice was significantly prolonged (Fig. 1 and ).

Fig. 1.

Deletion of the m6A RNA demethylase Alkbh5 sensitizes tumors to anti–PD-1 immunotherapy and alters immune cell recruitment. (A) Experimental design to investigate the role of m6A RNA methylation in anti–PD-1 therapy. Alkbh5 and Fto were deleted by CRISPR/Cas9 editing of B16 mouse melanoma cells and injected subcutaneously into C57BL/6 wild-type mice (5 × 105 per mouse). Control mice received NTC B16 cells. Because B16 cells are poorly immunogenic, all mice were injected subcutaneously with GVAX (irradiated B16 GM-CSF cells) on days 1 and 4 to elicit an anti-B16 immune response. Anti–PD-1 Ab (200 μg per mouse) was injected intraperitoneally on days 6, 9, and 12 (or as indicated for individual experiments). Similar experiments were performed for CT26 cells. The cells were inoculated in BALB/c mice and mice were treated with PD-1 Ab on days 11, 14, 17, 20, and 23. (B) Growth of NTC and Alkbh5-KO B16 tumors in C57BL/6 mice treated as described in A. Data are the mean ± SEM of the indicated total number of mice per group. For each gene, three B16 CRISPR cell lines with ≥4 mice per line were examined. (C) Kaplan–Meier survival curves for mice injected with NTC and Alkbh5-KO B16 cells and treated with GVAX and PD-1 Ab. NTC: n = 27; Alkbh5-KO: n = 28. Mice were killed and considered “dead” when the tumor size reached 2 cm at the longest axis. (D) Growth of NTC and Alkbh5-KO CT26 tumors in BALB/c mice treated with anti–PD-1 Ab. Data are the mean ± SEM of the indicated total number of mice per group. (E) Kaplan–Meier survival curves for mice injected with NTC and Alkbh5-KO CT26 cells and treated as described for D. NTC: n = 10; Alkbh5-KO: n = 10. Mice were killed and considered “dead” when the tumor size reached 2 cm at the longest axis. (F) FACS quantification of immune cells isolated from B16 NTC, Alkbh5-KO, and Fto-KO tumors as described in A. Tumor-infiltrating cells were analyzed using the gating strategies described in . CD4+ FoxP3+ (Treg), CD45+CD11b+Ly6G+Ly6CloF4/80−MHC-II− (PMN-MDSCs), and CD45+Ly6C−MHC-II+ CD24hi F4/80lo (DCs) were analyzed. Data are presented as the mean ± SEM. Points represent individual mice. See also . *P < 0.05, **P < 0.01, ***P < 0.001.

We then sought to determine whether the effects of Alkbh5 and Fto KO reflect a generalizable phenomenon during cancer immunotherapy. For this purpose, we employed a modestly immunogenic colorectal cancer line CT26, which responds to PD-1 Ab treatment (23–25). Similar to the B16 model, we found the tumor growth of Alkbh5 KO was significantly reduced compared with NTC in CT26 tumors treated with PD-1 Ab. However, Fto-KO tumors did not show significant changes although they grew slower than NTC (Fig. 1 and ). As observed in the B16 model, the survival of Alkbh5-deficient tumor-bearing mice were significantly prolonged in CT26 model (Fig. 1 and ). These data confirmed the role of Alkbh5-KO and Fto-KO tumors in immunotherapy independent of tumor types. Alkbh5 KO showed more dramatic effects than Fto KO in restricting tumor growth and prolonging mouse survival. In addition, there were no significant differences between the growth of NTC, Alkbh5-KO, and Fto-KO B16 cells either in vitro () or in vivo in untreated mice (), indicating that deletion of the m6A demethylases did not intrinsically impair their growth. Taken together, these data demonstrate that Alkbh5 expression is not required for their growth or survival in vitro or in vivo; however, the enzymes play a crucial role in the efficacy of anti–PD-1 therapy.

Deletion of Alkbh5 in Melanoma Cells Alters the Recruitment of Immune Cell Subpopulations during GVAX/Anti–PD-1 Treatment.

To examine the mechanism by which Alkbh5 modulates GVAX/anti–PD-1 therapy, we examined whether Alkbh5 and Fto deletion in tumor cells modulates immune cell recruitment during GVAX/anti–PD-1 therapy by flow cytometric analysis of tumor infiltrates on day 12 (). Compared with NTC B16 tumors, there is no significant difference in total number of tumor infiltrated lymphocytes (CD45+), CD4+, CD8+ cells in Alkbh5- and Fto-deficient mouse tumors, although a trend to higher abundance of granzyme B (GZMB)+ CD8, GZMB+ CD4 T cell, and NK cell numbers in Fto-null mice tumor (). However, the number of infiltrating regulatory T cells (Tregs) and polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs), but not myeloid (M)-MDSCs, was significantly decreased in Alkbh5-KO tumors compared with NTC tumors during GVAX/anti–PD-1 treatment (Fig. 1 and ). Interestingly, dendritic cells (DCs), but not macrophages, were also significantly elevated in Alkhb5-KO tumors compared with NTC tumors (Fig. 1 and ). In contrast, Fto-KO tumors did not show significant changes in MDSC, Tregs, or DC cell populations (Fig. 1 and ). In accordance with these observations, in Tcra-deficient mice, which lack the TCR-α chain and do not develop mature CD4+ and CD8+ T cells, the effects of Alkbh5 KO but not Fto KO on tumor growth were dampened, but not eliminated (Fig. 2 and ), suggesting that the effect of Alkbh5 in regulating GVAX/anti–PD-1 therapy was partially independent of the host T cell response. To verify the decrease in PMN-MDSCs, we performed immunohistochemical staining and found a marked reduction in the accumulation of MDSCs in Alkbh5-KO tumors compared with NTC tumors on day 12 (Fig. 2).

Fig. 2.

Alkbh5 regulates tumor infiltration of Treg and MDSCs and gene expression during GVAX/anti–PD-1 therapy. (A) As described for Fig. 1, except B16 cells were injected into B6.129S2-Tcratm1Mom/J (TCR-α–deficient) mice, which are devoid of mature CD8+ and CD4+ T cells. Data are presented as the mean ± SEM. *P < 0.05; n.s., not significant. (B) Immunohistochemical staining of Ly6G+ PMN-MDSCs in NTC or Alkbh5-KO tumors isolated from mice on day 12. Magnification: 50 µm. (C) Growth of NTC and Alkbh5-KO tumors in mice treated as described in Fig. 1 and additionally injected intraperitoneally with 10 mg/kg of control IgG or Treg-depleting anti-CD25 Ab on day 11. Data are presented as the mean ± SEM. *P < 0.05 vs. NTC control mice. (D) Growth of NTC and Alkbh5-KO tumors in mice treated as described in Fig. 1 and additionally injected intraperitoneally with 10 mg/kg of control IgG or MDSC-depleting anti-mouse Ly6G/Ly6C (Gr-1) Ab on day 10. Data are presented as the mean ± SEM. *P < 0.05, **P < 0.01 vs. NTC control mice. (E and F) GO analysis (E) and heatmap presentation (F) of DEGs in Alkbh5-KO tumors compared with NTC tumors. Genes satisfying the cut-off criteria of P < 0.05 and log fold-change >0.5 or <−0.5 are shown. See also .

Cross-talk between Tregs and other immune cells is an important contributor to tumor-induced immune suppression; for example, MDSCs can induce Treg amplification and decrease DC differentiation in the TME, and Tregs can greatly inhibit cytotoxic T cell function (27). To assess Treg function in GVAX/anti–PD-1 therapy of melanoma, we monitored the effect on tumor growth after injection of a Treg-depleting anti-CD25 Ab on day 11 of treatment (28, 29). We observed that Treg depletion in NTC tumors showed significant decrease in tumor growth (Fig. 2), while Alkbh5-KO tumors, which had lower numbers of Treg cells than NTC tumors (Fig. 1), did not show significant effects on tumor growth (Fig. 2). These data suggest that Treg cells played important roles in the effects of Alkbh5 KO to restrict tumor growth during therapy, since Treg depletion only worked in NTC tumors that had higher Treg cell numbers. Similarly, we also performed MDSC depletion to assess the tumor growth in NTC and Alkbh5-KO tumors during ICB therapy. Our results show that MDSC depletion had an effect similar to Treg depletion while growth kinetics may vary from these deletions in NTC tumors (Fig. 2). Collectively, these data demonstrate that tumor cell expression of Alkbh5 plays an important role in tumor growth by modulating the recruitment of immunosuppressive MDSCs and Tregs during GVAX/anti–PD-1 therapy.

m6A Demethylase Deletion Alters the Tumor Cell Transcriptome during GVAX/Anti–PD-1 Treatment.

To understand the regulatory role of Alkbh5 and Fto in tumor therapy at the molecular level, we performed RNA-sequencing (RNA-seq) to identify differentially expressed genes (DEGs) in NTC B16 tumors compared with Alkbh5-KO or Fto-KO tumors on day 12 of GVAX/anti–PD-1 treatment. Tumors were confirmed to be Alkbh5- or Fto-deficient before RNA-seq analysis (). Gene ontology (GO) analysis showed that the DEGs in Alkbh5-KO tumors were predominantly involved in metabolic processes, apoptosis, cell adhesion, transport, and hypoxia (Fig. 2 and ). Interestingly, however, DEGs in Fto-KO tumors were mostly immune response-associated genes (). Indeed, further analysis of GO pathways and heatmaps revealed that >80% of the DEGs differed between Alkbh5-KO and Fto-KO B16 tumors. Genes most affected by Alkbh5 KO were associated with regulation of tumor cell survival, adhesion, metastasis, and metabolism, such as Ralgps2, Mmp3, Epha4, Adgrg7, Reln, and Mct4/Slc16a3 (Fig. 2), whereas those most affected by Fto KO were associated with IFN-γ and chemokine signaling, including IRF1, IRF9, STAT2, Cxcl9, Ccl5, and Ccr5 (). To confirm this result, we exposed NTC, Alkbh5-KO, and Fto-KO B16 cells to IFN-γ in vitro and analyzed gene expression by qRT-PCR. As shown in , Fto-KO, but not Alkbh5-KO or NTC tumor cells showed increased expression of the IFN-γ pathway targets Pdl1 and Irf1 and the chemokines Cxcl9, Cxcl10, and Ccl5 after IFN-γ stimulation. These results suggest that, during anti–PD-1/GVAX therapy, Alkbh5 expression in B16 melanoma cells predominantly affects cell intrinsic changes and recruitment of immune cells to the TME, while Fto is involved in regulating IFN-γ and inflammatory chemokine pathways.

IFN-γ pathway activation has been shown to be an important indicator of the efficacy of PD-1 blockade in mouse model studies (22), whereas another study of melanoma patients identified associations between anti–PD-1 response and expression of genes involved in mesenchymal transition, inflammatory, wound healing, and angiogenesis, but not the IFN-γ pathway or other gene signatures indicative of sensitivity to ICB (30). Therefore, we analyzed a gene-expression dataset from 38 melanoma patients who did (n = 21) or did not (n = 17) respond to anti–PD-1 therapy, and searched for DEGs that were also identified here as DEGs in B16 tumors with Alkbh5 or Fto KO. This analysis identified 8 genes that were commonly down-regulated in Alkbh5-KO B16 tumors and responder melanoma patients, and 11 genes that were commonly down-regulated in Fto-KO B16 tumors and responder patients (). Fewer genes were commonly up-regulated between these groups (). These results suggest that the down-regulated genes conserved among mouse model and patients receiving PD-1 Ab treatment play important roles in regulating cancer immunotherapy response and are potential target genes of Alkbh5 and Fto.

Alkbh5 Deletion in Melanoma Cells Affects the m6A Epitranscriptome during GVAX/Anti–PD-1 Treatment.

Given the profound importance of m6A in regulating the function of target RNAs and gene expression (31, 32), we next examined how Alkbh5 affected m6A content in RNA by LC-MS/MS of B16 tumors on day 12 of GVAX/anti–PD-1 therapy (33–35). This analysis revealed that levels of m6A were significantly increased in Alkbh5-KO but not in Fto-KO tumors (Fig. 3). We then performed m6A RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-seq) to determine whether the altered gene expression observed in the KO tumors was a consequence of m6A/m6Am demethylation. To obtain the most robust data, we selected only m6A peaks identified by two independent peak calling algorithms and detected in tumors from all biological replicates per group (). In the NTC and Fto-KO B16 tumors, the majority of m6A peaks were detected in the coding sequence (CDS) and the 3′UTR and 5′UTR, which is consistent with previous studies (4, 5, 36). Notably, the density of m6A peaks in intronic regions was substantially higher in Alkhb5-KO tumors compared with NTC tumors during treatment (Fig. 3), and Alkbh5-KO tumors had more unique m6A peaks compared with NTC or Fto-KO tumors (Fig. 3 and ). Analysis of motifs in the m6A peaks showed that the canonical m6A motif DRACH (D = A, G, U; R = A, G; H = A, C, U) was the most common motif in all tumor groups. The putative m6Am motif BCA (B = C, U, or G; A* = methylatable A) was present in other enriched motifs. One motif enriched in Alkbh5-KO tumors contained the SAG core, which is reminiscent of the SRSF binding site motif known to affect gene splicing (Fig. 3 and ). These data suggest that Fto and Alkbh5 deletion had some common and some distinct effects on m6A/m6Am peaks in B16 tumors, which might contribute to the different mechanisms through which the two demethylases influence the efficacy of GVAX/anti–PD-1 therapy.

Fig. 3.

Alkbh5 regulates gene splicing during GVAX/anti–PD-1 immunotherapy. (A) LC-MS/MS quantification of m6A in ribosome-depleted total RNA isolated from NTC, Alkbh5-KO, and Fto-KO tumors. Data are presented as the mean ± SEM fold-change relative to the NTC control in four mice per group. *P < 0.05 vs. NTC control. (B) Genomic location of the conserved m6A peaks identified by MeRIP-seq in B16 tumors from mice treated as described in Fig. 1. Plot shows the proportion of m6A in the CDS, 5′ and 3′ UTRs, introns, transcription start site (TSS), transcription end site (TES), and intergenic regions. (C) Pie charts showing the proportions of common and unique m6A/m6Am peaks of NTC and Alkhb5-KO B16 tumors from mice treated as described in Fig. 1. (D) Top consensus motifs of MeRIP-seq peaks identified by MEME in NTC and Alkbh5-KO B16 tumors from mice treated as described in Fig. 1. (E) Genome browser tracks of NTC and Alkhb5-KO tumors after treatment were shown for Slc16A3/Mct4 with called m6A sites by MeRIP and corresponding inputs. Input was indicated by blue color in each track. Bed files of the called peaks were shown under the MeRIP track of each group. Scale of the peak density was set the same for all of the groups and shown in the corner. (F) MeRIP-qPCR of m6A enrichment in Mct4 gene for both peak 1 and peak 2 regions shown in E. **P < 0.01 vs. NTC control. (G) The density of m6A in the region of 100-nt exon regions from the 5′ splice site (“SS”) and the 3′ SS. The relative m6A peak density of a specific position in NTC and Alkbh5-deficient tumors was calculated as the scaled m6A peak density proportional to the average m6A peak density in the internal exonic regions. (H) Difference of PSI was calculated by MISO as NTC control minus either Alkbh5-KO or Fto-KO tumors. See also .

We next examined whether the down-regulation of the overlapped genes in Alkbh5-KO or Fto-KO tumors (responding better than NTC) and melanoma patients responding to immunotherapy was due to altered levels of m6A (). Five of eight common down-regulated genes had increased m6A peaks in Alkbh5-deficient mouse tumor (shown in red in ). While only 1 of a total of 11 common genes, Mex3d, had elevated m6A levels in Fto-deficient tumors (red in ). m6A peaks in Mex3d, common in both Alkbh5 and Fto down-regulated genes, increased compared to NTC (). Mct4/Slc16a3, found in only Alkbh5 down-regulated genes, had significantly increased m6A density in the Alkbh5-KO tumors compared to NTC (Fig. 3 ).

These results suggest that Alkbh5 KO increased m6A levels and reduced expression of certain genes involved in immunotherapy resistance. The overall levels of m6A in Fto-deficient tumors was not changed; however, it showed increased m6A at some genes, albeit the number of changed genes were much less than in Alkbh5-KO tumors (e.g., ).

m6A Density Is Increased Near Splice Sites and Leads to Aberrant RNA Splicing in Alkbh5-Deficient Tumors.

Although the regulatory role of m6A deposition in splicing is somewhat controversial (36, 37), Alkbh5 has been reported to affect splicing in an m6A demethylase-dependent manner (38). Our MeRIP-seq results showed that unique m6A peaks were more prevalent in Alkbh5-KO tumors compared with NTC or Fto-KO tumors during GVAX/anti–PD-1 treatment, and that one m6A motif enriched in Alkbh5-KO tumors had a sequence similar to the SRSF binding motif (Fig. 3 ) (36). GO analysis of mRNAs with unique m6A peaks in Alkbh5-KO tumors showed enrichment in splicing, cell cycle, and signaling pathway functions (), suggesting that Alkbh5 also regulates gene expression in B16 cells through effects on mRNA splicing. To test this hypothesis, we examined the location of m6A at 5′ or 3′ intron–exon splice junctions by positional assessment. Consistent with a previous study using m6A individual nucleotide-resolution cross-linking and immunoprecipitation (36, 37), we found that m6A deposition increased from both 5′ and 3′ splice sites to the internal exonic regions in NTC control tumors with immunotherapy (Fig. 3). Surprisingly, we found that in Alkbh5-deficient tumors, the m6A densities were elevated at the both 5′ and 3′ splice sites, with a dramatic increase at the proximal region to the 3′ splicing site (Fig. 3). In contrast, m6A deposition at splice sites in Fto-KO tumors was comparable to that in NTC tumors (), suggesting that Alkbh5 plays a role in gene splicing through depositing m6A modifications near the splicing sites.

Changes in m6Am by FTO have been reported to affect snRNA biogenesis and gene splicing (10), and we observed an increase in m6Am/m6A in U1, U2, and U3 snRNAs in Fto-KO tumors compared with NTC tumors (). To investigate this further, we analyzed our RNA-seq data using MISO to detect differences in RNA splicing. Although the global splicing profiles were unaffected by Alkbh5 or Fto deletion, the frequency of spliced-in transcripts (as reflected by the percent spliced-in index, PSI) in a subset of genes was increased by Alkbh5 deletion in tumors analyzed during GVAX/anti–PD-1 treatment (Fig. 3 and ). Categories of gene functions, where the PSI was changed in Alkbh5-KO tumors, included genes involved in important cellular processes, such as transcription, splicing, protein degradation, transport, translation, and cytokine-related pathways ().

To determine whether changes in m6A deposition were linked with mRNA splicing, we next asked whether the m6A density increased in mRNAs with higher spliced-in frequencies (i.e., higher PSI) in Alkbh5-KO compared with NTC tumors. Indeed, mRNA with high PSI due to Alkbh5 KO had higher m6A densities near intron–exon junctions compared with the same mRNAs in NTC tumors; these mRNAs included Usp15, Arid4b, and Eif4a2 (). Among the genes with altered PSI in Alkbh5-KO tumors after immunotherapy, Eif4a2 regulates gene translation, Arid4b regulates gene transcription, and Sema6d, Setd5, and Met regulate vasculature, the expression and secretion of vascular endothelial growth factor, and hepatocyte growth factor, both of which promote MDSC expansion (39–42). Usp15 affects signaling by transforming growth factor-β, which attracts and activates Tregs. Notably, Met and Usp15 are expressed as isoforms that have markedly different functions (43, 44), suggesting that gene-splicing changes may play a role in TME composition and eventually affecting the immunotherapy efficacy. Taken together, these data indicate that Alkbh5 regulates the density of m6A near spice sites in multiple mRNAs with functions potentially important during GVAX/anti–PD-1 therapy.

Alkbh5 Regulates Lactate and Vegfα Accumulation in the TME during GVAX/Anti–PD-1 Treatment.

Our findings above suggest that Alkbh5 KO regulates its targets by changing m6A levels, which leads to decreased gene expression or altered gene splicing. Some of these genes are involved in regulating cytokines (Vegfα and Tgfβ1) or metabolite (lactate) in TME, such as Mct4/Slc16a3, Usp15, Met, and Sema6d (Figs. 2 and 3 and ). Therefore, it is important to examine whether in Alkbh5-KO tumors, cytokines, or metabolites in the TME are altered that consequently modulate tumor infiltrated lymphocyte populations and immunotherapy efficacy (Figs. 1 and 2).

To address these questions, we quantified lactate, Vegfα, and Tgfβ1 concentrations in the tumor interstitial fluid (TIF), which contains proteins, metabolites, and other noncellular substances present in the TME (). Indeed, both the lactate concentration in TIF and the total lactate content in the TME were dramatically reduced in Alkbh5-KO tumors compared with NTC tumors (Fig. 4). Although the Vegfα concentration in TIF was comparable between NTC and Alkbh5-KO tumors, the total Vegfα content in the TME was reduced by Alkbh5 deletion (Fig. 4). In agreement with a previous study, we also found that Vegfα levels were much lower in plasma than in TIF (45), showing that our isolation of TIF was successful (). The lactate and Vegfα levels in plasma did not differ in mice bearing NTC vs. Alkbh5-KO tumors, suggesting that the effect of Alkbh5 deletion on lactate and Vegfα levels was restricted to the TME and was not systemic (). In contrast to lactate and Vegfα, we found that the concentration of Tgfβ1 in TIF was increased by Alkbh5 deletion, whereas the TME content of Tgfβ1 was reduced only in Alkbh5-deficient tumors (). Collectively, these results showed that Alkbh5 expression in melanoma modulates metabolite and cytokine content with the most significant change of lactate in TIF, suggesting another mechanism by which m6A demethylase could modulate the infiltration of immune cells during anti–PD-1/GVAX treatment.

Fig. 4.

Mct4/Slc16a3 is an Alkbh5 target gene and regulates lactate contents, Tregs and MDSC accumulation in the TME. (A) Lactate concentration and total content in TIF isolated from NTC or Alkbh5-KO tumors excised on day 12 from mice treated as described in Fig. 1. (Left) Absolute lactate concentration in TIF; (Right) lactate content per milligram tumor. Data are the presented as the mean ± SEM of five (NTC) or four (Alkbh5-KO) mice. (B) As for A, except Vegfα was analyzed. (C) mRNA decay analysis of Mct4/Slc16a3 in NTC and Alkbh5-KO B16 cells. NTC and Alkbh5-KO B16 cells were treated with actinomycin D (ActD) at concentration of 5 μg/mL and cells were collected for RNA extraction at indicated time points. Three independent experiments were performed and calculated. *P < 0.05. (D) Mct4 protein levels in NTC, Alkbh5-KO, and Alkbh5-KO cells expressing Mct4 (Alkbh5-KO + Mct4) B16 cells by Western blotting. (E) Extracellular lactate concentration in supernatants of NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 B16 cells. **P < 0.01, ***P < 0.001. (F) Growth of NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 B16 tumors in C57BL/6 mice treated as described in A. Data are the mean ± SEM of the indicated total number of mice per group. The mice number for each group was: NTC: n = 8; Alkbh5-KO: n = 8; Alkbh5-KO+Mct4: n = 10. (G) Lactate concentration and total content in TIF isolated from NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 B16 tumors excised on day 12 from mice treated as described in Fig. 1. Data are the presented as the mean ± SEM. Points represent individual mice. *P < 0.05. (H–I) FACS quantification of immune cells isolated from B16 NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 B16 tumors as described in Fig. 1. Treg cells (H) and PMN-MDSC cells (I) were analyzed. Data are presented as the mean ± SEM. Points represent individual mice. *P < 0.05. (J) PCR analysis of alternative splicing of Eif4a2 and Sema6d genes in NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 B16 cells are shown. (K) Growth of NTC, Alkbh5-KO, Alkbh5-KO + wild-type Alkbh5 (Alkbh5 KO+Alkbh5 Wt), Alkbh5 KO + catalytically mutant Alkbh5 (Alkbh5 KO+Alkbh5 Mut) tumors in C57BL/6 mice treated as described in Fig. 1. Data are the mean ± SEM of the indicated total number of mice per group. Growth of individual mice are shown in . See also .

Mct4/Slc16a3, an Alkbh5 Target Gene, Is Involved in Regulating Extracellular Lactate Concentration, Tregs, and MDSC Accumulation in the TME.

As shown above, we found lactate was the most dramatically decreased metabolite in Alkbh5-KO tumors compared with NTC tumors among all of the Alkbh5-related cytokines and metabolites we examined in the TME (Fig. 4 and ). In Alkbh5-KO tumors, Mct4/Slc16a3 mRNA level was decreased and m6A density was increased compared with NTC tumors during anti–PD-1/GVAX treatment (Figs. 2 and 3 ). Mct4 is a key enzyme catalyzing rapid transport across the plasma membrane of lactate. Lactate is the metabolite that directly affects MDSC and Treg recruitment in tumor sites (46, 47). Therefore, we hypothesized Mct4 is an Alkbh5 target gene in regulating lactate concentration and affecting Tregs and MDSC accumulation in TME during the treatment. To test this hypothesis, we first examined Mct4 expression and RNA stability in NTC and Alkbh5-deficient cells and tumors. We found that Mct4 mRNA levels were lower in Alkbh5-KO than in NTC cells in both B16 and CT26 mouse cell lines, as well as in two other human cell lines when compared ALKBH5 knockdown with control cells (). In mouse tumors under anti–PD-1/GVAX treatment, both mRNA and protein levels of Mct4 were decreased in Alkbh5-KO tumors compared with NTC (). Next, we performed an mRNA decay assay to determine Mct4 RNA stability in NTC and Alkbh5-KO cells. Our results showed that Mct4 mRNA stability was reduced in Alkbh5-KO cells compared with NTC in both B16 and CT26 cell lines (Fig. 4 and ). These results strongly suggest that Alkbh5 regulated Mct4 expression by changing its m6A levels and RNA stability.

To further delineate the role of Mct4 in Alkbh5 KO tumors during anti–PD-1/GVAX treatment, we constructed a stable cell line expressing Mct4 in Alkbh5-KO cells and examined the function of Mct4 in Alkbh5-KO cells in vitro and in vivo. First, we validated the cell lines by detecting the both mRNA and protein levels of Mct4, and performed in vitro proliferation assay for NTC, Alkbh5-KO, and Alkhb5-KO+Mct4 cells. Results of these analyses showed that there was no difference in cell proliferation of NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 cells in vitro (Fig. 4 and ). As Mct4 is a key enzyme mediating transport of lactate across the cell membrane, we examined extracellular lactate concentration in NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 B16 cells. As expected, we observed a reduction of lactate concentration in Alkhb5-KO cells, and an increased level of lactate in Alkhb5-KO+Mct4 cells in vitro (Fig. 4). Furthermore, we inoculated these cells to mice and treated them with anti–PD-1/GVAX and monitored tumor growth in vivo as described above (Fig. 1). We observed a significantly reduced tumor growth in Alkhb5-KO tumors but not Alkbh5-KO+Mct4 tumors compared with NTC, albeit Alkbh5-KO+Mct4 tumors also grew slower than NTC (Fig. 4). These results suggested that Mct4 was one of the major targets of Alkbh5 during anti–PD-1/GVAX treatment. Next, we isolated tumors and assayed the lactate concentration and amounts in TIF and found that lactate levels were significantly reduced in Alkbh5-KO but not Alkbh5-KO+Mct4 tumors, which was consistent with in vitro assay (Fig. 4). In accordance with that, flow cytometry analysis of the tumors showed that Tregs and PMN-MDSC populations were significantly decreased in in Alkbh5-KO but not Alkbh5-KO+Mct4 tumors (Fig. 4 ). Altogether, these results show that Mct4 is a key Alkhb5 target gene mediating reduced lactate levels, as well as Tregs and MDSC populations in Alkbh5-KO tumors during the anti–PD-1/GVAX treatment.

Although not significant, we observed a slower tumor growth in Alkbh5-KO+Mct4 than NTC cells-inoculated mice in vivo (Fig. 4). We speculate that other factors or events also play roles in Alkbh5-KO tumor, albeit not as significant as Mct4, during anti–PD-1/GVAX treatment. Therefore, we analyzed several genes whose splicing events were altered in Alkbh5-KO tumors, and compared their pattens in NTC, Alkbh5-KO, and Alkbh5-KO+Mct4 cells. The results showed that gene splicing changes remained the same in Alkbh5-KO and Alkbh5-KO+Mct4 cells, such as Eif4a2 and Sema6d, suggesting that besides Mct4, they may play roles in Alkbh5-KO tumors (Fig. 4 and ).

m6A mRNA Demethylase Activity of Alkbh5 Is Indispensable during GVAX/Anti–PD-1 Treatment.

Since Alkbh5 is an m6A RNA demethylase, we asked whether the m6A demethylase enzymatic activity is essential for the functions of Alkbh5. We constructed stable cell lines by expressing Alkbh5 CRISPR sgRNA-resistant wild-type or H205A/H267A catalytically inactive mutant [conserved enzymatic sites in human ALKBH5, H204H266 (8, 48)] of Alkhb5 in Alkbh5-KO cells and examined Mct4 expression. Our results showed that wild-type but not mutant Alkbh5 could rescue Mct4 mRNA and protein levels (). We also analyzed gene splicing of Eif4a2 and Sema6d genes, which had altered PSI in Alkbh5-KO tumors, in wild-type and mutant Alkbh5-expressed Alkbh5-KO cell lines. These results showed that Eif4a2 gene splicing was rescued in wild-type but not mutant Alkbh5-expressed Alkbh5-KO cells. MeRIP-seq also showed an increased signal of m6A in Alkhb5-KO cells compared with NTC in the exon–intron junction, which were involved in the alternative splicing of Eif4a2. On the other hand, gene splicing of Sema6d was not affected by the enzymatic activity of Alkbh5, and we did not observe m6A peaks around the spliced exons (). These results showed that enzymatic activity of Alkbh5 play important roles in regulating Mct4 RNA and protein expression, as well as certain genes with altered alternative splicing, which was directly affected by m6A and Alkbh5. Furthermore, we performed in vivo tumor growth experiments in mice treated with GVAX/anti–PD-1. As shown in Fig. 4 and , expressing wild-type but not catalytically inactive mutant Alkbh5 in Alkbh5-KO cells abolished the tumor restricting effects of Alkbh5 KO during GVAX/anti–PD-1 treatment. Altogether, these results demonstrate that the catalytic activity of Alkbh5 is indispensable for its effects on in vivo tumor growth during GVAX/anti–PD-1 treatment.

ALKBH5 and MCT4/SLC16A3 Levels in Melanoma Patients Correlate with the Response to Anti–PD-1 Therapy.

Our results thus far strongly suggest that ALKBH5 deletion enhances the efficacy of anti–PD-1 therapy. Therefore, we analyzed the Cancer Genome Atlas (TCGA) database to examine the correlation between expression level of ALKBH5 and survival time in metastatic melanoma patients. Consistent with our findings, low expression of ALKBH5 correlated with better patients’ survival (Fig. 5). Importantly, Treg cell numbers, as indicated by FOXP3/CD45 ratio, were significant lower in patients with less expression of ALKBH5 (Fig. 5). As described above, we found that Mct4/Slc16a3 was an important Alkbh5 target gene during immunotherapy, and Mct4 level decreased in Alkbh5-KO tumors during therapy. Mct4 is a key gene to mediate lactate secretion which led to reduced lactate in TIF contents and suppressive immune cell populations of Alkhb5-KO tumors during immunotherapy (Figs. 2, 3 , and 4 ). Therefore, we examined the gene expression of ALKBH5 and MCT4/SLC16A3 in the TCGA database. Consistent with our mouse tumor data, we found there is a positive correlation between ALKBH5 and MCT4/SLC16A3 expression in melanoma patients (Fig. 5). Consistent with our results (Fig. 2 and ), our analysis did not show any correlation of ALKBH5 expression with IFN pathway genes IRF1 and PDL1 (). We observed a negative correlation of PBRM1 and GZMB, which serves as a positive control for our analysis (49) ().

Fig. 5.

ALKBH5 expression influences the response of melanoma patients to anti–PD-1 therapy. (A) Kaplan–Meier survival rate analysis of TCGA metastasized melanoma patients grouped by ALKBH5 mRNA levels. Patients with follow-up history were included in the analysis; the mean ALKBH5 level for the entire group was used as the cut-off value. ALKBH5 low: n = 196; ALKBH5 high: n = 163. (B) FOXP3/CD45 expression ratio was calculated for metastatic melanoma patients grouped by ALKBH5 mRNA levels; the mean ALKBH5 level for the entire group was used as the cut-off value. ALKBH5 low: n = 196; ALKBH5 high: n = 163. *P < 0.05. (C) Pearson correlation of ALKHB5 and SLC16A3/MCT4 in melanoma patients from the TCGA database (n = 472). (D) Melanoma patients (n = 26) carrying low or high MCT4/SLC16A3 mRNA expression were treated with pembrolizumab or nivolumab anti–PD-1 Ab. Average expression was used as cut-off. The percentage with complete response (CR), partial response (PR), and progressive disease (PD) are shown. Data are from GSE78220. (E) Pearson correlation of ALKHB5 and MCT4/SLC16A3 in melanoma patients treated with pembrolizumab or nivolumab anti–PD-1 Ab (GSE78220). (F) Melanoma patients carrying wild-type (normal) or deleted/mutated ALKHB5 gene were treated with pembrolizumab or nivolumab. complete response (CR), partial response (PR), and progressive disease (PD) are shown (GSE78220). (G) scRNA-seq data presented as t-distributed stochastic neighbor embedding (t-SNE) plots. Cells were from a tumor biopsy collected from a melanoma patient who showed a response to anti–PD-1 therapy. Plots show the distribution of identified cells. (H) ALKBH5 expression in normal keratinocytes/melanocytes and melanoma tumor cells in melanoma patient receiving PD-1 therapy. See also .

MCT4/SLC16A3 was also found in the down-regulated gene list of 26 melanoma patients receiving pembrolizumab or nivolumab treatment (30); we then analyzed the percentage of patient response to PD-1 Ab in low- and high-expressed MCT4/SLC16A3 groups. Melanoma patients with low expression of MCT4/SLC16A3 has much higher complete or partial response rate than the high-expression group (Fig. 5). In the same cohort of melanoma patients receiving pembrolizumab or nivolumab treatment, we also observed a positive correlation of ALKBH5 and MCT4/SLC16A3 expression (Fig. 5). We next determined whether melanoma patients harboring ALKBH5 deletion/mutation were more sensitive to anti–PD-1 therapy than patients carrying wild-type ALKBH5. To this end, we examined the treatment response according to their ALKBH5 mutation and gene-expression status. As shown in Fig. 5, we found that more patients harboring deleted or mutated ALKBH5 achieved complete or partial responses to pembrolizumab or nivolumab therapy than did patients with wild-type ALKBH5.

Next, we performed single-cell RNA-seq (scRNA-seq) on tumor cells obtained from a patient with stage IV melanoma who had responded well to anti–PD-1 therapy. By using scRNA-seq, we were able to examine ALKBH5 expression in the resistant tumor cells in patients receiving PD-1 Ab. We identified 10 cell types in the tumor (Fig. 5), with substantial immune cell infiltration and very few residual melanoma cells, reflecting the response to therapy. We then examined ALKBH5 expression in the tumor cells and found that 16.7% of melanoma cells (16.7%) expressed ALKBH5 compared with only 6.6% of normal keratinocytes and melanocytes surrounding the tumor cells (Fig. 5). Taken together, these results indicate that tumor expression of ALKBH5 might be a predictive biomarker of patient’s survival and response to anti–PD-1 therapy, at least for melanoma patients.

A Small-Molecule Inhibitor of Alkbh5 Enhances the Efficacy of Anti–PD-1 Therapy.

Our results thus far indicate that loss of m6A demethylase Alkbh5, in B16 melanoma cells, potentiates the efficacy of GVAX/anti–PD-1 therapy. To identify clinically relevant pharmacological inhibitors of Alkbh5, we have identified a specific inhibitor of ALKBH5, named ALK-04, by in silico screening of compounds using the X-ray crystal structure of ALKBH5 (PDB ID code 4NRO) and by performing structure–activity relationship studies on a library of synthesized compounds. First, we tested the cytotoxicity of the inhibitor in vitro, and B16 cells proliferation was not significantly affected by inhibitor treatment (Fig. 6). Next, we compared tumor growth of control and inhibitor-treated mice during immunotherapy. Consistent with our previous findings of Alkbh5-KO tumor, mice treated with ALKBH5 inhibitor significantly reduced tumor growth compared to control (Fig. 6 and ). These results confirmed the function of Alkbh5 in restricting the efficacy of immunotherapy and provide a rational for future combinatorial therapy by using an ALKBH5 inhibitor.

Fig. 6.

ALKBH5 inhibitor enhances efficacy of immunotherapy in combination with GVAX and PD-1 Ab. (A) Proliferation assay of B16 cells treated with DMSO control and 10 µM, 30 µM, and 50 µM ALKBH5 inhibitor. (B) Treatment timeline and mouse B16 tumor growth of control and ALKHB5 inhibitor combined with PD-1 and GVAX immunotherapy. *P < 0.05. (C) Proposed model for ALKBH5-mediated regulation of immunotherapy. ALKBH5 influences anti–PD-1 therapy responses by modifying m6A levels and splicing of specific genes. Inhibition of ALKBH5 mRNA demethylation by CRISPR or a small molecule increased m6A on MCT4/SLC16A3, a lactate transporter, which reduced its mRNA levels leading to reduction of lactate in TIFs. Consequently, MDSC and Treg suppressive immune cell populations in the TME are decreased and immunotherapy therapy responses are enhanced. See Results section for details.

Discussion

A major challenge facing the future of ICB for cancer is to understand the mechanisms of resistance to ICB and to develop combination therapies that enhance antitumor immunity and durable responses. Using the poorly immunogenic B16 mouse model of melanoma, which is resistant to ICB, we discovered that genetic inactivation of the demethylases Alkbh5 and Fto in tumor cells rendered them more susceptible to anti–PD-1/GVAX therapy. The possibility that a similar approach could be employed for clinical applications is supported by the finding that Alkbh5 and Fto KO mice are viable (7, 8). This contrasts with m6A methyltransferases, which are known to be essential for embryonic development and stem cell differentiation (50, 51). Notably, a recent study showed that anti–PD-1 blockade responses were enhanced in FTO knockdown tumors (21). We also observed a similar trend with FTO-KO tumors during PD-1 Ab treatment, but it is not as robust as observed for Alkbh5-KO tumors (Fig. 1 and ). Therefore, Alklbh5 has more obvious effects on PD-1 Ab treatment alone or combined with GVAX compared to Fto (Fig. 1). Besides, it seems that the role of FTO in cell proliferation dominates the effects of FTO for in vivo tumor growth from the published report (21), which we did not observe in our experiments (). Overall, our data showed a more dramatic effects of Alkbh5 in regulating immunotherapy compared to Fto, and we further dissected the mechanisms of Alkbh5 in this process.

Tregs and MDSCs are the dominant immunosuppressive cell populations in antitumor immunity (27). In our study, we found that both cell populations were reduced in Alkbh5-KO tumors during GVAX/anti–PD-1 therapy, whereas the abundance of DCs increased. A decrease in tumor infiltration of MDSCs and Tregs was also observed in a mouse model of 4T1 tumors in response to the anti–PD-1/anti–CTLA-4 plus AZA/ENT treatment (23). Importantly, here we propose the link between m6A demethylase ALKBH5 and the altered tumor infiltrated lymphocytes composition during anti–PD-1/GVAX immunotherapy, providing a new target to regulate the mechanism of the TME and modulate of immunotherapy outcomes.

Our results showed that the function of Alkbh5 in regulating the TME and immunotherapy efficacy was not through the IFN-γ pathway, in accordance with the observation of unchanged infiltrated cytotoxic CD8 T cell population in Alkbh5-deficient tumors. Instead, Alkbh5 KO increased the m6A density in its targets and decreased mRNA expression or enhanced percentage of exon splice-in ratios. For example, Mex3d and Mct4/Slc16a3 mRNA expression was reduced in Alkhb5-KO tumors compared with NTC tumors during GVAX/anti–PD-1 therapy. Mex3d is an RNA-binding protein with putative roles in RNA turnover (52), and Mct4/Slc16a3 is important for pH maintenance, lactate secretion, and nonoxidative glucose metabolism in cancer cells (53). Reduced lactate concentration in the TME has been linked to impaired MDSC and Treg expansion and differentiation (46, 47). In this study, we found that Alkbh5 enzymatic activity is indispensable for regulating in vivo tumor growth during GVAX/anti–PD-1 therapy. Mct4/Slc16a3 was one of the major targets of Alkbh5 during this process. Alkbh5-KO B16 tumors displayed reductions in Mct4/Slc16a3 expression, lactate content in TIF, and MDSC and Treg abundance in the TME. Rescue experiments showed that Mct4/Slc16a3 was responsible for regulating lactate concentration and MDSC, Treg accumulation in Alkbh5-KO tumors during the GVAX/anti–PD-1 therapy. In addition, Mct4/Slc16a3 was reported to regulate VEGF expression in tumor cells (54). We also observed a reduction in the TME level of Vegfα in Alkbh5-KO tumors (Fig. 4 and ).

Except for Mct4, we also analyzed several genes with altered PSI in the Mct4-expressing Alkbh5-KO cells. Gene splicing did not change in the rescue cells compared with Alkbh5-KO cells (Fig. 4 and ); these results suggest that gene splicing may play a role independent of Mct4. Previous studies have shown that tumor-specific alternative splicing-derived neoepitopes were related to immunotherapy response (55). We examined the gene-mutation profiles of several of those genes with altered PSI in melanoma patients, and indeed we found that these genes harbored the mutations that affected gene splicing in patients (). The extract role and detailed mechanisms of gene splicing in Alkbh5-KO tumors during GVAX/anti–PD-1 therapy will need further investigations.

In summary, we have uncovered a previously unknown function for tumor-expressed Alkbh5 in regulating metabolite/cytokine content and filtration of immune cells in the TME during GVAX/anti–PD-1 therapy. Alkbh5-mediated alterations in the density of m6A was found to regulate the splicing and expression of mRNAs with potential roles in the control of tumor growth (Fig. 6). These findings highlight the importance of m6A demethylation in regulating the tumor response to immunotherapy and suggest that ALKBH5 could be a potential therapeutic target, alone or in combination with ICB, for cancer.

Materials and Methods

Tumor samples were obtained from a melanoma patient who had been treated with anti–PD-1 Ab. The procedures were approved by the University of California San Diego Institutional Review Board and the patient provided informed consent. Animal studies and procedures were approved by the University of California San Diego Institutional Animal Care and Use Committee. Details of materials regarding cell lines, mouse strains and human tumor specimens, antibodies, and reagents used for our study can be found in . Detailed methods of mouse models and treatments, CRISPR/Cas9-mediated generation of KO cell lines, flow cytometry analysis of tumor-infiltrating immune cells, qRT-PCR and RNA-seq, MeRIP-seq, MeRIP-seq data analysis, alternative splicing and splice junction analysis, scRNA-seq of human melanoma specimens, TIF isolation and analysis, IFN-γ stimulation of melanoma cells in vitro, cell proliferation assay, Western blot analysis, immunohistochemistry, and LC-MS/MS analysis of m6A RNA can also be found in .