Genomic Insights into CRISPR-Harboring Plasmids in the Klebsiella Genus: Distribution, Backbone Structures, Antibiotic Resistance, and Virulence Determinant Profiles

CRISPR systems are often encoded by many prokaryotes as adaptive defense against mobile genetic elements (MGEs), but several MGEs also recruit CRISPR components to perform additional biological functions. Type IV-A systems are identified in Klebsiella plasmids, yet the distribution, characterization, and role of these plasmids carrying CRISPR systems in the whole Klebsiella genus remain unclear. Here, we performed large-scale comparative analysis of these plasmids using publicly available plasmid genomes. CRISPR-harboring plasmids were mainly distributed in Klebsiella pneumoniae (9.09%), covering 19.23% of sequence types, but sparse in Klebsiella species outside Klebsiella pneumoniae (3.92%). Plasmid genome comparison reiterated that these plasmids often carried the cointegrates of IncFIB and IncHI1B replicons, occasionally linked to other replicons, such as IncFIA, IncFII, IncR, IncQ, and IncU. Comparative genome analysis showed that CRISPR-carrying Klebsiella plasmids shared a conserved pNDM-MAR-like conjugation module as their backbones and served as an important vector for the accretion of antibiotic resistance genes (ARGs) and even virulence genes (VGs). Moreover, compared with CRISPR-negative IncFIB/IncHIB plasmids, CRISPR-positive IncFIB/IncHIB plasmids displayed high divergences in terms of ARGs, VGs, GC content, plasmid length, and backbone structures, suggesting their divergent evolutionary paths. The network analysis revealed that CRISPR-positive plasmids yielded fierce competitions with other plasmid types, especially conjugative plasmids, thereby affecting the dynamics of plasmid transmission. Overall, our study provides valuable insights into the role of CRISPR-positive plasmids in the spread of ARGs and VGs in Klebsiella genus.