Comparative genomics reveals structural and functional features specific to the genome of a foodborne Escherichia coli O157:H7

BackgroundEscherichia coli O157:H7 (O157) has been linked to numerous foodborne disease outbreaks. The ability to rapidly sequence and analyze genomes is important for understanding epidemiology, virulence, survival, and evolution of outbreak strains. In the current study, we performed comparative genomics to determine structural and functional features of the genome of a foodborne O157 isolate NADC 6564 and infer its evolutionary relationship to other O157 strains.ResultsThe chromosome of NADC 6564 contained 5466 kb compared to reference strains Sakai (5498 kb) and EDL933 (5547 kb) and shared 41 of its 43 Linear Conserved Blocks (LCB) with the reference strains. However, 18 of 41 LCB had inverse orientation in NADC 6564 compared to the reference strains. NADC 6564 shared 18 of 19 bacteriophages with reference strains except that the chromosomal positioning of some of the phages differed among these strains. The additional phage (P19) of NADC 6564 was located on a 39-kb insertion element (IE) encoding several hypothetical proteins, an integrase, transposases, transcriptional regulators, an adhesin, and a phosphoethanolamine transferase (PEA). The complete homologs of the 39-kb IE were found in E. coli PCN061 of porcine origin. The IE-encoded PEA showed low homology (32–33%) to four other PEA in NADC 6564 and PEA linked to mobilizable colistin resistance in E. coli but was highly homologous (95%) to a PEA of uropathogenic, avian pathogenic, and enteroaggregative E. coli. NADC 6564 showed slightly higher minimum inhibitory concentration of colistin compared to the reference strains. The 39-kb IE also contained dndBCDE and dptFGH operons encoding DNA S-modification and a restriction pathway, linked to oxidative stress tolerance and self-defense against foreign DNA, respectively. Evolutionary tree analysis grouped NADC 6564 with lineage I O157 strains.ConclusionsThese results indicated that differential phage counts and different chromosomal positioning of many bacteriophages and genomic islands might have resulted in recombination events causing altered chromosomal organization in NADC 6564. Evolutionary analysis grouped NADC 6564 with lineage I strains and suggested its earlier divergence from these strains. The ability to perform S-DNA modification might affect tolerance of NADC 6564 to various stressors.

[1]  Samuel I. Miller,et al.  Analysis of the Genome of the Escherichia coli O157:H7 2006 Spinach-Associated Outbreak Isolate Indicates Candidate Genes That May Enhance Virulence , 2009, Infection and Immunity.

[2]  S. Huynh,et al.  RcsB Contributes to the Distinct Stress Fitness among Escherichia coli O157:H7 Curli Variants of the 1993 Hamburger-Associated Outbreak Strains , 2012, Applied and Environmental Microbiology.

[3]  V. Sharma,et al.  Complete Genome Sequences of Curli-Negative and Curli-Positive Isolates of Foodborne Escherichia coli O157:H7 Strain 86-24 , 2016, Genome Announcements.

[4]  David G. E. Smith,et al.  Escherichia coli O157 : H7 forms attaching and effacing lesions at the terminal rectum of cattle and colonization requires the LEE4 operon. , 2005, Microbiology.

[5]  T. Cebula,et al.  Genomic anatomy of Escherichia coli O157:H7 outbreaks , 2011, Proceedings of the National Academy of Sciences.

[6]  J. Keen,et al.  Mutations in the csgD Promoter Associated with Variations in Curli Expression in Certain Strains ofEscherichia coli O157:H7 , 2001, Applied and Environmental Microbiology.

[7]  Lukas Wagner,et al.  A Greedy Algorithm for Aligning DNA Sequences , 2000, J. Comput. Biol..

[8]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[9]  O. Gascuel,et al.  New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. , 2010, Systematic biology.

[10]  Nicola K. Petty,et al.  BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons , 2011, BMC Genomics.

[11]  J. Wells,et al.  A multistate outbreak of Escherichia coli O157:H7-associated bloody diarrhea and hemolytic uremic syndrome from hamburgers. The Washington experience. , 1994, JAMA.

[12]  B. Bosworth,et al.  Escherichia coli O157:H7 Requires Intimin for Enteropathogenicity in Calves , 1998, Infection and Immunity.

[13]  S. Huynh,et al.  Natural rpoS mutations contribute to population heterogeneity in Escherichia coli O157:H7 strains linked to the 2006 US spinach-associated outbreak. , 2014, Food microbiology.

[14]  Thomas K. F. Wong,et al.  ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates , 2017, Nature Methods.

[15]  Matthew R. Laird,et al.  IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets , 2017, Nucleic Acids Res..

[16]  Harry B. McGee,et al.  Hemorrhagic colitis associated with a rare Escherichia coli serotype , 1983 .

[17]  R. Hendriksen,et al.  Identification of a novel transposon-associated phosphoethanolamine transferase gene, mcr-5, conferring colistin resistance in d-tartrate fermenting Salmonella enterica subsp. enterica serovar Paratyphi B , 2017, The Journal of antimicrobial chemotherapy.

[18]  Alejandro A. Schäffer,et al.  Database indexing for production MegaBLAST searches , 2008, Bioinform..

[19]  Brian D. Ondov,et al.  The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes , 2014, Genome Biology.

[20]  T. Cebula,et al.  Genome Signatures of Escherichia coli O157:H7 Isolates from the Bovine Host Reservoir , 2011, Applied and Environmental Microbiology.

[21]  Zixin Deng,et al.  dndDB: A Database Focused on Phosphorothioation of the DNA Backbone , 2009, PloS one.

[22]  David G. E. Smith,et al.  Escherichia coli O 157 : H 7 forms attaching and effacing lesions at the terminal rectum of cattle and colonization requires the LEE 4 operon , 2005 .

[23]  V. Parreira,et al.  Phosphoethanolamine substitution in the lipid A of Escherichia coli O157 : H7 and its association with PmrC. , 2006, Microbiology.

[24]  Zixin Deng,et al.  A novel DNA modification by sulphur , 2005, Molecular microbiology.

[25]  E. Brzuszkiewicz,et al.  Genome sequence analyses of two isolates from the recent Escherichia coli outbreak in Germany reveal the emergence of a new pathotype: Entero-Aggregative-Haemorrhagic Escherichia coli (EAHEC) , 2011, Archives of Microbiology.

[26]  David S. Wishart,et al.  PHAST: A Fast Phage Search Tool , 2011, Nucleic Acids Res..

[27]  H. Moon,et al.  Escherichia coli O157:H7 Causes More-Severe Systemic Disease in Suckling Piglets than in Colostrum-Deprived Neonatal Piglets , 2000, Infection and Immunity.

[28]  P S Hayes,et al.  An outbreak of Escherichia coli O157:H7 infections associated with leaf lettuce consumption. , 1998, The Journal of infectious diseases.

[29]  V. Sperandio,et al.  Acyl-Homoserine-Lactone Autoinducer in the Gastrointesinal Tract of Feedlot Cattle and Correlation to Season, E. Coli O157:H7 Prevalence, and Diet , 2009, Current Microbiology.

[30]  T. Croley,et al.  Draft Genome Sequences of Two O104:H21 Escherichia coli Isolates Causing Hemorrhagic Colitis during a 1994 Montana Outbreak Provide Insight into Their Pathogenicity , 2013, Genome Announcements.

[31]  Eduardo N. Taboada,et al.  Genome evolution in major Escherichia coli O157:H7 lineages , 2007, BMC Genomics.

[32]  DNA Backbone Sulfur-Modification Expands Microbial Growth Range under Multiple Stresses by its anti-oxidation function , 2017, Scientific Reports.

[33]  A. von Haeseler,et al.  IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies , 2014, Molecular biology and evolution.

[34]  T. Strobaugh,et al.  Phage insertion in mlrA and variations in rpoS limit curli expression and biofilm formation in Escherichia coli serotype O157: H7. , 2013, Microbiology.

[35]  J. H. Green,et al.  The United States National Prospective Hemolytic Uremic Syndrome Study: microbiologic, serologic, clinical, and epidemiologic findings. , 2001, The Journal of infectious diseases.

[36]  Hong-Yu Ou,et al.  The dnd operon for DNA phosphorothioation modification system in Escherichia coli is located in diverse genomic islands , 2015, BMC Genomics.

[37]  A. Benson,et al.  Octamer-based genome scanning distinguishes a unique subpopulation of Escherichia coli O157:H7 strains in cattle. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[38]  H. Moon,et al.  Comparative Pathogenicity of Escherichia coli O157 and Intimin-Negative Non-O157 Shiga Toxin-Producing E. coli Strains in Neonatal Pigs , 2003, Infection and Immunity.

[39]  C. Parker,et al.  Distinct Acid Resistance and Survival Fitness Displayed by Curli Variants of Enterohemorrhagic Escherichia coli O157:H7 , 2011, Applied and Environmental Microbiology.

[40]  S. Salzberg,et al.  Versatile and open software for comparing large genomes , 2004, Genome Biology.

[41]  L. Gansheroff,et al.  Vaccination of Pregnant Dams with IntiminO157 Protects Suckling Piglets from Escherichia coli O157:H7 Infection , 2002, Infection and Immunity.

[42]  N. Perna,et al.  progressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement , 2010, PloS one.

[43]  T. Whittam,et al.  Clonal relationships among Escherichia coli strains that cause hemorrhagic colitis and infantile diarrhea , 1993, Infection and immunity.

[44]  V. Sperandio,et al.  SdiA Aids Enterohemorrhagic Escherichia coli Carriage by Cattle Fed a Forage or Grain Diet , 2013, Infection and Immunity.

[45]  M. Brandl,et al.  Curli fimbriae are conditionally required in Escherichia coli O157:H7 for initial attachment and biofilm formation. , 2016, Food microbiology.

[46]  C. Gyles.  Shiga toxin-producing Escherichia coli: an overview. , 2007, Journal of animal science.

[47]  W. Melchers,et al.  Characterization of verocytotoxin-producing Escherichia coli O157 isolates from patients with haemolytic uraemic syndrome in Western Europe , 1995, Epidemiology and Infection.

[48]  W. Elkhatib,et al.  Multilocus genotypic characterization of Escherichia coli O157:H7 recovered from food sources , 2014, Epidemiology and Infection.

[49]  T. Whittam,et al.  Variation in virulence among clades of Escherichia coli O157:H7 associated with disease outbreaks , 2008, Proceedings of the National Academy of Sciences.

[50]  Z. Deng,et al.  Site-specific degradation of Streptomyces lividans DNA during electrophoresis in buffers contaminated with ferrous iron. , 1988, Nucleic acids research.

[51]  Eric J. Alm,et al.  DNA phosphorothioation is widespread and quantized in bacterial genomes , 2011, Proceedings of the National Academy of Sciences.

[52]  Genome amplification and promoter mutation expand the range of csgD-dependent biofilm responses in an STEC population. , 2017, Microbiology.

[53]  L. Harrison,et al.  Genotyping Primers for Fully Automated Multilocus Variable-Number Tandem Repeat Analysis of Escherichia coli O157:H7 , 2004, Journal of Clinical Microbiology.

[54]  A. Benson,et al.  Identification of Common Subpopulations of Non-Sorbitol-Fermenting, β-Glucuronidase-Negative Escherichia coli O157:H7 from Bovine Production Environments and Human Clinical Samples , 2004, Applied and Environmental Microbiology.

[55]  V. Sharma,et al.  Role of curli and plant cultivation conditions on Escherichia coli O157:H7 internalization into spinach grown on hydroponics and in soil. , 2014, International journal of food microbiology.

[56]  T. Besser,et al.  Multilocus Genotype Analysis of Escherichia coli O157 Isolates from Australia and the United States Provides Evidence of Geographic Divergence , 2013, Applied and Environmental Microbiology.

[57]  A. Benson,et al.  Genotypic Characterization and Prevalence of Virulence Factors among Canadian Escherichia coli O157:H7 Strains , 2008, Applied and Environmental Microbiology.

[58]  R. Oakes,et al.  Hemolytic uremic syndrome; pathogenesis, treatment, and outcome , 2005, Current opinion in pediatrics.

[59]  L. Beutin,et al.  Derivation of Escherichia coli O157:H7 from Its O55:H7 Precursor , 2010, PloS one.

[60]  T. Besser,et al.  The control of VTEC in the animal reservoir. , 2001, International journal of food microbiology.

[61]  Eduardo N. Taboada,et al.  In silico genomic analyses reveal three distinct lineages of Escherichia coli O157:H7, one of which is associated with hyper-virulence , 2009, BMC Genomics.

[62]  M. Hattori,et al.  Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12. , 2001, DNA research : an international journal for rapid publication of reports on genes and genomes.

[63]  David L. Swerdlow,et al.  Epidemiology of Escherichia coli O157:H7 Outbreaks, United States, 1982–2002 , 2005, Emerging infectious diseases.

[64]  V. Sharma,et al.  Role of curli and cellulose expression in adherence of Escherichia coli O157:H7 to spinach leaves. , 2012, Foodborne pathogens and disease.

[65]  David G. E. Smith,et al.  Edinburgh Research Explorer An investigation of the expression and adhesin function of H7 flagella in the interaction of Escherichia coli O157 : H7 with bovine intestinal epithelium , 2022 .

[66]  V. Sharma,et al.  Contributions of EspA Filaments and Curli Fimbriae in Cellular Adherence and Biofilm Formation of Enterohemorrhagic Escherichia coli O157:H7 , 2016, PloS one.

[67]  Correlation of enterohemorrhagic Escherichia coli O157 prevalence in feces, hides, and carcasses of beef cattle during processing. , 2000 .

[68]  Minjun Yang,et al.  Genome analysis and in vivo virulence of porcine extraintestinal pathogenic Escherichia coli strain PCN033 , 2015, BMC Genomics.

[69]  J. Hadler,et al.  A multistate outbreak of Escherichia coli O157:H7 infections associated with consumption of mesclun lettuce. , 1999, Archives of internal medicine.

[70]  G. Remuzzi,et al.  Hemolytic uremic syndrome. , 2005, Journal of the American Society of Nephrology : JASN.

[71]  A. Bustamante,et al.  Molecular characterization of Verocytotoxigenic Escherichia coli O157:H7 isolates from Argentina by Multiple-Loci VNTR Analysis (MLVA). , 2009 .

[72]  Jianzhong Shen,et al.  Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. , 2015, The Lancet. Infectious diseases.

[73]  Nicole T. Perna,et al.  Molecular Evolution of a Pathogenicity Island from Enterohemorrhagic Escherichia coli O157:H7 , 1998, Infection and Immunity.

[74]  Dale D. Hancock,et al.  Carriage of stx2a Differentiates Clinical and Bovine-Biased Strains of Escherichia coli O157 , 2012, PloS one.

[75]  Masahira Hattori,et al.  Comparative genomics reveal the mechanism of the parallel evolution of O157 and non-O157 enterohemorrhagic Escherichia coli , 2009, Proceedings of the National Academy of Sciences.

[76]  T. Whittam,et al.  Probing genomic diversity and evolution of Escherichia coli O157 by single nucleotide polymorphisms. , 2006, Genome research.

[77]  L. Teel,et al.  Pathogenesis of Shiga-toxin producing escherichia coli. , 2012, Current topics in microbiology and immunology.

[78]  Jianzhong Shen,et al.  Novel Plasmid-Mediated Colistin Resistance Gene mcr-3 in Escherichia coli , 2017, mBio.

[79]  V. Sharma,et al.  Escherichia coli O157:H7 Lacking the qseBC-Encoded Quorum-Sensing System Outcompetes the Parental Strain in Colonization of Cattle Intestines , 2014, Applied and Environmental Microbiology.

[80]  A. von Haeseler,et al.  UFBoot2: Improving the Ultrafast Bootstrap Approximation , 2017, bioRxiv.

[81]  V. Sharma,et al.  Hha Controls Escherichia coli O157:H7 Biofilm Formation by Differential Regulation of Global Transcriptional Regulators FlhDC and CsgD , 2013, Applied and Environmental Microbiology.

[82]  Phillip I. Tarr,et al.  Escherichia coli O157:H7 Shiga Toxin-Encoding Bacteriophages: Integrations, Excisions, Truncations, and Evolutionary Implications , 2003, Journal of bacteriology.

[83]  Zhengwei Zhu,et al.  CD-HIT: accelerated for clustering the next-generation sequencing data , 2012, Bioinform..

[84]  N. W. Davis,et al.  Genome sequence of enterohaemorrhagic Escherichia coli O157:H7 , 2001, Nature.

[85]  C. Southward,et al.  Genomic Instability in Regions Adjacent to a Highly Conserved pch Prophage in Escherichia coli O157:H7 Generates Diversity in Expression Patterns of the LEE Pathogenicity Island , 2009, Journal of bacteriology.

[86]  D. Fleming,et al.  A prolonged outbreak of Escherichia coli O157:H7 infections caused by commercially distributed raw milk. , 1997, The Journal of infectious diseases.

[87]  D. Alt,et al.  Disruption of rcsB by a duplicated sequence in a curli-producing Escherichia coli O157:H7 results in differential gene expression in relation to biofilm formation, stress responses and metabolism , 2017, BMC Microbiology.

[88]  Adam Godzik,et al.  Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences , 2006, Bioinform..

[89]  Ying Xu,et al.  Prediction of pathogenicity islands in Enterohemorrhagic Escherichia coli O157:H7 using genomic barcodes , 2010, FEBS letters.

[90]  Dale D. Hancock,et al.  Evaluation of Pulsed-Field Gel Electrophoresis as a Tool for Determining the Degree of Genetic Relatedness between Strains of Escherichia coli O157:H7 , 2003, Journal of Clinical Microbiology.

[91]  Z. Deng,et al.  A Novel Target of IscS in Escherichia coli: Participating in DNA Phosphorothioation , 2012, PloS one.

[92]  J. Keen,et al.  Variations in the csgD Promoter of Escherichia coli O157:H7 Associated with Increased Virulence in Mice and Increased Invasion of HEp-2 Cells , 2002, Infection and Immunity.

[93]  T. Besser,et al.  Correlation between geographic distance and genetic similarity in an international collection of bovine faecal Escherichia coli O157[ratio ]H7 isolates , 2003, Epidemiology and Infection.

[94]  S. Huynh,et al.  Distinct Transcriptional Profiles and Phenotypes Exhibited by Escherichia coli O157:H7 Isolates Related to the 2006 Spinach-Associated Outbreak , 2011, Applied and Environmental Microbiology.

[95]  David S. Wishart,et al.  PHASTER: a better, faster version of the PHAST phage search tool , 2016, Nucleic Acids Res..

[96]  Martyn D Kirk,et al.  Global incidence of human Shiga toxin-producing Escherichia coli infections and deaths: a systematic review and knowledge synthesis. , 2014, Foodborne pathogens and disease.

[97]  Z. Deng,et al.  Analysis of a genomic island housing genes for DNA S‐modification system in Streptomyces lividans 66 and its counterparts in other distantly related bacteria , 2007, Molecular microbiology.

[98]  A. O’Brien,et al.  The Established Intimin Receptor Tir and the Putative Eucaryotic Intimin Receptors Nucleolin and β1 Integrin Localize at or near the Site of Enterohemorrhagic Escherichia coli O157:H7 Adherence to Enterocytes In Vivo , 2006, Infection and Immunity.

[99]  T. Whittam,et al.  Genotypic and phenotypic changes in the emergence of Escherichia coli O157:H7. , 1998, The Journal of infectious diseases.

[100]  K. Kurokawa,et al.  Diversification of Escherichia coli genomes: are bacteriophages the major contributors? , 2001, Trends in microbiology.

[101]  Zixin Deng,et al.  A novel host-specific restriction system associated with DNA backbone S-modification in Salmonella , 2010, Nucleic acids research.

[102]  J. Bono,et al.  Lineage and Genogroup-Defining Single Nucleotide Polymorphisms of Escherichia coli O157:H7 , 2013, Applied and Environmental Microbiology.