Characterization of an NDM-5 carbapenemase-producing Escherichia coli ST156 isolate from a poultry farm in Zhejiang, China
BMC Microbiologyvolume 19, Article number: 82 (2019)
The emergence of carbapenem-resistant Enterobacteriaceae strains has posed a severe threat to public health in recent years. The mobile elements carrying the New Delhi metallo-β-lactqtamase (NDM) gene have been regarded as the major mechanism leading to the rapid increase of carbapenem-resistant Enterobacteriaceae strains isolated from clinics and animals.
We describe an NDM-5-producing Escherichia coli strain, ECCRA-119 (sequence type 156 [ST156]), isolated from a poultry farm in Zhejiang, China. ECCRA-119 is a multidrug-resistant (MDR) isolate that exhibited resistance to 27 antimicrobial compounds, including imipenem and meropenem, as detected by antimicrobial susceptibility testing (AST). The complete genome sequence of the ECCRA-119 isolate was also obtained using the PacBio RS II platform. Eleven acquired resistance genes were identified in the chromosome; four were detected in plasmid pTB201, while six were detected in plasmid pTB202. Importantly, the carbapenem-resistant gene blaNDM-5 was detected in the IncX3 plasmid pTB203. In addition, seven virulence genes and one metal-resistance gene were also detected. The results of conjugation experiments and the transfer regions identification indicated that the blaNDM-5-harboring plasmid pTB203 could be transferred between E. coli strains.
The results reflected the severe bacterial resistance in a poultry farm in Zhejiang province and increased our understanding of the presence and transmission of the blaNDM-5 gene.
The overuse of antibiotics has led to the emergence of a large number of multidrug-resistant pathogens, which constitutes a serious threat to public health . Imipenem and meropenem are carbapenem antibiotics that have been used as last resorts in the treatment of infections caused by gram-negative bacteria, especially multidrug-resistant gram-negative pathogens . In 2008, a novel carbapenem resistance gene, New Delhi metallo-β-lactamase (NDM), was detected in Klebsiella pneumoniae isolated from a Swedish patient of Indian origin. This gene attracted international attention for the high level of resistance it confers to bacteria against most β-lactams, except aztreonam, and its spread to over 50 countries . The NDM variant NDM-5 was first reported in 2011 in Escherichia coli isolated from a patient in the United Kingdom who had received treatment in India . Subsequently, NDM-5 was reported in many other countries, including India , Algeria , Japan , South Korea , Australia , China , Denmark , Italy , America , Spain , Egypt , France , and New Zealand . In China, many pathogens carrying blaNDM-5 have been isolated from patients [18,19,20,21]. In addition, blaNDM-5 can also be isolated from pigs [22, 23], dairy cows  and vegetables . The complete sequences of blaNDM-5-harboring plasmids have been helpful for the study of the transmission of the blaNDM-5 gene, although not all of these plasmids have been reported.
In this study, we first describe the NDM-5-producing carbapenem-resistant E. coli strain, ECCRA-119, isolated from a layer hen farm in Zhejiang, China. We obtained the complete genome sequence, predicted the possible mechanism of multidrug resistance and assessed the transmission ability of the plasmid harboring blaNDM-5 from the ECCRA-119 isolate. These results increased our understanding of the diversity and complexity of the strains harboring blaNDM-5.
Two hundred nineteen of the samples studied tested positive for E. coli, and E. coli isolates from all of these samples were obtained and characterized by antimicrobial susceptibility testing (AST) using the VITEK® 2 COMPACT system (BioMérieux, France). The highest overall levels of resistance were observed toward ampicillin, with 74.43% of all isolates resistant to this antimicrobial. High rates of resistance were also observed toward trimethoprim (54.34%), with lower levels of resistance observed toward piperacillin (1.83%), amikacin (2.29%), and amoxicillin (0.91%). No strain was determined to be resistant to tigecycline. One hundred eighty isolates (82.2%) were resistant to at least one antimicrobial agent, and 92 isolates (42.01%) were resistant to three or more antimicrobial agents. Of the 219 E. coli isolates, a carbapenem-resistant strain was identified that showed resistance toward ertapenem and imipenem, which is rare in poultry.
The minimum inhibitory concentrations (MICs) of the ECCRA-119 isolate toward different antibiotics are shown in Table 1. The ECCRA-119 isolate was susceptible to colistin (MIC < 0.125 mg/L), polymyxin B (MIC 1 mg/L) and amikacin (MIC ≤4 mg/L), exhibited intermediate resistance toward gentamicin (MIC 8 mg/L), and was resistant to 27 different compounds from 7 antimicrobial classes that are frequently used in medical treatments, food animal feed and animal medicine (Table 1). In particular, this isolate was resistant to two carbapenems, imipenem (MIC 4 mg/L) and meropenem (MIC 8 mg/L). Therefore, we classified the ECCRA-119 isolate as a multidrug-resistant strain (MDR) due to its nonsusceptibility to many antimicrobial agents, including imipenem and meropenem.
Characterization of the genome sequence of strain ECCRA-119
The genome of the ECCRA-119 isolate consisted of a single circular chromosome and three circular plasmids (Table 2, Figs. 1 and 2b). The chromosome sequence of ECCRA-119 was determined to be 4,893,130 bp in length, have a GC content of 50.77% and encodes 5042 proteins that account for 90.96% of the genome. The average depth of coverage was 210.5×, and 22 rRNAs, 87 tRNAs, and 2 CRISPRs were detected. Three plasmids in the ECCRA-119 isolate were identified, pTB201, pTB202 and pTB203. The plasmid pTB201, which is a combination of IncFII- and IncFIB-type plasmid, was determined to be 146,268 bp in length and have an average GC content of 51.35%. The plasmid pTB202, a p0111-IncN-type plasmid, was determined to be 139,629 bp in length and have an average GC content of 49.13%. In addition, the blaNDM-5-harboring plasmid pTB203, an IncX3-type plasmid, was determined to be 46,161 bp in length and have an average GC content of 46.65%. Moreover, the three plasmids were characterized by S1-PFGE (Fig. 1a), the results of which were consistent with the whole genome sequencing analysis. Multilocus sequence typing (MLST) analysis classified E. coli ECCRA-119 as ST156, suggesting that E. coli ST156 strains have the potential to harbor blaNDM-5-like genes.
Twenty-two acquired resistance genes were identified in the ECCRA-119 genome that belong to eight antibiotic resistance categories (Table 3). Among these genes, 11 are located on the chromosome, four on plasmid pTB201, six on plasmid pTB202 and one on plasmid pTB203. In addition, several gene mutations were identified in the quinolone and fluoroquinolone resistance-determining region on the chromosome (Additional file 1: Table S1). Double gyrA mutations (giving rise to the amino acid substitutions S83 L and D87Y), parC mutation (giving rise to the amino acid substitution S80I) and parE mutations (giving rise to the amino acid substitution S458A) were also predicted in the ECRRA-119 isolate.
Seven virulence factors were detected in the whole genome sequence (Additional file 1: Table S2), four in the chromosome and three in the pTB201 plasmid, indicating the potential virulence of the ECCRA-119 isolate. These virulence factors are grouped into five classes (iss, gad, lpfA, iroN, and cma), which are related to serum survival, glutamate decarboxylase, long polar fimbriae, enterobactin siderophore receptor protein, and colicin M, respectively. In addition, one mercury resistance-related gene, merA, was identified on plasmid pTB201 (Additional file 1: Table S3).
Transferability of plasmids
Conjugation assays confirmed that blaNDM-5 could be transferred between E. coli strains, with an observed transfer frequency of (1.39 ± 0.12) × 10− 5. The antibiotic susceptibility testing results showed that the transconjugants, confirmed by PCR and sequencing, were resistant to meropenem (4 mg/L). The transfer regions of the three plasmids of strain ECCRA-119 were successfully identified (Figs. 1 and 2b) by oriTfinder, including the origin of transfer region (oriT), relaxase gene, bacterial type IV secretion system (T4SS) apparatus gene clusters and the type IV coupling protein (T4CP) gene. Plasmid pTB201 was observed to possess an oriT (52,884–52,969 bp in the plasmid), relaxase gene, T4CP and T4SSs, indicating a high potential for self-transferability . Plasmid pTB202 was observed to harbor a relaxase but lacked an oriT, T4CP and/or T4SS, indicating that it is not a mobilizable plasmid . Plasmid pTB203 possess a relaxase gene, T4CP and T4SSs, but lacked a typical oriT sequence, demonstrating its potential to be transferred to other bacteria , with its transfer ability having been confirmed experimentally.
Phylogenetic analysis of strain ECCRA-119 with other E. coli ST156 isolates
MLST analysis classified E. coli strain ECCRA-119 as ST156. Thus, we built a phylogenetic tree to determine its relationship among ST156 E. coli strains based on a SNP analysis (Additional file 1: Figure S1). We identified 52,076 SNPs from the 37 genome sequences available in GenBank. Of these, 17,953 and 34,123 were identified as core and noncore SNPs, respectively. We excluded the noncore SNPs for the further analysis and constructed a phylogenetic tree based on the genome-wide core SNPs. The core genome analysis identified 5 groups (Additional file 1: Figure S1). E. coli strain ECCRA-119 is grouped with the strains 174,900, SCEC020022 and VREC0575, which were isolated from Bangladesh, China and the United Kingdom, respectively. There were 7 group-specific core SNPs in this group. The number of strain-specific SNPs identified in strains ECCRA-119, 174,900, SCEC020022 and VREC0575 was 59, 71, 134 and 160, respectively. Interestingly, most isolates identified from the same region or source are not in the same lineage. Isolates from different countries were observed to be clustered together (strains 157–1949 and SE11). Similarly, strains isolated from different hosts (wild animals, livestock and poultry, and dog) clustered into the same branch (strains MOD1-EC5693, CVM N33633PS and MOD1-EC6498).
Phylogenetic and comparative analysis of pTB203 and other bla NDM-5-harboring IncX3 plasmids
An SNP-based phylogenetic analysis was conducted using the 52 complete sequences of blaNDM-5-harboring IncX3 plasmids available in GenBank (Fig. 2a). Among these sequences, 41 originated from bacterial strains from humans, 1 from a pig, 5 from geese, 1 from a vegetable, 1 from a layer hen, 1 from sewage, and 2 from unknown sources. Our results showed that the IncX3 plasmids have an extensive host range. Among these 52 plasmids, 43 were isolated in China and 33 were from E. coli. Five plasmids from geese became available in January 2019 but were not published. Among these plasmids, 9 published plasmids were selected and constructed by BRIG (Fig. 2b), including pVH1 (vegetable, China, 46,161 bp) , pNDM_MGR194 (human, India, 46,253 bp) , pECNDM101 (swine, China, 46,165 bp) , pEC463-NDM5 (human, China, 46,145 bp) , pBJ114–46 (human, China, 46,161 bp) , pEsco-5256cz (human, Czech, 46,161 bp) , pEc1929 (human, China, 46,164 bp) , pTB203 (layer hen, China, 46,161 bp, in this study), and pZSH6-blaNDM-5 (human, China, 46,161 bp) . The results of BLAST homology analyses showed that these plasmids had more than 99.9% identity and 99.8% query coverage with each other. The comparative analysis of 9 blaNDM-5-harboring IncX3 plasmids (~ 46 kb) revealed that these plasmids are highly similar to each other, possessing the same backbone that includes the IncX3 replication, blaNDM-5 gene and conjugation/type IV secretion components. This result was further confirmed by the comparative analysis of 52 blaNDM-5-harboring IncX3 plasmids, with the exception of pD2-NDM_1_1 (human, South Korea, 79,613 bp) (Additional file 1: Figure S2). The results of our analysis showed that blaNDM-5-harboring IncX3 plasmids with an ~ 46 kb backbone have extensive host adaptability in Enterobacteriaceae.
Complete sequences of plasmids harboring bla NDM variants from China
At present, 24 variant blaNDM sequences are available in GenBank, all of which were aligned by ClustalX (Additional file 1: Figure S3 and S4). These sequences are 813 bp in length, with the exception of blaNDM-18, and only 1–6 SNPs are observed among these sequences. In particular, the blaNDM-5 gene has the closest homology with blaNDM-17, blaNDM-20 and blaNDM-21 (Additional file 1: Figure S4). Relative to blaNDM-5, blaNDM-17, blaNDM-20 and blaNDM-21 contained point mutations at positions 508 (G → A), 809 (G → A), and 205 (G → A), generating amino acid substitutions Glu170Lys, Arg270His, and Gly69Ser, respectively. In China, 13 types of plasmids harbor blaNDM genes with complete sequence are reported in GenBank, including blaNDM-1, blaNDM-4, blaNDM-5, blaNDM-6, blaNDM-7, blaNDM-9, blaNDM-13, blaNDM-14, blaNDM-16, blaNDM-17, blaNDM-19, blaNDM-20, blaNDM-21 (Additional file 2: Table S4, Fig. 3). The blaNDM-1 and blaNDM-5 genes are the most prevalent blaNDM variants in China, with humans being the primary host source. Additionally, seven blaNDM-5-harboring plasmids have been detected in poultry and livestock in China, including 1 plasmid detected from swine in Sichuan in 2016, 1 plasmid detected from a layer hen in Zhejiang in 2017 (in this study), and 5 plasmids detected from geese in Jiangsu in 2018.
Comparative analysis of plasmids pTB201 and pTB202
We compared the plasmids pTB201 and pTB202 with the corresponding homologous plasmids from GenBank via BLAST analyses. The results showed that plasmid pTB201 shares homology with plasmid pSMS35_130 (CP000971), plasmid pJIE186_2 (JX077110) and the p300 iro gene cluster (AY205565) (Fig. 4a); plasmid pTB202 showed homology with plasmid p1079-IncFIB-N (MG825383) and part of plasmid pD90–3 (CP022453) (Fig. 4b). These comparisons revealed that these plasmids do not have full-length matching plasmids in the GenBank database, although they shared backbones with many other plasmids.
The extensive use of antibacterials has led to the emergence of drug resistance as an increasingly serious issue, which poses a great threat to public health. There have been widespread reports of the isolation of multidrug resistant E. coli from hospitals, poultry, livestock, food, and the environment . In this study, we identified a blaNDM-5-harboring E. coli isolate from a layer hen farm in Zhejiang, China, and we obtained detailed data through bioinformatics and experimental analyses. The AST results showed that the ECCRA-119 isolate is resistant to 27 different compounds used as therapeutics and food animal feeding, indicating its strong environmental adaptability under antibiotic selection pressure. There is no doubt that the multidrug resistance of this strain may present a serious risk to clinical and veterinary medicine. Apart from a few cases, acquired antimicrobial resistance genes and genomic mutations can largely explain drug resistance phenotypes. The identification of the antimicrobial resistance genes acquired by this strain show that this isolate may have a broad spectrum of drug resistance. For example, the presence of the fosA3 gene in the chromosome could result in fosfomycin resistance , but further experiments are needed to confirm this possibility.
The iss gene was detected on both the chromosome and on plasmid pTB201, indicating the potential virulence of the ECCRA-119 isolate. The protein encoded by the iss gene is part of an outer membrane protein and is involved in the anti-complement effect of bacteria, possibly enhancing the serum resistance of E. coli and enabling the strain to rapidly proliferate in the host. It is widely believed that the iss gene is closely associated with the virulence of avian E. coli .
The merA gene was detected in plasmid pTB201, which can confer resistance to mercury and increase the viability of the ECCRA-119 isolate. Furthermore, the results suggested that plasmid pTB201 had a high potential for self-transferability. Therefore, it is likely that the mercury resistance of the ECCRA-119 isolate may be transferred to other bacteria . Thus, the ECCRA-119 isolate has strong environmental resilience and a high potential to survival in a complicated breeding environment for a long time.
To the best of our knowledge, this is the first time a blaNDM-5-harboring plasmid has been reported in layer chickens. E. coli ST156 has not been a predominant multidrug-resistant clone observed worldwide in the past, but it is associated with the distribution of blaNDM-1 and blaCTX-M-15 in humans and poultry [37, 38]. The genes mcr-1 and blaNDM-5 have been reported to be detected in E. coli ST156 from Muscovy duck in China . E. coli ST156 has spread to many countries and can be isolated from many types of hosts, suggesting that E. coli ST156 has the potential to play an important role in the transmission of the blaNDM-5 gene. In this study, the blaNDM-5-harboring plasmid was first detected from E. coli ST156 in the feces of a layer hen in China, which may increase our understanding of the transmission of blaNDM-5.
IncX3 plasmids are narrow host range plasmids of Enterobacteriaceae and are believed to have a low prevalence . Since the first discovery of blaNDM-5 in China, this gene has been identified in a variety of Enterobacteriaceae [21, 31], with IncX3 being the primary type of Inc. to harbor blaNDM-5 . From our results, the IncX3 plasmids harboring blaNDM-5 were highly similar to each other in different countries and host sources, suggesting its ability to be an efficient vehicle for blaNDM-5 dissemination among humans, animals, food and the environment, potentially indicating its role in the rapid spread of blaNDM-5-harboring isolates [21, 28]. The BRIG analysis results showed that blaNDM-5-harboring IncX3 plasmids have a conserved backbone of ~ 46 kb, indicating that these plasmids had a common ancestor, and the conjugation/type IV secretion components in the backbone may be a factor promoting its transmission.
The blaNDM-5-harboring plasmids were initially detected from isolates from human [4, 27, 28]. However, they have also been detected in food, the environment and livestock and poultry sources in recent years. For example, plasmid pNDM5_025943 (unpublished) was detected in sewage, and plasmid pVH1 was detected from a cucumber . Carbapenem resistance is well known to be a universal phenomenon because of its frequent usage in clinics. Thus, it is interesting that the blaNDM-5-harboring plasmid has an increasing host range, which reflects the development of serious carbapenem resistance. In particular, the blaNDM-5 gene has been detected from livestock animals in recent years, such as swine  and dairy cows . In this study, the complete sequence of a blaNDM-5-harboring plasmid isolated from layer hen feces was first published, which is important evidence of blaNDM-5 transmission in poultry in China.
Materials and methods
Sample collection and antimicrobial susceptibility testing
The E. coli isolate recovered was designated ECCRA-119 and showed resistance to meropenem (8 mg/L) and imipenem (4 mg/L). This isolate was selected for AST using the broth dilution method with the Biofosun® Gram-negative panel (Fosun Diagnostics, Shanghai, China). The criteria from the Clinical and Laboratory Standards Institute (CLSI) were used to interpret the results, and the US National Antimicrobial Resistance Monitoring System (NARMS) and European Committee on Antimicrobial Susceptibility Testing (EUCAST) protocol was used when CLSI standards were not appropriate. The panel of antimicrobial compounds tested included ampicillin, amoxicillin/clavulanic acid, tetracycline, doxycycline, gentamicin, spectinomycin, sulfisoxazole, trimethoprim/sulfamethoxazole, ceftiofur, ceftazidime, florfenicol, enrofloxacin, ofloxacin, colistin, meropenem, and mequindox from among livestock antibiotics, and ampicillin, ampicillin/sulbactam, tetracycline, chloramphenicol, trimethoprim/sulfamethoxazole, cefazolin, cefotaxime, ceftazidime, cefoxitin, gentamicin, imipenem, nalidixic acid, azithromycin, sulfisoxazole, ciprofloxacin, amoxicillin/clavulanate, Cefotaxime/clavulanate, ceftazidime/clavulanate, colistin, polymyxin B, minocycline, amikacin, aztreonam, cefepime, meropenem, levofloxacin, doxycycline, kanamycin, streptomycin, and gemifloxacin from among medical antibiotics.
Whole genome sequencing, assembly and annotation
After genomic DNA extraction and quality checks, a 20-kb fragment library was constructed for the sample when the concentration and purity met the sequencing requirements. Whole-genome sequencing was performed using a PacBio RS II instrument . The assembly of the reads was performed following the Hierarchical Genome Assembly Process (HGAP) workflow . In this process, the Celera Assembler, following the OLC algorithm, was used to assemble the sequences , and Quiver was used to optimize the assembly results . The gene prediction and annotation of the genomes was performed using the NCBI Prokaryotic Genome Annotation Pipeline . The complete genome of the ECCRA-119 isolate was deposited in GenBank under the accession numbers CP029242 (chromosome), CP029243 (plasmid pTB201), CP029244 (plasmid pTB202) and CP029245 (plasmid pTB203).
CRISPRfinder (https://crispr.i2bc.paris-saclay.fr/Server/) was used to search for CRISPR loci in the genome of the ECCRA-119 isolate . MLST 2.0 (https://cge.cbs.dtu.dk/services/MLST/) was used to determine the ST . The plasmid replicon types were identified using PlasmidFinder-1.3 (https://cge.cbs.dtu.dk/services/PlasmidFinder/) . Acquired antimicrobial resistance genes were predicted using ResFinder (https://cge.cbs.dtu.dk/services/ResFinder/) . VirulenceFinder (https://cge.cbs.dtu.dk/services/VirulenceFinder/) was used to identify the virulence factors , and oriTfinder (http://184.108.40.206/oriTfinder/oriTfinder.html) was used to identify the origin of transfer in the genome . The genome was investigated for metal resistance genes using the Antibacterial Biocide and Metal Resistance Genes Database (BacMet) (http://bacmet.biomedicine.gu.se/) . Easyfig  and BIRG  were used in the comparative analysis of the plasmids. Phylogenetic analysis of genome and plasmids was performed by KSNP based on the maximum-likelihood method . Clustal X was used to perform the alignment analysis of blaNDM based on nucleotide sequences . The phylogenetic tree was generated used in MEGA X  and iTOL .
Plasmid conjugation experiments were performed on the ECCRA-119 isolate as described previously by Lin et al. [23, 61]. A rifamycin-resistant E. coli EC600 strain was used as the recipient in the plasmid conjugation assay to test the transferability of the carbapenem resistance gene and other resistance genes harbored by the ECCRA-119 isolate. Briefly, transconjugants were selected on LB agar plates (Landbridge., Beijing, China) supplemented with rifamycin (400 mg/L) (Sangon Biotech., Shanghai, China) and meropenem (4 mg/L) (J&K Chemical Ltd., Shanghai, China). The transfer frequencies were calculated by dividing the number of colony-forming units (CFUs) of transconjugants by the number of CFUs of the recipients. Genome DNA was extracted from the E. coli transconjugant using a bacterial DNA extract kit (Generay, Shanghai, China). The blaNDM-5 primers (F: 5′-GTCTGGCAGCACACTTCCTA-3′; R: 5′- TAGTGCTCAGTGTCGGCATC-3′) were used to confirm that the transconjugant harbored the plasmid.
S1-PFGE was performed according to a standard protocol using the contour-clamped homogeneous electric field (CHEF) technique with 0.5 × TBE buffer . Salmonella enterica serotype Braenderup H9812 was used as a size marker . The gels were run at 6 V/cm and 14 °C with an angle of 120°, and initial and final pulses were set at 2.16 and 63.8 s, respectively. The running time was 16 h using the CHEF apparatus (CHEF MAPPER XA; Bio-Rad, USA).
In this study, we reported the isolation and characterization of a carbapenem-resistant E. coli strain ST156 harboring the blaNDM-5 gene from a layer hen farm in Zhejiang province, China. Three plasmids in ECCRA-119 were identified based on whole genome sequencing and S1-PFGE. Twenty-two acquired resistance genes were identified, and this finding is consistent with the MDR phenotype of strain ECCRA-119. In particular, the blaNDM-5 gene has a high risk of spreading widely due to the potential transfer ability of the IncX3 plasmid pTB203 in this strain. The results of our study may reflect the level of antimicrobial resistance in poultry breeding in Zhejiang province and increase our knowledge of the presence and transmission of the blaNDM-5 gene.
Antimicrobial susceptibility testing
Antibacterial Biocide and Metal Resistance Genes Database
clamped homogeneous electric field
Clinical and Laboratory Standards Institute
European Committee on Antimicrobial Susceptibility Testing
Hierarchical Genome Assembly Process
Minimum inhibitory concentration
Multilocus sequence typing
The US National Antimicrobial Resistance Monitoring System
New Delhi metallo-β-lactamase
- oriT :
Origin of transfer region
Single nucleotide polymorphism
The type IV coupling protein
Bacterial type IV secretion system
Spellberg B, Guidos R, Gilbert D, Bradley J, Boucher HW, Scheld WM, Bartlett JG, Edwards J Jr. The epidemic of antibiotic-resistant infections: a call to action for the medical community from the Infectious Diseases Society of America. Clin Infect Dis. 2008;46(2):155–64. https://doi.org/10.1086/524891.
Yong D, Toleman MA, Giske CG, Cho HS, Sundman K, Lee K, Walsh TR. Characterization of a new metallo-beta-lactamase gene, bla NDM-1, and a novel erythromycin esterase gene carried on a unique genetic structure in Klebsiella pneumoniae sequence type 14 from India. Antimicrob Agents Chemother. 2009;53(12):5046–54. https://doi.org/10.1128/AAC.00774-09.
Zafer MM, Amin M, El Mahallawy H, Ashour MS, Al Agamy M. First report of NDM-1-producing Pseudomonas aeruginosa in Egypt. Int J Infect Dis. 2014;29:80–1. https://doi.org/10.1016/j.ijid.2014.07.008.
Hornsey M, Phee L, Wareham DW. A novel variant, NDM-5, of the New Delhi metallo-β-lactamase in a multidrug-resistant Escherichia coli ST648 isolate recovered from a patient in the United Kingdom. Antimicrob Agents Chemother. 2011;55(12):5952–4. https://doi.org/10.1128/AAC.05108-11.
Rahman M, Shukla SK, Prasad KN, Ovejero CM, Pati BK, Tripathi A, Singh A, Srivastava AK, Gonzalez-Zorn B. Prevalence and molecular characterisation of New Delhi metallo-β-lactamases NDM-1, NDM-5, NDM-6 and NDM-7 in multidrug-resistant Enterobacteriaceae from India. Int J Antimicrob Agents. 2014;44(1):30–7. https://doi.org/10.1016/j.ijantimicag.2014.03.003.
Sassi A, Loucif L, Gupta SK, Dekhil M, Chettibi H, Rolain JM. NDM-5 carbapenemase-encoding gene in multidrug-resistant clinical isolates of Escherichia coli from Algeria. Antimicrob Agents Chemother. 2014;58(9):5606–8. https://doi.org/10.1128/AAC.02818-13.
Nakano R, Nakano A, Hikosaka K, Kawakami S, Matsunaga N, Asahara M, Ishigaki S, Furukawa T, Suzuki M, Shibayama K, Ono Y. First report of metallo-β-lactamase NDM-5-producing Escherichia coli in Japan. Antimicrob Agents Chemother. 2014;58(12):7611–2. https://doi.org/10.1128/AAC.04265-14.
Cho SY, Huh HJ, Baek JY, Chung NY, Ryu JG, Ki CS, Chung DR, Lee NY, Song JH. Klebsiella pneumoniae co-producing NDM-5 and OXA-181 carbapenemases, South Korea. Emerg Infect Dis. 2015;21(6):1088–9. https://doi.org/10.3201/eid2106.150048.
Wailan AM, Paterson DL, Caffery M, Sowden D, Sidjabat HE. Draft genome sequence of NDM-5-producing Escherichia coli sequence type 648 and genetic context of bla NDM-5 in Australia. Genome Announc. 2015;3(2):e00194–15. https://doi.org/10.1128/genomeA.00194-15.
Zhang LP, Xue WC, Meng DY. First report of New Delhi metallo-β-lactamase 5 (NDM-5)-producing Escherichia coli from blood cultures of three leukemia patients. Int J Infect Dis. 2016;42:45–6. https://doi.org/10.1016/j.ijid.2015.10.006.
Bathoorn E, Rossen JW, Lokate M, Friedrich AW, Hammerum AM. Isolation of an NDM-5-producing ST16 Klebsiella pneumoniae from a Dutch patient without travel history abroad, August 2015. Euro Surveill. 2015;20(41):30040. https://doi.org/10.2807/1560-7917.ES.2015.20.41.30040.
Bitar I, Piazza A, Gaiarsa S, Villa L, Pedroni P, Oliva E, Nucleo E, Pagani L, Carattoli A, Migliavacca R. ST405 NDM-5 producing Escherichia coli in northern Italy: the first two clinical cases. Clin Microbiol Infect. 2017;23(7):489–90. https://doi.org/10.1016/j.cmi.2017.01.020.
Rojas LJ, Hujer AM, Rudin SD, Wright MS, Domitrovic TN, Marshall SH, Hujer KM, Richter SS, Cober E, Perez F, Adams MD, van Duin D, Bonomo RA. NDM-5 and OXA-181 beta-lactamases, a significant threat continues to spread in the Americas. Antimicro Agents Chemother. 2017;61(7):e00454–17. https://doi.org/10.1128/AAC.00454-17.
Pérez-Moreno MO, Ortega A, Pérez-Vázquez M, Centelles-Serrano MJ, Bautista V, Escrig-Monfort C, Oteo J. Simultaneous colonisation by ST340 Klebsiella pneumoniae producing NDM-5 and ST399 Escherichia coli producing NDM-7. Int J Antimicrob Agents. 2016;48(4):464–6. https://doi.org/10.1016/j.ijantimicag.2016.07.003.
Gamal D, Fernández-Martínez M, El-Defrawy I, Ocampo-Sosa AA, Martínez-Martínez L. First identification of NDM-5 associated with OXA-181 in Escherichia coli from Egypt. Emerg Microbes Infect. 2016;5:e30. https://doi.org/10.1038/emi.2016.24.
Almakki A, Maure A, Pantel A, Romano-Bertrand S, Masnou A, Marchandin H, Jumas-Bilak E, Licznar-Fajardo P. NDM-5-producing Escherichia coli in an urban river in Montpellier, France. Int J Antimicrob Agents. 2017;50(1):123–4. https://doi.org/10.1016/j.ijantimicag.2017.04.003.
Howard JC, Creighton J, Heffernan H, Werno A. Evidence of transmission of an NDM-5-producing Klebsiella pneumoniae in a healthcare facility in New Zealand. J Antimicrob Chemother. 2017;72(3):949–51. https://doi.org/10.1093/jac/dkw498.
Huang Y, Yu X, Xie M, Wang X, Liao K, Xue W, Chan EW, Zhang R, Chen S. Widespread dissemination of carbapenem-resistant Escherichia coli sequence type 167 strains harboring bla NDM-5 in clinical settings in China. Antimicrob Agents Chemother. 2016;60(7):4364–8. https://doi.org/10.1128/AAC.00859-16.
Li X, Jiang Y, Wu K, Zhou Y, Liu R, Cao Y, Wu A, Qiu Y. Whole-genome sequencing identification of a multidrug-resistant Salmonella enterica serovar typhimurium strain carrying bla NDM-5 from Guangdong, China. Infect Genet Evol. 2017;55:195–8. https://doi.org/10.1016/j.meegid.2017.09.005.
Zhu YQ, Zhao JY, Xu C, Zhao H, Jia N, Li YN. Identification of an NDM-5-producing Escherichia coli sequence type 167 in a neonatal patient in China. Sci Rep. 2016;6:29934. https://doi.org/10.1038/srep29934.
Zhang F, Xie L, Wang X, Han L, Guo X, Ni Y, Qu H, Sun J. Further spread of bla NDM-5 in Enterobacteriaceae via IncX3 plasmids in Shanghai, China. Front Microbiol. 2016;7:424. https://doi.org/10.3389/fmicb.2016.00424.
Ho PL, Wang Y, Liu MC, Lai EL, Law PY, Cao H, Chow KH. IncX3 epidemic plasmid carrying bla NDM-5 in Escherichia coli from swine in multiple geographic areas in China. Antimicrob Agents Chemother. 2018;62(3):e02295–17. https://doi.org/10.1128/AAC.02295-17.
Kong LH, Lei CW, Ma SZ, Jiang W, Liu BH, Wang YX, Guan R, Men S, Yuan QW, Cheng GY, Zhou WC, Wang HN. Various sequence types of Escherichia coli isolates coharboring bla NDM-5 and mcr-1 genes from a commercial swine farm in China. Antimicrob Agents Chemother. 2017;61(3):e02167–16. https://doi.org/10.1128/AAC.02167-16.
He T, Wang Y, Sun L, Pang M, Zhang L, Wang R. Occurrence and characterization of bla NDM-5-positive Klebsiella pneumoniae isolates from dairy cows in Jiangsu, China. J Antimicrob Chemother. 2017;72(1):90–4. https://doi.org/10.1093/jac/dkw357.
Liu BT, Zhang XY, Wan SW, Hao JJ, Jiang RD, Song FJ. Characteristics of carbapenem-resistant Enterobacteriaceae in ready-to-eat vegetables in China. Front Microbiol. 2018;9:1147. https://doi.org/10.3389/fmicb.2018.01147.
Smillie C, Garcillán-Barcia MP, Francia MV, Rocha EP, de la Cruz F. Mobility of plasmids. Microbiol Mol Biol Rev. 2010;74(3):434–52. https://doi.org/10.1128/MMBR.00020-10.
Krishnaraju M, Kamatchi C, Jha AK, Devasena N, Vennila R, Sumathi G, Vaidyanathan R. Complete sequencing of an IncX3 plasmid carrying bla NDM-5 allele reveals an early stage in the dissemination of the bla NDM gene. Indian J Med Microbiol. 2015;33(1):30–8. https://doi.org/10.4103/0255-0857.148373.
Li X, Fu Y, Shen M, Huang D, Du X, Hu Q, Zhou Y, Wang D, Yu Y. Dissemination of bla NDM-5 gene via an IncX3-type plasmid among non-clonal Escherichia coli in China. Antimicrob Resist Infect Control. 2018;7:59. https://doi.org/10.1186/s13756-018-0349-6.
Xie M, Li R, Liu Z, Chan EWC, Chen S. Recombination of plasmids in a carbapenem-resistant NDM-5-producing clinical Escherichia coli isolate. J Antimicrob Chemother. 2018;73(5):1230–4. https://doi.org/10.1093/jac/dkx540.
Paskova V, Medvecky M, Skalova A, Chudejova K, Bitar I, Jakubu V, Bergerova T, Zemlickova H, Papagiannitsis CC, Hrabak J. Characterization of NDM-encoding plasmids from Enterobacteriaceae recovered from Czech hospitals. Front Microbiol. 2018;9:1549. https://doi.org/10.3389/fmicb.2018.01549.
Chen D, Gong L, Walsh TR, Lan R, Wang T, Zhang J, Mai W, Ni N, Lu J, Xu J, Li J. Infection by and dissemination of NDM-5-producing Escherichia coli in China. J Antimicrob Chemother. 2016;71(2):563–5. https://doi.org/10.1093/jac/dkv352.
Fu L, Wang S, Zhang Z, Yan X, Yang X, Zhang L, Li Y, Wang G, Zhao K, Zhou Y. Co-carrying of KPC-2, NDM-5, CTX-M-3 and CTX-M-65 in three plasmids with serotype O89: H10 Escherichia coli strain belonging to the ST2 clone in China. Microb Pathog. 2019;128:1–6. https://doi.org/10.1016/j.micpath.2018.12.033.
Liu BT, Liao XP, Yue L, Chen XY, Li L, Yang SS, Sun J, Zhang S, Liao SD, Liu YH. Prevalence of β-lactamase and 16S rRNA methylase genes among clinical Escherichia coli isolates carrying plasmid-mediated quinolone resistance genes from animals. Microb Drug Resist. 2013;19(3):237–45. https://doi.org/10.1089/mdr.2012.0179.
Li G, Zhang Y, Bi D, Shen P, Ai F, Liu H, Tian Y, Ma Y, Wang B, Rajakumar K, Ou HY, Jiang X. First report of a clinical, multidrug-resistant Enterobacteriaceae isolate coharboring fosfomycin resistance gene fosA3 and carbapenemase gene bla KPC-2 on the same transposon, Tn1721. Antimicrob Agents Chemother. 2015;59(1):338–43. https://doi.org/10.1128/AAC.03061-14.
Johnson TJ, Wannemuehler YM, Nolan LK. Evolution of the iss gene in Escherichia coli. Appl Environ Microbiol. 2008;74(8):2360–9. https://doi.org/10.1128/AEM.02634-07.
Mcintosh D, Cunningham M, Ji B, Fekete FA, Parry EM, Clark SE, Zalinger ZB, Gilg IC, Danner GR, Johnson KA, Beattie M, Ritchie R. Transferable, multiple antibiotic and mercury resistance in Atlantic Canadian isolates of Aeromonas almonicida subsp. salmonicida is associated with carriage of an IncA/C plasmid similar to the Salmonella enterica plasmid pSN254. J Antimicrob Chemother. 2008;61(6):1221–8. https://doi.org/10.1093/jac/dkn123.
Mushtaq S, Irfan S, Sarma JB, Doumith M, Pike R, Pitout J, Livermore DM, Woodford N. Phylogenetic diversity of Escherichia coli strains producing NDM-type carbapenemases. J Antimicrob Chemother. 2011;66(9):2002–5. https://doi.org/10.1093/jac/dkr226.
Giufre M, Graziani C, Accogli M, Luzzi I, Busani L, Cerquetti M. Escherichia coli of human and avian origin: detection of clonal groups associated with fluoroquinolone and multidrug resistance in Italy. J Antimicrob Chemother. 2012;67(4):860–7. https://doi.org/10.1093/jac/dkr565.
Yang RS, Feng Y, Lv XY, Duan JH, Chen J, Fang LX, Xia J, Liao XP, Sun J, Liu YH. Emergence of NDM-5- and MCR-1-producing Escherichia coli clones ST648 and ST156 from a single muscovy duck (Cairina moschata). Antimicrob Agents Chemother. 2016;60(11):6899–902. https://doi.org/10.1128/AAC.01365-16.
Yang P, Xie Y, Feng P, Zong Z. bla NDM-5 carried by an IncX3 plasmid in Escherichia coli sequence type 167. Antimicrob Agents Chemother. 2014;58(12):7548–52. https://doi.org/10.1128/AAC.03911-14.
Shen ZQ, Hu YY, Sun QL, Hu FP, Zhou HW, Shu LB, Ma TF, Shen YB, Wang Y, Li J, Timothy RW, Zhang R, Wang SL. Emerging carriage of NDM-5 and MCR-1 in Escherichia coli from healthy people in multiple regions in China: a cross sectional observational study. EClinicalMedicine. 2018;6:11–20. https://doi.org/10.1016/j.eclinm.2018.11.003.
Leon-Velarde CG, Cai HY, Larkin C, Bell-Rogers P, Stevens RW, Odumeru JA. Evaluation of methods for the identification of Salmonella enterica serotype typhimurium DT104 from poultry environmental samples. J Microbiol Methods. 2004;58(1):79–86. https://doi.org/10.1016/j.mimet.2004.03.005.
Hassan ARHA, Salam HSH, Abdel-Latef GK. Serological identification and antimicrobial resistance of Salmonella isolates from broilers carcasses and human stools in Beni-Suef, Egypt. Beni-Suef University Journal of Basic and Applied Sciences. 2016;5(2):202–7. https://doi.org/10.1016/j.bjbas.2016.04.002.
Roberts RJ, Carneiro MO, Schatz MC. The advantages of SMRT sequencing. Genome Biol. 2013;14(7):405. https://doi.org/10.1186/gb-2013-14-6-405.
Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, Clum A, Copeland A, Huddleston J, Eichler EE, Turner SW, Korlach J. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10(6):563–9. https://doi.org/10.1038/nmeth.2474.
Miller JR, Koren S, Sutton G. Assembly algorithms for next-generation sequencing data. Genomics. 2010;95(6):315–27. https://doi.org/10.1016/j.ygeno.2010.03.001.
Tatusova T, Dicuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, Ostell J. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016;44(14):6614–24. https://doi.org/10.1093/nar/gkw569.
Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:52–7. https://doi.org/10.1093/nar/gkm360.
Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, Jelsbak L, Sicheritz-Ponten T, Ussery DW, Aarestrup FM, Lund O. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012;50(4):1355–61. https://doi.org/10.1128/JCM.06094-11.
Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O, Villa L, Moller Aarestrup F, Hasman H. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58(7):3895–903. https://doi.org/10.1128/AAC.02412-14.
Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, Aarestrup FM, Larsen MV. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67(11):2640–4. https://doi.org/10.1093/jac/dks261.
Joensen KG, Scheutz F, Lund O, Hasman H, Kaas RS, Nielsen EM, Aarestrup FM. Real-time whole-genome sequencing for routine typing, surveillance, and outbreak detection of verotoxigenic Escherichia coli. J Clin Microbiol. 2014;52(5):1501–10. https://doi.org/10.1128/JCM.03617-13.
Li X, Xie Y, Liu M, Tai C, Sun J, Deng Z, Ou HY. oriTfinder: a web-based tool for the identification of origin of transfers in DNA sequences of bacterial mobile genetic elements. Nucleic Acids Res. 2018;46(W1):W229–34. https://doi.org/10.1093/nar/gky352.
Pal C, Bengtsson-Palme J, Rensing C, Kristiansson E, Larsson DG. BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 2014;42:737–43. https://doi.org/10.1093/nar/gkt1252.
Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10. https://doi.org/10.1093/bioinformatics/btr039.
Alikhan NF, Petty NK, Ben Zakour NL, Beatson SA. BLAST ring image generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011;12:402. https://doi.org/10.1186/1471-2164-12-402.
Gardner SN, Hall BG. When whole-genome alignments just won't work: kSNP v2 software for alignment-free SNP discovery and phylogenetics of hundreds of microbial genomes. PLoS One. 2013;8(12):e81760. https://doi.org/10.1371/journal.pone.0081760.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Hiqqins DG. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8. https://doi.org/10.1093/bioinformatics/btm404.
Kumar S, Tamura K, Nei M. MEGA: molecular evolutionary genetics analysis software for microcomputers. Comput Appl Biosci. 1994;10(2):189–91. https://doi.org/10.1093/bioinformatics/10.2.189.
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:242–5. https://doi.org/10.1093/nar/gkw290.
Chen L, Chen ZL, Liu JH, Zeng ZL, Ma JY, Jiang HX. Emergence of RmtB methylase-producing Escherichia coli and Enterobacter cloacae isolates from pigs in China. J Antimicrob Chemother. 2007;59(5):880–5. https://doi.org/10.1093/jac/dkm065.
Gouby A, Neuwirth C, Bourg G, Bouziges N, Carles-Nurit MJ, Despaux E, Ramuz M. Epidemiological study by pulsed-field gel electrophoresis of an outbreak of extended-spectrum beta-lactamase-producing Klebsiella pneumoniae in a geriatric hospital. J Clin Microbiol. 1994;32(2):301–5.
Hunter SB, Vauterin P, Lambert-Fair MA, Van Duyne MS, Kubota K, Graves L, Wrigley D, Barrett T, Ribot E. Establishment of a universal size standard strain for use with the PulseNet standardized pulsed-field gel electrophoresis protocols: converting the national databases to the new size standard. J Clin Microbiol. 2005;43(3):1045–50. https://doi.org/10.1128/JCM.43.3.1045-1050.2005.
We would like to thank Shanghai Biofosun Biotech Co. Ltd. for providing antibiotic susceptibility test assistance.
This work was financially supported by the Public Welfare Technology Applied Research Projects in Zhejiang Province Project (2017C32021), the Natural Science Foundation for the Youth (31700007), the Youth Talent Training and International Cooperation program of Zhejiang Academy of Agricultural Sciences (2017R19R08E01), the Key research and development plan of Zhejiang Province (2015C02041) and the State Key Laboratory Breeding Base for Zhejiang Sustainable Pest and Disease Control (2011DS700124-ZZ1703). The funding bodies had no role in the design of the study, in the collection, analysis, interpretation of data, or in the writing of the manuscript.
Availability of data and materials
All of the data generated or analyzed during this study are included in this article and its supplementary information files. The genome sequence of Escherichia coli ECCRA-119 has been deposited in the GenBank database under accession numbers CP029242, CP029243, CP029244 and CP029245.
Ethics approval and consent to participate
For convenient sampling, fecal samples were collected for bacterial isolation after the farmer’s verbal consent. All of the activities in our study were approved by institutional review board of Zhejiang Academy of Agricultural Sciences.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The gene mutation on the chromosome of strain ECCRA-119. Table S2. Virulence factors of strain ECCRA-119. Table S3. Metal resistance genes of strain ECCRA-119. Figure S1. SNP tree of E. coli ST156 strains. Figure S2. The comparison analysis of 52 blaNDM-5-harboring IncX3 plasmid sequences. Figure S3. Comparison analysis of 24 blaNDM variant sequences based on nucleotide sequences. Figure S4. Phylogenetic relationships between blaNDM variants based on nucleotide sequences. (DOCX 2669 kb)
Table S4. The information for the blaNDM-harboring plasmids with complete sequences in China. (XLSX 20 kb)