Whole genome sequencing based typing and characterisation of Shiga-toxin producing Escherichia coli strains belonging to O157 and O26 serotypes and isolated in dairy farms.

In the present study, the genetic relationships as well as the virulome and resistome of newly sequenced O26 and O157 Shiga-toxin producing E. coli (STEC) isolates, collected from dairy farms in Italy, were investigated in comparison to publicly available genomes collected worldwide. The whole genome of Italian isolates was sequenced on Illumina MiSeq Platform. Reads quality control, de novo draft genome assembly, species confirmation and the 7- loci Multi-Locus Sequence Type assignment were performed using INNUca pipeline. Reference-based SNPs calling was performed on O157 and O26 genomes, separately, mapping contigs to high-quality finished genomes. Virulence and antimicrobial resistance determinants were detected in silico using the tool ABRicate. Phylogenetic reconstructions revealed that genomes clustered mainly based on their 7-loci MLST type. The virulome of tested genomes included 190 determinants. O157 genomes carried chu genes associated to heme mediated iron uptake, whereas O26 genomes harboured genes ybt associated to siderophore mediated iron uptake. Resistome analysis showed the presence of tet (34) on all but one O157 genomes and on only one O26 genomes. Only 4 genomes carried genes associated to multiresistance. In the present study, the genes chu and ybt were identified as potential biomarker for the differentiation of O157 and O26 serotypes.