Population genetics of Babesia vogeli based on the mitochondrial cytochrome b gene

The current study aimed at population genetic characterization of B. vogeli based on the cytochrome b (cyt b) gene sequences (≥ 685 bp) available in the GenBank. Phylogenetic trees placed all the sequences of B. vogeli in a single large monophyletic clade; however, it was further divided into two subclades (Bv1 and Bv2). Out of seven nucleotide variations observed between Bv1 and Bv2 subclades, four were synonymous (G92A, C170T, T488C and A659G), and three were non-synonymous (G324A, C438A and G465A) resulting in amino acid substitutions at three places (V108I, L146I and V155I). Within different B. vogeli populations, the nucleotide and haplotype diversities were low. The median-joining haplotype network revealed only two haplotypes (Hap_1 and Hap_2). A geographical sub-structuring was noticed in the B. vogeli populations, with moderate genetic differentiation (F_ST = 0.05000; P < 0.05) and a very high gene flow (Nm = 4.75) between Indian and Chinese populations. Neutrality tests and mismatch distributions for the Indian population and the overall dataset of B. vogeli indicated a constant population size. This study provides the first insight into the genetic characterization, population genetics and haplotype network of B. vogeli based on the cyt b gene.