Sequencing of whole microbial genomes has become a standard procedure for cluster detection, source tracking, outbreak investigation and surveillance of many microorganisms. An increasing number of laboratories are currently in a transition phase from classical methods towards next generation sequencing, generating unprecedented amounts of data. Since the precision of downstream analyses depends significantly on the quality of raw data generated on the sequencing instrument, a comprehensive, meaningful primary quality control is indispensable. Here, we present AQUAMIS, a Snakemake workflow for an extensive quality control and assembly of raw Illumina sequencing data, allowing laboratories to automatize the initial analysis of their microbial whole-genome sequencing data. AQUAMIS performs all steps of primary sequence analysis, consisting of read trimming, read quality control (QC), taxonomic classification, de-novo assembly, reference identification, assembly QC and contamination detection, both on the read and assembly level. The results are visualized in an interactive HTML report including species-specific QC thresholds, allowing non-bioinformaticians to assess the quality of sequencing experiments at a glance. All results are also available as a standard-compliant JSON file, facilitating easy downstream analyses and data exchange. We have applied AQUAMIS to analyze ~13,000 microbial isolates as well as ~1000 in-silico contaminated datasets, proving the workflow’s ability to perform in high throughput routine sequencing environments and reliably predict contaminations. We found that intergenus and intragenus contaminations can be detected most accurately using a combination of different QC metrics available within AQUAMIS.

Antimicrobial resistance (AMR) is a major threat to public health worldwide. Currently, AMR typing changes from phenotypic testing to whole-genome sequence (WGS)-based detection of resistance determinants for a better understanding of the isolate diversity and elements involved in gene transmission (e.g., plasmids, bacteriophages, transposons). However, the use of WGS data in monitoring purposes requires suitable techniques, standardized parameters and approved guidelines for reliable AMR gene detection and prediction of their association with mobile genetic elements (plasmids). In this study, different sequencing and assembly strategies were tested for their suitability in AMR monitoring in Escherichia coli in the routines of the German National Reference Laboratory for Antimicrobial Resistances. To assess the outcomes of the different approaches, results from in silico predictions were compared with conventional phenotypic- and genotypic-typing data. With the focus on (fluoro)quinolone-resistant E.coli, five qnrS-positive isolates with multiple extrachromosomal elements were subjected to WGS with NextSeq (Illumina), PacBio (Pacific BioSciences) and ONT (Oxford Nanopore) for in depth characterization of the qnrS1-carrying plasmids. Raw reads from short- and long-read sequencing were assembled individually by Unicycler or Flye or a combination of both (hybrid assembly). The generated contigs were subjected to bioinformatics analysis. Based on the generated data, assembly of long-read sequences are error prone and can yield in a loss of small plasmid genomes. In contrast, short-read sequencing was shown to be insufficient for the prediction of a linkage of AMR genes (e.g., qnrS1) to specific plasmid sequences. Furthermore, short-read sequencing failed to detect certain duplications and was unsuitable for genome finishing. Overall, the hybrid assembly led to the most comprehensive typing results, especially in predicting associations of AMR genes and mobile genetic elements. Thus, the use of different sequencing technologies and hybrid assemblies currently represents the best approach for reliable AMR typing and risk assessment.

In spring 2016, Greece reported an outbreak caused by a previously undescribed Salmonella enterica subsp. enterica serotype (antigenic formula 11:z41:e,n,z15) via the Epidemic Intelligence Information System for Food- and Waterborne Diseases and Zoonoses (EPIS-FWD), with epidemiological evidence for sesame products as presumptive vehicle. Subsequently, Germany, Czech Republic, Luxembourg and the United Kingdom (UK) reported infections with this novel serotype via EPIS-FWD. Concerned countries in collaboration with the European Centre for Disease Prevention and Control (ECDC) and European Food Safety Authority (EFSA) adopted a common outbreak case definition. An outbreak case was defined as a laboratory-confirmed notification of the novel Salmonella serotype. Between March 2016 and April 2017, 47 outbreak cases were notified (Greece: n = 22; Germany: n = 13; Czech Republic: n = 5; Luxembourg: n = 4; UK: n = 3). Whole genome sequencing revealed the very close genetic relatedness of isolates from all affected countries. Interviews focusing on sesame product consumption, suspicious food item testing and trace-back analysis following Salmonella spp. detection in food products identified a company in Greece where sesame seeds from different countries were processed. Through European collaboration, it was possible to identify and recall sesame spread as one contaminated food item serving as vehicle of infection and trace it back to its origin.

The aim of this study was to gain an overview of the genetic diversity of Salmonella found in wildlife in Germany. We were particularly interested in exploring whether wildlife acts as a reservoir of certain serovars/subtypes or antimicrobial resistance (AMR) genes. Moreover, we wanted to explore the potential of Salmonella in spreading from wildlife to livestock and humans. To answer these questions, we sequenced 260 Salmonella enterica subsp. enterica isolates sampled between 2002 and 2020 from wildlife across Germany, using short-read whole genome sequencing. We found, consistent with previous findings, that some Salmonella sequence types are associated with certain animal species, such as S. Choleraesuis ST145 with wild boar and S. Enteritidis ST183 with hedgehogs. Antibiotic resistance was detected in 14.2% of all isolates, with resistance against important WATCH group antibiotics present in a small number of isolates. We further found that wildlife isolates do not form separate phylogenetic clusters distant to isolates from domestic animals and foodstuff, thus indicating frequent transmission events between these reservoirs. Overall, our study shows that Salmonella in German wildlife are diverse, with a low AMR burden and close links to Salmonella populations of farm and food-production environments.

Non-O1 and non-O139 Vibrio cholerae (NOVC) can cause gastrointestinal infections in humans. Contaminated food, especially seafood, is an important source of human infections. In this study, the virulence potential of 63 NOVC strains isolated from retail seafood were characterized at the genotypic and phenotypic levels. Although no strain encoded the cholera toxin (CTX) and the toxin-coregulated pilus (TCP), several virulence factors, including the HlyA hemolysin, the cholix toxin ChxA, the heat-stable enterotoxin Stn, and genes coding for the type 3 and type 6 secretion systems, were detected. All strains showed hemolytic activity against human and sheep erythrocytes: 90% (n = 57) formed a strong biofilm, 52% (n = 33) were highly motile at 37 °C, and only 8% (n = 5) and 14% (n = 9) could resist ≥60% and ≥40% human serum, respectively. Biofilm formation and toxin regulation genes were also detected. cgMLST analysis demonstrated that NOVC strains from seafood cluster with clinical NOVC strains. Antimicrobial susceptibility testing (AST) results in the identification of five strains that developed non-wildtype phenotypes (medium and resistant) against the substances of the classes of beta-lactams (including penicillin, carbapenem, and cephalosporin), polymyxins, and sulphonamides. The phenotypic resistance pattern could be partially attributed to the acquired resistance determinants identified via in silico analysis. Our results showed differences in the virulence potential of the analyzed NOVC isolated from retail seafood products, which may be considered for further pathogenicity evaluation and the risk assessment of NOVC isolates in future seafood monitoring.

Non-O1/non-O139 Vibrio cholerae (NOVC) are ubiquitous in aquatic ecosystems. In rare cases, they can cause intestinal and extra-intestinal infections in human. This ability is associated with various virulence factors. The presence of NOVC in German North Sea and Baltic Sea was observed in previous studies. However, data on virulence characteristics are still scarce. Therefore, this work aimed to investigating the virulence potential of NOVC isolated in these two regions. In total, 31 NOVC strains were collected and subjected to whole genome sequencing. In silico analysis of the pathogenic potential was performed based on the detection of genes involved in colonization and virulence. Phenotypic assays, including biofilm formation, mobility and human serum resistance assays were applied for validation. Associated toxin genes (hlyA, rtxA, chxA and stn), pathogenicity islands (Vibrio pathogenicity island 2 (VPI-II) and Vibrio seventh pathogenicity island 2 (VSP-II)) and secretion systems (Type II, III and VI secretion system) were observed. A maximum likelihood analysis from shared core genes revealed a close relationship between clinical NOVCs published in NCBI and environmental strains from this study. NOVC strains are more mobile at 37 °C than at 25 °C, and 68% of the NOVC strains could form strong biofilms at both temperatures. All tested strains were able to lyse erythrocytes from both human and sheep blood. Additionally, one strain could survive up to 60% and seven strains up to 40% human serum at 37 °C. Overall, the genetic virulence profile as well as the phenotypic virulence characteristics of the investigated NOVC from the German North Sea and Baltic Sea suggest potential human pathogenicity.