Multicenter evaluation of a sequence-based protocol for subtyping Shiga toxins and standardizing Stx nomenclature

Abstract

When Shiga toxin-producing Escherichia coli (STEC) strains emerged as agents of human disease, two types of toxin were identified: Shiga toxin type 1 (Stx1) (almost identical to Shiga toxin produced by Shigella dysenteriae type 1) and the immunologically distinct type 2 (Stx2). Subsequently, numerous STEC strains have been characterized that express toxins with variations in amino acid sequence, some of which confer unique biological properties. These variants were grouped within the Stx1 or Stx2 type and often assigned names to indicate that they were not identical in sequence or phenotype to the main Stx1 or Stx2 type. A lack of specificity or consistency in toxin nomenclature has led to much confusion in the characterization of STEC strains. Because serious outcomes of infection have been attributed to certain Stx subtypes and less so with others, we sought to better define the toxin subtypes within the main Stx1 and Stx2 types. We compared the levels of relatedness of 285 valid sequence variants of Stx1 and Stx2 and identified common sequences characteristic of each of three Stx/Stx1 and seven Stx2 subtypes. A novel, simple PCR subtyping method was developed, independently tested on a battery of 48 prototypic STEC strains, and improved at six clinical and research centers to test the reproducibility, sensitivity, and specificity of the PCR. Using a consistent schema for nomenclature of the Stx toxins and stx genes by phylogenetic sequence-based relatedness of the holotoxin proteins, we developed a typing approach that should obviate the need to bioassay each newly described toxin and that predicts important biological characteristics.

Fig 1

Three clusters of Stx and Stx1: neighbor-joining cluster analysis of 13 unique sequences as described in the text and the proposed new designations.

Fig 2

Five clusters of Stx2: neighbor-joining cluster analysis of 93 unique sequences as described in the text and the proposed new designations.

Fig 3

Maximum parsimony tree of 57 aa sequences: Stx2a, Stx2c, and Stx2d analyzed separately.