BackgroundSET domain-containing histone lysine methyltransferases (HKMTs) and JmjC domain-containing histone demethylases (JHDMs) are essential for maintaining dynamic changes in histone methylation across parasite development and infection. However, information on the HKMTs and JHDMs in human pathogenic piroplasms, such as Babesia duncani and Babesia microti, and in veterinary important pathogens, including Babesia bigemina, Babesia bovis, Theileria annulata and Theileria parva, is limited.ResultsA total of 38 putative KMTs and eight JHDMs were identified using a comparative genomics approach. Phylogenetic analysis revealed that the putative KMTs can be divided into eight subgroups, while the JHDMs belong to the JARID subfamily, except for BdJmjC1 (BdWA1_000016) and TpJmjC1 (Tp Muguga_02g00471) which cluster with JmjC domain only subfamily members. The motifs of SET and JmjC domains are highly conserved among piroplasm species. Interspecies collinearity analysis provided insight into the evolutionary duplication events of some SET domain and JmjC domain gene families. Moreover, relative gene expression analysis by RT‒qPCR demonstrated that the putative KMT and JHDM gene families were differentially expressed in different intraerythrocytic developmental stages of B. duncani, suggesting their role in Apicomplexa parasite development.ConclusionsOur study provides a theoretical foundation and guidance for understanding the basic characteristics of several important piroplasm KMT and JHDM families and their biological roles in parasite differentiation.