# Next generation sequencing uncovers unexpected bacterial pathogens in ticks in western Europe

Muriel Vayssier Taussat, Sara Moutailler, Lorraine Michelet, Elodie Devillers, Sarah Bonnet, Justine Cheval, Charles Hébert, Marc Eloit

► To cite this version:

Muriel Vayssier Taussat, Sara Moutailler, Lorraine Michelet, Elodie Devillers, Sarah Bonnet, et al. Next generation sequencing uncovers unexpected bacterial pathogens in ticks in western Europe. PLoS ONE, 2013, 8 (11), 10.1371/journal.pone.0081439. [DOI: 10.1371/journal.pone.0081439](https://doi.org/10.1371/journal.pone.0081439)

HAL Id: hal-02649277

https://hal.inrae.fr/hal-02649277

Submitted on 29 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Next Generation Sequencing Uncovers Unexpected Bacterial Pathogens in Ticks in Western Europe

Muriel Vayssier-Taussat¹, Sara Moutailler¹, Lorraine Michelet¹, Elodie Devillers¹, Sarah Bonnet¹, Justine Cheval², Charles Hébert², Marc Eloit²,³,⁴*

1 USC Bipar, INRA, Anses, Maisons-Alfort, France, 2 PathoQuest SAS, Paris, France, 3 Ecole Nationale Vétérinaire d'Alfort, UMR 1161 Virologie ENVA, INRA, ANSES, Maisons-Alfort, France, 4 Institut Pasteur, Laboratory of Pathogen Discovery, Paris, France

## Abstract

**Background and Aims:** Ticks are highly susceptible to global environmental and socio-economical changes. Several tick-borne pathogens have been reported in new geographical regions while new species, strains or genetic variants of tick-borne microorganisms are continually being detected. However, tick-borne pathogens are still poorly understood, and it is estimated that half of all human tick-borne disease has an unknown origin. Therefore in order to prevent these diseases, more effort is required to identify unknown or unexpected tick-borne pathogens. *Ixodes ricinus* is the vector for a broad range of bacterial pathogens and the most prevalent tick in Europe. The aim of the present study was to evaluate the capability of Next Generation Sequencing (NGS) to extend the inventory of pathogenic bacteria carried by this species of tick in France.

**Methods:** RNA and DNA were extracted from 1450 *l. ricinus* questing nymphs collected by flagging in Alsace, France. RNA was pooled and used for NGS. Following *de novo* assembly, bacterial contigs were assigned to the closest known taxonomy. DNA was used for real time PCR to confirm taxonomic species assignment of NGS-derived contigs for the doubtful cases, and for determination of prevalence.

**Results:** We have generated a global in-depth picture of tick-borne bacteria. We identified RNA from the main pathogenic bacterial species known to be transmitted by *l. ricinus*. In addition we also identified unanticipated bacterial species for which we have estimated the prevalence within those ticks inhabiting the studied areas.

**Conclusions:** The data obtained from this study has proven that NGS has an enormous potential to detect the unexpected and provides the means to monitor pathogen occurrence.

**Citation:** Vayssier-Taussat M, Moutailler S, Michelet L, Devillers E, Bonnet S, et al. (2013) Next Generation Sequencing Uncovers Unexpected Bacterial Pathogens in Ticks in Western Europe. PLoS ONE 8(11): e81439. doi:10.1371/journal.pone.0081439

**Editor:** Brian Stevenson, University of Kentucky College of Medicine, United States of America

**Received** September 2, 2013; **Accepted** October 22, 2013; **Published** November 27, 2013

**Copyright:** © 2013 Vayssier-Taussat et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This study was partially funded by the EU grant FP7-261504 EDENext and is catalogued by the EDENext Steering Committee as EDENext 181 (http://www.edenext.eu). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** CH and JC are employees and ME is the chairman of PathoQuest, a spin-off of Institut Pasteur. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

* E-mail: marc.eloit@pasteur.fr

## Introduction

Ticks are highly susceptible to global environmental and socio-economical changes [1], such that the incidence of tick-borne disease is increasing worldwide [2], and that some tick-borne pathogens have been reported in new geographical regions [3,4]. *Ixodes ricinus* is the most common European tick species, and is the vector for a broad range of bacterial pathogens [2,5].

Due to recent advances in molecular biology, new species, strains or genetic variants of microorganisms are continually being detected in ticks, resulting in an ever-increasing list of pathogens capable of infecting livestock and humans. In Europeans, the most prevalent tick-borne human disease is Lyme borreliosis, caused by any of four different *Borrelia* species; *B. burgdorferi* sensu stricto, *B. afzelii*, *B. garini*, and *B. spielmanii* [6], with over 50 000 new cases each year in Europe [7]. Lyme disease has a high rate of diagnosis, mainly due to its distinctive symptoms and the existence of effective serological tests. Although, in addition to these typical

cases, patients with tick bites can also be infected by many other bacteria including *Anaplasma phagocytophilum*, *Rickettsia* spp., “Candidatus” Neoehrlichia mikurensis”, *Ehrlichia* spp., *Bartonella* spp. and *Francisella* spp. Some of these species have been identified in ticks decades prior to their association with human disease, whereas others have only been discovered very recently. Therefore we can confidently assume that many tick-associated pathogens are still unknown [3,4,8]. In cases of co-infection, symptoms induced by these microorganisms are often non-specific and can be aggravated by pre-existing or concurrent diseases [9] making diagnosis exceedingly difficult. It is estimated that up to half of human tick-borne disease has an unknown microbiological origin [10]. Consequently, in order to identify novel or unexpected tick-borne pathogens in their natural hosts, concerted efforts are required to develop efficient detection methods.

Up until now, most studies aiming to identify tick-borne pathogens in their vectors have only been able to assess minimal numbers of pathogens at a time, which is partly due to
technological limitations. Indeed, complete screens of microorganisms in natural host populations (both reservoirs and vectors) were out of reach using standard laboratory procedures. Identification of microorganisms in biological samples has been dominated by the use of culture-dependent methods and by conventional molecular approaches. These methodologies suffer from major limitations: (i) isolation in cell culture or in synthetic medium is difficult for many microorganisms and (ii) molecular identification was mostly based on the use of specific primers combined with real-time PCR, which are only capable of detecting a selected number of species. However, the recent rapid development of Next Generation Sequencing (NGS) methods combined with bioinformatics has revolutionized the research field of epidemiology and diagnosis of infectious diseases, and has effectively overcome those limitations. NGS could be effective tools in diagnosing pathogens for medical and veterinary purposes [11–18] and these techniques has recently been successfully used to identify the bacterial communities associated with *I. ricinus* ticks [19,20]. In these studies, the authors used the amplification and the sequencing of the hyper-variable region V6 of the 16S rRNA encoding genes. If this gene is well adapted to describe bacterial communities, the analysis of its sequence does not allow identifying bacteria at the species level. Ticks carry many bacterial species, including endogenous ones (symbionts and commensals). In ticks, symbionts are often phylogenetically close to pathogens and 16S rRNA diversity does not allow, for instance, distinguishing symbiotic *Rickettsia* species from pathogenic ones. Besides this problem, within a bacteria genus, many species exist. This is the case for *Bartonella* sp., for which some species are known to be pathogenic for humans (or animals) and other have never been associated to diseases. Thus identifying them at the species level is very important. This is also the case for the bacteria belonging to the complex of *Borrelia burgdorferi* sensu lato, responsible of Lyme diseases, for which the analysis of 16S rRNA diversity is not efficient to identify them at the species level.

In this context, we undertook a study to produce an inventory of predicted and/or unexpected pathogenic bacteria carried by ticks, in a geographical region with abundant ticks and a concomitant high risk of disease transmission. To overcome the limitation of deep sequencing of 16S RNA encoding genes, we sequenced the whole transcriptome of ticks and their resident microbiota.

# Materials and Methods

## Tick collection and extract preparation

No specific permits were required for the described field studies. These field studies did not involve endangered or protected species. A total of 1450 questing nymphs and 62 females of the *I. ricinus* tick were collected by flagging from three wooded habitats in Eastern France (Alsace Department, [Murbach (47°55'05"N, 7°8'46"E), Hohbuhl (48°27'33"N, 7°17'22"E) and Wasselonne (48°38'09"N, 7°21'45"E)]). All collected ticks were washed as previously described [21], nymphs were pooled into groups of 25 individuals (58 pools in total) and adults into groups of 2 individuals (31 pools in total) and crushed in 300 µl of Dulbecco's MEM (DMEM) medium supplemented by 10% fetal bovine serum. A pool of 15 *I. ricinus* nymphs from our pathogen-free colony was equivalently treated and used as a reference. This control colony originated from female ticks collected in Murbach (Alsace, France) and was bred as previously described [22]. The pathogen-free colony was obtained as follows; females from the field were engorged on rabbits and allowed to lay eggs. Female tick DNA samples were extracted and different PCRs were performed to test for the presence of *Borrelia* spp., *Bartonella* spp., *Anaplasma* spp.,

*Rickettsia* spp., *Francisella* spp. and *Coxiella* spp. Only larvae from “pathogen-free” female ticks were conserved and maintained in our colony.

## DNA and RNA extraction

Crushed tick pools were divided into two equal samples. Total RNA, used to identify viable and replicating microorganisms, and total DNA, for both pathogen taxonomic confirmation and prevalence studies, were then extracted from the nymphs using the Nucleospin RNA II kit (Macherey-Nagel, Duren, Germany), and the Wizard genomic DNA purification kit (Promega, Madison, WI, USA) following manufacturers’ instructions.

## High throughput sequencing and data analysis

High throughput sequencing was only performed on nymphs, and all RNA samples were pooled before sequencing. RNA was first retro-transcribed to cDNA, then randomly amplified using the Multiple Displacement Amplification (MDA) protocol with phi29 polymerase and random hexamers as described [17]. Library preparations and sequencing with an Illumina HiSeq2000 sequencer were outsourced to DNAVision (Belgium). Wild and pathogen-free samples were sequenced to a depth of 100 Mio and 62 Mio paired-end reads of 101 bp respectively. Raw sequence reads were trimmed according to their quality score. At the time of analysis, there was no published reference genome for *I. ricinus* available in public databases, so sequences from the ticks themselves, or from symbiotic or commensal bacteria naturally found in ticks, were removed from the analysis by subtracting sequences derived from the pathogen-free reference sample using the SOAP2 aligner tool. Finally, *de novo* assembly was performed on all remaining reads (7.7 Mio), producing 174 841 contigs. Taxonomic assignation of these contigs was achieved via successive sequence alignment using nucleotide and protein databases and the BLAST algorithm. Contigs were assigned the closest known taxonomy according to their identity percentage, and distant alignments were not considered. Of the assigned reads, 12.22% of the cDNA derived sequences were of bacterial origin.

## Specific Real-Time PCR for taxonomic species assignment of NGS-derived contigs confirmation

Real-time PCR on DNA samples was performed to confirm taxonomic species assignment of NGS-derived contigs for those cases in doubt: *Francisella tularensis*, *Coxiella burnettii*, *Ehrlichia canis*, *Rickettsia felis* and *R. canadensis*. All DNA samples were pooled as for RNA analysis. Primers and probes used for amplification are described in Table 1.

DNA pre-amplification was performed using the TaqMan PreAmp Master Mix kit (Applied Biosystems). Primers were pooled at a final concentration of 0.2X with equal primer concentrations (20 µM each). The reaction was performed in a final volume of 5 µl containing 2.5 µl TaqMan PreAmp Master Mix (2X), 1.2 µl of the pooled primer mix (0.2X) and 1.3 µl of DNA, with one cycle at 95°C for 10 min, 14 cycles at 95°C for 15 sec and 4 min at 60°C. After the pre-amplification PCR, products were diluted 1:20 and stored at -20°C prior to use in qPCR. All fluorogenic probes were synthesized with a 6-carboxy-fluorescein (FAM) reporter molecule attached to the 5' end and a Black Hole Quencher 1 (BHQ1) attached to the 3' end. Real-time Taqman PCR assays were performed in a final volume of 12 µl using the LightCycler® 480 Probe Master mix (Roche Applied Science, Germany) at a 1X final concentration, with primers and probes at 200 nM and pre-amplified DNA. Thermal cycling
<table><caption>**Table 1.** Primers and probes used for pathogen confirmation or prevalence estimation in ticks.</caption><thead><tr><th>Pathogen</th><th>Target</th><th>Primers and Probe Sequence (5' → 3')</th><th>Amplicon size (bp)</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="3">Borrelia miyamotoi</td><td rowspan="3">glpQ</td><td>CACGACCCAGAAATTGACACA</td><td rowspan="3">94</td><td rowspan="3">This study</td></tr><tr><td>GTGTGAAGTCAGTGGCGTAAT</td></tr><tr><td>TCGTCCGTTTTCTCTAGCTCGATTGGG</td></tr><tr><td rowspan="3">Coxiella burnetti</td><td rowspan="3">idc</td><td>AGGCCCGTCCGTTATTTTACG</td><td rowspan="3">74</td><td rowspan="3">Adapted from [47]</td></tr><tr><td>CGGAAAATCACCATATTCACCTT</td></tr><tr><td>TTCAGGCGTTTTGACCGGGCTTGGC</td></tr><tr><td rowspan="3">Ehrlichia canis</td><td rowspan="3">dsb</td><td>AATACTTGGTGAGCTTCACTCA</td><td rowspan="3">110</td><td rowspan="3">Adapted from [30]</td></tr><tr><td>GTTGCTTGTAAATGTAGTGCTGC</td></tr><tr><td>AAGTTGCCCAAGCAGCACTAGCTGTAC</td></tr><tr><td rowspan="3">Francisella tularensis</td><td rowspan="3">tul4</td><td>ACCCACAAGGAAGTGTAAGATTA</td><td rowspan="3">76</td><td rowspan="3">Adapted from [48]</td></tr><tr><td>GTAATTGGGAAGCTTGTATCATG</td></tr><tr><td>AATGGCAGGCTCCAGAAGGTTCTAAGT</td></tr><tr><td rowspan="3">Rickettsia felis</td><td rowspan="3">orfB</td><td>CCCTTTTCGTAACGCTTTGCT</td><td rowspan="3">149</td><td rowspan="3">[49]</td></tr><tr><td>GGGCTAAACCAGGGAAACCT</td></tr><tr><td>TGTTCCGGTTTTAACGGCAGATACCCA</td></tr><tr><td rowspan="3">Rickettsia canadensis</td><td rowspan="3">ompA</td><td>AGTGGGACATTGGGTGTTGC</td><td rowspan="3">129</td><td rowspan="3">This study</td></tr><tr><td>GGCAACACCTTGAGCATTATC</td></tr><tr><td>CTCCTGCCTGCGTTATACCATGCCA</td></tr><tr><td rowspan="3">Candidatus Neoehrlichia mikurensis</td><td rowspan="3">groEL</td><td>AGAGACATCATTCGCATTTTGGA</td><td rowspan="3">96</td><td rowspan="3">This study</td></tr><tr><td>TTCCGGTGACCATAAGGCTT</td></tr><tr><td>AGATGCTGTTGGATGTACTGCTGGACC</td></tr></tbody></table>

doi:10.1371/journal.pone.0081439.t001

conditions were as follows: 95°C for 5 min, 45 cycles at 95°C for 10 s and 60°C for 15 s and then 40°C for 10 s.

## Specific Real-Time PCR for prevalence of unexpected pathogens in ticks

Real-time PCR was performed to determine the prevalence of three unexpected tick-borne bacteria DNA in ticks: “Candidatus Neoehrlichia mikurensis”, *Borrelia miyamotoi*, and *R. felis*. We assessed 58 pools of 25 nymphs according to the collection location (corresponding to a total of 1450 nymphs) and 31 pools of 2 females (corresponding to a total of 62 females). The primers for amplification of “Candidatus Neoehrlichia mikurensis” and *Borrelia miyamotoi* were specifically designed for this study and are presented in Table 1. Real-time PCR was conducted as described above.

## Statistical analysis

For the pooled nymphs, we estimated prevalence at the individual level using the exact method of Hauck, assuming perfect sensitivity and specificity of our pathogen detection methods [23]. Hauck noted a one-to-one relationship between individual prevalence levels, $\pi$, and the prevalence of positive pools, P. A point estimate for the prevalence rate can thus be obtained from the pool positive rate by $\pi = 1-(1-P)^{1/k}$ where k is the number of nymphs or adults per pool. Exact 95% confidence intervals were then obtained by assuming a binomial distribution for the number of positive pools [24]. For the adults, as each pool only contained two individuals, we calculated the individual prevalence intervals corresponding to prevalence if 2/2 individuals

from the positive pools were positive and prevalence if 1/2 individuals from the positive pools were positive.

## Results

### Inventory of medically important bacterial RNA content in ticks

Bacterial content was organized into genera or putative species by comparison with the closest homolog (Figure 1). Species were only allocated unambiguous assignment to one specific species when either the NGS-derived contig identity percentage was ≥95% compared to a specific species sequence (and lower when compared to other species); or when NGS-derived contigs matched sequences from one single bacterial species. Following these criteria, all bacterial species unambiguously identified are listed in Table 2.

Global analysis of genera mainly associated with medically important tick-borne bacteria: *Anaplasma* spp., *Borrelia* spp., *Coxiella* spp., *Ehrlichia* spp., *Francisella* spp. and *Rickettsia* spp., indicated that among the alpha-proteobacteria family, *Rickettsia* spp. were the most highly represented. Species identified included *R. helvetica*, a member of the Spotted Fever Group of *Rickettsiae*, predicted to be associated with human illness in many European countries [25]; *R. canadensis* isolated from ticks belonging to *Haemaphysalis* spp. in Canada [26], never previously identified in Europe, and involved, on serological evidence of infection, in Rocky Mountain Spotted Fever-like disease and cerebral vasculitis in US patients; and *R. felis*, responsible for Spotted Fever in humans and whose only known vector is the flea [27]. As the number of reads was low for *R. canadensis* (16) and *R. felis* (180), and because the presence of
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Megan4-generated taxonomy assignment and read numbers based on NCBI reference sequences. Each node is labeled with a taxon, and the number of reads assigned to this taxon.
doi:10.1371/journal.pone.0081439.g001</figcaption></figure>
<table><caption>**Table 2.** List of bacteria unambiguously identified at the species level in ticks, including the number of NGS reads for each species, the corresponding list of identified genes allowing taxonomic affiliation (with% identity) to one specific bacteria, and the closest bacterial species to the unambiguously identified species (with% identity).</caption><thead><tr><th>Unambiguously identified bacterial species (read number)</th><th>Gene hit (Accession number) -% id</th><th>Closest species -% id</th></tr></thead><tbody><tr><td rowspan="2">Anaplasma phagocytophilum (5234)</td><td>16S RNA (JN990105) - 100%</td><td>A. ovis - 99%</td></tr><tr><td>23S RNA (AF416766) - 100%</td><td>A. marginale - 96%</td></tr><tr><td rowspan="7">Bartonella henselae (33)</td><td>dnaB (BH0270) - 100%</td><td>B. tribocorum - 94%</td></tr><tr><td>dnaJ (BH00660) - 100%</td><td>B. quintana - 91%</td></tr><tr><td>trkA (BH05860) - 100%</td><td>B. quintana - 93%</td></tr><tr><td>Phage related protein (BH09270) - 100%</td><td>B. grahamii - 84%</td></tr><tr><td>Hypothetical protein (BH022330) - 99%</td><td>B. tribocorum - 85%</td></tr><tr><td>Transketolase (BH15090) - 99%</td><td>B. quintana - 87%</td></tr><tr><td rowspan="2">Bartonella grahamii (39)</td><td>Phage tail protein encoding gene (ACS5068.1) - 90%</td><td>None</td></tr><tr><td>Adhesin (ACS50528.1) - 75%</td><td>None</td></tr><tr><td>Borrelia afzelii (1057)</td><td>Hypothetical protein (CP000395)</td><td>B. garinii - 96%</td></tr><tr><td rowspan="2">Borrelia garinii (723)</td><td>Hypothetical protein (CP001205) - 100%</td><td>B. afzelii - 98%</td></tr><tr><td>Hypothetical protein (CP003151) -100%</td><td>B. afzelii - 98%</td></tr><tr><td>Borrelia burgdorferi s.s. (475)</td><td>23S RNA (BORG23S) - 100%</td><td>B. garinii - 99%</td></tr><tr><td>Borrelia miyamotoi (2499)</td><td>16S RNA (JF951385) - 100%</td><td>Other Borrelia sp. &lt;98%</td></tr><tr><td>Candidatus "Neoehrlichia mikurensis" (2010)</td><td>16S RNA (JQ359061) - 100%</td><td>A. phagocytophilum &lt;97%</td></tr><tr><td rowspan="2">Ehrlichia canis* (660)</td><td>PhoT family ABC transporter (CP000107) - 75%</td><td>None</td></tr><tr><td>23SrRNA - 96%</td><td>E. chaffensis - 96%</td></tr><tr><td rowspan="3">Rickettsia canadensis (16)</td><td>MraW (CP000409) - 96%</td><td>R. felis - 93%</td></tr><tr><td>TraD (CP000409) - 97%</td><td>R. amblyommii - 85%</td></tr><tr><td>16S RNA methyltransferase: (CP003304) - 97%</td><td>R. montanensis - 88%</td></tr><tr><td rowspan="5">Rickettsia felis (180)</td><td>Plasmid pRF (GQ329881) - 98%</td><td>R. australis - 91%</td></tr><tr><td>Unknown protein (CP000053) - 95%</td><td>R. australis - 90%</td></tr><tr><td>AmpC (CP000053.1) - 99%</td><td>R. amblyoma - 97%</td></tr><tr><td>Sca13 (CP000053) - 96%</td><td>R. slovaca - 91%</td></tr><tr><td>Sca1 (AY355363) - 100%</td><td>R. felis - 91%</td></tr><tr><td rowspan="4">Rickettsia helvetica (2411)</td><td>ompB (JN036412.1) - 100%</td><td>R. asiatica - 94%</td></tr><tr><td>groEL (DQ442911) - 99%</td><td>R. canadensis - 97%</td></tr><tr><td>sca4 (FJ358501.1) - 100%</td><td>R. tamurae - 94%</td></tr><tr><td></td><td></td><td></td></tr></tbody></table>

doi:10.1371/journal.pone.0081439.t002

both pathogens in *I. ricinus* was not expected, we performed real-time PCR on DNA extracts corresponding to a single pool of 1450 nymphs, which confirmed the detection of *R. felis* DNA, but not *R. canadensis* indicating that for this bacteria the RNA-seq results is below the threshold of detection. We can thus consider that this bacteria is not present in our samples

Among the Anaplasmataceae family, we identified sequences belonging to A. phagocytophilum, the agent of human granulocytosis [25] and “*candidatus* Neoehrlichia mikurensis” a new tick-borne pathogen discovered in 2004. The first human case of “*candidatus* N. mikurensis” infection was reported in a Swedish patient in 2010, but has never been linked to human disease in France thus far [28]. The presence of DNA from this bacterium was confirmed by real-time PCR in the pool of 1450 nymphs. In 2012, we had previously identified this bacterium for the first time in France, in bank voles from North Eastern France [29]. We also identified species belonging to the Ehrlichia genus. NGS sequence analysis combined with specific real-time PCR confirmation facilitated the detection of E. canis, the agent of canine

ehrlichiosis, known to be transmitted by the dark tick Rhipicephalus sanguineus [30] (Table 2).

Among the *Bartonella* spp. sequences, the cat-specific *B. henselae* was clearly identified, as well as sequences related to the rodent-specific *B. grahamii* (Table 2). Both of these species are known to be mainly flea-borne zoonotic human pathogens but have also been previously detected in European ticks [31–33].

We also detected RNA sequences related to gamma-proteobacteria belonging to *Coxiella* spp. and *Francisella* spp. For these genera, the low number of reads prevented the identification of these bacteria at the species level. We failed to detect *Coxiella burnetti* or *Francisella tularensis* DNA within ticks by real-time PCR, these two species being the most common tick-borne pathogens of each genera.

As expected for ticks collected in Alsace, a French region with high rates of Lyme disease, we also identified *B. afzelii*, *B. garinii*, and *B. burgdorferi s.s.*, the three human pathogenic species known to circulate in France [34]. Interestingly, we also detected sequences from *B. miyamotoi*, which is responsible for relapsing fever in Russia, and which was recently detected in the USA [35]. The
presence of *B. miyamotoi* DNA was confirmed by real-time PCR. Even though sequences from this bacteria have previously been amplified from *I. ricinus* in Sweden, France and Slovakia [21,36], it has never thus far been associated with disease in Europe.

## Prevalence of unexpected human pathogen DNA in ticks

Considering their potential importance to human health and disease, we estimated the prevalence of “*candidatus* Neoehrlichia mikurensis”, *B. miyamotoi*, and *R. felis* microorganisms. Among the 58 tested pools of nymph DNA (each pool containing 25 nymphs), 19 were positive for “*candidatus* Neoehrlichia mikurensis” and 25 for *B. miyamotoi*. By estimating the individual prevalence of nymphs using the exact method of Hauck, we found a prevalence of 2.2% (1.4–3.3) for *Borrelia miyamotoi* and 1.6% (0.9–2.5) for “*candidatus* Neoehrlichia mikurensis”. For adults, we analyzed 31 pools, each containing two adults. We identified two positive pools for “*candidatus* Neoehrlichia mikurensis” which represented an individual prevalence of between 3.2% and 6.4%. All tested adult pools were negative for *B. miyamotoi*. For *R. felis*, we identified one positive pool of nymphs among the 58 tested (estimated prevalence of 0.09%), but none within the adult pools.

## Discussion

We have used a whole transcriptome analysis of ticks RNA to generate a global picture of tick-borne bacteria and validate an NGS strategy, by identifying bacteria at the RNA level. It has been recently shown that metagenomic 16S rDNA Illumina tags are a powerful alternative to 16S DNA amplicon sequencing[37]. This paper extends this approach by using wide transcriptome analysis, including but not limited to 16S and 23 S RNAs. We have identified bacterial species known to be transmitted by *I. ricinus* and causing disease in France or Europe: *Rickettsia helvetica*; *Anaplasma phagocytophilum*; *Borrelia garinii*, *B. afzelii*, *B. burgdorferi* s.s. and “*candidatus* Neoehrlichia mikurensis”. Considering “*candidatus* Neoehrlichia mikurensis”, has never been identified in ticks from France, we calculated an estimated prevalence of 1.6% in nymphs and from between 3.2 to 6.4% in adults, rates similar to those found elsewhere in Europe [3]. “*candidatus* Neoehrlichia mikurensis” infection has also been detected in humans from northeastern China, and field surveys showed that 1.6% of ticks and 3.8% of rodents collected from patient residences were also infected [38,39]. Although corresponding human or animal disease has not yet been reported in France, the prevalence of positive rodents and ticks is similar to levels found in China. In addition, the infected French bank vole (*Myodes glareolus*) and ticks are located in close proximity to human dwellings, suggesting that “*candidatus* Neoehrlichia mikurensis” is very likely to pose risks to public health in France. However, to date, the diagnosis of “*candidatus* Neoehrlichia mikurensis” infection has solely relied on PCR amplification of bacterial DNA. Therefore, the absence of serological tests combined with physicians’ lack of information relating to these bacteria, makes diagnosis particularly difficult. Thorough surveillance and improved diagnostic tools are required to provide more insight into the relevance of *Candidatus* N. mikurensis to public health.

In addition to the expected species, our study generated evidence of other unanticipated bacterial species. For example, *B. miyamotoi*, transmitted by *Ixodes* spp., and causing recurrent

fever, has not previously been associated with disease in Europe [35]. *B. miyamotoi* was isolated for the first time in Japan in 1995 from *I. persulcatus* ticks as well as from the blood of *Apodemus argentus* mice [40,41]. In Europe, *B. miyamotoi* was detected in ticks from Sweden, Germany, and very recently, from Estonia [36,42,43]. In addition to *A. argentus* mice, it has been shown that other rodents such as white-footed mice (*Peromyscus leucopus*) may serve as host reservoirs for *B. miyamotoi* [44]. In our study *B. miyamotoi* was estimated to occur at a mean prevalence rate of 2.2% in nymphs, a rate close to those obtained from different European countries, which fluctuates from 0.4% in Estonia to 3.5% in Germany [35,36,43,45]. Some symptoms caused by this particular *Borrelia* sp., such as recurring high fever, are very different from the skin rash and joint pain typical of Lyme disease. Identifying the RNA of this bacterium in French *I. ricinus* for the second time strongly suggests that the disease could exist in this country. However, in the absence of specific tests, it is impossible to know whether the disease is indeed present in France (or in Europe), and if it has the potential to cause serious long-term damage, as can result from untreated Lyme disease.

We also detected known pathogens which are either normally transmitted by other tick species, or by vectors other than ticks, such as for the *B. henselae*, *B. grahamii* and *R. felis* microorganisms. *B. henselae* and *R. felis* are typically transmitted from cat to cat via fleas, with human contamination arising from cat or flea bites. Establishing ticks as vectors of *B. henselae* has been fiercely debated for many years, even given the recent demonstration of the competence of *I. ricinus* to host *B. henselae* [46]. To our knowledge, our data present the first characterization of *R. felis* RNA in ticks. However even though we have confirmed the presence of RNA and DNA from *R. felis* in these ticks by real time-PCR with specific primers, the detection of *R. felis* DNA in only one pool of nymphs strongly suggests, if it exists, an anecdotic role for *I. ricinus* in *R. felis* circulation.

Detecting RNA of unexpected pathogens in *Ixodes* ticks indicates that these microorganisms are viable and actively replicating, however, it does not prove that transmission by *I. ricinus* actually occurs. Competence studies of *I. ricinus* for all detected pathogens should be experimentally clarified in the future before reclassifying them as tick-borne pathogens. Although further research is thus required to confirm the role of *I. ricinus* in transmitting these pathogens, the data obtained with NGS has proven that this method has enormous potential to detect the unexpected, and provides an excellent means for surveying pathogen occurrence in disease vectors from particular localities.

## Acknowledgments

We thank the group Tiques et Maladies à Tiques group (TMT) of Réseau Ecologie des Interactions Durables for stimulating discussion and support, the DIM Ile de France, the CoVetLab partner institutes, the ANSES, the Animal Health department of the INRA, the INRA Metaprogramme MEM (Patho-ID) and EU FP7-261504 EDENext for their support.

## Author Contributions

Conceived and designed the experiments: MVT ME. Performed the experiments: LM ED JC CH. Analyzed the data: MVT SM SB JC CH ME. Wrote the paper: MVT JC ME.

## References

1. Githeko AK, Lindsay SW, Confalonieri UE, Patz JA (2000) Climate change and vector-borne diseases: a regional analysis. Bull World Health Organ 78: 1136–1147.

2. Parola P, Raoult D (2001) Ticks and tickborne bacterial diseases in humans: an emerging infectious threat. Clin Infect Dis 32: 897–928.
3. Andersson M, Raberg L (2011) Wild rodents and novel human pathogen candidatus Neoehrlichia mikurensis, Southern Sweden. Emerg Infect Dis 17: 1716-1718.

4. Haselbarth K, Tenter AM, Brade V, Krieger G, Hunfeld KP (2007) First case of human babesiosis in Germany - Clinical presentation and molecular characterisation of the pathogen. Int J Med Microbiol 297: 197-204.

5. Parola P, Raoult D (2001) Tick-borne bacterial diseases emerging in Europe. Clin Microbiol Infect 7: 80-83.

6. Hubalek Z (2009) Epidemiology of lyme borreliosis. Curr Probl Dermatol 37: 31-50.

7. Piesman J, Eisen L (2008) Prevention of tick-borne diseases. Annu Rev Entomol 53: 323-343.

8. Kik M, Nijhof AM, Balk JA, Jongejan F (2011) Babesia sp. EU1 infection in a forest reindeer, The Netherlands. Emerg Infect Dis 17: 936-938.

9. Mayne PJ (2011) Emerging incidence of Lyme borreliosis, babesiosis, bartonellosis, and granulocyticehrlichiosis in Australia. Int J Gen Med 4: 845-852.

10. Feder HM, Jr., Johnson BJ, O'Connell S, Shapiro ED, Steere AC, et al. (2007) A critical appraisal of "chronic Lyme disease". N Engl J Med 357: 1422-1430.

11. Sauvage V, Cheval J, Foulongne V, Gouilh MA, Pariente K, et al. (2011) Identification of the first human Gyrovirus, a virus related to chicken anemia virus. J Virol.

12. Sauvage V, Foulongne V, Cheval J, Ar Gouilh M, Pariente K, et al. (2011) Human polyomavirus related to African green monkey lymphotropic polyomavirus. Emerg Infect Dis 17: 1364-1370.

13. Cox-Foster DL, Conlan S, Holmes EC, Palacios G, Evans JD, et al. (2007) A metagenomic survey of microbes in honey bee colony collapse disorder. Science 318: 283-287.

14. Buckwalter MR, Nga PT, Gouilh MA, Fiette L, Bureau JF, et al. (2011) Identification of a Novel Neuropathogenic Theiler's Murine Encephalomyelitis Virus. J Virol 85: 6893-6905.

15. Nakamura S, Maeda N, Miron IM, Yoh M, Izutsu K, et al. (2008) Metagenomic diagnosis of bacterial infections. Emerg Infect Dis 14: 1784-1786.

16. Nakamura S, Nakaya T, Iida T (2011) Metagenomic analysis of bacterial infections by means of high-throughput DNA sequencing. Exp Biol Med (Maywood) 236: 968-971.

17. Cheval J, Sauvage V, Frangeul L, Dacheux L, Guigon G, et al. (2011) Evaluation of high-throughput sequencing for identifying known and unknown viruses in biological samples. J Clin Microbiol 49: 3268-3275.

18. Hoffmann B, Scheuch M, Hoper D, Jungblut R, Holsteg M, et al. (2012) Novel orthobunyavirus in cattle, europe, 2011. Emerg Infect Dis 18: 469-472.

19. Carpi G, Cagnacci F, Wittekindt NE, Zhao F, Qi J, et al. (2011) Metagenomic profile of the bacterial communities associated with Ixodes ricinus ticks. PLoS One 6: e25604.

20. Nakao R, Abe T, Nijhof AM, Yamamoto S, Jongejan F, et al. (2013) A novel approach, based on BLSOMs (Batch Learning Self-Organizing Maps), to the microbiome analysis of ticks. ISME J 7: 1003-1015.

21. Reis C, Cote M, Paul RE, Bonnet S (2011) Questing ticks in suburban forest are infected by at least six tick-borne pathogens. Vector Borne Zoonotic Dis 11: 907-916.

22. Reis C, Cote M, Le Rhun D, Lecuelle B, Levin MI, et al. (2011) Vector competence of the tick Ixodes ricinus for transmission of Bartonella birtlesii. PLoS Negl Trop Dis 5: e1186.

23. Hauck WW (1991) Confidence intervals for seroprevalence determined from pooled sera. Ann Epidemiol 1: 277-281.

24. Cowling DW, Gardner IA, Johnson WO (1999) Comparison of methods for estimation of individual-level prevalence based on pooled samples. Prev Vet Med 39: 211-225.

25. Socolovschi C, Mediannikov O, Raoult D, Parola P (2009) Update on tick-borne bacterial diseases in Europe. Parasite 16: 259-273.

26. McKiel JA, Bell EJ, Lackman DB (1967) Rickettsia canada: a new member of the typhus group of rickettsiae isolated from Haemaphysalis leporispalustris ticks in Canada. Can J Microbiol 13: 503-510.

27. Brouqui P, Parola P, Fournier PE, Raoult D (2007) Spotted fever rickettsioses in southern and eastern Europe. FEMS Immunol Med Microbiol 49: 2-12.

28. Welinder-Olsson C, Kjellin E, Vaht K, Jacobsson S, Wenneras C (2010) First case of human "Candidatus Neoehrlichia mikurensis" infection in a febrile patient with chronic lymphocytic leukemia. J Clin Microbiol 48: 1956-1959.

29. Vayssier-Taussat M, Le Rhun D, Buffet JP, Maaoui N, Galan M, et al. (2012) Candidatus Neoehrlichia mikurensis in Bank Voles, France. Emerg Infect Dis 18: 2063-2065.

30. Doyle CK, Labruna MB, Breitschwerdt EB, Tang YW, Corstvet RE, et al. (2005) Detection of medically important Ehrlichia by quantitative multicolor TaqMan real-time polymerase chain reaction of the dsb gene. J Mol Diagn 7: 504-510.

31. Chomel BB, Boulouis HJ, Breitschwerdt EB, Kasten RW, Vayssier-Taussat M, et al. (2009) Ecological fitness and strategies of adaptation of Bartonella species to their hosts and vectors. Vet Res 40: 29.

32. Janecek E, Mietze A, Goethe R, Schnieder T, Strube C (2012) Bartonella spp. infection rate and B. grahamii in ticks. Emerg Infect Dis 18: 1689-1690.

33. Dietrich F, Schmidgen T, Maggi RG, Richter D, Matuschka FR, et al. (2010) Prevalence of Bartonella henselae and Borrelia burgdorferi sensu lato DNA in ixodes ricinus ticks in Europe. Appl Environ Microbiol 76: 1395-1398.

34. Marsot M, Sigaud M, Chapuis JL, Ferquel E, Cornet M, et al. (2011) Introduced Siberian chipmunks (Tamias sibiricus barberi) harbor more-diverse Borrelia burgdorferi sensu lato genospecies than native bank voles (Myodes glareolus). Appl Environ Microbiol 77: 5716-5721.

35. Platonov AE, Karan LS, Kolyasnikova NM, Makhneva NA, Toporkova MG, et al. (2011) Humans infected with relapsing fever spirochete Borrelia miyamotoi, Russia. Emerg Infect Dis 17: 1816-1823.

36. Geller J, Nazarova L, Katargina O, Jarvekulg L, Fomenko N, et al. (2012) Detection and genetic characterization of relapsing fever spirochete Borrelia miyamotoi in Estonian ticks. PLoS One 7: e51914.

37. Logares R, Sunagawa S, Salazar G, Cornejo-Castillo FM, Ferrera I, et al. (2013) Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities. Environ Microbiol.

38. Li H, Jiang J, Tang F, Sun Y, Li Z, et al. (2013) Wide distribution and genetic diversity of "Candidatus Neoehrlichia mikurensis" in rodents from China. Appl Environ Microbiol 79: 1024-1027.

39. Li H, Jiang JF, Liu W, Zheng YC, Huo QB, et al. (2012) Human infection with Candidatus Neoehrlichia mikurensis, China. Emerg Infect Dis 18: 1636-1639.

40. Fukunaga M, Koreki Y (1995) The flagellin gene of Borrelia miyamotoi sp. nov. and its phylogenetic relationship among Borrelia species. FEMS Microbiol Lett 134: 255-258.

41. Fukunaga M, Takahashi Y, Tsuruta Y, Matsushita O, Ralph D, et al. (1995) Genetic and phenotypic analysis of Borrelia miyamotoi sp. nov., isolated from the ixodid tick Ixodes persulcatus, the vector for Lyme disease in Japan. Int J Syst Bacteriol 45: 804-810.

42. Fraenkel CJ, Garpmo U, Berglund J (2002) Determination of novel Borrelia genospecies in Swedish Ixodes ricinus ticks. J Clin Microbiol 40: 3308-3312.

43. Richter D, Schlee DB, Matuschka FR (2003) Relapsing fever-like spirochetes infecting European vector tick of Lyme disease agent. Emerg Infect Dis 9: 697-701.

44. Scoles GA, Papero M, Beati L, Fish D (2001) A relapsing fever group spirochete transmitted by Ixodes scapularis ticks. Vector Borne Zoonotic Dis 1: 21-34.

45. Sarksyan DS, Platonov AE, Karan LS, Malinin IE, Khalitova LI, et al. (2012) [Clinical presentation of "new" tick-borne borreliosis caused by Borrelia miyamotoi]. Ter Arkh 84: 34-41.

46. Cotte V, Bonnet S, Le Rhun D, Le Naour E, Chauvin A, et al. (2008) Transmission of Bartonella henselae by Ixodes ricinus. Emerg Infect Dis 14: 1074-1080.

47. Klee SR, Tyczka J, Ellerbrok H, Franz T, Linke S, et al. (2006) Highly sensitive real-time PCR for specific detection and quantification of Coxiella burnetii. BMC Microbiol 6: 2.

48. Versage JL, Severin DD, Chu MC, Petersen JM (2003) Development of a multitarget real-time TaqMan PCR assay for enhanced detection of Francisella tularensis in complex specimens. J Clin Microbiol 41: 5492-5499.

49. Socolovschi C, Pages F, Raoult D (2012) Rickettsia felis in Aedes albopictus mosquitoes, Libreville, Gabon. Emerg Infect Dis 18: 1687-1689.