- Open Access
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
Parasites & Vectorsvolume 12, Article number: 1 (2019)
The geographical expansion of the tick Ixodes ricinus in northern Europe is a serious concern for animal and human health. The pathogen Anaplasma phagocytophilum is transmitted by ticks and causes emergences of tick-borne fever (anaplasmosis) in livestock. The transmission dynamics of the different ecotypes of A. phagocytophilum in the ecosystems is only partly determined. Red deer and roe deer contribute to circulation of different ecotypes of A. phagocytophilum in continental Europe, while the role of moose for circulation of different ecotypes is not fully established but an important issue in northern Europe.
We determined infection prevalence and ecotypes of A. phagocytophilum in moose (n = 111), red deer (n = 141), roe deer (n = 28) and questing ticks (n = 9241) in Norway.
As previously described, red deer was exclusively linked to circulation of ecotype I, while roe deer was exclusively linked to circulation of ecotype II. Surprisingly, we found 58% ecotype I (n = 19) and 42% of ecotype II (n = 14) in moose. Both ecotypes were found in questing ticks in areas with multiple cervid species present, while only ecotype I was found in ticks in a region with only red deer present. Hence, the geographical distribution of ecotypes in ticks followed the distribution of cervid species present in a given region and their link to ecotype I and II.
Moose probably function as reservoirs for both ecotype I and II, indicating that the ecotypes of A. phagocytophilum are not entirely host-specific and have overlapping niches. The disease hazard depends also on both host abundance and the number of immature ticks fed by each host. Our study provides novel insights in the northern distribution and expansion of tick-borne fever.
Pathogens transmitted by Ixodes ricinus cause disease emergence in northern areas of Europe with serious implications for animal and human health [1,2,3]. Changes in the distribution and local abundances of the tick Ixodes ricinus is one of the main factors for increases in disease incidences , but several other processes, such as changes in vertebrate communities, also affect disease incidences [5, 6]. Fundamental knowledge of disease hazard requires understanding of the role of vertebrate hosts in the intricate links between the tick life-cycle and the enzootic cycles of the pathogens they transmit . The distribution and abundance of tick-borne pathogens often depends on distribution of competent transmission hosts, not just distribution of the vector, and such insight provides another potential key to disease mitigation efforts.
Anaplasma phagocytophilum is the causative agent of human granulocytic anaplasmosis , and it also causes disease and economic losses in livestock . Tick-borne fever, caused by A. phagocytophilum, is considered a major problem for livestock production. Anaplasmosis leads to mortality and especially lambs are at risk . Infection with A. phagocytophilum also reduces body growth of lambs . Some tick-borne pathogens, such as Borrelia burgdorferi (sensu lato) and A. phagocytophilum, have a high genetic diversity. Due to the high genetic diversity, A. phagocytophilum was formerly even split in three different species . This genetic diversity has probably arisen as a response to evasion of immune defenses of phylogenetically different vertebrate host groups .
Elucidating the presence of different enzootic cycles for A. phagocytophilum is important to be able to determine which ones are pathogenic for livestock and which ones for humans, but the list of variants and transmission hosts are still incompletely understood [14, 15]. Common markers used for genotyping are 16S-rDNA, groEL, ankA and msp4 . One classification of the groEL gene has identified the circulation of four ecotypes among vertebrates of Europe . We define in this paper an ecotype as a cluster of genetically similar A. phagocytophilum isolates based on groEL sequences following Jahfari et al. . Ecotype I based on this classification is linked to a wide range of mammals including red deer (Cervus elaphus), livestock and humans, whereas ecotype II is linked mainly to roe deer (Capreolus capreolus). Both ecotypes are most probably vectored by I. ricinus. Ecotype III is associated with small mammals and expected to be vectored by Ixodes trianguliceps , while ecotype IV is linked to birds and potentially vectored by I. frontalis . In support of this scheme, the ecotype II found in roe deer was for example not the one found in livestock in France . However, due to scarcity of studies, this classification in ecotypes can be considered as a working hypothesis of circulation and transmission, rather than as a definitive answer.
Ixodes ricinus ticks are now expanding their distribution range towards northern latitudes [2, 18]. This also involves a shift into the distribution range of another community of mammals compared to continental Europe. In the case of A. phagocytophilum, potential cervid transmission hosts are no longer limited to roe deer and red deer, but also include moose (Alces alces). Moose is by far the most important game species in terms of meat and economy both in Sweden and Norway, with an annual harvest in 2016 of, respectively, 82,097  and 30,829  individuals. Moose is widespread in Fennoscandia, and it is important to assess its role in the propagation of ticks and the transmission of tick-borne diseases [21, 22]. Yet, we do not know much about the potential role of moose as a transmission host of A. phagocytophilum. Seroprevalence for A. phagocytophilum of moose in Norway [22, 23], and infection levels of moose in both Norway  and Sweden  was reported to be high. However, the ecotype of A. phagocytophilum has only been determined in two moose with no firm conclusion, as one type I and one type II was found (data from  reported in ). Further, the geographical distribution of red deer is quite rapidly expanding in both Norway and Sweden. This may introduce the assumed pathogenic variant (ecotype I) to new areas, currently with low incidence of anaplasmosis [3, 25]. Verification of presence of a pathogenic variant in cervids is therefore important.
We here report infection prevalence and determine ecotypes of A. phagocytophilum in moose (n = 111), red deer (n = 141), roe deer (n = 28) as well as in 9241 questing Ixodes ricinus ticks (nymphs, male and female adults) from three regions in Norway.
Samples come from different sites in the southern part of Norway (Fig. 1). Roe deer came from around Vestby municipality, Akershus county; moose came from Siljan municipality in Telemark county, Vennesla, Songdalen, Audnedal and Marnardal municipality in Vest-Agder county, and Andebu, Lardal and Re municipality in Vestfold county; red deer came from Florø municipality, Sogn & Fjordane county, Drangedal and Skien municipality in Telemark county and Lardal in Vestfold county; questing ticks came from counties Møre & Romsdal and Sogn og Fjordane in the west and Akershus and Østfold in the east.
Sample collection and DNA extraction
Questing ticks were sampled by aid of flagging along transects in Sogn og Fjordane, Møre og Romsdal and Akershus and Østfold counties, Norway. Further details on sampling protocols of ticks can be found elsewhere [25,26,27]. For ticks, as in our previous work [25, 28], DNA extraction was optimized by modifying the incubation step of Qiagen DNeasy® 96 Blood & Tissue kit as recommended . We incubated ticks, 2 mm zirconium oxide beads, 40 μl of proteinase K solution and 4 μl Antifoam-A (Sigma) at 56 °C overnight followed by 5 min of bead homogenization at 30 cycles per second using a Qiagen Tissue Lyser II. The samples were then frozen. This was done to more easily prepare high number of samples that later could be extracted. Later we added the remaining 160 μl of proteinase K solution and the homogenized mixture was incubated for 1 h at 56 °C. Centrifugation was carried out before removing seals at any step in the protocol to prevent cross contamination of samples. We then transferred the mixture to the DN easy plates.
Traffic killed roe deer were collected during spring 2016. The roe deer (n = 28) were skinned, frozen and then transported to CEES, UiO. Red deer (n = 141) and moose (n = 111, 29 yearling males, 19 yearling females, 45 adult males and 18 adult females) exterior pinna (ears) were collected by the Norwegian Food Safety Authority as part of the national monitoring for Chronic Wasting Disease during fall hunting season 2016 in Norway, and then sent as frozen samples to CEES, UiO.
A small piece (1/8 cm2) from the frozen exterior pinna was collected with a clean disposable scalpel and used for the DNA extraction. DNA was extracted in the same way as for ticks (except bead milling) with Qiagen DNeasy 96 Blood & Tissue kit according to the manufactures recommendations. A total of 94 samples were extracted at a time, leaving two spaces empty for controls. The DNA was stored at -80 °C for later use.
The extracted DNA was screened for A. phagocytophilum by realtime PCR (qPCR) in a duplex (with B. burgdorferi, not reported here) , following Allender et al.  and implemented at CEES, UiO as in previous work .
For A. phagocytophilum, the msp2 gene was targeted with the primers ApMSP2f (5'-ATG GAA GGT AGT GTT GGT TAT GGT ATT-3') and ApMSP2r (5'-TTG GTC TTG AAG CGC TCG TA-3') and the probe ApMSP2p-[HEX]TGG TGC CAG GGT TGA GCT TGA GAT TG[TAMRA]. Each duplex qPCR reaction was done in a total volume of 10 μl, with 1 μl of extracted DNA added to the mastermix as a template. The mastermix consisted of PCR buffer (5×) (TaqMan® polymerase, LightCycler®, Roche, Basel, Switzerland) and PCR water, 900 nM for both the ApMSP2f and ApMSP2r primers and 125 nM for the ApMSP2p probe, while for B. burgdorferi (s.l.) it was 700 nM for both primers and 175 nM for the probe. Each 96-well plate contained a negative and a positive control. A two-step program was used on the LightCycler® 96 System (Roche). It started with a pre incubation of 10 min at 95 °C, followed by 50 cycles of a two-step amplification with (i) 15 s at 95 °C and (ii) 60 s at 60 °C.
Sequencing of groEL from positive samples
The DNA from the samples that came up positive from the qPCR were amplified by conventional PCR, targeting a 530 bp fragment of the groEL gene of A. phagocytophilum. The Hotstart mastermix (Qiagen, Hilden, Germany) was mixed with primers and template. We used the primers EphplgroEL-A.phago-F (5'-ATG GTA TGC AGT TTG ATC GC-3') and EphgroEL-A.phago-R (5'-TTG AGT ACA GCA ACA CCA CCG GAA-3') for PCR amplification .
The PCR settings were a pre incubation of 95 °C for 15 min followed by 40 cycles of (i) 30 s at 94 °C, (ii) 30 s at 57 °C, (iii) 45 s at 72 °C, and a final extension of 10 min at 72 °C. To ensure that a product from the PCR was obtained, gel electrophoreses was utilized with 8 μl of DNA on a 1.5% agarose gel. The gel was stained with SYBR™ Gold Nucleic Acid Gel Stain (Invitrogen, Waltham, MA). If the PCR was successful, shown with a clear band on the gel, the PCR product was cleaned with ExoSAP-IT™ PCR Product Cleanup Reagent (Applied Biosystems, Waltham, MA) and sent for sequencing by the firm BaseClear (Leiden, the Netherlands). The chromatographs of the sequences were visually inspected and the primers sites were trimmed in Bionumerics version 7.6 (Applied Math, Belgium). Our sequences and those of known ecotypes were aligned and cluster analysis was performed as described in Jahfari et al. . The sequences from this study were deposited in the Genbank database under the accession numbers MK069678-MK069965.
Infection prevalence is binomial and therefore analysed with mixed-effects logistic regression models in R package lme4 . We used municipality as a random intercept to account for potential dependency of observations from the same region. For moose, we had data on sex and age class (yearling/adult) for all and carcass mass (kg) for 92 individuals. Carcass mass was used as a proxy for body mass . We used the Akaike information criterion to establish the most parsimonious model. For red deer, there were too few negative individuals to allow a meaningful analysis of infection pattern. The statistical modeling were done in R version 3.4.2.
We made a phylogenetic tree of the sequences to explore the link between A. phagocytophilum groEL sequences, host and region. The sequences were aligned and trimmed in MAFFT v.7.271. A maximum clade credibility tree was prepared with BEAST v1.8.4, and visualized in Figtree v1.4.3. When plotted in QGIS v. 2.18.7 , the coordinates of the 288 samples formed three clusters, which we have visualized in the geographical map that was extracted from the GADM database v. 3.4 . The southern cluster includes samples around the Oslofjord inlet. While Oslofjord forms a barrier, the scarceness of sampling on the east side made us chose to combine it into one location.
Anaplasma phagocytophilum DNA was found in 70% of the moose, in 82% of the roe deer and in 94% of the red deer (Table 1). The best model of infection prevalence in moose included both age and sex, but not the interaction term (ΔAIC = 1.592). Infection prevalence in moose was higher in males (76%) than females (59%) (Z = 2.175, df = 1, P = 0.030), but infection prevalence was not related to age class (Z = 1.553, df = 1, P = 0.120). For the subsample with known body mass, adding (log) body mass did not improve the model with age and sex included (ΔAIC = 1.958).
To distinguish between ecotype I and II of A. phagocytophilum, we amplified, sequenced and typed a 530 bp fragment of the groEL gene. We obtained successful sequences of the A. phagocytophilum positive samples in 42%, 57% and 57% of cases deriving from, respectively, moose, red deer and roe deer (Table 1). In moose, we found 58% ecotype I (n = 19) and 42% of ecotype II (n = 14) (Fig. 2). In red deer, all sequences grouped to ecotype I (n = 75). In roe deer, all sequences grouped to ecotype II (n = 13).
In questing ticks, we found ecotype I (n = 2) and ecotype II (n = 1) from a region with roe deer and moose (Akershus and Østfold county). Only ecotype I (n = 102) was found in Sogn & Fjordane county with only red deer present, and we found ecotype I (n = 38) and ecotype II (n = 5) in Møre & Romsdal county with mainly red deer and some roe deer populations (Fig. 2).
The phylogenetic analysis identified the two main clusters belonging to ecotypes I and II. The more fine scale genetic variation within ecotypes showed no clear geographical pattern or further association with host (Fig. 3).
Our study delineated ecotypes of A. phagocytophilum by standard molecular methods, based on the presence of DNA and bacterial groEL sequences originating from ear tissue of three species of cervids as well as in questing ticks coming from the northern distribution range of I. ricinus ticks in Europe. This study corroborates the preferential association of A. phagocytophilum ecotype I to red deer and ecotype II to roe deer, consistent with previous studies . Importantly and surprisingly, we found evidence of high infection prevalence of both ecotype I and II in moose. Infection prevalence is not direct evidence of transmission competence. However, the high infection rates in moose and the presence of ecotype I in questing ticks in a region without red deer, but with moose present, is at least suggestive that moose may play a role in transmission of both ecotype I and II of A. phagocytophilum.
Presence of bacterial DNA in tissue samples is not sufficient to document the viability or infectivity of a pathogen in the host. Culturing is the gold standard for proving the viability or infectiousness of microorganisms, but this is costly, time-consuming, difficult for many tick-borne pathogens, and generally suffers from a low sensitivity . PCR-array also has its limits, and there might be presence of bacteria below detectable concentrations . Also, vector competence requires larval (or nymphal) ticks to be infected from their blood meal and passing viable pathogens at the nymphal (or adult female) stage to a novel host. However, since many studies are indicating I. ricinus as the vector and cervids as potential hosts of A. phagocytophilum, DNA-based detection provides useful information when aiming to elucidate the different enzootic cycles. There are, however, several other difficult issues in determining circulation and enzootic cycles of pathogens. Adaptations of bacteria to specific host groups require mutations in coding genes, but many mutations are neutral and single mutation sequence dissimilarity is not sufficient in the meaning of ecotype characterization or bacterial speciation . That we found both ecotype I and II in moose, can be seen as evidence that ecotype are not entirely host-specific and have overlapping niches. The distinction between ecotype I and II are less marked than the distinction between ecotype I/II and the other two ecotypes, involving many different ruminant hosts. The ankA gene identifies some of the same correlation to host species as groEL . Both methods group sequences so that roe deer is probably not the major reservoir of granulocytic anaplasmosis in humans and domestic animals, and that red deer could be the reservoir for that ecotype. A limitation of our methodology is that co-infection of ecotypes would not be discovered. The high infection prevalence of both ecotypes in moose makes it likely that mixed infections can be common. This could be addressed with other methods such as Reverse Line Blot in future work.
Obtaining sequence data from positive samples of A. phagocytophilum required for determination of ecotype can be challenging. We had a fairly low success (42–57%) in obtaining sequence data, and the sequencing success was good when Ct values (i.e. the number of PCR cycles before getting a positive signal) were below 30. This highlights the need for large sample sizes. There are only two previous groEL sequences reported from A. phagocytophilum isolated from moose tissue. One sequence derived from a mainland moose in Sweden was an exact match of the sequence (GQ988754) from roe deer in Austria  clustering to ecotype II, while the sequence derived from a moose from the Island of Öland, Sweden clustered to ecotype I . There are no red deer on Öland, but still ecotype I was found there. Together with our discovery of frequent infection of ecotype I in moose, as well as in questing ticks in an area without red deer, this is indeed indicative of moose as a potential transmission host of ecotype I.
Our study provides one crucial step forward towards identifying assumed transmission hosts at the northern distribution ranges, and we have provided the first results potentially linking moose to both ecotype I and II of A. phagocytophilum. The importance of a given host for the disease hazard also depends on their abundance and the number of ticks feeding on them [40, 41]. Red deer can generally reach higher population densities than moose, while roe deer have highly variable densities that can be very high locally . We found a considerable amount of larvae attached to red deer close to the present western study site in Norway , and larvae were also found on roe deer in Germany . Lower proportion of larval ticks has been found on moose ears in Norway when compared to red deer and roe deer . We found much lower prevalence of A. phagocytophilum in questing ticks from the eastern study site with roe deer and moose compared to in the western study sites with red deer . Therefore, we still need more quantitative data to better understand the role different cervid species play in determining disease hazard.
Stuen S. Ticks and tick-borne infections in sheep in the Nordic countries. Small Ruminant Res. 2012;106:14–5.
Jore S, Viljugrein H, Hofshagen M, Brun-Hansen H, Kristoffersen A, Nygard K, et al. Multi-source analysis reveals latitudinal and altitudinal shifts in range of Ixodes ricinus at its northern distribution limit. Parasit Vectors. 2011;4:84.
Mysterud A, Jore S, Østerås O, Viljugrein H. Emergence of tick-borne diseases at northern latitudes in Europe: a comparative approach. Sci Rep. 2017;7:16316.
Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Peña A, George J-C, et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit Vectors. 2013;6:1.
Vuong HB, Chiu GS, Smouse PE, Fonseca DM, Brisson D, Morin PJ, et al. Influences of host community characteristics on Borrelia burgdorferi infection prevalence in blacklegged ticks. PloS One. 2017;12:e0167810.
Parham PE, Waldock J, Christophides GK, Hemming D, Agusto F, Evans KJ, et al. Climate, environmental and socio-economic change: weighing up the balance in vector-borne disease transmission. Phil Trans R Soc London Ser B. 2015;370:20130551.
Franke J, Hildebrandt A, Dorn W. Exploring gaps in our knowledge on Lyme borreliosis spirochaetes-updates on complex heterogeneity, ecology, and pathogenicity. Ticks Tick Borne Dis. 2013;4:11–25.
Dumler JS, Choi KS, Garcia-Garcia JC, Barat NS, Scorpio DG, Garyu JW, et al. Human granulocytic anaplasmosis and Anaplasma phagocytophilum. Emerg Infect Dis. 2005;11:1828–34.
Woldehiwet Z. Anaplasma phagocytophilum in ruminants in Europe. Ann NY Acad Sci. 2006;1078:446–60.
Stuen S, Bergström K. Serological investigation of granulocytic Ehrlichia infection in sheep in Norway. Acta Vet Scand. 2001;42:331–8.
Grøva L, Olesen I, Steinshamn H, Stuen S. Prevalence of Anaplasma phagocytophilium infection and effect on lamb growth. Acta Vet Scand. 2011;53:30.
Dumler JS, Barbet AF, Bekker CPJ, Dasch GA, Palmer GH, Ray SC, et al. Reorganization of genera in the families Rickettsiaceae and Anaplasmataceae in the order Rickettsiales: unification of some species of Ehrlichia with Anaplasma, Cowdria with Ehrlichia and Ehrlichia with Neorickettsia, descriptions of six new species combinations and designation of Ehrlichia equi and ‘HGE agent’ as subjective synonyms of Ehrlichia phagocytophila. Int J Syst Evol Microbiol. 2001;51:2145–65.
Estrada-Peña A, de la Fuente J, Ostfeld RS, Cabezas-Cruz A. Interactions between ticks and transmitted pathogens evolved to minimise competition through nested and coherent networks. Sci Rep. 2015;5:10361.
Scharf W, Schauer S, Freyburger F, Petrovec M, Schaarschimdt-Kiener D, Liebisch G, et al. Distinct host species correlate with Anaplasma phagocytophilum ankA gene clusters. J Clin Microbiol. 2011;49:790.
Jahfari S, Coipan EC, Fonville M, van Leeuwen AD, Hengeveld P, Heylen D, et al. Circulation of four Anaplasma phagocytophilum ecotypes in Europe. Parasit Vectors. 2014;7:365.
Bown KJ, Lambin X, Ogden NH, Begon M, Telford G, Woldehiwet Z, et al. Delineating Anaplasma phagocytophilum ecotypes in coexisting, discrete enzootic cyles. Emerg Infect Dis. 2009;15:1948–54.
Chastagner A, Pion A, Verheyden H, Lourtet B, Cargnelutti B, Picot 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.
Jaenson TGT, Eisen L, Comstedt P, Mejlon HA, Lindgren E, Bergström S, et al. Risk indicators for the tick Ixodes ricinus and Borrelia burgdorferi sensu lato in Sweden. Med Vet Entomol. 2009;23:226–37.
Älgdata. http://www.algdata.se/sv/pages/default.aspx . Länssyrelserna; 2018. Accessed 13 Nov 2018.
Statistics Norway. Official hunting statistics of Norway. http://www.ssb.no/jord-skog-jakt-og-fiskeri/statistikker/elgjakt. Oslo and Kongsvinger: Statistics Norway; 2018. Accessed 13 Nov 2018.
Malmsten J, Widén DG, Rydevik G, Yon L, Hutchings MR, Thulin CG, et al. Temporal and spatial variation in Anaplasma phagocytophilum infection in Swedish moose (Alces alces). Epidemiol Infect. 2014;142:1205–13.
Milner JM, van Beest F. Ecological correlates of a tick-borne disease, Anaplasma phagocytophilum, in moose in southern Norway. Eur J Wildl Res. 2013;59:399–406.
Stuen S, Åkerstedt J, Bergström K, Handeland K. Antibodies to granulocytic Ehrlichia in moose, red deer, and roe deer in Norway. J Wildl Dis. 2002;38:1–6.
Puraite I, Rosef O, Paulauskas A, Radzijevskaja J. Anaplasma phagocytophilum infection in moose (Alces alces) in Norway. Microb Inf. 2015;17:823–8.
Mysterud A, Stigum VM, Seland IV, Herland A, Easterday WJ, Jore S, et al. Tick abundance, pathogen prevalence, and disease incidence in two contrasting regions at the northern distribution range of Europe. Parasit Vectors. 2018;11:309.
Qviller L, Risnes-Olsen N, Bærum KM, Meisingset EL, Loe LE, Ytrehus B, et al. Landscape level variation in tick abundance relative to seasonal migration pattern of red deer. PloS One. 2013;8:e71299.
Mysterud A, Easterday WR, Stigum VM, Aas AB, Meisingset EL, Viljugrein H. Contrasting emergence of Lyme disease across ecosystems. Nature Comm. 2016;7:11882.
Mysterud A, Easterday WR, Qviller L, Viljugrein H, Ytrehus B. Spatial and seasonal variation in prevalence of Anaplasma phagocytophilum and Borrelia burgdorferi in Ixodes ricinus ticks in Norway. Parasit Vectors. 2013;6:187.
Allender CJ, Easterday WJ, Van Ert MN, Wagner DM, Keim P. High-throughput extraction of arthropod vector and pathogen DNA using bead milling. Biotechniques. 2004;37:730–4.
Courtney JW, Kostelnik LM, Zeidner NS, Massung RF. Multiplex real-time PCR for detection of Anaplasma phagocytophilum and Borrelia burgdorferi. J Clin Microbiol. 2004;42:3164–8.
Alberti A, Zobba R, Chessa B, Addis MF, Sparagano O, Pinna Parpaglia ML, et al. 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–22.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:45.
Langvatn R. Criteria of physical condition, growth and development in Cervidae, suitable for routine studies. Stockholm: Nordic Council for Wildlife Research; 1977.
. QGIS Geographic Information Team. QGIS Geographic Information System. Open Source Geospatial Foundation; 2018. https://qgis.org/en/site/. Accessed 13 Nov 2018.
GADM Database of Global Administrative Areas version 3.4. 2018. https://gadm.org/. Accessed 13 Nov 2018.
Tijsse-Klasen E, Koopmans MPG, Sprong H. Tick-borne pathogen - reversed and conventional discovery of disease. Front Publ Health. 2014;2:73.
Lachish S, Gopalaswamy AM, Knowles SCL, Sheldon BC. Site-occupancy modelling as a novel framework for assessing test sensitivity and estimating wildlife disease prevalence from imperfect diagnostic tests. Methods Ecol Evol. 2012;3:339–48.
Cohan FM, Kane M. Bacterial species and speciation. Syst Biol. 2001;50:513–24.
Silaghi C, Hamel D, Thiel C, Pfister K, Passos LMF, Rehbein S. Genetic variants of Anaplasma phagocytophilum in wild caprine and cervid ungulates from the Alps in Tyrol, Austria. Vector-Borne Zoonot. 2011;11:355–62.
Ostfeld RS, Canham CD, Oggenfuss K, Winchcombe RJ, Keesing F. Climate, deer, rodents, and acorns as determinants of variation in Lyme-disease risk. Plos Biol. 2006;4:1058–68.
Hofmeester TR, Coipan EC, Van Wieren SE, Prins HHT, Takken W, Sprong H. Few vertebrate species dominate the Borrelia burgdorferi s.l. life cycle. Environ Res Lett. 2016;11:043001.
Austrheim G, Solberg EJ, Mysterud A. Spatio-temporal variation in large herbivore pressure in Norway from 1949–1999: has decreased grazing by livestock been countered by increased browsing by cervids? Wildl Biol. 2011;17:286–98.
Mysterud A, Hatlegjerde IL, Sørensen OJ. Attachment site selection of life stages of Ixodes ricinus ticks on a main large host in Europe, the red deer (Cervus elaphus). Parasit Vectors. 2014;7:510.
Kiffner C, Lödige C, Alings M, Vor T, Rühe F. Abundance estimation of Ixodes ticks (Acari: Ixodidae) on roe deer (Capreolus capreolus). Exp Appl Acarol. 2010;52:73–84.
Handeland K, Qviller L, Vikøren T, Viljugrein H, Lillehaug A, Davidson RK. Ixodes ricinus infestation in free-ranging cervids in Norway - a study based upon ear examinations of hunted animals. Vet Parasitol. 2013;195:142–9.
We are grateful to Vidar Holthe for collection of traffic-killed roe deer, and to Katsuhisa Takumi for valuable help with the phylogenetic analysis.
We are grateful to the Research Council of Norway and the Norwegian Environment Agency for funding the #EcoTick project (254469), and to the Faculty of Mathematics and Natural Sciences, UiO, for giving a scholarship funding a sabbatical of AM of University of Antwerp facilitating this study.
Availability of data and materials
Data on proportions of ecotypes are provided in Table 1. The sequences from this study were deposited in the GenBank database under the accession numbers MK069678-MK069965.
Ethics approval and consent to participate
The paper contains only road killed and hunter killed animals, i.e. no particular permission is required for sampling.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.