Isolates collections
Firstly, eleven CREC isolates were collected from 6 hospitals during 2015 June-2015 December in Nanjing city. These strains were isolated from urine (n=4), sputum (n=3), blood (n=1), bile (n=1), secretion of uterus neck (n=1), and the source of one strain was unknown.
A total of 43 consecutive non-duplicate isolates collected in Nanjing Drum Tower hospital during 2013-2017 were further analyzed for isolates co-producing KPC-2 and NDM. Among them, 4 strains were isolated in 2013, 10 in 2014, 2 in 2015, and 11 from 2016 and 16 in 2017. The source of the samples was as follows: urine (n=18), blood (n=9), sputum (n=6), secretion (n=3), bile (n=3), abdominal dropsy (n=2), and pus (n=1).
All the CREC strains were identified by Vitek 2.0 (BioMérieux. Marcy l'Etoile, France) or ATB 32E Semi-auto identification machine (Bio-
Mérieux, France). Isolates resistant to at least one carbapenem (imipenem, meropenem, ertapenem) were included in the study.
DNA Extraction
The Ultraclean Microbial DNA Isolation Kit (MOBIO Laboratories, Carlsbad, CA, US) was used to extract genomic DNA. The NanoDrop 2000c spectrophotometer (Thermo Scientific, Waltham, MA, USA) was used for measuring the DNA concentration and purity for whole genome sequencing.
Whole Genome Sequencing, denovo Assembly, Scaffolding, and Annotation
The prepared pair-end DNA library was sequenced on the MiSeq (Illumina, SanDiego, CA, USA). Denovo assembly of the paired-end reads was performed by CLC Genomics Workbenchv7.0.4 (QIAGEN, Hilden, Germany) after quality trimming (Qs ≥ 20). Scaffolding was finished using SSPACE standard version 3.0 and the gaps within scaffolds were further closed by GapFiller (13, 14). Then genomes were then submitted to NCBI for annnotation.
Analysis of genomic epidemiology
The antimicrobial resistance determinants and virulence factors were identified using Resfinder v2.1 (http://cge.cbs.dtu.dk/services/ResFinder-2.1/) and Virulence Finder 2.0 (https://cge.cbs.dtu.dk/services/VirulenceFinder/); the 11 CREC were typed by multi-locus sequence typing (MLST) 2.0 (https://cge.cbs.dtu.dk/services/MLST/), Plasmid Finder 2.1 (https://cge.cbs.dtu.dk/services/PlasmidFinder/), Serotype Finder 2.0 (https://cge.cbs.dtu.dk/services/SerotypeFinder/), and FimTyper 1.0 (https://cge.cbs.dtu.dk/services/FimTyper/).
Phylogenetic relationship of 11 CREC
The core-genome phylogeny of the 11 CREC was constructed by using single-nucleotide polymorphisms (SNP)-sites detected from 1010 core genes (identity > 95%; coverage = 100% ) derived from 502 ST11 strains (15). A maximum-likelihood tree was calculated using RAxML version 8.2.8, with general time-reversible model and 100 bootstraps (16). Interactive Tree Of Life (https://itol.embl.de) was used to produce the phylogenetic tree (17).
The screening of Escherichia coli co-producing NDM and KPC-2 carbapenemases
A total of 43 CREC strains collected from our hospital during 2013-2017 were further analyzed for the prevalence of E. coli co-producing NDM-5 and KPC-2. Presence of genes encoding carbapenemases (blaKPC, blaAIM, blaSIM, blaSPM, blaVIM, blaIMP, blaOXA, blaDIM, blaNDM, blaGIM, and blaGES) was analyzed by PCR and DNA sequencing (18).
Pulsed-field gel electrophoresis
Isolates selected from the 43 CREC were further analyzed for genetic relatedness by PFGE, which was performed according to the protocol as previously described (19). Briefly, fresh colonies were mixed with proteinase K (Merck Sharp & Dohme Ltd, Germany) into plugs. After the plugs were digested by restriction endonuclease XbaI (Fermentas, ABI, Germany), the resultant DNA fragments were separated in a PFGE CHEF-DR III system (Bio-Rad Laboratories, Hercules, CA) in 0.5×Tris-borate-EDTA buffer at 120 V for 19 h. The pulse times ranged from 2.2 s to 54.2 s. Finally, the BioNumerics software (Applied Math, Sint-Maten-Latem, Belgium) was used to analyze the banding patterns.
Conjugation assay
For the 6 isolates co-carrying blaKPC and blaNDM, broth mating was performed in order to analyze the transferability of these genes according to the protocol prescribed previously(20). Azide resistant E. coli J53 was used as the recipient. Briefly, fresh colonies were inoculated into 5 ml LB broth and incubated at 37℃, 200 rpm. After 5 hours, 500 μl recipient cells and 100 ul donor were suspended in 5 ml LB broth for overnight culture at 37℃, 200 rpm, then 100 μl were plated onto the LB plates containing 30 mg/L cefoxitin and 100 mg/L sodium azide for E. coli J53. PCR (amplification for blaNDM and blaKPC) and Eric-PCR were used to verify conjugants.
Nucleotide Sequence Genbank Accession Numbers
The Whole Genome Shotgun BioProject for these CREC isolates has been deposited at GenBank, the accession number of each strain were as follows: JABEWI000000000(NJSETYY11), JABEWJ000000000 (NJSETYY14), JABEWL000000000 (NJLSYY40), JABEWM000000000 (NJSDYYY42), JABEWN000000000 (NJJYYY51), RZMM00000000 (NJGLYY3940), RZMI00000000 (NJGLYY3610), JABEWH000000000 (NJMJYY4), JABEWK000000000 (NJSFYBJY17).