Whole Genome Sequencing Applied to Pathogen Source Tracking in Food Industry: Key Considerations for Robust Bioinformatics Data Analysis and Reliable Results Interpretation

CamPype: an open-source workflow for automated bacterial whole-genome sequencing analysis focused on Campylobacter

BACKGROUND

The rapid expansion of Whole-Genome Sequencing has revolutionized the fields of clinical and food microbiology. However, its implementation as a routine laboratory technique remains challenging due to the growth of data at a faster rate than can be effectively analyzed and critical gaps in bioinformatics knowledge.

RESULTS

To address both issues, CamPype was developed as a new bioinformatics workflow for the genomics analysis of sequencing data of bacteria, especially Campylobacter, which is the main cause of gastroenteritis worldwide making a negative impact on the economy of the public health systems. CamPype allows fully customization of stages to run and tools to use, including read quality control filtering, read contamination, reads extension and assembly, bacterial typing, genome annotation, searching for antibiotic resistance genes, virulence genes and plasmids, pangenome construction and identification of nucleotide variants. All results are processed and resumed in an interactive HTML report for best data visualization and interpretation.

CONCLUSIONS

The minimal user intervention of CamPype makes of this workflow an attractive resource for microbiology laboratories with no expertise in bioinformatics as a first line method for bacterial typing and epidemiological analyses, that would help to reduce the costs of disease outbreaks, or for comparative genomic analyses. CamPype is publicly available at https://github.com/JoseBarbero/CamPype .

Serotype Screening of Salmonella enterica Subspecies I by Intergenic Sequence Ribotyping (ISR): Critical Updates

(1) Background: Foodborne illness from Salmonella enterica subspecies I is most associated with approximately 32 out of 1600 serotypes. While whole genome sequencing and other nucleic acid-based methods are preferred for serotyping, they require expertise in bioinformatics and often submission to an external agency. Intergenic Sequence Ribotyping (ISR) assigns serotype to Salmonella in coordination with information freely available at the National Center for Biotechnology Information. ISR requires updating because it was developed from 26 genomes while there are now currently 1804 genomes and 1685 plasmids. (2) Methods: Serotypes available for sequencing were analyzed by ISR to confirm primer efficacy and to identify any issues in application. Differences between the 2012 and 2022 ISR database were tabulated, nomenclature edited, and instances of multiple serotypes aligning to a single ISR were examined. (3) Results: The 2022 ISR database has 268 sequences and 40 of these were assigned new NCBI accession numbers that were not previously available. Extending boundaries of sequences resolved hdfR cross-alignment and reduced multiplicity of alignment for 37 ISRs. Comparison of gene cyaA sequences and some cell surface epitopes provided evidence that homologous recombination was potentially impacting results for this subset. There were 99 sequences that still had no match with an NCBI submission. (4) The 2022 ISR database is available for use as a serotype screening method for Salmonella enterica subspecies I. Finding that 36.9% of the sequences in the ISR database still have no match within the NCBI Salmonella enterica database suggests that there is more genomic heterogeneity yet to characterize.

Rapid Proteomic Characterization of Bacteriocin-Producing Enterococcus faecium Strains from Foodstuffs

Enterococcus belongs to a group of microorganisms known as lactic acid bacteria (LAB), which constitute a broad heterogeneous group of generally food-grade microorganisms historically used in food preservation. Enterococci live as commensals of the gastrointestinal tract of warm-blooded animals, although they also are present in food of animal origin (milk, cheese, fermented sausages), vegetables, and plant materials because of their ability to survive heat treatments and adverse environmental conditions. The biotechnological traits of enterococci can be applied in the food industry; however, the emergence of enterococci as a cause of nosocomial infections makes their food status uncertain. Recent advances in high-throughput sequencing allow the subtyping of bacterial pathogens, but it cannot reflect the temporal dynamics and functional activities of microbiomes or bacterial isolates. Moreover, genetic analysis is based on sequence homologies, inferring functions from databases. Here, we used an end-to-end proteomic workflow to rapidly characterize two bacteriocin-producing Enterococcus faecium (Efm) strains. The proteome analysis was performed with liquid chromatography coupled to a trapped ion mobility spectrometry-time-of-flight mass spectrometry instrument (TimsTOF) for high-throughput and high-resolution characterization of bacterial proteins. Thus, we identified almost half of the proteins predicted in the bacterial genomes (>1100 unique proteins per isolate), including quantifying proteins conferring resistance to antibiotics, heavy metals, virulence factors, and bacteriocins. The obtained proteomes were annotated according to function, resulting in 22 complete KEGG metabolic pathway modules for both strains. The workflow used here successfully characterized these bacterial isolates and showed great promise for determining and optimizing the bioengineering and biotechnology properties of other LAB strains in the food industry.

The power, potential, benefits, and challenges of implementing high-throughput sequencing in food safety systems

The development and application of modern sequencing technologies have led to many new improvements in food safety and public health. With unprecedented resolution and big data, high-throughput sequencing (HTS) has enabled food safety specialists to sequence marker genes, whole genomes, and transcriptomes of microorganisms almost in real-time. These data reveal not only the identity of a pathogen or an organism of interest in the food supply but its virulence potential and functional characteristics. HTS of amplicons, allow better characterization of the microbial communities associated with food and the environment. New and powerful bioinformatics tools, algorithms, and machine learning allow for development of new models to predict and tackle important events such as foodborne disease outbreaks. Despite its potential, the integration of HTS into current food safety systems is far from complete. Government agencies have embraced this new technology, and use it for disease diagnostics, food safety inspections, and outbreak investigations. However, adoption and application of HTS by the food industry have been comparatively slow, sporadic, and fragmented. Incorporation of HTS by food manufacturers in their food safety programs could reinforce the design and verification of effectiveness of control measures by providing greater insight into the characteristics, origin, relatedness, and evolution of microorganisms in our foods and environment. Here, we discuss this new technology, its power, and potential. A brief history of implementation by public health agencies is presented, as are the benefits and challenges for the food industry, and its future in the context of food safety.

Phylogenomic Analysis of Salmonella enterica subsp. enterica Serovar Bovismorbificans from Clinical and Food Samples Using Whole Genome Wide Core Genes and kmer Binning Methods to Identify Two Distinct Polyphyletic Genome Pathotypes

Salmonella enterica subsp. enterica serovar Bovismorbificans has caused multiple outbreaks involving the consumption of produce, hummus, and processed meat products worldwide. To elucidate the intra-serovar genomic structure of S. Bovismorbificans, a core-genome analysis with 2690 loci (based on 150 complete genomes representing Salmonella enterica serovars developed as part of this study) and a k-mer-binning based strategy were carried out on 95 whole genome sequencing (WGS) assemblies from Swiss, Canadian, and USA collections of S. Bovismorbificans strains from foodborne infections. Data mining of a digital DNA tiling array of legacy SARA and SARB strains was conducted to identify near-neighbors of S. Bovismorbificans. The core genome analysis and the k-mer-binning methods identified two polyphyletic clusters, each with emerging evolutionary properties. Four STs (2640, 142, 1499, and 377), which constituted the majority of the publicly available WGS datasets from >260 strains analyzed by k-mer-binning based strategy, contained a conserved core genome backbone with a different evolutionary lineage as compared to strains comprising the other cluster (ST150). In addition, the assortment of genotypic features contributing to pathogenesis and persistence, such as antimicrobial resistance, prophage, plasmid, and virulence factor genes, were assessed to understand the emerging characteristics of this serovar that are relevant clinically and for food safety concerns. The phylogenomic profiling of polyphyletic S. Bovismorbificans in this study corresponds to intra-serovar variations observed in S. Napoli and S. Newport serovars using similar high-resolution genomic profiling approaches and contributes to the understanding of the evolution and sequence divergence of foodborne Salmonellae. These intra-serovar differences may have to be thoroughly understood for the accurate classification of foodborne Salmonella strains needed for the uniform development of future food safety mitigation strategies.

Surveillance of Listeria monocytogenes: Early Detection, Population Dynamics, and Quasimetagenomic Sequencing during Selective Enrichment

In this study, we addressed different aspects regarding the implementation of quasimetagenomic sequencing as a hybrid surveillance method in combination with enrichment for early detection of Listeria monocytogenes in the food industry. Different experimental enrichment cultures were used, comprising seven L. monocytogenes strains of different sequence types (STs), with and without a background microbiota community. To assess whether the proportions of the different STs changed over time during enrichment, the growth and population dynamics were assessed using dapE colony sequencing and dapE and 16S rRNA amplicon sequencing. There was a tendency of some STs to have a higher relative abundance during the late stage of enrichment when L. monocytogenes was enriched without background microbiota. When coenriched with background microbiota, the population dynamics of the different STs was more consistent over time. To evaluate the earliest possible time point during enrichment that allows the detection of L. monocytogenes and at the same time the generation of genetic information that enables an estimation regarding the strain diversity in a sample, quasimetagenomic sequencing was performed early during enrichment in the presence of the background microbiota using Oxford Nanopore Technologies Flongle and Illumina MiSeq sequencing. The application of multiple displacement amplification (MDA) enabled detection of L. monocytogenes (and the background microbiota) after only 4 h of enrichment using both applied sequencing approaches. The MiSeq sequencing data additionally enabled the prediction of cooccurring L. monocytogenes strains in the samples.

IMPORTANCE We showed that a combination of a short primary enrichment combined with MDA and Nanopore sequencing can accelerate the traditional process of cultivation and identification of L. monocytogenes. The use of Illumina MiSeq sequencing additionally allowed us to predict the presence of cooccurring L. monocytogenes strains. Our results suggest quasimetagenomic sequencing is a valuable and promising hybrid surveillance tool for the food industry that enables faster identification of L. monocytogenes during early enrichment. Routine application of this approach could lead to more efficient and proactive actions in the food industry that prevent contamination and subsequent product recalls and food destruction, economic and reputational losses, and human listeriosis cases.

Interpretative labor and the bane of non-standardized metadata in public health surveillance and food safety

Open source DNA sequence databases have long been touted as beneficial to public health, including the facilitation of earlier detection and response to infectious disease outbreaks. Of critical importance to harnessing these benefits is the metadata which describes general and other domain specific attributes (e.g., collection location, isolate type, etc.) of a sample. Unlike the sequence data, the metadata is often incomplete and lacks adherence to an international standard. We describe the problem posed by such variable and incomplete metadata in terms of interpretative labor costs (the time and energy necessary to make sense of the signal in the genetic data), and the impact such metadata has on foodborne outbreak detection and response. Improving the quality of sequence-associated metadata would allow for earlier detection of emerging food safety hazards and allow faster response to foodborne outbreaks.