Comparative virulome analysis of four Staphylococcus epidermidis strains from human skin and platelet concentrates using whole genome sequencing.

Staphylococcus epidermidis is one of the predominant bacterial contaminants in platelet concentrates (PCs), a blood component used to treat bleeding disorders. PCs are a unique niche that triggers biofilm formation, the main pathomechanism of S. epidermidis infections. We performed whole genome sequencing of four S. epidermidis strains isolated from skin of healthy human volunteers (AZ22 and AZ39) and contaminated PCs (ST10002 and ST11003) to unravel phylogenetic relationships and decipher virulence mechanisms compared to 24 complete S. epidermidis genomes in GenBank. AZ39 and ST11003 formed a separate unique lineage with strains 14.1 .R1 and SE95, while AZ22 formed a cluster with 1457 and ST10002 closely grouped with FDAAGOS_161. The four isolates were assigned to sequence types ST1175, ST1174, ST73 and ST16, respectively. All four genomes exhibited biofilm-associated genes ebh, ebp, sdrG, sdrH and atl. Additionally, AZ22 had sdrF and aap, whereas ST10002 had aap and icaABCDR. Notably, AZ39 possesses truncated ebh and sdrG and harbours a toxin-encoding gene. All isolates carry multiple antibiotic resistance genes conferring resistance to fosfomycin (fosB), ß-lactams (blaZ) and fluoroquinolones (norA). This study reveales a unique lineage for S. epidermidis and provides insight into the genetic basis of virulence and antibiotic resistance in transfusion-associated S. epidermidis strains.

Complete genomes of Clostridium botulinum type B(F) isolates associated with a 1995 foodborne botulism outbreak from commercial pâté reveals a recombination event disrupting the ntnh gene.

Foodborne botulism is a neuroparalytic disease caused by ingestion of foods contaminated with botulinum neurotoxin (BoNT), produced by Clostridium botulinum. In 1995 a husband and wife from Québec, Canada, were hospitalized for several months with prolonged muscle paralysis after ingesting a commercial pâté de campagne. Examination of faecal samples from both patients and the pâté produced viable Group I (proteolytic) C. botulinum type B from each of the three samples. Whole genome sequencing revealed that all three isolates contain identical bont/B5 and bont/F2 genes encoded on a plasmid. Both faecal isolate genomes were identical in chromosome and plasmid length, as well as gene content. The genome of the pâté isolate was nearly identical to that of the faecal isolates with the notable difference of a missing 13-gene insertion on the bont/B5 cluster disrupting the ntnh gene. Examination of the insertion revealed several mobile genetic elements that participate in recombination.

Automated centrifugal microfluidic system for the preparation of adaptor-ligated sequencing libraries.

The intensive workload associated with the preparation of high-quality DNA libraries remains a key obstacle toward widespread deployment of sequencing technologies in remote and resource-limited areas. We describe the development of single-use microfluidic devices driven by an advanced pneumatic centrifugal microfluidic platform, the PowerBlade, to automate the preparation of Illumina-compatible libraries based on adaptor ligation methodology. The developed on-chip workflow includes enzymatic DNA fragmentation coupled to end-repair, adaptor ligation, first DNA cleanup, PCR amplification, and second DNA cleanup. This complex workflow was successfully integrated into simple thermoplastic microfluidic devices that are amenable to mass production with injection molding. The system was validated by preparing, on chip, libraries from a mixture of genomic DNA extracted from three common foodborne pathogens (Listeria monocytogenes, Escherichia coli and Salmonella enterica serovar Typhimurium) and comparing them with libraries made via a manual procedure. The two types of libraries were found to exhibit similar quality control metrics (including genome coverage, assembly, and relative abundances) and led to nearly uniform coverage independent of GC content. This microfluidic technology offers a time-saving and cost-effective alternative to manual procedures and robotic-based automation, making it suitable for deployment in remote environments where technical expertise and resources might be scarce. Specifically, it facilitates field practices that involve mid- to low-throughput sequencing, such as tasks related to foodborne pathogen detection, characterization, and microbial profiling.

Foodborne Botulism, Canada, 2006-20211.

During 2006-2021, Canada had 55 laboratory-confirmed outbreaks of foodborne botulism, involving 67 cases. The mean annual incidence was 0.01 case/100,000 population. Foodborne botulism in Indigenous communities accounted for 46% of all cases, which is down from 85% of all cases during 1990-2005. Among all cases, 52% were caused by botulinum neurotoxin type E, but types A (24%), B (16%), F (3%), and AB (1%) also occurred; 3% were caused by undetermined serotypes. Four outbreaks resulted from commercial products, including a 2006 international outbreak caused by carrot juice. Hospital data indicated that 78% of patients were transferred to special care units and 70% required mechanical ventilation; 7 deaths were reported. Botulinum neurotoxin type A was associated with much longer hospital stays and more time spent in special care than types B or E. Foodborne botulism often is misdiagnosed. Increased clinician awareness can improve diagnosis, which can aid epidemiologic investigations and patient treatment.

Targeted metagenomics using bait-capture to detect antibiotic resistance genes in retail meat and seafood.

Metagenomics analysis of foods has the potential to provide comprehensive data on the presence and prevalence of antimicrobial resistance (AMR) genes in the microbiome of foods. However, AMR genes are generally present in low abundance compared to other bacterial genes in the food microbiome and consequently require multiple rounds of in-depth sequencing for detection. Here, a metagenomics approach, using bait-capture probes targeting antimicrobial resistance and plasmid genes, is used to characterize the resistome and plasmidome of retail beef, chicken, oyster, shrimp, and veal enrichment cultures (n = 15). Compared to total shotgun metagenomics, bait-capture required approximately 40-fold fewer sequence reads to detect twice the number of AMR gene classes, AMR gene families, and plasmid genes across all sample types. For the detection of critically important extended spectrum beta-lactamase (ESBL) genes the bait capture method had a higher overall positivity rate (44%) compared to shotgun metagenomics (26%), and a culture-based method (29%). Overall, the results support the use of bait-capture for the identification of low abundance genes such as AMR genes from food samples.

Draft Genomes of 83 Vibrio Isolates from Fresh Canadian Mollusks.

A subset of Vibrio spp. isolated from fresh Canadian mollusks (2014 to 2018) were selected for sequencing based on antimicrobial resistance profiles. The resulting de novo draft genomes include 38 Vibrio alginolyticus, 32 V. diabolicus, 10 V. parahaemolyticus, 1 V. cholerae, 1 V. ordalii, and 1 Vibrio sp. isolate.

Draft Genomes of 75 Vibrio Isolates from Retail Mollusks Collected in Canada between 2014 and 2018.

Vibrio spp. isolated from fresh retail mollusk samples were selected for sequencing based on their antimicrobial resistance burden. The de novo genomes include those for Vibrio alginolyticus (n = 48), V. diabolicus (n = 15), V. parahaemolyticus (n = 3), V. cholerae (n = 2), V. metoecus (n = 1), V. vulnificus (n = 1), V. fluvialis (n = 1), and unidentified Vibrio spp. (n = 4).

Draft Genomes of 92 Vibrio Isolates from Warm-Water Shrimps Imported into Canada between 2009 and 2019.

Vibrio spp. were isolated from raw shrimps imported into Canada (2009 to 2019). A total of 92 isolates with various multidrug resistance profiles were sequenced, including 59 V. parahaemolyticus, 12 V. diabolicus, 10 V. cholerae, 7 V. alginolyticus, 1 V. campbellii, 1 V. harveyi, 1 V. owensii, and 1 V. vulnificus isolate.

Choice of DNA extraction method affects detection of bacterial taxa from retail chicken breast.

BACKGROUND: Sequence-based methods for the detection of bacteria such as 16S rRNA amplicon sequencing and metagenomics can provide a comprehensive view of the bacterial microbiome of food. These methods rely on the detection of gene sequences to indicate the presence of viable bacteria. This indirect form of detection can be prone to experimental artefacts. Sample handling and processing are key sources of variation that require standard approaches. Extracting sufficient quantities of high quality DNA from food matrices is challenging because target bacterial species are usually minor components of the microbiota and foods contain an array of compounds that are inhibitory to downstream DNA applications. Here, three DNA extraction methods are compared for their ability to extract high quality bacterial DNA from retail chicken breast rinses, with or without enrichment. Method performance was assessed by comparing ease of use, DNA yield, DNA quality, PCR amplicon yield, and the detection of bacterial taxa by 16S rRNA amplicon sequencing. RESULTS: All three DNA extraction methods yielded DNA of sufficient quantity and quality to perform quantitative PCR and 16S rRNA amplicon sequencing. The extraction methods differed in ease of use, with the two commercial kits (PowerFood, PowerSoil) offering considerable time and cost savings over a hybrid method that used laboratory reagents for lysis and commercial column based kits for further purification. Bacterial richness as determined by 16S rRNA amplicon sequencing was similar across the three DNA extraction methods. However, differences were noted in the relative abundance of bacterial taxa, with significantly higher abundance of Gram-positive genera detected in the DNA samples prepared using the PowerFood DNA extraction kit. CONCLUSION: The choice of DNA extraction method can affect the detection of bacterial taxa by 16S rRNA amplicon sequencing in chicken meat rinses. Investigators should be aware of this procedural bias and select methods that are fit for the purposes of their investigation.

Characterization of Atypical Shiga Toxin Gene Sequences and Description of Stx2j, a New Subtype.

Shiga toxin (Stx) is the definitive virulence factor of Shiga toxin-producing Escherichia coli (STEC). Stx variants are currently organized into a taxonomic system of three Stx1 (a, c, and d) and seven Stx2 (a, b, c, d, e, f, and g) subtypes. In this study, seven STEC isolates from food and clinical samples possessing stx2 sequences that do not fit current Shiga toxin taxonomy were identified. Genome assemblies of the STEC strains were created from Oxford Nanopore and Illumina sequence data. The presence of atypical stx2 sequences was confirmed by Sanger sequencing, as were Stx2 expression and cytotoxicity. A strain of O157:H7 was found to possess stx1a and a truncated stx2a, which were originally misidentified as an atypical stx2. Two strains possessed unreported variants of Stx2a (O8:H28) and Stx2b (O146:H21). In four of the strains, we found three Stx subtypes that are not included in the current taxonomy. Stx2h (O170:H18) was identified in a Canadian sprout isolate; this subtype has only previously been reported in STEC from Tibetan Marmots. Stx2o (O85:H1) was identified in a clinical isolate. Finally, Stx2j (O158:H23 and O33:H14) was found in lettuce and clinical isolates. The results of this study expand the number of known Stx subtypes, the range of STEC serotypes, and isolation sources in which they may be found. The presence of the Stx2j and Stx2o in clinical isolates of STEC indicates that strains carrying these variants are potential human pathogens.

Whole-Genome Sequences of Vibrio Species from Warm-Water Shrimps Imported into Canada: Detection of Genetic Elements Associated with Antimicrobial Resistance and Potential Mobilizing Capacities.

We present draft genome sequences of Vibrio species (Vibrio alginolyticus, Vibrio cholerae, and two Vibrio parahaemolyticus strains) that were isolated from warm-water shrimps imported into Canada. All four isolates harbor genetic elements associated with antimicrobial resistance (AMR), including mobile genetic elements that can promote horizontal transfer of AMR genes.

Complete Genome Sequence of Staphylococcus aureus CI/BAC/25/13/W, Isolated from Contaminated Platelet Concentrates in England.

We present the genome sequence of Staphylococcus aureus CI/BAC/25/13/W, which was isolated in 2013 as a contaminant of a platelet concentrate with abnormal clotting at the National Health Service Blood and Transplant. Assessment of the genome sequence showed the presence of one chromosome (2,719,347 bp) and one plasmid (1,533 bp).

Complete Genome Sequence of Staphylococcus aureus PS/BAC/169/17/W, Isolated from a Contaminated Platelet Concentrate in England.

We report the genome sequence of Staphylococcus aureus PS/BAC/169/17/W, which was isolated in 2017 from a contaminated platelet concentrate at the National Health Service Blood and Transplant. Assessment of the genome sequence of this strain showed the presence of a 2,753,746-bp chromosome and a plasmid of 2,762 bp.

Genome Sequence of Staphylococcus aureus Strain PS/BAC/317/16/W, Isolated from Contaminated Platelet Concentrates in England.

We present the genome sequence of Staphylococcus aureus strain PS/BAC/317/16/W, which was isolated from contaminated platelet concentrates by the National Health Service Blood and Transplant in England (2017). Genome sequence analysis revealed the presence of one chromosome (2,665,983 bp) and two plasmids (4,265 bp and 2,921 bp) in this strain.

Genome Sequence of Staphylococcus aureus Strain CBS2016-05, Isolated from Contaminated Platelet Concentrates in Canada.

We present the genome sequence of Staphylococcus aureus strain CBS2016-05, which was isolated from contaminated platelet concentrates by Canadian Blood Services in 2016. This strain caused a septic reaction in an acute leukemia patient. Genome sequence analysis revealed the presence of one chromosome (2,766,936 bp) and one plasmid (36,441 bp).

Complete Genome Sequence of Clostridium botulinum CJ0611A1, a Type A(B) Isolate Associated with an International Outbreak of Botulism from Commercial Carrot Juice.

We report the complete genome (3.9-Mb chromosome, 5.9-kb plasmid) of Clostridium botulinum CJ0611A1, a type A(B) strain isolated from carrot juice distributed in Canada and linked to an international 2006 foodborne botulism outbreak. This strain encodes a full-length bont/A1 gene and a truncated bont/B gene.

Genomic analysis of human noroviruses using combined Illumina-Nanopore data.

Whole-genome sequence analysis of noroviruses is routinely performed by employing a metagenomic approach. While this methodology has several advantages, such as allowing for the examination of co-infection, it has some limitations, such as the requirement of high viral load to achieve full-length or near full-length genomic sequences. In this study, we used a pre-amplification step to obtain full-length genomic amplicons from 39 Canadian GII isolates, followed by deep sequencing on Illumina and Oxford Nanopore platforms. This approach significantly reduced the required viral titre to obtain full-genome coverage. Herein, we compared the coverage and sequences obtained by both platforms and provided an in-depth genomic analysis of the obtained sequences, including the presence of single-nucleotide variants and recombination events.

Shared genome analyses of notable listeriosis outbreaks, highlighting the critical importance of epidemiological evidence, input datasets and interpretation criteria.

The persuasiveness of genomic evidence has pressured scientific agencies to supplement or replace well-established methodologies to inform public health and food safety decision-making. This study of 52 epidemiologically defined Listeria monocytogenes isolates, collected between 1981 and 2011, including nine outbreaks, was undertaken (1) to characterize their phylogenetic relationship at finished genome-level resolution, (2) to elucidate the underlying genetic diversity within an endemic subtype, CC8, and (3) to re-evaluate the genetic relationship and epidemiology of a CC8-delimited outbreak in Canada in 2008. Genomes representing Canadian Listeria outbreaks between 1981 and 2010 were closed and manually annotated. Single nucleotide variants (SNVs) and horizontally acquired traits were used to generate phylogenomic models. Phylogenomic relationships were congruent with classical subtyping and epidemiology, except for CC8 outbreaks, wherein the distribution of SNV and prophages revealed multiple co-evolving lineages. Chronophyletic reconstruction of CC8 evolution indicates that prophage-related genetic changes among CC8 strains manifest as PFGE subtype reversions, obscuring the relationship between CC8 isolates, and complicating the public health interpretation of subtyping data, even at maximum genome resolution. The size of the shared genome interrogated did not change the genetic relationship measured between highly related isolates near the tips of the phylogenetic tree, illustrating the robustness of these approaches for routine public health applications where the focus is recent ancestry. The possibility exists for temporally and epidemiologically distinct events to appear related even at maximum genome resolution, highlighting the continued importance of epidemiological evidence.

SNVPhyl: a single nucleotide variant phylogenomics pipeline for microbial genomic epidemiology.

The recent widespread application of whole-genome sequencing (WGS) for microbial disease investigations has spurred the development of new bioinformatics tools, including a notable proliferation of phylogenomics pipelines designed for infectious disease surveillance and outbreak investigation. Transitioning the use of WGS data out of the research laboratory and into the front lines of surveillance and outbreak response requires user-friendly, reproducible and scalable pipelines that have been well validated. Single Nucleotide Variant Phylogenomics (SNVPhyl) is a bioinformatics pipeline for identifying high-quality single-nucleotide variants (SNVs) and constructing a whole-genome phylogeny from a collection of WGS reads and a reference genome. Individual pipeline components are integrated into the Galaxy bioinformatics framework, enabling data analysis in a user-friendly, reproducible and scalable environment. We show that SNVPhyl can detect SNVs with high sensitivity and specificity, and identify and remove regions of high SNV density (indicative of recombination). SNVPhyl is able to correctly distinguish outbreak from non-outbreak isolates across a range of variant-calling settings, sequencing-coverage thresholds or in the presence of contamination. SNVPhyl is available as a Galaxy workflow, Docker and virtual machine images, and a Unix-based command-line application. SNVPhyl is released under the Apache 2.0 license and available at http://snvphyl.readthedocs.io/ or at https://github.com/phac-nml/snvphyl-galaxy.

Neptune: a bioinformatics tool for rapid discovery of genomic variation in bacterial populations.

The ready availability of vast amounts of genomic sequence data has created the need to rethink comparative genomics algorithms using 'big data' approaches. Neptune is an efficient system for rapidly locating differentially abundant genomic content in bacterial populations using an exact k-mer matching strategy, while accommodating k-mer mismatches. Neptune's loci discovery process identifies sequences that are sufficiently common to a group of target sequences and sufficiently absent from non-targets using probabilistic models. Neptune uses parallel computing to efficiently identify and extract these loci from draft genome assemblies without requiring multiple sequence alignments or other computationally expensive comparative sequence analyses. Tests on simulated and real datasets showed that Neptune rapidly identifies regions that are both sensitive and specific. We demonstrate that this system can identify trait-specific loci from different bacterial lineages. Neptune is broadly applicable for comparative bacterial analyses, yet will particularly benefit pathogenomic applications, owing to efficient and sensitive discovery of differentially abundant genomic loci. The software is available for download at: http://github.com/phac-nml/neptune.