Next Article in Journal
Mindfulness Training plus Nature Exposure for Veterans with Psychiatric and Substance Use Disorders: A Model Intervention
Next Article in Special Issue
Knowledge and Knowledge Needs about Lyme Disease among Occupational and Recreational Users of the Outdoors
Previous Article in Journal
Acrylamide and Thermal-Processing Indexes in Market-Purchased Food
Previous Article in Special Issue
Prolongation of Tick-Borne Encephalitis Cycles in Warmer Climatic Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Presence of Roe Deer Affects the Occurrence of Anaplasma phagocytophilum Ecotypes in Questing Ixodes ricinus in Different Habitat Types of Central Europe

1
Institute of Zoology, Slovak Academy of Sciences, Dúbravská cesta 9, 84506 Bratislava, Slovakia
2
Institute of Infectology, Friedrich-Loeffler-Institute, Südufer 10, 17493 Greifswald-Insel Riems, Germany
3
Laboratory for Zoonoses and Environmental Microbiology, National Institute for Public Health and Environment, Antonie van Leeuwenhoeklaan 9, P.O. Box 1, 3720 BA Bilthoven, The Netherlands
4
Institute of Vertebrate Biology, The Czech Academy of Sciences, Klášterní 2, 69142 Valtice, Czech Republic
5
Institute of Macromolecular Chemistry, The Czech Academy of Sciences, Heyrovského nám. 1888/2, 16200 Prague 6, Czech Republic
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(23), 4725; https://doi.org/10.3390/ijerph16234725
Submission received: 29 October 2019 / Revised: 22 November 2019 / Accepted: 25 November 2019 / Published: 27 November 2019
(This article belongs to the Special Issue Ticks and Tick Vectored Diseases—Biology to Society)

Abstract

:
The way in which European genetic variants of Anaplasma phagocytophilum circulate in their natural foci and which variants cause disease in humans or livestock remains thus far unclear. Red deer and roe deer are suggested to be reservoirs for some European A. phagocytophilum strains, and Ixodes ricinus is their principal vector. Based on groEL gene sequences, five A. phagocytophilum ecotypes have been identified. Ecotype I is associated with the broadest host range, including strains that cause disease in domestic animals and humans. Ecotype II is associated with roe deer and does not include zoonotic strains. In the present study, questing I. ricinus were collected in urban, pasture, and natural habitats in the Czech Republic, Germany, and Slovakia. A fragment of the msp2 gene of A. phagocytophilum was amplified by real-time PCR in DNA isolated from ticks. Positive samples were further analyzed by nested PCRs targeting fragments of the 16S rRNA and groEL genes, followed by sequencing. Samples were stratified according to the presence/absence of roe deer at the sampling sites. Geographic origin, habitat, and tick stage were also considered. The probability that A. phagocytophilum is a particular ecotype was estimated by a generalized linear model. Anaplasma phagocytophilum was identified by genetic typing in 274 I. ricinus ticks. The majority belonged to ecotype I (63.9%), 28.5% were ecotype II, and both ecotypes were identified in 7.7% of ticks. Ecotype II was more frequently identified in ticks originating from a site with presence of roe deer, whereas ecotype I was more frequent in adult ticks than in nymphs. Models taking into account the country-specific, site-specific, and habitat-specific aspects did not improve the goodness of the fit. Thus, roe deer presence in a certain site and the tick developmental stage are suggested to be the two factors consistently influencing the occurrence of a particular A. phagocytophilum ecotype in a positive I. ricinus tick.

1. Introduction

Anaplasma phagocytophilum is a tick-borne obligate intracellular bacterium causing anaplasmosis in humans and animals. Its ecology and epidemiology are very complex in Europe due to the occurrence of several ecotypes that circulate in distinct transmission cycles [1,2,3,4,5]. By using molecular markers, the presence of diverse A. phagocytophilum genetic variants has been reported in a wide range of free-living and domestic animals [6,7,8,9,10]. Cervids can be reservoirs for several A. phagocytophilum genetic variants transmitted by Ixodes ricinus. Certain variants associated with roe deer (Capreolus capreolus) are probably nonpathogenic to humans, domestic ruminants, horses, and dogs, while red deer (Cervus elaphus) are likely to be reservoirs for variants pathogenic to humans, domestic ruminants, and horses [3,11,12,13,14]. However, co-infection of roe deer with two–three distinct genetic variants, including those causing disease in domestic ruminants, has also been reported recently [15]. The potential reservoir role of wild boar (Sus scrofa) for the human granulocytic agent has been suggested [16,17], but generally the role of free-living animals in circulation of A. phagocytophilum strains of medical and veterinary importance is still incompletely understood.
Based on groEL gene sequences derived from vertebrates and ticks, four or five A. phagocytophilum ecotypes have been identified, each associated with different host ranges in Europe [3,18]. The majority of the analyzed sequences belonged to ecotype I and II. Members of ecotype I infect the broadest host range, including wild-living ungulates and I. ricinus ticks, and comprise also strains causing disease in domestic animals and humans. Ecotype II was found to be associated with roe deer and I. ricinus ticks and does not appear to include zoonotic strains. Ecotype III was identified only in rodents and the rodent-associated Ixodes trianguliceps ticks, and ecotype IV was found to be associated with birds and ticks feeding on birds (mainly Ixodes frontalis, but also I. ricinus) [3]. Ecotype V was identified in Ixodes ventalloi ticks from southern Europe, but its pathogenic potential is thus far unknown [18]. Although classification of ecotypes based on the groEL gene sequences was found to be more clear than classification based on the 16S rRNA gene sequences, still more research is required to elucidate the epidemiological significance of European A. phagocytophilum genetic linages and ecotypes [4,5,10].
The main objective of the present work was to investigate which factor(s) out of geographic location, habitat type, presence/absence of roe deer, or tick developmental stage could determine the presence of particular A. phagocytophilum ecotypes in questing I. ricinus ticks.

2. Materials and Methods

2.1. Study Sites

Questing I. ricinus ticks were collected in different habitats—urban, pasture, and natural, with regard to the occurrence of various potential reservoir hosts for A. phagocytophilum—in three countries of Central Europe: Czech Republic, Germany, and Slovakia (Table 1).
Six study areas were selected in Bavaria, southern Germany (see [19,20,21] for detailed descriptions of sites):
A1. 
The southern part of Englischer Garten is a city center park in Munich under strong anthropogenic influence, highly visited by citizens and dogs, maintained by gardening. Interaction between domestic and free-living animals is limited to mice, hedgehogs, foxes, and birds.
A2. 
Nymphenburger Schlosspark is a city park in Munich, surrounded by walls. Several forest-like structures with higher tree density, scrub, and free-living animals such as roe deer are present in this park. Grassy areas are extensively used and less frequently mowed.
A3. 
Dörnbergpark is situated in the city center of Regensburg. It is a small park, surrounded by walls, with strong anthropogenic influence, which is expressed by a high frequency of visitors spending their leisure time there. The site is a well-tended park with mostly grassy landscape and only a few old trees such as oaks and maples. Large free-living mammals like roe deer and wild boar are not present.
A4. 
Schlosspark Berg is a forest-like park on the eastern shore of Lake Starnberg with scrub and walkways. Deciduous trees and bushes, large free-living mammals such as roe deer and foxes exist and hunting is practiced.
A5. 
Kerschlach is an agricultural site formed by two fenced pastures, which are used for stock breeding of cattle, and is surrounded by a forest area. Contact between wild animals (also roe deer) and domestic cattle exists because animals from the forest can enter the pasture.
A6. 
Tussenhausen (Angelberger Forst) is a large mixed forest with low anthropogenic influence. This forest is mainly dominated by beech, oak, and spruce. Different free-living animal species are present, and the frequency of visitors is low.
Study sites in the Czech Republic were located in five areas and were characterized in detail in [22,23]:
B1. 
Valtice is an enclosed, fenced urban park in the castle grounds, where vegetation is periodically adjusted. The fauna is represented by birds, small mammals, and roe deer, the vegetation by broad-leaved deciduous trees and mown grass.
B2. 
Pohansko is a natural floodplain forest. The fauna consists of birds, various small and large mammals such as roe, red and fallow deer, and wild boar, the vegetation of deciduous forests.
B3. 
Suchovské mlýny is a pasture area; ticks were sampled along a sheep fence. The fauna is represented by birds, small mammals, and sheep, the vegetation by some trees and shrubs.
B4. 
Proskovice is a natural ecosystem outside the town of Ostrava. This mixed forest with dominant broad-leaved trees is rarely visited by people. The fauna consists of small and medium-sized mammals, roe deer, birds, and occasionally wild boar.
B5. 
Bělský les is an urban park in Ostrava. The local fauna is represented by birds and small mammals, and the vegetation by broad-leaved deciduous trees and grass. The park is surrounded by housing estates and used for leisure activities and dog-walking.
In Slovakia, four study sites were selected (for more details see [24,25]):
C1. 
The campus of the Slovak Academy of Sciences (SAS) is an enclosed urban area, within the built-up plots there are patches of broad-leaved forest vegetation, shrubby plots, mown and unmown grassy plots. The fauna consists mainly of birds, feral cats, hedgehogs, and roe deer. The density of rodents is very low.
C2. 
Mokrohájska street is an urban area situated along fences of gardens and cottages in the vicinity of the SAS campus. The fauna is represented by birds, small mammals, and lizards.
C3. 
Železná studienka is part of Bratislava Forest Park, located on the foothills of the Small Carpathians Mountains in the northern part of Bratislava. The vegetation is represented by broad-leaved trees, the fauna by birds, small and medium-sized mammals, deer, and wild boar. Železná studienka serves for relaxation of inhabitants of Bratislava and dog walking.
C4. 
Fúgelka is a natural forest site visited by tourists, rangers, and hunters. The vegetation is represented by a mixed forest, the fauna by various small, medium-sized, and large free-living mammals such as roe, red and fallow deer, mouflon, and wild boar.

2.2. Tick Sampling, DNA Extraction, and PCR Amplification

Questing I. ricinus ticks were collected from vegetation during previous studies carried out in Bavaria (Germany) in 2011–2012 [20], in eastern and south-eastern Czech Republic in 2010–2014 [22,23], and in south-western Slovakia in 2011–2013 [24].
In each laboratory, DNA was extracted from individual nymphal and adult ticks by commercial isolation kits. A 77-bp long fragment of the msp2 gene of A. phagocytophilum was amplified by real-time PCR, according to Courtney et al. [26] and as described in Overzier et al. [20], Venclikova et al. [22,23], and Svitálková et al. [24], for ticks from Germany, the Czech Republic, and Slovakia, respectively. To identify A. phagocytophilum variants, positive samples with Ct values ≤30 were further analyzed with a nested PCR targeting a 497-bp fragment of the 16S rRNA gene [27], as described in Overzier et al. [20], and with primers targeting a fragment of the groEL gene, according to Alberti et al. [28]. Primers and probes used in the PCR reactions are listed in Table 2. All sequences were confirmed by sequencing of both strands.

2.3. Sequence Analysis

16S rRNA gene sequences were analyzed by ChromasLite® (Technelysium Pty Ltd., South Brisbane, Australia) and by MEGA 6 [29] and aligned with MUSCLE [30]. Database searches and sequence comparisons were performed with the BLAST tool provided by the National Centre for Biotechnology (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The 16S rRNA gene variants were designated according to Schorn et al. [31] and Silaghi et al. [32]. GroEL gene sequences were analyzed in Bionumerics (Version 7.5, Applied Math, Sint-Martens-Latem, Belgium), after subtraction of the primer sequences.

2.4. Statistical Analyses

Anaplasma phagocytophilum ecotypes were assigned to each sequence, as described previously [3]. Samples were stratified according to the presence and absence of roe deer at the sampling locations. Only those samples for which amplicons could be obtained for both the partial 16S rRNA and the groEL gene were included in the analysis. Next to this, other factors such as geographic origin (country), landscape, and the different I. ricinus tick developmental stages were taken into account (Table 3). Each A. phagocytophilum ecotype identified in the I. ricinus ticks was recorded once into a database for this analysis. As genetic typing was performed at the 16S and groEL loci, sometimes both ecotype I and II were identified in a single I. ricinus tick. We consider these ticks to be co-infected with two different A. phagocytophilum ecotypes.
We estimated the probability that A. phagocytophilum is a particular ecotype by fitting a generalized linear model (glm) to the resulting database, assuming binomial distribution and a logit link. The best-fit model was identified by running the command “step” in forward and backward directions. Collinearity was identified by calculating the condition number of the model matrix. All computations were performed using R [33].
A chi-square (χ2) test of independence was performed by using PAST [34] to examine the relationship between presence/absence of roe deer and the occurrence of ticks infected with particular 16S rRNA gene variants.

3. Results

Anaplasma phagocytophilum ecotypes were identified by genetic typing in a total of 274 infected I. ricinus ticks that were collected in urban, suburban, park, pasture, and natural landscapes at 13 geographically distinct sites in three Central European countries (Table 1). In the majority of those samples, only one ecotype was identified: ecotype I (n = 175; 63.9%); ecotype II (n = 78; 28.5%). Both ecotype I and ecotype II were identified in 21 (7.7%) I. ricinus samples (Table 3; Supplementary file 1: Table S1).
A model taking into account roe deer presence in the sampling sites and nymphal and adult stages of ticks fitted best to the dataset, as judged by the minimal value of Bayesian Information Criterion. Ecotype II was more frequently identified in an I. ricinus tick when it was collected in a territory with presence of roe deer (p < 0.001). Ecotype I was more frequently identified in adult I. ricinus ticks compared with nymphs (p = 0.004). Alternative models taking additionally into account the country-specific, site-specific, and landscape-specific aspects in our dataset did not improve the goodness of the fit.
Akaike Information Criterion (AIC) is an alternative measure for the goodness of fit. Based on AIC, roe deer and tick’s developmental stage remained significant factors, but sampling site emerged as a new significant factor. None of the individual estimates for a specific site were significant. This led to an extra assessment of the new factor. Roe deer, tick stage, and site were collinear, as detected by the condition number (=62) of the model matrix. High collinearity and insignificant parameter estimates for individual sites led us to conclude that roe deer presence and tick’s developmental stage are the only two factors consistently influencing the occurrence of a particular ecotype in an A. phagocytophilum-positive I. ricinus tick.
By analyzing 16S rRNA gene sequences obtained from 291 ticks (the number of ticks for which amplification of the 16S rRNA fragment was successful was higher than for those for which both the 16S rRNA and the groEL fragments could be amplified), 13 A. phagocytophilum variants were identified. Their occurrence and frequency varied between the study sites (Table 4). Variant A was the most frequent (59.1%). The overall proportion of ticks infected with variant A was significantly related to the presence of roe deer (χ2 = 64.146, p < 0.001) and was higher in sites with absence of these ungulates (90.3% vs. 42.0% in sites with roe deer) (Table 4). This trend was similar in both nymphs and adult ticks (nymphs: χ2 = 8.579, p < 0.01; adults: χ2 = 54.307, p < 0.001). In sites where roe deer were present, ticks carrying variant X (total: χ2 = 21.228, p < 0.001; nymphs: χ2 = 8.811, p < 0.01; adults: χ2 = 10.896, p < 0.001) and Y (total: χ2 = 11.155, p < 0.01; nymphs: χ2 = 0.617, ns; adults: χ2 = 11.379, p < 0.001) occurred more frequently than in sites without deer (Table 4).

4. Discussion

Ixodes ricinus is a generalist three-host tick that is able to parasitize a wide range of vertebrate hosts and transmit a great variability of microorganisms, including A. phagocytophilum [35]. The ecology of A. phagocytophilum in Europe is very complex and involves specific host–microbe–tick associations and distinct transmission cycles in which zoonotic and nonzoonotic strains circulate [2]. In spite of existing knowledge on associations of particular European A. phagocytophilum genetic lineages to certain vertebrate hosts and vector ticks and on their pathogenicity to humans and domestic animals, the underlying mechanisms determining the genetic variability of the bacterium are still not clear. Ixodes ricinus, due to its complex life cycle and wide host range, but mainly vertebrate reservoirs, due to peculiarities of their immune system, seem to drive the evolution and distribution of A. phagocytophilum genetic lineages [5].
It has already been confirmed that roe deer are important feeding hosts for all developmental stages of I. ricinus and reservoirs for particular nonzoonotic genetic variants of A. phagocytophilum [3,12,32,36,37], whereby the overall prevalence of infection in roe deer in particular sites can be higher than 90% [36,37,38,39]. The main objective of this study was to find out if, and to what extent, presence/absence of roe deer in a particular site can determine the occurrence of A. phagocytophilum strains in questing I. ricinus ticks. For this purpose, 13 sites located in three European countries and different landscapes were selected for tick collection, and the A. phagocytophilum strains in ticks were characterized based on partial 16S rRNA and groEL gene sequences. Anaplasma phagocytophilum was previously detected by PCR in questing I. ricinus ticks from each of the study sites, with total prevalence ranging between 1.4% and 19.7% (see Table 5). However, presence/absence of roe deer did not seem to affect the overall infection rates, except for an urban fenced area in Bratislava (Slovakia) with a relatively high roe deer density and a 19.7% A. phagocytophilum prevalence (Table 5 and [24]). Because transovarial transmission has not been confirmed in ticks, different reservoir hosts are likely to maintain the circulation of the bacterium in sites without roe deer presence. It was previously indicated for particular regions that the distribution of A. phagocytophilum ecotypes in ticks was linked to the presence of certain cervid species [20,36,37,40]. The prevalence of the roe deer-associated ecotype II as well as of the 16S rRNA variants X and Y (ranked among ecotype II, see Supplementary file 1: Table S1) in questing I. ricinus nymphs and adults was indeed higher in our study sites with roe deer presence than in locations where roe deer were absent. However, infection rates in adult ticks and nymphs differed in those sites: ecotype II was detected in 54.2% of nymphs and 38.0% of adults, 16S rRNA variant X occurred in 34.6% of nymphs and 14.7% of adults, and 16S rRNA variant Y was detected in 11.5% of nymphs and 12.5% of adults. These differences are difficult to explain, but could be, for example, due to nymph mortality, dispersal to other locations, or diminishing of the infection during molting to the adult stage. Transstadial transmission of some A. phagocytophilum strains has been confirmed experimentally [41,42], but transmission rates for the roe deer-associated strains are unknown. It is also not clear if and to what extent the infection with a particular A. phagocytophilum strain perpetuates in adult ticks in case infected nymphs feed on uninfected or noncompetent hosts or feed on hosts infected with other A. phagocytophilum strains or other pathogen(s). For example, in an agricultural site in France with high cattle abundance and about 80% A. phagocytophilum prevalence in roe deer, cattle density was suggested to be the main factor affecting the bacterial prevalence and occurrence of genotypes in questing I. ricinus, because almost 80% of the examined ticks carried sequences associated with infections in cattle and humans, and the high infection prevalence in roe deer was probably due to superinfections [43].
Ticks harboring A. phagocytophilum ecotype I that comprises strains associated with a wide range of different mammal species including red deer, small free-living and domestic ruminants, wild boar, carnivores, hedgehogs, and humans, were found to prevail in both types of sites, that is, with presence/absence of roe deer, but the infection rates were higher in sites without roe deer. This suggests that in sites with presence of roe deer the transmission hosts are not limited only to these animals but include a wider spectrum of vertebrate species, whereas in sites without roe deer, the bacterium is maintained by other reservoir hosts. The occurrence of ticks infected with ecotype II in sites where roe deer were absent may be the result of transport of those ticks by birds or mammals with large home ranges.
In total, about 7% of ticks (nymphs as well as adults) from the two types of sites harbored both ecotype I and II, suggesting that either the ticks fed on hosts co-infected with different A. phagocytophilum strains or larvae and nymphs fed successively on hosts harboring various strains. Simultaneous infection of vertebrates with multiple A. phagocytophilum strains was found to exist and has recently been confirmed, for example, for roe deer from France that harbored 2–3 ankA gene variants [15] or cattle from Germany infected with two ankA variants and multiple MLST sequence types [44]. Co-infections may lead to bacterial recombination and emergence of new genetic variants [44,45]. However, it is not clear if multiple strains are co-transmitted by I. ricinus.
The rodent-associated ecotype III and the bird-associated ecotype IV were not found in this study, although presence of rodents and birds was observed in the majority of the study sites. The main reasons why these ecotypes were not detected are probably because (1) the representation of ecotypes III and IV is generally much lower than that of ecotypes I and II [3], and (2) only questing I. ricinus were analyzed, while the vector of ecotype III is the endophilic I. trianguliceps tick and ecotype IV can be found mainly in birds and bird-feeding I. frontalis and I. ricinus ticks.
Analysis of the variability of the partial 16S rRNA gene revealed the presence of 13 A. phagocytophilum variants in I. ricinus: nine, six, and ten from sites in Germany, the Czech Republic, and Slovakia, respectively (see Table 4). In addition to the roe deer-associated variants X and Y (discussed above), the occurrence and prevalence of the other variants also differed in the study sites, whereby their spectrum was generally wider in natural sites with occurrence of roe deer and numerous free-living vertebrate host species than in urban areas and pastures. Although a basic distinction of variants associated with free-living ruminants, domestic animals, or both groups can be achieved with 16S rRNA gene-based typing [32,37], this approach is not reliable to define A. phagocytophilum genotypes [12,14,36,44]. The occurrence of 16S rRNA variants in questing I. ricinus ticks and roe deer, their possible pathogenicity to domestic ruminants and humans, and their association with free-living reservoir hosts in the study sites in Germany have extensively been discussed in earlier publications [20,31,37]. However, the results obtained for sites in the Czech Republic and Slovakia are new. Based on analysis of a limited number of free-living ungulates and engorged ticks from the Small Carpathian Mountains in south-western Slovakia, where the forest park and natural tick collection site are located, variants Y and X (ecotype II) were previously detected in roe deer and roe deer attached ticks, respectively. Ecotype I, comprising variant B, was detected in wild boar and fallow deer attached ticks, variant S in fallow deer, variant S and W in red deer and fallow deer attached ticks, and variant W in mouflon [39]. These variants were also detected in questing ticks in our study and in addition, variants A, I, J, N, and V were identified, whereby variants I, J. N, V and some sequences of variant W were ranked among ecotype II (see Supplementary file 1: Table S1). Both variant A (ecotype I) and V (ecotype II) were previously found associated with, for example, questing I. ricinus ticks and hedgehogs, while variant A was also associated with red fox, cat, dog, horse, and a human case [36], and variant V was detected in roe deer, fallow deer, wild boar, rodents, and small ruminants [20,31,46,47,48]. Similar to Germany, variant A prevailed in the Slovak urban site without roe deer. Interestingly and in contrast to Germany, a relatively high proportion of ticks carrying variant A was found in Slovak suburban and natural sites with roe deer presence (see Table 4). This suggests that the host range of this variant may be wider than previously thought, and the spectrum and abundance of potential reservoir hosts differ between the sites in the two countries. Variant B, which is identical with the prototype human pathogenic variant of A. phagocytophilum, was detected only in one adult tick from the natural site, but was previously identified in wild boar from the same area [39] and also in horses and dogs with granulocytic anaplasmosis in Germany [47,48]. Variants I and J were identified, for example, in roe deer and red deer in Germany and Austria, respectively [32,37], and variant N in roe deer from Germany [37]. Variant S was found, for example, in red deer, roe deer, and ibex in Austria [32] and in red deer and fallow deer in the Slovak study area [39]. Variant W was found to be associated with clinical cases of tick-borne fewer in domestic ruminants, but was identified also, for example, in roe deer, red deer, chamois, and mouflon [32,37,39].
To our knowledge, there are no data on occurrence of 16S rRNA variants or ecotypes based on groEL sequences in vertebrate hosts from the study sites in the Czech Republic. Similar to Germany and Slovakia, I. ricinus ticks with variant A prevailed in the urban park (Bělský les) without roe deer. In the agricultural site Suchovské mlýny, ticks with variant W dominated, suggesting that they probably fed on infected sheep grazing in that area. Interestingly and in contrast to German and Slovak sites, ticks harboring variant A dominated in sites with roe deer presence, and only 16.7% of ticks carried roe deer-associated variants X and Y. Presence of variant B was confirmed in ticks from the urban park in Valtice and in the natural site Pohansko. The occurrence of variants B and W suggests that I. ricinus is also involved in the circulation of A. phagocytophilum variants pathogenic for domestic animals and humans in the Czech Republic.

5. Conclusions

In conclusion, the occurrence and prevalence of A. phagocytophilum ecotypes and partial 16S rRNA gene variants in questing I. ricinus ticks from different habitat types of Germany, the Czech Republic, and Slovakia that were investigated during our study depended mainly on the presence/absence of roe deer. In sites with presence of roe deer, the prevalence of ticks infected with nonpathogenic variants of the bacterium was significantly higher than in sites where roe deer were absent. The results suggest that presence of roe deer might, at least to a certain degree, decrease the risk of contracting granulocytic anaplasmosis in domestic animals and humans. However, to support this assumption, further research is needed to elucidate the evolution of A. phagocytophilum strains and their associations to vertebrate hosts.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/16/23/4725/s1, Table S1 List of analyzed I. ricinus ticks for 16S rRNA and GroEl genetic variants, including tick developmental stage (N: nymph, F: female, M: male), country of origin, area/site, type of habitat, sampling date, indication of presence/absence of roe deer, and the analyzed GroEL and 16S rRNA sequences.

Author Contributions

Conceptualization, study design: C.S., H.S., Z.H., M.K.; molecular analyses, primary data analysis: Z.H., K.G., I.R., C.S., H.S.; data analyses, statistical evaluation of the results: K.T., H.S., M.K.; drafting the manuscript: M.K., Z.H., K.T.; manuscript corrections and revisions: C.S., H.S., I.R.

Funding

The study 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. The work of Z.H., C.S., and M.K. was supported by the German Academic Exchange Service (DAAD), and the work of Z.H. was additionally supported by EurNegVec COST Action TD1303.

Acknowledgments

The authors would like to thank Ryanne Jaarsma (RIVM, Bilthoven, The Netherlands) and Claudia Thiel (Institute of Comparative Tropical Medicine and Parasitology, Faculty of Veterinary Medicine, LMU, Munich, Germany) for technical assistance in the laboratory.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Woldehiwet, Z. The natural history of Anaplasma phagocytophilum. Vet. Parasitol. 2010, 167, 108–122. [Google Scholar] [CrossRef] [PubMed]
  2. Stuen, S.; Granquist, E.G.; Silaghi, C. Anaplasma phagocytophilum—A widespread multihost pathogen with highly adaptive strategies. Front. Cell. Infect. Microbiol. 2013, 3, 31. [Google Scholar] [CrossRef] [PubMed]
  3. Jahfari, S.; Coipan, C.; Fonville, M.; van Leeuwen, A.D.; Hengeveld, P.; Heylen, D.; Heyman, P.; van Maanen, C.; Butler, C.M.; Földvári, G.; et al. Circulation of four Anaplasma phagocytophilum ecotypes in Europe. Parasit. Vectors 2014, 7, 365. [Google Scholar] [CrossRef] [PubMed]
  4. Dugat, T.; Lagrée, A.-C.; Maillard, R.; Boulouis, H.-J.; Haddad, N. Opening the black box of Anaplasma phagocytophilum diversity: Current situation and future perspectives. Front. Cell. Infect. Microbiol. 2015, 5, 61. [Google Scholar] [CrossRef]
  5. Jaarsma, R.I.; Sprong, H.; Takumi, K.; Kazimirova, M.; Silaghi, C.; Mysterud, A.; Rudolf, I.; Beck, R.; Földvári, G.; Tomassone, L.; et al. Anaplasma phagocytophilum evolves in geographical and biotic niches of vertebrates and ticks. Parasit. Vectors 2019, 12, 328. [Google Scholar] [CrossRef]
  6. Stuen, S. Anaplasma phagocytophilum - the most widespread tick-borne infection in animals in Europe. Vet. Res. Commun. 2007, 31 (Suppl. 1), 79–84. [Google Scholar] [CrossRef]
  7. Stuen, S.; Pettersen, K.S.; Granquist, E.G.; Bergström, K.; Bown, K.J.; Birtles, R.J. Anaplasma phagocytophilum variants in sympatric red deer (Cervus elaphus) and sheep in southern Norway. Ticks Tick Borne Dis. 2013, 4, 197–201. [Google Scholar] [CrossRef]
  8. Víchová, B.; Majláthová, V.; Nováková, M.; Stanko, M.; Hviščová, I.; Pangrácová, L.; Chrudimský, T.; Čurlík, J.; Peťko, B. Anaplasma infections in ticks and reservoir host in Slovakia. Infect. Genet. Evol. 2014, 22, 265–272. [Google Scholar] [CrossRef]
  9. Baráková, I.; Derdáková, M.; Carpi, G.; Rosso, F.; Collini, M.; Tagliapietra, V.; Ramponi, C.; Hauffe, H.-C.; Rizzoli, A. Genetic and ecologic variability among Anaplasma phagocytophilum strains in Northern Italy. Emerg. Infect. Dis. 2014, 20, 1082–1084. [Google Scholar] [CrossRef]
  10. Battilani, M.; De Arcangeli, S.; Balboni, A.; Dondi, F. Genetic diversity and molecular epidemiology of Anaplasma. Infect. Genet. Evol. 2017, 49, 195–211. [Google Scholar] [CrossRef]
  11. Petrovec, M.; Bidovec, A.; Sumner, J.W.; Nicholson, W.L.; Childs, J.E.; Avsic-Zupanc, T. Infection with Anaplasma phagocytophila in cervids from Slovenia: Evidence of two genotypic lineages. Wien. Klin. Wochenschr. 2002, 114, 641–647. [Google Scholar] [PubMed]
  12. Huhn, C.; Winter, C.; Wolfsperger, T.; Wüppenhorst, N.; Strašek Smrdel, K.; Skuballa, J.; Pfäffle, M.; Petney, T.; Silaghi, C.; Dyachenko, V.; et al. Analysis of the population structure of Anaplasma phagocytophilum using multilocus sequence typing. PLoS ONE 2014, 9, e93725. [Google Scholar] [CrossRef] [PubMed]
  13. Dugat, T.; Chastagner, A.; Lagrée, A.-C.; Petit, E.; Durand, B.; Thierry, S.; Corbière, F.; Verheyden, H.; Chabanne, L.; Bailly, X.; et al. A new multiple-locus variable-number tandem repeat analysis reveals different clusters for Anaplasma phagocytophilum circulating in domestic and wild ruminants. Parasit. Vectors 2014, 7, 439. [Google Scholar] [CrossRef] [PubMed]
  14. Langenwalder, D.B.; Schmidt, S.; Gilli, U.; Pantche, N.; Ganter, M.; Silaghi, C.; Aardema, M.L.; von Loewenich, F.D. Genetic characterization of Anaplasma phagocytophilum strains from goats (Capra aegagrus hircus) and water buffalo (Bubalus bubalis) by 16S rRNA gene, ankA gene and multilocus sequence typing. Ticks Tick Borne Dis. 2019, 10, 101267. [Google Scholar] [CrossRef] [PubMed]
  15. Jouglin, M.; Chagneau, S.; Faille, F.; Verheyden, H.; Bastian, S.; Malandrin, L. Detecting and characterizing mixed infections with genetic variants of Anaplasma phagocytophilum in roe deer (Capreolus capreolus) by developing an ankA cluster-specific nested PCR. Parasit. Vectors 2017, 10, 377. [Google Scholar] [CrossRef] [PubMed]
  16. Strasek Smrdel, K.; Bidovec, A.; Malovrh, T.; Petrovec, M.; Duh, D.; Avsic Zupanc, T. Detection of Anaplasma phagocytophilum in wild boar in Slovenia. Clin. Microbiol. Infect. 2009, 15 (Suppl. 2), 50–52. [Google Scholar] [CrossRef]
  17. Michalik, J.; Stańczak, J.; Cieniuch, S.; Racewicz, M.; Sikora, B.; Dabert, M. Wild boars as hosts of human-pathogenic Anaplasma phagocytophilum variants. Emerg. Infect. Dis. 2012, 18, 2094–2095. [Google Scholar] [CrossRef]
  18. Santos, A.S.; de Bruin, A.; Veloso, A.R.; Marques, C.; da Fonseca, I.P.; de Sousa, R.; Sprong, H.; Santos-Silva, M.M. Detection of Anaplasma phagocytophilum, Candidatus Neoehrlichia sp., Coxiella burnetii and Rickettsia spp. in questing ticks from a recreational park, Portugal. Ticks Tick Borne Dis. 2018, 9, 1555–1564. [Google Scholar] [CrossRef]
  19. Silaghi, C.; Woll, D.; Hamel, D.; Pfister, K.; Mahling, M.; Pfeffer, M. Babesia spp. and Anaplasma phagocytophilum in questing ticks, ticks parasitizing rodents and the parasitized rodents – Analyzing the host pathogen-vector interface in a metropolitan area. Parasit. Vectors 2012, 5, 191. [Google Scholar] [CrossRef]
  20. Overzier, E.; Pfister, K.; Thiel, C.; Herb, I.; Mahling, M.; Silaghi, C. Anaplasma phagocytophilum in questing Ixodes ricinus ticks: A comparison of the prevalence and partial 16S rRNA gene variants in urban, pasture and natural habitats. Appl. Environ. Microbiol. 2013, 79, 1730–1734. [Google Scholar] [CrossRef]
  21. Overzier, E.; Pfister, K.; Thiel, C.; Herb, I.; Mahling, M.; Silaghi, C. Diversity of Babesia and Rickettsia species in questing Ixodes ricinus: A longitudinal study in urban, pasture and natural habitats. Vector Borne Zoonot. Dis. 2013, 13, 559–564. [Google Scholar] [CrossRef] [PubMed]
  22. Venclikova, K.; Rudolf, I.; Mendel, J.; Betasova, L.; Hubalek, Z. Rickettsiae in questing Ixodes ricinus ticks in the Czech Republic. Ticks Tick Borne Dis. 2014, 5, 135–138. [Google Scholar] [CrossRef] [PubMed]
  23. Venclíková, K.; Mendel, J.; Betášová, L.; Blažejová, H.; Jedličková, P.; Straková, P.; Hubálek, Z.; Rudolf, I. Neglected tick-borne pathogens in the Czech Republic, 2011–2014. Ticks Tick Borne Dis. 2016, 7, 107–112. [Google Scholar] [CrossRef] [PubMed]
  24. 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]
  25. 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] [PubMed]
  26. Courtney, J.W.; Kostelnik, L.M.; Zeidner, N.S.; Massung, R.F. Multiplex real-time PCR for detection of Anaplasma phagocytophilum and Borrelia burgdorferi. J. Clin. Microbiol. 2004, 42, 3164–3168. [Google Scholar] [CrossRef]
  27. Massung, R.F.; Slater, K.; Owens, J.H.; Nicholson, W.L.; Mather, T.N.; Solberg, V.B.; Olson, J.G. Nested PCR assay for detection of granulocytic ehrlichiae. J. Clin. Microbiol. 1998, 36, 1090–1095. [Google Scholar]
  28. Alberti, A.; Zobba, R.; Chessa, B.; Addis, M.F.; Sparagano, O.; Pinna Parpaglia, M.L.; Cubeddu, T.; Pintori, G.; Pittau, M. Equine and canine Anaplasma phagocytophilum strains isolated on the island of Sardinia (Italy) are phylogenetically related to pathogenic strains from the United States. Appl. Environ. Microbiol. 2005, 71, 6418–6422. [Google Scholar] [CrossRef]
  29. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  30. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef]
  31. Schorn, S.; Pfister, K.; Reulen, H.; Mahling, M.; Manitz, J.; Thiel, C.; Silaghi, C. Prevalence of Anaplasma phagocytophilum in Ixodes ricinus in Bavarian public parks, Germany. Ticks Tick Borne Dis. 2011, 2, 196–203. [Google Scholar] [CrossRef] [PubMed]
  32. Silaghi, C.; Hamel, D.; Thiel, C.; Pfister, K.; Passos, L.M.; Rehbein, S. Genetic variants of Anaplasma phagocytophilum in wild caprine and cervid ungulates from the Alps in Tyrol, Austria. Vector Borne Zoonot. Dis. 2011, 11, 355–362. [Google Scholar] [CrossRef] [PubMed]
  33. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  34. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 9. [Google Scholar]
  35. 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. Public Health. 2014, 2, 251. [Google Scholar] [CrossRef]
  36. Scharf, W.; Schauer, S.; Freyburger, F.; Petrovec, M.; Schaarschmidt-Kiener, D.; Liebisch, G.; Runge, M.; Ganter, M.; Kehl, A.; Dumler, J.S.; et al. Distinct host species correlate with Anaplasma phagocytophilum ankA gene clusters. J. Clin. Microbiol. 2011, 49, 790–796. [Google Scholar] [CrossRef] [Green Version]
  37. Overzier, E.; Pfister, K.; Herb, I.; Mahling, M.; Bock, G. Jr.; Silaghi, C. Detection of tick-borne pathogens in roe deer (Capreolus capreolus), in questing ticks (Ixodes ricinus), and in ticks infesting roe deer in southern Germany. Ticks Tick Borne Dis. 2013, 4, 320–328. [Google Scholar] [CrossRef]
  38. Kauffmann, M.; Rehbein, S.; Hamel, D.; Lutz, W.; Heddergott, M.; Pfister, K.; Silaghi, C. Anaplasma phagocytophilum and Babesia spp. in roe deer (Capreolus capreolus), fallow deer (Dama dama) and mouflon (Ovis musimon) in Germany. Mol. Cell. Probes 2017, 31, 46–54. [Google Scholar] [CrossRef]
  39. Kazimírová, M.; Hamšíková, Z.; Špitalská, E.; Minichová, L.; Mahríková, L.; Caban, R.; Sprong, H.; Fonville, M.; Schnittger, L.; Kocianová, E. Diverse tick-borne microorganisms identified in free-living ungulates in Slovakia. Parasit. Vectors 2018, 11, 495. [Google Scholar] [CrossRef]
  40. Stigum, V.M.; Jaarsma, R.I.; Sprong, H.; Rolandsen, C.M.; Mysterud, A. Infection prevalence and ecotypes of Anaplasma phagocytophilum in moose Alces alces, red deer Cervus elaphus, roe deer Capreolus capreolus and Ixodes ricinus ticks from Norway. Parasit. Vectors 2019, 12, 1. [Google Scholar] [CrossRef]
  41. Telford, S.R.; Dawson, J.E.; Katavlos, P.; Warner, C.C.; Kolbert, K.P.; Persing, D.H. Perpetuation of the agent of human granulocytic ehrlichiosis in a deer tick-rodent cycle. Proc. Natl. Acad. Sci. USA 1996, 93, 6209–6214. [Google Scholar] [CrossRef] [Green Version]
  42. Ogden, N.H.; Casey, A.N.J.; Woldehiwet, Z.; French, N.P. Transmission of Anaplasma phagocytophilum to Ixodes ricinus ticks from sheep in the acute and post-acute phases of infection. Infect. Immun. 2003, 71, 2071–2078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Chastagner, A.; Pion, A.; Verheyden, H.; Lourtet, B.; Cargnelutti, B.; Picot, D.; Poux, V.; Bard, É.; Plantard, O.; McCoy, K.D.; et al. Host specificity, pathogen exposure, and superinfections impact the distribution of Anaplasma phagocytophilum genotypes in ticks, roe deer, and livestock in a fragmented agricultural landscape. Infect. Genet. Evol. 2017, 55, 31–44. [Google Scholar] [CrossRef] [PubMed]
  44. Tegtmeyer, P.; Ganter, M.; von Loewenich, F.D. Simultaneous infection of cattle with different Anaplasma phagocytophilum variants. Ticks Tick Borne Dis. 2019, 10, 1051–1056. [Google Scholar] [CrossRef] [PubMed]
  45. Majazki, J.; Wüppenhorst, N.; Hartelt, K.; Birtles, R.; von Loewenich, F.D. Anaplasma phagocytophilum strains from small mammals exhibit specific ankA gene sequences. BMC Vet. Res. 2013, 9, 235. [Google Scholar] [CrossRef] [Green Version]
  46. Silaghi, C.; Gilles, J.; Hohle, M.; Fingerle, V.; Just, F.T.; Pfister, K. Anaplasma phagocytophilum infection in Ixodes ricinus, Bavaria, Germany. Emerg. Infect. Dis. 2008, 14, 972–974. [Google Scholar] [CrossRef]
  47. Silaghi, C.; Kohn, B.; Chirek, A.; Thiel, C.; Nolte, I.; Liebisch, G.; Pfister, K. Relationship of molecular and clinical findings on Anaplasma phagocytophilum involved in natural infections of dogs. J. Clin. Microbiol. 2011, 49, 4413–4414. [Google Scholar] [CrossRef] [Green Version]
  48. Silaghi, C.; Liebisch, G.; Pfister, K. Genetic variants of Anaplasma phagocytophilum from 14 equine granulocytic anaplasmosis cases. Parasit. Vectors 2011, 4, 161. [Google Scholar] [CrossRef] [Green Version]
Table 1. Geographical location of tick sampling sites.
Table 1. Geographical location of tick sampling sites.
CountryAcronymName of the SiteLandscapeGeographic Coordinates
GermanyA1Englischer GartenUrban park48°09′01.73″ N11°35′24.19″ E
A2Nymphenburger SchlossparkUrban park48°09′38.93″ N11°29′33.31″ E
A3DörnbergparkUrban park49°00′55.72″ N12°05′08.89″ E
A4Schlosspark BergForest park47°57′43.85″ N11°20′53.35″ E
A5KerschlachAgricultural47°54′57.18″ N11°12′44.04″ E
A6TussenhausenNatural48°06′36.42″ N10°34′33.40″ E
Czech RepublicB1ValticeUrban park48°44′05.68″ N16°45′11.31″ E
B2PohanskoNatural48°43′37.68″ N16°53′08.35″ E
B3Suchovské mlýnyPastureland48°53′50.79″ N17°34′54.94″ E
B4ProskoviceNatural49°44′51.73″ N18°12′23.45″ E
B5Bělský lesUrban park49°47′5.65″ N18°14′28.29″ E
SlovakiaC1SAS campusUrban park48°10′14.58″ N17°04′1.15″ E
C2Mokrohájska streetUrban48°10′34.71″ N17°04′3.15″ E
C3Železná studienkaForest park48°12′14.16″ N17°05′47.05″ E
C4FúgelkaNatural48°22′44.14″ N17°18′52.86″ E
Table 2. Primers and probes used in PCR reactions for the amplification of Anaplasma phagocytophilum gene fragments.
Table 2. Primers and probes used in PCR reactions for the amplification of Anaplasma phagocytophilum gene fragments.
GenePrimers and Probes (P)Sequences (5′-3′)Ref.
msp2ApMSP2fATGGAAGGTAGTGTTGGTTATGGTATT[26]
ApMSP2rTTGGTCTTGAAGCGCTCGTA
ApMSP2p (P)HEX-TGGTGCCAGGGTTGAGCTTGAGATTG-TAMRA
16S rRNA1st amplification: [27]
ge3aCACATGCAAGTCGAACGGATTATTC
ge10rTTCCGTTAAGAAGGATCTAATCTCC
2nd amplification:
ge9fAACGGATTATTCTTTATAGCTTGCT
ge2GGCAGTATTAAAAGCAGCTCCAGG
groELEphplgroEL-A.phago-FATGGTATGCAGTTTGATCGC[28]
EphgroEL-A.phago-RTTGAGTACAGCAACACCACCGGAA
Table 3. Frequency of Anaplasma phagocytophilum ecotypes identified in questing Ixodes ricinus ticks based on 16S rRNA and groEL gene sequences. In the logit model describing the probability that a particular ecotype is present in an I. ricinus tick, the only significant factors are presence/absence of roe deer and tick developmental stage.
Table 3. Frequency of Anaplasma phagocytophilum ecotypes identified in questing Ixodes ricinus ticks based on 16S rRNA and groEL gene sequences. In the logit model describing the probability that a particular ecotype is present in an I. ricinus tick, the only significant factors are presence/absence of roe deer and tick developmental stage.
Roe DeerTick StageCountrySiteLandscapeEcotype IEcotype IIEcotype I & II
AbsentAGermanyA1UP14
A3UP47
CzechB3A 12
RepublicB5UP1 1
SlovakiaC2U121
NGermanyA1UP3
A3UP5
CzechB3A1 4
RepublicB5UP41
SUBTOTAL8737
%89.73.17.2
PresentAGermanyA2UP88
A4FP 5
A5A26
A6 N110
CzechB1UP5
RepublicB2N2 1
B4N5
SlovakiaC1U1361
C3FP1591
C4N1957
NGermanyA2UP1
A6N 2
CzechB1UP211
RepublicB4N32
SlovakiaC1U 7
C3FP792
C4N551
SUBTOTAL887514
%49.742.47.9
TOTAL 1757821
% 63.928.57.7
A = adult; N = nymph; for site acronyms see Table 1. UP = urban park, FP = forest park, A = agricultural, pastureland, N = natural forest, U = urban.
Table 4. Frequency of Anaplasma phagocytophilum variants identified in questing Ixodes ricinus ticks based on a fragment of the 16S rRNA sequence.
Table 4. Frequency of Anaplasma phagocytophilum variants identified in questing Ixodes ricinus ticks based on a fragment of the 16S rRNA sequence.
Roe DeerTick StageCountrySiteLandscape 16S rRNA Gene Variant
ABIJNOPSVWXYZ
AbsentAGermanyA1UP13 1
A3UP48
CzechB3A1 21
RepublicB5UP4
SlovakiaC2U13 1
NGermanyA1UP3
A3UP6
CzechB3A1 4
RepublicB5UP4 1
SUBTOTAL93 1 1611
%90.3 1.0 1.05.81.01.0
PresentAGermanyA2UP8 1 71
A4FP 32
A5A 2 123
A6 N 1 1 23 4
CzechB1UP5
RepublicB2N 1 1 1
B4N6
SlovakiaC1U14 1 5
C3FP15 11 261
C4N131 1 9173
NGermanyA2UP1
A6N 2
CzechB1UP21 1
RepublicB4N4 1
SlovakiaC1U 4 41
C3FP7 2 9
C4N4 1 241
SUBTOTAL7962211 1191538231
%42.03.21.11.10.50.5 5.84.88.020.212.20.5
TOTAL 17262211111102139241
% 59.12.10.70.70.30.30.33.83.47.213.48.20.3
A = adult; N = nymph; for site acronyms, see Table 1. Designation of 16S rRNA gene variants follows Schorn et al. [31] and Silaghi et al. [32]. UP = urban park, FP = forest park, A = agricultural, pastureland, N = natural forest, U = urban.
Table 5. Mean prevalence of Anaplasma phagocytophilum in questing Ixodes ricinus ticks per site for the study period.
Table 5. Mean prevalence of Anaplasma phagocytophilum in questing Ixodes ricinus ticks per site for the study period.
RoeAdultsNymphsTOTAL
CountrySiteLandscapeDeerPositive/Total%Positive/Total%Positive/Total%
GermanyA1UP-29/4806.06/2382.535/7184.9
A2UP+24/4914.92/2600.826/7513.5
A3UP-93/47519.69/2403.7102/71514.3
A4FP+8/3902.01/2260.49/6061.5
A5A+14/4593.00/260014/7191.9
A6N+24/3057.94/2400.828/5455.1
Czech B1UP+11/1786.26/2372.517/4154.1
RepublicB2N+3/378.10/17003/2071.4
B3A-5/865.87/4561.512/5422.2
B4N+11/8313.312/11141.123/11971.9
B5UP-13/9613.513/1807.226/2769.4
SlovakiaC1UP+139/39735.070/66310.6209/106019.7
C2U-23/19311.90/13023/20611.2
C3FP+47/34413.732/4047.979/74810.6
C4N+30/5205.820/13621.550/18822.7
Data are summarized from publications by Silaghi et al. [19], Overzier et al. [20], Venclikova et al. [22,23], and Svitálková et al. [24]. For site acronyms, see Table 1. UP = urban park, FP = forest park, A = agricultural, pastureland, N = natural forest, U = urban. +/- = presence/absence of roe deer.

Share and Cite

MDPI and ACS Style

Hamšíková, Z.; Silaghi, C.; Takumi, K.; Rudolf, I.; Gunár, K.; Sprong, H.; Kazimírová, M. Presence of Roe Deer Affects the Occurrence of Anaplasma phagocytophilum Ecotypes in Questing Ixodes ricinus in Different Habitat Types of Central Europe. Int. J. Environ. Res. Public Health 2019, 16, 4725. https://doi.org/10.3390/ijerph16234725

AMA Style

Hamšíková Z, Silaghi C, Takumi K, Rudolf I, Gunár K, Sprong H, Kazimírová M. Presence of Roe Deer Affects the Occurrence of Anaplasma phagocytophilum Ecotypes in Questing Ixodes ricinus in Different Habitat Types of Central Europe. International Journal of Environmental Research and Public Health. 2019; 16(23):4725. https://doi.org/10.3390/ijerph16234725

Chicago/Turabian Style

Hamšíková, Zuzana, Cornelia Silaghi, Katsuhisa Takumi, Ivo Rudolf, Kristyna Gunár, Hein Sprong, and Mária Kazimírová. 2019. "Presence of Roe Deer Affects the Occurrence of Anaplasma phagocytophilum Ecotypes in Questing Ixodes ricinus in Different Habitat Types of Central Europe" International Journal of Environmental Research and Public Health 16, no. 23: 4725. https://doi.org/10.3390/ijerph16234725

APA Style

Hamšíková, Z., Silaghi, C., Takumi, K., Rudolf, I., Gunár, K., Sprong, H., & Kazimírová, M. (2019). Presence of Roe Deer Affects the Occurrence of Anaplasma phagocytophilum Ecotypes in Questing Ixodes ricinus in Different Habitat Types of Central Europe. International Journal of Environmental Research and Public Health, 16(23), 4725. https://doi.org/10.3390/ijerph16234725

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