Comparative evaluation of genomic and laboratory approaches for determination of Shiga toxin subtypes in Escherichia coli

Abstract

The determination of Shiga toxin (ST) subtypes can be an important element in the risk characterization of foodborne STproducing Escherichia coli (STEC) isolates for making risk management decisions. ST subtyping methods include PCR techniques based on electrophoretic or pyrosequencing analysis of amplicons and in silico techniques based on whole genome sequence analysis using algorithms that can be readily incorporated into bioinformatics analysis pipelines for characterization of isolates by their genetic composition. The choice of technique will depend on the performance characteristics of the method and an individual laboratory's access to specialized equipment or personnel. We developed two whole genome sequence-based ST subtyping tools: (i) an in silico PCR algorithm requiring genome assembly to replicate a reference PCR-based method developed by the Statens Serum Institut (SSI) and (ii) an assembly-independent routine in which raw sequencing results are mapped to a database of known ST subtype sequence variants (V-Typer). These tools were evaluated alongside the SSI reference PCR method and a recently described PCR-based pyrosequencing technique. The V-Typer method results corresponded closely with the reference method in the analysis of 67 STEC cultures obtained from a World Health Organization National Reference Laboratory. In contrast, the in silico PCR method failed to detect ST subtypes in several cases, a result which we attribute to assemblyinduced errors typically encountered with repetitive gene sequences. The V-Typer can be readily integrated into bioinformatics protocols used in the identification and characterization of foodborne STEC isolates.