Application of quasimetagenomics methods to define microbial diversity and subtype Listeria monocytogenes in dairy and seafood production facilities

ABSTRACT Microorganisms frequently colonize surfaces within food production facilities. Detection of Listeria monocytogenes in this setting relies on culture-dependent methods, but the complex dynamics of bacterial interactions within these environments and their impact on pathogen detection remain largely unexplored. To address this challenge, we applied both 16S rRNA and shotgun quasimetagenomic (enriched microbiome) sequencing of swab culture enrichments from five seafood and seven dairy production environments. Utilizing 16S rRNA amplicon sequencing, we observed variability between 355 samples taken from these 12 production facilities and a distinctive microbiome for each environment. With shotgun quasimetagenomic sequencing, we were able to assemble L. monocytogenes metagenome-assembled genomes (MAGs) from 28 of the 32 culture-positive samples. We compared these MAGs to their corresponding whole-genome sequencing assemblies, which resulted in two polyphyletic clades consisting of L. monocytogenes lineages I and II with 13,195 and 25,556 single-nucleotide polymorphism sites, respectively. The remaining four MAGs did not produce sufficient genome coverage. To understand and establish limits for pathogen detection and subtyping using shotgun quasimetagenomics, these same data sets were downsampled in slilico to produce a titration series of abundances of L. monocytogenes and analyzed. Pathogen detection was achieved for all downsampled data sets, even those with only 3× genome coverage. This study contributes to the understanding of microbial diversity within food production environments and presents insights into the level of genome coverage needed in a metagenome sequencing data set to detect, subtype, and source track a foodborne pathogen IMPORTANCE In developed countries, the human diet is predominated by food commodities, which have been manufactured, processed, and stored in a food production facility. Little is known about the application of metagenomic sequencing approaches for detecting foodborne pathogens, such as L. monocytogenes, and characterizing microbial diversity in food production ecosystems. In this work, we investigated the utility of 16S rRNA amplicon and quasimetagenomic sequencing for the taxonomic and phylogenetic classification of Listeria culture enrichments of environmental swabs collected from dairy and seafood production facilities. We demonstrated that single-nucleotide polymorphism (SNP) analyses of L. monocytogenes metagenome-assembled genomes (MAGs) from quasimetagenomic data sets can achieve similar resolution as culture isolate whole-genome sequencing. To further understand the impact of genome coverage on MAG SNP cluster resolution, an in silico downsampling approach was employed to reduce the percentage of target pathogen sequence reads, providing an initial estimate of required MAG coverage for subtyping resolution of L. monocytogenes. In developed countries, the human diet is predominated by food commodities, which have been manufactured, processed, and stored in a food production facility. Little is known about the application of metagenomic sequencing approaches for detecting foodborne pathogens, such as L. monocytogenes, and characterizing microbial diversity in food production ecosystems. In this work, we investigated the utility of 16S rRNA amplicon and quasimetagenomic sequencing for the taxonomic and phylogenetic classification of Listeria culture enrichments of environmental swabs collected from dairy and seafood production facilities. We demonstrated that single-nucleotide polymorphism (SNP) analyses of L. monocytogenes metagenome-assembled genomes (MAGs) from quasimetagenomic data sets can achieve similar resolution as culture isolate whole-genome sequencing. To further understand the impact of genome coverage on MAG SNP cluster resolution, an in silico downsampling approach was employed to reduce the percentage of target pathogen sequence reads, providing an initial estimate of required MAG coverage for subtyping resolution of L. monocytogenes.