Dynamics of mcr-1 prevalence and mcr-1-positive Escherichia coli after the cessation of colistin use as a feed additive for animals in China: a prospective cross-sectional and whole genome sequencing based molecular epidemiological study

Background The global dissemination of colistin resistance encoded by mcr-1 has been attributed to extensive use of colistin in livestock, threatening colistin efficacy in medicine. The emergence of mcr-1 in common pathogens, such as Escherichia coli, is of particular concern. Therefore, China banned the use of colistin in animal feed from May 1ST 2017. We investigated subsequent changes in mcr-1 prevalence, and the genomic epidemiology of mcr-1-positive Escherichia coli (MCRPEC). Methods Sampling was conducted pre- (October-December 2016) and post-colistin ban (October-December, 2017 and 2018, respectively). 3675 non-duplicate pig fecal samples were collected from 14 provinces (66 farms) in China to determine intervention-related changes in mcr-1 prevalence. 15193 samples were collected from pigs, healthy human volunteers, colonized and infected hospital inpatients, food and the environment in Guangzhou, to characterize source-specific mcr-1 prevalence and the wider ecological impact of the ban. From these samples, 688 MCRPEC were analyzed with whole genome sequencing (WGS), plasmid conjugation and S1-PFGE/Southern blots to characterize associated genomic changes. Findings After the ban, mcr-1 prevalence decreased significantly in national pig farms, from 45·0% (308/684 samples) in 2016, to 19·4% (274/1416) in 2018 (p<0·0001). This trend was mirrored in samples from most sources in Guangzhou (overall 19·2% [959/5003 samples] in 2016; 5·3% [238/4489] in 2018; p<0·0001). The population structure of MCRPEC was diverse (23 sequence clusters [SCs]); ST10 clonal complex isolates were predominant (247/688 [36%]). MCRPEC causing infection in hospitalized inpatients were genetically more distinct and appeared less affected by the ban. mcr-1 was predominantly found on plasmids (632/688 [92%]). Common mcr-1 plasmid types included IncX4, IncI2 and IncHI2 (502/656 [76.5%]); significant increases in IncI2-associated mcr-1 and a distinct lineage of mcr-1-associated IncHI2 were observed post-ban. Changes in the frequency of mcr-1-associated flanking sequences (ISApl1-negative MCRPEC), 63 core genome SNPs and 30 accessory genes were also significantly different after the ban, consistent with rapid genetic adaptation in response to changing selection pressures. Interpretation A rapid, ecosystem-wide, decline in mcr-1 was observed after banning the use of colistin in animal feed, with associated genetic changes in MCRPEC. Genomic surveillance is key to assessing and monitoring stewardship interventions. Funding National Natural Science Foundation of China

[1]  Huanchun Chen,et al.  Characteristics of a Colistin-Resistant Escherichia coli ST695 Harboring the Chromosomally-Encoded mcr-1 Gene , 2019, Microorganisms.

[2]  T. V. Van Boeckel,et al.  Global trends in antimicrobial resistance in animals in low- and middle-income countries , 2019, Science.

[3]  T. Dagan,et al.  The Effect of Population Bottleneck Size and Selective Regime on Genetic Diversity and Evolvability in Bacteria , 2019, bioRxiv.

[4]  Yuan Liu,et al.  Chromosome-mediated mcr-1 in Escherichia coli strain L73 from a goose. , 2019, International Journal of Antimicrobial Agents.

[5]  J. Li,et al.  The rise and spread of mcr plasmid-mediated polymyxin resistance , 2019, Critical reviews in microbiology.

[6]  Yan Liu,et al.  Transmission of mcr-1-Producing Multidrug-resistant Enterobacteriaceae in Public Transportation in Guangzhou, China. , 2018, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[7]  Davide Heller,et al.  eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses , 2018, Nucleic Acids Res..

[8]  Jianzhong Shen,et al.  Anthropogenic and environmental factors associated with high prevalence of mcr-1 carriage in humans across China , 2018, Nature Microbiology.

[9]  S. Sheppard,et al.  Population genomics of bacterial host adaptation , 2018, Nature Reviews Genetics.

[10]  L. Leibovici,et al.  Colistin alone versus colistin plus meropenem for treatment of severe infections caused by carbapenem-resistant Gram-negative bacteria: an open-label, randomised controlled trial. , 2018, The Lancet. Infectious diseases.

[11]  E. Snesrud,et al.  The Birth and Demise of the ISApl1-mcr-1-ISApl1 Composite Transposon: the Vehicle for Transferable Colistin Resistance , 2018, mBio.

[12]  A. Walker,et al.  High Rates of Human Fecal Carriage of mcr-1–Positive Multidrug-Resistant Enterobacteriaceae Emerge in China in Association With Successful Plasmid Families , 2018, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[13]  F. Codjoe,et al.  Carbapenem Resistance: A Review , 2017, Medical sciences.

[14]  Z. Iqbal,et al.  The global distribution and spread of the mobilized colistin resistance gene mcr-1 , 2017, bioRxiv.

[15]  Heather Ganshorn,et al.  Restricting the use of antibiotics in food-producing animals and its associations with antibiotic resistance in food-producing animals and human beings: a systematic review and meta-analysis , 2017, The Lancet. Planetary health.

[16]  Miaomiao Xie,et al.  Widespread distribution of mcr-1-bearing bacteria in the ecosystem, 2015 to 2016 , 2017, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[17]  Jian Sun,et al.  Expanding landscapes of the diversified mcr-1-bearing plasmid reservoirs , 2017, Microbiome.

[18]  A. Letellier,et al.  Extended-spectrum β-lactamases, carbapenemases and the mcr-1 gene: is there a historical link? , 2017, International journal of antimicrobial agents.

[19]  Lonneke Scheffer,et al.  Erratum to: Rapid scoring of genes in microbial pan-genome-wide association studies with Scoary , 2016, Genome Biology.

[20]  Lonneke Scheffer,et al.  Rapid scoring of genes in microbial pan-genome-wide association studies with Scoary , 2016, Genome Biology.

[21]  Yongning Wu,et al.  China bans colistin as a feed additive for animals. , 2016, The Lancet. Infectious diseases.

[22]  E. Snesrud,et al.  A Model for Transposition of the Colistin Resistance Gene mcr-1 by ISApl1 , 2016, Antimicrobial Agents and Chemotherapy.

[23]  Jian-Hua Liu,et al.  Carbapenem-resistant and colistin-resistant Escherichia coli co-producing NDM-9 and MCR-1. , 2016, The Lancet. Infectious diseases.

[24]  Jianzhong Shen,et al.  Early emergence of mcr-1 in Escherichia coli from food-producing animals. , 2016, The Lancet. Infectious diseases.

[25]  Simon R. Harris,et al.  SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments , 2016, bioRxiv.

[26]  Daniel J. Wilson,et al.  Within-host evolution of bacterial pathogens , 2016, Nature Reviews Microbiology.

[27]  Jianzhong Shen,et al.  Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. , 2015, The Lancet. Infectious diseases.

[28]  Andrew J. Page,et al.  Roary: rapid large-scale prokaryote pan genome analysis , 2015, bioRxiv.

[29]  Justin Zobel,et al.  SRST2: Rapid genomic surveillance for public health and hospital microbiology labs , 2014, bioRxiv.

[30]  Torsten Seemann,et al.  Prokka: rapid prokaryotic genome annotation , 2014, Bioinform..

[31]  Sergey I. Nikolenko,et al.  SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing , 2012, J. Comput. Biol..

[32]  Thomas L. Madden,et al.  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. , 1997, Nucleic acids research.

[33]  B. Yangco,et al.  CDC definitions for nosocomial infections. , 1989, American journal of infection control.

[34]  J M Hughes,et al.  CDC definitions for nosocomial infections, 1988. , 1988, American journal of infection control.