Tick range expansion to higher elevations: does Borrelia burgdorferi sensu lato facilitate the colonisation of marginal habitats?
BMC Ecology and Evolution volume 22, Article number: 104 (2022)
Parasites can alter host and vector phenotype and thereby affect ecological processes in natural populations. Laboratory studies have suggested that Borrelia burgdorferi sensu lato, the causative agent of human Lyme borreliosis, may induce physiological and behavioural alterations in its main tick vector in Europe, Ixodes ricinus, which increase the tick’s mobility and survival under challenging conditions. These phenotypic alterations may allow I. ricinus to colonise marginal habitats (‘facilitation hypothesis’), thereby fuelling the ongoing range expansion of I. ricinus towards higher elevations and latitudes induced by climate change. To explore the potential for such an effect under natural conditions, we studied the prevalence of B. burgdorferi s.l. in questing I. ricinus and its variation with elevation in the Swiss Alps.
We screened for B. burgdorferi s.l. infection in questing nymphs of I. ricinus (N = 411) from 15 sites between 528 and 1774 m.a.s.l to test if B. burgdorferi s.l. prevalence is higher at high elevations (i.e. in marginal habitats). Opposite of what is predicted under the facilitation hypothesis, we found that B. burgdorferi s.l. prevalence in I. ricinus nymphs decreased with increasing elevation and that Borrelia prevalence was 12.6% lower in I. ricinus nymphs collected at the range margin compared to nymphs in the core range. But there was no association between Borrelia prevalence and elevation within the core range of I. ricinus. Therefore the observed pattern was more consistent with a sudden decrease in Borrelia prevalence above a certain elevation, rather than a gradual decline with increasing elevation across the entire tick range.
In conclusion, we found no evidence that B. burgdorferi s.l.-induced alterations of I. ricinus phenotype observed in laboratory studies facilitate the colonisation of marginal habitats in the wild. Rather, ticks in marginal habitats are substantially less likely to harbour the pathogen. These findings have implications for a better understanding of eco-evolutionary processes in natural host-parasite systems, as well as the assessment of Lyme borreliosis risk in regions where I. ricinus is newly emerging.
Parasites can alter the physiology, morphology and behaviour of their hosts [1,2,3]. Vectors can similarly be the target of parasite manipulation [4, 5]. Such manipulations can facilitate parasite transmission, and are thus adaptive from the parasite’s perspective . However, not all phenotypic alterations are the result of such adaptive manipulations [7, 8]. Phenotypic changes in infected hosts and vectors can also be a side effect of infection that is of no adaptive value for the parasite.
Surprisingly little is known about the ecological consequences of parasite-induced phenotypic alterations in hosts or vectors. Recent studies have suggested that parasite-induced changes may have profound effects on host population dynamics, trophic niche specialisation as well as interactions between competitors and predators at all trophic levels [9,10,11]. Consequently, parasites can have a strong ecological impact that becomes especially apparent when they are missing or introduced, e.g. during biological invasions [11, 12]. Invasive host populations often harbour fewer parasites than populations in their native range, which may enhance their population growth and competitive ability (i.e. ‘enemy release hypothesis’ ). Parasites may be ‘lost’ because of low host densities and founder effects during the three phases of an invasion: introduction, establishment and range expansion . However, empirical studies investigating the role of parasites during range expansion from core populations, i.e. without bottlenecks occurring during the introduction and establishment phases, are still rare , and both increasing and decreasing parasite loads at the host range margins have been reported (e.g. [16,17,18]).
The distribution of Ixodes spp. ticks is strongly influenced by abiotic factors, such as temperature and humidity, and there is accumulating evidence that in Europe and North America Ixodes spp. ticks have been expanding northwards as well as towards higher elevations during the last decades due to climate warming (reviewed in [19, 20]). In Europe, I. ricinus is the main vector of the bacterium Borrelia burgdorferi sensu lato, the causative agent of human Lyme borreliosis. Borrelia burgdorferi s.l. forms a complex comprising of at least 18 genospecies . A diverse host community, which includes rodents, insectivores, birds and reptiles, acts as reservoir hosts for B. burgdorferi s.l.. In Europe, B. afzelii and B. garinii are the most common genospecies . Borrelia afzelii is a rodent specialist, infecting mice (Apodemus sylvaticus and A. flavicollis) and voles (Myodes glareolus and Microtus agrestis), but also shrews (Sorex araneus and S. minutus; ), whereas B. garinii is a bird specialist, infecting a range of bird species such as Turdus spp., Sturnus vulgaris, Sylvia spp. and Parus major [21, 24].
Ixodes ricinus is a three host tick. Each life stage (i.e. larva, nymph and adult) seeks a host, feeds, and drops off to develop in the natural environment. Although I. ricinus is a host generalist, the host distribution differs across tick stages [25, 26]. The host selection is, among other things, driven by habitat type, host behaviour, microclimate conditions, and questing height of the different stages . Adults quest higher in the vegetation than nymphs, and nymphs higher than larvae. Larvae and nymphs feed on vertebrate hosts of all size classes but particularly on rodents and birds, while adult females usually feed on larger mammals . Unlike hematophagous insects, ticks typically feed for several days on their host, and the transmission of tick-borne pathogens such as B. burgdorferi s.l. to the host may start several hours after tick attachment .
Various studies have reported that bacteria and viruses may induce phenotypic alterations in Ixodes spp. [26, 27]. Specifically, laboratory studies have suggested that B. burgdorferi s.l.-infected I. ricinus are bigger, walk faster and longer, quest longer for a host, have higher energy reserves (for the same body mass) and survive better under challenging thermohygrometric conditions than uninfected I. ricinus. Yet to date, studies investigating such vector manipulation by pathogens remain scarce and results are heterogeneous ([26, 27], and references therein). These alterations of tick behaviour and physiology may benefit pathogen transmission by extending vector lifespan, or by increasing host-finding efficacy . Importantly, both may facilitate the colonisation of marginal habitats by ticks, and lead to particularly high prevalence of B. burgdorferi s.l. in I. ricinus populations at range margins . Particularly, higher fat contents in Borrelia-infected nymphs may provide an advantage under colder or more variable temperatures, such as at higher elevations, although laboratory data were not conclusive . Furthermore, increased host-finding efficacy may facilitate the establishment of tick populations in less diverse or more fragmented host communities. Finally by modulating questing behaviour, such as adjusting questing height or moving to shade, Borrelia infection may shape the interaction between ticks and host species [26, 27]. Specifically, questing higher up in the vegetation may allow ticks to attach to larger animals. Since larger animals have larger home ranges , questing at higher heights may allow ticks to disperse further. All these possible mechanisms may result in higher B. burgdorferi s.l. prevalences in I. ricinus populations at range margins. Assessing the role of B. burgdorferi s.l. in influencing vector range expansion to marginal habitats is crucial for our understanding of eco-evolutionary processes in natural host-parasite systems, but also for quantifying public health threats in regions where ticks are newly emerging. Yet, to our knowledge no study has tested this hypothesis in I. ricinus populations under natural conditions to date.
Although exclusion experiments (i.e. where a taxon is excluded) are the gold standard to assess the ecological role of parasites, such manipulations remain an enormous practical and logistic challenge in I. ricinus and its pathogens. Moreover, experiments in the wild testing whether B. burgdorferi s.l. enhances survival and reproduction of I. ricinus under harsh conditions (i.e. physiological alterations) would only provide a partial picture because B. burgdorferi s.l. might affect dispersal abilities of I. ricinus through its ability to find a host (i.e. behavioural alterations). Because of these practical challenges, we used a correlational approach in naturally infected tick populations in the Swiss Alps to evaluate the potential of B. burgdorferi s.l. to facilitate the colonisation of marginal habitats by I. ricinus. Previous studies monitored B. burgdorferi s.l. and I. ricinus within their core range [31,32,33,34,35,36], at the range margins without a direct comparison with core populations  or when B. burgdorferi s.l. prevalence was low [38,39,40] (but see ). In our study we monitored tick populations at different elevations, including core range and range margin populations, allowing for a direct comparison of Borrelia prevalence across the tick distribution range. We predict that if B. burgdorferi s.l. alters the phenotype of I. ricinus, thereby making them better colonizers of unfavourable habitats, B. burgdorferi s.l. prevalence will be disproportionally high at the tick range margin (i.e. at high elevations).
Questing I. ricinus abundance along elevation
We collected a total of 1138 questing nymphs and 270 questing adults across three 1 h-dragging sessions at 15 sites each (Table 1; Fig. 1; 0–109 nymphs and 0–44 adults per dragging event, Table 2). Overall, we collected at least four times more nymphs than adults. Patterns of tick abundance along elevation did not vary between tick stages (Interaction stage × elevation2: χ21 = 0.88, P = 0.349; Interaction stage × elevation: χ21 = 3.17, P = 0.075). Therefore, nymphs and adults were pooled for subsequent analyses.
Questing I. ricinus abundance decreased non-linearly with increasing elevation with an abundance peak at intermediate elevations (Elevation2: χ21 = 17.62, P < 0.001, β2 = − 1.057 ± 0.206; Elevation: χ21 = 17.06, P < 0.001, β1 = − 0.903 ± 0.170, Fig. 2a). No evidence for seasonal differences in tick abundance across elevations was found (Nymphs: Interaction session × elevation2: χ21 = 0.88, P = 0.645; Interaction session × elevation: χ21 = 2.95, P = 0.229; nymphs and adults pooled: Interaction session × elevation2: χ21 = 1.55, P = 0.461; Interaction session × elevation: χ21 = 3.28, P = 0.194).
Questing I. ricinus abundance was high up to an elevation of 1100 m.a.s.l., and then rapidly decreased (Table 1, Fig. 2a). No questing I. ricinus were found at the highest elevation site (VIL) at 1680 m.a.s.l.. We therefore defined sites below 1150 m.a.s.l. as the core range and sites between 1450 and 1680 m.a.s.l. to represent the elevational range margin of I. ricinus ticks in our study area (note that no sites between 1150–1450 m.a.s.l. were included in our study). The highest elevation site (VIL) was excluded from subsequent analyses because it is not yet colonised by ticks (i.e. outside the tick range). Questing tick abundance was higher in the core range than in the elevational range margin (Range: χ21 = 13.44, P < 0.001, excluding VIL, Fig. 2a).
Borrelia infection in I. ricinus nymphs across the tick range
Overall, 24.6% of the analysed I. ricinus nymphs (101 out of 411) were infected with B. burgdorferi s.l.. 83.2% of these infections were caused by B. afzelii (Table 1). Prevalence among sites varied between 0 and 47.1% for B. burgdorferi s.l., and 0–44.1% for B. afzelii (Table 1).
Borrelia burgdorferi s.l. prevalence decreased with increasing elevation (Elevation2: χ21 = 3.10, P = 0.078; Elevation: χ21 = 6.38, P = 0.011, β = − 0.664, Fig. 2b), whereas for B. afzelii the decrease with increasing elevation was marginally non-significant (Elevation2: χ21 = 3.19, P = 0.074; Elevation: χ21 = 3.39, P = 0.066; Fig. 2c). Borrelia prevalence was 12.6 and 9.8 times higher in I. ricinus nymphs collected within the core range than at the range margin for B. burgdorferi s.l. (χ21 = 9.55, P = 0.002, Fig. 2b) and B. afzelii (χ21 = 7.56, P = 0.006, Fig. 2c), respectively.
If abiotic factors underlie the decrease of Borrelia prevalence with increasing elevation across the entire tick range, we would predict to see a gradual decline of B. burgdorferi s.l. prevalence also within the core range of ticks. However, we found no evidence that Borrelia prevalence declined with elevation within the core range of ticks (B. burgdorferi s.l.: Elevation2: χ21 = 0.07, P = 0.786; Elevation: χ21 = 0.81, P = 0.369; B. afzelii: Elevation2: χ21 = 0.13, P = 0.718; Elevation: χ21 = 0.05, P = 0.823). Moreover, within the core range of ticks, no association between Borrelia prevalence and questing I. ricinus abundance was found for B. burgdorferi s.l. (Questing tick abundance2: χ21 = 2.19, P = 0.1439; Questing tick abundance: χ21 = 1.85, P = 0.174, Fig. 2b) or B. afzelii (Questing tick abundance2: χ21 = 0.57, P = 0.449; Questing tick abundance: χ21 = 0.47, P = 0.493, Fig. 2c).
Laboratory experiments have suggested that B. burgdorferi s.l. may increase tick tolerance to desiccation and boost host-finding efficacy ([26, 27], and references therein). Such vector manipulations by the pathogen may facilitate the range expansion of ticks to marginal habitats, such as higher elevations or higher latitudes . Here we used a correlational approach to investigate whether patterns of B. burgdorferi s.l. prevalence in questing I. ricinus nymphs, and its variation with elevation, provide support for this hypothesis.
Overall, questing I. ricinus abundance was high and B. burgdorferi s.l. infection common in the study area: I. ricinus were found at 14 out of 15 sites, and B. burgdorferi s.l. infection in questing I. ricinus nymphs was detected at 12 sites. Above 1400 m.a.s.l., however, I. ricinus abundance strongly decreased and no I. ricinus were found above 1700 m.a.s.l.. B. burgdorferi s.l. prevalence in questing I. ricinus nymphs across sites was similar to previous reports from central Europe (e.g. , and references therein, ), but variation across sites was high (0–47% compared to 1–20% across Europe according to ), and it exceeded 40% at two sites of contrasting elevations (630 and 1123 m.a.s.l.).
As expected, the most common B. burgdorferi s.l. genospecies was B. afzelii, accounting for 83.2% of I. ricinus infections in our study. This is higher than what has been reported previously (7–68% of I. ricinus infections caused by B. afzelii; ) across Europe but similar to what has recently been reported from the Italian Alps ( where 73.1% of the nymphs were infected by B. afzelii). Differences across studies and locations are likely due to differences in the composition of the local host community, i.e. the relative abundance of rodent, bird and reptile hosts. Indeed the different B. burgdorferi s.l. genospecies are associated with different reservoir hosts. For example, B. afzelii and B. bavariensis are rodent specialists, B. garinii and B. valaisiana are bird specialists, and B. burgdorferi s.s. is a generalist that can infect birds and rodents (reviewed in ). Alternatively, strong bottlenecks within sites could lead to the local dominance of specific genospecies that may vary over time (e.g. [43, 44]).
We found no evidence that B. burgdorferi s.l. prevalence is higher at the range margins of I. ricinus (i.e. at higher elevations), indicating that B. burgdorferi s.l.-mediated modifications of I. ricinus physiology and behaviour observed in the laboratory ([26, 27], and references therein) plays a minor role in the ongoing colonisation process of marginal habitats in the wild. Rather, I. ricinus nymphs at the range margin had substantially lower B. burgdorferi s.l. prevalence. This finding is in line with what was found in the Italian Alps, where higher Borrelia prevalence was observed at locations < 1400 m.a.s.l compared to > 1400 m.a.s.l (note that ticks were found up to > 1800 m.a.s.l at these study sites) . Our analytical approach allows us to describe changes in Borrelia prevalence with elevation in more detail and highlights that the patterns we observe are more consistent with a sudden decrease above a certain elevation, with reduced Borrelia prevalence specifically at the range margin, rather than a gradual decline with increasing elevation across the entire tick range. Indeed, analyses within the core range of ticks revealed no association between Borrelia prevalence and elevation.
The observed patterns may be due to different non-exclusive processes. The colder microclimate at high elevations may decrease the survival of B. burgdorferi s.l. in I. ricinus and slow down its multiplication and transmission. Effects of meteorological temperature on pathogen development (i.e. the development from the ingested infectious stage to the stage in the salivary glands) and transmission efficiency have been observed in malaria-mosquito [45, 46] and trypanosome-tsetse fly systems (references in ). Although protozoan development is not directly comparable to bacterial multiplication, many bacteria similarly utilize the environmental temperature as a signal to determine their location and to regulate expression of a large set of proteins necessary for survival, multiplication, migration and transmission [48, 49]. For example, the expressions of Outer Surface Protein A and C of B. burgdorferi s.l., which are involved in its survival in the tick midgut and its dissemination into the tick salivary glands or into its vertebrate host, appear to be mediated by differences in ambient temperature . To our knowledge, no study evaluated the effect of meteorological temperature on the expression of vector-specific or host-specific proteins, or the expression of virulence genes in B. burgdorferi s.l.. Although temperature and other abiotic factors are often involved to explain the lower pathogen prevalence at higher elevations , the lack of a gradual decline in Borrelia prevalence with increasing elevation within the core range of ticks makes it unlikely that such factors are the sole explanation for the low Borrelia prevalence at range margins.
Alternatively, the pattern may also be affected by other processes such as changes in host community, host species abundance, habitat quality or population structure of hosts/vectors in marginal habitats. Indeed, Patot et al. (2010) found no effect of temperature on filamentous virus transmission efficiency, but observed a clear relationship between virus prevalence and the density of its host, a parasitoid wasp. Similarly, B. burgdorferi s.l. transmission may be lower in I. ricinus populations at the range margin that are more fragmented , less dense  or/and less genetically diverse  than populations at the core range. Similar processes may act on the reservoir host populations of B. burgdorferi s.l. (e.g. rodents and birds). However detailed knowledge about the composition of the host community, host abundance, I. ricinus infestation of these host populations as well as B. burgdorferi s.l. prevalence would be necessary to link host abundance and B. burgdorferi s.l. prevalence in questing I. ricinus.
Finally, because edge populations are usually founded by only few individuals, stochastic processes may lead to the loss of pathogens in ticks at the range margins. Thus, the combination of stochastic events and host and vector population structure that hinder B. burgdorferi s.l. transmission may lead to local extinctions of B. burgdorferi s.l. at the range edges (e.g. ). It is often assume that vertebrate hosts dispersing ticks on wide ranges are mainly birds and wild cervids (e.g. red deers and chamois). One may argue that because the latter are incompetent hosts for Borrelia spp., it may dilute infection prevalence in ticks at the range margins. The lower Borrelia infection prevalence at the range margins supports this hypothesis. However a sudden drop in Borrelia prevalence above a certain elevation would also be consistent with a higher mortality of ticks driven by abiotic factors. Moreover the dominance of B. afzelii in this study, and in the Italian Alps , suggests that small rodents may play a particularly important role in the maintenance of Borrelia infection along elevational clines in the Alps. Though, further investigations are needed to evaluate these hypotheses in detail.
Our study is correlational and has limitations. First, fewer ticks were collected at the range margin, which will decrease the precision of the prevalence estimate at these sites. Second, a time-based approach (which is less reliable than an area-based approach) was used to estimate questing tick abundance because of the rugged mountain terrain of the study area. Third, vegetation structure varies substantially across the elevation range and could affect tick detection despite a consistent sampling effort. However, the decrease of I. ricinus abundance above 1150 m.a.s.l. was observed in different life-stages (Lemoine pers. obs.), which have different questing behaviours and biological needs. Therefore, we are confident that the decrease of I. ricinus abundance above 1150 m.a.s.l. is not a methodological artefact but describes an ecological pattern. Fourth, our approach does not allow to disentangle between range-expanding ticks bringing the pathogen with them (i.e., dispersal to a novel habitat) and the enzootic cycle being less likely to become established (i.e., establishment in a novel habitat, for example infecting less dense host populations). Fifth, whereas numerous studies in Europe and North America have documented a range expansion of ticks to higher elevations and latitudes due to climate change (reviewed in ), we do not have the long-term monitoring data to directly demonstrate this range expansion at our study sites. Finally, the presence of questing I. ricinus ticks in the vegetation does not prove that I. ricinus populations are established locally, or will survive through the winter. High elevation habitats might thus represent a ‘ragged edge’ rather than an expansion front. Ultimately, an experimental approach, in which the B. burgdorferi s.l. infection status of I. ricinus is manipulated and the colonisation of marginal habitats is monitored would be required to test conclusively whether B. burgdorferi s.l. facilitates tick range expansion to marginal habitats in the wild. But as stated previously, this would be an enormous practical challenge.
Overall, we found no evidence that B. burgdorferi s.l.-induced changes in I. ricinus behaviour or physiology that facilitate I. ricinus range expansion to higher elevations in the Swiss Alps. Rather, questing I. ricinus in marginal habitats are less likely to carry B. burgdorferi s.l.. These findings show that when bitten by a tick, the risk of human B. burgdorferi s.l. infection is lower, rather than higher, in regions where I. ricinus is newly emerging. Low B. burgdorferi s.l. prevalence at I. ricinus range margins may enhance population growth and competitive ability of hosts and vectors. Less infected, hosts may, for example, invest differently in immunity and reproduction than hosts in core populations [15, 54], which can affect host-parasite interactions when the parasite finally invades host populations at the range margins . A better understanding of eco-evolutionary processes between pathogens, vectors and hosts at range margins, and their effect on pathogen life-history and virulence evolution, will therefore be a fruitful next step (e.g. ), and will contribute to a better prediction of zoonotic disease risks in regions where vectors and pathogens are newly emerging due to climate change.
We sampled questing I. ricinus at fifteen sites across an elevational range from 528 to 1774 m.a.s.l. in the Swiss Alps (Kanton Graubünden, Fig. 1, Table 1). Each site was visited three times, once in June, July and August 2014 to obtain a measure of questing tick abundance per site. During each dragging event, four people dragged a white blanket (1 m × 1 m each) slowly over the ground vegetation during 15 min each. During these 15 min, each collector regularly checked for ticks on their blanket and transferred them to Eppendorf tubes containing 95% ethanol. The count of each collector was pooled to provide an estimate of questing tick abundance per dragging event per hour (Table 2). Although an area-based estimation is often more accurate than a time-based approach, it was not possible to apply such a method in our study because of the rugged mountain terrain of the Swiss Alps. Tick life stage and species were verified with a dissection microscope in the laboratory on the basis of morphologic features following .
Borrelia burgdorferi s.l. infection in questing Ixodes ricinus nymphs.
We randomly selected 34–35 questing I. ricinus nymphs per site for B. burgdorferi s.l. detection. When fewer than 35 nymphs were collected at a site (Table 1), additional questing nymphs were collected for Borrelia detection one day after estimating questing tick abundance. All collected nymphs were analysed when fewer than 35 nymphs were available (Table 1). DNA was extracted using the ‘HotShot’ method [57, 58] with slight modifications. Each nymph was incubated with 120 μL of alkaline lysis reagent (25 mM NaOH, 0.2 mM EDTA) and a metal bead (2 mm in diameter) at 95 °C for 15 min. After homogenisation with a Retsch TissueLyser (Haan, Germany) for 2 min, the samples were incubated at 95 °C for 15 min. After cooling, 60 μL 40 mM Tris–HCl was added to each tube.
To assess the B. burgdorferi s.l. infection status of I. ricinus nymphs, we used a combination of two complimentary approaches. First, we specifically targeted B. afzelii, the most common Borrelia genospecies at the study sites , using a highly sensitive quantitative real-time PCR (qPCR) assay using the B. afzelii-specific primers Fla5F 5′-CACCAGCATCACTTTCAGGA-3′ and Fla6R 5′-CTCCCTCACCAGCAAAAAGA-3′  on a StepOnePlus Real-Time PCR machine (Applied Biosystems). Each reaction contained 10 μL of SYBR Select Master Mix (Applied Biosystems), 0.8 μL of each primer (10 μM), 4.4 μL of water and 4 μL of template DNA in a final volume of 20 μL. Two series of four standards, two negative and two positive tick samples, as well as two no-template controls were included in each run. The PCR amplification protocol consisted of an initial denaturation step at 50 °C and 95 °C for 2 min each, followed by 42 cycles of 95 °C for 15 s, 59 °C for 30 s, and 72 °C for 30 s. The length of the amplicon was 129 bp and the melting temperature was between 77.3 °C and 77.9 °C. A 100% repeatability was obtained based on 40 samples (10 positives and 30 negatives) when amplification occurred before 34 cycles. Out of 85 B. afzelii-positive samples, only two samples amplified between 34 and 37 cycles. These samples were repeated and considered to be B. afzelii-infected only if they were found to be positive twice.
In a second step, all samples that were found to be B. afzelii-negative using the qPCR approach described above were screened for the presence of other B. burgdorferi s.l. genospecies using 16S-LD primers designed to target all Borrelia genospecies associated with Lyme borreliosis . Amplifications were performed in a total volume of 10 μL containing 0.2 μL JumpStart Taq DNA Polymerase (Sigma-Aldrich), 0.5 μL of each primer (300 nM; LD-F: 5′-ATGCACACTTGGTGTTAACTA-3′ and LD-R: 5′-GACTTATCACCGGCAGTCTTA-3′) and 2 μL of DNA template. The PCR amplification protocol consisted of an initial denaturation step at 94 °C for 1 min, followed by 42 cycles of denaturation at 94 °C for 30 s, annealing at 50 °C for 30 s, and extension at 72 °C for 90 s, with a final elongation step at 72 °C for 10 min. PCR products were visualized under UV light on 1% agarose gels that were stained with SYBR® Safe DNA gel stain (Thermo Scientific). The length of the amplicon was 351 bp. We used DNA from six reference genospecies B31 and Pbre (B. burgdorferi s.s), Phei (B. garinii), Pko and PVPM (B. afzelii) and Pbi (B. bavariensis) to verify that they are detected with the 16S approach. All references were successfully amplified.
Each run included B. afzelii-positive samples identified with the qPCR (see above) as controls. All of these samples were found to be B. burgdorferi s.l.-positive using the 16S method. The repeatability based on 39 samples (12 negative and 27 positive samples) was 85%. Out of the six samples that were not repeatable with the 16S method, three were previously quantified with the qPCR method and were found to have very low B. afzelii DNA concentrations (i.e. more than 30 cycles on the qPCR) suggesting that the sensitivity of the 16S approach decreases at low infection intensities. During the development of the methods, some fragments were sequenced and all were found to be B. burgdorferi s.l. suggesting that both methods were specific to the Lyme borreliosis causing genospecies.
Questing I. ricinus abundance (Table 2) was analysed using a generalized mixed effect model with the negative binomial distribution (with the parameterization σ = μ(1 + μ/k)) and site ID was included as a random effect. Preliminary analyses showed that the negative binomial distribution (with the parameterization σ = μ(1 + μ/k)) described the abundance of nymphs (or nymphs and adults combined) best among the following distributions: Poisson, negative binomial (with the parameterization σ = μ(1 + μ/k) or σ = φμ) with or without zero-inflation (results not shown). We first tested for an association between tick abundance across life stages (nymphs and adults) and elevation or season. To this end, we included elevation and elevation2 in interaction with tick stage, or in interaction with sampling session as fixed effects (i.e. as explanatory variables to quantify an average effect, β) in the mixed model. Finally, we defined the range margin based on questing tick abundance (see Results) and included range type (range core (C) vs range margin (M)) of each site as a fixed effect.
Borrelia infection (i.e. B. burgdorferi s.l. or B. afzelii infection) in I. ricinus nymphs (Table 2) was analysed using a generalized mixed effect model with a binomial error structure. The number of infected and non-infected nymphs per dragging event was used as the dependant variable and site ID was included as a random effect. In a first step, we tested for an association between B. burgdorferi s.l. or B. afzelii prevalence and elevation by including elevation and elevation2 as fixed effects (i.e. as explanatory variables to quantify an average effect, β) in the mixed model. In a second step, we tested if B. burgdorferi s.l. or B. afzelii prevalence was higher at the range margin of I. ricinus (i.e. at high elevations). To this end, we compared the B. burgdorferi s.l. or B. afzelii prevalence in I. ricinus at the range margin (sites > 1400 m.a.s.l.) to B. burgdorferi s.l. or B. afzelii prevalence in the core range of I. ricinus (sites < 1150 m.a.s.l.) by including range type (range core (C) vs range margin (M)) of each site as a fixed effect.
Finally, we focused on the patterns within the core range of I. ricinus by investigating the association between B. burgdorferi s.l. or B. afzelii prevalence and elevation, and the relationship between B. burgdorferi s.l. or B. afzelii prevalence and questing I. ricinus abundance within the core range of I. ricinus. The effect of questing I. ricinus abundance was investigated by including questing I. ricinus abundance and questing I. ricinus abundance2 as fixed effects. Site ID was again included as a random effect.
To ensure independence between main and quadratic terms, all continuous variables were standardized (by subtracting the population mean from each sample and dividing it by the standard deviation). The significance of explanatory variables was assessed using likelihood ratio tests. The significance of the linear terms was assessed after removing the non-significant quadratic terms. Statistical analyses were performed using the package glmmADMB [60, 61] in R 3.2.3 .
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article and its additional files.
Metres above mean sea level
- B. :
- I. :
Polymerase chain reaction
Outside of the range
Price PW, editor. Evolutionary biology of parasites. Princeton: Princeton Univ. Press; 1980.
Poulin R. The evolution of parasite manipulation of host behaviour: a theoretical analysis. Parasitology. 1994;109(S1):S109–18.
Thomas F, Schmidt-Rhaesa A, Martin G, Manu C, Durand P, Renaud F. Do hairworms (Nematomorpha) manipulate the water seeking behaviour of their terrestrial hosts? J Evol Biol. 2002;15(3):356–61.
Dobson AP. The population biology of parasite-induced changes in host behavior. Q Rev Biol. 1988;63(2):139–65.
Hurd H. Manipulation of medically important insect vectors by their parasites. Annu Rev Entomol. 2003;48:141–61.
Lefevre T, Thomas F. Behind the scene, something else is pulling the strings: emphasizing parasitic manipulation in vector-borne diseases. Infect, Genet Evol. 2008;8(4):504–19.
Cezilly F, Favrat A, Perrot-Minnot MJ. Multidimensionality in parasite-induced phenotypic alterations: ultimate versus proximate aspects. J Exp Biol. 2013;216(Pt 1):27–35.
Thomas F, Adamo S, Moore J. Parasitic manipulation: where are we and where should we go? Behav Processes. 2005;68(3):185–99.
Britton JR, Andreou D. Parasitism as a driver of trophic niche specialisation. Trends Parasitol. 2016;32(6):437–45.
Dunn AM, Torchin ME, Hatcher MJ, Kotanen PM, Blumenthal DM, Byers JE, et al. Indirect effects of parasites in invasions. Funct Ecol. 2012;26(6):1262–74.
Wood CL, Johnson PTJ. A world without parasites: exploring the hidden ecology of infection. Front Ecol Environ. 2015;13(8):425–34.
Hatcher MJ, Dick JTA, Dunn AM. Parasites that change predator or prey behaviour can have keystone effects on community composition. Biol Lett. 2014;10(1):20130879.
Heger T, Jeschke JM. The enemy release hypothesis as a hierarchy of hypotheses. Oikos. 2014;123(6):741–50.
Kolar CS, Lodge DM. Progress in invasion biology: predicting invaders. Trends Ecol Evol. 2001;16(4):199–204.
White TA, Perkins SE. The ecoimmunology of invasive species. Funct Ecol. 2012;26(6):1313–23.
Briers RA. Range limits and parasite prevalence in a freshwater snail. Proc R Soc Lond Ser B Biol Sci. 2003;270(Suppl 2):S178–80.
Patot S, Martinez J, Allemand R, Gandon S, Varaldi J, Fleury F. Prevalence of a virus inducing behavioural manipulation near species range border. Mol Ecol. 2010;19(14):2995–3007.
Phillips BL, Kelehear C, Pizzatto L, Brown GP, Barton D, Shine R. Parasites and pathogens lag behind their host during periods of host range advance. Ecology. 2010;91(3):872–81.
Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Pena A, George JC, et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasites Vectors. 2013;6.
Ostfeld RS, Brunner JL. Climate change and Ixodes tick-borne diseases of humans. Philos Trans R Soc Lond Ser B Biol Sci. 2015;370(1665):20140051.
Margos G, Vollmer SA, Ogden NH, Fish D. Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato. Infect, Genet Evol. 2011;11(7):1545–63.
Rauter C, Hartung T. Prevalence of Borrelia burgdorferi sensu lato genospecies in Ixodes ricinus ticks in Europe: a metaanalysis. Appl Environ Microbiol. 2005;71(11):7203–16.
Hellgren O, Andersson M, Raberg L. The genetic structure of Borrelia afzelii varies with geographic but not ecological sampling scale. J Evol Biol. 2011;24(1):159–67.
Comstedt P, Bergström S, Olsen B, Garpmo U, Marjavaara L, Mejlon H, et al. Migratory passerine birds as reservoirs of Lyme borreliosis in Europe. Emerg Infect Dis. 2006;12(7):1087.
McCoy KD, Léger E, Dietrich M. Host specialization in ticks and transmission of tick-borne diseases: a review. Front Cell Infect Microbiol. 2013;3:57.
Herrmann C, Gern L. Search for blood or water is influenced by Borrelia burgdorferi in Ixodes ricinus. Parasit Vectors. 2015;8:6.
Benelli G. Pathogens manipulating tick behavior-through a glass, darkly. Pathogens. 2020;9(8):664.
Gassner F, Hartemink N. Tick–Borrelia interactions: burden or benefit? Ecology of parasite-vector interactions: Springer; 2013. p. 141-54.
Herrmann C, Gern L. Survival of Ixodes ricinus (Acari: Ixodidae) nymphs under cold conditions is negatively influenced by frequent temperature variations. Ticks Tick Borne Dis. 2013;4(5):445–51.
Madhav NK, Brownstein JS, Tsao JI, Fish D. A dispersal model for the range expansion of blacklegged tick (Acari: Ixodidae). J Med Entomol. 2004;41(5):842–52.
Aeschlimann A, Chamot E, Gigon F, Jeanneret J-P, Kesseler D, Walther C. B. burgdorferi in Switzerland. Zentralblatt für Bakteriologie, Mikrobiologie und Hygiene. 1987;263(3):450–8.
Jouda F, Perret JL, Gern L. Ixodes ricinus density, and distribution and prevalence of Borrelia burgdorferi sensu lato infection along an altitudinal gradient. J Med Entomol. 2004;41(2):162–9.
Moran Cadenas F, Rais O, Jouda F, Douet V, Humair PF, Moret J, et al. Phenology of Ixodes ricinus and infection with Borrelia burgdorferi sensu lato along a north- and south-facing altitudinal gradient on Chaumont Mountain, Switzerland. J Med Entomol. 2007;44(4):683–93.
Gern L, Moran Cadenas F, Burri C. Influence of some climatic factors on Ixodes ricinus ticks studied along altitudinal gradients in two geographic regions in Switzerland. Int J Med Microbiol. 2008;298:55–9.
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(2):237–46.
Stünzner D, Hubálek Z, Halouzka J, Wendelin I, Sixl W, Marth E. Prevalence of Borrelia burgdorferi sensu lato in the tick Ixodes ricinus in the Styrian mountains of Austria. Wien Klin Wochenschr. 2006;118(21):682–5.
Ragagli C, Mannelli A, Ambrogi C, Bisanzio D, Ceballos L, Grego E, et al. Presence of host-seeking Ixodes ricinus and their infection with Borrelia burgdorferi sensu lato in the Northern Apennines. Italy Exp Appl Acarol. 2016;69(2):167–78.
Danielová V, Daniel M, Schwarzová L, Materna J, Rudenko N, Golovchenko M, et al. Integration of a tick-borne encephalitis virus and Borrelia burgdorferi sensu lato into mountain ecosystems, following a shift in the altitudinal limit of distribution of their vector, Ixodes ricinus (Krkonoše mountains, Czech Republic). Vector-Borne Zoonotic Dis. 2010;10(3):223–30.
Daniel M, Materna J, Hönig V, Metelka L, Danielová V, Harcarik J, et al. Vertical distribution of the tick Ixodes ricinus and tick-borne pathogens in the northern Moravian mountains correlated with climate warming (Jeseníky Mts., Czech Republic). Central Eur J Public Health. 2009;17(3):139.
Danielová V, Rudenko N, Daniel M, Holubová J, Materna J, Golovchenko M, et al. Extension of Ixodes ricinus ticks and agents of tick-borne diseases to mountain areas in the Czech Republic. Int J Med Microbiol. 2006;296:48–53.
Garcia-Vozmediano A, Krawczyk AI, Sprong H, Rossi L, Ramassa E, Tomassone L. Ticks climb the mountains: Ixodid tick infestation and infection by tick-borne pathogens in the Western Alps. Ticks Tick-borne Dis. 2020;11(5): 101489.
Herrmann C, Gern L, Voordouw MJ. Species co-occurrence patterns among Lyme Borreliosis pathogens in the tick vector Ixodes ricinus. Appl Environ Microbiol. 2013;79(23):7273–80.
Bruyndonckx N, Henry I, Christe P, Kerth G. Spatio-temporal population genetic structure of the parasitic mite Spinturnix bechsteini is shaped by its own demography and the social system of its bat host. Mol Ecol. 2009;18:3581–92.
Criscione CD, Blouin MS. Minimal selfing, few clones, and no among-host genetic structure in a hermaphroditic parasite with asexual larval propagation. Evolution. 2006;60(3):553–62.
Eling W, Hooghof J, van de Vegte-Bolmer M, Sauerwein R, Van Gemert G, editors. Tropical temperatures can inhibit development of the human malaria parasite Plasmodium falciparum in the mosquito. Proceedings of the Section Experimental and Applied Entomology of the Netherlands Entomological Society (N. E. V.) Amsterdam; 2001.
Noden B, Kent M, Beier JC. The impact of variations in temperature on early Plasmodium falciparum development in Anopheles stephensi. Parasitology. 1995;111(05):539–45.
Moore S, Shrestha S, Tomlinson KW, Vuong H. Predicting the effect of climate change on African trypanosomiasis: integrating epidemiology with parasite and vector biology. J R Soc Interface. 2012;9(70):817–30.
Miller JF, Mekalanos JJ, Falkow S. Coordinate regulation and sensory transduction in the control of bacterial virulence. Science. 1989;17(243):916–22.
Steinmann R, Dersch P. Thermosensing to adjust bacterial virulence in a fluctuating environment. Future Microbiol. 2013;8(1):85–105.
Schwan TG, Piesman J. Vector interactions and molecular adaptations of Lyme disease and relapsing fever spirochetes associated with transmission by ticks. Emerg Infect Dis. 2002;8(2):115.
Lafferty KD. The ecology of climate change and infectious diseases. Ecology. 2009;90(4):888–900.
Hanski I. Metapopulation ecology. New York: Oxford University Press; 1999.
Sexton JP, McIntyre PJ, Angert AL, Rice KJ. Evolution and ecology of species range limits. 2009.
Cornetti L, Hilfiker D, Lemoine M, Tschirren B. Small-scale spatial variation in infection risk shapes the evolution of a Borrelia resistance gene in wild rodents. Mol Ecol. 2018;27(17):3515–24.
Kelehear C, Brown GP, Shine R. Rapid evolution of parasite life history traits on an expanding range-edge. Ecol Lett. 2012;15(4):329–37.
Hillyard PD. Ticks of north-west Europe: Field Studies Council; 1996.
Råberg L. Infection intensity and infectivity of the tick-borne pathogen Borrelia afzelii. J Evol Biol. 2012;25:1448–53.
Truett GE, Heeger P, Mynatt RL, Truett AA, Walker JA, Warman ML. Preparation of PCR-quality mouse genomic DNA with hot sodium hydroxide and tris (HotSHOT). Biotechniques. 2000;29(1):52–4.
Marconi RT, Garon CF. Development of polymerase chain reaction primer sets for diagnosis of Lyme disease and for species-specific identification of Lyme disease isolates by 16S rRNA signature nucleotide analysis. J Clin Microbiol. 1992;30(11):2830–4.
Fournier DA, Skaug HJ, Ancheta J, Ianelli J, Magnusson A, Maunder M, et al. AD Model Builder: using automatic differentiation for statistical inference of highly parameterized complex nonlinear models. Optim Methods Softw. 2012;27:233–49.
Skaug HJ, Fournier DA, Bolker B, Magnusson A, Nielsen A. Generalized linear mixed models using AD model builder. R package version 0.8.0. 2014.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2013.
We thank the local authorities for their permission to collect ticks in their municipal areas, numerous students for their help with data collection, Reto Lienhard, ADMED Microbiologie, for providing the reference strains, and Kay Lucek and two anonymous reviewers for helpful comments on the manuscript.
The study was financially supported by the University of Zurich Research Priority Program ‘Evolution in Action: From Genomes to Ecosystems’ (to ML and BT), the Swiss National Science Foundation (PMPDP3_151361 and PMPDP3_161858 to ML; PP00P3_128386 and PP00P3_157455 to BT), the Baugarten Stiftung & Stiftung für wissenschaftliche Forschung an der Universität Zürich (STWF-17-027), and the Georges und Antoine Claraz-Schenkung.
Ethics approval and consent to participate
Consent for publication
The authors declare there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lemoine, M., Cornetti, L., Reeh, K. et al. Tick range expansion to higher elevations: does Borrelia burgdorferi sensu lato facilitate the colonisation of marginal habitats?. BMC Ecol Evo 22, 104 (2022). https://doi.org/10.1186/s12862-022-02058-x