Emergence of Babesia naoakii infection in Indonesian domestic cattle, a new host record in water buffaloes, and characterization of complete mitochondrial protein-coding genes.

Abstract

Babesia (B.) naoakii, previously referred to as Babesia sp. Mymensingh, is a recently characterized tick-borne haemoprotozoan parasite of cattle. In Indonesia, we first reported its presence in 2022 from clinically affected cattle in Central Java. To investigate the wider epidemiology of this neglected ruminant-associated Babesia species, we surveyed apparently healthy cattle (Bos indicus) and water buffaloes (Bubalus bubalis) across three districts of Java, Indonesia. A PCR assay targeting the B. naoakii-specific apical membrane antigen 1 (ama1) gene detected the parasite occurrence in 34.39% of assessed cattle (87/253; 95% CI: 28.80-40.44%) and 30.77% of water buffaloes (12/39; 95% CI: 18.47-46.52%). These results represent the first record of B. naoakii infection in water buffaloes in the country and confirm that the parasite circulates in subclinically infected bovine hosts. To characterise this apicomplexan parasite further at the molecular level, we assembled in full length the three mitochondrial protein-coding genes (PCGs): cytochrome c oxidase subunits 1 (cox1) and 3 (cox3), as well as cytochrome b (cytb). These genes were reconstructed by next-generation sequencing of blood DNA collected during the acute haemolytic-phase of B. naoakii infection, from calves that subsequently succumbed to the disease in the endemic area. Phylogenetic analyses of the concatenated amino-acid sequences of cox1, cox3, and cytb placed the Indonesian isolates within a well-supported monophyletic clade, distinct from all previously characterised ruminant-associated Babesia species and sister to the Babesia bigemina/Babesia ovata lineage. This placement confirmed species identity and reinforced the genetic distinctiveness of B. naoakii in Indonesia. Notably, although B. naoakii circulates in peripheral blood and mirrors the diagnostic behaviour of the mild pathogen B. bigemina, its clinical impact more closely resembles that of the severe pathogenic B. bovis, particularly in young animals. This diagnostic-clinical discordance highlights the need for B. naoakii-specific molecular surveillance and species-level differentiation in regions of co-endemicity. Given the high prevalence in subclinically B. naoakii-infected adults, the documented severity of babesiosis in calves, and the potential for substantial economic losses, broader epidemiological investigations and species-specific control measures for B. naoakii are urgently performed. The same holds true for future epizootiological investigations of underdiagnosed B. naoakii-infections possibly circulating in Indonesian endemic ruminant bovids such as the banteng (Bos javanicus), the lowland anoa (Bubalus depressicornis) and the tamaraw (Bubalus mindorensis).