# Broad range molecular detection methods identify only Borrelia spp. in erythema migrans biopsies and blood of tick-bitten patients

Philippe Pérot, Laura Tondeur, Sara Moutailler, Delphine Chrétien, Nicole Corre-Catelin, Muriel Vayssier-Taussat, Marc Eloit, Catherine Chirouze, Céline Cazorla

<figure><img src="image_1.png" /></figure>

## To cite this version:

Philippe Pérot, Laura Tondeur, Sara Moutailler, Delphine Chrétien, Nicole Corre-Catelin, et al. Broad range molecular detection methods identify only Borrelia spp. in erythema migrans biopsies and blood of tick-bitten patients. One Health, 2024, 19, pp.100886. 10.1016/j.onehlt.2024.100886. hal-04693283

HAL Id: hal-04693283

https://hal.inrae.fr/hal-04693283v1

Submitted on 12 Sep 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Broad range molecular detection methods identify only *Borrelia* spp. in erythema migrans biopsies and blood of tick-bitten patients

Philippe Pérot<sup>a,b,*</sup>, Laura Tondeur<sup>c</sup>, Sara Moutailler<sup>d</sup>, Delphine Chrétien<sup>a,b</sup>,
Nicole Corre-Catelin<sup>e</sup>, Muriel Vayssier-Taussat<sup>f</sup>, Marc Eloit<sup>a,b,g</sup>, Catherine Chirouze<sup>h,i</sup>,
Céline Cazorla<sup>j</sup>, On behalf of the OHTICKS Consortium

## ARTICLE INFO

### Keywords:

Ticks

Lyme disease

Erythema migrans

Metagenomics

RT-PCR

Pathogen discovery

## ABSTRACT

In this multicenter study conducted in France, we challenged the hypothesis of the transmission of pathogens other than *Borrelia* spp. in 22 patients developing erythema migrans following a tick bite. Using a combination of high-throughput microfluidic PCRs and agnostic metagenomics on skin biopsies and blood samples, no micro-organisms other than *Borrelia* spp. was found.

## 1. Introduction

Lyme borreliosis (LB) is the most common tick-borne disease in Europe, is caused by different genospecies of *Borrelia burgdorferi* sensu lato (Bbsl) complex [1,2] and its early clinical diagnosis is based on the presence of cutaneous erythema migrans (EM) [3]. *Borrelia burgdorferi* sensu lato bacteria, the causative agents of Lyme borreliosis, are transmitted in Europe mainly by the tick species *Ixodes ricinus*. Other tick species could be involved in maintaining the transmission cycle of this bacteria in wildlife but are rarely involved in the transmission to humans [4,5]. *Ixodes ricinus* can also transmit other infectious agents pathogenic for humans: bacteria such as *Anaplasma phagocytophilum*, virus such as Tick-borne Encephalitis, or parasites such as *Babesia* [6–10]. In western Europe, the estimated annual incidence of LB is 22/100000 inhabitants with a wide variation depending on the country [11]. However, incomplete surveillance, missed diagnoses and the use of insensitive

laboratory tests suggest significant underreporting. In 2022, approximately 63,000 cases of LD were reported to the CDC in the USA. Recent estimates using different methods suggest that around 476,000 people could be diagnosed with Lyme disease each year in the United States [12].

A minority of patients adequately treated by antibiotics against Lyme borreliosis complains about polymorphic, non-specific (asthenia, fever, myalgia) but persistent symptoms [13]. One hypothesis to explain this clinic persistency is an infection by unforeseen or novel pathogenic microorganisms transmitted by the tick bite. In this context, we proposed a multi-disciplinary project, OHTICKS, bringing together veterinarians, physicians, scientists and sociologists, to better characterize the tick-borne diseases present in humans and domestic animals following a tick bite, in a One Health approach. Four French university hospitals were involved in the OHTICKS project (Besançon, Saint-Etienne, Saint-Antoine (AP-HP) and Garches (AP-HP)). They were chosen either

<footnote><sup>a</sup> Pathogen Discovery Laboratory, Institut Pasteur, Université Paris Cité, F-75015 Paris, France</footnote>

<footnote><sup>b</sup> Institut Pasteur, Université Paris Cité, The WOAH (OIE) Collaborating Center for the Detection and Identification in Humans of Emerging Animal Pathogens, F-75015 Paris, France</footnote>

<footnote><sup>c</sup> Emerging Infectious Diseases Unit, Institut Pasteur, Paris, France</footnote>

<footnote><sup>d</sup> ANSES, INRAE, Ecole Nationale Vétérinaire d'Alfort, UMR BIPAR, Laboratoire de Santé Animale, Maisons-Alfort, F-94700, France</footnote>

<footnote><sup>e</sup> Institut Pasteur, Université Paris Cité, Clinical Investigation and Access to Research Bioresources (ICAReB) Platform, Paris, France</footnote>

<footnote><sup>f</sup> Animal Health Division, INRAE, 37380 Nouzilly, France</footnote>

<footnote><sup>g</sup> Ecole Nationale Vétérinaire d'Alfort, F-94700 Maisons-Alfort, France</footnote>

<footnote><sup>h</sup> Chrono-Environnement UMR6249, CNRS, Université Bourgogne Franche-Comté, F-25000 Besançon, France</footnote>

<footnote><sup>i</sup> Infectious Disease Department, University Hospital of Besançon, F-25000, France</footnote>

<footnote><sup>j</sup> Infectious Disease department, University Hospital of Saint Etienne, Avenue Albert-Raimond, 42055, Saint Etienne Cedex 02, France</footnote>

<footnote>* Corresponding author at: Institut Pasteur, Pathogen Discovery Laboratory, 25-28 rue du Dr. Roux, 75015 Paris, France.
E-mail address: philippe.perot@pasteur.fr (P. Pérot).</footnote>
because they were located in areas with a high incidence of Lyme disease, or because they had particular expertise in this field.

The main objective of OHTICKS consortium was to identify neglected, unsuspected or new micro-organisms in patients bitten by ticks, using two complementary approaches: (i) high-throughput microfluidic real-time PCR targeting 41 major tick-borne bacteria and parasites known to circulate in Europe; (ii) random deep sequencing (NGS) expanding the search to unexpected or novel tick-borne pathogens (bacteria, parasites and viruses). In this study, we used EM, which is specific enough for a clinical diagnosis of Lyme disease, as a marker of transmission of *Borrelia* spp. or other pathogens. In these patients, small in number but well characterized, we searched for *Borrelia* spp. and other infectious agents in one cutaneous biopsy, plus plasma and blood pellet sampled at Day 0 (concomitant with the biopsy) and Day 14, using these two methods.

## 2. Materials and methods

Twenty-six patients with suspected EM were included in the OHTICKS protocol. The diagnosis was made by infectiologists of the consortium, and the photos of all EM were blindly reviewed by an infectiologist from outside the inclusion centers. The lesions were classified as: "suggestive of EM", defined as skin lesion greater than 5 cm and centrifugal evolution with peripheral border darker than the center of the lesion (N = 12/26, 46 %); "not very typical", defined as skin lesion greater than 5 cm and homogeneous erythema (N = 11/26, 42 %); and "really not typical", defined as skin lesion less than 5 cm and homogeneous erythema (N = 3/26, 11 %, excluded from the analysis). The "really not typical" EM patients were initially enrolled in the study because they had reported a tick bite and had a little redness at the site of the bite. In conditions of arthropod bite, differential diagnoses were: hypersensitivity reactions at the site of the bite (non-infectious) or infections with pyogenes bacteria following bites; outside the context of the bite: fixed pigmented erythema, dermatophytosis, granuloma annulare and Morpheus patch. One skin biopsy was eventually performed from the active EM area of each of the 22 patients (N = 12/22 "suggestive", N = 10/22 "not very typical"). Among these 22 subjects, 13 (59 %) subjects were male and the median [IQR] age was 49 [40–64] years. Thirteen (59 %) patients had removed a tick from the site of EM (N = 6/13 "suggestive"; N = 7/13 "not very typical"). Sixteen (73 %) lesions were located on the lower limbs, 3 (14 %) on the upper limbs and 3 (14 %) on the trunk (Table 1). Demographic characteristics did not differ between the 2 groups, except for the variable "bite in the garden" (p = 0.03). Clinically, only lesion size was almost significantly greater in the "suggestive" group (median [IQR] = 13 [9–19]cm) than in the "not very typical" group (median [IQR] = 8 [6–15]cm) (p = 0.06) (Table 1), as expected by the classification criteria. The geographical distribution of tick bites of the 22 patients is given in Table 2.

### 2.1. Real-time microfluidic PCR

To search for the presence of 41 known tick-borne pathogens (bacteria and parasites among the genera *Borrelia*, *Anaplasma*, *Rickettsia*, *Ehrlichia*, *Bartonella*, *Coxiella*, *Francisella*, *Babesia*, *Theileria*), the 22 skin biopsies (DNA and RNA), 26 plasma samples (D0 and D14) and 26 blood pellet (D0 and D14) from patients were screened using a high-throughput real-time microfluidic PCRs [14]. Briefly, a 48.48 dynamic array were screened using the BioMark™ real-time PCR system (Standard Biotools, CA, USA). This chip can be used to perform 2304 real-time microfluidic PCR reactions thanks to 48 PCR mixes and 48 samples placed into individual wells, prior to transfer into individual chambers for the reaction. The thermal cycling conditions were 50 °C for 2 min, 95 °C for 10 mins, and 40 cycles at 95 °C for 15 s and 60 °C for 1 min. One negative water control, one inhibitory molecule control (*Escherichia coli* EDL933 strain), and one DNA extraction control (human targeting primers) were added to each chip. Then positive results were confirmed

<table><caption>**Table 1**

Clinical and biological characteristics of the 22 patients.</caption><thead><tr><th></th><th>"Suggestive EM<sup>+</sup>" group<br>N = 12</th><th>"Not typical EM<sup>+</sup>" group<br>N = 10</th><th>Total<br>N = 22</th><th>pvalue</th></tr></thead><tbody><tr><td>Gender, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Male</td><td>9 (75.0)</td><td>4 (40.0)</td><td>13<br>(59.1)</td><td>0.19</td></tr><tr><td>Female</td><td>3 (25.0)</td><td>6 (60.0)</td><td>9 (40.9)</td><td></td></tr><tr><td>Age (years), median (Q1-Q3)</td><td>52 [42–70]</td><td>47<br>[36–60]</td><td>49<br>[40–64]</td><td>0.58</td></tr><tr><td>Presence of a tick bite, n (%)</td><td>10 (83.3)</td><td>8 (80.0)</td><td>18<br>(81.8)</td><td>0.99</td></tr><tr><td>Place of suspected tick bite, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Forest</td><td>9 (75.0)</td><td>5 (50.0)</td><td>14<br>(63.6)</td><td>0.38</td></tr><tr><td>Meadow</td><td>2 (16.7)</td><td>3 (30.0)</td><td>5 (22.7)</td><td>0.62</td></tr><tr><td>Garden</td><td>0 (0)</td><td>4 (40.0)</td><td>4 (18.2)</td><td>0.03</td></tr><tr><td>Other</td><td>2 (16.7)</td><td>1 (10.0)</td><td>3 (13.6)</td><td>0.99</td></tr><tr><td>Appearance of ticks when bitten, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Large and white (swallowed)</td><td>0 (0)</td><td>2 (20.0)</td><td>2 (9.1)</td><td></td></tr><tr><td>Small and black (not swallowed)</td><td>6 (50.0)</td><td>5 (50.0)</td><td>11<br>(50.0)</td><td>0.50</td></tr><tr><td>Missing data</td><td>6 (50.0)</td><td>3 (30.0)</td><td>9 (40.9)</td><td></td></tr><tr><td>Time to onset of ME after tick bite, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Less than a week</td><td>2 (16.7)</td><td>2 (20.0)</td><td>4 (18.2)</td><td></td></tr><tr><td>Between 1 and 2 weeks</td><td>3 (25.0)</td><td>1 (10.0)</td><td>4 (18.2)</td><td></td></tr><tr><td>More than 2 weeks</td><td>6 (50.0)</td><td>6 (60.0)</td><td>12<br>(54.6)</td><td>0.91</td></tr><tr><td>Missing data</td><td>1 (8.3)</td><td>1 (10.0)</td><td>2 (9.1)</td><td></td></tr><tr><td>Anatomical location of ME, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Upper limb</td><td>1 (8.3)</td><td>2 (20.0)</td><td>3 (13.6)</td><td></td></tr><tr><td>Lower limb</td><td>8 (66.7)</td><td>8 (80.0)</td><td>16<br>(72.7)</td><td>0.39</td></tr><tr><td>Trunk</td><td>3 (25.0)</td><td>0 (0)</td><td>3 (13.6)</td><td></td></tr><tr><td>Size of lesion (cms)</td><td>13 [9–19]</td><td>8 [6–15]</td><td>10<br>[7–15]</td><td>0.06</td></tr><tr><td>Clinical signs at inclusion, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Asthenia</td><td>2 (16.7)</td><td>3 (30.0)</td><td>5 (22.7)</td><td>0.62</td></tr><tr><td>Neurological signs</td><td>2 (16.7)</td><td>4 (40.0)</td><td>6 (27.3)</td><td>0.35</td></tr><tr><td>Rheumatological signs</td><td>2 (16.7)</td><td>2 (20.0)</td><td>4 (18.2)</td><td>0.99</td></tr><tr><td>Sensory disorders</td><td>0 (0)</td><td>1 (10.0)</td><td>1 (4.6)</td><td>0.46</td></tr><tr><td>Other clinical signs</td><td>1 (8.3)</td><td>4 (40.0)</td><td>5 (22.7)</td><td>0.14</td></tr><tr><td>Prescription of antibiotic treatment following inclusion visit, n (%)</td><td>12 (100.0)</td><td>10 (100.0)</td><td>22<br>(100.0)</td><td></td></tr><tr><td>Detection of Borreliella, n (%)</td><td>2 (16.7)</td><td>3 (30.0)</td><td>5 (22.7)</td><td>0.62</td></tr></tbody></table>

by Nested PCR.

### 2.2. NGS

To complete the first analysis and look for the presence of common or unexpected pathogens transmitted by tick bite, we carried out a pathogen discovery approach based on random deep sequencing, as described previously in a context of routine analyses of patient's samples in a clinical setting [15]. Briefly, total RNA from skin biopsies taken from the active area of 22 patients presenting with erythema migrans were sequenced using the SMARTer Stranded Total RNA-Seq Kit v2 Pico Input Mammalian. RNA sequencing allows for identification of all types of microorganisms, as all express their genomes as RNA transcripts. Plasma samples from 26 patients (including the 22 patients who had skin biopsy) collected at the time of their first medical visit (D0, concomitant with skin biopsies) and two weeks later (D14) were treated

<footnote>† EM: Erythema Migrans.</footnote>
<table><caption>**Table 2**
Geographical distribution of tick bites of the 22 patients.</caption><thead><tr><th rowspan="2">French department</th><th rowspan="2">French region</th><th>"Suggestive EMj" group</th><th>"Not typical EMj" group</th><th rowspan="2">Total N = 22</th></tr><tr><th>N = 12</th><th>N = 10</th></tr><tr><th></th><th></th><th>n (%)</th><th>n (%)</th><th>n (%)</th></tr></thead><tbody><tr><td>25</td><td>Doubs</td><td>Bourgogne-Franche-Comté</td><td>1 (8.3)</td><td>1 (10.0)</td><td>2 (9.0)</td></tr><tr><td>28</td><td>Eure-et-Loir</td><td>Centre-Val de Loire</td><td>1 (8.3)</td><td>0 (0)</td><td>1 (4.6)</td></tr><tr><td>39</td><td>Jura</td><td>Bourgogne-Franche-Comté</td><td>0 (0)</td><td>2 (20.0)</td><td>2 (9.0)</td></tr><tr><td>42</td><td>Loire</td><td>Auvergne-Rhone-Alpes</td><td>7 (58.4)</td><td>4 (40.0)</td><td>11 (50.0)</td></tr><tr><td>70</td><td>Haute-Saône</td><td>Bourgogne-Franche-Comté</td><td>1 (8.3)</td><td>0 (0)</td><td>1 (4.6)</td></tr><tr><td>72</td><td>Sarthe</td><td>Pays de la Loire</td><td>0 (0)</td><td>1 (10.0)</td><td>1 (4.6)</td></tr><tr><td>74</td><td>Haute-Savoie</td><td>Auvergne-Rhone-Alpes</td><td>1 (8.3)</td><td>0 (0)</td><td>1 (4.6)</td></tr><tr><td>77</td><td>Seine et Marne</td><td>Ile de France</td><td>0 (0)</td><td>1 (10.0)</td><td>1 (4.6)</td></tr><tr><td>Missing data</td><td></td><td></td><td>1 (8.3)</td><td>1 (10.0)</td><td>2 (9.0)</td></tr></tbody></table>

† EM: Erythema Migrans.

by nucleases before total nucleic acids extraction to enrich for nucleic acids associated with bacterial bodies or virus particles, then sequenced using an adapted MALBAC protocol [16]. In parallel, total RNA extracted from blood pellets at D0 and D14 were sequenced using the SMARTer Stranded Total RNA-Seq Kit v2. Sequencing was done on Illumina NextSeq500 and NovaSeq instruments with average depth of 58 million clusters per sample. The search for microbial sequences was done using Kraken2 [17] on NCBI nt (version 2021-3-29) and a complementary search for viruses was achieved through Microseek [18].

### 3. Results

We first screened the 22 skin biopsies (RNA and DNA extracts), plus plasma (total nucleic acids) and blood pellets (RNA) collected at the time of the first medical visit (D0, concomitant with skin biopsies) and two weeks later (D14), for the presence of 41 tick-borne pathogens (bacteria and parasites) using a high-throughput real-time microfluidic PCR. On the 22 skin biopsies RNA extracts, one was positive and confirmed as *Borrelia afzelii*. Regarding the DNA analysis, although 14 samples were positive for *Borrelia* spp. using the microfluidic PCR, only 5 could be confirmed as *Borrelia afzelii* through Sanger sequencing (N = 2/5 "suggestive of EM"; N = 3/5 "not very typical"), including the sample that was positive following the RNA approach. All skin biopsies were negative for the other tick-borne pathogens tested. None of the 41 tick-borne pathogens were detected from plasma and blood pellets.

We then carried out an agnostic pathogen identification approach based on random deep sequencing. On the 22 skin biopsies sequenced (total RNA only), five (22%) were positive for *Borrelia burgdorferi* sensu lato, matching the 5 samples that were positive by the real-time microfluidic PCR from DNA, with 4 out of 5 being attributed to the species *Borrelia afzelii* (the last one attributed to the genus *Borrelia*). Other skin biopsies were negative for *Borrelia* spp., and no other viral, bacterial or fungal sequences were found. Plasma samples at D0 and D14 were all positive for Torque Teno Virus (TTV, *Anelloviridae*), a non-pathogenic virus whose load is proposed as a marker of immune status. TTV viral loads were therefore quantified using the TTV R-GENE real-time PCR assay (BioMérieux SA, Marcy l'Etoile, France), showing an average of 2.4 log₁₀ Cp/mL and a distribution of viral loads not associated with immune dysfunction or increased risk of infection [19,20], with no difference between D0 and D14. A few other viral reads were detected by NGS in certain plasma samples: STL Polyomavirus (N = 14/22, 63%), Merkel Cell Polyomavirus (N = 7/22, 31%) and Epstein-Barr Virus (6/22, 27%) but confirmatory specific PCR were all negative,

suggesting that STLPyV, MCPyV corresponded to skin flora and that EBV viremia was better detected by NGS than by PCR. No microbial sequences were detected from blood pellets.

### 4. Discussion

We challenged the hypothesis of the transmission of unexpected or novel pathogen other than *Borrelia* spp. in tick-bitten patients developing erythema migrans, using a combination of high-throughput real-time microfluidic PCRs and agnostic deep sequencing. No evidence of microorganisms other than *Borrelia* spp. was found in this small but well characterized cohort. The significance of these results is nevertheless limited by the fact that *Borrelia* spp. was identified only in a minority of samples, leaving open the possibility that other etiologies may also have escaped our investigation.

#### Ethics statement

This work was approved on December 18, 2017 by the Comité de Protection des Personnes (CPP) Est-II under protocol number CPP 17/567 (ANSM reference 2017-A02916-47). The CPP in France has a role similar to that of Ethics Committees or Institutional Review Boards (IRBs) in other countries. Volunteers received an information sheet detailing the study, and informed consent was obtained before their participation.

#### Funding

The OHTICKS study was funded by the French Agence Nationale de la Recherche (ANR) grant number ANR-16-CES35-0011.

#### CRediT authorship contribution statement

Philippe Pérot: Writing – review & editing, Writing – original draft, Validation, Methodology, Investigation. Laura Tondeur: Writing – review & editing, Writing – original draft, Resources, Methodology, Data curation. Sara Moutailler: Writing – review & editing, Writing – original draft, Validation, Methodology, Investigation, Conceptualization. Delphine Chrétien: Validation, Methodology, Investigation. Nicole Corre-Catelin: Resources, Project administration, Data curation. Muriel Vayssier-Taussat: Writing – review & editing, Supervision, Project administration, Methodology, Investigation, Funding acquisition, Conceptualization. Marc Eloit: Writing – review & editing, Writing – original draft, Supervision, Methodology, Investigation, Funding acquisition, Conceptualization. Catherine Chirouze: Writing – review & editing, Resources, Methodology, Investigation, Conceptualization. Céline Cazorla: Writing – review & editing, Writing – original draft, Resources, Methodology, Investigation, Conceptualization. Céline Cazorla: Writing – review & editing, Writing – original draft, Resources, Methodology, Investigation, Conceptualization.

#### Declaration of competing interest

The authors declare they have no conflicts of interest.

#### Data availability

Data will be made available on request.

#### Acknowledgments

We thank the Biomics platform of Institut Pasteur for providing access to the sequencers.
References

[1] A. Estrada-Peña, S. Cutler, A. Potkonjak, M. Vassier-Tussaut, W. Van Bortel, H. Zeller, N. Fernández-Ruiz, A.D. Mihalca, An updated meta-analysis of the distribution and prevalence of *Borrelia burgdorferi* s.l. in ticks in Europe, Int. J. Health Geogr. 17 (2018) 41, https://doi.org/10.1186/s12942-018-0163-7.

[2] M. Strnad, V. Hönig, D. Rüzek, L. Grubhoffer, R.O.M. Rego, Europe-wide Meta-analysis of *Borrelia burgdorferi* Sensu Lato prevalence in questing *Ixodes ricinus* ticks, Appl. Environ. Microbiol. 83 (2017) e00609-17, https://doi.org/10.1128/AEM.00609-17.

[3] Borréliose de Lyme et Autres Maladies Vectorielles à Tiques, Haute Autorité de Santé. https://www.has-sante.fr/jcms/c_2857558/fr/borreliose-de-lyme-et-autres-maladies-vectorielles-a-tiques, 2024 accessed June 24, 2024.

[4] N. Mencke, Future challenges for parasitology: vector control and “one health” in Europe: the veterinary medicinal view on CVBDs such as tick borreliosis, rickettsiosis and canine leishmaniosis, Vet. Parasitol. 195 (2013) 256–271, https://doi.org/10.1016/j.vetpar.2013.04.007.

[5] A. Mannelli, L. Bertolotti, L. Gern, J. Gray, Ecology of *Borrelia burgdorferi* sensu lato in Europe: transmission dynamics in multi-host systems, influence of molecular processes and effects of climate change, FEMS Microbiol. Rev. 36 (2012) 837–861, https://doi.org/10.1111/j.1574-6976.2011.00312.x.

[6] K.M. Hansford, B.W. Wheeler, B. Tschirren, J.M. Medlock, Questing *Ixodes ricinus* ticks and Borrelia spp. in urban green space across Europe: a review, Zoonoses Public Health 69 (2022) 153–166, https://doi.org/10.1111/zph.12913.

[7] I. Rochlin, A. Toledo, Emerging tick-borne pathogens of public health importance: a mini-review, J. Med. Microbiol. 69 (2020) 781–791, https://doi.org/10.1099/jmm.0.001206.

[8] A. Zintl, A. McManus, M. Galan, M. Diquattro, L. Giuffredi, N. Charbonnel, J. Gray, C. Holland, P. Stuart, Presence and identity of Babesia microti in Ireland, Ticks Tick Borne Dis. 14 (2023) 102221, https://doi.org/10.1016/j.ttbdis.2023.102221.

[9] J.C. Semenza, J.E. Suk, Vector-borne diseases and climate change: a European perspective, FEMS Microbiol. Lett. 365 (2018) fnx244, https://doi.org/10.1093/femsle/fnx244.

[10] S. Melis, G. Batisti Biffignandi, E. Olivieri, C. Galon, N. Vicari, P. Prati, S. Moutailler, D. Sassera, M. Castelli, High-throughput screening of pathogens in *Ixodes ricinus* removed from hosts in Lombardy, northern Italy, Ticks Tick Borne Dis. 15 (2024) 102285, https://doi.org/10.1016/j.ttbdis.2023.102285.

[11] R.A. Sykes, P. Makiello, An estimate of Lyme borreliosis incidence in Western Europe†, J. Public Health 39 (2017) 74–81, https://doi.org/10.1093/pubmed/fdw017.

[12] CDC, Lyme Disease Surveillance and Data, Lyme Disease. https://www.cdc.gov/lyme/data-research/facts-stats/index.html, 2024 accessed June 24, 2024.

[13] L. Geebelen, B. Devleesschauwer, T. Lernout, K. Tersago, Y. Parmentier, H. Van Oyen, N. Speybroeck, P. Beutels, Lyme borreliosis in Belgium: a cost-of-illness analysis, BMC Public Health 22 (2022) 2194, https://doi.org/10.1186/s12889-022-14380-6.

[14] L. Michelet, S. Delannoy, E. Devillers, G. Umhang, A. Aspan, M. Juremalm, J. Chirico, F.J. van der Wal, H. Sprong, T.P. Boye Pihl, K. Klitgaard, R. Bødker, P. Fach, S. Moutailler, High-throughput screening of tick-borne pathogens in Europe, Front. Cell. Infect. Microbiol. 4 (2014) 103, https://doi.org/10.3389/fcimb.2014.00103.

[15] J. Fourgeaud, B. Regnault, V. Ok, N. Da Rocha, É. Sitterlé, M. Mekouar, H. Faury, C. Milliancourt-Seels, F. Jagorel, D. Chrétien, T. Bigot, É. Troadec, I. Marques, A. Serris, D. Seilhean, B. Neven, P. Frange, A. Ferroni, M. Lecuit, X. Nassif, O. Lortholary, M. Leruez-Ville, P. Pérot, M. Eloit, A. Jamet, Performance of clinical metagenomics in France: a prospective observational study, Lancet Microbe (2023), https://doi.org/10.1016/S2666-5247(23)00244-6. S2666-5247(23) 00244-6.

[16] B. Regnault, T. Bigot, L. Ma, P. Pérot, S. Temmam, M. Eloit, Deep impact of random amplification and library construction methods on viral metagenomics results, Viruses 13 (2021) 253, https://doi.org/10.3390/v13020253.

[17] D.E. Wood, J. Lu, B. Langmead, Improved metagenomic analysis with kraken 2, Genome Biol. 20 (2019) 257, https://doi.org/10.1186/s13059-019-1891-0.

[18] P. Pérot, T. Bigot, S. Temmam, B. Regnault, M. Eloit, Microseek: a protein-based metagenomic pipeline for virus diagnostic and discovery, Viruses 14 (2022) 1990, https://doi.org/10.3390/v14091990.

[19] K. Doberer, F. Haupenthal, M. Nackenhorst, F. Bauernfeind, F. Dermuth, M. Eigenschink, M. Schiemann, J. Kläger, I. Görzer, F. Eskandary, R. Reindl-Schwaighofer, Ž. Kikić, G. Böhmig, R. Strassl, H. Regele, E. Puchhammer-Stöckl, G. Bond, Torque Teno virus load is associated with subclinical alloreactivity in kidney transplant recipients: a prospective observational trial, Transplantation 105 (2021) 2112–2118, https://doi.org/10.1097/TP.00000000000003619.

[20] M. Schiemann, E. Puchhammer-Stöckl, F. Eskandary, P. Kohlbeck, S. Rasoul-Rockenschaub, A. Heilos, N. Kozakowski, I. Görzer, Ž. Kikić, H. Herkner, G. A. Böhmig, G. Bond, Torque Teno virus load-inverse association with antibody-mediated rejection after kidney transplantation, Transplantation 101 (2017) 360–367, https://doi.org/10.1097/TP.00000000000001455.