- Open Access
Diversity of Borrelia burgdorferi sensu lato in ticks and small mammals from different habitats
Parasites & Vectors volume 15, Article number: 195 (2022)
Ixodid ticks are important vectors for zoonotic pathogens, with Ixodes ricinus being the most important in Europe. Rodents are hosts of immature life stages of I. ricinus ticks and are considered main reservoirs for tick-borne pathogens, e.g. Borrelia burgdorferi. The aim of this study was to analyse the prevalence as well as genospecies and sequence type (ST) diversity of Borrelia burgdorferi sensu lato in ticks and small mammals from central Germany and to elaborate on the influence of environmental and/or individual host and vector factors on Borrelia prevalence.
After species identification, 1167 small mammal skin samples and 1094 ticks from vegetation were screened by B. burgdorferi sensu lato real-time polymerase chain reaction, and positive samples were characterized by multilocus sequence typing. Generalized linear (mixed) models were used to estimate how seasonality, small mammal species/tick life stage and habitat affect individual infection status.
In total, 10 small mammal species and three tick species, Ixodes ricinus, Ixodes inopinatus (both considered members of the I. ricinus complex) and Dermacentor reticulatus, were investigated. Borrelia DNA was detected in eight host species, i.e. the striped field mouse (Apodemus agrarius), the yellow-necked field mouse (Apodemus flavicollis), the wood mouse (Apodemus sylvaticus), the water vole (Arvicola amphibius), the bank vole (Clethrionomys glareolus), the field vole (Microtus agrestis), the common vole (Microtus arvalis), and the common shrew (Sorex araneus). Two species were Borrelia negative, the greater white-toothed shrew (Crocidura russula) and the pygmy shrew (Sorex minutus). The average prevalence was 6.2%, with two genospecies detected, Borrelia afzelii and Borrelia garinii, and at least three STs that had not been previously reported in small mammals. Borrelia prevalence in small mammals did not differ between seasons. Six genospecies of Borrelia—Borrelia afzelii, Borrelia valaisiana, Borrelia garinii, Borrelia lusitaniae, Borrelia spielmanii, and Borrelia burgdorferi sensu stricto—and 25 STs of Borrelia, of which 12 have not been previously described at all and five have not been previously reported in Germany, were detected in 13% of I. ricinus complex ticks. Prevalence was highest in adult females (25.3%) and lowest in nymphs (11.4%). Prevalence was significantly higher in ticks from grassland (16.8%) compared to forests (11.4%).
The high level of small mammal diversity in this region of Germany seems to be reflected in a wide variety of genospecies and STs of B. burgdorferi.
Hard ticks (Ixodidae) are important vectors of zoonotic pathogens, including bacteria, protozoan parasites, and viruses. The most common tick species in Europe is the castor bean tick, Ixodes ricinus. Due to its wide spectrum of hosts and its environmental plasticity, it occurs in several types of habitats and it is important for the transmission of several zoonotic pathogens . Its immature life stages (larvae and nymphs) parasitize mostly small mammals and birds. The second most common tick species in temperate Europe is the ornate cow tick, Dermacentor reticulatus, whose juvenile developmental stages feed mostly on small mammals [2, 3]. Therefore, small mammals, predominantly rodents, are essential for the maintenance and distribution of ticks and tick-borne pathogens [4, 5].
In Europe, the most frequently reported tick-borne disease is Lyme disease, which is caused by Borrelia burgdorferi sensu lato (s.l.), a complex of gram-negative spirochetes. Borrelia burgdorferi s.l. is transmitted mostly by I. ricinus ticks, but has also occasionally been found in D. reticulatus [6,7,8].
To date, over 20 genospecies of the B. burgdorferi s.l. complex have been identified, with 11 of them occurring throughout Europe : Borrelia afzelii, Borrelia bavariensis, Borrelia bissettiae (previously known as Borrelia bissettii ), Borrelia burgdorferi sensu stricto (s.s.), Borrelia carolinensis, Borrelia finlandensis, Borrelia garinii, Borrelia kurtenbachii, Borrelia lusitaniae, Borrelia spielmanii, and Borrelia valaisiana. Most of them are known to be pathogenic for humans. I. ricinus is the main vector for the 11 European genospecies of Borrelia, but the reservoir host species range is much more diverse and genospecies specific [9, 11]. More than 40 different vertebrate species, including reptiles, birds, and mammals, may serve as reservoir hosts of Borrelia [4, 5, 12,13,14,15]. The adaptation of B. burgdorferi s.l. species to different reservoir hosts is an important driver for their spatial distribution and diversification in ecosystems . In a complex multi-host system, Borrelia burgdorferi s.l. circulation is difficult to monitor and forecast in the respective ecosystems . According to the biodiversity-driven dilution effect hypotheses, the balance between competent (reservoir) and incompetent (dilution, rescue) hosts is crucial for the risk of infection [18,19,20]. Feeding on a dilution host species may clear a pre-existing Borrelia infection in I. ricinus [21, 22]. Clearance might therefore occur more often in communities dominated by incompetent hosts. The rescue hosts mostly support the life cycle of ticks, but are capable of maintaining disease risk when reservoir host density is low.
Genotyping of Borrelia spp. is a helpful tool for the description of epidemiological and ecological patterns. Borrelia genospecies are usually further categorized into sequence types (STs) by multilocus sequence typing (MLST), which is based on eight housekeeping genes: nifS, pyrG, clpX, pepX, uvrA, rplB, cplA, and recG. These genes are characterized by slow evolutionary rates and an almost neutral variation resulting in robust phylogenetic relationships and the avoidance of potentially skewed clustering in single-locus sequence analyses [16, 23]. Factors such as habitat, geographical region, host, and vector species may have an impact on the genetic variability of certain Borrelia genospecies. The mechanisms of rapid adaptation through alterations in gene expression allow the spirochetes to adjust and survive in different host species [24,25,26]. However, there is scarce information on the variation of STs in mammals and ticks when considered on a small geographical scale, and the habitat associations of certain STs.
The aim of this study was to determine the prevalence and (geno)diversity of Borrelia burgdorferi s.l. in small mammals (21 sites) and ticks (17 sites) from different forest and grassland sites in central Germany, and to elaborate on the effects of habitat, season, and small mammal species richness on Borrelia prevalence.
The study sites (Fig. 1) were located in the Hainich Dün region surrounding the Hainich National Park, which is one of the largest continuous deciduous forests in Germany. Outside the national park the grasslands are grazed extensively, primarily by sheep, and the forests are managed by selection cutting. Due to its large area, which is protected and under a balanced system of local management, the Hainich Dün region has a high animal and plant biodiversity (https://www.biodiversity-exploratories.de/en/regions/hainich-duen/).
Small mammal and tick collection
Skin samples from small mammals trapped using snap traps at 21 sites in central Germany were available from a former study [27, 28] (Fig. 1). Trapping took place in spring, summer and autumn in 2017, 2018 and 2019. For this study, 1167 out of 1945 samples from all 21 sites, comprising 168 spring samples (range 24–146) (2018 and 2019) and 999 summer samples (range 72–639) (2017, 2018 and 2019), which represented 60% of the original material, were randomly selected.
Ticks were collected from vegetation at 17 of the 21 sites by the flagging method  once per season along 100-m2 transects in spring, summer and fall 2018 and 2019. The number of collection sites for ticks was limited to 17 due to financial restrictions. The ticks were stored in 70% ethanol and the small mammal skin samples in a freezer at − 20 °C until further examination.
At each site, small mammal trapping was performed in grasslands (I) and forests (II), and tick collection was conducted within the ecotone at the edge of forests, which comprised woodland and adjacent grassland (Fig. 2). Ticks were not collected from the grassland habitat where the small mammals were trapped (I), but were collected from the grassland (IIb) and woodland areas (IIa) of the forest habitat (II). Thus, grassland habitat of ticks refers to the grassland part of the forest habitat (IIb) and forest habitat of ticks refers to the woodland part of the forest habitat (IIa) (Fig. 2).
Sample preparation, DNA extraction and polymerase chain reaction methods
Small mammal species identification and dissection were conducted as described previously [27, 28]. Ticks were identified to species level using morphological characters in accordance with Siuda  and Estrada-Peña et al. . Before DNA extraction, the ticks were washed twice in distilled water and phosphate-buffered saline and then air dried. DNA extraction of ticks and small mammal skin samples was carried out for each sample individually using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) . The quality and quantity of the DNA samples were measured with a spectrophotometer (NanoDrop 2000c; Peqlab Biotechnologie, Erlangen, Germany). The extracted DNA was stored at − 20 °C until further examination.
For the detection of Borrelia spp. DNA, small mammal and tick samples were first screened by real-time polymerase chain reaction (PCR) targeting the p41 flagellin gene with an expected amplicon of 96 base pairs . In order to identify the Borrelia genospecies and STs of positive samples yielding a cycle threshold (Ct) value below 35, a MLST targeting the housekeeping genes nifS, pyrG, clpX, pepX, uvrA, rplB, cplA, and recG was performed, with slight modifications [14, 23]. In our experience, the amount of pathogen DNA in samples with Ct < 35 is too low for sufficient amplification by conventional PCR. For samples that did not yield PCR products for all eight MLST housekeeping genes and for which the identification of a specific ST was impossible, genospecies determination was based on recG gene sequences.
PCR products were visualized under ultraviolet (UV) light using the UVP GelSolo Simplified UV Gel Documentation System (Analytik Jena, Jena, Germany). For samples positive for all eight PCRs, all the housekeeping genes were sequenced (Eurofins Genomics, Ebersberg, Germany), with forward and reverse primers of each gene used for PCR amplification. The sequence results were analysed with Bionumerics software (version 7.6.1.; Applied Maths, Austin, TX), and compared with sequences published in GenBank using BLASTn (https://blast.ncbi.nlm.nih.gov/Blast.cgi). The aligned sequences of the eight housekeeping loci were compared with those in the online MLST database to assign allele and ST profiles. Novel STs were submitted to the curator, who allocated them consecutive numbers (IDs 3255–3267 with respective STs 986–998).
Confidence intervals (CIs; 95% CI) for the prevalences were determined by the Clopper and Pearson method using GraphPad software (Graph Pad Software, San Diego, CA). A chi-square test was used to test the independence of Borrelia prevalence for the factors Small mammal species, Tick life stage, Season, Habitat, and Collection site (GraphPad Software). To analyse B. burgdorferi s.l. prevalence in small mammals in relation to season, habitat, and small mammal species, we used a generalized linear mixed model (GLMM) with binomial error distribution using R software (version 4.1.2. for Windows; RStudio, Boston, MA) and the lme4 package . Infection status was used as a binary dependent variable (Borrelia spp. positive = 1; Borrelia spp. negative = 0). The GLMM was generated to estimate how (i) seasonality (independent binary variable: summer vs. spring), (ii) small mammal species (independent categorical variable), and (iii) habitat (independent binary variable: forest vs. grassland) affect individual infection status (dependent binary variable).
The computed similar approach was used for ticks with the dependent variable Borrelia prevalence and the independent variables Tick developmental stage (binary: adult vs. immature), Habitat (see above) and Season (three stages: autumn, spring and summer). Dermacentor reticulatus ticks were excluded from this GLMM as none of them was Borrelia positive and because the sample size was low.
For small mammals and ticks, the interaction term for the GLMM consisted of three variables with at least two levels each. Therefore, we computed marginal means using the emmeans package within R and a post hoc test to compare the effects of all independent variables separately . Additionally, a generalized linear model (lme4 package in R) was performed with a binary error distribution and a logit link function in order to determine if the number of small mammal species (species richness) per site  influences the probability of Borrelia infection in individual ticks and small mammals. The significance threshold was set at P ≤ 0.05.
Small mammal collection
The selected samples (n = 1167) were from 10 small mammal species, comprising seven rodent and three shrew species (Table 1). The most common species was the common vole (Microtus arvalis) (n = 407; 34.9%), followed by the bank vole (Clethrionomys glareolus) (n = 278; 23.8%) and the yellow-necked field mouse (Apodemus flavicollis) (n = 240; 20.6%); next were the wood mouse (Apodemus sylvaticus) (n = 108; 9.3%) and the striped field mouse (Apodemus agrarius) (n = 90; 7.7%). The common shrew (Sorex araneus) (n = 20; 1.7%), the field vole (Microtus agrestis) (n = 15; 1.3%), the water vole (Arvicola amphibius) (n = 4; 0.3%), the pygmy shrew (Sorex minutus) (n = 4; 0.3%), and the greater white-toothed shrew (Crocidura russula) (n = 1; 0.1%) were rare. The three least abundant species were excluded from the analyses due to their low numbers. Altogether, 504 (43.2%) of the small mammals were from grasslands and 663 (56.8%) from forests (Table 1). In the grasslands, the predominant species was M. arvalis (74.8%) and in the forests there were two dominant species, C. glareolus (40.1%) and A. flavicollis (33.2%). The majority of individuals were caught in the summer (n = 999; 85.6%); fewer of the small mammals were collected in the spring (n = 168; 14.4%).
In total, 1115 ticks were collected, but 21 individuals were excluded from further examination due to poor conservation of the material. The remaining 1094 ticks belonged to three species, the predominant I. ricinus (n = 998; 91.2%), followed by D. reticulatus (n = 80; 7.3%) and I. inopinatus (n = 16; 1.5%). For further analysis, data were combined for I. ricinus and I. inopinatus, and are henceforth referred to as those for ‘I. ricinus complex’ . In general, nymphs were the most abundant life stage found (Table 2). Overall, almost twice as many ticks were collected from forests (n = 717; 65.5%) compared to grasslands (n = 377; 34.5%). However, almost all D. reticulatus ticks were found in grasslands (74/80; 92.5%). Most of the ticks were collected in spring (n = 822, 75.1%; I. ricinus complex, n = 785, 77.4%; D. reticulatus, n = 37, 46.3%), followed by summer (n = 203, 18.6%; I. ricinus complex, n = 202, 19.9%; D. reticulatus, n = 1, 1.3%), and fall (n = 69, 6.3%; I. ricinus complex, n = 27, 2.7%; D. reticulatus, n = 42, 52.5%).
Borrelia spp. prevalence and genospecies in small mammals
Borrelia spp. DNA was detected in eight out of the 10 small mammal species, with an overall prevalence of 6.2% (n = 72; 95% CI 4.9–7.7) (Table 3). Among the small mammal genera, prevalence was lower in Apodemus compared to Microtus, Clethrionomys and Sorex (χ2 = 29.122, df = 3, P < 0.001). Borrelia prevalence did not differ between spring and summer (see GLMM in Additional file 1: Table S1). The individual infection probability was not affected by the small mammal species richness per site (Additional file 1: Table S2). The GLMM showed that there were no differences in prevalence between small mammals from grassland (I) and forest (II) (P = 0.795), and the bias in prevalence must have been confounded by infection rates per small mammal species, as the small mammal species composition differed between grassland (I) and forest (II), with M. arvalis having a high prevalence (9.6%) and being the most abundant small mammal species in the grassland (I), while A. flavicollis (prevalence 1.7%) was mostly present in the forest (II) (P = 0.478). In general, prevalence differed among sites, and ranged from 0 to 20.5% (χ2 = 68.392, df = 20, P < 0.001). This was due to differences in prevalence among grassland sites (I) (0–30.8%; χ2 = 59.381, df = 19, P < 0.001), while in forest sites (II) there was less difference between prevalences (0–13%; χ2 = 27.883, df = 20, P = 0.112).
Out of 72 positive samples identified by real-time PCR, eight yielded a Ct value < 35 and were sequenced (one from Cl glareolus, four from M. arvalis and three from S. araneus). The analysis revealed the presence of two genospecies (Additional file 1: Table S5; Table 4): B. afzelii in seven individuals of three host species, and B. garinii in one M. arvalis. The recG gene sequences showed 100% identity to existing GenBank entries (Table 4).
Borrelia spp. prevalence and genospecies in ticks
Borrelia spp. DNA was detected in 132 out of 1094 ticks (12.1%; 95% CI 10.3–14.1). All of the positive samples originated from I. ricinus complex ticks (n = 1014) (13%; 95% CI 11.1–15.2). The prevalence was higher in females (19/75, 25.3%; 95% CI 16.8–36.3) than in males (19/111, 17.1%; 95% CI 11.2–25.3), and lowest in nymphs (94/828, 11.4%; 95% CI 9.4–14.5) (χ2 = 13.721, df = 2, P = 0.001). Ixodes ricinus complex ticks collected in the grasslands (IIb) were significantly more often infected (51/303, 16.8%; 95% CI 13–21.5) than those from forests (IIa) (81/711, 11.4%; 95% CI 9.3–14) (see GLMM in Additional file 1: Table S3). More detailed analyses revealed that a higher prevalence in adults compared to nymphs was true for grasslands (P = 0.015) but not for forests (P = 0.849) (Additional file 1: Table S3). The GLMM also revealed that seasonal fluctuation of prevalence was a confounding factor (P = 0.323–0.989) and could be explained by a seasonal effect of the life stage composition of I. ricinus. Adults were significantly more often positive than nymphs (P = 0.0155), and in spring the ratio of nymphs to adults was 4:1, while in summer it was 13:1. The prevalence among sites varied significantly, from 0 to 28.6% (χ2 = 39.013, df = 16, P = 0.001); it was significant for grasslands (IIb), where the prevalence range was higher (0–75%; χ2 = 41.205, df = 16; p < 0.001), but not for forests (IIa) (0–15.8; χ2 = 23.952, df = 16, P = 0.09). The generalized linear model suggested that the probability of infection in ticks was positively correlated to the number of small mammal species per site (Additional file 1: Table S4).
Out of 132 positive tick samples, 69 recG gene PCR products were sequenced; this revealed the presence of six genospecies (Additional file 1: Table S5 and Table 4): B. afzelii (n = 37), B. garinii (n = 16), B. valaisiana (n = 9), B. lusitaniae (n = 3), B. burgdorferi s.s. (n = 2), and B. spielmanii (n = 2). All of the genospecies were detected in ticks collected from grasslands (IIb) and forests (IIa). However, the genospecies composition differed between habitats, with an almost equal genospecies distribution of B. afzelii and B. garinii (44% vs. 35%) in the grasslands and a 4:1 ratio of B. afzelii to B. garinii in forests (60% vs. 15%). The recG gene sequences showed very high (97.1–100%) identity to GenBank entries of the corresponding reference strains (Table 4).
MLST analyses were performed for 51 Borrelia-positive samples obtained from 46 I. ricinus complex ticks and five small mammals. The MLST results were ambiguous for 16 samples, and thus STs could not be assigned for these. In total, the remaining 35 samples revealed 28 STs, 25 in ticks and three in small mammals, indicating a genetically diverse Borrelia population (Additional file 1: Table S5). Only three STs were detected more than once, and all of these were from tick samples. The most common of these was ST 347, belonging to B. afzelii, which was found in five ticks from three different sites across both forest habitats (IIa and IIb), followed by ST 251 of B. garinii, which was detected in three ticks from two sites and habitats. The third ST that was found multiple times was B. garinii ST 187, but only in ticks from two different spots in forests (IIa). The three STs detected in common shrews and a common vole were not found in ticks (Additional file 1: Table S5). In grasslands, 10 different STs were found in ticks (IIb) and two in small mammals (I). In forest sites, 17 STs were identified in ticks (IIa) and one in small mammals (II). The highest diversity of STs was noted for B. afzelii (n = 16), followed by B. garinii (n = 9), B. valaisiana (n = 2) and B. burgdorferi s.s. (n = 1). No STs were assigned for B. spielmanii or B. lusitaniae.
Small mammals are important reservoir hosts for B. burgdorferi s.l. [4, 12, 13, 37, 38]. In our study, 10 small mammal species were examined, almost all of which are recognized reservoir hosts of B. burgdorferi [4, 5, 12, 13]. Borrelia spp. DNA was found in eight out of the 10 species, but not in C. russula and S. minutus, most probably due to their low sample size. However, these two shrew species were found to be Borrelia positive in previous studies [4, 39].
The mean prevalence in small mammals in this study was 6.2%, which is in line with the results of previous European studies that showed prevalences of 1.2–24% [40,41,42,43,44,45,46]. However, it was lower than in other studies with a lower number of small mammal species (31.2–49.1%) [14, 31]. Although Apodemus had moderate to high prevalence (e.g. 25.4–47.5%) in other studies [14, 31], in our study the lowest Borrelia spp. prevalence was in this genus. There was no seasonal pattern for Borrelia prevalence in small mammals, most probably because the samples derived only from spring and summer sampling. A similar absence of a seasonal pattern was also observed in a study from Spain .
Dermacentor reticulatus ticks have been found to carry Borrelia in Germany [47, 48], but Borrelia-negative D. reticulatus ticks have also been reported in several studies from Germany [14, 31, 49,50,51]. In our study, none of the D. reticulatus ticks (n = 80) was positive for B. burgdorferi s.l., although at two spots Borrelia DNA was detected in sympatric I. ricinus complex ticks. These findings support the view that D. reticulatus is of minor importance in the natural life cycle of the B. burgdorferi species complex .
Borrelia prevalence in the I. ricinus complex in the current study was 13%, which is in agreement with the results of European meta-analyses from the last decades, i.e. 13.7% (1984–2003) and 12.3% (2010–2016) [52, 53]. The level of infection was higher in adults (20.4%) than in nymphs (11.1%), as also found in other studies (18.6 vs. 10.1% ; and 14.9 vs. 11.8% ), and is a result of Borrelia transstadial transmission and a higher number of blood meals on reservoir hosts for adults. As seasonality was not an important explanatory variable in the GLMM analysis, the lack of significant differences in prevalence between the seasons could be explained by the different seasonal ratios and activity levels of immature and adult ticks.
Borrelia prevalence in ticks and small mammals differed significantly among grassland sites, which suggests that local factors may impact the maintenance of pathogen circulation. Such factors could be land use that influences the abundance of the potential hosts—not only small but also medium- and large-sized mammals, as well as birds. Interestingly, the prevalence was lower in ticks from forests (IIa) than from grasslands (IIb), possibly because of a higher diversity of wildlife in forests [54, 55].
Our analysis showed that Borrelia spp. prevalence in ticks is positively correlated to the number of small mammal species per site: the more reservoir host species, the higher the infection probability in ticks. This seems to be in line with the dilution effect hypothesis [56, 57], which assumes that reservoir competence should vary between host species. In addition, in our study, almost all of the tested animals were reservoir hosts (e.g. ) with a different prevalence of Borrelia; therefore, an amplification effect was observed, meaning the more reservoir species the higher the prevalence in vectors. The number of small mammal species per site did not influence Borrelia prevalence in the hosts themselves.
Six Borrelia genospecies were identified in the present study. B. afzelii was predominant in ticks (37/69) and small mammals (7/8), followed by B. garinii [which was found in ticks (16/69) and small mammals (1/8)]; similar results have been previously reported [31, 58] for Germany and elsewhere in Europe . The other genospecies—B. valaisiana (6/69), B. lusitaniae (3/69), B. burgdorferi s.s. (2/69) and B. spielmanii (2/69)—were found only in ticks, and have been described before in ticks from Germany at similar frequencies [59, 60]. The diversity of the detected genospecies was higher [14, 31, 58] or similar to that of other studies carried out in Germany, Belarus, and Romania [59,60,61,62]. The reason for this might be the association of certain genospecies with the abundance of their respective reservoirs. Rodents are considered the main reservoir hosts for following Borrelia genospecies: B. afzelii, B. lusitaniae, B. burgdorferi s.s. and B. spielmanii, while small birds for B. garinii, B. valaisiana as well as B. burgdorferi s.s. Some years ago, B. garinii OspA serotype 4 strains were reclassified as B. bavariensis, as they use rodents as reservoir hosts instead of birds [11, 63,64,65]. However, as we did not sequence ospA, which could have been used to discriminate the genospecies, we could only compare the sample from M. arvalis (KS19/1709) to GenBank data; our recG sequence was a 100% match to B. garinii (GenBank no. CP028861). Additionally, lizards are recognized reservoirs for B. lusitaniae, B. valaisiana and B. garinii, and larger mammals for B. burgdorferi s.s.
MLST has been performed for Borrelia spp. genotyping in ticks from 21 European countries, but for vertebrate hosts only for six European countries (pubmlst.org, accessed on 1 Feb 2022). Our MLST analyses for tick- and small mammal-derived samples revealed 15 different STs that have been previously published. Most of these STs were previously detected in humans and/or ticks from European countries. Six of these STs have been reported in ticks from eight European countries, but never from Germany. One ST (ST 467) that has not been previously reported in ticks was previously reported in a person from France . ST 338, which was detected in a small mammal in the current study, has been formerly described in ticks from Austria (pubmlst.org). The other ST found in the small mammals (ST 348) has been previously described from a person from Germany, but has never previously been reported for other animals. Additionally, we discovered 13 new STs that have never been reported before (STs 986–998): 12 from ticks and one from a shrew. The STs found in the ticks did not overlap with those detected in the shrews and voles, suggesting that the ticks acquired the pathogens from other reservoir hosts. However, it should be borne in mind that the number of STs determined in the small mammals was very low. This high variety of STs (at least 28) indicates a very diverse Borrelia population, even at individual sites (e.g. E1, E3, UH1). The quantity and variety of Borrelia STs reflected the generally high plant and animal biodiversity of this region, which has many potential host species [20, 57, 67,68,69]. However, as we tested only small mammals, other potential host groups should be investigated in future research.
This study was conducted in an area of high biodiversity with numerous species of animals and plants. Our results showed a high species richness of small mammals (n = 10), ticks (n = 3) and Borrelia genospecies (n = 6). Likewise, the number of detected Borrelia STs (≥ 28) was also very high. The genetic variation of B. afzelii in ticks was not aligned with that in shrews and voles, which points to other vertebrate host species contributing to the life cycle of this pathogen. The prevalence of Borrelia in ticks seemed to be influenced by habitat, as in grassland habitats adjacent to forest the prevalence was higher than in the forested habitats. Even though the prevalence of Borrelia in ticks and small mammals was moderate, all the genospecies detected have zoonotic potential and constitute a disease risk for humans. Therefore, the attention of public health authorities should be drawn to the results presented here for this region of Germany.
Availability of data and materials
The datasets generated and analysed during the current study are available from the corresponding author.
Generalized linear mixed model
Multilocus sequence typing
Polymerase chain reaction
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.
Drehmann M, Springer A, Lindau A, Fachet K, Mai S, Thoma D, et al. The spatial distribution of Dermacentor ticks (Ixodidae) in Germany-evidence of a continuing spread of Dermacentor reticulatus. Front Vet Sci. 2020;7:578220.
Földvári G, Široký P, Szekeres S, Majoros G, Sprong H. Dermacentor reticulatus: a vector on the rise. Parasit Vectors. 2016;9:314.
Gern L, Estrada-Peña A, Frandsen F, Gray JS, Jaenson T, Jongejan F, et al. European reservoir hosts of Borrelia burgdorferi sensu lato. Zentralbl Bakteriol. 1998;287:196–204.
Brunner JL, LoGiudice K, Ostfeld RS. Estimating reservoir competence of Borrelia burgdorferi hosts: prevalence and infectivity, sensitivity, andspecificity. J Med Entomol. 2008;45:139–47.
Rizzoli A, Hauffe H, Carpi G, Vourc HG, Neteler M, Rosa R. Lyme borreliosis in Europe. Euro Surveill. 2011. https://doi.org/10.2807/ese.16.27.19906-en .
Bonnet S, La De Fuente J, Nicollet P, Liu X, Madani N, Blanchard B, et al. Prevalence of tick-borne pathogens in adult Dermacentor spp. ticks from nine collection sites in France. Vector Borne Zoonotic Dis. 2013;13:226–36.
Mierzejewska EJ, Pawełczyk A, Radkowski M, Welc-Falęciak R, Bajer A. Pathogens vectored by the tick, Dermacentor reticulatus, in endemic regions and zones of expansion in Poland. Parasit Vectors. 2015;8:490.
Waindok P, Schicht S, Fingerle V, Strube C. Lyme borreliae prevalence and genospecies distribution in ticks removed from humans. Ticks Tick Borne Dis. 2017;8:709–14.
Margos G, Lane RS, Fedorova N, Koloczek J, Piesman J, Hojgaard A, et al. Borrelia bissettiae sp. nov. and Borrelia californiensis sp nov. prevail in diverse enzootic transmission cycles. Int J Syst Evol Microbiol. 2016;66:1447–52.
Kurtenbach K, de Michelis S, Etti S, Schäfer SM, Sewell H-S, Brade V, et al. Host association of Borrelia burgdorferi sensu lato—the key role of host complement. Trends Microbiol. 2002;10:74–9.
Humair PF, Turrian N, Aeschilimann A, Gern L. Borrelia burgdorferi in a focus of Lyme borreliosis: epizootiologic contribution of small mammals. Folia Parasitol (Praha). 1993;40:65–70.
Kurtenbach K, Dizij A, Seitz HM, Margos G, Moter SE, Kramer MD, et al. Differential immune responses to Borrelia burgdorferi in European wild rodent species influence spirochete transmission to Ixodes ricinus L. (Acari: Ixodidae). Infect Immun. 1994;62:5344–52.
Obiegala A, Król N, Oltersdorf C, Nader J, Pfeffer M. The enzootic life-cycle of Borrelia burgdorferi (sensu lato) and tick-borne rickettsiae: an epidemiological study on wild-living small mammals and their ticks from Saxony, Germany. Parasit Vectors. 2017;10:115.
Matuschka FR, Fischer P, Heiler M, Richter D, Spielman A. Capacity of European animals as reservoir hosts for the Lyme disease spirochete. J Infect Dis. 1992;165:479–83.
Margos G, Vollmer SA, Cornet M, Garnier M, Fingerle V, Wilske B, et al. A new Borrelia species defined by multilocus sequence analysis of housekeeping genes. Appl Environ Microbiol. 2009;75:5410–6.
Buhnerkempe MG, Roberts MG, Dobson AP, Heesterbeek H, Hudson PJ, Lloyd-Smith JO. Eight challenges in modelling disease ecology in multi-host, multi-agent systems. Epidemics. 2015;10:26–30.
Keesing F, Brunner J, Duerr S, Killilea M, Logiudice K, Schmidt K, et al. Hosts as ecological traps for the vector of Lyme disease. Proc Biol Sci. 2009;276:3911–9.
LoGiudice K, Duerr STK, Newhouse MJ, Schmidt KA, Killilea ME, Ostfeld RS. Impact of host community composition on Lyme disease risk. Ecology. 2008;89:2841–9.
LoGiudice K, Ostfeld RS, Schmidt KA, Keesing F. The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. Proc Natl Acad Sci USA. 2003;100:567–71.
Matuschka FR, Spielman A. Loss of Lyme disease spirochetes from Ixodes ricinus ticks feeding on European blackbirds. Exp Parasitol. 1992;74:151–8.
Richter D, Matuschka FR. Elimination of Lyme disease spirochetes from ticks feeding on domestic ruminants. Appl Environ Microbiol. 2010;76:7650–2.
Wang G, Liveris D, Mukherjee P, Jungnick S, Margos G, Schwartz I. Molecular typing of Borrelia burgdorferi. Curr Protoc Microbiol. 2014. https://doi.org/10.1002/9780471729259.mc12c05s34 .
Miller CL, Karna SLR, Seshu J. Borrelia host adaptation regulator (BadR) regulates rpoS to modulate host adaptation and virulence factors in Borrelia burgdorferi. Mol Microbiol. 2013;88:105–24.
Anguita J, Hedrick MN, Fikrig E. Adaptation of Borrelia burgdorferi in the tick and the mammalian host. FEMS Microbiol Rev. 2003;27:493–504.
Samuels DS. Gene regulation in Borrelia burgdorferi. Annu Rev Microbiol. 2011;65:479–99.
Jeske K, Weber S, Pfaff F, Imholt C, Jacob J, Beer M, et al. Molecular detection and characterization of the first cowpox virus isolate derived from a bank vole. Viruses. 2019. https://doi.org/10.3390/v11111075 .
Jeske K, Jacob J, Drewes S, Pfeffer M, Heckel G, Ulrich RG, et al. Hantavirus-Leptospira coinfections in small mammals from central Germany. Epidemiol Infect. 2021;149:e97.
Siuda K. Kleszcze polski (Acari: Ixodida): systematyka i rozmieszczenie. Warszawa: Polskie Towarzystwo Parazytologiczne; 1993. (in Polish).
Estrada-Peña A, Mihalca AD, Petney T, editors. Ticks of Europe and North Africa: a guide to species identification. Cham: Springer; 2017.
Galfsky D, Król N, Pfeffer M, Obiegala A. Long-term trends of tick-borne pathogens in regard to small mammal and tick populations from Saxony, Germany. Parasit Vectors. 2019;12:131.
Schwaiger M, Péter O, Cassinotti P. Routine diagnosis of Borrelia burgdorferi (sensu lato) infections using a real-time PCR assay. Clin Microbiol Infect. 2001;7:461–9.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015. https://doi.org/10.18637/jss.v067.i01 .
Lenth R, Singmann H, Love J, Buerkner P, Herve M. Emmeans: estimated marginal means, aka least-squares means: R package version; 2018.
Moore JC. Diversity, taxonomic versus functional. In: Moore JC, editor. Encyclopedia of biodiversity. 2nd ed. Amsterdam: Academic press; 2013. p. 648–56.
Estrada-Peña A, Nava S, Petney T. Description of all the stages of Ixodes inopinatus n. sp. (Acari: Ixodidae). Ticks Tick Borne Dis. 2014;5:734–43.
de Boer R, Hovius KE, Nohlmans M, Gray JS. The woodmouse (Apodemus sylvaticus) as a reservoir of tick-transmitted spirochetes (Borrelia burgdorferi) in the Netherlands. Zentralbl Bakteriol. 1993;279:404–16.
Kurtenbach K, Kampen H, Dizij A, Arndt S, Seitz HM, Schaible UE, et al. Infestation of rodents with larval Ixodes ricinus (Acari: Ixodidae) is an important factor in the transmission cycle of Borrelia burgdorferi s.l. in German woodlands. J Med Entomol. 1995;32:807–17.
Barandika JF, Hurtado A, García-Esteban C, Gil H, Escudero R, Barral M, et al. Tick-borne zoonotic bacteria in wild and domestic small mammals in northern Spain. Appl Environ Microbiol. 2007;73:6166–71.
Gryczyńska A, Gortat T, Kowalec M. Urban rodent reservoirs of Borrelia spp. in Warsaw, Poland. Epidemiol Infect. 2018;146:589–93.
Siński E, Pawełczyk A, Bajer A, Behnke J. Abundance of wild rodents, ticks and environmental risk of Lyme borreliosis: a longitudinal study in an area of Mazury Lakes district of Poland. Ann Agric Environ Med. 2006;13:295–300.
Kybicová K, Kurzová Z, Hulínská D. Molecular and serological evidence of Borrelia burgdorferi sensu lato in wild rodents in the Czech Republic. Vector Borne Zoonotic Dis. 2008;8:645–52.
Marsot M, Sigaud M, Chapuis JL, Ferquel E, Cornet M, Vourc’h G. Introduced Siberian chipmunks (Tamias sibiricus barberi) harbor more-diverse Borrelia burgdorferi sensu lato genospecies than native bank voles (Myodes glareolus). Appl Environ Microbiol. 2011;77:5716–21.
Pascucci I, Di Domenico M, Dall’Acqua F, Sozio G, Cammà C. Detection of Lyme disease and Q fever agents in wild rodents in Central Italy. Vector Borne Zoonotic Dis. 2015;15:404–11.
Gray JS, Kirstein F, Robertson JN, Stein J, Kahl O. Borrelia burgdorferi sensu lato in Ixodes ricinus ticks and rodents in a recreational park in south-western Ireland. Exp Appl Acarol. 1999;23:717–29.
Khanakah G, Kocianová E, Vyrosteková V, Rehácek J, Kundi M, Stanek G. Seasonal variations in detecting Borrelia burgdorferi sensu lato in rodents from north eastern Austria. Wien Klin Wochenschr. 2006;118:754–8.
Kahl O, Janetzki C, Gray JS, Stein J, Bauch RJ. Tick infection rates with Borrelia: Ixodes ricinus versus Haemaphysalis concinna and Dermacentor reticulatus in two locations in eastern Germany. Med Vet Entomol. 1992;6:363–6.
Kohn M, Krücken J, McKay-Demeler J, Pachnicke S, Krieger K, von Samson-Himmelstjerna G. Dermacentor reticulatus in Berlin/Brandenburg (Germany): activity patterns and associated pathogens. Ticks Tick Borne Dis. 2019;10:191–206.
Richter D, Kohn C, Matuschka FR. Absence of Borrelia spp., Candidatus Neoehrlichia mikurensis, and Anaplasma phagocytophilum in questing adult Dermacentor reticulatus ticks. Parasitol Res. 2013;112:107–11.
Opalińska P, Wierzbicka A, Asman M. The PCR and nested PCR detection of Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Babesia microti in Dermacentor reticulatus F. collected in a new location in Poland (Trzciel, western Poland). Acta Parasitol. 2016;61:849–54.
Kovryha N, Tsyhankova A, Zelenuchina O, Mashchak O, Terekhov R, Rogovskyy AS. Prevalence of Borrelia burgdorferi and Anaplasma phagocytophilum in ixodid ticks from southeastern Ukraine. Vector Borne Zoonotic Dis. 2021;21:242–6.
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, Hönig V, Růžek D, Grubhoffer L, Rego R. Europe-wide meta-analysis of Borrelia burgdorferi sensu lato prevalence in questing Ixodes ricinus ticks. Appl Environ Microbiol. 2017. https://doi.org/10.1128/AEM.00609-17 .
Burton V, Moseley D, Brown C, Metzger MJ, Bellamy P. Reviewing the evidence base for the effects of woodland expansion on biodiversity and ecosystem services in the United Kingdom. For Ecol Manage. 2018;430:366–79.
Ciccarese L, Mattsson A, Pettenella D. Ecosystem services from forest restoration: thinking ahead. New For. 2012;43:543–60.
Ostfeld RS, Keesing F, Eviner V. Infectious disease ecology: effects of ecosystems on disease and of disease on ecosystems. Princeton: Princeton University Press; 2008.
Ostfeld RS, Keesing F. Biodiversity series: the function of biodiversity in the ecology of vector-borne zoonotic diseases. Can J Zool. 2000;78:2061–78.
Răileanu C, Tauchmann O, Vasić A, Wöhnke E, Silaghi C. Borrelia miyamotoi and Borrelia burgdorferi (sensu lato) identification and survey of tick-borne encephalitis virus in ticks from north-eastern Germany. Parasit Vectors. 2020;13:106.
Springer A, Jordan D, Glass A, Kahl O, Fingerle V, Girl P, et al. Borrelia infections in ageing ticks: relationship with morphometric age ratio in field-collected Ixodes ricinus nymphs. Microorganisms. 2022. https://doi.org/10.3390/microorganisms10010166 .
Springer A, Raulf M-K, Fingerle V, Strube C. Borrelia prevalence and species distribution in ticks removed from humans in Germany, 2013–2017. Ticks Tick Borne Dis. 2020;11:101363.
Reye AL, Stegniy V, Mishaeva NP, Velhin S, Hübschen JM, Ignatyev G, et al. Prevalence of tick-borne pathogens in Ixodes ricinus and Dermacentor reticulatus ticks from different geographical locations in Belarus. PLoS ONE. 2013;8:e54476.
Raileanu C, Moutailler S, Pavel I, Porea D, Mihalca AD, Savuta G, et al. Borrelia diversity and co-infection with other tick borne pathogens in ticks. Front Cell Infect Microbiol. 2017;7:36.
Rudenko N, Golovchenko M, Grubhoffer L, Oliver JH. Updates on Borrelia burgdorferi sensu lato complex with respect to public health. Ticks Tick Borne Dis. 2011;2:123–8.
Huegli D, Hu CM, Humair P-F, Wilske B, Gern L. Apodemus species mice are reservoir hosts of Borrelia garinii OspA serotype 4 in Switzerland. J Clin Microbiol. 2002;40:4735–7.
Margos G, Fingerle V, Reynolds S. Borrelia bavariensis: vector switch, niche invasion, and geographical spread of a tick-borne bacterial parasite. Front Ecol Evol. 2019. https://doi.org/10.3389/fevo.2019.00401 .
Gallais F, de Martino SJ, Sauleau EA, Hansmann Y, Lipsker D, Lenormand C, et al. Multilocus sequence typing of clinical Borreliella afzelii strains: population structure and differential ability to disseminate in humans. Parasit Vectors. 2018;11:374.
Ostfeld RS, Keesing F. Biodiversity and disease risk: the case of Lyme disease. Conserv Biol. 2000;14:722–8.
Ostfeld RS, Keesing F. Effects of host diversity on infectious disease. Annu Rev Ecol Evol Syst. 2012;43:157–82.
Dantas-Torres F. Climate change, biodiversity, ticks and tick-borne diseases: the butterfly effect. Int J Parasitol Parasites Wildl. 2015;4:452–61.
The publication fee of this paper was covered by Leipzig University within the program of Open Access Publishing in the framework of the 16th CVBD World Forum Symposium. The authors would like to thank Dana Rüster for her excellent technical assistance and lab work, and Volker Fingerle for providing the positive controls.
The small mammal trapping was commissioned and funded by the Federal Environment Agency within the Environment Research Plan of the German Federal Ministry for the Environment, Nature Conservation, Building, and Nuclear Safety (to JJ; grant number 3716484310) and partially supported by the Bundesministerium für Bildung und Forschung through the research net Zoonotic Infectious Diseases (consortium RoBoPub; grant number 01KI1721A/01KI2004A to RGU and 01KI1721A/01KI2004E to JJ). This study was also partially funded by the Veterinarian Faculty of the University of Leipzig through the Young Scientists Award (to NK).
Ethics approval and consent to participate
Permits for small mammal trapping were approved by the Thuringia State Office of Consumer Protection (protocol code 22-2684-04-15-105/16).
Consent for publication
All the authors agreed to the publishing of the results.
The authors declare that they have no conflicts of interest. The funders were not involved in the study design, sample collection, or analyses.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Results of a GLMM with a binomial error distribution for effects of habitat, small mammal species and season on Borrelia spp. infection probability in small mammals. Table S2. Generalized linear model for Borrelia spp. infection probability in small mammals according to the number of small mammal species per site. Table S3. Results of a GLMM with a binomial error distribution for effects of habitat, tick developmental stage and seasonality on Borrelia spp. infection probability in ticks. Table S4. Generalized linear model results for Borrelia spp. infection probability in ticks according to the number of small mammal species detected per site. Table S5. Results of MLST analyses of Borrelia spp. samples from ticks and small mammals.
About this article
Cite this article
Król, N., Obiegala, A., Imholt, C. et al. Diversity of Borrelia burgdorferi sensu lato in ticks and small mammals from different habitats. Parasites Vectors 15, 195 (2022). https://doi.org/10.1186/s13071-022-05326-3
- Sequence type
- Multilocus sequence typing