Open access
Public Health
Announcement
5 January 2024

Draft whole-genome sequences of three Borrelia burgdorferi isolates from Western Canada

ABSTRACT

Whole-genome sequences are presented for three Borrelia burgdorferi, a causative agent of Lyme disease in North America, isolated from Ixodes pacificus ticks collected in British Columbia, Canada. Shotgun DNA libraries were prepared with Illumina DNA Prep and sequenced using the MiniSeq platform. Genome assemblies enabled multilocus sequence typing and ospC typing.

ANNOUNCEMENT

Borrelia burgdorferi sensu stricto is the predominant causative agent of Lyme disease in North America (1), the most common vector-borne disease in the region (13). In Canada, the highest disease risk is associated with the black-legged tick Ixodes scapularis in central to eastern provinces (47). B. burgdorferi is endemic but less prevalent in southern British Columbia, where it is spread primarily by the tick Ixodes pacificus (4). To investigate B. burgdorferi genome composition in Western Canada, we sequenced the whole genomes of isolates 21T 0066, 22T 0045, and 22T 0109.
The strains were cultured from three I. pacificus ticks found on humans. Strains 21T 0066 and 22T 0045 were isolated from ticks collected in Victoria (2021 and 2022, respectively), and strain 22T 0109 was isolated from a tick collected on Galiano Island (2022). Ticks were submitted to the British Columbia Center for Disease Control, surface disinfected using 10% hydrogen peroxide and 70% ethanol, and then dissected. The contents of their midgut were cultured at 34°C in modified Barbour-Stonner-Kelly medium containing kanamycin and rifampicin at concentrations of 10 µg/mL and 50 µg/mL, respectively. Cultures were checked weekly for spirochetes using dark-field microscopy, followed by quantitative PCR confirmation of B. burgdorferi using primers and probe for 23S rDNA, as previously described (8). Culturing was performed for periods up to 35 days, and positive cultures were then stored at −80°C before being shipped overnight on −80°C ice packs. Total genomic DNA was isolated from culture aliquots using the DNeasy UltraClean Microbial kit (Qiagen, Venlo, Netherlands). DNA extracts were quantified using the Qubit 2.0 fluorometer and dsDNA High Sensitivity assay kit (Thermo Fisher, Waltham, MA).
DNA sequencing libraries were prepared using Illumina DNA Prep (Illumina, CA, USA) and assumed to have an average fragment size of 600 bp. Final libraries were quantified using the NEBNext Library Quant Kit for Illumina (New England BioLabs; Ipswich, MA, USA) then normalized to 10 nM prior to pooling, diluting, and denaturing. Sequencing was performed on the Illumina MiniSeq platform using 300-cycle paired-end chemistry.
Sequence assembly was conducted entirely within BV-BRC Comprehensive Genome Analysis (9); default parameters were used at all stages of the workflow except where otherwise noted. Raw sequencing data were uploaded to BV-BRC in fastq format, and then Trim Galore (0.6.5dev) removed adapters and low-quality sequence per the “true” trim setting (10). Contigs were assembled by Unicycler (v0.4.8) (11) and then were filtered by default minimum thresholds of 300-bp length and fivefold coverage, followed by two polishing iterations using PILON (v1.23) (12). Assembly quality was assessed using QUAST (v5.2.0) (13). Genomes were subsequently annotated by the NCBI Prokaryotic Genome Annotation Pipeline (v6.5) (14); completeness assessed by CheckM (v1.2.2) (15) ranged from 99.4% to 99.7%.
Multilocus sequence typing (MLST) was performed using the PubMLST database (16), classifying isolate 21T 0066 as sequence type (ST) 13 and 22T 0045 as ST6. 22T 0109 matched ST6 at all MLST positions except for a single nucleotide variant 250G>A at the nifS locus, making 22T 0109 nearest ST6. The ospC gene encodes the primary surface antigen under immune selection pressure and so demonstrates high genetic diversity (17). Assemblies were queried against ospC reference sequences as in reference (18) using the BLASTn function in BLAST+ (v2.9.0) (19) and a >92.0% nucleotide identity threshold (20). 21T 0066 is ospC type M, while both 22T 0045 and 22T 0109 are ospC type H. Sequencing, assembly, genome quality, and annotation metrics for the three isolates are presented in Table 1.
TABLE 1
TABLE 1 Sequencing, assembly, quality, and annotation metrics for Borrelia burgdorferi isolates
Isolate nameNo. of paired readsNo. of basesAverage read lengthCoverage (×)No. of contigsN50 contig size (bp)Total size (bp)Assembly G + C content (%)No. of CDSNo. of RNAsNo. of pseudo genesMultilocus sequence typeospC typeGenbank accession no.SRA Accession no.
21T 0066504,431151,329,30015083.270894,8431,223,15128.61,210413213MJAUPGQ000000000SRR25420533
22T 0045427,022128,106,60015084.4124213,0281,249,51628.11,25138286HJAUPGP000000000SRR25420532
22T 0109230,53369,159,90015046.5124164,1931,234,77228.21,23340286aHJAUPGO000000000SRR25420531
a
Isolate 22T 0109 was nearest sequence type 6. Seven of the eight loci were perfect matches, but nifS differed from the nifS sequence type 8 at 250G>A.

ACKNOWLEDGMENTS

This project was supported in part by a Natural Sciences and Engineering Research Council of Canada (NSERC) Undergraduate Student Research Award (USRA) and by Public Health Agency of Canada grant number 1920-HQ-00069.

REFERENCES

1.
Marques AR, Strle F, Wormser GP. 2021. Comparison of Lyme disease in the United States and Europe. Emerg Infect Dis 27:2017–2024.
2.
Stanek G, Wormser GP, Gray J, Strle F. 2012. Lyme Borreliosis. The Lancet 379:461–473.
3.
Bouchard C, Leonard E, Koffi JK, Pelcat Y, Peregrine A, Chilton N, Rochon K, Lysyk T, Lindsay LR, Ogden NH. 2015. The increasing risk of Lyme disease in Canada. Can Vet J 56:693–699.
4.
Ogden NH, Lindsay LR, Morshed M, Sockett PN, Artsob H. 2009. The emergence of Lyme disease in Canada. CMAJ 180:1221–1224.
5.
Leighton PA, Koffi JK, Pelcat Y, Lindsay LR, Ogden NH. 2012. Predicting the speed of tick invasion: an empirical model of range expansion for the Lyme disease vector Ixodes scapularis in Canada. J. Appl. Ecol. 49:457–464.
6.
Ogden NH, Lindsay LR, Leighton PA, Sheppard A. 2013. Predicting the rate of invasion of the agent of Lyme disease Borrelia Burgdorferi. J. Appl. Ecol. 50:510–518.
7.
Ogden NH, Trudel L, Artsob H, Barker IK, Beauchamp G, Charron DF, Drebot MA, Galloway TD, O’Handley R, Thompson RA, Lindsay LR. 2006. Ixodes scapularis ticks collected by passive surveillance in Canada: analysis of geographic distribution and infection with Lyme borreliosis agent Borrelia burgdorferi. J Med Entomol 43:600–609.
8.
Courtney JW, Kostelnik LM, Zeidner NS, Massung RF. 2004. Multiplex real-time PCR for detection of Anaplasma phagocytophilum and Borrelia burgdorferi. J Clin Microbiol 42:3164–3168.
9.
Olson RD, Assaf R, Brettin T, Conrad N, Cucinell C, Davis JJ, Dempsey DM, Dickerman A, Dietrich EM, Kenyon RW, et al. 2023. Introducing the bacterial and viral bioinformatics resource center (BV-BRC): a resource combining PATRIC, IRD and ViPR. Nucleic Acids Res 51:D678–D689.
10.
Krueger, F. (2015). TrimGalore: a wrapper tool around Cutadapt and FastQC to consistently apply quality and adapter trimming to FastQ files. 516(517).
11.
Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLOS Comput Biol 13:e1005595.
12.
Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, Cuomo CA, Zeng Q, Wortman J, Young SK, Earl AM. 2014. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9:e112963.
13.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. 2013. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29:1072–1075.
14.
Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, Ostell J. 2016. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 44:6614–6624.
15.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. 2015. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 25:1043–1055.
16.
Jolley KA, Bray JE, Maiden MCJ. 2018. Open-access bacterial population genomics: BIGSdb software, the PubMLST. org website and their applications. Wellcome Open Res 3:124.
17.
Brisson D, Dykhuizen DE. 2004. ospC diversity in Borrelia burgdorferi: different hosts are different niches. Genetics 168:713–722.
18.
Ogden NH, Margos G, Aanensen DM, Drebot MA, Feil EJ, Hanincová K, Schwartz I, Tyler S, Lindsay LR. 2011. Investigation of genotypes of Borrelia burgdorferi in Ixodes scapularis ticks collected during surveillance in Canada. Appl Environ Microbiol 77:3244–3254.
19.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. 2009. BLAST+: architecture and applications. BMC Bioinformatics 10:421.
20.
Pérez D, Kneubühler Y, Rais O, Jouda F, Gern L. 2011. Borrelia afzelii ospC genotype diversity in Ixodes ricinus questing ticks and ticks from rodents in two Lyme borreliosis endemic areas: contribution of co-feeding ticks. Ticks Tick Borne Dis 2:137–142.

Information & Contributors

Information

Published In

cover image Microbiology Resource Announcements
Microbiology Resource Announcements
Volume 13Number 215 February 2024
eLocator: e00879-23
Editor: Irene L. G. Newton, Indiana University, Bloomington, Indiana, USA
PubMed: 38179914

History

Received: 16 September 2023
Accepted: 10 December 2023
Published online: 5 January 2024

Keywords

  1. Borrelia burgdorferi
  2. Lyme disease
  3. Ixodes pacificus

Data Availability

This whole-genome sequencing project has been documented as NCBI BioProject PRJNA998043. Under this BioProject, sequence assemblies have been deposited in GenBank, and raw sequence read files have been uploaded to the NCBI sequence Read Archive (SRA). Accession numbers to these submissions for each isolate can be found in Table 1.

Contributors

Authors

Institute for Microbial Systems and Society, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Department of Biology, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Author Contributions: Formal analysis, Investigation, Writing – original draft, and Writing – review and editing.
Min-Kuang Lee
Zoonotic Diseases & Emerging Pathogens, BC Centre for Disease Control-Public Health Laboratory, Vancouver, British Columbia, Canada
Author Contributions: Conceptualization and Investigation.
Institute for Microbial Systems and Society, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Department of Biology, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Author Contribution: Formal analysis.
Jennifer N. Russell
Institute for Microbial Systems and Society, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Department of Biology, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Author Contribution: Formal analysis.
Erin Fraser
BC Centre for Disease Control, Vancouver, British Columbia, Canada
School of Population & Public Health, University of British Columbia, Vancouver, British Columbia, Canada
Author Contributions: Conceptualization and Funding acquisition.
Zoonotic Diseases & Emerging Pathogens, BC Centre for Disease Control-Public Health Laboratory, Vancouver, British Columbia, Canada
Department of Pathology and Laboratory Medicine, University of British Columbia, Vancouver, British Columbia, Canada
Author Contributions: Conceptualization and Project administration.
Institute for Microbial Systems and Society, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Department of Biology, Faculty of Science, University of Regina, Regina, Saskatchewan, Canada
Author Contributions: Conceptualization, Funding acquisition, Project administration, Supervision, and Writing – review and editing.

Editor

Irene L. G. Newton
Editor
Indiana University, Bloomington, Indiana, USA

Notes

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures

Tables

Media

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy