Type 1, 2, and 1/2-Hybrid IncC Plasmids From China.

A collection of 11 IncC plasmids from China were fully sequenced herein and compared with reference plasmids pR148 and pR55. These 13 plasmids could be assigned into three different subgroups: type 1, type 2, and type 1/2 hybrid. Type 1/2-hybrid plasmids most likely emerged from homologous recombination between type 1 and type 2 plasmids. Different IncC plasmids had evolved to acquire quite different profiles of accessory modules and thus different collections of resistance genes. The accessory resistance modules included not only the bla CMY-carrying region, the ARI-A island, and the ARI-B island, but also various additional kinds of resistance islands such as the bla CTX-M-carrying regions and the MDR regions. Insertion of accessory modules was sometimes accompanied by deletion, inversion, and translocation of surrounding backbone regions. pR148 and pR55 were confirmed to have the most complete backbones for type 1 and type 2, respectively. This was the first report of a bla IMP- 8-carrying IncC plasmid, and that of three novel mobile elements: a Tn1696-derived unit transposon Tn6395, a class 2 integron In2-76, and an insertion sequence ISEcl10.

Introduction

Plasmids of IncA and IncC incompatibility groups have been combined into IncA/C since 1970s, but IncA and IncC groups are essentially compatible (Ambrose et al., 2018) and have significant nucleotide divergence across their backbones (Harmer and Hall, 2015) and, thus, IncA [reference plasmid pRA1 (Fricke et al., 2009)] and IncC should be separated from IncA/C (Harmer et al., 2017). IncC plasmids have two distinct types, namely, type 1 and type 2 (Harmer and Hall, 2014), which are represented by reference plasmids pR148 (Del Castillo et al., 2013), and pR55 (Benoît et al., 2012), respectively. These two types can be distinguished by two genes of substitution (orf1832 in type 1 or orf1847 in type 2, and rhs1 in type 1 or rhs2 in type 2) and two small intergenic insertions (i1 and i2 found in type 2 rather than type 1) (Harmer and Hall, 2014).

A large array of antibiotic resistance genes have been found in IncC plasmids and are clustered in various accessory modules inserted at various sites of the IncC backbones (Harmer and Hall, 2015). IncC plasmids contain at least three hotspots for integration of the blaCMY-carrying region and two antibiotic resistance islands designated ARI-A and ARI-B. The blaCMY-carrying region and ARI-A are frequently found in type 1 plasmids but not type 2; both type 1 and type 2 plasmids carry ARI-B, but this island is not always present in these plasmids (Harmer and Hall, 2015).

Our previous study (Ma et al., 2017) has shown considerable genomic variation in a collection of three type 1 IncC plasmids from China together with pR148. This follow-up study presents the complete nucleotide sequences of 11 new IncC plasmids from China and a further comprehensive genomic comparison of these 11 plasmids together with pR148 and pR55, providing a deeper insight into parallel evolution of IncC plasmids from China.

Materials and Methods

Bacterial Strains

The 11 isolates (Supplementary Table S1) were recovered from the sputum or blood specimens of 11 different patients with nosocomial infections in 10 Chinese public hospitals from 2012 to 2016. Bacterial species identification was performed using genome sequence-based average nucleotide identity analysis (Richter and Rossello-Mora, 2009).

Conjugal Transfer

Plasmid conjugal transfer experiments were carried out with rifampin-resistant Escherichia coli EC600 or sodium azide-resistant E. coli J53 being used as a recipient and each of the QD1501, A2293, and 205880 isolates as a donor. Three milliliters of overnight cultures of each of donor and recipient bacteria were mixed together, harvested, and resuspended in 80 μl of Brain Heart Infusion (BHI) broth (BD Biosciences). The mixture was spotted on a 1 cm2 hydrophilic nylon membrane filter with a 0.45 μm pore size (Millipore) that was placed on a BHI agar (BD Biosciences) plate and then incubated for mating at 37°C for 12 to 18 h. Bacteria were washed from filter membrane and spotted on Muller-Hinton (MH) agar (BD Biosciences) plates, for selecting an E. coli transconjugant carrying blaNDM (pQD1501-Ct1), blaCTX–M (pA2293-Ct2), or blaCMY (p205880-Ct1/2). 200 μg/ml sodium azide (for J53) or 1000 μg/ml rifampin (for EC600), together with 4 μg/ml imipenem (for blaNDM) or 200 μg/ml ampicillin (for blaCTX–M or blaCMY), was for E. coli transconjugant selection.

Sequencing and Sequence Assembly

Genomic DNA was isolated from each of the above 11 isolates using an UltraClean Microbial Kit (Qiagen, NW, Germany). The genomic DNA of strains 427113, T5282, and 397108 was sequenced from mate-pair libraries with an average insert size of 5 kb (ranging from 2 to 10 kb) using a MiSeq sequencer (Illumina, CA, United States) and DNA contigs were assembled based on their contig coverages using Newbler 2.6 (Nederbragt, 2014). Quality control, removing adapters, and low quality reads were performed using Trimmomatic 0.36 (Bolger et al., 2014). The filtered clean reads were then assembled using Newbler 2.6 (Nederbragt, 2014), followed by extraction of the consensus sequence with CLC Genomics Workbench 3.0 (Qiagen Bioinformatics). Gapfiller V1.11 (Boetzer and Pirovano, 2012) was used for gap closure. For all the other eight isolates, genome sequencing was performed with a sheared DNA library with an average size of 15 kb (ranging from 10 to 20 kb) on a PacBio RSII sequencer (Pacific Biosciences, CA, United States), as well as a paired-end library with an average insert size of 400 bp (ranging from 150 to 600 kb) on a HiSeq sequencer (Illumina, CA, United States). The paired-end short Illumina reads were used to correct the long PacBio reads utilizing proovread (Hackl et al., 2014), and then the corrected PacBio reads were assembled de novo utilizing SMARTdenovo[1].

Sequence Annotation and Comparison

Open reading frames (ORFs) and pseudogenes were predicted using RAST 2.0 (Brettin et al., 2015) combined with BLASTP/BLASTN (Boratyn et al., 2013) searches against the UniProtKB/Swiss-Prot database (Boutet et al., 2016) and the RefSeq database (O’Leary et al., 2016). Annotation of resistance genes, mobile elements, and other features was carried out using the online databases including CARD (Jia et al., 2017), ResFinder (Zankari et al., 2012), ISfinder (Siguier et al., 2006), INTEGRALL (Moura et al., 2009), and Tn Number Registry (Roberts et al., 2008). Multiple and pairwise sequence comparisons were performed using MUSCLE 3.8.31 (Edgar, 2004) and BLASTN, respectively. Gene organization diagrams were drawn in Inkscape 0.48.1[2]. Heatmaps were plotted with MeV 4.9.0 (Saeed et al., 2003).

Phylogenetic Analysis

The backbone regions of indicative plasmids were aligned using MUMmer 3.0 (Kurtz et al., 2004). Inference of homologous recombination was performed using ClonalFrameML (Didelot and Wilson, 2015) to remove recombination-associated single-nucleotide polymorphisms (SNPs). A maximum-likelihood tree was constructed from recombination-free SNPs using MEGA7 (Kumar et al., 2016) with a bootstrap iteration of 1000.

Antimicrobial Susceptibility Testing

Bacterial antimicrobial susceptibility was tested by BioMérieux VITEK 2 and interpreted as per the 2017 Clinical and Laboratory Standards Institute (CLSI) guidelines (CLSI, 2017).

Nucleotide Sequence Accession Numbers

The complete nucleotide sequences of plasmids p24845-Ct2, p205880-Ct1/2, pT5282-Ct2, pKpn47-Ct1/2, p11935-Ct1/2, pQD1501-Ct1, p12085-Ct1, pA2293-Ct2, p397108-Ct2, p427113-Ct1/2, and pA1763-Ct2 were submitted to GenBank under the accession numbers MF344572 to MF344574, MN310369, MN310375, MN310377, MN310378, MH917284, MH917285, and MG764552, respectively.

Results

Overview of Sequenced IncC Plasmids

Of the 11 IncC plasmids fully sequenced in this work (Table 1 and Figure 1), two and five could be assigned into type 1 and type 2, respectively, based on the presence or absence of orf1832/orf1847, rhs1/rhs2, i1, and i2, while the remaining four, containing orf1832 and rhs1 (characteristic of type 1) plus i2 (characteristic of type 2), were recognized as type 1/2 hybrid. A total of 597 core SNPs (among them 69 were recombination-free) were identified from the backbone regions of these 11 plasmids together with type 1 reference pR148 and type 2 reference pR55. As shown in the phylogenetic tree (Figure 1A) constructed from the 69 recombination-free SNPs, these 13 plasmids could be clustered into three separate clades, corresponding to the above three subtypes. A pairwise sequence comparison using BLASTN showed that these 13 plasmids displayed >99% nucleotide identity across ≥50% of their backbone sequences within each of subgroup IncC plasmids, and had >97% nucleotide identity across ≥47% of their backbone sequences between different subgroups (Figure 1B and Supplementary Table S2). The above results confirmed a parallel diversification and evolution of type 1, type 2, and type 1/2-hybrid IncC plasmids. These 13 plasmids varied in size from about 90 kb to nearly 206 kb with variation in the number of predicted ORFs from 118 to 246, and each plasmid was dissected into the IncC backbone regions and the accessory modules, which were defined as acquired DNA regions associated and bordered with mobile elements (Table 2 and Supplementary Figure S1). These 13 plasmids shared common IncC backbone genes repA, parAB, tra2, dcm2, yacC, int, kfrA, and uvrD.

TABLE 1

Type 1, 2, and 1/2-hybrid plasmids analyzed.

IncC subtype	Plasmid	Accession number	orf1832/orf1847	rhs1/rhs2	i1	i2	Host bacterium	Location	References
Type 1 reference	pR148	JX141473	orf1832	rhs1	−	−	Aeromonas hydrophila	Thailand	Del Castillo et al., 2013
Type 1	pQD1501-Ct1	MN310375	orf1832	Δrhs1	−	−	Klebsiella quasipneumoniae	China	This study
Type 1	p12085-Ct1	MN310377	-	Δrhs1	−	−	K. pneumoniae	China	This study
Type 2 reference	pR55	JQ010984	orf1847	rhs2	+	+	K. pneumoniae	France	Benoît et al., 2012
Type 2	pA2293-Ct2	MN310378	orf1847	Δrhs2	+	+	K. pneumoniae	China	This study
Type 2	pA1763-Ct2	MG764552	-	Δrhs2	+	+	K. pneumoniae	China	This study
Type 2	p24845-Ct2	MF344572	orf1847	Δrhs2	+	+	Enterobacter hormaechei	China	This study
Type 2	p397108-Ct2	MH917284	-	rhs2	+	−	K. pneumoniae	China	This study
Type 2	pT5282-Ct2	MF344574	orf1847	rhs2	−	+	E. hormaechei	China	This study
Type 1/2 hybrid	p205880-Ct1/2	MF344573	orf1832	rhs1	−	+	K. pneumoniae	China	This study
Type 1/2 hybrid	p427113-Ct1/2	MH917285	orf1832	rhs1	−	+	K. pneumoniae	China	This study
Type 1/2 hybrid	pKpn47-Ct1/2	MN310369	orf1832	rhs1	−	+	K. pneumoniae	China	This study
Type 1/2 hybrid	p11935-Ct1/2	MN310370	orf1832	rhs1	−	+	K. pneumoniae	China	This study

FIGURE 1

Evolutionary relationships of IncC plasmids. (A) A maximum-likelihood phylogenetic tree. The tree is constructed from recombination-free SNPs of core backbone regions. The degree of support (percentage) for each cluster of associated taxa, as determined by bootstrap analysis, is shown next to each branch. The bar corresponds to the scale of sequence divergence. Triangles indicate the references for type 1 and type 2 IncC plasmids. (B) A heatmap of pairwise comparison of plasmid backbone sequences. The original BLAST coverage and nucleotide identity values are shown in Supplementary Table S2.

TABLE 2

Major features of analyzed IncC plasmids.

Plasmid	Total length (bp)	Total number of ORFs	Mean G + C content (%)	Length of backbone (bp)	Mean G + C content of backbone (%)	Accessory modules
						Resistance				Non-resistance
						ARI-A	ARI-B	blaCMY-carrying region	Others
pR148	165,906	195	52.5	127,803	51.2	Tn6358	−	–	–	–
pQD1501–Ct1	132,407	161	51.9	113,444	51.2	+	−	Tn6538b	–	–
p12085–Ct1	90,252	122	51.7	79,675	51.1	−	−	blaCMY–6–blaNDM–1 region	−	−
pR55	170,810	203	53.0	129,210	51.2	−	+	−	Tn6187	−
pA2293–Ct2	181,726	205	51.0	124,622	51.1	−	+	−	−	ISKpn52:IS30, ISAs25, and IS100kyp
pA1763–Ct2	93,441	118	51.2	67,327	51.0	−	+	−	blaCTX–M–3 region	−
p24845–Ct2	167,619	192	49.9	89,202	51.0	−	−	−	MDR region	ΔTn6292–ΔIS26–IS26 region
p397108–Ct2	105,977	133	52.5	64,665	50.7	−	+	−	Tn6558	ISKpn18
pT5282–Ct2	152,215	173	52.8	94,867	51.2	−	+	−	blaCTX–M–14–blaIMP–8 region, and MDR region	ISEhe3, and ISEhe3
p205880–Ct1/2	153,373	182	52.3	117,240	51.2	Tn6395	+	Tn6538a	−	−
p427113–Ct1/2	205,674	246	53.3	116,457	51.3	+	+	Tn6538a	MDR region	−
pKpn47–Ct1/2	154,204	178	52.4	116,460	51.3	+	+	Tn6538c	−	ISKpn18, IS1R, and ISKpn18
p11935–Ct1/2	154,204	181	52.5	116,459	51.3	+	+	Tn6538c	−	ISKpn18, IS1R, and ISKpn18

Major Modular Differences in Backbone and Accessory Regions

At least seven major modular differences were identified across the whole plasmid sequences (Figure 2 and Supplementary Figure S2). First, two small regions (dsbA–orf429, and orf858–orf192/orf189–nuc) displayed <95% nucleotide identity between all the three type 1 plasmids and all the 10 type 1/2-hybrid and type 2 plasmids, for which a gene substitution (orf192 in the former three, and orf189 in the later ten) occurred. Second, the ARI-B island or the ΔTn6292–ΔIS26–IS26 region was inserted into the backbone dcm1 region of all the 10 type 1/2-hybrid and type 2 plasmids, leading to a 4.5- to 11.0-kb deletion within the dcm1 region of these 10 plasmids expect for pR55. Third, the three type 1 plasmids had a backbone gene orf489 showing <90% nucleotide identity to all the other plasmids; ISKpn18 was inserted at a site within the backbone gene orf501 in p397108-Ct2, resulting in a 60.1-kb deletion (containing the whole tra1 region and the dcm3 region). Fourth, compared to the prototype tra1 region as observed in pR148 and additional two plasmids, various insertions of accessory modules were found in the other plasmids: (i) Tn6538a or Tn6538b or Tn6538c or the blaCMY–6–blaNDM–1 region was inserted downstream of tivF1 of all the seven type 1 and type 1/2-hybrid plasmids except for pR148, resulting in a 48.1-kb deletion in p12085-Ct1; (ii) the insertion of ISEhe3 into orf672 of pT5282-Ct2 or that of blaCTX–M–3 region into rlx of pA1763-Ct2 led to a 24.5-kb or 57.4-kb deletion, respectively; and (iii) ISKpn18 plus IS1R was inserted into tivF6 of pKpn47-Ct1/2 or p11935-Ct1/2, while the MDR (multidrug resistance) region was inserted into the same gene of p427113-Ct1/2 or p24845-Ct2, leading to a 35.5-kb deletion in p24845-Ct2. Fifth, the backbone dcm3 region was completely or partially lost in five plasmids, which was caused by accessory module insertion or other unknown reasons. Sixth, the ARI-A island was inserted into the backbone gene orf240 of all the seven type 1 and type 1/2-hybrid plasmids except for p12085-Ct1, while ISKpn52:IS30 plus ISAs25 was inserted into a type 2 plasmid pA2293-Ct2. Seventh, Tn6187, IS100kyp and Tn6558 were inserted at different sites around the backbone ter gene of pR55, pA2293-Ct2, and p397108-Ct2, respectively; additionally, three separate insertions of the MDR region, the blaCTX–M–14–blaIMP–8 region, and the ISEhe3 element into the backbone ter–to–tivF3 region made its segmentation in pT5282-Ct2.

FIGURE 2

Major modular differences among IncC plasmids. Genes are denoted by arrows. Genes, mobile elements, and other features are colored based on function classification. Shading denotes regions of homology (≥95% nucleotide identity).

ARI-A and ARI-B Islands

As found in all the seven type 1 and type 1/2-hybrid plasmids except for p12085-Ct1 (type 1), the ARI-A islands (Figure 3) were identified as Tn1696 (Partridge et al., 2001) derivatives. The Tn1696 derivatives in pR148 and p205880-Ct1/2 were recognized as intact unit transposons Tn6358 (Ma et al., 2017) and Tn6395, respectively, which differed from Tn1696 mainly by insertion of In834 (Ma et al., 2017) and In1212 instead of In4 into a primary tnpAR–mer structure. The Tn1696 derivatives in the other four plasmids could not be discriminated as intact transposons due to 3′-terminal truncations, which were in some cases accompanied with insertion of additional resistance modules: blaNDM–1-carrying ΔTn125 (Poirel et al., 2012) in pQD1501-Ct1 and strAB-carrying ΔTn6029 (Cain et al., 2010) in p427113-Ct1/2.

FIGURE 3

Organization of ARI-A islands and comparison to related regions. Genes are denoted by arrows. Genes, mobile elements, and other features are colored based on their functional classification. Shading denotes regions of homology (nucleotide identity ≥95%). Numbers in brackets indicate nucleotide positions within corresponding plasmids. The accession numbers of Tn1696 (Partridge et al., 2001), Tn125 (Poirel et al., 2012), Tn1721 (Allmeier et al., 1992), and Tn6029 (Cain et al., 2010) for reference are U12338, JN872328, X61367, and HQ840942, respectively.

The ARI-B islands (Figure 4) were found in all the 10 type 2 and type 1/2-hybrid plasmids except for p24845-Ct2 (type 2), and all of them contained a sul2 gene that was originated from either GIsul2 (Nigro and Hall, 2011) or Tn6029 (Cain et al., 2010), because various sul2-carrying remnants of GIsul2 or Tn6029 were found in these ARI-B islands. In general, these ARI-B islands acquired dramatically different collection of resistance modules, which contained genes resistance to older (such as dfrA, strAB, tetA, and floR) and newer (such as blaCTX–M, armA, msr, and mph) antibiotics, and thus varied considerably in size with a highly mosaic nature. The ΔTn6292–ΔIS26–IS26 region (Figure 4), which did not contain any of resistance genes, was inserted at the site specific for ARI-B insertion.

FIGURE 4

Organization of ARI-B islands and comparison to related regions. Genes are denoted by arrows. Genes, mobile elements, and other features are colored based on their functional classification. Shading denotes regions of homology (nucleotide identity ≥95%). Numbers in brackets indicate nucleotide positions within corresponding plasmids. The accession numbers of Tn1696 (Partridge et al., 2001), Tn1548 (Galimand et al., 2005), aacC2–tmrB region (Partridge et al., 2012), Tn7 (see Zhan et al., 2018, for gene organization), Tn2 (Bailey et al., 2011), Tn6029 (Cain et al., 2010), GIsul2 (Nigro and Hall, 2011), Tn6292 (Wei et al., 2016), and Tn1722 (Allmeier et al., 1992) for reference are U12338, AF550415, JX101693, KX117211, HM749967, HQ840942, AE014073, KU886034, and X61367, respectively.

Other Accessory Resistance Modules

Four different blaCMY-carrying regions (Figure 5) were found in all the seven type 1 and type 1/2-hybrid plasmids except for pR148 and were manifested as three intact ISEcp1-based transposition units Tn6538a (p205880-Ct1/2 and p427113-Ct1/2), Tn6538b (pQD1501-Ct1), and Tn6538c (pKpn47-Ct1/2 and p11935-Ct1/2), and the blaCMY–6–blaNDM–1 region (p12085-Ct1). Tn6538a, Tn6538b, and Tn6538c were highly similar to one another but differed by substitution of blaCMY–2 with blaCMY–6 or by insertion of IS1R into ISEcp1. The blaCMY–6–blaNDM–1 region was a combination of ΔTn6538b carrying blaCMY–6, ISKpn14, and ΔTn125 containing blaNDM–1.

FIGURE 5

Organization of blaCMY-carrying regions and comparison to related regions. Genes are denoted by arrows. Genes, mobile elements and other features are colored based on their functional classification. Shading denotes regions of homology (nucleotide identity ≥95%). Numbers in brackets indicate nucleotide positions within corresponding plasmids. The accession number of Tn125 (Poirel et al., 2012) for reference is JN872328.

Four different blaCTX–M-carrying regions were found in four type 2 plasmids: Tn6558 (p397108-Ct2), the blaCTX–M–14–blaIMP–8 region (pT5282-Ct2), and the blaCTX–M–3 region (pA1763-Ct2) (Figure 6), as well as the ARI-B island (pA2293-Ct2) (Figure 4). Insertion of a truncated version of Tn6503, which was an ISEcp1-based, blaCTX–M–14-carrying transposition unit (Feng et al., 2015), into the cryptic unit transposon Tn1722 generated Tn6558, while the blaCTX–M–14–blaIMP–8 region was resulted from three separate insertions of Tn6503, In655 (carrying blaIMP–8), and IS26 into Tn1722. The blaCTX–M–3 region in pA1763-Ct2 was highly similar to a partial region of the ARI-B island in pA2293-Ct2, and both carried the ISEcp1-based blaCTX–M–3-carrying transposition unit Tn6502 (Liang et al., 2014).

FIGURE 6

Organization of blaCTX–M-carrying regions and comparison to related regions. Genes are denoted by arrows. Genes, mobile elements and other features are colored based on their functional classification. Shading denotes regions of homology (nucleotide identity ≥95%). Numbers in brackets indicate nucleotide positions within corresponding plasmids. The accession numbers of Tn1722 (Allmeier et al., 1992), Tn6292 (Wei et al., 2016) and Tn1696 (Partridge et al., 2001) for reference are X61367, KU886034 and U12338, respectively.

Four different MDR regions (Figure 7) were found in four plasmids: Tn6187 (Benoît et al., 2012) in pR55, and three MDR regions in p427113-Ct1/2, pT5282-Ct2, and p24845-Ct2. As similar to the above ARI-B islands, different MDR regions were assembled from different collections of resistance units [especially including ΔTn6296 carrying blaKPC–2, IS26–mph(A)–mrx–mphR(A)–IS6100 unit, In207, In384, and rmtB–tetA(G)–tetR(G)–floR region] via homologous or non-homologous recombination.

FIGURE 7

Organization of MDR regions and comparison to related regions. Genes are denoted by arrows. Genes, mobile elements and other features are colored based on their functional classification. Shading denotes regions of homology (nucleotide identity ≥95%). Numbers in brackets indicate nucleotide positions within corresponding plasmids. The accession number of Tn1696 (Partridge et al., 2001), Tn21 (Partridge et al., 2001), Tn1548 (Galimand et al., 2005), Tn2 (Bailey et al., 2011), aacC2–tmrB region (Partridge et al., 2012) and Tn6296 (Jiang et al., 2010) for reference are U12338, AF071413, AF550415, HM749967, JX101693 and FJ628167, respectively.

Transferability and Antimicrobial Susceptibility

Plasmids pQD1501-Ct1, pA2293-Ct2, and p205880-Ct1/2 were arbitrarily selected as the representatives of type 1, type 2, and type 1/2 hybrid, respectively, and could be transferred from the wild-type isolates into J53 or EC600 through conjugation, generating E. coli transconjugants QD1501-NDM-J53, A2293-CTXM-J53, and 205880-CMY-EC600, respectively. The self-transferable nature was consistent with the presence of two complete sets (tra1 and tra2) of conjugal transfer genes in each plasmid, although Tn6538a or Tn6538b was inserted at an intergenic site in tra1 of p205880-Ct1/2 or pQD1501-Ct1, respectively. All the three wild-type isolates QD1501, A2293, and 205880 and one transconjugant QD1501-NDM-J53 were resistant to imipenem with minimum inhibitory concentration (MIC) values ≥16 (due to production of NDM or KPC enzyme) but the other two transconjugants A2293-CTXM-J53 and 205880-CMY-EC600 remained susceptible to imipenem. All the above strains were resistant to cefuroxime with MIC values ≥32 (due to production of NDM or CTX-M or CMY enzyme).

Discussion

Comparative genomics analysis of the collection of IncC plasmids (n = 13 in this study) provides a deeper understanding of diversification of IncC plasmids. Different IncC plasmids have evolved to acquire very different profiles of accessory modules and thus different collections of resistance genes (Figure 8). Insertion of accessory modules is often accompanied by deletion of surrounding backbone regions. Type 1/2-hybrid IncC plasmids most likely emerged from homologous recombination between type 1 and type 2 plasmids. In addition to using the presence or absence of signature sequences orf1832/orf1847, rhs1/rhs2, i1, and i2, a phylogenetic analysis using recombination-free SNPs of core backbone regions can steadily discriminate type 1/2 hybrid from type 1 and type 2. As shown previously (Ma et al., 2017) and herein, additional resistance islands (e.g., the MDR regions and the blaCTX–M-carrying regions in this study) can be found in type 1, type 2, and type 1/2-hybrid plasmids.

FIGURE 8

A heatmap of prevalence of resistance genes in IncC plasmids. The original data are shown in Supplementary Table S3.

The IncC reference plasmids pR148 (166 kb in length, 195 ORFs) and pR55 (171 kb in length, 203 ORFs) had the most complete backbones for type 1 and type 2, respectively. By contrast, the type 1 plasmid p12085-Ct1 and the type 2 plasmids pA1763-Ct2 and p397108-Ct2 harbored extremely simple backbones due to deletion of large backbone regions. In addition to insertion of accessory modules and deletion of backbone regions, large inversion (such as those in p427113-Ct1/2 and pA1763-Ct2) and translocation (such as those in pT5282-Ct2) events occurred within backbone regions during genomic diversification of IncC plasmids (Ma et al., 2017; Supplementary Figure S2).

The carbapenemase genes blaNDM–1, blaKPC–2, and blaIMP–8 were identified in pQD1501-Ct1 and p12085-Ct1, p427113-Ct1/2, and pT5282-Ct2, respectively. This is the first report of a blaIMP–8-carrying IncC plasmid. This study presents three novel mobile elements [firstly identified in this study, including Tn6395 (p205880-Ct1/2), In2-76 (pA1763-Ct2), and ISEcl10 (p24845-Ct2)] and additionally six newly named mobile elements [firstly designated in this study but with previously determined sequences, including Tn6538a (p205880-Ct1/2 and p427113-Ct1/2), Tn6538b (pQD1501-Ct1), Tn6538c (pKpn47-Ct1/2 and p11935-Ct1/2), Tn6502 (pA2293-Ct2 and pA1763-Ct2), Tn6503 (pT5282-Ct2), and ISKpn52 (pA2293-Ct2)].

Data Availability Statement

The datasets generated for this study can be found in the complete nucleotide sequences of plasmids p24845-Ct2, p205880-Ct1/2, pT5282-Ct2, pKpn47-Ct1/2, p11935-Ct1/2, pQD1501-Ct1, p12085-Ct1, pA2293-Ct2, p397108-Ct2, p427113-Ct1/2, and pA1763-Ct2 were submitted to GenBank under the accession numbers MF344572 to MF344574, MN310369, MN310370, MN310375, MN310377, MN310378, MH917284, MH917285, and MG764552, respectively.

Ethics Statement

This study uses the clinical bacterial isolates obtained from the Chinese public hospitals as listed in Supplementary Table S1. The local legislation did not require the study to be reviewed or approved by an ethics committee, because the bacterial isolates involved in this study were part of the routine hospital laboratory procedures. The research involving biohazards and all related procedures were approved by the Biosafety Committee of the Beijing Institute of Microbiology and Epidemiology.

Author Contributions

ED and DZ conceived the study and designed the experimental procedures. QC, XJ, ZY, YX, LH, and HG performed the experiments. QC, WL, XJ, DZ, and ED analyzed the data. WY, HY, YZ, and XZ provided the reagents and materials. DZ, QC, and ED wrote the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.