Skip to main content

Effect of rodent density on tick and tick-borne pathogen populations: consequences for infectious disease risk

Abstract

Background

Rodents are considered to contribute strongly to the risk of tick-borne diseases by feeding Ixodes ricinus larvae and by acting as amplifying hosts for pathogens. Here, we tested to what extent these two processes depend on rodent density, and for which pathogen species rodents synergistically contribute to the local disease risk, i.e. the density of infected nymphs (DIN).

Methods

In a natural woodland, we manipulated rodent densities in plots of 2500 m2 by either supplementing a critical food source (acorns) or by removing rodents during two years. Untreated plots were used as controls. Collected nymphs and rodent ear biopsies were tested for the presence of seven tick-borne microorganisms. Linear models were used to capture associations between rodents, nymphs, and pathogens.

Results

Investigation of data from all plots, irrespective of the treatment, revealed a strong positive association between rodent density and nymphal density, nymphal infection prevalence (NIP) with Borrelia afzelii and Neoehrlichia mikurensis, and hence DIN’s of these pathogens in the following year. The NIP, but not the DIN, of the bird-associated Borrelia garinii, decreased with increasing rodent density. The NIPs of Borrelia miyamotoi and Rickettsia helvetica were independent of rodent density, and increasing rodent density moderately increased the DINs. In addition, NIPs of Babesia microti and Spiroplasma ixodetis decreased with increasing rodent density, which had a non-linear association with DINs of these microorganisms.

Conclusions

A positive density dependence for all rodent- and tick-associated tick-borne pathogens was found, despite the observation that some of them decreased in prevalence. The effects on the DINs were variable among microorganisms, more than likely due to contrasts in their biology (including transmission modes, host specificity and transmission efficiency). The strongest associations were found in rodent-associated pathogens that most heavily rely on horizontal transmission. Our results draw attention to the importance of considering transmission mode of a pathogen while developing preventative measures to successfully reduce the burden of disease.

Background

Lyme borreliosis is the most prevalent tick-borne disease in the northern hemisphere with increasing incidence and expanding endemic regions [1, 2]. The risk of acquiring Lyme borreliosis is partially determined by the density of questing ticks infected with its causative agent, Borrelia burgdorferi (sensu lato) [3, 4]. Particularly, the density of infected nymphs (DIN) is of interest, because humans are predominantly exposed to and infected with Lyme spirochetes, as well as other pathogens, by nymphs [5]. The density of infected questing ticks is a product of the density of questing ticks and infection prevalence of a pathogen, which both express high temporal variations, presumably attributed to changes in weather conditions and fluctuations in the abundance of vertebrate hosts [6,7,8]. The mechanisms underlying these variations are complex, as climatic conditions, vertebrate hosts and their food source, ticks, and tick-borne microorganisms form biological networks with multiple direct and indirect interactions [9]. Therefore, quantifying these interactions will help us to understand changes in the distribution and incidence of Lyme borreliosis and other tick-borne diseases.

The most common vectors of tick-borne diseases in the northern hemisphere are ticks of the Ixodes ricinus complex. Their survival primarily depends on their ability to find a vertebrate host, which may vary between life stages. In forested areas, larvae of the I. ricinus complex feed predominantly on rodents, nymphs on rodents and birds, and adults on ungulates, mostly deer [10, 11]. Although presence of deer is generally responsible for high abundance of ticks [12], variations in the density of nymphs (DON) has been associated with the density of rodents [13]. For instance, the density of host-seeking I. scapularis nymphs was correlated with abundance of white-footed mice in the previous year. White-footed mice are the main hosts for larval I. scapularis; high abundance of these mice provides more opportunities for larvae to feed successfully, and subsequently emerge as nymphs in the following year.

The abundance of rodent species is affected by many different factors, such as predation, vegetation cover, and food availability [10, 14, 15]. A key food supply for rodents is acorns and its seasonal availability has been shown to be responsible for the fluctuations in rodent densities between years and geographical locations [16,17,18]. In general, acorn availability increases the length of the breeding season and facilitates winter survival of forest rodents resulting in a higher rodent density in the following spring [9, 19,20,21,22]. As a consequence, in the temperate zone, an increased rodent density has been shown to cause upsurges in Puumala hantavirus disease in humans [23, 24]. In addition, several North American studies have suggested that acorns and rodents are good predictors for Lyme-disease risk because rodents are reservoir hosts of B. burgdorferi (s.l.) [9, 15, 25]. The causal relationship between rodent fluctuations and Lyme disease incidence, however, remains unresolved as this has not been investigated in experimental settings, enabling the exclusion of confounding factors.

In the Netherlands, the wood mice (Apodemus sylvaticus Linnaeus) and bank voles (Myodes glareolus (Schreber)) are amplifying hosts of several tick-borne pathogens including B. afzelii, B. miyamotoi, Babesia microti and Neoehrlichia mikurensis [26,27,28,29,30], and the most common hosts of larval I. ricinus [31]. Apart from rodent-borne pathogens mentioned above, I. ricinus carries many other microorganisms including B. garinii, Spiroplasma ixodetis and Rickettsia helvetica [32]. Most, if not all, of the pathogens are transmitted between ticks via a vertebrate host (horizontally), which can be broadly divided into co-feeding and systemic transmission (Table 1). Co-feeding relies on localized and temporal infection in the vertebrate skin and occurs when infected and uninfected ticks feed close to each other [33, 34]. Systemic transmission depends more on a persistent infection in a host, which can be local (e.g. skin) or systemic (e.g. blood) [34]. Amplifying hosts are responsible for producing infected ticks and therefore, for increasing risk of human exposure. In addition, ticks maintain microorganisms such as S. ixodetis via vertical transmission, with different efficiency (Table 1). Some bacteria such as R. helvetica and B. miyamotoi, can utilize both horizontal and vertical transmission routes [35, 36]. It is unclear how variations in rodent densities affect the disease risk of tick-borne pathogens with different transmission modes, particularly in the European setting.

Table 1 Transmission modes and amplification hosts of tick-borne microorganisms

The goal of the present study was to investigate how rodent densities, the density of I. ricinus nymphs and transmission dynamics of tick-borne pathogens interact in order to generate the density of infected ticks. To our knowledge, this is the first European study experimentally investigating these relationships in the field. In addition, no prior study has assessed the influence of rodent density on prevalence and density of tick-borne microorganisms other than rodent-borne. Our approach was to artificially manipulate the rodent densities by either acorn addition or rodent removal for two consecutive years in a natural habitat. We measured and quantified the rodent, nymph, and pathogen population responses to these treatments, as well as performed regression analysis. Using this approach, we aimed to learn whether rodent densities play a major role in shaping the density of questing ticks and transmission dynamics of tick-borne microorganisms, which in turn, will help assess and potentially predict disease risk and formulate possible intervention strategies.

Given that rodents are locally the most substantial hosts for larvae [31] and high rodent density results in high larval encounter rates, increase of the rodent density at a given yeart is expected to lead to a rise in density of nymphs in the following year (DONt+1). Along with the higher rodent densities, transmission events of tick-borne microorganisms are expected to increase. We anticipate that differences in the microorganisms’ modes of transmission as well as host amplification potential are main determinants in the change after manipulation. Our hypothesis is that the NIPt+1 (nymphal infection prevalence) of tick-borne pathogens, such as B. afzelii, N. mikurensis and B. microti, which are amplified by rodents, is dependent on the density of rodents. Consequently, we expect a synergistic effect of rodent densities on the density of infected nymphs one year later (DINt+1). Also, we hypothesise that rodent densities will not alter the NIPt+1 of tick-associated microorganisms, such as R. helvetica, B. miyamotoi and S. ixodetis, which predominantly rely on vertical transmission. Further, we expect that DINt+1 R. helvetica, DINt+1 B. miyamotoi and DINt+1 S. ixodetis will be only moderately affected by increasing rodent density. In case of B. garinii, a tick-borne pathogen amplified by birds [37, 38], we expect that increasing rodent density will increase the proportion of larvae feeding on rodents and, therefore, have a negative effect on NIPt+1 B. garinii. Lastly, we hypothesize that a higher rodent density will have no effect on DINt+1 B. garinii.

Methods

Study sites

The study was conducted at the forest reserves Planken Wambuis (52°01′45″N, 5°48′49″E) and Noord Ginkel (52°02′23″N, 5°45′09″E) near Wageningen, The Netherlands. Both forests are dominated by Scots pine (Pinus sylvestris) and harbour a diversity of bird and mammal species, including wood mice (A. sylvaticus), bank voles (M. glareolus), common shrews (Sorex araneus), wild boar (Sus scrofa), roe deer (Capreolus capreolus), red deer (Cervus elaphus) and a few free-ranging cattle and horses.

Manipulation of rodent density and estimation of nymphal density

In both forests, six plots of 50 × 50 m were selected with at least 350 m between plots (Additional file 1: Figure S1). Each plot was assigned to one of three treatments (rodent removal, control or acorn addition). In rodent removal plots, rodents were trapped for one night a month with Heslinga live traps (Heslinga Traps, Groningen, The Netherlands) in a 5 × 5 grid with 10 m inter-trap distance. Captured rodents were euthanized by cervical dislocation. The first rodent removal event was directly after the mark-recapture trapping in September 2012. Thereafter, rodents that accidentally found their way into the plots were removed monthly until December 2014 using the same grid with traps. Four control plots received no treatment. To increase rodent density, acorns were added to four plots [39]. Acorns were provided beneath feeding stations, which were made of 60 × 60 cm plates kept 5 cm above the ground to prevent acorn predation by birds and large mammals. In each plot, 16 of these feeding stations were placed in a 4 × 4 grid with 15 m between feeding stations. Each feeding station was provided with 6.25 kg of acorns in November and January of 2012 and 2013 (1600 kg in total). Control feeding stations without acorns were also placed in the control and rodent removal plots. A plastic screen, 40 cm high and dug 10 cm into the ground was placed as a barrier around the four rodent removal plots to prevent immigration of rodents (Additional file 1: Figure S1). To overcome a possible bias in large vertebrate community caused by a visual effect, screens were also placed around the control and acorn addition plots. However, the lowest 10 cm of these screens was left open to enable rodents to walk in and out freely.

Tick density was estimated monthly in each plot by blanket dragging over the vegetation. At each plot, a 1 m2 blanket was dragged over four transits of 50 m and inspected at 25 m intervals. All attached nymphs were counted. Dragging was performed in the afternoons (12:00–18:00 h CET) when the vegetation was dry. Given that nymphs have been shown to quest when the weekly mean daily maximum temperature exceeds 7 °C [40,41,42], we included temperature data from September 2012 to December 2015 to investigate the relationship between temperature and onset of tick activity. Daily measurements were collected from the nearest weather station (Deelen, KNMI, the Netherlands; Additional file 4: Table S1).

Rodent samples and nymph collection

Rodents were sampled at three-month intervals (March, June, September and December) from September 2012 until December 2014. At each plot, 25 Heslinga live traps were placed in a 5 × 5 grid. Traps were pre-baited with oats for 3 days, after which they were rebaited with grain, carrot and mealworms and set at 9:00 h CET. Traps were then inspected four times at 12-h intervals. Trapped rodents were marked by shaving a patch of fur from their side [43]. Rodent density was calculated per species according to the Schnabel method (multiple marking; [44, 45]. During the morning trappings, newly captured rodents were screened for ticks, and larvae were counted. A small ear biopsy was taken with sterile scissors from each newly captured rodent and stored in 70% ethanol at − 20 °C until further analysis.

Questing nymphs were collected during monthly density estimation. All nymphs attached to the blanket were collected and stored individually in 70% ethanol at − 20 °C until further analysis.

DNA extraction and pathogen detection

Ear biopsies and nymphs were analysed individually. DNA from a maximum of 40 nymphs per plot per month was extracted with ammonium hydroxide as described previously [46]. DNA from the ear biopsies was extracted using the Qiagen DNeasy Blood & Tissue Kit according to the manufacturer’s protocol (Qiagen, Venlo, The Netherlands). The lysates were stored at 4 °C. Samples were analysed with different (multiplex) real-time PCRs, based on various target genes depending on microorganism of interest such B. burgdorferi (s.l.) [47], B. miyamotoi [48], N. mikurensis [49], R. helvetica [50], B. microti and S. ixodetis (this study, Additional file 4: Text S1). A detailed description of the qPCR protocol is provided in Additional file 4: Text S1. Samples positive for B. burgdorferi (s.l.), were subjected to conventional PCR followed by sequencing to identify a genotype [47].

Data analysis and modelling

Data analysis and model building were performed in R version 3.5.1 [51] and RStudio [52]. To evaluate whether rodent removal and acorn addition treatments were successful, we compared means of densities of rodents (data from 2013 and 2014) between the treatments using the non-parametric Wilcoxon signed-rank test. The same test was used to evaluate whether the treatments influenced the DON, and density of rodent-associated pathogens, B. afzelii and N. mikurensis (data from 2014 and 2015). Since monthly at each plot 200 m2 were inspected for questing nymphs, we combined these measurements into a yearly DON per 2400 m2 (by summing up all nymphs from 12 months). The differences in prevalence of microorganisms and tick burdens between two rodent species were compared with the Chi-square test and non-parametric Wilcoxon signed-rank test, respectively.

To investigate how well the density of rodents from 2013 and 2014 predicts DONt+1, NIPt+1 and DINt+1, we performed regression analyses. Several linear models for DONt+1 (the annual median) were assessed with different interactions between rodent density, year, and treatment. For NIPt+1, binomial generalized linear models were assessed with different interactions between rodent density and year. Because NIP is represented by fraction data, we choose a binomial generalized linear model taking into account sample size with the logit link transform. For DINt+1, linear models were assessed with different interactions between rodent density and year. DIN data were calculated by multiplying DON and NIP, which are both potentially influenced by rodent density and therefore we have also included (rodent density)2 as a covariate.

Year (2013, 2014) and treatment (acorn addition, control and rodent removal) were categorical variables while DON, NIP and DIN were numerical variables. The ranges of DON, NIP and rodent density are provided in Additional file 4: Table S2 and Table S3, respectively. For all models, best-fitting models were compared on the basis of a likelihood ratio test, R2 (linear models) and AIC (linear and generalized models). Model selection was performed using histograms to visually evaluate normality of the residuals. If there clearly was no best model, the simpler model was selected.

Results

Effect of treatment on rodent density, DON, DINB. afzelii, and DINN. mikurensis

Rodent density was affected by treatment (Fig. 1). With bank voles the effects were apparent throughout the intervention period, while with wood mice addition of acorns led to a strong increase in density in the second year of the study. The removal of rodents led to a lower (P = 0.0031) rodent density and the addition of acorns led to a higher (P = 0.042) rodent density than in the control plots in years 2013 and 2014 (Fig. 1).

Fig. 1
figure 1

a Mean density of two rodent species, bank vole and wood mouse per plot. Solid arrows indicate events of acorn supplementation (November and January); dashed arrows indicate when monthly removal of rodents started (September 2012) and ended (December 2014). b Box plots of rodent density per plot for each treatment (data from 2013 and 2014). The lower and upper hinges correspond to the first and third quartiles (the 25th and 75th percentiles). The upper whisker shows the largest value no further than 1.5 * IQR from the hinge (where IQR is the inter-quartile range, or distance between the first and third quartiles) and the lower whisker shows the smallest value at most 1.5 * IQR of the hinge. The differences in the rodent density between the treatments was calculated based on the mean (black dot) with the Wilcoxon test and the overall difference is statistically significant (P < 0.0001). The diagram shows also the median observation (solid horizontal line)

The DON fluctuated over the years and was the highest from May until October (Fig. 2). We observed that the moment that nymphs started to quest was in the first month of the year with a mean temperature above 7 °C. The number of months with a mean temperature below 7 °C varied between the years (Fig. 2, Additional file 2: Figure S2). In 2013, five months had mean temperatures below 7 °C, whereas both 2014 and 2015 had three months with mean temperatures below 7 °C, but these were spread differently throughout the year. The mean DON of all plots in 2013, 2014 and 2015 were 581, 272 and 257 per 2400 m2 (200 m2 × 12 months), respectively. Mean nymphal density in 2014 and 2015 was significantly lower than in 2013 (P = 0.0083 and P = 0.013, respectively), whereas the mean nymphal densities of 2014 and 2015 were not significantly different (P = 0.63; not shown). There was no effect (P = 0.27) of acorn addition and a negative effect (P = 0.043) of rodent removal on the DON in the same (not shown) or following years. Nevertheless, when the density of nymphs from 2013 served as a baseline to measure the effect of a treatment on the DON in 2014 and 2015, there was no significant effect (Fig. 2).

Fig. 2
figure 2

a Mean density of questing nymphs (DON) per 200 m2. Horizontal solid lines just above the x-axis depict months with average temperature below 7 °C. In winter 2012/2013, the number of months with the mean temperature below 7 °C was five, while in both 2013/2014 and 2014/2015 was four, however different months. b Density of nymphs (DON) in 2014 and 2015 in all three treatments in comparison to 2013 (baseline year). c Differences in DON between the treatments in two separate years calculated with the Wilcoxon test with a correction for a baseline year (2013). The overall differences between the treatments were not significant either in 2014, or 2015 (P > 0.59 and P > 0.87, respectively)

To investigate the effect of treatment on the dynamics of tick-borne pathogens amplified by rodents, we compared the mean DINB. afzelii and DINN. mikurensis between the treatments in 2014 and 2015. Our analyses before and after a correction for a baseline DIN from 2013 showed that there was no effect of either acorn addition or rodent removal on the density of nymphs infected with B. afzelii and N. mikurensis in the following years (Fig. 3 and Additional file 3: Figure S3).

Fig. 3
figure 3

a Density of nymphs infected with B. afzelii (DIN B. afzelii) in 2014 and 2015 in all three treatments in comparison to 2013 (baseline year). b Differences in DIN B. afzelii between the treatments in two separate years calculated with the Wilcoxon test with a correction for a baseline year (2013). The overall differences between the treatments were not significant either in 2014, or 2015 (P = 0.69 and P = 0.53, respectively)

Rodent sample and nymph collection

A total of 2386 rodents was caught in the experiment. From these, 345 bank voles and 547 wood mice were inspected for ticks, from which 155 and 346 were infested with larvae, respectively. The average number of larvae found on wood mice (9.0; 95% CI: 7.6–10.4) was significantly higher (W = 118,520, P < 0.0001) than the average in bank voles (4.2; 95% CI: 3.0–5.4). None of the bank voles and 97 wood mice were infested with nymphs and the average nymphal burden was 0.2 (95% CI: − 0.2–0.6).

A total of 772 ear biopsies was taken (478 from wood mice and 294 from bank voles) and subjected to pathogen detection. In addition, 13,916 nymphs were collected by dragging, from which 7609 were tested for the presence of tick-borne pathogens. A detailed overview of rodent densities, number of analysed rodents, tick density and analysed ticks per treatment, month, and year are provided in Additional file 4: Table S2 and Table S3.

Pathogen detection

In the rodent ear biopsies and the collected questing nymphs we detected DNA of B. burgdorferi (s.l.), B. miyamotoi, N. mikurensis, B. microti, R. helvetica and S. ixodetis (Fig. 4). The sequencing success of qPCR-positive ticks (n = 1017) for B. burgdorferi (s.l.) was 64%, and four genospecies were identified: B. afzelii, B. garinii, B. valaisiana and B. burgdorferi (s.s.). Borrelia-positive rodent biopsies were not sequenced and were treated as B. afzelii in further analysis. A justification for this assumption derives from previous studies, which have shown that, in the Netherlands, more than 99% of the positive rodents infected with B. burgdorferi (s.l.) carried B. afzelii [53, 54]. The prevalence of B. afzelii as well as N. mikurensis was higher in bank voles than in wood mice (χ2 =  3.296, df = 1, P = 0.0694 and χ2 = 4.234, df = 1, P = 0.0396, respectively). Interestingly, S. ixodetis was almost exclusively detected in wood mice with prevalence significantly higher than in bank voles (χ2 = 14.264, df = 1, P = 0.0002), whereas B. microti was almost exclusively found in bank voles with prevalence significantly higher than in wood mice (χ2 = 27.012, df = 1, P < 0.0001). The prevalence of R. helvetica was not significantly different between two rodent species (χ2 = 0.803, df = 1, P = 0.3703). A complete overview of infection prevalence of all pathogens in ticks and rodent biopsies is provided in Additional file 4: Table S4.

Fig. 4
figure 4

Overview of tick-borne microorganism infections in rodents and nymphs. a Rodent infection prevalence separately for each collection month and rodent species. b Density of infected nymphs (DIN) separately for each collection month (data combined from 2013 and 2014)

The only pathogen consistently present in both rodent species throughout the year was B. afzelii (Fig. 4). The infection in rodents persisted despite infected nymphs not being active in months below 7 °C. Other pathogens, such as N. mikurensis, B. miyamotoi, R. helvetica and S. ixodetis, were present in the rodent population mostly when activity of rodents and (infected) nymphs overlapped (Fig. 4).

Rodent density versus DONt+1

Equations of all tested models investigating the association between rodent density and DONt+1, NIPt+1, and DINt+1, their R2, AIC values, and results of a likelihood test are provided in Additional file 4: Table S5. Full equations of the best-fitting models are provided in Additional file 4: Table S6, while Table 2 shows significant interactions incorporated in the models as well as the type of effect rodent density had on all responses (DONt+1, NIPt+1 and DINt+1). Because our treatments affected bank voles and wood mice simultaneously, rodent density data used in the models consist of rodent species added together.

Table 2 Best models for prediction of density of nymphs (DON), nymphal infection prevalence (NIP), and density of infected nymphs (DIN)

The model that fit the data best indicated that rodent density and DONt+1 were significantly positively associated (P = 0.000631). The best model was a linear model of rodent density incorporating year and treatment as covariates explaining 61% of the variance (Table 2, Eq. 1; Fig. 5).

Fig. 5
figure 5

Effect of rodent density on DONt+1. The plot shows the relationships between the number of rodents per plot in year t and DON (number per 200 m2 per plot) in the following year (t+1). Rodent density had significant positive effect on DON in all treatments and years

Rodent density versus rodent-associated pathogens

Regarding B. afzelii and N. mikurensis, there was a significant positive association between rodent density and NIPt+1 (P < 0.0001 and P < 0.0001), and rodent density and DINt+1 (P = 0.000187 and P < 0.0001; Fig. 6). The best model for both NIPt+1 B. afzelii and NIPt+1 N. mikurensis was a simple generalized linear model of rodent density (Table 2, Eq. 2 and Eq. 4). In the case of DINt+1, a simple linear model of rodent density was the best and explained 45% and 56% of the variance in DINt+1 B. afzelii and DINt+1 N. mikurensis, respectively (Table 2, Eq. 3 and Eq. 5). Regarding another pathogen amplified by rodents, B. microti, there was a negative effect (P < 0.0001) of rodent density on NIPt+1 and the best model was a simple generalized linear model of rodent density (Table 2, Eq. 8; Fig. 7). In the case of DINt+1 B. microti, the best model was a linear model of rodent density and (rodent density)2, Table 2, Eq. 9), which explained 20% of the variance. The model including a quadratic term allowed to reveal significant negative (P = 0.0141) non-linear association between rodent density and B. microti (Fig. 7).

Fig. 6
figure 6

Association between density of rodents and pathogens amplified by rodents. The graphs show the relationship between the number of rodents per plot in year t and NIP and DIN (number per 200 m2 per plot) in year t+1. a Effect of rodent density on NIPt+1 B. afzelii. Rodent density had significant positive effect on NIP. b Effect of rodent density on DINt+1 B. afzelii. Rodent density has significant positive effects on DIN. c Effect of rodent density on NIPt+1 N. mikurensis. Rodent density had significant positive effect on NIP. d Effect of rodent density on DONt+1 N. mikurensis. Rodent density had significant positive effect on DIN

Fig. 7
figure 7

Association between density of rodents and tick-associated microorganisms. The graphs show the relationship between the number of rodents per plot in year t and NIP and DIN (number per 200 m2 per plot) in year t+1. a Effect of rodent density on NIPt+1 S. ixodetis. Rodent density had significant negative effect on NIP. b Effect of rodent density on DONt+1 S. ixodetis. Rodent density had significant non-linear effect on DIN. c Effect of rodent density on NIPt+1 B. microti. Rodent density had significant negative effect on NIP. d Effect of rodent density on DINt+1 B. microti. Rodent density had significant non-linear effect on DIN

Rodent density versus a bird-associated pathogen

There was a significant negative association (P = 0.000149) between rodent density and NIPt+1 B. garinii and no association between rodent density and DINt+1 B. garinii (Fig. 8), which remained constant through the experiment. The best model for NIPt+1 B. garinii was a generalized linear model of rodent density and year (Table 2, Eq. 10), while none of the tested models for DINt+1 B. garinii was better than a null model (Table 2, Eq. 11).

Fig. 8
figure 8

Association between density of rodents and a pathogen amplified by birds. The graphs show the relationship between the number of rodents per plot in year t and NIP and DIN (number per 200 m2 per plot) in year t+1. a Effect of rodent density on NIPt+1 B. garinii. Rodent density had significant negative effect on NIP in both years. b Effect of rodent density on DINt+1 B. garinii. Rodent density had no effect on DIN

Rodent density versus vertically transmitted microorganisms

Rodent density had differential effect on NIPt+1 R. helvetica, and a significantly positive effect on DINt+1 R. helvetica (P = 0.0143; Fig. 9). In case of NIPt+1 R. helvetica, the best-fitting model was a generalized linear model taking into account the differences in association with respect to year (Table 2, Eq. 12). For DINt+1 R. helvetica, the best model was a simple linear model of rodent density explaining 21% of the variance (Table 2, Eq. 13).

Fig. 9
figure 9

Association between density of rodents and vertically transmitted pathogens. The graphs show the relationship between the number of rodents per plot in year t and NIP and DIN (number per 200 m2 per plot) in year t+1. a Effect of rodent density on NIPt+1 R. helvetica. Rodent density had inconsistent effect on NIP (no effect in 2013 and significant positive effect in 2014). b Effect of rodent density on DONt+1 R. helvetica. Rodent density had significant positive effect on DIN. c Effects of rodent density on NIPt+1 B. miyamotoi. Rodent density had inconsistent effect on NIP (negative but no significant effect in 2013 and significant positive effect in 2014). d Effects of rodent density on DONt+1 B. miyamotoi. Rodent density had significant positive effect on DIN

Rodent density had a differential association with NIPt+1 B. miyamotoi between the years (Fig. 9). In 2013, the association was negative but not significant (P = 0.15797) and in 2014, positive and significant (P = 0.00862). The association between rodent density and DINt+1 B. miyamotoi was significantly positive (P = 0.0119; Fig. 9). The best model for NIPt+1 B. miyamotoi was a generalized linear model of rodent density taking into account the differences in association with respect to year (Table 2, Eq. 6), and for DINt+1 B. miyamotoi, a simple linear model of rodent density explaining only 22% of the variance (Table 2, Eq. 7).

The association between rodent density and NIPt+1 S. ixodetis was significantly negative (P < 0.0001) and the best model was a simple generalized linear model of rodent density (Table 2, Eq. 14, Fig. 7). In case of DINt+1 S. ixodetis, the best model was a linear model of rodent density and (rodent density)2, which explained 45% of the variance (Table 2, Eq. 15). The model including a quadratic term allowed revealing significant negative (P = 0.005297) non-linear association between rodent density and S. ixodetis (Fig. 7).

Discussion

This study was designed to investigate the association between rodent density and I. ricinus nymphs and tick-borne microorganisms. We observed that the densities of rodents affect DON, NIP and DIN in the following year. We found positive associations between rodent density and DONt+1 regardless of the year and type of treatment (Fig. 5). The NIPt+1 and DINt+1 depending on tick-borne pathogens and microorganisms were associated with the rodent density to a different extent, determined by the infection dynamics of the microorganism species (Figs. 6, 7, 9). In addition, although the treatments affected rodent density in the following years, we did not observe an effect on either the DON (Fig. 2) or DINt+1 B. afzelii and DINt+1 N. mikurensis (Fig. 3, Additional file 3: Figure S3).

Rodent density versus DON

In all years and plots, we observed a positive association between the DONt+1 and rodent density, which as a predictor explained 61% of the variance (Fig. 5). Our findings are comparable to previous cross-sectional studies performed in the USA [13, 15] and support that rodents are the main hosts of larval ticks and consistently contribute to a new generation of nymphs in the following year [31]. Regarding the contribution of each rodent species in feeding ticks, wood mice were infested at significantly higher levels with larval ticks than bank voles, which has been reported before [28, 55,56,57,58]. The difference in larval tick burden between the two rodent species has been attributed to bank voles acquiring immunity to feeding ticks [59].

Rodent density versus pathogens amplified by rodents

As expected, transmission dynamics of B. afzelii and N. mikurensis were rodent density-dependent. A higher density of rodents increased the probability for larval ticks to feed on an infected rodent, and subsequently significantly increased the NIPt+1 B. afzelii and NIPt+1 N. mikurensis (Fig. 6). Since DON was also rodent density-dependent, there was a strong synergistic effect of rodent density on DINt+1 B. afzelii and DINt+1 N. mikurensis (Fig. 6). We observed a significantly higher NIPN. mikurensis than NIPB. afzelii. Possibly, N. mikurensis-infected rodents are more infectious than while being infected with B. afzelii, which may be due to different tissue tropism of these pathogens in the rodent [26, 55]. In addition, B. afzelii had a higher infection prevalence in bank voles than in wood mice, which has been reported previously [28, 53, 60,61,62]. Although one study has reported the opposite, these studies showed that infectivity of voles was much higher than that of mice [28, 55]. As mentioned above, the larval infestation was higher in wood mice whilst a greater proportion of bank voles was infected with B. afzelii. This indicates that these two rodent species play distinct but complementary roles in B. afzelii transmission dynamics.

There was a significantly negative association between rodent density and NIPt+1 B. microti and a non-linear association with DINt+1 B. microti (Fig. 7). We observed a positive association at low and a negative association at high densities of rodents. We detected the parasite almost exclusively in bank voles; thus, our results might be a consequence of increasing density of wood mouse, which probably is not an amplifying host of B. microti (Fig. 4). An alternative explanation for this non-linear association might be that I. ricinus is not a main vector of this parasite. Previous studies proposed I. trianguliceps, a nidicolous rodent tick species as the main vector [63,64,65]. It indicates that B. microti circulates in the, so called ‘cryptic cycle’ between specialist ticks and rodents, while I. ricinus sporadically becomes infected and perhaps acts as an occasional bridge vector to other host species [63].

Rodent density versus a pathogen amplified by birds

An increasing density of rodents negatively associated with NIPt+1 B. garinii (Fig. 8). This is probably due to the increased number of nymphs uninfected with B. garinii, which fed on the widely abundant rodents, B. garinii-incompetent hosts [66]. Our plots were not large enough to cover the territory of birds, B. garinii-amplifying hosts [35, 67], thus, we speculate that all (or the majority) of the collected B. garinii-infected nymphs were brought by birds from outside the experimental plots and that these events were more or less constant during the course of the study. The DINt+1 B. garinii remained unaltered which suggests that the increase in the DON eliminated the negative effect of rodents on NIPB. garinii (Fig. 8).

Rodent density versus vertically transmitted tick-borne pathogens and microorganisms

We observed a different association between rodent density and the NIPt+1 R. helvetica depending on the year of study (Fig. 9). Although several studies detected R. helvetica in rodent blood and skin samples, and in various ectoparasites feeding on rodents, to date, it is not clear which role rodents play in its transmission cycle [68,69,70]. Other vertebrates were suggested to be amplifying hosts, for instance, songbirds, which were shown to acquire bacteraemia [35, 36]. Here, we detected R. helvetica in rodent ears of both species; however, it is not possible to infer from our results whether rodents acquire systemic infection. On the other hand, we can speculate that R. helvetica causes short-term, localized infection in the skin, which is favourable for co-feeding transmission [34], and this transmission route has been attributed to R. helvetica on many occasions [35, 71, 72]. The infection prevalence in questing nymphs was significantly higher than in rodents, which indicates that ticks are the main amplification hosts of this bacterium. Rodent density was positively associated with DINt+1 R. helvetica, which is unexpected and requires further study (Fig. 9).

Rodent density was indifferently associated with NIPt+1 B. miyamotoi (Fig. 9). In general, the average NIP in questing ticks was only 3% (CI: 2.6–3.4%) and the fluctuations from year to year were small. Borrelia miyamotoi is a predominantly vertically transmitted bacterium, which means that a proportion of unfed larvae originating from an infected female tick, may also be infected [73, 74]. The efficiency of the transmission, in an experimental setting, was shown to vary between 6% and 73% [73]. Thus, it is surprising that despite this mode of transmission, the prevalence of B. miyamotoi in ticks was not higher. A possible explanation for this could be an inefficient horizontal transmission from infected amplification hosts to naïve ticks as it seems that B. miyamotoi does not cause a persistent infection in rodents [75, 76]. Nonetheless, the pathogen can sustain its widespread distribution with only small numbers of ticks being infected. Ultimately, we observed a significant positive effect of rodent density on DINt+1 B. miyamotoi, possibly related to the general increase in DON (Fig. 9).

An increasing density of rodents was significantly negatively associated with NIPt+1 S. ixodetis (Fig. 7). Decreasing NIPt+1 S. ixodetis and increasing DON along with the increasing rodent density resulted in non-linear association between rodents and DINt+1 S. ixodetis (Fig. 7). We observed a positive association at low and a negative association at high densities of rodents. Since from the two rodent species investigated in this study we detected S. ixodetis almost exclusively in wood mouse, a possible explanation is that increasing bank vole populations diluted the prevalence of this bacterium in ticks (Fig. 4).

Another tick-borne Spiroplasma species has been shown to amplify in rodents only in experimental settings [77] and have been reported to cause infections in humans [78, 79]. However, the role of vertebrates in the transmission cycle of S. ixodetis in natural conditions is largely unknown. Nevertheless, the detection of S. ixodetis in rodent ears indicates that these rodents may facilitate horizontal transfer of the bacterium to naïve ticks. Our findings are in line with a recent phylogenetic study, which has revealed that horizontal transmission is probably one of the drivers responsible for spreading of S. ixodetis across tick community [80]. This transmission mode is proposed in addition to the stable vertical transmission, for which spiroplasmas are known [81, 82].

Effect of treatments on rodent density, DON, and DIN

The variation in rodent density throughout the season was comparable with studies from other woodland areas [83]. The density of both rodent species was affected by food resource availability, here acorns, and with our treatment we succeeded to obtain study sites with significantly different densities of rodents (Fig. 1). It allowed us to study the effect of rodent density on tick population dynamics and associated pathogen infections. Our results with acorn addition are in accordance with previous findings, where bank vole and wood mouse populations increase after mast years [16, 19,20,21,22].

Nevertheless, the variation in tick density throughout the years did not follow fluctuations of rodent densities. As a consequence, DON was not affected by our treatment (Fig. 2). There was an effect of the treatment on rodents and an effect of rodents on DON; however, the effect of the treatments on rodents was apparently not enough to establish a significant change in DON. This can also be appreciated from considering the size of the confidence bands in Fig. 5. It is larger than the vertical distance between the treatment lines. In addition, there was no effect of either acorn addition or rodent removal on DINB. afzelii and DINN. mikurensis in the following years (Fig. 3, Additional file 3: Figure S3).

The discrepancy in effect of the treatment indicates that there are additional factors affecting nymphal densities, which expressed high natural variation despite experimental methods. This variation is probably affected by fluctuations in abundance of other vertebrates and/or meteorological conditions affecting seasonal activity of both rodents and ticks. Although in this study we did not assess abundance of other tick hosts, we observed that nymphal activity was affected by temperature, which has been noticed before ([84]; Fig. 2 and Additional file 2: Figure S2). The onset and annual duration of nymphal activity seemed to be related to a number of months with a mean temperature equal or below 7 °C.

Synchrony in activity of rodents and ticks and its influence on transmission dynamics of tick-borne microorganisms

In our study, rodent density had differential effects on NIP and DIN depending on the species of tick-borne microorganism, which indicates that there are additional factors playing a role in microorganism dynamics. Some of these factors might be timing of both activity and infection of rodents and ticks. In temperate European forests, there is a well-documented synchronization between questing larval ticks and rodents, which facilitates the transition of larvae to nymphs [33, 60, 85]. In addition to driving I. ricinus development, rodents contribute to the maintenance of vertically transmitted microorganisms. However, to propagate horizontally transmitted tick-borne pathogens, questing larvae have to be synchronized with infected rodents. Depending on the persistence of a pathogen in a rodent population, rodents may infect larvae directly at the onset of larval activity or after the pathogen has been introduced into the rodent population by infected nymphs. The former situation has been documented for B. afzelii, which causes infection in rodents for life, and therefore often persists over winter [60]. In this study, we observed that B. afzelii-infected rodents were, indeed, present throughout the year, also before the onset of ticks (Fig. 4).

The latter situation is probably applicable to N. mikurensis as the smaller proportion of rodents captured in March was infected with this pathogen, than in later months (Fig. 4). A possible explanation could be that N. mikurensis causes systemic blood infection and decreases the overwintering survival of infected rodents. This phenomenon was observed before in bank voles and Puumala virus (PUUV) despite the expectation that hantaviruses have become well adapted to their rodent hosts during co‐evolution [86, 87]. Thus, the most favorable scenario for N. mikurensis transmission is synchronization in activity of rodents and infected nymphs right before the onset of larvae [88, 89]. In the Netherlands, nymphs have been shown to start their seasonal activity at least one month before larvae [84], which seems to be advantageous for zoonotic pathogens overwintering in nymphs rather than in vertebrate hosts.

Study limitations

This study greatly enhanced our understanding about the role of rodents in the dynamics of tick populations and their associated microorganisms. However, we recognize that our semi-experimental approach has logistic limitations on the temporal and spatial extent that must be acknowledged.

First, our results on mechanisms driving the population of nymphs was measured at a relatively small temporal scale, which is only a transition from a larva to a nymph, and do not necessarily hold at a larger scale involving a complete tick life-cycle. A study of many years following all life stages would have added value and perhaps reveal the robustness of a rodent-tick relationship.

Secondly, the size of the plots was not large enough to cover territory of other vertebrate species, such as deer and birds, for which we had no data on density fluctuations. Since these vertebrates may substantially contribute to the tick and pathogen cycles, it is advisable to increase a plot size and obtain data on vertebrate abundance/arrival rate by, for instance, camera trapping [11].

In addition, increasing the plot size would be also beneficial for more accurate description of rodent population dynamics. It has been shown that along with growth and maturation, rodents change their home range, and therefore depending on population structure, they might have various effects on tick and pathogen populations [90].

Furthermore, in the first study year (2012), the experiment of acorn addition was already ongoing, thus we have no good baseline density of rodents to compare the effect of treatments to. It is advisable, in future field experiments, to have a longer monitoring period prior to implementation of the intervention, in order to have a solid baseline in place. This would also increase the statistical power to detect the effects of an intervention.

Lastly, it should be borne in mind that there was natural variation between plots, even within the experimental settings. Hence the data are obtained in a complex environment where rodent densities can vary by plot, year or treatment. Tick population and infection dynamics is intricately interwoven with the rodent dynamics, and we realize that a more involved modelling exercise is probably needed to fully understand the ecology. However, in the current approach our aim was to be ‘descriptiveʼ of responses of ticks and their infection, rather than finding the most appropriate mechanistic model.

Conclusions

We demonstrated experimentally that increase in rodent density positively affects populations of nymphal ticks in the following year. In addition, we show that prevalence and density of infected ticks with various tick-borne microorganisms are dependent on rodent density to a different extent. These differences probably arouse from varying transmission modes of tick-borne microorganisms and the strongest associations can be observed between rodent density and rodent-associated pathogens that rely on horizontal transmission. Nevertheless, it is not possible to predict disease risk solely on rodent density since we have shown that other factors, independent from our experiment, strongly affected tick density. Our results draw attention to the importance of considering transmission mode of a pathogen as well as other (spatial and temporal) factors while developing models to predict the tick-borne disease risk.

Availability of data and materials

Data are available from the corresponding author on request. Unique DNA sequences of B. burgdorferi (s.l.) were deposited to the GenBank database with the accession numbers MN515318-MN515341.

Abbreviations

qPCR:

quantitative polymerase chain reaction

DNA:

deoxyribonucleic acid

KNMI:

Koninklijk Nederlands Meteorologisch Instituut

References

  1. Rizzoli A, Hauffe H, Carpi G, Vourc HG, Neteler M, Rosa R. Lyme borreliosis in Europe. Euro Surveill. 2011;16:19906.

    PubMed  Google Scholar 

  2. Schwartz AM, Hinckley AF, Mead PS, Hook SA, Kugeler KJ. Surveillance for Lyme Disease—United States, 2008–2015. MMWR Surveill Summ. 2017;66:1–12.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Glass GE, Amerasinghe FP, Morgan JM III, Scott TW. Predicting Ixodes scapularis abundance on white-tailed deer using geographic information systems. Am J Trop Med Hyg. 1994;51:538–44.

    Article  CAS  PubMed  Google Scholar 

  4. Glass GE, Schwartz BS, Morgan JM III, Johnson DT, Noy PM, Israel E. Environmental risk factors for Lyme disease identified with geographic information systems. Am J Public Health. 1995;85:944–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Robertson JN, Gray JS, Stewart P. Tick bite and Lyme borreliosis risk at a recreational site in England. Eur J Epidemiol. 2000;16:647–52.

    Article  CAS  PubMed  Google Scholar 

  6. 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  PubMed  PubMed Central  Google Scholar 

  7. 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.

    Article  CAS  PubMed  Google Scholar 

  8. 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  PubMed  Google Scholar 

  9. Jones CG, Ostfeld RS, Richard MP, Schauber EM, Wolff JO. Chain reactions linking acorns to gypsy moth outbreaks and Lyme disease risk. Science. 1998;279:1023–6.

    Article  CAS  PubMed  Google Scholar 

  10. 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:20170453.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Takumi K, Sprong H, Hofmeester TR. Impact of vertebrate communities on Ixodes ricinus-borne disease risk in forest areas. Parasit Vectors. 2019;12:434.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 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  PubMed  PubMed Central  Google Scholar 

  13. Ostfeld RS, Schauber EM, Canham CD, Keesing F, Jones CG, Wolff JO. Effects of acorn production and mouse abundance on abundance and Borrelia burgdorferi infection prevalence of nymphal Ixodes scapularis ticks. Vector Borne Zoonotic Dis. 2001;1:55–63.

    Article  CAS  PubMed  Google Scholar 

  14. Bogdziewicz M, Zwolak R, Crone EE. How do vertebrates respond to mast seeding? Oikos. 2016;125:300–7.

    Article  Google Scholar 

  15. Ostfeld RS, Canham CD, Oggenfuss K, Winchcombe RJ, Keesing F. Climate, deer, rodents, and acorns as determinants of variation in Lyme-disease risk. PLoS Biol. 2006;4:e145.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Ostfeld RS, Jones CG, Wolff JO. Of mice and mast. BioScience. 1996;46:323–30.

    Article  Google Scholar 

  17. Ostfeld RS, Keesing F. Pulsed resources and community dynamics of consumers in terrestrial ecosystems. Trends Ecol Evol. 2000;15:232–7.

    Article  CAS  PubMed  Google Scholar 

  18. Wolff JO. Coexistence of white-footed mice and deer mice may be mediated by fluctuating environmental conditions. Oecologia. 1996;108:529–33.

    Article  PubMed  Google Scholar 

  19. Clotfelter ED, Pedersen AB, Cranford JA, Ram N, Snajdr EA, Nolan V Jr, et al. Acorn mast drives long-term dynamics of rodent and songbird populations. Oecologia. 2007;154:493–503.

    Article  PubMed  Google Scholar 

  20. Jensen TS. Seed production and outbreaks of non-cyclic rodent populations in deciduous forests. Oecologia. 1982;54:184–92.

    Article  PubMed  Google Scholar 

  21. McShea WJ. The influence of acorn crops on annual variation in rodent and bird populations. Ecology. 2000;81:228–38.

    Article  Google Scholar 

  22. Pucek Z, Jędrzejewski W, Jędrzejewska B, Pucek M. Rodent population dynamics in a primeval deciduous forest (Białowieża National Park) in relation to weather, seed crop, and predation. Acta Theriol. 1993;38:199–232.

    Article  Google Scholar 

  23. Swart A, Bekker DL, Maas M, de Vries A, Pijnacker R, Reusken CB, et al. Modelling human Puumala hantavirus infection in relation to bank vole abundance and masting intensity in the Netherlands. Infect Ecol Epidemiol. 2017;7:1287986.

    PubMed  PubMed Central  Google Scholar 

  24. Tersago K, Verhagen R, Servais A, Heyman P, Ducoffre G, Leirs H, et al. Hantavirus disease (nephropathia epidemica) in Belgium: effects of tree seed production and climate. Epidemiol Infect. 2009;137:250–6.

    Article  CAS  PubMed  Google Scholar 

  25. Schauber EM, Ostfeld RS, Evans AS Jr. What is the best predictor of annual Lyme disease incidence: weather, mice, or acorns? Ecol Appl. 2005;15:575–86.

    Article  Google Scholar 

  26. Burri C, Schumann O, Schumann C, Gern L. Are Apodemus spp. mice and Myodes glareolus reservoirs for Borrelia miyamotoi, Candidatus Neoehrlichia mikurensis, Rickettsia helvetica, R. monacensis and Anaplasma phagocytophilum? Ticks Tick Borne Dis. 2014;5:245–51.

    Article  CAS  PubMed  Google Scholar 

  27. Hu CM, Humair PF, Wallich R, Gern L. Apodemus sp. rodents, reservoir hosts for Borrelia afzelii in an endemic area in Switzerland. Zentralbl Bakteriol. 1997;285:558–64.

    Article  CAS  PubMed  Google Scholar 

  28. Gassner F, Takken W, van der Plas C, Kastelein P, Hoetmer AJ, Holdinga M, et al. Rodent species as natural reservoirs of Borrelia burgdorferi sensu lato in different habitats of Ixodes ricinus in The Netherlands. Ticks Tick Borne Dis. 2013;4:452–8.

    Article  PubMed  Google Scholar 

  29. Karbowiak G. Zoonotic reservoir of Babesia microti in Poland. Pol J Microbiol. 2004;53(Suppl.):61–5.

    PubMed  Google Scholar 

  30. 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  PubMed  Google Scholar 

  31. Hofmeester T, Coipan E, Van Wieren S, Prins H, Takken W, Sprong H. Few vertebrate species dominate the Borrelia burgdorferi s.l. life cycle. Environ Res Lett. 2016;11:043001.

    Article  Google Scholar 

  32. Stanek G. Pandoraʼs Box: pathogens in Ixodes ricinus ticks in Central Europe. Wien Klin Wochenschr. 2009;121:673–83.

    Article  PubMed  Google Scholar 

  33. Gern L, Rais O. Efficient transmission of Borrelia burgdorferi between cofeeding Ixodes ricinus ticks (Acari: Ixodidae). J Med Entomol. 1996;33:189–92.

    Article  CAS  PubMed  Google Scholar 

  34. Voordouw MJ. Co-feeding transmission in Lyme disease pathogens. Parasitology. 2015;142:290–302.

    Article  PubMed  Google Scholar 

  35. 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  PubMed  Google Scholar 

  36. Hornok S, Kovats D, Csorgo T, Meli ML, Gonczi E, Hadnagy Z, et al. Birds as potential reservoirs of tick-borne pathogens: first evidence of bacteraemia with Rickettsia helvetica. Parasit Vectors. 2014;7:128.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Kurtenbach K, Carey D, Hoodless AN, Nuttall PA, Randolph SE. Competence of pheasants as reservoirs for Lyme disease spirochetes. J Med Entomol. 1998;35:77–81.

    Article  CAS  PubMed  Google Scholar 

  38. Taragelʼova V, Koci J, Hanincova K, Kurtenbach K, Derdakova M, Ogden NH, et al. Blackbirds and song thrushes constitute a key reservoir of Borrelia garinii, the causative agent of borreliosis in central Europe. Appl Environ Microbiol. 2008;74:1289–93.

    Article  PubMed  CAS  Google Scholar 

  39. Sunyer P, Muñoz A, Mazerolle MJ, Bonal R, Espelta JM. Wood mouse population dynamics: interplay among seed abundance seasonality, shrub cover and wild boar interference. Mamm Bio. 2016;81:372–9.

    Article  Google Scholar 

  40. Perret JL, Guigoz E, Rais O, Gern L. Influence of saturation deficit and temperature on Ixodes ricinus tick questing activity in a Lyme borreliosis-endemic area (Switzerland). Parasitol Res. 2000;86:554–7.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  42. 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  PubMed  PubMed Central  Google Scholar 

  43. Panzacchi M, Linnell JD, Melis C, Odden M, Odden J, Gorini L, et al. Effect of land-use on small mammal abundance and diversity in a forest-farmland mosaic landscape in south-eastern Norway. Forest Ecol Manag. 2010;259:1536–45.

    Article  Google Scholar 

  44. Leirs H, Stenseth NC, Nichols JD, Hines JE, Verhagen R, Verheyen W. Stochastic seasonality and nonlinear density-dependent factors regulate population size in an African rodent. Nature. 1997;389:176.

    Article  CAS  PubMed  Google Scholar 

  45. Schnabel ZE. The estimation of the total fish population of a lake. Am Math Mon. 1938;45:348–52.

    Google Scholar 

  46. Wielinga PR, Gaasenbeek C, Fonville M, de Boer A, de Vries A, Dimmers W, et al. Longitudinal analysis of tick densities and Borrelia, Anaplasma, and Ehrlichia infections of Ixodes ricinus ticks in different habitat areas in The Netherlands. Appl Environ Microbiol. 2006;72:7594–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  Google Scholar 

  48. 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  PubMed  PubMed Central  Google Scholar 

  49. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Stenos J, Graves SR, Unsworth NB. A highly sensitive and specific real-time PCR assay for the detection of spotted fever and typhus group Rickettsiae. Am J Trop Med. 2005;73:1083–5.

    Article  CAS  Google Scholar 

  51. R Core Team. R: a language and environment for statistical computing. In: R version 3.5.1 (2018-07-02) “Feather Spray” edn. Vienna: R Foundation for Statistical Computing; 2018. https://www.R-project.org/

  52. Team R. RStudio: integrated development for R. In: 0.99.903 edn. Boston: RStudio, Inc.; 2015.

  53. Coipan CE, van Duijvendijk GLA, Hofmeester TR, Takumi K, Sprong H. The genetic diversity of Borrelia afzelii is not maintained by the diversity of the rodent hosts. Parasit Vectors. 2018;11:454.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Mysterud A, Stigum VM, Jaarsma RI, Sprong H. Genospecies of Borrelia burgdorferi sensu lato detected in 16 mammal species and questing ticks from northern Europe. Sci Rep. 2019;9:5088.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. 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.

    Article  CAS  PubMed  Google Scholar 

  56. Humair PF, Turrian N, Aeschilimann A, Gern L. Borrelia burgdorferi in a focus of Lyme borreliosis: epizootiologic contribution of small mammals. Folia Parasitol (Praha). 1993;40:65–70.

    CAS  PubMed  Google Scholar 

  57. Kurtenbach K, Kampen H, Dizij A, Arndt S, Seitz HM, Schaible UE, et al. Infestation of rodents with larval Ixodes ricinus (Acari: Ixodidae) is an important factor in the transmission cycle of Borrelia burgdorferi s.l. in German woodlands. J Med Entomol. 1995;32:807–17.

    Article  CAS  PubMed  Google Scholar 

  58. Nilsson A, Lundqvist L. Host selection and movements of Ixodes ricinus (Acari) larvae on small mammals. Oikos. 1978;31:313–22.

    Article  Google Scholar 

  59. Dizij A, Kurtenbach K. Clethrionomys glareolus, but not Apodemus flavicollis, acquires resistance to Ixodes ricinus L., the main European vector of Borrelia burgdorferi. Parasite Immunol. 1995;17:177–83.

    Article  CAS  PubMed  Google Scholar 

  60. Humair P, Rais O, Gern L. Transmission of Borrelia afzelii from Apodemus mice and Clethrionomys voles to Ixodes ricinus ticks: differential transmission pattern and overwintering maintenance. Parasitology. 1999;118:33–42.

    Article  PubMed  Google Scholar 

  61. Kybicova K, Kurzova Z, Hulinska D. Molecular and serological evidence of Borrelia burgdorferi sensu lato in wild rodents in the Czech Republic. Vector Borne Zoonotic Dis. 2008;8:645–52.

    Article  CAS  PubMed  Google Scholar 

  62. Talleklint L, Jaenson TG. Transmission of Borrelia burgdorferi s.l. from mammal reservoirs to the primary vector of Lyme borreliosis, Ixodes ricinus (Acari: Ixodidae), in Sweden. J Med Entomol. 1994;31:880–6.

    Article  CAS  PubMed  Google Scholar 

  63. Bown KJ, Lambin X, Telford GR, Ogden NH, Telfer S, Woldehiwet Z, et al. Relative importance of Ixodes ricinus and Ixodes trianguliceps as vectors for Anaplasma phagocytophilum and Babesia microti in field vole (Microtus agrestis) populations. Appl Environ Microbiol. 2008;74:7118–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Randolph SE. Quantifying parameters in the transmission of Babesia microti by the tick Ixodes trianguliceps amongst voles (Clethrionomys glareolus). Parasitology. 1995;110:287–95.

    Article  PubMed  Google Scholar 

  65. Cayol C, Jaaskelainen A, Koskela E, Kyrolainen S, Mappes T, Siukkola A, et al. Sympatric Ixodes-tick species: pattern of distribution and pathogen transmission within wild rodent populations. Sci Rep. 2018;8:16660.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Kurtenbach K, Peacey M, Rijpkema SG, Hoodless AN, Nuttall PA, Randolph SE. Differential transmission of the genospecies of Borrelia burgdorferi sensu lato by game birds and small rodents in England. Appl Environ Microbiol. 1998;64:1169–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Humair PF, Postic D, Wallich R, Gern L. An avian reservoir (Turdus merula) of the Lyme borreliosis spirochetes. Zentralbl Bakteriol. 1998;287:521–38.

    Article  CAS  PubMed  Google Scholar 

  68. Burri C, Dupasquier C, Bastic V, Gern L. Pathogens of emerging tick-borne diseases, Anaplasma phagocytophilum, Rickettsia spp., and Babesia spp., in ixodes ticks collected from rodents at four sites in Switzerland (Canton of Bern). Vector Borne Zoonotic Dis. 2011;11:939–44.

    Article  PubMed  Google Scholar 

  69. Miťková K, Berthová L, Kalúz S, Kazimírová M, Burdová L, Kocianová E. First detections of Rickettsia helvetica and R. monacensis in ectoparasitic mites (Laelapidae and Trombiculidae) infesting rodents in south-western Slovakia. Parasitol Res. 2015;114:2465–72.

    Article  PubMed  Google Scholar 

  70. Obiegala A, Oltersdorf C, Silaghi C, Kiefer D, Kiefer M, Woll D, et al. Rickettsia spp. in small mammals and their parasitizing ectoparasites from Saxony, Germany. Vet Parasitol. 2016;5:19–24.

    Google Scholar 

  71. Burgdorfer W, Aeschlimann A, Peter O, Hayes SF, Philip RN. Ixodes ricinus: vector of a hitherto undescribed spotted fever group agent in Switzerland. Acta Trop. 1979;36:357–67.

    CAS  PubMed  Google Scholar 

  72. Zemtsova G, Killmaster L, Mumcuoglu K, Levin MLA. Co-feeding as a route for transmission of Rickettsia conorii israelensis between Rhipicephalus sanguineus ticks. Exp Appl Acarol. 2010;52:383–92.

    Article  CAS  PubMed  Google Scholar 

  73. Scoles GA, Papero M, Beati L, Fish D. A relapsing fever group spirochete transmitted by Ixodes scapularis ticks. Vector Borne Zoonotic Dis. 2001;1:21–34.

    Article  CAS  PubMed  Google Scholar 

  74. van Duijvendijk G, Coipan C, Wagemakers A, Fonville M, Ersoz J, Oei A, et al. Larvae of Ixodes ricinus transmit Borrelia afzelii and B. miyamotoi to vertebrate hosts. Parasit Vectors. 2016;9:97.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  75. Taylor KR, Takano A, Konnai S, Shimozuru M, Kawabata H, Tsubota T. Borrelia miyamotoi infections among wild rodents show age and month independence and correlation with Ixodes persulcatus larval attachment in Hokkaido, Japan. Vector Borne Zoonotic Dis. 2013;13:92–7.

    Article  PubMed  Google Scholar 

  76. 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.

    Article  CAS  PubMed  Google Scholar 

  77. Tully JG, Rose DL, Yunker CE, Carle P, Bov JM, Williamson DL, et al. Spiroplasma ixodetis sp. nov., a new species from Ixodes pacificus ticks collected in Oregon. Int J Syst Evol Micr. 1995;45:23–8.

    CAS  Google Scholar 

  78. Aquilino A, Masiá M, López P, Galiana AJ, Tovar J, Andrés M, et al. First human systemic infection caused by Spiroplasma. 2015;53:719–21.

    Google Scholar 

  79. Lorenz B, Schroeder J, Reischl U. First evidence of an endogenous Spiroplasma sp. infection in humans manifesting as unilateral cataract associated with anterior uveitis in a premature baby. Graefes Arch Clin Exp Ophthalmol. 2002;240:348–53.

    Article  CAS  PubMed  Google Scholar 

  80. Binetruy F, Bailly X, Chevillon C, Martin OY, Bernasconi MV, Duron O. Phylogenetics of the Spiroplasma ixodetis endosymbiont reveals past transfers between ticks and other arthropods. Ticks Tick Borne Dis. 2019;10:575–84.

    Article  PubMed  Google Scholar 

  81. Fukatsu T, Tsuchida T, Nikoh N, Koga R. Spiroplasma symbiont of the pea aphid, Acyrthosiphon pisum (Insecta: Homoptera). Appl Environ Microb. 2001;67:1284–91.

    Article  CAS  Google Scholar 

  82. Hurst GD, Anbutsu H, Kutsukake M, Fukatsu T. Hidden from the host: Spiroplasma bacteria infecting Drosophila do not cause an immune response, but are suppressed by ectopic immune activation. Insect Mol Biol. 2003;12:93–7.

    Article  CAS  PubMed  Google Scholar 

  83. Montgomery WI. Population regulation in the wood mouse, Apodemus sylvaticus. II. Density dependence in spatial distribution and reproduction. J Anim Ecol. 1989;58:477–94.

    Article  Google Scholar 

  84. Hartemink N, van Vliet A, Sprong H, Jacobs F, Garcia-Martí I, Zurita-Milla R, et al. Temporal-spatial variation in questing tick activity in the Netherlands: the effect of climatic and habitat factors. Vector Borne Zoonotic Dis. 2019;19:494–505.

    Article  PubMed  Google Scholar 

  85. Randolph S. Quantitative ecology of ticks as a basis for transmission models of tick-borne pathogens. Vector Borne Zoonotic Dis. 2002;2:209–15.

    Article  PubMed  Google Scholar 

  86. Kallio ER, Voutilainen L, Vapalahti O, Vaheri A, Henttonen H, Koskela E, et al. Endemic hantavirus infection impairs the winter survival of its rodent host. Ecology. 2007;88:1911–6.

    Article  PubMed  Google Scholar 

  87. Nemirov K, Vaheri A, Pljusnin A. Hantaviruses: co-evolution with natural hosts. Rec Res Dev Virol. 2004;6:201–28.

    CAS  Google Scholar 

  88. De Boer R, Hovius K, Nohlmans M, Gray JS. The woodmouse (Apodemus sylvaticus) as a reservoir of tick-transmitted spirochetes (Borrelia burgdorferi) in The Netherlands. Zentralbl Bakteriol. 1993;279:404–16.

    Article  PubMed  Google Scholar 

  89. Talleklint L, Jaenson TG. Is the small mammal (Clethrionomys glareolus) or the tick vector (Ixodes ricinus) the primary overwintering reservoir for the Lyme borreliosis spirochete in Sweden? J Wildl Dis. 1995;31:537–40.

    Article  CAS  PubMed  Google Scholar 

  90. Korn H. Changes in home range size during growth and maturation of the wood mouse (Apodemus sylvaticus) and the bank vole (Clethrionomys glareolus). Oecologia. 1986;68:623–8.

    Article  PubMed  Google Scholar 

  91. Rudzinska MA, Spielman A, Riek RF, Lewengrub SJ, Plesman J. Intraerythrocytic ‛gametocytesʼ of Babesia microti and their maturation in ticks. Can J Zool. 1979;57:424–34.

    Article  CAS  PubMed  Google Scholar 

  92. Bell-Sakyi L, Palomar AM, Kazimirova M. Isolation and propagation of a Spiroplasma sp. from Slovakian Ixodes ricinus ticks in Ixodes spp. cell lines. Ticks Tick Borne Dis. 2015;6:601–6.

    Article  PubMed  PubMed Central  Google Scholar 

  93. Hersh MH, Tibbetts M, Strauss M, Ostfeld RS, Keesing F. Reservoir competence of wildlife host species for Babesia microti. Emerg Infect Dis. 2012;18:1951.

    Article  PubMed  PubMed Central  Google Scholar 

  94. Hanincová K, Taragelová V, Koci J, Schäfer SM, Hails R, Ullmann AJ, et al. Association of Borrelia garinii and B valaisiana with songbirds in Slovakia. Appl Environ Microbiol. 2003;69:2825–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgments

We are grateful to everyone who helped in sample collection during the experiment. We thank to all colleagues who provided valuable comments.

Funding

This study was financially supported by the Dutch Ministry of Health, Welfare and Sport (VWS).

Author information

Authors and Affiliations

Authors

Contributions

GLAD, HS. and WT jointly conceived and designed the study. GLAD and FHHJ performed the field experiment and collected samples. RIJ, MF and AIK performed DNA extraction and molecular detection of microorganisms. AIK and AS performed all statistical analyses. AIK wrote the first draft of the manuscript with the help of DH and HS. All authors contributed to the writing of the final manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Aleksandra I. Krawczyk or Hein Sprong.

Ethics declarations

Ethics approval and consent to participate

The experiment was approved by the animal ethical committee of Wageningen University & Research under permit numbers 2013017 and 2014064.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Figure S1.

An example of an experimental plot 50 × 50 m surrounded by screens to study the effect of rodent density on the density of infected nymphs.

Additional file 2: Figure S2.

Mean monthly temperature (in °C) from August 2012 to December 2015.

Additional file 3: Figure S3. a

Density of nymphs infected with N. mikurensis (DIN N. mikurensis) in 2014 and 2015 in all three treatments in comparison to 2013 (baseline year). b Differences in DIN N. mikurensis between the treatments in two separate years calculated with the Wilcoxon test with a correction for a baseline year (2013). The overall differences between the treatments were not significant either in 2014, or 2015 (P = 0.87 and P = 0.94, respectively).

Additional file 4: Table S1.

Daily temperature data from September 2012 to December 2015 collected from the nearest weather station (Deelen, KNMI, the Netherlands). Text S1. Primers and probes for detection of B. microti and S. ixodetis, and qPCR protocol. Table S2. Density of nymphs and number of analyzed nymphs for microorganisms per treatment, month and year. Table S3. Density of rodents and number of analyzed individuals for microorganisms per rodent species, treatment, month and year. Table S4. Prevalences of tick-borne microorganisms in rodents and nymphs. Table S5. All tested models for prediction of density of nymphs (DON), nymphal infection prevalence (NIP), and density of infected nymphs (DIN). Table S6. Full equations for the best fitting models for prediction of density of nymphs (DON), nymphal infection prevalence (NIP), and density of infected nymphs (DIN).

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Krawczyk, A.I., van Duijvendijk, G.L.A., Swart, A. et al. Effect of rodent density on tick and tick-borne pathogen populations: consequences for infectious disease risk. Parasites Vectors 13, 34 (2020). https://doi.org/10.1186/s13071-020-3902-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-020-3902-0

Keywords