Draft genome sequence of non-shiga toxin-producing Escherichia coli O157 NCCP15738

Construction of a Lectin-Glycan Interaction Network from Enterohemorrhagic Escherichia coli Strains by Multi-omics Analysis

Enterohemorrhagic Escherichia coli (EHEC) causes hemorrhagic colitis and hemolytic uremic syndrome. EHEC infection begins with bacterial adherence to the host intestine via lectin-like adhesins that bind to the intestinal wall. However, EHEC-related lectin–glycan interactions (LGIs) remain unknown. Here, we conducted a genome-wide investigation of putative adhesins to construct an LGI network. We performed microarray-based transcriptomic and proteomic analyses with E. coli EDL933. Using PSORTb-based analysis, potential outer-membrane-embedded adhesins were predicted from the annotated genes of 318 strains. Predicted proteins were classified using TMHMM v2.0, SignalP v5.0, and LipoP v1.0. Functional and protein–protein interaction analyses were performed using InterProScan and String databases, respectively. Structural information of lectin candidate proteins was predicted using Iterative Threading ASSEmbly Refinement (I-TASSER) and Spatial Epitope Prediction of Protein Antigens (SEPPA) tools based on 3D structure and B-cell epitopes. Pathway analysis returned 42,227 Gene Ontology terms; we then selected 2585 lectin candidate proteins by multi-omics analysis and performed homology modeling and B-cell epitope analysis. We predicted a total of 24,400 outer-membrane-embedded proteins from the genome of 318 strains and integrated multi-omics information into the genomic information of the proteins. Our integrated multi-omics data will provide a useful resource for the construction of LGI networks of E. coli.

Comparative genomics reveals differences in mobile virulence genes of Escherichia coli O103 pathotypes of bovine fecal origin

Escherichia coli O103, harbored in the hindgut and shed in the feces of cattle, can be enterohemorrhagic (EHEC), enteropathogenic (EPEC), or putative non-pathotype. The genetic diversity particularly that of virulence gene profiles within O103 serogroup is likely to be broad, considering the wide range in severity of illness. However, virulence descriptions of the E. coli O103 strains isolated from cattle feces have been primarily limited to major genes, such as Shiga toxin and intimin genes. Less is known about the frequency at which other virulence genes exist or about genes associated with the mobile genetic elements of E. coli O103 pathotypes. Our objective was to utilize whole genome sequencing (WGS) to identify and compare major and putative virulence genes of EHEC O103 (positive for Shiga toxin gene, stx1, and intimin gene, eae; n = 43), EPEC O103 (negative for stx1 and positive for eae; n = 13) and putative non-pathotype O103 strains (negative for stx and eae; n = 13) isolated from cattle feces. Six strains of EHEC O103 from human clinical cases were also included. All bovine EHEC strains (43/43) and a majority of EPEC (12/13) and putative non-pathotype strains (12/13) were O103:H2 serotype. Both bovine and human EHEC strains had significantly larger average genome sizes (P < 0.0001) and were positive for a higher number of adherence and toxin-based virulence genes and genes on mobile elements (prophages, transposable elements, and plasmids) than EPEC or putative non-pathotype strains. The genome size of the three pathotypes positively correlated (R² = 0.7) with the number of genes carried on mobile genetic elements. Bovine strains clustered phylogenetically by pathotypes, which differed in several key virulence genes. The diversity of E. coli O103 pathotypes shed in cattle feces is likely reflective of the acquisition or loss of virulence genes carried on mobile genetic elements.

Comparative genomic analysis of Shiga toxin-producing and non-Shiga toxin-producing Escherichia coli O157 isolated from outbreaks in Korea

Background

The Shiga toxin-producing Escherichia coli (STEC) O157 strain NCCP15739 and non-STEC O157 strain NCCP15738 were isolated from outbreaks in Korea. We characterized NCCP15739 and NCCP15738 by genome sequencing and a comparative genomic analysis using two additional strains, E. coli K-12 substr. MG1655 and O157:H7 EDL933.

Results

Using the Illumina HiSeq 2000 platform and the RAST server, the whole genomes of NCCP15739 and NCCP15738 were obtained and annotated. NCCP15739 and NCCP15738 clustered with different E. coli strains based on a whole-genome phylogeny and multi-locus sequence typing analysis. Functional annotation clustering indicated enrichment for virulence plasmid and hemolysis-related genes in NCCP15739 and conjugation- and flagellum-related genes in NCCP15738. Defense mechanism- and pathogenicity-related pathways were enriched in NCCP15739 and pathways related to the assimilation of energy sources were enriched in NCCP15738. We identified 66 and 18 virulence factors from the NCCP15739 and NCCP15738 genome, respectively. Five and eight antibiotic resistance genes were identified in the NCCP15739 and NCCP15738 genomes, respectively. Based on a comparative analysis of phage-associated regions, NCCP15739 and NCCP15738 had specific prophages. The prophages in NCCP15739 carried virulence factors, but those in NCCP15738 did not, and no antibiotic resistance genes were found in the phage-associated regions.

Conclusions

Our whole-genome sequencing and comparative genomic analysis revealed that NCCP15739 and NCCP15738 have specific genes and pathways. NCCP15739 had more genes (410), virulence factors (48), and phage-related regions (11) than NCCP15738. However, NCCP15738 had three more antibiotic resistance genes than NCCP15739. These differences may explain differences in pathogenicity and biological characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13099-017-0156-2) contains supplementary material, which is available to authorized users.