- Research Note
- Open Access
Isolation and multilocus sequence typing of Borrelia burgdorferi from Ixodes scapularis collected from dogs in Ontario, Canada
BMC Research Notes volume 16, Article number: 43 (2023)
To identify the multilocus sequence typing (MLST) sequence types of Borrelia burgdorferi from Ixodes scapularis in Ontario, Canada.
One hundred and eighty-five I. scapularis ticks were submitted from 134 dogs via participating clinics from April 1, 2019, to March 31, 2020. Seventeen MLST sequence types of B. burgdorferi were detected from fifty-eight cultured isolates from 21 ticks. The most common MLST sequence types were 12 and 16. Mixed infections of two MLST sequence types were detected in four ticks. Three sequence types (48, 317, 639) were new detections in Ontario.
Over the last two decades, there has been rapid northward range expansion of the blacklegged tick, Ixodes scapularis, in Ontario, Canada [1, 2]. This tick is the vector for several pathogens, including Borrelia burgdorferi sensu stricto (herein referred to as B. burgdorferi), the causative agent of Lyme disease in humans, dogs, and horses [3,4,5]. Lyme disease is the most common human vector-borne disease in Canada, with a notable increase in case counts from 144 reported in 2009 to 3147 reported in 2021 . Evidence of exposure to B. burgdorferi in dogs and horses has also been detected in areas of central and eastern Canada [7, 8].
Borrelia burgdorferi sensu stricto can be further characterized by genetic analyses. Restriction fragment length polymorphism (RFLP) analysis of the 16-23 S intragenic spacer (IGS) region and sequencing of the plasmid-encoded OspA and OspC surface proteins have been employed [9,10,11]. Multilocus sequence typing (MLST) uses the combined sequences of eight chromosomal bacterial housekeeping genes. Given these housekeeping genes are slow to evolve, this approach has been shown to have higher discriminatory power than previously utilized methods to characterize genetic diversity and is currently the recommended approach [9,10,11].
Associations have been detected between sequence types of B. burgdorferi and geographic distribution, reservoir host species, immune response in humans and severity of human clinical disease [12,13,14,15,16,17,18,19]. Approximately 5% of dogs exposed to B. burgdorferi will develop clinical disease, which is characterized by fever, anorexia and shifting lameness, and in rare cases, the disease can progress to a potentially fatal protein-losing nephropathy . Similarly in horses, only a subset of those exposed develop clinical disease, with three documented syndromes: cutaneous pseudolymphoma, uveitis, and a highly fatal neuroborreliosis . It is unknown if the variability of clinical disease presentation and severity in dogs and horses has any association with B. burgdorferi sequence type.
As emergence of I. scapularis and B. burgdorferi continue across Ontario, it is pertinent to expand our knowledge base on sequence types and the associated ecological, epidemiological, and clinical patterns. To date, a handful of studies have characterized MLST sequence types of B. burgdorferi in Canada [12,13,14,15,16,17]. This study builds on these efforts by describing the MLST sequence types of B. burgdorferi in ticks collected from dogs in southern Ontario.
A convenience sample of veterinary clinics from across Ontario was invited to participate in this study via email. Live I. scapularis ticks were collected from dogs (e.g., the tick could be attached or crawling on the animal) at participating clinics and shipped immediately to our laboratory during a one-year period (April 2019 – March 2020). Each tick submission was accompanied by a questionnaire that documented patient information, including use of antibiotics during the last two months, suspected geographic location of tick acquisition, and date of tick removal.
Tick samples were morphologically identified using a stereoscope and standard keys to confirm the species identification, life stage, and sex (for adult samples) . Ticks were surface disinfected by submerging and briefly vortexing each sample sequentially in the following reagents: sterile water, 0.5 − 1% sodium hypochlorite, 0.5% benzalkonium chloride, 3% hydrogen peroxide, 70% ethanol, sterile water. Each tick was macerated using a sterile scalpel blade, and a tick homogenate was prepared by adding 969 µl 1.5x complete BSK and 31 µl of an antibiotic mixture (phosphomycin (2 mg/mL), rifampicin (5 mg/mL), amphotericin B (0.25 mg/mL), kanamycin (10 µg/mL)). At this point, a 100 µL aliquot was removed from the sample for DNA extraction followed by real-time PCR targeting the 23S rRNA gene of B. burgdorferi sensu lato, as previously described . DNA was extracted using the DNeasy Blood and Tissue Kit [Qiagen, Toronto, ON] per manufacturer’s protocol, except ticks were incubated with the digestive reagent, proteinase K, overnight (vs. the prescribed 2 h). Real-time PCR runs included both negative (water) and positive controls. The tick homogenates of samples that were PCR positive for B. burgdorferi were plated using a reagent combination of 20 mL 1.5x complete BSK, 10 mL 2.1% agarose, and 930 µL antibiotic mixture. Plates were placed in a 37o Celsius incubator with 5% CO2 and monitored for at least 31 days. If colonies appeared, three colonies were sampled from each plate to prepare isolates. In cases where there were less than three colonies present, all colonies were taken. Each colony was extracted by removing the agarose plug and placing it in BSK-H with 10% glycerol for expansion. A 150 µL aliquot was removed for DNA extraction (as described above), followed by multilocus sequence typing (MLST). MLST was conducted for eight housekeeping genes (clpA, clpX, nifS, pepX, pyrG, recG, rplB and uvrA). Bidirectional automated Sanger sequencing was performed using previously published primers [9, 10]. Forward and reverse reads were assembled and trimmed using Geneious Prime (Version 2020.2.4, Biomatters Ltd., Auckland, New Zealand, https://www.geneious.com). Chromatograms were visualized to validate base calls. Sequences were then uploaded to PubMLST to determine locus identities and sequence types .
Spatial data were prepared using QGIS version 3.22.13 (www.qgis.org; 2022). Vector layers were accessed through the Scholars GeoPortal at the University of Guelph (http://geo2.scholarsportal.info).
One hundred and eighty-five I. scapularis ticks were submitted from 134 dogs from 19 participating clinics from April 1, 2019, to March 31, 2020 (Fig. 1). Nineteen of these ticks (10.3%) were submitted from 12 dogs that either had reported antibiotic use (5 dogs, 3.7%) or an unknown history (i.e., no information was provided, or owner did not know) (7 dogs, 5.2%). Thirty-eight (20.5%) ticks from 29 (21.6%) dogs were positive for B. burgdorferi by PCR. Four of these positive ticks were from dogs with reported or unknown antibiotic use. Borrelia burgdorferi was isolated from 21/38 (55.2%) PCR positive ticks. No culture growth was detected from the four positive ticks removed from dogs with reported or unknown antibiotic use. Fifty-eight B. burgdorferi isolates were subjected to MLST, but read quality was poor for three isolates. In total, 55 isolates were successfully typed. Seventeen B. burgdorferi MLST sequence types were detected, all of which have been previously described (Table 1). Four ticks had mixed infections, with isolates of different MLST sequence types identified concurrently (Table 1).
A diversity of B. burgdorferi MLST sequence types were isolated from I. scapularis collected in Ontario, Canada. Fourteen of the MLST sequence types isolated in this study have widespread distribution across northeastern and/or midwestern North America and have been previously identified from ticks and/or hosts in Ontario [12,13,14,15, 17]. However, three B. burgdorferi MLST sequence types identified in our study had only been reported in the provinces of Manitoba (48, 639) and Prince Edward Island (317), and thus represent new records for Ontario [14, 22]. Tick host movement, particularly that of migratory birds, has driven long distance dispersal of I. scapularis. However, this is predominately along a north-south corridor and these provinces are to the west (Manitoba) and east (Prince Edward Island) of Ontario . Movement of reservoir hosts, such as small mammals, via connected habitats, has been associated with geographic patterns of MLST sequence types . In this context, movement occurs on a local scale, not vast distances between these provinces. Additional research is needed to elucidate the ecological processes that drive MSLT sequence type diversity and dispersal. Moreover, as we develop a comprehensive understanding of geographic distribution of sequence types, further studies are warranted to explore the relevance of these sequence types to the epidemiology and clinical manifestation of Lyme disease, not only in humans, but dogs and horses as well.
The main limitations for this study were:
Specimen quality: To detect and subsequently culture B. burgdorferi from I. scapularis, high quality live specimens are required. Although collecting ticks via veterinary clinics was an efficient way to acquire samples, many were highly engorged, and this can impact DNA extraction and thus PCR analyses. Although all samples were requested to be shipped immediately, some ticks were dead upon arrival, which may have contributed to unsuccessful attempts to culture B. burgdorferi.
Geographic data: Location of tick acquisition was acquired through an owner-reported survey. Ticks can be attached for several days, over which time a dog can visit multiple locations, so it can be difficult to determine the precise location at which the tick was acquired.
Sample size: Given only a small number of isolates were available for MLST, there was insufficient power to explore any advanced spatial analyses (e.g., cluster analysis).
Genetic analysis: Although MLST has high discriminatory power, additional insight could have been gained by conducting whole genome analysis of isolates.
Borrelia spp. isolates from this study have been uploaded to PubMLST with identification numbers 3456 to 3510 (https://pubmlst.org).
Ogden NH, St-Onge L, Barker IK, Brazeau S, Bigras-Poulin M, Charron DF, et al. Risk maps for range expansion of the Lyme disease vector, Ixodes scapularis, in Canada now and with climate change. Int J Health Geogr. 2008;7:24.
Robinson EL, Jardine CM, Koffi JK, Russell C, Lindsay LR, Dibernardo A, et al. Range expansion of Ixodes scapularis and Borrelia burgdorferi in Ontario, Canada, from 2017 to 2019. Vector-Borne Zoonotic Dis. 2022;22(7):361–9.
Burgdorfer W, Barbour AG, Hayes SF, Benach JL, Davis JP. Lyme Disease - A tick-borne spirochetosis? Science. 1982;216(4552):1317–9.
Littman MP, Gerber B, Goldstein RE, Labato MA, Lappin MR, Moore GE. ACVIM consensus update on Lyme borreliosis in dogs and cats. J Vet Intern Med. 2018;32(3):887–903.
Divers TJ, Gardner RB, Madigan JE, Witonsky SG, Bertone JJ, Swinebroad EL, et al. Borrelia burgdorferi infection and Lyme disease in north american horses: a consensus statement. J Vet Intern Med. 2018;32(2):617–32.
Lyme disease: Surveillance. Government of Canada. 2023. https://www.canada.ca/en/public-health/services/diseases/lyme-disease/surveillance-lyme-disease.html. Accessed 21 Feb 2023.
Evason M, Stull JW, Pearl DL, Peregrine AS, Jardine C, Buch JS et al. Prevalence of Borrelia burgdorferi, Anaplasma spp., Ehrlichia spp. and Dirofilaria immitis in Canadian dogs, 2008 to 2015: a repeat cross-sectional study. Parasites Vectors. 2019;12:64.
Neely M, Arroyo LG, Jardine C, Moore A, Hazlett M, Clow KM, et al. Seroprevalence and evaluation of risk factors associated with seropositivity for Borrelia burgdorferi in Ontario horses. Equine Vet J. 2020;53:331–8.
Margos G, Gatewood AG, Aanensen DM, Hanincová K, Terekhova D, Vollmer SA, et al. MLST of housekeeping genes captures geographic population structure and suggests a european origin of Borrelia burgdorferi. Proc Natl Acad Sci. 2008;105(25):8730–5.
Hoen AG, Margos G, Bent SJ, Diuk-Wasser MA, Barbour A, Kurtenbach K, Fish D. Phylogeography of Borrelia burgdorferi in the eastern United States reflects multiple independent Lyme disease emergence events. Proc Natl Acad Sci USA. 2009;106(35):15013–8.
Ogden NH, Feil EJ, Leighton PA, Lindsay LR, Margos G, Mechai S, et al. Evolutionary aspects of emerging Lyme disease in Canada. Appl Environ Microbiol. 2015;81(21):7350–9.
Ogden NH, Bouchard C, Kurtenbach K, Margos G, Lindsay LR, Trudel L, et al. Active and passive surveillance and phylogenetic analysis of Borrelia burgdorferi elucidate the process of Lyme disease risk emergence in Canada. Environ Health Perspect. 2010;118(7):909–14.
Ogden NH, Margos G, Aanensen DM, Drebot MA, Feil EJ, Hanincova K, et al. Investigation of genotypes of Borrelia burgdorferi in Ixodes scapularis ticks collected during surveillance in Canada. Appl Environ Microbiol. 2011;77(10):3244–54.
Mechai S, Margos G, Feil EJ, Lindsay LR, Ogden NH. Complex population structure of Borrelia burgdorferi in southeastern and south-central Canada as revealed by phylogeographic analysis. Appl Environ Microbiol. 2015;81(4):1309–18.
Mechai S, Margos G, Feil EJ, Barairo N, Lindsay LR, Michel P et al. Evidence for host-genotype associations of Borrelia burgdorferi sensu stricto.PLoS One. 2016;1–26.
Ogden NH, Arsenault J, Hatchette TF, Mechai S, Lindsay LR. Antibody responses to Borrelia burgdorferi detected by western blot vary geographically in Canada. PLoS ONE. 2017;12(2):1–14.
Mechai S, Margos G, Feil EJ, Lindsay LR, Michel P, Kotchi SO, Ogden NH. Evidence for an effect of landscape connectivity on Borrelia burgdorferi sensu stricto dispersion in a zone of range expansion. Ticks Tickborne Dis. 2018;9:1407–15.
Hanincova K, Mukherjee P, Ogden N, Margos G, Wormser GP, Reed KD, et al. Multilocus sequence typing of Borrelia burgdorferi suggests existence of lineages of differential pathogenic properties in humans. PLoS ONE. 2013;8(9):e73066.
Khatchikian CE, Nadelman RB, Nowakowski J, Schwartz I, Wormser GP, Brisson D. Evidence for strain-specific immunity in patients treated for early Lyme disease. Infect Immun. 2014;82(4):1408–13.
Lindquist EE, Galloway TD, Artsob H, Lindsay LR, Drebot M, Wood H, et al. A handbook to the ticks of Canada (Ixodida: Ixodidae, Argasidae). Biological Survey of Canada; 2016.
Courtney JW, Kostelnik LM, Zeidner NS, Massung RF. Multiplex real-time PCR for detection of Anaplasma phagocytophilum and Borrelia burgdorferi. J Clin Microbiol. 2004;42(7):3164–8.
Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018;24(3):124. https://doi.org/10.12688/wellcomeopenres.
Ogden NH, Lindsay LR, Hanincová K, Barker IK, Bigras-Poulin M, Charron DF, et al. Role of migratory birds in introduction and range expansion of Ixodes scapularis ticks and of Borrelia burgdorferi and Anaplasma phagocytophilum in Canada. Appl Environ Microbiol. 2008;74(6):1780–90.
The authors would like to sincerely thank:
- All participating veterinary clinics, clients, and patients.
- Dr. Robbin Lindsay and Antonia Dibernardo (National Microbiology Laboratory, Public Health Agency Canada) for sharing their protocol for culturing B. burgdorferi from ticks.
- Joyce Rousseau (Centre for Public Health and Zoonoses, University of Guelph) for assistance with laboratory techniques.
- Drs. Nicole Ricker (Department of Pathobiology, University of Guelph), Nicholas Ogden (National Microbiology Laboratory, Public Health Agency of Canada) and Gabriele Margos (Bavarian Health and Food Safety Authority) for sharing their expertise on data analyses.
This study was funded by the Ontario Veterinary College Pet Trust Fund. GKN was supported by an Undergraduate Research Assistantship from the University of Guelph.
This study was approved by the Animal Care Committee at the University of Guelph (Animal Use Protocol #3981). All methods were performed in accordance with the approved protocol.
Consent for publication
The authors have no competing interests to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Nichol, G.K., Weese, J.S. & Clow, K.M. Isolation and multilocus sequence typing of Borrelia burgdorferi from Ixodes scapularis collected from dogs in Ontario, Canada. BMC Res Notes 16, 43 (2023). https://doi.org/10.1186/s13104-023-06315-0
- Borrelia burgdorferi
- Ixodes scapularis
- Multilocus sequence typing