- Genome Report
- Open Access
Genome analyses of blaNDM-4 carrying ST 315 Escherichia coli isolate from sewage water of one of the Indian hospitals
© The Author(s) 2018
- Received: 14 April 2018
- Accepted: 19 May 2018
- Published: 24 May 2018
Emergence of carbapenem resistant Escherichia coli pathovars and their environmental dissemination are alarming problems. E. coli isolated from sewage water of hospital setting conferred a high resistance towards β-lactams, particularly towards carbapenem. This prompted us to perform whole genome sequence analysis to investigate the antimicrobial determinants, pathogenicity status and mobile genetic elements associated with resistance genes.
To the best of our knowledge this is the first report of ST 315 carrying NDM-4 from India. The genome analysis has revealed the unknown characteristics associated with this sequence type (ST 315) like resistance and virulence factors. Based on virulence markers, its pathotype was identified as ExPEC. Furthermore, a mobile plasmid with multiple β-lactamases genes and clinically relevant resistance markers was detected. Phylogenetic analysis of Inc F plasmids sequences carrying ESBLs and NDM variants, revealed un-relatedness in these plasmids due to their varying size and backbone sequences.
Presence of carbapenem resistant E. coli ST 315 with high level antibiotic resistance, near hospital environment is an alarming situation in context to its spread. WGS based analyses have provided details on virulence and resistance status which could overcome the lack of information available on ST 315, globally. This could further help in its quick detection and control in clinical settings.
- ST 315
- Hospital setting
Escherichia coli share a commensal relationship with humans and animals. The extensive acquisition of virulence genes has potentiated E. coli to become pathogenic . The E. coli pathotypes are identified on the basis of virulence determinants present in the genome. A pathotype, ExPEC (extra intestinal pathogenic E. coli) has been reported in extra intestinal, neonatal meningitis and septicaemia infections. Over the years ExPECs are being increasingly recognised for plasmid-mediated carriage of extended spectrum β-lactamases and carbapenemases (metallo β-lactamases, MBLs) . The emergence of a carbapenemase, New Delhi-Metallo β-lactamase (NDM-1) has conferred resistance to last resort β-lactams which has further made the management of ExPECs difficult . A single amino acid variation (Met154Leu) in NDM-1 has resulted in the emergence of a novel NDM-4 which has extended and increased hydrolytic activity towards β-lactams, especially towards carbapenem .
The bacterial isolate AK-1 found in hospital sewage water was subjected to antibiotic susceptibility testing which revealed an exorbitantly elevated MIC values against β-lactams . This unusual resistance in AK-1 strain intrigued us to further explore other genetically predisposed features through whole genome sequencing.
E. coli isolation and characterisation
The isolate for whole genome sequence analysis was collected and identified as reported in our previous published study . It was identified and characterised as NDM-4 producing E. coli strain.
Whole genome sequencing, annotation and analysis
The bacterial DNA was isolated from the AK-1 strain by Qiagen’s QIAamp DNA mini kit and GE SimpliNano UV–Vis Spectrophotometer was used to measure the concentration and purity of the DNA. The genomic DNA was subjected to whole genome sequencing on Illumina NextSeq 500 platform using the paired-end 2 × 150 nt sequencing protocol. The raw sequence data was further analysed by FastQC tool for quality control purposes . SPAdes version 3.10.1 was used to create denovo assembly with genome coverage of 266.181× . Genome annotation was performed by NCBI Prokaryotic Genome Annotation Pipeline using Best-placed reference protein set and GeneMarkS+ methodology.
StringMLST tool  was used to determine sequence type and multi-locus sequence typing database of E. coli at PubMLST (http://ukmirror1.pubmlst.org/databases.shtml)  was used as reference. The serotype of the strain (fliC, wxy and wzx genes) was determined by SerotypeFinder . Resistance genes were identified using CARD and ARDB databases [11, 12]. Virulence factors were determined using the combination of UniProt and Virulence finder database (VFDB) . Further genome analyses was done using ISfinder  and RAC . Phylogenetic analysis was done on MEGA (version 7) using plasmid sequences retrieved from GenBank.
Escherichia coli ST 315 genomic DNA was extracted from a single colony and this strain was maintained. The 16s rRNA gene from the draft genome was confirmed to check for contamination. High quality of reads were confirmed by CLC bio Genomic Workbench version 9 (CGWB) and selected for assembly.
General features of AK-1 genome
6, 2, 5 (5S, 16S, 23S)
6, 1 (5S, 16S)
1, 5 (16S, 23S)
Pseudo genes (total)
AK-1 was found to carry a plasmid of 155,678 bp (Additional file 1: Figure S1). pMLST identified plasmid belonging to IncF group having alleles, repFIB and repFII.
The genetic environment of blaTEM1b was identified to be flanked by one truncated copy of IS26. The blaTEM1b gene was associated with partial Tn2, it was present distally to tnpR-encoding gene and proximally to an IR (Fig. 1b) .
A truncated Tn2 transposon unit was found to be associated with blaCTX-M-15. It was revealed that both DRL and DRR were bracketing this unit where ISEcp1 was distally located to blaCTX-M-15 gene (Fig. 1c). This arrangement is very common and ISEcp1 is reported to mediate the mobilization of blaCTX-M-15 [20, 21].
Other previously reported resistance genes and associated mobile genetic elements were found on AK-1 plasmid. (i) Class I integron carrying qacE∆ 1, sul 2, dfrA12 and aadA2 (Fig. 1d) , (ii) rmt B was found to be associated with partial sequence of IS 26 (Fig. 1e) and upstream to blaTEM1b (not shown) , (iii) sul1 gene upstream to strA and strB genes which are bracketed by Tn 5903 , (iv) tet (A) efflux protein and its regulator tet R (A) associated with Tn1721 mobile element .
AK-1 chromosome carries wide range of resistance markers towards major classes of antibiotics like flouroquinoles, macrolides, aminoglycosides, tetracycline, trimethoprim isoniazid, triclosan, elfamycin and β-lactams (Additional file 3: Table S2). Extensive numbers of genes conferring resistance towards β-lactams were found; four types of penicillin binding proteins and class C β-lactamases, blaCFE 1 and blaPEDO 2. Accumulative effect of these genes explains the high level of phenotypic resistance towards β-lactams in AK-1 .
Type 1 fimbriae NSBV01000003 (15138–21501) sfm fimbriae NSBV01000011 (7740–14015) yad fimbriae NSBV01000005 (66168–72121) MAT(ECP) fimbriae NSBV01000007 (119869–126709) FdeC NSBV01000007 (106372–110622) Curli NSBV01000004 (93069–97942)
ChuA NSBV01000028 (31777–40728) Enterobactin NSBV01000021 (27676–46582) Aerobactin NSBV01000005 (59212–52929)
HlyA NSBV01000011 (54096–55271) HlyE NSBV01000026 (46866–47777)
TypeIII secretion system (ETT2) NSBV01000002 (15607–31153) TypeVI secretion system NSBV01000032 (5165–28659)
ColV NSBV01000002 (84476–84964) ibeB NSBV01000021 (57575–58957) K1capsule NSBV01000002 (150300–158286) ompA NSBV01000004 (166414–167466)
Haemolytic toxins were also identified in AK-1 (i) Hemolysin α is associated with ExPEC virulence and attacks immune cell , (ii) membrane pore-forming toxin HlyE lyses mammalian cells and erythrocytes .
ExPECs cope up with low iron availability by secreting siderophores which retrieves sequestered iron from host proteins . AK-1 was identified to harbour aerobactin, enterobactin and chuA siderophores. Proectins/invasins like ibeB, ompA and K1 capsule have been reported in invasion of brain microvasular endothelial cells .
The genetic context of NDM-4 in AK-1 plasmid was similar to plasmids (pM109 FII, pMC-NDM, pGUE NDM, pCRKP-2297, pCRKP-1215, pM214 FII, and pNDM5-IBAC) carrying NDM variants. The genetic environment for blaNDM remained conserved in these plasmids. This suggests that alteration in bla gene originated new NDM variants. Furthermore, genetic context of blaTEM1b and blaCTX-M-15 were almost similar in these plasmids.
Presence of plasmid harbouring blaNDM-4 and other β-lactamase genes near the hospital setting environment is a serious concern in context to its circulation and spread in hospital settings. It is the first time NDM-4 producing ST 315 E. coli was detected in India. Moreover, no genome based information on ST 315 strain is yet available. Lack of information on mechanism of virulence, transmission sources and other genetic characteristics have become an impasse for alleviation of ExPECs infections. AK-1 genome based virulence profile provides cause of serious infections by ST 315 ExPEC, a common microflora of healthy individuals. Genome informed virulence and resistance mechanisms will definitely help in identification of this ExPECs in healthcare settings and controlling the clonal spread of carbapenamase carrying E. coli pathovars.
AZB, analysed the data and wrote the first draft of the MS. AUK, provided the NGS data, helped in interpretation of the data and checked first draft. Both authors read and approved the final manuscript.
We acknowledge the support of internal funds of Biotechnology Department.
The authors declare that they have no competing interests.
Availability of data and materials
The whole genome sequence project of AK-1 is deposited to GenBank under the accession id NSBV01000000 and consists of sequences NSBV01000001–NSBV01000129. The genome draft has also been published .
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Croxen MA, Finlay BB. Molecular mechanisms of Escherichia coli pathogenicity. Nat Rev Microbiol. 2010;8(1):26–38.View ArticlePubMedGoogle Scholar
- Pitout JDD. Extraintestinal pathogenic Escherichia coli: a combination of virulence with antibiotic resistance. Front Microbiol. 2012;3:9.View ArticlePubMedPubMed CentralGoogle Scholar
- Walsh TR, et al. Metallo-β-lactamases: the quiet before the storm? Clin Microbiol Rev. 2005;18(2):306–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Nordmann P, Boulanger AE, Poirel L. NDM-4 metallo-β-lactamase with increased carbapenemase activity from Escherichia coli. Antimicrob Agents Chemother. 2012;56(4):2184–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Khan AU, Parvez S. Detection of blaNDM-4 in Escherichia coli from hospital sewage. J Med Microbiol. 2014;63(10):1404–6.View ArticlePubMedGoogle Scholar
- Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc. Accessed Sept 2016.
- Bankevich A, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Gupta A, Jordan IK, Rishishwar L. stringMLST: a fast k-mer based tool for multilocus sequence typing. Bioinformatics. 2016;33(1):119–21.View ArticlePubMedGoogle Scholar
- Wirth T, et al. Sex and virulence in Escherichia coli: an evolutionary perspective. Mol Microbiol. 2006;60(5):1136–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Joensen KG, et al. Rapid and easy in silico serotyping of Escherichia coli using whole genome sequencing (WGS) data. J Clin Microbiol. 2015. https://doi.org/10.1128/JCM-00008-15.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu B, Pop M. ARDB—antibiotic resistance genes database. Nucleic Acids Res. 2008;37(suppl_1):D443–7.PubMedPubMed CentralGoogle Scholar
- McArthur AG, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57(7):3348–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Kleinheinz KA, Joensen KG, Larsen MV. Applying the ResFinder and VirulenceFinder web-services for easy identification of acquired antibiotic resistance and E. coli virulence genes in bacteriophage and prophage nucleotide sequences. Bacteriophage. 2014;4(2):e27943.View ArticlePubMedPubMed CentralGoogle Scholar
- Siguier P, et al. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006;34(suppl_1):D32–6.View ArticlePubMedGoogle Scholar
- Tsafnat G, Copty J, Partridge SR. RAC: repository of antibiotic resistance cassettes. Database. 2011. https://doi.org/10.1093/database/bar054.PubMedPubMed CentralView ArticleGoogle Scholar
- Peirano G, et al. Molecular epidemiology over an 11-year period (2000 to 2010) of extended-spectrum β-lactamase-producing Escherichia coli causing bacteremia in a centralized Canadian region. J Clin Microbiol. 2012;50(2):294–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Ewers C, et al. Extended-spectrum β-lactamase-producing and AmpC-producing Escherichia coli from livestock and companion animals, and their putative impact on public health: a global perspective. Clin Microbiol Infect. 2012;18(7):646–55.View ArticlePubMedGoogle Scholar
- Poirel L, et al. Genetic features of blaNDM-1-positive Enterobacteriaceae. Antimicrob Agents Chemother. 2011;55(11):5403–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, et al. Characterization of the genetic environment of blaESBL genes, integrons and toxin-antitoxin systems identified on large transferrable plasmids in multi-drug resistant Escherichia coli. Front Microbiol. 2015;5:716.View ArticlePubMedPubMed CentralGoogle Scholar
- Boyd DA, et al. Complete nucleotide sequence of a 92-kilobase plasmid harboring the CTX-M-15 extended-spectrum beta-lactamase involved in an outbreak in long-term-care facilities in Toronto, Canada. Antimicrob Agents Chemother. 2004;48(10):3758–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Poirel L, Decousser J-W, Nordmann P. Insertion sequence ISEcp1B is involved in expression and mobilization of a blaCTX-M β-lactamase gene. Antimicrob Agents Chemother. 2003;47(9):2938–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Gillings MR, Paulsen IT, Tetu SG. Genomics and the evolution of antibiotic resistance. Ann N Y Acad Sci. 2016;1388(1):92–107.View ArticlePubMedGoogle Scholar
- Yu FY, et al. High prevalence of plasmid-mediated 16S rRNA methylase gene rmtB among Escherichia coli clinical isolates from a Chinese teaching hospital. BMC Infect Dis. 2010;10(1):184.View ArticlePubMedPubMed CentralGoogle Scholar
- Sundin GW, Bender CL. Dissemination of the strA–strB streptomycin-resistance genes among commensal and pathogenic bacteria from humans, animals, and plants. Mol Ecol. 1996;5(1):133–43.View ArticlePubMedGoogle Scholar
- Wang W, et al. High-level tetracycline resistance mediated by efflux pumps Tet (A) and Tet (A)-1 with two start codons. J Med Microbiol. 2014;63(11):1454–9.View ArticlePubMedGoogle Scholar
- Mellata M. Human and avian extraintestinal pathogenic Escherichia coli: infections, zoonotic risks, and antibiotic resistance trends. Foodborne Pathog Dis. 2013;10(11):916–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Ideses D, et al. A degenerate type III secretion system from septicemic Escherichia coli contributes to pathogenesis. J Bacteriol. 2005;187(23):8164–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang S, et al. Escherichia coli type III secretion system 2 ATPase EivC is involved in the motility and virulence of avian pathogenic Escherichia coli. Front Microbiol. 2016;7:1387.PubMedPubMed CentralGoogle Scholar
- Bingle LEH, Bailey CM, Pallen MJ. Type VI secretion: a beginner’s guide. Curr Opin Microbiol. 2008;11(1):3–8.View ArticlePubMedGoogle Scholar
- Antao EM, Wieler LH, Ewers C. Adhesive threads of extraintestinal pathogenic Escherichia coli. Gut Pathog. 2009;1(1):22.View ArticlePubMedPubMed CentralGoogle Scholar
- Badouraly RA, et al. Escherichia coli K-12 possesses multiple cryptic but functional chaperone–usher fimbriae with distinct surface specificities. Environ Microbiol. 2010;12(7):1957–77.View ArticlePubMedGoogle Scholar
- Rendon MA, et al. Commensal and pathogenic Escherichia coli use a common pilus adherence factor for epithelial cell colonization. Proc Natl Acad Sci. 2007;104(25):10637–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Nesta B, et al. FdeC, a novel broadly conserved Escherichia coli adhesin eliciting protection against urinary tract infections. MBio. 2012;3(2):e00010–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas S, Holland IB, Schmitt L. The type 1 secretion pathway—the hemolysin system and beyond. Biochim Biophys Acta. 2013;1843(8):1629–41.View ArticlePubMedGoogle Scholar
- Hunt S, Green J, Artymiuk PJ. Hemolysin E (HlyE, ClyA, SheA) and related toxins. In: Proteins membrane binding and pore formation. 2010:116–26.Google Scholar
- Garenaux A, Caza M, Dozois CM. The Ins and Outs of siderophore mediated iron uptake by extra-intestinal pathogenic Escherichia coli. Vet Microbiol. 2011;153(1):89–98.View ArticlePubMedGoogle Scholar
- Wang S, et al. IbeB is involved in the invasion and pathogenicity of avian pathogenic Escherichia coli. Vet Microbiol. 2012;159(3):411–9.View ArticlePubMedGoogle Scholar
- Khan AU, Beg AZ, Verma PK. Draft genome sequence of the first NDM-4-producing Escherichia coli strain (AK1), isolated from sewage water of a North Indian hospital. Genome Announc. 2017;5(50):e01366-17.View ArticlePubMedPubMed CentralGoogle Scholar