Skip to content
1887

Abstract

The main agent of Lyme borreliosis (LB) in North America, the bacteria , is spreading in Canada following the northward expansion of its primary tick vectors, and . Despite the importance of this pathogen for human health, the precise geographical origin and genome structure of Canadian strains remain to be determined, and no complete genome sequence from this region is available. The complex genome structure of species makes their assembly challenging, but the latest long-read sequencing technologies and bioinformatics software now enable assembly of genomes with high efficacy. In this study, we sequenced and assembled the genomes of six Canadian strains to compare their content and structure to additional genomes from the USA and Europe. We successfully reconstructed the genome, comprising chromosomes and plasmids of the six Canadian strains. These genomes showed an overall similar structure compared to other strains. Phylogenetic inferences highlighted topological differences in the placement of strains between the chromosome and the cp26 and lp54 plasmids. Synteny analyses revealed important replicon sequence conservation across strains while highlighting a high proportion of shared gene sequences among the replicons of the same strain, especially for cp32 plasmids. We describe the first complete genomes of Canadian strains and present a strategy for the assembly, annotation, comparative analysis of plasmids and their evolution in the same bacterial genus. While the genome content and structure of Canadian strains are similar to other strains, the information in the plasmids and genes they harbour will be useful to elucidate the origins and evolution of LB in Canada.

Funding
This study was supported by the:
  • Caron Thorburn Institute
    • Principal Award Recipient: AnthonyPiot
  • G. Magnotta Foundation for Vector-Borne Diseases
    • Principal Award Recipient: AnthonyPiot
  • Génome Québec
    • Principal Award Recipient: AnthonyPiot
  • This is an open-access article distributed under the terms of the Creative Commons Attribution License.
Loading

Article metrics loading...

/content/journal/mgen/10.1099/mgen.0.001606
2026-01-08
2026-01-14

Metrics

Loading full text...

Full text loading...

/deliver/fulltext/mgen/12/1/mgen001606.html?itemId=/content/journal/mgen/10.1099/mgen.0.001606&mimeType=html&fmt=ahah

References

  1. Ogden NH, Trudel L, Artsob H, Barker IK, Beauchamp G et al. Ixodes scapularis ticks collected by passive surveillance in Canada: analysis of geographic distribution and infection with Lyme borreliosis agent Borrelia burgdorferi. J Med Entomol 2006; 43:600–609 [View Article] [PubMed]
    [Google Scholar]
  2. Ogden NH, St-Onge L, Barker IK, Brazeau S, Bigras-Poulin M 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 [View Article] [PubMed]
    [Google Scholar]
  3. Leighton PA, Koffi JK, Pelcat Y, Lindsay LR, Ogden NH. Predicting the speed of tick invasion: an empirical model of range expansion for the Lyme disease vector Ixodes scapularis in Canada. J Appl Ecol 2012; 49:457–464 [View Article]
    [Google Scholar]
  4. Li F, Mu Q, Ma D, Wu Q. Predicting the potential global distribution of Ixodes pacificus under climate change. PLoS One 2024; 19:e0309367 [View Article]
    [Google Scholar]
  5. Wilson C, Gasmi S, Bourgeois A-C, Badcock J, Carr J et al. Surveillance for Ixodes scapularis and Ixodes pacificus ticks and their associated pathogens in Canada, 2020. CCDR 2023; 49:288–298 [View Article]
    [Google Scholar]
  6. Public Health Agency of Canada Lyme disease: Monitoring; 2015 https://www.canada.ca/en/public-health/services/diseases/lyme-disease/surveillance-lyme-disease.html accessed 30 December 2024
  7. Tyler S, Tyson S, Dibernardo A, Drebot M, Feil EJ et al. Whole genome sequencing and phylogenetic analysis of strains of the agent of Lyme disease Borrelia burgdorferi from Canadian emergence zones. Sci Rep 2018; 8:10552 [View Article] [PubMed]
    [Google Scholar]
  8. Russell JN, Lee M-K, Uyaguari-Diaz MI, Sies AN, Suchan DM et al. Whole-genome sequencing of Western Canadian Borrelia spp. collected from diverse tick and animal hosts reveals short-lived local genotypes interspersed with longer-lived continental genotypes. Microbial Genomics 2024; 10:001276 [View Article]
    [Google Scholar]
  9. Margos G, Gofton A, Wibberg D, Dangel A, Marosevic D et al. The genus Borrelia reloaded. PLoS One 2018; 13:e0208432 [View Article] [PubMed]
    [Google Scholar]
  10. Margos G, Fingerle V, Cutler S, Gofton A, Stevenson B et al. Controversies in bacterial taxonomy: The example of the genus Borrelia. Ticks Tick Borne Dis 2020; 11:101335 [View Article] [PubMed]
    [Google Scholar]
  11. Wolcott KA, Margos G, Fingerle V, Becker NS. Host association of Borrelia burgdorferi sensu lato: a review. Ticks Tick Borne Dis 2021; 12:101766 [View Article] [PubMed]
    [Google Scholar]
  12. Barbour AG. Relapsing fever. In Tick-Borne Diseases of Humans John Wiley & Sons, Ltd; 2005 pp 268–291 [View Article]
    [Google Scholar]
  13. Barbour AG, Adeolu M, Gupta RS. Division of the genus Borrelia into two genera (corresponding to Lyme disease and relapsing fever groups) reflects their genetic and phenotypic distinctiveness and will lead to a better understanding of these two groups of microbes (Margos et al. (2016) There is inadequate evidence to support the division of the genus Borrelia. Int J Syst Evol Microbiol 2016; 67:2058–2067 [View Article]
    [Google Scholar]
  14. Gupta RS. Distinction between Borrelia and Borreliella is more robustly supported by molecular and phenotypic characteristics than all other neighbouring prokaryotic genera: response to Margos’ et al. “The genus Borrelia reloaded” (PLoS ONE 13(12): e0208432). PLoS One 2019; 14:e0221397 [View Article] [PubMed]
    [Google Scholar]
  15. Adeolu M, Gupta RS. A phylogenomic and molecular marker based proposal for the division of the genus Borrelia into two genera: the emended genus Borrelia containing only the members of the relapsing fever Borrelia, and the genus Borreliella gen. nov. containing the members of the Lyme disease Borrelia (Borrelia burgdorferi sensu lato complex). Antonie Van Leeuwenhoek 2014; 105:1049–1072 [View Article] [PubMed]
    [Google Scholar]
  16. Barbour AG, Gupta RS. The Family Borreliaceae (Spirochaetales), a diverse group in two genera of tick-borne spirochetes of mammals, birds, and reptiles. J Med Entomol 2021; 58:1513–1524 [View Article] [PubMed]
    [Google Scholar]
  17. Margos G, Marosevic D, Cutler S, Derdakova M, Diuk-Wasser M et al. There is inadequate evidence to support the division of the genus Borrelia. Int J Syst Evol Microbiol 2017; 67:1081–1084 [View Article] [PubMed]
    [Google Scholar]
  18. Margos G, Wormser GP, Schwartz I, Markowicz M, Henningsson AJ et al. Evidence of taxonomic bias in public databases: the example of the genus Borrelia. Ticks Tick Borne Dis 2022; 13:101994 [View Article] [PubMed]
    [Google Scholar]
  19. van Dam AP, Kuiper H, Vos K, Widjojokusumo A, de Jongh BM et al. Different genospecies of Borrelia burgdorferi are associated with distinct clinical manifestations of Lyme borreliosis. Clin Infect Dis 1993; 17:708–717 [View Article] [PubMed]
    [Google Scholar]
  20. Purser JE, Norris SJ. Correlation between plasmid content and infectivity in Borrelia burgdorferi. Proc Natl Acad Sci USA 2000; 97:13865–13870 [View Article] [PubMed]
    [Google Scholar]
  21. Tijsse-Klasen E, Pandak N, Hengeveld P, Takumi K, Koopmans MPG et al. Ability to cause erythema migrans differs between Borrelia burgdorferi sensu lato isolates. Parasit Vectors 2013; 6:23 [View Article] [PubMed]
    [Google Scholar]
  22. Cerar T, Strle F, Stupica D, Ruzic-Sabljic E, McHugh G et al. Differences in Genotype, Clinical Features, and inflammatory potential of Borrelia burgdorferi sensu stricto strains from Europe and the United States. Emerg Infect Dis 2016; 22:818–827 [View Article] [PubMed]
    [Google Scholar]
  23. Lemieux JE, Huang W, Hill N, Cerar T, Freimark L et al. Whole genome sequencing of human Borrelia burgdorferi isolates reveals linked blocks of accessory genome elements located on plasmids and associated with human dissemination. PLOS Pathog 2023; 19:e1011243 [View Article] [PubMed]
    [Google Scholar]
  24. Brisson D, Zhou W, Jutras BL, Casjens S, Stevenson B. Distribution of cp32 prophages among Lyme disease-causing spirochetes and natural diversity of their lipoprotein-encoding erp loci. Appl Environ Microbiol 2013; 79:4115–4128 [View Article] [PubMed]
    [Google Scholar]
  25. Norris SJ. vls Antigenic Variation Systems of Lyme Disease Borrelia: Eluding Host Immunity through both Random, Segmental Gene Conversion and Framework Heterogeneity. In Mobile DNA John Wiley & Sons, Ltd; 2015 pp 471–489 [View Article]
    [Google Scholar]
  26. Fraser CM, Casjens S, Huang WM, Sutton GG, Clayton R et al. Genomic sequence of a Lyme disease spirochaete, Borrelia burgdorferi. Nature 1997; 390:580–586 [View Article]
    [Google Scholar]
  27. Casjens SR, Gilcrease EB, Vujadinovic M, Mongodin EF, Luft BJ et al. Plasmid diversity and phylogenetic consistency in the Lyme disease agent Borrelia burgdorferi. BMC Genomics 2017; 18:165 [View Article]
    [Google Scholar]
  28. Margos G, Hepner S, Mang C, Marosevic D, Reynolds SE et al. Lost in plasmids: next generation sequencing and the complex genome of the tick-borne pathogen Borrelia burgdorferi. BMC Genomics 2017; 18:422 [View Article] [PubMed]
    [Google Scholar]
  29. Casjens SR, Di L, Akther S, Mongodin EF, Luft BJ et al. Primordial origin and diversification of plasmids in Lyme disease agent bacteria. BMC Genomics 2018; 19:218 [View Article]
    [Google Scholar]
  30. Casjens SR, Mongodin EF, Qiu W-G, Luft BJ, Schutzer SE et al. Genome Stability of Lyme disease spirochetes: comparative genomics of Borrelia burgdorferi plasmids. PLoS One 2012; 7:e33280 [View Article]
    [Google Scholar]
  31. Chaconas G, Norris SJ. Peaceful coexistence amongst Borrelia plasmids: getting by with a little help from their friends?. Plasmid 2013; 70:161–167 [View Article]
    [Google Scholar]
  32. Hinnebusch J, Barbour AG. Linear plasmids of Borrelia burgdorferi have a telomeric structure and sequence similar to those of a eukaryotic virus. J Bacteriol 1991; 173:7233–7239 [View Article]
    [Google Scholar]
  33. Chaconas G, Stewart PE, Tilly K, Bono JL, Rosa P. Telomere resolution in the Lyme disease spirochete. EMBO J 2001; 20:3229–3237 [View Article] [PubMed]
    [Google Scholar]
  34. Corona A, Schwartz I. Borrelia burgdorferi: carbon metabolism and the tick-mammal enzootic cycle. Microbiol Spectr 2015; 3: [View Article] [PubMed]
    [Google Scholar]
  35. Fraser CM, Eisen JA, Nelson KE, Paulsen IT, Salzberg SL. The value of complete microbial genome sequencing (You Get What You Pay For). J Bacteriol 2002; 184:6403–6405 [View Article]
    [Google Scholar]
  36. Wick RR, Judd LM, Holt KE. Assembling the perfect bacterial genome using oxford nanopore and illumina sequencing. PLoS Comput Biol 2023; 19:e1010905 [View Article] [PubMed]
    [Google Scholar]
  37. Grimm D, Elias AF, Tilly K, Rosa PA. Plasmid stability during in vitro propagation of Borrelia burgdorferi assessed at a clonal level. Infect Immun 2003; 71:3138–3145 [View Article] [PubMed]
    [Google Scholar]
  38. Kingry LC, Replogle A, Batra D, Rowe LA, Sexton C et al. Toward a complete North American Borrelia miyamotoi genome. Genome Announc 2017; 5:01557–16 [View Article] [PubMed]
    [Google Scholar]
  39. De Coster W, Rademakers R. NanoPack2: population-scale evaluation of long-read sequencing data. Bioinformatics 2023; 39:btad311 [View Article] [PubMed]
    [Google Scholar]
  40. De Coster W, D’Hert S, Schultz DT, Cruts M, Van Broeckhoven C. NanoPack: visualizing and processing long-read sequencing data. Bioinformatics 2018; 34:2666–2669 [View Article] [PubMed]
    [Google Scholar]
  41. Kolmogorov M, Yuan J, Lin Y, Pevzner PA. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol 2019; 37:540–546 [View Article]
    [Google Scholar]
  42. Kolmogorov M, Bickhart DM, Behsaz B, Gurevich A, Rayko M et al. metaFlye: scalable long-read metagenome assembly using repeat graphs. Nat Methods 2020; 17:1103–1110 [View Article]
    [Google Scholar]
  43. Andrews S, Krueger F, Seconds-Pichon A, Biggins F, FastQC WS. Andrews S, Krueger F, Seconds-Pichon A, Biggins F, Wingett S. FastQC. A quality control tool for high throughput sequence data. Babraham Institute; 2015 https://www.bioinformatics.babraham.ac.uk/projects/fastqc/
  44. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics 2014; 30:2114–2120 [View Article]
    [Google Scholar]
  45. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLOS Comput Biol 2017; 13:e1005595 [View Article] [PubMed]
    [Google Scholar]
  46. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 2018; 34:3094–3100 [View Article] [PubMed]
    [Google Scholar]
  47. Li H, Durbin R. Fast and accurate short read alignment with burrows–wheeler transform. Bioinformatics 2009; 25:1754–1760 [View Article]
    [Google Scholar]
  48. Wick RR, Holt KE. Polypolish: short-read polishing of long-read bacterial genome assemblies. PLOS Comput Biol 2022; 18:e1009802 [View Article] [PubMed]
    [Google Scholar]
  49. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 2014; 30:2068–2069 [View Article] [PubMed]
    [Google Scholar]
  50. Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015; 31:3210–3212 [View Article] [PubMed]
    [Google Scholar]
  51. Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 2015; 25:1043–1055 [View Article] [PubMed]
    [Google Scholar]
  52. Rhie A, Walenz BP, Koren S, Phillippy AM. Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies. Genome Biol 2020; 21:245 [View Article] [PubMed]
    [Google Scholar]
  53. Jain C, Rodriguez-R LM, Phillippy AM, Konstantinidis KT, Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun 2018; 9:5114 [View Article] [PubMed]
    [Google Scholar]
  54. Birolo G, Telatin A. BamToCov: an efficient toolkit for sequence coverage calculations. Bioinformatics 2022; 38:2617–2618 [View Article] [PubMed]
    [Google Scholar]
  55. Kuleshov KV, Margos G, Fingerle V, Koetsveld J, Goptar IA et al. Whole genome sequencing of Borrelia miyamotoi isolate Izh-4: reference for a complex bacterial genome. BMC Genomics 2020; 21:16 [View Article] [PubMed]
    [Google Scholar]
  56. Jones P, Binns D, Chang H-Y, Fraser M, Li W et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 2014; 30:1236–1240 [View Article] [PubMed]
    [Google Scholar]
  57. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 2013; 30:772–780 [View Article] [PubMed]
    [Google Scholar]
  58. Price MN, Dehal PS, Arkin AP. FastTree 2--approximately maximum-likelihood trees for large alignments. PLoS One 2010; 5:e9490 [View Article] [PubMed]
    [Google Scholar]
  59. Tonkin-Hill G, MacAlasdair N, Ruis C, Weimann A, Horesh G et al. Producing polished prokaryotic pangenomes with the panaroo pipeline. Genome Biol 2020; 21:180 [View Article] [PubMed]
    [Google Scholar]
  60. Margos G, Gatewood AG, Aanensen DM, Hanincová K, Terekhova D et al. MLST of housekeeping genes captures geographic population structure and suggests a European origin of Borrelia burgdorferi. Proc Natl Acad Sci USA 2008; 105:8730–8735 [View Article] [PubMed]
    [Google Scholar]
  61. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res 2015; 43:e15 [View Article] [PubMed]
    [Google Scholar]
  62. Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009; 25:1972–1973 [View Article] [PubMed]
    [Google Scholar]
  63. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD et al. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol Biol Evol 2020; 37:1530–1534 [View Article] [PubMed]
    [Google Scholar]
  64. Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res 2018; 3:124 [View Article] [PubMed]
    [Google Scholar]
  65. Barbour AG, Travinsky B. Evolution and distribution of the ospC Gene, a transferable serotype determinant of Borrelia burgdorferi. mBio 2010; 1:e00153-10 [View Article] [PubMed]
    [Google Scholar]
  66. Wang Y, Tang H, Debarry JD, Tan X, Li J et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res 2012; 40:e49 [View Article] [PubMed]
    [Google Scholar]
  67. Gu Z, Gu L, Eils R, Schlesner M, Brors B. circlize Implements and enhances circular visualization in R. Bioinformatics 2014; 30:2811–2812 [View Article] [PubMed]
    [Google Scholar]
  68. Lovell JT, Sreedasyam A, Schranz ME, Wilson M, Carlson JW et al. GENESPACE tracks regions of interest and gene copy number variation across multiple genomes. Elife 2022; 11:e78526 [View Article] [PubMed]
    [Google Scholar]
  69. Emms DM, Kelly S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol 2019; 20:238 [View Article] [PubMed]
    [Google Scholar]
  70. Hepner S, Kuleshov K, Tooming-Kunderud A, Alig N, Gofton A et al. A high fidelity approach to assembling the complex Borrelia genome. BMC Genomics 2023; 24:401 [View Article] [PubMed]
    [Google Scholar]
  71. Gofton AW, Margos G, Fingerle V, Hepner S, Loh S-M et al. Genome-wide analysis of Borrelia turcica and “Candidatus Borrelia tachyglossi” shows relapsing fever-like genomes with unique genomic links to Lyme disease Borrelia. Infect Genet Evol 2018; 66:72–81 [View Article] [PubMed]
    [Google Scholar]
  72. Margos G, Fedorova N, Becker NS, Kleinjan JE, Marosevic D et al. Borrelia maritima sp. nov., a novel species of the Borrelia burgdorferi sensu lato complex, occupying a basal position to North American species. Int J Syst Evol Microbiol 2020; 70:849–856 [View Article] [PubMed]
    [Google Scholar]
  73. Kneubehl AR, Krishnavajhala A, Leal SM, Replogle AJ, Kingry LC et al. Comparative genomics of the Western Hemisphere soft tick-borne relapsing fever borreliae highlights extensive plasmid diversity. BMC Genomics 2022; 23:410 [View Article] [PubMed]
    [Google Scholar]
  74. Kuleshov KV, Koetsveld J, Goptar IA, Markelov ML, Kolyasnikova NM et al. Whole-genome sequencing of six Borrelia miyamotoi clinical strains isolated in Russia. Genome Announc 2018; 6:01424–17 [View Article] [PubMed]
    [Google Scholar]
  75. Becker NS, Rollins RE, Nosenko K, Paulus A, Martin S et al. High conservation combined with high plasticity: genomics and evolution of Borrelia bavariensis. BMC Genomics 2020; 21:702 [View Article] [PubMed]
    [Google Scholar]
  76. Mongodin EF, Casjens SR, Bruno JF, Xu Y, Drabek E et al. Inter- and intra-specific pan-genomes of Borrelia burgdorferi sensu lato: genome stability and adaptive radiation. BMC Genomics 2013; 14:693 [View Article]
    [Google Scholar]
  77. Livey I, Gibbs CP, Schuster R, Dorner F. Evidence for lateral transfer and recombination in OspC variation in Lyme disease Borrelia. Mol Microbiol 1995; 18:257–269 [View Article] [PubMed]
    [Google Scholar]
  78. Eggers CH, Gray CM, Preisig AM, Glenn DM, Pereira J et al. Phage-mediated horizontal gene transfer of both prophage and heterologous DNA by ϕBB-1, a bacteriophage of Borrelia burgdorferi. Pathogens Dis 2016; 74:ftw107 [View Article]
    [Google Scholar]
  79. Iyer R, Caimano MJ, Luthra A, Axline D Jr, Corona A et al. Stage-specific global alterations in the transcriptomes of Lyme disease spirochetes during tick feeding and following mammalian host adaptation. Mol Microbiol 2015; 95:509–538 [View Article] [PubMed]
    [Google Scholar]
  80. Schwan TG, Burgdorfer W, Garon CF. Changes in infectivity and plasmid profile of the lyme disease spirochete, Borrelia burgdorferi, as a result of in vitro cultivation. Infect Immun 1988; 56:1831–1836
    [Google Scholar]
  81. Casjens S. Borrelia genomes in the year 2000. J Mol Microbiol Biotechnol 2000; 2:401–410 [PubMed]
    [Google Scholar]
  82. Biškup UG, Strle F, Ružić-Sabljić E. Loss of plasmids of Borrelia burgdorferi sensu lato during prolonged in vitro cultivation. Plasmid 2011; 66:1–6 [View Article] [PubMed]
    [Google Scholar]
  83. Margos G, Hepner S, Mang C, Sing A, Liebl B et al. Completed genome sequences of Borrelia burgdorferi Sensu Stricto B31(NRZ) and closely related patient isolates from Europe. Genome Announc 2017; 5:e00637-17 [View Article] [PubMed]
    [Google Scholar]
  84. Casjens S, van Vugt R, Tilly K, Rosa PA, Stevenson B. Homology throughout the multiple 32-kilobase circular plasmids present in Lyme disease spirochetes. J Bacteriol 1997; 179:217–227 [View Article] [PubMed]
    [Google Scholar]
/content/journal/mgen/10.1099/mgen.0.001606
Loading
/content/journal/mgen/10.1099/mgen.0.001606
Loading

Data & Media loading...

Supplements

Supplementary material 1

PDF

Supplementary material 2

EXCEL
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An error occurred
Approval was partially successful, following selected items could not be processed due to error