Skip to main content

Risk factors associated with tick infestations on equids in Khyber Pakhtunkhwa, Pakistan, with notes on Rickettsia massiliae detection



Studies on ticks infesting equids are lacking in various parts of the world, including Khyber Pakhtunkhwa (KP), Pakistan. The aim of this study was to investigate the diversity of ticks infesting equids, associated risk factors and rickettsial detection in ticks from equids in KP.


Inspection of 404 equid hosts from November 2018 to October 2019 resulted in the collection of 550 ticks. Data on tick-associated risk factors were collected from equid owners by means of a questionnaire. After morphological identification, partial DNA sequences of the tick mitochondrial 16S rRNA gene were used for taxonomic confirmation of species. Partial sequences of the gltA and ompA genes were used for Rickettsia detection in ticks.


A total of 550 tick specimens were collected on 324 (80.2%) of the equids inspected, of which 161 were horses (50%), 145 (45%) were donkeys and 18 were mules (5%). The ticks were identified as belonging to the following five species: Rhipicephalus microplus (341 specimens, 62% of the total ticks), Rh. haemaphysaloides (126, 23%), Rh. turanicus (39, 7%), Rh. sanguineus (s.l.) (33, 6%) and Hyalomma anatolicum (11, 2%). The most prevalent tick life stage was adult females (279, 51%) followed by adult males (186, 34%) and nymphs (85, 15%). Higher tick infestations were observed on male equids (relative risk [RR] 0.7432, P < 0.0005) and adult equids (RR 1.268, P < 0.0020). Ticks were frequently attached to the axial region of horses (55, 21%), sternum of donkeys (44, 21%) and belly of mules (19, 23%) (P < 0.04). Temporal patterns of tick infestation in association with temperature and humidity were highly significant (P < 0.05). Risk factors, such as animal housing (P < 0.0003), living management (P < 0.006), grazing type (P < 0.01) and location in hilly areas (P < 0.02), significantly enhanced the chances for tick infestation. Tick species analyzed in this study were phylogenetically related to species from Afghanistan, China, South Africa and Taiwan. Partial sequences of the gltA and ompA genes obtained from Rh. microplus and Rh. haemaphysaloides were 100% identical to the spotted fever group pathogen Rickettsia massiliae.


Equids exposed to significant risk factors were infected by one or more of at least five tick species in KP, Pakistan, and some of the ticks harbored the human pathogen R. massiliae.

Graphical abstract


Ticks cause substantial economic losses in various communities around the world, but especially in low-income livestock holders in tropical and subtropical regions, where approximately 80% of the world’s cattle population is at risk of infestation [1]. These hematophagous ectoparasites play a major role in the transmission of pathogens, including bacteria, such as intracellular coccobacilli of the genus Rickettsia, and several protozoans and viruses that cause disease and are a threat to human and veterinary health [2,3,4]. Only 10% of tick species have been identified as carriers of pathogens; the remaining species require further research to determine whether they are vectors of disease.

Both the persistence of various life stages of ticks and their distribution depend mainly on such factors as habitat, climatic conditions, host availability and anthropogenic activities. Both biotic (vegetation structure, host availability and management) and abiotic (rainfall regimen, temperature range) factors are most likely to influence tick distributions on a global scale [5,6,7]. The potential impacts of climate change and transboundary movements of humans and animals on the distribution of ticks have been widely studied [8, 9]. Understanding the niche of a tick species is a complex exercise due to the multiple variables operating on the different stages of the tick life-cycle, including regulation of the patterns of tick and host abundance and their spatial encounters [10]. Some tick species exhibit ecological flexibility and adapt easily to changing climate and novel habitats [10, 11]. Efforts have been made at global, regional and local scales to characterize the major environmental factors influencing tick distributions using statistical models and correlative approaches [12]. The documentation of risk factors associated with tick distribution contributes to the design of cost-effective control measures.

Livestock activities are an integral part of Pakistan’s economy and the backbone of rural income, as more than 70% of the population live in rural areas [13]. According to the Agricultural Census Organization, the total equid population in Pakistan was 4.8 million in 2006, which rose to 5.5 million as per the census report in 2012–2013. The Pakistan Economic Survey 2018–2019 reported that the horse population was 0.4 million, donkey population was 5.4 million and mule population was 0.2 million ( Khyber Pakhtunkhwa (KP) is a mountainous region of Pakistan that lacks sufficient access to transport infrastructure. Consequently, equids play important roles in providing transportation, but also in farming, pumping water and milling activities in elevated areas. Equids are also used in fairs and games, such as horse races and polo festivals. However, despite their economic importance, the welfare of equids has received far less attention than that given to other animals in Pakistan.

Equids are frequently employed in outdoor activities and are therefore exposed to tick infestations [14, 15]. Previous studies conducted in Pakistan were mostly confined to domestic animals but excluded equids, thereby limiting our understanding of the prevalence of and associated risk factors for tick infestations on these animals [13, 16, 17]. In addition, although Pakistan’s climatic conditions are suitable for the survival of tick species that parasitize domestic animals, investigations on the frequency and distribution of ticks infesting equids are lacking, as are data on Rickettsia associated to equid ticks that are pathogenic to humans. The present study was designed to investigate the diversity of ticks, assess the presence of Rickettsia and estimate the risk factors associated with ticks infesting equids in KP.


Study area

We selected nine districts in northern, southern and central KP (northwestern geographic state of Pakistan, previously known as the North-West Frontier Province) as study areas: Charsadda (71.669079°E, 34.209086°N), Mardan (72.079179°E, 34.161767°N), Nowshera (72.005571°E, 33.988264°N), Peshawar (72.378867°E, 34.093417°N), Buner (72.484021°E, 34.450940°N), Malakand (71.863614°E, 34.519290°N), Upper Dir (71.995056°E, 35.066089°N), Swabi (72.394640°E, 34.119706°N) and Bajaur (71.621760°E, 34.634959°N). Based on access, nine regions of each of these districts were selected for tick sampling. A global positioning system (GPS) was used to obtain the geographic coordinates data and loaded onto a Microsoft Excel (Microsoft Corp., Redmond, WA, USA) spreadsheet to develop a distribution map for the study area using ArcGIS 10.3.1 (Fig. 1).

Fig. 1
figure 1

Agro-ecological zones of the selected districts of the Khyber Pakhtunkhwa (KP) region of Pakistan showing the locations (blue circles) at which ticks were collected. The study districts were: Bajaur (BJR), Upper Dir (UDR), Malakand (MKD), Charsadda (CHD), Peshawar (PSR), Mardan (MDN), Nowshera (NHR), Buner (BNR), Swabi (SWB)

Ethical consent and information collection

Ethical consent was obtained from the advanced studies and research board of the Abdul Wali Khan University Mardan. Oral/written consent was obtained from the owner of equids for the collection of ticks. A questionnaire was provided to equid holders that assisted in collection of data such as: host-related information (gender and age), eco-geography and herd management.

Tick collection and morphological identification

Ticks were collected between November 2018 and October 2019 from horses, donkeys and mules at different sites (villages, towns) in tehsils (i. e. administrative divisions) within the nine districts (Fig. 1). Two visits per month were made to each area, and information related to climate was retrieved. In infested areas, ticks of all stages were removed from the base of the tail, ears, perianal area, sternum, scrotal and belly, using forceps without damaging the morphological features of the ticks. The collected ticks were placed in vials containing 100% ethanol and brought to the laboratory for identification.

Collected ticks were identified morphologically under a StereoZoom microscope (HT StereoZoom) and following taxonomic keys [18, 19]. The specimens were preserved in 100% ethanol for later molecular analysis.

DNA extraction and PCR

Tick specimens were rinsed in distilled water and 70% ethanol and dried prior to DNA extraction [17]. Specimens were cut into pieces with a sterilized scalpel inside Eppendorf tubes. Genomic DNA of 91 randomly selected ticks was extracted using the phenol–chloroform method [20, 21]. Sterile conditions were maintained throughout the process. The concentration of DNA in each sample was quantified using a spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). The extracted DNA was stored at − 20 °C.

Conventional PCR targeting a fragment (approx. 460 bp) of the tick mitochondrial 16S rRNA gene was carried out as previously described [22] to verify successful extractions (Table 1). A real-time PCR (StepOne; Applied Biosystems, Thermo Fisher Scientific, Waltham, MA USA) was performed to detect a short fragment (147 bp) of the Rickettsia citrate synthase-encoding gene (gltA) using primers CS-5 and CS-6 and an internal fluorogenic probe (6-carboxyfluorescein [6-FAM]) [2]. Positive samples were subjected to two conventional PCR assays, one to obtain a larger sequence (401 bp) of the rickettsial gltA gene, using primers CS-78 and CS-323 [2], and the other to amplify a 532-bp fragment of the rickettsial outer membrane protein A (ompA) gene [23]. Conventional PCR reactions consisted of a 25-µl mix (2 µl genomic DNA, 1 µl of each primer, 8.5 µl PCR water and 12.5 µl master mix; Thermo Fisher Scientific). PCR cycling parameters were as previously reported [2, 22, 24] (see Table 1). PCR products were electrophoresed in 1.5% agarose gels, and the results were captured with GelDoc (UVP BioDoc-It imaging system; Analytik Jena AG, Jena, Germany). Amplicons of the expected size were purified with ExoSAP-IT (Thermo Fisher Scientific) and sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Thermo Fisher Scientific) in an ABI 3500 Genetic Analyzer automated device (Applied Biosystems, Thermo Fisher Scientific).

Table 1 List of primers used for the amplification of target genes of ticks and Rickettsia spp.

Sequence and phylogenetic analyses

The obtained sequences were assembled and trimmed in SeqMan v 5.00 (DNASTAR, Madison, WI, USA) and subjected to BLASTn analyses ( to infer closest identities with the organisms available in GenBank [25]. Alignments with our sequences and GenBank-retrieved sequences were constructed with CLUSTAL W in BioEdit [26, 27] manually edited with GeneDoc [28]. A phylogenetic tree for the tick mitochondrial 16S rRNA gene was implemented in MEGA X [29] and constructed using the neighbor-joining method with 1000 bootstrap replicons [30]. For alignment of the ompA rickettsial gene, we constructed a phylogenetic tree in MrBayes using Bayesian statistics [31], with four independent Markov chain runs for 1,000,000 metropolis-coupled Markov chain Monte Carlo (MCMC) generations, sampling a tree every 100th generation. The first 25% of the trees represented burn-in, and the remaining trees were used to calculate the Bayesian posterior probability.

Statistical analysis

The recorded observations were assembled and arranged in the spreadsheets of Microsoft Excel V 2016. The chi-square test (χ2) was used for differences, and relative risk (RR) was analyzed using GraphPad Prism v. 5.00 (GraphPad Software Inc., San Diego, CA, USA), with the 95% confidence interval (CI). Significance was set at P < 0.05.


Demographic characteristics of the study population

A total of 404 animals were sampled, including 200 horses, 171 donkeys and 33 mules, of which 324 (80.2%) were found to be infested with ticks (mean infestation 110.3, standard deviation [SD] 80.79, 95% CI − 90.35–311). Of the equids inspected, horses were the most infested (161, 85.5%) followed by donkeys (145, 80%) and mules (18, 55%). The highest tick burdens were found in equids inspected in Peshawar (103, 19%), Mardan (101, 18%), Nowshera (100, 18%), Swabi (96, 12.5%) and Charsadda (65, 11%); the lowest burdens were observed in Upper Dir (11, 2%), Malakand (22, 4%), Bajaur (24, 4%) and Buner (27, 5%) (Fig. 1; Table 2). The geography of each collection point significantly influenced the pattern of tick infestation on the hosts (P < 0.002).

Table 2 Ticks collected from equids in selected districts of Khyber Pakhtunkhwa, Pakistan, between November 2018 and October 2019

The trend line of tick infestation increased from May (61, 11%) to August (77, 14%) of 2019, with the highest tick infestation found to be throughout the summer (P < 0.0354) and the lowest tick infestation during the winter between November (13, 2.3%) and January (6, 1%) (P = 0.1907) (Fig. 2). Linear regression analysis showed a significant association between tick infestation and high humidity (r2 = 0.9006, P < 0.0001) and temperature (r2 = 0.75827, P < 0.0002).

Fig. 2
figure 2

Temporal distribution of tick species collected during this study

Ticks and equid populations

A total of 550 collected ticks were morphologically identified into two genera and five species as follows: Rhipicephalus (539 specimens, 98% of all ticks) and Hyalomma (11, 2%). The most numerous tick species was Rhipicephalus microplus (341 specimens, 62%) followed by Rh. haemaphysaloides (126, 23%), Rh. turanicus (39, 7%), Rh. sanguineus (s.l.) (33, 6%) and Hyalomma anatolicum (11, 2%).

A total of 263 ticks (48% of all collected ticks) were collected from horses, of which 185 were Rh. microplus (70% of the ticks from horses), 44 were Rh. haemaphysaloides (17%), 18 were Rh. turanicus (0.06%), ten were Rh. sanguineus (s.l.) (0.03%) and six were H. anatolicum (0.03%). From donkeys, 206 ticks (37% of all collected ticks) were collected, of which 110 were Rh. microplus (53% of the ticks from donkeys), 62 were Rh. haemaphysaloides (30%), 13 were Rh. turanicus (06%), 18 were Rh. sanguineus (s.l.) (8%) and three were H. anatolicum (1%). Mules were infested by 81 tick specimens (15% of all collected ticks), of which 46 were Rh. microplus (57% of the ticks from mules), 20 were Rh. haemaphysaloides (24%), eight were Rh. turanicus (10%), five were Rh. sanguineus (s.l.) (6%) and two were H. anatolicum (2%) (Table 3).

Table 3 Tick stages and species collected across Khyber Pakhtunkhwa, Pakistan, from equids

Tick infestation on different body parts of the equid host

Different parts of the body of the equid hosts were examined for ticks. The highest infestations were observed on the axial region of horses (63, 24%), the sternum of donkeys (44, 21%) and the belly of mules (19, 23%). The lowest tick-infested body parts were the dewlap in horses (37, 14%), the axial region in donkeys (25, 12%) and the scrotum of mules (10, 12%). Statistical analysis of infestation of the body regions of the equid hosts revealed significant differences (χ2 = 18.66, 10; P < 0.045) (Table 4).

Table 4 Number of ticks collected from different body parts of the equid hosts

Assessment of risk factors associated with tick infestations

Several independent variables (Table 5) potentially associated with tick infestations on equids were analyzed and the risk factors calculated. Analysis of infestation according to sex revealed that there were more tick infestations on male equids than on female equids, with the difference being statistically significant in the tick acquisition. When age was considered, the highest tick infestation was found on adult hosts aged 4–6 years compared to hosts aged 1–3 years. Most of the equids examined were in the plain areas, but tick infestations were predominantly found on equids in hilly areas, and the association between a high risk of tick infestation and equids in hilly areas was significant. When animal housing was considered, the tick infestation rate was found to be higher on equids in muddy systems than those maintained in concrete and semi-concrete housing. The risk of tick infestation was also higher in those managed in mud housing than in those in open grazing systems. Equids raised in a herd were found to have a higher risk of tick infestation than those reared solely. Providing fresh or stored food to the equids had no impact on the tick infestations (Table 5).

Table 5 Risk factors associated with tick infestation on equid hosts

Molecular confirmation of taxonomic identification of ticks

Partial sequences of the tick mitochondrial 16S rRNA gene were successfully obtained from Rh. microplus (56 specimens), Rh. turanicus (15 specimens), Rh. sanguineus (s.l.) (2 specimens), Rh. haemaphysaloides (12 specimens) and H. anatolicum (5 specimens). BLAST of obtained 16S rRNA sequences for Rh. microplus, Rh. haemaphysaloides and H. anatolicum resulted in the screening of homologous sequences from Pakistan, India and China that showed 98–99% identity. The BLAST of Rh. turanicus showed closest identity (95–97%) with conspecific sequences reported from India, Afghanistan, Iran, South Africa and China, and the Rh. sanguineus (s.l.) sequences showed highest identities (96–98%) with conspecific sequences from Taiwan, USA and Singapore. The resultant consensus sequences after trimming were deposited in the GenBank for each species; Rh. microplus (accession MW071171), Rh. turanicus (MW074299), Rh. sanguineus (s.l.) (MW113239), Rh. haemaphysaloides (MW113238) and H. anatolicum (MW172215). The phylogenetic analysis of the collected five species revealed an evolutionary relatedness with homologous species reported in Pakistan’s neighboring countries. In a phylogenetic tree, the obtained sequences of prevalent tick species clustered with the same species reported in the bordering countries (Fig. 3).

Fig. 3
figure 3

Phylogenetic analysis based on 16S rRNA gene partial sequences of Rhipicephalus species and Hyalomma anatolicum from Pakistan using neighbor-joining method in Mega X. Accession numbers are followed by species and country names. The bootstrap values are shown at each node (1000). The bar represents 0.1 substitutions per site. Sequences obtained in the present study are denoted with black circles

Detection of rickettsial DNA in ticks

The DNA of 91 tick specimens were molecularly tested by real-time PCR, which detected rickettsial DNA in 23 specimens of ticks (Rh. microplus, Rh. turanicus, Rh. haemaphysaloides and H. anatolicum). Among these, ten speciments showed high Cq values (> 35) and were excluded from further analyses. The remaining 13 samples were found to be positive for Rickettsia through conventional PCR targeting the rickettsial gltA and ompA genes. These Rickettsia-positive ticks included Rh. microplus-infected ticks (10/91, 11%) from all sampled districts, and Rh. haemaphysaloides-infected ticks (3/91, 3%) from Bajaur, Mardan and Charsadda districts. No Rickettsia spp. were detected in Rh. turanicus, Rh. sanguineus (s.l.) and H. anatolicum through PCR. A consensus partial sequence (333 bp) of the highly conserved Rickettsia gltA gene, detected in Rh. microplus and Rh. haemaphysaloides, was 100% identical to several sequences of R. massiliae from GenBank (MF002497, KY069259, KT588058, MT309038, MT309037, JN043507, CP000683, U59719). From these same ticks, a consensus partial sequence (592 bp) of the highly polymorphic Rickettsia ompA gene was 100% identical to a sequence of R. massiliae from GenBank (CP003319). The obtained sequences of gltA and ompA were uploaded to GenBank with accessions numbers MW922872 and MW928497 respectively. The obtained sequence of the ompA gene for Rickettsia spp. clustered with sequences of R. massiliae in a phylogenetic tree (Fig. 4).

Fig. 4
figure 4

Bayesian phylogenetic tree inferred from partial sequences of the rickettsial outer membrane protein A (ompA) gene of Rickettsia spp. The tree is drawn to scale with the scale bar indicating nucleotide substitutions per site. Values of Bayesian posterior probabilities are indicated above or below each branch. The sequence of Rickettsia massiliae detected in ticks from Pakistan is highlighted in bold


Ecological changes and an increasingly evident increase in global temperature have widely impacted the distribution of ticks and tick-borne diseases, ultimately increasing their threat to public and veterinary health [9, 10]. Based on spatial distribution and molecular evidence, it has been suggested that global climate change has played a significant role in the spread of tick species [10]. Previously, we reported numerous tick species infesting diverse hosts, including equids (horse), in selected districts of KP, Pakistan [13, 17]. However, a knowledge gap remained on the identity and effect of various factors affecting ticks infesting equid hosts and their associated Rickettsia spp. To fill this gap, we designed this study to investigate equid hosts for tick infestation, associated risk factors and the molecular phylogeny of these ticks in selected districts of KP, Pakistan. The collected ticks from five selected districts were categorized into two genera which included five medically important tick species: Rh. microplus, Rh. haemaphysaloides, Rh. turanicus, Rh. sanguineus (s.l.) and H. anatolicum. These ticks were found to be infesting equid hosts (horses, mules, donkeys) throughout the year in extreme cool to extreme hot weather conditions. The bacterium Rickettsia massiliae was detected in all Rh. microplus ticks collected in all nine study districts and in Rh. haemaphysaloides collected in Bajaur, Mardan and Charsadda districts. Ecological conditions and biotopes drive the geographic distribution and the risks associated with the possible harmful effects of tick species. The favorable environmental conditions in the region extend tick adaptation and also result in infestation of novel hosts [13, 16, 17].

Risk factors, such as age, altitude, housing type, grazing, living management and food conditions, of equid hosts with respect to tick infestation were statistically significant, with the exception of food conditions. A strong correlation was observed between the intensity of tick infestation and temperature and relative humidity, with the highest infestation recorded in the summer (June–August) and lowest infestation recorded in the winter (December–January); these results agree with previous reports in the region [16, 17, 32]. The variable climatic conditions of the KP favor tick survival and reproduction, resulting in high infestations on various hosts [17, 33].

The younger age group (1–years) of equid hosts was found to have lower tick infestations, possibly correlated with a strong immune response, shorter exposure to questing ticks, shorter time in free grazing conditions and frequent grooming, all of which protect younger equids from tick infestation [34,35,36]. Higher tick infestations on adult hosts may be due to the larger surface area of the body, which enhances the chances of tick attachment, frequent use and ultimate resistance to acaricide and low immunity [32, 36,37,38]. The muddy-keeping system of hosts provided favorable conditions for the survival and reproduction of ticks, as compared to concrete and semi-concrete keeping systems. Host density in herds ultimately increases the chances of tick infestation, compared to animals kept alone, due to the availability of host-questing ticks detached in the herd area [39,40,41]. Pastures with long grasses in the hilly areas may provide favorable conditions for the questing ticks where they await and attach to open-grazing equids.

Ticks collected across the different study regions belonged to two genera and five species. The most prevalent tick was Rh. microplus, which has been reported to the dominant tick species in the region [13, 16, 17]. The highest numbers of ticks were collected in the months of June to August, while the lowest numbers were collected in the coldest months of November to January, which is in accordance with the results of a previous study [17]. It is important to mention that several Haemaphysalis and Hyalomma species other than H. anatolicum are also abundant tick species in Pakistan [13, 16, 17], but these latter species were not found to infest the equids inspected in this study throughout the study period, possibly due to the availability of other suitable and accessible hosts for blood feeding. Several markers have been utilized to delineate the accurate identity of a tick species at the molecular level. For example, 16S rRNA, COX1, 12S rRNA and ITS are well-known genetic markers that have been implied to separate closely related tick species [17, 42,43,44, 51]. Since the objective of this work was only to identify tick species—and not population structure—we used 16S rRNA to identify the collected tick species and to establish their evolutionary relationship [43,44,45,46]. The phylogenetic analysis of the collected ticks showed the closest relation with the same species reported from other countries, including Afghanistan, China, South Africa and Taiwan.

Tick-borne Rickettsia spp. have been collected from several Rhipicephalus ticks infesting different hosts in various regions of the world [13, 47,48,49,50,51]. In Pakistan, reports are not available on ticks infesting equids associated with Rickettsia infection. In this study, Rickettsia massiliae was detected in Rh. microplus and Rh. haemaphysaloides using gltA and ompA rickettsial markers. There have been several reports of spotted fever due to R. massiliae in Europe [49]. Noteworthy, this is the first report of R. massiliae associated with equid ticks in Pakistan. Our findings highlight that human illness due to this spotted fever group of pathogens could be underreported in Pakistan, where to date tick-borne rickettsioses have not been reported.


To our knowledge this is the first report of equid ticks from selected regions of KP. The risk factor analyses showed that temperature, humidity, age, keeping and grazing systems, sex and host species considerably affected the intensity of tick infestations. The phylogenetic tree of the collected ticks based on the tick mitochondrial 16S rRNA gene (2 genera and 5 species) revealed a close resemblance to tick species reported from Afghanistan, China, South Africa and Taiwan. The human pathogen R. massiliae, a tick-borne pathogenic spotted fever group Rickettsia species, was detected in Rh. microplus and Rh. haemaphysaloides. These results demonstrate that it is important to investigate the vector potential of ticks for other pathogens infesting equids in different regions of Pakistan to avoid the emergence of zoonotic infections. They will also provide valuable data for use in developing integrated control management strategies against ticks and tick-borne diseases in Pakistan.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article.


  1. Lew-Tabor AE, Rodriguez Valle MA. A review of reverse vaccinology approaches for the development of vaccines against ticks and tick borne diseases. Ticks Tick Borne Dis. 2016;7:573–85.

    Article  CAS  PubMed  Google Scholar 

  2. Labruna MB, Whitworth T, Bouyer DH, McBride J, Camargo LMA, Camargo EP, et al. Rickettsia bellii and Rickettsia amblyommii in Amblyomma ticks from the State of Rondonia. Western Amazon Brazil J Med Entomol. 2004;41:1073–81.

    PubMed  Google Scholar 

  3. Gondard M, Cabezas-Cruz A, Charles RA, Vayssier-Taussat M, Albina E, Moutailler S. Ticks and tick-borne pathogens of the Caribbean: current understanding and future directions for more comprehensive surveillance. Front Cell Infect Microbiol. 2017;7:1–16.

    Article  Google Scholar 

  4. Schorderet-Weber S, Noack S, Selzer PM, Kaminsky R. Blocking transmission of vector-borne diseases. Int J Parasitol-Drug. 2017;7:90–109.

    Article  Google Scholar 

  5. Estrada-Peña A, Ayllón N, De La Fuente J. Impact of climate trends on tick-borne pathogen transmission. Front Physiol. 2012;3:1–12.

    Article  Google Scholar 

  6. Szabó MPJ, Pinter A, Labruna MB. Ecology, biology and distribution of spotted-fever tick vectors in Brazil. Front Cell Infect Microbiol. 2013;3:1–9.

    Article  Google Scholar 

  7. Esser HJ, Mögling R, Cleton NB, Van Der Jeugd H, Sprong H, Stroo A, et al. Risk factors associated with sustained circulation of six zoonotic arboviruses: a systematic review for selection of surveillance sites in non-endemic areas. Parasit Vectors. 2019;12:1–17.

    Article  Google Scholar 

  8. Pietzsch ME, Medlock JM, Jones L, Avenell D, Abbott J, Harding P, et al. Distribution of Ixodes ricinus in the British Isles: investigation of historical records. Med Vet Entomol. 2005;19:306–14.

    Article  CAS  PubMed  Google Scholar 

  9. Leger E, Vourch G, Vial L, Chevillon C, McCoy KD. Changing distributions of ticks: causes and consequences. Exp Appl Acarol. 2013;59:219–44.

    Article  PubMed  Google Scholar 

  10. Estrada-Peña A. Climate, niche, ticks, and models: what they are and how we should interpret them. Parasitol Res. 2008;103:87–95.

    Article  Google Scholar 

  11. Egekwu N, Sonenshine DE, Bissinger BW, Roe RM. Transcriptome of the female synganglion of the black-legged tick Ixodes scapularis (Acari: Ixodidae) with comparison between Illumina and 454 systems. PLoS ONE. 2014;9(7):1–24.

    Article  CAS  Google Scholar 

  12. James AM, Burdett C, McCool MJ, Fox A, Riggs P. The geographic distribution and ecological preferences of the American dog tick, Dermacentor variabilis (Say), in the USA. Med Vet Entomol. 2015;29(2):178–88.

    Article  CAS  PubMed  Google Scholar 

  13. Karim S, Budachetri K, Mukherjee N, Williams J, Kausar A, Hassan MJ, et al. A study of ticks and tick-borne livestock pathogens in Pakistan. PLOS Negl Trop Dis. 2017;11(6):1–17.

    Article  Google Scholar 

  14. Laus F, Veronesi F, Passamonti F, Paggi E, Cerquetella M, Hyatt D, et al. Prevalence of tick-borne pathogens in horses from Italy. J Vet Med Sci. 2013;75(6):715–20.

    Article  PubMed  Google Scholar 

  15. Kamran K, Ali A, Villagra C, Siddiqui S, Alouffi AS, Iqbal A. A cross-sectional study of hard ticks (Acari: Ixodidae) on horse farms to assess the risk factors associated with tick-borne diseases. Zoonoses Public Health. 2021;68(3):247–62.

    Article  PubMed  Google Scholar 

  16. Rehman A, Nijhof AM, Sauter-Louism C, Schauer B, Staubach C, Conraths FJ. Distribution of ticks infesting ruminants and risk factors associated with high tick prevalence in livestock farms in the semi-arid and arid agro-ecological zones of Pakistan. Parasites Vectors. 2017;10(1):1–15.

    Article  Google Scholar 

  17. Ali A, Khan MA, Zahid H, Yaseen PM, Khan MQ, Nawab J, et al. Seasonal dynamics, record of ticks infesting humans, wild and domestic animals and molecular phylogeny of Rhipicephalus microplus in Khyber Pakhtunkhwa Pakistan. Front Physiol. 2019;10:1–15.

    Article  Google Scholar 

  18. Geevarghese G, Dhanda V. The Indian Hyalomma ticks (Ixodoidea: Ixodidae). New Delhi: ICAR Press; 1987. p. 119.

    Google Scholar 

  19. Walker JB, Keirans JE, Horak IG. The genus Rhipicephalus (Acari, Ixodidae): a guide to the brown ticks of the world. Cambridge: Cambridge University Press; 2005.

    Google Scholar 

  20. Sambrook J, Russell DW. Molecular cloning: a laboratory manual. New York: Cold Spring Harbor Laboratory Press; 2001. p. 23–44.

    Google Scholar 

  21. Sangioni LA, Horta MC, Vianna MC, Gennari SM, Soares RM, Galvão MA, et al. Rickettsial Infection in animals and Brazilian spotted fever endemicity. Emerg Infect Dis. 2005;11:265–70.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Mangold AJ, Bargues MD, Mas-Coma S. Mitochondrial 16S rDNA sequences and phylogenetic relationships of species of Rhipicephalus and other tick genera among Metastriata (Acari: Ixodidae). Parasitol Res. 1998;84:478–84.

    Article  CAS  PubMed  Google Scholar 

  23. Regnery RL, Spruill CL, Plikaytis BD. Genotypic identification of rickettsiae and estimation of intraspecies sequence divergence for portions of two rickettsial genes. J Bacteriol. 1991;173(5):1576–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Roux V, Fournier PE, Raoult D. Differentiation of spotted fever group rickettsiae by sequencing and analysis of restriction fragment length polymorphism of PCR-amplified DNA of the gene encoding the protein rOmpA. J Clin Microbiol. 1996;34:2058–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.

    Article  CAS  PubMed  Google Scholar 

  26. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22(22):4673–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Hall T, Biosciences I, Carlsbad C. BioEdit: an important software for molecular biology. GERF Bull Biosci. 2011;2:60–1.

    Google Scholar 

  28. Nicholas KB. GeneDoc: analysis and visualization of genetic variation. EMBNEW News. 1997;4:14.

    Google Scholar 

  29. Kumar S, Stecher G, Li M, Knyaz C, Tamura KMEGAX. molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5.

    Article  CAS  PubMed  Google Scholar 

  32. Kamran K, Ali A, Villagra AC, Bazai ZA, Iqbal A, Sajid MS. Hyalomma anatolicum resistance against ivermectin and fipronil is associated with indiscriminate use of acaricides in southwestern Balochistan. Pakistan Parasitol Res. 2020;120:1–25.

    Google Scholar 

  33. Zahid H, Munoz-Leal S, Khan MK, Alouffi A, Labruna MB, Ali A. Life cycle and genetic identification of Argas persicus infesting domestic fowl in Khyber Pakhtunkhwa. Pakistan Front Vet Sci. 2021;8:1–10.

    Google Scholar 

  34. Swai ES, Mbise AN, Kessy V, Kaaya E, Sanka P, Loomu PM. Farm constraints, cattle disease perception and tick management practices in pastoral Maasai community-Ngorongoro. Tanzania Livest Res Rural Dev. 2005;17(2):17–20.

    Google Scholar 

  35. Mooring B, Harte H. Testing the interspecific body size principle in ungulates: the smaller they come, the harder they groom. Anim Behav. 2000;60:35–45.

    Article  CAS  PubMed  Google Scholar 

  36. Lorusso V, Picozzi K, de Bronsvoort BM, Majekodunmi A, Dongkum C, Balak G, et al. Ixodid ticks of traditionally managed cattle in central Nigeria: where Rhipicephalus (Boophilus) microplus does not dare (yet?). Parasites Vectors. 2013;6:1–10.

    Article  Google Scholar 

  37. Kemal J, Tamerat N, Tuluka T. Infestation and identification of ixodid tick in cattle: the case of Arbegona District Southern Ethiopia. J Vet Sci. 2016;5:1–8.

    Google Scholar 

  38. Tabor AE, Ali A, Rehman G, Rocha Garcia G, Zangirolamo AF, Malardo T, et al. Cattle tick Rhipicephalus microplus-host interface: a review of resistant and susceptible host responses. Front Cell Infect Microbiol. 2017;7:1–18.

    Article  CAS  Google Scholar 

  39. Labruna MB, Kasai N, Ferreira F, Faccini JL, Gennari SM. Seasonal dynamics of ticks (Acari: Ixodidae) on horses in the state of São Paulo. Brazil Vet Parasitol. 2002;105(1):65–77.

    Article  PubMed  Google Scholar 

  40. Greenfield BPJ. Environmental parameters affecting tick (Ixodes ricinus) distribution during the summer season in Richmond Park. London Biosci Horiz. 2011;4(2):140–8.

    Article  Google Scholar 

  41. Gharbi M, Moussi N, Jedidi M, Mhadhbi M, Sassi L, Darghouth MA. Population dynamics of ticks infesting the one-humped camel (Camelus dromedarius) in central Tunisia. Ticks Tick Borne Dis. 2013;4(6):488–91.

    Article  PubMed  Google Scholar 

  42. Burger TD, Shao R, Barker SC. Phylogenetic analysis of mitochondrial genome sequences indicates that the cattle tick, Rhipicephalus (Boophilus) microplus, contains a cryptic species. Mol Phylogenet Evol. 2014;76:241–53.

    Article  PubMed  Google Scholar 

  43. Lv J, Wu S, Zhang Y, Chen Y, Feng C, Yuan X, et al. Assessment of four DNA fragments (COI, 16S rDNA, ITS2, 12S rDNA) for species identification of the Ixodida (Acari: Ixodida). Parasites Vectors. 2014;7(1):1–11.

    Article  CAS  Google Scholar 

  44. Low VL, Tay ST, Kho KL, Koh FX, Tan TK, Lim YAL, et al. Molecular characterization of the tick Rhipicephalus microplus in Malaysia: new insights into the cryptic diversity and distinct genetic assemblages throughout the world. Parasit Vectors. 2015;8(1):1–10.

    Article  CAS  Google Scholar 

  45. Ali A, Parizi LF, Ferreira BR, Junior V, Silva I. A revision of two distinct species of Rhipicephalus: R. microplus and R. australis. Cienc Rural. 2016;46:1240–8.

    Article  CAS  Google Scholar 

  46. Roy BC, Estrada-Peña A, Krücken J, Rehman A, Nijhof AM. Morphological and phylogenetic analyses of Rhipicephalus microplus ticks from Bangladesh, Pakistan and Myanmar. Ticks Tick Borne Dis. 2018;9(5):1069–79.

    Article  PubMed  Google Scholar 

  47. Eremeeva ME, Zambrano ML, Anaya L, Beati L, Karpathy SE, Santos-Silva MM, et al. Rickettsia rickettsii in ticks, Mexicali, Mexico. J Med Entomol. 2011;48(2):418–21.

    Article  PubMed  Google Scholar 

  48. Abarca K, López J, Acosta-Jamett G, Martínez-Valdebenito C. Rickettsia felis in Rhipicephalus sanguineus from two distant Chilean cities. Vector Borne Zoonotic Dis. 2013;13(8):607–9.

    Article  PubMed  Google Scholar 

  49. Parola P, Paddock CD, Socolovschi C, Labruna MB, Mediannikov O, Kernif T, et al. Update on tick-borne rickettsioses around the world: a geographic approach. Clin Microbiol Rev. 2013;26(4):657–702.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Wei QQ, Guo LP, Wang AD, Mu LM, Zhang K, Chen CF, et al. The first detection of Rickettsia aeschlimannii and Rickettsia massiliae in Rhipicephalus turanicus ticks, in northwest China. Parasites Vectors. 2015;8(1):1–4.

    Article  CAS  Google Scholar 

  51. Ortuño A, Sanfeliu I, Nogueras MM, Pons I, López-Claessens S, Castellà J, et al. Detection of Rickettsia massiliae/Bar29 and Rickettsia conorii in red foxes (Vulpes vulpes) and their Rhipicephalus sanguineus complex ticks. Ticks Tick Borne Dis. 2018;9(3):629–31.

    Article  PubMed  Google Scholar 

Download references


We acknowledge the partial financial support provided by the Pakistan Science Foundation and the Higher Education Commission of Pakistan.


Not applicable.

Author information

Authors and Affiliations



AA, SML, ASA and MBL designed the study and acquired financing. AA, HZ, IZ, MT, SK, MH, MT, BA and MH collected the samples. IZ, HZ and ASA carried out the statistical analysis. AA, HZ, IZ, MT, SK, MH, MT, BA and MH, SML and MBL conducted the experiments. AA, HZ, IZ, ASA, SML and MBL contributed to the revision of critically important intellectual contents and manuscript writing. All authors carried out critical revisions. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Abid Ali.

Ethics declarations

Ethics approval and consent to participate

Ethical consent was obtained from the advanced studies and research board of the Abdul Wali Khan University Mardan. Oral/written consent was obtained from the owner of equids for the collection of ticks and data such as host-related information (gender and age), eco-geography and herd management.

Consent for publication

Not applicable.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ali, A., Zahid, H., Zeb, I. et al. Risk factors associated with tick infestations on equids in Khyber Pakhtunkhwa, Pakistan, with notes on Rickettsia massiliae detection. Parasites Vectors 14, 363 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: