Next Article in Journal
Isolation of Bacillus subtilis and Bacillus pumilus with Anti-Vibrio parahaemolyticus Activity and Identification of the Anti-Vibrio parahaemolyticus Substance
Next Article in Special Issue
Study of Viral Coinfection of the Ixodes persulcatus Ticks Feeding on Humans in a Natural Focus of the South of the Far East
Previous Article in Journal
Potential Oral Microbial Markers for Differential Diagnosis of Crohn’s Disease and Ulcerative Colitis Using Machine Learning Models
Previous Article in Special Issue
The Effect of Rickettsia bellii on Anaplasma marginale Infection in Dermacentor andersoni Cell Culture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial and Temporal Variability in Prevalence Rates of Members of the Borrelia burgdorferi Species Complex in Ixodes ricinus Ticks in Urban, Agricultural and Sylvatic Habitats in Slovakia

by
Mária Kazimírová
1,*,
Lenka Mahríková
1,
Zuzana Hamšíková
1,
Michal Stanko
1,2,
Maryna Golovchenko
3 and
Natalie Rudenko
3
1
Institute of Zoology, Slovak Academy of Sciences, Dúbravská cesta 9, 84506 Bratislava, Slovakia
2
Institute of Parasitology, Slovak Academy of Sciences, Hlinkova 3, 04001 Košice, Slovakia
3
Biology Centre Czech Academy of Sciences, Institute of Parasitology, Branišovská 31, 37005 České Budějovice, Czech Republic
*
Author to whom correspondence should be addressed.
Microorganisms 2023, 11(7), 1666; https://doi.org/10.3390/microorganisms11071666
Submission received: 17 April 2023 / Revised: 19 June 2023 / Accepted: 25 June 2023 / Published: 27 June 2023
(This article belongs to the Special Issue Advanced Research on Ticks and Tick-Borne Diseases)

Abstract

:
Lyme borreliosis (LB) is the most prevalent tick-borne human infection in Europe, with increasing incidence during the latest decades. Abundant populations of Ixodes ricinus, the main vector of the causative agent, spirochetes from the Borrelia burgdorferi sensu lato (Bbsl) complex, have been observed in urban and suburban areas of Europe, in general, and Slovakia, particularly. Understanding the spread of infectious diseases is crucial for implementing effective control measures. Global changes affect contact rates of humans and animals with Borrelia-infected ticks and increase the risk of contracting LB. The aim of this study was to investigate spatial and temporal variation in prevalence of Bbsl and diversity of its species in questing I. ricinus from three sites representing urban/suburban, natural and agricultural habitat types in Slovakia. Ixodes ricinus nymphs and adults were collected by dragging the vegetation in green areas of Bratislava town (urban/suburban habitat), in the Small Carpathians Mountains (natural habitat) (south-western Slovakia) and in an agricultural habitat at Rozhanovce in eastern Slovakia. Borrelia presence in ticks was detected by PCR and Bbsl species were identified by restriction fragment length polymorphism (RFLP). Borrelia burgdorferi s.l. species in coinfected ticks were identified by reverse line blot. Significant spatial and temporal variability in prevalence of infected ticks was revealed in the explored habitats. The lowest total prevalence was detected in the urban/suburban habitat, whereas higher prevalence was found in the natural and agricultural habitat. Six Bbsl species were detected by RFLP in each habitat type —B. burgdorferi sensu stricto (s.s.), B. afzelii, B. garinii, B. valaisiana, B. lusitaniae and B. spielmanii. Coinfections accounted for 3% of the total infections, whereby B. kurtenbachii was identified by RLB and sequencing in mixed infection with B. burgdorferi s.s, B. garinii and B. valaisiana. This finding represents the first record of B. kurtenbachii in questing I. ricinus in Slovakia and Europe. Variations in the proportion of Bbsl species were found between nymphs and adults, between years and between habitat types. Spatial variations in prevalence patterns and proportion of Bbsl species were also confirmed between locations within a relatively short distance in the urban habitat. Habitat-related and spatial variations in Borrelia prevalence and distribution of Bbsl species are probably associated with the local environmental conditions and vertebrate host spectrum. Due to the presence of Borrelia species pathogenic to humans, all explored sites can be ranked as areas with high epidemiological risk.

1. Introduction

Lyme borreliosis (LB) is a multisystemic disorder with diverse clinical manifestations caused by spirochaetes of the Borrelia burgdorferi sensu lato (s.l.) (Bbsl) complex. In Europe, LB is the most prevalent tick-borne human infection, transmitted by Ixodes ricinus, the main vector of the disease. The Bbsl complex currently includes 23 recognised validated and proposed species, the majority of which occur in the temperate zones of the Northern Hemisphere. Nine species are present in Europe, viz., B. burgdorferi sensu stricto (s.s.), B. garinii, B. afzelii, B. valaisiana, B. lusitaniae, B. bavariensis, B. spielmanii, B. turdi and B. bissettiae [1,2,3], while the validity of B. finlandensis has recently been questioned [3]. Additional species such as B. carolinensis, B. californiensis, B. americana and B. lanei have sporadically been detected in European questing and host feeding Ixodes spp. and in wild animal tissues [4,5,6]. In Europe, pathogenicity to humans has been confirmed for eight species (B. burgdorferi s.s., B. garinii, B. afzelii, B. bavariensis, B. spielmanii, B. bissettiae, B. lusitaniae and B. valaisiana) [7]. Even though associations of a few other species with human cases have been anticipated [3,8], B. afzelii and B. garinii are responsible for the majority of contribution to human LB in Europe. Both species have been found to be the most prevalent among questing I. ricinus ticks in Central Europe and are associated with different clinical manifestations of the disease in humans [1,8,9,10].
The ecology of LB in Europe is very complex. Environmental factors affecting the phenology of the vector tick I. ricinus and vertebrate species abundance and diversity have been assumed to be the main drivers determining the presence and prevalence of Bbsl species in a particular site [3,11,12,13]. However, recently, it has been suggested that only a few widespread vertebrate species (rodents—Rodentia, thrushes—Turdidae, deer—Cervidae) maintain I. ricinus and Bbsl populations in an environment [14,15]. Ixodes ricinus is a generalist tick species with about 300 different vertebrate host species, some of which are competent reservoirs for Bbsl [16]. Such a diversity and quantity of hosts may contribute to tick dispersion and lead to intensive colonization of new areas and establishment of new enzootic LB foci. Studies on reservoir competence have indicated associations between Borrelia species and their vertebrate hosts [17]. Specific associations were identified, e.g., between small rodents such as mice, voles (Murinae, Arvicolinae) and B. afzelii [18]; squirrels (Sciuridae) and B. burgdorferi s.s. [19]; birds (Aves) and B. garinii and B. valaisiana [20,21]; lizards (Lacertidae) and B. lusitaniae [22,23]; or between dormice (Gliridae) and B. spielmanii [24]. However, the natural circulation and transmission patterns of different Bbsl species have still not been fully elucidated. The situation is even more complicated due to existence of multiple diverse strains within the individual Borrelia spp. [25]. Expansion of host specialized Borrelia spp. (e.g., B. afzelii or B. garinii) in different geographic regions results in geographic structuring of the pathogen and depends on the host distribution. In contrast, spirochete species without host specialization (e.g., B. burgdorferi s.s. or B. bissettiae) might infect phylogenetically distant hosts, sometimes covering as many as three orders, generating a spatially uniform population structure [26,27].
Ixodes ricinus has a significant impact on public health, being a major vector of multiple human diseases in Europe. Distribution and phenology of I. ricinus are influenced by various biotic (host species, host density and behaviour, vegetation structure) and abiotic factors (temperature, relative humidity, altitude, latitude), which are the main drivers determining the spatial and temporal changes in prevalence and genetic variability of Bbsl [3,9]. Some biotic and abiotic factors can be effectively managed to decrease the risk of tick bites [28]. Ixodes ricinus is an exophilic tick. It actively quests for hosts and its presence is strongly associated with specific biotopes. Ixodes ricinus mostly occurs in humid habitats, such as deciduous and mixed forests or forest edges with dense undergrowth and permanent leaf litter [29]. During the last decades, spread and increasing abundance of I. ricinus populations have been observed in urban and peri-urban areas of Europe [28,30]. In small urban parks and cemeteries, absence of large mammals serving as feeding hosts for all developmental stages of I. ricinus results in changes in tick–host associations as well as in the prevalence and distribution of some pathogens including Borrelia. Small- to medium-sized mammals, such as mice, voles, shrews (Soricidae), hedgehogs (Erinaceinae), squirrels, dormice, hares (Lagomorpha), foxes, domestic dogs (Canidae) and cats (Felidae), as well as birds (mostly competent reservoirs of Borrelia) become the main hosts for ticks. In contrast, high population densities of roe deer (Capreolus capreolus) in large city parks and peri-urban areas determine the establishment and maintenance of abundant tick populations. Deer (Cervidae) and other large mammals that are incompetent hosts for Borrelia were found to act as “dilution” hosts, which decrease infection prevalence of Bbsl in tick populations; however, their presence can still increase the risk for human populations to acquire LB [15,31].
Similar to the entirety of Europe, increasing incidence of human LB cases and their associations with B. afzelii and B. garinii have been reported in Slovakia during the last decades [32]. Borrelia burgdorferi s.l.-infected questing I. ricinus ticks have been found in different habitat types throughout Slovakia, with considerable spatial variations in prevalence and distribution of Bbsl species [18,20,21,23,33,34,35,36,37,38,39,40,41,42]. Presence of B. burgdorferi s.s., B. garinii, B. afzelii, B. valaisiana, B. lusitaniae, B. bavariensis, B. spielmanii and B. bissettiae has been confirmed [18,20,21,23,35,36,37,38,39,40,41,42,43,44]. Associations were found between B. afzelii, B. garinii and B. lusitaniae with small rodents [18], birds [20,21] and lizards [23], respectively, but more investigations are needed to unravel the factors that determine the differences in prevalence and distribution of Bbsl species and strains in habitats with various diversities of vertebrate hosts.
The main objectives of the present study were (i) assessment of prevalence of Bbsl-infected questing I. ricinus ticks in three habitat types in Slovakia with different anthropogenic impact; (ii) comparison of Bbsl prevalence in the three habitat types between three consecutive years; (iii) identification of Bbsl species and assessment of their diversity and spatial distribution.

2. Materials and Methods

2.1. Study Sites and Tick Collections

The research was carried out in the timeframe of the FP7 project EDENext (www.edenext.eu (accessed on 1 February 2011)), with the aim to investigate prevalence and spread of tick-borne diseases in various habitats in Europe. The methods used were standardised in the frame of the project [45].
Three study sites comprising different habitat types were selected (Figure 1; [45]). Two sites (urban/suburban and natural habitat) were located on the south-eastern slopes of the Small Carpathians Mountains (south-western Slovakia). The urban/suburban habitat (48.166667 N, 17.066667 E) comprised green patches in the north-western part of Bratislava town, including the campus of the Slovak Academy of Sciences (SAS) and Železná studnička (located in the Bratislava forest park recreational area), which are characterized by moderate-to-significant human intervention and high frequency of visitors. The site representing the natural habitat (48.366667 N, 17.300000 E) was located in a non-fragmented deciduous forest at Fúgelka, about 40 km away from Bratislava. The site was visited predominantly by forest workers, gamekeepers, hunters and hikers (for detailed characterisation of the sites, see [46,47]). The third site was located in eastern Slovakia at Rozhanovce village (Košická kotlina basin; 48.750000 N; 21.366667 E), within a game reserve, with a special facility for breeding mouflon (Ovis aries musimon) and fallow deer (Dama dama) along with a pheasantry. Rozhanovce is an agricultural habitat with patches of oak–hornbeam forest, i.e., ecotones between the forest and the cultivated fields and meadows (for details, see [48]), and is long-term monitored for diverse pathogen foci [49].
Questing ticks were collected by dragging of vegetation with a 1 m2 blanket along three 100 m transects set in each habitat. Random additional collections were carried out within close vicinity of the transects. Ticks were collected in monthly intervals during the season of the highest activity of I. ricinus (April–June) and in September–October of 2011–2013 (for details, see [46,47,48]). Ticks were stored in 70% ethanol and identified according to species and developmental stage according to Siuda [50] and Filippova [51]. Only I. ricinus were selected for molecular analyses.
In addition to the transects in Bratislava which were sampled regularly during three consecutive years (A, D, E), two transects (B, C) were selected for random tick collections from April to June 2011 (Figure S1).
The airline distances between the transects are as follows: 350 m between A and B, 420 m between B and C, 2.2 km between C and D and 1.5 km between D and E. The overall airline distance between transects A and E is about 4.5 km. For detailed characterization of the transects, see papers by Svitálková et al. [46] and Mydlová et al. [52].

2.2. DNA Isolation

Genomic DNA was isolated from individual ticks by using the Macherey-Nagel NucleoSpin® Tissue kit (Düren, Germany), following the manufacturer protocol. Quantity and quality of total DNA were measured with a Nanodrop 2000 c spectrophotometer and samples were stored at −20 °C for further analyses [46]. DNA isolation, subsequent PCR amplification and post-amplification analyses were set up in separate areas with all precautions against contamination.

2.3. PCR Amplification and Molecular Analysis

Molecular analyses were performed in all adult ticks and in randomly selected nymphs (whenever possible, in at least 30 individuals per transect and collection date). Total tick DNA was tested for the presence of B. burgdorferi s.l. DNA by PCR amplification of a 222–255-bp fragment of the 5S–23S rRNA intergenic spacer [53]. The primers used were IGSA 50-CGACCTTCTTCGCCTTAAAGC-30 and IGSB 50-AGCTCTTATTCGCTGATGGTA-30. The PCR reaction was carried out in a reaction volume of 25 µL in a BioRad T100TM Thermal Cycler (Hercules, CA, USA) by using the HotStarTaq PCR kit (Qiagen, Hilden, Germany). The PCR reaction mixture per sample contained 2.5 μL of Coral Load PCR buffer (containing 15 mM MgCl2), 1 μL of MgCl2 (25 mM), 0.125 µL of HotStarTaq Plus DNA Polymerase (5 U/µL; Qiagen, Hilden, Germany), 0.5 µL of each primer (10 µM), 0.5 µL of 10 mM dNTP (Thermo Fisher Scientific, Dreieich, Germany) and 14.875 µL of nuclease-free water. Then, 5 µL of total tick DNA was added to the reaction mixture as a template. Negative (nuclease-free water) and positive controls (tick DNA positive for B. burgdorferi s.l.) were included in each run. The PCR was set up according to the protocol by Derdáková et al. [53]. The PCR products were electrophoresed on a 1.5% agarose gel stained with GoodViewTM Nucleic Acid stain (SBS Genetech, Beijing, China) and visualised with a UV transilluminator. In PCR positive samples, borrelial species were identified by restriction fragment length polymorphism (RFLP) analysis [53]. PCR products were digested with FastDigest® Tru1I restriction enzyme (Fermentas, Thermo Fisher Scientific Biosciences, Leon-Rot, Germany), according to the instructions of the manufacturer. Restriction fragments were analysed by electrophoresis using precast gels (Spreadex gel EL300 mini, Elchrom Scientific, Cham, Switzerland) in the Origins submarine electrophoresis system (Elchrom Scientific, Cham, Switzerland). Gels were stained with SYBRGreen (Sigma-Aldrich, St. Louis, MO, USA) and visualised by UV transilluminator.

2.4. Identification of Borrelia Species in Samples with Mixed Infections

The presence of multiple species of LB spirochetes in tick samples was detected by initial PCR amplification. To identify the spirochete species, amplification of the partial genes encoding flagellin and outer surface protein C (ospC) was conducted, followed by sequencing of amplicons. To avoid the inhibition of the reaction by the excess of tick DNA, all PCR reactions were conducted in two steps—the first (spacer) and the second (nested) using the protocol previously described [54]. As a template for the second step, 5 µL of the resulting first-step reaction mixture was used. In case of overlap of sequence chromatograms that confirmed the presence of multiple infections in a tick, the cloning step of total PCR product into the TOPO TA cloning vector (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) was conducted as described below.
The fragment of ospC gene was amplified by spacer/nested PCR using the protocol of Bunikis et al. [55]. The first round of PCR was conducted using the primers F1, 5′-ATGAAAAAGAATACATTAAGTGC-3′ (positions 306–328 of U01894) and R1, 5′-ATTAATCTTATAATATTGATTTTAATTAAGG-3′ (963–933), targeting 657 bp-long amplicon. The second round of PCR was conducted using the primers F2, 5′-TATTAATGACTTTATTTTTATTTATATCT-3′ (331–359) and R2, 5′-TTGATTTTAATTAAGGTTTTTTTGG-3′ (948–924), targeting 617 bp-long fragment. The conditions of both reactions were the same: 30 cycles at 95 °C for 30 s, 52 °C for 30 s and 72 °C for 30 s, except for the annealing temperature, which was 50 and 52 °C for the first and the second rounds of PCR, respectively.
Primers described by Clark et al. [56] were used to amplify the fragment of flagellin gene—F1: 5′- AARGAATTGGCAGTTCAATC-3′ (271–290) and R1: 5′- GCATTTTCWATTTTAGCAAGTGATG-3′ (743–767), targeting 465 bp-long fragment in the first round of PCR, and F2: 5′- ACATATTCAGATGCAGACAGAGGTTCTA -3′ (301–328) and R2: 5′- GAAGGTGCTGTAGCAGGTGCTGGCTGT-3′ (663–689), targeting 388 bp-long fragment in the second round of PCR. The conditions of PCR reaction were exactly the same as in the case with ospC gene PCR amplification.
All PCR reactions were carried out in a final volume of 20 μL using 2x HotStarTaq Plus Master Mix (Qiagen, Hilden, Germany). Amplicons were visualized by electrophoresis in a 1% agarose gel (1×  TAE, pH 8.0). SYBR® Gold DNA stain (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) was included in the sample loading buffer before the loading of samples into the gel. In all cases, a reaction mixture with water instead of DNA template was used as negative control. Purified DNA from Bbsl was used as a positive control in all PCR reactions.

2.4.1. Sequencing

All PCR products of the expected sizes were excised from agarose gels, purified using Ultrafree-DA centrifugal filter units for DNA extraction from agarose (Millipore, Darmstadt, Germany) and sequenced in both directions using the same primers as for PCR. Sequence analysis was performed by SeqMe, (SEQme s.r.o., Dobříš, Czech Republic) and the sequences were compared to those available in the GenBank™ database by Basic Local Alignment Tool (BLAST, https://blast.ncbi.nlm.nih.gov/) (accessed on 1 April 2015) analysis.

2.4.2. Reverse Line Blotting (RLB)

DNA samples were subjected to a conventional PCR, targeting a fragment of the 5S–23S rRNA intergenic spacer (IGS) by use of biotin-linked forward B5S-Bor (biotin-5′-GAGTTCGCGGGAGAGTAGGTTATT-3′), 23SBor reverse (5′-TCAGGGTACTTAGATG GTTCACTT-3′) and biotin-linked forward BMiya-For (biotin-5′-TTAGGATTAATGATRT TKTTACC-3′) primers, as described by Alekseev et al. [57] and Blazejak et al. [58]. PCR cycling conditions included 45 cycles at 94 °C for 20 s, 52 °C for 30 s, 72 °C for 30 s and a final elongation step at 72 °C for 7 min. Each run included different Borrelia strains corresponding to the used probes as positive controls (Supplementary file): B. burgdorferi s.l. (-), B. burgdorferi s.s. (PAbe), B. garinii (PWudII), B. afzelii (PBas), B. valaisiana (VS116), B. lusitaniae (Poti B2), B. spielmanii (PHap), B. bissettiae (DN127), RF-like Borrelia (-) and B. miyamotoi (HT31). PCR products, as well as positive controls and a negative control, were hybridized to different oligonucleotide probes linked to the membrane (for details, see Supplementary file and [58,59,60]).

2.4.3. Cloning and Transformation of the Competent Cells

For separation of multiple species present as coinfection in one sample, total PCR products were cut from the gel, purified and cloned into pCR-4 TOPO plasmid, followed by the transformation of the competent cells and plating of them on LB/agar medium. TOPO™ TA Cloning™ Kit for Sequencing (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) was used. The detailed protocol is given in the Supplementary file. Next, 12 colonies were picked individually from the plates and grown overnight in 10 mL LB/amp liquid medium. The overnight culture was centrifuged and the collected cells were used for plasmid DNA isolation using the QIApred Miniprep Kit for plasmid DNA purification (Qiagen, Hilden, Germany). The plasmid DNAs were checked for the presence of the insert required using the PCR with the specific primers. The plasmids with the confirmed insert were sent for sequencing.

2.4.4. RFLP Analysis

RFLP analysis of Borrelia sequences was performed in silico using the free software available at http://insilico.ehu.es (accessed on 15 June 2015). The partial sequence of fla gene was digested with HapII, HhaI, HincII, CelII and DdeI. All RFLP patterns obtained were compared with patterns already published [61,62,63].

2.5. Statistical Analysis

Differences in the prevalence of Bbsl-infected ticks between habitats, years, developmental stages and sexes, as well as differences in the proportion of Borrelia species between habitats, were analysed by Chi-square test or Fisher’s exact test using the PAST programme [64]. Additionally, 95% confidence intervals (95% CI) were calculated individually for each proportion by using the online calculator at the website http://epitools.ausvet.com.au (accessed on 15 February 2023). Association between percentage proportion of nymphs in the local tick populations in Bratislava (transects A–E) and total Bbsl prevalence was analysed by Spearman rank correlation using the PAST programme. p values < 0.05 were considered significant in all tests.

3. Results

3.1. Prevalence of Borrelia burgdorferi s.l. in Questing Ixodes ricinus

In total, 9071 I. ricinus ticks (6954 nymphs, 976 females, 1141 males) were collected during 2011–2013. Of these, 4587 ticks (3305 nymphs, 590 females, 692 males) were collected in Bratislava, 2939 (2307 nymphs, 286 females, 346 males) in Fúgelka and 1545 (1342 nymphs, 100 females, 103 males) in Rozhanovce. A total of 6123 ticks (4006 nymphs, 976 females and 1141 males; 2816 ticks from Bratislava, 2251 from Fúgelka and 1056 from Rozhanovce) were examined, resulting in an overall Bbsl infection prevalence of 12.5% (765/6123; 95% CI: 11.7–13.3%). The total prevalence differed significantly among the three habitats (Table 1; Figure 2): Bratislava 9.4%; Fúgelka 15.4%; Rozhanovce 14.5% (χ2 = 44.813; p < 0.0001). Prevalence in Bratislava (urban/suburban habitat) was significantly lower than in Fúgelka (natural habitat) and Rozhanovce (agricultural habitat) (χ2 = 41.355; p < 0.0001 and χ2 = 20.236; p < 0.0001, respectively), while it did not differ between Fúgelka and Rozhanovce (χ2 = 0.437; p = 0.509). Similar to overall prevalence, statistical analysis revealed significant differences in the prevalence of Borrelia-infected nymphs among the three habitats (χ2 = 83.474; p < 0.0001) (Figure 2), between Bratislava and Fúgelka (χ2 = 82.781; p < 0.0001) and Bratislava and Rozhanovce (χ2 = 40.882; p < 0.0001). The difference between Fúgelka and Rozhanovce was not significant (χ2 = 2.745; p = 0.097). For adult ticks, infection prevalence did not differ among the habitats (χ2 = 2.923; p = 0.232) (Figure 2). Significant differences among habitats were also found for infection prevalence in nymphs in individual years (2011: χ2 = 46.474; p < 0.0001; 2012: χ2 = 8.022; p = 0.018; 2013: χ2 = 31.628; p < 0.0001), but not for adult ticks.
The total prevalence of infection in adult ticks (13.8%) was significantly higher than in nymphs (11.8%) (χ2 = 5.366; p = 0.020). With respect to habitat type, the difference between adults and nymphs was significant in Bratislava (χ2 = 46.850; p < 0.0001) but not in Fúgelka (χ2 = 3.383; p = 0.066) nor Rozhanovce (χ2 = 2.136; p = 0.144) (Figure 3).
Significantly higher total prevalence was determined in female I. ricinus than in males (χ2 = 29.364; p < 0.0001). The differences between sexes were also significant in Bratislava (χ2 = 17.987; p < 0.0001) and Fúgelka (χ2 = 13.346; p < 0.001), while in Rozhanovce the difference was not significant (χ2 = 0.216; p = 0.642) (Table 1).
By comparing the prevalence of infection in I. ricinus among the three years (2011–2013) separately in each habitat, significant differences were revealed for nymphs and for nymphs + adult ticks in both Bratislava and Fúgelka, as well as for adult ticks in Rozhanovce (Table 1).

3.2. Occurrence of Borrelia burgdorferi s.l. Species

Overall, six Bbsl species were detected by RFLP: B. burgdorferi s.s., B. afzelii, B. garinii (including B. bavariensis), B. valaisiana, B. lusitaniae and B. spielmanii. Coinfections with two or more species were detected in 23 (3%) out of the 765 Borrelia-infected ticks. The presence of the six species was confirmed in each habitat, in nymphs, in adults or both developmental stages and at least in one year (Table 2, Figure 4). Borrelia afzelii, B. garinii and B. valaisiana were present in all habitats and during all analysed years. Details on the numbers of ticks (nymphs, females, males) infected with the different Borrelia species per habitat type and year are presented in Table S1.
Differences in the proportion of ticks infected with individual Bbsl species among habitats were significant for B. burgdorferi s.s. (χ2 = 15.100; p = 0.0005), B. afzelii2 = 9.051; p = 0.011), B. garinii2 = 53.501; p < 0.0001), B. valaisiana2 = 14.803; p < 0.001) and B. spielmanii (Fisher’s exact test; p = 0.022), while they were not significant for B. lusitaniae2 = 1.162; p = 0.559). Paired comparisons between habitats revealed significant differences in the proportion of B. burgdorferi s.s. between Bratislava and Rozhanovce (χ2 = 13.167; p < 0.001) and between Fúgelka and Rozhanovce (χ2 = 7.752; p = 0.005), in the proportion of B. afzelii between Bratislava and Fúgelka (χ2 = 8.343; p = 0.004) and between Bratislava and Rozhanovce (χ2 = 3.959; p = 0.047), in the proportion of B. garinii between all habitats (Bratislava vs. Fúgelka χ2 = 53.223; p < 0.0001; Bratislava vs. Rozhanovce χ2 = 8.041; p = 0.004; Fúgelka vs. Rozhanovce χ2 = 7.911; p = 0.005), in the proportion of B. valaisiana between Bratislava and Rozhanovce (χ2 = 15.031; p < 0.001) and between Fúgelka and Rozhanovce (χ2 = 4.004; p = 0.045) and in the proportion of B. spielmanii between Bratislava and Fúgelka (Fisher’s exact test; p = 0.017). The proportion of ticks with coinfections differed significantly among the three habitats (Fisher’s exact test; p = 0.0001), whereas paired comparisons revealed significant differences between Bratislava and Fúgelka (Fisher’s exact test; p = 0.008) and between Bratislava and Rozhanovce (Fisher’s exact test; p < 0.0001).

3.3. Variability in Infection Prevalence and Distribution of Borrelia burgdorferi s.l. Species in Bratislava

The percentage share of nymphs in the total number of collected ticks varied in the studied transects between 6.7% (C) and 81.2% (A). Out of the 1694 examined ticks from this locality, the presence of Bbsl was confirmed in 192 individuals (11.3%; 95% CI: 9.8–12.8%). Significant differences among the transects were found for prevalence of Borrelia-infected nymphs and adults as well as for total prevalence (Table 3). A significantly negative association between the percentage proportion of nymphs in the tick populations in transects A–E and total Borrelia prevalence was observed (rs = −0.9, p = 0.037).
The proportion of the identified Borrelia species varied considerably between transects, even between those within a relatively short distance (Figure S2). Four species, B. afzelii, B. garinii, B. valaisiana and B. lusitaniae, were found in each of the five transects in nymphs, adults or both. Borrelia burgdorferi s.s. was present only in nymphs collected in transect B and B. spielmanii in nymphs in transect B, as well as adults in transects A, B and D. Prevalence of B. afzelii was the highest (37.5–67.4%) in all transects, except for transect E, where B. lusitaniae prevailed (50%). Borrelia garinii followed B. afzelii in transects A, C and D (31.2%, 25.6% and 42.9%, respectively). The highest proportion of B. valaisiana (18.5%) was determined in transect A, where, together with B. garinii, they accounted for 49.5% of infected ticks. Coinfections were detected in two ticks (B. garinii + B. valaisiana in a female from transect D, not identified in a female from transect C).

3.4. Borrelia burgdorferi s.l. Species in Mixed Infections

The presence of multiple spirochete species in ticks that feed few times during their lifespan, mostly on different host species, is not surprising. Identification of species in coinfected samples was performed by RLB, followed by a cloning step and sequencing. In general, by using these combined methods to separate multiple species detected in a single sample, the presence of nine spirochete species was confirmed: B. burgdorferi s.s., B. garinii, B. afzelii, B. valaisiana, B. bissettiae, B. spielmanii, B. lusitaniae, B. kurtenbachii and B. bavariensis. The majority of coinfections (52.2%) in the present study represented samples with paired host associated species, such as bird-associated B. garinii + B. valaisiana that occurred in each habitat type. Samples where bird-associated species were accompanied with rodent-associated B. afzelii were detected as well, though in a lower ratio (8.7%) (for detailed data, see Table S1). Considering the proportion of B. garinii + B. valaisiana, the difference among the three habitats was significant (Fisher’s exact test; p = 0.003), whereas paired comparisons revealed a significant difference only between Bratislava and Rozhanovce (Fisher’s exact test; p = 0.002). In 30.4% of tick DNA samples for which the RFLP profile indicated mixed infections, only B. garinii could be identified by RLB and sequencing, and no DNA was available for further analyses. It is worth mentioning that, in one nymph from Fúgelka (natural habitat), the presence of B. garinii + B. burgdorferi s.s. was detected and, in another nymph, three species (B. garinii, B. valaisiana, B. afzelii) were identified. In one single female from Rozhanovce (agricultural habitat), the presence of four species of LB spirochetes was detected: B. burgdorferi s.s., B. garinii, B. valaisiana and B. kurtenbachii. The presence of B. kurtenbachii was confirmed by the results of RLB and the sequence, which was identical to the type strain 25015. To the best of our knowledge, this is the first record of B. kurtenbachii in a natural I. ricinus population from Slovakia and Europe.

4. Discussion

Spatial and temporal variations in prevalence and genetic diversity of Bbsl are phenomena reported from a number of sites across Europe which differ in geographic position, habitat structure, biodiversity of vertebrate species and land-use patterns [65,66,67,68,69,70,71,72,73,74,75,76]. Due to the current distribution, diversity and invasiveness of Borrelia spp. in Europe [3,10] and an increasing overlap of human and spirochaete habitats [2,6,30], the risk to contract LB is higher than has been reported. Climate change, affecting abundance and distribution of ticks and vertebrates, prevalence of Borrelia infection in ticks, land use patterns and human activities, may have unexpectedly significant impact on risk prediction of LB [26].
In the present study, we investigated the prevalence of infection and distribution of Bbsl species in questing I. ricinus ticks in three habitat types in Slovakia during three consecutive years. We hypothesised that the proportion of infected ticks and diversity of Bbsl species differ between urban/suburban green areas, un-fragmented deciduous forests and agricultural areas due to variations in the structure of the habitats and differences in the diversity of vertebrate species serving as tick feeding hosts and reservoirs for the pathogens. Our results are largely consistent with previous observations. The overall prevalence of Bbsl in our study was lower or in agreement with the 13.7% and 12.3% mean prevalence reported in Europe-wide meta-analyses based on prevalence data in questing ticks from 1984–2003 [77] and 2010–2016 [10], respectively, but it was considerably lower than the mean prevalence (~20%) calculated for Central Europe [10,77]. The mean prevalence of infected nymphs in our study was higher or the same as the values calculated in the cited meta-analyses (10.1% and 11.8%, respectively), whereas, for adults, the prevalence was lower than the values 18.6% and 14.9% obtained by Rauter and Hartung [77] and Strnad et al. [10], respectively.
Generally, Bbsl prevalence in nymphs was found to be lower than in adults [10,77], which has been attributed to a higher probability to acquire the pathogen for adults than for nymphs due to transstadial transmission and an additional blood-meal from reservoir hosts. Results of our study follow this trend in the urban and agricultural habitat but, in the natural habitat, prevalence in nymphs was higher than in adults. This could be explained, in part, by the fact that, in the natural habitat, the diversity of feeding hosts which are also competent reservoirs is higher for larvae than in the other habitat types, while nymphs may feed more frequently on species (e.g., deer) that are not competent reservoirs for Bbsl. On the other hand, prevalence in adult ticks from the natural and urban habitats was almost the same, which is probably due to the presence of roe deer in some urban locations and connectivity of the other locations with areas where roe deer or other cervids are present. In the agricultural site where ticks were collected inside the enclosure with game (fallow deer, mouflon, pheasant Phasianus colchicus) and in the ecotone between fields and oak–hornbeam forest, the Bbsl prevalence in adult ticks was higher than in the other explored habitat types, probably due to presence of abundant populations of competent hosts such as rodents and ground-feeding birds [49]. In previous studies, circulation of Bbsl in Rozhanovce was confirmed based on presence of antibodies against Borrelia in 10.7% of mouflons, 45.8% of fallow deer [78] and in 9–40% of small mammals [49,79].
In a study by Rosá et al. [45], aimed at modelling density of questing I. ricinus nymphs (DON) and their infection rate with Bbsl, Anaplasma phagocytophilum and Rickettsia spp. in urban, agricultural and natural habitats of Italy, Germany, Czech Republic, Hungary and Slovakia by considering various climatic and environmental factors, the same dataset of nymphs (excluding randomly sampled ticks) from three 100 m-long transects per habitat and period from April to June of 2011–2013 was used, as that analysed in the present study. Interestingly, the highest DON was found in Slovakia (average 518.7 nymphs/year/300 m2). By summarising data from all countries, the highest DON was found in the natural habitat, followed by the urban and agricultural habitats, while in Slovakia, the highest DON was found in the urban habitat. The overall prevalence of Bbsl in nymphs based on data from Slovakia, Germany and Italy was 19.3% (i.e., a higher value than in the present study). The mentioned findings also suggest that, in the explored urban/suburban habitat of Slovakia, due to presence of tick feeding hosts that are not competent reservoirs of Borrelia (e.g., deer), the DON is high but the prevalence of Bbsl-infected nymphs is relatively low.
The proportion of Bbsl-infected ticks in previous studies from Slovakia was found to vary, depending on the location and habitat type, with values ranging from 1.0% to 27.8% in nymphs and from 0% to 49.4% in adults, and with a total prevalence ranging between 0% and 38.3% in urban/suburban sites, between 8.0% and 29.4% in natural sites and between 32.3% and 38.1% in agricultural sites and ecotones (Table S2, [80,81,82,83,84]). Thus, our prevalence values for the three studied habitat types are in the lower range limit of the prevalences reported from Slovakia.
Although the abovementioned meta-analyses [10,77] did not show any significant effect of different PCR-based detection methods on prevalence rates, we came to a different conclusion. In a previous study aimed at determination of mixed infections of Bbsl and B. miyamotoi, 136 randomly selected ticks from the same sample set as that analysed in the current study were screened for presence of Bbsl by qPCR targeting fragments of the ospA and flaB genes, and species were identified by sequencing [44]. A 21.3% prevalence of infection was determined and, in 13.2% of the samples, Bbsl species were identified, while, in the present study, the prevalence rate in the same selected ticks was 9.5%, but species could be identified in all positive samples.
Borrelia afzelii and B. garinii are the most frequent species in Europe, followed by B. burgdorferi s.s., B. valaisiana and B. lusitaniae; however, their representation was found to vary between regions and countries [9,10,77]. For Central Europe, the mean percentual proportions were as follows: B. afzelii (44%), B. garinii (23%), B. burgdorferi s.s. (15%), B. valaisiana (11%), B. lusitaniae (6%), B. spielmanii (2%), B. bissettiae (0.1%) and B. bavariensis (2%) [10]. The listed eight species were also reported from Slovakia, with B. afzelii prevailing in the majority of the explored sites (Table S2). Considering total prevalence of Bbsl species, our results partly agree with previous findings. Borrelia afzelii and B. garinii were most prevalent, while B. valaisiana and B. lusitaniae were more prevalent than B. burgdorferi s.s. and the proportion of B. spielmanii was 2.1%. By the PCR and RFLP method used in our study, identification of B. bavariensis was not possible (see above); thus, infections with B. garinii may include B. bavariensis. The assumption that B. bavariensis circulates in the three studied habitats is supported by its previous detection in a yellow-necked mouse (Apodemus flavicollis) specimen from Bratislava forest park and in a questing I. ricinus female from Fúgelka [44], as well as in I. ricinus nymphs (formerly B. garinii ospA serotype 4) from Rozhanovce [37].
Established foci of Bbsl spirochetes should possess an appropriate environment for vector survival, an adequate density of reservoirs and a temporal and spatial overlap between them. The same six Bbsl species as mentioned above were detected in the three studied habitats, suggesting that, in all of them, suitable conditions for their circulation exist. In the present study, the proportion patterns of the Bbsl species varied according to habitat and tick developmental stage, as well as between years. The proportion of B. afzelii was mostly higher in nymphs than in adult ticks, probably because larvae fed preferably on rodents while nymphs fed on other hosts. In the natural habitat, bird-associated species (B. garinii, B. valaisiana) prevailed. This may suggest that a higher portion of ticks fed on birds than in the urban/suburban habitat, although rodents were abundant there in 2012, while they declined in 2013 after a long and cold winter [46]. The marked temporal changes in the Bbsl species distribution patterns, with a decline in the proportion of B. afzelii in the agricultural habitat, are difficult to explain. Fluctuations in rodent populations, phenology of I. ricinus and presence of a pheasantry could partly account for the changing patterns in favour of bird-associated species.
Coinfections in tick vectors and reservoir hosts with multiple Borrelia species are common in Europe. The most prevalent (51%) are coinfections with the bird-associated species B. garinii + B. valaisiana, followed by B. afzelii and B. burgdorferi s.s. [77]; however, triple and quadruple infections with species associated with different hosts have also been detected in single ticks (e.g., [68,70,85]). Coinfections with B. garinii and B. valaisiana were found to prevail also in Slovakia, evidenced in [20,21,42] as well as in the present study. Recently, in neighbouring countries to Slovakia, the presence of “atypical” spirochete species, or their combinations, were detected in questing ticks in Poland [4] or in vertebrate hosts in the Czech Republic [6]. For the first time, the coinfection with B. bissettiae and B. carolinensis was shown in cadavers of red squirrels (Sciurus vulgaris) and European hedgehogs (Erinaceus europaeus). Using the multiple analyses of several tissues from cadavers, mixed infections of two or three species were shown rather often in northern white-breasted hedgehogs (Erinaceus roumanicus, including B. burgdorferi s.s. + B. afzelii + B. garinii in one animal) [6]. Previously considered as “rodent associated”, B. bavariensis and B. afzelii, were detected in tissues of common blackbirds (Turdus merula), together with the expected B. garinii, B. valaisiana and B. turdi [6]. In a single questing I. ricinus nymph from Slovakia, mixed infection with B. afzelii + B. garinii + B. bissettii was confirmed [20]. Detection of unexpected B. kurtenbachii as one of the species in a coinfected sample in the present study may be in line with our speculation about the role of migrating hosts in distribution of the causative agents of LB. Another explanation may be that B. kurtenbachii, formerly known as B. bissettii strain 25015 from human samples from Europe [86,87], has circulated in natural I. ricinus populations for a long time, but has remained unrecognized due to its low prevalence and insufficient detection methods.
Due to increasing exposure of humans and domestic animals to tick bites in urban parks, gardens, cemeteries and peri-urban recreational areas, especially during the last decades, these habitats are of special interest from an epidemiological point of view [30,88]. In a number of studies, the presence and prevalence of B. burgdorferi s.l. in I. ricinus ticks from urban and peri-urban areas in Europe have been reported. The findings from the period 1990–2021 were reviewed by Hansford et al. [88], who calculated a 17.3% mean prevalence of Borrelia infection—14.2% in nymphs and 21.1% in adult ticks. In our study, the total prevalence, as well as the prevalence in nymphs and adults in the urban/sub-urban habitat, were lower than the European average. Total prevalence in our study was also in the lower limit of previously published values from green areas of Bratislava and other urban areas in Slovakia [38,40,41] (Table S2).
Marked spatial variations in total prevalence of infected ticks have been reported from a number of European cities. For example, in Budapest, prevalence ranged from 31.3% in a cemetery to 50% in an urban forest [89] and values ranging from 1.8% to 60.5% were determined in 11 green urban locations in Prague [90]. In 10 different recreational areas in Hanover, Borrelia infection rates in I. ricinus varied between 14.8% and 40.5% [58,91]. In two green areas in Munich, 20.1% and 30.1% of ticks were infected [70], while, in another study, seven green plots were explored and the Bbsl prevalence varied from 17.4 to 59.5% in nymphs and from 0 to 61.5% in adults [92]. In two cities in Poland, Olsztyn and Warsaw, the proportion of infected ticks ranged between 18.2 and 25% [71] and 7.7 and 14% [93], respectively. Great spatial variability in the spectrum of Bbsl species and in their proportions, mostly with predominance of B. afzelii or B. garinii, were also reported for the mentioned cities.
Spatial variability in Borrelia prevalence has been reported also in Slovak cities. For example, 1.8–18% and 0–33.3% total prevalences were determined in different green locations of Košice and Bardejov (eastern Slovakia), respectively [38], while, in green areas of Bratislava, prevalence of infected ticks varied between 6.8 and 30.9% [40,41]. Thus, spatial differences in Borrelia prevalence and distribution of Bbsl species within the same cities are common and are likely affected by the local environmental conditions, microclimate, structure of the biotope, biodiversity of vertebrate hosts, connectivity of the locations to forested areas, etc. However, in some cases, the results may be biased due to low numbers of examined ticks. In general, our results on the spatial variability in Bbsl prevalence between the five locations (A–E) in Bratislava (see Figure S2) are in line with previous findings on the mosaic patterns of Bbsl infection rates in ticks from different patches of the same area [94,95]. The situation in the SAS campus where two sampling transects (A, B) were located at a distance of about 350 m deserves particular attention. The campus is a fenced area with low abundance of small rodents and a local roe deer population that obviously supports the high abundance of I. ricinus [46], which may act as a dilution factor for Borrelia. Significantly lower Bbls prevalence was found in transect A than B (see Table 3), which may suggest that roe deer concentrate mainly around transect A, while they avoid frequent encounters with people in the part of the campus where transect B is located. A third transect within the SAS campus with similar environmental conditions as transect A was sampled in the same period (2011–2012) by Chvostáč et al. [40], who determined a 6.8% prevalence by the same PCR method as applied in our study. The three locations in the SAS campus differed also in the proportion of Bbsl species. In the transect studied by Chvostáč et al. [40], bird-associated species predominated, similarly as in transect A, while, in transect B, prevalence of B. afzelii was higher than of B. garinii + B. valaisiana. Thus, even at low abundance of rodents, more than 30% of ticks carried B. afzelii. Overall, the findings in the SAS campus indicate that birds, rodents and other small- and medium-sized mammals (e.g., hedgehogs) are the important reservoir hosts for Borrelia, while roe deer maintain an abundant I. ricinus population. In contrast to the SAS campus, rodents were more abundant in transects D and E, located in Bratislava forest park [46], which was displayed in the high proportion of nymphs infected with B. afzelii in transect D. In contrast to transect D, considerably lower prevalence was determined in ticks from transect E (Table 3) (distance between the transects was 1.5 km). Interestingly, 50% of both nymphs and adult ticks collected in transect E were infected with B. lusitaniae, suggesting the presence of an abundant lizard population. In comparison with our results, in another location in the forest park at a distance of about 2 km from transect D, a total prevalence of 12.6% with predominance of B. afzelii (comparable with transect D) was found [41].

5. Conclusions

Presence of questing I. ricinus ticks infected with Bbsl was confirmed during a three-year study in three habitat types in Slovakia. Prevalence of infected ticks was found to differ between the habitats, locations, years and developmental stages (nymphs and adults). Presence of the same six Borrelia species was identified in each of the habitats with variable proportions, depending on the habitat type, year and tick developmental stage. The study confirmed that urban/suburban sites and parks can be focal points for circulation of LB spirochaetes. Due to the presence of human-pathogenic Borrelia species such as B. garinii, B. afzelii, B. burgdorferi s.s. and B. spielmanii, the explored habitats can be ranked among areas with high epidemiological risk. More research is needed to explain circulation patterns of individual Borrelia species in habitats with different diversity of reservoirs, as well as on the effect of migrating vertebrate hosts on Borrelia prevalence and establishment of new foci of pathogens with impact on public health in highly populated regions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11071666/s1, Figure S1: Localisation of transects in Bratislava (urban/suburban habitat) for tick collection in 2011; Supplementary file: Probes used for RLB and the Borrelia strains used as positive control, RLB experimental procedure, Cloning and transformation of the competent cells; Table S1: Numbers of Ixodes ricinus nymphs, females and males infected with species of Borrelia burgdorferi s.l. in three habitat types between 2011 and 2013; Figure S2: Proportion of Borrelia burgdorferi s.l. species (in %) in infected Ixodes ricinus per transects A–E in Bratislava (urban/suburban habitat) in 2011. Species: ss—sensu stricto, l—lusitaniae, a—afzelii, g—garinii (+ bavariensis), v—valaisiana, s—spielmanii, coinf—coinfections; Table S2: Prevalence of Borrelia burgdorferi s.l. and diversity of Borreia species reported from Slovakia.

Author Contributions

Conceptualisation, M.K.; methodology, L.M., Z.H., M.G. and N.R.; investigation, M.K., M.S. and N.R.; resources, M.K.; data curation: M.K., L.M. and Z.H.; writing—original draft preparation, M.K. and N.R.; writing—review and editing, M.K., N.R. and M.S.; visualisation, M.K.; supervision, M.K., M.S. and N.R.; project administration, M.K.; funding acquisition, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by EU grant FP7-261504 EDENext (http://www.edenext.eu). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission.

Data Availability Statement

All data generated and analysed during this study are included in the published article and its supplementary information files.

Acknowledgments

Partial analysis was carried out within the framework of the COST Action CA21170 “Prevention, anticipation and mitigation of tick-borne disease risk applying the DAMA protocol (PRAGMATICK)”. We thank Markéta Derdáková, Veronika Rusňáková Taragel’ová and Diana Selyemová for their help in methodology and Elena Kocianová, Mirko Slovák, Lenka Minichová, Martin Bona, Ladislav Mošanský, Jasna Kraljik and Lucia Blaňarová for help in fieldwork.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Rizzoli, A.; Hauffe, H.C.; Carpi, G.; Vourc’h, G.I.; Neteler, M.; Rosà, R. Lyme borreliosis in Europe. Euro Surveill. 2011, 16, 19906. [Google Scholar] [CrossRef] [PubMed]
  2. Rudenko, N.; Golovchenko, M.; Grubhoffer, L.; Oliver, J.H., Jr. Updates on Borrelia burgdorferi sensu lato complex with respect to public health. Ticks Tick Borne Dis. 2011, 2, 123–128. [Google Scholar] [CrossRef] [Green Version]
  3. Steinbrink, A.; Brugger, K.; Margos, G.; Kraiczy, P.; Klimpel, S. The evolving story of Borrelia burgdorferi sensu lato transmission in Europe. Parasitol. Res. 2022, 121, 781–803. [Google Scholar] [CrossRef] [PubMed]
  4. Wodecka, B.; Kolomiiets, V. Genetic diversity of Borreliaceae species detected in natural populations of Ixodes ricinus ticks in Northern Poland. Life 2023, 13, 972. [Google Scholar] [CrossRef] [PubMed]
  5. Wodecka, B.; Michalik, J.; Grochowalska, R. Red foxes (Vulpes vulpes) are exposed to high diversity of Borrelia burgdorferi sensu lato species infecting fox-derived Ixodes ticks in west-central Poland. Pathogens 2022, 11, 696. [Google Scholar] [CrossRef] [PubMed]
  6. Majerová, K.; Hönig, V.; Houda, M.; Papežík, P.; Fonville, M.; Sprong, H.; Rudenko, N.; Golovchenko, M.; Černá Bolfíková, B.; Hulva, P.; et al. Hedgehogs, squirrels, and blackbirds as sentinel hosts for active surveillance of Borrelia miyamotoi and Borrelia burgdorferi complex in urban and rural environments. Microorganisms 2020, 8, 1908. [Google Scholar] [CrossRef]
  7. Rudenko, N.; Golovchenko, M. Sexual transmission of Lyme borreliosis? The question that calls for an answer. Trop. Med. Infect. Dis. 2021, 6, 87. [Google Scholar] [CrossRef]
  8. Stanek, G.; Reiter, M. The expanding Lyme Borrelia complex-clinical significance of genomic species? Clin. Microbiol. Infect. 2011, 17, 487–493. [Google Scholar] [CrossRef] [Green Version]
  9. Estrada-Peña, A.; Ortega, C.; Sánchez, N.; Desimone, L.; Sudre, B.; Suk, J.E.; Semenza, J.C. Correlation of Borrelia burgdorferi sensu lato prevalence in questing Ixodes ricinus ticks with specific abiotic traits in the western Palearctic. Appl. Environ. Microbiol. 2011, 77, 3838–3845. [Google Scholar] [CrossRef] [Green Version]
  10. Strnad, M.; Honig, V.; Růžek, D.; Grubhoffer, L.; Rego, R.O.M. Europe-wide meta-analysis of Borrelia burgdorferi sensu lato prevalence in questing Ixodes ricinus ticks. Appl. Environ. Microbiol. 2017, 83, e00609-17. [Google Scholar] [CrossRef] [Green Version]
  11. Kurtenbach, K.; De Michelis, S.; Etti, S.; Schäfer, S.M.; Sewell, H.S.; Brade, V.; Kraiczy, P. Host association of Borrelia burgdorferi sensu lato–the key role of host complement. Trends Microbiol. 2002, 10, 74–79. [Google Scholar] [CrossRef] [PubMed]
  12. Kurtenbach, K.; Hanincová, K.; Tsao, J.I.; Margos, G.; Fish, D.; Ogden, N.H. Fundamental processes in the evolutionary ecology of Lyme borreliosis. Nat. Rev. Microbiol. 2006, 4, 660–669. [Google Scholar] [CrossRef]
  13. O’Keeffe, K.R.; Oppler, Z.J.; Brisson, D. Evolutionary ecology of Lyme Borrelia. Infect. Genet. Evol. 2020, 85, 104570. [Google Scholar] [CrossRef] [PubMed]
  14. Hofmeester, T.R.; Coipan, E.C.; van Wieren, S.E.; Prins, H.H.T.; Takken, W.; Sprong, H. Few vertebrate species dominate the Borrelia burgdorferi s.l. life cycle. Environ. Res. Lett. 2016, 11, 043001. [Google Scholar] [CrossRef] [Green Version]
  15. Heylen, D.; Lasters, R.; Adriaensen, F.; Fonville, M.; Sprong, H.; Matthysen, E. Ticks and tick-borne diseases in the city: Role of landscape connectivity and green space characteristics in a metropolitan area. Sci. Total Environ. 2019, 670, 941–949. [Google Scholar] [CrossRef]
  16. Gern, L. Borrelia burgdorferi sensu lato, the agent of Lyme borreliosis: Life in the wilds. Parasite 2008, 15, 244–247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Humair, P.F.; Gern, L. The wild hidden face of Lyme borreliosis in Europe. Microbes Infect. 2000, 2, 915–922. [Google Scholar] [CrossRef]
  18. Hanincová, K.; Schäfer, S.M.; Etti, S.; Sewell, H.S.; Tarageľová, V.; Žiak, D.; Labuda, M.; Kurtenbach, K. Association of Borrelia afzelii with rodents in Europe. Parasitology 2003, 126, 11–20. [Google Scholar] [CrossRef]
  19. Humair, P.F.; Gern, L. Relationship between Borrelia burgdorferi sensu lato species, red squirrels (Sciurus vulgaris) and Ixodes ricinus in enzootic areas in Switzerland. Acta Trop. 1998, 69, 213–227. [Google Scholar] [CrossRef] [Green Version]
  20. Hanincová, K.; Tarageľová, V.; Koči, J.; Schäfer, S.M.; Hails, R.; Ullmann, A.J.; Piesman, J.; Labuda, M.; Kurtenbach, K. Association of Borrelia garinii and B. valaisiana with songbirds in Slovakia. Appl. Environ. Microbiol. 2003, 69, 2825–2830. [Google Scholar] [CrossRef] [Green Version]
  21. Tarageľová, V.; Koči, J.; Hanincová, K.; Kurtenbach, K.; Derdáková, M.; Ogden, N.H.; Literák, I.; Kocianová, E.; Labuda, M. Blackbirds and song thrushes constitute a key reservoir of Borrelia garinii, the causative agent of borreliosis in central Europe. Appl. Environ. Microbiol. 2008, 74, 1289–1293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Richter, D.; Matuschka, F.R. Perpetuation of the Lyme disease spirochete Borrelia lusitaniae by lizards. Appl. Environ. Microbiol. 2006, 72, 4627–4632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Majláthová, V.; Majláth, I.; Víchová, B.; Derdáková, M.; Peťko, B. Borrelia lusitaniae and green lizards (Lacerta viridis), Karst Region, Slovakia. Emerg. Infect. Dis. 2006, 12, 1895–1901. [Google Scholar] [CrossRef] [PubMed]
  24. Richter, D.; Schlee, D.B.; Allgöwer, R.; Matuschka, F.-R. Relationships of a novel Lyme disease spirochete, Borrelia spielmanii sp. nov., with its hosts in Central Europe. Appl. Environ. Microbiol. 2004, 70, 6414–6419. [Google Scholar] [CrossRef] [Green Version]
  25. Durand, J.; Herrmann, C.; Genné, D.; Sarr, A.; Gern, L.; Voordouw, M.J. Multistrain infections with Lyme borreliosis pathogens in the tick vector. Appl. Environ. Microbiol. 2017, 83, e02552-16. [Google Scholar] [CrossRef] [Green Version]
  26. Rudenko, N.; Grubhoffer, L.; Golovchenko, M. Climate Impact on Lyme Borreliosis and Its Causative Agents. In Climate, Ticks and Disease; Nuttall, P., Ed.; University of Oxford: Oxford, UK, 2021; pp. 349–354. [Google Scholar]
  27. Hanincova, K.; Kurtenbach, K.; Diuk-Wasser, M.; Brei, B.; Fish, D. Epidemic spread of Lyme borreliosis, Northeastern United States. Emerg. Infect. Dis. 2006, 12, 604–611. [Google Scholar] [CrossRef]
  28. Černý, J.; Lynn, G.; Hrnková, J.; Golovchenko, M.; Rudenko, N.; Grubhoffer, L. Management options for Ixodes ricinus-associated pathogens: A review of prevention strategies. Int. J. Environ. Res. Public Health 2020, 17, 1830. [Google Scholar] [CrossRef] [Green Version]
  29. Kahl, O.; Gray, J.S. The biology of Ixodes ricinus with emphasis on its ecology. Ticks Tick Borne Dis. 2023, 14, 102114. [Google Scholar] [CrossRef]
  30. Rizzoli, A.; Silaghi, C.; Obiegala, A.; Rudolf, I.; Hubálek, Z.; Földvári, G.; Plantard, O.; Vayssier-Taussat, M.; Bonnet, S.; Špitalská, E.; et al. Ixodes ricinus and its transmitted pathogens in urban and peri-urban areas in Europe: New hazards and relevance for public health. Front. Publ. Health 2014, 2, 251. [Google Scholar] [CrossRef]
  31. Gandy, S.; Kilbride, E.; Biek, R.; Millins, C.; Gilbert, L. Experimental evidence for opposing effects of high deer density on tick-borne pathogen prevalence and hazard. Parasit. Vectors 2021, 14, 509. [Google Scholar] [CrossRef]
  32. Stanko, M.; Derdáková, M.; Špitalská, E.; Kazimírová, M. Ticks and their epidemiological role in Slovakia: From the past till present. Biologia 2022, 77, 575–1610. [Google Scholar] [CrossRef] [PubMed]
  33. Kmety, E.; Rehácek, J.; Výrosteková, V. Investigations of ticks for the presence of Borrelia in Czechoslovakia. Zentralbl. Bakteriol. Mikrobiol. Hyg. A 1987, 263, 468–470. [Google Scholar] [CrossRef] [PubMed]
  34. Drgoňová, M.; Řeháček, J. Prevalence of Lyme borrelia in ticks in Bratislava, Slovak Republic. Cent. Eur. J. Public Health 1995, 3, 134–137. [Google Scholar]
  35. Gern, L.; Hu, C.M.; Kocianova, E.; Vyrostekova, V.; Rehacek, J. Genetic diversity of Borrelia burgdorferi sensu lato isolates obtained from Ixodes ricinus ticks collected in Slovakia. Eur. J. Epidemiol. 1999, 15, 665–669. [Google Scholar] [CrossRef]
  36. Štepánová-Tresová, G.; Peťko, B.; Štefančíková, A.; Nadzamová, D. Occurrence of Borrelia burgdorferi sensu stricto, Borrelia garinii and Borrelia afzelii in the Ixodes ricinus ticks from Eastern Slovakia. Eur. J. Epidemiol. 2000, 16, 105–109. [Google Scholar] [CrossRef] [PubMed]
  37. Lenčáková, D.; Hizo-Teufel, C.; Petko, B.; Schulte-Spechtel, U.; Stanko, M.; Wilske, B.; Fingerle, V. Prevalence of Borrelia burgdorferi s.l. OspA types in Ixodes ricinus ticks from selected localities in Slovakia and Poland. Int. J. Med. Microbiol. 2006, 296 (Suppl. S40), 108–118. [Google Scholar] [CrossRef]
  38. Pangrácová, L.; Derdáková, M.; Pekárik, L.; Hviščová, I.; Víchová, B.; Stanko, M.; Hlavatá, H.; Peťko, B. Ixodes ricinus abundance and its infection with the tick-borne pathogens in urban and suburban areas of Eastern Slovakia. Parasit. Vectors 2013, 6, 238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Rusňáková Tarageľová, V.; Mahríková, L.; Selyemová, D.; Václav, R.; Derdáková, M. Natural foci of Borrelia lusitaniae in a mountain region of Central Europe. Ticks Tick Borne Dis. 2016, 7, 350–356. [Google Scholar] [CrossRef]
  40. Chvostáč, M.; Špitalská, E.; Václav, R.; Vaculová, T.; Minichová, L.; Derdáková, M. Seasonal patterns in the prevalence and diversity of tick-borne Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Rickettsia spp. in an urban temperate forest in south western Slovakia. Int. J. Environ. Res. Public Health 2018, 15, 994. [Google Scholar] [CrossRef] [Green Version]
  41. Vaculová, T.; Derdáková, M.; Špitalská, E.; Václav, R.; Chvostáč, M.; Rusňáková Tarageľová, V. Simultaneous occurrence of Borrelia miyamotoi, Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Rickettsia helvetica in Ixodes ricinus ticks in urban foci in Bratislava, Slovakia. Acta Parasitol. 2019, 64, 19–30. [Google Scholar] [CrossRef]
  42. Mtierová, Z.; Derdáková, M.; Chvostáč, M.; Didyk, Y.M.; Mangová, B.; Rusňáková Tarageľová, V.; Selyemová, D.; Šujanová, A.; Václav, R. Local population structure and seasonal variability of Borrelia garinii genotypes in Ixodes ricinus ticks, Slovakia. Int. J. Environ. Res. Public Health. 2020, 17, 3607. [Google Scholar] [CrossRef]
  43. Subramanian, G.; Sekeyova, Z.; Raoult, D.; Mediannikov, O. Multiple tick-associated bacteria in Ixodes ricinus from Slovakia. Ticks Tick Borne Dis. 2012, 3, 406–410. [Google Scholar] [CrossRef] [PubMed]
  44. Hamšíková, Z.; Coipan, C.; Mahríková, L.; Minichová, L.; Sprong, H.; Kazimírová, M. Borrelia miyamotoi and co-infection with Borrelia afzelii in Ixodes ricinus ticks and rodents from Slovakia. Microb. Ecol. 2017, 73, 1000–1008. [Google Scholar] [CrossRef] [PubMed]
  45. Rosà, R.; Andreo, V.; Tagliapietra, V.; Baráková, I.; Arnoldi, D.; Hauffe, H.C.; Manica, M.; Rosso, F.; Blaňarová, L.; Bona, M.; et al. Effect of climate and land use on the spatio-temporal variability of tick-borne bacteria in Europe. Int. J. Environ. Res. Public Health 2018, 15, 732. [Google Scholar] [CrossRef] [Green Version]
  46. Svitálková, Z.; Haruštiaková, D.; Mahríková, L.; Berthová, L.; Slovák, M.; Kocianová, E.; Kazimírová, M. Anaplasma phagocytophilum prevalence in ticks and rodents in an urban and natural habitat in South-Western Slovakia. Parasit. Vectors 2015, 8, 276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Kazimírová, M.; Hamšíková, Z.; Kocianová, E.; Marini, G.; Mojšová, M.; Mahríková, L.; Berthová, L.; Slovák, M.; Rosá, R. Relative density of host-seeking ticks in different habitat types of south-western Slovakia. Exp. Appl. Acarol. 2016, 69, 205–224. [Google Scholar] [CrossRef]
  48. Špitalská, E.; Stanko, M.; Mošanský, L.; Kraljik, J.; Miklisová, D.; Mahríková, L.; Bona, M.; Kazimírová, M. Sesonal analysis of Rickettsia species in ticks in an agricultural site of Slovakia. Exp. Appl. Acarol. 2016, 68, 315–324. [Google Scholar] [CrossRef]
  49. Stanko, M.; Mošanský, L.; Fričová, J.; Čisláková, L.; Schniererová, E. Long-Time Study of Natural Focus and Its Structure in Eastern Slovakia. In Stawonogi Roznorodnosc Form i Oddzialywan; Buczek, A., Blaszak, C., Eds.; Koliber: Lublin, Poland, 2005; pp. 113–119. [Google Scholar]
  50. Siuda, K. Kleszcze Polski (Acari: Ixodida). Część II. Systematyka i Rozmieszczenie. In Monografie Parazytologiczne; Polskie Towarzystwo Parazytologiczne: Warszawa, Poland, 1993; pp. 1–375. [Google Scholar]
  51. Filippova, N.A. Fauna of the USSR: Arachnida, Subfamily Ixodinae 4(4); Nauka Publisher: Leningrad, Russia, 1977. [Google Scholar]
  52. Mydlová, L.; Svitálková, Z.; Slovák, M.; Kazimírová, M. Variability in Prevalence and Distribution of Borrelia burgdorferi s.l. In Ixodes ricinus Ticks across Bratislava Forest Park and Adjacent Sites (SW Slovakia). In Zborník z Konferencie III. Labudove dni; Kocianová, E., Kazimírová, M., Lúčanová, V., Eds.; Virologický ústav SAV: Bratislava, Slovakia, 2013; pp. 26–29. [Google Scholar]
  53. Derdáková, M.; Beati, L.; Peťko, B.; Stanko, M.; Fish, D. Genetic variability within Borrelia burgdorferi sensu lato genospecies established by PCR-single-strand conformation polymorphism analysis of the RrfA-RrlB intergenic spacer in Ixodes ricinus ticks from the Czech Republic. Appl. Environ. Microbiol. 2003, 69, 509–516. [Google Scholar] [CrossRef] [Green Version]
  54. Rudenko, N.; Golovchenko, M.; Grubhoffer, L.; Oliver, J.H., Jr. Borrelia carolinensis sp. nov., a new (14th) member of the Borrelia burgdorferi sensu lato complex from the southeastern region of the United States. J. Clin. Microbiol. 2009, 47, 134–141. [Google Scholar] [CrossRef] [Green Version]
  55. Bunikis, J.; Garpmo, U.; Tsao, J.; Berglund, J.; Fish, D.; Barbour, A.G. Sequence typing reveals extensive strain diversity of the Lyme borreliosis agents Borrelia burgdorferi in North America and Borrelia afzelii in Europe. Microbiology 2004, 150, 1741–1755. [Google Scholar] [CrossRef] [Green Version]
  56. Clark, K.; Hendricks, A.; Burge, D. Molecular identification and analysis of Borrelia burgdorferi sensu lato in lizards in the southeastern United States. Appl. Environ. Microbiol. 2005, 71, 2616–2625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Alekseev, A.N.; Dubinina, H.V.; Van De Pol, I.; Schouls, L.M. Identification of Ehrlichia spp. and Borrelia burgdorferi in Ixodes ticks in the Baltic regions of Russia. J. Clin. Microbiol. 2001, 39, 2237–2242. [Google Scholar] [CrossRef] [Green Version]
  58. Blazejak, K.; Raulf, M.K.; Janecek, E.; Jordan, D.; Fingerle, V.; Strube, C. Shifts in Borrelia burgdorferi (s.l.) geno-species infections in Ixodes ricinus over a 10-year surveillance period in the city of Hanover (Germany) and Borrelia miyamotoi-specific Reverse Line Blot detection. Parasit. Vectors 2018, 11, 304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Gern, L.; Douet, V.; López, Z.; Rais, O.; Cadenas, F.M. Diversity of Borrelia genospecies in Ixodes ricinus ticks in a Lyme borreliosis endemic area in Switzerland identified by using new probes for reverse line blotting. Ticks Tick Borne Dis. 2010, 1, 23–29. [Google Scholar] [CrossRef] [PubMed]
  60. Springer, A.; Jordan, D.; Glass, A.; Kahl, O.; Fingerle, V.; Girl, P.; Chitimia-Dobler, L.; Strube, C. Borrelia infections in ageing ticks: Relationship with morphometric age ratio in field-collected Ixodes ricinus nymphs. Microorganisms 2022, 10, 166. [Google Scholar] [CrossRef] [PubMed]
  61. Fukunaga, M.; Okada, K.; Nakao, M.; Konishi, T.; Sato, Y. Phylogenetic analysis of Borrelia species based on flagellin gene sequences and its application for molecular typing of Lyme disease borreliae. Int. J. Syst. Bacteriol. 1996, 46, 898–905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Postic, D.; Assous, M.V.; Grimont, P.A.; Baranton, G. 1994. Diversity of Borrelia burgdorferi sensu lato evidenced by restriction fragment length polymorphism of rrf (5S)-rrl (23S) intergenic spacer amplicons. Int. J. Syst. Bacteriol. 1994, 44, 743–752. [Google Scholar] [CrossRef] [Green Version]
  63. Wang, G.; van Dam, A.P.; Schwartz, I.; Dankert, J. 1999. Molecular typing of Borrelia burgdorferi sensu lato: Taxonomic, epidemiological, and clinical implications. Clin. Microbiol. Rev. 1999, 2, 633–653. [Google Scholar] [CrossRef] [Green Version]
  64. Hammer, O.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontol. Electron. 2001, 4, 1–9. [Google Scholar]
  65. Jouda, F.; Perret, J.L.; Gern, L. Density of questing Ixodes ricinus nymphs and adults infected by Borrelia burgdorferi sensu lato in Switzerland: Spatio-temporal pattern at a regional scale. Vector Borne Zoonotic Dis. 2004, 4, 23–32. [Google Scholar] [CrossRef] [Green Version]
  66. Reye, A.L.; Hübschen, J.M.; Sausy, A.; Muller, C.P. Prevalence and seasonality of tick-borne pathogens in questing Ixodes ricinus ticks from Luxembourg. Appl. Environ. Microbiol. 2010, 76, 2923–2931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Venclíková, K.; Betášová, L.; Sikutová, S.; Jedličková, P.; Hubálek, Z.; Rudolf, I. Human pathogenic borreliae in Ixodes ricinus ticks in natural and urban ecosystem (Czech Republic). Acta Parasitol. 2014, 59, 717–720. [Google Scholar] [CrossRef] [Green Version]
  68. Hönig, V.; Svec, P.; Halas, P.; Vavruskova, Z.; Tykalova, H.; Kilian, P.; Vetiskova, V.; Dornakova, V.; Sterbova, J.; Simonova, Z.; et al. Ticks and tick-borne pathogens in South Bohemia (Czech Republic)—Spatial variability in Ixodes ricinus abundance, Borrelia burgdorferi and tick-borne encephalitis virus prevalence. Ticks Tick Borne Dis. 2015, 6, 559–567. [Google Scholar] [CrossRef]
  69. Honig, V.; Carolan, H.E.; Vavruskova, Z.; Massire, C.; Mosel, M.R.; Crowder, C.D.; Rounds, M.A.; Ecker, D.J.; Ruzek, D.; Grubhoffer, L.; et al. Broad-range survey of vector-borne pathogens and tick host identification of Ixodes ricinus from Southern Czech Republic. FEMS Microbiol. Ecol. 2017, 93, fix129. [Google Scholar] [CrossRef]
  70. Raileanu, C.; Moutailler, S.; Pavel, I.; Porea, D.; Mihalca, A.D.; Savuta, G.; Vayssier-Taussat, M. Borrelia diversity and co-infection with other tick borne pathogens in ticks. Front. Cell. Infect. Microbiol. 2017, 7, 36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Kubiak, K.; Dziekońska-Rynko, J.; Szymańska, H.; Kubiak, D.; Dmitryjuk, M.; Dzika, E. Questing Ixodes ricinus ticks (Acari, Ixodidae) as a vector of Borrelia burgdorferi sensu lato and Borrelia miyamotoi in an urban area of north-eastern Poland. Exp. Appl. Acarol. 2019, 78, 113–126. [Google Scholar] [CrossRef] [Green Version]
  72. Zubriková, D.; Wittmann, M.; Hönig, V.; Švec, P.; Víchová, B.; Essbauer, S.; Dobler, G.; Grubhoffer, L.; Pfister, K. Prevalence of tick-borne encephalitis virus and Borrelia burgdorferi sensu lato in Ixodes ricinus ticks in Lower Bavaria and Upper Palatinate, Germany. Ticks Tick Borne Dis. 2020, 11, 101375. [Google Scholar] [CrossRef] [PubMed]
  73. Hartemink, N.; van Vliet, A.J.H.; Gort, G.; Gassner, F.; Jacobs, F.; Fonville, M.; Takken, W.; Sprong, H. Seasonal patterns and spatial variation of Borrelia burgdorferi (sensu lato) infections in Ixodes ricinus in the Netherlands. Parasit. Vectors 2021, 14, 121. [Google Scholar] [CrossRef] [PubMed]
  74. Răileanu, C.; Silaghi, C.; Fingerle, V.; Margos, G.; Thiel, C.; Pfister, K.; Overzier, E. Borrelia burgdorferi sensu lato in questing and engorged ticks from different habitat types in southern Germany. Microorganisms 2021, 9, 1266. [Google Scholar] [CrossRef]
  75. Knoll, S.; Springer, A.; Hauck, D.; Schunack, B.; Pachnicke, S.; Fingerle, V.; Strube, C. Distribution of Borrelia burgdorferi s.l. and Borrelia miyamotoi in Ixodes tick populations in Northern Germany, co-infections with Rickettsiales and assessment of potential influencing factors. Med. Vet. Entomol. 2021, 35, 595–606. [Google Scholar] [CrossRef]
  76. Król, N.; Obiegala, A.; Imholt, C.; Arz, C.; Schmidt, E.; Jeske, K.; Ulrich, R.G.; Rentería-Solís, Z.; Jacob, J.; Pfeffer, M. Diversity of Borrelia burgdorferi sensu lato in ticks and small mammals from different habitats. Parasit. Vectors 2022, 15, 195. [Google Scholar] [CrossRef]
  77. 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–7216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Trávniček, M.; Guryčová, D.; Halánová, M.; Kožuch, O.; Nadzamová, D.; Miško, J. Presence of antibodies to some zoonoses in mouflons (Ovis musimon Pall.) and fallow deer (Dama dama L.) in eastern Slovakia. Vet. Med. 1999, 44, 215–219. [Google Scholar]
  79. Štefančíková, A.; Derdáková, M.; Lenčáková, D.; Ivanová, R.; Stanko, M.; Čisláková, L.; Peťko, B. Serological and molecular detection of Borrelia burgdorferi sensu lato and Anaplasmataceae in rodents. Folia Microbiol. 2008, 53, 493–499. [Google Scholar] [CrossRef] [PubMed]
  80. Nadzamová, D.; Petko, B.; Stefancíková, A.; Cisláková, L.; Stĕpánová, G.; Seba, M. Potential foci of Lyme borreliosis in towns. Cent. Eur. J. Public Health 2000, 8, 176–178. [Google Scholar] [PubMed]
  81. Peťko, B.; Štefančíková, A.; Nadzamová, D.; Stanko, M.; Siuda, K.; Karbowiak, G.; Winiarczyjk, S. 2001. Epizootiological Aspects of Lyme Borreliosis in the City Agglomerations of the Carpathian Regions of Slovakia and Poland and Their Peripheral Parts. In Stawonogi Pasożyty i Nosiciele; Buczek, A., Błaszak, C., Eds.; Wyd. KGM: Lublin, Poland, 2001; pp. 157–164. [Google Scholar]
  82. Derdakova, M.; Halanova, M.; Stanko, M.; Stefancikova, A.; Cislakova, L.; Pet’ko, B. Molecular evidence for Anaplasma phagocytophilum and Borrelia burgdorferi sensu lato in Ixodes ricinus ticks from eastern Slovakia. Ann. Agric. Environ. Med. 2003, 10, 269–271. [Google Scholar]
  83. Smetanová, K.; Schwarzová, K.; Kocianová, E. Detection of Anaplasma phagocytophilum, Coxiella burnetti, Rickettsia spp., and Borrelia burgdorferi s.l. in ticks and wild living animals in western and middle Slovakia. Ann. N. Y. Acad. Sci. 2006, 1078, 312–315. [Google Scholar] [CrossRef]
  84. Smetanová, K.; Burri, C.; Pérez, D.; Gern, L.; Kocianová, E. Detection and identification of Borrelia burgdorferi sensu lato genospecies in ticks from three different regions in Slovakia. Wien. Klin. Wochenschr. 2007, 119, 534–537. [Google Scholar] [CrossRef] [Green Version]
  85. Danielová, V.; Daniel, M.; Schwarzová, L.; Materna, J.; Rudenko, N.; Golovchenko, M.; Holubová, J.; Grubhoffer, L.; Kilián, P. Integration of a tick-borne encephalitis virus and Borrelia burgdorferi sensu lato into mountain ecosystems, following a shift in the altitudinal limit of distribution of their vector, Ixodes ricinus (Krkonose mountains, Czech Republic). Vector Borne Zoonotic Dis. 2010, 10, 223–230. [Google Scholar] [CrossRef]
  86. Picken, M.M.; Cheng, Y.; Strle, F.; Picken, R.N. Patient isolates of Borrelia burgdorferi sensu lato with genomic and phenotypic similarities of strain 25015. J. Infect. Dis. 1996, 174, 1112–1115. [Google Scholar] [CrossRef]
  87. Strle, F.; Picken, R.N.; Cheng, Y.; Cimperman, J.; Maraspin, V.; Lotric-Furlan, S.; Ruzic-Sabljic, E.; Picken, M.M. Clinical findings for patients with Lyme borreliosis caused by Borrelia burgdorferi sensu lato with genomic and phenotypic similarities of strain 25015. Clin. Infect. Dis. 1997, 25, 273–280. [Google Scholar] [PubMed] [Green Version]
  88. Hansford, K.M.; Wheeler, B.W.; Tschirren, B.; Medlock, J.M. Questing Ixodes ricinus ticks and Borrelia spp. in urban green space across Europe: A review. Zoonoses Public Health 2022, 69, 153–166. [Google Scholar] [CrossRef] [PubMed]
  89. Hornok, S.; Meli, M.L.; Gönczi, E.; Halász, E.; Takács, N.; Farkas, R.; Hofmann-Lehmann, R. Occurrence of ticks and prevalence of Anaplasma phagocytophilum and Borrelia burgdorferi s.l. in three types of urban biotopes: Forests, parks and cemeteries. Ticks Tick Borne Dis. 2014, 5, 785–789. [Google Scholar] [CrossRef] [PubMed]
  90. Richtrová, E.; Míchalová, P.; Lukavská, A.; Navrátil, J.; Kybicová, K. Borrelia burgdorferi sensu lato infection in Ixodes ricinus ticks in urban green areas in Prague. Ticks Tick Borne Dis. 2022, 13, 102053. [Google Scholar] [CrossRef] [PubMed]
  91. Glass, A.; Springer, A.; Raulf, M.K.; Fingerle, V.; Strube, C. 15-year Borrelia prevalence and species distribution monitoring in Ixodes ricinus/inopinatus populations in the city of Hanover, Germany. Ticks Tick Borne Dis. 2023, 14, 102074. [Google Scholar] [CrossRef]
  92. Rollins, R.E.; Yeyin, Z.; Wyczanska, M.; Alig, N.; Hepner, S.; Fingerle, V.; Margos, G.; Becker, N.S. Spatial variability in prevalence and genospecies distributions of Borrelia burgdorferi sensu lato from ixodid ticks collected in southern Germany. Ticks Tick Borne Dis. 2021, 12, 101589. [Google Scholar] [CrossRef]
  93. Kowalec, M.; Szewczyk, T.; Welc-Faleçiak, R.; Siński, E.; Karbowiak, G.; Bajer, A. Ticks and the city—Are there any differences between city parks and natural forests in terms of tick abundance and prevalence of spirochaetes? Parasit. Vectors 2017, 10, 573. [Google Scholar] [CrossRef]
  94. Gray, J.S.; Kahl, O.; Janetzki, C.; Stein, J. Studies on the ecology of Lyme disease in a deer forest in County Galway, Ireland. J. Med. Entomol. 1992, 29, 915–920. [Google Scholar] [CrossRef]
  95. Zeman, P.; Daniel, M. Mosaic pattern of Borrelia infection in a continuous population of the tick Ixodes ricinus (Acari: Ixodidae). Exp. Appl. Acarol. 1999, 23, 327–335. [Google Scholar] [CrossRef]
Figure 1. Tick collection sites in Slovakia. Bratislava—urban/suburban habitat, Fúgelka—natural habitat, Rozhanovce—agricultural habitat.
Figure 1. Tick collection sites in Slovakia. Bratislava—urban/suburban habitat, Fúgelka—natural habitat, Rozhanovce—agricultural habitat.
Microorganisms 11 01666 g001
Figure 2. Prevalence of Borrelia burgdorferi s.l. in questing Ixodes ricinus nymphs and adults per habitat type. Bars represent prevalence with 95% confidence intervals. Significance level of the χ2 test: *** p < 0.0001, ns—non-significant.
Figure 2. Prevalence of Borrelia burgdorferi s.l. in questing Ixodes ricinus nymphs and adults per habitat type. Bars represent prevalence with 95% confidence intervals. Significance level of the χ2 test: *** p < 0.0001, ns—non-significant.
Microorganisms 11 01666 g002
Figure 3. Differences in prevalence of Borrelia burgdorferi s.l. between questing Ixodes ricinus nymphs and adults in three habitat types. Bars represent prevalence with 95% confidence intervals. Significance level of the χ2 test: * p < 0.05, *** p < 0.0001, ns—non-significant.
Figure 3. Differences in prevalence of Borrelia burgdorferi s.l. between questing Ixodes ricinus nymphs and adults in three habitat types. Bars represent prevalence with 95% confidence intervals. Significance level of the χ2 test: * p < 0.05, *** p < 0.0001, ns—non-significant.
Microorganisms 11 01666 g003
Figure 4. Proportion of Borrelia burgdorferi s.l. species (in %) in infected Ixodes ricinus nymphs and adults per habitat and year. Species: ss—sensu stricto, l—lusitaniae, a—afzelii, g—garinii (+ bavariensis), v—valaisiana, s—spielmanii, coinf—coinfections.
Figure 4. Proportion of Borrelia burgdorferi s.l. species (in %) in infected Ixodes ricinus nymphs and adults per habitat and year. Species: ss—sensu stricto, l—lusitaniae, a—afzelii, g—garinii (+ bavariensis), v—valaisiana, s—spielmanii, coinf—coinfections.
Microorganisms 11 01666 g004
Table 1. Prevalence of Borrelia burgdorferi s.l. in Ixodes ricinus per habitat type in 2011–2013.
Table 1. Prevalence of Borrelia burgdorferi s.l. in Ixodes ricinus per habitat type in 2011–2013.
2011 2012 2013 Total
Site % (pos/ex)95% CI% (pos/ex)95% CI% (pos/ex)95% CIχ2p% (pos/ex)95% CI
BratislavaNymphs6.9 (61/887)5.4–8.711.6 (23/199)7.1–16.01.8 (8/448)0.6–3.026.227<0.00016.0 (92/1534)4.8–7.2
(urban/suburban habitat)Females18.6 (69/370)14.7–22.614.5 (9/62)5.7–23.317.7 (28/158)11.8–23.70.6240.73218.0 (106/590)14.9–21.1
Males9.9 (44/443)7.1–12.713.2 (9/68)5.2–21.38.29 (15/181)4.3–12.31.3810.5019.8 (68/692)7.6–12.0
Adults total13.9 (113/813)11.5–16.313.8 (18/130)7.9–19.812.3 (43/339)9.1–16.20.3100.85613.6 (174/1282)11.7–15.4
Total10.2 (174/1700)8.8–11.712.5 (41/329)8.9–16.06.5 (51/787)4.8–8.212.8290.0029.4 (266/2816)8.3–10.5
FúgelkaNymphs17.0 (184/1084)14.7–19.219.5 (58/298)15.0–24.08.9 (21/237)5.2–12.512.1900.00216.2 (263/1619)14.4–18.0
(natural habitat)Females17.0 (26/153)11.0–22.925.0 (15/60)14.0–36.016.4 (12/73)7.9–24.92.1140.34718.5 (53/286)14.0–23.0
Males7.2 (12/167)3.3–11.112.0 (10/83)5.0–19.08.3 (8/96)2.8–13.91.6740.4338.7 (30/346)5.7–11.6
Adults total11.9 (38/320)8.3–15.417.5 (25/143)11.3–23.711.8 (20/169)7.0–16.73.0650.21613.1 (83/632)10.5–15.8
Total15.8 (222/1404)13.9–17.718.8 (83/441)15.2–22.510.1 (41/406)7.2–13.012.9210.00215.4 (346/2251)13.9–16.9
RozhanovceNymphs15.4 (41/266)11.1–19.712.9 (56/433)9.8–16.113.0 (20/154)7.7–18.30.9410.62513.7 (117/853)11.4–16.0
(agricultural habitat)Females21.2 (7/33)7.3–35.227.3 (9/33)12.1–42.58.8 (3/34)0.0–18.4 0.119 F19.0 (19/100)11.3–26.7
Males13.3 (4/30)1.2–25.525.6 (11/43)12.5–38.66.7 (2/30)0.0–15.6 0.101 F16.5 (17/103)9.3–23.7
Adults total17.5 (11/63)8.1–26.826.3 (20/76)16.4–36.27.8 (5/64)1.2–14.48.1580.01717.7 (36/203)12.5–23.0
Total15.8 (52/329)11.9–19.714.9 (76/509)11.8–18.011.5 (25/218)7.2–15.72.1460.34214.5 (153/1056)12.4–16.6
(pos/ex)—number of positive/number of examined ticks; 95% CI—confidence interval; p—significance level of the χ2 test, F Fisher’s exact test, significant differences are indicated in bold.
Table 2. Representation of Borrelia burgdorferi s.l. species in Ixodes ricinus populations from three habitat types.
Table 2. Representation of Borrelia burgdorferi s.l. species in Ixodes ricinus populations from three habitat types.
Speciesss a g v l s coinf T
%95% CI%95% CI%95% CI%95% CI%95% CI%95% CI%95% CI
Bratislava
20110.10–0.34.53.5–5.43.12.3–3.90.60.2–0.91.10.6–1.60.80.3–1.20.060–0.21700
20120.60–1.455.53.0–7.92.10.6–3.73.31.4–5.390.90–1.90 0 329
20130.60–1.21.30.5–2.01.60.8–2.52.41.3–3.50.40–0.8 0.10–0.4787
Total0.30.1–0.53.73.0–4.42.62.0–3.21.41.0–1.90.90.5–1.20.50.2–0.70.10–0.22816
Fúgelka
20110.60.2–1.16.54.8–7.36.14.9–7.42.21.4–3.00.10–0.20.10–0.20.60.2–1.11404
20120 5.03.0–7.08.86.2–11.52.30.9–3.72.71.2–4.20 0 441
20130 3.21.5–4.95.43.2–7.60.70–1.60.20–0.70.20–0.70.20–0.7406
Total0.40.1–0.75.34.4–6.36.55.5–7.51.91.4–2.50.60.3–0.90.10–0.20.40.2–0.72251
Rozhanovce
20110.60–1.48.85.7–11.93.31.4–5.3 0.90–1.90.30–0.91.80.4–3.3329
20121.80.6–2.94.12.4–5.83.31.8–4.93.72.1–5.41.00.1–1.80 1.00.1–1.8509
20130.90–2.21.40–2.94.61.8–7.44.61.8–7.40 0 0 218
Total1.20.6–1.95.23.7–6.33.62.5–4.72.71.8–3.70.80.2–1.30.10–0.31.40.4–1.61056
TOTAL0.50.3–0.74.54.0–5.04.23.7–4.71.81.5–2.20.80.5–1.00.30.1–0.40.40.2–0.56123
Borrelia burgdorferi s.l. species: ss—sensu stricto, l—lusitaniae, a—afzelii, g—garinii (+ bavariensis), v—valaisiana, s—spielmanii, coinf—coinfections; 95% CI—confidence interval; T—total number of examined ticks.
Table 3. Prevalence of Borrelia burgdorferi s.l. in Ixodes ricinus per transects A–E in Bratislava in 2011.
Table 3. Prevalence of Borrelia burgdorferi s.l. in Ixodes ricinus per transects A–E in Bratislava in 2011.
TransectNymphs Adults Total %N
% (pos/ex)95% CI% (pos/ex)95% CI% (pos/ex)95% CI
A2.6 (10/381)1.0–4.26.8 (6/88)1.5–12.13.4 (16/469)1.8–5.081.2
B13.9 (25/180)8.8–18.98.6 (25/291)5.4–11.810.6 (50/471)7.8–13.438.2
C14.3 (2/14)0.0–32.621.0 (41/195)15.3–26.720.6 (43/209)15.1–26.06.7
D25.0 (19/76)15.3–34.725.4 (44/173)18.9–31.925.3 (63/249)19.9–30.730.5
E5.1 (12/233)2.3–8.012.7 (8/63)4.5–20.96.8 (20/296)3.9–9.678.7
χ258.554 33.947 101.798
p<0.0001 <0.0001 <0.0001
(pos/ex)—number of positive/number of examined ticks; 95% CI—confidence interval; %N—percentage share of nymphs in the total number of collected ticks; p—significance level of the χ2 test, significant differences are indicated in bold.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kazimírová, M.; Mahríková, L.; Hamšíková, Z.; Stanko, M.; Golovchenko, M.; Rudenko, N. Spatial and Temporal Variability in Prevalence Rates of Members of the Borrelia burgdorferi Species Complex in Ixodes ricinus Ticks in Urban, Agricultural and Sylvatic Habitats in Slovakia. Microorganisms 2023, 11, 1666. https://doi.org/10.3390/microorganisms11071666

AMA Style

Kazimírová M, Mahríková L, Hamšíková Z, Stanko M, Golovchenko M, Rudenko N. Spatial and Temporal Variability in Prevalence Rates of Members of the Borrelia burgdorferi Species Complex in Ixodes ricinus Ticks in Urban, Agricultural and Sylvatic Habitats in Slovakia. Microorganisms. 2023; 11(7):1666. https://doi.org/10.3390/microorganisms11071666

Chicago/Turabian Style

Kazimírová, Mária, Lenka Mahríková, Zuzana Hamšíková, Michal Stanko, Maryna Golovchenko, and Natalie Rudenko. 2023. "Spatial and Temporal Variability in Prevalence Rates of Members of the Borrelia burgdorferi Species Complex in Ixodes ricinus Ticks in Urban, Agricultural and Sylvatic Habitats in Slovakia" Microorganisms 11, no. 7: 1666. https://doi.org/10.3390/microorganisms11071666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop