Ticks transmit pathogens to animals and humans more often than any other arthropod vector. The rural economy of Pakistan mainly depends on livestock farming, and tick infestations cause severe problems in this sector. The present study aimed to molecularly characterize the Anaplasma spp. in hard ticks collected from six districts of Khyber Pakhtunkhwa, Pakistan. Ticks were collected from various livestock hosts, including cattle breeds (Holstein-Friesian, Jersey, Sahiwal, and Achai), Asian water buffaloes, sheep, and goats from March 2018 to February 2019. Collected ticks were morphologically identified and subjected to molecular screening of Anaplasma spp. by amplifying 16S rDNA sequences. Six hundred seventy-six ticks were collected from infested hosts (224/350, 64%). Among the nine morphologically identified tick species, the highest occurrence was noted for Rhipicephalus microplus (254, 37.6%), followed by Hyalomma anatolicum (136, 20.1%), Rhipicephalus haemaphysaloides (119, 17.6%), Rhipicephalus turanicus (116, 17.1%), Haemaphysalis montgomeryi (14, 2.1%), Hyalomma dromedarii (11, 1.6%), Haemaphysalis bispinosa (10, 1.5%), Hyalomma scupense (8, 1.2%), and Haemaphysalis kashmirensis (8, 1.2%). The occurrence of tick females was highest (260, 38.5%), followed by nymphs (246, 36.4%) and males (170, 25.1%). Overall, the highest occurrence of ticks was recorded in the Peshawar district (239, 35.3%), followed by Mardan (183, 27.1%), Charsadda (110, 16.3%), Swat (52, 7.7%), Shangla (48, 7.1%), and Chitral (44, 6.5%). Among these ticks, Anaplasma marginale was detected in R. microplus, R. turanicus, and R. haemaphysaloides. The 16S rDNA sequences showed high identity (98–100%) with A. marginale reported from Australia, China, Japan, Pakistan, Thailand, Uganda, and the USA. In phylogenetic analysis, the sequence of A. marginale clustered with the same species reported from Australia, China, Pakistan, Thailand, Uruguay, and the USA. Further molecular work regarding the diversity of tick species and associated pathogens is essential across the country.

Hard ticks (Ixodida: Ixodidae) are medically important ectoparasites that feed on all classes of terrestrial vertebrates. Recently, we molecularly characterized hard ticks and associated Anaplasma spp. in the northern and central regions of Khyber Pakhtunkhwa (KP), Pakistan; however, this knowledge was missing in the southern regions. This study aimed to investigate tick prevalence, host range, genetic diversity, and molecular survey of Anaplasma spp. in a wide range of tick species in two distinct physiographic regions of southern KP. A total of 1873 hard ticks were randomly collected from 443/837 hosts (cattle, Asian water buffaloes, horses, goats, sheep, dogs, and camels) in Lakki Marwat, Bannu, and Orakzai districts of KP. Overall, 12 tick species were morphologically identified, among which Hyalomma dromedarii was the most prevalent species (390/1873, 20.9%), followed by Hy. anatolicum (294, 15.7%), Rhipicephalus microplus (262, 14%), Hy. scupense (207, 11.1%), R. sanguineus (136, 7.3%), R. turanicus (121, 6.5%), Haemaphysalis cornupunctata (107, 5.7%), R. haemaphysaloides (110, 5.9%), Ha. montgomeryi (87, 4.6%), Hy. isaaci (58, 3.1%), Ha. bispinosa (54, 2.9%), and Ha. sulcata (47, 2.5%). The extracted DNA from a subset of each tick species was subjected to PCR to amplify cox1 or 16S rRNA sequences of ticks and 16S rRNA sequences of Anaplasma spp. The tick cox1 sequences showed 99–100% identities with the sequences of the same species, whereas 16S rRNA sequences of R. turanicus, Ha. montgomeryi and Ha. sulcata showed 97–100% identities with the corresponding species. The 16S rRNA sequence of Ha. cornupunctata showed 92% identity with the species from the same subgenus, such as Ha. punctata. The 16S rRNA sequence of Anaplasma spp. showed 100% identity with Anaplasma marginale. Moreover, 54 ticks were found positive for A. marginale with a total infection rate of 17.2%. The highest infection rate was recorded in Hy. dromedarii (31.1%) and the lowest in each R. haemaphysaloides and R. sanguineus (20%). All the cox1 or 16S rRNA sequences in phylogenetic trees clustered with the same species, except Ha. cornupunctata, which clustered with the Ha. (Aboimisalis) punctata. In this study, Ha. cornupunctata was reported for the first time at the molecular level. The genetic characterization of ixodid ticks and molecular detection of associated A. marginale will assist in the epidemiological surveillance of these parasites in the region.

Ticks are ectoparasites that act as vectors for transmission of various pathogens to wild and domesticated animals and pose a serious threat to human health. Because of the hot and humid conditions in different agro-ecological zones of Pakistan, ticks are abundant and parasitize a variety of animals. The aim of this study was to identify different tick species and distribution on different hosts especially livestock, such as sheep, goat, cattle, buffalo, and camel, and livestock associated canines and equines, such as horse, donkey, and dog, across different agro-ecological zones of Pakistan. The ticks samples were collected and morphologically identified at genus and species level using morphological keys under stereomicroscope. A total of 2,846 animals were examined for the tick infestation, and 408 animals were tick-infested. Eleven tick species belonging to 4 genera were identified: Hyalomma anatolicum, Hyalomma scupense, Hyalomma dromedarii, Hyalomma isaaci, Rhipicephalus microplus, Rhipicephalus haemaphysaloides, Rhipicephalus turanicus, Haemaphysalis cornupunctata, Haemaphysalis montgomeryi, Haemaphysalis bispinosa, and Ixodes kashmiricus. The overall tick prevalence was 14.3%; host-wise infestation rate was 12.2% in sheep; 12.6%, goat; 11.7%, buffalo; 11.7%, cattle; 19.6%, camel; 27.4%, donkey; 23.5%, horse; and 24.3%, dog. Tick infestation of different animals differed on the basis of the zones. Camels showed the highest tick infestation rate in zones 1 and 2 (21.4 and 26.7%, respectively), whereas donkeys showed the highest infestation rate in zones 3, 4, 6, and 7 (25, 39.3, 3.3, and 21.4%, respectively). The infestation rates of Hyalomma and Rhipicephalus were the highest in zone 2 (71.4 and 52.9%, respectively). The infestation rate of Hyalomma was the highest (47.4%) in sheep; Haemaphysalis (46.9%), goat; Rhipicephalus (69.7%), buffalo; Rhipicephalus (62.3%), cattle; Hyalomma (70%), camel; Ixodes (60.9%), donkey; Ixodes (75%), horse; and Rhipicephalus (61.1%), dog. This study showed the diversity and infestation rate of different ticks with respect to their hosts and agro-ecological zones of Pakistan. High tick burdens and infestation rates are responsible for the spread of different tick-borne infections, resulting in loss of animal productivity and posing a threat to animal and human health. Understanding different tick species and their distribution across different zones will be helpful for developing efficient control strategies against different tick born infections.

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.

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.