ABSTRACT
The incidence of ciprofloxacin resistance in Salmonella has increased dramatically in the past decade. To track the evolutionary trend of ciprofloxacin resistance-encoding genetic elements during this period, we surveyed the prevalence of Salmonella in food products in Shenzhen, China, during the period of 2012 to 2017 and performed whole-genome sequencing and genetic analysis of 566 ciprofloxacin-resistant clinical Salmonella strains collected during this survey. We observed that target gene mutations have become much less common, with single gyrA mutation currently detectable in Salmonella enterica serovar Typhimurium only. Multiple plasmid-mediated quinolone resistance (PMQR) genes located in the chromosome and plasmids are now frequently detectable in ciprofloxacin-resistant Salmonella strains of various serotypes. Among them, the qnrS1 gene was often harbored by multiple plasmids, with p10k-like plasmids being the most dominant. Importantly, p10k-like plasmids initially were not conjugative but became transmissible with the help of a helper plasmid. Ciprofloxacin resistance due to combined effect of carriage of the qnrS1 gene and other resistance mechanisms is common. In S. Typhimurium, carriage of qnrS1 is often associated with a single gyrA mutation; in other serotypes, combination of qnrS1 and other PMQR genes located in the chromosomal fragment or plasmid is observed. Another major mechanism of ciprofloxacin resistance, mainly observable in S. Derby, involves a chromosomal fragment harboring the qnrS2–aac(6′)lb-cr–oqxAB elements. Intriguingly, this chromosomal fragment, flanked by IS26, could form a circular intermediate and became transferrable. To conclude, the increase in the incidence of various PMQR mobile genetic elements and their interactions with other resistance mechanism contribute to a sharp increase in the prevalence of ciprofloxacin-resistant clinical Salmonella strains in recent years.
IMPORTANCE Resistance of nontyphoidal Salmonella to fluoroquinolones such as ciprofloxacin is known to be mediated by target mutations. This study surveyed the prevalence of Salmonella strains recovered from 2,989 food products in Shenzhen, China, during the period 2012 to 2017 and characterized the genetic features of several PMQR gene-bearing plasmids and ciprofloxacin resistance-encoding DNA fragments. The emergence of such genetic elements has caused a shift in the genetic location of ciprofloxacin resistance determinants from the chromosomal mutations to various mobile genetic elements. The distribution of these PMQR plasmids showed that they exhibited high serotype specificity, except for the p10k-like plasmids, which can be widely detected and efficiently transmitted among Salmonella strains of various serotypes by fusing to a new conjugative helper plasmid. The sharp increase in the prevalence of ciprofloxacin resistance in recent years may cause a predisposition to the emergence of multidrug-resistant Salmonella strains and pose huge challenges to public health and infection control efforts.
INTRODUCTION
Salmonella spp. remain one of the most important bacterial pathogens that cause foodborne diseases. Although the majority of Salmonella infections are self-limiting, they may occasionally cause systemic infections, for which the mortality rate is particularly high among immunocompromised and elderly patients. In such cases, antimicrobial treatment is required (1, 2). Three antibiotics have been approved by the FDA in the United States to treat infections caused by Salmonella: ciprofloxacin, ceftriaxone, and azithromycin (3–6). In recent years, however, resistance to these antibiotics has been increasingly common, with a particularly high rate being reported in Asia. Importantly, this rapid increase in the rate of resistance to ciprofloxacin is observable among all serotypes of Salmonella.
In the past, mutational changes in target genes, which often involved double mutations in the gyrA gene and a single mutation in the parC gene, were the primary mechanisms of resistance to ciprofloxacin (7, 8). Since 2000, plasmid-mediated quinolone resistance (PMQR) determinants, such as qnrA, qnrB, qnrC, qnrD, qnrS, aac(6′)lb-cr, and oqxAB, are being increasingly reported in Salmonella (9–14). PMQR genes normally mediate expression of quinolone resistance and reduction in susceptibility to ciprofloxacin. However, single mutation in gyrA, along with carriage of a PMQR gene, was found to mediate expression of ciprofloxacin resistance in Salmonella (15). In addition, ciprofloxacin resistance mediated by multiple PMQR genes is not uncommon (16–18). Nevertheless, these new mechanisms normally mediate a lower level of resistance to ciprofloxacin, with a MIC of >8 μg/ml being commonly reported in strains of various serotypes (14–16, 18). However, one feature of these new resistance mechanisms is that they are encoded by mobile genetic elements, thereby enabling Salmonella to acquire ciprofloxacin resistance rapidly without having to pay the fitness cost associated with mutational changes in the drug target genes. Most of the currently known PMQR-encoding plasmids are nonconjugative, limiting the rate by which they are transmitted among Salmonella strains. However, several studies have recently reported the discovery of conjugative plasmids that harbor PMQR genes in Salmonella and other Gram-negative bacteria (18). Worse still, conjugative helper plasmids that can be fused with nonconjugative PMQR-encoding plasmids and mediate their transmission among Salmonella strains have also been reported (17). However, these studies did not show how genetic elements that encode ciprofloxacin resistance evolved in Salmonella. In this work, we conducted surveillance of Salmonella strains recovered in various food samples in Shenzhen, China, during the period 2012 to 2017 and witnessed a rapid increase in the incidence of ciprofloxacin resistance in foodborne Salmonella isolates. Whole-genome sequencing and genetic analysis enable us to depict the molecular mechanisms underlying the rapid evolution of ciprofloxacin resistance-encoding genetic elements in different serotypes of Salmonella and provide valuable insights into the key genetic elements concerned, facilitating the development of new strategies to control the dissemination of resistance-encoding elements in Salmonella.
RESULTS
Prevalence of Salmonella isolates in food products.A total of 1,116 nonrepeated Salmonella strains were collected from 2,989 food samples purchased during the period of 2012 to 2017 in Shenzhen, China, as part of the food safety surveillance program organized by the Shenzhen Government. Among the 2,989 food samples, 1,558 were purchased from wet markets and 1,131 were purchased from supermarkets. These food samples contained 1,459 pork, 230 beef, 543 chicken, and 342 shrimp samples. Among the 1,116 Salmonella isolates, 82 (26%) were isolated from 317 food samples in 2013, 157 (25%) from 440 food samples in 2014, 287 (23%) from 754 food samples in 2015, 445 (23%) from 1,107 food samples in 2016, and 145 (20%) from 371 food samples in 2017 (Table 1). Among these Salmonella isolates, 746 were isolated from the 1,459 pork samples (30%), 282 were recovered from the 543 chicken samples (28%), 65 were from the 230 beef samples (17%), and 22 were from the 342 shrimp samples (5%) (Table 1).
Prevalence of Salmonella in different food products in Shenzhen, 2013 to 2017 (n = 2,689)
Antimicrobial susceptibility of Salmonella food isolates.These Salmonella isolates were found to be resistant to most of the antibiotics tested, with 51%, 13%, and 5% being resistant to ciprofloxacin, ceftriaxone, and azithromycin, respectively. These three antibiotics are the current choices for treatment for Salmonella infections. The rate of resistance to other agents is 76% to tetracycline, 60% to sulfamethoxazole-trimethoprim, 60% to ampicillin, 52% to chloramphenicol, and 42% to nalidixic acid. On the other hand, these strains were mostly susceptible to meropenem (100%), amikacin (97%), and kanamycin (75%). It is worrisome that a much higher proportion of Salmonella isolates that were resistant to ciprofloxacin, ceftriaxone, and azithromycin were also resistant to other antibiotics, unlike isolates that were susceptible to these three antibiotics. The rate of resistance to ciprofloxacin was also found to have increased significantly over the years, from 39% (33/82) in 2012 to 2013 to 77% (112/145) in 2017. Cefotaxime resistance rate, on the other hand, decreased significantly, from 17% (14/82) in 2012 to 2013 to 8% (12/157) in 2014 and then to 3% (8/290) in 2015, yet the rate increased again and reached 16% (70/445) in 2016 and further climbed to 25% (36/145) in 2017 (see Table S1 in the supplemental material).
TABLE S1
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Ciprofloxacin resistance in different serotypes of Salmonella.The 1,116 Salmonella isolates were found to belong to 54 serotypes, of which Salmonella enterica serovar Derby, Salmonella Typhimurium, and Salmonella Rissen were the most common, accounting for 24%, 17%, and 9% of all test strains, respectively. Each of several serotypes of Salmonella isolates, such as Salmonella Corvallis, Salmonella London, and Salmonella Agona, accounted for approximately 6% of the test strains. These serotypes of Salmonella were also prevalent among those that were resistant to ciprofloxacin. Although the distribution patterns of these serotypes, including Salmonella Corvallis strains, which were recovered after 2015, were highly similar over the past few years, the ciprofloxacin resistance rate of these serotypes has increased sharply, reaching the highest rate in the year 2017. Among the 269 S. Derby isolates tested, resistance to ciprofloxacin increased dramatically over the years, with a rate of 40%, 32%, 54%, 50%, and 86% being recorded in years 2012 to 2013, 2014, 2015, 2016, and 2017, respectively (Table 2). Likewise, ciprofloxacin-resistant strains accounted for 64%, 59%, 57%, 48%, and 84% of the Salmonella Typhimurium isolates collected in these 5 years (a total of 198 isolates), respectively, suggesting that the high resistance rate recorded throughout the study period increased further in 2017. For the 104 Salmonella strains that belonged to the serotype Rissen, none of the isolates recovered in 2012 to 2013 and 2014 was resistant to ciprofloxacin, yet 7%, 29%, and 57% of isolates were found to be ciprofloxacin resistant in 2015, 2016, and 2017, respectively. A similar trend of ciprofloxacin resistance was observed in other serotypes of Salmonella, with some exhibiting over 90% resistance (S. Corvallis, S. London, S. Kentucky, and S. Give) (Table 2).
Prevalence rate and mechanism of ciprofloxacin resistance in various Salmonella serotypes
Evolution of genetic elements that encode ciprofloxacin resistance in Salmonella isolates.To further investigate the genetic basis of variation in the prevalence of ciprofloxacin resistance in Salmonella during the study period, phylogenetic and genomic analyses of 816 isolates, including 250 ciprofloxacin-sensitive (Cips) and 566 ciprofloxacin-resistant (Cipr) isolates, were performed to study the pattern of distribution of PMQR genes and the types of target mutations among the resistant strains as well as the mechanism of transfer of resistance-encoding genetic elements between resistant and susceptible Salmonella strains. Ciprofloxacin resistance was serotype dependent, with some serotypes exhibiting a very high rate, whereas a low rate was recorded in the others. No significant genetic differences were observed between Cipr and Cips Salmonella strains within the same serotype (Fig. 1). Ciprofloxacin resistance was mainly mediated by PMQR genes, with 95% (539/566) of Cipr isolates being found to harbor PMQR genes; importantly, 78% (445/566) of these PMQR gene-bearing Cipr strains did not harbor any mutations in the quinolone resistance-determining regions (QRDR) of the target genes. The carriage rate of PMQR genes in Salmonella isolates increased from 20% in 2012 to 2013 to 65% in 2017. These genes were commonly detected in serotypes such as S. Derby, S. London, S. Typhimurium, S. Risen, S. Corvallis, and S. Agona (Fig. 1). Target mutations were detected in 122 of the 566 Cipr isolates and were most commonly detected in strains of the following serotypes: S. Typhimurium (74/122), S. Infantis (3/122), S. Kentucky (7/122). S. Albany (9/122), and S. Indiana (15/122). However, the rate of carriage of target gene mutations among the Cipr Salmonella strains was found to decline steadily from 39% (13/33) in 2013 to 15% (17/145) in 2017, despite a marked increase in the prevalence of ciprofloxacin-resistant strains during this period. Among the Cipr strains with target mutations, S. Typhimurium, S. Infantis, and S. Albany contained only a single mutation in the gyrA gene (S83F, S83G, S83L, S83N, S83Y, D87N, D87Y, or D87G), whereas S. Indiana and S. Kentucky harbored the double mutation S83L and D87Y in the gyrA gene and a single mutation, S80I/S80R, in the parC gene. Importantly, different PMQR genes, including qnrB6-qnrB4, qnrS1, oqxAB, and aac(6′)Ib-cr, were concurrently detected in isolates that contained target gene mutations, especially isolates collected after 2014. Overall, among the 566 Cipr Salmonella isolates, 8 isolates carried only single gyrA mutations with a CIP MIC of 1 μg/ml; 91 isolates carried single gyrA mutations and PMQR genes with CIP MICs ranging from 1 to ∼8 μg/ml; 22 isolates carried double gyrA and single parC mutations and exhibited CIP MICs of >16 μg/ml with and without PMQR genes; and 445 isolates only carried different PMQR genes, with CIP MICs ranging from 1 to ∼16 μg/ml.
Phylogenetic tree of different serotypes and sequence types among 816 Salmonella isolates collected during the period of 2012 to 2017. Circle 1 depicts strain isolation time; circle 2 denotes distribution of serotypes, using various colors shown in the data set legend (the main serotypes are also highlighted in clades of different colors); circle 3, of blue color, denotes ciprofloxacin-resistant isolates; and circle 4 and circle 5, of purple and deep red color, denote the distribution of PMQR genes and target mutations, respectively.
Detailed analysis of the genome sequences of these Salmonella isolates revealed that various combinations of PMQR genes were present and that some of the resistance gene combinations appeared to be closely associated with specific serotypes. The qnrS1-oqxAB-qnrS1 combination was the most common PMQR gene cluster and was found in 67% (375/566) of the ciprofloxacin-resistant Salmonella strains of various serotypes. The combination of qnrS2 and oqxAB–aac(6′)-Ib-cr (112/566) was the second most common but was limited to S. Derby. The qnrB6-qnrB4–aac(6′)-Ib-cr (79/566) combination was commonly detectable in S. London and S. Kentucky. The aac(6′)lb-cr–oqxAB (80/566) combination was most commonly observed in S. Typhimurium (Table 3). In summary, the rapid increase in prevalence of ciprofloxacin resistance in Salmonella might be attributed to the following mechanisms: (i) horizontal transfer of PMQR gene-bearing mobile genetic elements among Salmonella strains and (ii) carriage of a single target mutation and acquisition of exogenous PMQR genes (14).
Plasmids and TUs that harbored four major types of PMQR gene-bearing genetic elements in strains of different Salmonella serotypesa
Transmission of PMQR genetic mobile elements in ciprofloxacin-resistant Salmonella.It has been known that most of the PMQR genes in Salmonella are located in plasmids (19, 20). Therefore, the Cipr Salmonella strains were subjected to screening of plasmids and gene cassettes harboring these PMQR mobile elements. A total of nine different types of plasmids, one chromosomal fragment, and five transposable units (TUs) were found to carry PMQR genes (Table 3). Six types of plasmids, namely, p10k-like plasmid (Fig. S1), pSA1892-248kb (Fig. S2a), pSa21-28kb (Fig. S2b), pSa4-237kb (Fig. S3a), pSH01 (Fig. S3b), and pSa27-186kb (Fig. S4a), and two types of TUs, including TU_4kb and TU_13kb, were found to harbor a qnrS1-bearing mobile element (Table 3, Fig. 2). Among these plasmids and TUs, plasmid p10k (GenBank accession no. CP025337), which belonged to the ColRNA1 plasmid type, exhibited a low-level ciprofloxacin resistance in Salmonella (13, 21–23). The second most common type of mobile element was those carrying the qnrS2–aac(6′)lb-cr–oqxAB resistance gene combination (Fig. S4b), which was often found in a chromosomal DNA fragment in S. Derby, S. Typhimurim, and S. Corvallis, with S. Derby being the most dominant (108/112). The third most common type was mobile elements carrying the qnrB6-qnrB4–aac(6′)lb-cr genes. Two kinds of plasmids, pSa76-CIP (Fig. S5a) and pLA_64kb (Fig. S5b), three kinds of TUs, TU_100kb, TU_30kb (Fig. S6a), and TU_44kb (Fig. S6b), and a class I integron (Fig. S7a) were involved in the transmission of this PMQR gene-bearing element. Most of these elements were located in pSa76-CIP, a conjugative plasmid. A few of them were located in nonconjugative plasmids, such as pLA-64, and TUs with unknown transmission potential. Unlike the plasmids and TUs, class I integron carrying the qnrB4 gene was found in contigs of both 16,876 and 186,40 bp from S. Thompson, in which ISCR1 and genes encoding permeases (sapA, sapB, and sapC) and phage shock proteins (pspA, pspB, pspC, and pspD), as well as the AmpC-encoding blaDHA-1 gene, were also observed. This same region has been reported in different plasmids of various bacterial species (24). The fourth type of mobile element, carrying the aac(6′)lb-cr–oqxAB genes and located in a nonconjugative plasmid, pCFSA244-1, was mainly transmitted among S. Typhimurium isolates, with 73 isolates harboring such an element being detected (Fig. S7b). The transmission patterns of these PMQR genes in the major serotypes of ciprofloxacin-resistant Salmonella are described below (Table 3).
Structure alignment of different transposable units (TUs) carrying the qnrS1– or qnrB6–aac(6′)lb-cr genes, and comparison with plasmids with identical gene environments. (a) Alignment of TUs and plasmids carrying qnrS1 using Easyfig. Representative TUs, including TU_4kb, TU_13kb, and p10-like plasmid pSA1423-CIP, were aligned. Key genetic loci in TUs and plasmids were labeled. TU sequence was generated from Illumina sequencing data. Red denotes the resistance genes, and yellow denotes ISs. (b) Alignment of TUs and plasmid carrying qnrB6–aac(6′)lb-cr by Easyfig. Plasmid pSa76-CIP was mainly prevalent in strains of the S. London serotype. It shares antibiotic resistance gene structure identical to that of two other transposable units, which could only be detected in S. Kentucky and S. Corvallis, respectively. Red denotes the resistance genes, yellow denotes ISs, and the blue region depicts the presence of an integron, In27.
FIG S1
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
FIG S2
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
FIG S3
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
FIG S4
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
FIG S5
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
FIG S6
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
FIG S7
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Conjugative transmission of PMQR-bearing plasmids in Salmonella.Transmission of PMQR genes in Salmonella can be mediated by conjugative plasmids or conjugative helper plasmids (17, 18). Unlike the pSa1852-248kb, pSa76-CIP, and pSa4-237kb plasmids that were self-conjugative, the conjugation process of pSa21-CIP and pSa27-CIP needs to be mediated by a conjugative helper plasmid. Accordingly, the IncI1 conjugative helper plasmids pSa21-HP and pSa27-HP could be detected in most (8/64) of the Salmonella strains that carried pSa76-CIP or pSa21-CIP (Table S1). Nonconjugative plasmids such as pLA-64 and pSH01 were reported in only a few isolates, with the exception of pCFSA244-1, which harbored the aac(6′)lb-cr–oqxAB element, which was detectable in 73 isolates of S. Typhimurium (Fig. S7b). The molecular mechanism underlying the conjugation of this plasmid was not fully understood, and more work is required to investigate how these plasmids are transmitted among S. Typhimurium strains. A chromosomal MDR fragment carrying the PMQR genes qnrS2–aac(6′)lb-cr–oqxAB was mainly detected in S. Derby, although it was occasionally detectable in other serotypes, such as S. Typhimurium (n = 3) and S. Corvallis (n = 1), suggesting that this chromosomal fragment can be transmitted among strains of different serotypes. Our data confirmed that this chromosomal fragment could form a circular intermediate, which may facilitate its transmission between different Salmonella strains (Fig. S7c).
Importantly, the most predominant plasmid, p10k, was found to be widely distributed in multiple serotypes of Salmonella (Fig. S1, Table S2). This plasmid was first detected in 2014 and was found to be disseminated extensively among clinical Salmonella strains collected in 2015 to 2016. This phenomenon, which suggests that the transmission of p10k-like plasmids in Salmonella was not due to clonal transmission, is not consistent with the nonconjugative nature of this plasmid. To better understand the transmission mechanisms of this nonconjugative ColRNA1 plasmid, 50 randomly selected Salmonella strains carrying this plasmid were subjected to conjugation experiments using Escherichia coli J53 as the recipient strain. A total of 10 isolates could transfer the p10k-like plasmid to J53, suggesting that some currently unknown genetic mechanisms mediate the conjugation of p10k-like plasmid. One Salmonella strain, SA1423, and its transconjugant, SA1423-TC, were selected for complete plasmid sequencing by using both Illumina and Nanopore platforms. Our data showed that the p10k-like plasmid in SA1423 was a 10,218-bp plasmid that contained the qnrS1 gene, with 50.6% GC content. It exhibited 99% identity at 99% coverage with plasmid p10k and was designated pSa1423-CIP. pSa1423-CIP was able to form a hybrid plasmid with a 49,279-bp conjugative plasmid and was designated pSa1423-HP, and it exhibited 71% identity at 99% coverage with the plasmid pVQS1 (GenBank accession no. JQ609357.1), carried by a Salmonella strain originating from a traveler in Germany, through homologous recombination at the ISKpn19 site (Fig. 3). Formation of the hybrid plasmid pSa1423-CIP-HP facilitated the transmission of pSa1423-CIP among different Salmonella isolates. Screening of pSa1423-HP in all Salmonella isolates identified 18 isolates (Table S3) that carried this helper plasmid, further confirming that this helper plasmid was involved in the transmission of p10k-like plasmids in Salmonella isolates.
Mechanisms of plasmid recombination. (A) Structure alignment of three plasmids. Antimicrobial resistance genes, including qnrS1 and tet(A) genes, are highlighted in red; yellow depicts IS elements, including IS26, IS2, and ISKpn19; and the tra genes responsible for conjugation are shown by a green arrow. (B) Mechanism of plasmid fusion. Proposed IS element-mediated plasmid fusion through homologous recombination in Salmonella Sa1423. pSa1423-HP acted as the target plasmid, which attacked the hot spot in the donor plasmid, pSa1423-CIP, to facilitate plasmid fusion.
TABLE S2
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TABLE S3
Copyright © 2020 Chen et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
DISCUSSION
The genetic features of foodborne Salmonella strains have changed significantly in recent years as ciprofloxacin-resistant Salmonella strains have become prevalent. Strains of some serotypes, such as S. Corvallis, S. Kentucky, and S. Stanley, that were highly resistant to ciprofloxacin have emerged. We observed that the mechanisms of ciprofloxacin resistance delineated in Salmonella strains have evolved significantly during this period. A major observation is that target gene mutations have become much less common in ciprofloxacin-resistant Salmonella strains. Currently, target gene mutations, in particular single mutation in the gyrA gene, may be detected in S. Typhimurium but is rarely seen in other serotypes of Salmonella. Double mutations in gyrA and single mutation in parC that mediated high-level resistance were seen only in S. Indiana, S. Kentucky, and S. Typhimurium strains but not in other serotypes (7, 25–27). Instead, PMQR genes, often in the form of multiple PMQR genes located in the chromosome and plasmids, now are commonly detectable in ciprofloxacin-resistant Salmonella strains of various serotypes. This shift in genetic location of ciprofloxacin resistance determinants from the chromosome to plasmids is responsible for the sharp increase in the prevalence of ciprofloxacin-resistant clinical Salmonella strains.
Mobile genetic elements carrying multiple PMQR genes are increasingly being detected in clinical Salmonella isolates and are responsible for the transmission of these ciprofloxacin resistance-encoding genes among such strains (16–18). The qnrS1 gene was found to be harbored by multiple plasmids, with the p10k-like plasmids being the most dominant among a wide range of serotypes. The p10k-like plasmids were not conjugative but have become transmissible with the help of a helper plasmid, which could be lost upon being successfully transmitted to a new Salmonella strain. Our data showed that the widespread dissemination of p10k-like plasmids in Salmonella was closely associated with the emergence of the helper plasmid. Before 2015, only two Salmonella strains carrying p10k-like plasmids were detected, and no helper plasmid was found. In 2015, with the detection of pSA1423-HP, transmission of p10k in Salmonella accelerated, with as many as 119 isolates being found to carry this plasmid. The rate of detection of this plasmid continued to increase in 2016, with 133 isolates being found to carry this plasmid in that year. Meanwhile, the helper plasmid became detectable in a large number of isolates. In 2017, however, the helper plasmid was not detected in any strain, and the number of p10k-like plasmids in Salmonella was also found to have reduced dramatically. Other plasmids and TUs that carry the qnrS1 gene have become dominant. This is an excellent example of dynamic evolution of ciprofloxacin resistance mechanisms in Salmonella, in which one dominant mechanism is not only gradually being replaced by another but also the transferrable nature of the new resistance mechanism resulted in a spike of the ciprofloxacin resistance rate of Salmonella to as high as 77% in 2017. In various serotypes of Salmonella, the qnrS1 gene product was found to act together with other mechanisms to mediate the expression of the ciprofloxacin resistance phenotype. For example, in S. Typhimurium, it is often paired with a single gyrA mutation to cause ciprofloxacin resistance; in other serotypes of Salmonella, a combination of qnrS1 and other PMQR genes located in chromosomal fragments or plasmids is often observed. In addition to being harbored by the p10k-like plasmid, the qnrS1 gene can also be found in various conjugative and nonconjugative plasmids that are often serotype specific.
Another major mechanism of ciprofloxacin resistance involves the chromosomal fragment harboring the qnrS2–aac(6′)lb-cr–oqxAB elements. This chromosomal fragment, which was mainly detected in S. Derby, could produce a CIP MIC of ≥8 μg/ml. S. Derby isolates carrying this fragment were increasingly detected in recent years as a result of clonal spread. Our work showed that this chromosomal DNA fragment, flanked by IS26, was able to form a circular intermediate that facilitates the transmission of this DNA fragment between different Salmonella strains. The qnrB6-qnrB4 and aac(6′)lb-cr genes were commonly detected in different plasmids and TUs that are uniquely present in different serotypes of Salmonella. One limitation of this study is that we are not able to depict the genetic structure of some TU-bearing plasmids and the mechanisms underlying the transmission of the pCFSA244-1 plasmid that contains the aac(6′)lb-cr–oqxAB elements in 80 isolates of S. Typhimurium. Since these isolates of S. Typhimurium are not derived from the same clone, clonal transmission of this plasmid should not be the sole mechanism of transmission of this ciprofloxacin resistance-encoding genetic fragment. Further studies are required to elucidate the mechanism of transmission of this nonconjugative plasmid.
Data from this study and others have suggested the underlying mechanisms of the rapid development of fluoroquinolone resistance in Salmonella spp. Salmonella strains used to be susceptible to fluoroquinolone antibiotics due to its low rate of development of mutations in their target genes, namely, double gryA mutations and single parC mutation, which cause high-level fluoroquinolone resistance (7, 8). Since their emergence around 2010, the PMQR genes have been found in the Salmonella chromosome or nonconjugative plasmid. The products of such genes produce low-level fluoroquinolone resistance when acting in combination with a single gyrA mutation, speeding up the development of fluoroquinolone resistance in Salmonella (15). The number and types of PMQR genes have increased rapidly in Salmonella due to transmission of these PMQR genes by conjugative plasmids and helper plasmids that could help conjugate a nonconjugative PMQR gene-coding plasmid to other Salmonella strains (16–18). The efficient transmission of PMQR genes in Salmonella enables one Salmonella strain to acquire multiple types of PMQR genes that could mediate the expression of low-level fluoroquinolone resistance in Salmonella without mutation in the target genes. The efficient transmission of PMQR genes resulted in rapid development of fluoroquinolone resistance in Salmonella in recent years. In conclusion, this study has provided comprehensive insight into the rapid evolution of ciprofloxacin resistance in Salmonella in the past 6 years. Continuous surveillance is warranted to depict the future trend of development of ciprofloxacin resistance in Salmonella.
MATERIALS AND METHODS
Bacterial isolates.A total of 1,116 nonduplicate Salmonella isolates were collected from food samples during the period from 2012 to 2017 in Shenzhen, China. The species identities of these isolates were confirmed by detection of the Salmonella-specific gene invA and by using the matrix-assisted laser desorption ionization–time of flight mass spectrometry biotyper system (Bruker, Germany). Serotyping of the isolates then was performed according to the Kauffmann-White scheme by using a commercial antiserum (Difco, Detroit, MI), followed by testing of susceptibility to antimicrobial drugs according to CLSI guidelines (28). Ciprofloxacin-resistant isolates then were subjected to molecular screening for qnrA, qnrB, qnrC, qnrD, qnrS, aac(6′)lb-cr, and oqxAB as previously described (15). A total of 566 ciprofloxacin-resistant (Cipr) and 250 ciprofloxacin-sensitive (Cips) isolates were selected for further investigation.
Filter mating assay.The transmission potential of PMQR genes was assessed by performing conjugation experiments, in which the filter-mating method was used as previously described (29). Transconjugants were selected on eosin methylene blue agar containing ciprofloxacin (0.5 mg/liter) and sodium azide (100 mg/liter). The transconjugants and their parental strains both were tested for their susceptibility to ciprofloxacin.
Plasmid sequencing and analysis.Conjugative plasmids were collected from the test strains and their corresponding transconjugants using the Qiagen Plasmid Midi kit (Qiagen, Valencia, CA). The Illumina platform and Nanopore MinION long-read sequencing platform were used to draft whole-plasmid maps. The Illumina paired-end libraries were prepared by using the NEBNext Ultra DNA library prep kit for Illumina (NEB) and then sequenced on an Illumina NextSeq 500 platform. De novo assemblies of MinION long reads and Illumina reads were performed with SPAdes 3.12.1 (30) and the CLC Genomics Workbench (CLC bio, Denmark), respectively. Long assembled contigs obtained from MinION long reads were used to align and join the contigs obtained from the Illumina assembly results. The completed plasmid sequence was annotated by the RAST tool (31) and the NCBI Prokaryotic Genome Annotation Pipeline (PGAP). All plasmids were sequenced using both Illumina and MinION long reads, and only high-quality data were used for further analysis. Screening and alignment of plasmids with similar structures were produced by BLAST Ring Image Generator (BRIG), version 0.95.22. (32).
Whole-genome sequencing and genomic analysis.A Pure-Link genomic DNA minikit (Invitrogen, USA) was used to extract DNA. DNA libraries were constructed by using the Nextera XT DNA sample preparation kit (Illumina, San Diego, CA, USA). Samples were multiplexed and sequenced on an Illumina Hiseq X for 300 cycles (250-bp paired end). Raw reads were trimmed and quality filtered using Trimmomatic v0.36 (33). Draft genomes were acquired by using SPAdes version 3.10.1 (34). Species identity and serotypes of the test strains were confirmed by SISTR (35). Whole-genome phylogenetic trees were created containing reference isolates for identification of S. Derby 14-Sa79. To detect PMQR genes and assess the distribution patterns of such genes in plasmids, the draft genome searches were conducted by BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi), ResFinder (36), PlasmidFinder (37), and the CLC Genomics Workbench (CLC bio, Denmark) (Table 2). Integrons identified in the genomes were categorized according to INTEGRALL (http://integrall.bio.ua.pt).
Phylogenetic analysis.Trimmed and quality-filtered Illumina sequencing reads obtained from Salmonella strains were mapped to a ciprofloxacin-resistant S. Derby strain, 14-Sa79, which was sequenced using Nanopore and was confirmed to carry the qnrS2–aac(6′)lb-cr–oqxAB genes in its chromosome. Single-nucleotide polymorphisms (SNPs) were called using Snippy v3.1 with default settings (38), which used BWA-MEM v0.7.12 for short-read mapping. Snippy generates a core SNP alignment as well as a whole-genome SNP alignment. The whole-genome alignment was used to infer the ML phylogenies using Fasttree v2.1.10 with default parameters (39). The phylogenetic tree was visualized by iTOL version 3 (40). Eight hundred sixteen of the 1,116 Salmonella isolates, including 566 Cipr and 250 Cips isolates, were included in this study. SISTR analysis, which identified 47 serotypes from 816 genomes, including 34 serotypes from 566 strains that exhibited the Cipr phenotype, was also performed in the clustering analysis (41). Six serotypes of Salmonella that were predominantly prevalent among the foodborne strains and exhibited resistance to ciprofloxacin were also subjected to maximum likelihood phylogeny analysis.
Data availability.The sequencing data of the chromosome have been deposited in GenBank under BioProject ID PRJNA682289, and the sequences of the plasmids were deposited under the following accession numbers: pSA1852-248kb (MT513102), pSA1423-CIP (MK356559), pSA1423-HP (MK356560), pSA1423-TC-CIP (MK356561), and Sa79-chr (SGWG00000000).
ACKNOWLEDGMENTS
We gratefully appreciate financial support of the project by the National Key Research and Development Plan (2018YFD0500300), Basic Research Fund of Shenzhen City (20170410160041091), and Collaborative Research Fund of Hong Kong Research Grant Council (C5026-16G).
We have no conflicts of interest to declare.
K.C.C. and Y.C. performed the experiments; D.N., M.M.X., and L.W.Y. helped with sequencing and data analysis; E.W.C.C. and S.C. contributed to study design and manuscript writing; S.C. supervised the whole project.
FOOTNOTES
- Received November 30, 2020.
- Accepted December 7, 2020.
- Copyright © 2020 Chen et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.