
==== Front
Pol J Microbiol
Pol J Microbiol
pjm
pjm
Polish Journal of Microbiology
1733-1331
2544-4646
Sciendo

39268958
pjm-2024-032
10.33073/pjm-2024-032
Original Paper
Molecular Epidemiology of mcr-1-Positive Polymyxin B-Resistant Escherichia coli Producing Extended-Spectrum β-Lactamase (ESBL) in a Tertiary Hospital in Shandong, China
Liu Yue 1#
Wang Qian 1#
Qi Ting 2
Zhang Meng 3
Chen Ran 1
Si Zaifeng 4
Li Jinmei 45
Jin Yan 1
Xu Qingbing 6
Li Ping 715168886362@163.com

Hao Yingying 1haoyingying@sdfmu.edu.cn

1 Department of Clinical Laboratory, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, China
2 Department of Clinical Laboratory, Jinan Gangcheng District People’s Hospital, Jinan, China
3 Department of Clinical Laboratory, Liaocheng Second People’s Hospital, Liaocheng, China
4 Department of Clinical Laboratory, Shandong Provincial Hospital, Cheeloo College of Medicine, Shandong University, Jinan, China
5 Department of Clinical Laboratory, Jinan Seventh People’s Hospital, Jinan, China
6 Central Laboratory of Liaocheng People’s Hospital, Liaocheng, China
7 Department of Joint Surgery, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, China
# Yue Liu and Qian Wang contributed equally to this work.

13 9 2024
9 2024
73 3 363375
15 5 2024
23 7 2024
© 2024 Yue Liu et al., published by Sciendo
2024
Yue Liu et al., published by Sciendo
https://creativecommons.org/licenses/by-nc-nd/4.0/ This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.
Abstract

Escherichia coli, a rod-shaped Gram-negative bacterium, is a significant causative agent of severe clinical bacterial infections. This study aimed to analyze the epidemiology of extended-spectrum β-lactamase (ESBL)-producing mcr-1 -positive E. coli in Shandong, China. We collected 668 non-duplicate ESBL-producing E. coli strains from clinical samples at Shandong Provincial Hospital between January and December 2018, and estimated their minimum inhibitory concentrations (MICs) using a VITEK® 2 compact system and broth microdilution. Next-generation sequencing and bioinformatic analyses identified the mcr-1 gene and other resistance genes in the polymyxin B-resistant strains. The conjugation experiment assessed the horizontal transfer capacity of the mcr-1 gene. Of the strains collected, 24 polymyxin B-resistant strains were isolated with a positivity rate of 3.59% and among the 668 strains, 19 clinical strains carried the mobile colistin resistance gene mcr-1, with a positivity rate of approximately 2.8%. All 19 clinical strains were resistant to ampicillin, cefazolin, ceftriaxone, ciprofloxacin, levofloxacin, and polymyxin B. Seventeen strains successfully transferred the mcr-1 gene into E. coli J53. All transconjugants were resistant to polymyxin B, and carried the drug resistance gene mcr-1. The 19 clinical strains had 14 sequence types (STs), with ST155 (n = 4) being the most common. The whole-genome sequencing results of pECO-POL-29_mcr1 revealed that no ISApl1 insertion sequences were found on either side of the mcr-1 gene. Our study uncovered the molecular epidemiology of mcr-1-carrying ESBL-producing E. coli in the region and suggested horizontal transmission mediated by plasmids as the main mode of mcr-1 transmission.

Keywords

Escherichia coli
extended-spectrum β-lactamase (ESBL)
mcr-1
resistance gene
insertion sequence ISAPl1
The study was supported by grants from National Natural Science Foundation of China (81902119) and Shandong Province Natural Science Foundation (ZR2021MH214).
==== Body
pmcIntroduction

Escherichia coli is a Gram-negative rod-shaped bacterium that is one of the most responsible for hospital-acquired infections, such as gastrointestinal infections, urinary tract infections, bloodstream infections, and neonatal meningitis (Riley 2020). The increasing prevalence of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Gram-negative bacteria is indeed a concerning issue and has become a significant public health concern worldwide (Shen et al. 2018). Many Gram-negative bacteria produce extended-spectrum β-lactamase (ESBL) that can inactivate cephalosporins, monobactams, and penicillins. The production of ESBLs by bacteria poses a significant challenge in treating infections caused by these organisms since it renders many commonly used antibiotics less effective or even ineffective (Onduru et al. 2021). Meanwhile, carbapenemase-producing E. coli isolated from various countries contain unique carbapenemase genes, and their prevalence exhibits considerable global variation. In one study, eight carbapenemases, including KPC, IMP, NDM, VIM, OXA, FRI, GES, and IMI, were identified among the carbapenem-resistant E. coli (CRECs) analyzed, with NDM being the most predominant (52.15%), followed by OXA (30.09%) and KPC (14.72%) (Li et al. 2024b).

The resurgence of polymyxins as a last-resort treatment in the early 21st century was driven by the rising prevalence of MDR bacteria and the lack of new antibiotics (Rodríguez-Santiago et al. 2021). As a cationic polypeptide, it possesses a positive charge that enables it to bind to lipid A’s negatively charged phosphate group in the lipopolysaccharide (LPS) of bacterial outer membranes. This interaction leads to membrane destabilization, causing leakage of cytoplasmic contents and eventual cell lysis (Ling et al. 2020). However, as with other antibiotics, bacteria have developed resistance to polymyxins. Numerous studies have reported colistin-resistant bacteria’s emergence and global spread (Karki et al. 2021).

In 2015, Liu et al. (2016) discovered a plasmid-borne mcr-1 gene that confers polymyxin resistance. An initial report identified mcr-1 in E. coli and Klebsiella pneumoniae isolates from China’s food, animals, and humans (Rodríguez-Santiago et al. 2021). An initial report identified mcr-1 in E. coli and K. pneumoniae isolates from food, animals, and humans in China (Rodríguez-Santiago et al. 2021). The mobile colistin resistance (mcr) genes carried by plasmids have become increasingly widespread worldwide. As of now, 10 mcr variants (mcr-1 to mcr-10) have been identified in over 60 countries (Yin et al. 2022), with mcr-1, mcr-3, mcr-5, mcr-9, and mcr-10 variants being the most prevalent in mcr-positive E. coli (Dadashi et al. 2022). E. coli is the most common host for mcr gene, which are frequently carried by a variety of plasmids, including IncI2, IncHI2 and IncX4 (Shen et al. 2018; Xie et al. 2022). Research has indicated that the ease of horizontal transfer of mcr-family genes contributes to their global dissemination, increasing colistin tolerance in various bacterial species (Ling et al. 2020). Many studies have shown mcr-1-har-bouring E. coli belong to different sequence types, with a high prevalence of ST10, ST131, and ST155 (Wang et al. 2020; Shi et al. 2023; De La Cadena et al. 2023).

Countries across the world are now paying attention to the presence of mcr-1-positive E. coli producing ESBL, focusing on the environment, animals, water sources, and humans (Chotinantakul et al. 2022; Jamil et al. 2022; Lima et al. 2022; Shafiq et al. 2022; Bastidas-Caldes et al. 2023). However, very few studies have focused on the prevalence and characterization of ESBL-producing mcr-1-harbouring E. coli in clinical patients in Shandong Province; only a few studies have been focused on individuals living in rural areas (Bi et al. 2017; Berglund et al. 2018).

The current study aimed to investigate molecular epidemiology of polymyxin B resistance among ESBL-producing E. coli in the region and the presence of the mcr-1 gene. The findings could provide a theoretical evidence for the prevention and control of mcr-1 transmission.

Experimental

Materials and Methods

Clinical strains identification. In this study, 668 non-duplicate ESBL-producing E. coli strains isolated from clinical samples at Shandong Provincial Hospital between January and December 2018 were collected for analysis. VITEK® 2 compact system (bioMérieux, France) and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF/MS) were used to identify the isolates. Twenty-four polymyxin B-resistant clinical strains were obtained by colistin agar test (agar plates containing 4 μg/ml of colistin). Then, polymyxin B’s minimum inhibitory concentration (MIC) was determined using the broth microdilution method and interpreted following the guidelines of the Clinical and Laboratory Standards Institute (CLSI 2023). The Ethics Committee of Shandong Provincial Hospital approved the clinical data collection of patients.

Antimicrobial susceptibility testing (AST). Antimicrobial susceptibility testing was conducted utilizing the VITEK® 2 compact system for ampicillin (AMP), ampicillin/sulbactam (SAM), piperacillin/tazobactam (TZP), cefazolin (CZO), cefotetan (CTT), cefatriaxone (CRO), ceftazidime (CAZ), cefepime (FEP), aztreonam (ATM), ertapenem (ETP), ertapenem (IPM), amikacin (AK), gentamicin (CN), tobramycin (TOB), ciprofloxacin (CIP), ciprofloxacin (LEV), nitrofurantoin (NIT), and trimethoprim/sulfamethoxazole (STX), except for polymyxin B (POL), which was analyzed by broth microdilution. E. coli ATCC® 25922™ was used as the reference strain. Results of the susceptibility assays were interpreted using the breakpoints provided by the Clinical and Laboratory Standards Institute (CLSI 2023).

Conjugation assay and confirmation of resistance genes. Sodium azide-resistant E. coli J53 was utilized as the recipient bacterium while mcr-1-harboring strains were employed as the donor strains. Transconjugants carrying the mcr-1 resistance genes were selectively grown on China-blue agar plates supplemented with 4 μg/ml polymyxin B and 100 μg/ml sodium azide. To confirm the transfer of the mcr-1 gene, polymerase chain reaction (PCR) was performed using primers mcr-1-F (5’-AGTCCGTTTGTTCTTGTGGC-3’) and mcr-1-R (5’-AGATCCTTGGTCTCGGCTTG-3’).

Pulsed-field gel electrophoresis (PFGE). The protocols were as described previously (Ma et al. 2020; Tang et al. 2022). 19 clinical strains and 17 transconjugants were selected and analyzed using S1-PFGE to reveal the number and approximate size of the plasmids. The XbaI-digested fragments of Salmonella enterica serovar Braenderup (H9812) were used as size markers. The XbaI-digested genomic DNA samples were subjected to pulsed-field gel electrophoresis (PFGE) using a CHEF Mapper® XA apparatus manufactured by Bio-Rad (USA) to assess the genetic relatedness across the clinical strains. Gel-J software, version 2.0 (Heras et al. 2015), was used to compare the pulsed-field gel electrophoresis (PFGE) patterns. The pulsotypes were categorized into clusters based on a similarity threshold of 80% (Chen et al. 2019).

Whole-genome sequencing (WGS) and analyses. Clinical strains and transconjugants were enriched overnight in Luria-Bertani broth at 37°C. Genomic DNA was extracted from the enriched cultures using a commercially available bacterial DNA extraction kit (Aidlab Inc., China). Both concentration and purity of the extracted genomic DNA were assessed using a NanoDrop™ 1000 Spectrophotometer (Thermo Fisher Scientific, Inc., USA), and the genomic DNA was subjected to next-generation sequencing using an Illumina Hiseq platform at the Novogene (Novogene Co., Ltd., China). Illumina sequences were assembled de novo using the SPAdes v3.10 (Prjibelski et al. 2020).

Multilocus sequence typing (MLST) was analyzed in silico by using Abricate software. FimH typing was performed using the Center for Genomic Epidemiology website tool. Antibiotic resistance genes were identified using the ResFinder database (Zankari et al. 2012). Plasmid replicon types were identified using the Plasmid Finder database (Carattoli et al. 2014).

Genomic DNA of ECO-POL-29 with the most common sequence type among these clinical strains, was sequenced using the Illumina HiSeq platform (Novogene Co., Ltd., China) and a PacBio RSII sequencer (Biozeron Biological Technology Co., Ltd., China). Paired-end short Illumina reads were utilized to correct long PacBio reads using the proovread software. Subsequently, the corrected PacBio reads were assembled de novo using SMARTdenovo software.

Gene annotation of the sequence was initially completed using the Rapid Annotation using Subsystem Technology (RAST) online annotation platform for prokaryotic genomes. The gene annotation results were further refined by manual curation using the BLAST tool. The mobile element annotation was performed using the online IS finder database (Siguier et al. 2006).

Phylogenetic analysis of isolates. Single-nucleotide polymorphism (SNP) calling was conducted with Snippy v3.1, and recombinant variants were filtered out with ClonalFrameML v1.0 (Lu et al. 2019). Recombination-free SNPs were used to construct maximumlikelihood phylogenetic trees with the RAxML software (Stamatakis et al. 2014). The tree file was visualized and the annotated information was edited using iTOL v6 (Letunic and Bork 2024). Serotyping was performed using the EcoSP website tool (Dong et al. 2023).

Analysis of the plasmid harboring mcr-1. The sequences of five plasmids with high sequence similarity (93–97% query coverage and 98.35–99.14% identity) to pECO-POL-29_mcr1 were obtained through the BLASTn analysis. Then these five plasmids were compared with pECO-POL-29_mcr1 by a BLAST Ring Image Generator tool (Alikhan et al. 2011), including pSLy21 (swine, United States) (GenBank accession no. CP016405), pAH01-2 (chicken, China) (GenBank accession no. CP055253), pRIVM_C029515_2 (Homo sapiens, Netherlands) (GenBank accession no. CP081339), pmcr1_IncI2 (Homo sapiens, China) (GenBank accession no. KU761326), and pEC26-mcr-1.9 (Homo sapiens, China) (GenBank accession no. MG946761).

Results

Clinical isolates and antimicrobial susceptibility. Twenty-four polymyxin B-resistant strains were isolated from 668 non-duplicate ESBL-producing E. coli strains, with a positivity rate of 3.59%. Of the 668, 19 clinical strains carried the mobile colistin resistance gene mcr-1, with a positivity rate of approximately 2.8%. Clinical data of patients infected with these strains were shown in Table I.

Table I Clinical data of patients infected with mcr-1 -positive Escherichia coli isolates producing extended-spectrum β-lactamase.

Clinical isolates	Gender	Department	Specimen	Clinical diagnosis	Antimicrobial therapy	
ECO-POL-4	male	Intensive Care Unit	urine	aspiration pneumonia, ARDS	cefoperazone sodium/sulbactam sodium	
ECO-POL-6	female	Orthopedic Oncology and Joint Surgery Department	urine	recurrence of intradural teratoma	levofloxacin, levofloxacin sodium chloride, ceftazidime	
ECO-POL-7	female	Gastroenterology Department	bile	postoperative recurrence of pancreatic cancer	piperacillin/tazobactam	
ECO-POL-8	female	Emergency Department	secretion	acute purulent appendicitis	meropenem sodium/sulbactam sodium	
ECO-POL-9	female	Urology Department	urine	ureteral stone	–	
ECO-POL-10	female	Minimally Invasive Urology Department	urine	left renal abscess	cefotaxime	
ECO-POL-15	male	Minimally Invasive Urology Department	urine	post-prostatectomy for prostate cancer	cefotaxime	
ECO-POL-17	male	General Surgery Department	puncture fluid	post-pancreaticoduodenectomy	cefoperazone sodium/sulbactam sodium, imipenem/cilastatin	
ECO-POL-18	female	Urology Department	urine	bilateral ureteral stones	meropenem	
ECO-POL-19	female	Minimally Invasive Urology Department	urine	urinary retention	ceftriaxone sodium	
ECO-POL-21	male	Respiratory Department	sputum	severe pneumonia	fluconazole, cefoperazone sodium/sulbactam sodium, carbapenem	
ECO-POL-22	male	Intensive Care Unit	drainage fluid	open pelvic fracture with ARDS	imipenem/cilastatin	
ECO-POL-23	male	Intensive Care Unit	bronchoalveolar lavage fluid	septic shock	carbapenem, imipenem/cilastatin, daptomycin, meropenem + polymyxin B, vancomycin + fluconazole	
ECO-POL-24	male	Trauma and Orthopedics	secretion	traumatic hemorrhagic shock	cefotaxime, clindamycin hydrochloride	
ECO-POL-25	female	Spinal Surgery Department	puncture fluid	postoperative infection following lumbar spine surgery	penicillin, cefepime + levofloxacin, imipenem/cilastatin	
ECO-POL-27	male	Thoracic Surgery Department	sputum	post-chemotherapy for right upper lobe lung cancer	cefuroxime	
ECO-POL-28	male	Intensive Care Unit	pus	gastrointestinal bleeding, gastric ulcer	cefoperazone sodium/sulbactam sodium, linezolid	
ECO-POL-29	male	General Surgery Department	ascitic puncture fluid	post-appendectomy for appendiceal cancer, colonic lesion	imipenem/cilastatin	
ECO-POL-30	male	Trauma and Orthopedics	secretion	hemorrhagic shock	cefotaxime	

All 19 clinical strains were resistant to ampicillin, cefazolin, ceftriaxone, ciprofloxacin, levofloxacin, and polymyxin B. However, there was high sensitivity to piperacillin/tazobactam, with a sensitivity rate of 89.5% (17/19). The rate of aztreonam resistance was 84.2% (16/19). One of the 19 strains was resistant to ertapenem and imipenem. The rate of amikacin resistance was 26.3% (5/19). The rate of resistance to gentamicin was 73.7% (14/19). The resistance rates to nitrofurantoin and trimethoprim/sulfamethoxazole were 5.3% (1/19) and 73.7% (14/19), respectively.

Prior to conjugation, the recipient strain, J53, was sensitive to all the tested drugs. All transconjugants were sensitive to cefotetan, ertapenem, imipenem, amikacin, ciprofloxacin, levofloxacin, and nitrofurantoin but resistant to polymyxin B. Of the 17 transconjugants, seven were resistant to ampicillin, ampicillin/sulbactam, cefazolin, and ceftriaxone, and six were resistant to ceftazidime and cefepime. Sensitivity to piperacillin/tazobactam and tobramycin was high, with only one strain being intermediate or resistant out of the 17. Resistance rates to gentamicin and trimethoprim/sulfamethoxazole were both 17.6%. The results of the antibiotic susceptibility tests are presented in Table II.

Table II Antibiotic susceptibility testing of Escherichia coli clinical strains, transconjugants, and the recipient bacterium J53AziR.

Isolates	MIC (μg/ml)	
AMP	SAM	TZP	CZO	CTT	CRO	CAZ	FEP	ATM	ETP	IPM	AK	CN	TOB	CIP	LEV	NIT	STX	POL	
ECO-POL-4	≥ 32	≥ 32	≥ 128	≥ 64	≥ 64	≥ 64	≥ 64	≥ 64	16	≥ 8	≥ 16	≤ 2	≤ 1	≤ 1	≥ 4	≥ 8	32	≥ 320	8	
ECO-POL-6	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	32	≥ 64	≤ 1	16	≤ 0.5	≤ 1	≤ 2	≥ 16	≥ 16	≥ 4	≥ 8	64	40	4	
ECO-POL-7	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	16	≥ 64	16	≤ 0.5	≤ 1	≤ 2	≥ 16	≥ 16	≥ 4	≥ 8	≤ 16	≥ 320	4	
ECO-POL-8	≥ 32	8	≤ 4	≥ 64	≤ 4	≥ 64	16	2	16	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≥ 4	≥ 8	32	≥ 320	4	
ECO-POL-9	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	4	2	16	≤ 0.5	≤ 1	≤ 2	≥ 16	8	≥ 4	≥ 8	≤ 16	≥ 320	4	
ECO-POL-10	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	16	≥ 64	16	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≥ 4	4	32	≥ 320	4	
ECO-POL-15	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	16	≥ 64	≥ 64	≤ 0.5	≤ 1	≥ 64	≥ 16	≥ 16	≥ 4	≥ 8	64	≤ 20	8	
ECO-POL-17	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	16	≥ 64	≥ 64	≤ 0.5	≤ 1	≥ 64	≥ 16	≥ 16	≥ 4	≥ 8	32	≤ 20	4	
ECO-POL-18	≥ 32	≥ 32	16	≥ 64	≤ 4	≥ 64	16	32	≥ 64	≤ 0.5	≤ 1	16	≥ 16	≥ 16	≥ 4	≥ 8	≤ 16	≥ 320	≥ 64	
ECO-POL-19	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	16	32	32	≤ 0.5	≤ 1	≥ 64	≥ 16	≥ 16	≥ 4	≥ 8	64	≥ 320	4	
ECO-POL-21	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≥ 16	≥ 16	≥ 4	≥ 8	64	≥ 320	4	
ECO-POL-22	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	≤ 1	2	16	≤ 0.5	≤ 1	≥ 64	≥ 16	≥ 16	≥ 4	≥ 8	64	≥ 320	4	
ECO-POL-23	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≥ 16	≥ 16	≥ 4	≥ 8	64	≥ 320	4	
ECO-POL-24	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	16	32	16	≤ 0.5	≤ 1	≤ 2	≥ 16	8	≥ 4	≥ 8	≤ 16	≥ 320	4	
ECO-POL-25	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	4	≥ 64	32	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≥ 4	≥ 8	128	≤ 20	4	
ECO-POL-27	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	4	2	16	≤ 0.5	≤ 1	≤ 2	≥ 16	8	≥ 4	≥ 8	64	≤ 20	4	
ECO-POL-28	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	≤ 1	2	4	≤ 0.5	≤ 1	≥ 64	≥ 16	≥ 16	≥ 4	≥ 8	64	≥ 320	4	
ECO-POL-29	≥ 32	8	≤ 4	≥ 64	≤ 4	≥ 64	4	2	16	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≥ 4	≥ 8	32	≥ 320	4	
ECO-POL-30	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	16	32	16	≤ 0.5	≤ 1	≤ 2	≥ 16	8	≥ 4	≥ 8	≤ 16	≥ 320	4	
J4	≥ 32	≥ 32	64	≥ 64	≤ 4	≥ 64	16	≥ 64	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	8	
J6	8	8	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J7	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	16	≥ 64	≥ 64	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≥ 320	4	
J8	4	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J9	4	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	8	
J10	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	16	≥ 64	≥ 64	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J15	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	4	2	16	≤ 0.5	≤ 1	≤ 2	8I	≥ 16	≤ 0.25	≤ 0.25	32	≤ 20	16	
J18	8	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	8	
J19	8	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J21	8	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J22	8	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J23	8	4	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J24	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	≥ 64	8	≥ 64	≤ 0.5	≤ 1	≤ 2	≥ 16	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≥ 320	8	
J25	≥ 32	16	≤ 4	≥ 64	≤ 4	≥ 64	16	≥ 64	≥ 64	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	4	
J27	8	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	8	
J29	4	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	8	
J30	≥ 32	≥ 32	≤ 4	≥ 64	≤ 4	≥ 64	≥ 64	8	≥ 64	≤ 0.5	≤ 1	≤ 2	≥ 16	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≥ 320	8	
J53AziR	8	≤ 2	≤ 4	≤ 4	≤ 4	≤ 1	≤ 1	≤ 1	≤ 1	≤ 0.5	≤ 1	≤ 2	≤ 1	≤ 1	≤ 0.25	≤ 0.25	≤ 16	≤ 20	0.125	
1 MIC – minimal inhibitory concentrations, AMP – ampicillin, SAM – ampicillin/sulbactam, TZP – piperacillin/tazobactam,

1 CZO – cefazolin, CTT – cefotetan, CRO – ceftriaxone, CAZ – ceftazidime, FEP – cefepime, ATM – aztreonam, ETP – ertapenem,

1 IPM – imipenem, AK – amikacin, CN – gentamicin, TOB – tobramycin, CIP – ciprofloxacin, LEV – levofloxacin,

1 NIT – nitrofurantoin, STX – trimethoprim/sulfamethoxazole, POL – polymyxin B

Plasmid carriage status of the clinical strains and transconjugants. S1-PFGE was used to analyze the number and approximate size of plasmids. Among the 19 strains of mcr-1-positive polymyxin B-resistant ESBL-producing E. coli, there were different numbers of plasmids, ranging from 2 to 6, with sizes roughly estimated between 20 kb and 300 kb. Seventeen transcon-jugants were obtained using conjugation assays. The transconjugants produced 1–4 plasmids with plasmid sizes ranging from approximately 20 kb to 300 kb (Table III).

Table III The number, approximate size and replicon types of plasmids in 19 strains of mcr-1-positive polymyxin B-resistant Escherichia coli producing extended-spectrum β-lactamase and 17 transconjugants.

Isolates	Plasmid numbers	Size (kbp)	Replicon types	
ECO-POL-4	3	~ 110, ~ 65, ~ 45	ColRNAI, IncFIA, IncFIC, IncFIC, IncX3	
ECO-POL-6	6	~ 210, ~ 140, ~ 110, ~ 104.5, ~ 78.2, ~ 65	Col, ColRNAI, ColpVC, IncFIB, IncFII, IncHI1B, IncHI2A,	
			IncHI2, IncI1, IncI2, RepA, p0111	
ECO-POL-7	5	~ 220, ~ 138, ~ 100, ~ 65, ~ 54	IncFIB, IncFIC, IncFII, IncHI2A, IncHI2, IncI2a,	
			RepA_1_pKPC-CAV1321	
ECO-POL-8	4	~ 140, ~ 100, ~ 33, ~ 21	IncFIB, IncFIC, IncX4	
ECO-POL-9	3	~ 140, ~ 100, < 20	Col156, Col440II, ColRNAI, IncB/O/K/Z, IncFIA, IncFIB1,	
			IncFIC, IncI2, IncX1, IncY	
ECO-POL-10	4	~ 140, ~ 65, ~ 34, ~ 30	Col, Col156, Col8282, IncFIB, IncFII, IncI2, IncX4	
ECO-POL-15	3	~ 265, ~ 65, ~ 45	Col, IncFII, IncHI2A, IncHI2, IncN, IncX1, RepA	
ECO-POL-17	3	~ 104.5, ~ 65, ~ 48	Col440I, ColRNAI, IncI1, IncI2, IncN, IncR, IncX1	
ECO-POL-18	3	~ 110, ~ 104.5, ~ 65	IncFIB, IncFIC, IncHI2A, IncHI2, IncI2, RepA	
ECO-POL-19	3	~ 265, ~ 140, ~ 48	Col156, Col440I, IncFIB, IncFIC, IncFII, IncHI2A, IncHI2, RepA	
ECO-POL-21	2	~ 140, ~ 65	Col156, IncFIB, IncFIC, IncI2	
ECO-POL-22	3	~ 180, ~ 110, ~ 65	Col440II, ColRNAI, IncFIB, IncFIC, IncI2, IncI	
ECO-POL-23	2	~ 140, ~ 65	Col156, IncFIB, IncFIC, IncI2	
ECO-POL-24	2	~ 140, ~ 65	IncFIB, IncI2	
ECO-POL-25	3	~ 140, ~ 95, ~ 65	IncB/O/K/Z, IncFIB, IncI2, p0111	
ECO-POL-27	6	~ 140, ~ 130, ~ 95, ~ 70, ~ 65, ~ 43	Col(MG828), ColE10, ColRNAI, IncB/O/K/Z, IncFIB, IncFIC,	
			IncFII, IncI2, IncX1, p0111	
ECO-POL-28	2	~ 90, ~ 28	Col, Col156, Col440I, ColRNAI, IncFIA, IncFIB, IncI, IncX1	
ECO-POL-29	4	128.293, 83.176, 62.869, 30.619	IncFIB/IncFIC, IncFIC, IncI2, IncX4	
ECO-POL-30	2	~ 130, ~ 60	IncFIB, IncI2	
J4	1	~ 65	IncI2	
J6	1	~ 65	IncI2	
J7	4	~ 103, ~ 100, ~ 65, ~ 60	IncI2, IncFIB, IncFII, IncFIC	
J8	2	~ 100, ~ 33	IncFIB, IncX4	
J9	1	~ 65	IncI2	
J10	3	~ 65, ~ 33, ~ 30	IncI2, IncX4	
J15	1	~ 255	IncHI2, RepA, IncHI2A, IncN	
J18	1	~ 65	IncI2	
J19	1	~ 255	IncHI2A, IncHI2, RepA	
J21	1	~ 65	IncI2	
J22	1	~ 65	IncI2, ColRNAI	
J23	1	~ 65	IncI2	
J24	2	~ 140, ~ 65	IncI2, IncFIB	
J25	1	~ 65	IncI2	
J27	1	~ 65	IncI2	
J29	2	~ 65, ~ 33	IncI2, IncX4	
J30	2	~ 140, ~ 65	IncI2, IncFIB	

Molecular typing and genotyping of the phylogenetic groups. Of the 19 strains of mcr-1 -positive polymyxin B-resistant E. coli-producing ESBL, there were 14 sequence types (STs). The most common ST was ST155 (n = 4, 21%), followed by ST162 (n = 2), ST393 (n = 2), ST457, ST88, ST450, ST617, ST744, ST354, ST359, ST93, ST224, ST69, and ST167. ECO-POL-21 and ECO-POL-23 shared the same serotype (O29:H12) and ST. Similarly, ECO-POL-24 and ECO-POL-30 had the same serotype (O15: H1) and ST (Fig. 1).

Fig. 1. A dendrogram was constructed using the similarity of pulsed-field gel electrophoresis (PFGE) patterns across the 19 strains of mcr-1 -positive polymyxin-B-resistant extended-spectrum β-lactamase-producing Escherichia coli. Names of the isolates, sequence types, serotypes, and fimH types are displayed on the right side.

Based on the PFGE profiles, only three distinct clusters, identified as clone groups A to C, were observed, indicating a lower homology among the 19 isolates. According to the PFGE profile and MLST analysis, there was homology between strains 24 and 30, both of which belong to ST393. The ST162 isolates showed a high degree of homology. The 19 strains were distributed across different branches, suggesting considerable genetic diversity among the strains, with variations in the overall SNPs. According to branch length, ECO-POL-24 and ECO-POL-30 exhibited the most closely related evolutionary lineage (Fig. 2). The predominant replicon type of transconjugants was IncI2 (n = 14) followed by IncFIB (n = 4) (Table III).

Fig. 2. Phylogenetic tree distribution of 19 strains of mcr-1-positive polymyxin-B-resistant extended-spectrum β-lactamase-producing Escherichia coli.

Distribution of drug-resistance genes in clinical strains and transconjugants. Next-generation sequencing and analysis revealed that all 19 clinical strains harbored mcr-1. No other mcr variants were detected. Additionally, the drug resistance genes blactx-m and blatem-1 were detected. Variants detected for blactx-m resistance gene included blactx-m-104, blactx-m-132, blactx-m-14, blactx-m-27, blactx-m-55, blactx-m-64 and blactx-m-65. In this study, blactx-m was the most common ESBL resistance gene. The carriage rate of blactx-m in 19 clinical strains was as high as 100% (19/19), with the majority being blactx-m-55 (seven strains). Notably, ECO-POL-15 and ECO-POL-18 both harbored blactx-m-55 and blactx-m-65, while ECO-POL-10 harbored blactx-m-27 and blactx-m-132.

The carriage rate of blatem-1 was 36.8% (7/19) and blashv was not detected. Four strains carried blactx-m-64, blatem-1, and mcr-1 simultaneously.

Regarding aminoglycoside resistance, 14 isolates harbored APH(6)-Id, followed by APH(3’’)-Ib (13/19), AAC(3)-IV (8/19), and aadA2 (8/19). ECO-POL-19 carried the most aminoglycoside resistance genes, totaling to 8. Blandm-5 resistance gene was detected in the imipenem-resistant isolate (5.3%, 1/19). Notably, the plasmid-encoded fosfomycin resistance gene fosA3 was detected in ten isolates (52.6%, 10/19), whereas the fluoroquinolone gene QnrS1 was detected in three. In sulfonamide resistance genes, sul2 was the most prevalent. All of the 19 clinical strains and 17 transconjugants carried the β-lactamase resistance gene ampC.

DNA sequencing and PCR results showed that all 17 transconjugants carried the drug resistance gene mcr-1. ESBL genes such as blactx-m, were identified, with blactx-m-64 being the most prevalent, compared to other types such as blactx-m-55, blactx-m-65 and blactx-m-132. Nine transconjugants lacked any blactx-m. Five transconjugants carried blatem-1. Transconjugants co-harboring blactx-m, blatem-1, and mcr-1 comprised 29.4% of all the transconjugants, and blactx-m was included blactx-m-55, blactx-m-64 and blactx-m-65. The aminoglycoside resistance genes AAC(3)-IId and APH(6)-Id were relatively common among the 17 transconjugants. None of the transconjugants harbouring the carbapenemase or fos-fomycin genes was acquired. The results are collectively shown in Fig. 3.

Fig. 3. a) Detection of resistance genes in 19 mcr-1-positive clinical strains, b) detection of resistance genes in 17 transconjugants.

Analysis of the plasmid harboring mcr-1. The clinical strain ECO-POL-29 was subjected to wholegenome sequencing (WGS). One chromosome and four plasmids were present in the ECO-POL-29. The plasmid pECO-POL-29_mcr1 was an IncI2-type plasmid with a GC content of 42.7% and 88 open reading frames and the only resistance gene in this plasmid was mcr-1. There was no insertion sequence near mcr-1. The Type IV secretion system (T4SS) was detected in pECO-POL-29_mcr1. A transconjugant harboring pECO-POL-29_mcr1 was successfully obtained using the conjugation assay. Based on the BLASTn results for pECO-POL-29_mcr1, five plasmids from different countries, regions, and hosts were identified and selected. All five plasmids originated from E. coli and shared the same plasmid type IncI2. As indicated in Fig. 4, the main differences across the plasmids were distributed in region 1 (ranging from 7.32 kb to 10.19 kb) and region 2 (ranging from 57.8 kb to 58.9 kb). These regions contain genes associated with proteins of unknown function (Fig. 4).

Fig. 4. The external ring is the schematic representation of plasmid pECO-POL-29_mcr1 (GenBank accession number PP721192). Genes are marked with different colors based on their functional annotations. The internal five rings show a comparative analysis of mcr-1-harboring plasmids with pECO-POL-29_mcr1, including pSLy21 (sky blue), pAH01-2 (light purple), pRIVM_C029515_2 (light blue), pmcr1_IncI2 (light purple pink), and pEC26-mcr-1.9 (purple) (constructed using BRIG).

Nucleotide sequence accession numbers. The raw reads of all 19 isolates have been deposited in GenBank (BioProject PRJNA1104118). The complete sequence of pECO-POL-29_mcr1 has been deposited in GenBank under accession number PP721192.

Discussion

E. coli is a significant causative agent of severe clinical bacterial infections with increasing levels of antibiotic resistance(Arafi et al. 2023; Ibrahim et al. 2023). Polymyxins have re-entered our field of view as the last line of defense against infections caused by multidrugresistant (MDR) Gram-negative bacteria (Rhouma et al. 2023). Since the first detection of the mcr-1 gene in 2015, ten mcr variants have been identified to date. Among those, mcr-1 remains the most dominant in clinical isolates (Nang et al. 2022; Rodríguez-Santiago et al. 2022).

Among the 668 strains of ESBL-producing E. coli collected in this study, 24 were multidrug-resistant strains, resistant to polymyxin B, with a resistance rate of 3.59%. According to a review published in 2022, the global prevalence of polymyxins resistance was 1.6% (Binsker et al. 2022). The result of this study was higher than that of the previous review. Among the 668 strains, 19 carried the colistin resistance gene mcr-1, with a positivity rate of approximately 2.8%. According to a tudy in 2016, mcr-1 was detected in only two isolates (4.2%) of 48 ESBL-producing E. coli from a hospital in France (Robin et al. 2017). In a study on healthy children in Central South China, the prevalence of mcr-1 was 0.84%, with one positive sample of over 118 ESBL-producing E. coli (Liu et al. 2022). The 2.8% positivity rate was significantly lower than that of environmental strains (5.43%) (Fan et al. 2020).

It would be worth noting that the mcr gene was associated with other antimicrobial resistance genes, such as ESBLs, carbapenemases, fluoroquinolons, and fosfomycins, making the treatment of infections caused by these bacteria even more difficult (Tian et al. 2020; Rodríguez-Santiago et al. 2022; Di Francesco et al. 2023; Szmolka et al. 2023). In the present study, all 19 mcr-1-positive clinical strains carried blactx-m, with blactx-m-55 being the most predominant. Fosfomycin, an old bactericidal antibiotic, was re-evaluated for its potential in the combined therapy of multidrugresistant (MDR) Enterobacterale infections, particularly those caused by ESBL-producing strains (Loras et al. 2020). The positive rate of fosA3 in ESBL-producing E. coli was 3.0% (5/165) (Lee et al. 2012) to 11% (51/465) (Cao et al. 2017). In this study, fosA3 was detected in 52.6% of the 19 isolates, which was higher than that reported in previous studies. Despite detecting the co-existence of blandm-5 and mcr-1 genes in only one isolate in our study, it is noteworthy that such co-producing isolates have been increasingly reported. (Han et al. 2020; Chen et al. 2023).

Many studies have indicated that the STs of E. coli isolates carrying mcr-1 gene exhibit a high degree of diversity (Liu et al. 2018; Zhang et al. 2019; Vu Thi Ngoc et al. 2022). The STs in this study were also characterized by such diversity. According to a previous study, the most common ST among the isolates carrying the mcr-1 gene was ST10, also commonly found in E. coli isolates carrying ESBL from humans and animals (Shen et al. 2018). In our study, ST155 was the most common sequence type, which was consistent with the description of a previous study (Wang et al. 2020).

The mcr gene carried by plasmids is present in various types of plasmids, such as IncI2, IncHI2, IncX4, IncP, IncF, and uncommon IncY (Li et al. 2022; Protonotariou et al. 2022; Rodríguez-Santiago et al. 2022). The IncI2 plasmid is considered the most successful vector for horizontal dissemination of the mcr-1 gene, according to previous studies (Al Mana et al. 2022; He et al. 2022). However, in the Shanghai region, between 2012 and 2015, the predominant plasmid type among mcr-1-positive E. coli was IncX4, followed by IncI2. After 2015, IncI2 gradually surpassed IncX4 (Xie et al. 2022). Nevertheless, this result differed from those of the main types detected in Hong Kong (IncI1, IncX4, and IncN) (Chan et al. 2018). In the present study, the most common replicon types of transconjugants associated with mcr-1 were IncI2 and IncFIB.

Whole-genome sequencing of ECO-POL-29 revealed that the mcr-1 gene was located on an IncI2-type plasmid approximately 62.8 kb in size, and the plasmid did not carry any other resistance gene. The transconjugant carried only the mcr-1 gene. Additionally, the blactx-m-55 gene was located on IncFIB/IncFIC-type plasmids. Furthermore, there was no insertion sequence surrounding the mcr-1 gene. The size of the five plasmids included in the comparison analysis with pECO-POL-29_mcr1 ranged from 62 kb to 65 kb. Five plasmids and pECO-POL-29_mcr1 shared the same plasmid type IncI2 and high similarity backbone structure. The characteristics of plasmids from different countries, regions, and sources were found to be conserved. IncI2 type plasmids harboring mcr-1 gene played an important role in the rapid and efficient dissemination of mcr-1 gene. Thus, more attention should be paid to monitoring and controlling the horizontal transmission of mcr-1 mediated by IncI2 type plasmids.

The genetic background of mcr-1 was found to be highly diverse. Many studies have indicated that the mobilization of the mcr-1 gene is primarily associated with transposition mediated by ISApl1 (Shen et al. 2018; Ding et al. 2021; Zhou et al. 2021). ISApl1 is considered an element that can mobilize genes on plasmids and integrate them into the chromosome (Al Mana et al. 2022). It was first described in Acineto-bacillus pleuropneumoniae as a part of the IS30 family (Tegetmeyer et al. 2008). Previous studies have shown that the ISApl1-mcr-1-pap2-ISApl1 and ISApl1-mcr-1-pap2 structures carrying mcr-1 are involved in ISApl1 (Xie et al. 2022; Li et al. 2024a). Shen et al. showed that most mcr-1 genes are not adjacent to IS elements (Shen et al. 2018). A previous study had reported that 68.3% of mcr-1-positive E. coli strains completely lost ISApl1 (Huang et al. 2021). After analyzing our sequencing results, we found that the mcr-1 gene was not adjacent to ISApl1. Our findings collectively suggested that ISApl1 is not an essential mobile element for the transfer of mcr-1.

According to the data compiled in one study, ISApl1 is not commonly found in IncI2 and IncX4 plasmids (Li et al. 2024a). This description aligned with the genetic environment of pECO-POL-29_mcr1, since the plasmid type of pECO-POL-29_mcr1 is IncI2. Based on the result of pECO-POL-29_mcr1 genomic analysis, there was no mobile element surrounding mcr-1, suggesting that the mechanism by which the plasmid acquired mcr-1 would require further study. The T4SS on the pECO-POL-29_mcr1plasmid plays a crucial role in the spread of the mcr-1 gene, enabling its transmission as a plasmid.

Conclusions

In conclusion, our study uncovered the presence and spread of mcr-1-carrying ESBL-producing E. coli in the region and provided their genomic characteristics. As we increasingly recognise the importance of “One Health” taking more measures to monitor and stem the prevalence and transmission of mcr-1 would be essential.

Acknowledgments

We thank Mingju Hao from Department of Clinical Laboratory Medicine, The First Affiliated Hospital of Shandong First Medical University, Jinan, Shandong, China and Zhen Liu from Novogene(Beijing, China) for the help and support.

Availability of data and material

The data presented in this study are available on reasonable request from the corresponding author.

Author contributions

Yingying Hao and Ping Li contributed to experiment conception and design. Yue Liu and Qian Wang conducted bioinformatics analysis and the wrote the paper. Ting Qi and Meng Zhang performed data analysis. Ran Chen and Zaifeng Si carried out the bacteria identification. Jinmei Li, Yan Jin and Qingbing Xu prepared the Tables and Figures. The corresponding author is responsible for submitting a competing interests statement on behalf of all authors of the paper.

Conflict of interest

The authors do not report any financial or personal connections with other persons or organizations, which might negatively affect the contents of this publication and/or claim authorship rights to this publication.
==== Refs
Literature

Al Mana H Johar AA Kassem II Eltai NO. Transmissibility and persistence of the plasmid-borne mobile colistin resistance gene, mcr-1, harbored in poultry-associated E. coli . Antibiotics . 2022 NaN ; 11 ( 6 ): 774 . 10.3390/antibiotics11060774 35740180
Alikhan NF Petty NK Ben Zakour NL Beatson SA. BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons . BMC Genomics . 2011 NaN ; 12 : 402 . 10.1186/1471-2164-12-402 21824423
Arafi V Hasani A Sadeghi J Varshochi M Poortahmasebi V Hasani A Hasani R. Uropathogenic Escherichia coli endeavors: An insight into the characteristic features, resistance mechanism, and treatment choice . Arch Microbiol . 2023 NaN ; 205 ( 6 ): 226 . 10.1007/s00203-023-03553-5 37156886
Bastidas-Caldes C Cisneros-Vásquez E Zambrano A Mosquera-Maza A Calero-Cáceres W Rey J Yamamoto Y Yamamoto M Calvopiña M de Waard JH. Co-harboring of beta-lactamases and mcr-1 genes in Escherichia coli and Klebsiella pneumoniae from healthy carriers and backyard animals in rural communities in Ecuador . Antibiotics . 2023 NaN ; 12 ( 5 ): 856 . 10.3390/antibiotics12050856 37237759
Berglund B Chen B Tärnberg M Sun Q Xu L Welander J Li Y Bi Z Nilsson M Nilsson LE. Characterization of extended-spectrum β-lactamase-producing Escherichia coli harboring mcr-1 and toxin genes from human fecal samples from China . Future Microbiol . 2018 NaN ; 13 : 1647 – 1655 . 10.2217/fmb-2018-0242 30489158
Bi Z Berglund B Sun Q Nilsson M Chen B Tärnberg M Ding L Stålsby Lundborg C Bi Z Prevalence of the mcr-1 colistin resistance gene in extended-spectrum ß-lactamase-producing Escherichia coli from human faecal samples collected in 2012 in rural villages in Shandong Province, China . Int J Antimicrob Agents . 2017 NaN ; 49 ( 4 ): 493 – 497 . 10.1016/j.ijantimicag.2016.12.018 28263896
Binsker U Käsbohrer A Hammerl JA. Global colistin use: A review of the emergence of resistant Enterobacterales and the impact on their genetic basis . FEMS Microbiol Rev . 2022 NaN ; 46 ( 1 ): fuab049 . 10.1093/femsre/fuab049 34612488
Cao XL Shen H Xu YY Xu XJ Zhang ZF Cheng L Chen JH Arakawa Y. High prevalence of fosfomycin resistance gene fosA3 in blactx-m-harbouring Escherichia coli from urine in a Chinese tertiary hospital during 2010–2014 . Epidemiol Infect . 2017 NaN ; 145 ( 4 ): 818 – 824 . 10.1017/s0950268816002879 27938421
Carattoli A Zankari E García-Fernández A Voldby Larsen M Lund O Villa L Møller Aarestrup F Hasman H. In silico detection and typing of plasmids using Plasmid Finder and plasmid multilocus sequence typing . Antimicrob Agents Chemother . 2014 NaN ; 58 ( 7 ): 3895 – 3903 . 10.1128/aac.02412-14 24777092
Chan WS Au CH Ho DN Chan TL Ma ES Tang BS. Prospective study on human fecal carriage of Enterobacteriaceae possessing mcr-1 and mcr-2 genes in a regional hospital in Hong Kong . BMC Infect Dis . 2018 NaN ; 18 ( 1 ): 81 . 10.1186/s12879-018-2987-y 29439654
Chen J Wang D Ding Y Zhang L Li X. Molecular epidemiology of plasmid-mediated fosfomycin resistance gene determinants in Klebsiella pneumoniae carbapenemase-producing Klebsiella pneumoniae isolates in China . Microb Drug Resist . 2019 NaN ; 25 ( 2 ): 251 – 257 . 10.1089/mdr.2018.0137 30113251
Chen R Wang G Wang Q Zhang M Wang Y Wan Z Si Z Bai Y Song Z Lu X Antimicrobial resistance and molecular epidemiology of carbapenem-resistant Escherichia coli from urinary tract infections in Shandong, China . Int Microbiol . 2023 NaN ; 26 ( 4 ): 1157 – 1166 . 10.1007/s10123-023-00369-7 37145384
Chotinantakul K Chusri P Okada S. Detection and characterization of ESBL-producing Escherichia coli and additional co-existence with mcr genes from river water in northern Thailand . PeerJ . 2022 NaN ; 10 : e14408 . 10.7717/peerj.14408 36405029
CLSI . Performance standards for antimicrobial susceptibility testing . 33th CLSI supplement M100 . Wayne (USA) : Clinical and Laboratory Standards Institute ; 2023 .
Dadashi M Sameni F Bostanshirin N Yaslianifard S Khosravi-Dehaghi N Nasiri MJ Goudarzi M Hashemi A Hajikhani B. Global prevalence and molecular epidemiology of mcr-mediated colistin resistance in Escherichia coli clinical isolates: A systematic review . J Glob Antimicrob Resist . 2022 NaN ; 29 : 444 – 461 . 10.1016/j.jgar.2021.10.022 34788692
De La Cadena E Mahecha M Velandia AM García-Betancur JC Rojas LJ Porras J Pallares C Villegas MV. Identification of mcr-1 genes and characterization of resistance mechanisms to colistin in Escherichia coli isolates from Colombian hospitals . Antibiotics . 2023 NaN ; 12 ( 3 ): 488 . 10.3390/antibiotics12030488 36978355
Di Francesco A Salvatore D Sakhria S Bertelloni F Catelli E Ben Yahia S Tlatli A. Colistin resistance genes in broiler chickens in Tunisia . Animals . 2023 NaN ; 13 ( 8 ): 1409 . 10.3390/ani13081409 37106971
Ding Y Saw WY Tan LWL Moong DKN Nagarajan N Teo YY Seedorf H. Extended-spectrum β-lactamase-producing and mcr-1-positive Escherichia coli from the gut microbiota of healthy Singaporeans . Appl Environ Microbiol . 2021 NaN ; 87 ( 20 ): e0048821 . 10.1128/aem.00488-21 34347523
Dong A Liu C Hua X Yu Y Guo Y Wang D Liu X Chen H Wang H Zhu L. Bioinformatic analysis of structures and encoding genes of Escherichia coli surface polysaccharides sheds light on the heterologous biosynthesis of glycans . BMC Genomics . 2023 NaN ; 24 ( 1 ): 168 . 10.1186/s12864-023-09269-6 37016299
Fan R Li C Duan R Qin S Liang J Xiao M Lv D Jing H Wang X. Retrospective screening and analysis of mcr-1 and blandm in Gramnegative bacteria in China, 2010–2019 . Front Microbiol . 2020 NaN ; 11 : 121 . 10.3389/fmicb.2020.00121 32117144
Han H Liu W Cui X Cheng X Jiang X. Co-existence of mcr-1 and bla NDM-5 in an Escherichia coli strain isolated from the pharmaceutical industry, WWTP . Infect Drug Resist . 2020 NaN ; 13 : 851 – 854 . 10.2147/idr.S245047 32214831
He K Li W Zhao B Xu H Pan Y He D Hu G Wu H Yuan L. Spreading advantages of coresident plasmids blactx-m-bearing IncFII and mcr-1-bearing IncI2 in Escherichia coli . Microbiol Spectr . 2022 NaN ; 10 ( 1 ): e0170621 . 10.1128/spectrum.01706-21 35171014
Heras J Domínguez C Mata E Pascual V Lozano C Torres C Zarazaga M. GelJ--a tool for analyzing DNA fingerprint gel images . BMC Bioinformatics . 2015 NaN ; 16 : 270 . 10.1186/s12859-015-0703-0 26307353
Huang S Wang S Li Y Fang M Kou Z Chen B Xu L Bi Z Xu H Chi X Prevalence and transmission of mobilized colistin resistance (mcr-1) gene positive Escherichia coli in healthy rural residents in Shandong province, China . Microbiol Res . 2021 NaN ; 253 : 126881 . 10.1016/j.micres.2021.126881 34592562
Ibrahim DR Dodd CER Stekel DJ Meshioye RT Diggle M Lister M Hobman JL. Multidrug-resistant ESBL-producing E. coli in clinical samples from the UK . Antibiotics . 2023 NaN ; 12 ( 1 ): 169 . 10.3390/antibiotics12010169 36671370
Jamil A Zahoor MA Nawaz Z Siddique AB Khurshid M. Genetic diversity of Escherichia coli coharboring mcr-1 and extended spectrum beta-lactamases from poultry . Biomed Res Int . 2022 NaN ; 2022 : 8224883 . 10.1155/2022/8224883 36246985
Karki D Dhungel B Bhandari S Kunwar A Joshi PR Shrestha B Rijal KR Ghimire P Banjara MR. Antibiotic resistance and detection of plasmid mediated colistin resistance mcr-1 gene among Escherichia coli and Klebsiella pneumoniae isolated from clinical samples . Gut Pathog . 2021 NaN ; 13 ( 1 ): 45 . 10.1186/s13099-021-00441-5 34225805
Lee SY Park YJ Yu JK Jung S Kim Y Jeong SH Arakawa Y. Prevalence of acquired fosfomycin resistance among extended-spectrum β-lactamase-producing Escherichia coli and Klebsiella pneumoniae clinical isolates in Korea and IS26-composite transposon surrounding fosA3 . J Antimicrob Chemother . 2012 NaN ; 67 ( 12 ): 2843 – 2847 . 10.1093/jac/dks319 22893681
Letunic I Bork P. Interactive Tree of Life (iTOL) v6: Recent updates to the phylogenetic tree display and annotation tool . Nucleic Acids Res . 2024 NaN ; 52 ( Wl ): W78 – W82 . 10.1093/nar/gkae268 38613393
Li Q Qian C Zhang X Zhu T Shi W Gao M Feng C Xu M Lin H Lin L Colistin resistance and molecular characterization of the genomes of mcr-1-positive Escherichia coli clinical isolates . Front Cell Infect Microbiol . 2022 NaN ; 12 : 854534 . 10.3389/fcimb.2022.854534 35601104
Li W He Z Di W Xu W Li Y Sun B. Transposition mechanism of ISApl1-the determinant of colistin resistance dissemination . Anti-microb Agents Chemother . 2024a NaN ; 68 ( 3 ): e0123123 . 10.1128/aac.01231-23
Li Y Sun X Dong N Wang Z Li R. Global distribution and genomic characteristics of carbapenemase-producing Escherichia coli among humans, 2005–2023 . Drug Resist Updat . 2024b NaN ; 72 : 101031 . 10.1016/j.drup.2023.101031 38071860
Lima T Loureiro D Henriques A Ramos F Pomba C Domingues S da Silva GJ. Occurrence and biological cost of mcr-1-carrying plasmids co-harbouring beta-lactamase resistance genes in zoonotic pathogens from intensive animal production . Antibiotics . 2022 NaN ; 11 ( 10 ): 1356 . 10.3390/antibiotics11101356 36290014
Ling Z Yin W Shen Z Wang Y Shen J Walsh TR. Epidemiology of mobile colistin resistance genes mcr-1 to mcr-9 . J Antimicrob Chemother . 2020 NaN ; 75 ( 11 ): 3087 – 3095 . 10.1093/jac/dkaa205 32514524
Liu X Li X Yang AW Tang B Jian ZJ Zhong YM Li HL Li YM Yan Q Liang XH Community fecal carriage and molecular epidemiology of extended-spectrum β-lactamase-and carbapene-mase-producing Escherichia coli from healthy children in the Central South China . Infect Drug Resist . 2022 NaN ; 15 : 1601 – 1611 . 10.2147/idr.S357090 35418762
Liu X Wang Y Cui L Li Y Xue F Liu J Wang L Shen Y Lv Y. A retrospective study on mcr-1 in clinical Escherichia coli and Klebsiella pneumoniae isolates in China from 2007 to 2016 . J Antimicrob Chemother . 2018 NaN ; 73 ( 7 ): 1786 – 1790 . 10.1093/jac/dky092 29912359
Liu YY Wang Y Walsh TR Yi LX Zhang R Spencer J Doi Y Tian G Dong B Huang X Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: A microbiological and molecular biological study . Lancet Infect Dis . 2016 NaN ; 16 ( 2 ): 161 – 168 . 10.1016/s1473-3099(15)00424-7 26603172
Loras C Mendes AC Peixe L Novais Â Alós JI. Escherichia coli resistant to fosfomycin from urinary tract infections: Detection of the fosA3 gene in Spain . J Glob Antimicrob Resist . 2020 NaN ; 21 : 414 – 416 . 10.1016/j.jgar.2020.01.023 32061811
Lu X Zeng M Xu J Zhou H Gu B Li Z Jin H Wang X Zhang W Hu Y Epidemiologic and genomic insights on mcr-1-harbour-ing Salmonella from diarrhoeal outpatients in Shanghai, China, 2006–2016 . EBioMedicine . 2019 NaN ; 42 : 133 – 144 . 10.1016/j.ebiom.2019.03.006 30905850
Ma J Wang J Feng J Liu Y Yang B Li R Bai L He T Wang X Yang Z. Characterization of three porcine Acinetobacter towneri strains coharboring tet(X3) and blaoxa-58 . Front Cell Infect Microbiol . 2020 NaN ; 10 : 586507 . 10.3389/fcimb.2020.586507 33363052
Nang SC Li M Harper M Mandela E Bergen PJ Rolain JM Zhu Y Velkov T Li J. Polymyxin causes cell envelope remodelling and stress responses in mcr-1-harbouring Escherichia coli . Int J Antimicrob Agents . 2022 NaN ; 59 ( 2 ): 106505 . 10.1016/j.ijantimicag.2021.106505 34954369
Onduru OG Mkakosya RS Aboud S Rumisha SF. Genetic determinants of resistance among ESBL-producing Enterobacteriaceae in Community and Hospital Settings in East, Central, and Southern Africa: A systematic review and meta-analysis of prevalence . Can J Infect Dis Med Microbiol . 2021 NaN ; 2021 : 5153237 . 10.1155/2021/5153237 34122680
Prjibelski A Antipov D Meleshko D Lapidus A Korobeynikov A. Using SPAdes De Novo Assembler . Curr Protoc Bioinformatics . 2020 NaN ; 70 ( 1 ): e102 . 10.1002/cpbi.102 32559359
Protonotariou E Meletis G Malousi A Kotzamanidis C Tychala A Mantzana P Theodoridou K Ioannidou M Hatzipantelis E Tsakris A First detection of mcr-1-producing Escherichia coli in Greece . J Glob Antimicrob Resist . 2022 NaN ; 31 : 252 – 255 . 10.1016/j.jgar.2022.10.008 36265802
Rhouma M Madec JY Laxminarayan R. Colistin: From the shadows to a One Health approach for addressing antimicrobial resistance . Int J Antimicrob Agents . 2023 NaN ; 61 ( 2 ): 106713 . 10.1016/j.ijantimicag.2023.106713 36640846
Riley LW. Distinguishing pathovars from nonpathovars: Escherichia coli . Microbiol Spectr . 2020 NaN ; 8 ( 4 2024032 ): 10.1128/microbiolspec.ame-0014-2020
Robin F Beyrouthy R Colot J Saint-Sardos P Berger-Carbonne A Dalmasso G Delmas J Bonnet R. MCR-1 in ESBL-producing Escherichia coli responsible for human infections in New Caledonia . J Antimicrob Chemother . 2017 NaN ; 72 ( 3 ): 946 – 947 . 10.1093/jac/dkw508 28031275
Rodríguez-Santiago J Cornejo-Juárez P Silva-Sánchez J Garza-Ramos U. Polymyxin resistance in Enterobacterales: Overview and epidemiology in the Americas . Int J Antimicrob Agents . 2021 NaN ; 58 ( 5 ): 106426 . 10.1016/j.ijantimicag.2021.106426 34419579
Rodríguez-Santiago J Rodríguez-Medina N Tamayo-Legorreta EM Silva-Sánchez J Téllez-Sosa J Duran-Bedolla J Aguilar-Vera A Lecona-Valera AN Garza-Ramos U Alpuche-Aranda C. Molecular and genomic insights of mcr-1-producing Escherichia coli isolates from piglets . Antibiotics . 2022 NaN ; 11 ( 2 ): 157 . 10.3390/antibiotics11020157 35203760
Shafiq M Rahman SU Bilal H Ullah A Noman SM Zeng M Yuan Y Xie Q Li X Jiao X. Incidence and molecular characterization of ESBL-producing and colistin-resistant Escherichia coli isolates recovered from healthy food-producing animals in Pakistan . J Appl Microbiol . 2022 NaN ; 133 ( 3 ): 1169 – 1182 . 10.1111/jam.15469 35094463
Shen Y Wu Z Wang Y Zhang R Zhou HW Wang S Lei L Li M Cai J Tyrrell J Heterogeneous and flexible transmission of mcr-1 in hospital-associated Escherichia coli . mBio . 2018 NaN ; 9 ( 4 ): e00943 - 18 . 10.1128/mBio.00943-18 29970465
Shi J Zhu H Liu C Xie H Li C Cao X Shen H. Epidemiological and genomic characteristics of global mcr-positive Escherichia coli isolates . Front Microbiol . 2023 NaN ; 13 : 1105401 . 10.3389/fmicb.2022.1105401 36741897
Siguier P Perochon J Lestrade L Mahillon J Chandler M. IS finder: The reference centre for bacterial insertion sequences . Nucleic Acids Res . 2006 NaN ; 34 ( Suppl_1 ): D32 – D36 . 10.1093/nar/gkj014 16381877
Stamatakis A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies . Bioinformatics . 2014 NaN ; 30 ( 9 ): 1312 – 1313 . 10.1093/bioinformatics/btu033 24451623
Szmolka A Gellért Á Szemerits D Rapcsák F Spisák S Adorján A. Emergence and genomic features of a mcr-1 Escherichia coli from duck in Hungary . Antibiotics . 2023 NaN ; 12 ( 10 ): 1519 . 10.3390/antibiotics12101519 37887221
Tang B Chang J Chen Y Lin J Xiao X Xia X Lin J Yang H Zhao G. Escherichia fergusonii, an underrated repository for antimicrobial resistance in food animals . Microbiol Spectr . 2022 NaN ; 10 ( 1 ): e0161721 . 10.1128/spectrum.01617-21 35138151
Tegetmeyer HE Jones SC Langford PR Baltes N. ISApl1, a novel insertion element of Actinobacillus pleuropneumoniae, prevents ApxIV-based serological detection of serotype 7 strain AP76 . Vet Microbiol . 2008 NaN ; 128 ( 3–4 ): 342 – 353 . 10.1016/j.vetmic.2007.10.025 18065168
Tian X Fang R Wu Q Zheng X Zhao Y Dong G Wang C Zhou T Cao J. Emergence of a multidrug-resistant ST 27 Escherichia coli co-harboring blandm-1, mcr-1, and fosA3 from a patient in China . J Antibiot (Tokyo) . 2020 NaN ; 73 ( 9 ): 636 – 641 . 10.1038/s41429-020-0306-5 32341508
Vu Thi Ngoc B Le Viet T Nguyen Thi Tuyet M Nguyen Thi Hong T Nguyen Thi Ngoc D Le Van D Chu Thi L Tran Huy H Penders J Wertheim H Characterization of genetic elements carrying mcr-1 gene in Escherichia coli from the community and hospital settings in Vietnam . Microbiol Spectr . 2022 NaN ; 10 ( 1 ): e0135621 . 10.1128/spectrum.01356-21 35138158
Wang Y Xu C Zhang R Chen Y Shen Y Hu F Liu D Lu J Guo Y Xia X Changes in colistin resistance and mcr-1 abundance in Escherichia coli of animal and human origins following the ban of colistin-positive additives in China: An epidemiological comparative study . Lancet Infect Dis . 2020 NaN ; 20 ( 10 ): 1161 – 1171 . 10.1016/s1473-3099(20)30149-3 32505232
Xie J Liang B Xu X Yang L Li H Li P Qiu S Song H. Identification of mcr-1-positive multidrug-resistant Escherichia coli isolates from clinical samples in Shanghai, China . J Glob Antimicrob Resist . 2022 NaN ; 29 : 88 – 96 . 10.1016/j.jgar.2022.02.008 35182776
Yin Y Qiu L Wang G Guo Z Wang Z Qiu J Li R. Emergence and transmission of plasmid-mediated mobile colistin resistance gene mcr-10 in humans and companion animals . Microbiol Spectr . 2022 NaN ; 10 ( 5 ): e0209722 . 10.1128/spectrum.02097-22 36000890
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 NaN ; 67 ( 11 ): 2640 – 2644 . 10.1093/jac/dks261 22782487
Zhang P Wang J Wang X Bai X Ma J Dang R Xiong Y Fanning S Bai L Yang Z. Characterization of five Escherichia coli isolates coexpressing ESBL and mcr-1 resistance mechanisms from different origins in China . Front Microbiol . 2019 NaN ; 10 : 1994 . 10.3389/fmicb.2019.01994 31555232
Zhou Y Ai W Cao Y Guo Y Wu X Wang B Rao L Xu Y Zhao H Wang X The co-occurrence of NDM-5, mcr-1, and fosa3-encoding plasmids contributed to the generation of extensively drug-resistant Klebsiella pneumoniae . Front Microbiol . 2022 NaN ; 12 : 811263 . 10.3389/fmicb.2021.811263 35046925
