- Open Access
Landscape structure affects the prevalence and distribution of a tick-borne zoonotic pathogen
© The Author(s). 2018
- Received: 18 July 2018
- Accepted: 13 November 2018
- Published: 4 December 2018
Landscape structure can affect pathogen prevalence and persistence with consequences for human and animal health. Few studies have examined how reservoir host species traits may interact with landscape structure to alter pathogen communities and dynamics. Using a landscape of islands and mainland sites we investigated how natural landscape fragmentation affects the prevalence and persistence of the zoonotic tick-borne pathogen complex Borrelia burgdorferi (sensu lato), which causes Lyme borreliosis. We hypothesized that the prevalence of B. burgdorferi (s.l.) would be lower on islands compared to the mainland and B. afzelii, a small mammal specialist genospecies, would be more affected by isolation than bird-associated B. garinii and B. valaisiana and the generalist B. burgdorferi (sensu stricto).
Questing (host-seeking) nymphal I. ricinus ticks (n = 6567) were collected from 12 island and 6 mainland sites in 2011, 2013 and 2015 and tested for B. burgdorferi (s.l.). Deer abundance was estimated using dung transects.
The prevalence of B. burgdorferi (s.l.) was significantly higher on the mainland (2.5%, 47/1891) compared to island sites (0.9%, 44/4673) (P < 0.01). While all four genospecies of B. burgdorferi (s.l.) were detected on the mainland, bird-associated species B. garinii and B. valaisiana and the generalist genospecies B. burgdorferi (s.s.) predominated on islands.
We found that landscape structure influenced the prevalence of a zoonotic pathogen, with a lower prevalence detected among island sites compared to the mainland. This was mainly due to the significantly lower prevalence of small mammal-associated B. afzelii. Deer abundance was not related to pathogen prevalence, suggesting that the structure and dynamics of the reservoir host community underpins the observed prevalence patterns, with the higher mobility of bird hosts compared to small mammal hosts leading to a relative predominance of the bird-associated genospecies B. garinii and generalist genospecies B. burgdorferi (s.s.) on islands. In contrast, the lower prevalence of B. afzelii on islands may be due to small mammal populations there exhibiting lower densities, less immigration and stronger population fluctuations. This study suggests that landscape fragmentation can influence the prevalence of a zoonotic pathogen, dependent on the biology of the reservoir host.
- Habitat fragmentation
- Pathogen persistence
- Host community
- Borrelia burgdorferi (sensu lato)
Landscape structure can influence the persistence of a pathogen by affecting the movement of hosts and vectors, and thus the potential for transmission. As host populations become reduced within smaller habitat patches, stochastic fadeout of transmission is expected to become more common. Pathogen persistence will then depend on the degree of isolation of habitat patches, and a balance between colonization and extinction events [1, 2]. Different biological processes are predicted to affect colonization and extinction events. For example, colonization of a pathogen into new habitat patches is likely to be affected by the mobility and ecology of the host species. In contrast, host population dynamics, particularly large seasonal or interannual fluctuations, can contribute to stochastic fadeout of a pathogen within patches .
Environmental change can affect disease risk to humans or livestock, with vector-borne pathogens especially sensitive to changes in habitat and climate . Following historic reduction and fragmentation of woodlands in Europe, there are now widespread efforts to increase the amount of woodland, as well as connectivity between patches, to promote biodiversity and restore ecological processes . This has the potential to affect the distribution and prevalence of vector-borne pathogens associated with forest ecosystems. In general, woodlands are favourable habitats for the tick vector Ixodes ricinus, as they support higher densities of hosts that can provide blood meals for the tick compared to open habitats, and the ground leaf litter layer provides humid conditions for off-host tick development and survival . Ixodes ricinus is a generalist tick vector that transmits many pathogens of importance for human and animal health, including the bacteria of the Borrelia burgdorferi (sensu lato) (s.l.) complex, which are the cause of Lyme borreliosis, an emerging disease in the northern hemisphere [7, 8].
Host community composition is an important factor determining persistence of tick-borne pathogens including B. burgdorferi (s.l.) and one mechanism by which environmental change can affect disease risk [9–12]. While many vertebrate species provide blood meals for the tick I. ricinus, not all species are competent to transmit B. burgdorferi (s.l.) At least nine genospecies of B. burgdorferi (s.l.) with differing pathogenicities, clinical signs and reservoir host associations are known to be circulating in Europe [13, 14]. Three of these together, B. afzelii which is adapted to small mammals , B. garinii which is adapted to bird reservoir hosts  and B. burgdorferi (s.s.), a generalist genospecies, are collectively responsible for most cases of human disease . Deer are generally considered non-competent hosts for genospecies of B. burgdorferi (s.l.) and feed all life-stages of the tick I. ricinus [18–20]. As deer feed significant numbers of adult female ticks which go onto produce immature larval stages they are sometimes termed ‘tick reproduction hosts’ . Deer density is frequently positively linked with nymphal tick density, e.g. [22, 23]; however no significant relationship has also been found in other studies . As incompetent hosts for B. burgdorferi (s.l.) and as hosts of ticks, deer can affect transmission, with both amplification and dilution of transmission predicted from theoretical models depending on how tick vectors are distributed amongst the host community .
Little is known about the effects of landscape structure and fragmentation on the risk of Lyme borreliosis in Europe. Previous studies suggested that key small mammal reservoir hosts such as bank voles (Myodes glareolus) may decrease in density in fragmented woodland [26, 27]. This could reduce persistence of small mammal-associated genospecies of B. burgdorferi (s.l.) within isolated patches of habitat due to stochastic extinction events. Also, as birds are unable to efficiently transmit small mammal-associated B. afzelii, this limits the reintroduction of this pathogen among suitable areas of habitat in fragmented landscapes . Therefore, we expect landscape fragmentation to have different effects on the prevalence and persistence of the genospecies of B. burgdorferi (s.l.) in Europe, with small mammal-associated B. afzelii being less likely to persist and present at lower prevalence in fragmented landscapes compared to bird-associated and generalist genospecies.
Here we use a natural island system located in a large freshwater lake in Scotland as a model to investigate the effect of landscape fragmentation on B. burgdorferi (s.l.) prevalence. Four genospecies of B. burgdorferi (s.l.) are known to occur in Scotland, B. afzelii, B. garinii and B. burgdorferi (s.s.) are present and also B. valaisiana which is associated with birds and is uncommonly associated with human disease. The two species of small mammal which are widely distributed and present at highest abundance in woodland habitats are the bank vole and the wood mouse (Apodemus sylvaticus). Other mammal species that may be significant reservoir hosts in the area include the grey squirrel (Sciurus carolinensis) , the red squirrel (Sciurus vulgaris)  and the common (Sorex araneus) and pygmy shrews (Sorex minutus) . Passerine birds are also hosts for ticks in UK woodlands , with particularly high tick burdens on ground foraging species such as the common blackbird (Turdus merula) and the song thrush (Turdus philomelos). Fallow deer (Dama dama) are present in the study area.
We first hypothesized that islands would be less likely to support endemic circulation of genospecies relying on small mammal hosts compared to continuous mainland habitat. We predicted that islands would have a lower prevalence of B. afzelii, due to existing evidence suggesting that there are lower densities of small mammals in fragmented habitat [26, 27] and that islands experience local extinction of B. afzelii more frequently. In contrast, no such restriction is expected for bird specialist and generalist genospecies since birds traveling between the mainland and islands can reintroduce infected ticks following local extinctions. Secondly, we hypothesized that there would be decreasing likelihood of B. burgdorferi (s.l.) presence with increasing isolation from the mainland, due to less frequent movement between the mainland and isolated habitat patches and reduced opportunities for recolonization. According to this, smaller and more distant islands are predicted to be less likely to support pathogen persistence due to lower reservoir host density and increased effects of stochastic events. From this we predict that (i) islands overall have a lower B. burgdorferi (s.l.) prevalence than mainland sites; (ii) B. garinii would be more prevalent than B. afzelii on islands due to greater dispersal by birds compared to small mammals; and (iii) smaller more isolated islands would have a lower prevalence than larger islands located closer to the mainland.
Study sites on Loch Lomond mainland and islands, with vegetation, two indices of deer density based on deer counts and deer dung, and island size and distance to the mainland
Mean deer density (per km2)b
Deer density (dung transects)c
Island area (km2)
Distance to the mainland (km)
Scots pine, oak and birch, understorey of grass, moss and bracken.
Oak with understorey of dense areas of bramble and bracken.
Oak with mixed understory of blueberry, grass and moss
Oak with birch and holly. Understorey relatively open with blueberry, bramble and grasses.
Oak with understorey of dense areas of bramble and bracken.
Oak/managed pine with understorey of bramble, blueberry and grasses
Scots pine, birch and oak. Understorey of rhododendrons, bracken and grass species
Oak, birch, holly and alder. Understorey dense with bracken and grass species.
Oak, relatively open understorey with blueberry and bracken.
Oak, grass and bracken understorey.
Oak with blueberry and moss understorey
Birch, oak and alder. Understoreys relatively open with blueberry, bracken and moss species.
Oak, with a sparse understorey of bracken and bramble.
Oak, birch, and yew trees. Understorey dense with bracken.
Birch and alder. Understorey of rhododendrons, blueberry, gorse and bog myrtle.
Oak, birch and Scots pine. An understorey of blueberry, bracken and bramble.
Oak with blueberry and moss understorey
Oak and birch. Understorey dense with bluberry and bracken.
Field collection of ticks
Questing ticks were collected from the vegetation by blanket dragging during the peak tick questing period, between April and August. Blanket dragging was carried out between 10:00 and 16:00 h, when the vegetation was dry or at least two hours after rainfall. Ticks were collected by dragging a white blanket slowly across the surface of the vegetation until 200 nymphs were collected, or for a minimum of 3 h per site in order to collect a representative sample of ticks to estimate the prevalence of B. burgdorferi (s.l.). Nymphs attached to the blanket were counted and stored in 70% ethanol after each drag. Standardized drags to estimate nymph density were also carried out, the methodology and results are provided in Additional file 1: Text S2 and Table S3. We focussed on nymphs as they are considered to be the most important life-stage in B. burgdorferi (s.l.) transmission to humans .
Host and environmental predictors of tick abundance and B. burgdorferi (s.l.) prevalence
Both island and mainland sites were known prior to the study to have populations of I. ricinus and woodland mammals and birds typical of the area  (Additional file 1: Table S1). While it was not practical to estimate the density of all potential bird and small mammal reservoir hosts across the multiple sites and years of this study, deer density was estimated. As both amplification and dilution effects of deer density on B. burgdorferi (s.l.) prevalence have been predicted from theoretical models , this was considered an important parameter to include in statistical models to investigate the effect of landscape fragmentation on B. burgdorferi (s.l.) prevalence . Deer density was quantified on both island and mainland sites in 2015 using a deer dung transect method (, Additional file 1: Text S1, Table S2). There was no significant difference between the mean deer density estimated on island and mainland sites (Kruskal-Wallis χ2 = 0.11, df = 1, P = 0.74). Data from two deer count surveys of the Loch Lomond Islands, carried out by park staff in 2008 and 2012 (Additional file 1: Text S1, Table S2) indicated that deer density remained consistent through time (Spearman’s rank correlation coefficient = 0.65, P = 0.02). There was also no difference between the mean deer density estimated in the deer surveys in 2008, 2012 and the deer dung survey in 2015 (Kruskal-Wallis χ2 = 0.78, df = 2, P = 0.68). Therefore, there was strong reason to assume that densities estimated from the deer dung transect survey as the only method for which data from both islands and mainland were available, could reasonably be used as a measure of deer density across all years of this study.
DNA extraction, B. burgdorferi (s.l.) detection and genospecies determination
Tick DNA extraction, PCR preparation and genospecies determination using PCR products were carried out in separate laboratories. Each nymph was placed in a separate Eppendorf tube and DNA extracted using a basic ammonia extraction technique . One DNA extraction control which included reagents but no tick was included with every 11 samples. Detection of B. burgdorferi (s.l.) was either by a nested PCR  which targets the 5S-23S intergenic spacer (IGS) region or by a real time PCR  which detects a fragment of the 23S ribosomal RNA gene. The PCR protocol and cycling conditions were followed as previously described . A negative and a positive control were included in each PCR run. Genospecies determination was carried out on all PCR-positive nymphs, either by reverse line blotting as previously described [37, 38] and/or by sequencing the PCR product from the 5S-23S intergenic spacer region . Where sequencing was performed, PCR products were sequenced with the forward and reverse primers using Sanger sequencing. Sequences were trimmed in Geneious version 7.0.6 (Biomatters Ltd). To determine the genospecies each sequenced and trimmed IGS PCR product was subjected to a BLAST search against the National Centre for Biotechnology Nucleotide BLAST database. Sequences were also examined for polymorphisms characteristic of each genospecies as described previously .
All GLMMs were fitted using the lme4 package , starting with a maximal global model, which included all fixed effects and interaction terms where specified. Model selection was performed by backward stepwise elimination, based on minimising the Akaike information criterion (AICc) [43, 44]. The AICcmodavg package  was used to calculate the second-order Akaike’s information criterion (AICc) to account for small sample size. The delta AICc, defined as the rise in AICc after removing each variable from the selected model was calculated for each explanatory variable in turn. We gauged the explanatory power of the models by partitioning the total variance into three components: variation explained by the fixed effects; unexplained variation between sites; and unexplained variation between observations/drags within sites. The proportion of the sum of these three variances explained by the fixed effects is defined as the latent-scale marginal R2, R2GLMM(m)*. The 95% confidence intervals (CI) for proportions were calculated using the prop.test( ) function in ‘R’.
Testing for association between habitat fragmentation and B. burgdorferi (s.l.) prevalence
B. burgdorferi (s.l.) infection in nymphs (infected or uninfected) was modelled with a logit-binomial general linear mixed model (GLMM) as a function of the following explanatory variables: fragmentation class (island or mainland), year (2011, 2013 or 2015), deer density (estimated from dung transects) and an interaction between fragmentation class (island/mainland) and year. To control for repeated sampling of sites and overdispersion, random effects of sampling site and observation  were included. This model was repeated with the individual genospecies B. afzelii, B. garinii and B. burgdorferi (s.s.) infection in nymphs (infected or not infected) as the outcome variable. There were insufficient data to repeat the model with B. valaisiana. To test if effects of landscape structure on B. burgdorferi (s.l.) prevalence were a result of differences in B. afzelii prevalence, the model was repeated with infected nymphs of the three genospecies: B. garinii, B. burgdorferi (s.s.) and B. valaisiana as the outcome variable. As deer density was not maintained in the models during model selection, they were re-run with data from all site-year combinations (n = 42 compared to n = 40, as one site did not have a deer density estimate).
Testing for associations between island size, island isolation and B. burgdorferi (s.l.) prevalence
Data collected from all island sites in 2011, 2013 and 2015 were used for this model. Borrelia burgdorferi (s.l.) infection in nymphs (infected or uninfected) was modelled with a logit-binomial GLMM as a function of the following explanatory variables: island size, distance to the mainland, deer density and year with random effects of site and observation. There were insufficient data to model individual genospecies separately.
There was variation in genospecies distribution over the study period, with many apparent introductions and losses of genospecies from individual sites (Fig. 2). Whereas 10 out of 11 mainland site collections included infective nymphs, this was only true for 13 out of 31 collections on the islands. When comparing only collections in which infected nymphs were detected, the prevalence on island and mainland sites was comparable (2.4% and 2.75%). Using our overall prevalence estimate for the mainland as the expected prevalence, we calculated the probability of not detecting any infected ticks in our sample for each of the island site collections. For the majority of island site collections without B. burgdorferi (s.l.) detection (14 out of 18, Additional file 1: Table S3), the probability of not detecting the pathogen if it was present was low (< 0.03). For the remaining 4 collections, the probability exceeded 0.05, indicating that absence of the pathogen could not be inferred with confidence for these collections due to limited sample size.
Borrelia burgdorferi (s.l.) prevalence on island and mainland sites
Best-fit models explaining variation in Borrelia burgdorferi (s.l.), B. afzelii and B. garinii prevalence in questing nymphs. Data from all sites and years are included, no measure of deer abundance is included. Results for each of the best-fit models are shown in each case. Delta AICc indicates the change in AICc after removing each variable from the best-fit model
B. burgdorferi (s.l.) prevalence
< 2 × 10-16
Location (Mainland vs Island)
Year (2013 vs 2011)
Year (2015 vs 2011)
Model to predict B. afzelii prevalence
1.2 × 10-13
Location (Mainland vs Island)
Model to predict B. garinii prevalence
< 2 × 10-16
Model to predict B. burgdorferi (s.s.) prevalence
4.2 × 10-9
Model to predict combined B. burgdorferi (s.s.), B. garinii and B. valaisiana prevalence
< 2 × 10-16
When considering the effect of island size, distance to the mainland, deer density and year of sampling on B. burgdorferi (s.l.) prevalence, the best-fit model was the intercept only model.
Borrelia burgdorferi (s.l.) genospecies prevalence on island and mainland sites
The best-fit model for B. afzelii prevalence included a fixed effect of fragmentation class (island/mainland) (Table 2). Mainland sites were significantly more likely to harbour B. afzelii infected nymphs compared to islands, OR = 22.5 (95% CI: 2.3–238.2) (delta AICc = 5.0). The proportion of variance explained by the fixed effect R2GLMM(m)* was 54% with all the remaining variation accounted for by the observation level random effect. The best-fit model for B. garinii and B. burgdorferi (s.s.) prevalence did not include a significant effect of year or fragmentation class (island/mainland) (Table 2). The best-fit model for B. burgdorferi (s.l.) prevalence with B. afzelii infected nymphs removed (combined B. garinii, B. burgdorferi (s.s.) and B. valaisiana prevalence) also did not include a significant effect of year or fragmentation class (island/mainland).
This study investigated whether landscape fragmentation (islands versus mainland) affects the distribution and persistence of the tick-borne pathogen B. burgdorferi (s.l.) and three genospecies B. afzelii, B. garinii and B. burgdorferi (s.s.). Consistent with our prediction for fragmented landscapes, we found a significantly lower prevalence of B. burgdorferi (s.l.) among island sites compared to mainland sites, due to an absence of B. burgdorferi (s.l.) from over half of the sample collections at island sites. The prevalence of B. burgdorferi (s.l.) on islands with infected ticks was similar to that on the mainland, suggesting that islands were capable of supporting similar levels of transmission as the mainland. For most island site collections, the sampling effort was sufficient to detect the pathogen with reasonable confidence, had it been present. This suggests local extinctions could contribute to the high proportion of apparent absences and the low overall prevalence of B. burgdorferi (s.l.) on the islands. In general, only the bird-associated genospecies B. garinii and B. valaisiana or the generalist genospecies B. burgdorferi (s.s.) were found on islands, though their detection there was also more sporadic than on mainland sites.
As expected, the prevalence and frequency of detection of small mammal-associated B. afzelii was significantly lower in fragmented habitat. There are several ecological factors which may affect the persistence of B. afzelii on islands. A previous study reported lower densities of competent small mammal reservoir hosts on the islands compared to the mainland . Low rodent densities might reduce the proportion of blood meals on competent versus non-competent transmission hosts below a threshold necessary for pathogen persistence . In addition, seasonal and inter-annual population fluctuations, which are common in small mammal populations, could lead to stochastic extinction of the pathogen. While these effects might also affect B. afzelii dynamics on mainland sites, their isolated nature makes island sites less likely to be recolonised by infected rodent hosts.
In contrast, the patterns for bird-associated B. garinii and B. valaisiana and generalist genospecies B. burgdorferi (s.s.) suggested more frequent introductions to islands, likely reflecting the mobility of the bird hosts of these genospecies. However, the high mobility of bird hosts did not necessarily result in pathogen persistence on islands, and despite recording similar communities of competent bird reservoir hosts (Additional file 1: Table S1), bird-associated or generalist genospecies of B. burgdorferi (s.l.) were not detected on three islands (CA, MO, TO) and were absent from other islands during some years (Fig. 2).
Questing nymph densities were not significantly different between island and mainland sites (Additional file 1: Text S2, Table S4) and were, in general, equivalent or higher than elsewhere in Scotland where B. burgdorferi (s.l.) transmission has been detected (e.g. 0.6–11.5 nymphs/10 m2, [24, 47]). Previous work in our study system had found higher levels of tick infestation on rodents on islands compared to mainland sites . Therefore, vector abundance was not the limiting factor for B. burgdorferi (s.l.) presence on islands . Likewise deer density which had been found to influence B. burgdorferi (s.l.) prevalence in some studies (e.g. [25, 47]) had no detectable effect in our study (Table 2). Instead, we suggest that the observed variability in B. burgdorferi (s.l.) prevalence and distribution among islands and between island and mainland sites is most likely to be due to effects of habitat fragmentation on reservoir host movements and community composition with resulting effects on pathogen persistence.
Contrary to our predictions, and results of studies on other pathogen systems [49, 50], we did not find an association between island size, or spatial isolation and B. burgdorferi (s.l.) prevalence. Within the spatial scale of our study, these predictors might be less relevant to the persistence of bird-associated genospecies that were found to dominate among island sites. In addition, the low prevalence of B. burgdorferi (s.l.) among island sites may have resulted in reduced power to detect these effects.
More generally, the low overall prevalence of B. burgdorferi (s.l.) in our system may have affected our power to detect some biological effects. For some islands, the numbers of ticks available for testing were insufficient to discriminate with confidence between B. burgdorferi (s.l.) absence and non-detection for any given year and location. Moreover, we cannot rule out persistence on islands at a lower prevalence than that found on the mainland, which could lead us to overestimate the frequency of B. burgdorferi (s.l.) extinction from islands. However, these considerations would not alter our main findings regarding the association between landscape fragmentation and decreased B. burgdorferi (s.l.) prevalence or the effect of fragmentation on the distribution and prevalence of B. afzelii.
Woodland cover in the UK has been reduced over the millennia from historic levels to among the lowest proportion in Europe . As a result, in many areas, suitable woodland habitat for I. ricinus and vertebrate hosts is fragmented into patches separated by agricultural land or islands of green space within urban areas. It is possible that processes similar to those discussed above affect the persistence of B. burgdorferi (s.l.) within these fragmented mainland areas. However, a difference from our island study system and fragmented woodland surrounded by land is that ecotonal habitat surrounding fragments of woodland may be an important habitat for reservoir hosts in the latter . For example, a recent study has found that ecotonal habitat and edge density affect B. burgdorferi (s.l.) prevalence in European forest fragments . Surveys from England so far suggest that the dominant genospecies present are the bird-associated genospecies B. garinii and B. valaisiana with B. afzelii and B. burgdorferi (s.s.) being apparently absent from many sites [52–54]. Our results suggest that local extinction and limited (re)colonization opportunities for B. afzelii might help to explain these patterns. Policies to increase woodland area and contiguity could potentially increase the prevalence of B. burgdorferi (s.l.) by facilitating movements of reservoir host species and tick vectors between patches of suitable habitat.
Our results contrast with those from eastern North America, where fragmented woodland habitats have been associated with increased prevalence of B. burgdorferi [55, 56] likely due to increased densities of small mammals . However, in North America, birds are still considered important in transporting infected ticks with resulting invasion of the pathogen in new areas . The presence of a single generalist genospecies, B. burgdorferi (s.s.) in North America with a broad host niche that includes birds and small mammals may facilitate colonisation and persistence in fragmented habitat . Our findings also disagree those of a European study which found an increased prevalence of B. burgdorferi (s.l.) in fragmented woodland but did not identify the genospecies . Further studies are required to clarify the effect of habitat fragmentation on the distribution and abundance of B. burgdorferi (s.l.) in Europe and to identify commonalities and differences to North American systems.
We found significantly lower prevalence of B. burgdorferi (s.l.) among island sites compared to the nearby mainland in a naturally fragmented landscape. Effects of landscape structure on pathogen prevalence was not associated with either nymph or deer abundance and was mainly a result of a significantly lower detection of pathogen presence on island sites and a lower prevalence of small mammal-associated B. afzelii on islands in particular. Host biology, in particular the higher mobility and less extreme population fluctuations of bird hosts compared to small mammals, is likely to play a role in explaining the wider distribution and greater prevalence of the bird-associated genospecies of B. garinii and the generalist genospecies B. burgdorferi (s.s.). Seasonal population fluctuations of small mammal hosts may contribute to pathogen extinction, as has been seen in other disease systems , while limited immigration reduces opportunities for colonisation events. Pathogens maintained by highly mobile hosts such as birds which can migrate between habitat patches are predicted to be more successful in fragmented habitats.
The authors thank the staff at SCENE for assistance with access to field sites in this study.
This work was supported by a Biotechnology and Biological Sciences Research Council (BBSRC) Doctoral Training grant (BB/F016786/1) to CM and by a Natural Environment Research Council grant NE/G000786/1 to RBk. ED was supported by a bursary from the Scottish Centre for Ecology and the Natural Environment (SCENE).
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional file.
CM participated in the study design, data collection, statistical analysis and drafted the manuscript. ED participated in the study design, data analysis, laboratory work and field sampling. PI, MJ, VP, EK, FT, AW, RBi carried out field and/or laboratory work. LG provided comments on the study design and manuscript. PJ advised on the statistical analysis. RBk participated in the study design and provided comments on drafts of the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Grenfell BT, Harwood J. (Meta)population dynamics of infectious diseases. Trends Ecol Evol. 1997;12:395–9.View ArticlePubMedGoogle Scholar
- Lloyd-Smith JO, Cross PC, Briggs CJ, Daugherty M, Getz WM, Latto J, et al. Should we expect population thresholds for wildlife disease? Trends Ecol Evol. 2005;20:511–9.View ArticlePubMedGoogle Scholar
- Guzzetta G, Tagliapietra V, Perkins SE, Hauffe HC, Poletti P, Merler S, et al. Population dynamics of wild rodents induce stochastic fadeouts of a zoonotic pathogen. J Anim Ecol. 2017;86:451–9.View ArticlePubMedGoogle Scholar
- Kilpatrick AM, Randolph SE. Drivers, dynamics, and control of emerging vector-borne zoonotic diseases. Lancet. 2012;380:1946–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Quine CP, Cahalan C, Hester A, Humphrey J, Kirby K, Moffat A, et al. Woodlands. In: UK National Ecosystem Assessment: Technical Report. Cambridge: UNEP-WCMC; 2011. p. 241–95.Google Scholar
- Medlock JM, Shuttleworth H, Copley V, Hansford K, Leach S. Woodland biodiversity management as a tool for reducing human exposure to Ixodes ricinus ticks: a preliminary study in an English woodland. J Vector Ecol. 2012;37:307–15.View ArticlePubMedGoogle Scholar
- Dennis D, Hayes E. Epidemiology of Lyme Borreliosis. In: Gray J, Kahl O, Lane R, Stanek G, editors. Lyme Borreliosis, Biology, Epidemiology and Control. New York: CABI; 2002. p. 251–80.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Gilbert L, Norman R, Laurenson KM, Reid HW, Hudson PJ. Disease persistence and apparent competition in a three-host community: an empirical and analytical study of large-scale, wild populations. J Anim Ecol. 2001;70:1053–61.View ArticleGoogle Scholar
- LoGiudice K, Ostfeld RS, Schmidt K a, Keesing F. The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. Proc Natl Acad Sci U S A. 2003;100:567–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosà R, Pugliese A, Norman R, Hudson PJ. Thresholds for disease persistence in models for tick-borne infections including non-viraemic transmission, extended feeding and tick aggregation. J Theor Biol. 2003;224:359–76.View ArticlePubMedGoogle Scholar
- Hartemink N, Randolph SE, Davis S, Heesterbeek JAP. The basic reproduction number for complex disease systems: defining R0 for tick-borne infections. Am Nat. 2008;171:743–54.View ArticlePubMedGoogle Scholar
- Rauter C, Hartung T. Prevalence of Borrelia burgdorferi sensu lato genospecies in Ixodes ricinus ticks in Europe: a metaanalysis. Appl Environ Microbiol. 2005;71:7203–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Strnad M, Honig V, Ruzek D, Grubhoffer L, Rego ROM. Europe-wide meta-analysis of Borrelia burgdorferi sensu lato prevalence in questing Ixodes ricinus ticks. Appl Environ Microbiol. 2017;83:e00609-17.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- Taragel’ová V, Koči J, Hanincová K, Kurtenbach K, Derdáková M, Ogden NH, et al. Blackbirds and song thrushes constitute a key reservoir of Borrelia garinii, the causative agent of borreliosis in central Europe. Appl Environ Microbiol. 2008;74:1289–93.View ArticlePubMedGoogle Scholar
- Stanek G, Reiter M. The expanding Lyme Borrelia complex-clinical significance of genomic species? Clin Microbiol Infect. 2011;17:487–93.View ArticlePubMedGoogle Scholar
- Telford SR, Mather TN, Moore SI, Wilson ML, Spielman A. Incompetence of deer as reservoirs of the Lyme disease spirochete. Am J Trop Med Hyg. 1988;39:105–9.View ArticlePubMedGoogle Scholar
- Jaenson TGT, Tälleklint L. Incompetence of roe deer (Capreolus capreolus) as reservoirs of the Lyme disease spirochete. J Med Entomol. 1992;29:813–7.View ArticlePubMedGoogle Scholar
- Matuschka FR, Heiler M, Eiffert H, Fischer P, Lotter H, Spielman A. Diversionary role of hoofed game in the transmission of Lyme disease spirochetes. Am J Trop Med Hyg. 1993;48:693–9.View ArticlePubMedGoogle Scholar
- Gray JS, Robertson JN, Key S. Limited role of rodents as reservoirs of Borrelia burgdorferi sensu lato in Ireland. Eur J Epidemiol. 2000;16:101–3.View ArticlePubMedGoogle Scholar
- Gilbert L. Altitudinal patterns of tick and host abundance: a potential role for climate change in regulating tick-borne diseases? Oecologia. 2010;162:217–25.View ArticlePubMedGoogle Scholar
- James MC, Gilbert L, Bowman AS, Forbes KJ. The heterogeneity, distribution, and environmental associations of Borrelia burgdorferi sensu lato, the agent of Lyme borreliosis, in Scotland. Front Public Health. 2014;2:129.View ArticlePubMedPubMed CentralGoogle Scholar
- Millins C, Gilbert L, Johnson P, James M, Kilbride E, Birtles R, et al. Spatial and temporal heterogeneity in the abundance and distribution of Ixodes ricinus and Borrelia burgdorferi (sensu lato) in Scotland: implications for risk prediction. Parasit Vectors. 2016;9:595.View ArticlePubMedPubMed CentralGoogle Scholar
- Mannelli A, Bertolotti L, Gern L, Gray J. Ecology of Borrelia burgdorferi sensu lato in Europe: transmission dynamics in multi-host systems, influence of molecular processes and effects of climate change. FEMS Microbiol Rev. 2012;36:837–61.View ArticlePubMedGoogle Scholar
- Rajska-Jurgiel E. Demography of woodland rodents in fragmented habitat. Environ Microbiol. 1992;37:73–90.Google Scholar
- Paterson VL. Population dynamics of rodents and their parasite communities in a naturally fragmented landscape. PhD Thesis, University of Glasgow; 2012.Google Scholar
- Heylen D, Matthysen E, Fonville M, Sprong H. Songbirds as general transmitters but selective amplifiers of Borrelia burgdorferi sensu lato genotypes in Ixodes ricinus ticks. Environ Microbiol. 2014;16:2859–68.View ArticlePubMedGoogle Scholar
- Craine NG, Nuttall PA, Marriott AC, Randolph SE. Role of grey squirrels and pheasants in the transmission of Borrelia burgdorferi sensu lato, the Lyme disease spirochaete, in the UK. Folia Parasitol. 1997;44:155–60.PubMedGoogle Scholar
- Gern L, Estrada-Peña A, Frandsen F, Gray JS, Jaenson TGT, Jongejan F, et al. European reservoir hosts of Borrelia burgdorferi sensu lato. Zentralblatt für Bakteriol. 1998;287:196–204.View ArticleGoogle Scholar
- James MC, Furness RW, Bowman AS, Forbes KJ, Gilbert L. The importance of passerine birds as tick hosts and in the transmission of Borrelia burgdorferi, the agent of Lyme disease: a case study from Scotland. Ibis. 2011;153:293–302.View ArticleGoogle Scholar
- Robertson JN, Gray JS, Stewart P. Tick bite and Lyme borreliosis risk at a recreational site in England. Eur J Epidemiol. 2000;16:647–52.View ArticlePubMedGoogle Scholar
- Marques FFC, Buckland ST, Goffin D, Dixon CE, Borchers DL, Mayle BA, et al. Estimating deer abundance from line transect surveys of dung: sika deer in southern Scotland. J Appl Ecol. 2001;38:349–63.View ArticleGoogle Scholar
- Gern L, Douet V, López Z, Rais O, Cadenas FM. Diversity of Borrelia genospecies in Ixodes ricinus ticks in a Lyme borreliosis endemic area in Switzerland identified by using new probes for reverse line blotting. Ticks Tick Borne Dis. 2010;1:23–9.View ArticlePubMedGoogle Scholar
- Rijpkema SG, Molkenboer MJ, Schouls LM, Jongejan F, Schellekens J. Simultaneous detection and genotyping of three genomic groups of Borrelia burgdorferi sensu lato in Dutch Ixodes ricinus ticks by characterization of the amplified intergenic spacer region between 5S and 23S rRNA genes. J Clin Microbiol. 1995;33:3091–5.PubMedPubMed CentralGoogle Scholar
- Courtney JW, Kostelnik LM, Zeidner NS, Massung RF. Multiplex real-time PCR for detection of Anaplasma phagocytophilum and Borrelia burgdorferi. J Clin Microbiol. 2004;42:3164–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Alekseev AN, Dubinina HV, Van De PI, Schouls LM. Identification of Ehrlichia spp. and Borrelia burgdorferi in Ixodes ticks in the Baltic Regions of Russia. J Clin Microbiol. 2001;39:2237–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Poupon M, Lommano E, Douet V, Rais O, Schaad M, Jenni L, et al. Prevalence of Borrelia burgdorferi sensu lato in ticks collected from migratory birds in Switzerland. Appl Environ Microbiol. 2006;72:976–9.View ArticlePubMed CentralGoogle Scholar
- Rijpkema SG, Herbes RG, Verbeek-De Kruif N, Schellekens JF. Detection of four species of Borrelia burgdorferi sensu lato in Ixodes ricinus ticks collected from roe deer (Capreolus capreolus) in the Netherlands. Epidemiol Infect. 1996;117:563–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Millins C, Magierecka A, Gilbert L, Edoff A, Brereton A, Kilbride E, et al. An invasive mammal (gray squirrel, Sciurus carolinensis) commonly hosts diverse and atypical genotypes of the zoonotic pathogen Borrelia burgdorferi sensu lato. Appl Environ Microbiol. 2015;81:109–15.View ArticleGoogle Scholar
- R Development Core Team. R version 3.4.1: Available at: https://www.r-project.org
- Bates D, Maechler M, Bolker B, Walker S. lme4: Linear mixed-effects , models using Eigen and S4. 2014. R package version 1.1-6. http://cran.r-project.org/package=lme4.
- Mazerolle M. Improving data analysis in herpetology: using Akaike’s Information Criterion (AIC) to assess the strength of biological hypotheses. Amphibia-Reptilia. 2006;27:169–80.View ArticleGoogle Scholar
- Mazerolle MJ. Model Selection and Multimodel Inference Based on (Q)AIC(c). 2017. R Package version 1.1-1. https://cran.r-project.org/web/packages/AICcmodavg/AICcmodavg.pdf
- Nakagawa S, Johnson PCD, Schielzeth H. Coefficient of determination R2 and intra-class correlation coefficient ICC from generalized linear mixed-effects models. J R Soc Interface. 2017;14:20170213.View ArticlePubMedPubMed CentralGoogle Scholar
- Harrison XA. A comparison of observation-level random effect and beta-binomial models for modelling overdispersion in binomial data in ecology & evolution. PeerJ. 2015;3:e1114.View ArticlePubMedPubMed CentralGoogle Scholar
- James MC, Bowman AS, Forbes KJ, Lewis F, McLeod JE, Gilbert L. Environmental determinants of Ixodes ricinus ticks and the incidence of Borrelia burgdorferi sensu lato, the agent of Lyme borreliosis, in Scotland. Parasitology. 2013;140:237–46.View ArticlePubMedGoogle Scholar
- Norman R, Bowers RG, Begon M, Hudson PJ. Persistence of tick-borne virus in the presence of multiple host species: tick reservoirs and parasite-mediated competition. J Theor Biol. 1999;200:111–8.View ArticlePubMedGoogle Scholar
- Lindström KM, Foufopoulos J, Pärn H, Wikelski M. Immunological investments reflect parasite abundance in island populations of Darwin’s finches. Proc R Soc B Biol Sci. 2004;271:1513–9.View ArticleGoogle Scholar
- Spurgin GL, Illera JC, Padilla DP, Richardson DS. Biogeographical patterns and co-occurrence of pathogenic infection across island populations of Berthelot’s pipit (Anthus berthelotii). Oecologia. 2012;168:691–701.View ArticlePubMedGoogle Scholar
- Ehrmann S, Ruyts SC, Scherer-Lorenzen M, Bauhus J, Brunet J, Cousins SAO, et al. Habitat properties are key drivers of Borrelia burgdorferi (s.l.) prevalence in Ixodes ricinus populations of deciduous forest fragments. Parasit Vectors. 2018;11:23.View ArticlePubMedPubMed CentralGoogle Scholar
- Bettridge J, Renard M, Zhao F, Bown KJ, Birtles RJ. Distribution of Borrelia burgdorferi sensu lato in Ixodes ricinus populations across central Britain. Vector-Borne Zoonotic Dis. 2013;13:139–46.View ArticlePubMedGoogle Scholar
- Vollmer SA, Bormane A, Dinnis RE, Seelig F, Dobson ADM, Aanensen DM, et al. Host migration impacts on the phylogeography of Lyme borreliosis spirochaete species in Europe. Environ Microbiol. 2011;13:184–92.View ArticlePubMedGoogle Scholar
- Hansford KM, Fonville M, Jahfari S, Sprong H, Medlock JM. Borrelia miyamotoi in host-seeking Ixodes ricinus ticks in England. Epidemiol Infect. 2014;143:1079–87.View ArticlePubMedGoogle Scholar
- Allan B, Keesing F, Ostfeld R. Effect of forest fragmentation on Lyme disease risk. Conserv Biol. 2003;17:267.View ArticleGoogle Scholar
- Brownstein JS, Skelly DK, Holford TR, Fish D. Forest fragmentation predicts local scale heterogeneity of Lyme disease risk. Oecologia. 2005;146:469–75.View ArticlePubMedGoogle Scholar
- Ostfeld R, Keesing F. Biodiversity and disease risk: the case of Lyme disease. Conserv Biol. 2000;14:722–8.View ArticleGoogle Scholar
- Ogden NH, Lindsay LR, Leighton PA. Predicting the rate of invasion of the agent of Lyme disease Borrelia burgdorferi. J Appl Ecol. 2013;50:510–8.View ArticleGoogle Scholar
- Hanincova K, Kurtenbach K, Diuk-Wasser M, Brei B, Fish D. Epidemic spread of Lyme borreliosis, northeastern United States. Emerg Infect Dis. 2006;12:604–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Halos L, Bord S, Cotté V, Gasqui P, Abrial D, Barnouin J, et al. Ecological factors characterizing the prevalence of bacterial tick-borne pathogens in Ixodes ricinus ticks in pastures and woodlands. Appl Environ Microbiol. 2010;76:4413–20.View ArticlePubMedPubMed CentralGoogle Scholar