New insights into developmental biology of Eimeria tenella revealed by comparative analysis of mRNA N6-methyladenosine modification between unsporulated oocysts and sporulated oocysts

Evidence showed that N6-methyladenosine(m^(6)A)modification plays a pivotal role in influencing RNA fate and is strongly associated with cell growth and developmental processes in many species.However,no information regarding m^(6)A modification in Eimeria tenella is currently available.In the present study,we surveyed the transcriptome-wide prevalence of m^(6)A in sporulated oocysts and unsporulated oocysts of E.tenella.Methylated RNA immunoprecipitation sequencing(MeRIP-seq)analysis showed that m^(6)A modification was most abundant in the coding sequences,followed by stop codon.There were 3,903 hypermethylated and 3,178 hypomethylated mRNAs in sporulated oocysts compared with unsporulated oocysts.Further joint analysis suggested that m^(6)A modification of the majority of genes was positively correlated with mRNA expression.The mRNA relative expression and m^(6)A level of the selected genes were confirmed by quantitative reverse transcription PCR(RT-qPCR)and MeRIP-qPCR.GO and KEGG analysis indicated that differentially m^(6)A methylated genes(DMMGs)with significant differences in mRNA expression were closely related to processes such as regulation of gene expression,epigenetic,microtubule,autophagy-other and TOR signaling.Moreover,a total of 96 DMMGs without significant differences in mRNA expression showed significant differences at protein level.GO and pathway enrichment analysis of the 96 genes showed that RNA methylation may be involved in cell biosynthesis and metabolism of E.tenella.We firstly present a map of RNA m^(6)A modification in E.tenella,which provides significant insights into developmental biology of E.tenella.