Comparative genomics reveals unique features of two Babesia motasi subspecies: Babesia motasi lintanensis and Babesia motasi hebeiensis.

Parasites of the Babesia genus are prevalent worldwide and infect a wide diversity of domestic animals and humans. Herein, using Oxford Nanopore Technology and Illumina sequencing technologies, we sequenced two Babesia subspecies, Babesia motasi lintanensis and Babesia motasi hebeiensis. We identified 3,815 one-to-one ortholog genes that are specific to ovine Babesia spp. Phylogenetic analysis reveals that the two B. motasi subspecies form a distinct clade from other piroplasmas. Consistent with their phylogenetic position, comparative genomic analysis reveals that these two ovine Babesia spp. share higher colinearity with Babesia bovis than with Babesia microti. Concerning the speciation date, B. m. lintanensis split from B. m. hebeiensis approximately 17 million years ago. Genes correlated to transcription, translation, protein modification and degradation, as well as differential/specialized gene family expansions in these two subspecies may favor adaptation to vertebrate and tick hosts. The close relationship between B. m. lintanensis and B. m. hebeiensis is underlined by a high degree of genomic synteny. Compositions of most invasion, virulence, development, and gene transcript regulation-related multigene families, including spherical body protein, variant erythrocyte surface antigen, glycosylphosphatidylinositol anchored proteins, and transcription factor Apetala 2 genes, is largely conserved, but in contrast to this conserved situation, we observe major differences in species-specific genes that may be involved in multiple functions in parasite biology. For the first time in Babesia spp., we find abundant fragments of long terminal repeat-retrotransposons in these two species. We provide fundamental information to characterize the genomes of B. m. lintanensis and B. m. hebeiensis, providing insights into the evolution of B. motasi group parasites.