Skip to main content

Diversity of Borrelia burgdorferi sensu lato in ticks and small mammals from different habitats



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.

Graphical abstract


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 [1]. 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 [9]: Borrelia afzelii, Borrelia bavariensis, Borrelia bissettiae (previously known as Borrelia bissettii [10]), 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 [16]. In a complex multi-host system, Borrelia burgdorferi s.l. circulation is difficult to monitor and forecast in the respective ecosystems [17]. 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.


Study area

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 (

Fig. 1
figure 1

a Location of the study area (inset) in Germany (, and b locations of the 21 collection sites within the study area (source, with the authors’ own modifications)

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 [28] 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).

Fig. 2
figure 2

Photograph of study sites showing paired forest and grassland habitats [grassland habitat of small mammals (I), forest habitat of small mammals (II), forest habitat of ticks (IIa), grassland habitat of ticks (IIb)] (source Google Earth with the authors’ own modifications)

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 [29] and Estrada-Peña et al. [30]. 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) [31]. 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 [32]. 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 ( 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).

Statistical analysis

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 [33]. 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 [34]. 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 [35] 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%).

Table 1 Numbers of tested individuals per small mammal species and habitat, and in total

Tick collection

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’ [36]. 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%).

Table 2 Collected ticks

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 Sorex2 = 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).

Table 3 Borrelia prevalence in small mammals per habitat type and in ticks per grassland and forested parts of the forest habitat, and in total

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).

Table 4 Sequencing results based on the recG gene in tick and small mammal samples that did not yield multilocus sequence typing products

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 analysis

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 [39].

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 [3].

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% [52]; and 14.9 vs. 11.8% [53]), 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. [4]) 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 [53]. 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 (, 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 [66]. ST 338, which was detected in a small mammal in the current study, has been formerly described in ticks from Austria ( 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.



Confidence interval


Cycle threshold


Generalized linear mixed model


Multilocus sequence typing


Polymerase chain reaction


Sensu lato


Sensu stricto


Sequence types


  1. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Földvári G, Široký P, Szekeres S, Majoros G, Sprong H. Dermacentor reticulatus: a vector on the rise. Parasit Vectors. 2016;9:314.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 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.

    Article  CAS  PubMed  Google Scholar 

  5. 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.

    Article  PubMed  Google Scholar 

  6. Rizzoli A, Hauffe H, Carpi G, Vourc HG, Neteler M, Rosa R. Lyme borreliosis in Europe. Euro Surveill. 2011. .

    Article  PubMed  Google Scholar 

  7. 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.

    Article  CAS  PubMed  Google Scholar 

  8. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 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.

    Article  PubMed  Google Scholar 

  10. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  CAS  PubMed  Google Scholar 

  12. 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.

    CAS  Google Scholar 

  13. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  CAS  PubMed  Google Scholar 

  16. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. 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.

    Article  PubMed  Google Scholar 

  18. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. 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.

    Article  PubMed  Google Scholar 

  20. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Matuschka FR, Spielman A. Loss of Lyme disease spirochetes from Ixodes ricinus ticks feeding on European blackbirds. Exp Parasitol. 1992;74:151–8.

    Article  CAS  PubMed  Google Scholar 

  22. Richter D, Matuschka FR. Elimination of Lyme disease spirochetes from ticks feeding on domestic ruminants. Appl Environ Microbiol. 2010;76:7650–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Wang G, Liveris D, Mukherjee P, Jungnick S, Margos G, Schwartz I. Molecular typing of Borrelia burgdorferi. Curr Protoc Microbiol. 2014. .

    Article  PubMed  PubMed Central  Google Scholar 

  24. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Anguita J, Hedrick MN, Fikrig E. Adaptation of Borrelia burgdorferi in the tick and the mammalian host. FEMS Microbiol Rev. 2003;27:493–504.

    Article  CAS  PubMed  Google Scholar 

  26. Samuels DS. Gene regulation in Borrelia burgdorferi. Annu Rev Microbiol. 2011;65:479–99.

    Article  CAS  PubMed  Google Scholar 

  27. 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. .

    Article  PubMed  PubMed Central  Google Scholar 

  28. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Siuda K. Kleszcze polski (Acari: Ixodida): systematyka i rozmieszczenie. Warszawa: Polskie Towarzystwo Parazytologiczne; 1993. (in Polish).

    Google Scholar 

  30. Estrada-Peña A, Mihalca AD, Petney T, editors. Ticks of Europe and North Africa: a guide to species identification. Cham: Springer; 2017.

    Google Scholar 

  31. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 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.

    Article  CAS  PubMed  Google Scholar 

  33. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015. .

    Article  Google Scholar 

  34. Lenth R, Singmann H, Love J, Buerkner P, Herve M. Emmeans: estimated marginal means, aka least-squares means: R package version; 2018.

    Google Scholar 

  35. Moore JC. Diversity, taxonomic versus functional. In: Moore JC, editor. Encyclopedia of biodiversity. 2nd ed. Amsterdam: Academic press; 2013. p. 648–56.

    Chapter  Google Scholar 

  36. 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.

    Article  PubMed  Google Scholar 

  37. 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.

    Article  PubMed  Google Scholar 

  38. 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.

    Article  CAS  PubMed  Google Scholar 

  39. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Gryczyńska A, Gortat T, Kowalec M. Urban rodent reservoirs of Borrelia spp. in Warsaw, Poland. Epidemiol Infect. 2018;146:589–93.

    Article  CAS  PubMed  Google Scholar 

  41. 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.

    PubMed  Google Scholar 

  42. 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.

    Article  PubMed  Google Scholar 

  43. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 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.

    Article  CAS  PubMed  Google Scholar 

  46. 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.

    Article  PubMed  Google Scholar 

  47. 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.

    Article  CAS  PubMed  Google Scholar 

  48. 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.

    Article  PubMed  Google Scholar 

  49. 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.

    Article  PubMed  Google Scholar 

  50. 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.

    Article  CAS  PubMed  Google Scholar 

  51. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  52. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. 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. .

    Article  PubMed  PubMed Central  Google Scholar 

  54. 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.

    Article  Google Scholar 

  55. Ciccarese L, Mattsson A, Pettenella D. Ecosystem services from forest restoration: thinking ahead. New For. 2012;43:543–60.

    Article  Google Scholar 

  56. Ostfeld RS, Keesing F, Eviner V. Infectious disease ecology: effects of ecosystems on disease and of disease on ecosystems. Princeton: Princeton University Press; 2008.

    Google Scholar 

  57. 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.

    Article  Google Scholar 

  58. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. 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. .

    Article  PubMed  PubMed Central  Google Scholar 

  60. 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.

    Article  PubMed  Google Scholar 

  61. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  64. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. 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. .

    Article  Google Scholar 

  66. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Ostfeld RS, Keesing F. Biodiversity and disease risk: the case of Lyme disease. Conserv Biol. 2000;14:722–8.

    Article  Google Scholar 

  68. Ostfeld RS, Keesing F. Effects of host diversity on infectious disease. Annu Rev Ecol Evol Syst. 2012;43:157–82.

    Article  Google Scholar 

  69. Dantas-Torres F. Climate change, biodiversity, ticks and tick-borne diseases: the butterfly effect. Int J Parasitol Parasites Wildl. 2015;4:452–61.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


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).

Author information

Authors and Affiliations



Conceptualization, CI, JJ, RGU, MP, AO, NK; methodology, CI, AO, NK; sample collection, CI, KJ, AO, NK; validation, NK, AO, CI; formal analysis, NK, AO, ES; investigation, NK, AO, CA, ES, ZRS; data curation, NK; writing—original draft preparation, NK and AO; writing—review and editing, CI, JJ, RGU, KJ, ZRS, MP; visualization, NK and AO. All the authors have read and approved the final version of the manuscript.

Corresponding author

Correspondence to Martin Pfeffer.

Ethics declarations

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.

Competing interests

The authors declare that they have no conflicts of interest. The funders were not involved in the study design, sample collection, or analyses.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

. 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.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: