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

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.

Development of a targeted amplicon sequencing method for genotyping Cyclospora cayetanensis from fresh produce and clinical samples with enhanced genomic resolution and sensitivity.

Outbreaks of cyclosporiasis, an enteric illness caused by the parasite Cyclospora cayetanensis, have been associated with consumption of various types of fresh produce. Although a method is in use for genotyping C. cayetanensis from clinical specimens, the very low abundance of C. cayetanensis in food and environmental samples presents a greater challenge. To complement epidemiological investigations, a molecular surveillance tool is needed for use in genetic linkage of food vehicles to cyclosporiasis illnesses, estimation of the scope of outbreaks or clusters of illness, and determination of geographical areas involved. We developed a targeted amplicon sequencing (TAS) assay that incorporates a further enrichment step to gain the requisite sensitivity for genotyping C. cayetanensis contaminating fresh produce samples. The TAS assay targets 52 loci, 49 of which are located in the nuclear genome, and encompasses 396 currently known SNP sites. The performance of the TAS assay was evaluated using lettuce, basil, cilantro, salad mix, and blackberries inoculated with C. cayetanensis oocysts. A minimum of 24 markers were haplotyped even at low contamination levels of 10 oocysts in 25 g leafy greens. The artificially contaminated fresh produce samples were included in a genetic distance analysis based on haplotype presence/absence with publicly available C. cayetanensis whole genome sequence assemblies. Oocysts from two different sources were used for inoculation, and samples receiving the same oocyst preparation clustered together, but separately from the other group, demonstrating the utility of the assay for genetically linking samples. Clinical fecal samples with low parasite loads were also successfully genotyped. This work represents a significant advance in the ability to genotype C. cayetanensis contaminating fresh produce along with greatly expanding the genomic diversity included for genetic clustering of clinical specimens.

Weather factors, soil microbiome, and bacteria-fungi interactions as drivers of the epiphytic phyllosphere communities of romaine lettuce.

Lettuce is associated with seasonal outbreaks of Shiga toxin-producing Escherichia coli (STEC) infections. Little is known about how various biotic and abiotic factors affect the lettuce microbiome, which in turn impacts STEC colonization. We characterized the lettuce phyllosphere and surface soil bacterial, fungal, and oomycete communities at harvest in late-spring and -fall in California using metagenomics. Harvest season and field type, but not cultivar, significantly influenced the microbiome composition of leaves and surface soil near plants. Phyllosphere and soil microbiome compositions were correlated with specific weather factors. The relative abundance of Enterobacteriaceae, but not E. coli, was enriched on leaves (5.2%) compared to soil (0.4%) and correlated positively with minimum air temperature and wind speed. Co-occurrence networks revealed seasonal trends in fungi-bacteria interactions on leaves. These associations represented 39%-44% of the correlations between species. All significant E. coli co-occurrences with fungi were positive, while all negative associations were with bacteria. A large proportion of the leaf bacterial species was shared with those in soil, indicating microbiome transmission from the soil surface to the canopy. Our findings provide new insight into factors that shape lettuce microbial communities and the microbial context of foodborne pathogen immigration events in the lettuce phyllosphere.

Development and Single Laboratory Evaluation of a Refined and specific Real-time PCR Detection Method, Using Mitochondrial Primers (Mit1C), for the Detection of Cyclospora cayetanensis in Produce.

Regulatory methods for detection of the foodborne protozoan parasite Cyclospora cayetanensis must be specific and sensitive. To that end, we designed and evaluated (in a single laboratory validation) a novel and improved primer/probe combination (Mit1C) for real-time PCR detection of C. cayetanensis in produce. The newly developed primer/probe combination targets a conserved region of the mitochondrial genome of C. cayetanensis that varies in other closely related organisms. The primer/probe combination was evaluated both in silico and using several real-time PCR kits and polymerases against an inclusivity/exclusivity panel comprised of a variety of C. cayetanensis oocysts, as well as DNA from other related Cyclospora spp. and closely related parasites. The new primer/probe combination amplified only C. cayetanensis, thus demonstrating specificity. Sensitivity was evaluated by artificially contaminating cilantro, raspberries, and romaine lettuce with variable numbers (200 and 5) of C. cayetanensis oocysts. As few as 5 oocysts were detected in 75%, 67.7%, and 50% of the spiked produce samples (cilantro, raspberries, and romaine lettuce), respectively, all uninoculated samples and no-template real-time PCR controls were negative. The improved primer/probe combination should prove an effective analytical tool for the specific detection of C. cayetanensis in produce.

Restriction Endonuclease-Based Modification-Dependent Enrichment (REMoDE) of DNA for Metagenomic Sequencing.

Metagenomic sequencing is a swift and powerful tool to ascertain the presence of an organism of interest in a sample. However, sequencing coverage of the organism of interest can be insufficient due to an inundation of reads from irrelevant organisms in the sample. Here, we report a nuclease-based approach to rapidly enrich for DNA from certain organisms, including enterobacteria, based on their differential endogenous modification patterns. We exploit the ability of taxon-specific methylated motifs to resist the action of cognate methylation-sensitive restriction endonucleases that thereby digest unwanted, unmethylated DNA. Subsequently, we use a distributive exonuclease or electrophoretic separation to deplete or exclude the digested fragments, thus enriching for undigested DNA from the organism of interest. As a proof of concept, we apply this method to enrich for the enterobacteria Escherichia coli and Salmonella enterica by 11- to 142-fold from mock metagenomic samples and validate this approach as a versatile means to enrich for genomes of interest in metagenomic samples. IMPORTANCE Pathogens that contaminate the food supply or spread through other means can cause outbreaks that bring devastating repercussions to the health of a populace. Investigations to trace the source of these outbreaks are initiated rapidly but can be drawn out due to the labored methods of pathogen isolation. Metagenomic sequencing can alleviate this hurdle but is often insufficiently sensitive. The approach and implementations detailed here provide a rapid means to enrich for many pathogens involved in foodborne outbreaks, thereby improving the utility of metagenomic sequencing as a tool in outbreak investigations. Additionally, this approach provides a means to broadly enrich for otherwise minute levels of modified DNA, which may escape unnoticed in metagenomic samples.

Seasonality, shelf life and storage atmosphere are main drivers of the microbiome and E. coli O157:H7 colonization of post-harvest lettuce cultivated in a major production area in California.

BACKGROUND: Lettuce is linked to recurrent outbreaks of Shiga toxin-producing Escherichia coli (STEC) infections, the seasonality of which remains unresolved. Infections have occurred largely from processed lettuce, which undergoes substantial physiological changes during storage. We investigated the microbiome and STEC O157:H7 (EcO157) colonization of fresh-cut lettuce of two cultivars with long and short shelf life harvested in the spring and fall in California and stored in modified atmosphere packaging (MAP) at cold and warm temperatures. RESULTS: Inoculated EcO157 declined significantly less on the cold-stored cultivar with short shelf life, while multiplying rapidly at 24 °C independently of cultivar. Metagenomic sequencing of the lettuce microbiome revealed that the pre-storage bacterial community was variable but dominated by species in the Erwiniaceae and Pseudomonadaceae. After cold storage, the microbiome composition differed between cultivars, with a greater relative abundance (RA) of Erwiniaceae and Yersiniaceae on the cultivar with short shelf life. Storage at 24 °C shifted the microbiome to higher RAs of Erwiniaceae and Enterobacteriaceae and lower RA of Pseudomonadaceae compared with 6 °C. Fall harvest followed by lettuce deterioration were identified by recursive partitioning as important factors associated with high EcO157 survival at 6 °C, whereas elevated package CO2 levels correlated with high EcO157 multiplication at 24 °C. EcO157 population change correlated with the lettuce microbiome during 6 °C storage, with fall microbiomes supporting the greatest EcO157 survival on both cultivars. Fall and spring microbiomes differed before and during storage at both temperatures. High representation of Pantoea agglomerans was a predictor of fall microbiomes, lettuce deterioration, and enhanced EcO157 survival at 6 °C. In contrast, higher RAs of Erwinia persicina, Rahnella aquatilis, and Serratia liquefaciens were biomarkers of spring microbiomes and lower EcO157 survival. CONCLUSIONS: The microbiome of processed MAP lettuce evolves extensively during storage. Under temperature abuse, high CO2 promotes a lettuce microbiome enriched in taxa with anaerobic capability and EcO157 multiplication. In cold storage, our results strongly support a role for season and lettuce deterioration in EcO157 survival and microbiome composition, suggesting that the physiology and microbiomes of fall- and spring-harvested lettuce may contribute to the seasonality of STEC outbreaks associated with lettuce grown in coastal California.

Comparative Transcriptomics of Shiga Toxin-Producing and Commensal Escherichia coli and Cytokine Responses in Colonic Epithelial Cell Culture Infections.

Ingestion of Shiga toxin-producing Escherichia coli (STEC) can result in a range of illness severity from asymptomatic to hemorrhagic colitis and death; thus risk assessment of STEC strains for human pathogenicity is important in the area of food safety. Illness severity depends in part on the combination of virulence genes carried in the genome, which can vary between strains even of identical serotype. To better understand how core genes are regulated differently among strains and to identify possible novel STEC virulence gene candidates that could be added to the risk assessment repertoire, we used comparative transcriptomics to investigate global gene expression differences between two STEC strains associated with severe illness and a commensal E. coli strain during in vitro intestinal epithelial cell (IEC) infections. Additionally, we compared a wide array of concomitant cytokine levels produced by the IECs. The cytokine expression levels were examined for a pattern representing STEC pathogenicity; however, while one STEC strain appeared to elicit a proinflammatory response, infection by the other strain produced a pattern comparable to the commensal E. coli. This result may be explained by the significant differences in gene content and expression observed between the STEC strains. RNA-Seq analysis revealed considerable disparity in expression of genes in the arginine and tryptophan biosynthesis/import pathways between the STEC strains and the commensal E. coli strain, highlighting the important role some amino acids play in STEC colonization and survival. Contrasting differential expression patterns were observed for genes involved in respiration among the three strains suggesting that metabolic diversity is a strategy utilized to compete with resident microflora for successful colonization. Similar temporal expression results for known and putative virulence genes were observed in the STEC strains, revealing strategies used for survival prior to and after initial adherence to IECs. Additionally, three genes encoding hypothetical proteins located in mobile genetic elements were, after interrogation of a large set of E. coli genomes, determined to likely represent novel STEC virulence factors.

Draft Genome Sequences of Escherichia albertii, Escherichia fergusonii, and Strains Belonging to Six Cryptic Lineages of Escherichia spp.

We report here the genome sequences of 55 strains belonging to the genus Escherichia from multiple animal and environmental sources. These strains include representatives of Escherichia albertii, Escherichia fergusonii, and six additional genetically distinct lineages of Escherichia spp., one of which is newly discovered and is being reported for the first time here.

Species-Wide Collection of Escherichia coli Isolates for Examination of Genomic Diversity.

Pathogenic and nonpathogenic Escherichia coli strains present a vast genomic diversity. We report the genome sequences of 2,244 E. coli isolates from multiple animal and environmental sources. Their phylogenetic relationships and potential risk to human health were examined.

Strain-Level Discrimination of Shiga Toxin-Producing Escherichia coli in Spinach Using Metagenomic Sequencing.

Consumption of fresh bagged spinach contaminated with Shiga toxin-producing Escherichia coli (STEC) has led to severe illness and death; however current culture-based methods to detect foodborne STEC are time consuming. Since not all STEC strains are considered pathogenic to humans, it is crucial to incorporate virulence characterization of STEC in the detection method. In this study, we assess the comprehensiveness of utilizing a shotgun metagenomics approach for detection and strain-level identification by spiking spinach with a variety of genomically disparate STEC strains at a low contamination level of 0.1 CFU/g. Molecular serotyping, virulence gene characterization, microbial community analysis, and E. coli core gene single nucleotide polymorphism (SNP) analysis were performed on metagenomic sequence data from enriched samples. It was determined from bacterial community analysis that E. coli, which was classified at the phylogroup level, was a major component of the population in most samples. However, in over half the samples, molecular serotyping revealed the presence of indigenous E. coli which also contributed to the percent abundance of E. coli. Despite the presence of additional E. coli strains, the serotype and virulence genes of the spiked STEC, including correct Shiga toxin subtype, were detected in 94% of the samples with a total number of reads per sample averaging 2.4 million. Variation in STEC abundance and/or detection was observed in replicate spiked samples, indicating an effect from the indigenous microbiota during enrichment. SNP analysis of the metagenomic data correctly placed the spiked STEC in a phylogeny of related strains in cases where the indigenous E. coli did not predominate in the enriched sample. Also, for these samples, our analysis demonstrates that strain-level phylogenetic resolution is possible using shotgun metagenomic data for determining the genomic relatedness of a contaminating STEC strain to other closely related E. coli.

Towards next-generation sequencing analytics for foodborne RNA viruses: Examining the effect of RNA input quantity and viral RNA purity.

Detection and identification of viruses in food samples are technically challenging due largely to the low viral copy number in contaminated food items, and the lack of effective culture enrichment methods that are amenable to regulatory applications for many of the common foodborne viruses. Using an Illumina MiSeq platform and two hepatitis A virus (HAV) cell-culture adapted strains as a representative enteric virus species, this study examined the limits of single-stranded RNA (ssRNA) viral detection following next-generation sequencing without pre-amplification of the viral genome. Complete viral genome sequences were obtained from HAV samples of varying purities and with an input as low as 2ng total RNA containing 1.4×10(5) copies of viral RNA. In addition, single nucleotide variations were reproducibly detected over the range of concentrations examined, and their identity confirmed by alternate sequencing technology. In summary, next-generation sequencing technology has the potential for sensitive detection/identification of a viral genome at a low copy number. This study provides a benchmark for metagenomic sequencing application as is required for virus detection in complex food matrices using a culture-independent diagnostic approach.

Investigating the Relatedness of Enteroinvasive Escherichia coli to Other E. coli and Shigella Isolates by Using Comparative Genomics.

Enteroinvasive Escherichia coli (EIEC) is a unique pathovar that has a pathogenic mechanism nearly indistinguishable from that of Shigella species. In contrast to isolates of the four Shigella species, which are widespread and can be frequent causes of human illness, EIEC causes far fewer reported illnesses each year. In this study, we analyzed the genome sequences of 20 EIEC isolates, including 14 first described in this study. Phylogenomic analysis of the EIEC genomes demonstrated that 17 of the isolates are present in three distinct lineages that contained only EIEC genomes, compared to reference genomes from each of the E. coli pathovars and Shigella species. Comparative genomic analysis identified genes that were unique to each of the three identified EIEC lineages. While many of the EIEC lineage-specific genes have unknown functions, those with predicted functions included a colicin and putative proteins involved in transcriptional regulation or carbohydrate metabolism. In silico detection of the Shigella virulence plasmid (pINV), which is essential for the invasion of host cells, demonstrated that a form of pINV was present in nearly all EIEC genomes, but the Mxi-Spa-Ipa region of the plasmid that encodes the invasion-associated proteins was absent from several of the EIEC isolates. The comparative genomic findings in this study support the hypothesis that multiple EIEC lineages have evolved independently from multiple distinct lineages of E. coli via the acquisition of the Shigella virulence plasmid and, in some cases, the Shigella pathogenicity islands.

Hybrid Shiga Toxin-Producing and Enterotoxigenic Escherichia sp. Cryptic Lineage 1 Strain 7v Harbors a Hybrid Plasmid.

UNLABELLED: Hybrid isolates of Shiga toxin-producing Escherichia coli (STEC) and enterotoxigenic E. coli (ETEC) encoding heat-stable enterotoxin (ST) are being reported with increasing frequency from a variety of sources. However, information regarding the plasmids that these strains harbor is scarce. In this study, we sequence and characterize a plasmid, p7v, from the STEC/ETEC hybrid strain 7v. Whole-genome phylogenetic analyses of STEC/ETEC hybrid strains and prototype E. coli isolates of other pathotypes placed 7v in the Escherichia sp. cryptic lineage 1 (CL1) clade. The complete plasmid, p7v, was determined to be 229,275 bp and encodes putative virulence factors that are typically carried on STEC plasmids as well as those often carried on ETEC plasmids, indicating that the hybrid nature of the strain extends beyond merely encoding the two toxins. Plasmid p7v carries two copies of sta with identical sequences, which were discovered to be divergent from the sta sequences found in the prototype human ETEC strains. Using a nomenclature scheme based on a phylogeny constructed from sta and stb sequences, the sta encoded on p7v is designated STa4. In silico analysis determined that p7v also encodes the K88 fimbria, a colonization factor usually associated with porcine ETEC plasmids. The p7v sequence and the presence of plasmid-encoded virulence factors are compared to those of other STEC/ETEC CL1 hybrid genomes and reveal gene acquisition/loss at the strain level. In addition, the interrogation of 24 STEC/ETEC hybrid genomes for identification of plasmid replicons, colonization factors, Stx and ST subtypes, and other plasmid-encoded virulence genes highlights the diversity of these hybrid strains. IMPORTANCE: Hybrid Shiga toxin-producing Escherichia coli/enterotoxigenic Escherichia coli (STEC/ETEC) strains, which have been isolated from environmental, animal, and human clinical samples, may represent an emerging threat as food-borne pathogens. Characterization of these strains is important for assessing virulence potential, aiding in the development of pathogen detection methods, and understanding how the hybrid strains evolve to potentially have a greater impact on public health. This study represents, to our knowledge, both the first characterization of a closed plasmid sequence from a STEC/ETEC hybrid strain and the most comprehensive phylogenetic analysis of available STEC/ETEC hybrid genomes to date. The results demonstrate how the mobility of plasmid-associated virulence genes has resulted in the creation of a diverse plasmid repertoire within the STEC/ETEC hybrid strains.

Stx-Producing Shigella Species From Patients in Haiti: An Emerging Pathogen With the Potential for Global Spread.

Shiga toxins (Stx) are commonly produced by Shigella dysenteriae serotype 1 and Stx-producing Escherichia coli. However, the toxin genes have been detected in additional Shigella species. We recently reported the emergence of Stx-producing Shigella in travelers in the United States and France who had recently visited Hispaniola (Haiti and the Dominican Republic). In this study, we confirm this epidemiological link by identifying Stx-producing Shigella from Haitian patients attending clinics near Port-au-Prince. We also demonstrate that the bacteriophage encoding Stx is capable of dissemination to stx-negative Shigella species found in Haiti, suggesting that Stx-producing Shigella may become more widespread within that region.

Application of metagenomic sequencing to food safety: detection of Shiga Toxin-producing Escherichia coli on fresh bagged spinach.

Culture-independent diagnostics reduce the reliance on traditional (and slower) culture-based methodologies. Here we capitalize on advances in next-generation sequencing (NGS) to apply this approach to food pathogen detection utilizing NGS as an analytical tool. In this study, spiking spinach with Shiga toxin-producing Escherichia coli (STEC) following an established FDA culture-based protocol was used in conjunction with shotgun metagenomic sequencing to determine the limits of detection, sensitivity, and specificity levels and to obtain information on the microbiology of the protocol. We show that an expected level of contamination (∼10 CFU/100 g) could be adequately detected (including key virulence determinants and strain-level specificity) within 8 h of enrichment at a sequencing depth of 10,000,000 reads. We also rationalize the relative benefit of static versus shaking culture conditions and the addition of selected antimicrobial agents, thereby validating the long-standing culture-based parameters behind such protocols. Moreover, the shotgun metagenomic approach was informative regarding the dynamics of microbial communities during the enrichment process, including initial surveys of the microbial loads associated with bagged spinach; the microbes found included key genera such as Pseudomonas, Pantoea, and Exiguobacterium. Collectively, our metagenomic study highlights and considers various parameters required for transitioning to such sequencing-based diagnostics for food safety and the potential to develop better enrichment processes in a high-throughput manner not previously possible. Future studies will investigate new species-specific DNA signature target regimens, rational design of medium components in concert with judicious use of additives, such as antibiotics, and alterations in the sample processing protocol to enhance detection.

Acquisition of the lac operon by Salmonella enterica.

BACKGROUND: Classical bacteriological characteristics of Salmonella enterica indicate that the members of this species are unable to utilize lactose as a carbon source. However, lactose-fermenting (Lac+) strains of several Salmonella serovars have been isolated from different foodborne outbreaks as well as different geographical regions worldwide. In the present study, we sequenced the genomes of 13 Lac + S. enterica isolates and characterized the lac region, comparing it to the lac region in other enteric bacterial species. RESULTS: Genetic analysis of the lac operons in the S. enterica genomes revealed that they all contain intact lacI, lacZ, and lacY genes. However, lacA was truncated in all of the S. enterica subsp. enterica isolates, encoding a 56 amino acid peptide rather than the full length 220 amino acid LacA protein. Molecular analyses of the 13 isolates revealed that the lac operon resided on a plasmid in some strains and in others was integrated into the bacterial chromosome. In most cases, an insertion sequence flanked at least one end of the operon. Interestingly, the S. enterica Montevideo and S. enterica Senftenberg isolates were found to harbor a plasmid with a high degree of sequence similarity to a plasmid from Klebsiella pneumoniae strain NK29 that also harbors the lac operon. In addition, two S. enterica Tennessee isolates carried two copies of the lac operon. Phylogenetic analysis based on lacIZY gene sequences determines distinct clusters, and reveals a greater correlation between lacIZY sequence and flanking organization than with either bacterial species or genomic location. CONCLUSIONS: Our results indicate that the lac region is highly mobile among Enterobacteriaceae and demonstrate that the Lac + S. enterica subsp. enterica serovars acquired the lac region through parallel events. The acquisition of the lac operon by several S. enterica serovars may be indicative of environmental adaptation by these bacteria.

Draft Genome Sequences of the Enteroinvasive Escherichia coli Strains M4163 and 4608-58.

We report here the draft genome sequences of enteroinvasive Escherichia coli (EIEC) O124:H30 strain M4163 isolated from imported French cheese and EIEC O143:H26 strain 4608-58. The assembled data determined that both strains contain multiple copies of the ipaH gene, as well as the pINV A form of the invasion plasmid.

Draft Genome Sequence of the Multidrug-Resistant Escherichia coli Strain LR09, Isolated from a Wastewater Treatment Plant.

We report the draft genome sequence of Escherichia coli O1:H6 strain LR09, which was isolated from a wastewater treatment plant and displays high resistance to five fluoroquinolone antimicrobials. The assembled data determine that the strain clusters with E. coli phylogroup F and harbors a plasmid conferring resistance to a broad spectrum of antibiotics.