Molecular identification of badger-associated Babesia sp. DNA in dogs: updated phylogeny of piroplasms infecting Caniformia

Background Piroplasms are unicellular, tick-borne parasites. Among them, during the past decade, an increasing diversity of Babesia spp. has been reported from wild carnivores. On the other hand, despite the known contact of domestic and wild carnivores (e.g. during hunting), and a number of ixodid tick species they share, data on the infection of dogs with babesiae from other families of carnivores are rare. Methods In this study blood samples were collected from 90 dogs and five road-killed badgers. Ticks were also removed from these animals. The DNA was extracted from all blood samples, and from 33 ticks of badgers, followed by molecular analysis for piroplasms with PCR and sequencing, as well as by phylogenetic comparison of detected genotypes with piroplasms infecting carnivores. Results Eleven of 90 blood DNA extracts from dogs, and all five samples from badgers were PCR-positive for piroplasms. In addition to the presence of B. canis DNA in five dogs, sequencing identified the DNA of badger-associated “Babesia sp. Meles-Hu1” in six dogs and in all five badgers. The DNA of “Babesia sp. Meles-Hu1” occurred significantly more frequently in dogs often taken to forests (i.e. the preferred habitat of badgers in Hungary), than in dogs without this characteristic. Moreover, detection of DNA from this Babesia sp. was significantly associated with hunting dogs in comparison with dogs not used for hunting. Two PCR-positive dogs (in one of which the DNA of the badger-associated Babesia sp. was identified, whereas in the other the DNA of B. canis was present) showed clinical signs of babesiosis. Engorged specimens of both I. canisuga and I. hexagonus were collected from badgers with parasitaemia, but only I. canisuga contained the DNA of “Babesia sp. Meles-Hu1”. This means a significant association of the DNA from “Babesia sp. Meles-Hu1” with I. canisuga. Phylogenetically, “Babesia sp. Meles-Hu1” belonged to the “B. microti” group. Conclusions This is the first detection of the DNA from a badger-associated Babesia sp. in dogs, one of which also showed relevant clinical signs. Based on the number of dogs with blood samples containing the DNA of “Babesia sp. Meles-Hu1” in this study (i.e. exceeding the number of B. canis-positives), these findings should not be regarded as isolated cases. It is assumed that dogs, which are used for hunting or frequently visit forests, are more likely to be exposed to this piroplasm, probably as a consequence of infestation with I. canisuga from badgers or from the burrows of badgers. The above results suggest that “Babesia sp. Meles-Hu1” should be added to the range of piroplasms, which are naturally capable of infecting hosts from different families of Caniformia.


Background
Piroplasms (Apicomplexa: Piroplasmida) are unicellular, tick-borne parasites, which infect red and white blood cells of their vertebrate hosts [1]. Among them, Babesia spp. and Theileria spp. are the most significant from the veterinary point of view, in terms of both their geographical distribution and the range of affected host species.
In general, babesiae have long been regarded as host species-specific [2], but in the era of molecular biological methods the reported ranges of hosts susceptible to infection with a particular Babesia species have started to expand. When the currently known diversity of Babesia spp., infecting domestic carnivores as their typical hosts, is concerned, the number of relevant disease agents appears to be stable, despite the fact that during the past years the taxonomic rank and/or the naming of certain variants has changed or is still unresolved. Thus, in a worldwide context, dogs are susceptible to three species of the so-called large babesiae and three species of small babesiae [3].
These data highlight carnivores as a group important from the point of studying piroplasms, especially at the interface of domestic and wild carnivores, for which mutually infective Babesia spp. have been reported [12]. In the above context, the present study aimed to molecularly investigate piroplasms in a region of Hungary, where emerging tick species and tick-borne pathogens have been reported [13,14]. To achieve this, blood samples and ticks were collected from dogs and badgers, consequently analyzed with PCR and sequencing.

Methods
Blood and tick samples were collected from dogs and European badgers (Meles meles) originating from 20 locations (17 and 3 for dogs and badgers, respectively) in southwestern Hungary (Somogy county) between March and October (dogs) and January to April (badgers) of 2017. From 90 dogs, EDTA-anticoagulated blood samples were drawn under clinical conditions (after skin surface sterilization, with sterile instruments) by cephalic venipuncture. The majority of sampled dogs (n = 80) were selected randomly (from regular patients of the Veterinary Clinic in Csurgó) and appeared to be healthy, but ten showed at least one clinical sign relevant to babesiosis. Animal data (owner, sex, age, mode of keeping) were recorded. From five freshly road-killed badgers EDTA-anticoagulated blood samples were collected (with sterile needle and syringe) from the heart, shortly after being reported to the Veterinary Authority. EDTA blood samples were frozen at -20°C until further processing. In addition, ixodid ticks were removed from all animals with pointed tweezers then transferred into 96% ethanol for storage in separate vials according to host individuals. Tick species were identified according to standard keys [15], and ticks of the subgenus Pholeoixodes according to [16].
DNA was obtained using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) following the manufacturer's instruction and using extraction controls to monitor cross contamination of samples in each run. DNA was extracted from 200 μl blood of 90 dogs and 5 badgers, as well as from 33 ticks collected from badgers. These ticks belonged to two species, i.e. Ixodes canisuga and I. hexagonus (subgenus Pholeoixodes), selected on the basis of literature data supporting their role in the transmission of Babesia spp. of the B. microti group [17,18]. DNA was extracted from the ticks individually, with the same method as from the blood samples, but including an overnight digestion in tissue lysis buffer and proteinase-K at 56°C, and incubation in lysis buffer for 10 min at 70°C. DNA samples were finally taken up in 130 μl elution buffer.
Each PCR was run with positive and negative controls (i.e. sequence-verified DNA of Babesia canis, and nontemplate reaction mixture, respectively). PCR products were visualized in 1.5% agarose gel. Negative controls and extraction controls remained PCR negative in all tests. Purification and sequencing (directly from the PCR product, using the forward primer and generating the sequence twice per sample) were performed from all piroplasm PCRpositive samples at Biomi Inc. (Gödöllő, Hungary). Sequences were aligned and compared to reference Gen-Bank sequences by nucleotide BLASTn program (https:// blast.ncbi.nlm.nih.gov). All blood samples, which yielded unusual results (i.e. the presence of DNA from a Babesia sp. not yet reported in the relevant host) were re-tested in duplicates, by repeating the whole procedure from the beginning (DNA extraction, PCR and sequencing). Representative sequences were submitted to the GenBank database under the accession numbers MG778912-MG778915. Phylogenetic analyses were performed with the Maximum Likelihood method and Tamura-Nei model, using MEGA 6.0.
Rates of PCR-positivity were compared with the Fisher's exact test and differences were considered significant when P < 0.05.

Results
Among blood DNA extracts, 11 of 90 samples from dogs, and all five samples from badgers were PCR-positive for piroplasms. Sequencing identified two Babesia spp. in dogs. In six dog blood samples the DNA of a Babesia sp. was present, which was formerly reported from European badgers (Meles meles) in Hungary and designated as "Babesia sp. Meles-Hu1" (GenBank: KX218234). This genotype also corresponded to "Babesia sp. badger type-A" reported from Spain (GenBank: KT223484) and the UK (GenBank: KX528553), with 100% (472/472 bp) and 99.8% (471/472 bp) identity, respectively. In addition, the DNA of B. canis was detected in five blood samples from dogs: in one case corresponding to "genotype A" (430/430 bp; 100% identity with KP835549), and in four samples to "genotype B" (430/430 bp; 100% identity with KP835550) formerly reported in Hungary. In the blood of all badgers the DNA of "Babesia sp. Meles-Hu1" was present, which had 100% (472/472 bp) identity with the above dog isolate.
Concerning the distribution of Babesia species in dogs according to their mode of keeping, the DNA of "Babesia sp. Meles-Hu1" occurred significantly (P = 0.0001) more frequently in dogs often taken to forests (5 of 12) than in dogs without this characteristic (1 of 78) ( Table 1). Moreover, the presence of DNA from this Babesia sp. was significantly (P = 0.00008) associated with hunting dogs (4 of 6 were infected) in comparison with dogs not used for hunting (2 of 84) ( Table 1). Taking into account the sampling time, the dog samples containing the DNA of different Babesia spp. also showed an uneven seasonal distribution: while the DNA of "Babesia sp. Meles-Hu1" was detected in samples collected in March (n = 3) and August (n = 3), the DNA of B. canis was present in samples obtained in May (n = 2) or September (n = 3). Importantly, two PCR-positive dogs (one yielding the sequence of the badger-associated Babesia sp., the other that of B. canis) showed clinical signs of babesiosis (fatigue, renal failure and/or anaemia, icterus).
In the molecularly analyzed 33 Pholeoixodes tick samples, only the DNA of "Babesia sp. Meles-Hu1" was found. Interestingly, while engorged specimens of both I. canisuga (n = 27) and I. hexagonus (n = 6) were collected from badgers, which had this Babesia sp. in their blood (as shown above), none of the I. hexagonus specimens were PCR-positive, but 66.7% (18 of 27) of I. canisuga DNA extracts contained "Babesia sp. Meles-Hu1" ( Table 1). This was a significant association of this piroplasm with I. canisuga (P = 0.0045).
Phylogenetically, sequences of B. canis genotypes amplified from dogs in the present study clustered with other conspecific isolates of the clade "Babesia (sensu stricto)" (Fig. 1). On the other hand, sequences of "Babesia sp. Meles-Hu1" obtained from both dogs and badgers here belonged to the "B. microti phylogenetic group" (Fig. 1).

Discussion
In this study, blood and tick samples collected from dogs and European badgers were molecularly analyzed for the presence of piroplasm DNA. The DNA from both B. canis and "Babesia sp. Meles-Hu1" were found. Based on literature reviews, in European badgers "Babesia sp. badger type-A" and " Babesia sp. badger type-B", as well as "B. annae" were reported [18]. In dogs as typical hosts, in a worldwide context, three small Babesia spp., i.e. B. gibsoni, "B. annae" (syn. "T. annae", B. cf. microti)  (Fig. 1). Rare accounts also attest the occasional finding of the DNA of piroplasms in dogs, for which the typical hosts are outside the order Carnivora (e.g. B. caballi [10]). However, this is the first molecular identification of the DNA from a badger-associated Babesia sp. in canine hosts. Taking into account the simultaneous occurrence of DNA from "Babesia sp. Meles-Hu1" in six dogs, one of which showed typical clinical signs of babesiosis, results of the present study suggest that the occurrence of "Babesia sp. Meles-Hu1" is not exceptional in dogs, and it might even have a pathogenic role in this host species. The badger-associated Babesia genotype identified here corresponded to "Babesia sp. Meles-Hu1", which was formerly identified in a badger in eastern Hungary (GenBank: KX218234), suggesting a widespread occurrence of this species in the country. Babesia canis genotypes "A" and "B" (GenBank: KP835549 and KP835550, Fig. 1 Maximum Likelihood tree of piroplasm DNA sequences reported from carnivores, shown according to their families. Sequences from this study are highlighted with red color. Piroplasm categories are named after [1] (except for Theileria equi, Cytauxzoon spp. and Theileria spp. united in one group). For Hyaenidae, the sequence KF270649 had a short coverage and therefore is not included (but in a sequence comparison it aligned with B. lengau). Vertical red columns mark adjacent families of Caniformia. The scale-bar indicates the number of substitutions per site respectively) were also reported formerly in Hungary from bats [21] and from questing D. reticulatus ticks [22].
The significantly higher number of "Babesia sp. Meles-Hu1"-positive samples among dogs frequently taken to forests (including hunting dogs) highlights the importance of the preferred habitat of badgers (i.e. forest in Hungary [23]) in the epidemiology of this infection. In particular, dogs in such environments, especially when getting into contact with badgers during hunting, may acquire infestation with I. canisuga from badgers or from the burrows of badgers.
On the other hand, all B. canis-positive dogs were identified in the sample group not associated with forests, most likely reflecting that the vector of this piroplasm (D. reticulatus) is an "open country tick species" [24] (in Hungary [25]). In addition, the time shift between PCR-positive samples containing either B. canis or "Babesia sp. Meles-Hu1" is also most likely a consequence of differences in the activity periods of their vectors.
Ixodes canisuga was proposed to be the vector of the badger-associated piroplasm, B. missirolli [17], and this vector potential was confirmed by the present results from the molecular analysis of badger ticks (I. canisuga vs I. hexagonus). Dogs were formerly reported to harbor I. canisuga in Hungary [26]. Such occasions would allow tick-borne badger-to-dog transmission. On the other hand, taking into account the usage of dogs during badger hunting in the evaluated region of Hungary, oral infection of relevant dogs with "Babesia sp. Meles-Hu1" cannot be completely ruled out. In support of this, (i) blood-borne transmission of babesiae was reported in fighting dogs (B. gibsoni [27]), and (ii) oral infection is possible in the case of the type-species of the phylogenetic group where "Babesia sp. Meles-Hu1" belongs (B. microti [28]).
Phylogenetic analysis showed that "Babesia sp. Meles-Hu1" identified in dogs belongs to the B. microti group, where piroplasms of the Mustelidae, Canidae and Procyonidae (suborder Caniformia) cluster together and are separated (with 100% support) from a phylogenetic group containing piroplasms of Feliformia (Fig. 1). Several examples attest that different species of carnivores within a family (e.g. Canidae) are susceptible to the same piroplasm (e.g. B. canis), while Babesia spp. are also known to infect carnivorous hosts from different families of Caniformia, as exemplified by high prevalences of B. microti-like infections in both Canidae and Procyonidae [12].

Conclusions
This is the first detection of the DNA from a badgerassociated Babesia sp. in dogs, one of which also showed relevant clinical signs. In this study, the number of dogs with blood samples containing the DNA of "Babesia sp.
Meles-Hu1" exceeded the number of B. canis-positive samples, therefore these findings should not be regarded as isolated cases. It is assumed that dogs, which are used for hunting or frequently visit forests, are more likely to be exposed to this piroplasm, probably as a consequence of infestation with I. canisuga from badgers or from the burrows of badgers. These results suggest that "Babesia sp. Meles-Hu1" should be added to the range of piroplasms, which are naturally capable of infecting hosts from different families of Caniformia.