
==== Front
One Health
One Health
One Health
2352-7714
Elsevier

S2352-7714(24)00212-X
10.1016/j.onehlt.2024.100886
100886
Short Communication
Broad range molecular detection methods identify only Borrelia spp. in erythema migrans biopsies and blood of tick-bitten patients
Pérot Philippe philippe.perot@pasteur.fr
ab⁎
Tondeur Laura c
Moutailler Sara d
Chrétien Delphine ab
Corre-Catelin Nicole e
Vayssier-Taussat Muriel f
Eloit Marc abg
Chirouze Catherine hi
Cazorla Céline j
OHTICKS ConsortiumArowas Laurence
Botelho-Nevers Elisabeth
Cazorla Céline
Chirouze Catherine
Chrétien Delphine
Corre-Catelin Nicole
Eloit Marc
Frey-Klett Pascale
Fontanet Arnaud
Galon Clémence
Lacombe Karine
Moutailler Sara
Pérot Philippe
Perronne Véronique
Piquard Valentine
Préau Marie
Puppo Costanza
Tondeur Laura
Ungeheuer Marie-Noelle
Vayssier-Taussat Muriel
Zayoud Ayla

a Pathogen Discovery Laboratory, Institut Pasteur, Université Paris Cité, F-75015 Paris, France
b 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
c Emerging Infectious Diseases Unit, Institut Pasteur, Paris, France
d ANSES, INRAE, Ecole Nationale Vétérinaire d'Alfort, UMR BIPAR, Laboratoire de Santé Animale, Maisons-Alfort, F-94700, France
e Institut Pasteur, Université Paris Cité, Clinical Investigation and Access to Research Bioresources (ICAReB) Platform, Paris, France
f Animal Health Division, INRAE, 37380 Nouzilly, France
g Ecole Nationale Vétérinaire d'Alfort, F-94700 Maisons-Alfort, France
h Chrono-Environnement UMR6249, CNRS, Université Bourgogne Franche-Comté, F-25000 Besançon, France
i Infectious Disease Department, University Hospital of Besançon, F-25000, France
j Infectious Disease department, University Hospital of Saint Etienne, Avenue Albert-Raimond, 42055, Saint Etienne Cedex 02, France
⁎ Corresponding author at: Institut Pasteur, Pathogen Discovery Laboratory, 25-28 rue du Dr. Roux, 75015 Paris, France. philippe.perot@pasteur.fr
30 8 2024
12 2024
30 8 2024
19 10088626 6 2024
27 8 2024
27 8 2024
© 2024 The Authors
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
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 microorganisms other than Borrelia spp. was found.

Keywords

Ticks
Lyme disease
Erythema migrans
Metagenomics
RT-PCR
Pathogen discovery
==== Body
pmc1 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], [7], [8], [9], [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 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.Table 1 Clinical and biological characteristics of the 22 patients.

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

Table 2 Geographical distribution of tick bites of the 22 patients.

Table 2French department	French region	“Suggestive EM†” group
N = 12

n (%)	“Not typical EM†” group
N = 10

n (%)	Total
N = 22

n (%)	
25	Doubs	Bourgogne-Franche-Comté	1 (8.3)	1 (10.0)	2 (9.0)	
28	Eure-et-Loir	Centre-Val de Loire	1 (8.3)	0 (0)	1 (4.6)	
39	Jura	Bourgogne-Franche-Comté	0 (0)	2 (20.0)	2 (9.0)	
42	Loire	Auvergne- Rhone-Alpes	7 (58.4)	4 (40.0)	11 (50.0)	
70	Haute-Saône	Bourgogne-Franche-Comté	1 (8.3)	0 (0)	1 (4.6)	
72	Sarthe	Pays de la Loire	0 (0)	1 (10.0)	1 (4.6)	
74	Haute-Savoie	Auvergne- Rhone-Alpes	1 (8.3)	0 (0)	1 (4.6)	
77	Seine et Marne	Ile de France	0 (0)	1 (10.0)	1 (4.6)	
Missing data	1 (8.3)	1 (10.0)	2 (9.0)	
† EM: Erythema Migrans.

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 mins, 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 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 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 log10 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.
==== Refs
References

1 Estrada-Peña A. Cutler S. Potkonjak A. Vassier-Tussaut M. Van Bortel W. Zeller H. Fernández-Ruiz N. Mihalca A.D. 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 10.1186/s12942-018-0163-7 30514310
2 Strnad M. Hönig V. Růžek D. Grubhoffer L. Rego R.O.M. Europe-wide Meta-analysis of Borrelia burgdorferi Sensu Lato prevalence in questing Ixodes ricinus ticks Appl. Environ. Microbiol. 83 2017 e00609-17 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 Mencke N. 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 10.1016/j.vetpar.2013.04.007 23680539
5 Mannelli A. Bertolotti L. Gern L. Gray J. 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 10.1111/j.1574-6976.2011.00312.x 22091928
6 Hansford K.M. Wheeler B.W. Tschirren B. Medlock J.M. Questing Ixodes ricinus ticks and Borrelia spp. in urban green space across Europe: a review Zoonoses Public Health 69 2022 153 166 10.1111/zph.12913 35122422
7 Rochlin I. Toledo A. Emerging tick-borne pathogens of public health importance: a mini-review J. Med. Microbiol. 69 2020 781 791 10.1099/jmm.0.001206 32478654
8 Zintl A. McManus A. Galan M. Diquattro M. Giuffredi L. Charbonnel N. Gray J. Holland C. Stuart P. Presence and identity of Babesia microti in Ireland Ticks Tick Borne Dis. 14 2023 102221 10.1016/j.ttbdis.2023.102221
9 Semenza J.C. Suk J.E. Vector-borne diseases and climate change: a European perspective FEMS Microbiol. Lett. 365 2018 fnx244 10.1093/femsle/fnx244
10 Melis S. Batisti Biffignandi G. Olivieri E. Galon C. Vicari N. Prati P. Moutailler S. Sassera D. Castelli M. High-throughput screening of pathogens in Ixodes ricinus removed from hosts in Lombardy, northern Italy Ticks Tick Borne Dis. 15 2024 102285 10.1016/j.ttbdis.2023.102285
11 Sykes R.A. Makiello P. An estimate of Lyme borreliosis incidence in Western Europe† J. Public Health 39 2017 74 81 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 Geebelen L. Devleesschauwer B. Lernout T. Tersago K. Parmentier Y. Van Oyen H. Speybroeck N. Beutels P. Lyme borreliosis in Belgium: a cost-of-illness analysis BMC Public Health 22 2022 2194 10.1186/s12889-022-14380-6 36443755
14 Michelet L. Delannoy S. Devillers E. Umhang G. Aspan A. Juremalm M. Chirico J. van der Wal F.J. Sprong H. Boye Pihl T.P. Klitgaard K. Bødker R. Fach P. Moutailler S. High-throughput screening of tick-borne pathogens in Europe Front. Cell. Infect. Microbiol. 4 2014 103 10.3389/fcimb.2014.00103 25120960
15 Fourgeaud J. Regnault B. Ok V. Da Rocha N. Sitterlé É. Mekouar M. Faury H. Milliancourt-Seels C. Jagorel F. Chrétien D. Bigot T. Troadec É. Marques I. Serris A. Seilhean D. Neven B. Frange P. Ferroni A. Lecuit M. Nassif X. Lortholary O. Leruez-Ville M. Pérot P. Eloit M. Jamet A. Performance of clinical metagenomics in France: a prospective observational study Lancet Microbe 2023 10.1016/S2666-5247(23)00244-6 S2666-5247(23)00244–6
16 Regnault B. Bigot T. Ma L. Pérot P. Temmam S. Eloit M. Deep impact of random amplification and library construction methods on viral metagenomics results Viruses 13 2021 253 10.3390/v13020253 33562285
17 Wood D.E. Lu J. Langmead B. Improved metagenomic analysis with kraken 2 Genome Biol. 20 2019 257 10.1186/s13059-019-1891-0 31779668
18 Pérot P. Bigot T. Temmam S. Regnault B. Eloit M. Microseek: a protein-based metagenomic pipeline for virus diagnostic and discovery Viruses 14 2022 1990 10.3390/v14091990 36146797
19 Doberer K. Haupenthal F. Nackenhorst M. Bauernfeind F. Dermuth F. Eigenschink M. Schiemann M. Kläger J. Görzer I. Eskandary F. Reindl-Schwaighofer R. Kikić Ž. Böhmig G. Strassl R. Regele H. Puchhammer-Stöckl E. Bond G. Torque Teno virus load is associated with subclinical alloreactivity in kidney transplant recipients: a prospective observational trial Transplantation 105 2021 2112 2118 10.1097/TP.0000000000003619 33587432
20 Schiemann M. Puchhammer-Stöckl E. Eskandary F. Kohlbeck P. Rasoul-Rockenschaub S. Heilos A. Kozakowski N. Görzer I. Kikić Ž. Herkner H. Böhmig G.A. Bond G. Torque Teno virus load-inverse association with antibody-mediated rejection after kidney transplantation Transplantation 101 2017 360 367 10.1097/TP.0000000000001455 27525643
