Prevalence of Neoehrlichia mikurensis in ticks and rodents from North-west Europe

1The most prevalent tick-borne infection of humans in the Northern hemisphere is Lyme The same tick species transmitting the etiologic agents of Lyme disease also serve as the vector of pathogens causing tick-borne encephalitis and several forms of rickettsioses, anaplasmoses and ehrlichioses. Members of the family Anaplasmataceae are obligatory intracellular bacteria that reside within membrane-enclosed vacuoles. Human ehrlichiosis and anaplasmosis are two closely related diseases caused by various members of the genera Ehrlichia and Anaplasma. A major difference between these two members is their cellular tropism. Ehrlichia chaffeensis, the etiologic agent of human monocytotropic ehrlichiosis (HME), is an emerging zoonosis that causes clinical manifestations ranging from a mild febrile illness to a fulminant disease characterized by multi-organ system failure. Anaplasma phagocytophilum causes human granulocytotropic anaplasmosis (HGA), previously known as human granulocytotropic ehrlichiosis. Despite the presence of Anaplasma phagocytophilum in questing Ixodes ricinus ticks in the Netherlands, only one human case has been reported. Seropositivity against anaplasmosis was observed in risk groups, such as foresters and suspected Lyme disease patients, but not in control groups. Still, the incidence of these tick-borne diseases and the associated public health risks remain largely unknown.

2A novel candidate species in the family of Anaplasmataceae, called Candidatus Neoehrlichia mikurensis (N. mikurensis), was first isolated from wild rats and was also found in I. ovatus in Japan. Neoehrlichia mikurensis can be distinguished from other genera based on sequence analysis of 16S rDNA, citrate synthase (gltA) and heat shock protein GroEL genes. This recently identified bacterium is detected in several tick species and rodents in different parts of the world under different names. The N. mikurensis found in I. ricinus ticks in Italy has been referred to as Candidatus Ehrlichia walkerii and the Ehrlichia species isolated from a rat in China was called “Rattus strain”. Furthermore, a N. mikurensis has been described in I. persulcatus in Russia and I. ovatus from China and Japan. In the US, an Ehrlichia-like organism, closely related to N. mikurensis, was previously detected in raccoons. This variant is called Candidatus Neoehrlichia lotoris. The Asian N. mikurensis isolates showed a 99% similarity based on the 16S rDNA to the Ehrlichia schotti. Ehrlichia schotti was first described in 1999 in I. ricinus in the Netherlands by Leo Schouls and was named after his technician. Later this species was reported in I. ricinus in Russia and subsequently in Germany and Slovakia. These findings raised the question whether Ehrlichia schotti is the same as N. mikurensis.

3It is unclear whether N. mikurensis poses a risk to public health. Until recently, there were no human infections reported. In 2010, the first case of human N. mikurensis infection was reported in a patient from Sweden. In the same year, five other human infections were described in Germany, Switzerland and the Czech Republic. More recently, a canine infection was reported in Germany. The symptoms described in all of these cases were generally non-specific and usually seen in any other ordinary inflammation reaction (Table 1). These reported cases of human infections imply that re-evaluation is needed regarding the pathogenesis of this species. All but one case that have been described so far have occurred in patients who were immuno-compromised. The non-specificity of the reported symptoms, poor diagnostic tools and the lack of awareness of public health professionals could explain the absence of (reported) patients.

4In this study we aim to investigate (i) whether Ehrlichia schotti is similar to the described N. mikurensis family, (ii) the distribution and prevalence of N. mikurensis in the Netherlands, Belgium and the UK, (iii) possible transmission routes of N. mikurensis in non-experimental settings and (vi) its putative mammalian hosts.

5Questing I. ricinus from all stages and Dermacentor reticulatus adults were collected in 2009 and 2010 by flagging the vegetation at geographically different locations in the Netherlands and Belgium. Ticks collected in the UK and Vrouwenpolder (NL) have been described before. For global geographic location, see Additional file 1: Figure S1. Questing I. arboricola were collected from bird nests in two different areas in Belgium. Ixodes hexagonus feeding on hedgehogs were collected in a hedgehog-shelter in 2010. Ixodes ricinus feeding on red deer (Cervus elaphus), European mouflon (Ovis orientalis musimon), wild boar (Sus scrofa) and sheep (Ovis aries) were collected. All the collected ticks were immersed in 70% alcohol and stored at −20°C until the DNA extraction. Based on morphological criteria, tick species and stages were identified to species level, with stage and sex recorded. In doubtful cases, sequencing of tick mitochondrial 16S rDNA confirmed the tick-species. DNA from questing ticks was extracted by alkaline lysis. DNA from engorged ticks was extracted using the Qiagen DNeasy Blood & Tissue Kit according to the manufacturer’s manual (Qiagen, 2006, Hilden; Germany) following the manufacturer’s protocol for the purification of total DNA from ticks.

6Longworth traps (Bolton Inc., UK), baited with hay, apple, carrot, oatmeal and mealworm were used to capture different species of rodents and insectivores at 7 different locations in the Netherlands between 2007 and 2010. Animals were anaesthetized with isoflurane and euthanized by cardiac puncture. Serum was collected and stored at – 20°C. Spleen, liver, kidney, brain and other organs were collected and frozen at −80°C. DNA was extracted using the Qiagen DNeasy Blood & Tissue Kit according to the manufacturer’s manual (Qiagen, 2006, Hilden; Germany). All animals were handled in compliance with Dutch laws on animal handling and welfare (RIVM/DEC permits).

7Polymerase chain reaction (PCR) amplifications were performed in a Px2 Thermal Cycler (Thermo Electron Corporation, Waltham, Massachusetts, USA). The presence of Ehrlichia schotti in questing I. ricinus was studied by Reverse Line Blotting as described. Fragments of the 16S rDNA, citrate synthase gene gltA, and the chaperonin GroEL of ehrlicial species were amplified from tick lysates and rodent tissue samples using novel primers and primers that were previously described (Table 2). Amplification of gltA and GroEL were both done in 50 μl reaction volumes containing 5 μl template DNA. GltA DNA was amplified using a final concentration of 800 nM of each primer, NMik fo-gltA and NMik re-gltA with the following PCR program, 15 min at 95°C, 40 cycles each consisting of 30 sec at 94°C, 25 sec at 53°C, and 10 min at 72°C. GroEL DNA was amplified using, 500 nM of each primer NMik fo-groEL and NMik re-groEL. The PCR program used is as followed: 15 min at 95°C, 40 cycles each consisting of 30 sec at 94°C, 75 sec at 49°C, and 10 min at 72°C. The nested reaction was carried out at the same temperature as the first reactions; only 25 cycles were carried out with 1 μL of the first amplification product. The HotStarTaq Polymerase Kit (Qiagen) was used for all PCR experiments. PCR products were detected by electrophoresis in a 1.5% agarose gel stained with SYBR gold (invitrogen).

8Oligonucleotide primer and probe sequences were designed to be specific for the N. mikurensis GroEL gene using Visual OMP DNA (Software, Inc., Ann Arbor, USA). Primer sequences for the N. mikurensis GroEL gene were NMikGroEL-F2a and NMikGroEL-R2b and generated a 99-bp fragment which was detected with the NMikGroEL-P2a TaqMan probe (Table 2). Sequences were evaluated on the basis of the following criteria: predicted cross-reactivity with closely related organisms, internal primer binding properties for hairpin and primer-dimer potential, length of the desired amplicon, G-C content, and melting temperatures (Tms) of probes and primers. The specificity of the N. mikurensis GroEL primers for N. mikurensis in the multiplex real-time PCR assay was tested with DNA extracted from the following microorganisms: Rickettsia rickettsiiAnaplasma phagocytophilumR. helveticaBartonella henselaeEhrlichia canisB. afzeliiB. garinii, B. sensu stricto, Babesia microti, Candidatus Midichloria mitochondrii and tick lysates containing Wolbachia species None were amplified. Random samples of tick lysates which were N. mikurensis-positive in the Q-PCR were routinely confirmed by conventional PCR using NMik fo-gltA and NMik re-gltA primers, followed by DNA sequencing.

9PCR was performed in a multiplex format with a reaction volume of 20 μl, using the iQ Multiplex Powermix (Bio-Rad Laboratories, Hercules, USA), in the LightCycler 480 Real-Time PCR System (F. Hoffmann-La Roche, Basel, Switzerland). Final PCR reaction concentrations were 1x iQ Powermix, primers ApMSP2F and ApMSP2R at 250 nM each, probe ApMSP2P-FAM at 125 nM, primers NMikGroEL-F2a and NMikGroEL-R2b at 250 nM each, probe NMikGroEL-P2a-RED at 250 nM, and 3 μl of template DNA. Cycling conditions were: 95°C for 5 min, followed by 60 cycles of a 5 sec denaturation at 95°C followed by a 35 sec annealing-extension step at 60°C. Ticks lysates were considered positive if the Ct-value of a proper sigmoid curve was maximally three cycles more than the highest dilution of the positive control sample.

10For each PCR and real-time multiplex PCR, positive, negative controls and blank samples were included. A 10-3 to 10-5 dilution of a mixture of sequencing-confirmed N. mikurensis-positive tick lysates were used as positive controls. In order to minimize contamination, the reagent setup, the extraction and sample addition, and the real-time PCR as well as sample analysis were performed in three separate rooms, of which the first two rooms were kept at positive pressure and had airlocks.

11PCR amplicons were sequenced using the described primers (Table 2) and the BigDye Terminator Cycle sequencing Ready Reaction kit (Perkin Elmer, Applied Biosystems). All sequences were confirmed by sequencing both strands. Sequences were compared with sequences in Genbank using BLAST after subtraction of the primer sequences (http://www.ncbi.nlm.nih.gov/genbank/). The collected sequences were assembled, edited, and analysed with BioNumerics version 6.5 (Applied Maths NV, Sint-Martens-Latem, Belgium). Resulting sequences were aligned with those from related organisms in Genbank. Phylogenetic analyses of the sequences and related organisms were conducted using the BioNumerics program using the neighbour-joining algorithm with Kimura's two-parameter model. Bootstrap proportions were calculated by the analysis of 1000 replicates for neighbour-joining trees. DNA sequences are available upon request.

12Twenty-three tick lysates that were previously tested positive for the presence of Ehrlichia schotti by PCR and Reverse Line Blotting were amplified by PCR on the three loci 16S rDNA, gltA and GroEL using primers specific for N. mikurensis (Table 2). Amplicons of all three partial genes were obtained in 21 cases. None of these three loci were successfully amplified in 15 Ehrlichia schotti-negative ticks. The PCR products of all three loci were sequenced and compared with each other and with N. mikurensis sequences available in Genbank. All the Ehrlichia schotti sequences were identical to each other on the three loci, 16S rDNA, gltA as well as the GroEL. The 1740 base pairs of the 16S rDNA sequences from the Ehrlichia schotti were 99.6% to 100% similar to the N. mikurensis sequences and the Candidatus Ehrlichia walkerii sequence in Genbank (Table 3). The 233 base-pair fragment of the gltA sequences from the Ehrlichia schotti were identical to the Candidatus Ehrlichia walkerii gltA sequence (Table 3). The 1238 base-pairs of the GroEL isolates amplified from the tick lysates showed a 94.3% and 95.5%, 98.7% and 100% (AB084583 and AB074461, EF633745 and FJ966365) match with the N. mikurensis GroEL sequences in Genbank, respectively. Phylogenetic analyses of the gltA and GroEL sequences showed that the Ehrlichia schotti clustered with N. mikurensis isolates, but not with A. phagocytophilum or any of the Ehrlichia species present in Genbank (Figure 1).

13In order to estimate the prevalence and distribution of N. mikurensis in North-West Europe, questing I. ricinus nymphs (~88%) and adults (~12%) were tested using a Q-PCR for the simultaneous detection of N. mikurensis and A. phagocytophilum. In all 12 study-areas in the Netherlands, N. mikurensis was detected with a prevalence varying from 1% to 16% (Table 4). Ticks from one study area (Duin en Kruidberg) were tested in two consecutive years. In 2009, 16% of the questing nymphs and adults I. ricinus were infected with N. mikurensis. The prevalence of N. mikurensis in questing ticks decreased to 8% in 2010. Neoehrlichia mikurensis-positive I. ricinus ticks were found in two out of three regions in Belgium. A fraction of the ticks from the N. mikurensis-negative area (Brussels) were positive for A. phagocytophilum, which was comparable to other regions (data not shown), indicating that the processing and testing of ticks from this area was not affecting the outcome of the results. The results for the A. phagocytophilum will be published elsewhere. To determine whether N. mikurensis is present in the UK, 338 I. ricinus and 63 D. reticularis ticks from a previous study were tested. These ticks were collected at 7 dispersed study areas in the UK and were partially caught by blanket dragging and removed from wildlife, pets and humans. Anaplasma phagocytophilum, but not N. mikurensis, was detected in these tick lysates.

14Transovarial (vertical) transmission has been implicated for Rickettsia and Anaplasma, but not for Ehrlichia species. Whether N. mikurensis is transmitted transovarially in I. ricinus has not been investigated so far. The prevalence of N. mikurensis was determined in 55 pools of 5 questing I. ricinus larvae from Vrouwenpolder, where nymphal and adult ticks were found to be positive for N. mikurensis (Table 4). None of the 55 pools were N. mikurensis-positive (Table 5). Some of the pools were positive for A. phagocytophilum, approving the used methodology. The prevalence of N. mikurensis in questing I. ricinus nymphs was ~7%, whereas the prevalence in adult ticks was ~ 11% (Table 5). No significant differences were observed in the prevalence between questing male and female I. ricinus ticks. To investigate the role of other tick species in the transmission of N. mikurensis: Dermacentor reticulatusI. hexagonus and I. arboricola were analysed for the presence of N. mikurensis (in the multiplex real-time PCR). None were found positive (Table 6). Again, some were found positive for the A. phagocytophilum msp2 gene (data not shown), indicating that there is no significant inhibition within these samples.

15To investigate the possible mammalian hosts for N. mikurensis, 79 spleen samples of different wild small mammals were tested by (nested)-PCR for the presence of gltA and GroEL (Table 7). PCR-positive samples were sequenced to confirm the presence of N. mikurensis. Both the GroEL and gltA sequences isolated from spleen were identical to the N. mikurensis sequences found in the questing ticks in the Netherlands (Figure 1). Spleen samples from Apodemus sylvaticus, Microtus arvalis and Myodes glareolus were N. mikurensis-positive. After the spleen was found positive, other organs (kidney, liver and brain) were also tested for N. mikurensis. All the tested organs were positive.

16Whether other mammals in the Netherlands are reservoir hosts is difficult to address, due to the protective status of these animals. An animal can be considered a potential reservoir host when the prevalence of N. mikurensis in ticks feeding on this animal is significantly higher than the prevalence in questing ticks. This is for example the case for Anaplasma phagocytophilum. I. ricinus feeding on red deer (Cervus elaphus), European mouflon (Ovis orientalis musimon), wild boar (Sus scrofa) and sheep (Ovis aries) were tested by multiplex real-time PCR. The prevalence of N. mikurensis in feeding ticks was comparable to the prevalence in questing ticks (Table 8).

17Recently, six human and one canine case of N. mikurensis infection were reported in different locations in Europe. These reports advocate a re-assessment of the occurrence of this microorganism in questing ticks. Schouls and colleagues described an Ehrlichia-like organism (Ehrlichia schotti) in Dutch ticks. In our study, the three genetic markers 16S rDNA, gltA and GroEL of Ehrlichia schotti-positive field isolates turned out to be similar and identical to DNA sequences available from N. mikurensis. Thus, Ehrlichia schotti and N. mikurensis are most likely one and the same species. Previous findings on E. schotti can be interpreted as findings on N. mikurensis. Thus, N. mikurensis has already been present in the Netherlands in 1999. Furthermore, 11% of 289 engorged I. ricinus removed from humans were N. mikurensis-positive, indicating that the Dutch population is being exposed to ticks infected with N. mikurensis. Remarkably, human and animal cases of N. mikurensis infection in the Netherlands have not yet been described.

18The development of a Q-PCR specific for N. mikurensis allowed us to test significant numbers of ticks without having to perform the labour-intensive Reverse-Line blotting. These analyses showed that the N. mikurensis is present in vegetation ticks throughout the Netherlands and Belgium. No N. mikurensis-positive ticks were found in one location in Belgium. One possible explanation is that this location in the Brussels-area is exceptional due to its reduced fauna and flora caused by human interference. This forest in the Brussels-area is also highly fragmented because of a railroad and several major motorways that run through the forest. Several parts of it can be ecologically considered ‘islands’, which could -through isolation of mammal and tick populations- explain the absence of the pathogen in this forest. More ticks of this unique area need to be tested in order to address this hypothesis. Neoehrlichia mikurensis was also not detected in ticks from the UK. This could indicate that these species have not (yet) been established on this island.

19The overall prevalence of N. mikurensis in questing nymphs and adults is approximately 7%. From the public health point of view, it indicates that a significant proportion of people contracting a tick bite are exposed to N. mikurensis. Transmission of the N. mikurensis in ticks appears to occur horizontally rather than vertically. None of the tested larvae were found positive, even though the prevalence of nymphs is approximately 7% and 11% for adults. Other tick species, with more restricted host preference than I. ricinus, were also tested for the presence of N. mikurensis. Dermacentor reticulatus, I. hexagonus and questing I. arboricola were found negative. The data indicate that these tick species probably play insignificant roles in the transmission of N. mikurensis. In contrast, I. ricinus can be considered as its main vector in the Netherlands and Belgium.

20A potential group of reservoir hosts for N. mikurensis are wild rodents. Indeed, spleen samples and other organs (kidney, liver and brain) of some rodent species turned out to be N. mikurensis-positive, which indicates a systemic infection of these rodents with N. mikurensis. The N. mikurensis isolates from ticks and wild rodents (Table 7) were genetically identical, indicating that rodents are potential reservoir hosts. However, the reservoir potential of rodents can only be by xenodiagnosis or experimental infection. The prevalence of N. mikurensis in I. ricinus ticks feeding on red deer, European mouflon, wild boar and sheep were comparable to the prevalence in questing ticks. From these prevalence data, it was not possible to infer the role of these animals in the transmission of N. mikurensis. However, it is clear that these animals are being exposed to the N. mikurensis through tick bites. Further experiments are necessary to determine whether there are other mammalian reservoirs than wild rodents.

21Although human infection of the N. mikurensis has not been reported in the Netherlands, it is unclear to what extent N. mikurensis poses risks to public health. The symptoms described in all of the N. mikurensis infection cases were generally non-specific and usually seen in any other ordinary inflammatory reaction. What’s more, most of the Ehrlichia infections are known to be either asymptomatic or mild, self-limiting diseases. In other words, infection can occur without causing disease. So far, diagnosis has relied only on PCR amplification of the N. mikurensis. The lack of serological tests makes diagnosis particularly difficult. Against these backdrops, the actual incidence of human infection with Ehrlichia is likely to be much higher than currently reported in Europe. Thorough surveillance and improvement of diagnostic tools will probably increase the number of identified human cases, and consequently provide more insight in the public health relevance of N. mikurensis.

22The authors declare that they have no competing interests.

23CR, EJS, WT, PH, DH, JK and JM organized and participated in the fieldwork for the collection and determination of wildlife samples and ectoparasites. SJ and MF developed the laboratory tests for the simultaneous detection of N. mikurensis and A. phagocytophilum. SJ, MF and MF carried out the laboratory preparation and testing of all the samples. SJ performed the sequencing and sequence analyses. SJ and HS drafted the manuscript and wrote the final version. HS organized and planned the study. All authors read and approved the final manuscript.