New records of Amblyomma gervaisi from Pakistan, with detection of a reptile-associated Borrelia sp

Leading report on the molecular prevalence of emerging pathogens Hepatozoon sp. and Lankesterella sp. in the blood samples of seven wild lizard species

Despite of having rich reptilian fauna, the wild lizards from Pakistan remained unexplored for the presence of blood borne parasites. The present study was designed to report the molecular epidemiology and phylogenetic evaluation of various blood borne pathogens (Hepatozoon sp., Schellackia spp., Toxoplasma gondii, Plasmodium spp., Haemoproteus spp. and Leucocytozoon spp.) in blood samples of wild lizards (N = 84), trapped during March 2022 till June 2023 from district Karak and Dir in Khyber Pakhtunkhwa (Pakistan). Seven lizard species were identified during present study. Our results revealed that 11 out of 84 (13%) lizards were Haemogregarines infected. Infected lizards included Laudakia (L.) tuberculata (1/4, 25% prevalence), L. pakistanica (3/15, 20%) and L. agrorensis (7/53, 13%). DNA sequencing and BLAST analysis confirmed the presence of Hepatozoon sp. and Lankesterella sp. While the lizards were negative for the remaining screened pathogens. Phylogenetic analysis of both pathogens revealed genetic diversity among the Pakistani isolates and they clustered with isolates detected in reptiles, birds and rodents from different countries. For L. tuberculata, Haemogregarines prevalence significantly varied between the sample collection sites. In conclusion, this is the first report from Pakistan documenting a relatively higher Haemogregarines infection rate in wild Pakistani lizards. Further, comprehensive and large-scale studies must be conducted in unexplored geo-climatic regions of Pakistan to report the actual prevalence of Haemogregarines among the wild lizards as well as in other wildlife species. These findings will add to our knowledge regarding the genetic diversity and the interactions of these parasites with their hosts that will lead towards parasite control.

Molecular detection and phylogeny of Hepatozoon ophisauri and Toxoplasma gondii in wild lizards from Khyber Pakhtunkhwa, Pakistan

Despite Pakistan's rich reptilian diversity, wild lizards have largely been unexplored for vector-borne parasites. This study reports the prevalence and phylogenetic assessment of Hepatozoon ophisauri and Toxoplasma gondii DNA in the blood samples of 101 wild lizards from Khyber Pakhtunkhwa, Pakistan, captured at altitudes of 1200 to 2250 m above sea level between March 2022 and June 2023. PCR-based molecular analysis identified H. ophisauri DNA in three (3%) lizards of all Laudakia tuberculata. Additionally, 45 (44.5%) lizards tested positive for T. gondii DNA, including L. agrorensis, L. pakistanica, L. tuberculata, and Abblepharus pannonicus. Sequence analysis and BLAST confirmed the presence of H. ophisauri and T. gondii. Phylogenetic analysis showed genetic diversity among the H. ophisauri and T. gondii nucleotide sequences, clustering with reference sequences from reptiles, birds, and ticks. Toxoplasma gondii prevalence varied among species, with the highest rates in L. agrorensis and A. pannonicus (67%), followed by L. pakistanica (45%) and L. tuberculata (43%). This study is the first from Pakistan to report high T. gondii prevalence and low H. ophisauri prevalence in wild lizards. Similar and larger-scale studies in unexplored regions of Pakistan are needed to enhance understanding of these pathogens' genetic diversity and host-parasite interactions.

Unveiling hidden threats: Molecular surveillance of bacterial and protozoan infections in Asian water monitors (Varanus salvator) at Thailands Khao-zon Wildlife Breeding Station

Asian water monitors (Varanus salvator) are often exposed to harmful microorganisms such as species of Hepatozoon, Anaplasma, Rickettsia, and Borrelia. These pathogens pose significant risks to wildlife and public health. In this study, we aimed to investigate the presence and genetic diversity of pathogenic microorganisms in V. salvator from the Khao-zon Wildlife Breeding Station, Ratchaburi, Thailand. Sixteen V. salvator were collected, and thin blood smears and polymerase chain reaction (PCR) assays of the DNA isolated from blood were used to identify infections. Blood smear analysis revealed two positive samples (2/16; 12.5%) for Hepatozoon sp. PCR results confirmed the presence of Hepatozoon sp. (6/16; 37.5%), Anaplasma sp. (3/16; 18.75%), Rickettsia sp. (2/16; 12.5%), Borrelia sp. (4/16; 25.0%) and Burkholderia sp. (1/16; 6.25%). Phylogenetic analysis revealed 99.86% similarity of Hepatozoon sp. sequences with Hepatozoon ophisauri (GenBank: MN723845) in Pseudopus apodus from Iran, whereas the gltA gene of Rickettsia sp. was closely related to Rickettsia conorii raoultii (GenBank: MF002515) in Dermacentor nuttalli from China. This study represents the first detection of the five pathogens in V. salvator from Thailand and provides valuable insights into the genetic diversity of these microorganisms. Our findings suggest that V. salvator may serve as reservoir for multiple pathogens, posing potential threats to both wildlife and humans. The presence of zoonotic agents such as Rickettsia spp., Anaplasma spp., and Borrelia spp. underscores the importance of continued surveillance in wildlife populations to mitigate the risk of emerging infectious diseases.

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.

Description of a new Ornithodoros (Pavlovskyella) (Ixodida: Argasidae) tick species from Pakistan

The genus Ornithodoros is notably diverse within the family Argasidae, comprising approximately 134 species distributed among 4 subgenera, 1 of which is the subgenus Pavlovskyella. In an earlier study, we identified distinct soft ticks as Ornithodoros (Pavlovskyella) sp., which were collected from animal shelters in Khyber Pakhtunkhwa, Pakistan. Providing additional collections from that same locality and a comprehensive analysis involving detailed morphological and mitogenome-based comparisons with closely related species, this study formally designates a novel species for these specimens. Adults and late-instar nymphs of the new species display a dorsoventral groove, small cheeks not covering the capitulum, 5 small even humps on tarsus I and a transverse postanal groove intersecting the median postanal groove perpendicularly. It also lacks a tuft of setae on the ventral surface of the hood which separates the novel species from Ornithodoros papillipes. Ventral chaetotaxy of tarsus IV indicates 4–7 setal pairs in nymphs and 5–7 pairs in adults that separate the new species from Ornithodoros tholozani sensu stricto and Ornithodoros crossi, 2 morphologically closely related species that occur in geographical proximity. Phylogenetic analyses of the full-length mitochondrial genome and the 18S and 28S ribosomal RNA genes, combined with pairwise nucleotide comparisons of cox1, cox2, atp8, atp6, cox3, nad3, nad5, nad4, nad4L, nad6, cytb, nad1, nad2, 12S rDNA, 16S rDNA, 18S rDNA and 28S rDNA further support that the new species belongs to the Pavlovskyella subgenus, clustering with O. tholozani, Ornithodoros verrucosus and Ornithodoros tartakovskyi.

First report of Anaplasma spp., Ehrlichia spp., and Rickettsia spp. in Amblyomma gervaisi ticks infesting monitor lizards (Varanus begalensis) of Pakistan

Ticks pose significant health risks to both wildlife and humans due to their role as vectors for various pathogens. In this study, we investigated tick infestation patterns, tick-associated pathogens, and genetic relationships within the tick species Amblyomma gervaisi, focusing on its prevalence in monitor lizards (Varanus bengalensis) across different districts in Pakistan. We examined 85 monitor lizards and identified an overall mean intensity of 19.59 ticks per infested lizard and an overall mean abundance of 11.98 ticks per examined lizard. All collected ticks (n = 1019) were morphologically identified as A. gervaisi, including 387 males, 258 females, 353 nymphs, and 21 larvae. The highest tick prevalence was observed in the Buner district, followed by Torghar and Shangla, with the lowest prevalence in Chitral. Lizard captures primarily occurred from May to October, correlating with the period of higher tick infestations. Molecular analysis was conducted on tick DNA, revealing genetic similarities among A. gervaisi ticks based on 16S rDNA and ITS2 sequences. Notably, we found the absence of A. gervaisi ITS2 sequences in the NCBI GenBank, highlighting a gap in existing genetic data. Moreover, our study identified the presence of pathogenic microorganisms, including Ehrlichia sp., Candidatus Ehrlichia dumleri, Anaplasma sp., Francisella sp., Rickettsia sp., and Coxiella sp., in these ticks. BLAST analysis revealed significant similarities between these pathogenic sequences and known strains, emphasizing the potential role of these ticks as vectors for zoonotic diseases. Phylogenetic analyses based on nuclear ITS2 and mitochondrial 16S rDNA genes illustrated the genetic relationships of A. gervaisi ticks from Pakistan with other Amblyomma species, providing insights into their evolutionary history. These findings contribute to our understanding of tick infestation patterns, and tick-borne pathogens in monitor lizards, which has implications for wildlife health, zoonotic disease transmission, and future conservation efforts. Further research in this area is crucial for a comprehensive assessment of the risks associated with tick-borne diseases in both wildlife and humans.

First record of Alectorobius coniceps (Ixodoidea: Argasidae) and Dermacentor sp. (Ixodoidea: Ixodidae) in Pakistan

Alectorobius species are soft ticks primarily infesting birds, such as swallows, while Dermacentor species are hard ticks mainly infesting mammals, such as small ruminants. This study for the first time reported on the morphological and molecular bases of two tick species, namely A. coniceps and a Dermacentor sp. in Pakistan. The former species was examined in swallows' nests in Khyber Pakhtunkhwa province, while the latter species was examined in small ruminants in Balochistan province. In total, 25 ticks were collected, with 14 ticks morphologically identified as A. coniceps (males = 9 and females = 5) and 11 ticks identified as Dermacentor sp. (males = 7 and females = 4). Following morphological identification, molecular identification was gained by obtaining 16S rDNA and cox1 sequences for these ticks. The BLAST results for the 16S rDNA and cox1 sequences from A. coniceps shared a maximum identity of 97.46% and 96.49% with the same species from Malta. The BLAST analysis of the 16S rDNA and cox1 sequences from Dermacentor sp. showed maximum identities of 98.42% and 97.45% with Dermacentor pavlovskyi from China. The phylogenetic analysis based on 16S rDNA and cox1 of A. coniceps showed a close evolutionary relationship with the same species. The case of Dermacentor sp., based on 16S DNA and cox1, indicated a close evolutionary relationship with Dermacentor pavlovskyi from China.

First report of spotted fever group Rickettsia aeschlimannii in Hyalomma turanicum, Haemaphysalis bispinosa, and Haemaphysalis montgomeryi infesting domestic animals: updates on the epidemiology of tick-borne Rickettsia aeschlimannii

Tick-borne Rickettsia spp. have long been known as causative agents for zoonotic diseases. We have previously characterized Rickettsia spp. in different ticks infesting a broad range of hosts in Pakistan; however, knowledge regarding Rickettsia aeschlimannii in Haemaphysalis and Hyalomma ticks is missing. This study aimed to obtain a better understanding about R. aeschlimannii in Pakistan and update the knowledge about its worldwide epidemiology. Among 369 examined domestic animals, 247 (66%) were infested by 872 ticks. Collected ticks were morphologically delineated into three genera, namely, Rhipicephalus, Hyalomma, and Haemaphysalis. Adult females were the most prevalent (number ₌ 376, 43.1%), followed by nymphs (303, 34.74%) and males (193, 22.13%). Overall, genomic DNA samples of 223 tick were isolated and screened for Rickettsia spp. by the amplification of rickettsial gltA, ompA, and ompB partial genes using conventional PCR. Rickettsial DNA was detected in 8 of 223 (3.58%) ticks including nymphs (5 of 122, 4.0%) and adult females (3 of 86, 3.48%). The rickettsial gltA, ompA, and ompB sequences were detected in Hyalomma turanicum (2 nymphs and 1 adult female), Haemaphysalis bispinosa (1 nymph and 1 adult female), and Haemaphysalis montgomeryi (2 nymphs and 1 adult female). These rickettsial sequences showed 99.71-100% identity with R. aeschlimannii and phylogenetically clustered with the same species. None of the tested Rhipicephalus microplus, Hyalomma isaaci, Hyalomma scupense, Rhipicephalus turanicus, Hyalomma anatolicum, Rhipicephalus haemaphysaloides, Rhipicephalus sanguineus, Haemaphysalis cornupunctata, and Haemaphysalis sulcata ticks were found positive for rickettsial DNA. Comprehensive surveillance studies should be adopted to update the knowledge regarding tick-borne zoonotic Rickettsia species, evaluate their risks to humans and livestock, and investigate the unexamined cases of illness after tick bite among livestock holders in the country.

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.

Molecular evidence of Borrelia theileri and closely related Borrelia spp. in hard ticks infesting domestic animals

Ticks pose significant threats to hosts by transmitting Borrelia spp., which are grouped into Lyme borreliae, relapsing fever borreliae (RF), and reptiles- and monotremes-associated borreliae. The RF borreliae encompass a group of Borrelia species predominantly transmitted by soft ticks, but some of its members can also be transmitted by hard ticks. Information on the detection and genetic characterization of tick-borne RF borreliae, including Borrelia theileri, is notably rare in Asia, particularly in Pakistan. Herein, we employed molecular techniques to detect borreliae in hard ticks collected from domestic animals in Khyber Pakhtunkhwa, Pakistan. Ticks were subjected to morphological analysis, followed by DNA extraction and PCR amplification of partial fragments of borrelial 16S rRNA and flaB genes. A total of 729 ticks were collected from 264 hosts, with Haemaphysalis cornupunctata (12.9%; 94/729) being the most prevalent, followed by Hyalomma anatolicum (11.7%; 85/729), Rhipicephalus microplus (10.0%; 73/729), Haemaphysalis kashmirensis (9.1%; 66/729), Haemaphysalis bispinosa (8.5%; 62/729), Rhipicephalus sanguineus (8%; 58/729), Haemaphysalis montgomeryi (6.2%; 45/729), Rhipicephalus turanicus (5.5%; 40/729), Hyalomma dromedarii and Ixodes kashmirensis (4.4%; 32/729 each), Rhipicephalus haemaphysaloides (4.1%; 30/729), Haemaphysalis sulcata and Hyalomma scupense (3.8%; 28/729 each), Haemaphysalis danieli (2.9%; 21/729), Hyalomma kumari (2.6%; 19/729), and Hyalomma isaaci (2.2%; 16/729). Based on 16S rRNA detection of Borrelia spp., only R. turanicus yielded positive results, resulting in an overall infection rate of 0.3% (2/160), while using flaB-based detection, four tick species including R. microplus, R. turanicus, Ha. sulcata, and Ha. cornupunctata showed positive results, yielding an overall infection rate of 6.9% (11/160). The amplified DNA fragments of borrelial 16S rRNA and flaB in R. turanicus from goats shared maximum identities of 100 and 99.40% with Borrelia theileri, respectively. Amplified borrelial flaB fragments in R. microplus from cows and sheep displayed 100% identity with B. theileri, while flaB fragments in Ha. cornupunctata and Ha. sulcata from goats revealed identities of 99.32 and 99.75% with undetermined RF Borrelia spp., respectively. Phylogenetic analysis revealed clustering of B. theileri from R. microplus and R. turanicus with the same species, while Borrelia spp. from Ha. cornupunctata and Ha. sulcata with undetermined RF Borrelia spp. Notably, this research marks the first documentation of B. theileri in R. turanicus and the identification of RF Borrelia spp. in Ha. cornupunctata and Ha. sulcata.

Molecular Survey of Rickettsia raoultii in Ticks Infesting Livestock from Pakistan with Notes on Pathogen Distribution in Palearctic and Oriental Regions

Ticks are hematophagous ectoparasites that transmit different pathogens such as Rickettsia spp. to domestic and wild animals as well as humans. Genetic characterizations of Rickettsia spp. from different regions of Pakistan are mostly based on one or two genetic markers and are confined to small sampling areas and limited host ranges. Therefore, this study aimed to molecularly screen and genetically characterize Rickettsia spp. in various tick species infesting camels, sheep, and goats. All the collected tick specimens were morphologically identified, and randomly selected tick species (148) were screened molecularly for the detection of Rickettsia spp. by amplifying three rickettsial DNA fragments, namely, the citrate-synthase gene (gltA), outer-membrane protein A (ompA), and outer-membrane protein B (ompB). After examining 261 hosts, 161 (61.7%) hosts were found infested by 564 ticks, including 287 (50.9%) nymphs, 171 (30.3%) females, and 106 (18.8%) males in five districts (Kohat, Dera Ismail Khan, Lower Dir, Bajaur, and Mansehra). The highest occurrence was noted for Hyalomma dromedarii (number = 72, 12.8%), followed by Haemaphysalis sulcata (n = 70, 12.4%), Rhipicephalus turanicus (n = 64, 11.3%), Rhipicephalus microplus (n = 55, 9.7%), Haemaphysalis cornupunctata (n = 49, 8.7%), Hyalomma turanicum (n = 48, 8.5%), Hyalomma isaaci (n = 45, 8.0%), Haemaphysalis montgomeryi (n = 44, 7.8%), Hyalomma anatolicum (n = 42, 7.5%), Haemaphysalis bispinosa (n = 38, 6.7%), and Rhipicephalus haemaphysaloides (n = 37, 6.6%). A subset of 148 ticks were tested, in which eight (5.4%) ticks, including four Hy. turanicum, two Ha. cornupunctata, one Ha. montgomeryi, and one Ha. bispinosa, were found positive for Rickettsia sp. The gltA, ompA, and ompB sequences revealed 100% identity and were phylogenetically clustered with Rickettsia raoultii reported in China, Russia, USA, Turkey, Denmark, Austria, Italy, and France. Additionally, various reports on R. raoultii from Palearctic and Oriental regions were summarized in this study. To the best of our knowledge, this is the first report regarding genetic characterization and phylogenetic analysis of R. raoultii from Pakistan. Further studies to investigate the association between Rickettsia spp. and ticks should be encouraged to apprise effective management of zoonotic consequences.

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.

Molecular Detection of Rickettsia hoogstraalii in Hyalomma anatolicum and Haemaphysalis sulcata: Updated Knowledge on the Epidemiology of Tick-Borne Rickettsia hoogstraalii

Ticks are hematophagous ectoparasites that transmit pathogens to animals and humans. Updated knowledge regarding the global epidemiology of tick-borne Rickettsia hoogstraalii is dispersed, and its molecular detection and genetic characterization are missing in Pakistan. The current study objectives were to molecularly detect and genetically characterize Rickettsia species, especially R. hoogstraalii, in hard ticks infesting livestock in Pakistan, and to provide updated knowledge regarding their global epidemiology. Ticks were collected from livestock, including goats, sheep, and cattle, in six districts of Khyber Pakhtunkhwa (KP) Pakistan. Overall, 183 hosts were examined, of which 134 (73.2%), including goats (number = 39/54, 72.2%), sheep (23/40, 57.5%), and cattle (71/89, 80%) were infested by 823 ticks. The most prevalent tick species was Rhipicephalus microplus (number = 283, 34.3%), followed by Hyalomma anatolicum (223, 27.0%), Rhipicephalus turanicus (122, 14.8%), Haemaphysalis sulcata (104, 12.6%), Haemaphysalis montgomeryi (66, 8.0%), and Haemaphysalis bispinosa (25, 3.03%). A subset of 210 ticks was selected and screened for Rickettsia spp. using PCR-based amplification and subsequent sequencing of rickettsial gltA and ompB fragments. The overall occurrence rate of R. hoogstraalii was 4.3% (number = 9/210). The DNA of Rickettsia was detected in Hy. anatolicum (3/35, 8.5%) and Ha. sulcata (6/49, 12.2%). However, no rickettsial DNA was detected in Rh. microplus (35), Rh. turanicus (35), Ha. montgomeryi (42), and Ha. bispinosa (14). The gltA and ompB fragments showed 99-100% identity with R. hoogstraalii and clustered phylogenetically with the corresponding species from Pakistan, Italy, Georgia, and China. R. hoogstraalii was genetically characterized for the first time in Pakistan and Hy. anatolicum globally. Further studies should be encouraged to determine the role of ticks in the maintenance and transmission of R. hoogstraalii in different hosts.

First report on detection of Hepatozoon ayorgbor in Rhipicephalus haemaphysaloides and Hepatozoon colubri in Haemaphysalis sulcata and Hyalomma anatolicum: risks of spillover of Hepatozoon spp. from wildlife to domestic animals

This study aimed to detect Hepatozoon spp. in ticks infesting asymptomatic domestic animals and to provide insight into their potential spillover from wild to domestic animals. In total, 537 tick specimens were collected in Khyber Pakhtunkhwa, Pakistan, and morphologically identified. The most prevalent tick species was Haemaphysalis cornupunctata (69; 12.8%), followed by Haemaphysalis kashmirensis (62; 11.5%), Rhipicephalus microplus (58; 10.8%), Haemaphysalis montgomeryi (51; 9.5%), Rhipicephalus sanguineus (49; 9.1%), each Haemaphysalis bispinosa and Haemaphysalis sulcata (43; 8.0%), each Hyalomma anatolicum and Rhipicephalus turanicus (37; 6.9%), Rhipicephalus haemaphysaloides (33; 6.1%) Hyalomma scupense (30; 5.6%), and Hyalomma isaaci (25; 4.7%). The extracted DNA from a subset of each tick species was subjected to PCR to amplify 18S rRNA fragments of Hepatozoon spp. By BLAST analysis, the Hepatozoon sp. detected in Hy. anatolicum infesting cows and in Ha. sulcata infesting sheep showed 99.7% maximum identity with Hepatozoon colubri. Similarly, the Hepatozoon sp. detected in R. haemaphysaloides infesting goats shared 99.49% maximum identity with Hepatozoon ayorgbor, and the Hepatozoon sp. detected in R. sanguineus infesting dogs exhibited 99.7% identity with Hepatozoon canis. Having an overall infection rate (9.3%; 16/172), the highest infection rate was recorded for each H. canis, and H. colubri (3.5%; 6/172), followed by H. ayorgbor (2.3%; 4/172). In the phylogenetic tree, H. colubri clustered with corresponding species from Iran, H. ayorgbor clustered with the same species from Croatia, Ghana, and Portugal, and H. canis clustered with the conspecifics from Iran, Israel, Romania, and Zambia. Regarding the potential spillover of Hepatozoon spp. from wildlife through ticks, free ranging animals was at higher risk compared to confined animals (RR = 3.05), animals consuming food from wildlife habitats were at higher risk compared to those consuming domestic food (RR = 3.06), and animals residing in farm buildings located in wildlife habitats were at higher risk compared to those residing in farm buildings located in villages (RR = 3.28). In addition to the first report on H. canis in R. sanguineus in Pakistan, this is the earliest data showing H. ayorgbor in R. haemaphysaloides and H. colubri in Ha. sulcata and Hy. anatolicum. These preliminary findings suggest a potential spillover of Hepatozoon spp. from wild to domestic animals via ticks under certain risk factors.

First report on tick-borne pathogens detected in ticks infesting stray dogs near butcher shops

Public health is a major concern for several developing countries due to infectious agents transmitted by hematophagous arthropods such as ticks. Health risks due to infectious agents transmitted by ticks infesting butcher-associated stray dogs (BASDs) in urban and peri-urban regions have been neglected in several developing countries. To the best of the authors' knowledge, this is the first study assessing public health risks due to ticks infesting BASDs in Pakistan's urban and peri-urban areas. A total of 575 ticks (390 from symptomatic and 183 from asymptomatic BASDs) were collected from 117 BASDs (63 symptomatic and 54 asymptomatic); the ticks belonged to 4 hard tick species. A subset of each tick species' extracted DNA was subjected to polymerase chain reaction (PCR) to amplify the 16S rDNA and cox1 sequences of the reported tick species, as well as bacterial and protozoal agents. The ticks' 16S rDNA and cox1 sequences showed 99-100% identities, and they were clustered with the sequence of corresponding species from Pakistan and other countries in phylogenetic trees. Among the screened 271 ticks' DNA samples, Anaplasma spp. was detected in 54/271 (19.92%) samples, followed by Ehrlichia spp. (n = 40/271, 14.76%), Rickettsia spp. (n = 33/271, 12.17%), Coxiella spp. (n = 23/271, 4.48%), and Hepatozoon canis (n = 9/271, 3.32%). The obtained sequences and phylogenetic analyzes revealed that the pathogens detected in ticks were Ehrlichia minasensis, Ehrlichia sp., Hepatozoon canis, Coxiella burnetii, Coxiella sp., Anaplasma capra, Anaplasma platys, Anaplasma sp., Rickettsia massiliae, "Candidatus Rickettsia shennongii" and Rickettsia aeschlimannii. Tick-borne pathogens such as E. minasensis, H. canis, A. capra, A. platys, and R. aeschlimannii, were detected based on the DNA for the first time in Pakistan. This is the first report on public health risks due to ticks infesting BASDs. These results not only provided insights into the occurrence of novel tick-borne pathogens in the region but also revealed initial evidence of zoonotic threats to both public health and domestic life.

Characterisation and comparative genomics of three new Varanus-associated Borrelia spp. from Indonesia and Australia

BACKGROUND

Borrelia are important disease-causing tick- and louse-borne spirochaetes than can infect a wide variety of vertebrates, including humans and reptiles. Reptile-associated (REP) Borrelia, once considered a peculiarity, are now recognised as a distinct and important evolutionary lineage, and are increasingly being discovered worldwide in association with novel hosts. Numerous novel Borrelia spp. associated with monitor lizards (Varanus spp.) have been recently identified throughout the Indo-Pacific region; however, there is a lack of genomic data on these Borrelia.

METHODS

We used metagenomic techniques to sequence almost complete genomes of novel Borrelia spp. from Varanus varius and Varanus giganteus from Australia, and used long- and short-read technologies to sequence the complete genomes of two strains of a novel Borrelia sp. previously isolated from ticks infesting Varanus salvator from Indonesia. We investigated intra- and interspecies genomic diversity, including plasmid diversity and relatedness, among Varanus-associated Borrelia and other available REP Borrelia and, based on 712 whole genome orthologues, produced the most complete phylogenetic analysis, to the best of our knowledge, of REP Borrelia to date.

RESULTS

The genomic architecture of Varanus-associated Borrelia spp. is similar to that of Borrelia spp. that cause relapsing fever (RF), and includes a highly conserved megaplasmid and numerous smaller linear and circular plasmids that lack structural consistency between species. Analysis of PF32 and PF57/62 plasmid partitioning genes indicated that REP Borrelia plasmids fall into at least six distinct plasmid families, some of which are related to previously defined Borrelia plasmid families, whereas the others appear to be unique. REP Borrelia contain immunogenic variable major proteins that are homologous to those found in Borrelia spp. that cause RF, although they are limited in copy number and variability and have low sequence identities to RF variable major proteins. Phylogenetic analyses based on single marker genes and 712 single copy orthologs also definitively demonstrated the monophyly of REP Borrelia as a unique lineage.

CONCLUSIONS

In this work we present four new genomes from three novel Borrelia, and thus double the number of REP Borrelia genomes publicly available. The genomic characterisation of these Borrelia clearly demonstrates their distinctiveness as species, and we propose the names Borrelia salvatorii, 'Candidatus Borrelia undatumii', and 'Candidatus Borrelia rubricentralis' for them.

Association of SFG Rickettsia massiliae and Candidatus Rickettsia shennongii with Different Hard Ticks Infesting Livestock Hosts

Ixodid ticks are responsible for the transmission of various intracellular bacteria, such as the Rickettsia species. Little Information is available about the genetic characterization and epidemiology of Rickettsia spp. The current study was designed to assess the tick species infesting various livestock hosts and the associated Rickettsia spp. in Pakistan. Ticks were collected from different livestock hosts (equids, cattle, buffaloes, sheep, goats, and camels); morphologically identified; and screened for the genetic characterization of Rickettsia spp. by the amplification of partial fragments of the gltA, ompA and ompB genes. Altogether, 707 ticks were collected from 373 infested hosts out of 575 observed hosts. The infested hosts comprised 105 cattle, 71 buffaloes, 70 sheep, 60 goats, 34 camels, and 33 equids. The overall occurrence of Rickettsia spp. was 7.6% (25/330) in the tested ticks. Rickettsia DNA was detected in Rhipicephalus haemaphysaloides (9/50, 18.0%), followed by Rhipicephalus turanicus (13/99, 13.1%), Haemaphysalis cornupunctata (1/18, 5.5%), and Rhipicephalus microplus (2/49, 4.1%); however, no rickettsial DNA was detected in Hyalomma anatolicum (71), Hyalomma dromedarii (35), and Haemaphysalis sulcata (8). Two Rickettsia agents were identified based on partial gltA, ompA, and ompB DNA sequences. The Rickettsia species detected in Rh. haemaphysaloides, Rh. turanicus, and Rh. microplus showed 99-100% identity with Rickettsia sp. and Candidatus Rickettsia shennongii, and in the phylogenetic trees clustered with the corresponding Rickettsia spp. The Rickettsia species detected in Rh. haemaphysaloides, Rh. turanicus, Rh. microplus, and Ha. cornupunctata showed 100% identity with R. massiliae, and in the phylogenetic trees it was clustered with the same species. Candidatus R. shennongii was characterized for the first time in Rh. haemaphysaloides, Rh. turanicus, and Rh. microplus. The presence of SFG Rickettsia spp., including the human pathogen R. massiliae, indicates a zoonotic risk in the study region, thus stressing the need for regular surveillance.

Molecular characterization of Amblyomma geoemydae using CO1 mitochondrial gene to validate phenotypic taxonomical evaluation

Animal ectoparasites are linked to the spread of serious medical and veterinary important pathogens. Our research intends to close the knowledge gap concerning the numerous ectoparasites that inhabit animals in Wayanad. Ectoparasites in animals brought to the veterinary dispensaries in Wayanad were retrieved and identified morphologically and molecularly. Using a high-quality stereomicroscope, the taxonomic features of the four following species were examined and identified: Haemaphysalis bispinosa, Rhipicephalus annulatus, R. microplus, and Amblyomma geoemydae. The important disease vector A. geoemydae was reported for the first time in Kerala. The important phenotypic characters of the highlighted species A. geoemydae are the edge of the basis capituli is circular without cornua, and the hypostomal dental formula is 2/2. The taxonomically identified four species were subjected to CO1 gene sequence analysis. The evolutionary relationship was inspected through the neighbour-joining method, and the phylogenetic tree was built through the Maximum Likelihood method. The present study has also estimated the diversity index of R. microplus, R. annulatus, H. bispinosa, and A. geoemydae. Among them, R. microplus 0.36638 have reported with the maximum diversity index score. The significance of the study is the presence of Lyme disease vector A. geoemydae, in the Wayanad District of Kerala, and it is the first report of the species from where an outbreak of Lyme disease occurred in 2013.

Assessment of cypermethrin and amitraz resistance and molecular profiling of voltage-gated sodium channel and octopamine tyramine genes of Rhipicephalus microplus

Control of ticks and tick-borne pathogens is a priority for human and animal health. Livestock-holders extensively rely on acaricide applications for tick control. Different groups of acaricides including cypermethrin and amitraz have been consistently used in Pakistan. There has been a gap in understanding the susceptibility or resistance of Rhipicephalus microplus, the most prevalent tick in Pakistan, to acaricides. The present study aimed to molecularly characterize cypermethrin and amitraz targeted genes such as voltage-gated sodium channel (VGSC) and octopamine tyramine (OCT/Tyr) of R. microplus ticks in Khyber Pakhtunkhwa (KP), Pakistan to monitor the acaricides resistance. Tick specimens were collected from cattle and buffaloes in northern (Chitral, Shangla, Swat, Dir, and Buner), central (Peshawar, Mardan, Charsadda, Swabi, and Nowshera), and southern districts (Kohat, Karak, Lakki Marwat, Tank, and Dera Ismail Khan) of KP, Pakistan. Different concentrations of commercially available cypermethrin (10%) and amitraz (12.5%) were prepared for in vitro larval immersion tests (LIT). In LIT, the average mortality rate of immersed larvae was recorded that was increased gradually with an increase in the concentration of specific acaricide. The larvae's highest mortality rates (94.5% and 79.5%) were observed at 100-ppm of cypermethrin and amitraz, respectively. A subset of 82 R. microplus ticks was subjected to extract genomic DNA, followed by PCR to amplify partial fragments of VGSC (domain-II) and OCT/Tyr genes. The BLAST results of the consensus sequence of VGSC gene (domain-II) showed 100% identity with the acaricides susceptible tick sequence from the United States (reference sequence). Obtained identical sequences of OCT/Tyr genes showed maximum identity (94-100%) with the identical sequences reported from Australia (reference sequence), India, Brazil, Philippines, USA, South Africa, and China. Thirteen single nucleotide polymorphisms (10 synonymous and three non-synonymous) were observed at various positions of partial OCT/Tyr gene fragments. The SNP at position A-22-C (T-8-P) in OCT/Tyr gene has been linked to amitraz resistance in R. microplus ticks. Molecular analysis and LIT bioassay's findings indicate the availability of resistant R. microplus ticks in the KP region. To our understanding, this is the first preliminary study to monitor cypermethrin and amitraz resistance via molecular profiling of cypermethrin and amitraz targeted genes (VGSC and OCT/Tyr) in combination with in vitro bioassays (LIT) in R. microplus ticks from Pakistan.

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.

Phylogenetic Position of Haemaphysalis kashmirensis and Haemaphysalis cornupunctata, with Notes on Rickettsia spp

Despite high diversity in the Oriental region, ticks of the genus Haemaphysalis have been neglected regarding their genetic data and vector potential. This study aimed to genetically characterize three species of the genus Haemaphysalis: Haemaphysalis cornupunctata, Haemaphysalis kashmirensis, and Haemaphysalis montgomeryi infesting goats and sheep, and Rickettsia spp. associated with these tick species in the Hindu Kush Himalayan range of Pakistan. Altogether, 834 ticks were collected by examining 120 hosts including goats (64/120, 53.3%) and sheep (56/120, 46.6%), in which 86 (71.6%) hosts were found to be tick-infested. The morphologically identified ticks were subjected to DNA extraction and PCR for the amplification of partial 16S rDNA and cox fragments. Rickettsia spp. associated with the collected ticks were detected through the amplification of gltA, ompA and ompB partial fragments. The 16S rDNA of H. cornupunctata and H. montgomeryi showed a maximum identity of 100% with the sequences of the same species, whereas the 16S rDNA of H. kashmirensis showed the highest identity of 93-95% with Haemaphysalis sulcata. The cox sequence of H. montgomeryi displayed 100% identity with the same species. In comparison, the cox sequences of H. cornupunctata and H. kashmirensis showed maximum identities of 87.65-89.22% with Haemaphysalis punctata and 89.34% with H. sulcata, respectively. The gltA sequence of Rickettsia sp. from H. kashmirensis showed the highest identity of 97.89% with Rickettsia conorii subsp. raoultii, while the ompA and ompB fragments from the same DNA samples revealed 100% and 98.16% identity with Rickettsia sp. and "Candidatus Rickettsia longicornii", respectively. Another gltA sequence amplified from H. montgomeryi ticks showed 100% identity with Rickettsia hoogstraalii, while the attempts to amplify ompA and ompB for R. hoogstraalii were unsuccessful. In the phylogenetic tree, the 16S rDNA of H. cornupunctata clustered with the corresponding species while its cox clustered with H. punctata. Both 16S rDNA and cox sequences of H. kashmirensis clustered with H. sulcata. The gltA sequence of Rickettsia sp. was clustered individually in the spotted fever (SF) group of Rickettsia, while the gltA sequence of R. hoogstraalii was clustered with the same species in the transition group of Rickettsia. In the SF group, the rickettsial ompA and ompB sequence clustered with undetermined Rickettsia sp. and "Candidatus Rickettsia longicornii", respectively. This is the earliest study regarding the genetic characterization of H. kashmirensis. This study indicated that ticks belong to the genus Haemaphysalis have the potential of harboring and/or transmitting Rickettsia spp. in the region.