Nicholas Petronella

Navigating Microbiological Food Safety in the Era of Whole-Genome Sequencing

SUMMARY The epidemiological investigation of a foodborne outbreak, including identification of related cases, source attribution, and development of intervention strategies, relies heavily on the ability to subtype the etiological agent at a high enough resolution to differentiate related from nonrelated cases. Historically, several different molecular subtyping methods have been used for this purpose; however, emerging techniques, such as single nucleotide polymorphism (SNP)-based techniques, that use whole-genome sequencing (WGS) offer a resolution that was previously not possible. With WGS, unlike traditional subtyping methods that lack complete information, data can be used to elucidate phylogenetic relationships and disease-causing lineages can be tracked and monitored over time. The subtyping resolution and evolutionary context provided by WGS data allow investigators to connect related illnesses that would be missed by traditional techniques. The added advantage of data generated by WGS is that these data can also be used for secondary analyses, such as virulence gene detection, antibiotic resistance gene profiling, synteny comparisons, mobile genetic element identification, and geographic attribution. In addition, several software packages are now available to generate in silico results for traditional molecular subtyping methods from the whole-genome sequence, allowing for efficient comparison with historical databases. Metagenomic approaches using next-generation sequencing have also been successful in the detection of nonculturable foodborne pathogens. This review addresses state-of-the-art techniques in microbial WGS and analysis and then discusses how this technology can be used to help support food safety investigations. Retrospective outbreak investigations using WGS are presented to provide organism-specific examples of the benefits, and challenges, associated with WGS in comparison to traditional molecular subtyping techniques.

Choice of Reference Sequence and Assembler for Alignment of Listeria monocytogenes Short-Read Sequence Data Greatly Influences Rates of Error in SNP Analyses

The wide availability of whole-genome sequencing (WGS) and an abundance of open-source software have made detection of single-nucleotide polymorphisms (SNPs) in bacterial genomes an increasingly accessible and effective tool for comparative analyses. Thus, ensuring that real nucleotide differences between genomes (i.e., true SNPs) are detected at high rates and that the influences of errors (such as false positive SNPs, ambiguously called sites, and gaps) are mitigated is of utmost importance. The choices researchers make regarding the generation and analysis of WGS data can greatly influence the accuracy of short-read sequence alignments and, therefore, the efficacy of such experiments. We studied the effects of some of these choices, including: i) depth of sequencing coverage, ii) choice of reference-guided short-read sequence assembler, iii) choice of reference genome, and iv) whether to perform read-quality filtering and trimming, on our ability to detect true SNPs and on the frequencies of errors. We performed benchmarking experiments, during which we assembled simulated and real Listeria monocytogenes strain 08-5578 short-read sequence datasets of varying quality with four commonly used assemblers (BWA, MOSAIK, Novoalign, and SMALT), using reference genomes of varying genetic distances, and with or without read pre-processing (i.e., quality filtering and trimming). We found that assemblies of at least 50-fold coverage provided the most accurate results. In addition, MOSAIK yielded the fewest errors when reads were aligned to a nearly identical reference genome, while using SMALT to align reads against a reference sequence that is ∼0.82% distant from 08-5578 at the nucleotide level resulted in the detection of the greatest numbers of true SNPs and the fewest errors. Finally, we show that whether read pre-processing improves SNP detection depends upon the choice of reference sequence and assembler. In total, this study demonstrates that researchers should test a variety of conditions to achieve optimal results.

The Listeria monocytogenes Core-Genome Sequence Typer (LmCGST): a bioinformatic pipeline for molecular characterization with next-generation sequence data

BackgroundNext-generation sequencing provides a powerful means of molecular characterization. However, methods such as single-nucleotide polymorphism detection or whole-chromosome sequence analysis are computationally expensive, prone to errors, and are still less accessible than traditional typing methods. Here, we present the Listeria monocytogenes core-genome sequence typing method for molecular characterization. This method uses a high-confidence core (HCC) genome, calculated to ensure accurate identification of orthologs. We also developed an evolutionarily relevant nomenclature based upon phylogenetic analysis of HCC genomes. Finally, we created a pipeline (LmCGST; https://sourceforge.net/projects/lmcgst/files/) that takes in raw next-generation sequencing reads, calculates a subject HCC profile, compares it to an expandable database, assigns a sequence type, and performs a phylogenetic analysis.ResultsWe analyzed 29 high-quality, closed Listeria monocytogenes chromosome sequences and identified loci that are reliable targets for automated molecular characterization methods. We identified 1013 open-reading frames that comprise our high-confidence core (HCC) genome. We then populated a database with HCC profiles from 114 taxa. We sequenced 84 randomly selected isolates from the Listeriosis Reference Service for Canada’s collection and analysed them with the LmCGST pipeline. In addition, we generated pulsed-field gel electrophoresis, ribotyping, and in silico multi-locus sequence typing (MLST) data for the 84 isolates and compared the results to those obtained using the CGST method. We found that all of the methods yielded results that are generally congruent. However, due to the increased numbers of categories, the CGST method provides much greater discriminatory power than the other methods tested here.ConclusionsWe show that the CGST method provides increased discriminatory power relative to typing methods such as pulsed-field gel electrophoresis, ribotyping, and multi-locus sequence typing while it addresses several shortcomings of other methods of molecular characterization with next-generation sequence data. It uses discrete, well-defined groupings (types) of organisms that are phylogenetically relevant and easily interpreted. In addition, the CGST scheme can be expanded to include additional loci and HCC profiles in the future. In total, the CGST method provides an approach to the molecular characterization of Listeria monocytogenes with next-generation sequence data that is highly reproducible, easily standardized, portable, and accessible.

Draft Genome Sequences of Four Vibrio parahaemolyticus Isolates from Clinical Cases in Canada

ABSTRACT Vibrio parahaemolyticus is a leading cause of bacterial gastroenteritis following ingestion of contaminated seafood. This report presents the draft genome sequences of four clinical strains of V. parahaemolyticus isolated in Canada. All four strains lack traditional pathogenic markers and possess uniquely individual characteristics identified using other typing criteria.

Draft Genome Sequences of Two Clostridium botulinum Group II (Nonproteolytic) Type B Strains (DB-2 and KAPB-3)

ABSTRACT Clostridium botulinum is important for food safety and studies of neurotoxins associated with human botulism. We present the draft genome sequences of two strains belonging to group II type B: one collected from Pacific Ocean sediments (DB-2) and another obtained during a botulism outbreak (KAPB-3).

Characterization of the Genomic Diversity of Norovirus in Linked Patients Using a Metagenomic Deep Sequencing Approach.

Norovirus (NoV) is the leading cause of gastroenteritis worldwide. A robust cell culture system does not exist for NoV and therefore detailed characterization of outbreak and sporadic strains relies on molecular techniques. In this study, we employed a metagenomic approach that uses non-specific amplification followed by next-generation sequencing to whole genome sequence NoV genomes directly from clinical samples obtained from 8 linked patients. Enough sequencing depth was obtained for each sample to use a de novo assembly of near-complete genome sequences. The resultant consensus sequences were then used to identify inter-host nucleotide variations that occur after direct transmission, analyze amino acid variations in the major capsid protein, and provide evidence of recombination events. The analysis of intra-host quasispecies diversity was possible due to high coverage-depth. We also observed a linear relationship between NoV viral load in the clinical sample and the number of sequence reads that could be attributed to NoV. The method demonstrated here has the potential for future use in whole genome sequence analyses of other RNA viruses isolated from clinical, environmental, and food specimens.

Draft Whole-Genome Sequences of 14 Vibrio parahaemolyticus Clinical Isolates with an Ambiguous K Serogroup

ABSTRACT Vibrio parahaemolyticus is a bacterial pathogen responsible for mild to severe gastroenteritis, wound infections, and septicemia resulting from the ingestion or handling of raw or undercooked contaminated seafood. Here, we report the draft whole-genome sequences and annotations of 14 Canadian V. parahaemolyticus clinical isolates that were serologically identified as K group II using polyvalent antisera but were not specifically K serogrouped using monovalent antisera.

Draft Genome Sequences of Two Salmonella enterica Strains Isolated from Sprouted Chia and Flax Seed Powders.

ABSTRACT A 2014 foodborne salmonellosis outbreak in Canada and the United States implicated, for the first time, sprouted chia seed powder as the vehicle of transmission. Here, we report the draft whole genome sequences of two Salmonella enterica strains isolated from sprouted powders related to the aforementioned outbreak.

Draft Genome Sequence of Bacillus megaterium Type Strain ATCC 14581

ABSTRACT Bacillus megaterium is a Gram-positive, rod-shaped, spore-forming bacterium of biotechnological importance. Here, we report a 5.7-Mbp draft genome sequence of B. megaterium ATCC 14581, which is the type strain of the species.

The mechanisms that regulate Vibrio parahaemolyticus virulence gene expression differ between pathotypes

Most Vibrio parahaemolyticus isolates found in marine environments are non-pathogenic; however, certain lineages have acquired genomic pathogenicity islands (PAIs) that enable these isolates to cause human illness. The V. parahaemolyticus PAI contains one or both of two toxins: thermostable direct haemolysin (TDH) or TDH-related haemolysin (TRH) and type III secretion system 2 (T3SS2). Recently, a few V. parahaemolyticus isolates that do not have this PAI were obtained from clinical samples, and there has been interest in determining whether these isolates possess novel virulence factors. In this investigation, we have selected four V. parahaemolyticus isolates: a canonical pathogenic strain containing TDH, TRH and T3SS2; two strains from clinical cases which do not contain a PAI; and an environmental isolate which also does not contain a PAI. For each isolate, we analyzed differential gene expression after crude bile exposure. Several enteric bacterial pathogens are known to use bile as a signal to enhance virulence gene expression. We have shown that in the tdh-positive trh-positive pathotype gene virulence gene expression was not up-regulated in response to crude bile, strongly indicating that the current dogma of virulence gene regulation in V. parahaemolyticus needs to be revisited and separately investigated for each pathotype. In addition, we have created a list of genes of interest that were up-regulated in the non-canonical pathotypes which may contribute to virulence in these isolates.