- Short report
- Open Access
Genomic patterns and characterizations of chromosomally-encoded mcr-1 in Escherichia coli populations
Gut Pathogens volume 12, Article number: 55 (2020)
The emergence and transmission of the mobile colistin resistance gene (mcr-1) threatened the extensive use of polymyxin antimicrobials. Accumulated evidence showed that the banning of colistin additive in livestock feed efficiently reduce mcr-1 prevalence, not only in animals but also in humans and environments. However, our previous study has revealed that a small proportion of Escherichia coli could continually carry chromosomally-encoded mcr-1. The chromosomally-encoded events, indicated the existence of stabilized heritage of mcr-1 and revealed a potential threat in the antimicrobial stewardship interventions, are yet to be investigated. In this study, we systematically investigated the genetic basis of chromosomally-encoded mcr-1 in prevalence and potential mechanisms of lineage, plasmid, insertion sequence, and phage. Our results demonstrated that the emergence of chromosomally-encoded mcr-1 could originate from multiple mechanisms, but mainly derived through the recombination of ISApl1/Tn6330. We reported a specific transmission mechanism, which is a phage-like region without lysogenic components, could associate with the emergence and stabilization of chromosomally-encoded mcr-1. These results highlighted the potential origin and risks of chromosomally-encoded mcr-1, which could be a heritable repository and thrive again when confronted with new selective pressures. To the best of our knowledge, this is the first study to systematically reveal the genomic basis of chromosomally-encoded mcr-1, and report a specific transmission pattern involved in phage-like region. Overall, we demonstrate the origin mechanisms and risks of chromosomally-encoded mcr-1. It highlights the need of public attention on chromosome-encoded mcr-1 to prevent from its reemergence.
The emergence and rapid dissemination of plasmid-mediated mobile colistin resistance gene (mcr-1) have become a severe threat to public health . The predominant carriers of mcr-1 were IncX4, IncI2, and IncHI2 plasmids, which are transferable and adaptive plasmid types with broad host range and contributed to the spread of mcr-1 among various sources and bacterial species [2,3,4]. Besides, recombination of transposons, especially Tn6330 (ISApl1-mcr-1-pap2-ISApl1), the primary vehicle for transmission of mcr-1, and phage-like sequences enable mcr-1 to transfer across plasmids and isolates. Such contributed factors facilitated high mcr-1 prevalence in several sources around the world, pushing local governments in Europe, Brazil and China to prohibit the use of colistin as growth promoter additive for livestock [5,6,7,8].
Accumulated evidence showed that banning of colistin in animal feed efficiently restricted mcr-1 prevalence, not only in animals but also in humans and the whole ecosystem in China [2,3,4]. However, our previous study showed that a low proportion of Escherichia coli carrying chromosomally-encoded mcr-1 continually existed in the ecosystem , which was sporadically reported by other studies as well [9,10,11]. On account of the plasmid that could be lost under certain circumstances due to instability, the chromosomally-encoded events could stabilize the heritage of mcr-1, threatening the intervention of colistin stewardship. In current study, we systematically investigate the epidemiological and genomic characterizations of E. coli population with chromosomally-encoded mcr-1.
Based on our previous large-scale epidemiological study from 2016 to 2018 in Guangzhou, China , we identified 24 (3.5%) out of 688 mcr-1-positive E. coli isolates with the chromosomally-encoded mcr-1 (Table 1). The prevalence of chromosomally-encoded mcr-1-positive E. coli was from 0 to 9.8% for each source and from 2.2 to 4.8% for each epoch, indicating that the chromosomally-encoded mcr-1 was at a low prevalence state in different dimensions (Table 1). Additionally, the comparison of prevalence for chromosomally-encoded mcr-1 between different niches or epochs showed no significant difference (Fisher’s exact test, p > 0.05 for each comparison), suggesting that the emergence of chromosomally-encoded mcr-1 was sporadic without temporal or source-specific signals.
To systematically illustrate the genomic basis of chromosomally-encoded mcr-1-positive E. coli population, we collected other 30 E. coli genomes with chromosomally-encoded mcr-1 from published literature for subsequent analysis (Additional file 1: Table S1). Through in silico multilocus sequence typing (MLST) assignment, 32 different sequence types (STs) within 10 ST complexes were determined (Fig. 1). The most common ST among chromosomally-encoded mcr-1-positive E. coli isolates was ST10 (n = 10, 18.5%), which is consistent with the main host for plasmid-mediated mcr-1 on E. coli species [3, 4, 12]. The phylogeny demonstrated two sequence clusters (SCs), except for two isolates which were distinct from two SCs as the outgroup (Fig. 1). The sources and serotypes of these genomes were scattered on the phylogeny, suggesting that the emergence of chromosomally-encoded mcr-1 was random without source- or lineage-based specificity (Fig. 1). Since most of the chromosomally-encoded mcr-1-positive E. coli isolates have been identified in China (n = 40, 74.1%), which was attributed to the extensive screening of mcr-1 in China, the associations between locations and SCs was ambiguous (SC1 [11/16] vs SC2 [29/36], Fisher’s exact test, p = 0.49).
The mcr-1 gene was initially found on plasmids in Enterobacteriaceae and on a transposon Tn6330, prompting that the chromosomally-encoded mcr-1 could come from recombination of plasmid segments or transposition of Tn6330 [13,14,15]. Therefore, we investigated the plasmidome of 54 genomes to illustrate the potential origin of chromosomally-encoded mcr-1. We identified 33 plasmid Inc types among all isolates, and the results showed that the most common Inc type was IncFIB(K) (45.8%, n = 22), followed by IncColRNAI (43.8%, n = 21), IncHI1 (33.3%, n = 16), IncX1 (31.3%, n = 15), IncFIB (AP001918) (27.1%, n = 13), and IncY (20.8%, n = 10). Remarkably, the common Inc types of mcr-1-harboring plasmids, such as IncX4, IncI2, IncHI2, and IncpO111 [1, 3, 4, 12], were rarely detected among these isolates (Fig. 1), indicating that the chromosomally-encoded mcr-1 may derive from ISApl1/Tn6330 through transposition, but not from the plasmid.
We subsequently analyzed the genetic context of mcr-1 for each isolate to investigate the genetic model of chromosomally-encoded mcr-1, except seven isolates were excluded due to short mcr-1-harboring contigs. We found that most of the mcr-1 genes (93.6%, 44/47) were flanked by ISApl1, comprising 24 isolates harboring upstream ISApl1 and 20 isolates carrying composite Tn6330, which complied with the hypothesis of transposition-mediated chromosome insertion.
By mapping the insertion site onto the chromosome of E. coli MG1655, we noted that the distribution of chromosomally-encoded mcr-1 insertion sites was sporadic (Fig. 2a). Thirty-seven clusters of mcr-1-harboring segments were generated based on sequence clustering analysis (Fig. 2a), which included three clusters involving more than one isolates (Fig. 2b) and 34 clusters only containing a single isolate (Additional file 2: Figure S1). The most common genetic pattern of chromosomally-encoded mcr-1 (19.1%, 9/47) involves in an insertion segment in size of ~ 25.7 kb, containing an incomplete phage-like region (score = 40 for phage Vibrio 12B8 [NC_021073] by PHASTER) and a truncated Tn6330 (ISApl1-mcr-1-pap2), which was inserted into the E. coli genome between lysN and hicB (toxin-antitoxin system) loci (Fig. 2b). The incomplete phage-like region only contains head, tail, and fiber protein, and lacks some necessary functional components (Fig. 2b), which seems unfunctional under current conditions. We used BLASTn to search this phage-like sequence in NCBI non-redundant nucleotide database, and the results showed that only five sequences, which are located on E. coli chromosome, were identified with ≥ 60% coverage and ≥ 90% identity, indicating the correlation between chromosomally-encoded mcr-1 and such phage-like region. Collectively, we heuristically concluded that such a phage-like region could mediate the emergence of chromosomally-encoded mcr-1, and then the phage may lose the lysogenic components, stabilization the genetic inheritance of chromosomally-encoded mcr-1. Additionally, the mcr-1 of two isolates showed the insertion of mcr-1 located on an integrative element region and a plasmid segment respectively, suggesting that chromosomally-encoded mcr-1 could be derived from the integration of the integrative region and plasmid segment (Fig. 2c).
In conclusion, our study comprehensively investigated the genetic basis of chromosomally-encoded mcr-1 in prevalence and potential mechanisms of lineage, plasmid, insertion sequence, and phage. Our results showed that chromosomally-encoded mcr-1 was mainly derived from ISApl1 insertion in genomic locations sporadically. Notably, we reported a new transmission mechanism, a phage-like region without functional components, could associate with the emergence and stabilization of chromosomally-encoded mcr-1. The chromosomally-encoded mcr-1 in current situations seems not a severe threat for public health, however, it could be a heritable repository and thrive again if the new selective pressure emerges, because the chromosome-mediated antimicrobial resistance genes (ARGs) might be conferred with genetic sustainability. In-depth investigations are needed to illustrate the genomic and epidemiological dynamics of chromosomally-encoded mcr-1, which may be changed after the approval of colistin in human clinical therapeutics in China .
We searched PubMed using the terms of “mcr-1” [MeSH]/[All Fields] AND “chromosome” [MeSH]/[All Fields] AND “Escherichia coli” [MeSH]/[All Fields] for articles published before 1th October 2020, and identified 20 publications, including 30 available E. coli genomes with chromosome-mediated mcr-1 (Additional file 3: Figure S2).
Antimicrobial resistance genes screening, plasmid incompatibility typing and serotype identification were performed by Center for Genomic Epidemiology (http://www.genomicepidemiology.org/). Multilocus sequence typing (MLST) was assigned using Enterobase (http://enterobase.warwick.ac.uk/). Prophage prediction was implemented by PHASTER . The phylogeny was constructed using RAxML v8.2 with GTR+G model and 1000 bootstrap  based on core genome single-nucleotide polymorphisms (cgSNPs) produced by Roary v3.11.2 and snp-site v2.4.1 . Population structure was assessed using cgSNPs with hierBAPS . The chromosome map was drawn by BRIG v0.95 and marked with insertion pattern manually by Easyfig v2.2.2 [21, 22]. The sequence clustering was performed by CD-HIT-EST .
The significance of prevalence variation of chromosomally-encoded mcr-1 between niches and epochs were tested by Fisher’s exact test using Statistical Package for the Social Sciences (SPSS), version 20.0.
Availability of data and materials
The datasets generated and analysed during the current study are available in the NCBI GenBank repository. The accession number for each genome can be obtained in Additional file 4: Appendix material.
- mcr-1 :
Mobile colistin resistance gene
Multilocus sequence typing
Core genome single-nucleotide polymorphisms
Antimicrobial resistance gene
Liu YY, Wang Y, Walsh TR, Yi LX, Zhang R, Spencer J, Doi Y, Tian GB, Dong BL, Huang XH, Yu LF, Gu DX, Ren HW, Chen XJ, Lv LC, He DD, Zhou HW, Liang ZS, Liu JH, Shen JZ. 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;16:161–8.
Wang Y, Xu C, Zhang R, Chen Y, Shen Y, Hu F, Liu D, Lu J, Guo Y, Xia X, Jiang J, Wang X, Fu Y, Yang L, Wang J, Li J, Cai C, Yin D, Che J, Fan R, Wang Y, Qing Y, Li Y, Liao K, Chen H, Zou M, Liang L, Tang J, Shen Z, Wang S, Yang X, Wu C, Xu S, Walsh TR, Shen J. 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. https://doi.org/10.1016/S1473-3099(20)30149-3.
Jiang Y, Zhang Y, Lu J, Wang Q, Cui Y, Wang Y, Quan J, Zhao D, Du X, Liu H, Li X, Wu X, Hua X, Feng Y, Yu Y. Clinical relevance and plasmid dynamics of mcr-1-positive Escherichia coli in China: a multicentre case-control and molecular epidemiological study. Lancet Microbe. 2020;1:e24–33.
Shen C, Zhong L-L, Yang Y, Doi Y, Paterson DL, Stoesser N, Ma F, El-Sayed Ahmed MAE-G, Feng S, Huang S, Li H-Y, Huang X, Wen X, Zhao Z, Lin M, Chen G, Liang W, Liang Y, Xia Y, Dai M, Chen D-Q, Zhang L, Liao K, Tian GB. 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. Lancet Microbe. 2020;1:e34–43.
Walsh TR, Wu Y. China bans colistin as a feed additive for animals. Lancet Infect Dis. 2016;16:1102–3.
Miguela-Villoldo P, Hernandez M, Moreno MA, Rodriguez-Lazaro D, Quesada A, Dominguez L, Ugarte-Ruiz M. National colistin sales versus colistin resistance in Spanish pig production. Res Vet Sci. 2019;123:141–3.
Monte DF, Mem A, Fernandes MR, Cerdeira L, Esposito F, Galvao JA, Franco B, Lincopan N, Landgraf M. Chicken meat as a reservoir of colistin-resistant Escherichia coli strains carrying mcr-1 genes in South America. Antimicrob Agents Chemother. 2017;61:e02718.
Duggett NA, Randall LP, Horton RA, Lemma F, Kirchner M, Nunez-Garcia J, Brena C, Williamson SM, Teale C, Anjum MF. Molecular epidemiology of isolates with multiple mcr plasmids from a pig farm in Great Britain: the effects of colistin withdrawal in the short and long term. J Antimicrob Chemother. 2018;73:3025–33.
Li R, Yu H, Xie M, Chen K, Dong N, Lin D, Chan EW, Chen S. Genetic basis of chromosomally-encoded mcr-1 gene. Int J Antimicrob Agents. 2018;51:578–85.
Lu XY, Xiao X, Liu Y, Li Y, Li RC, Wang ZQ. Chromosome-mediated mcr-1 in Escherichia coli strain L73 from a goose. Int J Antimicrob Agents. 2019;54:99–101.
Peng Z, Hu ZZ, Li ZG, Li XS, Jia CY, Zhang XX, Wu B, Chen HC, Wang XR. Characteristics of a colistin-resistant Escherichia coli ST695 harboring the chromosomally-encoded mcr-1 gene. Microorganisms. 2019;7:558.
Wang S, Shen J. Active surveillance of the spread of mcr-1-positive E coli. The Lancet Microbe. 2020;1:e4–5.
Snesrud E, He S, Chandler M, Dekker JP, Hickman AB, McGann P, Dyda F. A model for transposition of the colistin resistance gene mcr-1 by ISApl1. Antimicrob Agents Chemother. 2016;60:6973–6.
Li R, Chen K, Chan EW, Chen S. Characterization of the stability and dynamics of Tn6330 in an Escherichia coli strain by nanopore long reads. J Antimicrob Chemother. 2019;74:1807–11.
Snesrud E, McGann P, Chandler M. The birth and demise of the ISApl1-mcr-1-ISApl1 composite transposon: the vehicle for transferable colistin resistance. Mbio. 2018;9:e02381.
Huang H, Dong N, Shu L, Lu J, Sun Q, Chan EW, Chen S, Zhang R. Colistin-resistance gene mcr in clinical carbapenem-resistant Enterobacteriaceae strains in China, 2014-2019. Emerg Microbes Infect. 2020;9:237–45.
Arndt D, Marcu A, Liang Y, Wishart DS. PHAST, PHASTER and PHASTEST: tools for finding prophage in bacterial genomes. Brief Bioinform. 2019;20:1560–7.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.
Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, Fookes M, Falush D, Keane JA, Parkhill J. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31:3691–3.
Tonkin-Hill G, Lees JA, Bentley SD, Frost SDW, Corander J. RhierBAPS: an R implementation of the population clustering algorithm hierBAPS. Wellcome Open Res. 2018;3:93.
Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27:1009–10.
Alikhan NF, Petty NK, Ben Zakour NL, Beatson SA. BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011;12:402.
Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28:3150–2.
We acknowledge Ms. Lujie Liang (Department of Microbiology, Zhongshan School of Medicine, Sun Yat-sen University) for assistance in writing of the manuscript.
This work was supported by the National Natural Science Foundation of China (Grant Numbers 81830103, 81722030), China Postdoctoral Science Foundation (BX20200394), National Key Research and Development Program (Grant Number 2017ZX10302301), Guangdong Natural Science Foundation (Grant Number 2017A030306012), Project of high-level health teams of Zhuhai at 2018 (The Innovation Team for Antimicrobial Resistance and Clinical Infection), 111 Project (Grant Number B12003), Open project of Key Laboratory of Tropical Disease Control (Sun Yat-sen University), Ministry of Education (Grant Number 2020kfkt04, 2020kfkt07), Innovative Development Program for Outstanding Graduate Students (No. 19ykyj557), and The Science and Technology Planning Project of Guangdong (2017A020215017).
Ethics approval and consent to participate
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.
The genetic structure of chromosomally-encoded mcr-1 patterns which included only one isolate. The number for each pattern was identical to Fig. 2a.
Flow diagram of the study selection process.
About this article
Cite this article
Shen, C., Zhong, LL., Ma, F. et al. Genomic patterns and characterizations of chromosomally-encoded mcr-1 in Escherichia coli populations. Gut Pathog 12, 55 (2020). https://doi.org/10.1186/s13099-020-00393-2
- Antimicrobial resistance
- Genomic pattern
- Insertion sequence