Genomic Diversity, Virulence Gene, and Prophage Arrays of Bovine and Human Shiga Toxigenic and Enteropathogenic Escherichia coli Strains Isolated in Hungary.

Escherichia coli belonging to the enterohemorrhagic (EHEC), Shiga toxin-producing (STEC) and atypical enteropathogenic (aEPEC) pathotypes are significant foodborne zoonotic pathogens posing serious health risks, with healthy cattle as their main reservoir. A representative sampling of Hungarian cattle farms during 2017-2018 yielded a prevalence of 6.5 and 5.8% for STEC and aEPEC out of 309 samples. The draft genomes of twelve STEC (of them 9 EHEC) and four aEPEC of bovine origin were determined. For comparative purposes, we also included 3 EHEC and 2 aEPEC strains of human origin, as well four commensal isolates and one extraintestinal pathogenic E. coli (ExPEC) obtained from animals in a final set of 26 strains for a WGS-based analysis. Apart from key virulence genes, these isolates harbored several additional virulence genes with arrays characteristic for the site of isolation. The most frequent insertion site of Shiga toxin (stx) encoding prophages was yehV for the Stx1 prophage and wrbA and sbcB for Stx2. For O157:H7 strains, the locus of enterocyte effacement pathogenicity island was present at the selC site, with integration at pheV for other serotypes, and pheU in the case of O26:H11 strains. Several LEE-negative STEC and aEPEC as well as commensal isolates carried additional prophages, with an average of ten prophage regions per isolate. Comparative phylogenomic analysis showed no clear separation between bovine and human lineages among the isolates characterized in the current study. Similarities in virulence gene arrays and close phylogenetic relations of bovine and human isolates underline the zoonotic potential of bovine aEPEC and STEC and emphasize the need for frequent monitoring of these pathogens in livestock.

Multiplex polymerase chain reaction typing scheme based on Escherichia coli O157:H7 Sakai prophage (Sp)-associated genes.

OBJECTIVES: Escherichia coli strains of the O157 serogroup include significant foodborne pathogens: enterohemorrhagic E. coli (EHEC) and enteropathogenic E. coli, which are responsible for a considerable number of hospitalizations and deaths worldwide each year. There is a constant need for rapid, reliable, and easy-to-use methods for their identification, typing, and phylogenetic classification. In this study, we proposed a new multiplex polymerase chain reaction (PCR)-based typing system for pathogenic E. coli, focusing on the O157 serogroup. METHODS: We designed primers targeting 12 lambdoid prophage regions carried by the prototypic polylysogenic strain of EHEC, the O157:H7 Sakai strain. The reactions were tested in vitro as well as in silico with the PubMLST database. RESULTS: The PCR assays can be grouped into four multiplex reactions, and their results can be given as a four-digit code. In vitro and in silico testing showed that these Sakai prophage regions are prevalent not only in E. coli O157 strains but also in Shiga toxigenic E. coli non-O157 strains and the method provides appropriate resolution. CONCLUSIONS: The proposed method could be a valuable tool in epidemiologic tracing and preliminary phylogenetic grouping of this diverse group of pathogens.

Characterisation of new anti-O157 bacteriophages of bovine origin representing three genera.

Shiga-toxin-producing Escherichia coli (STEC) strains of the serogroup O157 are foodborne pathogens associated with severe clinical disease. As antibiotics are counter-indicated for treatment of these infections, they represent prime candidates for targeted application of bacteriophages to reduce infection burden. In this study, we characterised lytic bacteriophages representing three phage genera for activity against E. coli O157 strains. The phages vb_EcoM_bov9_1 (Tequatrovirus), vb_EcoM_bov11CS3 (Vequintavirus), and vb_EcoS_bov25_1D (Dhillonvirus) showed effective lysis of enterohaemorrhagic E. coli EHEC O157:H7 strains, while also exhibiting activity against other strains of the O157 serogroup, as well as of the 'big six' (STEC) serogroups, albeit with lower efficiency. They had a burst size of 293, 127 and 18 per cell and a latent period of 35, 5 and 30 min, respectively. In situ challenge experiments using the O157 Sakai strain on minced beef showed a reduction by 2-3-fold when treated with phages at a 0.1 MOI (multiplicity of infection), and approximately 1 log reduction when exposed to MOI values of 10 and 100. A cocktail of the phages, applied at 10 × and 100 × MOI showed 2 to 3 log reduction when samples were treated at room temperature, and all treatments at 37 °C with 100 × MOI resulted in a 5 to 6 log reduction in cell count. Our results indicate that the phages vb_EcoM_bov9_1 and vb_EcoM_bov11CS3, which have higher burst sizes, are promising candidates for biocontrol experiments aimed at the eradication of E. coli O157 strains in animals or foodstuff.

Novel prophage-like sequences in Mycoplasma anserisalpingitidis.

Mycoplasma anserisalpingitidis is a bacterial waterfowl pathogen. In these days of growing antibiotic resistance, it is necessary to search for alternative methods of defense against Mycoplasma impacts in flocks. In order to identify prophage-like sequences, three established bioinformatics tools (PHASTER, PhiSpy, Prophage Hunter) were used in this study for the in silico screening of 82 M. anserisalpingitidis whole genomes. The VIBRANT software was used as a novel approach to further investigate the possibility of prophages in the sequences. The commonly used softwares found prophage-like sequences in the strains, but the results were inconclusive. The VIBRANT search resulted in multiple hits, and many of them were over 10,000 base pairs (bp). These putative prophages are comparable in size to the few described mycoplasma phages. The translated coding DNA sequences of the putative prophages were checked with protein BLAST. The functions of the proteins found by the BLASTP search are common among bacteriophages. The BLASTN search of the sequences found that many of these were more similar to the M. anatis NCTC 10156 strain, rather than the available M. anserisalpingitidis strains. The initial screening pointed at the presence of novel bacteriophages in the M. anserisalpingitidis and M. anatis strains. The VIBRANT search results were very similar to each other and none of these sequences were part of the core genome of M. anserisalpingitidis, with a few exceptions. The VIBRANT analysis explored presumably intact, novel prophages.

Complete Genome Sequences of Novel Bovine T4, rv5-Like, and Dhillonviruses Effective against Escherichia coli O157.

Enterohemorrhagic Escherichia coli (EHEC) isolates of serotype O157:H7 are serious foodborne zoonotic pathogens and prime targets for biocontrol using bacteriophages. We report on the complete genome sequences of 11 novel lytic bacteriophages, representing three viral genera, isolated from cattle in Hungary that target E. coli O157 strains.

Genome Analysis of a Historical Shigella dysenteriae Serotype 1 Strain Carrying a Conserved Stx Prophage Region.

Shigella dysenteriae are significant agents of bacillary dysentery, accounting for a considerable number of illnesses with high morbidity worldwide. The Shiga toxin (Stx) encoded by a defective prophage is the key virulence factor of S. dysenteriae type 1 (SD1) strains. Here we present the full genome sequence of an SD1 strain HNCMB 20080 isolated in 1954, compare it to other sequenced SD1 genomes, and assess the diversity of Stx-prophages harbored by previously sequenced SD1 strains. The genome of HNCMB 20080 consists of a chromosome sized 4,393,622 bp containing 5,183 CDSs, as well as two small plasmids. Comparative genomic analysis revealed a high degree of uniformity among SD1 genomes, including the structure of Stx prophage regions, which we found to form two subgroups termed PT-I and PT-II. All PT-I strains are members of the sequence type (ST) 146 or ST260, while the only PT-II harboring strain, Sd1617 proved to be ST untypeable. In accordance with data from previous reports, the Stx1 prophage could not be induced from HNCMB 20080. Our cumulative data do not support the notion that stx-harboring phages in STEC are derived from historical SD1 isolates.

R18C is a new viable P2-like bacteriophage of rabbit origin infecting Citrobacter rodentium and Shigella sonnei strains.

Here, we report a novel virulent P2-like bacteriophage, R18C, isolated from rabbit faeces, which, in addition to Escherichia coli K-12 strains, was able to be propagated on Citrobacter rodentium strain ICC169 and a range of Shigella sonnei strains with high efficiency of plating (EOP). It represents the first lytic bacteriophage originating from rabbit and the first infectious P2-like phage of animal origin. In the three characteristic moron-containing regions of P2-like phages, R18C contains genes with unknown function that have so far only been found in cryptic P2-like prophages.

Complete genome sequence of C130_2, a novel myovirus infecting pathogenic Escherichia coli and Shigella strains.

The genome sequence of a novel virulent bacteriophage, named " C130_2", that is morphologically a member of the family Myoviridae is reported. The 41,775-base-pair double-stranded DNA genome of C130_2 contains 59 ORFs but exhibits overall low sequence similarity to bacteriophage genomes for which sequences are publicly available. Phylogenetic analysis indicated that C130_2 represents a new phage type. C130_2 could be propagated well on enterohemorrhagic Escherichia coli (EHEC) O157:H7 and other pathogenic E. coli strains, as well as on strains of various Shigella species.

Identification and characterization of new broad host-range rV5-like coliphages C203 and P206 directed against enterobacteria.

We isolated and characterized two novel rV5-like lytic bacteriophages from independently collected food samples. Nucleotide sequence analysis revealed that these phages have linear double-stranded DNA genomes comprising 138,073 bp with 213 CDS and 5 tRNA genes. The two genomes contain completely identical nucleotide sequence, albeit there is a 10,718 bp-long shift in the sequence. The GC content of the phage genomes was 43.7% and they showed high general homology to rV5-like phages. The new phages were termed C203 and P206. The genome of both phages contains a unique ORF that encodes for a putative phage homing endonuclease. The phage produced clear plaques with a burst size of approx. 1000 viral particles and a latent period of 60 min. Morphological investigation indicated that the new phages are members of the family Myoviridae with an approximate head length of 85 nm, tail length of 75 nm, and a head width of 96 nm. C203 and P206 exhibit a broad and uniform host range, which included enterohemorrhagic Escherichia coli strains of serogroup O157, multi drug resistant (MDR) E. coli strains of various sero- and pathotypes, and both Shigella sonnei and S. dysenteriae strains. C203 and P206 both effectively reduced the number of living EHEC O157:H7 Sakai in experimentally inoculated minced meat. The same broad host range, the lack of any virulence related genes, the stability and its short latent period suggest that these newly found phages could be suitable candidates as a bio-control agents against food-borne pathogenic Enterobacteria.

Identification and Characterization of T5-Like Bacteriophages Representing Two Novel Subgroups from Food Products.

During recent years, interest in the use of bacteriophages as biocontrol agents against foodborne pathogens has increased, particularly for members of the family Enterobacteriaceae, with pathogenic Escherichia coli, Shigella, and Salmonella strains among them. Here, we report the isolation and characterisation of 12 novel T5-like bacteriophages from confiscated food samples. All bacterophages effectively lysed E. coli K-12 strains and were able to infect pathogenic E. coli strains representing enterohaemorrhagic (EHEC), enteropathogenic (EPEC), enterotoxigenic (ETEC), and enteroinvasive (EIEC) pathotypes, Shigella dysenteriae, S. sonnei strains, as well as multidrug-resistant (MDR) E. coli and multiple strains representing different Salmonella enterica serovars. All the bacteriophages exhibited Siphoviridae morphology. Whole genome sequencing of the novel T5-like bacteriophages showed that they represent two distinct groups, with the genome-based grouping correlating to the different host spectra. As these bacteriophages are of food origin, their stability and lack of any virulence genes, as well as their broad and mutually complementary host spectrum makes these new T5-like bacteriophages valuable candidates for use as biocontrol agents against foodborne pathogenic enterobacteria.

Comparative Genomic and Phylogenetic Analysis of a Shiga Toxin Producing Shigella sonnei (STSS) Strain.

Shigella strains are important agents of bacillary dysentery, and in recent years Shigella sonnei has emerged as the leading cause of shigellosis in industrialized and rapidly developing countries. More recently, several S. sonnei and Shigella flexneri strains producing Shiga toxin (Stx) have been reported from sporadic cases and from an outbreak in America. In the present study we aimed to shed light on the evolution of a recently identified Shiga toxin producing S. sonnei (STSS) isolated in Europe. Here we report the first completely assembled whole genome sequence of a multidrug resistant (MDR) Stx-producing S. sonnei (STSS) clinical strain and reveal its phylogenetic relations. STSS 75/02 proved to be resistant to ampicillin, streptomycin, tetracycline, chloramphenicol, thrimetoprim, and sulfomethoxazol. The genome of STSS 75/02 contains a 4,891,717 nt chromosome and seven plasmids including the 214 kb invasion plasmid (pInv) harboring type III secretion system genes and associated effectors. The chromosome harbors 23 prophage regions including the Stx1 converting prophage. The genome carries all virulence determinants necessary for an enteroinvasive lifestyle, as well as the Stx1 encoding gene cluster within an earlier described inducible converting prophage. In silico SNP genotyping of the assembled genome as well as 438 complete or draft S. sonnei genomes downloaded from NCBI GenBank revealed that S. sonnei 75/02 belongs to the more recently diverged global MDR lineage (IIIc). Targeted screening of 1131 next-generation sequencing projects taken from NCBI Short Read Archive of confirms that only a few S. sonnei isolates are Stx positive. Our results suggest that the acquisition of Stx phages could have occurred in different environments as independent events and that multiple horizontal transfers are responsible for the appearance of Stx phages in S. sonnei strains.

Cytolethal distending toxin producing Escherichia coli O157:H43 strain T22 represents a novel evolutionary lineage within the O157 serogroup.

Enterohemorrhagic Escherichia coli (EHEC) O157:H7/NM strains are significant foodborne pathogens intensively studied, while other sero- and pathotypes of the O157 serogroup only began to receive more attention. Here we report the first genome sequence of a cytolethal distending toxin (CDT-V) producing E. coli O157:H43 strain (T22) isolated from cattle. The genome consists of a 4.9Mb chromosome assembled into three contigs and one plasmid of 82.4kb. Comparative genomic investigations conducted with the core genomes of representative E. coli strains in GenBank (n=62) confirmed the separation of T22 from the EHEC and enteropathogenic (EPEC) O157 lineages. Gene content based pangenome analysis revealed as many as 261 T22-specific coding sequences without orthologs in EDL933 EHEC O157 prototypic and two phylogenetically related commensal E. coli strains. The genome sequence revealed 10 prophage-like regions which harbor several virulence-associated genes including cdt and heat-labile enterotoxin (LT-II) encoding operons. Our results indicate that the evolutionary path of T22 is largely independent from that of EHEC and EPEC O157:H7/NM strains. Thus, the CDT-producing T22 E. coli O157:H43 strain represents a unique lineage of E. coli O157.

Comparative analysis of the Shiga toxin converting bacteriophage first detected in Shigella sonnei.

Here we report the first complete nucleotide sequence of a Shiga toxin (Stx) converting phage from a Shigella sonnei clinical isolate that harbors stx1 operon, first identified in the chromosome of Shigella dysenteriae type 1. The phage named Shigella phage 75/02 Stx displayed Podoviridae morphology. It proved to be transferable to Escherichia coli K-12 strains, and cytotoxicity of the lysogenized strains was demonstrated in Vero cell cultures. Genomic analysis revealed that the prophage genome is circular and its size is 60,875 nt that corresponds to 76 ORFs. The genome of Shigella phage 75/02 Stx shows a great degree of mosaic structure and its architecture is related to lambdoid phages. All the deduced proteins, including the 37 hypothetical proteins showed significant homologies to Stx phage proteins present in databases. The phage uniformly inserted into the ynfG oxidoreductase gene framed by phage integrase and antirepressor genes in parental S. sonnei and in the three lysogenized K-12 strains C600, DH5α and MG1655. The Stx1 prophage proved to be stable in its bacterial hosts and remained inducible.

Cytolethal distending toxin A, B and C subunit proteins are necessary for the genotoxic effect of Escherichia coli CDT-V.

Cytolethal distending toxins (CDT) are considered the prototype of inhibitory cyclomodulins, and are produced by a wide range of Gram-negative pathogenic bacteria, including Escherichia coli strains of various sero- and pathotypes. CDT is a heterotripartite toxin consisting of three protein subunits, CdtA, CdtB and CdtC. The active subunit, CdtB has DNase activity and causes DNA damage and cell cycle arrest in the target cell. However, several studies have highlighted different roles for CdtA and CdtC subunits. In order to reveal the necessity of CdtA and CdtC subunit proteins in the CDT-specific phenotype, expression clones containing the cdt-V subunit genes were constructed. Using cell culture assays, we demonstrated that clones expressing only the CdtB subunit or in combination with only CdtA or CdtC were unable to trigger the specific cell cycle arrest and changes in cell morphology in HeLa cells. At the same time, the recombinant clone harbouring the whole cdt-V operon caused all the CDT-associated characteristic phenotypes. All these results verify that all the three CDT subunit proteins are necessary for the genotoxic effect caused by CDT-V.

A novel transducible chimeric phage from Escherichia coli O157:H7 Sakai strain encoding Stx1 production.

Shiga toxin-producing Escherichia coli (STEC), and especially enterohaemorrhagic E. coli (EHEC) are important, highly virulent zoonotic and food-borne pathogens. The genes encoding their key virulence factors, the Shiga toxins, are distributed by converting bacteriophages, the Stx phages. In this study we isolated a new type of inducible Stx phage carrying the stx1 gene cluster from the prototypic EHEC O157:H7 Sakai strain. The phage showed Podoviridae morphology, and was capable of converting the E. coli K-12 MG1655 strain to Shiga toxin-producing phenotype. The majority of the phage genes originate from the stx2-encoding Sakai prophage Sp5, with major rearrangements in its genome. Beside certain minor recombinations, the genomic region originally containing the stx2 genes in Sp5 was replaced by a region containing six open reading frames from prophage Sp15 including stx1 genes. The rearranged genome, together with the carriage of stx1 genes, the morphology and the capability of lysogenic conversion represent a new type of recombinant Stx1 converting phage from the Sakai strain.

Sequence variability of P2-like prophage genomes carrying the cytolethal distending toxin V operon in Escherichia coli O157.

Cytolethal distending toxins (CDT) are potent cytotoxins of several Gram-negative pathogenic bacteria, including Escherichia coli, in which five types (CDT-I to CDT-V) have been identified so far. CDT-V is frequently associated with Shiga-toxigenic E. coli (STEC), enterohemorrhagic E. coli (EHEC) O157 strains, and strains not fitting any established pathotypes. In this study, we were the first to sequence and annotate a 31.2-kb-long, noninducible P2-like prophage carrying the cdt-V operon from an stx- and eae-negative E. coli O157:H43 strain of bovine origin. The cdt-V operon is integrated in the place of the tin and old phage immunity genes (termed the TO region) of the prophage, and the prophage itself is integrated into the bacterial chromosome between the housekeeping genes cpxP and fieF. The presence of P2-like genes (n = 20) was investigated in a further five CDT-V-positive bovine E. coli O157 strains of various serotypes, three EHEC O157:NM strains, four strains expressing other variants of CDT, and eight CDT-negative strains. All but one CDT-V-positive atypical O157 strain uniformly carried all the investigated genomic regions of P2-like phages, while the EHEC O157 strains missed three regions and the CDT-V-negative strains carried only a few P2-like sequences. Our results suggest that P2-like phages play a role in the dissemination of cdt-V between E. coli O157 strains and that after integration into the bacterial chromosome, they adapted to the respective hosts and became temperate.

Draft Genome Sequence of an Escherichia coli O157:H43 Strain Isolated from Cattle.

Here we report the draft genome sequence of an Escherichia coli O157:H43 strain, designated T22, with an atypical virulence gene profile and isolated from healthy cattle. T22 produces cytolethal distending toxin V (CDT-V) and belongs to phylogenetic group B1 and sequence type 155 (ST155).

The long polar fimbriae operon and its flanking regions in bovine Escherichia coli O157:H43 and STEC O136:H12 strains.

Long polar fimbriae (Lpf) are intestinal adhesins and important virulence factors of pathogenic Escherichia coli strains. We cloned and sequenced the lpf2-1 operon (lpf2ABCD) and its flanking regions of an intimin- and Shiga toxin-negative atypical E. coli O157:H43 strain of bovine origin and also sequenced the lpf2-1 operon of six additional atypical O157 bovine E. coli strains of various serotypes. Nucleotide sequence comparison of these lpf operons showed sequence conservation as they contained only four polymorphic nucleotide positions. Investigation of these O157 strains as well as 13 E. coli Reference Collection (ECOR) strains carrying the lpf2-1 allele revealed high degree of sequence conservation in the lpf2 flanking regions. The lpf2-1 allele is also present in a bovine Shiga toxin-producing E. coli STEC O136:H12 strain, and in vitro adherence assays revealed that the absence of lpf2-1 in this strain did not affect its host cell-binding properties. Our data indicate that lpf2 loci are highly conserved in E. coli isolates; however, its role in adherence might be masked by other uncharacterized adhesins.

Allelic types of long polar fimbriae in bovine and human Escherichia coli O157 strains.

Long polar fimbriae (Lpf) are recently discovered adhesins and increasingly important genetic markers of pathogenic Escherichia coli strains. The presence and genotype diversity of Lpf operons was screened in a collection of 97 Escherichia coli O157 strains representing different pathotypes, isolated from healthy cattle (n = 43) and human patients (n = 54) in several countries. Individual structural genes of Lpf were scanned by PCR, and allelic variants were detected with a recently developed typing scheme. Ninety-five strains carried at least one whole Lpf operon (genes lpfABCD and/or lpfABCDE). The 64 enterohaemorrhagic (EHEC) and 24 enteropathogenic (EPEC) strains all carried two Lpf operons, allele 3 of lpfA1 and allele 2 of lpfA2, a combination characteristic of the O157:H7/NM serotype. Out of the 9 bovine atypical (AT; stx-, eae-) strains, 7 carried one complete Lpf operon, allele 1 of lpfA2. The atypical strains belonged to main phylogenetic groups A and B1, while the EHEC and EPEC strains were from group D. Lpf variants carried by the 72 strains of the Escherichia coli Reference Collection (ECOR) were determined with the same typing scheme. Alleles were detected in 25 strains, of which 6 were found negative for the respective Lpf operons in earlier studies. The marker value of the Lpf allelic combination for the O157:H7/NM serotype was confirmed, and further evidence was given for the presence of at least two different genetic lineages of atypical bovine E. coli O157 strains.

Preferential ligation during TA-cloning of multitemplate PCR products--a factor causing bias in microbial community structure analysis.

The description of microbial community structure is always biased by the selectivity of the methods applied. Although TA cloning of PCR amplified community DNA is one of the most widely used techniques in bacterial community analysis, no thorough comparative testing has been carried out on different TA cloning systems. In this study, we measured and compared the selectivity of two widely used TA-cloning kits in experimental setups where the length heterogeneity of the inserts modeled the natural length variation of the 16S rRNA gene and the 16S-23S intergenic spacer region. Both TOPO TA (Invitrogen, CA USA) and pGem-T vector system (Promega, WI USA) cloning kits showed significant and reproducible insert size related selectivity. The effect of ligation time and temperature was also studied in case of the pGem-T vector system. We compared the performance of the two cloning kits on an environmental sample, along with a semiquantitative community fingerprinting method to gain reference data free of cloning bias. The two clone libraries showed significantly different compositions, and were also differing from the community structure revealed by length heterogeneity PCR.