Phylogenetic analysis, genetic diversity and geographical distribution of Babesia caballi based on 18S rRNA gene

Molecular Diagnosis and Identification of Equine Piroplasms: Challenges and Insights from a Study in Northern Italy

Equine piroplasmosis is a tick-borne disease caused by Babesia and Theileria species. Despite its presence in Europe, no laboratory testing is required for animal movement, even though some countries remain free of this disease. Differentiating between species and genotypes is crucial to determine the most effective treatment, as dosage, active compounds, and duration vary. However, diagnosis is often challenging due to genetic variability and the limited sensitivity of molecular methods. The aims of this study were to compare the performances of different molecular diagnostic tests to identify the most effective assay for piroplasm DNA detection and to genetically characterize the piroplasms identified in horses in Northern Italy from 2016 to 2022. Among 63 horses tested, 24 (38.1%) were positive in at least one of the tests used. Four horses tested positive for Babesia caballi with identical nucleotide sequences, while 22 horses tested positive for genetically different Theileria species, including Theileria equi, Theileria haneyi-like species, and Theileria sp. Africa. Two horses were coinfected by Babesia caballi and Theileria haneyi-like species. The best diagnostic approach to avoid false negative results was a combination of different assays. Further studies will be necessary to better assess the prevalence and genetic diversity of piroplasms involved in equine piroplasmosis.

Genetic diversity and molecular evolution of the internal transcribed spacer regions (ITS1-5.8S-ITS2) of Babesia vogeli

Canine babesiosis, a severe haemoparasitic disease caused by Babesia species, has a significant global presence and can be fatal if left untreated. The current study was aimed to perform the population genetic characterization of B. vogeli on the basis of the internal transcribed spacer regions (ITS1-5.8S-ITS2). A maximum likelihood tree constructed with the Hasegawa-Kishino-Yano model grouped all sequences into a single major clade (BvG1), with the exception of a Taiwanese isolate (EF186914), which branched separately. This Taiwanese isolate represented a novel genotype (BvG2) identified in the present study. Nucleotide sequences (n = 62) exhibited 92.5-100 % nucleotide identity among themselves. However, the BvG1 and BvG2 genotypes shared a lower identity of 92.5-93.8 % between them. Notably, the newly generated Indian sequences (n = 21) demonstrated a high degree of homology, with 98.3-100 % identity. Alignment of the nucleotide sequences revealed 58 variations across the dataset. Additionally, 32 sites exhibited variation within the BvG1 genotype, while 56 sites differed between BvG1 and BvG2 genotypes. Within different B. vogeli populations, the nucleotide diversity (π) was low, but the haplotype diversity (Hd) was high. The haplotype diversity of the Indian population, BvG1 genotype, and the combined dataset was ∼0.8 suggesting a high haplotype diversity. The median-joining haplotype network displayed a total of 21 haplotypes, out of which six haplotypes consisted of more than one sequence (2-25 sequences). Haplotype distribution showed significant geographical structuring, with most haplotypes confined to a single country. Only two haplotypes (9.52 %; Hap_1 and Hap_4) were shared between countries, whereas 19 haplotypes (90.48 %) were country-specific. Hap_1, Hap_6, and Hap_4 were the most representative haplotypes, comprising 25, 10, and four sequences, respectively. India exhibited the highest number of haplotypes (h = 13) followed by China (h = 4), the United States of America (h = 3), Taiwan and Tunisia (h = 2), and Thailand (h = 1). Both location-wise and genotype-wise median joining haplotype networks clustered the haplotypes in two groups, representing two distinct genotypes (BvG1 and BvG2). The B. vogeli populations between Thailand and Tunisia exhibited the highest genetic differentiation (FST = 0.80) with a low gene flow (Nm = 0.125) between them. Results of AMOVA revealed a higher genetic variation within populations (69.43 %) as compared to the variation between them (30.57 %). Neutrality indices and the mismatch distributions of the Indian population and the overall dataset of B. vogeli indicated a constant population size to population expansion and population expansion, respectively, with the presence of two distinct genotypes. These data provide information about parasite population genetics and highlight the importance of starting a long-term molecular surveillance program. In conclusion, a high genetic diversity along with the presence of two distinct genotypes of B. vogeli were observed on the basis of internal transcribed spacer regions (ITS1-5.8S-ITS2).

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 (FST = 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.

The first study of the prevalence and genetic diversity of Theileria equi and Babesia caballi in horses in Russia

Equine piroplasmosis (EP) is a global worldwide infection, which can lead to the death of animals. Despite the causative agents of EP being well studied, there are no data on the distribution and genetic characteristics of EP agents in any region of Russia. In this study, blood samples from 750 horses from Novosibirsk province, Irkutsk province, and Altai region of Russian Siberia were examined for the presence of EP agents. Theileria equi and Babesia caballi were detected in all examined regions, with mean prevalence rates of 60.4% and 7.2%, respectively. The identified pathogens were genetically characterized by the 18S rRNA gene. The determined T. equi sequences were highly conserved and belonged to genotypes A and E, with genotype E being found in 88.6% of genotyped samples. In contrast to T. equi, B. caballi sequences were genetically diverse. Seven sequence variants of B. caballi were identified, and only two of them matched known sequences from the GenBank database. The determined B. caballi sequences belonged to four distinct branches within genotype A. Mixed infections with several variants of B. caballi or with T. equi and B. caballi were common. The conducted phylogenetic analysis based on all available B. caballi sequences of the 18S rRNA gene (> 900 bp) from GenBank and from this study first demonstrated the presence of five monophyletic clusters within genotype A and three clusters within genotype B. Thus, the genetic study of B. caballi from Siberia has significantly expanded the data on the genetic diversity of this pathogen.

Revisiting the genotypes of Theileria equi based on the V4 hypervariable region of the 18S rRNA gene

Introduction

Equine theileriosis, an economically important disease that affects horses and other equids worldwide, is caused by a tick-borne intracellular apicomplexan protozoa Theileria equi. Genotyping of T. equi based on the 18S rRNA gene revealed the presence of two, three, four or five genotypes. In previous published reports, these genotypes have been labelled either alphabetically or numerically, and there is no uniformity in naming of these genotypes. The present study was aimed to revisit the phylogeny, genetic diversity and geographical distribution of T. equi based on the nucleotide sequences of the V4 hypervariable region of the 18S rRNA gene available in the nucleotide databases.

Methods

Out of 14792 nucleotide sequences of T. equi available in the GenBank™, only 736 sequences of T. equi containing the complete V4 hypervariable region of the 18S rRNA gene (>207 bp) were used in multiple sequence alignment. Subsequently, a maximum likelihood phylogenetic tree was constructed based on the Kimura 2-parameter model (K2+I).

Results

The phylogenetic tree placed all the sequences into four distinct clades with high bootstrap values which were designated as T. equi clades/ genotypes A, B, C and D. Our results indicated that the genotype B of Nagore et al. and genotype E of Qablan et al. together formed the clade B with a high bootstrap value (95%). Furthermore, all the genotypes probably originated from clade B, which was the most dominant genotype (52.85%) followed by clades A (27.58%), and C (9.78%) and D (9.78%). Genotype C manifested a comparatively higher genetic diversity (91.0-100% identity) followed by genotypes A (93.2-99.5%), and B and D (95.7-100%). The alignment report of the consensus nucleotide sequences of the V4 hypervariable region of the 18S rRNA gene of four T. equi genotypes (A-D) revealed significant variations in one region, between nucleotide positions 113-183, and 41 molecular signatures were recognized. As far as geographical distribution is concerned, genotypes A and C exhibited far-extending geographical distribution involving 31 and 13 countries of the Asian, African, European, North American and South American continents, respectively. On the contrary, the genotypes B and D exemplified limited distribution with confinement to 21 and 12 countries of Asian, African and European continents, respectively. Interestingly, genotypes A and C have been reported from only two continents, viz., North and South America. It was observed that genotypes A and C, and B and D exhibit similar geographical distribution.

Discussion

The present study indicated the presence of only four previously described T. equi genotypes (A, B, C and D) after performing the molecular analyses of all available sequences of the complete V4 hypervariable region of the 18S rRNA gene of T. equi isolates in the GenBank™.

Population genetic characterization of Theileria annulata based on the cytochrome b gene, with genetic insights into buparvaquone susceptibility in Haryana (India)

The present investigation was aimed at population genetic characterization of Theileria annulata on the basis of the cytochrome b (cyt b) gene along with the evaluation of status of buparvaquone resistance in Haryana (India). The sequences originating from China, Egypt, India, Iran, Iraq, Tunisia, Turkey and Sudan were included in the analysis. The maximum likelihood tree based on the Tamura-Nei (TN93+G) model placed all the sequences of T. annulata into a single clade. The median-joining haplotype network exemplified geographical clustering between T. annulata haplotypes originating from each country. Only five haplotypes (7.81 %) were shared between any two countries, while the remaining 59 haplotypes (92.19 %) were singleton and unique to one country. The values of pairwise genetic distance (FST) between all the populations indicated huge genetic differentiation (> 0.25) between different T. annulata populations, barring the FST value between Iraq and Turkey (0.14454) which suggested a moderate differentiation. Contrary to the FST index, the values of gene flow (Nm) between T. annulata populations were very low. The neutrality indices and mismatch distributions indicated a population expansion in the Indian T. annulata population. Furthermore, the secondary structure and homology modeling of the partial cyt b protein is also reported. The molecular analysis of newly generated sequences for buparvaquone resistance revealed that all the isolates were susceptible to buparvaquone treatment. However, two novel mutations at positions V203I and V219I in between the Q01 and Q02 drug-binding regions of the cyt b gene were observed for the first time.

Phylogenetic analysis and haplotype networking of Hepatozoon felis infecting wild animals in Gir National Park, Gujarat, India

The present study was aimed to determine the phylogenetic relationship, haplotype network, and demographic dynamics of H. felis infecting the endangered Asiatic lions in Gir National Park, Gujarat, India, on the basis of partial 18S rRNA gene. The phylogenetic analysis based on the partial 18S rRNA gene sequences of H. felis exhibited the presence of two distinct genotypes of H. felis (HfG1 and HfG2) infecting the Indian wild felids, viz., the Asiatic lion, Royal Bengal tiger, and Indian leopard. The HfG1 and HfG2 genotypes exhibited 97.6-100% and 99.7-100%, and 96.9-98.7% nucleotide identity within and between themselves, respectively. The HfG1 genotype exhibited a higher genetic diversity as compared to HfG2. A total of 22 molecular signatures were identified in the 18S rRNA gene between these genotypes. Further, analysis of a total 67 sequences of H. felis (13 different host species from 13 countries of Africa, South America, Asia, and Europe) that were downloaded from GenBankTM, generated 30 haplotypes. Among all the haplotypes, Hap_17 (h=12) was the most frequent followed by Hap_12 (h=09) and Hap_4 (h=05). Out of 13 location-wise populations, India (h=12) shared the highest number of haplotypes followed by Japan (h=08), and the least number of haplotypes were found in Hungary (h=02). Population dynamics study involving neutrality tests and mismatch distribution, and genetic differentiation indices, revealed the presence of phylogeographic population structure and a constant population size indicating a uniform gene flow among the populations worldwide. In conclusion, a high genetic diversity along with the presence of two distinct genotypes of H. felis were observed on the basis of 18S rRNA gene sequence analysis.

Microfluidic PCR and network analysis reveals complex tick-borne pathogen interactions in the tropics

BACKGROUND

Ixodid ticks, particularly Rhipicephalus sanguineus s.l., are important vectors of various disease-causing agents in dogs and humans in Cuba. However, our understading of interactions among tick-borne pathogens (TBPs) in infected dogs or the vector R. sanguineus s.l. remains limited. This study integrates microfluidic-based high-throughput real-time PCR data, Yule's Q statistic, and network analysis to elucidate pathogen-pathogen interactions in dogs and ticks in tropical western Cuba.

METHODS

A cross-sectional study involving 46 client-owned dogs was conducted. Blood samples were collected from these dogs, and ticks infesting the same dogs were morphologically and molecularly identified. Nucleic acids were extracted from both canine blood and tick samples. Microfluidic-based high-throughput real-time PCR was employed to detect 25 bacterial species, 10 parasite species, 6 bacterial genera, and 4 parasite taxa, as well as to confirm the identity of the collected ticks. Validation was performed through end-point PCR assays and DNA sequencing analysis. Yule's Q statistic and network analysis were used to analyse the associations between different TBP species based on binary presence-absence data.

RESULTS

The study revealed a high prevalence of TBPs in both dogs and R. sanguineus s.l., the only tick species found on the dogs. Hepatozoon canis and Ehrlichia canis were among the most common pathogens detected. Co-infections were observed, notably between E. canis and H. canis. Significant correlations were found between the presence of Anaplasma platys and H. canis in both dogs and ticks. A complex co-occurrence network among haemoparasite species was identified, highlighting potential facilitative and inhibitory roles. Notably, H. canis was found as a highly interconnected node, exhibiting significant positive associations with various taxa, including A. platys, and E. canis, suggesting facilitative interactions among these pathogens. Phylogenetic analysis showed genetic diversity in the detected TBPs.

CONCLUSIONS

Overall, this research enhances our understanding of TBPs in Cuba, providing insights into their prevalence, associations, and genetic diversity, with implications for disease surveillance and management.

Chromosome-level genome assembly of Babesia caballi reveals diversity of multigene families among Babesia species

BACKGROUND

Babesia caballi is an intraerythrocytic parasite from the phylum Apicomplexa, capable of infecting equids and causing equine piroplasmosis. However, since there is limited genome information available on B. caballi, molecular mechanisms involved in host specificity and pathogenicity of this species have not been fully elucidated yet.

RESULTS

Genomic DNA from a B. caballi subclone was purified and sequenced using both Illumina and Nanopore technologies. The resulting assembled sequence consisted of nine contigs with a size of 12.9 Mbp, rendering a total of 5,910 protein-coding genes. The phylogenetic tree of Apicomplexan species was reconstructed using 263 orthologous genes. We identified 481 ves1-like genes and named "ves1c". In contrast, expansion of the major facilitator superfamily (mfs) observed in closely related B. bigemina and B. ovata species was not found in B. caballi. A set of repetitive units containing an open reading frame with a size of 297 bp was also identified.

CONCLUSIONS

We present a chromosome-level genome assembly of B. caballi. Our genomic data may contribute to estimating gene expansion events involving multigene families and exploring the evolution of species from this genus.

Analysis of genetic diversity and population structure of Babesia gibsoni

Babesia gibsoni is a tick-borne apicomplexan protozoan causing canine babesiosis. This parasite has diploid sexual reproduction in ticks, during which genetic exchanges can occur leading to increased genetic diversity, which is an important factor in adapting to environmental changes. Exploring the genetic variation of B. gibsoni population can provide a foundation for understanding the patterns of disease transmission and developing babesiosis control strategies. Partial 18S rRNA fragment sequences were obtained from 11 B. gibsoni isolates collected from different regions in China and 117 publicly available sequences were from 12 geographical areas including China. The genetic variation, demographic expansion and population structure were examined. A total of 34 haplotypes were identified among B. gibsoni populations. Analysis of molecular variance, pairwise Fst and structure analysis showed that high genetic variation within populations, low genetic differentiation and obvious mixture haplotype were apparent in a single continent, but higher genetic differentiation was detected across different continents. Neutrality tests implied that B. gibsoni populations had experienced population extension. These findings will contribute to understand the genetics and evolution of B. gibsoni and will be useful for formulating effective management strategies to prevent and control this parasite.

An insight into misidentification of the small-subunit ribosomal RNA (18S rRNA) gene sequences of Theileria spp. as Theileria annulata

BACKGROUND

There had been isolated reports of the presence of novel Theileria annulata genotypes based on the 18S rRNA gene sequence data from India, Pakistan and Saudi Arabia; but, these studies were restricted to limited field samples. Additionally, no comparative study has been conducted on all the isolates of this parasite from different countries whose sequences are available in the nucleotide databases. Therefore, we aimed to study the genetic diversity of T. annulata based on all available nearly complete 18S rRNA gene sequences in the GenBank™. Out of a total of 312 gene sequences of T. annulata available in the NCBI database, only 70 nearly complete sequences (> 1527 bp) were used for multiple sequence alignment.

RESULTS

The maximum likelihood tree obtained using TN93 + G + I model manifested two major clades. All the valid host-cell transforming Theileria species clustered in one clade. The T. annulata designated sequences occupying this clade clustered together, excluding two isolates (DQ287944 and EU083799), and represented the true T. annulata sequences (n = 54). DQ287944 and EU083799 exhibited close association with Theileria lestoquardi. In addition, 14 Indian sequences formed a large monophyletic group with published Theileria orientalis sequences. The broad range of sequence identity (95.8-100%) of T. annulata designated sequences indicated the presence of different Theileria spp. A closer analysis revealed the presence of three Theileria spp., namely, T. annulata, T. orientalis, and two isolates (DQ287944 and EU083799) closely related to T. lestoquardi. The true T. annulata sequences manifested 98.8-100% nucleotide identity within them. EU083799 and 14 misidentified Indian T. annulata sequences exhibited the highest similarity with T. lestoquardi (98.6-98.8%) and T. orientalis (98.0-99.9%) in comparison with the other Theileria spp. of domestic and wild ruminants.

CONCLUSION

In the course of analyzing the genetic diversity of T. annulata, we identified the nearly complete 18S rRNA gene sequences of other Theileria spp. that have not only been misidentified as T. annulata in the GenBank™, but are also published as T. annulata. Moreover, a high level of sequence conservation was noticed in the 18S rRNA gene of true T. annulata and T. orientalis sequences.

Genetic diversity in the Tams1 gene of Theileria annulata (Duschunkowsky and Luhs, 1904) infecting cattle

The present study describes the genetic diversity in the Tams1 gene (733 bp) of Theileria annulata along with the sequence, phylogenetic and haplotype analyses of the Indian isolates. The phylogenetic analyses displayed distinct clustering of the Indian isolates into three groups suggesting the presence of three genotypes, hitherto designated as T. annulata genotypes 1-3 (G1-G3). Genotype 3 seems to be novel containing only two newly generated sequences. Indian isolates displayed 88.4-100% and 82.2-100% similarity with each other at nucleotide (nt) and amino acid (aa) levels, respectively. However, the newly generated sequences (n = 36) showed 90.5-100% and 84.3-100% identity between them at nt and aa levels, respectively. The most diverse and heterogeneous genotype, G1, exhibited the highest number of polymorphic sites (S = 148), haplotypes (h = 16) and nucleotide differences (k = 43.23) besides haplotype (Hd = 0.903 ± 0.031) and nucleotide (π = 0.059 ± 0.005) diversities. Neutrality indices suggested a respective decrease and increase in population sizes of G1 and G2 genotypes in India. The nucleotide sequence analyses indicated the presence of extensive sequence variations between nucleotide positions 1-124, 194-257 and 396-494. The N-terminus of Tams1 protein displayed a considerable sequence variability with extensive variations in two regions, between amino acid positions 1-39 and 127-172, as compared to the conserved carboxyl terminus.