Characterization and Genetic Diversity of Bacillus cereus Strains Isolated from Baby Wipes.

Bacillus cereus, a ubiquitous environmental microorganism known to cause foodborne illness, was isolated from samples taken from imported baby wipes from two different countries. These strains were characterized using a comprehensive molecular approach involving endpoint PCR, whole genome sequencing (WGS), comparative genomics, and biochemical analyses. A multiplex endpoint PCR assay was used to identify the enterotoxins: hemolysin BL, nonhemolytic enterotoxin, cytotoxin K, and enterotoxin FM toxin genes. Phylogenetically, the strains clustered into two major groups according to sequence type (ST) and singleton. We used the Center for Food Safety and Applied Nutrition (CFSAN) GalaxyTrakr BTyper computational tool to characterize the strains further. As an additional means of characterization, we investigated the possible role of carbohydrate transport systems and their role in nutrient uptake by performing a BLAST analysis of the 40 B. cereus genomes recovered from baby wipes. This study outlines a multifaceted workflow that uses the analysis of enterotoxigenic potential, bioinformatics, genomic diversity, genotype, phenotype, and carbohydrate utilization as a comprehensive strategy to characterize these B. cereus strains isolated from baby wipes and further our understanding of the phylogenetic relatedness of strains associated with baby wipe production facilities that could potentially pose an infection risk to a vulnerable infant population.

Characterization of Cronobacter sakazakii Strains Originating from Plant-Origin Foods Using Comparative Genomic Analyses and Zebrafish Infectivity Studies.

Cronobacter sakazakii continues to be isolated from ready-to-eat fresh and frozen produce, flours, dairy powders, cereals, nuts, and spices, in addition to the conventional sources of powdered infant formulae (PIF) and PIF production environments. To understand the sequence diversity, phylogenetic relationship, and virulence of C. sakazakii originating from plant-origin foods, comparative molecular and genomic analyses, and zebrafish infection (ZI) studies were applied to 88 strains. Whole genome sequences of the strains were generated for detailed bioinformatic analysis. PCR analysis showed that all strains possessed a pESA3-like virulence plasmid similar to reference C. sakazakii clinical strain BAA-894. Core genome analysis confirmed a shared genomic backbone with other C. sakazakii strains from food, clinical and environmental strains. Emerging nucleotide diversity in these plant-origin strains was highlighted using single nucleotide polymorphic alleles in 2000 core genes. DNA hybridization analyses using a pan-genomic microarray showed that these strains clustered according to sequence types (STs) identified by multi-locus sequence typing (MLST). PHASTER analysis identified 185 intact prophage gene clusters encompassing 22 different prophages, including three intact Cronobacter prophages: ENT47670, ENT39118, and phiES15. AMRFinderPlus analysis identified the CSA family class C ß-lactamase gene in all strains and a plasmid-borne mcr-9.1 gene was identified in three strains. ZI studies showed that some plant-origin C. sakazakii display virulence comparable to clinical strains. Finding virulent plant-origin C. sakazakii possessing significant genomic features of clinically relevant STs suggests that these foods can serve as potential transmission vehicles and supports widening the scope of continued surveillance for this important foodborne pathogen.

Analysis of the Molecular Diversity Among Cronobacter Species Isolated From Filth Flies Using Targeted PCR, Pan Genomic DNA Microarray, and Whole Genome Sequencing Analyses.

Cronobacter species are opportunistic pathogens capable of causing life-threatening infections in humans, with serious complications arising in neonates, infants, immuno-compromised individuals, and elderly adults. The genus is comprised of seven species: Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite a multiplicity of genomic data for the genus, little is known about likely transmission vectors. Using DNA microarray analysis, in parallel with whole genome sequencing, and targeted PCR analyses, the total gene content of two C. malonaticus, three C. turicensis, and 14 C. sakazaki isolated from various filth flies was assessed. Phylogenetic relatedness among these and other strains obtained during surveillance and outbreak investigations were comparatively assessed. Specifically, microarray analysis (MA) demonstrated its utility to cluster strains according to species-specific and sequence type (ST) phylogenetic relatedness, and that the fly strains clustered among strains obtained from clinical, food and environmental sources from United States, Europe, and Southeast Asia. This combinatorial approach was useful in data mining for virulence factor genes, and phage genes and gene clusters. In addition, results of plasmidotyping were in agreement with the species identity for each strain as determined by species-specific PCR assays, MA, and whole genome sequencing. Microarray and BLAST analyses of Cronobacter fly sequence datasets were corroborative and showed that the presence and absence of virulence factors followed species and ST evolutionary lines even though such genes were orthologous. Additionally, zebrafish infectivity studies showed that these pathotypes were as virulent to zebrafish embryos as other clinical strains. In summary, these findings support a striking phylogeny amongst fly, clinical, and surveillance strains isolated during 2010-2015, suggesting that flies are capable vectors for transmission of virulent Cronobacter spp.; they continue to circulate among United States and European populations, environments, and that this "pattern of circulation" has continued over decades.

Prevalence, Distribution, and Phylogeny of Type Two Toxin-Antitoxin Genes Possessed by Cronobacter Species where C. sakazakii Homologs Follow Sequence Type Lineages.

Cronobacter species are a group of foodborne pathogenic bacteria that cause both intestinal and systemic human disease in individuals of all age groups. Little is known about the mechanisms that Cronobacter employ to survive and persist in foods and other environments. Toxin-antitoxin (TA) genes are thought to play a role in bacterial stress physiology, as well as in the stabilization of horizontally-acquired re-combinatorial elements such as plasmids, phage, and transposons. TA systems have been implicated in the formation of a persistence phenotype in some bacterial species including Escherichia coli and Salmonella. This project's goal was to understand the phylogenetic relatedness among TA genes present in Cronobacter. Preliminary studies showed that two typical toxin genes, fic and hipA followed species evolutionary lines. A local database of 22 TA homologs was created for Cronobacter sakazakii and a Python version 3 shell script was generated to extract TA FASTA sequences present in 234 C. sakazakii genomes previously sequenced as part of Center for Food Safety and Applied Nutrition's (CFSAN) GenomeTrakr project. BLAST analysis showed that not every C. sakazakii strain possessed all twenty-two TA loci. Interestingly, some strains contained either a toxin or an antitoxin component, but not both. Five common toxin genes: ESA_00258 (parDE toxin-antitoxin family), ESA_00804 (relBE family), ESA_01887 (relBE family), ESA_03838 (relBE family), and ESA_04273 (YhfG-Fic family) were selected for PCR analysis and the primers were designed to detect these genes. PCR analysis showed that 55 of 63 strains possessed three of these genes Sequence analysis identified homologs of the target genes and some of the strains were PCR-negative for one or more of the genes, pointing to potential nucleotide polymorphisms in those loci or that these toxin genes were absent. Phylogenetic studies using a Cronobacter pan genomic microarray showed that for the most part TAs follow species evolutionary lines except for a few toxin genes possessed by some C. malonaticus and C. universalis strains; this demonstrates that some TA orthologues share a common phylogeny. Within the C. sakazakii strains, the prevalence and distribution of these TA homologs by C. sakazakii strain BAA-894 (a powdered infant formula isolate) followed sequence-type evolutionary lineages. Understanding the phylogeny of TAs among the Cronobacter species is essential to design future studies to realize the physiological mechanisms and roles for TAs in stress adaptation and persistence of Cronobacter within food matrices and food processing environments.

Draft genomes of Cronobacter sakazakii strains isolated from dried spices bring unique insights into the diversity of plant-associated strains.

Cronobacter sakazakii is a Gram-negative opportunistic pathogen that causes life- threatening infantile infections, such as meningitis, septicemia, and necrotizing enterocolitis, as well as pneumonia, septicemia, and urinary tract and wound infections in adults. Here, we report 26 draft genome sequences of C. sakazakii, which were obtained from dried spices from the USA, the Middle East, China, and the Republic of Korea. The average genome size of the C. sakazakii genomes was 4393 kb, with an average of 4055 protein coding genes, and an average genome G + C content of 56.9%. The genomes contained genes related to carbohydrate transport and metabolism, amino acid transport and metabolism, and cell wall/membrane biogenesis. In addition, we identified genes encoding proteins involved in osmotic responses such as DnaJ, Aquaproin Z, ProQ, and TreF, as well as virulence-related and heat shock-related proteins. Interestingly, a metabolic island comprised of a variably-sized xylose utilization operon was found within the spice-associated C. sakazakii genomes, which supports the hypothesis that plants may serve as transmission vectors or alternative hosts for Cronobacter species. The presence of the genes identified in this study can support the remarkable phenotypic traits of C. sakazakii such as the organism's capabilities of adaptation and survival in response to adverse growth environmental conditions (e.g. osmotic and desiccative stresses). Accordingly, the genome analyses provided insights into many aspects of physiology and evolutionary history of this important foodborne pathogen.

Analysis of enterotoxigenic Bacillus cereus strains from dried foods using whole genome sequencing, multi-locus sequence analysis and toxin gene prevalence and distribution using endpoint PCR analysis.

Bacillus cereus strains were isolated from dried foods, which included international brands of spices from South East Asia, Mexico and India purchased from several retail stores, samples of powdered infant formula (PIF), medicated fish feed and dietary supplements. The genetic diversity of 64 strains from spices and PIF was determined using a multiplex endpoint PCR assay designed to identify hemolysin BL, nonhemolytic enterotoxin, cytotoxin K, and enterotoxin FM toxin genes. Thirteen different B. cereus toxigenic gene patterns or profiles were identified among the strains. Randomly selected B. cereus strains were sequenced and compared with reference Genomic Groups from National Center Biotechnology Information using bioinformatics tools. A comprehensive multi-loci sequence analysis (MLSA) was designed using alleles from 25 known MLST genes specifically tailored for use with whole genome assemblies. A cohort of representative genomes of strains from a few FDA regulated commodities like dry foods and medicated fish feed was used to demonstrate the utility of the 25-MLSA approach for rapid clustering and identification of Genome Groups. The analysis clustered the strains from medicated fish feed, dry foods, and dietary supplements into phylogenetically-related groups. 25-MLSA also pointed to a greater diversity of B. cereus strains from foods and feed than previously recognized. Our integrated approach of toxin gene PCR, and to our knowledge, whole genome sequencing (WGS) based sequence analysis, may be the first of its kind that demonstrates enterotoxigenic potential and genomic diversity in parallel.

Genomic characterization of malonate positive Cronobacter sakazakii serotype O:2, sequence type 64 strains, isolated from clinical, food, and environment samples.

BACKGROUND: Malonate utilization, an important differential trait, well recognized as being possessed by six of the seven Cronobacter species is thought to be largely absent in Cronobacter sakazakii (Csak). The current study provides experimental evidence that confirms the presence of a malonate utilization operon in 24 strains of sequence type (ST) 64, obtained from Europe, Middle East, China, and USA; it offers explanations regarding the genomic diversity and phylogenetic relatedness among these strains, and that of other C. sakazakii strains. RESULTS: In this study, the presence of a malonate utilization operon in these strains was initially identified by DNA microarray analysis (MA) out of a pool of 347 strains obtained from various surveillance studies involving clinical, spices, milk powder sources and powdered infant formula production facilities in Ireland and Germany, and dried dairy powder manufacturing facilities in the USA. All ST64 C. sakazakii strains tested could utilize malonate. Zebrafish embryo infection studies showed that C. sakazakii ST64 strains are as virulent as other Cronobacter species. Parallel whole genome sequencing (WGS) and MA showed that the strains phylogenetically grouped as a separate clade among the Csak species cluster. Additionally, these strains possessed the Csak O:2 serotype. The nine-gene, ~ 7.7 kbp malonate utilization operon was located in these strains between two conserved flanking genes, gyrB and katG. Plasmidotyping results showed that these strains possessed the virulence plasmid pESA3, but in contrast to the USA ST64 Csak strains, ST64 Csak strains isolated from sources in Europe and the Middle East, did not possess the type six secretion system effector vgrG gene. CONCLUSIONS: Until this investigation, the presence of malonate-positive Csak strains, which are associated with foods and clinical cases, was under appreciated. If this trait was used solely to identify Cronobacter strains, many strains would likely be misidentified. Parallel WGS and MA were useful in characterizing the total genome content of these Csak O:2, ST64, malonate-positive strains and further provides an understanding of their phylogenetic relatedness among other virulent C. sakazakii strains.

Draft Genome Sequence of Cronobacter sakazakii GP1999, Sequence Type 145, an Epiphytic Isolate Obtained from the Tomato's Rhizoplane/Rhizosphere Continuum.

We present here the draft genome of Cronobacter sakazakii GP1999, a sequence type 145 strain isolated from the rhizosphere of tomato plants. Assembly and annotation of the genome resulted in a genome of 4,504,670 bp in size, with 4,148 coding sequences, and a GC content of 56.8%.

Comparative Genomic Characterization of the Highly Persistent and Potentially Virulent Cronobacter sakazakii ST83, CC65 Strain H322 and Other ST83 Strains.

Cronobacter (C.) sakazakii is an opportunistic pathogen and has been associated with serious infections with high mortality rates predominantly in pre-term, low-birth weight and/or immune compromised neonates and infants. Infections have been epidemiologically linked to consumption of intrinsically and extrinsically contaminated lots of reconstituted powdered infant formula (PIF), thus contamination of such products is a challenging task for the PIF producing industry. We present the draft genome of C. sakazakii H322, a highly persistent sequence type (ST) 83, clonal complex (CC) 65, serotype O:7 strain obtained from a batch of non-released contaminated PIF product. The presence of this strain in the production environment was traced back more than 4 years. Whole genome sequencing (WGS) of this strain together with four more ST83 strains (PIF production environment-associated) confirmed a high degree of sequence homology among four of the five strains. Phylogenetic analysis using microarray (MA) and WGS data showed that the ST83 strains were highly phylogenetically related and MA showed that between 5 and 38 genes differed from one another in these strains. All strains possessed the pESA3-like virulence plasmid and one strain possessed a pESA2-like plasmid. In addition, a pCS1-like plasmid was also found. In order to assess the potential in vivo pathogenicity of the ST83 strains, each strain was subjected to infection studies using the recently developed zebrafish embryo model. Our results showed a high (90-100%) zebrafish mortality rate for all of these strains, suggesting a high risk for infections and illness in neonates potentially exposed to PIF contaminated with ST83 C. sakazakii strains. In summary, virulent ST83, CC65, serotype CsakO:7 strains, though rarely found intrinsically in PIF, can persist within a PIF manufacturing facility for years and potentially pose significant quality assurance challenges to the PIF manufacturing industry.

Draft Genome Sequences of Enterotoxigenic Bacillus cereus Strains Obtained from Powdered Infant Formula.

We introduce the draft genome sequences of five enterotoxigenic Bacillus cereus strains: Bc 12, Bc 67, Bc 111, Bc 112, and Bc 113, which were obtained from powdered infant formula. The genome sizes of the strains ranged from 5.5 to 5.8 Mb, and the G+C contents were ~35.2%.

Analysis and Characterization of Proteins Associated with Outer Membrane Vesicles Secreted by Cronobacter spp.

Little is known about secretion of outer membrane vesicles (OMVs) by Cronobacter. In this study, OMVs isolated from Cronobacter sakazakii, Cronobacter turicensis, and Cronobacter malonaticus were examined by electron microscopy (EM) and their associated outer membrane proteins (OMP) and genes were analyzed by SDS-PAGE, protein sequencing, BLAST, PCR, and DNA microarray. EM of stained cells revealed that the OMVs are secreted as pleomorphic micro-vesicles which cascade from the cell's surface. SDS-PAGE analysis identified protein bands with molecular weights of 18 kDa to >100 kDa which had homologies to OMPs such as GroEL; OmpA, C, E, F, and X; MipA proteins; conjugative plasmid transfer protein; and an outer membrane auto-transporter protein (OMATP). PCR analyses showed that most of the OMP genes were present in all seven Cronobacter species while a few genes (OMATP gene, groEL, ompC, mipA, ctp, and ompX) were absent in some phylogenetically-related species. Microarray analysis demonstrated sequence divergence among the OMP genes that was not captured by PCR. These results support previous findings that OmpA and OmpX may be involved in virulence of Cronobacter, and are packaged within secreted OMVs. These results also suggest that other OMV-packaged OMPs may be involved in roles such as stress response, cell wall and plasmid maintenance, and extracellular transport.

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.

Genome Sequence of Cronobacter sakazakii Serogroup O:4, Sequence Type 4 Strain CDC 2009-03746, Isolated from a Fatal Case of Infantile Meningitis.

We report the draft genome sequence of a Cronobacter sakazakii serogroup O:4, sequence type 4 strain, CDC 2009-03746 (=NM1240=2009-06-01), isolated from a fatal case of infantile meningitis. The draft genome has a size of 4,492,904 bp and a G+C% content of 56.7.

A proposed harmonized LPS molecular-subtyping scheme for Cronobacter species.

Cronobacter are opportunistic pathogens, which cause infections in all age groups. To aid the characterization of Cronobacter in foods and environments a harmonized LPS identification scheme for molecular serotyping is needed. To this end, we studied 409 Cronobacter isolates representing the seven Cronobacter species using two previously reported molecular serotyping schemes, described here as Mullane-Jarvis (M-J) and Sun schemes. PCR analysis revealed many overlapping results that were obtained when independently applying the two serotyping schemes. There were complete agreements between the two PCR schemes for Cronobacter sakazakii (Csak) O:1, Csak O:3, and Csak O:7 serotypes. However, only thirty-five of 41 Csak O:4 strains, identified using the M-J scheme, were PCR-positive with the Sun scheme primers. Also the Sun scheme Csak O:5 primers failed to identify this serotype in any of the C. sakazakii strains tested, but did recognize seven Cronobacter turicensis strains, which were identified as Ctur O:3 using the M-J scheme. Similarly, the Sun scheme Csak O:6 primers recognized 30 Cronobacter malonaticus O:2 strains identified with the M-J scheme, but failed to identify this serotype in any C. sakazakii strain investigated. In this report, these findings are summarized and a harmonized molecular-serotyping scheme is proposed which is predicated on the correct identification of Cronobacter species, prior to serotype determination. In summary, fourteen serotypes were identified using the combined protocol, which consists of Csak O:1-O:4, and Csak O:7; Cmal O:1-O:2; Cdub O:1-O:2, Cmuy O:1-O:2, Cuni O:1, as well as Ctur O:1 and Ctur O:3.

Re-examination of the taxonomic status of Enterobacter helveticus, Enterobacter pulveris and Enterobacter turicensis as members of the genus Cronobacter and their reclassification in the genera Franconibacter gen. nov. and Siccibacter gen. nov. as Franconibacter helveticus comb. nov., Franconibacter pulveris comb. nov. and Siccibacter turicensis comb. nov., respectively.

Recently, a taxonomical re-evaluation of the genus Enterobacter, based on multi-locus sequence typing (MLST) analysis, has led to the proposal that the species Enterobacter pulveris, Enterobacter helveticus and Enterobacter turicensis should be reclassified as novel species of the genus Cronobacter. In the present work, new genome-scale analyses, including average nucleotide identity, genome-scale phylogeny and k-mer analysis, coupled with previously reported DNA-DNA hybridization values and biochemical characterization strongly indicate that these three species of the genus Enterobacter are not members of the genus Cronobacter, nor do they belong to the re-evaluated genus Enterobacter. Furthermore, data from this polyphasic study indicated that all three species constitute two new genera. We propose reclassifying Enterobacter pulveris and Enterobacter helveticus in the genus Franconibacter gen. nov. as Franconibacter pulveris comb. nov. (type strain 601/05(T) = LMG 24057(T) = DSM 19144(T)) and Franconibacter helveticus comb. nov. (type strain 513/05(T) = LMG 23732(T) = DSM 18396(T)), respectively, and Enterobacter turicensis in the genus Siccibacter gen. nov. as Siccibacter turicensis comb. nov. (type strain 508/05(T) = LMG 23730(T) = DSM 18397(T)).

Genome Sequences of Two Enterobacter pulveris Strains, 601/05T (=LMG 24057T =DSM 19144T) and 1160/04 (=LMG 24058 =DSM 19146), Isolated from Fruit Powder.

We report the draft genome sequences of the Enterobacter pulveris strains 601/05(T) (=LMG24057(T) =DSM19144(T)) and 1160/04 (=LMG24058 =DSM19146), isolated from fruit powder. The genome assemblies for the E. pulveris type strain, LMG24057, and strain LMG24058 have sizes of 4,708,624 and 4,811,103 bp and G+C contents of 56.6% and 56.5%, respectively.

Genome Sequence of Enterobacter turicensis Strain 610/05 (LMG 23731), Isolated from Fruit Powder.

We report the draft genome sequence of Enterobacter turicensis strain 610/05 (LMG 23731), isolated from fruit powder. The draft genome has a size of 4,182,790 bp and a G+C% content of 58.0.

Genome Sequence of an Enterobacter helveticus Strain, 1159/04 (LMG 23733), Isolated from Fruit Powder.

We report the draft genome sequence of Enterobacter helveticus strain LMG 23733, isolated from fruit powder. The draft genome assembly for E. helveticus strain LMG 23733 has a size of 4,635,476 bp and a G+C content of 55.9%.