Genetic Diversity of Merozoite Surface Antigens in Global Babesia bovis Populations

Cattle can be severely infected with the tick-borne protozoa Babesia bovis, giving rise to serious economic losses. Invasion of the host’s RBCs by the parasite merozoite/sporozoites depends largely on the MSA (merozoite surface antigens) gene family, which comprises various fragments, e.g., MSA-1, MSA-2a1, MSA-2a2, MSA-2b and MSA-2c, highlighting the importance of these antigens as vaccine candidates. However, experimental trials documented the failure of some developed MSA-based vaccines to fully protect animals from B. bovis infection. One reason for this failure may be related to the genetic structure of the parasite. In the present study, all MSA-sequenced B. bovis isolates on the GenBank were collected and subjected to various analyses to evaluate their genetic diversity and population structure. The analyses were conducted on 199 MSA-1, 24 MSA-2a1, 193 MSA-2b and 148 MSA-2c isolates from geographically diverse regions. All these fragments displayed high nucleotide and haplotype diversities, but the MSA-1 was the most hypervariable and had the lowest inter- and intra-population gene flow values. This fragment also displayed a strong positive selection when testing its isolates for the natural selection, which suggests the potential occurrence of more genetic variations. On the contrary, the MSA-2c was the most conserved in comparison to the other fragments, and displayed the highest inter- and intra-population gene flow values, which was evidenced by a significantly negative selection and negative neutrality indices (Fu’s Fs and Tajima’s D). The majority of the MSA-2c tested isolates had two conserved amino acid repeats, and earlier reports have found these repeats to be highly immunogenic, which underlines the importance of this fragment in developing vaccines against B. bovis. Results of the MSA-2a1 analyses were also promising, but many more MSA-2a1 sequenced isolates are required to validating this assumption. The genetic analyses conducted for the MSA-2b fragment displayed borderline values when compared to the other fragments.