Molecular Characterization and Assessment of Risk Factors Associated with Theileria annulata Infection

Monitoring of Coxiella burnetii in the Iberian lynx (Lynx pardinus)

Coxiella burnetii is a multi-host bacterium of major public and animal health concern. This pathogen circulates among several wild species in the Iberian Peninsula, however, the role of the Iberian lynx (Lynx pardinus) in the epidemiology of this emerging pathogen is still unknown. The objective of this work was to assess the circulation of C. burnetii in Iberian lynx populations from the Iberian Peninsula and to study the molecular characterisation of this pathogen in lynxes and their feeding ticks. A total of 922 lynxes, including free-ranging and captive individuals, were sampled between 2010 and 2022 for the collection of sera (n = 543), spleen samples (n = 390) and ticks (n = 357 from 61 lynxes). The overall seroprevalence was 7.7 % (42/543; 95 %CI: 5.5-10.0 %), with age being significantly associated with the C. burnetii exposure in free-ranging lynxes. A longitudinal study was also carried out to assess the dynamics of the circulation of C. burnetii in this wild host, revealing that 7 of the 37 longitudinally surveyed individuals seroconverted during the study period. The PCR prevalence was 4.4 % (17/390, 95 %CI: 2.3-6.4 %) for spleen samples and 1.1 % (4/357; 95 % CI: 0.0-2.2) in ticks. This is the first study to evaluate the circulation of C. burnetii in the Iberian lynx and to confirm the infection in this felid. The results obtained show a moderate, wide, homogeneous, and endemic circulation of this bacterium in the Iberian lynx populations.

Babesia bigemina and Theileria annulata infections in cattle: molecular detection, phylogenetic analysis, and assessment of risk factors

Babesia bigemina and Theileria annulata are tick-borne protozoans that cause piroplasmosis in cattle, resulting in huge damages to the livestock industry. The prevalence of these infections depends on various intrinsic and extrinsic risk factors. In Pakistan, there is no information regarding the molecular characterization of Babesia bigemina and the risk factors associated with piroplasmosis. This study aimed to molecularly characterize Babesia spp. and Theileria spp. infecting various cattle breeds in Khyber Pakhtunkhwa, Pakistan, and to shed light on risk factors associated with these infections. Altogether, 219 blood samples were collected from various symptomatic cattle breeds, including Holstein Friesian (65.3%; 143/219), Jersey (21.5%; 47/219) and Sahiwal (13.2%; 29/219). Isolated genomic DNA from these blood samples was used in PCR for the amplification of the 18S rRNA fragment of apicomplexan parasites. Obtained 18S rDNA sequences from cattle hosts showed 99.5% identity with B. bigemina, or 100% with T. annulata. Having an overall infection rate of 61.6% (135/219), the highest infection rate was recorded for T. annulata (43.8%; 95/219), followed by B. bigemina (18.3%; 40/219). Phylogenetic analysis of 18S rDNA sequences revealed that B. bigemina clustered with corresponding species reported from Bolivia, and South Africa, while T. annulata grouped with same species from Italy, India, and Turkey. Among the different risk factors, the breed, season, and tick infestation were found to have a significant (P < 0.05) association with the piroplasmid infections. The information obtained in this study can be employed for effective surveillance and control of babesiosis and theileriosis in Pakistan. In addition to confirming our previous molecular detection of T. annulata in cattle, this study provides the first molecular surveillance and phylogenetic position of B. bigemina and associated risk factors in the study region.

Microscopic and molecular detection of piroplasms among sheep in Upper Egypt

Introduction

Blood parasites pose a significant threat to livestock production in southern Egypt, yet there is a scarcity of information regarding their circulation and epidemiology in sheep in this region. This study aimed to investigate the seroprevalence of blood parasite infections in sheep in Assiut governorate, Upper Egypt.

Methods

A total of 400 blood samples were collected from sheep of varying ages and genders. The preliminary screening for the presence of piroplasms, mainly Babesia and Theileria spp., via microscopic examination, followed by investigation of the potential risk factors linked with the exposure to infection. Moreover, molecular identification of both parasites on some of positive samples was performed using PCR targeting Babesia 18S rRNA and Theileria annulata Tams1 gene.

Results

The microscopic examination revealed that among the examined sheep, there was an overall prevalence of blood parasites at 44% (176 out of 400), with Babesia spp. observed in 14% (56 out of 400) and Theileria spp. in 30% (120 out of 400). Furthermore, the infection rate was non-significantly higher in young animals (50%) compared to adults (38.5%) (P = 0.246). Male sheep exhibited a significantly higher vulnerability to both parasites' infection (63.3%) compared to females (35.7%) (P = 0.011). Interestingly, the prevalence of both blood parasites was significantly higher during the cold season (66.1%) compared to the hot season (15.9%) (P = < 0.001). The molecular analysis identified the presence of Babesia ovis and Theileria annulata among a subsample of the positive sheep's bloods films. The identified species were recorded in the GenBank™ databases and assigned specific accession numbers (OQ360720 and OQ360719 for B. ovis), and (OP991838 for T. annulata).

Conclusions

Taken together, this study confirms a high prevalence of piroplasmosis and offers epidemiological and molecular insights into blood parasites in sheep from Upper Egypt, highlighting the importance of detecting these parasites in various hosts and their competent vectors (ticks).

Seroprevalence and risk factors associated with Theileria parva infection among calves in Narok County, Kenya

Background and Aim

East Coast fever (ECF), caused by Theileria parva, is a devastating disease that causes significant economic losses to cattle production in sub-Saharan Africa. Prevention and control of ECF are challenging in pastoral settings due to inadequate epidemiological information. This study aimed to estimate the seroprevalence and risk factors associated with T. parva infection among calves in different production systems to help design appropriate control interventions.

Materials and Methods

Blood samples were collected from 318 calves and tested using an indirect enzyme-linked immunosorbent assay targeting antibodies against polymorphic immunodominant molecules found on the surface of T. parva. Information on calf characteristics and management practices was also collected during sampling. Descriptive statistics and logistic regression were used to analyze potential risk factors, such as age and acaricide application, where p < 0.05 was considered significant.

Results

Of the 318 calves sampled, 41 (12.89%) were positive for T. parva, with a higher proportion in pastoral systems (36.58%) than in mixed farming systems (34.10%) and agropastoral systems (29.27%). From univariate analysis, calf age (p = 0.002), body weight (p = 0.001), suckling status (p = 0.026), rectal temperature (p = 0.06), calves on pasture (p = 0.022), other feeds (p = 0.004), feed grown within the farm (p = 0.004), acaricide application (p = 0.001), and acaricide application frequency (p = 0.001) were significantly associated with seropositivity. However, calf age (odds ratio [OR], 0.96; 95% confidence interval [CI], 0.91-0.99; p = 0.04), other feeds (OR, 8.82; 95% CI, 1.74-44.63; p = 0.009), and suckling status (OR, 0.38; 95% CI, 0.15-0.99; p = 0.05) were significantly associated with T. parva infection in the multivariable mixed logistic model.

Conclusion

T. parva is circulating in young calves in the study area (and possibly in cattle populations due to maternal transfer of antibodies to the calves). There is a need for molecular surveillance to determine the presence and burden of T. parva infection.

First Molecular-Based Confirmation of Dermacentor marginatus and Associated Rickettsia raoultii and Anaplasma marginale in the Hindu Kush Mountain Range

Ticks of the genus Dermacentor Koch, 1844 (Acari: Ixodidae) are poorly known systematically due to their habitation in harsh topographic environments and high mountains. Dermacentor ticks are diversely distributed in the Palearctic, Nearctic, and Oriental regions. There is no available information on the occurrence of Dermacentor marginatus in Pakistan; thus, the current investigation aimed the first morphological and molecular confirmation of this species and associated Anaplasma marginale and Rickettsia raoultii. Ticks were collected from goats (Capra hircus) and morphologically identified. Genomic DNA was extracted from 18/26 (69.23%) tick specimens, including 11 males and 7 females (1 unfed and 6 fed females). Extracted DNA was subjected to PCR for the amplification of genetic markers like 16S rDNA and cox1 for ticks, 16S rDNA for Anaplasma spp., and gltA and ompB for Rickettsia spp. A total of 26 D. marginatus ticks composed of 19 males (73.07%) and 7 females (26.9%) [1 (3.84%) unfed and 6 (23.07%) fed females] were collected from goats. According to amplicons via BLAST analysis, the 16S rDNA sequence showed 97.28-98.85% identity and the cox1 sequence showed 95.82-98.03% identity with D. marginatus. Additionally, the 16S rDNA sequence for Anaplasma sp. was detected in D. marginatus that showed 100% identity with Anaplasma marginale. Rickettsial gltA and ompB sequences for Rickettsia sp. showed 100% identity with Rickettsia raoultii. In phylogenetic analysis, ticks' 16S rDNA and cox1 sequences clustered with the same species. In phylogenetic analysis, A. marginale based on 16 rDNA clustered with A. marginale, while gltA and ompB sequences clustered with R. raoultii. This is the first study on the genetic characterization of D. marginatus and associated A. marginale and R. raoultii in Pakistan. The northern areas of Pakistan, which need to be explored in terms of ticks and associated pathogens due to their zoonotic threats, have been neglected due to the inaccessible climatic conditions.

Argas persicus and Carios vespertilionis Ticks Infesting Ducks, Domestic Fowls and Bats in Pakistan: First Report on Molecular Survey and Phylogenetic Position of Borrelia anserina

Argasid ticks have the vectorial potential for transmitting disease-causing pathogens to avian hosts, resulting in economic losses that may not be fully estimated. Borrelia species are the responsible agents of borreliosis in poultry, animals and humans. Our previous studies have reported a high prevalence of Argas persicus infesting domestic fowls in Khyber Pakhtunkhwa (KP), Pakistan. However, molecular screening and genetic characterization of Borrelia spp. in A. persicus have been neglected in Pakistan. In this study, we focused on the molecular epidemiology and genetic characterization of Borrelia spp. associated with A. persicus ticks infesting domestic fowls and ducks, and Carios vespertilionis infesting bats in selected districts of KP. Overall, 1818 ticks, including females (415; 23%), males (345; 19%), nymphs (475; 26%) and larvae (583; 32%), were collected from 27 locations in nine districts (Peshawar, Mardan, Swabi, Charsadda, Chitral, Lakki Marwat, Bannu, Bajaur and Hangu) from domestic fowls, ducks and their shelters, and bats. A subset of 197 ticks was selected for DNA extraction and PCR to amplify fragments of the cytochrome c oxidase (cox) gene for ticks and flagellin B (flaB) for the detection and genetic characterization of associated Borrelia spp. Among these, only Borrelia anserina DNA was detected in 40 ticks (27.2%) of different life stages, where highest prevalence was found in female ticks (18; 45%), followed by nymphs (12; 30%), larvae (7; 17.5%) and males (3; 7.5%). Tick infestation in shelters (1081; 77%) was higher than on hosts (323; 23%). The resultant cox amplicons of A. persicus showed 100% identity with the same species reported from Pakistan, China, Iran, Kenya, Kazakhstan, Algeria and Egypt and C. vespertilionis show 100% identity with the species reported from Pakistan, China, Japan, Kenya, Vietnam, Spain, Netherlands, the United Kingdom and Hungry, and clustered with the aforementioned species in the phylogenetic tree. The obtained Borrelia sequences showed 100% identity with B. anserina and revealed a close resemblance to the relapsing fever group and clustered in a monophyletic clade with B. anserina from India, Iran and Brazil in a phylogenetic tree. These results establish the first molecular characterization of B. anserina in A. persicus infesting domestic fowls and ducks in the region, as well as their shelters. To effectively control zoonotic consequences, country-wide surveillance research should be encouraged to screen soft ticks infesting various birds for associated pathogens.

First Detection of Theileria sinensis-like and Anaplasma capra in Ixodes kashmiricus: With Notes on cox1-Based Phylogenetic Position and New Locality Records

Ixodes ticks transmit Theileria and Anaplasma species to a wide range of animals. The spreading of ticks and tick-borne pathogens has been attributed to transhumant herds, and research on these uninvestigated issues has been neglected in many countries, including Pakistan. Recently, we used internal transcribed spacer (ITS) and 16S ribosomal DNA partial sequences to genetically characterize Ixodes kashmiricus ticks and their associated Rickettsia spp. However, the data on its cox1 sequence and associated Theileria spp. and Anaplasma spp. are missing. This study aimed to genetically characterize I. kashmiricus based on the cox1 sequence and their associated Theileria spp. and Anaplasma spp. The I. kashmiricus ticks were collected from small ruminants: sheep (Ovis aries) and goats (Capra hircus) of transhumant herds in district Shangla, Dir Upper and Chitral, Khyber Pakhtunkhwa (KP), Pakistan. Out of 129 examined hosts, 94 (72.87%) (56 sheep and 38 goats) were infested by 352 ticks, including adult females (175; 49.7%) followed by nymphs (115; 32.7%) and males (62; 17.6%). For molecular analyses, 121 ticks were subjected to DNA isolation and PCR for the amplification of the cox1 sequence for I. kashmiricus, 18S rDNA for Theileria spp. and 16S rDNA sequences for Anaplasma spp. The obtained cox1 sequence showed 89.29%, 88.78%, and 88.71% identity with Ixodes scapularis, Ixodes gibbosus, and Ixodes apronophorus, respectively. Phylogenetically, the present cox1 sequence clustered with the Ixodes ricinus complex. Additionally, the 18S rDNA sequence showed 98.11% maximum identity with Theileria cf. sinensis and 97.99% identity with Theileria sinensis. Phylogenetically, Theileria spp. clustered with the T. cf. sinensis and T. sinensis. In the case of Anaplasma spp., the 16S rDNA sequence showed 100% identity with Anaplasma capra and phylogenetically clustered with the A. capra. PCR-based DNA detection targeting the amplification of groEL and flaB sequences of Coxiella spp. and Borrelia spp., respectively, was unsuccessful. This is the first phylogenetic report based on cox1 and new locality records of I. kashmiricus, and the associated T. sinensis-like and A. capra. Significant tick surveillance studies are needed in order to determine the epidemiology of Ixodes ticks and their associated pathogens.

Molecular prevalence and associated risk factors of Theileria annulata infections in dromedary camels in Egypt

Theileria annulata (T. annulata) is an intracellular protozoan, transmitted by ixodid ticks of the genus Hyalomma and affects camels. There are few epidemiological data on T. annulata infection and its associated risk factors in Egyptian camels. Thus, the present study aimed to determine the prevalence of T. annulata in camels using PCR and assess the associated risk factors for infection. A total of 380 blood samples were collected from camels raising in three Egyptian governorates and examined by PCR assay targeting 30-kDa gene to detect the presence of T. annulata infection, beside statistical analysis of associated factors. The results revealed presence of T. annulata with overall prevalence of 21.1%. In addition, the univariate analysis revealed significant (P<0.05) association between prevalence of T. annulata in camels and locality, age, sex, tick infestation, and application of acaricides. Whereas the prevalence of T. annulata was higher in camels of age group >6 to 10 years (38%), females (25.7%) and in infested camels with ticks (29%) and in case of absence of acaricides application (25%). In contrast, the body condition of camels had not significant effect on prevalence of theileriosis in camels. The current study concluded that T. annulata is prevalent in Egyptian camels and that a tick control program is required to reduce the risk of infection.

Genomic Characterization of Theileria luwenshuni Strain Cheeloo

Theileria, a tick-borne intracellular protozoan, can cause infections of various livestock and wildlife around the world, posing a threat to veterinary health. Although more and more Theileria species have been identified, genomes have been available only from four Theileria species to date. Here, we assembled a whole genome of Theileria luwenshuni, an emerging Theileria, through next-generation sequencing of purified erythrocytes from the blood of a naturally infected goat. We designated it T. luwenshuni str. Cheeloo because its genome was assembled by the researchers at Cheeloo College of Medicine, Shandong University, China. The genome of T. lunwenshuni str. Cheeloo was the smallest in comparison with the other four Theileria species. T. luwenshuni str. Cheeloo possessed the fewest gene gains and gene family expansion. The protein count of each category was always comparable between T. luwenshuni str. Cheeloo and T. orientalis str. Shintoku in the Eukaryote Orthologs annotation, though there were remarkable differences in genome size. T. luwenshuni str. Cheeloo had lower counts than the other four Theileria species in most categories at level 3 of Gene Ontology annotation. Kyoto Encyclopedia of Genes and Genomes annotation revealed a loss of the c-Myb in T. luwenshuni str. Cheeloo. The infection rate of T. luwenshuni str. Cheeloo was up to 81.5% in a total of 54 goats from three flocks. The phylogenetic analyses based on both 18S rRNA and cox1 genes indicated that T. luwenshuni had relatively low diversity. The first characterization of the T. luwenshuni genome will promote better understanding of the emerging Theileria. IMPORTANCE Theileria has led to substantial economic losses in animal husbandry. Whole-genome sequencing data of the genus Theileria are currently limited, which has prohibited us from further understanding their molecular features. This work depicted whole-genome sequences of T. luwenshuni str. Cheeloo, an emerging Theileria species, and reported a high prevalence of T. luwenshuni str. Cheeloo infection in goats. The first assembly and characterization of T. luwenshuni genome will benefit exploring the infective and pathogenic mechanisms of the emerging Theileria to provide scientific basis for future control strategies of theileriosis.

First Report of Rickettsia conorii in Hyalomma kumari Ticks

As a vector of wide range of pathogenic agents, ticks pose health threats to wild and domestic animals, and humans. Information is unavailable about the prevalence and spatial survey of Hyalomma kumari ticks and associated Rickettsia spp. in Pakistan. Concerning this knowledge gap, the present study aimed to molecularly detect Rickettsia species associated with H. kumari infesting small ruminants in Khyber Pakhtunkhwa (KP), Pakistan. A total of 409 H. kumari ticks were collected from 163/295 infested hosts with an infestation rate of 55.25%. A total of 204 females, 158 males, and 47 nymphs were collected. Goats were heavily infested by 224 ticks having an infestation rate of 58.33% (98/168), whereas sheep were infested by 185 ticks having a lesser infestation rate of 51.18% (65/127). Genomic DNA extracted from ticks was used for the amplification of tick (cox I, 16S rRNA, ITS-2) species and Rickettsia (gltA, ompA, and ompB) partial genes. Eighty-three ticks were subjected to PCR, and 8/83 (9.6%) were found positive for rickettsial agents. The cox I and 16S rRNA sequences of H. kumari showed 98.90-99.74% identity with H. kumari sequences reported from Pakistan, and phylogenetically clustered to the corresponding species reported from Pakistan and India. The obtained rickettsial gltA, ompA, and ompB sequences showed 100% identity with Rickettsia sp. of the Rickettsia conorii reported from Pakistan. In the phylogenetic trees, rickettsial sequences clustered with uncharacterized Rickettsia sp. from Pakistan and R. conorii from Israel, Russia, South Africa, and India. The present molecular based detection of H. kumari-associated R. conorii will facilitate effective surveillance in the region.

Low Genetic Polymorphism in the Immunogenic Sequences of Rhipicephalus microplus Clade C

Rhipicephalus microplus tick highly affects the veterinary sector throughout the world. Different tick control methods have been adopted, and the identification of tick-derived highly immunogenic sequences for the development of an anti-tick vaccine has emerged as a successful alternate. This study aimed to characterize immunogenic sequences from R. microplus ticks prevalent in Pakistan. Ticks collected in the field were morphologically identified and subjected to DNA and RNA extraction. Ticks were molecularly identified based on the partial mitochondrial cytochrome C oxidase subunit (cox) sequence and screened for piroplasms (Theileria/Babesia spp.), Rickettsia spp., and Anaplasma spp. PCR-based pathogens-free R. microplus-derived cDNA was used for the amplification of full-length cysteine protease inhibitor (cystatin 2b), cathepsin L-like cysteine proteinase (cathepsin-L), glutathione S-transferase (GST), ferritin 1, 60S acidic ribosomal protein (P0), aquaporin 2, ATAQ, and R. microplus 05 antigen (Rm05Uy) coding sequences. The cox sequence revealed 100% identity with the nucleotide sequences of Pakistan's formerly reported R. microplus, and full-length immunogenic sequences revealed maximum identities to the most similar sequences reported from India, China, Cuba, USA, Brazil, Egypt, Mexico, Israel, and Uruguay. Low nonsynonymous polymorphisms were observed in ATAQ (1.5%), cathepsin-L (0.6%), and aquaporin 2 (0.4%) sequences compared to the homologous sequences from Mexico, India, and the USA, respectively. Based on the cox sequence, R. microplus was phylogenetically assembled in clade C, which includes R. microplus from Pakistan, Myanmar, Malaysia, Thailand, Bangladesh, and India. In the phylogenetic trees, the cystatin 2b, cathepsin-L, ferritin 1, and aquaporin 2 sequences were clustered with the most similar available sequences of R. microplus, P0 with R. microplus, R. sanguineus and R. haemaphysaloides, and GST, ATAQ, and Rm05Uy with R. microplus and R. annulatus. This is the first report on the molecular characterization of clade C R. microplus-derived immunogenic sequences.