Tracing Origins of the Salmonella Bareilly Strain Causing a Food-borne Outbreak in the United States.

BACKGROUND: Using a novel combination of whole-genome sequencing (WGS) analysis and geographic metadata, we traced the origins of Salmonella Bareilly isolates collected in 2012 during a widespread food-borne outbreak in the United States associated with scraped tuna imported from India. METHODS: Using next-generation sequencing, we sequenced the complete genome of 100 Salmonella Bareilly isolates obtained from patients who consumed contaminated product, from natural sources, and from unrelated historically and geographically disparate foods. Pathogen genomes were linked to geography by projecting the phylogeny on a virtual globe and produced a transmission network. RESULTS: Phylogenetic analysis of WGS data revealed a common origin for outbreak strains, indicating that patients in Maryland and New York were infected from sources originating at a facility in India. CONCLUSIONS: These data represent the first report fully integrating WGS analysis with geographic mapping and a novel use of transmission networks. Results showed that WGS vastly improves our ability to delimit the scope and source of bacterial food-borne contamination events. Furthermore, these findings reinforce the extraordinary utility that WGS brings to global outbreak investigation as a greatly enhanced approach to protecting the human food supply chain as well as public health in general.

A Nonautochthonous U.S. Strain of Vibrio parahaemolyticus Isolated from Chesapeake Bay Oysters Caused the Outbreak in Maryland in 2010.

UNLABELLED: In the summer of 2010, Vibrio parahaemolyticus caused an outbreak in Maryland linked to the consumption of oysters. Strains isolated from both stool and oyster samples were indistinguishable by pulsed-field gel electrophoresis (PFGE). However, the oysters contained other potentially pathogenic V. parahaemolyticus strains exhibiting different PFGE patterns. In order to assess the identity, genetic makeup, relatedness, and potential pathogenicity of the V. parahaemolyticus strains, we sequenced 11 such strains (2 clinical strains and 9 oyster strains). We analyzed these genomes by in silico multilocus sequence typing (MLST) and determined their phylogeny using a whole-genome MLST (wgMLST) analysis. Our in silico MLST analysis identified six different sequence types (STs) (ST8, ST676, ST810, ST811, ST34, and ST768), with both of the clinical and four of the oyster strains being identified as belonging to ST8. Using wgMLST, we showed that the ST8 strains from clinical and oyster samples were nearly indistinguishable and belonged to the same outbreak, confirming that local oysters were the source of the infections. The remaining oyster strains were genetically diverse, differing in >3,000 loci from the Maryland ST8 strains. eBURST analysis comparing these strains with strains of other STs available at the V. parahaemolyticus MLST website showed that the Maryland ST8 strains belonged to a clonal complex endemic to Asia. This indicates that the ST8 isolates from clinical and oyster sources were likely not endemic to Maryland. Finally, this study demonstrates the utility of whole-genome sequencing (WGS) and associated analyses for source-tracking investigations. IMPORTANCE: Vibrio parahaemolyticus is an important foodborne pathogen and the leading cause of bacterial infections in the United States associated with the consumption of seafood. In the summer of 2010, Vibrio parahaemolyticus caused an outbreak in Maryland linked to oyster consumption. Strains isolated from stool and oyster samples were indistinguishable by pulsed-field gel electrophoresis (PFGE). The oysters also contained other potentially pathogenic V. parahaemolyticus strains with different PFGE patterns. Since their identity, genetic makeup, relatedness, and potential pathogenicity were unknown, their genomes were determined by using next-generation sequencing. Whole-genome sequencing (WGS) analysis by whole-genome multilocus sequence typing (wgMLST) allowed (i) identification of clinical and oyster strains with matching PFGE profiles as belonging to ST8, (ii) determination of oyster strain diversity, and (iii) identification of the clinical strains as belonging to a clonal complex (CC) described only in Asia. Finally, WGS and associated analyses demonstrated their utility for trace-back investigations.

Closed Genome Sequences of Salmonella enterica Serovar Javiana Isolates from Various Sources.

Salmonella enterica serovar Javiana is a major Salmonella serovar that causes human Salmonella infection in the United States. The complete genomic sequences of 9 S. Javiana isolates collected from food, environmental, and kratom sources in the United States were determined by hybrid assembly using Nanopore long-read sequencing and MiSeq short-read sequencing.

Changes in the genomes and methylomes of three Salmonella enterica serovars after long-term storage in ground black pepper.

Low moisture foods (LMFs) have traditionally been recognized as safe for consumption, as most bacteria require higher water content to grow. However, outbreaks due to LMF foods are increasing, and the microbial pathogen Salmonella enterica is frequently implicated. S. enterica can survive in LMFs for years, but few serovars have been studied, and the mechanisms which underlie this longevity are not well understood. Here, we determine that S. enterica serovars S. Tennessee, S. Anatum, and S. Reading but not S. Oranienburg can survive in the ground black pepper for 6 years. S. Reading was not previously associated with any LMF. Using both Illumina and Pacific Biosciences sequencing technologies, we also document changes in the genomes and methylomes of the surviving serovars over this 6-year period. The three serovars acquired a small number of single nucleotide polymorphisms (SNPs) including seven substitutions (four synonymous, two non-synonymous, and one substitution in a non-coding region), and two insertion-deletions. Nine distinct N6-methyladenine (m6A) methylated motifs across the three serovars were identified including five which were previously known, G m6ATC, CAG m6AG, BATGC m6AT, CRT m6AYN6CTC, and CC m6AN7TGAG, and four novel serovar-specific motifs, GRT m6AN8TTYG, GA m6ACN7GTA, GAA m 6A CY, and CAA m6ANCC. Interestingly, the BATGCAT motif was incompletely methylated (35-64% sites across the genome methylated), suggesting a possible role in gene regulation. Furthermore, the number of methylated BATGC m6AT motifs increased after storage in ground black pepper for 6 years from 475 to 657 (S. Tennessee), 366 to 608 (S. Anatum), and 525 to 570 (S. Reading), thus warranting further study as an adaptive mechanism. This is the first long-term assessment of genomic changes in S. enterica in a low moisture environment, and the first study to examine the methylome of any bacteria over a period of years, to our knowledge. These data contribute to our understanding of S. enterica survival in LMFs, and coupled with further studies, will provide the information necessary to design effective interventions which reduce S. enterica in LMFs and maintain a healthy, safe food supply.

Choice of library preparation affects sequence quality, genome assembly, and precise in silico prediction of virulence genes in shiga toxin-producing Escherichia coli.

Whole genome sequencing (WGS) provides essential public health information and is used worldwide for pathogen surveillance, epidemiology, and source tracking. Foodborne pathogens are often sequenced using rapid library preparation chemistries based on transposon technology; however, this method may miss random segments of genomes that can be important for accurate downstream analyses. As new technologies become available, it may become possible to achieve better overall coverage. Here we compare the sequence quality obtained using libraries prepared from the Nextera XT and Nextera DNA Prep (Illumina, San Diego, CA) chemistries for 31 Shiga toxin-producing Escherichia coli (STEC) O121:H19 strains, which had been isolated from flour during a 2016 outbreak. The Nextera DNA Prep gave superior performance metrics including sequence quality, assembly quality, uniformity of genome coverage, and virulence gene identification, among other metrics. Comprehensive detection of virulence genes is essential for making educated assessments of STECs virulence potential. The phylogenetic SNP analysis did not show any differences in the variants detected by either library preparation method which allows isolates prepared from either library method to be analysed together. Our comprehensive comparison of these chemistries should assist researchers wishing to improve their sequencing workflow for STECs and other genomic risk assessments.

Virulence and resistance genes profiles and clonal relationships of non-typhoidal food-borne Salmonella strains isolated in Tunisia by whole genome sequencing.

Whole genome sequencing (WGS) has made impressive progress in the field of molecular biology. Its most common application for public health is in the area of surveillance of food-borne diseases. WGS has the potential for providing a large amount of information, such as the identification of the strain type, the characterization of antibiotic resistance and virulence, and phylogeny. In our study, thirty-nine non-typhoidal Salmonella strains were isolated from diverse sources in Tunisia. Non-typhoidal Salmonella are among the most common pathogens contaminating food animals. The presence of virulence and antimicrobial resistance determinants in those strains were investigated using whole genome sequencing (WGS) and appropriate data analysis. The genomes were screened for several Salmonella virulence genes using the Virulence Factor Database VFDB. Twelve different virulence profiles, which correspond to the 12 identified serovars, were recognized. Several antimicrobial resistance genes were also detected: aac (6')-Iaa, sul1, tetA, bla-TEM and qnrS genes. Phylogenetic relationships among the strains were further assessed by a cgMLST analysis. The resulting phylogenetic tree consisted of several clusters consistently with the in silico multilocus sequence typing (MLST) and serotyping. Our findings demonstrated that WGS and subsequent data analysis provided an accurate tool for genetic characterization of bacterial strains compared to usual molecular typing techniques. To the best of our knowledge, this is the first report of an application of WGS for genetic characterization of food-borne Tunisian strains.

Genomic evidence of environmental and resident Salmonella Senftenberg and Montevideo contamination in the pistachio supply-chain.

Pistachios have been implicated in two salmonellosis outbreaks and multiple recalls in the U.S. This study performed an in-depth retrospective data analysis of Salmonella associated with pistachios as well as a storage study to evaluate the survivability of Salmonella on inoculated inshell pistachios to further understand the genetics and microbiological dynamics of this commodity-pathogen pair. The retrospective data analysis on isolates associated with pistachios was performed utilizing short-read and long-read sequencing technologies. The sequence data were analyzed using two methods: the FDA's Center for Food Safety and Applied Nutrition Single Nucleotide Polymorphism (SNP) analysis and Whole Genome Multilocus Sequence Typing (wgMLST). The year-long storage study evaluated the survival of five strains of Salmonella on pistachios stored at 25 °C at 35% and 54% relative humidity (RH). Our results demonstrate: i) evidence of persistent Salmonella Senftenberg and Salmonella Montevideo strains in pistachio environments, some of which may be due to clonal resident strains and some of which may be due to preharvest contamination; ii) presence of the Copper Homeostasis and Silver Resistance Island (CHASRI) in Salmonella Senftenberg and Montevideo strains in the pistachio supply chain; and iii) the use of metagenomic analysis is a novel tool for determining the composition of serovar survival in a cocktail inoculated storage study.

AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence.

Antimicrobial resistance (AMR) is a significant public health threat. With the rise of affordable whole genome sequencing, in silico approaches to assessing AMR gene content can be used to detect known resistance mechanisms and potentially identify novel mechanisms. To enable accurate assessment of AMR gene content, as part of a multi-agency collaboration, NCBI developed a comprehensive AMR gene database, the Bacterial Antimicrobial Resistance Reference Gene Database and the AMR gene detection tool AMRFinder. Here, we describe the expansion of the Reference Gene Database, now called the Reference Gene Catalog, to include putative acid, biocide, metal, stress resistance genes, in addition to virulence genes and species-specific point mutations. Genes and point mutations are classified by broad functions, as well as more detailed functions. As we have expanded both the functional repertoire of identified genes and functionality, NCBI released a new version of AMRFinder, known as AMRFinderPlus. This new tool allows users the option to utilize only the core set of AMR elements, or include stress response and virulence genes, too. AMRFinderPlus can detect acquired genes and point mutations in both protein and nucleotide sequence. In addition, the evidence used to identify the gene has been expanded to include whether nucleotide or protein sequence was used, its location in the contig, and presence of an internal stop codon. These database improvements and functional expansions will enable increased precision in identifying AMR genes, linking AMR genotypes and phenotypes, and determining possible relationships between AMR, virulence, and stress response.

Five Closed Salmonella enterica Genome Sequences from a 2017-2018 Multistrain, Multistate Kratom Outbreak.

We report here the closed genomes of Salmonella enterica strains from the 2017-2018 multistrain, multistate kratom outbreak using single-molecule real-time DNA sequencing. Four of the genomes consist of one circular chromosome, and the fifth has a circular chromosome and a single plasmid.

Complete Genome Sequences of Four Salmonella enterica Strains Associated with Pistachios Assembled Using a Combination of Short- and Long-Read Sequencing.

Here, we report the genomes of two Salmonella enterica subsp. enterica serovar Montevideo strains (CFSAN005645 and FCC0123) and two Salmonella enterica subsp. enterica serovar Senftenberg strains (FSW0104 and CFSAN087304) isolated from pistachios. The genomes were closed using a hybrid assembly method using short- and long-read sequencing technology.

Erratum for Haendiges et al., "Complete Genome Sequences of Five Salmonella enterica Strains Used in Inoculation Cocktails in Low-Moisture Food Storage Studies".

[This corrects the article DOI: 10.1128/MRA.01588-18.].

Complete Genome Sequences of Five Salmonella enterica Strains Used in Inoculation Cocktails in Low-Moisture Food Storage Studies.

Survival kinetics of Salmonella enterica have been previously studied using an inoculum cocktail composed of different strains that have been associated with low-moisture foods. Here, we report the closed genome sequences of five strains of Salmonella enterica that are commonly used in these storage studies.

Whole-Genome Analysis of Halomonas sp. Soap Lake #7 Reveals It Possesses Putative Mrp Antiporter Operon Groups 1 and 2.

The genus Halomonas possesses bacteria that are halophilic or halotolerant and exhibit a wide range of pH tolerance. The genome of Halomonas sp. Soap Lake #7 was sequenced to provide a better understanding of the mechanisms for salt and pH tolerance in this genus. The bacterium's genome was found to possess two complete multiple resistance and pH antiporter systems, Group 1 and Group 2. This is the first report of both multiple resistance and pH antiporter Groups 1 and 2 in the genome of a haloalkaliphilic bacterium.

Salmonella enterica serovar Kentucky recovered from human clinical cases in Maryland, USA (2011-2015).

Salmonella Kentucky is among the most frequently isolated S. enterica serovars from food animals in the United States. Recent research on isolates recovered from these animals suggests there may be geographic and host specificity signatures associated with S. Kentucky strains. However, the sources and genomic features of human clinical S. Kentucky isolated in the United States remain poorly described. To investigate the characteristics of clinical S. Kentucky and the possible sources of these infections, the genomes of all S. Kentucky isolates recovered from human clinical cases in the State of Maryland between 2011 and 2015 (n = 12) were sequenced and compared to a database of 525 previously sequenced S. Kentucky genomes representing 12 sequence types (ST) collected from multiple sources on several continents. Of the 12 human clinical S. Kentucky isolates from Maryland, nine were ST198, two were ST152, and one was ST314. Forty-one per cent of isolates were recovered from patients reporting recent international travel and 58% of isolates encoded genomic characteristics similar to those originating outside of the United States. Of the five isolates not associated with international travel, three encoded antibiotic resistance genes conferring resistance to tetracycline or aminoglycosides, while two others only encoded the cryptic aac(6')-Iaa gene. Five isolates recovered from individuals with international travel histories (ST198) and two for which travel was not recorded (ST198) encoded genes conferring resistance to between 4 and 7 classes of antibiotics. Seven ST198 genomes encoded the Salmonella Genomic Island 1 and substitutions in the gyrA and parC genes known to confer resistance to ciprofloxacin. Case report data on food consumption and travel were, for the most part, consistent with the inferred S. Kentucky phylogeny. Results of this study indicate that the majority of S. Kentucky infections in Maryland are caused by ST198 which may originate outside of North America.

First Complete Genome Sequence of Salmonella enterica subsp. enterica Serovar Worthington Strain CFSAN051295, Isolated from Pistachio.

We report here, using third-generation, single-molecule, real-time DNA sequencing, the first complete genome sequence of Salmonella enterica serovar Worthington CFSAN051295, isolated from pistachios in the United States. The genome consists of a single 4.9-Mb chromosome.

Closed Genome Sequences of Two Clostridium botulinum Strains Obtained by Nanopore Sequencing.

Here we report the genome sequences of two toxin-producing Clostridium botulinum strains, one environmental sample (83F) and one clinical sample (CDC51232). The genomes were closed by a combination of long-read and short-read sequencing. The strains belong to C. botulinum sequence type 4 (ST4) and ST7, respectively.

Complete Genome Sequences of Four Salmonella enterica subsp. enterica Serovar Senftenberg and Montevideo Isolates Associated with a 2016 Multistate Outbreak in the United States.

A multistate outbreak of 11 Salmonella infections linked to pistachio nuts occurred in 2016. In this announcement, we report the complete genome sequences of four Salmonella enterica subsp. enterica serovar Senftenberg and S. enterica subsp. enterica serovar Montevideo isolates from pistachios collected during the 2016 outbreak investigation.

Genetic diversity and virulence profiles of Listeria monocytogenes recovered from bulk tank milk, milk filters, and milking equipment from dairies in the United States (2002 to 2014).

Unpasteurized dairy products are known to occasionally harbor Listeria monocytogenes and have been implicated in recent listeriosis outbreaks and numerous sporadic cases of listeriosis. However, the diversity and virulence profiles of L. monocytogenes isolates recovered from these products have not been fully described. Here we report a genomic analysis of 121 L. monocytogenes isolates recovered from milk, milk filters, and milking equipment collected from bovine dairy farms in 19 states over a 12-year period. In a multi-virulence-locus sequence typing (MVLST) analysis, 59 Virulence Types (VT) were identified, of which 25% were Epidemic Clones I, II, V, VI, VII, VIII, IX, or X, and 31 were novel VT. In a multi-locus sequence typing (MLST) analysis, 60 Sequence Types (ST) of 56 Clonal Complexes (CC) were identified. Within lineage I, CC5 and CC1 were among the most abundant, and within lineage II, CC7 and CC37 were the most abundant. Multiple CCs previously associated with central nervous system and maternal-neonatal infections were identified. A genomic analysis identified variable distribution of virulence markers, Listeria pathogenicity islands (LIPI) -1, -3, and -4, and stress survival island-1 (SSI-1). Of these, 14 virulence markers, including LIPI-3 and -4 were more frequently detected in one lineage (I or II) than the other. LIPI-3 and LIPI-4 were identified in 68% and 28% of lineage I CCs, respectively. Results of this analysis indicate that there is a high level of genetic diversity among the L. monocytogenes present in bulk tank milk in the United States with some strains being more frequently detected than others, and some being similar to those that have been isolated from previous non-dairy related outbreaks. Results of this study also demonstrate significant number of strains isolated from dairy farms encode virulence markers associated with severe human disease.