Research | Open | Published:
Temporal dynamics of the tick Ixodes ricinus in northern Europe: epidemiological implications
Parasites & Vectorsvolume 10, Article number: 166 (2017)
Tick-borne pathogens pose an increasing threat to human and veterinary health across the northern hemisphere. While the seasonal activity of ticks is largely determined by climatic conditions, host-population dynamics are also likely to affect tick abundance. Consequently, abundance fluctuations of rodents in northern Europe are expected to be translated into tick dynamics, and can hence potentially affect the circulation of tick-borne pathogens. We quantified and explained the temporal dynamics of the tick Ixodes ricinus in the northernmost part of its European geographical range, by estimating (i) abundance in vegetation and (ii) infestation load in the most common rodent species in the study area, the bank vole Myodes glareolus.
Ixodes ricinus nymphs and adult females, the life stages responsible for the most of tick bites in humans, peaked in May-June and August-September. Larvae and nymphs were simultaneously active in June and abundance of questing larvae and nymphs in the vegetation showed a positive association with bank vole abundance. Moreover, infesting larvae and nymphs were aggregated on bank voles, and the infestation of bank voles with I. ricinus larvae and nymphs was positively associated with bank vole abundance.
Our results indicate early summer and early autumn as periods of increased risk for humans to encounter I. ricinus ticks in boreal urban forests and suggest a 2 years life-cycle for I. ricinus with two cohorts of ticks during the same year. Moreover, we identified a simultaneous activity of larvae and nymphs which allows co-feeding on the rodent host, which in turn supports the transmission of several important zoonotic tick-borne pathogens. Finally, we showed that a high density of the rodent host may enhance the risk that ticks and, potentially, tick-borne pathogens pose to human health.
Tick-borne pathogens are a growing burden for European public health policies [1–3]. The current observed increase in tick-borne disease incidence in Europe may be explained by the geographical expansion of Ixodes ricinus, the growing share of space between humans and wild animals, and the improvement of diagnostics tools [4–6]. The epidemiology of tick-borne zoonoses, such as Lyme borreliosis, anaplasmosis or tick-borne encephalitis (TBE), depends on tick abundance and population dynamics, infection prevalence within the tick population, and land use that may affect human exposure to ticks [7, 8]. In order to predict the risks that tick-borne diseases pose to humans, an assessment of factors underlying the temporal variation of tick abundance is necessary.
The abundance of I. ricinus varies in time and space and is highly dependent on environmental conditions, including habitat quality, host availability, and abiotic conditions [9–12]. In northern Fennoscandia, at the northernmost part of the European range of I. ricinus, abiotic conditions undergo extreme seasonal variation; there are 145 to 160 days of snow cover with short day-lengths, during which ticks are not active. This is followed by a quick elevation in temperature leading to a short summer with long day-lengths . In these conditions, I. ricinus activity is likely to show distinctive seasonal patterns, which have not been characterized to date (but see  for southwest Finland).
Ixodes ricinus is dependent on vertebrate hosts to complete its life-cycle. Larvae typically feed on small vertebrates, such as rodents; nymphs, the more common biting stage for humans, parasitize mostly medium-sized mammals; and adults feed mainly on large hosts, such as deer [15, 16]. The population dynamics of ticks and rodents are expected to be linked: some studies have indicated delayed density dependence of questing nymphs on rodent abundance, suggesting that high rodent abundance provides augmented opportunities for successful larvae feeding and nymph development [17, 18]. The bank vole (Myodes glareolus) is a common rodent species throughout Europe ; this species is commonly infested by immature I. ricinus [16, 20, 21]. In northern Europe, vole population abundance shows both seasonality, driven by seasonal breeding, and multiannual density fluctuations shaped by predation, food availability and food quality [22–24]. These seasonal and multiannual density fluctuations are likely to be translated into the dynamics of ticks, and consequently, into the epidemiology of tick-borne pathogens. To date, there are few studies that have investigated the association between the dynamics of cyclic small rodents and ticks .
The bank vole is also an important reservoir host for many tick-borne pathogens, such as Borrelia afzelii, tick-borne encephalitis virus (TBEV) and Babesia microti [26, 27]. Typically, tick larvae acquire infections from an infected rodent host that has become infected while feeding infected nymph(s) . Alternatively, larvae acquire infections via simultaneous feeding with infected nymphs without systemic infection of the host [28, 29]. Infectivity is transstadially maintained in the tick to the following life stage .
Here, we present results from a 4 years of longitudinal bank vole monitoring and tick sampling in central Finland at the northernmost part of the European range of I. ricinus, where abiotic conditions undergo extreme seasonal variation. Our primary aim is to characterize temporal dynamics and quantify the importance of host related factors and abiotic conditions on temporal dynamics of I. ricinus. We also aim to identify seasonal patterns that are relevant for tick-borne pathogen circulation, with the ultimate goal of providing information concerning the risk of tick-borne diseases in our study area.
Sampling took place monthly from May to October in 2012 − 2015 in four periurban forests in the Jyväskylä area in Central Finland: (Kylmänoro (62°13′36.220″, 25°45′1.739″); Jyskänlaakso (62°13′55.398″, 25°49′34.269″); Hämeenlahti (62°12′40.119″, 25°47′11.052″); and Sippulanniemi (62°11′9.019″, 25°44′58.147″) . One trapping period within a month will be referred to as “session” in the following paragraphs. Forests were dominated by Scots pine (Pinus sylvestris) and silver birch (Betula pendula) or by spruce (Picea abies). The herbaceous stratum was typically composed of Vaccinium myrtillus, V. vitis-idaea, Maianthemum bifolium, Linnaea borealis and Oxalis acetosella.
Monthly tick dragging was performed during or within a few days of the vole trapping, using a 1 × 1 m cotton flannel flag sewed to a wooden rod . The fabric was randomly dragged over the vegetation for 300 − 500 m per site around the rodent trapping transects and checked every 20−25 m for ticks, which were removed with tweezers and stored in alcohol at -20 °C. No dragging was performed during rain. In October 2014, due to early snow cover, dragging was not performed. Due to the duration, coverage and interval of the dragging (less than 30 min, 300–500 m2 once a month in each site) it is unlikely that the flag dragging affected the overall tick population abundance and it should not have interfered with the ticks parasitizing rodents in the area.
Vole trapping and tick infestation on voles
As the active tick population consists in parasitizing, questing and resting ticks, sampling targeted questing ticks and parasitizing ticks on their rodent host. This latter buffers the effects of microclimate changes and rodent sampling, in particular, also buffers the effect of the patchy distribution of larvae .
Vole trapping was carried out with two lines of 10 Ugglan Special multiple-capture live traps (Grahnab Company, Sweden), positioned 10 − 15 m apart, located near to rodent burrows. Traps were prebaited for 1−3 nights with sunflower seeds (Helianthus annuus), after which traps were set with sunflower seeds (for food) and a piece of potato (for water) for two consecutive nights. Wood shavings were provided as bedding in wet or cold weather. Traps were checked once per day and trapped voles were handled and sampled before release close to their capture site. Bycatch of species other than voles, as well as recapture of the same individual during the same session, were released immediately on site.
All trapped voles were marked individually with electronic identification chips (microchip Trovan Unique™), which were injected subcutaneously at their first capture. During each capture, voles were identified, body mass was measured as a proxy for age (as in e.g. ), and sex and reproductive condition were recorded. The presence of fleas was recorded and all voles were examined for ticks, with special attention to the area around the ears and face. All ticks were removed with tweezers and stored in alcohol at -20 °C until further identification. All ticks - both those removed from rodents and those collected from vegetation - were identified to species level and life stage under a dissection microscope using morphological identification keys [34–36]. Species identification of seven ticks identified as I. ricinus and three as I. trianguliceps was further confirmed with PCR following a method described elsewhere . Briefly, PCR targeted the mitochondrial 16S rRNA gene and the amplicons obtained were successfully sequenced for eight of the ten ticks. Thereafter sequence identity was determined by BLAST search against the NCBI Nucleotide database and the obtained sequences confirmed our morphological tick identification.
We assessed the overall bank vole population abundance by computing the overall minimum number of voles alive (MNA) at a given trapping session (t) as follows: total number of individuals caught at a given trapping session (t) summed with the total number of individuals marked when caught during subsequent sessions, but not caught at (t) .
We trapped 658 bank voles, an average of 1.53 times (range 1–6), for a total of 1007observations for which all variables described above were available. The minimum number of voles alive per session varied from 5 (in May 2013) to 120 individuals (in September 2014). Three other rodent species were bycaught, consisting of 52 observations of yellow-necked mouse (Apodemus flavicollis), one observation of field vole (Microtus agrestis), and two observations of house mouse (Mus musculus) (Additional file 1: Figure S3).
Ticks in vegetation
We characterized the temporal activity of I. ricinus in the vegetation (i.e. collected by flagging), by examining tick questing activity separately for each life stage, i.e. larvae, nymphs, adults (males and females), in relation to the following variables: year (2012 − 2015), month (May-October), estimated bank vole abundance per given session (MNA), abundance of other life stages present during the same session (number of ticks/100 m2), and the abundance of previous tick life stages collected during the previous session (for larvae: adult, for nymph: larvae, for adults: nymphs). To further identify the effect of current climatic conditions on tick activity, we computed the mean daily saturation deficit (SatDef, in millimetres of mercury) during tick flagging days, based on daily average humidity (in percent) and daily average temperature (in °C) [9, 39, 40] recorded at the meteorological station of Nenäinniemi in Jyväskylä, located 0.72−3.7 km from the study sites (http://www.jyv-weather.info/index.php) (Additional file 1: Figure S1). SatDef was used as an explanatory variable rather than month, with which it showed collinearity. Thus, the second set of models included SatDef and its second-degree polynomial term SatDef2, MNA, current and previous tick abundances as described above. Furthermore, the abundance of nymphs and females pooled together was also modelled with two sets of models: the first one included vole abundance, month and year and the second one included year, vole abundance, SatDef and SatDef2.
Models were fitted using generalized linear mixed models (GLMM) with a negative binomial error distribution (with log-link function) and site was included as a random effect to control for potential pseudoreplication . To take into account the variation in the distance flags were dragged, an offset term (log(distance flagged)/100) was introduced in the models. The model selection (provided in Additional file 1: Table S2) was an automated selection process starting from the full model and based on AICc (Akaike Information Criteria corrected for small sample size ), using dredge function in R software. We kept the most parsimonious model that lay within 2AICc difference from the best model fitted  (Additional file 1: Tables S1 and S2).
Ticks infesting voles
Ixodes ricinus infestation load on bank voles was examined separately for larvae and nymphs. We assessed whether tick infestation showed seasonality and/or between year variation and whether it was affected by individual host characteristics or by concomitant parasitism (by other tick stages, other tick species or fleas). For that purpose, we fitted a GLMM with a negative binomial error distribution to test the fixed effects of month, year, vole sex, body mass (centred value) and its second order polynomial term, presence of fleas, presence of other life stages of I. ricinus and I. trianguliceps, body mass*vole abundance (MNA) interaction term and body mass*sex interaction term. ‘Trapping site’ and ‘vole individual nested in the trapping site’ were included as random effects in the models. Model selection was performed as described before except that we utilized the function drop1 in R software (Additional file 1: Tables S4 and S5).
All statistical analyses were performed with R version 3.2.3 (2015, The R Foundation for Statistical Computing), and using the packages stats (http://www.R-project.org/), MASS (https://cran.r-project.org/web/packages/MASS/index.html), glmmADMB (http://glmmadmb.r-forge.r-project.org/) and MuMIn (https://cran.r-project.org/web/packages/MuMIn/index.html).
Ticks in vegetation
We sampled and identified 943 I. ricinus larvae, 867 nymphs, 239 (adult) females and 294 males from the vegetation. The mean abundance of I. ricinus per session and per area varied from 0 to 22.7 ticks/100 m2 when considering all tick life stages and from 0 to 6.25 ticks/100 m2 when taking into account only female adults and nymphs. Overall, the density of questing ticks collected from vegetation was 7.1/100 m2 (Additional file 1: Figure S2). The ratio between I. ricinus larvae, nymphs and adults was 3.5:3.3:2.0. In addition, one I. trianguliceps nymph was identified.
Models revealed unimodal questing patterns for larvae, which were mostly found in June. Conversely, a bimodal questing pattern was found for nymphs as well as nymph and females modelled together, with the highest abundances found in May-June and September, which therefore appears as the higher risk period for tick bites on humans. Questing adults (males and females) were more abundant in May-June and August-September, and their abundance varied with year, with the highest abundance found in 2015 (Table 1, Fig. 1).
For any given month, the abundances of questing I. ricinus larvae and nymphs showed positive associations with vole abundance (Tables 1, 2; Fig. 2; Additional file 1: Figure S4). For each addition of one individual to the bank vole population, an increase of larvae abundance by approximately 3% and of nymph abundance by 1% was predicted (Fig. 2; Additional file 1: Figure S4).
The abundance of questing larvae was positively associated with saturation deficit, while the abundance of questing adults showed a non-linear relationship with saturation deficit. The number of adults found in vegetation was positively associated with saturation deficit until an optimal value (3.16 mm Hg), after which the abundance of adult ticks was negatively affected by any further increase in saturation deficit. Nymph abundance was not associated with saturation deficit in the best model selected (Table 2).
We found a positive relationship between the number of questing larvae and the abundance of adults observed in the vegetation one session before. Nymph abundance increased with adult abundance during the same flagging session whereas adult abundance showed a negative relationship with nymph abundance during the previous session (Table 2).
Ticks infesting voles
From bank voles, two tick species were identified: I. trianguliceps, the vole tick and I. ricinus. The proportion of infestation with either of these tick species was 75.8%. The ratio of I. ricinus larvae to nymphs found feeding on bank voles was 13:1. The total number of ticks sampled from voles was 3564, out of which 14 ticks could not be identified due to poor condition (Additional file 1: Table S3 and Figure S5).
Models revealed a clear seasonal pattern in the infestation burden of I. ricinus larvae on bank voles (Table 3, Fig. 4): larval infestation underwent seasonality, with a peak in June and a trough in August-October. The highest infestation level was in 2013 and the lowest in 2014. In addition, bank vole infestation load with I. ricinus nymphs underwent seasonal variation, with a peak in May, but was stable between years (Table 4, Fig. 3).
For any given month, nymph infestation on voles was positively associated with bank vole abundance (Table 4). Similarly, larval infestation level increased with bank vole abundance, but the increase was more pronounced among female bank voles than among males (Table 3, Fig. 4). Moreover, infestation with larvae was positively associated with the amount of questing larvae observed in the environment (Table 3), whereas the bank vole infestation load with nymphs was not associated with the amount of questing nymph (i.e. the abundance of questing nymphs was not selected in the best model, Additional file 1: Table S5).
Tick infestation intensity on a host varied with individual characteristics such as age, sex and co-infestation. The oldest male bank voles (i.e. those with highest body mass) were the most intensely infested with larvae (Table 3). Moreover, bank vole infestation load with I. ricinus larvae was positively associated with co-infesting I. trianguliceps females and nymphs and I. ricinus nymphs (Table 3), whereas the infestation with I. ricinus nymphs increased with the presence of I. trianguliceps larvae and females (Table 4). In addition, the infestation load with nymphs showed a non-linear relationship with body mass: infestation load increased until voles reached 32.4 g, whereupon any further increase in body mass led to a reduction of the infestation burden (Table 4).
In this study, we characterized the temporal dynamics of I. ricinus by assessing its abundance in the vegetation and its infestation load in one of its main host in the northernmost part of its range. We focused on bank voles, which parasitism provides insightful information concerning the local immature tick communities. Moreover, we identified risk periods - when humans are likely to encounter tick bites - in boreal forests and seasonal patterns that might be relevant for tick-borne pathogen circulation.
We identified that the highest tick abundance was in early summer (May-June) and early autumn (August-September), which are consequently the periods of increased risk for humans to encounter I. ricinus ticks in boreal forests. The same pattern of bimodal questing activity was previously found in southern Finland for nymph and adult ticks in coniferous and deciduous forests, whereas larvae showed a bimodal occurrence with a larger peak in September than in June . Overall, two types of tick questing activity patterns have been described in Europe: in highly seasonal climates, such as those in central Europe, a bimodal questing activity with early spring and autumn peaks has been described for all life stages of I. ricinus . However, in milder climates, with less climatic variation between seasons, only one peak of activity was observed for all life stages; in either spring or early summer . In the present study, nymphs and adults showed bimodal activity, whereas larvae showed a unimodal activity pattern. This unimodal activity pattern could arise from egg production during the preceding year, the product of which overwintered as eggs or as larvae  or from egg production during the same spring. It could be argued that the inclusion of a year*month interaction term in the model would have captured between year seasonal variations suggested by the raw data (Additional file 1: Figure S1), and would have revealed both unimodal and bimodal activity patterns for larvae. However, data from a longer time series would be needed in order to clarify this point.
Our results seem to indicate the coexistence of two age cohorts of ticks during the same year. Larvae detected on bank voles and in the vegetation in early summer become nymphs in autumn, which can exhibit two different behaviors: immediate questing behavior in autumn; or activity postponed until the next spring after a behavioral diapause [45–47]. In our study, the largest peak of nymphal activity was observed in spring, suggesting that the second behavior was predominant . In addition, we found a peak in questing adults 2 years after the largest larval infestation, indicating a probable 2-year period between larvae blood meal and adults, and suggesting a 2 to 3 year life-cycle from egg to adult for I. ricinus in our study area. Furthermore, in an additional model, the amount of questing nymphs was explained, amongst others explanatory variables, by the total amount of larvae that fed on bank vole the year before (GLMM negative binomial: estimate (± SE) = 0.005 ± 0.002, P = 0.0071, see Additional file 1: Tables S6 and S7). This model confirmed firstly, that the variation in bank vole larval infestation was translated into nymph abundance and secondly, a 1 year delayed relationship between bank vole larval infestation and questing nymphs.
We observed an effect of saturation deficit on larval and adult questing behavior, but not on nymphs, as described in other studies [40, 45]. Ticks respond to microclimate, but climatic variations measured in this study presumably reflect only roughly microclimatic variations and could explain the lack of association between saturation deficit and nymph activity found in our study. On the other hand, nymphs might also be acclimatized to local conditions and therefore their questing behavior may vary compared to nymphs studied in other locations [48, 49]. This is further supported by the optimal saturation deficit value of 3.16 mm Hg over which the adult questing activity decreased, when an optimum of 4.4 mm Hg has been previously noted elsewhere .
Our ratio larvae:nymph:adult was 3.5:3.3:2.0 and was therefore different from theoretical biological expectation (100:10:2) , indicating a possible underestimation of nymphs, and particularly of larvae, which show an important patchiness in distribution. Indeed, the blanket dragging technique is limited by variability in sampling efficiency given the nature of the substrate, the wind speed during sampling, and the height, type and growth stage of the vegetation [50, 51]. Moreover, the total tick population is not accessible by flagging given that diapausing ticks, parasitizing ticks, quiescent ticks or rehydrating individual are not questing in vegetation. Associating bank vole screening to the blanket dragging provided a broader view of immature ticks’ population by beneficiating from the buffer effect that the host offers against the larvae patchiness and the drop of activity in case of unfavourable microclimate . The two approaches are complementary and had led to similar results supporting further the idea that bank voles play an important role as host for immature ticks in the area.
Dynamics of tick in immature stages and bank vole population are related
Questing and parasitizing abundances of larvae and nymphs showed positive associations with bank vole abundance during a given session. Regarding nymphs, this positive relationship might arise from better engorgement success for larvae in high bank vole abundance. However, regarding larvae, this correlation does not imply a causative relationship since bank voles do not contribute to larvae production, which relies on large mammal availability . Consequently, this positive relationship between the abundance of larvae and bank voles might reflect large-mammal density variations or might reveal a functional response: larvae may increase questing behavior in response to increased chemical signals produced by large bank vole populations . This hypothesis requires further attention and needs to be experimentally quantified. Additionally, abundance of other species known to host adult stages needs to be quantified.
The largest burden of nymph parasitism in voles was observed in May, whereas peaks of questing nymphs in the vegetation were observed in May-June and September. In May, vole populations are mainly composed of overwintered sexually active adults; highly mobile males exhibit large home ranges in their search for receptive females . Therefore, the probability of encountering questing nymphs present at a low level in the recovering spring vegetation is increased . In September, bank vole contact rate with nymphs might be lower due to taller vegetation, which allows nymphs to quest higher on plants, where they can contact larger mammal hosts . Moreover, the bank vole develops an acquired resistance to ticks, leading to a significant reduction of infestation success after the first infestation [56, 57], which could lead to poor infestation success during the second nymph peak in September. However, our data (Table 3) provide little support for this hypothesis as regards larvae infestation that increases with animal weight, which is used here as a proxy for age, when a decrease in the relationship was expected under acquired immunity hypothesis. As a consequence, the main period for larval and nymph co-infestation on bank voles is in early summer. The epidemiological consequences of these co-infestations are discussed below. Concerning larvae, we identified an infestation peak in June, which is in accordance with the peak of larvae questing activity and in accordance with previous surveys .
Male bank voles were more commonly infested with nymphs than females and the infestation increased with bank vole abundance. This sex-specific infestation load has been described previously  and may not only be due to the immunosuppressing role of testosterone [60–62], but also to sex-specific behavioral differences, e.g. in home range sizes . Surprisingly, we found larvae infestation differs with population density; females carried more larvae at high population density, whereas males carried more nymphs at any population density. A different use of vertical space by bank vole males and females in high population densities can be hypothesized, leading males to come into contact with more nymphs that quest higher in vegetation, whereas females, which exhibit aggressive defensive behavior against intruders during the reproductive season , would stay close to the ground, i.e. at larvae level. More attention should be paid to the use of vertical space by bank voles in order to clarify the potential role of vertical space use causing differences between individuals in their tick infestation load.
Our data show a concomitant early summer questing activity between larvae and nymphs, and a parasitic aggregation between larvae and nymphs of I. ricinus on bank voles, which are relevant from an epidemiological point of view. The simultaneous activity of larvae and potentially infected nymphs occurs when rapidly rising temperatures in spring allow the simultaneous emergence of larvae and nymphs from overwintering diapause. In these conditions, pathogen transmission from infected nymphs to susceptible larvae can occur via simultaneous feeding on the same host, even without systemic infection of the host. This co-feeding transmission pathway is important for several zoonotic tick-borne pathogens, especially those with short-lived or non-systemic infections in the rodent host, such as Anaplasma phagocytophilum or tick-borne encephalitis virus (TBEV), respectively [65–68].
Synchronous infestations on bank voles
In addition to aggregation between I. ricinus life stages on bank voles, we found a significant aggregation between tick species, with I. ricinus infestation load increasing with the presence of I. trianguliceps. I. trianguliceps is a nidicolous species associated with rodents and insectivores, which does not quest in the vegetation and hence does not come into contact with humans . Even if it is not involved in zoonotic transmission, I. trianguliceps is responsible for maintaining the enzootic cycle of potential zoonotic pathogens such as Anaplasma phagocytophilum [70, 71] or Babesia microti [72, 73]. Both of these pathogens have been identified in Finnish bank voles . Ixodes trianguliceps could contribute to the sylvatic cycle of pathogens that the generalist I. ricinus could transmit to humans, who are considered as dead-end hosts. Hence, the between-species ectoparasite aggregation is also relevant from an epidemiological point of view.
In northern European urban forests, population dynamics of bank voles and questing I. ricinus larvae and nymphs are related, suggesting higher tick abundance and consequently higher risk of tick-borne pathogens for human during the rodent population peak. Larvae and nymphs showed synchronous activity, which increases the transmission opportunity for several pathogens and which are the prerequisite conditions for the maintenance of some pathogens such as TBEV. Further studies should focus on assessing the prevalence of tick-borne pathogens in the bank vole and in questing ticks in order to specify the zoonotic risk. Recent models demonstrate a dampening of vole population cycles in northern Europe , which could therefore be translated into the population dynamics of ticks.
Heyman P, Cochez C, Hofhuis A, van der Giessen J, Sprong H, Porter SR, et al. A clear and present danger: tick-borne diseases in Europe. Expert Rev Anti Infect Ther. 2010;8:33–50.
Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Peña A, George J-C, et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit Vectors. 2013;6:1.
Lindgren E, Jaenson TGT. Lyme borreliosis in Europe: influences of climate and climate change, epidemiology, ecology and adaptation measures. Copenhagen, Denmark: WHO Regional Office for Europe; 2006.
Jaenson TG, Jaenson DG, Eisen L, Petersson E, Lindgren E. Changes in the geographical distribution and abundance of the tick Ixodes ricinus during the past 30 years in Sweden. Parasit Vectors. 2012;5:8.
Rizzoli A, Hauffe HC, Carpi G, Vourc’h GI, Neteler M, Rosà R. Lyme borreliosis in Europe. Euro Surveill. 2011;16:ppi = 19906.
Sormunen J, Penttinen R, Klemola T, Hänninen J, Vuorinen I, Laaksonen M, et al. Tick-borne bacterial pathogens in southwestern Finland. Parasit Vectors. 2016;9:168.
Estrada-Peña A, Jongejan F. Ticks feeding on humans: A review of records on human-biting Ixodoidea with special reference to pathogen transmission. Exp Appl Acarol. 1999;23:685–715.
Vanwambeke SO, Sumilo D, Bormane A, Lambin EF, Randolph SE. Landscape predictors of tick-borne encephalitis in Latvia: land cover, land use, and land ownership. Vector Borne Zoonotic Dis. 2010;10:497–506.
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.
Pfäffle M, Littwin N, Muders SV, Petney TN. The ecology of tick-borne diseases. Int J Parasitol. 2013;43:1059–77.
Randolph SE. The shifting landscape of tick-borne zoonoses: tick-borne encephalitis and Lyme borreliosis in Europe. Philos Trans R Soc Lond B Biol Sci. 2001;356:1045–56.
Sonenshine D. Ecological Dynamics of Tick-Borne Zoonoses. Oxford: Oxford University Press; 1994.
Snow statistics - Finnish Meteorological Institute [Internet]. [cited 2016 Aug 12]. Available from: http://en.ilmatieteenlaitos.fi/snow-statistics
Sormunen J, Klemola T, Vesterinen E, Vuorinen I, Hytönen J, Hänninen J, et al. Assessing the abundance, seasonal questing activity, and Borrelia and tick-borne encephalitis virus (TBEV) prevalence of Ixodes ricinus ticks in a Lyme borreliosis endemic area in Southwest Finland. Ticks Tick Borne Dis. 2016;7:208–15.
Wilhelmsson P, Lindblom P, Fryland L, Nyman D, Jaenson TGT, Forsberg P, et al. Ixodes ricinus ticks removed from humans in northern Europe: seasonal pattern of infestation, attachment sites and duration of feeding. Parasit Vectors. 2013;6:362.
Tälleklint L, Jaenson TG. Infestation of mammals by Ixodes ricinus ticks (Acari: Ixodidae) in south-central Sweden. Exp Appl Acarol. 1997;21:755–71.
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:1058–68.
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.
Amori G, Hutterer R, Kryštufek B, Yigit N, Mitsain G, Palomo LJ, et al. Myodes glareolus. The IUCN Red List of Threatened Species 2008: e.T4973A11104409. [Internet]. 2008 [cited 2016 Nov 13]. Available from: http://dx.doi.org/10.2305/IUCN.UK.2008.RLTS.T4973A11104409.en
Gray JS, Kirstein F, Robertson JN, Stein J, Kahl O. Borrelia burgdorferi sensu lato in Ixodes ricinus ticks and rodents in a recreational park in south-western Ireland. Exp Appl Acarol. 1999;23:717–29.
Hanincová K, Schäfer SM, Etti S, Sewell HS, Taragelová V, Ziak D, et al. Association of Borrelia afzelii with rodents in Europe. Parasitology. 2003;126:11–20.
Kallio ER, Begon M, Henttonen H, Koskela E, Mappes T, Vaheri A, et al. Cyclic hantavirus epidemics in humans - Predicted by rodent host dynamics. Epidemics. 2009;1:101–7.
Hanski I, Hansson L, Henttonen H. Specialist predators, generalist predators, and the microtine rodent cycle. J Anim Ecol. 1991;60:353–67.
Massey FP, Smith MJ, Lambin X, Hartley SE. Are silica defences in grasses driving vole population cycles? Biol Lett. 2008;4:419–22.
Rosà R, Pugliese A. Effects of tick population dynamics and host densities on the persistence of tick-borne infections. Math Biosci. 2007;208:216–40.
Gern L, Estrada-Peña A, Frandsen F, Gray JS, Jaenson TG, Jongejan F, et al. European reservoir hosts of Borrelia burgdorferi sensu lato. Zentralbl Bakteriol. 1998;287:196–204.
Meerburg BG, Singleton GR, Kijlstra A. Rodent-borne diseases and their risks for public health. Crit Rev Microbiol. 2009;35:221–70.
Rais O, Gern L. Efficient transmission of Borrelia burgdorferi between co-feeding Ixodes ricinus ticks (Acari: Ixodidae). J Med Entomol. 1996;33:189–92.
Labuda M, Kozuch O, Zuffová E, Elecková E, Hails RS, Nuttall PA. Tick-borne encephalitis virus transmission between ticks cofeeding on specific immune natural rodent hosts. Virology. 1997;235:138–43.
Randolph SE. Ticks are not insects: Consequences of contrasting vector biology for transmission potential. Parasitol Today. 1998;14:186–92.
Siukkola A. Seasonality of Ixodes ricinus and Ixodes trianguliceps tick on the bank vole (Myodes glareolus) and on vegetation in Central Finland. Master’s thesis, University of Jyväskylä; 2014.
Nilsson A. Seasonal occurrence of Ixodes ricinus (Acari) in vegetation and on small mammals in southern Sweden. Ecography. 1988;11:161–5.
Kallio ER, Begon M, Henttonen H, Koskela E, Mappes T, Vaheri A, et al. Hantavirus infections in fluctuating host populations: the role of maternal antibodies. Proc R Soc Lond B. 2010;277:3783–91.
Snow KR. Identification of larval ticks found on small mammals in Britain. Berkshire: Mammal Society; 1978.
Filippova NA. Arachnida class: ixodid ticks of the subfamily Ixodinae. Fauna SSSR. Leningdrad: Nauka; 1977
Arthur DR. British ticks. ix. London: Butterworths; 1963.
Caporale DA, Rich SM, Spielman A, Telford SR, Kocher TD. Discriminating between Ixodes ticks by means of mitochondrial DNA sequences. Mol Phylogenet Evol. 1995;4:361–5.
Krebs CJ, Park T, Station RF. Demographic changes in fluctuating populations of Microtus californicus. Ecol Monogr. 1966;36:239–73.
Randolph SE, Storey K. Impact of microclimate on immature tick-rodent host interactions (Acari: Ixodidae): implications for parasite transmission. J Med Entomol. 1999;36:741–8.
Tagliapietra V, Rosà R, Arnoldi D, Cagnacci F, Capelli G, Montarsi F, et al. Saturation deficit and deer density affect questing activity and local abundance of Ixodes ricinus (Acari, Ixodidae) in Italy. Vet Parasitol Elsevier BV. 2011;183:114–24.
Paterson S, Lello J. Mixed models: Getting the best use of parasitological data. Trends Parasitol. 2003;19:370–5.
Akaike H. A New look at the statistical model identification. Autom Control IEEE Trans. 1974;19:716–23.
Kurtenbach K, Hanincová K, Tsao JI, Margos G, Fish D, Ogden NH. Fundamental processes in the evolutionary ecology of Lyme borreliosis. Nat Rev Microbiol. 2006;4:660–9.
Hamer SA, Hickling GJ, Sidge JL, Walker ED, Tsao JI. Synchronous phenology of juvenile Ixodes scapularis, vertebrate host relationships, and associated patterns of Borrelia burgdorferi ribotypes in the midwestern United States. Ticks Tick Borne Dis. 2012;3:65–74.
Dobson ADM, Finnie TJR, Randolph SE. A modified matrix model to describe the seasonal population ecology of the European tick Ixodes ricinus. J Appl Ecol. 2011;48:1017–28.
Talleklint L, Jaenson T. Seasonal variations in density of questing Ixodes ricinus (Acari: Ixodidae) nymphs and prevalence of infection with B. burgdorferi s.l. in south central Sweden. J Med Entomol. 1996;33:592–7.
Randolph SE. Tick ecology: processes and patterns behind the epidemiological risk posed by ixodid ticks as vectors. Parasitology. 2004;129(Suppl):S37–65.
Arsnoe IM, Hickling GJ, Ginsberg HS, McElreath R, Tsao JI. Different populations of blacklegged tick nymphs exhibit differences in questing behavior that have implications for human Lyme disease risk. PLoS One. 2015;10:e0127450.
Gilbert L, Aungier J, Tomkins JL. Climate of origin affects tick (Ixodes ricinus) host-seeking behavior in response to temperature: Implications for resilience to climate change? Ecol Evol. 2014;4:1186–98.
Dobson ADM. Ticks in the wrong boxes: assessing error in blanket-drag studies due to occasional sampling. Parasit Vectors. 2013;6:344.
Dobson ADM, 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:67–74.
Wilson ML, Adler GH, Spielman A. Correlation between abundance of deer and that of the deer tick, Ixodes dammini (Acari: Ixodidae). Ann Entomol Soc Am. 1985;78:172–6.
van Duijvendijk G, Sprong H, Takken W. Multi-trophic interactions driving the transmission cycle of Borrelia afzelii between Ixodes ricinus and rodents: a review. Parasit Vectors. 2015;8:643.
Kozakiewicz M, Chołuj A, Kozakiewicz A. Long-distance movements of individuals in a free-living bank vole population: an important element of male breeding strategy. Acta Theriol (Warsz). 2007;52:339–48.
Boyer N, Reale D, Marmet J, Pisanu B, Chapuis JL. Personality, space use and tick load in an introduced population of Siberian chipmunks Tamias sibiricus. J Anim Ecol. 2010;79:538–47.
Dizij A, Kurtenbach K. Clethrionomys glareolus, but not Apodemus flavicollis, acquires resistance to lxodes ricinus L., the main European vector of Borrelia burgdorferi. Parasite Immunol. 1995;17:177–83.
Wikel SK, Bergman D. Tick-host immunology: Significant advances and challenging opportunities. Parasitol Today. 1997;13:383–9.
Paulauskas A, Rosef O, Radzijevskaja J, Turcinaviciene J, Ambrasiene D. Infestation of mice and voles with Ixodes ricinus ticks in Lithuania and Norway. Est J Ecol. 2009;58:112–25.
Perkins SE, Cattadori IM, Tagliapietra V, Rizzoli AP, Hudson PJ. Empirical evidence for key hosts in persistence of a tick-borne disease. Int J Parasitol. 2003;33:909–17.
Hughes VL, Randolph SE. Testosterone depresses innate and acquired resistance to ticks in natural rodent hosts: a force for aggregated distributions of parasites. J Parasitol. 2001;87:49–54.
Mills S, Grapputo A, Jokinen I, Koskela E, Mappes T, Oksanen TA, et al. Testosterone-mediated effects on fitness related phenotypic traits and fitness. Am Nat. 2009;173:475–87.
Mills S, Grapputo A, Jokinen I, Koskela E, Mappes T, Poikonen T. Fitness trade-offs mediated by immunosuppression costs in a small mammal. Evolution. 2010;64:166–79.
Ims RA. Male spacing systems in microtine rodents. Am Nat. 1987;130:475–84.
Koskela E, Mappes T, Ylonen H. Territorial behaviour and reproductive success of bank vole Clethrionomys glareolus females. J Anim Ecol. 1997;66:341–9.
Harrison A, Bennett NC. The importance of the aggregation of ticks on small mammal hosts for the establishment and persistence of tick-borne pathogens: an investigation using the R 0 model. Parasitology. 2012;139:1605–13.
Randolph SE, Gern L, Nuttall P. Co-feeding ticks: Epidemiological significance for tick-borne pathogen transmission. Parasitol Today. 1996;12:472–9.
Randolph SE. Tick-borne encephalitis incidence in central and eastern Europe: consequences of political transition. Microbes Infect. 2008;10:209–16.
Randolph SE, Green RM, Peacey MF, Rogers DJ. Seasonal synchrony: the key to tick-borne encephalitis foci identified by satellite data. Parasitology. 2000;121:15–23.
Cotton MJ, Watts CH. The ecology of the tick Ixodes trianguliceps Birula (Arachnida; Acarina; Ixodoidea). Parasitology. 1967;57:525–31.
Bown KJ, Begon M, Bennett M, Birtles RJ, Burthe S, Lambin X, et al. Sympatric Ixodes trianguliceps and Ixodes ricinus ticks feeding on field voles (Microtus agrestis): Potential for increased risk of Anaplasma phagocytophilum in the United Kingdom? Vector-Borne Zoonotic Dis. 2006;6:404–10.
Bown KJ, Begon M, Bennett M, Woldehiwet Z, Ogden NH. Seasonal dynamics of Anaplasma phagocytophila in a rodent-tick (Ixodes trianguliceps) system, United Kingdom. Emerg Infect Dis. 2003;9:63–70.
Nefedova VV, Korenberg EI, Kovalevskii YV, Samokhvalov MV, Gorelova NB. The role of Ixodes trianguliceps tick larvae in circulation of Babesia microti in the Middle Urals. Entomol Rev. 2013;93:258–66.
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.
Kallio ER, Begon M, Birtles RJ, Bown KJ, Koskela E, Mappes T, et al. First report of Anaplasma phagocytophilum and Babesia microti in rodents in Finland. Vector Borne Zoonotic Dis. 2014;14:389–93.
Cornulier T, Yoccoz NG, Bretagnolle V, Brommer JE, Butet A, Ecke F, et al. Europe-wide dampening of population cycles in keystone herbivores. Science. 2013;340:63–6.
We thank all field and laboratory assistants involved in this project; especially Susanne Varjola, but also Meeri Väätäinen, Taru Niittynen, Tuuli-Maria Kailio, Anniina Runtuvuori, Otso Mappes, Angela Sims, Heikki Helle, Zbyszek Boratyński, Juha Ahonen, Anna Giermek, Risto Siekkinen and Sami Kyröläinen. We thank two anonymous reviewers for improving the quality of the manuscript.
This project was supported by Kone Foundation and The Academy of Finland (Eva Kallio 250524, Esa Koskela 257340 and Tapio Mappes 132190, 268670).
Availability of data and materials
The datasets analysed during the current study are available in the institutional repository of the University of Jyväskylä https://jyx.jyu.fi/dspace/handle/123456789/53330
ERK and TM designed the monitoring. CC, AS and ERK collected the field data. CC, ERK, EK and TM performed the statistical analysis. CC drafted the manuscript. ERK, EK and TM critically revised the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
The trapping methods applied in this study were approved by the Finnish Animal Experiment Board and the Finnish Ministry of the Environment, under the authorization ESAVI/3834/04.10.03/2011 and ESAVI/7256/04.10.07/2014.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Average monthly saturation deficit and temperature during the monitoring years. Figure S2. Observed mean abundance of ticks in vegetation per session, from May 2012 to October 2015. Figure S3. Mean number of vole captured per trap-night at each session and in each site, from May 2012 to October 2015. Table S1. Selection table for models explaining the abundance of ticks questing in the vegetation. Figure S4. Predicted number of larvae, nymphs and pooled nymphs and females per 100 m2 of vegetation explained by bank vole abundance. Table S2. Selection table for models explaining the abundance of ticks questing in the vegetation. Table S3. Total number of ticks (per species and stage) collected on voles. Figure S5. Vole infestation per session with I. ricinus larvae and nymphs from May 2012 to October 2015. Table S4. Selection table for models explaining the abundance of infesting larvae. Table S5. Selection table for models explaining the abundance of infesting nymphs. Table S6. Additional model for the abundance of nymphs questing in the vegetation. Table S7. Selection table for the additional model explaining the abundance of questing nymphs. (DOCX 3706 kb)