Landscape structure affects the prevalence and distribution of a tick-borne zoonotic pathogen
Parasites & Vectors volume 11, Article number: 621 (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.
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,10,11,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,19,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.
To investigate I. ricinus abundance, B. burgdorferi (s.l.) prevalence, and environmental and host associations in relation to landscape fragmentation, we sampled from among six mainland sites and twelve islands sites in 2011, 2013 and 2015 within Loch Lomond and Trossachs National Park (Fig. 1, Table 1). Sampling was conducted every other year to enable patterns of pathogen prevalence and persistence to be studied over an extended temporal time frame. The park is in the west of Scotland and contains the largest natural freshwater area in the United Kingdom. The islands are situated within the southern part of this freshwater lake. The Loch Lomond shoreline is surrounded by mature oak (Quercus spp.) and birch (Betula spp.) woodland with areas of managed coniferous forest. The ground vegetation is a mixture of grasses, mosses, bracken (Pteridium aquilinum), bilberry (Vaccinium myrtillus), bramble (Rubus fruticosus), and ericaceous species (Calluna vulgaris and Rhododendron spp.). The islands sampled within Loch Lomond varied in size from 0.03 km2 to 1.2 km2 and are situated between 0.3 km and 2.4 km from the mainland (Table 1). The predominant woodland cover on the islands, like the mainland sites, is also oak and birch woodland. Mainland sites were chosen to have similar habitat characteristics to the islands, and were located within 0.5 km of the Loch shoreline. Both island and mainland sites were known prior to this study to have populations of I. ricinus, fallow deer (Dama dama), and other woodland mammal and bird species typical of the area (Additional file 1: Table S1) .
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 statistical analyses were carried out in R version 3.4.1 . We calculated the probability p of not detecting B. burgdorferi (s.l.) given a sample of n tested nymphs and an expected prevalence P as
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.
A total of 6567 tick nymphs were collected across all sites. The overall mean prevalence of B. burgdorferi (s.l.) in questing nymphs was 1.4% (91/6567; 95% CI: 1.1–1.7%; range 0–21.2%) (Additional file 1: Table S3). The mean prevalence in questing nymphs on the six mainland sites was 2.5% (47/1891; 95% CI: 1.9–3.3%; range 0–4.6%). The prevalence in questing nymphs on the twelve islands was 0.9% (44/4673; 95% CI: 0.7–1.2%; range 0–21.2%). All four genospecies of B. burgdorferi (s.l.) which are known to occur in the UK [24, 47] were detected at mainland and island sites. Of the 47 mainland nymphs that tested positive for B. burgdorferi (s.l.) 40.4% (19/47) were infected with B. garinii, 31.9% (15/47) with B. burgdorferi (s.s.), 19.1% (9/47) with B. afzelii, and 8.5% (4/47) with B. valaisiana. On island sites, of the 44 nymphs that tested positive for B. burgdorferi (s.l.), 63.6% (28/44) were infected with B. garinii, 22.7% (10/44) with B. burgdorferi (s.s.), 11.4% with B. valaisiana (5/44) and 2.3% with B. afzelii (1/44) (Fig. 2, Additional file 1: Table S3).
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
Data from all sites and years were used to model B. burgdorferi (s.l.) prevalence in questing nymphs (n = 42 site/year sampling events). The best-fit model for B. burgdorferi (s.l.) prevalence included year and fragmentation class (island/mainland) (Table 2), an interaction between these variables was not supported. The proportion of variance explained by the fixed effects R2GLMM(m)* was 34% while the site level random effect accounted for only 3.4% (unexplained variation between sites), and the remaining 62.6% was attributed to the observation level random effect (unexplained variation within sites). Mainland sites were significantly more likely to harbour infected nymphs compared to islands, odds ratio (OR) = 4.8 (95% CI: 1.5–15.5) (delta AICc = 3.4). Year of sampling did not contribute significantly to explaining variation in prevalence (delta AICc for year = 0.75).
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,53,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.
Grenfell BT, Harwood J. (Meta)population dynamics of infectious diseases. Trends Ecol Evol. 1997;12:395–9.
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.
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.
Kilpatrick AM, Randolph SE. Drivers, dynamics, and control of emerging vector-borne zoonotic diseases. Lancet. 2012;380:1946–55.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Stanek G, Reiter M. The expanding Lyme Borrelia complex-clinical significance of genomic species? Clin Microbiol Infect. 2011;17:487–93.
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.
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.
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.
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.
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.
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.
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.
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.
Rajska-Jurgiel E. Demography of woodland rodents in fragmented habitat. Environ Microbiol. 1992;37:73–90.
Paterson VL. Population dynamics of rodents and their parasite communities in a naturally fragmented landscape. PhD Thesis, University of Glasgow; 2012.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Allan B, Keesing F, Ostfeld R. Effect of forest fragmentation on Lyme disease risk. Conserv Biol. 2003;17:267.
Brownstein JS, Skelly DK, Holford TR, Fish D. Forest fragmentation predicts local scale heterogeneity of Lyme disease risk. Oecologia. 2005;146:469–75.
Ostfeld R, Keesing F. Biodiversity and disease risk: the case of Lyme disease. Conserv Biol. 2000;14:722–8.
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.
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.
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.
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.
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.
Text S1. Description of deer dung transect and deer survey methodology used to estimate deer density at Loch Lomond and the Trossachs National Park. Text S2. Methodology to collect nymph density data and to carry out analysis of data collected on islands and mainland sites at Loch Lomond and the Trossachs National Park. Table S1. Counts of selected bird species at study sites from a point transect study carried out in summer 2015. Table S2. Estimates of deer density at mainland and island sites (Locations of sites shown in Fig. 1) at Loch Lomond and the Trossachs National Park. Deer density was estimated using two methods described in Text S1. These methods were: counts of deer carried out on island sites only (Deer survey 2008 & 2012), an estimate of density calculated from deer dung counts along transects spaced at 200 m through each of the study sites carried out in 2015 (Dung Transects). Table S3. Numbers of nymphs tested for Borrelia burgdorferi (sensu lato) at each site, density of nymphs (Nymphs/10 m2), overall prevalence (Prev %) of B. burgdorferi (s.l.) and 95% CI, and the prevalence of each genospecies: B. garinii (B.g); B. afzelii (B.a), Borrelia valaisiana (B.v); Borrelia burgdorferi (sensu stricto) (B.ss) and the number of infected nymphs (n). p(non-detect) represents the probability of failing to detect infected ticks in a given island sample with an estimated prevalence of 0%, assuming an expected B. burgdorferi (s.l.) prevalence of 2.5% (as estimated for the mainland). Asterisks indicate cases for which the calculated probability was lower than 0.05. Table S4. Best model explaining questing nymphal tick variation among eleven island and 5 mainland sites in 2013 using a Poisson Generalised Linear Mixed Model. The best-fit model included vegetation type at the site of the blanket drag. Delta AICc indicates the change in AICc after removing each variable from the best-fit model. (DOCX 53 kb)
About this article
Cite this article
Millins, C., Dickinson, E.R., Isakovic, P. et al. Landscape structure affects the prevalence and distribution of a tick-borne zoonotic pathogen. Parasites Vectors 11, 621 (2018). https://doi.org/10.1186/s13071-018-3200-2