O157 Genome Research Articles

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.

Virulence, Antimicrobial Resistance Properties and Phylogenetic Background of Non-H7 Enteropathogenic Escherichia coli O157.

Escherichia coli (E.coli) O157 that do not produce Shiga toxin and do not possess flagellar antigen H7 are of diverse H serotypes. In this study, the antibiotic resistance properties, genotype of a set of virulence associated genes and the phylogenetic background of E. coli O157:non-H7 groups were compared. Whole genome sequencing was performed on fourteen O157:non-H7 isolates collected in the STEC-ID-net study. The genomes were compared with E. coli O157 genomes and a typical Enteropathogenic E. coli (tEPEC) genome downloaded from NCBI. Twenty-six (86%) of the analyzed genomes had the intimin encoding gene eae but of different types mostly correlating with their H types, e.g., H16, H26, H39, and H45 carried intimin type ε, β, κ, and α, respectively. They belonged to several E. coli phylogenetic groups, i.e., to phylogenetic group A, B1, B2, and D. Seven (50%) of our collected O157:non-H7 isolates were resistant to two or more antibiotics. Several mobile genetic elements, such as plasmids, insertion elements, and pathogenicity islands, carrying a set of virulence and resistance genes were found in the E. coli O157:non-H7 isolates. Core genome phylogenetic analysis showed that O157:non-H7 isolates probably evolved from different phylogenetic lineages and were distantly related to the E. coli O157:H7 lineage. We hypothesize that independent acquisition of mobile genetic elements by isolates of different lineages have contributed to the different molecular features of the O157:non-H7 strains. Although distantly related to the STEC O157, E. coli O157:non-H7 isolates from multiple genetic background could be considered as pathogen of concern for their diverse virulence and antibiotic resistance properties.

Putative classification of clades of enterohemorrhagic Escherichia coli O157 using an IS-printing system.

Enterohemorrhagic Escherichia coli O157 (O157) strains can be classified in clades by single nucleotide polymorphisms (SNPs), but this analysis requires significant laboratory effort. As the distribution of insertion sequence (IS) 629 insertions has been reported to be biased among different clades, O157 isolates can be putatively classified in clades by comparison with an IS629 distribution database. A database of the IS629 distribution in O157 strains isolated in Chiba Prefecture and their classification in clades was determined by SNP analysis and IS-printing, an easy and quick analytical tool for IS629 in the O157 genome. The IS629 distribution in O157 strains isolated in Fukuoka and Yamagata Prefectures was determined by IS-printing. These strains were putatively classified in clades by Relative Likelihood calculations that compared the IS-printing data and the IS629 distribution database. Concordance Ratios were calculated, which compared the number of strains putatively classified in a clade by Relative Likelihood to the number of strains classified in that clade by SNP analysis. For the Fukuoka and Yamagata strains, the Concordance Ratios for clades 3, 6 and 8 were 97-100%, for clade 7 about 88%, and for clades 2 and 12 over 90%. In conclusion, O157 clade 2, 3, 6, 7, 8 and 12 strains could be putatively classified by IS-printing. This study demonstrated that enterohemorrhagic E.coli O157 (O157) strains could be putatively classified in clades using an IS-printing system. IS-printing was previously developed as a relatively quick and easy tool for analysis of insertion sequence 629 in the O157 genome. Since most local government public health institutes in Japan carry out IS-printing for early detection of O157 outbreaks, these data should be useful for putative classification of O157 strains in each area.

Rapid Detection of Viable Escherichia coli O157 by Coupling Propidium Monoazide with Loop-Mediated Isothermal Amplification

Conventional molecular detection methods cannot distinguish between viable and dead Escherichia coli O157 cells. In this study, the loop-mediated isothermal amplification (LAMP) method combined with propidium monoazide (PMA) treatment was developed to selectively detect viable E. coli O157 cells. Four primers, including outer primers and inner primers, were specially designed for the recognition of six distinct sequences on the serogroups (O157) of the specific rfbE gene of the E. coli O157 genome. PMA selectively penetrated through the compromised cell membranes and intercalated into DNA. Amplification of DNA from dead cells was completely inhibited by 3.0 μg/ml PMA, whereas the DNA derived from viable cells was amplified remarkably within 1 h by PMA-LAMP. Exhibiting high sensitivity and specificity, PMA-LAMP is a suitable method for evaluating the inactivation efficacy of slightly acidic electrolyzed water in broth. PMA-LAMP can selectively detect viable E. coli O157 cells. This study offers a novel molecular detection method to distinguish between viable and dead E. coli O157 cells.

Polymorphic Amplified Typing Sequences and Pulsed-Field Gel Electrophoresis Yield Comparable Results in the Strain Typing of a Diverse Set of BovineEscherichia coliO157:H7 Isolates

Polymorphic amplified typing sequences (PATS), a PCR-based Escherichia coli O157:H7 (O157) strain typing system, targets insertions-deletions and single nucleotide polymorphisms at XbaI and AvrII restriction enzyme sites, respectively, and the virulence genes (stx1, stx2, eae, hlyA) in the O157 genome. In this study, the ability of PATS to discriminate O157 isolates associated with cattle was evaluated. An in-depth comparison of 25 bovine O157 isolates, from different geographic locations across Northwest United States, showed that about 85% of these isolates shared the same dendogram clade by PATS and pulsed-field gel electrophoresis (PFGE), irrespective of the restriction enzyme sites targeted. The Pearson's correlation coefficient, r, calculated at about 0.4, 0.3, and 0.4 for XbaI-based, AvrII-based and combined-enzymes PATS and PFGE similarities, respectively, indicating that these profiles shared a good but not high correlation, an expected inference given that the two techniques discriminate differently. Isolates that grouped differently were better matched to their locations using PATS. Overall, PATS discriminated the bovine O157 isolates without interpretive biases or sophisticated analytical software, and effectively complemented while not duplicating PFGE. With its quick turnaround time, PATS has excellent potential as a convenient tool for early epidemiological or food safety investigations, enabling rapid notification/implementation of quarantine measures.

Insertion sequence-excision enhancer removes transposable elements from bacterial genomes and induces various genomic deletions

Insertion sequences (ISs) are the simplest transposable elements and are widely distributed in bacteria. It has long been thought that IS excision rarely occurs in bacterial cells because most bacteria exhibit no end-joining activity to regenerate donor DNA after IS excision. Recently, however, we found that excision of IS629, an IS3 family member, occurs frequently in Escherichia coli O157. In this paper, we describe a protein IS-excision enhancer (IEE) that promotes IS629 excision from the O157 genome in an IS transposase-dependent manner. Various types of genomic deletions are also generated on IEE-mediated IS excision, and IEE promotes the excision of other IS3 family members and ISs from several other IS families. These data and the phylogeny of IEE homologues found in a broad range of bacteria suggest that IEE proteins have coevolved with IS elements and have pivotal roles in bacterial genome evolution by inducing IS removal and genomic deletion.

Chromosomal instability in enterohaemorrhagic Escherichia coli O157:H7: impact on adherence, tellurite resistance and colony phenotype

Tellurite (Tel) resistant enterohaemorrhagic Escherichia coli (EHEC) O157:H7 is a global pathogen. In strain EDL933 Tel resistance (TelR) is encoded by duplicate ter cluster in O islands (OI) 43 and 48, which also harbour iha, encoding the adhesin and siderophore receptor Iha. We identified five EHEC O157:H7 strains that differentiate into large (L) colonies and small (S) colonies with high and low Tel minimal inhibitory concentrations (MICs) respectively. S colonies (Tel-MICs ≤ 4 µg ml−1) sustained large internal deletions within the TelR OIs via homologous recombination between IS elements and lost ter and iha. Moreover, complete excision of the islands occurred by site-specific recombination between flanking direct repeats. Complete excision of OI 43 and OI 48 occurred in 1.81 × 10−3 and 1.97 × 10−4 cells in culture, respectively; internal deletion of OI 48 was more frequent (9.7 × 10−1 cells). Under iron limitation that promotes iha transcription, iha-negative derivatives adhered less well to human intestinal epithelial cells and grew slower than did their iha-positive counterparts. Experiments utilizing iha deletion and complementation mutants identified Iha as the major factor responsible for these phenotypic differences. Spontaneous deletions affecting TelR OIs contribute to EHEC O157 genome plasticity and might impair virulence and/or fitness.

Inference of the impact of insertion sequence (IS) elements on bacterial genome diversification through analysis of small-size structural polymorphisms in Escherichia coli O157 genomes

Mobile genetic elements play important roles in the evolution and diversification of bacterial genomes. In enterohemorrhagic Escherichia coli O157, a major factor that affects genomic diversity is prophages, which generate most of the large-size structural polymorphisms (LSSPs) observed in O157 genomes. Here, we describe the results of a systematic analysis of numerous small-size structural polymorphisms (SSSPs) that were detected by comparing the genomes of eight clinical isolates with a sequenced strain, O157 Sakai. Most of the SSSPs were generated by genetic events associated with only two insertion sequence (IS) elements, IS629 and ISEc8, and a number of genes that were inactivated or deleted by these events were identified. Simple excisions of IS629 and small deletions (footprints) formed by the excision of IS629, both of which are rarely described in bacteria, were also detected. In addition, the distribution of IS elements was highly biased toward prophages, prophage-like integrative elements, and plasmids. Based on these and our previous results, we conclude that, in addition to prophages, these two IS elements are major contributors to the genomic diversification of O157 strains and that LSSPs have been generated mainly by bacteriophages and SSSPs by IS elements. We also suggest that IS elements possibly play a role in the inactivation and immobilization of incoming phages and plasmids. Taken together, our results reveal the true impact of IS elements on the diversification of bacterial genomes and highlight their novel role in genome evolution.

Probing genomic diversity and evolution of Escherichia coli O157 by single nucleotide polymorphisms

Infections by Shiga toxin-producing Escherichia coli O157:H7 (STEC O157) are the predominant cause of bloody diarrhea and hemolytic uremic syndrome in the United States. In silico comparison of the two complete STEC O157 genomes (Sakai and EDL933) revealed a strikingly high level of sequence identity in orthologous protein-coding genes, limiting the use of nucleotide sequences to study the evolution and epidemiology of this bacterial pathogen. To systematically examine single nucleotide polymorphisms (SNPs) at a genome scale, we designed comparative genome sequencing microarrays and analyzed 1199 chromosomal genes (a total of 1,167,948 bp) and 92,721 bp of the large virulence plasmid (pO157) of eleven outbreak-associated STEC O157 strains. We discovered 906 SNPs in 523 chromosomal genes and observed a high level of DNA polymorphisms among the pO157 plasmids. Based on a uniform rate of synonymous substitution for Escherichia coli and Salmonella enterica (4.7x10(-9) per site per year), we estimate that the most recent common ancestor of the contemporary beta-glucuronidase-negative, non-sorbitolfermenting STEC O157 strains existed ca. 40 thousand years ago. The phylogeny of the STEC O157 strains based on the informative synonymous SNPs was compared to the maximum parsimony trees inferred from pulsed-field gel electrophoresis and multilocus variable numbers of tandem repeats analysis. The topological discrepancies indicate that, in contrast to the synonymous mutations, parts of STEC O157 genomes have evolved through different mechanisms with highly variable divergence rates. The SNP loci reported here will provide useful genetic markers for developing high-throughput methods for fine-resolution genotyping of STEC O157. Functional characterization of nucleotide polymorphisms should shed new insights on the evolution, epidemiology, and pathogenesis of STEC O157 and related pathogens.

Complexity of the Genomic Diversity in Enterohemorrhagic Escherichia coli O157 Revealed by the Combinational Use of the O157 Sakai OligoDNA Microarray and the Whole Genome PCR scanning

Escherichia coli O157, an etiological agent of hemorrhagic colitis and hemolytic uremic syndrome, is one of the leading worldwide public health threats. Genome sequencing of two O157 strains have revealed that the chromosome is comprised of a 4.1 Mb backbone shared by K-12 and a total of 1.4 Mb O157-specific sequences. Most of the large O157-specific sequences are prophages and prophage-like elements, which have carried many virulence genes into the O157 genome. This suggests that bacteriophages have played the key roles in the emergence of O157. The Whole Genome PCR Scanning (WGPScanning) analysis of O157 strains, on the other hand, revealed a high level of genomic diversity in O157. Variation of prophages has also been suggested as a major factor generating such diversity. In this study, we analyzed the gene content of O157 strains, by an oligoDNA microarray, using the same set of strains as examined by the WGPScanning method. Although most of the strains were typical O157 : H7, they differed remarkably in gene composition, particularly in those on prophages, and we identified more than 400 'variably absent or present' genes which included virulence-related genes. This confirms the role of prophages in generating the genomic diversity, and raises a possibility that some level of variation in potential virulence is present among O157 strains. Fine comparison of the two datasets obtained by microarray and WGPScanning provided much further details on the O157 genome diversity than illustrated by each method alone, indicating the usefulness of this combinational approach in the genomic comparison of closely related strains.

Evolution of Genomic Content in the Stepwise Emergence ofEscherichia coliO157:H7

Genome comparisons have demonstrated that dramatic genetic change often underlies the emergence of new bacterial pathogens. Evolutionary analysis of Escherichia coli O157:H7, a pathogen that has emerged as a worldwide public health threat in the past two decades, has posited that this toxin-producing pathogen evolved in a series of steps from O55:H7, a recent ancestor of a nontoxigenic pathogenic clone associated with infantile diarrhea. We used comparative genomic hybridization with 50-mer oligonucleotide microarrays containing probes from both pathogenic and nonpathogenic genomes to infer when genes were acquired and lost. Many ancillary virulence genes identified in the O157 genome were already present in an O55:H7-like progenitor, with 27 of 33 genomic islands of >5 kb and specific for O157:H7 (O islands) that were acquired intact before the split from this immediate ancestor. Most (85%) of variably absent or present genes are part of prophages or phage-like elements. Divergence in gene content among these closely related strains was approximately 140 times greater than divergence at the nucleotide sequence level. A >100-kb region around the O-antigen gene cluster contained highly divergent sequences and also appears to be duplicated in its entirety in one lineage, suggesting that the whole region was cotransferred in the antigenic shift from O55 to O157. The beta-glucuronidase-positive O157 variants, although phylogenetically closest to the Sakai strain, were divergent for multiple adherence factors. These observations suggest that, in addition to gains and losses of phage elements, O157:H7 genomes are rapidly diverging and radiating into new niches as the pathogen disseminates.

Insertions, deletions, and single-nucleotide polymorphisms at rare restriction enzyme sites enhance discriminatory power of polymorphic amplified typing sequences, a novel strain typing system for Escherichia coli O157:H7.

Polymorphic amplified typing sequences (PATS) for Escherichia coli O157:H7 (O157) was previously based on indels containing XbaI restriction enzyme sites occurring in O-island sequences of the O157 genome. This strain-typing system, referred to as XbaI-based PATS, typed every O157 isolate tested in a reproducible, rapid, straightforward, and easy-to-interpret manner and had technical advantages over pulsed-field gel electrophoresis (PFGE). However, the system was less discriminatory than PFGE and was unable to differentiate fully between unrelated isolates. To overcome this drawback, we enhanced PATS by using another infrequently cutting restriction enzyme, AvrII (also known as BlnI), to identify additional polymorphic regions that could increase the discriminatory ability of PATS typing. Referred to as AvrII-based PATS, the system identified seven new polymorphic regions in the O157 genome. Unlike XbaI, polymorphisms involving AvrII sites were caused by both indels and single-nucleotide polymorphisms occurring in O-island and backbone sequences of the O157 genome. AvrII-based PATS by itself provided poor discrimination of the O157 isolates tested. However, when primer pairs amplifying the seven polymorphic AvrII sites were combined with those amplifying the eight polymorphic XbaI sites (combined PATS), the discriminatory power of PATS was enhanced. Combined PATS matched related O157 isolates better than PFGE while differentiating between unrelated isolates. PATS typed every O157 isolate tested and directly targeted polymorphic sequences responsible for differences in the restriction digest patterns of O157 genomic DNA, utilizing PCR rather than relying on gel electrophoresis. This enabled PATS to resolve the ambiguity in PFGE typing, including that arising from the "more distantly related" and "untypeable" profiles.

Multiple-Locus Variable-Number Tandem-Repeats Analysis of Escherichia coli O157 using PCR multiplexing and multi-colored capillary electrophoresis

The Multiple-Locus Variable-Number Tandem-Repeats Analysis (MLVA) method is currently being used as the primary typing tool for Shiga-toxin-producing Escherichia coli (STEC) O157 isolates in our laboratory. The initial assay was performed using a single fluorescent dye and the different patterns were assigned using a gel image. Here, we present a significantly improved assay using multiple dye colors and enhanced PCR multiplexing to increase speed, and ease the interpretation of the results. The different MLVA patterns are now based on allele sizes entered as character values, thus removing the uncertainties introduced when analyzing band patterns from the gel image. We additionally propose an easy numbering scheme for the identification of separate isolates that will facilitate exchange of typing data. Seventy-two human and animal strains of Shiga-toxin-producing E. coli O157 were used for the development of the improved MLVA assay. The method is based on capillary separation of multiplexed PCR products of VNTR loci in the E. coli O157 genome labeled with multiple fluorescent dyes. The different alleles at each locus were then assigned to allele numbers, which were used for strain comparison.

A mosaic pathogenicity island made up of the locus of enterocyte effacement and a pathogenicity island of Escherichia coli O157:H7 is frequently present in attaching and effacing E. coli.

Enteropathogenic Escherichia coli (EPEC) and enterohemorragic E. coli (EHEC) possess a pathogenicity island (PAI), termed the locus of enterocyte effacement (LEE), which confers the capability to cause the characteristic attaching and effacing lesions of the brush border. Due to this common property, these organisms are also termed attaching and effacing E. coli (AEEC). Sequencing of the EHEC O157 genome recently revealed the presence of other putative PAIs in the chromosome of this organism. In this article, we report on the presence of four of those PAIs in a panel of 133 E. coli strains belonging to different pathogroups and serotypes. One of these PAIs, termed O122 in strain EDL 933 and SpLE3 in strain Sakai, was observed in most of the AEEC strains examined but not in the other groups of E. coli. It was also found to contain the virulence-associated gene efa1/lifA. In EHEC O157, PAI O122 is located 0.7 Mb away from the LEE. Conversely, we demonstrated that in many EHEC non-O157 strains and EPEC strains belonging to eight serogroups, PAI O122 and the LEE are physically linked to form a cointegrated structure. This structure can be considered a mosaic PAI that could have been acquired originally by AEEC. In some clones, such as EHEC O157, the LEE-O122 mosaic PAI might have undergone recombinational events, resulting in the insertion of the portion referred to as PAI O122 in a different location.

Genomic diversity of enterohemorrhagic Escherichia coli O157 revealed by whole genome PCR scanning.

Enterohemorrhagic Escherichia coli O157 is one of the leading worldwide public health concerns, causing large outbreaks of hemorrhagic colitis as well as numerous small outbreaks and sporadic cases. The variability of restriction enzyme-digestion patterns of O157 genomes, which is widely used to distinguish strains in the molecular epidemiology of O157 infections, suggests the presence of some genomic diversity among the strains. Based on the complete genome sequence of O157 Sakai, we analyzed the whole genome structures of eight O157 strains displaying diverse XbaI-digestion patterns by a systematic PCR analysis that we have named whole genome PCR scanning. This analysis identified not only the O157-specific sequences that are highly conserved among the strains, but also revealed an unexpectedly high degree of genomic diversity. In particular, prophages, including Shiga toxin-transducing phages, exhibited extensive structural and positional diversity, implying that variation of bacteriophages is a major factor in generating genomic diversity among the O157 lineage.

Analysis of codon usage diversity of bacterial genes with a self-organizing map (SOM): characterization of horizontally transferred genes with emphasis on the E. coli O157 genome

With increases in the amounts of available DNA sequence data, it has become increasingly important to develop tools for comprehensive systematic analysis and comparison of species-specific characteristics of protein-coding sequences for a wide variety of genomes. In the present study, we used a novel neural-network algorithm, a self-organizing map (SOM), to efficiently and comprehensively analyze codon usage in approximately 60,000 genes from 29 bacterial species simultaneously. This SOM makes it possible to cluster and visualize genes of individual species separately at a much higher resolution than can be obtained with principal component analysis. The organization of the SOM can be explained by the genome G+C% and tRNA compositions of the individual species. We used SOM to examine codon usage heterogeneity in the E. coli O157 genome, which contains ‘O157-unique segments’ (O-islands), and showed that SOM is a powerful tool for characterization of horizontally transferred genes.