Purpose: In this study, a diverse group of 357 Cronobacter strains obtained during environmental sampling of five powdered infant formula production facilities in Ireland and Germany, and 93 milk powder and cheese manufacturing facilities in the USA were analyzed. This study describes the genomic diversity and phylogenetic relationships among these strains.
Methods: Genomic DNA was isolated using a QIAcube workstation and hybridizations were performed using a previously described pan-genome DNA microarray. Whole genome sequencing (WGS) was carried out on a MiSeq platform using Nextera XT chemistry.
Results: Microarray analysis (MA) determined that 22 phylogenetically-related Csak strains possessed genes of the C. malonaticus-like malonate utilization operon; and all, but one strain, could utilize malonate. The strains phylogenetically grouped as a separate clade among the Csak cluster. Additionally, these strains possessed the Csak O:2 serotype and 21 of the strains were identified as sequence type (ST) ST64. The malonate-negative strain was shown by MA and WGS to lack the transcriptional regulator, mdcR; and its ST was ST1. WGS confirmed that the nine gene, ~7.7 kbp malonate utilization operon was located in these strains between two flanking genes, gyrB and katG.
Significance: Until this investigation, the presence of malonate-positive Csak strains that are associated with foods was underappreciated, possibly leading to misidentification when relying on phenotypic identification alone. A custom designed pan genome microarray was useful in characterizing the total genome content of these CsakO:2, ST64/ST1, malonate-positive strains; and showed for the first time that these Csak strains were highly related in total gene content, serotype and ST. This study further establishes MA as a powerful platform for genomics research of Cronobacter.