The correlation between CRISPR-Cas systems and plasmid-mediated bacterial antibiotic resistance is increasingly growing attention. However, currently no reports exist on the relationship between the CRISPR-Cas systems and the carriage of blaNDM or plasmids in E. coli. Here, molecular characterization and phylogenetic analysis of 639 E. coli isolated from humans in China were carried out. Depending on similarity in sequence, the type I-E CRISPR-Cas systems in E. coli can be grouped into two distinct clades, which we refer to for descriptive purposes within this study as the type I-E-S1 and I-E-S2, whereas the type I-E-S2 CRISPR-Cas system is further divided into I-E-S2a and I-E-S2b systems based on the presence of cas8e and cas11. ST167 (phylogroup A) and ST410 (phylogroup C) E. coli were observed bearing the type I-E-S1 and I-E-S2b systems, respectively. Compared with strains carrying the I-E-S1 type CRISPR-Cas system, the blaNDM carrying rate, the positive rate of IncX3 plasmid, and the positive rate of IncF plasmid of strains with the I-E-S2a type CRISPR-Cas system were evidently lower (p < 0.05); the blaNDM carrying rate and the positive rate of IncF plasmid of strains with the I-E-S2b type CRISPR-Cas system were evidently higher (p < 0.05). The blaNDM positive rate and IncF plasmid positive rate of strains carrying the I-E-S2a type CRISPR-Cas system were significantly lower than those of strains carrying the I-E-S2b type CRISPR-Cas system (p < 0.001). It proves that the I-E-S1, I-E-S2a, and I-E-S2b type CRISPR-Cas systems are beneficial for spreading blaNDM and IncX3 plasmids. We found significant differences in the cas gene sequences of the I-E-S1 and I-E-S2 type CRISPR loci. The type I-E CRISPR-Cas systems in E. coli isolated from Chinese sources are classified further for the first time, revealing their high correlation with blaNDM, phylogenetic groups, and multilocus sequence typing. This work paves the way for a deeper understanding of the role that CRISPR-Cas systems play in the rise of resistant E. coli ST167 and ST410.
Get full access to this article
View all access options for this article.
References
1.
MiottoM, OssaiSA, MeredithJE, et al.Genotypic and phenotypic characterization of Escherichia coli isolated from mollusks in Brazil and the United States. Microbiologyopen, 2019; 8(5):e00738; doi: 10.1002/mbo3.738
2.
FritzenwankerM, ImirzaliogluC, HeroldS, et al.Treatment options for carbapenem- resistant gram-negative infections. Dtsch Arztebl Int, 2018; 115(20–21):345–352; doi: 10.3238/arztebl.2018.0345
3.
ClermontO, ChristensonJK, DenamurE, et al.The Clermont Escherichia coli phylo-typing method revisited: Improvement of specificity and detection of new phylo-groups. Environ Microbiol Rep, 2013; 5(1):58–65; doi: 10.1111/1758-2229.12019
4.
TouchonM, CharpentierS, PognardD, et al.Antibiotic resistance plasmids spread among natural isolates of Escherichia coli in spite of CRISPR elements. Microbiology (Reading), 2012; 158(12):2997–3004; doi: 10.1099/mic.0.060814-0
5.
AydinS, PersonneY, NewireE, et al.Presence of Type I-F CRISPR/Cas systems is associated with antimicrobial susceptibility in Escherichia coli. J Antimicrob Chemother, 2017; 72(8):2213–2218; doi: 10.1093/jac/dkx137
6.
Al-FarsiHM, CamporealeA, IninbergsK, et al.Clinical and molecular characteristics of carbapenem non-susceptible Escherichia coli: A nationwide survey from Oman. PLoS One, 2020; 15(10):e0239924; doi: 10.1371/journal.pone.0239924
7.
WeissmanSJ, HansenNI, Zaterka-BaxterK, et al.Emergence of antibiotic resistance-associated clones among Escherichia coli recovered from newborns with early-onset sepsis and meningitis in the United States, 2008-2009. J Pediatric Infect Dis Soc, 2016; 5(3):269–276; doi: 10.1093/jpids/piv013
8.
MaidenMC. Multilocus sequence typing of bacteria. Annu Rev Microbiol, 2006; 60:561–588; doi: 10.1146/annurev.micro.59.030804.121325
9.
FengY, LiuL, LinJ, et al.Key evolutionary events in the emergence of a globally disseminated, carbapenem resistant clone in the Escherichia coli ST410 lineage. Commun Biol, 2019; 2:322; doi: 10.1038/s42003-019-0569-1
10.
WuW, FengY, TangG, et al.NDM Metallo-β-lactamases and their bacterial producers in health care settings. Clin Microbiol Rev, 2019; 32(2):e00115–e00118; doi: 10.1128/CMR.00115-18
11.
ZhangR, LiuL, ZhouH, et al.Nationwide Surveillance of Clinical Carbapenem-resistant Enterobacteriaceae (CRE) Strains in China. EBioMedicine, 2017; 19:98–106; doi: 10.1016/j.ebiom.2017.04.032
12.
ZongZ, FennS, ConnorC, et al.Complete genomic characterization of two Escherichia coli lineages responsible for a cluster of carbapenem-resistant infections in a Chinese hospital. J Antimicrob Chemother, 2018; 73(9):2340–2346; doi: 10.1093/jac/dky210
13.
YongD, TolemanMA, GiskeCG, et al.Characterization of a new metallo-beta-lactamase gene, bla(NDM-1), and a novel erythromycin esterase gene carried on a unique genetic structure in Klebsiella pneumoniae sequence type 14 from India. Antimicrob Agents Chemother, 2009; 53(12):5046–5054; doi: 10.1128/AAC.00774-09
14.
IslamMA, IslamM, HasanR, et al.Environmental Spread of New Delhi Metallo-β-Lactamase-1-Producing Multidrug-Resistant Bacteria in Dhaka, Bangladesh. Appl Environ Microbiol, 2017; 83(15):e00793–e00717; doi: 10.1128/AEM.00793-17
15.
DingY, WangY, HsiaY, et al.Systematic review of carbapenem-resistant Enterobacteriaceae causing neonatal sepsis in China. Ann Clin Microbiol Antimicrob, 2019; 18(1):36; doi: 10.1186/s12941-019-0334-9
16.
SugawaraY, AkedaY, SakamotoN, et al.Genetic characterization of blaNDM-harboring plasmids in carbapenem-resistant Escherichia coli from Myanmar.PLoS One, 2017; 12(9):e0184720; doi: 10.1371/journal.pone.0184720
17.
GarneauJE, DupuisM-È, VillionM, et al.The CRISPR/Cas bacterial immune system cleaves bacteriophage and plasmid DNA. Nature, 2010; 468(7320):67–71; doi: 10.1038/nature09523
18.
BarrangouR, FremauxC, DeveauH, et al.CRISPR provides acquired resistance against viruses in prokaryotes. Science, 2007; 315(5819):1709–1712; doi: 10.1126/science.1138140
19.
JinekM, ChylinskiK, FonfaraI, et al.A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science, 2012; 337(6096):816–821; doi: 10.1126/science.1225829
20.
LinTL, PanYJ, HsiehPF, et al.Imipenem represses CRISPR-Cas interference of DNA acquisition through H-NS stimulation in Klebsiella pneumoniae. Sci Rep, 2016; 6:31644; doi: 10.1038/srep31644
21.
ZhouY, TangY, FuP, et al.The type I-E CRISPR-Cas system influences the acquisition of blaKPC-IncF plasmid in Klebsiella pneumonia. Emerg Microbes Infect, 2020; 9(1):1011–1022; doi: 10.1080/22221751.2020.1763209
22.
TangY, FuP, ZhouY, et al.Absence of the type I-E CRISPR-Cas system in Klebsiella pneumoniae clonal complex 258 is associated with dissemination of IncF epidemic resistance plasmids in this clonal complex. J Antimicrob Chemother, 2020; 75(4):890–895; doi: 10.1093/jac/dkz538
AhatorSD, SagarS, ZhuM, et al.Nutrient availability and phage exposure alter the Quorum-Sensing and CRISPR-Cas-controlled population dynamics of Pseudomonas aeruginosa. mSystems, 2022; 7(4):e0009222; doi: 10.1128/msystems.00092-22
26.
MakarovaKS, ShmakovSA, WolfYI, et al.An updated evolutionary classification of CRISPR-Cas systems including rare variants. Nat Microbiol, 2025; 10(12):3346–3361; doi: 10.1038/s41564-025-02180-8
27.
JacksonSA, McKenzieRE, FagerlundRD, et al.CRISPR-Cas: Adapting to change. Science, 2017; 356(6333):eaal5056; doi: 10.1126/science.aal5056
28.
RoerL, Overballe-PetersenS, HansenF, et al.Escherichia coli sequence type 410 is causing new international high-risk clones. mSphere, 2018; 3(4):e00337–e00318; doi: 10.1128/mSphere.00337-18
29.
AlmendrosC, GuzmánNM, Díez-VillaseñorC, et al.Target motifs affecting natural immunity by a constitutive CRISPR-Cas system in Escherichia coli. PLoS One, 2012; 7(11):e50797; doi: 10.1371/journal.pone.0050797
30.
Díez-VillaseñorC, AlmendrosC, García-MartínezJ, et al.Diversity of CRISPR loci in Escherichia coli. Microbiology (Reading), 2010; 156(5):1351–1361; doi: 10.1099/mic.0.036046-0
31.
LongJ, XuY, OuL, et al.Polymorphism of Type I-F CRISPR/Cas system in Escherichia coli of phylogenetic group B2 and its application in genotyping. Infect Genet Evol., 2019; 74:103916; doi: 10.1016/j.meegid.2019.103916
32.
TrinettaV, MagossiG, AllardMW, et al.Characterization of Salmonella enterica Isolates from Selected U.S. Swine Feed Mills by Whole-Genome Sequencing. Foodborne Pathog Dis, 2020; 17(2):126–136; doi: 10.1089/fpd.2019.2701
FeldgardenM, BroverV, HaftDH, et al.Erratum for Feldgarden et al., “Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates”. Antimicrob Agents Chemother, 2020; 64(4):e00361-20; doi: 10.1128/AAC.00361-20
36.
BeghainJ, Bridier-NahmiasA, Le NagardH, et al.ClermonTyping: An easy-to-use and accurate in silico method for Escherichia genus strain phylotyping. Microb Genom, 2018; 4(7):e000192; doi: 10.1099/mgen.0.000192
37.
GrissaI, VergnaudG, PourcelC. CRISPRFinder: A web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res, 2007; 35(Web Server issue)issue:W52–W57; doi: 10.1093/nar/gkm360
38.
AbbySS, NéronB, MénagerH, et al.MacSyFinder: A program to mine genomes for molecular systems with an application to CRISPR-Cas systems. PLoS One, 2014; 9(10):e110726; doi: 10.1371/journal.pone.0110726
39.
CouvinD, BernheimA, Toffano-NiocheC, et al.CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids Res, 2018; 46(W1):W246–W251; doi: 10.1093/nar/gky425
40.
OlawoyeIB, FrostSDW, HappiCT. The Bacteria Genome Pipeline (BAGEP): An automated, scalable workflow for bacteria genomes with Snakemake. PeerJ, 2020; 8:e10121; doi: 10.7717/peerj.10121
41.
CroucherNJ, PageAJ, ConnorTR, et al.Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res, 2015; 43(3):e15; doi: 10.1093/nar/gku1196
42.
PriceMN, DehalPS, ArkinAP. FastTree: Computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol., 2009; 26(7):1641–1650; doi: 10.1093/molbev/msp077
43.
PriceMN, DehalPS, ArkinAP. FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One, 2010; 5(3):e9490; doi: 10.1371/journal.pone.0009490
44.
LetunicI, BorkP. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res, 2016; 44(W1):W242–W245; doi: 10.1093/nar/gkw290
45.
EdgarRC. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res, 2004; 32(5):1792–1797; doi: 10.1093/nar/gkh340
46.
BiswasA, GagnonJN, BrounsSJ, et al.CRISPRTarget: Bioinformatic prediction and analysis of crRNA targets. RNA Biol, 2013; 10(5):817–827; doi: 10.4161/rna.24046
47.
TouchonM, CharpentierS, ClermontO, et al.CRISPR distribution within the Escherichia coli species is not suggestive of immunity-associated diversifying selection. J Bacteriol, 2011; 193(10):2460–2467; doi: 10.1128/JB.01307-10
48.
ColpanA, JohnstonB, PorterS, et al.; VICTORY (Veterans Influence of Clonal Types on Resistance: Year 2011) Investigators. Escherichia coli sequence type 131 (ST131) subclone H30 as an emergent multidrug-resistant pathogen among US veterans. Clin Infect Dis, 2013; 57(9):1256–1265; doi: 10.1093/cid/cit503
49.
ShenJ, LvL, WangX, et al.Comparative analysis of CRISPR-Cas systems in Klebsiella genomes. J Basic Microbiol, 2017; 57(4):325–336; doi: 10.1002/jobm.201600589
50.
LiHY, KaoCY, LinWH, et al.Characterization of CRISPR-Cas systems in clinical Klebsiella pneumoniae isolates uncovers its potential association with antibiotic susceptibility. Front Microbiol, 2018; 9:1595; doi: 10.3389/fmicb.2018.01595
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
