YTHDF1 is associated with reduced IFN-γ-related gene signatures and poor prognosis in CRC

By analysing copy number variations of 21 known m⁶A regulators using TCGA data, we identified that YTHDF1 is the top m⁶A regulator displaying copy number gain or amplification in nearly 80% of CRC tumours (online supplemental figure S1A). Such an increase of copy number is also concordant with upregulation of mRNA expression (online supplemental figure S1A) and protein expression (online supplemental figure S1C), implicating that YTHDF1 functions in promoting CRC. To establish a link between antitumour immunity and m⁶A regulators, we then performed GSEA to analyse the correlation between m⁶A regulators with IFN-γ response gene signature. We found that YTHDF1 showed a striking negative correlation (q<0.0001) with the IFN-γ response pathway (online supplemental figure S1B). Of note, IFN-γ response is associated with the induction of antitumour immunity10 11 and responsiveness to ICB therapy.12 13 Consistently, YTHDF1 expression strongly anticorrelated with an 18 IFN-γ-related gene signature (figure 1A), which predicts anti-PD1 responsiveness in multiple cancer types.14 15 Furthermore, expression of CD8A or CD8⁺ T cell signature16 is negatively correlated with YTHDF1 expression (online supplemental figure S1D). In agreement with TCGA data, IHC staining of TMA from our CRC TMA cohorts showed that high protein expression of YTHDF1 was correlated with low infiltration of CD8⁺ T cells in cohort I (p<0.001, r=−0.248, n=206) (figure 1B) and cohort II (p<0.0001, r=−0.269, n=202) (figure 1C). These data strongly imply that m⁶A reader YTHDF1 is associated with impaired antitumour immunity and reduced ICB treatment efficacy.

• Download figure
• Open in new tab
• Download powerpoint

Figure 1

YTH N⁶-methyladenosine RNA binding protein 1 (YTHDF1) is associated with an immunosuppressive microenvironment in patients with colorectal cancer (CRC) and its validation by scRNA-seq in immunocompetent mice. (A) Spearman’s correlation in mRNA expression between interferon gamma (IFN-γ)-related genes and m⁶A regulators in The Cancer Genome Atlas dataset. (B) Correlation between protein levels of YTHDF1 and CD8⁺ T cell infiltration was determined by immunohistochemistry staining in cohort I (p<0.0001, r=−0.2686, n=202), and (C) validation in cohort II (p<0.001, r=−0.2477, n=206). (D) Experimental design for single-cell analysis of CD45⁺ cells sorted from the tumours with or without Ythdf1 by FACS, followed by droplet-enabled scRNA-seq. NC: cells with control sgRNA. Ythdf1-KO: cells with CRISPR knockout of Ythdf1. (E) Representative image (left), tumour volume (middle), and weight (right) of MC38 syngeneic tumours with or without Ythdf1 knockout injected in C57BL/6 mice. (F) Upper panel: tSNE plot showing the components of immune cells in MC38 syngeneic tumours with or without Ythdf1 knockout. Each dot represents a single cell. The same cell type was colour coded. Lower panel: subset analysis of T cell and NK cell clusters shown in tSNE projection regions. (G) Feature plots of characteristic markers of all cell types showing expression levels with low expression in dark red to high expression in bright red. (H) Gating strategies for flow cytometry. Identification of myeloid-derived suppressor cells (MDSC), M-MDSC, and G-MDSC (left). Identification of CD8⁺ T and CD4⁺ T cells (middle). Representative images showing increase of IFN-γ⁺ CD8⁺ T cells in Ythdf1-KO tumours in (D) (Right). (I) Flow cytometry analysis performed with tumours derived from the indicated cells in (D) (**p<0.01; ****p<0.0001) (linear regression (B, C), two tailed t-test (E, H, I), analysis of variance test (E)). NC: cells with control sgRNA. KO1: cells with Ythdf1 sgRNA #1. KO2: cells with Ythdf1 sgRNA #2. GrB: granzyme B.

In agreement with the notion that low CD8⁺ T cell infiltration is associated with poor prognosis,17 high expression of YTHDF1 (54.3%, 100/184) predicted poor survival of patients with CRC (p<0.01, log-rank test) (online supplemental figure S1E,F). Multivariate Cox regression analysis validated that YTHDF1 was an independent prognostic factor for CRC in cohort II (HR, 1.764; 95% CI 1.058 to 2.939; p<0.05) (online supplemental figure S1F). These findings were further verified in cohort I by multivariate Cox regression analysis (online supplemental figure S1E).

Single-cell transcriptomics reveal YTHDF1-induced immunosuppression

To investigate the role of YTHDF1 in regulating antitumour immunity, we used the CRISPR-Cas9 system to knockout Ythdf1 (Ythdf1-KO) in MC38 murine MSI-H CRC cells and injected the cells to syngeneic C57BL6 mice (figure 1D). We found that both tumour volume and weight were reduced by knockout of Ythdf1 compared with control (NC) (figure 1E). To ask if Ythdf1-KO impact TIME, we isolated CD45⁺ immune cells from the tumours and performed single cell RNA-seq (scRNA-seq) (NC: 1480 cells; Ythdf1-KO: 1816 cells). Tumours with Ythdf1-KO exhibited strong reduction of granulocytic myeloid-derived suppressor cells (G-MDSCs, cluster 3) and neutrophils compared with NC group (figure 1F and online supplemental figure S2A). In contrast, T cells and NK cells were largely increased in Ythdf1-KO tumours (figure 1F and online supplemental figure S2A). We further reclustered T and NK cells into CD4⁺ T, CD8⁺ T, NKT and NK cell subsets, and identified that they were simultaneously increased in Ythdf1-KO tumours compared with controls (figure 1F). We thus speculated that YTHDF1 can suppress the antitumour immunity through inducing MDSC accumulation. Functional markers of MDSCs, Il1b, Arg2, Cxcr2 and Ccr2 were examined, and these genes were mainly enriched in MDSC clusters (cluster 3 and cluster 4), especially from tumours without Ythdf1 knockout (figure 1G). Thus, data from this syngeneic model supports an immunosuppressive function of YTHDF1 in CRC.

Ythdf1 knockout reduces MDSCs but increases cytotoxic T cell infiltration

To validate our findings in scRNA-seq analyses, the composition of tumour-infiltrating immune cells in MC38 syngeneic mice was determined by flow cytometry. We confirmed that Ythdf1 knockout significantly suppressed tumour weight and volume (figure 1C), and flow cytometry revealed that Ythdf1 knockout decreased MDSCs, but increased CD8⁺ T and CD4⁺ T cells in tumours (figure 1H,I). Of MDSCs, G-MDSCs were the predominant subset, and knockout of Ythdf1 led to the remarkable reduction of G-MDSCs (figure 1I). In line with the immunosuppressive function of MDSCs, we observed significant increase of functional T cells, including IFN-γ⁺ CD8⁺ T cells, Granzyme B⁺ CD8⁺ T cells, and IFN-γ⁺ CD4⁺ T cells in Ythdf1-KO group (figure 1H,J).

We next performed the experiment in CT26 (MSS-CRC) with Ythdf1-KO to validate the role of YTHDF1 in regulating antitumour immunity. As expected, Ythdf1-KO led to reduced tumour volume and weight (figure 2A,B) in CT26 syngeneic mice, together with reduced functional T cells and accumulation of MDSCs (figure 2C,D). Immunofluorescence staining confirmed the decreased infiltration of MDSCs (CD11b⁺Gr-1⁺) in both MC38 and CT26 syngeneic tumours on Ythdf1 knockout (figure 2E and online supplemental figure S2B). Collectively, Ythdf1 depletion in CRC cells reduced MDSCs and increased functional T cell infiltration. These findings are in line with clinical data demonstrating YTHDF1 anticorrelated with CD8⁺ T cell and IFN-γ-related signatures (figure 1A–C, online supplemental figureS1B,D).

• Download figure
• Open in new tab
• Download powerpoint

Figure 2

YTH N⁶-methyladenosine RNA binding protein 1 (Ythdf1) knockout induces antitumour immunity by reduction of myeloid-derived suppressor cell (MDSC) and increase of functional T cells in syngeneic tumours, an effect reversed by CD8⁺ T cell depletion. (A) Western blot validated Ythdf1 knockout in CT26 cells. (B) Representative image of CT26 syngeneic tumours with or without Ythdf1 Knockout (left). Knockout of Ythdf1 in CT26 cells inhibits tumour growth (middle) and tumour weight (right) in BALB/c mice. (C) MDSC, G-MDSC, M-MDSC, CD4⁺ T cells and CD8⁺ T cells from tumours in (B) were analysed by flow cytometry. (D) Flow cytometry analysis assessing the percentage of T cell functional markers interferon gamma (IFN-γ) and granzyme B (GrB) from tumours in (B). (E) Immunofluorescence identifying MDSCs in subcutaneous tumours from BALB/c injected with the indicated CT26 cells (n=5 each group). (F) MC38 NC or Ythdf1 knockout cells were implanted in C57BL/6 (n=6). Isotype control (IgG) or anti-mouse CD8 antibody (αCD8) were given at 200 µg/mouse on days 4, 7 and 9 post cell injection. Representative images of tumours from each group were shown. (G) Tumour volume (left) and weight (right) from tumours in (F). (H) CT26 NC or Ythdf1 knockout cells were implanted in BALB/c (n=6). Isotype control (IgG) or αCD8 were given at 200 µg/mouse on days 4, 7 and 9 post cell injection. Representative images of tumours from each group were shown. (I) Tumour volume (left) and weight (right) from tumours in (H) (*p<0.05; **p<0.01; ***p<0.001; ****p<0.0001) (two tailed t-test (B, C, D, E, G, H), two-way analysis of variance test (D, G, I)). NC: cells with control sgRNA. KO1: cells with YTHDF1 sgRNA #1. KO2: cells with YTHDF1 sgRNA #2.

We asked if the attenuated tumour formation in Ythdf1-KO was dependent on CD8⁺ T cell antitumour immunity. To address this, we depleted CD8⁺ T cells with an anti-CD8 antibody in the MC38 syngeneic model. Consistent with our hypothesis, depletion of CD8⁺ T cells restored growth of Ythdf1-KO tumours (figure 2F,G), demonstrating that tumour-suppressing function of Ythdf1-KO was dependent, at least partially, on CD8⁺ T cells. This was confirmed in CT26 syngeneic mice showing that anti-CD8 antibody treatment rescued arrested tumour growth in Ythdf1-KO group (figure 2H,I). Depletion of CD8⁺ T cells by anti-CD8 antibody was confirmed by flow cytometry (online supplemental figure S3A,B). Together, Ythdf1 knockout suppresses CRC growth through the induction of CD8⁺ T cell-dependent antitumour immunity.

• Download figure
• Open in new tab
• Download powerpoint

Figure 3

Intestine-specific YTH N⁶-methyladenosine RNA binding protein 1 (Ythdf1) knockin promotes colorectal tumourigenesis and inhibits antitumour immunity in mice. (A) Scheme for azoxymethane-dextran sulphate-sodium (AOM/DSS) mouse model (left). Representative images of the colon at sacrifice (right). (B) Scheme for Apc^Min/+ mouse model (left). Representative images of the colon at sacrifice (right). (C) Tumour number, and tumour size in WT littermates (n=16) and Ythdf1 Ki mice (n=16) mice treated with AOM-DSS (left). Tumour number, and tumour size in Apc^Min/+ Ythdf1^loxp/loxp (n=10) and Apc^Min/+ Ythdf1^loxp/loxpCDX2-cre mice (n=5) mice (right). (D) Infiltration of the indicated immune cells in tumours derived from AOM/DSS mice assessed by flow cytometry. (E) Composition of functional T cells in the tumours derived from AOM/DSS mice assessed by flow cytometry. (F) Infiltration of the indicated immune cells in the tumours derived from Apc^Min/+ mice assessed by flow cytometry. (G) Composition of functional T cells in the tumours derived from Apc^Min/+ mice assessed by flow cytometry. Two tailed t-test (C, E, F, G, H). WT: Ythdf1 wild type; cKI: conditional Ythdf1 knockin.

Intestine-specific Ythdf1 knockin promotes colorectal tumourigenesis and suppresses antitumour immunity in mice

To verify the role of YTHDF1 in spontaneous colorectal tumourigenesis, we generated intestine-specific Ythdf1 knockin mice (Ythdf1^loxp/loxpCDX2-cre) and initiated CRC in these mice by AOM/DSS treatment (figure 3A). We found that the overexpression of Ythdf1 resulted in increased colon tumour number and size in the AOM/DSS model (figure 3C). Flow cytometry revealed increased MDSC infiltration together with reduced NK, CD4⁺ T and CD8⁺ T cells in colon tumours of Ythdf1 knockin mice compared with wildtype mice (figure 3D). Moreover, we found that Ythdf1 knockin reduced the proportion of functional T cells as identified by granzyme B, INF-γ and TNF-α expression (figure 3E).

We next sought to validate these results in Apc^Min/+ -driven spontaneous CRC (figure 3B) by establishing Apc^Min/+Ythdf1^loxp/loxpCDX2-cre mice. Consistently, Apc^Min/+ mice with intestine-specific knockin of Ythdf1 developed significantly bigger colon tumours than their wild-type littermates (figure 3C). Analysis of the tumour-infiltrating immune cells demonstrated remarkably decreased infiltration of NK, CD4⁺ T and CD8⁺ T cells in tumours with Apc^Min/+ with knockin of Ythdf1, together with the induction of MDSCs (figure 3F). Furthermore, Ythdf1 knockin reduced granzyme B⁺, INF-γ⁺, or TNF-α⁺ T cells in Apc^Min/+ mice (figure 3G). Collectively, these results support that YTHDF1 fosters an immunosuppressive microenvironment promoting spontaneous CRC.

Targeting YTHDF1 via VNP-siYTHDF1 boosts antitumour immunity in CD34⁺ humanised mice

To confirm the role of YTHDF1 in modulating human antitumour immune response, we established the CD34⁺ humanised mouse model.18 Mice with peripheral blood mononuclear cells (PBMCs) comprising >20% human CD45⁺ cells were used (figure 4A). To target YTHDF1 in vivo, we developed VNPs19 to carry siRNA against YTHDF1. The humanised NSG mice bearing human CRC HCT116 xenografts were treated with VNP-siNC or -siYTHDF1 after tumours reached 50–100 mm³ (figure 4B). VNP-siYTHDF1 significantly inhibited tumour volume and weight compared with VNP-siNC (figure 4B,C). We also performed flow cytometry to analyse the TIME (figure 4D). VNP-siYTHDF1 decreased MDSC infiltration, but increased CD4⁺ T cell, CD8⁺ T cell and NK cell accumulation (figure 4E). In addition, more IFN-γ⁺, TNF-α⁺ and granzyme B⁺ CD8⁺ T cells were identified in tumours receiving VNP-siYTHDF1 (figure 4E). We also determined the safety of VNP-siYTHDF1 by measuring serum markers of liver (alanine aminotransferase and aspartate transaminase) and kidney function (creatinine and blood urea nitrogen). Mice treated with VNP-siNC or VNP-siYTHDF1 did not demonstrate abnormal liver or kidney function indicators (online supplemental figure S4), showing that VNP treatment is well tolerated. Thus, the targeting of YTHDF1 using VNP-siYTHDF1 is a safe and effective means to enhance antitumour immunity in humanised mice.

• Download figure
• Open in new tab
• Download powerpoint

Figure 4

VNP-siYTHDF1 boosted antitumour immunity in CD34⁺ humanised mice. (A) Workflow for establishment of CD34⁺ humanised mice (left). Percentage of human CD45⁺ cells in the CD34⁺ humanised mice was identified by flow cytometry (right). (B) Design for establishment of HCT116 xenografts and treatment with VNP-siNC or VNP-siYTHDF1 in CD34⁺ humanised mice (left). Representative image of xenografts in different groups (right). (C) Tumour volume (left) and tumour weight (right) in mice in (B). (D) Gating strategies of immune cells in CD34⁺ humanised mice. (E) Infiltration of CD4⁺ T, CD8⁺ T, NK cells, MDSCs and functional CD8 T cells was assessed by flow cytometry in the tumours from (B) (two tailed t-test (C, E), analysis of variance test (C)). VNP, vesicle-like nanoparticles.

YTHDF1 promotes p65 translation to activate TNF/NF-κB signalling

To identify the molecular mechanism by which YTHDF1 elicits immunosuppression, we performed RNA-seq and Ribo-seq in CRC cells with or without knockout of YTHDF1. By RNA-seq analysis, differentially expressed genes between MC38-NC and MC38-YTHDF1-KO cells were enriched in TNF and NF-κB signalling pathways (online supplemental figure S5A). Consistent result was obtained in another CRC cell line CT26 showing that YTHDF1 regulated TNF/NF-κB signalling pathway (online supplemental figure S5B). In support of these, qPCR validated that Ythdf1-KO reduced mRNA expression of TNF/NF-κB targets (online supplemental figure S6A). Furthermore, Ribo-seq data revealed that loss of YTHDF1 was significantly associated with inactivation of TNF signalling (figure 5A). Accordingly, YTHDF1 knockout reduced ribosome protected fragment abundance of genes involved in TNF signalling (figure 5A). Thus, YTHDF1 could regulate TNF/NF-κB signalling by promoting protein translation. Since YTHDF1 functions as a m⁶A reader, we next performed m⁶A immunoprecipitation sequencing (MeRIP-seq) to pinpoint m⁶A-modified transcripts. By screening m⁶A peaks in mRNAs involved in TNF signalling, we identified two m⁶A sites close to the stop codon of p65 mRNA (figure 5B), which was validated by MeRIP-qPCR (figure 5C). Importantly, we identified a direct interaction between YTHDF1 and p65 mRNA by RNA immunoprecipitation (RIP) sequencing and RIP-qPCR with anti-YTHDF1 antibody (figure 5D). Thus, p65 mRNA is a direct target of YTHDF1. We found that knockout of Ythdf1 attenuated p65 protein expression, especially nuclear p65 expression, in both CT26 and MC38 cells, without affecting the expression of other regulators of NF-κB pathways such as IKKα and IκBα (figure 5E). Notably, mRNA expression of p65 was unchanged in Ythdf1-KO cells (online supplemental figure S6B), supporting that YTHDF1 regulates p65 primarily at the level of protein translation, which is in line with the reported YTHDF1 function of facilitating translation of its targets. Consistent results were obtained in human CRC cells showing that overexpression of wild-type YTHDF1, but not dysfunctional mutant,5 6 elevated p65 protein expression; conversely, YTHDF1 knockdown attenuated p65 protein in human CRC cells (figure 5F). RIP-qPCR using anti-YTHDF1 antibody also confirmed a direct interaction between YTHDF1 and p65 mRNA in human CRC cells (figure 5G). We next sought to validate the association of YTHDF1 and p65 in vivo. In Ythdf1 knockin mice, both p65 and phospho-p65 protein were increased in colon tumours (figure 5H). Collectively, YTHDF1 promotes p65 protein expression to activate TNF and NF-κB signalling in vitro and in vivo.

• Download figure
• Open in new tab
• Download powerpoint

Figure 5

YTH N⁶-methyladenosine RNA binding protein 1 (YTHDF1) promotes TNF/NF-kB signalling in CRC through promoting RELA (p65) mRNA translation. (A) Differentially expressed genes between CT26 cells with and without knockout of Ythdf1 were enriched in TNF signalling pathway identifying by Ribo-seq (left). Heatmap of genes on TNF signalling pathway (right). (B) Methylated RNA immunoprecipitation (MeRIP) sequencing on CT26 cells, showing m⁶A modifications on last exon or 3’UTR of RELA (p65) mRNA. (C) Scheme showing the design of primers for MeRIP-qPCR to validate m⁶A modifications on p65 mRNA. Potential m⁶A sites were highlighted in red (upper). MeRIP-qPCR primers are indicated by arrows. MeRIP-qPCR validated m⁶A modification with primers #2 and #3 (lower). (D) RNA immunoprecipitation sequencing with anti-YTHDF1 antibody (YTHDF1-RIP) showing the enrichment of p65 mRNA compared with IgG control (left). YTHDF1-RIP-qPCR with primers specific to murine p65 mRNA (right). (E) Western blot of key effectors of TNF/NF-κB signalling pathway on knockout of Ythdf1. (F) Western blot analysis performed with indicated human CRC cells overexpressing wild-type (WT) or mutant (mut) YTHDF1, or YTHDF1 knockdown with shRNA. Non-targeted shRNA (shNC) were used as control. (G) YTHDF1-RIP-qPCR with primers specific to human p65 mRNA. (H) Expression of p65 and phospho-p65 was determined by western blot in colon tumours from the intestine-specific Ythdf1 knockin mice and wildtype littermates. cKI: conditional knockin (two tailed t-test (C, D, G)).

• Download figure
• Open in new tab
• Download powerpoint

Figure 6

Loss of YTH N⁶-methyladenosine RNA binding protein 1 (Ythdf1) promotes reduction of myeloid-derived suppressor cells (MDSCs) by decreasing CXCL1 secretion. (A) Multiplex mouse cytokine immunoassay assessing levels of 23 cytokines in cell culture supernatants from MC38 and CT26 cell culture, and tumour lysates and sera from MC38 and CT26 syngeneic tumour models. (B) ELISA validation of levels of CXCL1 in the indicated materials as in (A) (n=4 each group). (C) MDSC migration was assessed in vitro with transwell migration assay using conditional medium derived from the indicated cells. CXCR2 inhibitor SB265610 was used at 5 µM. (D) RT-qPCR determination of Cxcl1 mRNA expression in the indicated cells. (E) Flow chart of T cell suppression assay. (F) Representative images of flow cytometry assessing the proliferation of CFSE-labelled T cells cocultured with MDSCs isolated as (E). (G) Quantification of (F) (two tailed t-test (C, D, E, G)).

YTHDF1 promotes MDSC migration via p65-CXCL1 axis

To understand the link between YTHDF1-induced p65 and its immunosuppression, we performed a cytokine multiplex immunoassay for the detection of 23 different mouse cytokines in conditioned medium of CRC cells, tumour lysates, and serum from mice bearing syngeneic tumours. Among these cytokines, CXCL1 was consistently reduced by Ythdf1-KO (figure 6A), and this reduction of CXCL1 was confirmed by ELISA assay (figure 6B). CXCL1 has been reported as a transcriptional target of NF-κB signalling,20 21 and it promotes MDSC chemotaxis via interaction with its receptor CXCR2.22–24 Given that Ythdf1-KO reduced MDSC infiltration in CRC TIME, we asked whether YTHDF1 modulates MDSC migration. As such, in vitro MDSC migration assay was conducted. We found that conditioned medium from wild-type CRC cells enhanced MDSC migration, which was impaired by the knockout of Ythdf1 (figure 6C). Blocking CXCL1-CXCR2 interaction by CXCR2 inhibitor SB265610 eliminated the difference between control and Ythdf1-KO culture supernatant in mediating MDSC migration (figure 6C). Thus, YTHDF1 promotes MDSC migration through CXCL1/CXCR2 axis. We next asked whether YTHDF1 could regulate Cxcl1 mRNA expression. As expected, Ythdf1 knockout significantly inhibited Cxcl1 mRNA levels in murine (figure 6D) and human CRC cell lines (online supplemental figure S7A). Conversely, the overexpression of wild-type YTHDF1, but not its dysfunctional mutant, elevated CXCL1 transcription and protein levels in human CRC cells (online supplemental figure S7A). NF-κB activators, such as TNF-α and IL-1β, have been reported to induce CXCL1 expression.25 We thus treated MC38 cells with TNF-α and measured CXCL1 mRNA and secretion. TNF-α stimulated Cxcl1 mRNA and secretion, an effect abrogated by the loss of Ythdf1 (online supplemental figure S7B). To confirm the link between YTHDF1 and CXCL1 in human CRC, we examined the association between YTHDF1 and CXCL1 expression in TCGA CRC cohort. Consistently, YTHDF1-high CRC demonstrated high CXCL1 expression (online supplemental figure S7C). Besides CXCL1, CXCL2 also positively correlated with YTHDF1 (online supplemental figure S7C). In contrast, CXCL5 and CXCL8 expression negatively correlated with YTHDF1 (online supplemental figure S7C). Considering the role of YTHDF1 in promoting MSDC infiltration, we next performed CRC TMA of CD33, a marker for G-MDSCs.26 27 We found that YTHDF1 protein levels positively correlated with the proportion of intratumoural CD33⁺ cells (online supplemental figure S8). Our findings thus support a YTHDF1-p65-CXCL1/CXCR2 axis mediating MDSC migration in CRC.

• Download figure
• Open in new tab
• Download powerpoint

Figure 7

Targeting YTH N⁶-methyladenosine RNA binding protein 1 (YTHDF1) augmented anti-PD1 blockade therapy in both microsatellite instability-high (MSI-H) and microsatellite stable colorectal cancer (CRC). (A) MC38 cells with or without Ythdf1 knockout were implanted into C57BL/6 mice and treated with IgG or anti-PD1 (αPD1; n=5 per group). Tumour volume (left) and mice survival (right) were presented. Euthanasia was applied when tumour burden was greater than or equal to 1500 mm³ or when the mouse was moribund. (B) Growth curve of subcutaneous tumours from C57BL/6 injected with MC38 and treated with the indicated treatments. Vesicle-like nanoparticle (VNP)-siRNAs were given intratumorally. And αPD1 was given with intraperitoneal injection. (C) Images and weight of the tumours from (B). (D) The indicated CT26 cells were subcutaneously injected in BALB/c mice. The mice were treated with control or αPD1 (n=8 each group). Tumours removed from the mice at sacrifice were shown. (E) Tumour growth curve (left) and tumour weight (right) in (D). (F) Flow cytometry performed with the tumours in (D). (G) Top panel: flow cytometry assessing the proportion of MDSCs, CD4⁺ T cells, and CD8⁺ T cells in total CD45⁺ cells in tumours in (B). Lower panel: flow cytometry assessing functional CD8⁺ T cells in tumours in (B). (H) Schematic of the m⁶A-YTHDF1-p65-CXCL1 axis inducing MDSC accumulation and hence suppressing T cell function in CRC (two tailed t-test (C, D, F, G, H); analysis of variance test (A, C, F); long-rank test (A)). siY1: siYTHDF1. GrB: granzyme B.

We next investigated if YTHDF1 affects MDSC function. CD11b⁺Gr-1⁺ MDSCs were isolated from MC38 syngeneic tumours and then cocultured with T cells in vitro (figure 6E). MDSCs isolated from control tumours suppressed T cell proliferation; however, MDSCs from Ythdf1-KO tumours exhibited significantly less suppressive activity against the proliferation of both CD8⁺ T cells and CD4⁺ T cells compared with MDSCs from the control (figure 6F,G). In support of this, MDSCs isolated from Ythdf1-KO tumours have reduced expression of MDSC functional markers Nos2, Arg1, and Cd274 compared with Ythdf1-WT tumours (online supplemental figure S9A). Flow cytometry confirmed that iNOS⁺ MDSCs was attenuated Ythdf1-KO tumours (online supplemental figure S9B). These data validate the reported immunosuppressive function of MDSCs on key effector cells, including CD8⁺ T cells and CD4⁺ T cells,28 and support that YTHDF1-expressing CRC recruits functional MDSCs.

YTHDF1 is a potential therapeutic target for CRC immunotherapy

Considering that reduced infiltration of MDSCs has been reported to correlate with enhanced immunotherapeutic efficacy in various cancer types,22 29 30 we sought to test if targeting YTHDF1 augments anti-PD1 therapy in CRC. As expected, we found that knockout of Ythdf1 boosted efficacy of anti-PD1 in MC38 (MSI-H) syngeneic tumours and prolonged survival of the tumour-bearing mice (figure 7A). We further exploited the VNPs system to deliver specific Ythdf1-siRNA into tumours. When MC38 syngeneic tumours reached 50~100 mm³, we treated the mice with VNP-siYthdf1 (or VNP-siNC) and anti-PD1 treatment (or IgG). VNP-siYthdf1 significantly suppressed MC38 tumour growth compared with VNP-siNC (figure 7B,C). Remarkably, the combination of VNP-siYthdf1 plus anti-PD1 exerted the strongest inhibitory effect against tumour growth (figure 7B,C).

We further addressed if targeting YTHDF1 can overcome anti-PD1 resistance in MSS CRC based on syngeneic CT26 (MSS CRC) tumour model. CT26 cells with Ythdf1 knockout were thus injected to the syngeneic mice and treated with anti-PD1. We found that knockout of Ythdf1 significantly enhanced anti-PD1 treatment efficacy in CT26 syngeneic tumours, which otherwise were non-responsive to ICB therapy (figure 7D,E).

Flow cytometry analysis further revealed that combination of Ythdf1 silencing and anti-PD1 remarkably increased tumour-infiltrating functional CD8⁺ T cells, including IFN-γ⁺ CD8⁺ T cells and Granzyme B⁺ CD8⁺ T cells in both CT26 and MC38 syngeneic models (figure 7F,G). Furthermore, combinational treatment strongly reduced accumulation of MDSCs, whereas CD4⁺ T cells and CD8⁺ T cells were induced (figure 7G). Thus, targeting of YTHDF1 not only potentiates ICB therapeutic efficacy in MSI-H CRC, but also overcomes ICB resistance in MSS CRC by suppressing recruitment of MDSCs and improving functionality of CD8⁺ T cells.