Precision long-read metagenomics sequencing for food safety by detection and assembly of Shiga toxin-producing Escherichia coli in irrigation water

Shiga toxin-producing Escherichia coli (STEC) contamination of agricultural water might be an important factor to recent foodborne illness and outbreaks involving leafy greens. Whole genome sequencing generation of closed bacterial genomes plays an important role in source tracking. We aimed to determine the limits of detection and classification of STECs by qPCR and nanopore sequencing using enriched irrigation water artificially contaminated with E. coli O157:H7 (EDL933). We determined the limit of STEC detection by qPCR to be 30 CFU/reaction, which is equivalent to 105 CFU/ml in the enrichment. By using Oxford Nanopores EPI2ME WIMP workflow and de novo assembly with Flye followed by taxon classification with a k-mer analysis software (Kraken), E. coli O157:H7 could be detected at 103 CFU/ml (68 reads) and a complete fragmented E. coli O157:H7 metagenome-assembled genome (MAG) was obtained at 105-108 CFU/ml. Using a custom script to extract the E. coli reads, a completely closed MAG was obtained at 107-108 CFU/ml and a complete, fragmented MAG was obtained at 105-106 CFU/ml. In silico virulence detection for E. coli MAGs for 105-108 CFU/ml showed that the virulotype was indistinguishable from the spiked E. coli O157:H7 strain. We further identified the bacterial species in the un-spiked enrichment, including antimicrobial resistance genes, which could have important implications to food safety. We propose this workflow could be used for detection and complete genomic characterization of STEC from a complex microbial sample and could be applied to determine the limit of detection and assembly of other foodborne bacterial pathogens. IMPORTANCEFoodborne illness caused by Shiga toxin-producing E. coli (STEC) ranges in severity from diarrhea to hemolytic uremic syndrome and produce-related incidence is increasing. The pervasive nature of E. coli requires not only detection, but also a complete genome to determine potential pathogenicity based on stx and eae genes, serotype, and other virulence factors. We have developed a pipeline to determine the limits of nanopore sequencing for STECs in a metagenomic sample. By utilizing the current qPCR in the FDA Bacteriological Analytical Manual (BAM) Chapter 4A, we can quantify the amount of STEC in the enrichment and then sequence and classify the STEC in less than half the time as current protocols that require a single isolate. These methods have wide implications for food safety, including decreased time to STEC identification during outbreaks, characterization of the microbial community, and the potential to use these methods to determine the limits for other foodborne pathogens.