Comprehensive analysis of m6A modifications in oral squamous cell carcinoma by MeRIP sequencing

N6-methyladenosine (m6A) modifications are the most abundant internal modifications of mRNA and have a significant role in various cancers; however, the m6A methylome profile of oral squamous cell carcinoma (OSCC) in the mRNA-wide remains unknown. In this study, we examined the relationship between m6A and OSCC. Four pairs of OSCC and adjacent normal tissues were compared by Methylated RNA immunoprecipitation sequencing (MeRIP-seq). Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Ingenuity Pathway Analysis (IPA) analyses were used to further analyze the MeRIP-seq data. A total of 2,348 different m6A peaks were identified in the OSCC group, including 85 m6A upregulated peaks and 2,263 m6A downregulated peaks. Differentially methylated m6A binding sites were enriched in the coding sequence in proximity to the stop codon of both groups. KEGG analysis revealed genes with upregulated m6A-modified sites in the OSCC group, which were prominently associated with the forkhead box O (FOXO) signaling pathway. Genes containing downregulated m6A-modified sites were significantly correlated with the PI3K/Akt signaling pathway, spliceosome, protein processing in the endoplasmic reticulum, and endocytosis. IPA analysis indicated that several genes with differential methylation peaks form networks with m6A regulators. Overall, this study established the mRNA-wide m6A map for human OSCC and indicated the potential links between OSCC and N6-methyladenosine modification.