Search

For successful elucidation of a food-borne infection chain, the availability of high-quality sequencing data from suspected microbial contaminants is a prerequisite. Commonly, those investigations are a joint effort undertaken by different laboratories and institutes. To analyze the extent of variability introduced by differing wet-lab procedures on the quality of the sequence data we conducted an interlaboratory study, involving four bacterial pathogens, which account for the majority of food-related bacterial infections: Campylobacter spp., Shiga toxin-producing Escherichia coli , Listeria monocytogenes , and Salmonella enterica . The participants, ranging from German federal research institutes, federal state laboratories to universities and companies, were asked to follow their routine in-house protocols for short-read sequencing of 10 cultures and one isolated bacterial DNA per species. Sequence and assembly quality were then analyzed centrally. Variations within isolate samples were detected with SNP and cgMLST calling. Overall, we found that the quality of Illumina raw sequence data was high with little overall variability, with one exception, attributed to a specific library preparation kit. The variability of Ion Torrent data was higher, independent of the investigated species. For cgMLST and SNP analysis results, we found that technological sequencing artefacts could be reduced by the use of filters, and that SNP analysis was more suited than cgMLST to compare data of different contributors. Regarding the four species, a minority of Campylobacter isolate data showed the in comparison highest divergence with regard to sequence type and cgMLST analysis. We additionally compared the assembler SPAdes and SKESA for their performance on the Illumina data sets of the different species and library preparation methods and found overall similar assembly quality metrics and cgMLST statistics., Competing Interests: ClS was employed by the company Labor Kneißler GmbH & Co. KG. TH was emplyoed by StarSEQ GmbH. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Forth, Brinks, Denay, Fawzy, Fiedler, Fuchs, Geuthner, Hankeln, Hiller, Murr, Petersen, Reiting, Schäfers, Schwab, Szabo, Thürmer, Wöhlke, Fischer, Lüth, Projahn, Stingl, Borowiak, Deneke, Malorny and Uelze.)

Salmonella enterica subsp. enterica serovar Agona has a history of causing food-borne outbreaks and any emergence of multidrug-resistant (MDR) isolates in novel food products is of concern. Particularly, in food products frequently consumed without sufficient heating prior to consumption. Here, we report about the MDR isolate, 18-SA00377, which had been isolated from a dietary supplement in Germany in 2018 and submitted to the German National Reference Laboratory for Salmonella . WGS-based comparative genetic analyses were conducted to find a potential reservoir of the isolate itself or mobile genetic elements associated with MDR. As a phylogenetic analysis did not yield any closely related S. Agona isolates, either globally or from Germany, a detailed analysis of the largest plasmid (295,499 bp) was performed as it is the main carrier of resistances. A combined approach of long-read and short-read sequencing enabled the assembly of the isolate's chromosome and its four plasmids. Their characterization revealed the presence of 23 different antibiotic resistance genes (ARGs), conferring resistance to 12 different antibiotic drug classes, as well as genes conferring resistance to six different heavy metals. The largest plasmid, pSE18-SA00377-1, belongs to the IncHI2 plasmid family and carries 16 ARGs, that are organized as two distinct clusters, with each ARG associated with putative composite transposons. Through a two-pronged approach, highly similar plasmids to pSE18-SA00377-1 were identified in the NCBI database and a search for Salmonella isolates with a highly similar ARG resistance profile was conducted. Mapping and structural comparisons between pSE18-SA00377-1 and these plasmids and Salmonella isolates showed that both the plasmid backbone and identical or similar ARG clusters can be found not only in Salmonella isolates, originating mostly from a wide variety of livestock, but also in a diverse range of bacterial genera of varying geographical origins and isolation sources. Thus, it can be speculated that the host range of pSE18-SA00377-1 is not restricted to Salmonella and its spread already occurred in different bacterial populations. Overall, this hints at a complex history for pSE18-SA00377-1 and highlights the importance of surveilling multidrug-resistant S. enterica isolates, especially in novel food items that are not yet heavily regulated., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Bartsch, Borowiak, Deneke, Gruetzke, Hammerl, Malorny, Szabo, Alter, Nguyen and Fischer.)

Caenorhabditis ( C .) elegans has become a popular toxicological and biological test organism in the last two decades. Furthermore, the role of C. elegans as an alternative for replacing or reducing animal experiments is continuously discussed and investigated. In the current study, we investigated whether C. elegans survival assays can help in determining differences in the virulence of Salmonella enterica strains and to what extent C. elegans assays could replace animal experiments for this purpose. We focused on three currently discussed examples where we compared the longevity of C. elegans when fed (i) with S . enterica serovar Enteritidis vaccination or wild-type strains, (ii) with lipopolysaccharide (LPS) deficient rough or LPS forming smooth S . enterica serovar Enteritidis, and (iii) with an S. enterica subsp. diarizonae strain in the presence or absence of the typical pSASd plasmid encoding a bundle of putative virulence factors. We found that the C. elegans survival assay could indicate differences in the longevity of C. elegans when fed with the compared strain pairs to a certain extent. Putatively higher virulent S. enterica strains reduced the lifespan of C. elegans to a greater extent than putatively less virulent strains. The C. elegans survival assay is an effective and relatively easy method for classifying the virulence of different bacterial isolates in vivo , but it has some limitations. The assay cannot replace animal experiments designed to determine differences in the virulence of Salmonella enterica strains. Instead, we recommend using the described method for pre-screening bacterial strains of interest to select the most promising candidates for further animal experiments. The C. elegans assay possesses the potential to reduce the number of animal experiments. Further development of the C. elegans assay in conjunction with omics technologies, such as transcriptomics, could refine results relating to the estimation of the virulent potential of test organisms., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Burkhardt, Salzinger, Fischer, Malorny, Fischer and Szabo.)

Whole-genome sequencing (WGS)-based outbreak investigation has proven to be a valuable method for the surveillance of bacterial pathogens. Its utility has been successfully demonstrated using both gene-by-gene (cgMLST or wgMLST) and single-nucleotide polymorphism (SNP)-based approaches. Among the obstacles of implementing a WGS-based routine surveillance is the need for an exchange of large volumes of sequencing data, as well as a widespread reluctance to share sequence and metadata in public repositories, together with a lacking standardization of suitable bioinformatic tools and workflows. To address these issues, we present chewieSnake , an intuitive and simple-to-use cgMLST workflow. ChewieSnake builds on the allele calling software chewBBACA and extends it by the concept of allele hashing. The resulting hashed allele profiles can be readily compared between laboratories without the need of a central allele nomenclature. The workflow fully automates the computation of the allele distance matrix, cluster membership, and phylogeny and summarizes all important findings in an interactive HTML report. Furthermore, chewieSnake can join allele profiles generated at different laboratories and identify shared clusters, including a stable and intercommunicable cluster nomenclature, thus facilitating a joint outbreak investigation. We demonstrate the feasibility of the proposed approach with a thorough method comparison using publically available sequencing data for Salmonella enterica . However, chewieSnake is readily applicable to all bacterial taxa, provided that a suitable cgMLST scheme is available. The workflow is freely available as an open-source tool and can be easily installed via conda or docker., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Deneke, Uelze, Brendebach, Tausch and Malorny.)

Despite extensive monitoring programs and preventative measures, Salmonella spp. continue to cause tens of thousands human infections per year, as well as many regional and international food-borne outbreaks, that are of great importance for public health and cause significant socio-economic costs. In Germany, salmonellosis is the second most common cause of bacterial diarrhea in humans and is associated with high hospitalization rates. Whole-genome sequencing (WGS) combined with data analysis is a high throughput technology with an unprecedented discriminatory power, which is particularly well suited for targeted pathogen monitoring, rapid cluster detection and assignment of possible infection sources. However, an effective implementation of WGS methods for large-scale microbial pathogen detection and surveillance has been hampered by the lack of standardized methods, uniform quality criteria and strategies for data sharing, all of which are essential for a successful interpretation of sequencing data from different sources. To overcome these challenges, the national GenoSalmSurv project aims to establish a working model for an integrated genome-based surveillance system of Salmonella spp. in Germany, based on a decentralized data analysis. Backbone of the model is the harmonization of laboratory procedures and sequencing protocols, the implementation of open-source bioinformatics tools for data analysis at each institution and the establishment of routine practices for cross-sectoral data sharing for a uniform result interpretation. With this model, we present a working solution for cross-sector interpretation of sequencing data from different sources (such as human, veterinarian, food, feed and environmental) and outline how a decentralized data analysis can contribute to a uniform cluster detection and facilitate outbreak investigations., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Uelze, Becker, Borowiak, Busch, Dangel, Deneke, Fischer, Flieger, Hepner, Huber, Methner, Linde, Pietsch, Simon, Sing, Tausch, Szabo and Malorny.)

We compared the consistency, accuracy and reproducibility of next-generation short read sequencing between ten laboratories involved in food safety (research institutes, state laboratories, universities and companies) from Germany and Austria. Participants were asked to sequence six DNA samples of three bacterial species ( Campylobacter jejuni , Listeria monocytogenes and Salmonella enterica ) in duplicate, according to their routine in-house sequencing protocol. Four different types of Illumina sequencing platforms (MiSeq, NextSeq, iSeq, NovaSeq) and one Ion Torrent sequencing instrument (S5) were involved in the study. Sequence quality parameters were determined for all data sets and centrally compared between laboratories. SNP and cgMLST calling were performed to assess the reproducibility of sequence data collected for individual samples. Overall, we found Illumina short read data to be more accurate (higher base calling accuracy, fewer miss-assemblies) and consistent (little variability between independent sequencing runs within a laboratory) than Ion Torrent sequence data, with little variation between the different Illumina instruments. Two laboratories with Illumina instruments submitted sequence data with lower quality, probably due to the use of a library preparation kit, which shows difficulty in sequencing low GC genome regions. Differences in data quality were more evident after assembling short reads into genome assemblies, with Ion Torrent assemblies featuring a great number of allele differences to Illumina assemblies. Clonality of samples was confirmed through SNP calling, which proved to be a more suitable method for an integrated data analysis of Illumina and Ion Torrent data sets in this study., (Copyright © 2020 Uelze, Borowiak, Bönn, Brinks, Deneke, Hankeln, Kleta, Murr, Stingl, Szabo, Tausch, Wöhlke and Malorny.)

The polymyxin antibiotic colistin has been used in decades for treatment and prevention of infectious diseases in livestock. Nowadays, it is even considered as last-line treatment option for severe human infections caused by multidrug- and carbapenem-resistant Gram-negative bacteria. Therefore, the discovery of plasmid-mediated mobile colistin resistance ( mcr ) genes raised major public health concern. The aim of our study was to analyze colistin-resistant Salmonella enterica strains from animals, food, feed and the environment collected at the National Reference Laboratory for Salmonella in Germany on the presence of mcr-1 to mcr-9 genes. Altogether 407 colistin-resistant (MIC >2 mg/L) Salmonella isolates received between 2011 and 2018 were selected and screened by PCR using a published mcr-1 to mcr-5 as well as a newly developed mcr-6 to mcr-9 multiplex PCR protocol. 254 of 407 (62.4%) isolates harbored either mcr-1 ( n = 175), mcr-4 ( n = 53), mcr-5 ( n = 18) or mcr-1 and mcr-9 ( n = 8). The number of mcr -positive isolates ranged from 19 (2017) to 64 (2012) per year. WGS revealed that none of our isolates harbored the mcr-9.1 gene. Instead, two novel mcr-9 variants were observed, which both were affected by frameshift mutations and are probably non-functional. The mcr -harboring isolates were mainly derived from animals (77.2%) or food (20.1%) and could be assigned to ten different Salmonella serovars. Many of the isolates were multidrug-resistant. Co-occurrence of mcr-1 and AmpC or ESBL genes was observed in eight isolates. Our findings suggest that mcr genes are widely spread among colistin-resistant Salmonella isolates from livestock and food in Germany. Potential transfer of mcr -harboring isolates along the food chain has to be considered critically., (Copyright © 2020 Borowiak, Baumann, Fischer, Thomas, Deneke, Hammerl, Szabo and Malorny.)

In food safety the detection of food contaminations with pathogenic microorganisms is a race against time and often outpaced by error-prone epidemiological approaches. For evidence-based outbreak investigations fast and reliable techniques and procedures are required to identify the source of infection. Metagenomics has the potential to become a powerful tool in the field of modern food safety, since it allows the detection, identification and characterization of a broad range of pathogens in a single experiment without pre-cultivation within a couple of days. Nevertheless, sample handling, sequencing and data analysis are challenging and can introduce errors and biases into the analysis. In order to evaluate the potential of metagenomics in food safety, we generated a mock community containing DNA of foodborne bacteria. Herewith, we compare the aptitude of the two prevalent approaches - 16S rDNA amplicon sequencing and whole genome shotgun sequencing - for the detection of foodborne bacteria using different parameters during sample preparation, sequencing and data analysis. 16S rDNA sequencing did not only result in high deviations from the expected sample composition on genus and species level, but more importantly lacked the detection of several pathogenic species. While shotgun sequencing is more suitable for species detection, abundance estimation, genome assembly and species characterization, the performance can vary depending on the library preparation kit, which was confirmed for a naturally Francisella tularensis contaminated game meat sample. The application of the Nextera XT DNA Library Preparation Kit for shotgun sequencing did not only result in lower reference genome recovery and coverage, but also in distortions of the mock community composition. For data analysis, we propose a publicly available workflow for pathogen detection and characterization and demonstrate its benefits on the usability of metagenomic sequencing in food safety by analyzing an authentic metagenomic sample.

Salmonella enterica serovar 4,[5],12:i:- is a monophasic variant of S. Typhimurium incapable of expressing the second-phase flagellar antigen (fljAB operon), and it is recognized to be one of the most prevalent serovars causing human infections. A clonal lineage characterized by phage type DT193, PulseNet PFGE profile STYMXB.0131 and multidrug resistance to ampicillin, streptomycin, sulphonamides and tetracycline (R-type ASSuT) is commonly circulating in Europe. In this study we determined the deletions affecting the fljAB operon and the resistance region responsible for the R-type ASSuT in a strain of Salmonella enterica serovar 4,5,12:i:- DT193/STYMXB.0131, through an approach based on PCRs and Southern blot hybridization of genomic DNA. Using a set of nine specific PCRs, the prevalence of the resistance region was assessed in a collection of 144 S. enterica serovar 4,[5],12:i:-/ASSuT/STYMXB.0131 strains isolated from Germany, Switzerland and Italy. A 28 kb-region is embedded between the loci STM2759 and iroB, replacing the DNA located in between, including the fljAB operon. It encompasses the genes bla TEM-1, strA-strB, sul2 and tet(B) responsible for the R-type ASSuT together with genes involved in plasmid replication and orfs of unknown function characteristically located on IncH1 plasmids. Its location and internal structure is fairly conserved in S. enterica serovar 4,[5],12:i:-/ASSuT/STYMXB.0131 strains regardless of the isolation source or country. Hence, in the S. enterica serovar 4,[5],12:i:-/ASSuT/STYMXB.0131 clonal lineage widespread in Germany, Switzerland and Italy, a resistance region derived from IncH1 plasmids has replaced the chromosomal region encoding the second flagellar phase and is an example of the stabilization of new plasmid-derived genetic material due to integration into the bacterial chromosome.