Abstract
We used whole genome sequencing (WGS) as an epidemiologic surveillance tool to elucidate the transmission dynamics of Shiga toxin-producing Escherichia coli (STEC) strains along the beef production chain in South Africa. Isolates were obtained from a cattle farm, abattoirs and retail outlets. Isolates were analysed using WGS on a MiSeq platform (Illumina, San Diego, CA, USA) and phylogenetic analysis was carried out. Of the 85 isolates, 39% (33) carried the stx gene and 61% (52) had lost the stx gene. The prevalence of stx subtypes was as follows; stx1a 55% (18/33), stx1b 52% (17/33), stx2a 55% (18/33), stx2b 27% (9/33), stx2dB 30% (10/33) and stx2d1A 15% (5/33). Thirty-five different serogenotypes were detected, of which 65% (56) were flagellar H-antigens and 34% (29) were both O-antigens and flagellar H-antigens. We identified 50 different sequence types (STs), and only nine of the isolates were assigned to three different clonal complexes. Core genome phylogenetic analysis revealed genetic relatedness, and isolates clustered mainly according to their STs and serogenotypes regardless of stx subtypes. This study provides evidence of horizontal transmission and recirculation of STEC strains in Gauteng province and demonstrates that every stage of the beef production chain plays a significant role in STEC entry into the food chain.