A Vegetable Fermentation Facility Hosts Distinct Microbiomes Reflecting the Production Environment.

Fermented vegetables are highly popular internationally in part due to their enhanced nutritional properties, cultural history, and desirable sensorial properties. In some instances, fermented foods provide a rich source of the beneficial microbial communities that could promote gastrointestinal health. The indigenous microbiota that colonize fermentation facilities may impact food quality, food safety, and spoilage risks and maintain the nutritive value of the product. Here, microbiomes within sauerkraut production facilities were profiled to characterize variance across surfaces and to determine the sources of these bacteria. Accordingly, we used high-throughput sequencing of the 16S rRNA gene in combination with whole-genome shotgun analyses to explore biogeographical patterns of microbial diversity and assembly within the production facility. Our results indicate that raw cabbage and vegetable handling surfaces exhibit more similar microbiomes relative to the fermentation room, processing area, and dry storage surfaces. We identified biomarker bacterial phyla and families that are likely to originate from the raw cabbage and vegetable handling surfaces. Raw cabbage was identified as the main source of bacteria to seed the facility, with human handling contributing a minor source of inoculation. Leuconostoc and Lactobacillaceae dominated all surfaces where spontaneous fermentation occurs, as these taxa are associated with the process. Wall, floor, ceiling, and barrel surfaces host unique microbial signatures. This study demonstrates that diverse bacterial communities are widely distributed within the production facility and that these communities assemble nonrandomly, depending on the surface type.IMPORTANCE Fermented vegetables play a major role in global food systems and are widely consumed by various global cultures. In this study, we investigated an industrial facility that produces spontaneous fermented sauerkraut without the aid of starter cultures. This provides a unique system to explore and track the origins of an "in-house" microbiome in an industrial environment. Raw vegetables and the surfaces on which they are handled were identified as the likely source of bacterial communities rather than human contamination. As fermented vegetables increase in popularity on a global scale, understanding their production environment may help maintain quality and safety goals.

Genome sequencing and comparative genomics of enterohemorrhagic Escherichia coli O145:H25 and O145:H28 reveal distinct evolutionary paths and marked variations in traits associated with virulence & colonization.

BACKGROUND: Enterohemorrhagic Escherichia coli (EHEC) O145 are among the top non-O157 serogroups associated with severe human disease worldwide. Two serotypes, O145:H25 and O145:H28 have been isolated from human patients but little information is available regarding the virulence repertoire, origin and evolutionary relatedness of O145:H25. Hence, we sequenced the complete genome of two O145:H25 strains associated with hemolytic uremic syndrome (HUS) and compared the genomes with those of previously sequenced O145:H28 and other EHEC strains. RESULTS: The genomes of the two O145:H25 strains were 5.3 Mbp in size; slightly smaller than those of O145:H28 and other EHEC strains. Both strains contained three nearly identical plasmids and several prophages and integrative elements, many of which differed significantly in size, gene content and organization as compared to those present in O145:H28 and other EHECs. Furthermore, notable variations were observed in several fimbrial gene cluster and intimin types possessed by O145:H25 and O145:H28 indicating potential adaptation to distinct areas of host colonization. Comparative genomics further revealed that O145:H25 are genetically more similar to other non-O157 EHEC strains than to O145:H28. CONCLUSION: Phylogenetic analysis accompanied by comparative genomics revealed that O145:H25 and O145:H28 evolved from two separate clonal lineages and that horizontal gene transfer and gene loss played a major role in the divergence of these EHEC serotypes. The data provide further evidence that ruminants might be a possible reservoir for O145:H25 but that they might be impaired in their ability to establish a persistent colonization as compared to other EHEC strains.

Pan-genome analysis of the emerging foodborne pathogen Cronobacter spp. suggests a species-level bidirectional divergence driven by niche adaptation.

BACKGROUND: Members of the genus Cronobacter are causes of rare but severe illness in neonates and preterm infants following the ingestion of contaminated infant formula. Seven species have been described and two of the species genomes were subsequently published. In this study, we performed comparative genomics on eight strains of Cronobacter, including six that we sequenced (representing six of the seven species) and two previously published, closed genomes. RESULTS: We identified and characterized the features associated with the core and pan genome of the genus Cronobacter in an attempt to understand the evolution of these bacteria and the genetic content of each species. We identified 84 genomic regions that are present in two or more Cronobacter genomes, along with 45 unique genomic regions. Many potentially horizontally transferred genes, such as lysogenic prophages, were also identified. Most notable among these were several type six secretion system gene clusters, transposons that carried tellurium, copper and/or silver resistance genes, and a novel integrative conjugative element. CONCLUSIONS: Cronobacter have diverged into two clusters, one consisting of C. dublinensis and C. muytjensii (Cdub-Cmuy) and the other comprised of C. sakazakii, C. malonaticus, C. universalis, and C. turicensis, (Csak-Cmal-Cuni-Ctur) from the most recent common ancestral species. While several genetic determinants for plant-association and human virulence could be found in the core genome of Cronobacter, the four Cdub-Cmuy clade genomes contained several accessory genomic regions important for survival in a plant-associated environmental niche, while the Csak-Cmal-Cuni-Ctur clade genomes harbored numerous virulence-related genetic traits.

Whole-genome sequencing of gentamicin-resistant Campylobacter coli isolated from U.S. retail meats reveals novel plasmid-mediated aminoglycoside resistance genes.

Aminoglycoside resistance in Campylobacter has been routinely monitored in the United States in clinical isolates since 1996 and in retail meats since 2002. Gentamicin resistance first appeared in a single human isolate of Campylobacter coli in 2000 and in a single chicken meat isolate in 2007, after which it increased rapidly to account for 11.3% of human isolates and 12.5% of retail isolates in 2010. Pulsed-field gel electrophoresis analysis indicated that gentamicin-resistant C. coli isolates from retail meat were clonal. We sequenced the genomes of two strains of this clone using a next-generation sequencing technique in order to investigate the genetic basis for the resistance. The gaps of one strain were closed using optical mapping and Sanger sequencing, and this is the first completed genome of C. coli. The two genomes are highly similar to each other. A self-transmissible plasmid carrying multiple antibiotic resistance genes was revealed within both genomes, carrying genes encoding resistance to gentamicin, kanamycin, streptomycin, streptothricin, and tetracycline. Bioinformatics analysis and experimental results showed that gentamicin resistance was due to a phosphotransferase gene, aph(2")-Ig, not described previously. The phylogenetic relationship of this newly emerged clone to other Campylobacter spp. was determined by whole-genome single nucleotide polymorphisms (SNPs), which showed that it clustered with the other poultry isolates and was separated from isolates from livestock.

Identification and characterization of five new molecular serogroups of Cronobacter spp.

Cronobacter spp. (formerly Enterobacter sakazakii) is an emerging foodborne pathogen consisting of seven species including C. sakazakii, C. malonaticus, C. muytjensii, C. turicensis, C. dublinensis (with three subspecies, dublinensis, lausannensis, and lactaridi), C. universalis, and C. condimenti. To date, 12 Cronobacter serogroups have been identified. In this study, MboII restriction fragment length polymorphism patterns and DNA sequences of O-antigen gene clusters were used to identify novel serogroups of Cronobacter spp. Sequence analysis of the O-antigen regions, located between galF and gnd, of strains with distinct restriction fragment length polymorphism patterns revealed five unique gene clusters. These new O-antigen gene clusters were species specific and were termed C. turicensis O3, C. muytjensii O2, C. dublinensis O1, C. dublinensis O2, and C. universalis O1. Polymerase chain reaction assays were developed using primers specific to O-antigen processing genes and used to screen a collection of Cronobacter strains to determine the frequency of these newly identified serotypes.

Rapid genomic-scale analysis of Escherichia coli O104:H4 by using high-resolution alternative methods to next-generation sequencing.

Two technologies, involving DNA microarray and optical mapping, were used to quickly assess gene content and genomic architecture of recent emergent Escherichia coli O104:H4 and related strains. In real-time outbreak investigations, these technologies can provide congruent perspectives on strain, serotype, and pathotype relationships. Our data demonstrated clear discrimination between clinically, temporally, and geographically distinct O104:H4 isolates and rapid characterization of strain differences.

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.

Complete Genome Sequences of Four Enterohemolysin-Positive (ehxA) Enterocyte Effacement-Negative Shiga Toxin-Producing Escherichia coli Strains.

Shiga toxin-producing Escherichia coli (STEC) strains are important foodborne pathogens associated with human disease. Most disease-associated STEC strains carry the locus of enterocyte effacement (LEE); however, regularly LEE-negative STEC strains are recovered from ill patients. Few reference sequences are available for these isolate types. Here, we report here the complete genome sequences for four LEE-negative STEC strains.

Genomic variability among enteric pathogens: the case of the mutS-rpoS intergenic region.

The mutS-rpoS intergenic region of enteric bacteria ranges in size from 88 bp in Yersinia to > 12000 bp in Salmonella. We interpret this expansion as the result of the horizontal transfer of segments of DNA from diverse origins. Both comparative genomic analysis and selective sequencing of a variety of Escherichia coli pathogens have provided additional evidence for reassortment of segments within this region.

Development of a Custom-Designed, Pan Genomic DNA Microarray to Characterize Strain-Level Diversity among Cronobacter spp.

Cronobacter species cause infections in all age groups; however neonates are at highest risk and remain the most susceptible age group for life-threatening invasive disease. The genus contains seven species:Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite an abundance of published genomes of these species, genomics-based epidemiology of the genus is not well established. The gene content of a diverse group of 126 unique Cronobacter and taxonomically related isolates was determined using a pan genomic-based DNA microarray as a genotyping tool and as a means to identify outbreak isolates for food safety, environmental, and clinical surveillance purposes. The microarray constitutes 19,287 independent genes representing 15 Cronobacter genomes and 18 plasmids and 2,371 virulence factor genes of phylogenetically related Gram-negative bacteria. The Cronobacter microarray was able to distinguish the seven Cronobacter species from one another and from non-Cronobacter species; and within each species, strains grouped into distinct clusters based on their genomic diversity. These results also support the phylogenic divergence of the genus and clearly highlight the genomic diversity among each member of the genus. The current study establishes a powerful platform for further genomics research of this diverse genus, an important prerequisite toward the development of future countermeasures against this foodborne pathogen in the food safety and clinical arenas.

Optical mapping and 454 sequencing of Escherichia coli O157 : H7 isolates linked to the US 2006 spinach-associated outbreak.

Optical maps for five representative clinical, food-borne and bovine-derived isolates from the 2006 Escherichia coli O157 : H7 outbreak linked to fresh spinach in the United States showed a common set of 14 distinct chromosomal markers that define the outbreak strain. Partial 454 DNA sequencing was used to characterize the optically mapped chromosomal markers. The markers included insertions, deletions, substitutions and a simple single nucleotide polymorphism creating a BamHI site. The Shiga toxin gene profile of the spinach-associated outbreak isolates (stx1(-) stx2(+) stx2c(+)) correlated with prophage insertions different from those in the prototypical EDL933 and Sakai reference strains (stx1(+) stx2(+) stx2c(-)). The prophage occupying the yehV chromosomal position in the spinach-associated outbreak isolates was similar to the stx1(+) EDL933 cryptic prophage V, but it lacked the stx1 gene. In EDL933, the stx2 genes are within prophage BP933-W at the wrbA chromosomal locus; this locus was unoccupied in the spinach outbreak isolates. Instead, the stx2 genes were found within a chimeric BP933-W-like prophage with a different integrase, inserted at the argW locus in the outbreak isolates. An extra set of Shiga toxin genes, stx2c, was found in the outbreak isolates within a prophage integrated at the sbcB locus. The optical maps of two additional clinical isolates from the outbreak showed a single, different prophage variation in each, suggesting that changes occurred in the source strain during the course of this widespread, multi-state outbreak.

Optical maps distinguish individual strains of Escherichia coli O157 : H7.

Optical maps of 11 Escherichia coli O157 : H7 strains have been generated by the assembly of contiguous sets of restriction fragments across their entire 5.3 to 5.6 Mbp chromosomes. Each strain showed a distinct, highly individual configuration of 500-700 BamHI fragments, yielding a map resembling a DNA 'bar code'. The accuracy of optical mapping was assessed by comparing directly the in silico restriction maps of two wholly sequenced reference genomes of E. coli O157 : H7, i.e. EDL933 and the Sakai isolate (RIMD 0509952), with the optical maps of the same strains. The optical maps of nine other E. coli O157 : H7 strains were compared similarly, using the sequence-based maps of the Sakai and EDL933 strains as references. A total of 91 changes at 28 loci were positioned and sized; these included complex chromosomal inversions, insertions, deletions, substitutions, as well as a number of simple RFLPs. The optical maps defined unique genome landmarks in each of the strains and demonstrated the ability of optical mapping to distinguish and differentiate, at the individual level, strains of this important pathogen.

Detection of recombination among Salmonella enterica strains using the incongruence length difference test.

Particular serovars of Salmonella enterica have emerged as significant foodborne pathogens in humans. At the chromosomal level, discrete regions in the Salmonella genome have been identified that are known to play important roles in the maintenance, survival, and virulence of S. enterica within the host. Interestingly, several of these loci appear to have been acquired by horizontal transfer of DNA among and between bacterial species. The profound importance of recombination in pathogen emergence is just now being realized, perhaps explaining the sudden interest in developing novel and facile ways for detecting putative horizontal transfer events in bacteria. The incongruence length difference (ILD) test offers one such means. ILD uses phylogeny to trace sequences that may have been acquired promiscuously by exchange of DNA during chromosome evolution. We show here that the ILD test readily detects recombinations that have taken place in several housekeeping genes in Salmonella as well as genes composing the type 1 pilin complex (14 min) and the inv-spa invasion gene complex (63 min). Moreover, the ILD test indicated that the mutS gene (64 min), whose product helps protect the bacterial genome from invasion by foreign DNA, appears to have undergone intragenic recombination within S. enterica subspecies I. ILD findings were supported using additional tests known to be independent of the ILD approach (e.g., split decomposition analysis and compatibility of sites). Taken together, these data affirm the application of the ILD test as one approach for identifying recombined sequences in the Salmonella chromosome. Furthermore, horizontally acquired sequences within mutS support a model whereby evolutionarily important recombinants of S. enterica are rescued from strains carrying defective mutS alleles via horizontal transfer.

Evolution of multi-gene segments in the mutS-rpoS intergenic region of Salmonella enterica serovar Typhimurium LT2.

The nucleotide sequence of the 12.6 kb region between the mutS and rpoS genes of Salmonella enterica serovar Typhimurium LT2 (S. typhimurium) was compared to other enteric bacterial mutS-rpoS intergenic regions. The mutS-rpoS region is composed of three distinct segments, designated HK, O and S, as defined by sequence similarities to contiguous ORFs in other bacteria. Inverted chromosomal orientations of each of these segments are found between the mutS and rpoS genes in related ENTEROBACTERIACEAE: The HK segment is distantly related to a cluster of seven ORFs found in Haemophilus influenzae and a cluster of five ORFs found between the mutS and rpoS genes in Escherichia coli K-12. The O segment is related to the mutS-rpoS intergenic region found in E. coli O157:H7 and Shigella dysenteriae type 1. The third segment, S, is common to diverse Salmonella species, but is absent from E. coli. Despite the extensive collinearity and conservation of the overall genetic maps of S. typhimurium and E. coli K-12, the insertions, deletions and inversions in the mutS-rpoS region provide evidence that this region of the chromosome is an active site for horizontal gene transfer and rearrangement.