- Open Access
Low genetic diversity of Ehrlichia canis associated with high co-infection rates in Rhipicephalus sanguineus (s.l.)
- Alejandro Cabezas-Cruz†1,
- Eleonore Allain†1,
- Abdullah S. Ahmad†2, 3,
- Muhammad A. Saeed2,
- Imran Rashid3,
- Kamran Ashraf3,
- Lena Yousfi1,
- Wasim Shehzad4,
- Lea Indjein2,
- Manuel Rodriguez-Valle5,
- Agustin Estrada-Peña6,
- Dasiel Obregón7,
- Abdul Jabbar2Email author and
- Sara Moutailler1Email author
© The Author(s). 2019
- Received: 5 October 2018
- Accepted: 9 November 2018
- Published: 7 January 2019
Rhipicephalus sanguineus sensu lato (s.l.) is the most widely distributed ixodid tick and is a vector of major canine and human pathogens. High-throughput technologies have revealed that individual ticks carry a high diversity of pathogens, including bacteria, protozoa and viruses. Currently, it is accepted that co-infections (multiple pathogen species within an individual) are very common in ticks and influence pathogen acquisition and transmission as well as host infection risk. However, little is known on the impact of the genetic diversity of pathogens on the incidence of co-infections. Herein, we studied the frequency of co-infections in R. sanguineus (s.l.) and their association with the genetic diversity of Ehrlichia canis.
Rhipicephalus sanguineus (s.l.) female ticks (n = 235) were collected from healthy farm dogs in three districts of Pakistan. Microfluidic real-time PCR, a powerful nanotechnology for high-throughput molecular detection of pathogens, was used to test the presence of 25 bacterial and seven parasitic species in individual ticks. The genetic diversity of E. canis was evaluated by characterizing the trp36 gene.
A total of 204 ticks were infected with at least one pathogen and 109 co-infected with two (80%) or three (20%) pathogens. Rickettsia massiliae (human pathogen) and E. canis (zoonotic dog pathogen) were the most common pathogens co-infecting (30.4%) ticks. Furthermore, all identified co-infections included R. massiliae and/or E. canis. Multiple correspondence analysis (MCA) revealed that single infections did not show clear regional association whereas some co-infections were restricted to certain geographical regions. The sequence analysis of trp36 in representative samples allowed the identification of three E. canis strains with low genetic diversity, and the strain found in Muzaffargarh district appeared to be more adapted to co-infection with R. massiliae.
Rhipicephalus sanguineus (s.l.) harbors multiple co-infections with human and dog pathogens of zoonotic potential. Findings of this study suggest that genetic diversity of E. canis may favor co-infections with different pathogens.
- Rhipicephalus sanguineus (s.l.)
- Ehrlichia canis
- Rickettsia massiliae
- Genetic diversity
Ticks (Acari: Ixodidae) are important ectoparasites of animals and humans that cause mechanical damage to their hosts and serve as vectors for the transmission of various pathogens of medical and veterinary importance [1–3]. Although tick-borne pathogens are maintained in stable natural cycles involving ticks and animals (domestic and/or wild), humans may serve as accidental hosts . In the past few years, the distribution and abundance of many tick species have risen in canids worldwide, resulting in an increased prevalence of tick-borne pathogens in these animals [2, 5–7]. Given that domestic animals (such as dogs) live in close proximity to humans and they may act as reservoir hosts for pathogens that can infect humans, the distribution and dissemination of ticks in these animals could be a significant public health concern [8–10].
Rhipicephalus sanguineus (sensu lato), the brown dog tick, is the most widely distributed ixodid tick that infests human as well as canine dwellings . Being a three-host tick, three different developmental stages (larvae, nymphs and adults) of R. sanguineus (s.l.) usually feed on dogs [1, 11]. Rhipicephalus sanguineus (s.l.) is known to transmit a number of pathogens that produce disease in dogs, including babesiosis (caused by Babesia canis, B. gibsoni and B. vogeli), ehrlichiosis (Ehrlichia canis and Anaplasma platys), hepatozoonosis (Hepatozoon canis), and rickettsioses/spotted fever (Rickettsia rickettsii, R. massiliae and R. conorii), especially in tropical and subtropical regions of the world [1, 11–15]. There are also reports of R. typhi, a rickettsia of the Typhus group mainly transmitted by fleas, detected in R. sanguineus (s.l.) . Other important canine ticks (and pathogens they transmit) include Amblyomma spp. (Hepatozoon spp. and E. canis), Dermacentor spp. (Babesia spp., Ehrlichia spp. and Rickettsia spp.), Haemaphysalis spp. (Babesia spp. and H. canis), Hyalomma spp. (Theileria annulata), Ixodes spp. (Borrelia spp., Ehrlichia spp. and Anaplasma spp.) and Rhipicephalus microplus (Anaplasma spp.) in various parts of the world [3, 16–18]. Importantly, many of these tick-borne pathogens in dogs could be of zoonotic importance, which warrants for studies assessing tick-borne pathogens in ticks collected from dogs. The risk of zoonotic tick-borne diseases in dogs is especially important in geographical regions which offer conducive environments for the increasing of tick abundance [2, 6, 9].
Pakistan is an important agricultural country, located in South Asia (30°0' N, 70°0'E) and exhibits 10 different agro-ecological zones [6, 19, 20]. The vast majority of the country consists of sub-tropical and partially temperate regions which extend from Himalayas in the north to the Arabian Sea in the south . Due to such favourable climatic conditions, a high prevalence of ticks and tick-borne diseases could be expected in these regions. Although limited data are available on canine ticks from Pakistan, R. sanguineus (s.l.) appears to be the most widespread canine tick species in this country, with prevalence rates as high as 98% in dogs [21–23]. Other ticks, including Dermacentor spp., Haemaphysalis spp. and Hyalomma spp. have also been reported in dogs from Pakistan [21, 23]. Likewise, a high prevalence of tick-borne pathogens such as Anaplasma spp. and Babesia spp. have been observed in conventional diagnostic methods (e.g. blood smear examination) from dogs in Pakistan [21, 24–26]. Recent studies from Pakistan utilized molecular (18S rRNA gene amplification) approaches and reported a high prevalence (46% and 12%) of H. canis from farm and pet dogs, respectively [6, 27]. Overall, these studies indicate that multiple tick species and a variety of tick-borne pathogens are prevalent in dogs from Pakistan. However, there is paucity of information on the occurrence and magnitude of tick-borne pathogens and their co-infections in ticks collected from dogs in Pakistan. Co-infections of tick-borne pathogens have been reported to be common in several tick species in different geographical areas [28–30], including R. sanguineus (s.l.) [31–33], and dogs were found infected with more than one tick-borne pathogen [34–37]. This highlights the importance of considering co-infections in tick-borne pathogen surveys.
Given the spectrum of pathogens transmitted by ticks, some organisms may not be identified using artificial culture methods; hence, molecular approaches are indispensable for thorough detection of multiple tick-borne pathogens [2, 3, 38]. However, conventional molecular approaches (e.g. PCR) could amplify and detect only known target pathogens. Furthermore, only a restricted number of target pathogens could be amplified and tested due to limiting factors (e.g. quantity of DNA) in PCR. To address these issues, a high throughput epidemiological surveillance method (based on a microfluidic system) was developed [38, 39]. The system holds the capacity to perform parallel real-time PCR using a small volume/quantity of DNA and can process up to 9216 individual reactions simultaneously [38, 39]. This high throughput system has been used successfully for the detection of four tick-species and 37 tick-borne pathogens in Europe .
Very little is known about the tick-borne pathogens in canine ticks from Pakistan. Therefore, this study aimed to use the microfluidic high-throughput system to determine tick-borne pathogens of zoonotic importance in ticks collected from dogs in three different agro-ecological zones of Punjab, Pakistan.
Study area and tick samples
DNA extraction and PCR pre-amplification
Prior to DNA extraction, ethanol was removed, and individual ticks were washed three times (30 min) in distilled H2O, and then ground using a plastic mortar. DNA was extracted using a DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol. DNA was stored at -20 °C until further use.
For better detection of pathogen DNAs, total DNAs were pre-amplified with the Perfecta Preamp Supermix (Quanta Biosciences, Beverly, USA) according to the manufacturer’s instructions. Primers (targeted all pathogens) were pooled combining equal volume of primers (200 nM final each). The reaction was performed in a final volume of 5 μl containing 1 μl Perfecta Preamp 5X, 1.25 μl pooled primers mix, 1.5 μl distilled water and 1.25 μl DNA, with one cycle at 95 °C for 2 min, 14 cycles at 95 °C for 10 s and 60 °C for 3 min. At the end of the cycling program the reactions were diluted as 1:10. Pre-amplified DNAs were stored at -20 °C until further use.
Microfluidic real-time PCR
To detect major tick-borne pathogens (see below), the BioMarkTM real-time PCR system (Fluidigm, California, USA) was used for high-throughput microfluidic real-time PCR amplification using the 48.48 dynamic arrays (Fluidigm). These chips dispensed 48 PCR mixes and 48 samples into individual wells, after which on-chip microfluidics assembled real-time PCR reactions in individual chambers prior to thermal cycling resulting in 2304 individual reactions. Targeted pathogens (and markers) included: Borrelia spp. (23S), B. burgdorferi (rpoB), B. garinii (rpoB), B. afzelii (fla), B. valaisiana (ospA), B. lusitaniae (rpoB), B. spielmanii (fla), B. bissettii (rpoB), B. bissettii (glpQ), Anaplasma spp. (16S), A. marginale (msp1), A. platys (groEL), A. phagoctyophilum (msp2), A. ovis (msp4), A. centrale (groEL), A. bovis (groEL), Ehrlichia spp. (16S), Ehrlichia canis (gltA), Neorickettsia mikurensis (groEL), Rickettsia spp. (gltA), R. conorii (ITS), R. slovaca (ITS), R. massiliae (ITS), R. helvetica (ITS), R. helvetica (ITS), R. felis (orfB), Bartonella spp. (ssrA), B. henselae (pap31), Francisella tularensis (tul4 and fopA), Coxiella endosymbiont (IS1111 and icd), Apicomplexa (18S), Babesia microti (CCTeta), B. canis (hsp70), B. ovis (18S), B. bovis (CCTeta), B. caballi (rap1), Babesia str. EU1 (18S), B. divergens (hsp70), Mycoplasma spp. (16S), Theileria spp. (18S) and Hepatozoon spp. (18S).
Briefly, amplifications were performed using 6-carboxyfluorescein (FAM)- and black hole quencher (BHQ1)-labeled TaqMan probes with TaqMan Gene expression master mix in accordance with manufacturer’s instructions (Applied Biosystems, Massachusetts, USA) . PCR cycling comprised of a denaturation step at 95 °C for 5 min followed by 45 cycles at 95 °C for 10 s, 60 °C for 15 s and 40 °C for 10 s. One negative water control was included per chip. To assess inhibitory molecules (that could inhibit PCR) in samples, DNA from Escherichia coli (EDL933 strain) was added to each sample as an internal inhibition control, and primers and probe specific for the eae gene of E. coli were used.
Validation of microfluidic real-time PCR results
Primers used for microfluidic real-time PCR
Forward primer (5'-3')
Reverse primer (5'-3')
Babesia spp., Theileria spp., Hepatozoon spp.
Anaplasma spp., Ehrlichia spp.
Analysis of Ehrlichia canis trp36 gene and encoded amino acid sequence
To assess the genetic diversity of E. canis, the full-length of the tandem repeat protein 36 (trp36) gene was amplified and sequenced. This gene codes for a major immunogenic protein and commonly used as a genetic marker to characterize the genetic diversity of this bacterium [41–43]. The amplification of trp36 was performed using primers and following a protocol previously published . The tandem repeats finder (TRF) database  was used to predict the presence of tandem repeats in the trp36. For sequence analysis, the predicted amino acid sequence of TRP36 was divided into three regions as previously reported [42, 43]. Region I, the 5' end pre-repeat region composed of 426–429 bp/142–143 aa at the N terminus of the encoded protein; Region II, the tandem repeat region (variable number of the 27 bp/9 aa repeat units depending on the strain); and Region III the 3' end post-repeat region (48, 84 and 90 bp/16, 28 and 30 aa) at the C terminus of the encoded protein.
Phylogenetic analysis using trp36
Nucleotide sequences of the pre-repeat region of E. canis trp36 were used for phylogenetic tree reconstruction as previously reported [41, 43]. Two sequences of the trp36 ortholog in E. chaffeensis (trp47) were used as the outgroup. Sequences available in the GenBank database were collected and aligned on MAFFT software , configured for the highest accuracy. The best-fit model of sequence evolution was selected based on Corrected Akaike Information Criterion (cAIC) and Bayesian Information Criterion (BIC) as implemented in MEGA v.6.00 . Tamura’s 3-parameter (T92)  model of amino acid evolution, which had the lowest values of cAIC and BIC, was chosen for tree reconstruction. The phylogenetic tree was reconstructed using the Maximum Likelihood (ML) method in MEGA. Initial tree for the heuristic search was obtained automatically by applying Nearest-Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using Maximum Composite Likelihood (MCL), and then selecting the topology with superior log-likelihood value. The analysis involved 30 nucleotide sequences. All positions containing gaps and missing data were eliminated. Reliability of internal branches was assessed using the bootstrapping method (1000 replicates). Graphical representation and editing of the phylogenetic trees were performed in MEGA.
Chi-square test was used to evaluate the co-occurrence of bacteria and to compare prevalence among different regions. Student’s t-test was used to assess significant differences of multiple co-infections in ticks according to sites of collection. Data were analyzed using GraphPad 5 Prism program (GraphPad Software Inc., La Jolla, CA, USA). Multiple correspondence analysis (MCA) was used to analyze the pattern of relationships between single and co-infections and tick origin. To this aim, a standard ‘Burt matrix’ analysis was performed. All possible pairwise pathogen combinations were calculated as (a)(a-1)/2 and all possible triple pathogen combinations were calculated as (a)(a-2)/3, where ‘a’ is the number of individual pathogen(s).
Prevalence of pathogens and Coxiella endosymbiont
Prevalence of pathogens in ticks collected from farm dogs in three districts of Pakistan
R. massiliae + E. canis
R. massiliae + Coxiella-like
R. massiliae + A. marginale
R. massiliae + A. phagocytophilum
E. canis + B. canis
R. massiliae + A. centrale
R. massiliae + B. canis
R. massiliae + A. platys
E. canis + A. phagocytophilum
E. canis + A. centrale
R. massiliae + A. ovis
R. massiliae + E. canis + Coxiella-like
R. massiliae + E. canis + A. marginale
R. massiliae + E. canis + B. canis
R. massiliae + E. canis + A. phagocytophilum
R. massiliae + E. canis + . platys
R. massiliae + E. canis + A. ovis
R. massiliae + E. canis + A. centrale
E. canis + A. centrale + Coxiella-like
Co-infections of ticks
Among the 204 infected ticks, 108 (53%) were found to be co-infected with two (80%) or three (20%) pathogens. Some microorganisms (e.g. Coxiella-like, B. canis, A. ovis and A. centrale) were found only in mixed infections (Table 2). The most frequent co-infection was due to R. massiliae and E. canis with a prevalence of 26.4% (P < 0.005). All other co-infections had prevalence values lower than 5% which were not statistically significant.
Of 21 possible triple pathogen combinations included in the assay, only eight were detected in this study although all without any significant associations (Table 2). All triple co-infections included both R. massiliae and E. canis except for one co-infection (E. canis, A. centrale and Coxiella-like) with a very low prevalence (Table 2). A significant association (P < 0.05) was observed between E. canis and R. massiliae in both Muzaffargarh and Rawalpindi districts. The prevalence of R. massiliae + E. canis co-infection was significantly different between Muzaffargarh and Kasur (χ2 = 38.1, P < 0.0001) and between Rawalpindi and Muzaffargarh (χ2 = 17.01, P < 0.0003).
Genetic variability of Ehrlichia canis
The full coding region of the gene trp36 was amplified and sequenced in four E. canis-positive samples per region (n = 12). Following theoretical translation of all trp36 nucleotide sequences into proteins, three different sequences were identified and named Kasur (GenBank: MH608289), Muzaffargarh (GenBank: MH608290) and Rawalpindi (GenBank: MH549197). To test the genetic diversity of E. canis in ticks from Pakistan, we performed three types of analyses based on trp36 nucleotide and protein sequences as described previously , including (i) nucleotide and amino acid sequences identity; (ii) molecular signatures of the pre-repeat region of TRP36 (i.e. presence/absence of glycine (G) 117 and sequence of putative glycosylation sequon at asparagine (N) 125; (iii) tandem repeat composition of TRP36; (iv) length of post-repeat region; and (v) phylogenetic analysis using pre-repeat region of the gene trp36.
Percentage of identity of trp36 nucleotide (above the diagonal) and putative amino acid sequences (below the diagonal) between Ehrlichia canis strains
E. canis strain
Comparison of Ehrlichia canis strains prevalence in single infections and co-infections with Rickettsia massiliae
E. canis strain
Single infection n (%)
Co-infection n (%)
Kasur vs Muzaffargarh
χ2 = 45.2, P < 0.0001
Kasur vs Rawalpindi
χ2 = 5.1, P = 0.07
Muzaffargarh vs Rawalpindi
χ2 = 73.3, P < 0.0001
To our knowledge, this is the first comprehensive study on genetic characterization of pathogens in R. sanguineus (s.l.) ticks collected from farm dogs in three different agro-climatic regions of Punjab, Pakistan. Rhipicephalus sanguineus (s.l.) is known to transmit a number of pathogens that produce diseases in dogs and humans, especially in tropical and subtropical regions of the world [1, 11–15]. However, there is paucity of information on the prevalence of pathogens in R. sanguineus (s.l.) from dogs in Asia. Although, standard PCR is most frequently used tool for the detection of tick-borne pathogens, it could be biased as pathogen detection is strongly influenced by particular research interests . Microfluidic real-time PCR (a high-throughput state-of-the-art technology) presents an alternative solution as it has the ability to detect a diverse array of tick-borne pathogens [28, 38]. Therefore, we used microfluidic real-time PCR approach to test the presence of 25 bacterial and seven parasite species in individual ticks. Furthermore, the genetic diversity of the most commonly detected pathogen (E. canis) was assessed using the trp36 gene sequencing.
Rhipicephalus sanguineus (s.l.) harbored multiple co-infections with human and dog pathogens of zoonotic potential. Importantly, all the ticks used for this study were partially or fully engorged and therefore the detection of pathogen DNA either represents the most recent blood meal and/or a prior infection of the tick. Once acquired by a tick, tick-borne pathogens are transmitted transstadially and thus, each blood meal increases the chance of acquiring new pathogens which in turn increases the likelihood of co-infections within the tick. This is supported by several field surveys in which co-infections have been reported in several tick species, including Ixodes ricinus [28, 49], I. scapularis [29, 50], Dermacentor reticulatus , Haemaphysalis longicornis  and R. sanguineus (s.l.) [31–33].
An important finding in this study is the detection of a wide array of pathogen combinations in R. sanguineus (s.l.) ticks from Pakistan, potentially due to the use of high-throughput approach (Table 2). For instance, only eight different combinations of tick-borne pathogens in I. ricinus have been reported in previous studies using standard serological and molecular (PCR) approaches [49, 51–55]. However, 31 different tick-borne pathogen combinations were detected in the same tick species using microfluidic real-time PCR, indicating that this technique could produce superior results than the standard PCR . We identified 19 different pathogen combinations occurring in R. sanguineus (s.l.) from Pakistan (Table 2). However, only nine pathogen combinations have been reported in R. sanguineus (s.l.) in previous studies, including Ehrlichia spp., Anaplasma spp., Hepatozoon spp., Babesia spp., Leishmania spp. and Cercopithifilaria spp. [31–33]. Importantly, our assay included not only pathogens commonly identified in R. sanguineus (s.l.) (e.g. E. canis)  but also pathogens rarely found in this tick species (e.g. Borrelia spp.) .
Ehrlichia canis was found in vast majority (12/19; 63%) of co-infections in R. sanguineus (s.l.), which is in agreement with previous studies [31, 33]. In addition, our results also showed a ubiquitous presence (15/19; 79%) of R. massiliae in the identified co-infections. Whether E. canis and/or R. massiliae facilitate the acquisition and transmission of other pathogens in R. sanguineus (s.l.) remains elusive, although at least one of these pathogens was found in all co-infections. Further studies are required in this context as co-infections could considerably impact the epidemiology of tick-borne pathogens. For instance, recent studies have indicated that Babesia microti (the causative agent of human babesiosis) is emerging in areas endemic for Borrelia burgdorferi (the causative agent of Lyme disease) in the USA . The emergence of Babesia microti has become difficult to explain because this pathogen has a low ecological fitness characterized by poor transmission from Peromyscus leucopus mice to larval ticks and poor transstadial transmission from larvae to nymphs . The current hypothesis, supported by empirical data, is that B. burgdorferi increases B. microti transmission from Peromyscus leucopus mice, which act as reservoirs for both B. microti and B. burgdorferi (s.l.) to ticks . Interestingly, a recent study using Pseudomonas aeruginosa and five different phage virus parasites as a model showed that co-infections can accelerate host adaptation and diversification . Considering the findings in , it can be hypothesized that infection with multiple bacteria species may accelerate ecological innovation in ticks with a potential impact in tick fitness and pathogen transmission .
Since, the genetic diversity of tick-borne bacterial pathogen impacts pathogenicity, virulence, host specificity, prevalence and transmission [48, 60, 61], we hypothesize that it could also influence co-infections. The effect of strain diversity on tick-borne pathogen co-infection remains poorly studied, although co-infection with multiple strains of A. marginale is well-documented in ticks and hosts [60, 61]. Our study did not specifically evaluate the effect of genetic diversity of E. canis on co-infection within individual ticks. However, all the samples from the same region contained an identical E. canis trp36 sequence which is suggestive of local adaptation of E. canis strains. Remarkably, the strain Muzaffargarh appeared to be more adapted or permissive to co-infection with R. massiliae. Some strains of Flavobacterium columnare, an environmental opportunistic bacterium, were found to be more permissible to co-infections in Danio rerio . In co-infection systems, interacting pathogens can compete, cooperate or coexist . The high rate of co-infection between E. canis strain Muzaffargarh and R. massiliae and the low rate of single infection of E. canis strain Muzaffargarh suggest that this strain has very low fitness in single infections and/or that it cooperates effectively with R. massiliae to infect ticks.
Previous studies indicated that trp36 sequences of E. canis strains were the most variable among other immunodominant protein-encoding genes (e.g. gp200, gp140, gp19) sequenced . Thus, small nucleotide variation in trp36 indicates low genetic diversity in E. canis. We found low genetic diversity of E. canis in the three regions of Pakistan. Similar results were reported in Taiwan where all identified strains were highly similar between them . Low genetic diversity of E. canis appears to be associated with high prevalence of this pathogen in ticks which contrasts with results obtained with A. marginale where high prevalence was associated with high genetic diversity of these bacteria in cattle. Despite low genetic diversity, trp36 sequences differences allowed the identification of three distinct E. canis strains (one from each region) which grouped under the same clade. Intrinsic transmission efficiency of specific E. canis strains may explain strain predominance in regions of high prevalence. Collectively, these data suggest a minor strain variation and vast geographical spread of the bacterial parasite in Pakistan. It is possible that these genotypes might have diverged within the country and have not been introduced as such from other countries. However, intercontinental movement of dogs could also be a contributory factor in the spread of R. sanguineus (s.l.) and associated pathogens [23, 65].
This study reports that tick-borne pathogen co-infections are very common in R. sanguineus (s.l.) ticks from Pakistan. The high prevalence of co-infection with E. canis and R. massiliae in R. sanguineus (s.l.) suggests a synergism between the two bacterial pathogens which in turn potentially increases the likelihood of acquiring a third pathogen. Low genetic diversity of E. canis was associated with high prevalence of this bacterium in R. sanguineus (s.l.) of Pakistan. The strain E. canis Muzaffargarh seems to be more adapted than others to co-infection with R. massiliae in R. sanguineus (s.l.)
We are grateful to farmers for allowing us the collection of ticks from their dogs for this study. We would like to thank Dr Muhammad Luqman Sohail (The Islamia University, Bahawalpur, Pakistan) for his assistance in sample collection.
ASA is grateful to the Higher Education Commission of Pakistan for the financial assistance under the International Research Support Initiative Programme (IRSIP).
Availability of data and materials
All data generated or analysed during this study are included in this published article. Nucleotide sequences have been submitted to the GenBank database under the accession nos. MH608289, MH608290 and MH549197.
ACC, AJ and MRV conceived the idea. AS performed field work and extracted DNA from ticks. EA and LI performed molecular detection of pathogens and gene amplification. EA, ACC, SM, AEP and DO performed data analysis. ACC, AS, MAS and AJ drafted the manuscript with critical input from IR, KA, WS, LI and MRV. All authors read and approved the final manuscript.
Ethics approval and consent to participate
This study was approved by the Animal Ethics Committee of the University of Veterinary and Animal Sciences, Lahore, Pakistan.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Dantas-Torres F. Biology and ecology of the brown dog tick, Rhipicephalus sanguineus. Parasit Vectors. 2010;3:26.View ArticleGoogle Scholar
- Maia C, Ferreira A, Nunes M, Vieira ML, Campino L, Cardoso L. Molecular detection of bacterial and parasitic pathogens in hard ticks from Portugal. Ticks Tick Borne Dis. 2014;5:409–14.View ArticleGoogle Scholar
- Shaw SE, Day MJ, Birtles RJ, Breitschwerdt EB. Tick-borne infectious diseases of dogs. Trends Parasitol. 2001;17:74–80.View ArticleGoogle Scholar
- de la Fuente J, Estrada-Pena A, Venzal JM, Kocan KM, Sonenshine DE. Overview: ticks as vectors of pathogens that cause disease in humans and animals. Front Biosci. 2008;13:6938–46.View ArticleGoogle Scholar
- Abd Rani PA, Irwin PJ, Coleman GT, Gatne M, Traub RJ. A survey of canine tick-borne diseases in India. Parasit Vectors. 2011;4:141.View ArticleGoogle Scholar
- Ahmad AS, Saeed MA, Rashid I, Ashraf K, Shehzad W, Traub RJ, et al. Molecular characterization of Hepatozoon canis from farm dogs in Pakistan. Parasitol Res. 2018;117:1131–8.View ArticleGoogle Scholar
- Dantas-Torres F, Latrofa MS, Annoscia G, Giannelli A, Parisi A, Otranto D. Morphological and genetic diversity of Rhipicephalus sanguineus sensu lato from the New and Old Worlds. Parasit Vectors. 2013;6:213.View ArticleGoogle Scholar
- Beugnet F, Marié JL. Emerging arthropod-borne diseases of companion animals in Europe. Vet Parasitol. 2009;163:298–305.View ArticleGoogle Scholar
- Dantas-Torres F, Chomel BB, Otranto D. Ticks and tick-borne diseases: a One Health perspective. Trends Parasitol. 2012;28:437–46.View ArticleGoogle Scholar
- Wolfe ND, Dunavan CP, Diamond J. Origins of major human infectious diseases. Nature. 2007;447:279–83.View ArticleGoogle Scholar
- Gray J, Dantas-Torres F, Estrada-Peña A, Levin M. Systematics and ecology of the brown dog tick, Rhipicephalus sanguineus. Ticks Tick Borne Dis. 2013;4:171–80.View ArticleGoogle Scholar
- Chisu V, Masala G, Foxi C, Socolovschi C, Raoult D, Parola P. Rickettsia conorii israelensis in Rhipicephalus sanguineus ticks, Sardinia, Italy. Ticks Tick Borne Dis. 2014;5:446–8.View ArticleGoogle Scholar
- Cicuttin GL, Brambati DF, Rodriguez Eugui JI, Gonzalez Lebrero C, De Salvo MN, Beltran FJ, et al. Molecular characterization of Rickettsia massiliae and Anaplasma platys infecting Rhipicephalus sanguineus ticks and domestic dogs, Buenos Aires (Argentina). Ticks Tick Borne Dis. 2014;5:484–8.View ArticleGoogle Scholar
- Dzul-Rosado K, Lugo-Caballero C, Tello-Martin R, López-Avila K, Zavala-Castro J. Direct evidence of Rickettsia typhi infection in Rhipicephalus sanguineus ticks and their canine hosts. Open Vet J. 2017;7:165–9.View ArticleGoogle Scholar
- Ramos R, Latrofa MS, Giannelli A, Lacasella V, Campbell BE, Dantas-Torres F, Otranto D. Detection of Anaplasma platys in dogs and Rhipicephalus sanguineus group ticks by a quantitative real-time PCR. Vet Parasitol. 2014;205:285–8.View ArticleGoogle Scholar
- Chomel B. Tick-borne infections in dogs-an emerging infectious threat. Vet Parasitol. 2011;179:294–301.View ArticleGoogle Scholar
- de Miranda RL, de Castro JR, Olegário MM, Beletti ME, Mundim AV, O’Dwyer LH, et al. Oocysts of Hepatozoon canis in Rhipicephalus (Boophilus) microplus collected from a naturally infected dog. Vet Parasitol. 2011;177:392–6.Google Scholar
- Solano-Gallego L, Sainz Á, Roura X, Estrada-Peña A, Miró G. A review of canine babesiosis: the European perspective. Parasit Vectors. 2016;9:336.View ArticleGoogle Scholar
- Farooqi AB, Khan AH, Mir H. Climate change perspective in Pakistan. Pakistan J Meteorol. 2005;2:3.Google Scholar
- Maps. Maps of World. 2017. https://www.mapsofworld.com/lat_long/pakistan-lat-long.html. Accessed 18 Sept 2017.
- Bashir I, Chaudhry Z, Ahmed S, Saeed M. Epidemiological and vector identification studies on canine babesiosis. Pak Vet J. 2009;29:51–4.Google Scholar
- 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:e0005681.View ArticleGoogle Scholar
- Ul-Hasan M, Abubakar M, Muhammad G, Khan MN, Hussain M. Prevalence of tick infestation (Rhipicephalus sanguineus and Hyalomma anatolicum anatolicum) in dogs in Punjab, Pakistan. Vet Ital. 2012;48:95–8.PubMedGoogle Scholar
- Ahmad S, Khan M, Khan M. Prevalence of canine babesiosis in Lahore, Pakistan. J Anim Plant Sci. 2007;17:1–2.Google Scholar
- Ahmad S, Khan M, Khan M. Epidemilogy and seasonal abundance of canine babesiosis in Lahore, Pakistan. J Anim Plant Sci. 2011;21:351–3.Google Scholar
- Shah SSA, Khan MI, Rafiullah KMA, Khan H, Ali A, et al. Tick-borne diseases-possible threat to humans-dog interspecie bond. Adv Anim Vet Sci. 2017;5:115–20.View ArticleGoogle Scholar
- Qamar M, Malik MI, Latif M, Ain QU, Aktas M, Shaikh RS, Iqbal F. Molecular detection and prevalence of Hepatozoon canis in dogs from Punjab (Pakistan) and hematological profile of infected dogs. Vector Borne Zoonotic Dis. 2017;17:179–84.View ArticleGoogle Scholar
- Moutailler S, Valiente Moro C, Vaumourin E, Michelet L, Tran FH, Devillers E, et al. Co-infection of ticks: the rule rather than the exception. PLoS Negl Trop Dis. 2016;10:e0004539.View ArticleGoogle Scholar
- Prusinski M, Kokas J, Hukey K, Kogut S, Lee J, Backenson P. Prevalence of Borrelia burgdorferi (Spirochaetales: Spirochaetaceae), Anaplasma phagocytophilum (Rickettsiales: Anaplasmataceae), and Babesia microti (Piroplasmida: Babesiidae) in Ixodes scapularis (Acari: Ixodidae) collected from recreational lands in the Hudson Valley region, New York State. J Med Entomol. 2014;51:226–36.Google Scholar
- Zając V, Wójcik-Fatla A, Sawczyn A, Cisak E, Sroka J, Kloc A, et al. Prevalence of infections and co-infections with 6 pathogens in Dermacentor reticulatus ticks collected in eastern Poland. Ann Agric Environ Med. 2017;24:26–32.View ArticleGoogle Scholar
- Gonçalves LR, Filgueira KD, Mendes Ahid SM, Pereira JS, do Vale AM, Machalo RZ, André MR. Study on coinfecting vector-borne pathogens in dogs and ticks in Rio Grande do Norte, Brazil. Rev Bras Parasitol Vet. 2014;23:407–12.View ArticleGoogle Scholar
- Latrofa MS, Dantas-Torres F, Giannelli A, Otranto D. Molecular detection of tick-borne pathogens in Rhipicephalus sanguineus group ticks. Ticks Tick Borne Dis. 2014;5:943–6.View ArticleGoogle Scholar
- Zhang J, Liu Q, Wang D, Li W, Beugnet F, Zhou J. Epidemiological survey of ticks and tick-borne pathogens in pet dogs in south-eastern China. Parasite. 2017;24:35.View ArticleGoogle Scholar
- Adamu M, Troskie M, Oshadu DO, Malatji DP, Penzhorn BL, Matjila PT. Occurrence of tick-transmitted pathogens in dogs in Jos, Plateau State, Nigeria. Parasit Vectors. 2014;7:119.View ArticleGoogle Scholar
- Attipa C, Hicks CAE, Barker EN, Christodoulou V, Neofytou K, Mylonakis ME, et al. Canine tick-borne pathogens in Cyprus and a unique canine case of multiple co-infections. Ticks Tick Borne Dis. 2017;8:341–6.View ArticleGoogle Scholar
- Kordick SK, Breitschwerdt EB, Hegarty BC, Southwick KL, Colitz CM, Hancock SI, et al. Coinfection with multiple tick-borne pathogens in a Walker Hound kennel in North Carolina. J Clin Microbiol. 1999;37:2631–8.PubMedPubMed CentralGoogle Scholar
- Starkey LA, Newton K, Brunker J, Crowdis K, Edourad EJP, Meneus P, Little SE. Prevalence of vector-borne pathogens in dogs from Haiti. Vet Parasitol. 2016;224:7–12.View ArticleGoogle Scholar
- Michelet L, Delannoy S, Devillers E, Umhang G, Aspan A, Juremalm M, et al. High-throughput screening of tick-borne pathogens in Europe. Front Cell Infect Microbiol. 2014;4:103.View ArticleGoogle Scholar
- Liu J, Hansen C, Quake SR. Solving the “world-to-chip” interface problem with a microfluidic matrix. Anal Chem. 2003;75:4718–23.View ArticleGoogle Scholar
- Pérez-Eid C. Les tiques. Identification, biologie, importance medicale et veterinaire. Paris: Editions Tec & Doc; 2007.Google Scholar
- Cabezas-Cruz A, Valdés JJ, de la Fuente J. The glycoprotein TRP36 of Ehrlichia sp. UFMG-EV and related cattle pathogen Ehrlichia sp. UFMT-BV evolved from a highly variable clade of E. canis under adaptive diversifying selection. Parasit Vectors. 2014;7:584.View ArticleGoogle Scholar
- Hsieh Y-C, Lee C-C, Tsang C-L, Chung Y-T. Detection and characterization of four novel genotypes of Ehrlichia canis from dogs. Vet Microbiol. 2010;146:70–5.View ArticleGoogle Scholar
- Zweygarth E, Cabezas-Cruz A, Josemans AI, Oosthuizen MC, Matjila PT, Lis K, et al. In vitro culture and structural differences in the major immunoreactive protein gp36 of geographically distant Ehrlichia canis isolates. Ticks Tick Borne Dis. 2014;5:423–31.View ArticleGoogle Scholar
- Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27:573–80.View ArticleGoogle Scholar
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.View ArticleGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.View ArticleGoogle Scholar
- Tamura K. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+ C-content biases. Mol Biol Evol. 1992;9:678–87.PubMedGoogle Scholar
- Cabezas-Cruz A, Vayssier-Taussat M, Greub G. Tick-borne pathogen detection: what’s new? Microbes Infect. 2018;S1286-4579:30004–2.Google Scholar
- Raulf M-K, Jordan D, Fingerle V, Strube C. Association of Borrelia and Rickettsia spp. and bacterial loads in Ixodes ricinus ticks. Ticks Tick Borne Dis. 2018;9:18–24.View ArticleGoogle Scholar
- Xu G, Mather TN, Hollingsworth CS, Rich SM. Passive surveillance of Ixodes scapularis (Say), their biting activity, and associated pathogens in Massachusetts. Vector Borne Zoonotic Dis. 2016;16:520–7.View ArticleGoogle Scholar
- Asman M, Solarz K, Cuber P, Gąsior T, Szilman P, Szilman E, et al. Detection of protozoans Babesia microti and Toxoplasma gondii and their co-existence in ticks (Acari: Ixodida) collected in Tarnogórski district (Upper Silesia, Poland). Ann Agric Environ Med. 2015;22:80–3.View ArticleGoogle Scholar
- Kjelland V, Paulsen KM, Rollum R, Jenkins A, Stuen S, Soleng A, et al. Tick-borne encephalitis virus, Borrelia burgdorferi sensu lato, Borrelia miyamotoi, Anaplasma phagocytophilum and Candidatus Neoehrlichia mikurensis in Ixodes ricinus ticks collected from recreational islands in southern Norway. Ticks Tick Borne Dis. 2018;9:1098–102.Google Scholar
- Koetsveld J, Tijsse-Klasen E, Herremans T, Hovius JW, Sprong H. Serological and molecular evidence for spotted fever group Rickettsia and Borrelia burgdorferi sensu lato co-infections in The Netherlands. Ticks Tick Borne Dis. 2016;7:371–7.View ArticleGoogle Scholar
- Sytykiewicz H, Karbowiak G, Chorostowska-Wynimko J, Szpechciński A, Supergan-Marwicz M, Horbowicz M, et al. Coexistence of Borrelia burgdorferi s.l. genospecies within Ixodes ricinus ticks from central and eastern Poland. Acta Parasitol. 2015;60:654–61.View ArticleGoogle Scholar
- Tappe J, Jordan D, Janecek E, Fingerle V, Revisited SC. Borrelia burgdorferi sensu lato infections in hard ticks (Ixodes ricinus) in the city of Hanover (Germany). Parasit Vectors. 2014;7:441.View ArticleGoogle Scholar
- Dantas-Torres F. The brown dog tick, Rhipicephalus sanguineus (Latreille, 1806) (Acari: Ixodidae): from taxonomy to control. Vet Parasitol. 2008;152:173–85.View ArticleGoogle Scholar
- Diuk-Wasser MA, Vannier E, Krause PJ. Coinfection by Ixodes tick-borne pathogens: ecological, epidemiological, and clinical consequences. Trends Parasitol. 2016;32:30–42.View ArticleGoogle Scholar
- Betts A, Gray C, Zelek M, MacLean R, King K. High parasite diversity accelerates host adaptation and diversification. Science. 2018;360:907–11.View ArticleGoogle Scholar
- Cabezas-Cruz A, Estrada-Peña A, Rego RO, de la Fuente J. Tick-pathogen ensembles: do molecular interactions lead ecological innovation? Front Cell Infect Microbiol. 2017;7:74.PubMedPubMed CentralGoogle Scholar
- Castañeda-Ortiz EJ, Ueti MW, Camacho-Nuez M, Mosqueda JJ, Mousel MR, Johnson WC, Palmer GH. Association of Anaplasma marginale strain superinfection with infection prevalence within tropical regions. PLoS One. 2015;10:e0120748.View ArticleGoogle Scholar
- Noh SM, Dark MJ, Reif KE, Ueti MW, Kappmeyer LS, Scoles GA, et al. Superinfection exclusion of the ruminant pathogen Anaplasma marginale in the tick vector is dependent on time between exposures to the strains. Appl Environ Microbiol. 2016;82:3217–24.View ArticleGoogle Scholar
- Kinnula H, Mappes J, Sundberg LR. Coinfection outcome in an opportunistic pathogen depends on the inter-strain interactions. BMC Evol Biol. 2017;17:77.View ArticleGoogle Scholar
- Abdullah AS, Moffat CS, Lopez-Ruiz FJ, Gibberd MR, Hamblin J, Zerihun A. Host-multi-pathogen warfare: pathogen interactions in co-infected plants. Front Plant Sci. 2017;8:1806.View ArticleGoogle Scholar
- Zhang X, Zhang X, Luo T, Keysary A, Baneth G, Miyashiro S, et al. Genetic and antigenic diversities of major immunoreactive proteins in globally distributed Ehrlichia canis strains. Clin Vaccine Immunol. 2008;15:1080–8.View ArticleGoogle Scholar
- Stich RW, Blagburn BL, Bowman DD, Carpenter C, Cortinas MR, Ewing SA, et al. Quantitative factors proposed to influence the prevalence of canine tick-borne disease agents in the United States. Parasit Vectors. 2014;7:417.View ArticleGoogle Scholar
- Masatani T, Hayashi K, Andoh M, Tateno M, Endo Y, Asada M, et al. Detection and molecular characterization of Babesia, Theileria, and Hepatozoon species in hard ticks collected from Kagoshima, the southern region in Japan. Ticks Tick Borne Dis. 2017;8:581–7.View ArticleGoogle Scholar
- Loh SM, Gofton AW, Lo N, Gillett A, Ryan UM, Irwin PJ, Oskam CL. Novel Borrelia species detected in echidna ticks, Bothriocroton concolor, in Australia. Parasit Vectors. 2016;9:339.View ArticleGoogle Scholar
- Rar VA, Fomenko NV, Dobrotvorsky AK, Livanova NN, Rudakova SA, Fedorov EG, et al. Tickborne pathogen detection, western Siberia, Russia. Emerg Infect Dis. 2005;11:1708–15.View ArticleGoogle Scholar
- Duron O, Jourdain E, McCoy KD. Diversity and global distribution of the Coxiella intracellular bacterium in seabird ticks. Ticks Tick Borne Dis. 2014;5:557–63.View ArticleGoogle Scholar
- Regnery RL, Spruill CL, Plikaytis B. Genotypic identification of rickettsiae and estimation of intraspecies sequence divergence for portions of two rickettsial genes. J Bacteriol. 1991;173:1576–89.View ArticleGoogle Scholar