Skip to main content

Role of mustelids in the life-cycle of ixodid ticks and transmission cycles of four tick-borne pathogens



Elucidating which wildlife species significantly contribute to the maintenance of Ixodes ricinus populations and the enzootic cycles of the pathogens they transmit is imperative in understanding the driving forces behind the emergence of tick-borne diseases. Here, we aimed to quantify the relative contribution of four mustelid species in the life-cycles of I. ricinus and Borrelia burgdorferi (sensu lato) in forested areas and to investigate their role in the transmission of other tick-borne pathogens. Road-killed badgers, pine martens, stone martens and polecats were collected in Belgium and the Netherlands. Their organs and feeding ticks were tested for the presence of tick-borne pathogens.


Ixodes hexagonus and I. ricinus were found on half of the screened animals (n = 637). Pine martens had the highest I. ricinus burden, whereas polecats had the highest I. hexagonus burden. We detected DNA from B. burgdorferi (s.l.) and Anaplasma phagocytophilum in organs of all four mustelid species (n = 789), and Neoehrlichia mikurensis DNA was detected in all species, except badgers. DNA from B. miyamotoi was not detected in any of the investigated mustelids. From the 15 larvae of I. ricinus feeding on pine martens (n = 44), only one was positive for B. miyamotoi DNA, and all tested negative for B. burgdorferi (s.l.), N. mikurensis and A. phagocytophilum. The two feeding larvae from the investigated polecats (n = 364) and stone martens (n = 39) were negative for all four pathogens. The infection rate of N. mikurensis was higher in feeding nymphs collected from mustelids compared to questing nymphs, but not for B. burgdorferi (s.l.), B. miyamotoi or A. phagocytophilum.


Although all stages of I. ricinus can be found on badgers, polecats, pine and stone martens, their relative contribution to the life-cycle of I. ricinus in forested areas is less than 1%. Consequently, the relative contribution of mustelids to the enzootic cycles of I. ricinus-borne pathogens is negligible, despite the presence of these pathogens in organs and feeding ticks. Interestingly, all four mustelid species carried all stages of I. hexagonus, potentially maintaining enzootic cycles of this tick species apart from the cycle involving hedgehogs as main host species.


Tick-borne diseases such as Lyme borreliosis (LB) and tick-borne encephalitis pose serious health concerns in Europe [1, 2]. Further geographical spread and long-lasting increases in the incidence of these two diseases have been observed in several European countries. Furthermore, diseases involving other tick-borne pathogens such as Anaplasma phagocytophilum, Borrelia miyamotoi and Neoehrlichia mikurensis are emerging or being (re)discovered [3]. Some of these pathogens are also of veterinary relevance, particularly A. phagocytophilum, which causes disease in dogs, horses and domesticated ruminants [4,5,6,7,8]. Understanding which factors drive population densities of ticks and the transmission cycles of these pathogens are important steps in assessing disease risk and formulating possible intervention strategies. One of these factors is the vertebrate community: wildlife and free ranging domestic animals act as feeding and propagation hosts of I. ricinus and as reservoir of tick-borne pathogens in nature areas [9, 10].

The causative agents of Lyme borreliosis are spirochetes belonging to the Borrelia burgdorferi (sensu lato) complex. At least five genospecies of B. burgdorferi (s.l.) have been shown to be pathogenic: B. afzelii, B. garinii, B. burgdorferi (sensu stricto), B. bavariensis and B. spielmanii [11, 12]. These genospecies are maintained through distinct enzootic cycles involving ticks and vertebrates acting as amplifying hosts [13]. For example, B. afzelii is mainly transmitted by small mammals, while B. garinii is mainly transmitted by birds [14, 15]. Even within genospecies, host species differ in their ability to transmit B. burgdorferi (s.l.) [16]. In Europe, these pathogens are predominantly transmitted by Ixodes ricinus. This generalist tick species has a three-host life-cycle; it depends on three vertebrate hosts for blood meals. Although I. ricinus utilizes a multitude of host species, these host species differ considerably in the number of ticks and the different life stages they feed [17]. Larvae mainly feed on small rodents and nymphs on multiple host species, whereas deer, most often roe deer (Capreolus capreolus), are the main hosts for adults [18, 19]. Finding out which host species contribute most in maintaining I. ricinus populations and transmitting B. burgdorferi (s.l.) is imperative in understanding the driving forces behind the enzootic cycles of this bacterial complex.

Both the number of ticks that a host successfully feeds and its ability to maintain and transmit a particular B. burgdorferi (s.l.) genospecies (i.e. reservoir competence) are important host characteristics [19,20,21]. The relative contribution to maintain B. burgdorferi (s.l.) in enzootic cycles, therefore, also depends on the density of the various vertebrate host species in the environment. Thus, the resulting distribution of B. burgdorferi (s.l.) in the questing ticks is a function of the densities of different host species, their capacity to feed ticks and to transmit the bacteria to those ticks. We previously performed an extensive systematic review to quantify the contribution of various vertebrate host species to feeding the different I. ricinus life stages, and transmitting B. burgdorferi (s.l.) genospecies to feeding larvae [18]. As the knowledge of the vertebrate assemblage in each environment is essential to perform these quantifications, we used data from the literature to approximate the vertebrate assemblage, including the most widespread vertebrate species occurring in most European forests. An important limitation of that review was that no, or only little, information was available for some widespread host species. In this study, we aimed to collect part of the missing data, namely that of several species of mustelids.

The family of mustelids (Mustelidae) consists of medium-sized to small carnivores, such as the Eurasian badger (Meles meles), the Eurasian otter (Lutra lutra), the European pine marten (Martes martes), the stone marten (Martes foina), the European polecat (Mustela putorius) and the common weasel (Mustela nivalis). Some mustelid species are closely bound to specific environments, such as the Eurasian otter to water, while others are habitat generalists occurring in most habitat types, such as the common weasel and the European polecat [23]. Of all European mustelid species, there are several that commonly occur in forested habitat and their ecotones and could therefore be important hosts for I. ricinus, namely the Eurasian badger, European pine marten, European polecat and stone marten [24]. These species are hereinafter referred to as badger, pine marten, polecat and stone marten.

Mustelids are often infested with I. ricinus [18, 25, 26], but can also carry Ixodes hexagonus [27] and other tick species [28]. Ixodes hexagonus is a nest-dwelling, three-host tick species. It parasitizes hedgehogs and carnivorous mammals from the families Canidae, Mustelidae and Felidae [29]. This tick species does not show any host preference between the adult and immature stages [27] and has been reported to bite humans as well [27, 30]. The vector capacity of I. hexagonus has been experimentally demonstrated for B. burgdorferi (s.l.) and tick-borne encephalitis virus [31,32,33], and has been suspected for other I. ricinus-borne pathogens as well, but without experimental proof [34,35,36].

In this study, we aimed to quantify the relative importance of four mustelid species hosting the different life stages of I. ricinus. We estimated the tick burden of these mustelids for other tick species and investigated their role in the transmission cycles of four tick-borne pathogens: B. burgdorferi (s.l.), A. phagocytophilum, B. miyamotoi and N. mikurensis. For B. burgdorferi (s.l.), we aimed to quantify the realized reservoir competence of these mustelid species, which is defined as the proportion of blood-fed larvae that become infected with a pathogen [37].


Sample collection in the Netherlands

Mustelids killed by traffic or otherwise found dead were collected by volunteers throughout the Netherlands between 1995 and 2015. Animals were collected throughout the year in a range of habitats, from forest to agricultural fields and urbanized areas. Only animals that were perceived as recently killed and intact were taken into account. These animals were transported in a sealed plastic bag, ensuring that the ectoparasites that were still present on the animals were preserved (Figure 1). The animals were transferred to the Dutch Wildlife Health Centre (DWHC; 6 badgers), the Dutch Mammal Society or Wageningen Environmental Research (formerly Alterra). The animals sent to the DWHC were stored at 4 °C overnight and checked for ticks and dissected the following morning. All other animals were stored at -20 °C until further analysis. The animals were screened for ticks by thoroughly inspecting the skin and hairs of the head, neck, axillae, inguinal region and legs. The rest of the body was roughly inspected for ectoparasites, as was the contents of the transportation plastic bags. Ectoparasites were collected in a single tube with 70% ethanol per individual. These tubes were stored at -20 °C until further analysis. After screening for ectoparasites, tissue samples were taken from a subset of animals. Small parts of the liver and spleen were collected in 2 ml tubes when these organs were still intact after the impact of the vehicle. Ear biopsies were collected only from some of the animals. The dissection material and table were cleaned thoroughly between animals.

Fig. 1
figure 1

Road-killed polecat with fed ticks from collection bag. Inset: Feeding I. hexagonus female on pine marten

Sample collection in Belgium

Analogous to the Netherlands, road-killed mustelids were collected since 1995 by a voluntary network (Marternetwerk; throughout Flanders, i.e. the northern part of Belgium, hereinafter referred to as Belgium. Systematic attention for ticks was given in the period until 2005. As pine martens were rare in Belgium until recent years [38], no specimens of this species were available for this study. Animals were collected throughout the year in a similar range of habitats as in the Netherlands. Dead animals were transferred in a sealed plastic bag to a local network freezer at -20 °C and later autopsied at the Research Institute for Nature and Forest (INBO). The screening of dead animals for ticks was performed as described above. The ectoparasites were collected in a small sealed plastic bag and stored at -20 °C. During the autopsy, liver samples were taken if possible and stored at -20 °C (no spleen samples or ear biopsies were taken). Not all animals from which a liver sample was taken were also screened for ticks, resulting in different sample sizes for tick burden and infection rate in organs. The dissection material and table were cleaned thoroughly between animals.

DNA extraction, qPCR assays and sequencing procedures

DNA from ticks and tissues were extracted using the Qiagen DNeasy Blood and Tissue Kit (Qiagen, Venlo, the Netherlands) according to the manufacturer’s protocol [39]. For the detection of B. burgdorferi (s.l.) DNA, a duplex qPCR was used, based on the detection of fragments of the outer surface protein A (ospA) and flagellin genes [40]. A conventional PCR assay, targeting the 5S-23S intergenic spacer region (IGS), was performed for B. burgdorferi (s.l.) genospecies identification [41]. Conventional PCR assays were carried out in a Px2 thermal cycler (Thermo Electron Corporation, Breda, the Netherlands) and visualized on a 2% agarose gel. Both strands of PCR products were sequenced by BaseClear (Leiden, the Netherlands), according to the company’s protocol and using the same forward and reverse primers as in the conventional PCR [42]. BLAST analyses and in-house molecular epidemiological databases (Bionumerics 7.1 Applied Math, Sint-Martens-Latem, Belgium) were used to identify B. burgdorferi (s.l.) genospecies. These databases contain all our DNA sequences from (field) isolates, together with (reference) sequences from GenBank [11, 41]. For detection of B. miyamotoi, a qPCR assay was used that targets a region of the flagellin gene, specific for B. miyamotoi [43]. For detection of A. phagocytophilum and N. mikurensis DNA, a duplex qPCR assay was used, as described by Jahfari et al. [44, 45]. This qPCR assay targets specific regions of the major surface protein 2 gene (msp2) for A. phagocytophilum and a heat-shock protein gene (groEL) for N. mikurensis. All qPCR runs were carried out in a final volume of 20 μl, containing IQ Multiplex Powermix (Bio-Rad Laboratories, Hercules, CA, USA), 400 nM of primers and hydrolysis probes and 3 μl of DNA template. Conditions for PCR amplification were as follows: 95 °C for 5 min, 60 cycles at 95 °C for 5 s and 60 °C for 35 s, followed by a final incubation step at 37 °C for 20 s. We carried out qPCR assays on a LightCycler 480 instrument (Roche Diagnostics Nederland B.V, Almere, the Netherlands) and analysis was performed by the instrument’s software (release Quantification cycle (Cq) values were calculated using the second derivative method.

Quantifying the relative importance of mustelids as hosts for I. ricinus and B. burgdorferi (s.l.)

We quantified the relative importance of the four mustelid species using the method described by Hofmeester et al. [18]. Briefly, we modeled a simplified host community based on species presence and densities in European forests, using the same species composition as [18] with the addition of the four mustelid species described in this study. For each of these species, we calculated the relative importance of the species in feeding I. ricinus ticks of the different life stages (\( {RI}_{l_i} \), \( {RI}_{n_i} \), and \( {RI}_{a_i} \) for the larval, nymphal and adult stage, respectively) as,

$$ {RI}_{l_i}=\frac{B_{l_i}{D}_i}{\sum_{j=1}^n{B}_{l_j}{D}_j} $$

where \( {B}_{l_i} \) is the mean larval burden (no. of larvae per host individual) of species i, and Di is the average density (km-2) in which species i occurs. To get a relative proportion, the total number of larvae fed by all species in the community (\( {\sum}_{j=1}^n{B}_{l_j}{D}_j \)) is used. \( {B}_{l_i} \) can be substituted by \( {B}_{n_i} \) (the mean nymphal burden) or \( {B}_{a_i} \) (the mean adult burden) to calculate \( {RI}_{n_i} \) and \( {RI}_{a_i} \), respectively.

Following Hofmeester et al. [18], we added the realized reservoir competence for B. burgdorferi (s.l.) (\( {RC}_{bb_i} \)) to Equation 1 to estimate the relative importance of each species in infecting I. ricinus larvae with B. burgdorferi (s.l.) (\( {RI}_{bb_i} \)).

As our study contained animals that were collected from a range of habitat types and we wanted to model a host community in a forested habitat, we estimated I. ricinus burdens for the four species of mustelid described in this study by modelling tick burden as a function of forest cover. First, we used the Pan-European forest map for 2000 [46] and estimated the percentage forest cover within a buffer with a 564.2 m radius (~1 km2: the average home range size of the studied species [23]). We did this for a subset of the roadkill samples with information on the exact location where the animal was found and used the map of 2000 as this year overlapped with both the sampling in the Netherlands and Belgium. Secondly, we modelled individual I. ricinus burden as a function of species and percentage forest cover using a generalized linear model (GLM) with a negative binomial distribution and a log link function and used the model fit to predict the tick burden for each of the four species in a site with 100% forest cover. We combined estimates of female and male I. ricinus to estimate the adult I. ricinus burden. In addition to these predicted tick burdens, we used average densities for badger, polecat and pine marten as published in the PanTHERIA database [47]. Stone marten density was not available from PanTHERIA and there are, to our knowledge, no recent estimates of average stone marten densities in forested habitats. We therefore used the same estimate as for the related pine marten.

Statistical analyses

We tested for differences between the species and the two countries in which animals were collected assuming that biases due to habitat and season were similar between species and countries. We did this to be able to use the full dataset, as only a subset contained full information on location and date when the animal was collected. We used GLMs with a negative binomial distribution and a log link function to test for differences in tick burden between the different mustelid species while accounting for overdispersion in the data. We also took into account potential differences between the two countries as a factor, resulting in GLMs with two factors as independent variables: species and country. Each tick species and life stage was tested separately except for males of both tick species, due to inadequate sample size. Similarly, we used GLMs with a binomial distribution and a log link function to test for differences in infection rate between the different mustelid species. Again, we included country as a factor in all models, resulting in GLMs with two factors as independent variables, species and country, and ran the same model for each pathogen. We defined animals as being infected when at least one tissue sample was found positive for the pathogen. All models were run in R 3.4.4. We performed a Tukey HSD post-hoc test to test for differences between the species using the multcomp package [48].

In a second part of the analysis, we used GLMs with a binomial distribution and a logit link function to test for differences in infection rate for the different pathogens in feeding ticks. First, we tested for a difference between feeding nymphs from the combined pool of mustelids and questing nymphs, for which we used data on the infection rate of questing I. ricinus nymphs from three recent studies in Belgium [49,50,51] and one from the Netherlands [52]. We performed this analysis to test for a potential increase in infection rate in feeding nymphs, as an indication that the mustelid species might be able to transmit the different pathogens. Secondly, we tested for differences between feeding nymphs and feeding adults from the combined pool of mustelids, again as an indication that the mustelid species might be able to transmit the different pathogens.


Tick burden

In total, 637 animals were inspected for the presence of ticks. Only half of the animals carried ticks (n = 308), and less than 20% of the animals carried ~80% of the ticks (Table 1). We collected a total of 642 I. ricinus and 2621 I. hexagonus and did not find any other ixodid species, such as Ixodes rugicollis or I. canisuga. In general, pine martens had the highest I. ricinus burden (Table 2 and Additional file 1: Table S1). For larvae, pine martens had a higher burden than badgers (β = 4.7, Z = 3.2, P = 0.006) and stone martens (β = 4.1, Z = 4.7, P < 0.0001). For nymphs, pine martens had a higher burden than stone martens (β = 2.7, Z = 3.5, P = 0.003). For females, stone martens had a higher burden than polecats (β = 1.2, Z = 2.6, P = 0.049). We found a consistent difference between the two countries, where mustelids carried more larvae (β = 4.0, Z = 2.3, P = 0.02), nymphs (β = 5.4, Z = 4.1, P < 0.0001) and females (β = 2.4, Z = 5.3, P < 0.0001) in the Netherlands compared to Belgium.

Table 1 Tick burden of the different mustelid species
Table 2 Tick burden of the different stages of I. ricinus and I. hexagonus. Tick burden is expressed as the mean number of ticks per host individual. Superscript letters indicate significant differences between species (P < 0.05), where estimates increase from a-c, as described in Methods

Ixodes hexagonus showed a different pattern, where polecats generally had the highest burden (Table 2 and Additional file 1: Table S1). For larvae, polecats had a higher burden than badgers (β = 4.5, Z = 4.9, P < 0.0001) and stone martens (β = 2.4, Z = 3.3, P = 0.005). For nymphs, polecats had a higher burden than badgers (β = 1.1, Z = 3.8, P < 0.0001), pine martens (β = 2.3, Z = 4.0, P < 0.0001) and stone martens (β = 2.2, Z = 4.8, P < 0.0001). For females, polecats had a higher burden than badgers (β = 1.1, Z = 3.3, P = 0.004), pine martens (β = 2.9, Z = 4.0, P < 0.001) and stone martens (β = 2.6, Z = 4.5, P < 0.001), and badgers had a higher burden than pine martens (β = 1.8, Z = 2.6, P = 0.04) and stone martens (β = 1.5, Z = 2.6, P = 0.04). Again, we found a consistent difference between the two countries, where mustelids carried more larvae (β = 3.9, Z = 5.2, P < 0.0001), nymphs (β = 2.6, Z = 5.7, P < 0.0001) and females (β = 1.5, Z = 2.8, P = 0.005) in the Netherlands compared to Belgium.

Infection rate in animals

In total, organs of 789 animals were inspected for pathogens. Based on the qPCR analyses of spleen, liver and ear samples, we detected infection with A. phagocytophilum in 27 polecats (n = 556), 11 pine martens (n = 51), two badgers (n = 114) and one stone marten (n = 68) (Table 3). There were no differences in infection rate between the species (Additional file 1: Table S2), but infection rate was higher in the animals from the Netherlands compared to Belgium (odds ratio: 7.2, Z = 2.7, P = 0.006). Based on conventional PCR followed by sequencing on the Anaplasma-positive samples, liver samples from three polecats could be typed further as ecotype I [44]. DNA of B. burgdorferi (s.l.) was found in organs of two stone martens, two pine martens, one polecat and one badger. Spleen, liver and ear samples were found positive for B. burgdorferi (s.l.). Again, there was no difference in infection rate between the species (Additional file 1: Table S2), but the infection rate was higher in the samples from the Netherlands compared to Belgium (odds ratio: 44.9, Z = 2.3, P = 0.02). Despite several attempts, only two B. burgdorferi (s.l.)-positive samples could be typed to the genospecies level. Borrelia afzelii was identified in two organs from two different stone martens (Table 3). Two pine martens, one stone marten and one polecat were positive for N. mikurensis, whereas we did not detect B. miyamotoi DNA in any of the 789 animals tested. DNA of N. mikurensis was detected in spleen, liver and ear samples. There was no difference in infection rate for N. mikurensis between the species (Additional file 1: Table S2) or between the countries (odds ratio: 14.4, Z = 1.1 P = 0.27).

Table 3 Presence of tick-borne pathogens in mustelid tissues. DNA lysates from mustelid liver (L), spleen (S) and ear (E) biopsies were tested by qPCR for tick-borne pathogens. Number of positive animals and infection rates of each pathogen are shown

Infection rate in feeding ticks

In total, 1828 feeding ticks, both I. ricinus (n = 215) and I. hexagonus (n = 1613), were removed from the mustelids and tested for the presence of four tick-borne pathogens (Table 4). From the 15 I. ricinus larvae feeding on pine martens, only one tested positive for B. miyamotoi DNA, and all tested negative for B. burgdorferi (s.l.), N. mikurensis and A. phagocytophilum. The two feeding larvae from a polecat (n = 1) and a stone marten (n = 1) were negative for all four pathogens (Table 4). All four pathogens were found in feeding I. ricinus nymphs and adults, with nymphal infection rates of 7.5% [B. burgdorferi (s.l.)], 2.5% (B. miyamotoi), 23% (N. mikurensis) and 10% (A. phagocytophilum). We compared these to infection rates of questing I. ricinus nymphs (n = 10,628) based on previous studies [41, 42], which were 13.3% for B. burgdorferi (s.l.), 2.5% for B. miyamotoi, 3.7% for N. mikurensis and 4.0% for A. phagocytophilum. The infection rate with N. mikurensis was higher in feeding I. ricinus nymphs collected from mustelids compared to questing nymphs (odds ratio = 7.6, Z = 5.3, P < 0.0001), but not for B. burgdorferi (s.l.) (odds ratio = 0.53, Z = 1.1, P = 0.29), B. miyamotoi (odds ratio = 1.0, Z = 0.0, P = 1.00) or A. phagocytophilum (odds ratio = 2.7, Z = 1.9, P = 0.06). We found adult I. ricinus infection rates of 13% [B. burgdorferi (s.l.)], 3.8% (B. miyamotoi), 6.3% (N. mikurensis) and 21% (A. phagocytophilum). The infection rate was lower in feeding I. ricinus adults compared to nymphs for N. mikurensis (odds ratio = 4.3, Z = 2.9, P = 0.004), but not for the other pathogens (Additional file 1: Table S3).

Table 4 Presence of tick-borne pathogens in ticks feeding on mustelids. DNA lysates from I. ricinus (I. ric) and I. hexagonus (I. hex) feeding on mustelids (n) were tested by qPCR for tick-borne pathogens. Number of pathogen-positive ticks are shown, as well as the number of animals with positive ticks. Infection rate is only shown when more than 10 ticks were tested

None of the 1613 I. hexagonus ticks feeding on mustelids tested positive for B. miyamotoi (Table 4). Two I. hexagonus larvae from polecats tested positive for A. phagocytophilum (1%). One larvae from a stone marten tested positive for N. mikurensis (6%). The three other pathogens were found in feeding I. hexagonus nymphs (n = 836) and adults (n = 533), with nymphal infection rates of 0.7% [B. burgdorferi (s.l.)], 0.6% (N. mikurensis) and 4.3% (A. phagocytophilum). The infection rates in I. hexagonus adults were higher than in I. hexagonus nymphs for B. burgdorferi (s.l.) (2.3%: odds ratio = 3.2, Z = 2.3, P = 0.02), and A. phagocytophilum (9.8%: odds ratio = 2.4, Z = 3.9, P < 0.0001), but not for N. mikurensis (0.8%: odds ratio = 1.3, Z = 0.3, P = 0.73; Additional file 1: Table S3). The infection rates in I. hexagonus nymphs were lower than the infection rates in I. ricinus nymphs for B. burgdorferi (s.l.) (odds ratio = 11.2, Z = 3.3, P = 0.0009) and N. mikurensis (odds ratio = 48.3, Z = 6.6, P < 0.0001), but not for A. phagocytophilum (odds ratio = 2.5, Z = 1.6, P = 0.10; Additional file 1: Table S3). Infection rates in I. hexagonus adults were consistently lower than infection rates in I. ricinus adults (B. burgdorferi (s.l.): odds ratio = 6.7, Z = 5.1, P < 0.0001; A. phagocytophilum: odds ratio = 2.4, Z = 3.7, P = 0.0003; N. mikurensis: odds ratio = 8.9, Z = 3.7, P = 0.0003).

Relative importance of mustelids in maintaining I. ricinus and B. burgdorferi (s.l.)

Because we found few I. ricinus larvae feeding on any of the mustelid species, the realized reservoir competences of these hosts could not be determined. Therefore, only the relative importance of the four host species in terms of feeding the different stages of I. ricinus was determined (Table 5). Modelled I. ricinus burdens were based on a subset of 110 badgers, 52 pine martens, 307 polecats and 69 stone martens for which we could estimate the percentage forest cover around the location where they were found. This percentage did not differ between species or countries (GLM with binomial distribution and logit link, all P > 0.2). In this subset, I. ricinus larval burden increased with forest cover (β = 3.7, Z = 3.1, P = 0.002), while nymphal burden (β = -0.4, Z = -0.4, P = 0.70) and female burden (β = 1.0, Z = 1.5, P = 0.13) did not. Using fitted GLMs to predict I. ricinus burden at 100% forest cover for each of the mustelid species resulted in the estimated tick burdens in Table 5. Due to their relatively low burden and low density, estimated relative importance for all four mustelid species was low, with the highest estimate of 0.5% of the adult I. ricinus fed by pine martens.

Table 5 Relative importance of mustelids for feeding the different stages of I. ricinus and for infecting I. ricinus larvae with B. burgdorferi (s.l.). Vertebrate species considered in our model host-community and the density that was used in these calculations. Mustelid species are shown in bold


Unravelling which host species maintain tick and pathogen populations is important for understanding pathogen dynamics and for developing measures to reduce disease burden [1]. In this study, we investigated the tick burden and infection rate with four tick-borne pathogens of four species of mustelid (badger, pine marten, polecat and stone marten) that generally occur in European forests [24]. Overall, we found relatively low I. ricinus and I. hexagonus burdens on all four mustelid species and low infection rates with three pathogens Anaplasma phagocytophilum, Borrelia burgdorferi (s.l.) and Neoehrlichia mikurensis. Given the low densities in which these four host species generally occur, this leads us to conclude that these common mustelids play a limited role in maintaining I. ricinus populations and in infecting I. ricinus with tick-borne pathogens in forested areas.

The distribution of I. ricinus on mustelids was highly overdispersed, like many other vertebrate species [53], and summarizing this distribution by a mean value might not result in the most accurate parameter [54]. However, we still chose to present the mean (Table 2) to make our findings comparable with other studies. Subsequently, we did take overdispersion into account in our analyses of tick burden. Pine martens had the highest I. ricinus burden for all stages. This could be due to the fact that pine martens were only collected in the Netherlands and mustelids in the Netherlands had a higher I. ricinus burden. However, it is more likely due to the habitat preference of both tick and host. Ixodes ricinus is mainly found in forested areas [18, 55,56,57,58], which is also the habitat that is preferred by the pine marten [59]. This was further supported by a positive correlation of I. ricinus burden with the percentage forest cover around the roadkill site for some tick stages. The three other mustelid species are regularly found in forested areas, but also spend a considerable amount of time in other, more open, habitats [23], which could explain their lower I. ricinus burden.

As only a few previous studies have been performed on the burden of I. ricinus on mustelids, it is difficult to put our results into perspective. We found a higher prevalence of infestation and average I. ricinus burden on pine martens, stone martens and badgers than in previous studies [18, 25, 26, 60]. Two other studies [22, 23] used a similar methodology, using road-killed animals. The risk with using road-killed animals is that ticks might have left the host before collection. We tried to collect animals as soon as possible after they died, but the exact timing of death could not be determined. Therefore, our estimates of tick burden might be underestimations due to loss of ticks post-mortem. However, as Christian [60] found similar tick burdens on live-trapped animals, we think our estimates are representative. Furthermore, several studies on other host species also used road-killed specimens, increasing the comparability between datasets [25, 61].

We found a consistently lower number of I. ricinus of all stages on badgers, polecats and stone martens collected in Belgium compared to the Netherlands. This could be caused by a difference in sampling between the two countries. For example, almost all Belgian (i.e. Flemish) animals checked for ticks were collected before 2006, while most Dutch specimens were collected and screened after 2006 and up to 2015. Therefore, the difference in tick burden might be caused by spatio-temporal variations in abundance of ticks and of their infection rate in the past decades [62]. Another explanation might be a difference between the two countries in habitat where the mustelids were collected. This is, however, unlikely as there was no difference in estimated forest cover around the road-kill site between species or countries. There might have been a difference in forest composition between the countries, which could result in differences in questing tick densities and infection rate [63]. A more elaborate analysis, including detailed information on the landscape surrounding the place where the animals were found, is needed to further investigate this relationship.

All four mustelid species were infested with I. hexagonus of all stages, suggesting that they play an important role in maintaining I. hexagonus populations. Little is known about the population dynamics of I. hexagonus, although previous studies have described the tick species to be parasitizing several members of the Mustelidae family: the common weasel, stoat (Mustela erminea), polecat, badger and otter [27, 29]. Furthermore, hedgehogs (Erinaceus europaeus and Erinaceus roumanicus) have often been found to host I. hexagonus [34, 35]. As I. hexagonus is a nidiculous tick species (searching for a host in the nest or burrow) and as all these vertebrate species use different burrows to nest in, it is possible that I. hexagonus forms many cryptic subpopulations with only occasional genetic exchange between subpopulations, potentially giving rise to host races as was found for I. uriae [64]. Based on our results, it seems possible that polecats, and to some extent all the studied mustelids, can play a role in maintaining I. hexagonus subpopulations.

Borrelia burgdorferi (s.l.) and A. phagocytophilum DNA were detected in tissue samples of European badgers, polecats, pine martens and stone martens (Table 3), suggesting that all four mustelid species can be infected and consequently act as reservoir hosts for these two pathogens. Neoehrlichia mikurensis DNA was present in organs of polecats, stone martens and pine martens, but not in badgers. One explanation might be that badgers are not a reservoir host for N. mikurensis. Recent studies provide evidence that badgers and foxes are reservoir hosts for a related Neoehrlichia species [65]. Borrelia miyamotoi was not detected in any of the 867 investigated organ samples (Table 3), indicating that these animals do not play a major role in the transmission of B. miyamotoi. The higher infection rate in mustelids collected in the Netherlands compared to the animals collected in Belgium could simply be a response to the difference in tick burden we found between the countries as the two parameters are highly correlated [18].

This study assessed, by means of molecular methods, the circulation of four pathogens in their enzootic cycle. The presence of the DNA of a pathogen in internal organs such as liver or spleen is an indication that the host is infected [66, 67]. Although culturing is considered the most reliable method in determining the presence of viable or infectious microorganisms, it suffers from low sensitivity. Furthermore, it is costly, time consuming and difficult, if not impossible for some pathogens [3]. Our investigation only tested for the presence of DNA from four pathogens and not their viability or infectivity. However, previous studies implicate I. ricinus as their vector and the investigated animals as potential, but poorly studied, hosts [35, 68, 69]. Therefore, the inability of a DNA-based detection method to assess infectiousness in host species was expected to be a minor issue. On the other hand, the absence of detectable DNA from these microorganisms in tick lysates or animal tissues does not imply complete absence of these infectious agents either as pathogens might be present in tissues in concentrations lower than the detection limit of the PCR-array. Furthermore, some pathogens might be missed. For example B. miyamotoi, might give rise to short-term, limited infections [70, 71], whereas the tissue tropism for B. afzelii is most likely skin, rather than internal organs [72, 73]. As a result, one should be careful when interpreting results based on only PCR methods as in this study.

Although we detected B. burgdorferi (s.l.) DNA in organs of all mustelid species, we did not find any feeding larvae that were infected with B. burgdorferi (s.l.). Only one study previously investigated the infection rate of B. burgdorferi (s.l.) in organs for any of the species in this study. Gern & Sell [69] found a far higher infection rate compared to our findings. This could be due to the larger number of skin samples tested by these authors, which is the preferred tissue for B. afzelii. The high infection rate in their study could also, however, have been a result of the low sample size (6 animals). Infection rate in feeding I. ricinus nymphs was similar to that in questing nymphs. These results suggested that none of the studied animal species are important reservoirs for B. burgdorferi (s.l.) Unfortunately, the number of feeding larvae that we found on the mustelids was too low to be able to calculate a representative realized reservoir competence [18].

The infection rate of B. burgdorferi (s.l.), A. phagocytophilum and N. mikurensis were higher in feeding I. hexagonus adults compared to feeding I. hexagonus nymphs, a pattern that would be expected when a host species amplifies a pathogen [35, 44]. Furthermore, two I. hexagonus larvae feeding on polecats were infected with A. phagocytophilum, while one I. hexagonus larva feeding on a stone marten was infected with N. mikurensis. However, these pathogens are not transovarially transmitted [74]. Therefore, these findings suggest that the studied mustelids can transmit these pathogens to I. hexagonus ticks. In theory, I. hexagonus can function as a bridging vector for some tick-borne pathogens, maintaining enzootic cycles with occasional spillover to I. ricinus. However, the low densities of mustelids and relatively low I. ricinus burden plus the low number of people that are bitten by I. hexagonus [30, 75] indicate that this enzootic cycle does not pose a public health risk. The absence of B. miyamotoi DNA in the investigated I. hexagonus (n = 1613), and tissues from mustelids (n = 789) indicate that this tick species and these animal species do not play a significant role in the enzootic transmission cycle of B. miyamotoi.

Given the relatively low densities of badgers, pine martens, polecats and stone martens in forest areas (0.44–2.52 km-2), and the low to moderate abundance of I. ricinus feeding on them (larvae: 0.1–47.5 individual-1, nymphs: 0–1 individual-1, adults: 0.1–3.6 individual-1; Table 5), their relative importance in feeding any stage of I. ricinus is low. Species with a high relative importance for maintaining the different stages of I. ricinus all have either high average I. ricinus burdens (roe deer: 25.3 adults individual-1), high density (small rodents: 1000–1200 km-2) or moderate I. ricinus burden and density (blackbird Turdus merula: 4.3 nymphs individual-1 and 200 km-2) [18]. Although we could not quantify the realized reservoir competence for B. burgdorferi for any of the mustelid species, we infer that the relative importance of mustelids in the transmission cycle of B. burgdorferi (s.l.), and other tick-borne pathogens, is low due to their low population density, low I. ricinus burden and low infection rate with pathogens (0.2–3.9%). Hofmeester et al. [18] showed that transmission cycles can exist when hosts reach either high densities (small rodents) or high infection rates (blackbirds: 86%), neither of which is the case for the four species of mustelid. Quantitative (meta) analyses on the relative importance of vertebrate hosts in the transmission cycle of B. miyamotoi, A. phagocytophilum and N. mikurensis have never been performed and are highly needed to understand which animals contribute most to the disease risk of these tick-borne pathogens [1, 63]. As all four mustelid species occur in (peri-)urban areas, they might act as propagation hosts for adult I. ricinus [76, 77], and might be involved in maintaining relatively small populations of I. ricinus in these areas.


Although all stages of I. ricinus were found on mustelids, the relative contribution of mustelids to the life-cycle of this tick species in forested areas is low due to low host densities and tick burdens. Consequently, we conclude that their relative contribution to the enzootic cycle of tick-borne pathogens is low. Moreover, a recent study showed that (meso-)carnivores, such as foxes and stone martens, can lower the number of ticks feeding on rodents, and therefore may have cascading negative effects on tick-borne disease risk [78]. Therefore, the role of mustelids in I. ricinus-borne pathogen dynamics is more likely indirect through species interactions with their prey than direct as hosts for ticks and pathogens. In contrast, all four mustelid species carried all stages of I. hexagonus, potentially maintaining an enzootic cycle of this species, apart from the cycle involving hedgehogs as main host species. Similarly, as badgers, martens and polecats carry adult I. ricinus, they might be able to act as propagation hosts in sites where deer are absent [56]. As such, they may contribute to the maintenance of I. ricinus populations in (peri)-urban areas. Finally, our study shows the importance of filling knowledge gaps in terms of the tick burden and infection rate of little-studied host species that occur in many European forests [1, 18].


Cq :

Quantification cycle


Generalized linear models


5S-23S intergenic spacer region


Lyme borreliosis


tick-borne pathogens


  1. 1.

    Sprong H, Azagi T, Hoornstra D, Nijhof AM, Knorr S, Baarsma ME, et al. Control of Lyme borreliosis and other Ixodes ricinus-borne diseases. Parasit Vectors. 2018;11:145.

    Article  Google Scholar 

  2. 2.

    Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Pena A, George JC, et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit Vectors. 2013;6:1.

    Article  Google Scholar 

  3. 3.

    Tijsse-Klasen E, Koopmans MP, Sprong H. Tick-borne pathogen - reversed and conventional discovery of disease. Front Public Health. 2014;2:73.

    Article  Google Scholar 

  4. 4.

    Leschnik MW, Kirtz GC, Thalhammer JG. Tick-borne encephalitis (TBE) in dogs. Int J Med Microbiol. 2002;291(Suppl. 33):66–9.

    Article  Google Scholar 

  5. 5.

    Zintl A, Mulcahy G, Skerrett HE, Taylor SM, Gray JS. Babesia divergens, a bovine blood parasite of veterinary and zoonotic importance. Clin Microbiol Rev. 2003;16:622–36.

    Article  Google Scholar 

  6. 6.

    Pantchev N, Pluta S, Huisinga E, Nather S, Scheufelen M, Vrhovec MG, et al. Tick-borne diseases (borreliosis, anaplasmosis, babesiosis) in German and Austrian dogs: status quo and review of distribution, transmission, clinical findings, diagnostics and prophylaxis. Parasitol Res. 2015;114(Suppl. 1):S19–54.

    Article  Google Scholar 

  7. 7.

    Gilbert L. Louping ill virus in the UK: a review of the hosts, transmission and ecological consequences of control. Exp Appl Acarol. 2016;68:363–74.

    Article  Google Scholar 

  8. 8.

    Hovius E, de Bruin A, Schouls L, Hovius J, Dekker N, Sprong H. A lifelong study of a pack Rhodesian ridgeback dogs reveals subclinical and clinical tick-borne Anaplasma phagocytophilum infections with possible reinfection or persistence. Parasit Vectors. 2018;11:238.

    Article  Google Scholar 

  9. 9.

    Levi T, Keesing F, Holt RD, Barfield M, Ostfeld RS. Quantifying dilution and amplification in a community of hosts for tick-borne pathogens. Ecol Appl. 2016;26:484–98.

    Article  Google Scholar 

  10. 10.

    Millins C, Gilbert L, Medlock J, Hansford K, Thompson DB, Biek R. Effects of conservation management of landscapes and vertebrate communities on Lyme borreliosis risk in the United Kingdom. Philos Trans R Soc Lond B Biol Sci. 2017;372(1722).

    Article  Google Scholar 

  11. 11.

    Coipan EC, Jahfari S, Fonville M, Oei GA, Spanjaard L, Takumi K, et al. Imbalanced presence of Borrelia burgdorferi s.l. multilocus sequence types in clinical manifestations of Lyme borreliosis. Infect Genet Evol. 2016;42:66–76.

    CAS  Article  Google Scholar 

  12. 12.

    Stanek G, Wormser GP, Gray J, Strle F. Lyme borreliosis. Lancet. 2012;379:461–73.

    Article  Google Scholar 

  13. 13.

    Jahfari S, Krawczyk A, Coipan EC, Fonville M, Hovius JW, Sprong H, et al. Enzootic origins for clinical manifestations of Lyme borreliosis. Infect Genet Evol. 2017;49:48–54.

    Article  Google Scholar 

  14. 14.

    Hamsikova Z, Coipan C, Mahrikova L, Minichova L, Sprong H, Kazimirova M. Borrelia miyamotoi and co-infection with Borrelia afzelii in Ixodes ricinus ticks and rodents from Slovakia. Microb Ecol. 2017;73:1000–8.

    CAS  Article  Google Scholar 

  15. 15.

    Hanincova K, Schafer SM, Etti S, Sewell HS, Taragelova V, Ziak D, et al. Association of Borrelia afzelii with rodents in Europe. Parasitology. 2003;126:11–20.

    CAS  Article  Google Scholar 

  16. 16.

    Kurtenbach K, Dizij A, Seitz HM, Margos G, Moter SE, Kramer MD, et al. Differential immune responses to Borrelia burgdorferi in European wild rodent species influence spirochete transmission to Ixodes ricinus L. (Acari: Ixodidae). Infect Immun. 1994;62:5344–52.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Jaenson TG, Talleklint L, Lundqvist L, Olsen B, Chirico J, Mejlon H. Geographical distribution, host associations, and vector roles of ticks (Acari: Ixodidae, Argasidae) in Sweden. J Med Entomol. 1994;31:240–56.

    CAS  Article  Google Scholar 

  18. 18.

    Hofmeester TR, Coipan EC, van Wieren SE, Prins HH, Takken W, Sprong H. Few vertebrate species dominate the Borrelia burgdorferi s.l. life cycle. Environ Res Lett. 2016;11:043001.

    Article  Google Scholar 

  19. 19.

    Mannelli A, Bertolotti L, Gern L, Gray J. Ecology of Borrelia burgdorferi sensu lato in Europe: transmission dynamics in multi-host systems, influence of molecular processes and effects of climate change. FEMS Microbiol Rev. 2012;36:837–61.

    CAS  Article  Google Scholar 

  20. 20.

    Barbour AG, Bunikis J, Fish D, Hanincova K. Association between body size and reservoir competence of mammals bearing Borrelia burgdorferi at an endemic site in the northeastern United States. Parasit Vectors. 2015;8:299.

    Article  Google Scholar 

  21. 21.

    Huang ZY, de Boer WF, van Langevelde F, Olson V, Blackburn TM, Prins HH. Species’ life-history traits explain interspecific variation in reservoir competence: a possible mechanism underlying the dilution effect. PLoS One. 2013;8:e54341.

    CAS  Article  Google Scholar 

  22. 22.

    Marsot M, Henry PY, Vourc'h G, Gasqui P, Ferquel E, Laignel J, et al. Which forest bird species are the main hosts of the tick, Ixodes ricinus, the vector of Borrelia burgdorferi sensu lato, during the breeding season? Int J Parasitol. 2012;42:781–8.

    CAS  Article  Google Scholar 

  23. 23.

    Stubbe M, Krapp F. Band 5 teil II: Raubsäuger - Carnivora (Fissipedia), Mustelidae, Viverridae, Herpestidae, Felidae, vol. 5. AULA-Verlag: Wiesbaden; 1993.

    Google Scholar 

  24. 24.

    Hofmeester TR, Rowcliffe JM, Jansen PA. Quantifying the availability of vertebrate hosts to ticks: a camera-trapping approach. Front Vet Sci. 2017;4:115.

  25. 25.

    Lorusso V, Lia RP, Dantas-Torres F, Mallia E, Ravagnan S, Capelli G, et al. Ixodid ticks of road-killed wildlife species in southern Italy: new tick-host associations and locality records. Exp Appl Acarol. 2011;55:293–300.

    Article  Google Scholar 

  26. 26.

    Millan J, Ruiz-Fons F, Marquez FJ, Viota M, Lopez-Bao JV, Paz Martin-Mateo M. Ectoparasites of the endangered Iberian lynx Lynx pardinus and sympatric wild and domestic carnivores in Spain. Med Vet Entomol. 2007;21:248–54.

    CAS  Article  Google Scholar 

  27. 27.

    Arthur DR. The host relationships of Ixodes hexagonus leach in Britain. Parasitology. 1953;43:227–38.

    CAS  Article  Google Scholar 

  28. 28.

    Hornok S, Sandor AD, Beck R, Farkas R, Beati L, Kontschan J, et al. Contributions to the phylogeny of Ixodes (Pholeoixodes) canisuga, I. (Ph.) kaiseri, I. (Ph.) hexagonus and a simple pictorial key for the identification of their females. Parasit Vectors. 2017;10:545.

    Article  Google Scholar 

  29. 29.

    Walker MD. The hedgehog tick, Ixodes hexagonus (Leach, 1815) (Acari: Ixodidae); The natural history and ecology of a nest ectoparasite. Syst Appl Acarol. 2018;23:680–714.

    Article  Google Scholar 

  30. 30.

    Estrada-Pena A, Jongejan F. Ticks feeding on humans: a review of records on human-biting Ixodoidea with special reference to pathogen transmission. Exp Appl Acarol. 1999;23:685–715.

    CAS  Article  Google Scholar 

  31. 31.

    Gern L, Rouvinez E, Toutoungi LN, Godfroid E. Transmission cycles of Borrelia burgdorferi sensu lato involving Ixodes ricinus and/or I. hexagonus ticks and the European hedgehog, Erinaceus europaeus, in suburban and urban areas in Switzerland. Folia Parasitol (Praha). 1997;44:309–14.

    CAS  PubMed  Google Scholar 

  32. 32.

    Gern L, Toutoungi LN, Hu CM, Aeschlimann A. Ixodes (Pholeoixodes) hexagonus, an efficient vector of Borrelia burgdorferi in the laboratory. Med Vet Entomol. 1991;5:431–5.

    CAS  Article  Google Scholar 

  33. 33.

    Streissle G. Studies on the transmission of the virus of early summer meningoencephalitis by the tick Ixodes hexagonus Leach. Zentralbl Bakteriol. 1960;179:289–97 (In German).

    CAS  PubMed  Google Scholar 

  34. 34.

    Foldvari G, Jahfari S, Rigo K, Jablonszky M, Szekeres S, Majoros G, et al. Candidatus Neoehrlichia mikurensis and Anaplasma phagocytophilum in urban hedgehogs. Emerg Infect Dis. 2014;20:496–8.

    Article  Google Scholar 

  35. 35.

    Jahfari S, Ruyts SC, Frazer-Mendelewska E, Jaarsma R, Verheyen K, Sprong H. Melting pot of tick-borne zoonoses: the European hedgehog contributes to the maintenance of various tick-borne diseases in natural cycles urban and suburban areas. Parasit Vectors. 2017;10:134.

    Article  Google Scholar 

  36. 36.

    Skuballa J, Petney T, Pfaffle M, Taraschewski H. Molecular detection of Anaplasma phagocytophilum in the European hedgehog (Erinaceus europaeus) and its ticks. Vector Borne Zoonotic Dis. 2010;10:1055–7.

    Article  Google Scholar 

  37. 37.

    LoGiudice K, Ostfeld RS, Schmidt KA, Keesing F. The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. Proc Natl Acad Sci USA. 2003;100:567–71.

    CAS  Article  Google Scholar 

  38. 38.

    Van Den Berge K, Broekhuizen S, Müskens GJDM. Voorkomen van de boommarter Martes martes in Vlaanderen en het zuiden van Nederland. Lutra. 2000;43:125–36.

    Google Scholar 

  39. 39.

    Halos L, Jamal T, Vial L, Maillard R, Suau A, Le Menach A, et al. Determination of an efficient and reliable method for DNA extraction from ticks. Vet Res. 2004;35:709–13.

    CAS  Article  Google Scholar 

  40. 40.

    Heylen D, Tijsse E, Fonville M, Matthysen E, Sprong H. Transmission dynamics of Borrelia burgdorferi s.l. in a bird tick community. Environ Microbiol. 2013;15:663–73.

    Article  Google Scholar 

  41. 41.

    Coipan EC, Fonville M, Tijsse-Klasen E, van der Giessen JW, Takken W, Sprong H, et al. Geodemographic analysis of Borrelia burgdorferi sensu lato using the 5S-23S rDNA spacer region. Infect Genet Evol. 2013;17:216–22.

    CAS  Article  Google Scholar 

  42. 42.

    Zimmermann J, Voss H, Schwager C, Stegemann J, Ansorge W. Automated Sanger dideoxy sequencing reaction protocol. FEBS Lett. 1988;233:432–6.

    CAS  Article  Google Scholar 

  43. 43.

    Hovius JW, de Wever B, Sohne M, Brouwer MC, Coumou J, Wagemakers A, et al. A case of meningoencephalitis by the relapsing fever spirochaete Borrelia miyamotoi in Europe. Lancet. 2013;382:658.

    Article  Google Scholar 

  44. 44.

    Jahfari S, Coipan EC, Fonville M, van Leeuwen AD, Hengeveld P, Heylen D, et al. Circulation of four Anaplasma phagocytophilum ecotypes in Europe. Parasit Vectors. 2014;7:365.

    Article  Google Scholar 

  45. 45.

    Jahfari S, Fonville M, Hengeveld P, Reusken C, Scholte EJ, Takken W, et al. Prevalence of Neoehrlichia mikurensis in ticks and rodents from north-west Europe. Parasit Vectors. 2012;5:74.

    CAS  Article  Google Scholar 

  46. 46.

    Pekkarinen A, Reithmaier L, Strobl P. Pan-European forest/non-forest mapping with Landsat ETM+ and CORINE Land Cover 2000 data. ISPRS J Photogramm Remote Sens. 2009;64:171–83.

    Article  Google Scholar 

  47. 47.

    Jones KE, Bielby J, Cardillo M, Fritz SA, O'Dell J, Orme CDL, et al. PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals. Ecology. 2648;2009:90.

    Google Scholar 

  48. 48.

    Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom J. 2008;50:346–63.

    Article  Google Scholar 

  49. 49.

    Heylen D, Fonville M, van Leeuwen AD, Sprong H. Co-infections and transmission dynamics in a tick-borne bacterium community exposed to songbirds. Environ Microbiol. 2016;18:988–96.

    Article  Google Scholar 

  50. 50.

    Heylen D, Matthysen E, Fonville M, Sprong H. Songbirds as general transmitters but selective amplifiers of Borrelia burgdorferi sensu lato genotypes in Ixodes rinicus ticks. Environ Microbiol. 2014;16:2859–68.

    Article  Google Scholar 

  51. 51.

    Ruyts SC, Tack W, Ampoorter E, Coipan EC, Matthysen E, Heylen D, et al. Year-to-year variation in the density of Ixodes ricinus ticks and the prevalence of the rodent-associated human pathogens Borrelia afzelii and B. miyamotoi in different forest types. Ticks Tick Borne Dis. 2018;9:141–5.

    Article  Google Scholar 

  52. 52.

    Coipan EC, Jahfari S, Fonville M, Maassen CB, van der Giessen J, Takken W, et al. Spatiotemporal dynamics of emerging pathogens in questing Ixodes ricinus. Front Cell Infect Microbiol. 2013;3:36.

    Article  Google Scholar 

  53. 53.

    Randolph SE. Tick ecology: processes and patterns behind the epidemiological risk posed by ixodid ticks as vectors. Parasitology. 2004;129(Suppl.):S37–S65.

    Article  Google Scholar 

  54. 54.

    Rozsa L, Reiczigel J, Majoros G. Quantifying parasites in samples of hosts. J Parasitol. 2000;86:228–32.

    CAS  Article  Google Scholar 

  55. 55.

    Braks M, van Ginkel R, Wint W, Sedda L, Sprong H. Climate change and public health policy: translating the science. Int J Environ Res Public Health. 2013;11:13–29.

    Article  Google Scholar 

  56. 56.

    Hofmeester TR, Sprong H, Jansen PA, Prins HHT, van Wieren SE. Deer presence rather than abundance determines the population density of the sheep tick, Ixodes ricinus, in Dutch forests. Parasit Vectors. 2017;10:433.

    Article  Google Scholar 

  57. 57.

    Sprong H, Hofhuis A, Gassner F, Takken W, Jacobs F, van Vliet AJ, et al. Circumstantial evidence for an increase in the total number and activity of Borrelia-infected Ixodes ricinus in the Netherlands. Parasit Vectors. 2012;5:294.

    Article  Google Scholar 

  58. 58.

    Swart A, Ibanez-Justicia A, Buijs J, van Wieren SE, Hofmeester TR, Sprong H, et al. Predicting tick presence by environmental risk mapping. Front Public Health. 2014;2:238.

    Article  Google Scholar 

  59. 59.

    Brainerd SM, Helldin J-O, Lindström ER, Rolstad E, Rolstad J, Storch I. Pine marten (Martes martes) selection of resting and denning sites in Scandinavian managed forests. Ann Zool Fenn. 1995;32:151–7.

    Google Scholar 

  60. 60.

    Christian A. Zeckenbefall am baummarder in Mecklenburg. Abh Ber Naturkundemus Görlitz. 2002;74:15–20.

    Google Scholar 

  61. 61.

    Ruyts SC, Frazer-Mendelewska E, van Den Berge K, verheyen K, Sprong H. Molecular detection of tick-borne pathogens Borrelia afzelii, Borrelia miyamotoi and Anaplasma phagocytophilum in Eurasian red squirrels (Sciurus vulgaris). Eur J Wildl Res. 2017;63:43.

    Article  Google Scholar 

  62. 62.

    Takken W, van Vliet AJ, Verhulst NO, Jacobs FH, Gassner F, Hartemink N, et al. Acarological risk of Borrelia burgdorferi sensu lato infections across space and time in the Netherlands. Vector Borne Zoonotic Dis. 2017;17:99–107.

    Article  Google Scholar 

  63. 63.

    Ruyts SC, Ampoorter E, Coipan EC, Baeten L, Heylen D, Sprong H, et al. Diversifying forest communities may change Lyme disease risk: extra dimension to the dilution effect in Europe. Parasitology. 2016;143:1310–9.

    CAS  Article  Google Scholar 

  64. 64.

    McCoy KD, Boulinier T, Tirard C, Michalakis Y. Host-dependent genetic structure of parasite populations: differential dispersal of seabird tick host races. Evolution. 2003;57:288–96.

    Article  Google Scholar 

  65. 65.

    Hornok S, Trauttwein K, Takacs N, Hodzic A, Duscher GG, Kontschan J. Molecular analysis of Ixodes rugicollis, Candidatus Neoehrlichia sp. (FU98) and a novel Babesia genotype from a European badger (Meles meles). Ticks Tick Borne Dis. 2017;8:41–4.

    Article  Google Scholar 

  66. 66.

    Persing DH. Diagnostic molecular microbiology. Current challenges and future directions. Diagn Microbiol Infect Dis. 1993;16:159–63.

    CAS  Article  Google Scholar 

  67. 67.

    Vaneechoutte M, Van Eldere J. The possibilities and limitations of nucleic acid amplification technology in diagnostic microbiology. J Med Microbiol. 1997;46:188–94.

    CAS  Article  Google Scholar 

  68. 68.

    Garcia-Perez AL, Oporto B, Espi A, del Cerro A, Barral M, Povedano I, et al. Anaplasmataceae in wild ungulates and carnivores in northern Spain. Ticks Tick Borne Dis. 2016;7:264–9.

    CAS  Article  Google Scholar 

  69. 69.

    Gern L, Sell K. Isolation of Borrelia burgdorferi sensu lato from the skin of the European badger (Meles meles) in Switzerland. Vector Borne Zoonotic Dis. 2009;9:207–8.

    Article  Google Scholar 

  70. 70.

    Wagemakers A, Koetsveld J, Narasimhan S, Wickel M, Deponte K, Bleijlevens B, et al. Variable major proteins as targets for specific antibodies against Borrelia miyamotoi. J Immunol. 2016;196:4185–95.

    CAS  Article  Google Scholar 

  71. 71.

    Wagemakers A, Staarink PJ, Sprong H, Hovius JW. Borrelia miyamotoi: a widespread tick-borne relapsing fever spirochete. Trends Parasitol. 2015;31:260–9.

    Article  Google Scholar 

  72. 72.

    Jahfari S, Hofhuis A, Fonville M, van der Giessen J, van Pelt W, Sprong H. Molecular detection of tick-borne pathogens in humans with tick bites and erythema migrans, in the Netherlands. PLoS Negl Trop Dis. 2016;10:e0005042.

    Article  Google Scholar 

  73. 73.

    Lin YP, Benoit V, Yang X, Martinez-Herranz R, Pal U, Leong JM. Strain-specific variation of the decorin-binding adhesin DbpA influences the tissue tropism of the lyme disease spirochete. PLoS Pathog. 2014;10:e1004238.

    Article  Google Scholar 

  74. 74.

    Silaghi C, Beck R, Oteo JA, Pfeffer M, Sprong H. Neoehrlichiosis: an emerging tick-borne zoonosis caused by Candidatus Neoehrlichia mikurensis. Exp Appl Acarol. 2016;68:279–97.

    Article  Google Scholar 

  75. 75.

    Fonville M, Friesema IH, Hengeveld PD, Docters van Leeuwen A, Jahfari S, Harms MG, et al. Human exposure to tickborne relapsing fever spirochete Borrelia miyamotoi, the Netherlands. Emerg Infect Dis. 2014;20:1244–5.

    Article  Google Scholar 

  76. 76.

    Tomassone L, Berriatua E, De Sousa R, Duscher GG, Mihalca AD, Silaghi C, et al. Neglected vector-borne zoonoses in Europe: into the wild. Vet Parasitol. 2018;251:17–26.

    Article  Google Scholar 

  77. 77.

    Szekeres S, Majláthová V, Majláth I, Földvári G. Neglected hosts: the role of lacertid lizards and medium-sized mammals in the eco-epidemiology of Lyme borreliosis. In: Braks M, van Wieren S, Takken W, Sprong H, editors. Ecology and Prevention of Lyme Borreliosis, vol. 4. Wageningen: Wageningen Academic Publishers; 2016. p. 103–26.

    Chapter  Google Scholar 

  78. 78.

    Hofmeester TR, Jansen PA, Wijnen HJ, Coipan EC, Fonville M, Prins HHT, et al. Cascading effects of predator activity on tick-borne disease risk. Proc Biol Sci. 2017;284(1859).

    Article  Google Scholar 

Download references


We are grateful to Hugh Jansman, Sim Broekhuizen, Gerard Müskens and the rest of the staff and volunteers at Wageningen Environmental Research; the staff and volunteers of the Dutch Mammal Society; the staff of the DWHC; as well as to Han ten Seldam, Marije Lageschaar, Filip Berlengee, Dirk Vansevenan, students at INBO and volunteers of the Flemish “marternetwerk” for help in collecting and dissecting animals. We would like to thank the two anonymous reviewers for their valuable and constructive comments on a previous version of this manuscript.


This study was financially supported by the Dutch Ministry of Health, Welfare and Sport (VWS). The funders had no role in study design, data collection, interpretation and analysis, decision to publish or preparation of the manuscript.

Availability of data and materials

The data supporting the conclusions of this article are included within the article and its additional files. The datasets used and/or analyzed during the present study are available from the corresponding author upon reasonable request.

Author information




TH, KB, SR, KH and HS conceived the study; TH, KB, MM and JG collected samples; AK, ADL, MF and HS performed DNA extractions and molecular analyses; TH performed statistical analyses; HS, AK, KB and TH wrote the parts of the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Tim R. Hofmeester or Hein Sprong.

Ethics declarations

Ethics approval

The study was carried out according to national animal welfare regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Table S1. Estimated differences from outcomes of models. Model outcomes of generalized linear models with a negative binomial distribution and log link function testing for differences in tick burden between the mustelid species while correcting for differences between countries. One model per tick species and life stage combination. Estimated differences in tick burden (on a log scale) and standard error (between brackets) are given for each combination of species. The species in the first column is the “base” value, and P-values, based on a Tukey post-hoc test, are represented by · for P < 0.1, * for P < 0.05, ** for P < 0.01 and *** for P < 0.001. Table S2. Estimated odds ratio from outcomes of models. Model outcomes of generalized linear models with a binomial distribution and logit link function testing for differences in infection prevalence between the mustelid species while correcting for differences between countries. One model per microorganism. Estimated ln(odds ratio) and standard error are given for each combination of species. The species in the first column is the “base” value. P-values based on a Tukey post-hoc test were all > 0.1 and are not presented in the table. Table S3. Estimated odds ratio from outcomes models. Model outcomes of generalized linear models with a binomial distribution and logit link function testing for differences in infection prevalence between feeding ticks of different stages found on all mustelids. Each cell gives the outcome of a single model. Estimated ln(odds ratio) and standard error are given for each combination. P-values are represented by · for P < 0.1, * for P < 0.05, ** for P < 0.01 and *** for P < 0.001. Table S4. DNA sequences of B. afzelii (IGS) and A. phagocytophilum (GroEL) from tissue samples. (DOCX 28 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hofmeester, T.R., Krawczyk, A.I., van Leeuwen, A.D. et al. Role of mustelids in the life-cycle of ixodid ticks and transmission cycles of four tick-borne pathogens. Parasites Vectors 11, 600 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Meles meles
  • Mustela putorius
  • Martes martes
  • Martes foina
  • Ixodes ricinus
  • Ixodes hexagonus
  • Anaplasma phagocytophilum
  • Borrelia burgdorferi (s.l.)
  • Neoehrlichia mikurensis
  • Borrelia miyamotoi