INTRODUCTION
Escherichia coli is a facultative anaerobic bacterium which is commonly spread on biosphere. E. coli normally colonizes the lower intestine of animals and humans (1). However, Some of the serotypes such as Enterohemorrhagic E. coli (EHEC), Enterotoxigenic E. coli (ETEC), Enteropathogenic E. coli (EPEC) and Shiga toxin-producing E. coli (STEC) can cause foodborne illnesses in people.
E. coli K_EC180 was isolated from swine feces that were collected from a livestock farm in Haenamgun, Jeollanam-do, Korea. E. coli K_EC180 was streaked to Luria-Bertani (LB) agar and incubated at 37℃ for 24 h. The suspected colony in LB agar was inoculated into LB broth and incubated at 37℃ for 24 h. To analyze the complete genome, the E. coli K_EC180 genome was sequenced by PacBio RS II (Pacific Biosciences, Menlo Park, CA, USA) at Insilicogen (Yongin, Korea) and Illumina NextSeq 500 (Illumina, San Diego, CA, USA) platform at LabGenomics (Seongnam, Korea). The genomic DNA of E. coli K_EC180 for PacBio and Illumina sequencing was extracted using the MagAttract HMW DNA Kit (QIAGEN), and NucleoSpin® Microbial DNA kit (TAKARA) according to the manufacturer’s instructions. Library preparation was conducted using SMRTbell™ Template Prep Kit 1.0 for Pacbio (Pacific Biosciences) and TruSeq DNA Sample Preparation Kit for Illumina (Illumina) according to the manufacturer’s instructions. PacBio sequencing yielded 1,131,537,370 base pairs and 145,423 long reads after filtering, and 9, 199, 306 paired-end reads with 1,389,095,206 bp were obtained with Illumina sequencing. De novo assembly was conducted using the hierarchical genome assembly process (HGAP v2.3.0) workflow (Chin et al., 2013) and polished using Quiver. Subsequently, Illumina NextSeq reads were aligned to the PacBio RSII assembly using Burrows-Wheeler Aligner (BWA)-MEM v0.7.17-r1188, and the errors were corrected by using Pilon version 1.23 (2, 3). The quality of genome assembly and the validaty of the final genome were assessed using Quality Assessment Tool for Genome Assemblies (QUAST) v5.0.2 and Benchmarking Universal Single-Copy Orthologs (BUSCO) v3.0.2 (4, 5). Open reading frames (ORFs) and RNA genes of E. coli K_EC180 were predicted and functionally annotated through rapid prokaryotic genome annotation (PROKKA) v1.14.5 (6) and Rapid Annotation using Subsystem Technology (RAST) v2.0 (7). The functional categorization and classification of all predicted ORFs were conducted using the RAST server-based SEED viewer and Clusters of Orthologous Groups (COG) – based EggNOG. The putative virulence factors and Antimicrobial resistance were described using BLAST according to the Virulence Factor Database (VFDB) (8). The whole genome of E. coli K_EC180 is composed of one circular chromosome (5,017,281 bp) with 50.4% of G+C content, 4, 935 of coding sequence (CDS), 88 of tRNA, and 22 of rRNA genes.
The complete genome of E. coli K_EC180 contains the toxin genes encoding shiga-like toxin (stx2e subunit A and stx2e subunit B), which may cause diseases in humans by damaging small blood vessels in places such as the digestive tract, kidneys and central nervous system (9, 10). E. coli K-EC180 also possessed essC, escV, escR, escS, escV, and escJ genes which involved in a type III secretion system. In addition, there were fim (A to H) genes encoding Type I fimbriae. We summarized the general properties of the E. coli K_EC180’s complete genome in the Fig. 1 and Table 1.
Fig. 1. Genome map of Escherichia coli K_EC180. The outer circle denotes the locations of all annotated ORFs, and the inner circle with the red denotes GC content. Pink, and green peaks denote GC skew. The orange arrows denote rRNAs, and the sky blue arrows denote the tRNA operons. All annotated ORFs are colored differently based on the COG assignments. ORFs, open reading frames; G, guanine; C, cytosine; COG, clusters of orthologous groups.
Table 1. Genome features of Escherichia coli K_EC180
DATA AVAILABILITY
The complete genome sequences of E. coli K_EC180 were deposited in GeneBank under the accession numbers CP062203. The BioSample accession number is SAMN16277032, and BioProject accession number is PRJNA666028.
References
- Guevarra RB, Lee JH, Lee SH, Seok MJ, Kim DW, Kang BN, et al. Piglet gut microbial shifts early in life: causes and effects. J Anim Sci Biotechnol. 2019;10:1. https://doi.org/10.1186/s40104-018-0308-3
- Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv. https://arxiv.org/abs/1303.3997
- Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLOS ONE. 2014;9:e112963. https://doi.org/10.1371/journal.pone.0112963
- Simao FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31:3210-2. https://doi.org/10.1093/bioinformatics/btv351
- Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29:1072-5. https://doi.org/10.1093/bioinformatics/btt086
- Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068-9. https://doi.org/10.1093/bioinformatics/btu153
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75. https://doi.org/10.1186/1471-2164-9-75
- Liu B, Zheng D, Jin Q, Chen L, Yang J. VFDB 2019: a comparative pathogenomic platform with an interactive web interface. Nucl Acids Res. 2019;47:D687-92. https://doi.org/10.1093/nar/gky1080
- Tesh VL, O'Brien AD. The pathogenic mechanisms of Shiga toxin and the Shiga-like toxins. Mol Microbiol. 1991;5:1817-22. https://doi.org/10.1111/j.1365-2958.1991.tb00805.x
- O'Brien AD, Holmes RK. Shiga and Shiga-like toxins. Microbiol Rev. 1987;51:206-20. https://doi.org/10.1128/MR.51.2.206-220.1987