Whole-Genome Sequencing Analysis of a stx-Negative Escherichia coli O63:H6 Isolate Associated with Hemolytic Uremic Syndrome

Shiga toxin-encoding genes (stx) of enterohemorrhagic Escherichia coli (EHEC) can be lost during infection or in vitro cultivation, and in clinical practice, it is difficult to distinguish EHEC that have lost stx (EHEC-LST) from enteropathogenic E. coli (EPEC), as both are stx-negative and eae-positive. In this study, we performed whole-genome sequencing (WGS) of a stx-negative, eae-positive E. coli O63:H6 isolate from a child with hemolytic uremic syndrome and compared its genome with those of nine E. coli O63:H6 strains in public databases. Virulence gene profiles were analyzed and core-genome multilocus sequence typing (cgMLST) was conducted. The virulence gene profile of our isolate was consistent with EHEC, except for the absence of stx, and the isolate clustered with seven EHEC strains but was distant from two EPEC strains in cgMLST. In genome alignment, our isolate exhibited a high nucleotide identity with EHEC strain 377323_2f but displayed a gap corresponding to the stx-harboring prophage sequence. Overall, our isolate was genetically closely related to EHEC strains, consistent with this being an EHEC-LST strain. As EHEC-LST may be misdiagnosed as EPEC in routine laboratories, comparative genomic analysis using WGS can be useful to determine whether stx-negative and eae-positive isolates are EHEC-LST or EPEC.