- Research article
- Open Access
Quasimetagenomic source tracking of Listeria monocytogenes from naturally contaminated ice cream
BMC Infectious Diseases volume 20, Article number: 83 (2020)
The more quickly bacterial pathogens responsible for foodborne illness outbreaks can be linked to a vehicle of transmission or a source, the more illnesses can be prevented. Whole genome sequencing (WGS) based approaches to source tracking have greatly increased the speed and resolution with which public health response can pinpoint the vehicle and source of outbreaks. Traditionally, WGS approaches have focused on the culture of an individual isolate before proceeding to DNA extraction and sequencing. For Listeria monocytogenes (Lm), generation of an individual isolate for sequencing typically takes about 6 days. Here we demonstrate that a hybrid, “quasimetagenomic” approach ie; direct sequencing of microbiological enrichments (first step in pathogen detection and recovery) can provide high resolution source tracking sequence data, 5 days earlier than response that focuses on culture and sequencing of an individual isolate. This expedited approach could save lives, prevent illnesses and potentially minimize unnecessary destruction of food.
Naturally contaminated ice cream (from a 2015 outbreak) was enriched to recover Listeria monocytogenes following protocols outlined in the Bacteriological Analytic Manual (BAM). DNA from enriching microbiota was extracted and sequenced at incremental time-points during the first 48 h of pre-enrichment using the Illumina MiSeq platform (2 by 250), to evaluate genomic coverage of target pathogen, Listeria monocytogenes.
Quasimetagenomic sequence data acquired from hour 20 were sufficient to discern whether or not Lm strain/s were part of the ongoing outbreak or not. Genomic data from hours 24, 28, 32, 36, 40, 44, and 48 of pre-enrichments all provided identical phylogenetic source tracking utility to the WGS of individual isolates (which require an additional 5 days to culture).
The speed of this approach (more than twice as fast as current methods) has the potential to reduce the number of illnesses associated with any given outbreak by as many as 75% percent of total cases and potentially with continued optimization of the entire chain of response, contribute to minimized food waste.
Next generation sequencing (NGS) technologies have revolutionized our ability to source track pathogens such as Salmonella enterica , Listeria monocytogenes , Escherichia coli , and many others by facilitating the whole genome sequencing (WGS) of these small bacterial genomes in as little as 20h hours. Sequence data is organized using pipelines that generate matrices of single nucleotide polymorphisms (SNPs) that can be used with a variety of algorithms to infer phylogenetic relationships. This approach can differentiate highly clonal strains within serovars that differ by as few as 1 to 29 nucleotides . GenomeTrakr is the first of its kind network, organized by the U.S. Food and Drug Agency (FDA), the National Center for Biotechnology Information (NCBI) and the Centers for Disease Control (CDC), to standardize and coordinate the collection and sharing of WGS data with the goal of rapidly identifying origins of pathogens associated with illness outbreaks . The biggest time constraint for WGS source tracking is simply the culturing of the bacteria from which DNA is to be extracted and sequenced. Typical recovery of Lm from a food source, as described in the U.S. Food and Drug Administration (FDA) Bacteriological Analytical Manual (BAM) begins with 48 h of pre-enrichment, followed by incubation on selective agars (another 48 h), growth on nutrient agars (24 h) and then confirmation (3 to 24 h) of individually selected colonies . Time commitment for the full protocol is between120 h and 144 h (~ 6 days). Culture independent approaches are an alternative but still require financially and computationally impractical amounts of sequence data to achieve sufficient genomic coverage of low abundant targets that occur amidst a complex microbial ecology.
To evaluate hybrid options for improving the speed of source tracking, we sequenced DNA from enrichment microbiomes (every 4 hours) using short read (Illumina Miseq) sequencing. Using naturally contaminated ice cream previously linked to an outbreak of listeriosis that resulted in three deaths , shotgun data from incremental time-points over the first 48 h of pre-enrichment were evaluated to identify when Lm genomic data with source tracking utility could be obtained from the quasimetagenomes. Because the word “metagenomics” refers to true culture independent (CI) genomics , we use the term “quasimetagenomics”  to describe sequence data from pre-enrichments, which by their nature, are only partial selections of the environmental pan-genome.
Enrichment, DNA extraction and sequencing
Homogenized ice cream scoop samples were added to Buffered Listeria Enrichment Broth (BLEB) according to the specifications for Lm recovery outlined in the BAM . Four replicates of negative (no ice-cream) and positive controls (Lm cells) were assayed at all time-points.
DNA was extracted using DNeasy Blood and Tissue kit (Qiagen) following the protocol for Gram-positive bacteria with minor modifications: 1.5 ml of the culture was pelleted (5000×g, 15 min) and the pellet resuspended in 200 μL of enzymatic lysis buffer containing 20 mM Tris-HCl (pH -8.0), 2 mM Sodium EDTA, 1.2% Triton X- 100, 20 mg/ml of lysozyme. The samples were incubated for 60 min at 37 °C.
Libraries were prepared with Nextera XT (Illumina) according to the manufacturers specifications. For the first and second enrichments, approximately 16 libraries (independent replicates) were multiplexed per run (total of 3 runs) using Illumina MiSeq V2 (2 by 250).
WGS analyses were performed using the CFSAN SNP pipeline 0.6.1 [9, 10]. Raw reads from each replicate from all hours were mapped to the complete genome of CFSAN029793, (BioSample SAMN03386937) using default settings within Bowtie2 v2.2.2 . A resulting BAM (xyz) file was sorted using Samtools v1.3.1  and a pileup file for each isolate was produced. These files were then processed using VarScan2 v2.3.9  to identify high quality variant sites, using the mpileup2snp option. Additional information about these procedures, e.g. codes and instructions, is available at https://github.com/CFSAN-Biostatistics/snp-pipeline. Phylogenetic trees were created using the Genetic Algorithm for Rapid Likelihood Inference (GARLI) . Taxonomic profiles of enrichment microbiomes were assigned using CosmosID’s bacterial database and inhouse – k-mer based bacterial databases and pipelines.
All data is available at NCBI associated with BioProject PRJNA370011.
Evaluating genomic coverage of Listeria monocytogenes for trace-back utility
At h 20, from an average of 2.4 million sequences per replicate, an average of 1.68x coverage was obtained for Lm genome/s from ice-cream enrichments (Table 1). This coverage was sufficient to identify whether or not the candidate strain/s of Lm were part of the outbreak cluster or not. Determination of inclusion or exclusion after only 20 h (plus sequencing time) in cases where a reference genome is available facilitates a greatly expedited ability to remove contaminated commodities from the food supply and protect safe products from costly holds. High quality draft Lm genomes were achieved with data from h 24 through h 48 with coverage (12.74 to > 52.75x) -shown for hours 20–40 in Table 1. Lm genomes from enrichments of all time-points from h 20 through h 48 were used with WGS genomes of environmental and clinical isolates from the full outbreak cluster (available at NCBI) (Fig. 1a). An abbreviated representational set of clinical, food and environmental isolates is also shown in Fig. 1b with quasimetagenomes from enrichments of hours 20, 36 and 40. A list of genomes comprising the complete outbreak is available in Additional file 1.
Relative abundance of Listeria monocytogenes and other co-occurring genera
As microbiota from the naturally contaminated ice cream grew in the FDA BAM delineated media (BLEB), the following predominant bacterial genera were observed; Thermus (Thermales, Deinococcus-Thermus), Anoxybacillus, and Geobacillus (Bacillales, Firmicutes), and Lactococcus, Enterococcus and Streptococcus (Lactobacillales, Firmicutes) (Fig. 2). Lm remained at very low levels until h 20, introducing the possibility that the co-culturing Bacilli genera (Anoxybacillus and Geobacillus) play a role in inhibiting Listeria monocytogenes until it proliferates in h 20. This same phenomenon was also noted in previous work that focused on the description of the total ecology of enriching microbiota from naturally contaminated ice-cream using the three most commonly used enrichment protocols by; the International Organization for Standardization (ISO), the United States Department of Agriculture (USDA) and the FDA .
Metagenomic approaches have been accidentally and intentionally employed for decades to describe “black boxes” of human and environmental microbiota. These inquiries coupled with rapidly advancing NGS technologies have contributed to a renaissance in our understanding of numerous complex, previously poorly understood microbial ecologies. Focusing this approach on the “black box” of culture based pathogen recovery methods, many of which have not changed in 30+ years, provides novel insight to further expedite and optimize state of the art pathogen recovery methods. This applies both to the microbiological dynamics, ie; can we identify chemical, temperature or antibiotic inputs that would optimize recovery of targets by inhibiting competitors. Or, as presented here, by direct capture of genomic data of enriching targets from within the enrichment microbial community (quasimetagenome).
A true CI metagenomic approach was used to successfully identify strains of Shigatoxigenic E. coli O104: H4 linked to a 2011 outbreak that resulted in over 40 deaths . This powerful demonstration relied on extremely high throughput sequencing and extensive bioinformatic analyses. The authors describe the exciting potential, while acknowledging the significant challenges that remain, such as “speeding up and simplifying workflows, reducing costs, and improving diagnostic sensitivity”. Thus, the application of metagenomics for source tracking is valuable for certain inquiries but is not yet practical as a rapid, inexpensive approach that can be implemented in field and public health laboratories. The high costs of generating and analyzing sufficient genomic data to distinguish between the highly clonal organisms responsible for many foodborne illnesses such as E. coli, Salmonella enterica and Listeria monocytogenes, when they occur as low abundant members of highly diverse ecologies, are still prohibitive for most operations. Thus, validation of the quasimetagenomic approach (culture and shotgun sequencing) presented here, offers a protocol that can be easily integrated into public health laboratories at no additional cost to current programs.
The benefits achieved by the quasimetagenomic (qMGS) approach presented here, efficiently address previously described challenges of cost and speed, while maintaining diagnostic sensitivity. The source tracking phylogenetic success achieved by this experiment is particularly valuable because it validates the hybrid approach using naturally contaminated food samples. With modeled (spiked in) experiments, it is difficult to understand if results represent “real life” contamination dynamics or not. Quasimetagenomic sequencing can be applied to all culturable pathogens from all types of matrices (food or other). Ultimately this approach might have incredibly significant utility in hospital situations for rapid identification of strains causing septicemia. Temporal parameters, enrichment schema, and quantity of data necessary for optimal genomic coverage will be different for every pathogen, substrate and matrix. As shown here for source tracking Listeria monocytogenes from ice-cream the qMGS approach provides the exact same answer as WGS in a timeframe that when leveraged efficiently, could potentially eliminate the majority of illnesses associated with many outbreaks and potentially reduce unnecessary destruction of food.
Availability of data and materials
All data is available at NCBI associated with BioProject PRJNA370011 https://www.ncbi.nlm.nih.gov/bioproject/PRJNA370011.
Bacteriolgical Analytical Manual
Buffered Listeria enrichment broth
Centers for Disease Control
Center for Food Safety and Applied Nutrition
Food and Drug Administration
Genetic Algorithm for Rapid Likelihood Inference
Hydrogen chloride, Hydrochloric acid
International Organization for Standardization
- Lm :
National Center for Biotechnology Information
Next generation sequencing
Single nucleotide polymorphism
United States Department of Agriculture
Whole genome sequencing
Lienau EK, et al. Identification of a salmonellosis outbreak by means of molecular sequencing. N Engl J Med. 2011;364:981–2.
Jensen AK, et al. Whole-genome sequencing used to investigate a nationwide outbreak of listeriosis caused by ready-to-eat delicatessen meat, Denmark, 2014. Clin Infect Dis. 2016;63:64–70.
Joensen KG, et al. Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J Clin Microbiol. 2014;52:1501–10.
Chen Y, et al. Assessing the genome level diversity of Listeria monocytogenes from contaminated ice cream and environmental samples linked to a listeriosis outbreak in the United States. PLoS One. 2017;12:e0171389.
Allard MW, et al. The practical value of food pathogen traceability through building a whole-genome sequencing network and database. J Clin Microbiol. 2016;54:1975–83.
Hitchens A, Jinneman K, Chen Y. Detection and enumeration of Listeria monocyogenes in foods. Bacteriol Anal Manual. 2016;10.
Handelsman J. Metagenomics: application of genomics to uncultured microorganisms. Microbiol Mol Biol Rev. 2004;68:669–85.
Ottesen A, Ramachandran P. Food Microbiomes: A New Paradigm for Food and Food Ecology. Food Microbiology: Fundamentals and Frontiers 5th Ed. ASM Press. 2019;37:963–70.
Davis S, et al. CFSAN SNP pipeline: an automated method for constructing SNP matrices from next-generation sequence data. Peer J Comput Sci. 2015;1:e20.
Pettengill JB, et al. An evaluation of alternative methods for constructing phylogenies from whole genome sequence data: a case study with Salmonella. Peer J. 2014;2:e620.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:1.
Li H, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Koboldt DC, et al. VarScan: variant detection in massively parallel sequencing of individual and pooled samples. Bioinformatics. 2009;25:2283–5.
Gutell RR, Jansen RK. Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion; 2006.
Ottesen A, Ramachandran P, Reed E, White J, Chen Y. Enrichment dynamics of Listeria monocytogenes and the associated microbiome from naturally contaminated ice cream linked to a listeriosis outbreak. BMC Microbiol. 2016;16:275.
Loman NJ, et al. A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104: H4. JAMA. 2013;309:1502–10.
We would like to acknowledge Dr. John Callahan and Dr. Vincent Bunning for their tireless efforts in keeping the science of the Office of Regulatory Science at CFSAN moving forward.
We would like to acknowledge the Division of Microbiology, Office of Regulatory Science, Center for Food Safety and Applied Nutrition (CFSAN), FDA for all funding.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ottesen, A., Ramachandran, P., Chen, Y. et al. Quasimetagenomic source tracking of Listeria monocytogenes from naturally contaminated ice cream. BMC Infect Dis 20, 83 (2020). https://doi.org/10.1186/s12879-019-4747-z
- Source tracking
- Listeria monocytogenes
- Ice cream
- Microbiological enrichment