Differences in N6-methyladenosine (m6A) methylation among the three major clonal lineages of Toxoplasma gondii tachyzoites

Toxoplasma gondii is an important zoonotic parasite which has over 200 genotypes worldwide. N6-methyladenosine (m6A) methylation is a common epigenetic modification in messenger RNAs (mRNAs), and has been implicated in many aspects of mRNA biology. However, little is known about the difference in m6A methylation among different genotypes of T. gondii. In the present study, we employed methylated RNA immunoprecipitation sequencing (MeRIP-seq) technology to identify key genes exhibiting m6A methylation in the three major clonal lineages (Types I, II and III) of T. gondii tachyzoites. A total of 7650, 8359 and 7264 m6A peaks were identified in 5211, 5607 and 4974 genes in tachyzoites of RH strain (Type I), ME49 strain (Type II) and VEG strain (Type III), respectively. By comparing RH vs. ME49, RH vs. VEG, and ME49 vs. VEG, 735, 192 and 615 differentially methylated peaks (DMPs) were identified in 676, 168 and 553 genes, respectively. A combined MeRIP-seq and RNA-seq analysis revealed 172, 41 and 153 differentially methylated genes (DMGs) at both the m6A methylation and transcriptional level. Gene ontology term enrichment analysis of the DMPs identified differences related to Golgi apparatus, plasma membrane, signal transduction, RNA processing and catalytic step 2 spliceosome. KEGG pathway enrichment analysis showed that the DMGs are mainly involved in endocytosis, systemic lupus erythematosus and mTOR signaling pathway. These findings reveal genotype-specific differences in m6A methylation, which provide new resources for further investigations of the role of m6A in the pathobiology of T. gondii.