Choice of library preparation affects sequence quality, genome assembly, and precise in silico prediction of virulence genes in shiga toxin-producing Escherichia coli

Precision metagenomics sequencing for food safety: hybrid assembly of Shiga toxin-producing Escherichia coli in enriched agricultural water

Culture-independent metagenomic sequencing of enriched agricultural water could expedite the detection and virulotyping of Shiga toxin-producing Escherichia coli (STEC). We previously determined the limits of a complete, closed metagenome-assembled genome (MAG) assembly and of a complete, fragmented MAG assembly for O157:H7 in enriched agricultural water using long reads (Oxford Nanopore Technologies, Oxford), which were 107 and 105 CFU/ml, respectively. However, the nanopore assemblies did not have enough accuracy to be used in Single Nucleotide Polymorphism (SNP) phylogenies and cannot be used for the precise identification of an outbreak STEC strain. The present study aimed to determine the limits of detection and assembly for STECs in enriched agricultural water by Illumina MiSeq sequencing technology alone, followed by establishing the limit of hybrid assembly with nanopore long-read sequencing using three different hybrid assemblers (SPAdes, Unicycler, and OPERA-MS). We also aimed to generate a genome with enough accuracy to be used in a SNP phylogeny. The classification of MiSeq and nanopore sequencing identified the same highly abundant species. Using the totality of the MiSeq output and a precision metagenomics approach in which the E. coli reads are binned before assembly, the limit of detection and assembly of STECs by MiSeq were determined to be 105 and 107 CFU/ml, respectively. While a complete, closed MAG could not be generated at any concentration, a complete, fragmented MAG was produced using the SPAdes assembler with an STEC concentration of at least 107 CFU/ml. At this concentration, hybrid assembled contigs aligned to the nanopore-assembled genome could be accurately placed in a neighbor-joining tree. The MiSeq limit of detection and assembly was less sensitive than nanopore sequencing, which was likely due to factors including the small starting material (50 vs. 1 μg) and the dilution of the library loaded on the cartridge. This pilot study demonstrates that MiSeq sequencing requires higher coverage in precision metagenomic samples; however, with sufficient concentration, STECs can be characterized and phylogeny can be accurately determined.

Results of the 2020 Genomic Proficiency Test for the network of European Union Reference Laboratory for Antimicrobial Resistance assessing whole-genome-sequencing capacities

The global surveillance and outbreak investigation of antimicrobial resistance (AMR) is amidst a paradigm shift from traditional biology to bioinformatics. This is due to developments in whole-genome-sequencing (WGS) technologies, bioinformatics tools, and reduced costs. The increased use of WGS is accompanied by challenges such as standardization, quality control (QC), and data sharing. Thus, there is global need for inter-laboratory WGS proficiency test (PT) schemes to evaluate laboratories' capacity to produce reliable genomic data. Here, we present the results of the first iteration of the Genomic PT (GPT) organized by the Global Capacity Building Group at the Technical University of Denmark in 2020. Participating laboratories sequenced two isolates and corresponding DNA of Salmonella enterica, Escherichia coli and Campylobacter coli, using WGS methodologies routinely employed at their laboratories. The participants' ability to obtain consistently good-quality WGS data was assessed based on several QC WGS metrics. A total of 21 laboratories from 21 European countries submitted WGS and meta-data. Most delivered high-quality sequence data with only two laboratories identified as overall underperforming. The QC metrics, N50 and number of contigs, were identified as good indicators for high-sequencing quality. We propose QC thresholds for N50 greater than 20 000 and 25 000 for Campylobacter coli and Escherichia coli, respectively, and number of contigs >200 bp greater than 225, 265 and 100 for Salmonella enterica, Escherichia coli and Campylobacter coli, respectively. The GPT2020 results confirm the importance of systematic QC procedures, ensuring the submission of reliable WGS data for surveillance and outbreak investigation to meet the requirements of the paradigm shift in methodology.

Sequencing of Enteric Bacteria: Library Preparation Procedure Matters for Accurate Identification and Characterization

Enzymatic library preparation kits are increasingly used for bacterial whole genome sequencing. While they offer a rapid workflow, the transposases used in the kits are recognized to be somewhat biased. The aim of this study was to optimize and validate a protocol for the Illumina DNA Prep kit (formerly Nextera DNA Flex) for sequencing enteric pathogens and compare its performance against the Nextera XT kit. One hundred forty-three strains of Campylobacter, Escherichia, Listeria, Salmonella, Shigella, and Vibrio were prepared with both methods and sequenced on the Illumina MiSeq using 300 and/or 500 cycle chemistries. Sequences were compared using core genome multilocus sequence typing (cgMLST), 7-gene multilocus sequence typing (MLST), and detection of markers encoding serotype, virulence, and antimicrobial resistance. Sequences for one Escherichia strain were downsampled to determine the minimum coverage required for the analyses. While organism-specific differences were observed, the Prep libraries generated longer average read lengths and less fragmented assemblies compared to the XT libraries. In downstream analysis, the most notable difference between the kits was observed for Escherichia, particularly for the 300 cycle sequences. The O group was not predicted in 32% and 4% of XT sequences when using blast and kmer algorithms, respectively, while the O group was predicted from all Prep sequences regardless of the algorithm. In addition, the ehxA gene was not detected in 6% of XT sequences and 34% were missing one or more of the type III secretion systems and/or plasmid-associated genes, which were detected in the Prep sequences. The coverage downsampling revealed that acceptable assembly quality and allele detection was achieved at 30 × coverage with the Prep libraries, whereas 40-50 × coverage was required for the XT libraries. The better performance of the Prep libraries was attributed to more even coverage, particularly in genome regions low in GC content.

Desiccation Survival of Salmonella enterica,Escherichia coli, and Enterococcus faecium Related to Initial Cell Level and Cellular Components

ABSTRACT

Salmonella enterica is well known for its ability to survive and persist in low-moisture environments. Previous studies have indicated a link between the initial cell level and the population of Salmonella that survives after desiccation and subsequent storage; however, how the initial cell concentration affects survival is unknown. This study was conducted to examine this phenomenon and to determine whether it occurs in other microorganisms, specifically Shiga toxigenic Escherichia coli (STEC) and Enterococcus faecium. Salmonella, STEC, and E. faecium were grown as sessile cells on Trypticase soy agar with yeast extract (TSAYE) and harvested in buffered peptone water (BPW). To determine recovery at different initial cell levels, cultures were diluted to 9, 7, and 5 log CFU/mL and applied to filters. Filters were dried for 24 h and then stored for 28 days at 25°C and 33% relative humidity. During storage, cells were recovered from filters with BPW and cultivated on TSAYE. Recovery of both Salmonella and E. coli, but not E. faecium, was nonproportional. Lower initial populations were less viable after 24 h of desiccation; ≥10 log CFU/mL was recovered when 11 log CFU/mL was desiccated, but <3 log CFU/mL was recovered when 5 log CFU/mL was desiccated. Once dried, persistence did not appear affected by initial cell concentration. When inactivated (heat-treated) cells were added to the diluent, recovery of Salmonella was proportional with respect to the initial cell level. To further examine the response to desiccation, Salmonella was diluted in BPW containing 1 of 11 test cell components related to quorum sensing or known to affect desiccation resistance to assess recovery and persistence. Of the 11 additions, only cell debris fractions, cell-free extract, and peptidoglycan improved recovery of Salmonella. Desiccation survival appears related to cell wall components; however, the exact mechanism affecting survival remains unknown.

HIGHLIGHTS

Evaluating whole-genome sequencing quality metrics for enteric pathogen outbreaks

Background

Whole genome sequencing (WGS) has gained increasing importance in responses to enteric bacterial outbreaks. Common analysis procedures for WGS, single nucleotide polymorphisms (SNPs) and genome assembly, are highly dependent upon WGS data quality.

Methods

Raw, unprocessed WGS reads from Escherichia coli, Salmonella enterica, and Shigella sonnei outbreak clusters were characterized for four quality metrics: PHRED score, read length, library insert size, and ambiguous nucleotide composition. PHRED scores were strongly correlated with improved SNPs analysis results in E. coli and S. enterica clusters.

Results

Assembly quality showed only moderate correlations with PHRED scores and library insert size, and then only for Salmonella. To improve SNP analyses and assemblies, we compared seven read-healing pipelines to improve these four quality metrics and to see how well they improved SNP analysis and genome assembly. The most effective read healing pipelines for SNPs analysis incorporated quality-based trimming, fixed-width trimming, or both. The Lyve-SET SNPs pipeline showed a more marked improvement than the CFSAN SNP Pipeline, but the latter performed better on raw, unhealed reads. For genome assembly, SPAdes enabled significant improvements in healed E. coli reads only, while Skesa yielded no significant improvements on healed reads.

Conclusions

PHRED scores will continue to be a crucial quality metric albeit not of equal impact across all types of analyses for all enteric bacteria. While trimming-based read healing performed well for SNPs analyses, different read healing approaches are likely needed for genome assembly or other, emerging WGS analysis methodologies.

Comparison of four enzymatic library preparation kits for sequencing Shiga toxin-producing Escherichia coli for surveillance and outbreak detection

Four enzymatic DNA library preparation kits were compared for sequencing Shiga toxin-producing E. coli. All kits produced high quality sequence data which performed equally well in the downstream analyses for surveillance and outbreak detection. Important differences were noted in the workflow user-friendliness and per sample cost.