Comprehensive analysis of the transcriptome-wide m6A methylome in colorectal cancer by MeRIP sequencing

Accumulating evidence has demonstrated that N6-methyladenosine (m6A) plays important roles in various cancers, making it essential to profile m6A modifications at a transcriptome-wide scale in colorectal cancer (CRC). In the present study, we performed high-throughput sequencing to determine the m6A methylome in CRC. We obtained six pairs of CRC samples and tumour-adjacent normal tissues from Peking University People's Hospital. We used MeRIP-seq to determine that compared to the tumour-adjacent normal tissues, the CRC samples had 1343 dysregulated m6A peaks, and 625 m6A peaks were significantly upregulated and 718 m6A peaks were significantly downregulated. Genes with altered m6A peaks play critical roles in regulating glucose metabolism, RNA metabolism, and cancer stem cells. Furthermore, we identified 297 hypermethylated m6A peaks and 328 hypomethylated m6A peaks in mRNAs through conjoint analyses of MeRIP-seq and RNA-seq data. After analysing these genes with differentially methylated m6A peaks and synchronously differential expression, we identified four genes (WDR72, SPTBN2, MORC2, and PARM1) that were associated with prognosis of colorectal cancer patients by searching The Cancer Genome Atlas (TCGA). Our study suggests that m6A modifications play important roles in tumour progression and survival of CRC patients. The results also indicate that modulating m6A modifications may represent an alternative strategy to predict the survival of cancer patients and interfere with tumour progression in the future.