In this study, Babesia screening was conducted in 55 rodents and 160 tick samples collected from primary forests and an oil palm plantation in Sarawak, Malaysian Borneo. PCR targeting the 18S ribosomal DNA revealed the presence of Babesia spp. DNA detected in two questing male Haemaphysalis shimoga ticks collected from the oil palm plantation. Sequence analysis revealed that both sequences were identical and had 98.6% identity to a Babesia macropus sequence obtained from Eastern grey kangaroos (Macropus giganteus) in Australia. Phylogenetic tree revealed clustering with marsupial-associated Babesia spp. in the Babesia sensu stricto clade. Whether or not H. shimoga is the competent vector and the importance of the Babesia sp. detected in this study warrants more investigation.
Tick-borne diseases have recently been considered a potential emerging public health threat in Malaysia; however, fundamental studies into tick-borne pathogens and microbiome appear limited. In this study, six tick species (Ixodes granulatus, Haemaphysalis hystricis, Haemaphysalis shimoga, Dermacentor compactus, Dermacentor steini and Dermacentor atrosignatus) collected from two primary forests and an oil palm plantation in Sarawak, Malaysian Borneo, were used for microbiome analysis targeting bacterial 16S rDNA using next-generation sequencing (NGS). In addition, bacterial species were further characterized in conventional PCRs to identify potential pathogens. Sequences generated from NGS were first filtered with the Decontam package in R before subsequent microbial diversity analyses. Alpha and beta analyses revealed that the genus Dermacentor had the highest microbial diversity, and H. shimoga significantly differed in microbial composition from other tick species. Alpha and beta diversities were also significantly different between developmental stages of H. shimoga. Furthermore, we observed that some bacterial groups were significantly more abundant in certain tick species and developmental stages of H. shimoga. We tested the relative abundances using pairwise linear discriminant analysis effect size (LEfSe), which also revealed significant microbial composition differences between Borrelia-positive and Borrelia-negative I. granulatus ticks. Finally, pathogenic and potentially pathogenic bacteria circulating in different tick species, such as Rickettsia heilongjiangensis, Ehrlichia sp., Anaplasma sp. and Bartonella spp. were characterized by PCR and sequencing. Moreover, Coxiella and Francisella-like potential symbionts were identified from H. shimoga and D. steini, respectively. More studies are required to unravel the factors associated with the variations observed in this study.