To comprehensively profile the landscape of the mRNA N-methyladenosine (mA) modification in human colorectal cancer (CRC). Methylated RNA immunoprecipitation sequencing (MeRIP-seq) was explored to compare the difference in mRNA N-methyladenosine (mA) methylation between CRC tissues and adjacent normal control (NC) tissue. RNA-sequencing (RNA-seq) was performed to transcribe differentially expressed mRNAs. Conjoint analysis of MeRIP-seq and RNA-seq data was conducted to predict RNA-binding proteins (RBPs). MeRIP-seq identified 1110 differentially mA methylated sites (DMMSs) and 980 differentially mA methylated genes (DMMGs) in CRC, with 50.13% of all modified genes showing unique mA-modified peaks in CRC. RNA-seq showed 915 upregulated genes and 1463 downregulated genes in CRC. QRT-PCR verified the RNA-seq results by detecting the expression of some mRNAs. Conjoint analysis of MeRIP-seq and RNA-seq identified 400 differentially mA methylated and expressed genes (DEGs), and pathway analysis detected that DMMGs and DEGs were closely related to cancer. After analyzing these DMMGs and DEGs through the GEPIA database, we found that the expression of B3GNT6, DKC1, SRPK1, and RIMKLB were associated with prognosis, and the expression of B3GNT6 and RIMKLB were associated with clinical stage. 17 RBPs were identified based on the DMMGs and DEGs, among which FXR1, FXR2, FMR1, IGF2BP2, IGF2BP3, and SRSF1 were obviously highly expressed in CRC, and FMR1, IGF2BP2, and IGF2BP3 were closely related to methylation, and might be involved in the development of CRC. This study comprehensively profiled mA modification of mRNAs in CRC, which revealed possible mechanisms of mA-mediated gene expression regulation.
