Open Access Research

Circumstantial evidence for an increase in the total number and activity of borrelia-infected ixodes ricinus in the Netherlands

Hein Sprong1*, Agnetha Hofhuis1, Fedor Gassner12, Willem Takken2, Frans Jacobs2, Arnold JH van Vliet3, Marijn van Ballegooijen1, Joke van der Giessen1 and Katsuhisa Takumi1

Author Affiliations

1 Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands

2 Laboratory of Entomology, Wageningen University, Wageningen, The Netherlands

3 Environmental Systems Analyses Group, Wageningen University, Wageningen, The Netherlands

For all author emails, please log on.

Parasites & Vectors 2012, 5:294  doi:10.1186/1756-3305-5-294


The electronic version of this article is the complete one and can be found online at: http://www.parasitesandvectors.com/content/5/1/294


Received:24 September 2012
Accepted:9 December 2012
Published:17 December 2012

© 2012 Sprong et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Between 1994 and 2009, a threefold increase has been observed in consultations of general practitioners for tick bites and Lyme disease in The Netherlands. The objective of this study was to determine whether an increase in the number of questing ticks infected with B. burgdorferi sensu lato is a potential cause of the rise in Lyme disease incidence.

Methods

Historic data on land usage, temperature and wildlife populations were collected and analyzed together with data from two longitudinal field studies on density of questing ticks. Effective population sizes of Borrelia burgdorferi s.l. were calculated.

Results

Long-term trend analyses indicated that the length of the annual tick questing season increased as well as the surface area of tick-suitable habitats in The Netherlands. The overall abundances of feeding and reproductive hosts also increased. Mathematical analysis of the data from the field studies demonstrated an increase in mean densities/activities of questing ticks, particularly of larvae between 2006 and 2009. No increase in infection rate of ticks with Borrelia burgdorferi sensu lato was found. Population genetic analysis of the collected Borrelia species points to an increase in B. afzelii and B. garinii populations.

Conclusions

Together, these findings indicate an increase in the total number of Borrelia-infected ticks, providing circumstantial evidence for an increase in the risk of acquiring a bite of a tick infected with B. burgdorferi s.l. Due to the high spatiotemporal variation of tick densities/activities, long-term longitudinal studies on population dynamics of I. ricinus are necessary to observe significant trends.

Keywords:
Borrelia burgdorferi sensu lato; Ixodes ricinus; Population dynamics; Lyme disease; The Netherlands

Background

Lyme disease, or Lyme borreliosis, is the most prevalent tick-borne infection of humans. In the early stages, Lyme disease is clinically manifested by an erythema migrans, an expanding skin lesion occurring after several days or weeks at the site of the tick bite. Late and more serious Lyme disease is a multi-system disorder with skin, neurological, cardiac and musculoskeletal manifestations [1,2]. Over the last decade, the incidence of Lyme disease has increased in at least nine European countries [3]. In The Netherlands, a long-term retrospective study among all general practitioners (GP) has shown a strong and continuing increase in GP consultations for erythema migrans during the past fifteen years from 4 per 10,000 inhabitants in 1994 up to 13 per 10,000 inhabitants in 2009 [4,5]. A similar magnitude of increase was observed in the number of people that reported a tick bite: The incidence of GP consultations for tick bites was 19 per 10,000 inhabitants in 1994, and increased to 56 tick bites per 10,000 inhabitants in 2009 [4,5]. The most straightforward explanation for the increase in Lyme disease is therefore the increase in the exposure to infected ticks. Despite considerable efforts during the past decades in education of the public, aiming to reduce human exposure to ticks and promote timely removal of ticks from the human skin, the rise in the incidence of Lyme disease continued.

The risk of acquiring Lyme disease depends on many different biological, environmental, and societal factors [6-14]. In short, however, it depends on two main factors: the abundance of questing Ixodes ricinus ticks infected with the causal agent B. burgdorferi sensu lato (s.l.), and the level of human exposure to ticks. In this study, we assessed whether the number of questing ticks infected with B. burgdorferi s.l. in The Netherlands has increased during the past decades. Besides direct methods such as tick sampling from the field, we explore several indirect measures to investigate longitudinal trends in the number of questing Borrelia-infected nymphs.

The number of questing ticks in an area is determined by the tick density and their level of activity, which, in turn, is determined by a complex interplay between vegetation, climatic conditions and the presence of blood hosts. For example, vegetation provides questing sites for ticks, but it also affects micro-climatic conditions, such as humidity and temperature, which determines tick survival and activity [15-18]. Vegetation also influences utilization by host animals and affects questing times by providing different degrees of shelter for ticks. An increase in tick bites may be due to an increase in tick density or activity but could also be due to an increase of tick suitable areas or exposure.

Temporal weather and climate conditions appear to be predictors for the tick activity, and to a lesser extent tick density [19-22]. Nymphal and adult ticks tend to quest for a blood meal once the weekly mean daily maximum temperature exceeds 7°C [23-25]. The development of ticks depends upon the consumption of vertebrate blood. Therefore, the abundance of feeding hosts can affect the abundance of ticks. All mobile life stages of Ixodes ricinus can feed on a broad range of warm- and cold-blooded vertebrate hosts [12,21,26,27]. Ixodes ricinus larvae infest small mammals, but also feed on larger animals such as roe deer. Nymphs and adults usually feed on medium-sized and large mammals [28]. These differences are probably due to the differential vertical distribution of instars [29]. To investigate trends in the relative abundance of blood hosts during the past decades, we analyzed readily available data on the abundance of roe deer and fallow deer, birds, rodents, and proxies for the abundance of rodents such as birds of prey [30,31].

The best available estimate for the tick density is the activity of questing ticks measured by standard blanket dragging, despite the limitation that moulting, resting and feeding ticks are not caught [32,33]. Furthermore, this technique is not equally efficient in different vegetation types, and therefore will not provide absolute numbers of questing ticks from a given area. Measuring the number of ticks searching for a blood meal, and testing these ticks for infection with B. burgdorferi s.l. provides an estimate of local and temporal variation in public health risk. To search for longitudinal trends in the number of questing ticks and tick infection with B. burgdorferi s.l., tick sampling data from two field studies was analyzed. Our analyses involved a long-term field study at a single location [18,34,35], and an ongoing longitudinal field study at thirteen locations geographically spread throughout The Netherlands [18].

Finally, we applied a population genetic approach as a complementary indirect exploration for changes in the abundance of B. burgdorferi s.l. infected ticks that were collected in these field studies. For this purpose, we used genetic information of Borrelia-DNA found in ticks to construct phylogenetic trees, mapping the genetic divergence within and between B. burgdorferi genospecies. Using this reconstruction of the recent demographic history of B. burgdorferi, we looked for indications of changes in the genetic diversity over time, which may suggest changes in the sylvatic transmission of B. burgdorferi s.l. in The Netherlands.

Methods

Surface of tick suitable areas

A database of land usage in The Netherlands (LGN) was used as a primary data source (Alterra, Wageningen). The total area of each land type during the periods 1999–2000, 2003 – 2004 and 2007–2008 was extracted from LGN using ArcGIS 9.3 (ESRI: Redlands CA). Resolution of the map was 25m by 25m. The country was partitioned into 42 land types in LGN. To each of the 42 land types, we assigned one of three new labels indicating habitats with probabilities for a high, low, or zero-tick density (see Additional file 1: Table S1). These values are partially based on several (ongoing) nationwide surveys [18,34,36]. Forests, dunes, and natural grasslands were assigned the category ‘high tick density area’. Pastures, orchards, tree(s) and grass patches within (sub) urban areas, bog and peat area, and heather were assigned the categories ‘low tick density area’. Agricultural area, (sub) urban areas, motor and railways were assigned the category ‘zero tick density area’. The total area for each of the three habitats was calculated by adding the areas of corresponding land types.

Additional file 1 Table S1. Land usage types were assigned to three categories indicating habitats with probabilities for a high, low, or zero-tick density.

Format: DOC Size: 68KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Length of tick questing activity season

Mean daily temperatures at the weather station in De Bilt during 1985–2010 were extracted from the online database (Wolfram Research, Inc., Champaign, IL). The number of days above 7°C was counted annually and a linear trend was assessed using a generalized linear model for Poisson counts. Data points are shown in Additional file 2: Figure S1. The analysis was repeated using only the days above 7°C of the first half year (January until June), to focus more on the onset of the annual tick questing activity and to avoid the inhibitory effects of the shortening day length in the autumn on the questing activity of ticks [24].

Additional file 2 Figure S1. A: Number of days above 7 degrees Celsius since 1985. B: Number of days above 7 degrees Celsius since 1985 measured in the first half years (January – June).

Format: PPT Size: 78KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

Abundance of vertebrate hosts: birds and small rodents

We analyzed a database containing a total of 215 distinct species of birds that were monitored across the country between 1990 and 2008 (Netwerk Ecologische Monitoring, SOVON&CBS, http://www.sovon.nl webcite). Sixty-five species of winter birds and twenty species of water birds were excluded from the analysis because they were deemed less relevant for the sylvatic transmission of Borrelia[37]. Population index values for the remaining group of 130 bird species were described by the negative binomial distribution, where the annual mean index value is an exponential function in time with the rate r. The presence (r > 0) or the absence (r = 0) of trends in the population index values in time was assessed by the likelihood ratio test. Lacking direct measurements of rodent populations, we selected bird of prey species as indicator species for a trend in the rodent population densities on a nationwide scale [30,31]. Barn owl (Tyto alba) and common buzzard (Buteo buteo) are birds of prey that specialize on mice and other small rodent species, and they are distributed widely in The Netherlands. A trend in these bird population indexes was assessed as described for other bird species.

Abundance of vertebrate hosts: roe deer and fallow deer

Roe deer (Capreolus capreolus) are especially important for maintaining (and determining the level of) tick populations because they feed reproductive females [22]. After hatching and diapause, the larvae become active the following spring [38]. In Western Europe, roe deer is one of the most abundant propagation hosts of I. ricinus[39,40]. To detect a possible association between roe deer population density and tick abundance in the following year, the trend in national population density of roe deer in 1980, 2003, 2004, and 2008 (Dutch Royal Hunting Society, KNJV) was analyzed using a generalized linear model. For the spatial correlation analysis, location-specific roe deer population density estimates (in 2008) at the tick-sampling sites were extracted from the roe deer database and paired with the estimated numbers of questing ticks of all stages in 2009. Data on roe deer and fallow deer (Dama dama) abundance from Duin and Kruidberg, a coastal dune area, were available from the wildlife management unit in Noord Holland [41]. Correlation between roe deer and questing ticks was assessed using generalized linear models with a gamma distribution and the negative binomial link.

Data on tick densities/activities and infection rate

Data on tick densities/activities were derived from two large field studies conducted in The Netherlands [18,34,35]. In one study, ticks were collected by standard blanket dragging at thirteen locations geographically spread throughout the country from 2006 to 2009 at monthly intervals [18]. Criteria for selection of these study areas were based on the likelihood of finding any I. ricinus ticks, and the availability of a team of trained volunteers to make monthly blanket drag collections [18]. In a separate longitudinal study, a single area in the coastal dunes, Duin and Kruidberg, between 2000 and 2009 was selected at that time because of its unusually high tick density/activity. This area is rich in vegetation and open to the public for recreation [18,34,35]. Monthly tick collections were performed from April to October. For the distribution of the study areas in The Netherlands, see Additional file 3: Figure S2.

Additional file 3 Figure S2. Map of the Netherlands showing the 14 study sites where ticks were sampled.

Format: PPT Size: 391KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

Seven genospecies of Borrelia burgdorferi s.l. known to circulate in Europe have previously been identified in The Netherlands by reverse line blot analysis [34]. The prevalence of Borrelia genospecies was determined by PCR followed by reverse line blotting. DNA sequences of the variable 5S-23S intergenic spacer region were determined from 211 randomly selected B. burgdorferi sensu lato isolates from infected nymphs (98%) and adults (2%), as described by Wielinga et al.[34]. These genospecies were categorized by previously defined reference sequences [42]. Borrelia bavariensis, a recently defined genospecies, was also found by sequencing some B. garinii isolates. As our Reverse Line Blot cannot distinguish between B. garinii and B. bavariensis, little can be inferred for the prevalence of this newly named genospecies in The Netherlands.

A model for the numbers of questing ticks

Numbers of questing ticks are highly variable measurements. Hence, a model describing a tick sampling process is necessary to make inference about a multi-year trend in tick densities. The variation in questing tick densities was greater than Poisson. Thus, the probability that i-ticks are collected at a particular study site on a particular sampling time was empirically shown by the negative binomial distribution [43,44].

<a onClick="popup('http://www.parasitesandvectors.com/content/5/1/294/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/5/1/294/mathml/M1">View MathML</a>

The parameter m is the mean number of questing ticks per study site. The parameter k is a measure of aggregation where a highly aggregated number of questing ticks is indicated by a small value of k. Seasonality and other location-specific conditions are likely sources of variations in the numbers of questing ticks. Thus, we wrote the mean (m) as a product of the seasonal variation (s) and a multi-year trend (u)

<a onClick="popup('http://www.parasitesandvectors.com/content/5/1/294/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/5/1/294/mathml/M2">View MathML</a>

The time (t) is expressed as a fraction of the sampling period where one year was scaled to a period of 2π. Seasonality is a trigonometric function

<a onClick="popup('http://www.parasitesandvectors.com/content/5/1/294/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/5/1/294/mathml/M3">View MathML</a>

The mean number of questing ticks is therefore periodic and it reaches the maximum questing activity when t = π + τ. The multi-year trend is the linear function

<a onClick="popup('http://www.parasitesandvectors.com/content/5/1/294/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/5/1/294/mathml/M4">View MathML</a>

The parameter a, is the mean tick density/activity in the baseline year per drag area. The parameter b is the change in tick density/activity per time unit.

Parameter estimates

Numbers of ticks captured at different times at the study sites were used to estimate the parameters of the model; annual mean number of questing ticks (a), multi-year trend (b), the moment of peak questing activity (τ), and the parameter k. The set of parameters was estimated specifically for each study site by maximizing the log-likelihood of the questing tick counts. Probability (P-values) that the null hypothesis of a constant tick density/activity is the correct model (i.e. whether b = 0) was calculated for each study site by likelihood ratio. P-values less than or equal to 0.05 were considered to be a significant support for the model describing the multi-year trend.

The annual total estimated abundance of questing ticks

The total number of questing ticks in one year was calculated by integrating the estimated tick activity function over one sampling year for each site, and subsequently averaging the resulting total numbers over the study sites. The overall change during the sampling period was calculated by dividing the total annual tick activity for 2009 by the total annual tick activity for 2006.

Trends in the prevalence of ticks with B. burgdorferi sensu lato

Nymphal and adult ticks, collected from the field studies, were tested for the presence of B. burgdorferi by PCR followed by reverse line blotting. Numbers of tested ticks were highly variable. To remove effects of variable sampling sizes on the estimate for the prevalence, the numbers of questing infected ticks were described by the binomial distribution. The probability that a tick is positive for B. burgdorferi was modeled by the logit function.

<a onClick="popup('http://www.parasitesandvectors.com/content/5/1/294/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.parasitesandvectors.com/content/5/1/294/mathml/M5">View MathML</a>

To test seasonality, we set z = x + yCos(t − θ) where time t is expressed as a fraction of the sampling period where one year was scaled to 2π. To test for site-specific prevalence, we determined two clusters of the study sites using the squared Euclidean distance of the raw prevalence values. The parameters x, y and θ were estimated by maximizing the likelihood of observed numbers of ticks with B. burgdorferi. Using the likelihood ratio we assessed whether all study sites are characterized by the same seasonal oscillation or one cluster of study sites was characterized by higher prevalence values than the other cluster.

Population dynamics of Borrelia burgdorferi s.l

The DNA sequences of the intergenic 5s-23S spacer region from Borrelia isolates derived from questing ticks were used to estimate changes in the effective population size of infected ticks. This analysis was limited to the two main Borrelia genospecies for which a sufficient number of samples were available: B. afzelii (n=86) and B. garinii (n=60). Note that in these two cases the number of samples is relatively small for coalescent analyses [45]. We excluded the so-called B. ruskii subclade from the B. afzelii samples, because of the unusual sequence divergence between the -ruskii and the other -afzelii samples (see Additional file 4: Figure S4).

Additional file 4 Figure S4. Genetic divergence within and between Borrelia burgdorferi s.l. species based on the IGS-sequence.

Format: PPT Size: 161KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

We tested changes of genetic diversity of Borrelia over time by calculating Fu’s Fs statistics using Arlequin v.3.5 [46]. A large negative value of Fs indicates an excess of rare alleles, as would be expected from a recent population expansion or from genetic hitchhiking [47]. We also used the software Beast for Bayesian analysis of molecular sequences [48] to estimate coalescent effective population size with skyline plots [45]. Skyline plots represent a nonparametric flexible method to estimate changes in the coalescent effective population size, and here it is used as a measure for changes in genetic diversity, which is complementary to Fu’s Fs. We used the general time reversible model of DNA evolution with site specific mutation rates (gamma distribution) and default priors. Simulations were run for 30 million updates after discarding burn-in.

Results

Increase in surface of tick suitable areas

Between 2000 and 2008, habitats suitable for high-tick densities increased in surface area by 19% (+847 km2, Table 1), predominantly by expansion of forest areas. Low-tick density habitats increased only marginally in area (Table 1). Zero-tick density habitats decreased by 8% (−1192 km2, Table 1).

Table 1. Areas in The Netherlands in 2000 – 2008 categorized according to three levels of probable tick densities

Length of tick questing activity

The mean number of days above 7°C (the temperature above which nymphal and adult ticks become active) has increased linearly with the year since 1985 from 235 to 262 days in 2009 (P-value = 0.01, see also Additional file 2: Figure S1A). To exclude the adverse effects of the shortening day length in the autumn on the questing activity of ticks [24], the analysis was repeated using only the days above 7°C of the first half year (January until June). Between the years 1985 and 2009, the number of days in the first half of the year when the average daily temperature exceeded 7 degrees Celsius increased significantly (p = 0.002) from 96 days in 1985 to 118 days in 2009 (see also Additional file 2: Figure S1B). Changing the threshold temperature to 10°C (the temperature above which larval ticks also become active), the mean number of days above 10°C has increased linearly with the year since 1985 from 177 to 204 days in 2009 (P-value = 0.003). This proxy for the tick questing season length has been increasing since 1985. Trends in reported tick bites and erythema migrans have only been available since 1994. In the shorter time window of 1994 – 2009, only the mean number of days above 10°C significantly increased.

Abundance of birds and birds of prey as proxy for density of small rodents

The analysis of the longitudinal bird records revealed that the pooled populations of breeding birds increased significantly during the period 1990 – 2008 (P-value < 0.001, Figure 1). In 2008, there were 1.4 birds for each breeding bird in the baseline year, 1990, i.e. a 40% increase (Figure 1). This could be an indication that the tick host availability, and therefore feeding capacity, has increased during the past twenty years.

thumbnailFigure 1. Bird population densities 1980–2008. Dots represent population index values of each bird species per year. The bird population index was calculated by dividing the number of a bird species each year by the number of that bird species in the baseline year in 1990 and multiplying the result by 100. Lines represent the best-fit negative binomial distribution with the mean (solid) and the 95% confidence interval (dashed).

Since 1990, the barn owl population index in the country has increased at the rate of 0.084 birds per year (P-Value < 0.0001) (Figure 2). The common buzzard population in the country has also increased during 1990–2007 at the rate of 0.033 birds per year (P-value < 0.001). As populations of birds of prey are assumed to be correlated with rodent abundance [30,31], these data suggest that the abundance of rodents increased between 1990 and 2007.

thumbnailFigure 2. Change in the number of Barn owl (Tyto alba) 1990 – 2006. The change in the number of Barn owl was calculated by dividing the number of barn owls each year by the number in the baseline year 1990 and multiplying the result by 100. The solid line is the best fit exponential function.

Abundance of roe deer

The estimated nationwide population size of roe deer in The Netherlands has increased linearly between 1980 and 2008 (P-value = 0. 0009), growing each year on average by 1820 animals to ca. 70,000 animals in 2008. Roe deer population densities at the fourteen tick-sampling sites in 2008 were associated significantly (P-value = 0.005) with the numbers of questing ticks of all stages at the corresponding sites in the following year (Figure 3). When we tested the same association in the 2-year period between 2004 and 2006 at the fourteen sites, the association was significant as well (P-value = 0.008). An increase over low deer densities (0–3) has rather little effect on tick abundance until the upper end of (>5) the range of deer density (Figure 3). Thus, the increasing availability of blood meal sources is positively associated with tick reproduction (natality), host finding or tick survival.

thumbnailFigure 3. Roe deer population densities and questing tick abundance. Horizontal axis is the local population density of Roe deer at all fourteen study sites in 2008. Vertical axis is the peak number of questing ticks at the corresponding study sites in 2009. The solid line is the best generalized-linear-model fit.

Tick density/activity in the 13 study areas

All ticks that were collected at the thirteen study sites (Additional file 3: Figure S2) in monthly intervals from 2006 to 2009, were identified to be I. ricinus[18]. By applying an empirical model-based statistical analysis to each study site, we detected significant three-year trends in the density/activity of questing ticks at six study sites (Table 2). The cumulative number of questing larvae, nymphs, and adult ticks increased in density/activity at five of these study sites and decreased at one other site. A three-year trend was absent at the remaining seven sites (Table 2). The estimated overall density/activity of questing ticks increased by 21% during the three-year period (2006 – 2009). The number of nymphs and adult ticks increased significantly only at two study sites and decreased at one site (Table 3). At the remaining ten study sites, significant trends in tick density/activity were not observed. The estimated overall density/activity of questing nymphs and adults increased by 4% during the same period. Hence, the increase in the total tick density/activity was mainly due to questing larval ticks.

Table 2. Estimates for the parameters in the model describing the densities of questing larval, nymphal and adult Ixodes ricinus 2006 – 2009

Table 3. Estimates for the parameters in the model describing the densities of questing nymphal plus adult Ixodes ricinus 2006 – 2009

Case study: longitudinal population sampling of ticks in duin and kruidberg

In Duin and Kruidberg, the density of ticks (sum of instars) increased during the last decade (P-value = 0.055). The abundance of larvae increased significantly at this location (P-value < 0.001, Figure 4A). In contrast, the density of nymphs and adult ticks did not increase (P-value = 0.9, Figure 4B). The abundance of small rodents was assessed in the vicinity of the tick-sampling sites (2000 – 2008). Four species were identified among the sampled animals: Sorex araneus, Apodemus sylvaticus, Microtus arvalis, and Myodes glareolus. The total numbers of these rodents per year were 1 (2000), 35 (2001), 34 (2002), 11 (2006), 6 (2007), and 45 (2008; see also: [49]). No significant change in the rodent population was observed in the sampling period. Abundance data of birds from this study area were not available.

thumbnailFigure 4. Numbers of questing ticks at Duin and Kruidberg 2000 – 2009. The population of larvae is shown in panel A. The population of nymphs plus adults is shown in panel B. A dot represents the number of ticks collected per drag (100 m2). Shaded area represents the best-fit for negative binomial distributions with the 95% confidence interval.

Both roe deer and fallow deer (Dama dama) are present in Duin and Kruidberg. The population of roe deer has tripled from 200 in 1985 to 600 in 1997, and it has been stable since that time, varying between 400 – 600 roe deer per year. In contrast, in the same time period the population of fallow deer increased significantly. Only a few fallow deer were present in 1990 whereas the population of fallow deer in 2003 exceeded 800. Thus, since the start of our tick surveillance in Duin and Kruidberg in 2000, the total propagation capacity for ticks has steadily increased.

Trends in the prevalence of ticks with B. burgdorferi sensu lato

A total of 372 out of 3104 nymphs and adults, collected between 2000 and 2009, tested positive for the presence of B. burgdorferi s.l. The prevalence showed seasonal variation (P-value < 0.001; see Additional file 5: Figure S3). A spatial variation in Borrelia prevalence was also detected (P-value < 0.001). The Borrelia prevalences of several study sites (Bilthoven, Ede, Eijsden, Gieten, Kwade Hoek, Montferland) were higher (annual mean 20%) compared to the other cluster of the remaining study sites (annual mean 8%). Only from one location (Duin and Kruidberg), long-term prevalence data were available: Despite seasonal and yearly variation, no significant increase or decline in Borrelia infection in questing ticks was observed (data not shown).

Additional file 5 Figure S3. Seasonal infection rates (fraction) of Borrelia Borrelia burgdorferi s.l. The study areas were separated into two groups: One group (black line) had an annual mean of 20% (Bilthoven, Ede, Eijsden, Gieten, Kwade Hoek, Montferland). The other group (dotted line) had an annual mean of 8% (Appelscha, Duin and Kruidberg, Hoog Baarlo, Schiermonnikoog, Twiske, Vaals, Veldhoven, Wassenaar).

Format: PPT Size: 147KB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

Population dynamics of B. burgdorferi s.l

The most prevalent genospecies was B. afzelii (~65%), whereas B. garinii, B. burgdorferi sensu stricto, B. spielmanii and B. valaisiana were found less frequently (~10% each). Borrelia lusitaniae was identified for the first time in 2009 at one location.

A change in the genetic diversity of Borrelia over time is an indication that the sylvatic transmission of Borrelia in The Netherlands at large is changing. Negative estimates of Fs statistics indicate an excess of rare sequences due to a recent population expansion or genetic hitch hiking of B. afzelii (Fu’s Fs = −9.56, P-value = 0.02) and B. garinii (Fu’s Fs = −5.55, P-value = 0.03). The skyline plot appears to be more compatible with population growth within the last decade for B. afzelii and B. garinii (Figure 5) than an alternative explanation involving genetic hitchhiking.

thumbnailFigure 5. Skyline plot of B. afzelii and B. garinii populations. Changes in coalescent population sizes over years are shown in solid lines indicating B. afzelii (panel A) and B. garinii (panel B). Dashed lines enclose 95% credible interval.

Discussion

Several independent analyses on various aspects of the tick cycle were performed in this study, which provide circumstantial evidence for a rise in abundance of questing infected ticks. In the period 2000 – 2008 the total surface area of habitats suitable for ticks increased by 19%. Although the actual tick density/activity was not measured in any of the newly arisen “high-risk” areas, it is expected that these areas are readily inhabited by ticks and their vertebrate hosts as they are generally adjacent to established tick-suitable areas.

The length of tick questing activity, based on the number of days above 7°C and 10°C per year since 1985, has increased. The onset of tick activity, as measured by the number of days above 7°C in the first half of the year, became significantly earlier over time. Here, only very simple and limited models were used. Other relevant climatic factors, such as saturation deficit, were not taken into account. Furthermore, climatic effects on the tick cycle as such, the vegetation, vertebrate hosts, diapause and human (leisure) activities were not taken into account either.

The long-term dynamics of vertebrate hosts is an indirect method for measuring changes in tick densities [21,26]. As the densities of tick hosts increase, it is to be expected that tick abundance will increase concomitantly. At least two principal sources of vertebrate blood for I. ricinus, namely the population of birds and deer have increased significantly in The Netherlands since the end of the last century. Also, a simultaneous increase of the predatory bird species barn owls as well as common buzzards was observed. Barn owl and common buzzard prey on small rodents and in this study are assumed to be indicator species for the nationwide rodent density [30,31]. Other factors, such as intensified nature conservation activities, may have resulted in the increase of birds of prey as well. The effect of host abundance on tick density is partially supported by our findings in the Duin and Kruidberg area. The density of larval ticks has increased significantly since 2000, which may be due to an increase in local populations of roe deer and fallow deer [41]. Remarkably, no simultaneous increase in nymphal and adult tick density/activity was observed in this area. Deer not only act as blood hosts for larvae and nymphs, they are the primary host for the adult deer tick and are key to the reproductive success of the tick [12,22,50]. Perhaps, an increasing abundance of deer resulted in a (much) better reproductive success of ticks, but only a minor/insignificant increase in the total number of feeding hosts.

Analysis of the total numbers of I. ricinus ticks collected in the 13 field sites between 2006 and 2009 showed an overall increase of 21% in the total tick-density/activity similar to what was observed in Duin and Kruidberg. The larval population attributed the most to the observed change in density, similar to that observed in Duin and Kruidberg. The 4% change in densities of nymphs and adult ticks may not be more than change variation. Densities of nymphs and adults may be more subject to density-dependence than larvae, which probably brought the densities back towards the long-term equilibrium. The massive reproductive potential of ticks allows larval numbers to vary widely, but become stabilized by density-dependent mortality at later stages [24].

As the majority of people acquire Lyme disease from infected questing nymphs and adult ticks, the observed limited increase in nymphal plus adult density/activity alone cannot satisfactorily explain the increase in tick bites and patients with erythema migrans. These direct measurements on tick density/activity seem to contradict the indirect measurements, such as the abundance of feeding and propagation hosts. The high spatiotemporal variation in tick densities/activities in combination with the relatively short-term measurement (~ 1 tick cycle) may be one explanation. Collecting data for longer periods may result in a more robust long-term trend analysis. Furthermore, the strongest increase in tick density/activity would have occurred in forest areas that arose between 2000 and 2008. Unfortunately, none of these areas were included in the longitudinal studies.

Another indirect method for measuring the human risk of acquiring Lyme disease is measuring the population dynamics of Borrelia in questing ticks. Using a population genetic approach, strong indications are found for recent population expansions of the two most abundant Borrelia species, B. afzelii and B. garinii in questing ticks, which are known to cause Lyme disease. These findings are explained best by an increase in tick abundance, as no indications were found that the fraction of ticks infected with B. burgdorferi s.l. altered in the Duin and Kruidberg region between 2000 and 2009. The three-year prevalence of B. burgdorferi s.l. has not increased at other locations in The Netherlands either [9], supporting our findings. Overall, increased tick abundance and increased Borrelia populations, without increase in Borrelia infection rate, can be explained by increase in the surface of tick- suitable areas.

The expansion of suitable habitats for I. ricinus concomitant with an uncontrolled increase in the roe deer population underscores the necessity to investigate how a sustainable natural environment can be maintained while minimizing the threats of tick-borne diseases in densely populated areas such as The Netherlands. These investigations should not only focus on Lyme disease, but also on other tick-borne pathogens [9,35,36,51-54]. Measures to reduce the size of the tick population, particularly in recreational areas, might be one sensible option to reduce the human risk for Lyme disease. Intervention measures may include habitat management, deer management and the introduction of biological control [55]. Alternatively, reduction of tick exposure may be facilitated by for example (temporal) disclosure of high-risk areas for the general public and by continued public information efforts focusing at prompt removal of attached ticks.

Conclusions

Several biotic and abiotic factors that contribute to the overall abundance of questing infected ticks in The Netherlands have concomitantly increased in the past decades. Demographic analysis on genetic data from Borrelia in questing ticks supports these findings. Together, it supports the idea of an overall increase in abundance of questing infected ticks, and provides one explanation for the increase in reported tick bites and erythema migrans. The expected increase in the overall abundance of questing infected ticks was only marginally supported by longitudinal studies in 14 sites across the country between 2006 and 2009.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

FG, AV, FJ, WT, JG organized and participated in the fieldwork for the collection and tick data. KT, AH, and HS collected and analyzed data on incidence of tick bites and Erythema migrans, land use, climate, and data on vertebrate abundance. MB, KT and HS collected and analyzed data on Borrelia infection rate and performed population genetic analyses. HS and KT drafted the manuscript and wrote the final version. All authors read and approved the final manuscript.

Acknowledgements

We are grateful to the many volunteers who collected a substantial part of the ticks used for the model validation. We thank Margriet Montizaan of the Royal Dutch Hunting Association (Koninklijke Nederlandse Jagers Vereniging) for allowing us to access the database containing regional roe deer population densities. Moreover, we are grateful to Cor Gaasenbeek and Fred Borgsteede (Animal Science Group, Lelystad), for their collaboration in the Duin and Kruidberg study. We thank Marieta Braks and the two anonymous reviewers for critically reading the manuscript, their helpful suggestions and highly constructive criticism.

This study was financially supported by the Dutch Food and Consumer Product Safety Authority (NVWA) and by the Ministry of Health, Welfare and Sport (VWS).

References

  1. Girschick HJ, Morbach H, Tappe D: Treatment of Lyme borreliosis.

    Arthritis Res Ther 2009, 11(6):258. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  2. Stanek G, Fingerle V, Hunfeld KP, Jaulhac B, Kaiser R, Krause A, Kristoferitsch W, O'Connell S, Ornstein K, Strle F, et al.: Lyme borreliosis: clinical case definitions for diagnosis and management in Europe.

    Clin Microbiol Infect 2012, 17(1):69-79. OpenURL

  3. Hubalek Z: Epidemiology of lyme borreliosis.

    Curr Probl Dermatol 2009, 37:31-50. PubMed Abstract OpenURL

  4. Hofhuis A, van der Giessen JW, Borgsteede FH, Wielinga PR, Notermans DW, van Pelt W: Lyme borreliosis in the Netherlands: strong increase in GP consultations and hospital admissions in past 10 years.

    Euro Surveill 2006, 11(6):E060622.

    060622

    PubMed Abstract OpenURL

  5. Hofhuis A, Harms MG, Van der Giessen JWB, Sprong H, Notermans DW, Van Pelt W: Ziekte van Lyme in Nederland 1994–2009: aantal huisartsconsulten blijft toenemen. Is voorlichting en curatief beleid genoeg?

    Infectieziekten, Bulletin 2010, 3(21):84-87. OpenURL

  6. Lommano E, Bertaiola L, Dupasquier C, Gern L: Infections and coinfections of questing ixodes ricinus ticks by emerging zoonotic pathogens in western switzerland.

    Appl Environ Microbiol 2012, 78(13):4606-4612. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Heyman P, Cochez C, Hofhuis A, van der Giessen J, Sprong H, Porter SR, Losson B, Saegerman C, Donoso-Mantke O, Niedrig M, et al.: A clear and present danger: tick-borne diseases in Europe.

    Expert Rev Anti Infect Ther 2010, 8(1):33-50. PubMed Abstract | Publisher Full Text OpenURL

  8. Huegli D, Moret J, Rais O, Moosmann Y, Erard P, Malinverni R, Gern L: Prospective study on the incidence of infection by borrelia burgdorferi sensu lato after a tick bite in a highly endemic area of Switzerland.

    Ticks Tick Borne Dis 2011, 2(3):129-136. PubMed Abstract | Publisher Full Text OpenURL

  9. Tijsse-Klasen E, Jacobs JJ, Swart A, Fonville M, Reimerink JH, Brandenburg AH, van der Giessen JW, Hofhuis A, Sprong H: Small risk of developing symptomatic tick-borne diseases following a tick bite in the Netherlands.

    Parasit Vectors 2011, 4:17. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  10. Godfrey ER, Randolph SE: Economic downturn results in tick-borne disease upsurge.

    Parasit Vectors 2011, 4:35. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  11. Stefanoff P, Rosinska M, Samuels S, White DJ, Morse DL, Randolph SE: A national case–control study identifies human socio-economic status and activities as risk factors for tick-borne encephalitis in poland.

    PLoS One 2012, 7(9):e45511. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Gray JS: The ecology of ticks transmitting Lyme borreliosis.

    Exp Appl Acarology 1998, 22:249-258. Publisher Full Text OpenURL

  13. Gray JS, Dautel H, Estrada-Pena A, Kahl O, Lindgren E: Effects of climate change on ticks and tick-borne diseases in europe.

    Interdiscip Perspect Infect Dis 2009 2009:593232. OpenURL

  14. Gray JS, Stanek G: Risk assessment in Lyme borreliosis.

    Zentralbl Bakteriol 1998, 287(3):175-176. PubMed Abstract OpenURL

  15. Randolph SE, Storey K: Impact of microclimate on immature tick-rodent host interactions (acari: ixodidae): implications for parasite transmission.

    J Med Entomol 1999, 36(6):741-748. PubMed Abstract OpenURL

  16. Bisanzio D, Amore G, Ragagli C, Tomassone L, Bertolotti L, Mannelli A: Temporal variations in the usefulness of normalized difference vegetation index as a predictor for ixodes ricinus (acari: ixodidae) in a borrelia lusitaniae focus in Tuscany, central Italy.

    J Med Entomol 2008, 45(3):547-555. PubMed Abstract | Publisher Full Text OpenURL

  17. Dister SW, Fish D, Bros SM, Frank DH, Wood BL: Landscape characterization of peridomestic risk for Lyme disease using satellite imagery.

    AmJTrop Med Hyg 1997, 57(6):687-692. OpenURL

  18. Gassner F, van Vliet AJ, Burgers SL, Jacobs F, Verbaarschot P, Hovius EK, Mulder S, Verhulst NO, van Overbeek LS, Takken W: Geographic and temporal variations in population dynamics of ixodes ricinus and associated borrelia infections in the Netherlands.

    Vector Borne Zoonotic Dis 2011, 11(5):523-532. PubMed Abstract | Publisher Full Text OpenURL

  19. Schulze TL, Jordan RA: Influence of meso- and microscale habitat structure on focal distribution of sympatric ixodes scapularis and amblyomma americanum (acari: ixodidae).

    J Med Entomol 2005, 42(3):285-294. PubMed Abstract | Publisher Full Text OpenURL

  20. Schulze TL, Jordan RA, Schulze CJ, Hung RW: Precipitation and temperature as predictors of the local abundance of ixodes scapularis (acari: ixodidae) nymphs.

    J Med Entomol 2009, 46(5):1025-1029. PubMed Abstract | Publisher Full Text OpenURL

  21. 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(6):e145. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Dobson AD, Taylor JL, Randolph SE: Tick (ixodes ricinus) abundance and seasonality at recreational sites in the UK: hazards in relation to fine-scale habitat types revealed by complementary sampling methods.

    Ticks Tick Borne Dis 2011, 2(2):67-74. PubMed Abstract | Publisher Full Text OpenURL

  23. 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(7):554-557. PubMed Abstract | Publisher Full Text OpenURL

  24. Randolph SE: Tick ecology: processes and patterns behind the epidemiological risk posed by ixodid ticks as vectors.

    Parasitology 2004, 129(Suppl):S37-65. PubMed Abstract OpenURL

  25. Macleod J: Ixodes ricinus in relation to its physical environment. II. The factors governing survival and activity.

    Parasitology 1935, 27:123-144. Publisher Full Text OpenURL

  26. Rizzoli A, Hauffe HC, Tagliapietra V, Neteler M, Rosa R: Forest structure and roe deer abundance predict tick-borne encephalitis risk in Italy.

    PLoS One 2009, 4(2):e4336. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. 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(4):837-861. PubMed Abstract | Publisher Full Text OpenURL

  28. 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(6):880-886. PubMed Abstract OpenURL

  29. Mejlon HA, Jaenson TGT: Questing behaviour of ixodes ricinus ticks (acari: ixodidae).

    Expl appl acarology 1997, 21:747-754. Publisher Full Text OpenURL

  30. Southern HN: The natural control of a population of tawny owls (strix aluco).

    J Zool (London) 1970, 162:197-285. OpenURL

  31. Southern HN, Lowe VPW: Predation by tawny owls (srrir aluco) on bank voles (clethrionomys glareolus) and wood mice (apodemus sylvaticus).

    J Anim Ecol 1982, 198:83-102. OpenURL

  32. Daniels TJ, Falco RC, Fish D: Estimating population size and drag sampling efficiency for the blacklegged tick (acari: ixodidae).

    J Med Entomol 2000, 37(3):357-363. PubMed Abstract | Publisher Full Text OpenURL

  33. Talleklint-Eisen L, Lane RS: Efficiency of drag sampling for estimating population sizes of ixodes pacificus (acari: ixodidae) nymphs in leaf litter.

    J Med Entomol 2000, 37(3):484-487. PubMed Abstract | Publisher Full Text OpenURL

  34. Wielinga PR, Gaasenbeek C, Fonville M, de Boer A, de Vries A, Dimmers W, Akkerhuis Op Jagers G, Schouls LM, Borgsteede F, van der Giessen JW: 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(12):7594-7601. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Sprong H, Wielinga PR, Fonville M, Reusken C, Brandenburg AH, Borgsteede F, Gaasenbeek C, van der Giessen JW: Ixodes ricinus ticks are reservoir hosts for rickettsia Helvetica and potentially carry flea-borne rickettsia species.

    Parasit Vectors 2009, 2(1):41. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  36. Tijsse-Klasen E, Fonville M, Reimerink JH, Spitzen-van der Sluijs A, Sprong H: Role of sand lizards in the ecology of Lyme and other tick-borne diseases in the Netherlands.

    Parasit Vectors 2010, 3:42. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  37. Humair PF: Birds and borrelia.

    Int J Med Microbiol 2002, 291(Suppl 33):70-74. PubMed Abstract OpenURL

  38. Randolph SE, Green RM, Hoodless AN, Peacey MF: An empirical quantitative framework for the seasonal population dynamics of the tick ixodes ricinus.

    Int J Parasitol 2002, 32(8):979-989. PubMed Abstract | Publisher Full Text OpenURL

  39. Kiffner C, Lodige C, Alings M, Vor T, Ruhe F: Abundance estimation of ixodes ticks (acari: ixodidae) on roe deer (capreolus capreolus).

    Exp Appl Acarology 2010, 52(1):73-84. Publisher Full Text OpenURL

  40. Vor T, Kiffner C, Hagedorn P, Niedrig M, Ruhe F: Tick burden on European roe deer (capreolus capreolus).

    Exp Appl Acarology 2010, 51(4):405-417. Publisher Full Text OpenURL

  41. Faunabeheerplan Noord-Holland.

    [http://www.fbenoordholland.nl webcite]

    OpenURL

  42. Postic D, Assous MV, Grimont PA, Baranton G: Diversity of Borrelia burgdorferi sensu lato evidenced by restriction fragment length polymorphism of rrf (5S)-rrl (23S) intergenic spacer amplicons.

    Int J Syst Bacteriol 1994, 44(4):743-752. PubMed Abstract | Publisher Full Text OpenURL

  43. Anderson RM, May RM: Infectious Diseases of Humans: Dynamics and Control. Oxford University Press, Oxford; 1992. OpenURL

  44. Wilson K, Grenfell BT: Generalized linear modelling for parasitologists.

    Parasitol Today 1997, 13(1):33-38. PubMed Abstract | Publisher Full Text OpenURL

  45. Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences.

    Mol Biol Evol 2005, 22(5):1185-1192. PubMed Abstract | Publisher Full Text OpenURL

  46. Excoffier L, Lischer HE: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows.

    Mol Ecol Resour 2010, 10(3):564-567. PubMed Abstract | Publisher Full Text OpenURL

  47. Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection.

    Genetics 1997, 147(2):915-925. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees.

    BMC Evol Biol 2007, 7:214. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  49. Akkerhuis Jo, G.A.J.M , Dimmers WJ, van Kats RJM, van der Hoef M, Gaasenbeek CPH, Schouls LM, can de Pol I, Borgsteede FHM:

    Observations on lyme disease, tick abundance and the diversity of small mammals in different habitats in the Netherlands and implications for the 'inflation' and 'dilution' effect approaches. 2003, 31-38. [Annual meeting of entomologists in the Netherlands: 2003] OpenURL

  50. Gray JS, Kahl O, Janetzki C, Stein J: Studies on the ecology of Lyme disease in a deer forest in county Galway, Ireland.

    J Med Entomol 1992, 29(6):915-920. PubMed Abstract OpenURL

  51. Jahfari S, Fonville M, Hengeveld P, Reusken C, Scholte EJ, Takken W, Heyman P, Medlock J, Heylen D, Kleve J, et al.: Prevalence of neoehrlichia mikurensis in ticks and rodents from north-west Europe.

    Parasit Vectors 2012, 5(1):74. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  52. Tijsse-Klasen E, Fonville M, Gassner F, Nijhof AM, Hovius EK, Jongejan F, Takken W, Reimerink JR, Overgaauw PA, Sprong H: Absence of zoonotic bartonella species in questing ticks: first detection of bartonella clarridgeiae and rickettsia felis in cat fleas in the Netherlands.

    Parasit Vectors 2011, 4(1):61. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  53. Wielinga PR, Fonville M, Sprong H, Gaasenbeek C, Borgsteede F, Giessen JW: Persistent detection of babesia EU1 and babesia microti in ixodes ricinus in the Netherlands during a 5-year surveillance: 2003–2007.

    Vector Borne Zoonotic Dis 2009, 9(1):119-122. PubMed Abstract | Publisher Full Text OpenURL

  54. Sprong H, Tijsse-Klasen E, Langelaar M, De Bruin A, Fonville M, Gassner F, Takken W, Van Wieren S, Nijhof A, Jongejan F, et al.: Prevalence of coxiella burnetii in ticks after a large outbreak of Q fever.

    Zoonoses Public Health 2012, 59:69-75. PubMed Abstract | Publisher Full Text OpenURL

  55. Piesman J, Eisen L: Prevention of tick-borne diseases.

    Annu Rev Entomol 2008, 53:323-343. PubMed Abstract | Publisher Full Text OpenURL