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

N6-methyladenosine (m6A) modifications are the most abundant internal modi-fications 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 Meth-ylated 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 upregu-lated 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-meth-yladenosine modification.

Automated summary

This study identified m6A modifications in oral squamous cell carcinoma (OSCC) using MeRIP-seq and revealed potential connections between m6A modifications and OSCC. The findings showed that genes with m6A modifications are involved in multiple signaling pathways in OSCC.