quasi-metagenomic-analysis-salmonella-from-food-environmental

Quasi-metagenomic Analysis of <em>Salmonella</em> from Food and Environmental Samples

Quasi-metagenomics sequencing refers to the sequencing-based analysis of modified microbiomes of food and environmental samples. In this protocol, ...

Center for Food Safety

Salmonella enterica serotype Enteritidis is one of the most commonly reported causes of human salmonellosis. Its low genetic diversity, measured by fingerprinting methods, has made subtyping a challenge. We used whole-genome sequencing to characterize 125 S. enterica Enteritidis and 3 S. enterica serotype Nitra strains. Single-nucleotide polymorphisms were filtered to identify 4,887 reliable loci that distinguished all isolates from each other. Our whole-genome single-nucleotide polymorphism typing approach was robust for S. enterica Enteritidis subtyping with combined data for different strains from 2 different sequencing platforms. Five major genetic lineages were recognized, which revealed possible patterns of geographic and epidemiologic distribution. Analyses on the population dynamics and evolutionary history estimated that major lineages emerged during the 17th-18th centuries and diversified during the 1920s and 1950s. 

Genomic epidemiology of Salmonella enterica serotype Enteritidis based on population structure of prevalent lineages.

A retrospective investigation was performed to evaluate whole-genome sequencing as a benchmark for comparing molecular subtyping methods for Salmonella enterica serotype Enteritidis and survey the population structure of commonly encountered S. enterica serotype Enteritidis outbreak isolates in the United States. A total of 52 S. enterica serotype Enteritidis isolates representing 16 major outbreaks and three sporadic cases collected between 2001 and 2012 were sequenced and subjected to subtyping by four different methods: (i) whole-genome single-nucleotide-polymorphism typing (WGST), (ii) multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA), (iii) clustered regularly interspaced short palindromic repeats combined with multi-virulence-locus sequence typing (CRISPR-MVLST), and (iv) pulsed-field gel electrophoresis (PFGE). WGST resolved all outbreak clusters and provided useful robust phylogenetic inference results with high epidemiological correlation. While both MLVA and CRISPR-MVLST yielded higher discriminatory power than PFGE, MLVA outperformed the other methods in delineating outbreak clusters whereas CRISPR-MVLST showed the potential to trace major lineages and ecological origins of S. enterica serotype Enteritidis. Our results suggested that whole-genome sequencing makes a viable platform for the evaluation and benchmarking of molecular subtyping methods. 

Comparative analysis of subtyping methods against a whole-genome-sequencing standard for Salmonella enterica serotype Enteritidis.

Serotyping forms the basis of national and international surveillance networks for Salmonella, one of the most prevalent foodborne pathogens worldwide (1-3). Public health microbiology is currently being transformed by whole-genome sequencing (WGS), which opens the door to serotype determination using WGS data. SeqSero (www.denglab.info/SeqSero) is a novel Web-based tool for determining Salmonella serotypes using high-throughput genome sequencing data. SeqSero is based on curated databases of Salmonella serotype determinants (rfb gene cluster, fliC and fljB alleles) and is predicted to determine serotype rapidly and accurately for nearly the full spectrum of Salmonella serotypes (more than 2,300 serotypes), from both raw sequencing reads and genome assemblies. The performance of SeqSero was evaluated by testing (i) raw reads from genomes of 308 Salmonella isolates of known serotype; (ii) raw reads from genomes of 3,306 Salmonella isolates sequenced and made publicly available by GenomeTrakr, a U.S. national monitoring network operated by the Food and Drug Administration; and (iii) 354 other publicly available draft or complete Salmonella genomes. We also demonstrated Salmonella serotype determination from raw sequencing reads of fecal metagenomes from mice orally infected with this pathogen. SeqSero can help to maintain the well-established utility of Salmonella serotyping when integrated into a platform of WGS-based pathogen subtyping and characterization. 

Salmonella serotype determination utilizing high-throughput genome sequencing data.

During the past century, the microbiological safety of the US food supply has improved; however, many foodborne illnesses and outbreaks occur annually. Hence, opportunities for the food industry to improve the safety of both domestic and imported food exist through the adoption of risk-based preventive measures. Challenging food safety issues that are on the horizon include demographic changes to a population whose immune system is more susceptible to foodborne and opportunistic pathogens, climate changes that will shift where food is produced, and consumers' preferences for raw and minimally processed foods. Increased environmental and product testing and anonymous data sharing by the food industry with the public health community would aid in identifying system weaknesses and enabling more targeted corrective and preventive actions. Clinicians will continue to play a major role in reducing foodborne illnesses by diagnosing and reporting cases and in helping to educate the consumer about food safety practices. 

The food industry's current and future role in preventing microbial foodborne illness within the United States.

As we are approaching the twentieth anniversary of PulseNet, a network of public health and regulatory laboratories that has changed the landscape of foodborne illness surveillance through molecular subtyping, public health microbiology is undergoing another transformation brought about by so-called next-generation sequencing (NGS) technologies that have made whole-genome sequencing (WGS) of foodborne bacterial pathogens a realistic and superior alternative to traditional subtyping methods. Routine, real-time, and widespread application of WGS in food safety and public health is on the horizon. Technological, operational, and policy challenges are still present and being addressed by an international and multidisciplinary community of researchers, public health practitioners, and other stakeholders. 

Genomic Epidemiology: Whole-Genome-Sequencing-Powered Surveillance and Outbreak Investigation of Foodborne Bacterial Pathogens.

An epidemiological paradox surrounds Salmonella enterica serovar Enteritidis. In high-income settings, it has been responsible for an epidemic of poultry-associated, self-limiting enterocolitis, whereas in sub-Saharan Africa it is a major cause of invasive nontyphoidal Salmonella disease, associated with high case fatality. By whole-genome sequence analysis of 675 isolates of S. Enteritidis from 45 countries, we show the existence of a global epidemic clade and two new clades of S. Enteritidis that are geographically restricted to distinct regions of Africa. The African isolates display genomic degradation, a novel prophage repertoire, and an expanded multidrug resistance plasmid. S. Enteritidis is a further example of a Salmonella serotype that displays niche plasticity, with distinct clades that enable it to become a prominent cause of gastroenteritis in association with the industrial production of eggs and of multidrug-resistant, bloodstream-invasive infection in Africa.

Distinct Salmonella Enteritidis lineages associated with enterocolitis in high-income settings and invasive disease in low-income settings.

We presented the first attempt to combine immunomagnetic separation (IMS), whole genome amplification by multiple displacement amplification (MDA) and real-time PCR for detecting a bacterial pathogen in a food sample. This method was effective in enabling real-time PCR detection of low levels of Salmonella enterica Serotype Enteritidis (SE) (∼10 CFU/g) in raw chicken breast without culture enrichment. In addition, it was able to detect refrigeration-stressed SE cells at lower concentrations (∼0.1 CFU/g) in raw chicken breast after a 4-h culture enrichment, shortening the detection process from days to hours and displaying no statistical difference in detection rate in comparison with a culture-based detection method. By substantially improving performance in SE detection over conventional real-time PCR, we demonstrated the potential of IMS-MDA real-time PCR as a rapid, sensitive and affordable method for detecting Salmonella in food.

Rapid detection of Salmonella in raw chicken breast using real-time PCR combined with immunomagnetic separation and whole genome amplification.

Modern epidemiology of foodborne bacterial pathogens in industrialized countries relies increasingly on whole genome sequencing (WGS) techniques. As opposed to profiling techniques such as pulsed-field gel electrophoresis, WGS requires a variety of computational methods. Since 2013, United States agencies responsible for food safety including the CDC, FDA, and USDA, have been performing whole-genome sequencing (WGS) on all  found in clinical, food, and environmental samples. Each year, more genomes of other foodborne pathogens such as , and  are being sequenced. Comparing thousands of genomes across an entire species requires a fast method with coarse resolution; however, capturing the fine details of highly related isolates requires a computationally heavy and sophisticated algorithm. Most  investigations employing WGS depend on being able to identify an outbreak clade whose inter-genomic distances are less than an empirically determined threshold. When the difference between a few single nucleotide polymorphisms (SNPs) can help distinguish between genomes that are likely outbreak-associated and those that are less likely to be associated, we require a fine-resolution method. To achieve this level of resolution, we have developed Lyve-SET, a high-quality SNP pipeline. We evaluated Lyve-SET by retrospectively investigating 12 outbreak data sets along with four other SNP pipelines that have been used in outbreak investigation or similar scenarios. To compare these pipelines, several distance and phylogeny-based comparison methods were applied, which collectively showed that multiple pipelines were able to identify most outbreak clusters and strains. Currently in the US PulseNet system, whole genome multi-locus sequence typing (wgMLST) is the preferred primary method for foodborne WGS cluster detection and outbreak investigation due to its ability to name standardized genomic profiles, its central database, and its ability to be run in a graphical user interface. However, creating a functional wgMLST scheme requires extended up-front development and subject-matter expertise. When a scheme does not exist or when the highest resolution is needed, SNP analysis is used. Using three  outbreak data sets, we demonstrated the concordance between Lyve-SET SNP typing and wgMLST. : Lyve-SET can be found at https://github.com/lskatz/Lyve-SET.

A Comparative Analysis of the Lyve-SET Phylogenomics Pipeline for Genomic Epidemiology of Foodborne Pathogens.

Salmonella survival was characterized and modeled during thermal dehydration of fresh garlic and storage of dehydrated garlic products. In our experiments that simulated commercial dehydration processing at 80±5°C, moderate level of Salmonella contamination (4-5logCFU/g) on fresh garlic was reduced below the enumeration limit (1.7logCFU/g) after 4.5h of dehydration and not detectable by culture enrichment after 7h. With high level of contamination (7-8logCFU/g), the Salmonella population persisted at 3.6logCFU/g after 8h of processing. By increasing the dehydration temperature to 90±5°C, the moderate and high levels of initial Salmonella load on fresh garlic dropped below the enumeration limit after 1.5 and 3.75h of processing and became undetectable by culture enrichment after 2.5 and 6h, respectively. During the storage of dried garlic products, Salmonella was not able to grow under all tested combinations of temperature (25 and 35°C) and water activity (0.56-0.98) levels, suggesting active inhibition. Storage temperature played a primary role in determining Salmonella survival on dehydrated garlic flakes. Under a typical storage condition at 25°C and ambient relative humidity, Salmonella could persist over months with the population gradually declining (4.3 log reduction over 88days). Granular size of dehydrated garlic had an impact on Salmonella survival, with better survival of the pathogen observed in bigger granules. At the early stage of dehydrated garlic storage (until 7days), rising water activity appeared to initially promote but then inhibited Salmonella survival, resulting in a water activity threshold at 0.73 where Salmonella displayed strongest persistence. However, this phenomenon was less apparent during extended storage (after 14days).

Salmonella survival during thermal dehydration of fresh garlic and storage of dehydrated garlic products.

Metagenomics analysis of food samples promises isolation-independent detection and subtyping of foodborne bacterial pathogens in a single workflow. Selective concentration of  genomic DNA through immunomagnetic separation (IMS) and multiple displacement amplification (MDA) were shown to shorten culture enrichment of -spiked raw chicken breast samples by over 12 hours while permitting serotyping and high-fidelity single nucleotide polymorphisms (SNP) typing of the pathogen using short shotgun sequencing reads. The herein termed quasi-metagenomics approach was evaluated on -spiked lettuce and black peppercorn samples as well as retail chicken parts naturally contaminated with different serotypes of  Between 8 and 24 h culture enrichment was required for detecting and subtyping naturally occurring  from unspiked chicken parts compared with 4 to 12 h culture enrichment when -spiked food samples were analyzed, indicating the likely need for longer culture enrichment to revive low levels of stressed or injured  cells in food. Further acceleration of the workflow was achieved by real-time nanopore sequencing. After 1.5 hours of analysis on a potable sequencer, sufficient data were generated from sequencing IMS-MDA product of a cultured-enriched lettuce sample to allow serotyping and robust phylogenetic placement of the inoculated isolate. Both culture enrichment and next-generation sequencing remain to be time-consuming processes for food testing where rapid methods for pathogen detection are widely available. Our study demonstrated substantial acceleration of the respective process through IMS-MDA and real-time nanopore sequencing. In one example, the combined use of the two methods delivered a less than 24 h turnaround time from a -contaminated lettuce sample to phylogenetic identification of the pathogen. Improved efficiency like this is important for further expanding the use of whole genome and metagenomics sequencing in microbial analysis of food. Our results suggest the potential of the quasi-metagenomics approach in areas where rapid detection and subtyping of foodborne pathogens is important, such as foodborne outbreak response and precision tracking and monitoring of foodborne pathogens in production environments and supply chains.

Xiangyu Deng

Center for Food Safety

Quasi-metagenomic Analysis of Salmonella from Food and Environmental Samples