Skip to content
1887

Microbial Genomics logo Microbial Genomics

Volume 9, Issue 2

Alteration of a Shiga toxin-encoding phage associated with a change in toxin production level and disease severity in Open Access

Abstract

Among the nine clades of Shiga toxin (Stx)-producing O157:H7, clade 8 is thought to be highly pathogenic, as it causes severe disease more often than other clades. Two subclades have been proposed, but there are conflicting reports on intersubclade differences in Stx2 levels, although Stx2 production is a risk factor for severe disease development. The global population structure of clade 8 has also yet to be fully elucidated. Here, we present genome analyses of a global clade 8 strain set (=510), including 147 Japanese strains sequenced in this study. The complete genome sequences of 18 of the 147 strains were determined to perform detailed clade-wide genome analyses together with 17 publicly available closed genomes. Intraclade variations in Stx2 production level and disease severity were also re-evaluated within the phylogenetic context. Based on phylogenomic analysis, clade 8 was divided into four lineages corresponding to the previously proposed SNP genotypes (SGs): SG8_30, SG8_31A, SG8_31B and SG8_32. SG8_30 and the common ancestor of the other SGs were first separated, with SG8_31A and SG8_31B emerging from the latter and SG8_32 emerging from SG8_31B. Comparison of 35 closed genomes revealed the overall structure of chromosomes and pO157 virulence plasmids and the prophage contents to be well conserved. However, Stx2a phages exhibit notable genomic diversity, even though all are integrated into the locus, indicating that subtype changes in Stx2a phage occurred from the γ subtype to its variant (γ_v1) in SG8_31A and from γ to δ in SG8_31B and SG8_32 via replacement of parts or almost entire phage genomes, respectively. We further show that SG8_30 strains (all carrying γ Stx2a phages) produce significantly higher levels of Stx2 and cause severe disease more frequently than SG8_32 strains (all carrying δ Stx2a phages). Clear conclusions on SG8_31A and SG8_31B cannot be made due to the small number of strains available, but as SG8_31A (carrying γ_v1 Stx2a phages) contains strains that produce much more Stx2 than SG8_30 strains, attention should also be paid to this SG.

  • Received:
  • Accepted:
  • Published Online:
Keyword(s): disease severity , genome analysis , population structure , STEC O157:H7 clade 8 , Stx2 production level and Stx2a phage subtype
Funding
This study was supported by the:
  • AMED (Award 22fk0108611h0502)
    • Principle Award Recipient: TetsuyaHayashi
© 2023 The Authors
  • This is an open-access article distributed under the terms of the Creative Commons Attribution License.
Loading

Article metrics loading...

/content/journal/mgen/10.1099/mgen.0.000935
2023-02-23
2025-04-27
Loading full text...

Full text loading...

/deliver/fulltext/mgen/9/2/mgen000935.html?itemId=/content/journal/mgen/10.1099/mgen.0.000935&mimeType=html&fmt=ahah

References

  1. Hayashi T, Miyata T, Taniguchi I, Nakamura K, Gotoh Y et al. Alteration of a Shiga toxin-encoding phage associated with a change in toxin production level and disease severity in Escherichia coli. Figshare 2022 [View Article]
     [Google Scholar]
  2. Croxen MA, Law RJ, Scholz R, Keeney KM, Wlodarska M et al. Recent advances in understanding enteric pathogenic Escherichia coli. Clin Microbiol Rev 2013; 26:822–880 [View Article]
     [Google Scholar]
  3. Kaper JB, Nataro JP, Mobley HL. Pathogenic Escherichia coli. Nat Rev Microbiol 2004; 2:123–140 [View Article]
     [Google Scholar]
  4. Croxen MA, Finlay BB. Molecular mechanisms of Escherichia coli pathogenicity. Nat Rev Microbiol 2010; 8:26–38 [View Article]
     [Google Scholar]
  5. Tobe T, Beatson SA, Taniguchi H, Abe H, Bailey CM et al. An extensive repertoire of type III secretion effectors in Escherichia coli O157 and the role of lambdoid phages in their dissemination. Proc Natl Acad Sci USA 2006; 103:14941–14946 [View Article]
     [Google Scholar]
  6. Schmidt MA. LEEways: tales of EPEC, ATEC and EHEC. Cell Microbiol 2010; 12:1544–1552 [View Article]
     [Google Scholar]
  7. Trachtman H, Austin C, Lewinski M, Stahl RAK. Renal and neurological involvement in typical Shiga toxin-associated HUS. Nat Rev Nephrol 2012; 8:658–669 [View Article] [PubMed]
     [Google Scholar]
  8. Karpman D. Management of Shiga toxin-associated Escherichia coli-induced haemolytic uraemic syndrome: randomized clinical trials are needed. Nephrol Dial Transplant 2012; 27:3669–3674 [View Article] [PubMed]
     [Google Scholar]
  9. Johannes L, Römer W. Shiga toxins – from cell biology to biomedical applications. Nat Rev Microbiol 2010; 8:105–116 [View Article]
     [Google Scholar]
  10. Scheutz F, Teel LD, Beutin L, Piérard D, Buvens G et al. Multicenter evaluation of a sequence-based protocol for subtyping Shiga toxins and standardizing Stx nomenclature. J Clin Microbiol 2012; 50:2951–2963 [View Article] [PubMed]
     [Google Scholar]
  11. Bai X, Fu S, Zhang J, Fan R, Xu Y et al. Identification and pathogenomic analysis of an Escherichia coli strain producing a novel Shiga toxin 2 subtype. Sci Rep 2018; 8:6756 [View Article]
     [Google Scholar]
  12. Lacher DW, Gangiredla J, Patel I, Elkins CA, Feng PCH. Use of the Escherichia coli identification microarray for characterizing the health risks of Shiga toxin-producing Escherichia coli isolated from foods. J Food Prot 2016; 79:1656–1662 [View Article]
     [Google Scholar]
  13. Yang X, Bai X, Zhang J, Sun H, Fu S et al. Escherichia coli strains producing a novel Shiga toxin 2 subtype circulate in China. Int J Med Microbiol 2020; 310:151377 [View Article]
     [Google Scholar]
  14. De Rauw K, Buyl R, Jacquinet S, Piérard D. Risk determinants for the development of typical haemolytic uremic syndrome in Belgium and proposition of a new virulence typing algorithm for Shiga toxin-producing Escherichia coli. Epidemiol Infect 2018; 147:e6 [View Article]
     [Google Scholar]
  15. Ylinen E, Salmenlinna S, Halkilahti J, Jahnukainen T, Korhonen L et al. Hemolytic uremic syndrome caused by Shiga toxin-producing Escherichia coli in children: incidence, risk factors, and clinical outcome. Pediatr Nephrol 2020; 35:1749–1759 [View Article] [PubMed]
     [Google Scholar]
  16. Brandal LT, Wester AL, Lange H, Løbersli I, Lindstedt B-A et al. Shiga toxin-producing Escherichia coli infections in Norway, 1992–2012: characterization of isolates and identification of risk factors for haemolytic uremic syndrome. BMC Infect Dis 2015; 15:324 [View Article]
     [Google Scholar]
  17. Manning SD, Motiwala AS, Springman AC, Qi W, Lacher DW et al. Variation in virulence among clades of Escherichia coli O157:H7 associated with disease outbreaks. Proc Natl Acad Sci USA 2008; 105:4868–4873 [View Article]
     [Google Scholar]
  18. Iyoda S, Manning SD, Seto K, Kimata K, Isobe J et al. Phylogenetic clades 6 and 8 of enterohemorrhagic Escherichia coli O157:H7 with particular stx subtypes are more frequently found in isolates from hemolytic uremic syndrome patients than from asymptomatic carriers. Open Forum Infect Dis 2014; 1:ofu061 [View Article]
     [Google Scholar]
  19. Ogura Y, Mondal SI, Islam MR, Mako T, Arisawa K et al. The Shiga toxin 2 production level in enterohemorrhagic Escherichia coli O157:H7 is correlated with the subtypes of toxin-encoding phage. Sci Rep 2015; 5:16663 [View Article] [PubMed]
     [Google Scholar]
  20. Kulasekara BR, Jacobs M, Zhou Y, Wu Z, Sims E et al. Analysis of the genome of the Escherichia coli O157:H7 2006 spinach-associated outbreak isolate indicates candidate genes that may enhance virulence. Infect Immun 2009; 77:3713–3721 [View Article] [PubMed]
     [Google Scholar]
  21. Eppinger M, Mammel MK, Leclerc JE, Ravel J, Cebula TA. Genomic anatomy of Escherichia coli O157:H7 outbreaks. Proc Natl Acad Sci USA 2011; 108:20142–20147 [View Article]
     [Google Scholar]
  22. Hirai S, Yokoyama E, Wakui T, Ishige T, Nakamura M. Correction: enterohemorrhagic Escherichia coli O157 subclade 8b strains in Chiba Prefecture, Japan, produced larger amounts of Shiga toxin 2 than strains in subclade 8a and other clades. PLoS One 2018; 13:e0196027 [View Article]
     [Google Scholar]
  23. Riordan JT, Viswanath SB, Manning SD, Whittam TS. Genetic differentiation of Escherichia coli O157:H7 clades associated with human disease by real-time PCR. J Clin Microbiol 2008; 46:2070–2073 [View Article] [PubMed]
     [Google Scholar]
  24. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol 1990; 215:403–410 [View Article] [PubMed]
     [Google Scholar]
  25. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009; 25:1754–1760 [View Article] [PubMed]
     [Google Scholar]
  26. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J et al. The sequence alignment/map format and SAMtools. Bioinformatics 2009; 25:2078–2079 [View Article]
     [Google Scholar]
  27. Kajitani R, Toshimoto K, Noguchi H, Toyoda A, Ogura Y et al. Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res 2014; 24:1384–1395 [View Article] [PubMed]
     [Google Scholar]
  28. Kajitani R, Yoshimura D, Ogura Y, Gotoh Y, Hayashi T et al. Platanus_B: an accurate de novo assembler for bacterial genomes using an iterative error-removal process. DNA Res 2020; 27:dsaa014 [View Article]
     [Google Scholar]
  29. De Coster W, D’Hert S, Schultz DT, Cruts M, Van Broeckhoven C. NanoPack: visualizing and processing long-read sequencing data. Bioinformatics 2018; 34:2666–2669 [View Article] [PubMed]
     [Google Scholar]
  30. Murigneux V, Roberts LW, Forde BM, Phan M-D, Nhu NTK et al. MicroPIPE: validating an end-to-end workflow for high-quality complete bacterial genome construction. BMC Genomics 2021; 22:474 [View Article]
     [Google Scholar]
  31. Kolmogorov M, Yuan J, Lin Y, Pevzner PA. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol 2019; 37:540–546 [View Article] [PubMed]
     [Google Scholar]
  32. Vaser R, Sović I, Nagarajan N, Šikić M. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res 2017; 27:737–746 [View Article] [PubMed]
     [Google Scholar]
  33. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 2018; 34:3094–3100 [View Article] [PubMed]
     [Google Scholar]
  34. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES et al. Integrative genomics viewer. Nat Biotechnol 2011; 29:24–26 [View Article] [PubMed]
     [Google Scholar]
  35. Ohtsubo Y, Ikeda-Ohtsubo W, Nagata Y, Tsuda M. GenomeMatcher: a graphical user interface for DNA sequence comparison. BMC Bioinformatics 2008; 9:376 [View Article] [PubMed]
     [Google Scholar]
  36. Yoshimura D, Kajitani R, Gotoh Y, Katahira K, Okuno M et al. Evaluation of SNP calling methods for closely related bacterial isolates and a novel high-accuracy pipeline: BactSNP. Microb Genom 2019; 5:e000261 [View Article]
     [Google Scholar]
  37. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res 2015; 43:e15 [View Article]
     [Google Scholar]
  38. Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006; 22:2688–2690 [View Article] [PubMed]
     [Google Scholar]
  39. Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res 2016; 44:W242–W245 [View Article]
     [Google Scholar]
  40. Tonkin-Hill G, Lees JA, Bentley SD, Frost SDW, Corander J. Fast hierarchical Bayesian analysis of population structure. Nucleic Acids Res 2019; 47:5539–5549 [View Article] [PubMed]
     [Google Scholar]
  41. Nishida R, Nakamura K, Taniguchi I, Murase K, Ooka T et al. The global population structure and evolutionary history of the acquisition of major virulence factor-encoding genetic elements in Shiga toxin-producing Escherichia coli O121:H19. Microb Genom 2021; 7:000716 [View Article]
     [Google Scholar]
  42. Menardo F, Loiseau C, Brites D, Coscolla M, Gygli SM et al. Treemmer: a tool to reduce large phylogenetic datasets with minimal loss of diversity. BMC Bioinformatics 2018; 19:164 [View Article]
     [Google Scholar]
  43. Rambaut A, Lam TT, Max Carvalho L, Pybus OG. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol 2016; 2:vew007 [View Article] [PubMed]
     [Google Scholar]
  44. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 2012; 29:1969–1973 [View Article]
     [Google Scholar]
  45. Barnett DW, Garrison EK, Quinlan AR, Stromberg MP, Marth GT. BamTools: a C++ API and toolkit for analyzing and managing BAM files. Bioinformatics 2011; 27:1691–1692 [View Article] [PubMed]
     [Google Scholar]
  46. Tanizawa Y, Fujisawa T, Kaminuma E, Nakamura Y, Arita M. DFAST and DAGA: web-based integrated genome annotation tools and resources. Biosci Microbiota Food Health 2016; 35:173–184 [View Article] [PubMed]
     [Google Scholar]
  47. Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: a fast phage search tool. Nucleic Acids Res 2011; 39:W347–W352 [View Article]
     [Google Scholar]
  48. Arndt D, Grant JR, Marcu A, Sajed T, Pon A et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res 2016; 44:W16–W21 [View Article]
     [Google Scholar]
  49. Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res 2006; 34:D32–D36 [View Article]
     [Google Scholar]
  50. Ondov BD, Treangen TJ, Melsted P, Mallonee AB, Bergman NH et al. Mash: fast genome and metagenome distance estimation using MinHash. Genome Biol 2016; 17:132 [View Article]
     [Google Scholar]
  51. Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother 2014; 58:3895–3903 [View Article] [PubMed]
     [Google Scholar]
  52. Hofinger C, Karch H, Schmidt H. Structure and function of plasmid pColD157 of enterohemorrhagic Escherichia coli O157 and its distribution among strains from patients with diarrhea and hemolytic-uremic syndrome. J Clin Microbiol 1998; 36:24–29 [View Article] [PubMed]
     [Google Scholar]
  53. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother 2012; 67:2640–2644 [View Article] [PubMed]
     [Google Scholar]
  54. Nakamura K, Tokuda C, Arimitsu H, Etoh Y, Hamasaki M et al. Development of a homogeneous time-resolved FRET (HTRF) assay for the quantification of Shiga toxin 2 produced by E. coli. PeerJ 2021; 9:e11871 [View Article]
     [Google Scholar]
  55. Dallman TJ, Ashton PM, Byrne L, Perry NT, Petrovska L et al. Applying phylogenomics to understand the emergence of Shiga-toxin-producing Escherichia coli O157:H7 strains causing severe human disease in the UK. Microb Genom 2015; 1:e000029 [View Article]
     [Google Scholar]
  56. Rusconi B, Sanjar F, Koenig SSK, Mammel MK, Tarr PI et al. Whole genome sequencing for genomics-guided investigations of Escherichia coli O157:H7 outbreaks. Front Microbiol 2016; 7:985 [View Article]
     [Google Scholar]
  57. Hauser JR, Atitkar RR, Petro CD, Lindsey RL, Strockbine N et al. The virulence of Escherichia coli O157:H7 isolates in mice depends on Shiga toxin type 2a (Stx2a)-induction and high levels of Stx2a in stool. Front Cell Infect Microbiol 2020; 10:62 [View Article]
     [Google Scholar]
  58. Hosein AM, Breidt F, Smith CE. Modeling the effects of sodium chloride, acetic acid, and intracellular pH on survival of Escherichia coli O157:H7. Appl Environ Microbiol 2011; 77:889–895 [View Article]
     [Google Scholar]
  59. Lindsey RL, Batra D, Smith P, Patel PN, Tagg KA et al. PacBio genome sequences of Escherichia coli serotype O157:H7, diffusely adherent E. coli, and Salmonella enterica strains, all carrying plasmids with an mcr-1 resistance gene. Microbiol Resour Announc 2018; 7:14 [View Article]
     [Google Scholar]
  60. Hwang D, Kim SM, Kim HJ. Transcriptome changes and polymyxin resistance of acid-adapted Escherichia coli O157:H7 ATCC 43889. Gut Pathog 2020; 12:52 [View Article]
     [Google Scholar]
  61. Timme RE, Lafon PC, Balkey M, Adams JK, Wagner D et al. Gen-FS coordinated proficiency test data for genomic foodborne pathogen surveillance, 2017 and 2018 exercises. Sci Data 2020; 7:402 [View Article]
     [Google Scholar]
  62. Shaaban S, Cowley LA, McAteer SP, Jenkins C, Dallman TJ et al. Evolution of a zoonotic pathogen: investigating prophage diversity in enterohaemorrhagic Escherichia coli O157 by long-read sequencing. Microb Genom 2016; 2:000096 [View Article]
     [Google Scholar]
  63. Katani R, Cote R, Kudva IT, DebRoy C, Arthur TM et al. Comparative genomics of two super-shedder isolates of Escherichia coli O157:H7. PLoS One 2017; 12:e0182940 [View Article]
     [Google Scholar]
  64. Cote R, Katani R, Moreau MR, Kudva IT, Arthur TM et al. Comparative analysis of super-shedder strains of Escherichia coli O157:H7 reveals distinctive genomic features and a strongly aggregative adherent phenotype on bovine rectoanal junction squamous epithelial cells. PLoS One 2015; 10:e0116743 [View Article]
     [Google Scholar]
  65. Teng L, Lee S, Park D, Jeong KC. Genetic and functional analyses of virulence potential of an Escherichia coli O157:H7 strain isolated from super-shedder cattle. Front Cell Infect Microbiol 2020; 10:271 [View Article]
     [Google Scholar]
  66. Fitzgerald SF, Lupolova N, Shaaban S, Dallman TJ, Greig D et al. Genome structural variation in Escherichia coli O157:H7. Microb Genom 2021; 7:000682 [View Article]
     [Google Scholar]
  67. Poole SJ, Diner EJ, Aoki SK, Braaten BA, t’Kint de Roodenbeke C et al. Identification of functional toxin/immunity genes linked to contact-dependent growth inhibition (CDI) and rearrangement hotspot (Rhs) systems. PLoS Genet 2011; 7:e1002217 [View Article]
     [Google Scholar]
  68. Eppinger M, Mammel MK, Leclerc JE, Ravel J, Cebula TA. Genome signatures of Escherichia coli O157:H7 isolates from the bovine host reservoir. Appl Environ Microbiol 2011; 77:2916–2925 [View Article] [PubMed]
     [Google Scholar]
  69. Greig DR, Mikhail AFW, Dallman TJ, Jenkins C. Analysis Shiga toxin-encoding bacteriophage in a rare strain of Shiga toxin-producing Escherichia coli O157:H7 stx2a/stx2c. Front Microbiol 2020; 11:577658 [View Article]
     [Google Scholar]
  70. Zuppi M, Tozzoli R, Chiani P, Quiros P, Martinez-Velazquez A et al. Investigation on the evolution of Shiga toxin-converting phages based on whole genome sequencing. Front Microbiol 2020; 11:1472 [View Article]
     [Google Scholar]
  71. Chen C, Lewis CR, Goswami K, Roberts EL, DebRoy C et al. Identification and characterization of spontaneous deletions within the Sp11-Sp12 prophage region of Escherichia coli O157:H7 Sakai. Appl Environ Microbiol 2013; 79:1934–1941 [View Article] [PubMed]
     [Google Scholar]
  72. Nakamura K, Murase K, Sato MP, Toyoda A, Itoh T et al. Differential dynamics and impacts of prophages and plasmids on the pangenome and virulence factor repertoires of Shiga toxin-producing Escherichia coli O145:H28. Microb Genom 2020; 6:e000323 [View Article]
     [Google Scholar]
  73. Makino K, Yokoyama K, Kubota Y, Yutsudo CH, Kimura S et al. Complete nucleotide sequence of the prophage VT2-Sakai carrying the verotoxin 2 genes of the enterohemorrhagic Escherichia coli O157:H7 derived from the Sakai outbreak. Genes Genet Syst 1999; 74:227–239 [View Article] [PubMed]
     [Google Scholar]
  74. Plunkett G, Rose DJ, Durfee TJ, Blattner FR. Sequence of Shiga toxin 2 phage 933W from Escherichia coli O157:H7: Shiga toxin as a phage late-gene product. J Bacteriol 1999; 181:1767–1778 [View Article]
     [Google Scholar]
  75. Mondal SI, Islam MR, Sawaguchi A, Asadulghani M, Ooka T et al. Genes essential for the morphogenesis of the Shiga toxin 2-transducing phage from Escherichia coli O157:H7. Sci Rep 2016; 6:39036 [View Article] [PubMed]
     [Google Scholar]
  76. Court DL, Oppenheim AB, Adhya SL. A new look at bacteriophage lambda genetic networks. J Bacteriol 2007; 189:298–304 [View Article] [PubMed]
     [Google Scholar]
/content/journal/mgen/10.1099/mgen.0.000935
Loading
Alteration of a Shiga toxin-encoding phage associated with a change in toxin production level and disease severity in Escherichia coli
M Gen 9, 000935 (2023); https://doi.org/10.1099/mgen.0.000935
/content/journal/mgen/10.1099/mgen.0.000935
/content/journal/mgen/10.1099/mgen.0.000935
Loading

Data & Media loading...

Supplements

Loading data from figshare Loading data from figshare

Most cited Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error