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.

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.

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.

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.

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.