Multilocus Sequence Typing of Cronobacter Strains Isolated from Retail Foods and Environmental Samples

Microbiological survey and genomic analysis of Cronobacter sakazakii strains isolated from US households and retail foods

Cronobacter species are opportunistic pathogens that are capable of causing morbidity and mortality, particularly in infants. Although the transmission dynamics involved in Cronobacter infections remain largely unknown, contaminated powdered infant formula (PIF) has been linked to 30% of Cronobacter sakazakii cases involving invasive illness in infants. As several lines of evidence have implicated the domestic environment in PIF contamination, we undertook a microbiological survey of homes (N = 263) across the US. Cronobacter spp. and C. sakazakii were isolated from 36.1% and 24.7% of US homes, respectively, with higher recovery rates observed for floor and kitchen surfaces. Multi-locus sequence typing indicated that the dominant strain was C. sakazakii ST4, the sequence type most commonly associated with neonatal meningitis. For comparison purposes, retail foods (N = 4,009) were also surveyed, with the highest contamination frequencies (10.1%-26.3%) seen for nut products, seeds, and grains/baked goods/flours. The sequence type profile of isolates recovered from homes mirrored that of isolates recovered from retail foods, with increased representation of ST1, ST4, ST13, ST17, and ST40. Analysis of 386 whole genomic sequences revealed significant diversity. Redundancies were only observed for isolates recovered from within the same domicile, and there were no identical matches with sequences archived at the NCBI pathogen database. Genes coding for putative virulence and antibiotic resistance factors did not segregate with clinically significant sequence types. Collectively, these findings support the possibility that contamination events occurring within the home should not be overlooked as a contributor to community-onset Cronobacter infections.

IMPORTANCE

Cronobacter sakazakii is an opportunistic pathogen that can cause significant morbidity and mortality in neonates. Its transmission dynamics are poorly understood, though powered infant formula (PIF) is thought to be the major transmission vehicle. How the PIF becomes contaminated remains unknown. Our survey shows that roughly 1/4 of US homes are contaminated with Cronobacter sakazakii, particularly in the kitchen setting. Our analyses suggest that the domestic environment may contribute to contamination of PIF and provides insights into mitigating the risk of transmission.

Detection of Genus and Three Important Species of Cronobacter Using Novel Genus- and Species-Specific Genes Identified by Large-Scale Comparative Genomic Analysis

The genus Cronobacter includes seven species; however, the strains of Cronobacter sakazakii, Cronobacter malonaticus, and Cronobacter turicensis were highly correlated with clinical infections. Rapid and reliable identification of these three species of Cronobacter is important in monitoring and controlling diseases caused by these bacteria. Here, we identified four pairs of novel marker genes for the Cronobacter genus, C. sakazakii, C. malonaticus, and C. turicensis based on large-scale comparative genomic analysis from 799 Cronobacter and 136,146 non-Cronobacter genomes, including 10 Franconibacter and eight Siccibacter, which are close relatives of Cronobacter. Duplex and multiplex PCR methods were established based on these newly identified marker genes. The reliability of duplex and multiplex PCR methods was validated with 74 Cronobacter and 90 non-Cronobacter strains. Strains of C. sakazakii, C. malonaticus, and C. turicensis could be detected accurately at both the genus and species level. Moreover, the newly developed methods enable us to detect 2.5 × 10³ CFU/ml in pure culture. These data indicate that the accurate and sensitive established methods for Cronobacter can serve as valuable tools for the identification of these strains recovered from food, environmental, and clinical samples.

Role of the multiple efflux pump protein TolC on growth, morphology, and biofilm formation under nitric oxide stress in Cronobacter malonaticus

Nitric oxide (NO) is a biological signal molecule that can control and prevent the growth of most pathogens. Cronobacter species are a group of gram-negative foodborne pathogens that cause severe diseases, including neonatal meningitis, septicemia, and necrotizing enterocolitis, especially among newborns and infants consuming contaminated powdered infant formula. Cronobacter species might be tolerant to NO, resulting in severe infections. However, the specific mechanism of tolerance to NO in Cronobacter species is unclear. Here, we explore the effects of a key component, the protein TolC, of a multiple efflux pump on the growth, morphological changes, and biofilm formation of Cronobacter malonaticus under NO stress. We found that deletion of tolC resulted in a decreased growth rate under 100 mM sodium nitroprusside (NO donor) and led to more disruptive morphological injury to the bacterial cells. Furthermore, C. malonaticus lacking the TolC protein (ΔtolC mutant) showed weaker biofilm formation than the wild-type strain under normal or NO stress conditions. We have proved that TolC plays an important role in cell growth and biofilm formation of C. malonaticus. Therefore, our results may provide valuable theoretical basis for formulating clinical guidelines for treatment of disease caused by C. malonaticus and ensuring food safety.

Diversity of Cronobacter genus isolated between 1970 and 2019 on the American continent and genotyped using multi-locus sequence typing

This study aimed to evaluate the Cronobacter spp. strains isolated on the American continent and characterized using multi-locus sequence typing (MLST) available in the PubMLST database and current literature. From 465 Cronobacter spp. strains, the majority (n = 267, 57.4%) was from North America, mainly from USA (n = 234) and 198 (42.6%) were from South America, mainly from Brazil (n = 196). A total of 232 (49.9%) were isolated from foods, 102 (21.9%) from environmental, 87 (18.7%) from clinical, 27 (5.8%) from PIF, one from water (0.2%) and 16 (3.5%) from unknown sources. A total of five species were represented: Cronobacter sakazakii (374, 80.4%), Cronobacter malonaticus (41, 8.8%), Cronobacter dublinensis (29, 6.2%), Cronobacter turicensis (16, 3.5%) and Cronobacter muytjensii (5, 1.1%). The strains with complete MLST profile (n = 345) were assigned to 98 STs, a ratio of 3.5 strain by ST found and the calculated Simpson`s index was 0.93. The strains showed a high diversity and after eBURST analysis, 30 STs (n = 189) formed 12 single and/or double-locus variant clonal complexes (CC). A total of 38 STs (38.7%) were associated with clinical cases of infection, including well established C. sakazakii CC 1, 4, 8 and 83; C. malonaticus ST60, 307, 394 and 440; and C. sakazakii ST 12 and 494.

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.

Pan-genome diversification and recombination in Cronobacter sakazakii, an opportunistic pathogen in neonates, and insights to its xerotolerant lifestyle

Background

Cronobacter sakazakii is an emerging opportunistic bacterial pathogen known to cause neonatal and pediatric infections, including meningitis, necrotizing enterocolitis, and bacteremia. Multiple disease outbreaks of C. sakazakii have been documented in the past few decades, yet little is known of its genomic diversity, adaptation, and evolution. Here, we analyzed the pan-genome characteristics and phylogenetic relationships of 237 genomes of C. sakazakii and 48 genomes of related Cronobacter species isolated from diverse sources.

Results

The C. sakazakii pan-genome contains 17,158 orthologous gene clusters, and approximately 19.5% of these constitute the core genome. Phylogenetic analyses reveal the presence of at least ten deep branching monophyletic lineages indicative of ancestral diversification. We detected enrichment of functions involved in proton transport and rotational mechanism in accessory genes exclusively found in human-derived strains. In environment-exclusive accessory genes, we detected enrichment for those involved in tryptophan biosynthesis and indole metabolism. However, we did not find significantly enriched gene functions for those genes exclusively found in food strains. The most frequently detected virulence genes are those that encode proteins associated with chemotaxis, enterobactin synthesis, ferrienterobactin transporter, type VI secretion system, galactose metabolism, and mannose metabolism. The genes fos which encodes resistance against fosfomycin, a broad-spectrum cell wall synthesis inhibitor, and mdf(A) which encodes a multidrug efflux transporter were found in nearly all genomes. We found that a total of 2991 genes in the pan-genome have had a history of recombination. Many of the most frequently recombined genes are associated with nutrient acquisition, metabolism and toxin production.

Conclusions

Overall, our results indicate that the presence of a large accessory gene pool, ability to switch between ecological niches, a diverse suite of antibiotic resistance, virulence and niche-specific genes, and frequent recombination partly explain the remarkable adaptability of C. sakazakii within and outside the human host. These findings provide critical insights that can help define the development of effective disease surveillance and control strategies for Cronobacter-related diseases.

The occurrence and distribution characteristics of Cronobacter in diverse cereal kernels, flour, and flour-based products

Cronobacter was positive in cereals at a relatively high rate. In the present study, we investigated the occurrence and characteristics of this pathogen systematically in diverse cereals. All sampled food (N = 467) contained Cronobacter with a high positive rate of 54.0%. The enumeration experiment showed the concentration ranged from 0.3 to more than 110 MPN/100 g, and 87.9% of 127 samples were less than 10 MPN/100 g. There were significant differences (p < 0.05) in positive rates for Cronobacter among cereal kernels (40.2%), cereal flour (66.7%), cereal products made from raw cereal flour (87.6%), and cereal products made from flour (ready-to-eat) (17.4%). The dominant Cronobacter species was C. sakazakii and C. dublinensis, followed by C. malonaticus and C. turicensis. Two interesting clusters with more than 90% similarities were identified by pulsed-field gel electrophoresis (PFGE). The C1 cluster (four isolates) indicated these strains were derived from a common source and persisted in the food production environment for an extended time. The C2 cluster (six isolates) indicated the pathogen could be transmitted via cereal processing. Our research provided baseline data for Cronobacter in diverse cereals and was helpful for understanding Cronobacter transmission. The results also indicate that additional control measures should be developed to reduce the risk of infection by these opportunistic pathogenic bacteria.

Cutibacterium avidum is phylogenetically diverse with a subpopulation being adapted to the infant gut

The infant gut harbors a diverse microbial community consisting of several taxa whose persistence depends on adaptation to the ecosystem. In healthy breast-fed infants, the gut microbiota is dominated by Bifidobacterium spp.. Cutibacterium avidum is among the initial colonizers, however, the phylogenetic relationship of infant fecal isolates to isolates from other body sites, and C. avidum carbon utilization related to the infant gut ecosystem have been little investigated. In this study, we investigated the phylogenetic and phenotypic diversity of 28 C. avidum strains, including 16 strains isolated from feces of healthy infants. We investigated the in vitro capacity of C. avidum infant isolates to degrade and consume carbon sources present in the infant gut, and metabolic interactions of C. avidum with infant associated Bifidobacterium longum subsp. infantis and Bifidobacterium bifidum. Isolates of C. avidum showed genetic heterogeneity. C. avidum consumed d- and l-lactate, glycerol, glucose, galactose, N-acetyl-d-glucosamine and maltodextrins. Alpha-galactosidase- and β-glucuronidase activity were a trait of a group of non-hemolytic strains, which were mostly isolated from infant feces. Beta-glucuronidase activity correlated with the ability to ferment glucuronic acid. Co-cultivation with B. infantis and B. bifidum enhanced C. avidum growth and production of propionate, confirming metabolic cross-feeding. This study highlights the phylogenetic and functional diversity of C. avidum, their role as secondary glycan degraders and propionate producers, and suggests adaptation of a subpopulation to the infant gut.

Updates on the Cronobacter Genus

There has been considerable concern related to Cronobacter spp. in foods, especially due to their highlighted association with neonatal infections through the ingestion of reconstituted powdered infant formula (PIF). This concern resulted in improved microbiological criteria recommendations by the Codex Alimentarius Commission and revised WHO advice on the preparation of infant feeds. In recent years, the diversity of the genus has been well described, and various detection and typing methods have been developed. This review considers our current knowledge of the genus and how DNA-sequence-based methods have contributed considerably to research into improved detection methods and more reliable identification procedures, genotyping schemes, and genomic analysis. The broader occurrence of Cronobacter in food ingredients, finished products, and food manufacturing environments is covered. This review also highlights the significance of clonal lineages in microbial source tracking and the use of CRISPR-cas array profiling.

DNA-Sequence Based Typing of the Cronobacter Genus Using MLST, CRISPR-cas Array and Capsular Profiling

The Cronobacter genus is composed of seven species, within which a number of pathovars have been described. The most notable infections by Cronobacter spp. are of infants through the consumption of contaminated infant formula. The description of the genus has greatly improved in recent years through DNA sequencing techniques, and this has led to a robust means of identification. However some species are highly clonal and this limits the ability to discriminate between unrelated strains by some methods of genotyping. This article updates the application of three genotyping methods across the Cronobacter genus. The three genotyping methods were multilocus sequence typing (MLST), capsular profiling of the K-antigen and colanic acid (CA) biosynthesis regions, and CRISPR-cas array profiling. A total of 1654 MLST profiled and 286 whole genome sequenced strains, available by open access at the PubMLST Cronobacter database, were used this analysis. The predominance of C. sakazakii and C. malonaticus in clinical infections was confirmed. The majority of clinical strains being in the C. sakazakii clonal complexes (CC) 1 and 4, sequence types (ST) 8 and 12 and C. malonaticus ST7. The capsular profile K2:CA2, previously proposed as being strongly associated with C. sakazakii and C. malonaticus isolates from severe neonatal infections, was also found in C. turicensis, C. dublinensis and C. universalis. The majority of CRISPR-cas types across the genus was the I-E (Ecoli) type. Some strains of C. dublinensis and C. muytjensii encoded the I-F (Ypseudo) type, and others lacked the cas gene loci. The significance of the expanding profiling will be of benefit to researchers as well as governmental and industrial risk assessors.

Molecular Surveillance of Cronobacter spp. Isolated from a Wide Variety of Foods from 44 Different Countries by Sequence Typing of 16S rRNA, rpoB and O-Antigen Genes

Cronobacter spp. are emerging infectious bacteria that can cause acute meningitis and necrotizing enterocolitis in neonatal and immunocompromised individuals. Although this opportunistic human-pathogenic microorganism has been isolated from a wide variety of food and environmental samples, it has been primarily linked to foodborne outbreaks associated with powdered infant formula. The U.S. Food and Drug Administration use the presence of these microbes as one of the criteria to assess food adulteration and to implement regulatory actions. In this study, we have examined 195 aliquots of enrichments from the nine major categories of foods (including baby and medical food, dairy products, dried food, frozen food, pet food, produce, ready-to-eat snacks, seafood, and spices) from 44 countries using conventional microbiological and molecular techniques. The typical colonies of Cronobacter were then identified by VITEK2 and real-time PCR. Subsequently, sequence typing was performed on the 51 recovered Cronobacter isolates at the 16S rRNA, rpoB and seven O-antigen loci for species identification in order to accomplish an effective surveillance program for the control and prevention of foodborne illnesses.

Genetic Characterization of Cronobacter sakazakii Recovered from the Environmental Surveillance Samples During a Sporadic Case Investigation of Foodborne Illness

Cronobacter sakazakii is an opportunistic human-pathogenic bacterium known to cause acute meningitis and necrotizing enterocolitis in neonates and immunocompromised individuals. This human-pathogenic microorganism has been isolated from a variety of food and environmental samples, and has been also linked to foodborne outbreaks associated with powdered infant formula (PIF). The U.S. Food and Drug Administration have a policy of zero tolerance of these organisms in PIF. Thus, this agency utilizes the presence of these microorganisms as one of the criteria in implementing regulatory actions and assessing adulteration of food products of public health importance. In this study, we recovered two isolates of Cronobacter from the 91 environmental swab samples during an investigation of sporadic case of foodborne illness following conventional microbiological protocols. The isolated typical colonies were identified using VITEK2 and real-time PCR protocols. The recovered Cronobacter isolates were then characterized for species identification by sequencing the 16S rRNA locus. Further, multilocus sequence typing (MLST) was accomplished characterizing seven known C. sakazakii-specific MLST loci (atpD, fusA, glnS, gltB, gyrB, infB, and pps). Results of this study confirmed all of the recovered Cronobacter isolates from the environmental swab samples to be C. sakazakii. The MLST profile matched with the published profile of the complex 31 of C. sakazakii. Thus, rRNA and 7-loci MLST-based sequencing protocols are robust techniques for rapid detection and differentiation of Cronobacter species, and these molecular diagnostic tools can be used in implementing successful surveillance program and in the control and prevention of foodborne illness.