Fig. 3

SOGA1 is identified as the direct target of METTL16. A Heatmap showing the expression profile of differentially methylated genes after METTL16 knockdown in SW620 cells. B A total of 11 differential genes were classified according to the level of mRNA and m6A peak in SW620 cells with METTL16 knockdown. C 7 candidate target genes of METTL16 came from the intersection of RNA-sequencing and MeRIP-sequencing. D Volcano plot showed the 7 candidate target genes of METTL16. E-F The 7 candidate target genes mRNA expression in SW620 cells with METTL16 knockdown or overexpression were detected by qRT-PCR. G, H The SOGA1 protein expression in SW620 and HCT15 cells with METTL16 knockdown (G) and in SW620 and HCT116 cells with METTL16 overexpression (H) were detected by western blotting. I The m6A peak visualization of m6A-seq in SOGA1 transcripts in SW620 cells with or without METTL16 knockdown was shown. J Relative m6A enrichment of SOGA1 mRNA in SW620 cells with or without METTL16 knockdown was analyzed and normalized to input by using MeRIP-qPCR. K Representative images of SOGA1 IHC staining (left) and scores (right) in CRC tissues and adjacent tissues. L, M SOGA1 expression in the TCGA CRC cohort. N-P Association of SOGA1 mRNA expression with tumor size (N), lymph node metastasis (O), and distant metastasis (P) in CRC patients in TCGA database. R, S IHC (METTL16 and SOGA1)-stained paraffin-embedded sections obtained from xenografts based on SW620 cells with METTL16 knockdown or overexpression. **P < 0.01