Strain-Level Metagenomic Analysis of the Fermented Dairy Beverage Nunu Highlights Potential Food Safety Risks

ABSTRACT The rapid detection of pathogenic strains in food products is essential for the prevention of disease outbreaks. It has already been demonstrated that whole-metagenome shotgun sequencing can be used to detect pathogens in food but, until recently, strain-level detection of pathogens has relied on whole-metagenome assembly, which is a computationally demanding process. Here we demonstrated that three short-read-alignment-based methods, i.e., MetaMLST, PanPhlAn, and StrainPhlAn, could accurately and rapidly identify pathogenic strains in spinach metagenomes that had been intentionally spiked with Shiga toxin-producing Escherichia coli in a previous study. Subsequently, we employed the methods, in combination with other metagenomics approaches, to assess the safety of nunu, a traditional Ghanaian fermented milk product that is produced by the spontaneous fermentation of raw cow milk. We showed that nunu samples were frequently contaminated with bacteria associated with the bovine gut and, worryingly, we detected putatively pathogenic E. coli and Klebsiella pneumoniae strains in a subset of nunu samples. Ultimately, our work establishes that short-read-alignment-based bioinformatics approaches are suitable food safety tools, and we describe a real-life example of their utilization. IMPORTANCE Foodborne pathogens are responsible for millions of illnesses each year. Here we demonstrate that short-read-alignment-based bioinformatics tools can accurately and rapidly detect pathogenic strains in food products by using shotgun metagenomics data. The methods used here are considerably faster than both traditional culturing methods and alternative bioinformatics approaches that rely on metagenome assembly; therefore, they can potentially be used for more high-throughput food safety testing. Overall, our results suggest that whole-metagenome sequencing can be used as a practical food safety tool to prevent diseases or to link outbreaks to specific food products.

[1]  Pelin Yilmaz,et al.  The SILVA ribosomal RNA gene database project: improved data processing and web-based tools , 2012, Nucleic Acids Res..

[2]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[3]  Claudio Donati,et al.  MetaMLST: multi-locus strain-level bacterial typing from metagenomic samples , 2016, Nucleic acids research.

[4]  Richard J Ellis,et al.  Whole-genome sequencing for national surveillance of Shiga toxin-producing Escherichia coli O157. , 2015, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[5]  M. Wiedmann,et al.  Genomics tools in microbial food safety , 2015 .

[6]  Siu-Ming Yiu,et al.  IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth , 2012, Bioinform..

[7]  N. Loman,et al.  A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4. , 2013, JAMA.

[8]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[9]  Duy Tin Truong,et al.  MetaPhlAn2 for enhanced metagenomic taxonomic profiling , 2015, Nature Methods.

[10]  P. Griffin,et al.  An assessment of the human health impact of seven leading foodborne pathogens in the United States using disability adjusted life years , 2015, Epidemiology and Infection.

[11]  Susan R. Leonard,et al.  Strain-Level Discrimination of Shiga Toxin-Producing Escherichia coli in Spinach Using Metagenomic Sequencing , 2016, PloS one.

[12]  Paul D. Cotter,et al.  Thermus and the Pink Discoloration Defect in Cheese , 2016, mSystems.

[13]  Paul D. Cotter,et al.  Impacts of Seasonal Housing and Teat Preparation on Raw Milk Microbiota: a High-Throughput Sequencing Study , 2016, Applied and Environmental Microbiology.

[14]  Robin Patel,et al.  Reevaluation of Streptococcus bovis Endocarditis Cases from 1975 to 1985 by 16S Ribosomal DNA Sequence Analysis , 2002, Journal of Clinical Microbiology.

[15]  X. Zhao,et al.  Comparative genomics Lactobacillus reuteri from sourdough reveals adaptation of an intestinal symbiont to food fermentations , 2015, Scientific Reports.

[16]  Torsten Seemann,et al.  Prospective Whole-Genome Sequencing Enhances National Surveillance of Listeria monocytogenes , 2015, Journal of Clinical Microbiology.

[17]  P. Pevzner,et al.  metaSPAdes: a new versatile metagenomic assembler. , 2017, Genome research.

[18]  A. van Dorsselaer,et al.  Carcinogenic properties of proteins with pro-inflammatory activity from Streptococcus infantarius (formerly S.bovis). , 2004, Carcinogenesis.

[19]  K. Kakinuma,et al.  Phylogenetic Analysis of Salmonella, Shigella, and Escherichia coli Strains on the Basis of the gyrB Gene Sequence , 2002, Journal of Clinical Microbiology.

[20]  Bernard Henrissat,et al.  Metabolic Reconstruction for Metagenomic Data and Its Application to the Human Microbiome , 2012, PLoS Comput. Biol..

[21]  Bas E. Dutilh,et al.  SUPER-FOCUS: a tool for agile functional analysis of shotgun metagenomic data , 2015, Bioinform..

[22]  Susan R. Leonard,et al.  Application of Metagenomic Sequencing to Food Safety: Detection of Shiga Toxin-Producing Escherichia coli on Fresh Bagged Spinach , 2015, Applied and Environmental Microbiology.

[23]  R. P. Ross,et al.  Fermented beverages with health-promoting potential: Past and future perspectives , 2014 .

[24]  Uday Deshpande,et al.  Draft Genome Sequence of Enterotoxigenic Escherichia coli Strain E24377A, Obtained from a Tribal Drinking Water Source in India , 2015, Genome Announcements.

[25]  Jack A. Gilbert,et al.  Metatranscriptomics reveals temperature-driven functional changes in microbiome impacting cheese maturation rate , 2016, Scientific Reports.

[26]  F. Akabanda Microbiological Characteristics of Ghanaian Traditional Fermented Milk Product, Nunu , 2010 .

[27]  D. Ercolini,et al.  A Few Pseudomonas Oligotypes Dominate in the Meat and Dairy Processing Environment , 2017, Front. Microbiol..

[28]  Christina Boucher,et al.  Use of Metagenomic Shotgun Sequencing Technology To Detect Foodborne Pathogens within the Microbiome of the Beef Production Chain , 2016, Applied and Environmental Microbiology.

[29]  D. Nielsen,et al.  A traditional Sudanese fermented camel's milk product, Gariss, as a habitat of Streptococcus infantarius subsp. infantarius. , 2008, International journal of food microbiology.

[30]  Paul P. Gardner,et al.  An evaluation of the accuracy and speed of metagenome analysis tools , 2015, Scientific Reports.

[31]  William A. Walters,et al.  QIIME allows analysis of high-throughput community sequencing data , 2010, Nature Methods.

[32]  G. Funke,et al.  Comprehensive Study of Strains Previously Designated Streptococcus bovis Consecutively Isolated from Human Blood Cultures and Emended Description of Streptococcus gallolyticus and Streptococcus infantarius subsp. coli , 2008, Journal of Clinical Microbiology.

[33]  Marcus J. Claesson,et al.  Comparing Apples and Oranges?: Next Generation Sequencing and Its Impact on Microbiome Analysis , 2016, PloS one.

[34]  Wenjun Liu,et al.  Expanding the biotechnology potential of lactobacilli through comparative genomics of 213 strains and associated genera , 2015, Nature Communications.

[35]  Sylvie Mireille Kouamé-Sina,et al.  Prevalence and comparison of Streptococcus infantarius subsp. infantarius and Streptococcus gallolyticus subsp. macedonicus in raw and fermented dairy products from East and West Africa. , 2013, International journal of food microbiology.

[36]  J. Rothberg,et al.  Prospective Genomic Characterization of the German Enterohemorrhagic Escherichia coli O104:H4 Outbreak by Rapid Next Generation Sequencing Technology , 2011, PloS one.

[37]  F. Ryan,et al.  SPINGO: a rapid species-classifier for microbial amplicon sequences , 2015, BMC Bioinformatics.

[38]  Robert A. Edwards,et al.  Quality control and preprocessing of metagenomic datasets , 2011, Bioinform..

[39]  L. Jespersen,et al.  Taxonomic and molecular characterization of lactic acid bacteria and yeasts in nunu, a Ghanaian fermented milk product. , 2013, Food microbiology.

[40]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[41]  Rob Knight,et al.  PyNAST: a flexible tool for aligning sequences to a template alignment , 2009, Bioinform..

[42]  Susan Holmes,et al.  phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data , 2013, PloS one.

[43]  Duy Tin Truong,et al.  Strain-level microbial epidemiology and population genomics from shotgun metagenomics , 2016, Nature Methods.

[44]  Marcus J Claesson,et al.  Translating Omics to Food Microbiology. , 2017, Annual review of food science and technology.

[45]  J. Tiedje,et al.  Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy , 2007, Applied and Environmental Microbiology.

[46]  Robert C. Edgar,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2001 .

[47]  Marcus J. Claesson,et al.  Microbial Succession and Flavor Production in the Fermented Dairy Beverage Kefir , 2016, mSystems.

[48]  Timothy L. Tickle,et al.  Compact graphical representation of phylogenetic data and metadata with GraPhlAn , 2015, PeerJ.

[49]  Sharon L. Grim,et al.  Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data , 2013, Methods in ecology and evolution.

[50]  Duy Tin Truong,et al.  Studying Vertical Microbiome Transmission from Mothers to Infants by Strain-Level Metagenomic Profiling , 2016, mSystems.