Abstract
BACKGROUND: Growing evidence indicates that N(7)-methylguanosine (m(7)G) modification plays critical roles in epigenetic regulation. However, no data regarding m(7)G modification are currently available in Eimeria tenella, a highly virulent species causing coccidiosis in chickens. METHODS: In the present study, we explore the distribution of internal messenger RNA (mRNA) m(7)G modification in sporulated and unsporulated oocysts of E. tenella as well as its potential biological functions during oocyst development using methylated RNA immunoprecipitation sequencing (MeRIP-seq) and mRNA sequencing (mRNA-seq), and the mRNA-seq and MeRIP-seq data were verified by the quantitative reverse transcription polymerase chain reaction (RT-qPCR) and MeRIP-qPCR, respectively. RESULTS: Our data showed that m(7)G peaks were detected throughout the whole mRNA body, and the coding DNA sequence (CDS) region displayed the most methylation modification. Compared with unsporulated oocysts, 7799 hypermethylated peaks and 1945 hypomethylated peaks were identified in sporulated oocysts. Further combined analysis of differentially methylated genes (DMGs) and differentially expressed genes (DEGs) showed that there was a generally positive correlation between m(7)G modification levels and gene transcript abundance. Unsurprisingly, the mRNA-seq and MeRIP-seq data showed good consistency with the results of the RT-qPCR and MeRIP-qPCR, respectively. Gene Ontology (GO) and pathway enrichment analysis of DEGs with altered m(7)G-methylated peaks were involved in diverse biological functions and pathways, including DNA replication, RNA transport, spliceosome, autophagy-yeast, and cAMP signaling pathway. CONCLUSIONS: Altogether, our findings revealed the potential significance of internal m(7)G modification in E. tenella oocysts, providing some directions and clues for later in-depth research.