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

Molecular occurrence and genetic identification of Babesia spp. and Theileria spp. in naturally infected cattle from Thailand

Piroplasm including Babesia spp. and Theileria spp. in cattle can cause illness that affects livestock productivity, resulting in significant production losses, especially in tropical and subtropical regions such as Thailand. This study aimed to investigate the prevalence of bovine piroplasms and to identify these blood parasites based on the 18S ribosomal RNA gene in cattle in the northeastern part of Thailand. Piroplasmid infections among beef and dairy cattle were examined using nested PCR. Furthermore, amplicon DNA was sequenced and analyzed, and a phylogenetic tree was constructed to determine the genetic diversity and relationships of the parasite in each area. A total of 141 out of 215 (65.6%) cattle were positive for infection with Babesia or Theileria. DNA analysis revealed that infection by Babesia bigemina, Babesia bovis, Theileria orientalis, Theileria sinensis, and Theileria sp. were common piroplasms in cattle in this region, with a high sequence shared identity and similarity with each other and clustered with isolates from other countries. This study provides information on the molecular epidemiology and genetic identification of Babesia spp. and Theileria spp. in beef and dairy cattle to provide a better understanding of piroplasm infection in cattle in this region, which will help control these blood parasites. Moreover, this is the first report identifying T. sinensis circulating among Thai cattle.

First Insight into the Phylogenetic Diversity of Bovicola caprae Infesting Goats of Different Agro-climatic Locations in India

Bovicola caprae is an important obligate ectoparasite of goats worldwide including India. The present study aimed at the molecular confirmation, phylogenetics and population structure analyses of B. caprae infesting goats of three different agro-climatic locations in India, by targeting the mitochondrial cytochrome C oxidase subunit 1 (cox1) genetic marker. The phylogenetic tree exhibited the presence of two different lineages of B. caprae. The sequences generated herein clustered in lineage 2 along with the GenBank™ archived sequences from China and Iran. The sequences generated herein also showed the circulation of sub-lineages of B. caprae in India based on the analysis of pairwise genetic distances between sequences and median-joining haplotype network. The population structure analyses revealed low nucleotide (0.00353 ± 0.00291 and 0.02694 ± 0.00363) and high haplotype (0.667 ± 0.314 and 0.618 ± 0.104) diversities for the present study isolates as well as for the complete dataset, respectively, which evinced a recent demographic expansion. High genetic differentiation (FST value = 0.97826) and low gene flow (Nm = 0.00556) were also recorded in the different lineages/populations. In conclusion, the present study addressed the research gap and provided the first insight into the phylogenetics of the goat louse B. caprae and highlighted the circulation of sub-lineages of the ectoparasite in India.

Cladistics of Rhipicephalus microplus and laboratory assessment of acaricidal efficacy of different Cassia fistula extracts against R. microplus larvae

The present study aimed to analyze the cladistics and population structure analysis of Rhipicephalus microplus ticks infesting buffaloes in Haryana, India, as well as the assessment of the anti-tick efficacy of the ethanolic extracts of Cassia fistula (bark, pod pulp, and flowers) against R. microplus larvae. The molecular characterization and population structure analysis were performed by targeting the amplification of the partial mitochondrial cytochrome C oxidase subunit 1 (cox1) gene, whereas anti-tick efficacy was evaluated using a larval packet test. The sequences generated herein revealed a homology of 98.26-100% to the GenBank-archived R. microplus sequences. In population structure analysis, high haplotype (0.500 ± 0.265) and low nucleotide (0.002 ± 0.001) diversities were recorded for the sequences generated in this study. Significantly negative neutrality indices were recorded for the overall dataset. The extracts were found to significantly affect mortality rates in a dose-dependent manner, and the ethanolic extracts of the bark, pod pulp, and flowers of C. fistula exhibited median lethal concentration (LC50) values of 27.989, 40.457, and 49.43 mg/mL, respectively. The LC50 value recorded for the combination of the ethanolic extracts of the bark, flower, and pod pulp of C. fistula was 19.724 mg/mL, whereas the synthetic acaricide ivermectin had an LC50 value of 351.56 mg/mL. In conclusion, R. microplus populations infesting cattle and buffalo hosts in India exhibited negligible genetic differentiation and high gene flow between them. Additionally, the combination of all C. fistula extracts could serve as a potential substitute for the synthetic acaricide.

Population structure and haplotype network analyses of Hyalomma anatolicum based on the large subunit ribosomal RNA (16S rRNA) gene

Hyalomma anatolicum, an Anatolian hard tick is a well-recognized vector involved in the transmission of various pathogens to animals and humans. The present study elucidated the population structure and haplotype network of H. anatolicum based on the mitochondrial large subunit ribosomal RNA (16S rRNA) gene sequence. The population structure and haplotype network analysis of 75 sequences archived in the GenBank, including the 15 sequences generated herein, yielded 24 haplotypes. Haplotype 1 (Hap_1) was the predominant haplotype consisting of 45 sequences from India, China, Pakistan, Turkey, Egypt, Iraq, and Tajikistan. The complete haplotype network exhibited a stellate conformation, highlighting a recent population expansion. The overall dataset, together with the sequences corresponding to India, China, and Pakistan, showed a high haplotype (0.638 ± 0.065, 0.671 ± 0.103, 0.753 ± 0.099, and 0.854 ± 0.061, respectively) and low nucleotide (0.00407 ± 0.00090, 0.00525 ± 0.00196, 0.00680 ± 0.00233, and 0.00453 ± 0.00056, respectively) diversity, further emphasized a recent population expansion. The neutrality indices including Tajima's D, Fu and Li's D, and Fu and Li's F for the complete dataset (- 2.661, - 6.008, and - 5.649, respectively) as well as for the sequences from India (- 2.223, - 3.414, and - 3.567, respectively) were negative, suggesting deviation from neutrality and a recent population expansion. The present study provided novel insights into the population structure and haplotype networks of H. anatolicum based on the mitochondrial 16S rRNA gene, and the different tests inferred a low genetic differentiation and suggested a recent population expansion of this economically important tick species.

Unveiling genotypic diversity of Theileria orientalis in lethal outbreaks among bovines in Karnataka, India

Theileria orientalis, the causal agent of oriental theileriosis, is known to cause mild disease in cattle and buffalo across the world. Recently, different genotypes of T. orientalis have emerged as pathogenic, causing high reported morbidity in cattle. This study focuses on investigating three suspected outbreaks of oriental theileriosis that resulted in fatalities among crossbred and indigenous bulls in Karnataka, India. Examination of blood smears revealed the presence of T. orientalis piroplasms within erythrocytes. The genetic characterization of T. orientalis was conducted by targeting specific markers, including the mpsp gene, p23 gene, and ribosomal DNA markers (18S rRNA gene, ITS-1, and ITS-2). Analysis based on the 18S rRNA gene unveiled the presence of both Type A and Type E genotypes of T. orientalis in the outbreaks. The mpsp gene-based analysis identified genotype 7 of T. orientalis in crossbred cows, whereas genotype 1 (Chitose B) was found to be present in indigenous bulls. Haplotype network analysis based on the mpsp gene revealed the presence of 39 distinct haplotypes within the 12 defined genotypes of T. orientalis with a high haplotype diversity of 0.9545 ± 0.017. Hematological and biochemical analysis revealed a decrease in calcium, hemoglobin levels, red blood cell counts, and phosphorus. This study constitutes the initial documentation of a clinical outbreak of oriental theileriosis in indigenous bulls with genotype 1 (Chitose 1B). Substantial epidemiological investigations are imperative to gain a comprehensive understanding of the geographical distribution of distinct genotypes and the diverse clinical manifestations of the disease across various hosts.

Understanding the role of pigs in the transmission of zoonotic Echinococcus ortleppi in Haryana, India

The etiological agents of zoonotic cystic echinococcosis comprise the Echinococcus granulosus sensu lato (s.l.) species complex. The present study was aimed at investigating the zoonotic genotypes of Echinococcus granulosus s.l. circulating in the pig population of Haryana, India. Out of 253 slaughtered pigs screened, 5 showed the presence of hydatid cysts. The amplification of the partial mitochondrial NADH dehydrogenase subunit 1 (nad1) gene for the molecular confirmation and phylogenetics of the retrieved metacestodes (n = 2) revealed the presence of E. ortleppi. The sequences generated herein exhibited 99.80% homology to the GenBank archived E. ortleppi sequences. Cladistics targeting genetic diversity and haplotype network analysis involved 37 E. granulosus s.l. GenBank archived sequences from India corresponding to different hosts (large and small ruminants and humans) along with the sequences (n = 2) generated in the present study. Overall, 14 haplotypes with high haplotype (0.780 ± 0.059) and low nucleotide (0.033 ± 0.010) diversities were recorded for the overall data set, which evinced a population expansion. The median-joining haplotype network revealed a stellate shape of E. granulosus sensu stricto (s.s.) sequences, which was indicative of rapid population expansion. High genetic differentiation (FST = 0.840 - 0.983) and low gene flow (Nm = 0.003 - 0.047) were recorded between the pig intermediate hosts infected with E. ortleppi and other hosts infected with E. granulosus s.s. The findings are of paramount significance for the formulation of effective control strategies considering the public health and economic impact of cystic echinococcosis.

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.

Cladistics of Echinococcus granulosus Sensu Stricto Genotypes Infecting the Slaughtered Pigs

PURPOSE

Echinococcus granulosus sensu stricto (s. s.) consists of the most widespread genotypes (G1, G3) implicated in human cystic echinococcosis worldwide. The present study aimed to evaluate the role of pigs in the transmission dynamics of E. granulosus s. s. genotypes, including the phylogenetics, evolutionary divergence and haplotype network analyses of north Indian pig isolates along with GenBank™ archived sequences.

METHODS

In totality, 223 slaughtered pigs were thoroughly screened for the presence of hydatid cysts. The amplification of the partial mitochondrial cytochrome C oxidase subunit 1 gene was performed for the molecular confirmation and phylogenetics of the retrieved metacestodes.

RESULTS

The molecularly confirmed and sequenced present study isolates (n = 2) clustered with the E. granulosus genotype G3. The very low evolutionary divergence values (< 0.0027) of the present study isolates with E. ganulosus genotype G3 indicated the absence of diverged lineages. The significantly negative values of Tajima's D (- 2.598) and Fu and Li's D (- 4.054) of the overall dataset and for the Asian sequences signified an expansion of population size. The overall dataset exhibited low nucleotide (0.067 ± 0.055) and high haplotype (0.574 ± 0.015) diversities, also alluding to demographic expansion. The haplotype network showed that the pig isolates from South America and Europe constituted the predominant haplotype, Hap_2 along with Hap_3 and Hap_6, primarily associated to E. granulosus genotype G1; whereas, the Indian isolates formed different haplotypes (Hap_1 and Hap_5) belonging to genotype G3.

CONCLUSIONS

The present study highlighted the important role of pigs in the transmission of E. granulosus s. s., which is of paramount significance given the public health and economic impact of cystic echinococcosis.

Eimeria zuernii (Eimeriidae: Coccidia): mitochondrial genome and genetic diversity in the Chinese yak

BACKGROUND

Coccidiosis caused by Eimeria zuernii (Eimeriidae: Coccidia) represents a significant economic threat to the bovine industry. Understanding the evolutionary and genetic biology of E. zuernii can assist in new interaction developments for the prevention and control of this protozoosis.

METHODS

We defined the evolutionary and genetic characteristics of E. zuernii by sequencing the complete mitogenome and analyzing the genetic diversity and population structure of 51 isolates collected from eight yak breeding parks in China.

RESULTS

The 6176-bp mitogenome of E. zuernii was linear and encoded typical mitochondrial contents of apicomplexan parasites, including three protein-coding genes [PCGs; cytochrome c oxidase subunits I and III (cox1 and cox3), and cytochrome b (cytb)], seven fragmented small subunit (SSU) and 12 fragmented large subunit (LSU) rRNAs. Genome-wide comparative and evolutionary analyses showed cytb and cox3 to be the most and least conserved Eimeria PCGs, respectively, and placed E. zuernii more closely related to Eimeria mephitidis than other Eimeria species. Furthermore, cox1-based genetic structure defined 24 haplotypes of E. zuernii with high haplotype diversities and low nucleotide diversities across eight geographic populations, supporting a low genetic structure and rapid evolutionary rate as well as a previous expansion event among E. zuernii populations.

CONCLUSIONS

To our knowledge, this is the first study presenting the phylogeny, genetic diversity, and population structure of the yak E. zuernii, and such information, together with its mitogenomic data, should contribute to a better understanding of the genetic and evolutionary biological studies of apicomplexan parasites in bovines.

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.

Phylogeography and demographic dynamics of Rhipicephalus microplus from North India

Rhipicephalus microplus, a hematophagous vector prevalent in the tropics and subtropics, is responsible for huge economic losses throughout the globe. However, the taxonomy of the tick species, especially prevalent in north India and south China has been challenged in the recent past. The present study attempted to assess the cryptic status of R. microplus ticks of north India based on two mitochondrial markers; the 16S rRNA and cox1 gene sequences. The phylogenetic tree corresponding to both markers demonstrated the presence of three distinct genetic assemblages/ clades of R. microplus. The present study isolates (n = five and seven for the cox1 and 16S rRNA gene sequences, respectively) from north India along with other isolates from India assorted in the R. microplus clade C sensu. Based on the median joining network analysis corresponding to the 16S rRNA gene sequence, 18 haplotypes were recorded, exhibiting a stellate shape, which was indicative of rapid population expansion. For the cox1 gene, the haplotypes corresponding to clades A, B and C were distantly placed with two exceptions. While performing the population structure analysis, low nucleotide (0.04745 ± 0.00416 and 0.01021 ± 0.00146) and high haplotype diversities (0.913 ± 0.032 and 0.794 ± 0.058) were recorded for the different clades of R. microplus based on the cox1 and 16S rRNA mitochondrial markers, respectively. Eventually, high genetic differentiation and low gene flow were recorded among the different clades. A negative value for the neutrality indices (Tajima's D = -1.44125, Fu's Fs = -4.879, Fu and Li's D = -2.78031 and Fu and Li's F = -2.75229) corresponding to the 16S rRNA gene for the overall dataset evinced an expansion of population size. Based on the detailed studies, it was inferred that the R. microplus tick species circulating in north India belonged to clade C sensu, similar to that of the species prevalent in the other parts of the country as well as in the Indian subcontinent.

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 and population diversity of Toxocara cati (Schrank, 1788) Brumpt, 1927, on the basis of the internal transcribed spacer (ITS) region

The present investigation was aimed to study the sequence, phylogenetic and haplotype analyses of Toxocara cati based on the ITS region, along with the genetic diversity, demographic history and population-genetic structure. The maximum likelihood tree based on Kimura 2-parameter model was constructed using the complete ITS region of all the nucleotide sequences (n = 57) of Toxocara spp. and other related ascarid worms available in the GenBank™. It placed all the sequences of T. cati into four major clades designated as T. cati genotypes 1-4 (TcG1-G4). A total of 66 signature nucleotides were identified in the ITS region between genotypes. The median-joining haplotype network displayed a total of 24 haplotypes, with China exhibiting the highest number of haplotypes (h = 20) followed by India (h = 4), and Japan and Russia (h = 1). It indicated a clear distinction between all the four genotypes. The pairwise F_ST values between all the genotypes indicated huge genetic differentiation (> 0.25) between different T. cati genotypes. Moreover, the gene flow (Nm) between T. cati genotypes was very low. Results of AMOVA revealed higher genetic variation between genotypes (92.82%) as compared to the variation within genotypes (7.18%). The neutrality indices and mismatch distributions for the G1-G4 genotypes, Indian isolates and the overall dataset of T. cati indicated either a constant population size or a slight population increase. The geographical distribution of all the genotypes of T. cati is also reported. This is the first report of genotyping of T. cati on the basis of the ITS region.

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.