Strain-level metagenomic data analysis of enriched in vitro and in silico spiked food samples : paving the way towards a culture-free foodborne outbreak investigation using STEC as a case study

Culture-independent diagnostics, such as metagenomic shotgun sequencing of food samples, could not only reduce the turnaround time of samples in an outbreak investigation, but also allow the detection of multi-species and multi-strain outbreaks. For successful foodborne outbreak investigation using a metagenomic approach, it is, however, necessary to bioinformatically separate the genomes of individual strains, including strains belonging to the same species, present in a microbial community, which has up until now not been demonstrated for this application. The current work shows the feasibility of strain-level metagenomics of enriched food matrix samples making use of data analysis tools that classify reads against a sequence database. It includes a brief comparison of two database-based read classification tools, Sigma and Sparse, using a mock community obtained by in vitro spiking minced meat with a Shiga toxin-producing Escherichia coli (STEC) isolate originating from a described outbreak. The more optimal tool Sigma was further evaluated using in silico simulated metagenomic data to explore the possibilities and limitations of this data analysis approach. The performed analysis allowed us to link the pathogenic strains from food samples to human isolates previously collected during the same outbreak, demonstrating that the metagenomic approach could be applied for the rapid source tracking of foodborne outbreaks. To our knowledge, this is the first study demonstrating a data analysis approach for detailed characterization and phylogenetic placement of multiple bacterial strains of one species from shotgun metagenomic WGS data of an enriched food sample

A practical method to implement strain-level metagenomics-based foodborne outbreak investigation and source tracking in routine

The management of a foodborne outbreak depends on the rapid and accurate identification of the responsible food source. Conventional methods based on isolation of the pathogen from the food matrix and target-specific real-time polymerase chain reactions (qPCRs) are used in routine. In recent years, the use of whole genome sequencing (WGS) of bacterial isolates has proven its value to collect relevant information for strain characterization as well as tracing the origin of the contamination by linking the food isolate with the patient’s isolate with high resolution. However, the isolation of a bacterial pathogen from food matrices is often time-consuming and not always successful. Therefore, we aimed to improve outbreak investigation by developing a method that can be implemented in reference laboratories to characterize the pathogen in the food vehicle without its prior isolation and link it back to human cases. We tested and validated a shotgun metagenomics approach by spiking food pathogens in specific food matrices using the Shiga toxin-producing Escherichia coli (STEC) as a case study. Different DNA extraction kits and enrichment procedures were investigated to obtain the most practical workflow. We demonstrated the feasibility of shotgun metagenomics to obtain the same information as in ISO/TS 13136:2012 and WGS of the isolate in parallel by inferring the genome of the contaminant and characterizing it in a shorter timeframe. This was achieved in food samples containing different E. coli strains, including a combination of different STEC strains. For the first time, we also managed to link individual strains from a food product to isolates from human cases, demonstrating the power of shotgun metagenomics for rapid outbreak investigation and source tracking

Perceived utility and feasibility of pathogen genomics for public health practice : a survey among public health professionals working in the field of infectious diseases, Belgium, 2019

Background Pathogen genomics is increasingly being translated from the research setting into the activities of public health professionals operating at different levels. This survey aims to appraise the literacy level and gather the opinions of public health experts and allied professionals working in the field of infectious diseases in Belgium concerning the implementation of next-generation sequencing (NGS) in public health practice. Methods In May 2019, Belgian public health and healthcare professionals were invited to complete an online survey containing eight main topics including background questions, general attitude towards pathogen genomics for public health practice and main concerns, genomic literacy, current and planned NGS activities, place of NGS in diagnostic microbiology pathways, data sharing obstacles, end-user requirements, and key drivers for the implementation of NGS. Descriptive statistics were used to report on the frequency distribution of multiple choice responses whereas thematic analysis was used to analyze free text responses. A multivariable logistic regression model was constructed to identify important predictors for a positive attitude towards the implementation of pathogen genomics in public health practice. Results 146 out of the 753 invited public health professionals completed the survey. 63% of respondents indicated that public health agencies should be using genomics to understand and control infectious diseases. Having a high level of expertise in the field of pathogen genomics was the strongest predictor of a positive attitude (OR = 4.04, 95% CI = 1.11 – 17.23). A significantly higher proportion of data providers indicated to have followed training in the field of pathogen genomics compared to data end-users (p < 0.001). Overall, 79% of participants expressed interest in receiving further training. Main concerns were related to the cost of sequencing technologies, data sharing, data integration, interdisciplinary working, and bioinformatics expertise. Conclusions Belgian health professionals expressed favorable views about implementation of pathogen genomics in their work activities related to infectious disease surveillance and control. They expressed the need for suitable training initiatives to strengthen their competences in the field. Their perception of the utility and feasibility of pathogen genomics for public health purposes will be a key driver for its further implementation

Combining short and long read sequencing to characterize antimicrobial resistance genes on plasmids applied to an unauthorized genetically modified Bacillus

Antimicrobial resistance (AMR) is a major public health threat. Plasmids are able to transfer AMR genes among bacterial isolates. Whole genome sequencing (WGS) is a powerful tool to monitor AMR determinants. However, plasmids are difficult to reconstruct from WGS data. This study aimed to improve the characterization, including the localization of AMR genes using short and long read WGS strategies. We used a genetically modified (GM) Bacillus subtilis isolated as unexpected contamination in a feed additive, and therefore considered unauthorized (RASFF 2014.1249), as a case study. In GM organisms, AMR genes are used as selection markers. Because of the concern of spread of these AMR genes when present on mobile genetic elements, it is crucial to characterize their location. Our approach resulted in an assembly of one chromosome and one plasmid, each with several AMR determinants of which five are against critically important antibiotics. Interestingly, we found several plasmids, containing AMR genes, integrated in the chromosome in a repetitive region of at least 53 kb. Our findings would have been impossible using short reads only. We illustrated the added value of long read sequencing in addressing the challenges of plasmid reconstruction within the context of evaluating the risk of AMR spread

DISTILLER: a data integration framework to reveal condition dependency of complex regulons in Escherichia coli

DISTILLER, a data integration framework for the inference of transcriptional module networks, is presented and used to investigate the condition dependency and modularity in Escherichia coli networks

Application of a strain- level shotgun metagenomics approach on food samples : resolution of the source of a Salmonella food-borne outbreak

Food- borne outbreak investigation currently relies on the time- consuming and challenging bacterial isolation from food, to be able to link food- derived strains to more easily obtained isolates from infected people. When no food isolate can be obtained, the source of the outbreak cannot be unambiguously determined. Shotgun metagenomics approaches applied to the food samples could circumvent this need for isolation from the suspected source, but require downstream strain- level data analysis to be able to accurately link to the human isolate. Until now, this approach has not yet been applied outside research settings to analyse real food- borne outbreak samples. In September 2019, a Salmonella outbreak occurred in a hotel school in Bruges, Belgium, affecting over 200 students and teachers. Following standard procedures, the Belgian National Reference Center for human salmonellosis and the National Reference Laboratory for Salmonella in food and feed used conventional analysis based on isolation, serotyping and MLVA (multilocus variable number tandem repeat analysis) comparison, followed by wholegenome sequencing, to confirm the source of the contamination over 2 weeks after receipt of the sample, which was freshly prepared tartar sauce in a meal cooked at the school. Our team used this outbreak as a case study to deliver a proof of concept for a short- read strain- level shotgun metagenomics approach for source tracking. We received two suspect food samples: the full meal and some freshly made tartar sauce served with this meal, requiring the use of raw eggs. After analysis, we could prove, without isolation, that Salmonella was present in both samples, and we obtained an inferred genome of a Salmonella enterica subsp. enterica serovar Enteritidis that could be linked back to the human isolates of the outbreak in a phylogenetic tree. These metagenomics- derived outbreak strains were separated from sporadic cases as well as from another outbreak circulating in Europe at the same time period. This is, to our knowledge, the first Salmonella food- borne outbreak investigation uniquely linking the food source using a metagenomics approach and this in a fast time frame

Validation strategy of a bioinformatics whole genome sequencing workflow for Shiga toxin-producing Escherichia coli using a reference collection extensively characterized with conventional methods

Whole genome sequencing (WGS) enables complete characterization of bacterial pathogenic isolates at single nucleotide resolution, making it the ultimate tool for routine surveillance and outbreak investigation. The lack of standardization, and the variation regarding bioinformatics workflows and parameters, however, complicates interoperability among (inter)national laboratories. We present a validation strategy applied to a bioinformatics workflow for Illumina data that performs complete characterization of Shiga toxin-producing Escherichia coli (STEC) isolates including antimicrobial resistance prediction, virulence gene detection, serotype prediction, plasmid replicon detection and sequence typing. The workflow supports three commonly used bioinformatics approaches for the detection of genes and alleles: alignment with blast+, kmer-based read mapping with KMA, and direct read mapping with SRST2. A collection of 131 STEC isolates collected from food and human sources, extensively characterized with conventional molecular methods, was used as a validation dataset. Using a validation strategy specifically adopted to WGS, we demonstrated high performance with repeatability, reproducibility, accuracy, precision, sensitivity and specificity above 95 % for the majority of all assays. The WGS workflow is publicly available as a ‘push-button’ pipeline at https://galaxy.sciensano.be. Our validation strategy and accompanying reference dataset consisting of both conventional and WGS data can be used for characterizing the performance of various bioinformatics workflows and assays, facilitating interoperability between laboratories with different WGS and bioinformatics set-ups

Evaluation of WGS performance for bacterial pathogen characterization with the Illumina technology optimized for time-critical situations

Whole genome sequencing (WGS) has become the reference standard for bacterial outbreak investigation and pathogen typing, providing a resolution unattainable with conventional molecular methods. Data generated with Illumina sequencers can however only be analysed after the sequencing run has finished, thereby losing valuable time during emergency situations. We evaluated both the effect of decreasing overall run time, and also a protocol to transfer and convert intermediary files generated by Illumina sequencers enabling real-time data analysis for multiple samples part of the same ongoing sequencing run, as soon as the forward reads have been sequenced. To facilitate implementation for laboratories operating under strict quality systems, extensive validation of several bioinformatics assays (16S rRNA species confirmation, gene detection against virulence factor and antimicrobial resistance databases, SNP-based antimicrobial resistance detection, serotype determination, and core genome multilocus sequence typing) for three bacterial pathogens (Mycobacterium tuberculosis, Neisseria meningitidis, and Shiga-toxin producing Escherichia coli) was performed by evaluating performance in function of the two most critical sequencing parameters, i.e. read length and coverage. For the majority of evaluated bioinformatics assays, actionable results could be obtained between 14 and 22 h of sequencing, decreasing the overall sequencing-to- results time by more than half. This study aids in reducing the turn-around time of WGS analysis by facilitating a faster response in time-critical scenarios and provides recommendations for time-optimized WGS with respect to required read length and coverage to achieve a minimum level of performance for the considered bioinformatics assay(s), which can also be used to maximize the cost-effectiveness of routine surveillance sequencing when response time is not essential.The Belgian Federal Public Service of Health, Food Chain Safety and Environment and Sciensano RP-PJ - Belgium.https://www.microbiologyresearch.org/content/journal/mgenam2022Genetic

Transforming Shiga toxin-producing Escherichia coli surveillance through whole genome sequencing in food safety practices

IntroductionShiga toxin-producing Escherichia coli (STEC) is a gastrointestinal pathogen causing foodborne outbreaks. Whole Genome Sequencing (WGS) in STEC surveillance holds promise in outbreak prevention and confinement, in broadening STEC epidemiology and in contributing to risk assessment and source attribution. However, despite international recommendations, WGS is often restricted to assist outbreak investigation and is not yet fully implemented in food safety surveillance across all European countries, in contrast to for example in the United States.MethodsIn this study, WGS was retrospectively applied to isolates collected within the context of Belgian food safety surveillance and combined with data from clinical isolates to evaluate its benefits. A cross-sector WGS-based collection of 754 strains from 1998 to 2020 was analyzed.ResultsWe confirmed that WGS in food safety surveillance allows accurate detection of genomic relationships between human cases and strains isolated from food samples, including those dispersed over time and geographical locations. Identifying these links can reveal new insights into outbreaks and direct epidemiological investigations to facilitate outbreak management. Complete WGS-based isolate characterization enabled expanding epidemiological insights related to circulating serotypes, virulence genes and antimicrobial resistance across different reservoirs. Moreover, associations between virulence genes and severe disease were determined by incorporating human metadata into the data analysis. Gaps in the surveillance system were identified and suggestions for optimization related to sample centralization, harmonizing isolation methods, and expanding sampling strategies were formulated.DiscussionThis study contributes to developing a representative WGS-based collection of circulating STEC strains and by illustrating its benefits, it aims to incite policymakers to support WGS uptake in food safety surveillance

Whole-genome sequencing-based antimicrobial resistance characterization and phylogenomic investigation of 19 multidrug-resistant and extended-spectrum beta-lactamase-positive Escherichia coli strains collected from hospital patients in Benin in 2019

The increasing worldwide prevalence of extended-spectrum beta-lactamase (ESBL) producing Escherichia coli constitutes a serious threat to global public health. Surgical site infections are associated with high morbidity and mortality rates in developing countries, fueled by the limited availability of effective antibiotics. We used whole-genome sequencing (WGS) to evaluate antimicrobial resistance and the phylogenomic relationships of 19 ESBL-positive E. coli isolates collected from surgical site infections in patients across public hospitals in Benin in 2019. Isolates were identified by MALDI-TOF mass spectrometry and phenotypically tested for susceptibility to 16 antibiotics. Core-genome multi-locus sequence typing and single-nucleotide polymorphism-based phylogenomic methods were used to investigate the relatedness between samples. The broader phylogenetic context was characterized through the inclusion of publicly available genome data. Among the 19 isolates, 13 different sequence types (STs) were observed, including ST131 (n = 2), ST38 (n = 2), ST410 (n = 2), ST405 (n = 2), ST617 (n = 2), and ST1193 (n = 2). The blaCTX-M-15 gene encoding ESBL resistance was found in 15 isolates (78.9%), as well as other genes associated with ESBL, such as blaOXA-1 (n = 14) and blaTEM-1 (n = 9). Additionally, we frequently observed genes encoding resistance against aminoglycosides [aac-(6')-Ib-cr, n = 14], quinolones (qnrS1, n = 4), tetracyclines [tet(B), n = 14], sulfonamides (sul2, n = 14), and trimethoprim (dfrA17, n = 13). Nonsynonymous chromosomal mutations in the housekeeping genes parC and gyrA associated with resistance to fluoroquinolones were also detected in multiple isolates. Although the phylogenomic investigation did not reveal evidence of hospital-acquired transmissions, we observed two very similar strains collected from patients in different hospitals. By characterizing a set of multidrug-resistant isolates collected from a largely unexplored environment, this study highlights the added value for WGS as an effective early warning system for emerging pathogens and antimicrobial resistance.The ARES (Académie de la Recherche pour l’Enseignement Supérieur), Belgium.http://www.frontiersin.org/Microbiologyam2022Genetic