Comprehensive analysis of circRNAs for N7-methylguanosine methylation modification in human oral squamous cell carcinoma

Abstract

N7-methylguanosine (m7G) modification is closely related to the occurrence of tumors. However, the m7G modification of circRNAs in oral squamous cell carcinoma (OSCC) remains to be investigated. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) was used to measure the methylation levels of m7G and identify m7G sites in circRNAs in human OSCC and normal tissues. The host genes of differentially methylated and differentially expressed circRNAs were analyzed by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, and circRNA–miRNA–mRNA networks were predicted using the miRanda and miRDB databases. The analysis identified 2348 m7G peaks in 624 circRNAs in OSCC tissues. In addition, the source of m7G-methylated circRNAs in OSCC was mainly the sense overlap region compared with normal tissues. The most conserved m7G motif in OSCC tissues was CCUGU, whereas the most conserved motif in normal tissues was RCCUG (R = G/A). Importantly, GO enrichment and KEGG pathway analysis showed that the host genes of differentially methylated and differentially expressed circRNAs were involved in many cellular biological functions. Furthermore, the significantly differentially expressed circRNAs were analyzed to predict the circRNA–miRNA–mRNA networks. This study revealed the whole profile of circRNAs of differential m7G methylation in OSCC and suggests that m7G-modified circRNAs may impact the development of OSCC.