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

Characterization and identification of novel polymorphisms in the OR51H1 gene associated with resistance/tolerance to natural Theileria infection in Vechur (Bos indicus) and crossbred (B. taurus × B. indicus) cattle in a Theileria-endemic region

Theileriosis is a haemo protozoan disease affecting livestock, causing high mortality and production loss. Although therapeutic and vector control measures are available, selecting Theileria-resistant/tolerant animals is considered an alternate control measure for which the identification of genetic markers is essential. The present study screened dwarf Vechur (Bos indicus) and crossbred (B. taurus × B. indicus) cattle for natural infection to Theileria using peripheral blood smear and clinical and haematology evaluation. Polymorphisms in a previously identified candidate gene were screened in olfactory receptor family 51 sub-family H member-1 (OR51H1, single exon, 948 bp) associated with Theileria tolerance. The association of Theileria resistance/tolerance with different genotypes was also evaluated. Theileria was detected in 32.45 % of crossbred and 6 % of Vechur. Resistance and tolerance to Theileria infection varied significantly across breeds, and Vechur was found resistant (no parasites detected in the blood smear) compared to crossbred (Odds ratio-0.14, 95 % CI- 0.03-0.42, p < 0.01). The OR51H1 gene was characterised by sequencing in Vechur cattle, and putative single nucleotide variations (SNVs) were identified by pooled sample sequencing and genotyped using PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism). The five SNVs, T271C, C272T, T362C, G180A and T181G, exhibited polymorphic patterns. SNVs T271C and C272T were tolerant, while T362C, G180A and T181G were non-tolerant. The population was in Hardy-Weinberg equilibrium for SNV G180A in crossbred cattle and for SNVs T271C, C272T, and G180A in Vechur cattle. The SNVs T271C in crossbred and T362C, G180A and T181G in Vechur were significantly associated with Theileria tolerance. The identified SNVs also had moderate to high heterozygosity and polymorphic information content. The findings imply that these SNVs can be used as potential markers to select Theileria-resistant Vechur and tolerant crossbred cattle using marker-assisted selection.

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.

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.

Comprehensive genetic diversity and molecular evolutionary analysis of Theileria annulata isolates based on TAMS 1 gene

Molecular epidemiological studies related to the phylogenetic characterization of Theileria annulata are important in delineating the evolutionary history of the parasite. In the current study, the Theileria annulata (T. annulata) merozoite surface antigen 1 (TAMS 1) gene from 14 bovine isolates of T. annulata originating from semi-arid zone of northern India were amplified and sequenced. TAMS 1 gene sequences (n= 337) reported from 16 countries were subsequently analyzed for haplotype network along with genetic diversity. A total of five haplotypes out of the 14 sequenced isolates and 92 haplotypes out of 337 worldwide sequences are documented in this study. Phylogenetic and molecular evolutionary analyses based on TAMS 1 gene sequences showed that T. annulata is dissipated across different countries and numerous strains are closely linked, even though they belong to different geographical locations. The nucleotide homology between 14 isolates from northern India varied between 91.3 and 100%, whereas it was between 31.5 and 100% when sequences across the globe were compared. Haplotype 14 was recognized as most widely distributed haplotype, with 46 isolates circulating in 10 countries. Globally, negligible genetic distance (FST˂0.15) and very high gene flow (Nm˃1) was found in the five populations of the world (South Asia, East Asia, West Asia, Europe and Africa), supporting the absence of clearly defined subgroups in the phylogenetic analysis. Significant negative values of neutrality tests; Tajima's D (D) and Fu and Li's F (F) provided evidence for recent population expansion through positive selection of advantageous variations.

Molecular investigation of bovine tropical theileriosis outbreak in an organized dairy cattle farm in Madhya Pradesh, India

The present study aimed to investigate an outbreak of Theileria annulata (T. annulata) infection in an organized dairy cattle farm in Madhya Pradesh, India, using clinical and molecular techniques. Following the deaths of two crossbred cattle in March 2021, 43 blood samples were collected from infected and apparently healthy animals and examined by blood smear and polymerase chain reaction (PCR) techniques. The blood smear examination showed that 23.25% of samples were positive for Theileria organisms, while conventional PCR targeting the 18S ribosomal RNA (18S rRNA) and T. annulata merozoite surface antigen-1 (TAMS-1) genes revealed that 32.55% of samples were positive for T. annulata. PCR targeting cytochrome b (Cytb) gene showed 46.51% of samples were positive for T. annulata. Haematological analysis confirmed clinical signs of infection in affected animals, which were treated with buparvaquone @ 2.5 mg/kg body weight intramuscularly along with supportive medicine. Two 18S rRNA gene amplicons were sequenced and analysed in a phylogenetic tree and haplotype network with 54 Indian and 38 foreign sequences. The phylogenetic tree revealed two groups with a high posterior probability and bootstrap value, while the haplotype network revealed 35 haplotypes, with haplotype 1 (H1) being the most abundant and several single haplotypes clustering around it, indicating fast and widespread expansion. Genetic diversity indices and neutrality tests confirmed that the population was expanding. These studies highlight the significance of prompt and precise diagnosis and management of T. annulata outbreaks and provide insights into its evolutionary history and population dynamics of T. annulata in India, which could aid improving disease preventive and control strategies.

Simultaneous molecular detection of Anaplasma marginale and Theileria annulata in cattle blood samples collected from Pakistan-Afghanistan boarder region

Theileria annulata (T. annulata) and Anaplasma marginale (A. marginale) are among the most extensively reported tick borne pathogens and are associated with huge economic losses worldwide. A total of 298 cattle blood samples were screened to report the presence of these two pathogens. The samples were collected from apparently healthy cattle (Achai, n = 155, Jersy, n = 88 and crossbred, n = 55) in Bajaur district of Khyber Pakhtunkhwa (KPK) during June and July of 2022. A total of 31 out of 298 cattle (10.4%) were found infected with T. annulata as PCR amplified a 156 base pair fragment from Tams-1 gene of T. annulata from their blood. While 16/298 animals (5.4%) were found infected with A. marginale as they amplified a 382 base pair fragment specific for msp5 gene of this bacterium. Three animals (1%) were found co infected. Cattle susceptibility to T. annulata infection was significantly higher than A. marginale infection (P < 0.001). Phylogenetic analysis revealed that Pakistani isolates of both detected pathogen clustered together and were closely related isolates from worldwide countries. Prevalence of T. annulata varied significantly among the sampling sites (P = 0.05) while no such association was observed for A. marginale among the tested cattle. Epidemiological data analysis revealed that none of the studied risk factors was found associated either with the prevalence of T. annulata or A. marginale (P > 0.05) among enrolled cattle. In conclusion, our study has revealed a relatively higher prevalence of T. annulata than A. marginale in cattle from the Bajaur district in KPK. This information is important for improving the productivity of the livestock sector, which is one of the main sources of income in the country. It is recommended that this data be taken into account for the development and implementation of effective tick control programs, as well as for the improvement of livestock management practices to prevent and manage TBDs in Pakistan.

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.

Metabolomic profiling of bovine leucocytes transformed by Theileria annulata under BW720c treatment

Background

When Theileria annulata infects host cells, it undertakes unlimited proliferation as tumor cells. Although the transformed cells will recover their limited reproductive characteristics and enter the apoptosis process after treatment with buparvaquone (BW720c), the metabolites and metabolic pathways involved are not clear.

Methods

The transformed cells of T. annulata were used as experimental materials, and the buparvaquone treatment group and DMSO control group were used. Qualitative and quantitative analysis was undertaken of 36 cell samples based on the LC–QTOF platform in positive and negative ion modes. The metabolites of the cell samples after 72 h of drug treatment were analyzed, as were the different metabolites and metabolic pathways involved in the BW720c treatment. Finally, the differential metabolites and metabolic pathways in the transformed cells were found.

Results

A total of 1425 metabolites were detected in the negative ion mode and 1298 metabolites were detected in the positive ion mode. After drug treatment for 24 h, 48 h, and 72 h, there were 56, 162, and 243 differential metabolites in negative ion mode, and 35, 121, and 177 differential metabolites in positive ion mode, respectively. These differential metabolites are mainly concentrated on various essential amino acids.

Conclusion

BW720c treatment induces metabolic disturbances in T. annulata-infected cells by regulating the metabolism of leucine, arginine, and l-carnitine, and induces host cell apoptosis.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05450-0.

Molecular Epidemiology and Species Diversity of Tick-Borne Pathogens of Animals in Egypt: A Systematic Review and Meta-Analysis

Ticks and tick-borne pathogens (TTBPs) are listed among the most serious concerns harming Egyptian livestock’s productivity. Several reports on tick-borne pathogens (TBPs) from various geographical regions in the country were published. However, data on the molecular characterization of TBPs are the most beneficial for understanding the epidemiology of this important group of pathogens. In this study, we present the first meta-analysis on the molecular epidemiology and species diversity of TBPs infecting animals in Egypt. All published studies on TBPs were systematically collected from various databases (PubMed, Scopus, ScienceDirect, the Egyptian Knowledge Bank, and Google Scholar). Data from eligible papers were extracted and subjected to various analyses. Seventy-eight studies were found to be eligible for inclusion. Furthermore, ticks infesting animals that were molecularly screened for their associated pathogens were also included in this study to display high species diversity and underline the high infection risk to animals. Theileria annulata was used as parasite model of TBPs to study the genetic diversity and transmission dynamics across different governorates of Egypt. This study extends cross-comparisons between all published molecular data on TBPs in Egypt and provides resources from Egyptian data in order to better understand parasite epidemiology, species diversity, and disease outcome as well as the development and implementation of prevention and control methods for public health, veterinary care practitioners, and animal owners all over the country.

Molecular insights into the population structure and haplotype network of Theileria annulata based on the small-subunit ribosomal RNA (18S rRNA) gene

The present study was conducted to elucidate the population genetic diversity and haplotype network of Theileria annulata based on all available nearly complete 18S rRNA gene sequences in the GenBank™. In total, 52 sequences of the nuclear 18S rRNA gene used to assess the relationship of T. annulata with their country of origin identified 34 haplotypes. Haplotype 4 was widespread, occurring in India, China, Turkey and Iran, while the remaining haplotypes were singleton and unique to one country. Haplotype 4 displayed numerous single haplotypes around it and the stellate shape of the network suggested a rapid population expansion. India exhibited the largest number of haplotypes (h = 25) followed by Turkey (h = 6), China (h = 4), and Iran and Italy (h = 1). No geographical clustering of haplotypes was recorded. Nucleotide diversity was the highest in the Turkish followed by the Indian and Chinese populations. Similarly, haplotype diversity was the highest in China followed by Turkey, and the lowest in India. The overall dataset exhibited a low nucleotide diversity (0.00253 ± 0.00035), but high haplotype diversity (0.917 ± 0.034). It suggested the presence of only minor differences (01-11 nucleotide) between haplotypes which was also evident from the haplotype network. A high level of genetic diversity was documented within the Indian, Chinese and Turkish populations of T. annulata, whereas little genetic differentiation was noticed among these populations with a very high level of gene flow. Analysis of molecular variance (AMOVA) of T. annulata sequences revealed higher genetic variation within countries (83.58%) as compared to the variation among countries (16.42%). Neutrality indices, viz., Tajima's D, Fu and Li's F, Fu's F_s, and R_2, along with the unimodal mismatch distributions demonstrated a recent population expansion of T. annulata in India and the overall dataset. However, the non-significant values of Tajima's D, Fu and Li's F, and Fu's F_s for the Chinese population along with a bimodal mismatch distribution signified a constant population size. For the Turkish population, the neutrality and mismatch distribution tests either indicated a constant or a slight increase in population size. The present study provides novel insights into the population genetics and haplotype network of T. annulata based on the 18S rRNA gene for the first time.