Clinical data collection
Of 654 unique individuals with CNS-GNB samples, 549 (83.9%, 549/654) were from Southern China, mainly Guangdong province (97.3%, 534/549). The highest proportion was from Guangzhou in Guangdong (Supplementary Figure 2); 71 (10.9%, 71/654) were from Central China, 16 (2.4%, 16/654) from Eastern China; and 18 (2.8%, 18/654) from other regions in China. A total of 415 (63.5%, 415/654) were men and 239 (36.5%, 239/654) were women and the median age was 55 years (range, < 1–91 years; Table 3). The samples were collected from almost all the wards across the hospitals, including ICU (34.3%, 224/654) and nearly 40 wards (65.7%, 430/654) of non-ICU departments. A total of 653 (99.8%, 653/654) had at least one underlying comorbid condition with a median Charlson comorbidity index of 2 (range, 0–4) and 1 (0.2%, 1/654) had no documented underlying condition. The most commonly reported conditions included diabetes (15.1%, 99/654) and neurological disorders (28.7%, 188/654); the ratio of malignancy (15.9%, 104/654) and chronic pulmonary disease (13.5%, 88/654) also surpassed other general underlying diseases. Bacterial infections were prominent, and most of the cases had a history of antibiotic use.
Identification of isolates
A total of 745 samples harbored CNS-GNB in 10,000 fecal samples, and 680 CNS-GNB were identified in 654 samples excluding 91 samples with Stenotrophomonas maltophilia only (Supplementary Figure 1). A total of 22 genera were identified, including nine types of 504 (74.1%, 504/680) carbapenem non-susceptible non-fermenting bacillus (CNS-NFB), eight types of 156 (22.9%, 156/680) carbapenem non-susceptible Enterobacteria (CNSE), and four types of 20 (2.9%, 20/680) other CNS-GNB. The genera isolated in the highest proportion were Acinetobacter spp. (33.8%, 230/680), Pseudomonas spp. (25.1%, 171/680), and Enterobacter spp. (9.0%, 61/680) (Table 1).
Antimicrobial susceptibility testing
Most of the isolates were resistant to cephalosporin, but they were sensitive to colistin. Amikacin resistance was identified in 147 isolates (21.6%, 147/680). Most of the CNS-GNB had a sensitivity rate of more than 80%, and sensitivity to amikacin of Acinetobacter baumannii was 54.8% (23/42). Most of CNSE strains had a resistance rate of less than 10% to amikacin; E.coli had a resistance rate of 25.9% (7/27). A total of 490 (72.1%, 490/680) isolates were ciprofloxacin resistant, 391 (57.5%, 391/680) isolates were levofloxacin resistant, and the susceptibility of different CNSE species to ciprofloxacin and levofloxacin varied. CNSE species were more sensitive to trimethoprim/sulfamethoxazole than non-fermenting bacillus, and 527 (77.5%, 527/680) isolates showed resistance to trimethoprim/sulfamethoxazole. A total of 354 (52.1%, 354/680) isolates were resistant to tetracycline, and non-fermenting bacillus species were more sensitive to tetracycline than CNSE (Table 2).
Investigation of carbapenemase and colistin resistance genes
A total of 307 (45.1%, 307/680) isolates harbored carbapenemase genes. Among them, blaNDM was detected in 206 (67.1%, 206/307) isolates, including 188 (61.2%, 188/307) blaNDM-1, 2 (0.7%, 2/307) blaNDM-4, 15 (4.9%, 15/307) blaNDM-5, and one (0.3%, 1/307) blaNDM-7. In addition, 51 (16.6%, 51/307) isolates harbored the blaVIM-2 gene, and blaIMP and blaKPC-2 were observed in 48 (15.6%, 48/307) and seven (2.3%, 7/307) isolates, respectively. Four isolates were found to co-exist with the blaIMP and blaKPC genes. Seven isolates were found to harbor the blaAFM-1 gene, a novel metallo-β-lactamase gene; these isolates include Alcaligenes faecalis, Bordetella trematum, Comamonas testosterone, and Comamonas aquatica. Five strains of Acinetobacter and one strain of K. pneumoniae were positive for the colistin-resistant gene mcr-1 and co-existed with the blaNDM gene. No other carbapenemase-related and colistin resistance genes were identified (Table 1). Acinetobacter was the most abundant carbapenemase-producing species (87.8%, 115/131), followed by Escherichia coli (70.4%, 19/27), Providencia rettgeri (64.7%, 11/17), Citrobacter freundii (64.3%, 9/14), and Klebsiella pneumoniae (58.5%, 24/41). Pseudomonas aeruginosa and Klebsiella pneumoniae had the most carbapenemase gene types; blaNDM-1 was the most common carbapenemase gene in Acinetobacter spp. (100.0%, 115/115), blaNDM-5 was the most common type in E. coli (68.4%, 13/19), and blaKPC-2 was the most common in K. pneumoniae.
Conjugation
The transconjugants were selected on screening plates. PCR and sequencing by Sangon Biotech Co. Ltd. (Shanghai, China) revealed that the transconjugant NFYY0065-N contained the blaNDM-5 gene, and the transconjugant NFYY0065-M contained the mcr-1 gene. Antimicrobial susceptibility testing showed that NFYY0065-N was resistant to imipenem, meropenem, ceftazidime, and aztreonam, whereas NFYY0065-M was resistant to colistin (Table 4).
Genomic analysis of K. pneumoniae co-producing blaNDM-5 and mcr-1
The chromosome was 5,235,159 bp with 57.55% GC content; it harbored 5186 genes, including five drug resistance genes encoding a class A broad-spectrum beta-lactamase, multidrug efflux transporter transcriptional repressor, tetracycline repressor protein, multidrug efflux transporter, and fosfomycin resistance protein (Table 5).
pAN65-MCR was 276,370 bp, possessed a typical IncHI2-type backbone, and contained 319 predicted open reading frames. It harbored the mcr-1 gene, which was located between the IS30-like element ISApl1 and a PAP2-like encoding gene (Figure 1). The plasmid co-carried the resistance gene blaCTX-M-14. BLAST analysis of the complete nucleotide sequence against the NCBI database indicated the highest similarity (95% query cover and 99% nucleotide identity) to plasmid pMCR_WCHEC050613 (CP019214.2), which was from E. coli isolated from sewage in Chengdou. The genetic structure of mcr-1 in plasmid pMCR_WCHEC050613 was also the same as pAN65-MCR. BLAST homology analyses revealed that pAN65-MCR has 93% identity and 99% query coverage with pHNSHP45 [6], which was isolated from pigs. Comparative analysis of the genetic structures of mcr-1 in the reference plasmids pHN6DS2 and pMCR1-PA revealed there were identical to pAN65-1, and pHNSHP45 had lost the PAP2 gene downstream of mcr (Figure 2).
The blaNDM-5-bearing plasmid pAN65-NDM was 46,161 bp with a GC content of 46.65% and belonged to the IncX3 group (Figure 3). BLAST analysis revealed high similarity (100% query cover and 99% nucleotide identity) to previously reported blaNDM-5-positive plasmids, including the blaNDM-5-positive plasmids of various sources, such as MH234502 (human, China) [8], CP027204 (sewage, China) [8], CP029245 (laying hen, China) [8], CP028705 (vegetable, China) [8], KX507346 (swine, China) [8], and CP034744 (goose, China) [8]. Comparative analysis revealed that these plasmids are highly similar, possessing the same backbone, which includes IncX3 replication. The genetic structure of the blaNDM-5 gene was identical to the plasmids pNDM5_020038, pK516_NDM5, and tig00001251; the mobile genetic elements IS30 and IS5 were upstream of blaNDM-5 and the genes bleMBL and trpF were downstream (Figure 4).