IncHI2 Plasmids Are Predominant in Antibiotic-Resistant Salmonella Isolates

The wide usage of antibiotics contributes to the increase in the prevalence of antibiotic-resistant Salmonella. Plasmids play a critical role in horizontal transfer of antibiotic resistance markers in Salmonella. This study aimed to screen and characterize plasmid profiles responsible for antibiotic resistance in Salmonella and ultimately to clarify the molecular mechanism of transferable plasmid-mediated antibiotic resistance. A total of 226 Salmonella isolates were examined for antimicrobial susceptibility by a disk diffusion method. Thirty-two isolates (14.2%) were resistant to at least one antibiotic. The presence of plasmid-mediated quinolone resistance (PMQR) genes and β-lactamase genes were established by PCR amplification. PCR-based replicon typing revealed that these 32 isolates represented seven plasmid incompatibility groups (IncP, HI2, A/C, FIIs, FIA, FIB, and I1), and the IncHI2 (59.4%) was predominant. Antibiotic resistance markers located on plasmids were identified through plasmid curing. Fifteen phenotypic variants were obtained with the curing efficiency of 46.9% (15/32). The cured plasmids mainly belong to the HI2 incompatibility group. The elimination of IncHI2 plasmids correlated with the loss of β-lactamase genes (blaOXA-1 and blaTEM-1) and PMQR genes (qnrA and aac(6′)-Ib-cr). Both IncHI2 and IncI1 plasmids in a S. enterica serovar Indiana isolate SJTUF 10584 were lost by curing. The blaCMY -2-carrying plasmid pS10584 from SJTUF 10584 was fully sequenced. Sequence analysis revealed that it possessed a plasmid scaffold typical for IncI1 plasmids with the unique genetic arrangement of IS1294-ΔISEcp1-blaCMY -2-blc-sugE-ΔecnR inserted into the colicin gene cia. These data suggested that IncHI2 was the major plasmid lineage contributing to the dissemination of antibiotic resistance in Salmonella and the activity of multiple mobile genetic elements may contribute to antibiotic resistance evolution and dissemination between different plasmid replicons.

The wide usage of antibiotics contributes to the increase in the prevalence of antibioticresistant Salmonella. Plasmids play a critical role in horizontal transfer of antibiotic resistance markers in Salmonella. This study aimed to screen and characterize plasmid profiles responsible for antibiotic resistance in Salmonella and ultimately to clarify the molecular mechanism of transferable plasmid-mediated antibiotic resistance. A total of 226 Salmonella isolates were examined for antimicrobial susceptibility by a disk diffusion method. Thirty-two isolates (14.2%) were resistant to at least one antibiotic. The presence of plasmid-mediated quinolone resistance (PMQR) genes and β-lactamase genes were established by PCR amplification. PCR-based replicon typing revealed that these 32 isolates represented seven plasmid incompatibility groups (IncP, HI2, A/C, FIIs, FIA, FIB, and I1), and the IncHI2 (59.4%) was predominant. Antibiotic resistance markers located on plasmids were identified through plasmid curing. Fifteen phenotypic variants were obtained with the curing efficiency of 46.9% (15/32). The cured plasmids mainly belong to the HI2 incompatibility group. The elimination of IncHI2 plasmids correlated with the loss of β-lactamase genes (bla OXA and bla ) and PMQR genes (qnrA

INTRODUCTION
Salmonella is recognized as the predominant pathogen implicated in bacterial foodborne diseases worldwide (Yang et al., 2016). The emergence of multidrug-resistant (MDR) Salmonella strains due to the prolonged and extensive usage of antibiotics has become a public health issue. Such resistance toward medically significant antimicrobial agents including fluoroquinolones and extended-spectrum cephalosporins (ESCs) that are regarded as primary treatment options for bacterial infections will make it difficult for antibiotic therapy due to the reduced efficiency of empirical strategies as well as limited choice of treatment (Wong et al., 2014;Folster et al., 2015). The dissemination of undesirable antibiotic resistance in Gram-negative bacteria has been largely attributed to the acquisition of multiple plasmid-located antibiotic resistance genes by horizontal gene transfer (HGT; Carattoli, 2013;Wang et al., 2015). The emergence of plasmid-mediated quinolone resistance (PMQR) predominantly encoded by qnr variants, aac(6 )-Ib-cr, qepA, and oqxAB genes in Salmonella, has further facilitated selection of high-level chromosomal quinolone resistance via mutations in topoisomerase genes, as well as the wide spread of quinolone resistance via HGT (Jiang et al., 2014;Skov et al., 2015). In addition, the apparent correlation between the PMQR genes (such as qnr or aac(6 )-Ib-cr genes) and extended-spectrum β-lactamase (ESBL) genes or ampC genes on the same plasmid and their association with other antibiotic resistance genes have raised particular global concern, since their transmission could be driven by other mobile genetic elements located on plasmids among various bacterial species Jiang et al., 2014). Particular plasmid families associate with the emergence and dissemination of specific antibiotic resistance traits with differential prevalence and distribution (Johnson et al., 2007;Carattoli, 2013). Thus, characterization of plasmids based on PCR-based replicon typing (PBRT) is an indispensable part of plasmid epidemiological surveillance enhancing discrimination between Salmonella strains as well as tracing the spread and evolution of antibiotic resistance genes . Certain replicon types were found to be associated with MDR as well as with bacterial disease outbreaks (Huang et al., 2012). It has been reported that plasmids with IncA/C, B/O, HI1, HI2, I1, N, F, and P replicons are often associated with MDR in Salmonella, while many of them are found to be co-resident in some MDR Salmonella strains (Poole et al., 2009;Glenn et al., 2013).
Transmissible plasmids contain a full set of conjugationencoding genes facilitating their spread over large taxonomic distances, and they harbor drug-resistance determinants, virulence factors and addiction systems promoting their stability and maintenance in bacterial hosts under different environmental conditions. It is essential to clarify the roles of these genetic elements and the corresponding surrounding genetic structures by plasmid sequencing in order to trace the transmission and persistence of antibiotic resistance (Carattoli, 2013;Wang et al., 2015).
The aim of this study was to screen and characterize plasmid replicon profiles containing PMQR and β-lactamase genes in Salmonella so as to clarify molecular mechanism of transferable plasmid-mediated antibiotic resistance.

Salmonella Isolates
A total of 226 Salmonella isolates were used in this study, and all the non-susceptible strains are listed in Supplementary   Table S1, of which 39 were food isolates and 39 were clinical isolates. Among these Salmonella isolates, clinical isolates were collected by Shanghai Municipal Center for Disease Control and Prevention and Wuhan Municipal Center for Disease Control and Prevention, while food isolates were collected from beef, poultry, pork, shrimp, vegetables, fresh juice, and shellfish. PCR-based serotyping was performed as previously described (Ranieri et al., 2013). Two different primer pairs and an additional primer pair were used in this study (Supplementary Table  S2). The detailed information of these isolates is listed in Supplementary Table S1.

Plasmid Replicon Typing
Plasmid incompatibility (Inc) groups were assigned by PBRT using genomic DNA of Salmonella isolates as template. Amplification by PCR was performed with eighteen specific primer pairs designed for FIA, FIB, FIC, HI1, HI2, I1, L/M, N, P, W, T, A/C, K, B/O, X, Y, F, and FIIA basic replicons, and the corresponding protocol was previously described . Resulting amplicons were sequenced by Shanghai Majorbio Bio-pharm Technology Co., Ltd. Comparative analysis of nucleotide sequences was performed using the BLAST program at the National Center for Biotechnology Information (NCBI) site 1 .

Plasmid Curing
To determine if the antibiotic resistance observed in this study was plasmid-mediated, the identified antibiotic-resistant Salmonella isolates were subjected to a curing treatment using sodium dodecyl sulphate (SDS), according to Tomoeda et al. (1968) with modifications. Briefly, the overnight culture in Luria-Bertani (LB) broth was inoculated (1:10 4 dilution) in LB broth containing 5% SDS and incubated for 48 h at 44.5 • C with constant agitation. Aliquots of the cultures were then appropriately diluted and spread on LB agar plates without antibiotics and incubated at 37 • C for 16-18 h. The resulting single clones were streaked onto two LB agar plates with or without ampicillin. Because ampicillin was the most prevalent among antibiotic-resistant isolates, for the screening convenience of those plasmid-loss mutants, ampicillin was used as the priority selective marker. The LB agar plates were incubated at 37 • C overnight, and ampicillin-susceptible clones were considered positive for loss of the incompatible plasmid. After treatment with the curing agent, the putative cured derivatives were examined for the antibiotic susceptibility profile using the same set of antibiotics, and for the presence of related antibiotic resistance genes as well as plasmid replicon type by PCR method as described above.

Plasmid Sequencing and Analysis
The Salmonella isolate SJTUF 10584 identified in this study was chosen as the plasmid donor for its complete loss of plasmids and MDR after one single trial of plasmid curing. The plasmid was obtained through conjugation experiments by the liquid mating assay using rifampin-resistant E. coli Nk5449 as the recipient strain as previously described (Zurfluh et al., 2014). Transconjugants were selected on LB agar plates containing ampicillin (100 µg/ml) and rifampin (200 µg/ml). Antimicrobial susceptibility tests and antibiotic resistance gene detection were performed to confirm the plasmid transfer, followed by PBRT to test which plasmid and antibiotic resistance markers were transferred.
Plasmid DNA was purified from the transconjugant by the Qiagen Plasmid Midi Kit (Qiagen, Germany) according to the manufacturer's protocol. Plasmid sequencing was performed on the IIIumina Hiseq sequencing platform at the Shanghai Biotechnology Corporation (Shanghai, China). The complete sequence was annotated using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP), followed by manual inspection. Each predicted protein was further confirmed against the NCBI all-protein database using BLASTP 1 . Insertion sequences and repetitive elements were identified using IS finder 2 . DNA sequence comparisons and alignments were performed using BLAST 1 and Geneious software (version 9.1.2; Kearse et al., 2012), and a schematic plasmid map was constructed with WinPlas 2.7 software. 1 http://blast.ncbi.nlm.nih.gov/Blast.cgi 2 https://www-is.biotoul.fr//

Nucleotide Sequence Accession Number
The complete sequence of the plasmid pS10584 has been deposited in Genbank/EMBL/DDBJ with accession number KX058576.

Plasmid Curing and Antibiotic Resistance Localization
All the 32 antibiotic-resistant Salmonella isolates underwent plasmid curing using 5% SDS as well as elevated temperature (44.5 • C) for 48 h to ensure that the antibiotic resistance profile was related to the plasmid. Fifteen phenotypic variants were obtained after up to 10 trials of curing (Figure 2), and the curing efficiency was 46.9% (15/32). The cured plasmids mainly belonged to IncHI2 group (7/8, 87.5%), followed by FIGURE 2 | A heat-map showing the changes of Salmonella isolates after plasmid curing based on antibiotic resistance profiles as well as plasmid replicon type(s). "C" denoted the cured variant. Black and white squares denotes the presence and absence of a particular feature, respectively. The antibiotics listed are abbreviated as follows: AMP, ampicillin; TZP, piperacillin/tazobactam; SAM, ampicillin/sulbactam; CRO, ceftriaxone; CAZ, ceftazidime; ATM, aztreonam; CZO, cephazolin; CIP, ciprofloxacin; AMK, amikacin; GEN, gentamicin; TOB, tobramycin; LEV, levofloxacin; SXT, sulfamethoxazole/trimethoprim. IncF group (3/3, 100%), and IncI1 group (1/1, 100%); however, no strain showed curing of IncP group plasmids (0/5, 0%). When attributing the loss of resistance phenotypes to the loss of plasmids, we could infer that the antibiotic resistant phenotypes AMP-SAM-TOB of SJTUF 10236, AMP-SAM-GEN-TOB-LEV-SXT of SJTUF 10476, AMP-SAM of SJTUF 10479, AMP-SAM-GEN-TOB-CZO of SJTUF 10565, AMP of SJTUF 10577, and AMP-SAM-TOB-SXT of SJTUF 10580 were all related to IncHI2 plasmids while AMP-SAM-CZO-CAZ-CRO-CIP-SXT of SJTUF 10584 was determined together by IncI1 and IncHI2 plasmids. And the resistant phenotype AMP-SAM-CZO-GEN-CIP-LEV of SJTUF 10713 and AMP-GEN-TOB of SJTUF 10718 were both determined by IncF plasmids. In addition, the elimination of IncHI2 plasmids correlated with the loss of β-lactamase genes (bla OXA−1 and bla TEM−1 ) as well as PMQR genes (qnrA and aac(6 )-Ib-cr), and the loss of IncFIA and FIB plasmids resulted in the absence of the aac(6 )-Ib-cr gene. Furthermore, four cured variants (SJTUF 10211 C , 10700 C , 10702 C and 10705 C ) lost antibiotic resistance genes as well as drug resistance phenotypes but still retained the identified plasmids (Figure 2). This may be due to the loss of divergent or novel replicons which were untypeable by the current PBRT scheme. The cured SJTUF 10577 C lost the AMP resistance phenotype, bla genes, and the IncHI2 plasmid simultaneously, even though it retained the IncP plasmid. The function of IncP plasmids in antibiotic resistance needs further study.
Particular plasmids belonging to the Inc families usually carry multiple physically linked genetic determinants, playing a major role in the diffusion of specific antibiotic resistance genes among different bacteria (Carattoli, 2013;Wang et al., 2013). In this study, we characterized IncHI2 plasmids as the main vehicle in horizontal transfer of β-lactamase genes (bla OXA−1 and bla TEM−1 ) as well as PMQR genes (qnrA and acc(6 )-ib-cr). The co-existence of aac(6 )-Ib-cr with bla OXA−1 , with or without bla TEM−1 was also identified. IncHI2 plasmids are known to play a significant role in the acquisition of antibiotic resistance, and they have recently been implicated in the dissemination of genes encoding extended-spectrum β-lactamase (Cain and Hall, 2012). In contrast, IncA/C plasmids are prevalent in several MDR Salmonella serovars and closely linked to the expansion of MDR in the United States (Folster et al., 2010(Folster et al., , 2015. We also discovered the high prevalence of IncHI2 plasmids in S. Typhimurium, similar to the report from Li et al. (2013). The identical genetic context of antibiotic resistance genes on IncHI2 plasmids was also detected in S. Indiana in China  as well as S. Typhimurium in Europe (Campos et al., 2015), suggesting a similar evolutionary origin and highlighting the potentially global spread of IncHI2 plasmids among Salmonella. Furthermore, a similar genetic context of aac(6 )-Ib-cr with bla OXA−1 was also identified on an IncR plasmid of a K. oxytoca strain in Spain (Ruiz et al., 2011), and an IncN plasmid of an E. coli isolate in Hong Kong (Ho et al., 2013). The similar resistance gene module infers that the transfer of resistanceassociated module was driven by mobile genetic elements like integrons, insertion sequences or transposons between different plasmid replicons, and was probably mediated by IS26 in terms of the above cases (Miriagou et al., 2005;Li et al., 2013;He et al., 2015). Thus, further active surveillance is needed to minimize the spread of particular plasmids such as IncHI2 group to control the dissemination of antibiotic resistance. As mobile genetic elements play a major role in the acquisition and dissemination of antibiotic resistance genes, further study may be required to determine the distinct genetic environment of antibiotic resistance genes on plasmids.

Analysis of pS10584
The MDR Salmonella enterica serovar Indiana isolate SJTUF 10584 identified in this study was chosen as the plasmid donor. This isolate contains IncHI2 and IncI1 plasmids simultaneously, and these two types of plasmids were easily eliminated after one trial of curing. The plasmid designated as pS10584 was acquired through conjugation experiments recovered from ampicillin and rifampin resistance selection. Replicon typing showed that the only transferred plasmid pS10584 belonged to the IncI1 group and PCR analysis confirmed that it was positive for bla CMY−2 gene but not for bla OXA−1 or aac(6 )-Ib-cr genes. All of the transconjugants shared the same antibiotic resistance pattern of AMP-SAM-CZO-CAZ-CRO.
Complete DNA sequence of pS10584 revealed that it was a circular molecule of 94,697-bp with an average G+C content of 50.0% and 110 putative CDSs. In general, pS10584 possessed a plasmid scaffold typical for IncI1 plasmids, with important IncI1-associated genetic modules including traABCD regulatory cluster, traLMNOPQRSTUVWXY transfer region, nikB-trbABC region, ssb-psiAB region, and type IV thin pilus formation region pilIJKLMNOPQRSTUV (Figure 3). Besides the approximate 89,866-bp "core genome" region encoding characteristic IncI1 plasmid replication, transfer, maintenance and stability functions, a 4,831-bp accessory module encoding the antibiotic resistance gene as well as flanking insertion sequences was also discovered. BLASTN comparison revealed that the overall genetic organization of pS10584 was very similar to that of another fully sequenced IncI1 plasmid pJIE512b (GenBank accession no. HG970648; Tagg et al., 2014) with high sequence identity (>99% at nucleotide level) except an inversion within the pilV shufflon and a 2,307-bp deletion (nucleotide 6,685 th to 8,991 nd of pS10584, containing a truncated colicin gene cia and an intact putative DNA-binding transcriptional regulator gene yagA) in pJIE5112b.
In pS10584, the 4,831-bp antibiotic resistance-associated module comprised 161-bp of the right end of ISEcp1 truncated by IS1294 in the reverse orientation and a 2,823-bp region originating from the Citrobacter freundii chromosome including a bla CMY−2 gene, with blc (an outer membrane lipoprotein, lipocalin), sugE (a small MDR transporter), and ecnR (LuxR family transcriptional regulator) genes immediately downstream, followed by a 159-bp fragment of IncA/C ending in GTTC which matched the last 4-bp of terIS of IS1294 (Figure 4). In addition, the 4831-bp bla CMY−2 segment inserted within the cia gene at the 698-bp from the 5 end without causing adjacent deletion, followed by the remaining 1,184-bp segment of the cia gene.
A similar genetic context surrounding the bla CMY−2 gene was identified in an IncI1 plasmid p2735 (KP017243; Yassine et al., 2014) isolated from Klebsiella pneumoniae as well as six IncI1 plasmids isolated from E. coli including two completely sequenced plasmids pJIE512b (HG970648; Tagg et al., 2014), and pR7AC (KF434766; Minoia et al., 2012), and four partially sequenced plasmids pTN38148 (FM246883), pTN10, pTN13, and pTN386 (Verdet et al., 2009) (Figure 4). Compared to pS10584, the bla CMY−2 insertion was within the yagA gene 16-bp from the 5 end, removing the rest of the yagA gene, as well as a large part FIGURE 3 | An overview of the bla CMY−2 -carrying plasmid pS10584. Genes are color-coded as follows: dark brown, replication-associated genes; green, insertion sequences; red, antibiotic resistance gene; blackish green, partition-or stability-associated genes; army green, transfer-leading region associated genes; pink, tra locus; blue, type IV pilus region; purple, trb locus; and gray, hypothetical protein. of the cia gene, leaving the same part of the cia gene in the 3 end followed by the same 159-bp IncA/C backbone ending in GTTC in pJIE512b. Whereas, the flanking fragments of bla CMY−2 gene in other six IncI1 plasmids (pTN10, pTN13, pTN386, pTN38148, p2735, and pR7AC) were much more different from that in pS10584, including an IS1294-like element designated as IS1294b (Yassine et al., 2014), followed by a 372-bp truncated ISEcp1, the identical 2,823-bp region of the C. freundii chromosome that contained bla CMY−2 and the additional 1395-bp IncA/C backbone also ending in GTTC (Figure 4). This bla CMY−2 insertion was downstream of the repZ gene, removing a larger adjacent part and truncating the 3 end of yafA followed by yafB.
Although S. Indiana frequently exhibited MDR phenotypes, very little research has been done to clarify the molecular mechanisms conferring its MDR or the transferability of the resistance properties . Among ampC-like genes, the bla CMY−2 gene which probably originated from C. freundii, is the most prevalent gene worldwide conferring resistance to clinically important β-lactam-based compounds particularly the third-generation cephalosporins. (Verdet et al., 2009;Call et al., 2010;Karczmarczyk et al., 2014;Tagg et al., 2014). In this study, two bla CMY−2 -positive isolates both showed resistance phenotypes against the third-generation cephalosporins (CRO, CAZ, and CZO) while the other one was susceptible to them, probably due to the low-level expression or silence of this gene. In previous research, it was observed that resistance to a second generation cephalosporin, cefoxitin in Salmonella isolates harboring bla CMY−2 (Giles et al., 2004;Martin et al., 2012;Noda et al., 2015). However, none of the three bla CMY−2 -positive isolates was resistant to cefotetan, which is a second-generation cephalosporin like cefoxitin. This was also found in a previous study (Noda et al., 2015). This may also be attributed to the low-level expression or silence of bla CMY−2 gene.
Among various plasmid types associated with the carriage of the bla CMY−2 gene, IncI1 and IncA/C plasmids are often reported as the predominant carriers (Mata et al., 2012;Bortolaia et al., 2014). Additional research discovered that bla CMY−2 carrying IncA/C plasmids always conferred additional antimicrobial resistance phenotypes such as chloramphenicol and sulfisoxazole, whereas in this study, bla CMY−2 carrying IncI1 plasmids only conferred a bla CMY−2 -associated phenotype, which was also observed in Folster et al. (2011). The bla CMY−2 genetic context on the IncI1 plasmid identified in this study was IS1294-ISEcp1bla CMY−2 -blc-sugE-ecnR followed by 159-bp IncA/C backbone, which was identical to that of pJIE512b from an E. coli isolate (Tagg et al., 2014). The genetic arrangement ISEcp1-bla CMY−2blc-sugE-ecnR was conserved on IncA/C plasmids (Verdet et al., 2009;Ye et al., 2016). ISEcp1 was presumed to be responsible for the mobilization of bla CMY−2 (Verdet et al., 2009) or bla CTX−M−2 (Lartigue et al., 2006) as well as its adjacent region from the bacterial chromosome to the IncA/C plasmids by a one-ended mobilization mechanism. On the basis of sequence alignments, it was suggested that mobilization of bla CMY−2 from the IncA/C plasmid to the IncI1 replicon occurred through IS1294-mediated transposition activity (Tagg et al., 2014). The bla CMY−2 insertion region of pS10584 started in CTTG and ended in GTTC, both of which matched the last 4-bp of oriIS and terIS of IS1294 (Tavakoli et al., 2000;Yassine et al., 2014), implying that IS1294 inserted into ISEcp1 in an IncA/C plasmid and subsequently mobilized the adjacent region, including the bla CMY−2 segment plus a 159bp IncA/C backbone (Tagg et al., 2014). A similar bla CMY−2 context was also observed in six reported IncI1 plasmids p2735 (Yassine et al., 2014), pR7AC (Minoia et al., 2012), pTN38148, TN10, TN13, and TN386 (Verdet et al., 2009); however, in these plasmids the bla CMY−2 insertion segment varied in the length of IncA/C backbone and the position that ISEcp1 interrupted by an IS1294-like element designated as IS1294b. And the integration of the bla CMY−2 gene by IS1294b between IncA/C and IncI1 replicons was experimentally confirmed by Yassine et al. (2014). By comparing the integration sites of bla CMY−2 insertion segments on these reported IncI1 plasmids, we inferred that the downstream of repZ gene could represent a "hot spot" for the integration of IS1294b element in the IncI1 plasmid scaffold but determination of the "hot spot" for IS1294 may require further study.
Besides the resistance determinants IncI1 plasmids encode, additional factors may be responsible for their prevalent diffusion. The shufflon multiple inversion system on IncI1 plasmids involves four segments separated by seven inverted repeats and a shufflon-specific recombinase (Rci) which acts to catalyze recombination between these repeats (Sampei et al., 2010). This shufflon system is responsible for generating variation in the C terminus of the pilV tip adhesins of the pili, leading to the specificity of recipient cells in liquid mating (Tagg et al., 2014). The resulting type IV pilus is also a known virulence factor, contributing to adhesion and invasion of Shiga toxinproducing E. coli (Kim and Komano, 1997), which accounts for the prevalent existence of IncI1 plasmids in pathogenic E. coli (Johnson et al., 2007). Thus, the association of epidemic ability as well as resistance determinants may favor the wide dissemination of plasmids belonging to the IncI1 group (Carattoli, 2009;Folster et al., 2010).

CONCLUSION
We screened and characterized antibiotic resistance plasmids from clinical and food Salmonella isolates in this study and discovered that IncHI2 was the major plasmid lineage contributing to the dissemination of antibiotic resistance in Salmonella. We also reported the complete sequence of an IncI1 plasmid harboring an antibiotic resistance-encoding region with the genetic arrangement of IS1294-ISEcp1-bla CMY−2 -blc-sugE-ecnR, interrupting the cia gene from a MDR S. Indiana isolate. To our knowledge, this is the first report of this bla CMY−2 context on the IncI1 plasmid in Salmonella. These results revealed that plasmids represent a potential threat for the dissemination of antibiotic resistance since they carry relevant resistance determinants and complete transfer-associated genes. Thus, in addition to phenotypic monitoring of antimicrobial resistance in Salmonella, further active surveillance is needed to minimize the spread of particular plasmids such as those belonging to the IncHI2 group. In addition, investigations of mobile genetic elements that are located on plasmids, such as insertion sequences, integrons, and transposons, may aid in understanding the source and spread of antibiotic resistance.

AUTHOR CONTRIBUTIONS
WC completed the antimicrobial susceptibility test, PMQR and β-lactamase gene detection, plasmid replicon typing, plasmid curing, plasmid sequencing, and prepare the manuscript; TF completed the isolate collection, antimicrobial susceptibility test, PMQR and β-lactamase gene detection, plasmid replicon typing, and plasmid curing; XZ helped to finish the isolate collection and antimicrobial susceptibility test; DZ helped to finish the plasmid sequencing and analysis; XS helped to design the project and give some valuable suggestions; and CS designed the project, completed the data analysis, and prepared the manuscript.