- Open Access
Low genetic diversity of Ehrlichia canis associated with high co-infection rates in Rhipicephalus sanguineus (s.l.)
Parasites & Vectorsvolume 12, Article number: 12 (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.
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,2,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,6,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,9,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,12,13,14,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,17,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,22,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,25,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,29,30], including R. sanguineus (s.l.) [31,32,33], and dogs were found infected with more than one tick-borne pathogen [34,35,36,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
Single engorged adult female ticks (n = 235) were collected from clinically healthy farm dogs of both sexes (male, n = 70; female, n = 165) between June and October 2016 from three different districts in the province of Punjab, including Kasur (31°12'21"N, 74°45'81"E; n = 87), Muzaffargarh (30°07'36"N, 71°18'05"E; n = 75) and Rawalpindi (33°59'84"N, 73°04'41"E; n = 73) (Fig. 1). Following collection, each tick specimen was stored individually in 70% ethanol until used. For morphological identification, each tick was examined using a dissecting microscope (Olympus, Tokyo, Japan). Ticks were identified using the keys as described previously .
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
Conventional and nested PCR using different primers than those of the BioMarkTM system were used to confirm the presence of the detected pathogens in the samples (Table 1). Amplicons were sequenced by Eurofins MWG Operon (Ebersberg, Germany) and assembled using the BioEdit software (Ibis Biosciences, Carlsbad, USA). An online BLAST (National Center for Biotechnology Information) was used to identify the sequenced organisms.
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,42,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
A total of 204 (of 235, 87%) ticks were infected with at least one pathogen. Rickettsia massiliae was the most commonly (68%) detected pathogen followed by E. canis (51%) and Anaplasma spp. (17%), including A. phagocytophilum (6%), A. marginale (4%), A. centrale (3%), A. platys (3%) and A. ovis (1%). Furthermore, Coxiella-like endosymbionts (hereafter Coxiella-like) and B. canis were detected in 5% and 3% of samples, respectively. Only a few pathogens were identified as single infections, including R. massiliae (24.3%), E. canis (12.8%), A. platys (0.9%), A. phagocytophilum (1.7%) and A. marginale (0.9%), respectively (Table 2).
Based on location, the highest prevalence (of mono-specific infection) was found in Muzaffargarh (98%) followed by Rawalpindi (84%) and Kasur (78%). The prevalence was significantly different between Muzaffargarh and Kasur (χ2 = 18.9, P < 0.0001), and Rawalpindi and Muzaffargarh (χ2 = 12.0, P < 0.003). The region with the lowest microorganism diversity was Muzaffargarh where only five of the nine microorganisms were detected (Fig. 2). Rickettsia massiliae, E. canis and A. platys were detected in all three regions studied. However, Coxiella-like was detected only in Muzaffargarh and Rawalpindi. Likewise, A. phagocytophilum, A. marginale and A. centrale were found in Kasur and Rawalpindi only (Fig. 2). The prevalence of E. canis was significantly different between Muzaffargarh and Kasur (χ2 = 4.7, P < 0.05), Rawalpindi and Muzaffargarh (χ2 = 34.9, P < 0.0001), and Rawalpindi and Kasur (χ2 = 19.9, P < 0.0001), although no differences were found for R. massiliae among regions.
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).
Multiple correspondence analysis was performed to further test correlation between infections and geographic regions. When only a single pathogen (regardless of being part of a co-infection) was considered, the MCA confirmed the variation in the presence and prevalence of single pathogens among various agroclimatic regions (Fig. 3a). When double and triple co-infections were considered, the MCA revealed associations between regions and co-infections. Based on double co-infection analyses, Kasur and Rawalpindi clustered together with co-infections, including R. massiliae + A. platys and R. massiliae + A. marginale, respectively (Fig. 3b). Muzaffargarh clustered separately and harbored co-infection of R. massiliae + E. canis. In addition, two pairs of co-infections, including A. phagocytophilum (i.e. R. massiliae + A. phagocytophilum and E. canis + A. phagocytophilum) and B. canis (i.e. E. canis + B. canis and R. massiliae + B. canis) clustered irrespective of any particular region (Fig. 3b). When triple infections were considered, Kasur, Rawalpindi and Muzaffargarh lied in the center of the multidimensional matrix, and around these points the clustered triple infections were of R. massiliae + E. canis + B. canis, and R. massiliae + E. canis + A. phagocytophilum, E. canis + A. centrale + Coxiella-like, and R. massiliae + E. canis + A. centrale, indicating an independent distribution of triple infections in three regions (Fig. 3c). In addition, a strong correspondence was observed between the triple infections R. massiliae + E. canis + A. ovis and R. massiliae + Coxiella-like + A. ovis, which also showed the highest values of inertia (less independence).
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.
The E. canis trp36 nucleotide and putative protein sequences reported in this study shared more than 99% and 97% identity at the nucleotide and amino acid levels, respectively (Table 3). The amino acid sequence variation of Pakistan strains relative to those from the USA (Louisiana) and Taiwan (TWN17) is displayed in Fig. 4a. These were included as representatives of low and high genetic diversity strains, respectively [41, 43]. The pre-repeat region of TRP36 amino acid sequences from Pakistan had molecular signatures of previously reported E. canis strains of low genetic diversity (Fig. 4a). Particularly, the presence of G at position 117 and the sequence of the putative glycosylation sequon at N125 was NPS (arrows in Fig. 4a). The tandem repeat region of the three TRP36 sequences obtained in this study contained 9 repetitions of the sequence ‘TEDSVSAPA’ (highlighted in yellow in Fig. 4a). Unexpectedly, the trp36 sequences from Muzaffargarh and Kasur were 100% identical at the nucleotide level but shared only 97.9% identity at the amino acid level (Table 3). Multiple sequence alignment revealed that the Muzaffargarh trp36 sequence had an insertion of thymine (t) at position 707 (Fig. 4b). This insertion shifted the open reading frame of trp36 which resulted in amino acid changes in the post-repeat region (Fig. 4a) and an early stop codon ‘TGA’ (Fig. 4b). As a result, while the post-repeat region in sequences from Kasur and Rawalpindi was 30 amino acids long, the post-repeat region in that from Muzaffargarh was only 16 amino acids long. Variable sizes (16, 28 and 30 amino acids) in the post-repeat regions of TRP36 have been described previously [41, 43].
A previous study reported that trp36 grouped E. canis strains in two well-defined phylogenetic clusters, i.e. one containing highly variable strains and the other with similar ones . A maximum likelihood phylogenetic analysis was then performed to test the phylogenetic position of E. canis strains from Pakistan. In agreement with the molecular signatures described above, Fig. 5 shows that the E. canis strains from Pakistan clustered in the low genetic diversity clade.
Because each region had one specific E. canis strain, and ticks in each region had different prevalence of single infection with E. canis and co-infection with R. massiliae and E. canis, we hypothesized that different E. canis strain may have different fitness for single infection or co-infection with R. massiliae between regions. To test this hypothesis, the prevalence of E. canis strains in single infections and co-infections with R. massiliae was compared between regions using Chi-square test (Table 4). Results showed that the prevalence of co-infection with E. canis and R. massiliae in Muzaffargarh was significantly higher than mono-infection of E. canis in Kasur and Rawalpindi.
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,12,13,14,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,32,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,52,53,54,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,32,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.)
Dantas-Torres F. Biology and ecology of the brown dog tick, Rhipicephalus sanguineus. Parasit Vectors. 2010;3:26.
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.
Shaw SE, Day MJ, Birtles RJ, Breitschwerdt EB. Tick-borne infectious diseases of dogs. Trends Parasitol. 2001;17:74–80.
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.
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.
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.
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.
Beugnet F, Marié JL. Emerging arthropod-borne diseases of companion animals in Europe. Vet Parasitol. 2009;163:298–305.
Dantas-Torres F, Chomel BB, Otranto D. Ticks and tick-borne diseases: a One Health perspective. Trends Parasitol. 2012;28:437–46.
Wolfe ND, Dunavan CP, Diamond J. Origins of major human infectious diseases. Nature. 2007;447:279–83.
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.
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.
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.
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.
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.
Chomel B. Tick-borne infections in dogs-an emerging infectious threat. Vet Parasitol. 2011;179:294–301.
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.
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.
Farooqi AB, Khan AH, Mir H. Climate change perspective in Pakistan. Pakistan J Meteorol. 2005;2:3.
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.
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.
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.
Ahmad S, Khan M, Khan M. Prevalence of canine babesiosis in Lahore, Pakistan. J Anim Plant Sci. 2007;17:1–2.
Ahmad S, Khan M, Khan M. Epidemilogy and seasonal abundance of canine babesiosis in Lahore, Pakistan. J Anim Plant Sci. 2011;21:351–3.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Liu J, Hansen C, Quake SR. Solving the “world-to-chip” interface problem with a microfluidic matrix. Anal Chem. 2003;75:4718–23.
Pérez-Eid C. Les tiques. Identification, biologie, importance medicale et veterinaire. Paris: Editions Tec & Doc; 2007.
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.
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.
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.
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27:573–80.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
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.
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.
Cabezas-Cruz A, Vayssier-Taussat M, Greub G. Tick-borne pathogen detection: what’s new? Microbes Infect. 2018;S1286-4579:30004–2.
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.
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.
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.
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.
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.
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.
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.
Dantas-Torres F. The brown dog tick, Rhipicephalus sanguineus (Latreille, 1806) (Acari: Ixodidae): from taxonomy to control. Vet Parasitol. 2008;152:173–85.
Diuk-Wasser MA, Vannier E, Krause PJ. Coinfection by Ixodes tick-borne pathogens: ecological, epidemiological, and clinical consequences. Trends Parasitol. 2016;32:30–42.
Betts A, Gray C, Zelek M, MacLean R, King K. High parasite diversity accelerates host adaptation and diversification. Science. 2018;360:907–11.
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.
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.
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.
Kinnula H, Mappes J, Sundberg LR. Coinfection outcome in an opportunistic pathogen depends on the inter-strain interactions. BMC Evol Biol. 2017;17:77.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.