Ixodes ricinus abundance and its infection with the tick-borne pathogens in urban and suburban areas of Eastern Slovakia
© Pangrácová et al.; licensee BioMed Central Ltd. 2013
Received: 17 April 2013
Accepted: 14 August 2013
Published: 16 August 2013
Raising abundance of ticks and tick-borne diseases in Europe is the result of multiple factors including climate changes and human activities. Herein, we investigated the presence and seasonal activity of Ixodes ricinus ticks from 10 urban and suburban sites in two different geographical areas of southeastern and northeastern Slovakia during 2008–2010. Our aim was to study the abundance of ticks in correlation with the environmental factors and their infection with Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Neoehrlichia mikurensis.
Questing I. ricinus ticks were collected from ten urban and suburban sites in Eastern Slovakia. A total of 670 ticks were further analysed for the presence of B. burgdorferi s.l., A. phagocytophilum and N. mikurensis by molecular methods. Tick site and environmental relations were analysed using General Linear Models (LM). The differences between the number of Lyme borreliosis cases between the Košice and Bardejov regions during a ten-year period were tested by Wilcoxon matched pairs test.
In total, 2921 (1913 nymphs, 1008 adults) I. ricinus ticks were collected from 10 study sites during the main questing season. Tick activity and relative abundance differed between locations and months. Temperature and humidity were the main factors affecting the tick abundance and questing activity. Out of 670 examined ticks, 10.15% were infected with spirochetes from B. burgdorferi s.l. complex (represented by B. afzelii, B. garinii, B.valaisiana and B. burgdorferi s.s.), 2.69% with the A. phagocytophilum and 2.39% with N. mikurensis. The number of Lyme borreliosis cases per 100,000 inhabitants in the Bardejov region was significantly higher than in the Košice region.
Our data indicate that the risk of infection with tick-borne pathogens in Eastern Slovakia is common since 15.2% of ticks were infected at least with one of the tested microorganisms. Even though the abundance of ticks was affected by the microclimatic conditions and the prevalence of pathogens differed between the habitats, the infection risk for humans is also affected by human activities leading to an increased contact with infected ticks.
KeywordsIxodes ricinus Borrelia burgdorferi sensu lato Anaplasma phagocytophilum Neoehrlichia mikurensis PCR-RFLP Lyme borreliosis Anaplasmosis
In Europe, the changing climate and human activities in the environment have caused the changes in tick abundance together with the spread of ticks into the northern regions, urban and suburban areas as well as higher altitudes. This phenomenon is associated with the spread of tick-borne pathogens and new foci in the areas, previously free of the tick-borne diseases, which have been established [1–5]. Abiotic factors of the microclimate such as temperature, humidity, saturation or vapour pressure deficit and wind influence the survival of ticks and their questing behaviour in the habitat [6–8]. Structure of the habitat and the host availability for ticks also largely influence their phenology. The epidemiologically most important tick in Europe, Ixodes ricinus, transmits viral, bacterial as well as protozoan pathogens to humans and animals. The most commonly occurring and the most serious bacterial agents transmitted by this tick in Europe are spirochetes from the Borrelia burgdorferi sensu lato complex. They are causative agents of Lyme borreliosis, the multisystemic disorder that is maintained in natural foci in a wide spectrum of vertebrate reservoir hosts . Currently 19 different genospecies belong to this complex out of which at least 9 are present in Europe . The specific associations of different genospecies with the reservoir hosts as well as clinical symptomatics have been assigned , however, this association is not strict and differences have been observed . The occurrence of Lyme borreliosis has been reported from various habitats of Europe between 35 to 60°N with the focal distribution even within small countries and thus following the occurrence of ticks . The highest yearly incidence is in Central Europe, namely in Austria and Slovenia, with 130 and 136 cases per 100 000 inhabitants . Another bacterial zoonotic disease, transmitted by I. ricinus, is granulocytic anaplasmosis caused by Anaplasma phagocytophilum. Anaplasmosis is the common tick-borne bacterial disease of domestic animals often causing tick-borne fever in ruminants on pastures . In Europe, human cases are less common than in US [15, 16]. The prevalence of A. phagocytophilum in questing ticks in Slovakia varies from 1.1 to 7.8% . Neoehrlichia mikurensis is another tick-borne pathogen from the family Anaplasmataceae that attracts the attention of public health professionals in Europe. It was detected in questing ticks throughout Europe [18–21]. Rodents have been proposed as potential reservoir hosts since it was detected in blood and endothelial cells of their spleens and livers [22, 23]. Recently its pathogenicity in humans was reported as it was detected in patients with septicaemia and immunosuppressed patients [24–27]. Moreover, it was detected in a chronically neutropenic dog from Germany .
The main aim of this study was to investigate the abundance and activity of I. ricinus ticks in urban and suburban areas of two cities in Eastern Slovakia in relation to the tick habitat and environmental conditions. The infection rates with the three most important tick-borne bacterial pathogens (B. burgdoferi s.l., A. phagocytophilum and N. mikurensis) were investigated as well. Furthermore, the occurrence of ticks and the presence of Borrelia were analysed in conjunction with the incidence of human cases reported to the State Health Institute during the last 10 years from both studied regions.
Collection of ticks
Description of tick collection sites
425 m a.s.l.
-Birch-beech dry forest with hornbean shrubs
345 m a.s.l.
-Maple-oak forest with shruby humid vegetation around pathways, close to the agricultural land
310 m a.s.l.
-Beech-oak dry forest,
332 m a.s.l.
-Suburban beech forest with shruby vegetation, in close proximity of cattle pastures
5 Bardejovské kúpele
283 m a.s.l.
-Urban beech-oak forest with park recultivation in some areas
321 m a.s.l.
- Hornbeam suburban forest with shrubby vegetation
200 m a.s.l.
- Urban park with large open areas without trees
8 Botanická záhrada
208 m a.s.l.
- Urban hornbeam-oak park with shrubs
9 Verejný cintorín
200 m a.s.l.
- Urban park at the cemetery
192 m a.s.l.
- Urban hornbeam forest
Tick site and environmental relations were analysed using General Linear Models (LM). Basic models (with only fixed terms) were extended by random terms to deal with replications in the samples (repeated samples per sites). The response variable, tick relative abundance was (log + 1) transformed to meet the requirements of normal error distribution. To assess the influence of selected variables on tick relative abundance, a model containing fixed independent terms, humidity, temperature, saturation deficit, vapour pressure deficit, elevation and biotope type and their quadratic forms and random term site identity was constructed. To select the best model, the least significant variables were removed from the model and new models were refitted. In the case that the fitting method based on restricted maximum likelihood (REML) did not allow model comparison, models were refitted by the maximum likelihood method. These procedures were repeated until the model changed significantly and until the Akaike information criterion (AIC) was decreasing. To assess the differences between sites in tick relative abundance a model containing fixed term site and random term month was constructed. If the results of contrasts (default treatment contrast) showed similar estimates and standard errors for particular sites, these sites were pooled and the model was refitted again. The best model selection followed the steps described above. All modelling procedures followed the approach described in Zuur et al.. The differences between the number of Lyme borreliosis cases between the Košice and Bardejov regions during a ten-year period were tested by Wilcoxon matched pairs test. Data on Lyme borreliosis case incidence in studied regions (Košice and Bardejov) were obtained from the Epidemiological Information System of Slovakia . Confidence intervals (CI) for infection rates were obtained by Exact Binomial test. All statistical analyses were performed in R statistical software environment  and R package nlme .
Molecular identification of tick-borne pathogens
Genomic DNA from each tick was isolated after its removal from ethanol and drying on filter paper by the method of alkaline-hydrolysis using 1.25% of ammonium solution, according to the previously described protocol . Each tick was homogenized with a sterile pestle and negative extraction controls containing only ammonium solution were prepared for each set of DNA extraction to monitor the possible contamination. Extracted DNA from ticks was further analysed for the presence of B. burgdorferi s.l. complex by amplification of a 250 bp long fragment of 5S-23S (rrfA-rrlB) rDNA intergenic spacer using primers IgsA (5′CGACCTTCTTCGCCTTAAAGC′3) and IgsB (5′AGCTCTTATTCGCTGATGTA′3) . For the molecular detection of A. phagocytophilum, PCR amplification of a 849 bp long fragment of msp4 gene was used, with primers MAP4Ap5 (5′ATGAATTACAGAGAATTGCTTGTAGG′3) and MSP4Ap3 (5′TTAATTGAAAGCAAA TCTTGCTCCTATG′3) . To detect N. mikurensis a 560 bp long fragment of 16S rRNA gene was amplified using IS58-594r (5′CTATCCTCTCTCGATCTCTAGT′3) and IS58-62f (5′GGAATAGCTGTTAGAAATGAC′3) primers . MasterTaq DNA polymerase kit (Eppendorf AG, Hamburg, Germany) was used for PCR amplifications. A total volume of 25 μl of reaction mixture consisted of: 2.5 μl template DNA (sample), 7.6 μl of nuclease free water, 12.5 μl of PCR Master mix, and 1.2 μl of each primer (10 pmole/μl). Positive and negative controls were used in each PCR reaction. The PCR products were electrophoresed on 2% agarose gels stained with GoldView Nucleic Acid Stain (Beijing SBS Genetech, Beijing, China). Amplified fragments were visualised in a transilluminator under UV light.
In the case of Borrelia positive ticks, samples were further assigned to the different genospecies by RFLP method using Tru1 restriction endonuclease (Fermentas, Vilnius, Lithuania) as described before .
Selected positive PCR products of the 16S rDNA fragment of N. mikurensis were purified by using a QIAquick PCR purification kit (Qiagen, Hilden, Germany) and sequenced in both directions with the same primers as for the PCR amplifications. Sequencing was performed at the University of Veterinary Medicine and Pharmacy in Košice; Department of Microbiology and Immunology. The complementary strands of each sequenced product were manually assembled. Sequences were compared to GenBank entries by Blast N2.2.13 . The GenBank accession number for the nucleotide sequences of partial 16S rDNA of N. mikurensis is JN378917.
A relative density (RD) of ticks per one hour of collection of total sampling at model sites, a percentage proportion of developmental stage and sex
RD of ticks
% of nymphs
% of females
% of males
5 B. kúpele
8 Botanická záhrada
9 Verejný cintorín
In the northeastern suburban forest habitat of Bardejov (site no. 2), seasonal activity of ticks also had a unimodal pattern with the peak in May (Figure 3).
Infection rate and 95% confidence interval (CI) of N. mikurensis, A. phagocytophilum and B. burgdorferi s.l. in I. ricinus ticks from sampling sites in Slovakia
No. of ticks tested
B. burgdorferi s.l.
5 B. kúpele
8 Botanical garden
9 Verejný cintorín
I. ricinus ticks are widely distributed in moderate climatic regions of Europe in both natural and urban habitats. The occurrence and recent expansion of ticks into new areas are limited by temperature and saturation deficit [5, 7, 38–40]. The abundance of ticks in our sites correlated to the humidity and temperature. Neither saturation deficit nor vapour pressure deficit was significant. This might be due to the differences between the microclimatic conditions at our sites and data obtained from the meteorological stations as previously observed . The seasonal activity of ticks can be unimodal with one maximum peak usually in late spring or early summer or bimodal with maximum peaks in spring or summer [7, 40]. We have observed unimodal patterns for all of our sites. The tick activity in northeastern sites had a maximum peak in May, one month later than for southeastern sites. This correlates with the lower temperature increase in the north as one of the significant environmental factors affecting the tick abundance observed in our models. In neighbouring Hungary, Egyed et al. reported bimodal activity for all their sites. In the statistical model our tick sites grouped into three “site groups”– A, B, C that represented the sites with the similar tick abundance and seasonal activities. Tick group “A” had the highest abundance of ticks and only one site belonged to this group – a dense suburban forest with shrubby vegetation in southeastern Slovakia. Group “B” consisted of sites represented by urban parks in southeastern Slovakia where the vegetation was more fragmented and forested sites from the northeastern Slovakia. Third group “C” grouped together sites with the least favourable conditions for tick abundace –maintained urban park with large open spaces in southeastern Slovakia and dry suburban and urban forest in the northeast. Grouping of sites into three categories according to tick abundance showed that even small differences in the latitude (southern site vs. northern sites) with the lower daily mean temperature (Figure 1) can affect the tick abundance. Generally, less ticks were found in northeastern Slovakia in appropriate tick habitats as opposed to south.
Interestingly, the number of Lyme borreliosis cases per 100 000 inhabitants were higher for the Bardejov region in northeastern Slovakia than for Košice in southeastern Slovakia. This is probably due to larger rural areas and different outdoor human behaviour patterns in the district of Bardejov, even though Košice is the second largest city in Slovakia. The link between human activities and incidence of tick-borne diseases has already been highlighted in previous studies [41, 42]. Positive correlation between the abundance of ticks and seroprevalence against borrelia and TBE was observed among farmers in neighbouring Poland .
Questing ticks in our study were infected with all tested zoonotic bacteria with the dominance of B. burgdorferi s.l. as it was detected in 10.15% (CI: 7.95-12.69) of ticks. This is in agreement with the data from Eastern Slovakia obtained by Lenčáková et al.  where 11% of ticks were Borrelia positive. Similar infectious rates were detected in I. ricinus ticks from Estonia . Prevalence of Borrelia in neighbouring countries in Hungary  and Poland  was slightly lower. In our dataset, the highest infection rate (18%) was detected in a suburban forest in Košice, southeastern Slovakia, where the highest abundance of ticks was also recorded. Moreover, at this locality the highest diversity of Borrelia species was observed; probably due to a higher availability of hosts than in urban parks within the area. In the European countries, the infection rate of A. phagocytophilum infection in ticks is generally low. The results from the study in 11 sites in Switzerland showed 1.5% infection rate and patchy distribution . We obtained similar results with 2.69% (CI: 1.6-4.2%) infection rates and it was detected at four out of ten sites. In contrast to Norway, at the areas with the higher density of the red deer, the prevalence of A. phagocytophilum was more consistent and higher (8.8%) . N. mikurensis, the recently emerging pathogen, was detected in 2.23% (CI:1.37-3.85%) of ticks. Its distribution was, however, more homogenous than that for A. phagocytophilum, as it was detected in all northeastern sites in Bardejov. Recent studies show that N. mikurensis is common and frequently infects I. ricinus that is widely distributed in Europe [20, 21].
Our data indicate that the risk of infection with tick-borne pathogens in Eastern Slovakia is common since 15.2% of ticks were infected with at least with one of the tested microorganism. Even though the abundance of ticks was affected by the microclimatic conditions and the prevalence of pathogens differed between the habitats, the infection risk for humans is also affected by human activities leading to an increased contact with infected ticks.
We thank Slavka Barláková for the English editing of the manuscript, Viktória Dandárová and Lucia Bajáčková for their help with collecting ticks. The study was supported by the projects VEGA 2/0055/11, VEGA 2/0137/10, by the Slovak Research and Development Agency under contract No. APVV – 0267–10. This contribution is partially the result of the project implementation: Development of the diagnostic methods for the detection of tick-borne pathogens and the techniques for the preparation of the vaccine development (code ITMS: 26240220044), supported by the Research & Development Operational Programme funded by the ERDF.
- Lukáň M, Bullová E, Peťko B: Climate warming and tick-borne encephalitis, Slovakia. Emerg Infect Dis. 2010, 16: 524-526. 10.3201/eid1603.081364.PubMed CentralView ArticlePubMedGoogle Scholar
- Bullová E, Lukáň M, Stanko M, Peťko B: Spatial distribution of Dermacentor reticulatus tick in Slovakia in the beginning of the 21st century. Vet Parasitol. 2009, 165: 357-360. 10.1016/j.vetpar.2009.07.023.View ArticlePubMedGoogle Scholar
- Materna J, Daniel M, Danielová V: Altitudinal distribution limit of the tick Ixodes ricinus shifted considerably towards higher altitudes in central Europe: results of three years monitoring in the Krkonose Mts. (Czech Republic). Cent Eur J Public Health. 2005, 13: 24-28.PubMedGoogle Scholar
- Jaenson TGT, Jaenson DGE, Eisen L, Petersson E, Lindgren E: Changes in the geographical distribution and abundance of the tick Ixodes ricinus during the past 30 years in Sweden. Parasit Vectors. 2012, 5: 1-15. 10.1186/1756-3305-5-1.View ArticleGoogle Scholar
- Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Peña A, George JC, Golovljova I, Jaenson TG, Jensen JK, Jensen PM, Kazimirova M, Oteo JA, Papa A, Pfister K, Plantard O, Randolph SE, Rizzoli A, Santos-Silva MM, Sprong H, Vial L, Hendrickx G, Zeller H, Van Bortel W: Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit Vectors. 2013, 6: 1-11. 10.1186/1756-3305-6-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Randolph SE, Storey K: Impact of microclimate on immature tick-rodent host interactions (Acari: Ixodidae): implications for parasite transmission. J Med Entomol. 1999, 36: 741-748.View ArticlePubMedGoogle Scholar
- Tagliapietra V, Rosa R, Arnoldi D, Cagnacci F, Capelli G, Montarsi F, Hauffe HC, Rizzoli A: Saturation deficit and deer density affect questing activity and local abundance of Ixodes ricinus (Acari, Ixodidae) in Italy. Vet Parasitol. 2011, 183: 114-124. 10.1016/j.vetpar.2011.07.022.View ArticlePubMedGoogle Scholar
- Li S, Heyman P, Cochez C, Simons L, Vanwambeke SO: A multi-level analysis of the relationship between environmental factors and questing Ixodes ricinus dynamics in Belgium. Parasit Vectors. 2012, 5: 149-10.1186/1756-3305-5-149.PubMed CentralView ArticlePubMedGoogle Scholar
- Gern L: Borrelia burgdorferi sensu lato, the agent of lyme borreliosis: life in the wilds. Parasite. 2008, 15: 244-247. 10.1051/parasite/2008153244.View ArticlePubMedGoogle Scholar
- Margos G, Vollmer SA, Ogden NH, Fish D: Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato. Infect Genet Evol. 2011, 11: 1545-1563. 10.1016/j.meegid.2011.07.022.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Dam AP, Kuiper H, Vos K, Widjojokusumo A, De Jongh BM, Spanjaard L, Ramselaar AC, Kramer MD, Dankert J: Different genospecies of Borrelia burgdorferi are associated with distinct clinical manifestations of Lyme borreliosis. Clin Infect Dis. 1993, 17: 708-717. 10.1093/clinids/17.4.708.View ArticlePubMedGoogle Scholar
- Bazovská S, Ďurovská J, Derdáková M, Tarageľová V, Pancák J, Záborská M, Traubner P: The genospecies Borrelia burgdorferi s.l., isolated from ticks and from neurological patients with suspected Lyme borreliosis. Neuro Endocrinol Lett. 2011, 32: 491-495.PubMedGoogle Scholar
- Hubálek Z: Epidemiology of Lyme Borreliosis. Curr Probl Dermatol. 2009, 37: 31-50.View ArticlePubMedGoogle Scholar
- Woldehivet Z: The natural history of Anaplasma phagocytophilum. Vet Parasitol. 2010, 167: 108-122. 10.1016/j.vetpar.2009.09.013.View ArticleGoogle Scholar
- Doudier B, Olano J, Parola P, Brouqui P: Factors contributing to emergence of Ehrlichia and Anaplasma spp. as human pathogens. Vet Parasitol. 2010, 167: 149-154. 10.1016/j.vetpar.2009.09.016.View ArticlePubMedGoogle Scholar
- Nováková M, Víchová B, Majláthová V, Lesňáková A, Pochybová M, Peťko B: First case of human granulocytic anaplasmosis from Slovakia. Ann Agric Environ Med. 2010, 17: 173-175.PubMedGoogle Scholar
- Derdáková M, Štefančíková A, Špitálska E, Tarageľová V, Košťálová T, Hrkľová G, Kybicová K, Schánilec P, Majláthová V, Várady M, Peťko B: Emergence and genetic variability of Anaplasma species in small ruminants and ticks from Central Europe. Vet Microbiol. 2011, 153: 293-298. 10.1016/j.vetmic.2011.05.044.View ArticlePubMedGoogle Scholar
- Schouls LM, Van de Pol I, Rijpkema SG, Schot CS: Detection and identification of Ehrlichia, Borrelia burgdorferi sensu lato, and Bartonella species in Dutch Ixodes ricinus ticks. J Clin Microbiol. 1999, 37: 2215-2222.PubMed CentralPubMedGoogle Scholar
- Špitálska E, Literak I, Sparagano OAE, Golovchenko M, Kocianova E: Ticks (Ixodidae) from passerine birds in the Carpathian region. Wien Klin Wochenschr. 2006, 118: 759-764. 10.1007/s00508-006-0729-4.View ArticlePubMedGoogle Scholar
- Lommano E, Bertaiola L, Dupasquier C, Gern L: Infections and coinfections of questing Ixodes ricinus ticks by emerging zoonotic pathogens in western Switzerland. App Environ Microbiol. 2012, 78: 4606-4612. 10.1128/AEM.07961-11.View ArticleGoogle Scholar
- Jahfari S, Fonville M, Hengeveld P, Reusken C, Scholte EJ, Takken W, Heyman P, Medlock P, Heylen D, Kleve J, Sprong H: Prevalence of Neoehrlichia mikurensis in ticks and rodents from North-west Europe. Parasit Vectors. 2012, 5: 74-10.1186/1756-3305-5-74.PubMed CentralView ArticlePubMedGoogle Scholar
- Kawahara M, Rikihisa Y, Isogai E, Takahashi M, Misumi H, Suto C, Shibata S, Zhang C, Tsuji M: Ultrastructure and phylogenetic analysis of ‘Candidatus Neoehrlichia mikurensis’ in the family Anaplasmataceae, isolated from wild rats and found in Ixodes ovatus ticks. Int J Syst Evol Microbiol. 2004, 54: 1837-1843. 10.1099/ijs.0.63260-0.View ArticlePubMedGoogle Scholar
- Vayssier-Taussat M, Le Rhun D, Buffet JP, Maaoui N, Galan M, Guivier E, Charbonnel N, Cosson JF: Candidatus Neoehrlichia mikurensis in bank voles, France. Emerg Infect Dis. 2012, 18: 2063-2065. 10.3201/eid1812.120846.PubMed CentralView ArticlePubMedGoogle Scholar
- Fehr JS, Bloemberg GV, Ritter C, Hombach M, Luscher TF, Weber R, Keller PM: Bacterial pathogen Candidatus Neoehrlichia mikurensis. Emerg Infect Dis. 2010, 16: 1127-1129. 10.3201/eid1607.091907.PubMed CentralView ArticlePubMedGoogle Scholar
- Welinder-Olsson C, Kjellin E, Vaht K, Jacobsson S, Wennerås C: First case of human Neoehrlichia mikurensis infection in a febrile patient with chronic lymphocytic leukemia. J Clin Microbiol. 2010, 48: 1956-1959. 10.1128/JCM.02423-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Von Loewenich FD, Geissdörfer W, Disqué C, Matten J, Schett G, Sakka SG, Bogdan C: Detection of Candidatus Neoehrlichia mikurensis in two patients with severe fibrile illnesses: evidence for a European sequence variant. J Clin Microbiol. 2010, 48: 2630-2635. 10.1128/JCM.00588-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Pekova S, Vydra J, Kabickova H: Candidatus Neoehrlichia mikurensis infection identified in 2 hematooncologic patients: benefit of molecular techniques for rare pathogen detection. Diag Microbiol Infect Dis. 2011, 69: 266-270. 10.1016/j.diagmicrobio.2010.10.004.View ArticleGoogle Scholar
- Diniz PP, Schulz BS, Hartmann K, Breitschwerdt EB: Candidatus Neoehrlichia mikurensis infection in a dog from Germany. J Clin Microbiol. 2011, 49: 2059-2062. 10.1128/JCM.02327-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Derdáková M, Halánová M, Stanko M, Štefančíková A, Čisláková L, Peť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.PubMedGoogle Scholar
- Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith G: Mixed effects models and extensions in ecology with R. 2009, New York: SpringerView ArticleGoogle Scholar
- Epidemiological Information System. http://www.epis.sk,
- R Core Team: A language and environment for statistical computing. http://www.R-project.org,
- Pinheiro J, Bates D, DebRoy S, Sarkar D: R package version 3.1-108. The R development core team nlme: linear and nonlinear mixed effect models. 2013, http://www.R-project.org,Google Scholar
- Guy EC, Stanek G: Detection of Borrelia burgdorferi in patients with Lyme disease by the polymerase chain-reaction. J Clin Pathol. 1991, 44: 610-611. 10.1136/jcp.44.7.610.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1128/AEM.69.1.509-516.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- De La Fuente J, Massung RF, Wong SJ, Chu FK, Lutz H, Meli M, Von Loewenich FD, Grzeszczuk A, Torina A, Caracappa S, Mangold AJ, Naranjo V, Stuen S, Kocan KM: Sequence analysis of the msp4 g ene of Anaplasma phagocytophilum s trains. J Clin Microbiol. 2005, 43: 1309-1317. 10.1128/JCM.43.3.1309-1317.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSIBLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Perez JL, Guigoz E, Rais O, Gern L: Influence of saturation deficit and temperature on Ixodes ricinus tick questing activity in a Lyme borreliosis-endemic area (Switzerland). Parasitol Res. 2000, 86: 554-557. 10.1007/s004360000209.View ArticleGoogle Scholar
- Knap K, Durmiši E, Saksida A, Korva M, Petrovec M, Avšič-Županc T: Influence of climatic factors on dynamics of questing Ixodes ricinus ticks in Slovenia. Vet Parasitol. 2009, 164: 275-281. 10.1016/j.vetpar.2009.06.001.View ArticlePubMedGoogle Scholar
- Egyed L, Élö P, Sreter-Lancz Z, Szell Z, Balogh Z, Sréter T: Seasonal activity and tick-borne pathogen infection rates of Ixodes ricinus ticks in Hungary. Ticks Tick Borne Dis. 2012, 3: 90-94. 10.1016/j.ttbdis.2012.01.002.View ArticlePubMedGoogle Scholar
- Randolph SE: Human activities predominate in determining changing incidence of tick-borne encephalitis in Europe. Euro Surveill. 2010, 15: 24-21.PubMedGoogle Scholar
- Stefanoff P, Rosinska M, Samuels S, White DJ, Morse DL, Randolph SE: A national case control study identifies human socio-economic status and activities as risk factors for tick-borne encephalitis in Poland. PLoS One. 2012, 7: e4551-View ArticleGoogle Scholar
- Cisak E, Wójcik-Fatla A, Stojek NM, Chmielewska-Badora J, Zwoliński J, Buczek A, Dutkiewicz J: Prevalence of Borrelia burgdorferi genospecies in Ixodes ricinus ticks from Lublin region (Eastern Poland). Ann Agric Environ Med. 2006, 13: 301-306.PubMedGoogle Scholar
- Lenčáková D, Hizo-Teufel C, Pet”ko 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, 40: 108-118.View ArticleGoogle Scholar
- Geller J, Nazarova L, Katargina O, Golovljova I: Borrelia burgdorferi sensu lato prevalence in tick populations in Estonia. Parasit Vectors. 2013, 6: 202-10.1186/1756-3305-6-202.PubMed CentralView ArticlePubMedGoogle Scholar
- Mysterud A, Easterday WR, Qviller L, Viljugrein H, Ytrehus B: Spatial and seasonal variation in the prevalence of Anaplasma phagocytophilum and Borrelia burgdorferi sensu lato in questing Ixodes ricinus ticks in Norway. Parasit Vectors. 2013, 6: 187-10.1186/1756-3305-6-187.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.