Contribution of Different Mechanisms to Ciprofloxacin Resistance in Salmonella spp.

Development of fluoroquinolone resistance can involve several mechanisms that include chromosomal mutations in genes (gyrAB and parCE) encoding the target bacterial topoisomerase enzymes, increased expression of the AcrAB-TolC efflux system, and acquisition of transmissible quinolone-resistance genes. In this study, 176 Salmonella isolates from animals with a broad range of ciprofloxacin MICs were collected to analyze the contribution of these different mechanisms to different phenotypes. All isolates were classified according to their ciprofloxacin susceptibility pattern into five groups as follows: highly resistant (HR), resistant (R), intermediate (I), reduced susceptibility (RS), and susceptible (S). We found that the ParC T57S substitution was common in strains exhibiting lowest MICs of ciprofloxacin while increased MICs depended on the type of GyrA mutation. The ParC T57S substitution appeared to incur little cost to bacterial fitness on its own. The presence of PMQR genes represented an route for resistance development in the absence of target-site mutations. Switching of the plasmid-mediated quinolone resistance (PMQR) gene location from a plasmid to the chromosome was observed and resulted in decreased ciprofloxacin susceptibility; this also correlated with increased fitness and a stable resistance phenotype. The overexpression of AcrAB-TolC played an important role in isolates with small decreases in susceptibility and expression was upregulated by MarA more often than by RamA. This study increases our understanding of the relative importance of several resistance mechanisms in the development of fluoroquinolone resistance in Salmonella from the food chain.

Introduction

Salmonella species are important foodborne pathogens and associated infections can be life-threatening, particularly in elderly and immunocompromised patients (Majowicz et al., 2010). Fluoroquinolones (FQ) are the primary treatment option for life-threatening Salmonella infections in clinical practice (Hohmann, 2001). However, the emergence of FQ resistance in Salmonella spp. has limited therapeutic options and the World Health Organization categorized FQ-resistant Salmonella as a high priority pathogen for the research and development of new antibiotics in 2017 (WHO, 2017).

Mechanisms of FQ resistance include chromosomal mutations in the target enzymes, DNA gyrase and topoisomerase IV, down-regulation of outer membrane porin expression coupled with increases in active drug efflux, as well as acquisition of transmissible quinolone-resistance genes. The development of FQ resistance is often mediated by the accumulation of multiple mutations in a stepwise process and the interplay between multiple resistance mechanisms can result in the development of highly resistant mutants (Redgrave et al., 2014; Hooper and Jacoby, 2015). It has been suggested that, mutation of chromosomally located primary target genes represents the first step in the development of quinolone resistance. For Gram-negative organisms such as Escherichia coli, the primary target is typically the GyrA subunit of gyrase. The most common mutation site in the quinolone resistance determining region (QRDR) is serine 83 of GyrA. Highly resistant organisms typically carry a combination of mutations within gyrA and parC (Machuca et al., 2014; Redgrave et al., 2014; Hooper and Jacoby, 2015). However, whilst similar mutations are often recovered in Salmonella spp. as in E. coli after quinolone exposure (Marcusson et al., 2009; Machuca et al., 2014), the phenotypic impact resulting from target site mutations seem to have a smaller impact (conferring smaller MIC increases) than in E. coli (Heisig, 1996; Lin et al., 2015).

The chromosomal multidrug efflux pump AcrAB-TolC is capable of actively removing FQs and other drugs from the bacterial cell and is known to play an important role in the development of high level FQ resistance (Redgrave et al., 2014),overexpression of acrAB has been suggested to be a first step that facilitates high-level resistance development following acquisition of target site mutations (Giraud et al., 2000; Singh et al., 2012).

Transmissible quinolone-resistance mechanisms are often plasmid associated and so are known as “plasmid-mediated quinolone resistance” (PMQR) determinants. These typically confer a low-level of decreased susceptibility to quinolones. As with efflux over-expression, carriage of PMQRs can also facilitate the selection of mutants with higher levels of quinolone resistance through additional chromosomally encoded mechanisms (Hooper and Jacoby, 2015).

There is evidence for clinical quinolone resistance to emerge without target site mutation, a fully resistant E. coli strain with a ciprofloxacin MIC of 4 mg/L but without topoisomerase mutations has been reported and was shown to harbor plasmid-mediated qnrS1 and oqxAB, as well as overexpressing acrAB and genes encoding for other efflux pumps (Sato et al., 2013). Similarly, an experimentally derived mutant carrying five plasmid copies of qnrA1 elevated the ciprofloxacin MIC from 0.25 to 2 mg/L; this MIC value exceeds the CLSI breakpoint for resistance (Vinue et al., 2019). In a separate study, several PMQR genes either alone or in combination were reported to mediate ciprofloxacin resistance development in Salmonella isolates that did not contain target gene mutations (Lin et al., 2015).

The acquisition of antibiotic resistance in bacteria via either chromosomal or plasmid mechanisms is often accompanied by a fitness cost. But this is not always the case in quinolone resistance development. A study (Marcusson et al., 2009) investigating the fitness effect of various QRDR mutations and efflux activity in E. coli revealed that GyrA 83 mutations conferred a fitness gain on the isolates while GyrA 87 substitutions were associated with fitness costs. Acquisition of a triplet of GyrA 83, GyrA 87, and ParC 80 substitutions resulted in fitness gain when efflux was not involved. De-repressed efflux activity was, however, associated with fitness costs. Resistant mutants harboring the same triplet mutations, which was obtained by exposure of a quinolone-susceptible Salmonella Typhimurium clinical isolate to increasing concentrations of ciprofloxacin also showed a diminished ability to grow (Fabrega et al., 2014). Huseby et al. (2017) demonstrated that the even low-levels of ciprofloxacin strongly favor growth of topoisomerase mutants and even genotypes (gyrA 83 and parC 80) which carry a cost in drug free media are favored in sub-inhibitory concentrations of quinolone.

Most of our understanding of the contributions and fitness costs of different mechanisms to quinolone resistance is from laboratory studies or analysis of small numbers of mutants. Given the importance of quinolone resistance as a phenotype, a thorough understanding of the impact of the different FQ resistance mechanisms in the real world is of crucial importance. In the present study, 176 Salmonella isolates with a broad range of ciprofloxacin MIC distributions were collected. These isolates were analyzed for correlations between ciprofloxacin resistance phenotypes and: target site mutations, expression levels of acrAB, and prevalence of PMQR genes. This research assessed the relative contributions of defined mutant combinations toward fitness in different FQ concentrations and informs how FQ resistance evolves in the food chain.

Materials and Methods

Bacterial Strains and Susceptibility Testing

Salmonella isolates (126) used in this study were isolated from pork samples in a large-scale slaughterhouse between 2013 and 2014 (Yang et al., 2017). The ciprofloxacin MICs ranged from <0.015 to 64 mg/L although none of these isolates were inhibited by ciprofloxacin concentrations between 0.125 and 2 mg/L. In order to generate a consecutive ciprofloxacin MIC distribution range for each MIC value, we selected 50 additional strains from food-producing animals and raw meat (chicken and pork) as described elsewhere (Zhang et al., 2016; Zhang C.Z. et al., 2019; Yang et al., 2017). This included 11 isolates from swine, 9 from chickens, 17 from ducks and 13 from meat samples (10 chicken and 3 pork) (Table 1). Although efforts were made to screen for strains with a ciprofloxacin MIC between 0.25 and 2 mg/L, these strains were rare. We only identified 3, 3, and 4 isolates in the entire strain collection that possessed ciprofloxacin MICs of 0.5, 1, and 2 mg/L, respectively.

TABLE 1

Table 1. Sources and ciprofloxacin MIC values for Salmonella isolates used in this study.

The majority of the strains isolated from pork samples were serovars Derby, Rissen and Indiana as previously described (Yang et al., 2017). The 50 additional isolates included Typhimurium (15), Indiana (12), Enteritidis (7), Derby (5), and 11 others.

Salmonella enterica subsp. enterica serovar Typhimurium SL1344 (MIC_CIP = 0.015 mg/L) was kindly provided by Prof. Laura Piddock (Institute of Microbiology and Infection, University of Birmingham) and used as a standard reference strain in this study. The MICs of ciprofloxacin were determined by the agar dilution method following CLSI guidelines (CLSI, 2019).

Detection of Target Gene Mutations and the Presence of PMQR Genes

Identification of mutations in the quinolone resistance determining regions (QRDR) of gyrA, gyrB, parC, and parE and identification of qnr, aac(6′)-Ib-cr, qepA and oqxABR genes were performed using PCR amplification, primers used are listed in Supplementary Table 1. All PCR amplicons were sequenced by the BGI (Shenzhen, China), and the wild-type Salmonella Typhimurium LT2 strain was used as a comparison to identify mutations.

Detection of the oqxABR Circular Intermediates

All oqxABR-positive Salmonella isolates were screened for the presence of the IS26 insertion sequence (IS), and cyclization of the IS26-oqxABR region was determined using primers oqx-IF and oqx-IR. Sequence identity of all PCR amplicons was confirmed by DNA sequence analysis.

Location of PMQR Genes

The location of PMQR genes and the sizes of PMQR gene-positive plasmids were estimated by S1-PFGE and/or I-CeuI-PFGE followed by Southern blot hybridization with digoxigenin-dUTP labeled probes of oqxAB, qnrS, and/or 23S rDNA genes. Primers used are listed in Supplementary Table 1. Briefly, agarose-embedded DNA for each Salmonella strain was digested with 1 U S1 nuclease (Takara Biotechnology, Dalian, China) at 37°C for 45 min or I-CeuI (NEB, Ipswich, MA, United States) for 3 h. The restriction fragments were separated by electrophoresis in 0.5 × Tris-borate-EDTA buffer at 14°C for 19 h with pulse times of 2.2 to 54.2 s using a CHEF-MAPPER System (Bio-Rad Laboratories, Hercules, CA, United States). Agarose-embedded DNA from strain H9812 that had been in situ digested with XbaI (Takara Biotechnology, Dalian, China) at 37°C for 3 h was used as a DNA size marker. The gel was stained with ethidium bromide and DNA bands were visualized under UV light. DNA was transferred to a Hybond-N + membrane for Southern hybridization using conditions suggested by the manufacturer (GE Healthcare, Little Chalfont, United Kingdom). Blots were probed with PMQR genes containing digoxigenin labels using a commercial DIG High Prime DNA Labeling and Detection Starter Kit I (Roche Applied Science, Mannheim, Germany).

Gene Expression Analysis by qRT-PCR

Bacterial strains were grown in minimal media to an OD₆₀₀ nm of 0.6 and total RNA was isolated using the RNAiso Plus Kit (Takara); mRNAs were reverse transcribed using PrimeScript RT with gDNA Eraser (Takara). The cDNA was quantified using an iQ5 multicolor real-time PCR system (Bio-Rad) using gene-specific primers (Supplementary Table 1), and an iQ SYBR Green Supermix (Bio-Rad). The 16S rRNA gene was used as an internal control gene and relative expression levels of each gene were calculated using the ΔΔCt method and the software provided with the instrument. Data are presented as mean ± SD from three independent assays, in which each RNA sample was tested in triplicate.

Competitive Fitness Measurements

To directly measure fitness impacts of different resistance mutations we used competition assays between pairs of strains selected to carry different combinations of resistance mechanisms, and with different susceptibility levels to ciprofloxacin. The selected pairs competed against each other were; S(−)/RS(−), RS(c +)/I(p +), I(p +)/R(c +), and I(p +)/HR(p +). In addition, to eradicate any possible strain and serovar specific bias, competition assays were also completed between selected serovar Indiana strains with different resistance mechanisms, RS(−)/HR(p +) (strains HB137/K46) and I(p +)/HR(p +) (strains SP80/K46, and SP80/CL108).

Fitness was analyzed as previously described with some modifications (Zhang C.Z. et al., 2017). Growth competition was determined by pyrosequencing the single nucleotide variations in ParC which distinguished each pair of competed strains. Briefly, the competitive paired strains were mix and co-cultured (1:1 ratio) in antibiotic-free LB broth for 18 h and DNA was extracted using a standard phenol-chloroform-heat method. The DNA from the competitive growth assays was amplified by PCR in triplicate using biotinylated primer pairs targeting the region containing the single nucleotide polymorphism to distinguish the two organisms in the assay by mutations in ParC. All PCR amplifications were visualized on 1% agarose gels prior to pyrosequencing. The purified PCR products were pyrosequenced at BGI. Competition coefficients were defined as the ratio of the strain pairs subjected to different incubation conditions and calculated by measuring percentage yield of the single nucleotide mutations in parC. A competition coefficient <1 indicated that the less susceptible strain was more abundant after competition.

Results

Assembly and Definition of Study Panel of Isolates

A primary aim of our study was to identify Salmonella isolates with a broad range of ciprofloxacin MICs. We selected 176 Salmonella strains from 1,280 isolates, amongst which 91 were susceptible to ciprofloxacin, 20 were deemed intermediate and 65 were resistant according to CLSI breakpoints. However, epidemiologically there were clearly more than three levels of susceptibility in the populations, likely to reflect distinct combinations of resistance mechanisms. Therefore, we applied a 5-level interpretation system using “epidemiological breakpoints.” This system relies on overlapping normal distributions of MIC values (Kronvall, 2010; Zayed et al., 2015). The resulting categories were designated as: highly resistant (HR), resistant (R), intermediate (I), reduced susceptibility (RS), and susceptible (S) (Figure 1). We then randomly selected a total of 32 strains, 5∼7 strains from each group, for detailed study of gene expression and carriage of specific mutations (Table 2 and Supplementary Figure 1).

FIGURE 1

Figure 1. Ciprofloxacin MIC distributions in 176 Salmonella strains used in this study.

TABLE 2

Table 2. Characteristics of Salmonella enterica isolates from food-producing animals in China.

Correlation of QRDR Target Gene Mutations With MIC Category

Among the isolates that were analyzed, none contained QRDR mutations in gyrB and parE. In contrast, parC mutations that resulted in the ParC T57S substitution were present in S and RS strains with an MIC_CIP of 0.06 mg/L. ParC T57S was not seen in RS strains when the MIC_CIP increased to 0.25 and 0.5 mg/L. Isolates with these MICs carried single GyrA S83Y or D87Y substitutions. The “I” strains tested only possessed GyrA D87N single substitutions; interestingly, the ParC T57S substitution was again seen in some R strains, but these lacked GyrA mutations (although had other mechanisms as outlined below). The HR strains all possessed multiple substitutions in both GyrA (S83F/D87G) and ParC (T57S/S80R) (Table 2). The identification of these mutations was positively correlated with their assignment into the 5-level interpretation system.

The Presence and Location of PMQR Genes

Reduced susceptibility strains (ParC T57S) that possessed MIC_CIP values of 0.06 mg/L also harbored chromosomal copies of oqxABR. PMQR genes were not present in the RS strains that contained the single GyrA substitutions S83Y or D87Y. When the MIC_CIP levels reached 1–2 mg/L (I strains), the presence of plasmid-borne copies of both oqxAB and aac(6′)-Ib-cr were detected on ∼180-kb plasmids. In the R strains, oqxAB, aac(6′)-Ib-cr and qnrS were detected and all were located on the chromosome. In HR strains, aac(6′)-Ib-cr was found on either ∼240-kb plasmids or was present with oqxAB and qnrS2 on the chromosome (Table 2).

Expression of AcrAB-TolC and Regulators

We next analyzed whether increases in efflux pump expression correlated with the MIC categories of our isolates. The expression of efflux complex (acrA, acrB, and tolC) and regulatory (ramA, marA, and soxS) genes for the subset of 32 Salmonella isolates were measured (Table 2 and Supplementary Figure 2). Compared with S. Typhimurium SL1344, the RS strains exhibited higher expression of acrAB-tolC and the regulator marA but not ramA. Furthermore, the majority of I strains exhibited increased expression of efflux pump genes with concurrent overexpression of marA. Surprisingly, the efflux pump genes were not expressed more in R strains compared with S. Typhimurium SL1344, whereas there was increased expression of acrA, acrB, ramA, and marA in HR strains compared with S. Typhimurium SL1344. A two-way ANOVA test of expression values for the regulator genes indicated the differences between groups were not likely to be observed by chance (p 0.007).

oqxABR Circular Intermediates

All 16 oqxABR-positive strains possessed IS26 sequences flanking the oqxABR operon (Figure 2A). In addition, all these strains, regardless of whether the genes were located on a plasmid or chromosome, were positive for the presence of circular intermediates containing IS26 and the oqxABR genes (Figures 2B,C). Sequencing of this circular intermediate confirmed that the length of the intermediate was 6,026 bp and that oqxABR harbored a complete copy of IS26. This intermediate structure was like the oqxABR operon present in the widespread composite transposon Tn6010 found in plasmid pOLA52 from E. coli (Norman et al., 2008). These results indicate that the Tn6010 element was unstable and prone to excision in our Salmonella isolates.

FIGURE 2

Figure 2. Formation of a circular intermediate by oqxABR operon. (A) Comparison between genetic structure of transposon Tn6010 located in pOLA52 (accession no. NC_010378) and the oqxABR operon positive strains. (B) The circular form of Tn6010 and the approximate locations of the reverse primers oqx-IF and oqx-IR (line arrows). (C) Gel electrophoresis of PCR amplicons corresponding to the circular intermediate of Tn6010, detectable in oqxABR-positive strains using the reverse primers oqx-IF and oqx-IR. Primers are listed in Supplementary Table 1. NC, Negative Control.

Competitive Fitness

We compared the competitive fitness of closely related pairs of isolates selected to represent the different combinations of resistance mechanisms and susceptibility levels observed in the whole panel. This included the presence of target gene mutations and presence of PMQR genes on the chromosome (c +) or on plasmids (p +). The competitive growth experiments included the pairs S(−)/RS(−), RS(c +)/I(p +), I(p +)/R(c +), and I(p +)/HR(p +) and Indiana strains RS(−)/HR(p +) (strain HB137/K46) and I(p +)/HR(p +) (strain SP80/K46, and SP80/CL108). All strains used in competition assays are shown in Table 3.

TABLE 3

Table 3. Characteristics of strains used in growth competition assays.

The competition coefficient for the S (−)/RS (−) pair was 1.13, suggesting that the presence of a substitution of GyrA D87Y presented a burden to strain HB137 (RS), while ParC T57S did not impair strain M36 (S). The competition coefficient for RS(c +)/I(p +) counterparts was 1.53, indicating that carriage of both GyrA D87N and the two plasmid PMQR genes produced a burden on strain SP79 (I). The I(p +)/R(c +) competition yielded a competition coefficient of 0.76, indicating the three chromosomal PMQR genes in strain F64 (R) produced less of a burden than when these genes were plasmid-borne. The HR strain K46 carried 2 mutations in GyrA and 2 in ParC. The presence of the plasmid-borne PMQR gene aac(6′)-Ib-cr resulted in a competitive disadvantage when compared with the I strain SP79 that possessed a GyrA single mutation; this resulted in a competition coefficient of 1.95 for I(p +)/HR(p +) (Figure 3). The same patterns were seen from the corresponding experiments using identical Indianan strains [I(p +)/HR(p +)], SP80/K46, and SP80/CL108. However, RS strain HB137, harboring a substitution of GyrA D87Y and overexpressing AcrAB, presented a competitive disadvantage when compared with HR strain K46.

FIGURE 3

Figure 3. In vitro bacterial competition assays. Competition coefficient values were obtained from each independent experiment as indicated. Broken horizontal line, competition coefficient of 1; mean values, short continuous horizontal lines. –, not detected; C+, chromosomally located; P+, plasmid located.

Discussion

The evolution of antimicrobial resistance in bacteria is driven by the pressures exerted following exposure to antimicrobials. However, the impacts of different resistance mechanisms are still not well understood. We analyzed isolates of different Salmonella serovars with a broad range of ciprofloxacin MIC distributions in this study. From studying this curated panel of isolates, the current study resulted in the following conclusions:

ParC T57S Was Present in Isolates With Low Ciprofloxacin MICs

The development of FQ resistance is understood to often require accumulation of multiple mutations in a stepwise process, with mutations that alter GyrA at codon 83 often being considered the “first step” to selection of high-level resistance. Here, we observed a high prevalence of the ParC T57S substitution in strains with lowest MICs of ciprofloxacin (S and RS strains). Whether the ParC T57S substitution contributes to quinolone resistance, is still controversial (Baucheron et al., 2002; Ling et al., 2003; Eaves et al., 2004; Gunell et al., 2009; Zhang W.H. et al., 2017). We found that the presence of ParC T57S was altered as the ciprofloxacin MIC increased; this mutation only appeared alone in strains with MIC_CIP ≤ 0.06 mg/L. Therefore, we suggest that any impact of this mutation alone of quinolone susceptibility must be low and can only provide a low degree of protection against ciprofloxacin. This change alone may carry little fitness cost which would explain its prevalence in the S and RS strains, however, it may incur a different impact on fitness when associated with GyrA substitutions as no strain carried a single substitution within ParC and any GyrA change in any of the groups. Most resistant strains did carry multiple substitutions within both GyrA and ParC and the specific nature of these changes is likely to impact fitness.

In the intermediate isolates where the MIC_CIP was higher (I; 4–8 mg/L), strains carried a QRDR mutation, such as S83Y or D87Y/N but the ParC T57S mutation was not present (suggesting this is not a precursor required for development of resistance). ParC T57S was seen in R and HR strains, along with multiple substitutions in both GyrA and ParC, this demonstrates the ParC T57S substitution is viable when in combination with multiple other topoisomerase mutations.

Expression of acrAB-tolC was elevated in strains carrying the ParC T57S substitution, but only those where the MIC_CIP was below the resistance breakpoint. This again suggests that a combination of the ParC T57S change with increased active efflux activity may be important in the early stages of resistance development (where MICs are relatively low). The growth competition assays revealed that the ParC T57S mutation was associated with enhanced bacterial fitness again supporting a low-impact but high-fitness role for this substitution.

PMQR Genes Were Associated With Ciprofloxacin Resistance Development

Plasmid-mediated quinolone resistance determinants cause decreased susceptibility to FQs while also facilitating the selection of higher levels of FQ resistance (Rodriguez-Martinez et al., 2016). In the present study, oqxAB was the most prevalent PMQR gene detected. It was not present in S strains but was first observed in RS strains with an MIC_CIP 0.06 mg/L and the gene was chromosomally located in all the positive RS strains. The occurrence of oqxAB in combination with existing mechanisms, such as ParC T57S and elevated expression of acrAB-tolC, were linked to significant increases in MIC_CIP, compared to strains without oqxAB but with these other changes. The chromosomal location of oqxAB did not appear to incur a fitness cost in strains which also increased acrAB-tolC expression and may represent a successful combination of these two mechanisms in terms of susceptibility and fitness. Isolates with an MIC_CIP of 0.25–0.5 mg/L, also carried mutations in gyrA, which imposed a fitness cost.

Some isolates where the MIC_CIP achieved clinical resistance levels (although without additional mutations in the target genes) carried >2 PMQR genes (oqxAB, aac(6′)-Ib-cr, and qnrS) combined with mutation in gyrA. These genes were found in plasmid locations in intermediate strains (MIC_CIP 1–2 mg/L) while being observed at chromosomal locations in R strains (MIC_CIP 4–8 mg/L). Interestingly, the R group did not exhibit any known chromosomal mechanisms, such as target gene mutations or overexpression of acrAB-tolC for increased resistance, but did carry three plasmid-borne PMQR genes. This phenomenon is consistent with previous observations where the presence of 2–3 PMQR genes was observed in ciprofloxacin resistant Salmonella isolates from pork samples that lacked any target gene mutations (Lin et al., 2015). Our results confirmed that the presence of multiple PMQR genes or multiple copies of PMQR genes can confer resistance to ciprofloxacin.

The Switch of PMQR Gene Location Between Chromosome and Plasmid Mediated by IS26-oqxAB Correlated With Fitness and Stable Development of Resistance

OqxAB is an efflux pump from the RND family that mediates resistance to olaquindox, chloramphenicol and nalidixic acid. The presence of this pump elevates MICs for other antimicrobial reagents including ampicillin and gentamicin (Hansen et al., 2007). The pump is encoded by an operon (oqxABR) flanked by IS26 sequences, constituting a composite transposon (Tn6010) (Norman et al., 2008).

The oqxABR operon can form a circular structure to facilitate mobilization (He et al., 2015). In this study, we observed that oqxABR shifted from chromosomal to plasmid locations in isolates with relatively high MIC_CIP. We performed inverse PCR to determine whether an IS26-oqxAB circular intermediate was involved in this process. We found that all 16 oqxAB-positive strains possessed the circular intermediate regardless of genomic location. This indicated that the circular intermediate promotes mobilization of this PMQR between chromosome and plasmid. Interestingly, we detected qnrS in R strains located on the chromosome with aac(6′)-Ib-cr and oqxAB. This is consistent with other several studies, where qnrS co-existing with aac(6′)-Ib-cr and oqxAB were frequently detected in Salmonella spp. of food animal origin in China (Zhang et al., 2016; Wang et al., 2017; Kuang et al., 2018), and qnrS were also detected on the chromosome along with other PMQR genes (Lin et al., 2015). Co-existence of qnrS in the chromosome with other genes may be a strategy that confers a relatively high degree of tolerance to FQs without a prohibitive fitness cost which may explain the frequent observation of this combination.

The switch from plasmid to chromosomal location might represent a bacterial strategy for maintaining a resistance trait while adapting to antimicrobial stress (Machuca et al., 2014). Our growth competition assays confirmed that strains carrying PMQR genes on the chromosome had a growth advantage and could outcompete strains possessing PMQR genes on a plasmid. Interestingly, in the majority of our HR strains, only aac(6′)-Ib-cr was present. The latter gene was located on a plasmid and its presence coincided with multiple target gene mutations and overexpression of acrAB. aac(6′)-Ib-cr was associated with a fitness cost, in consistent with Machuca et al. (2014) observation in E. coli. Growth competition assays demonstrated that transfer of the aac(6′)-Ib-cr gene back to a plasmid location resulted in a fitness cost. Whole genome sequencing data identified a broad host range plasmid IncQ1 carrying no resistance genes in this HR strain. Thus, it is likely that the presence of this plasmid played a role in aac(6′)-Ib-cr gene movement from chromosome to plasmid, which has been previously suggested (Loftie-Eaton and Rawlings, 2012; Barry et al., 2019; Chen et al., 2019).

Increased Expression of acrAB-tolC Plays a Role in Isolates With Low Ciprofloxacin MICs via marA Rather Than ramA

Increased drug efflux primarily occurs via overexpression of AcrAB-TolC, the primary efflux pump, in Enterobacteriaceae and drug efflux is a primary mechanism for quinolone resistance in Salmonella (Baucheron et al., 2004; Blair et al., 2009, 2014). In accordance with studies that analyzed experimentally derived ciprofloxacin mutants (Redgrave et al., 2014; Hooper and Jacoby, 2015), we observed elevated expression of AcrAB in the early stages of resistance development prior to the occurrence of gyrA mutations (Vidovic et al., 2019). Significant levels of acrAB-tolC were observed in the RS strains.

Multidrug efflux pumps are expressed under precise transcriptional control. In Salmonella, RamA is often considered the master regulator of acrAB (Ricci et al., 2014). We did not detect the overexpression of ramA in strains that expressed high levels of acrAB but had relatively low ciprofloxacin MICs (Groups S and RS). Instead, this overexpression was accompanied by increases in the expression of marA. A similar observation was made with laboratory selected ciprofloxacin-resistant mutants where mutation of ramR was mainly seen in mutants with very high ciprofloxacin MICs (Vidovic et al., 2019). Interestingly, efflux pump genes were not expressed more in R strains compared with S. Typhimurium SL1344, in agreement with earlier observations made with MDR Klebsiella and MDR E. coli that strains capable of evolving multiple energetically favorable QRDR mutations will use less efflux than isolates with fewer favorable QRDR alterations (Toth et al., 2014; Johnson et al., 2015). The overexpression of ramA was only detected in highly resistant strains. Our data suggests that whilst both MarA and RamA can up-regulate acrAB, there is a preference for utilizing MarA in less resistant strains and RamA in the most resistant strains. Both MarA and RamA are pleiotropic regulators and the overexpression of ramA may play additional roles as well as regulating acrAB in highly resistant strains more compatible with this phenotype than over-expression of marA. This is in line with the previous research in which a deletion in the RamR-binding sites in the ramA promoter correlated with highest expression levels of efflux pumps (Fabrega et al., 2016).

In conclusion, we observed a complex interplay between a pool of quinolone resistance mechanisms, their contribution to susceptibility and fitness costs (Figure 4). In general, we found no significant differences between different serovars of Salmonella suggesting common impacts from these genotypes across the genus although only S. Indiana isolates were found in the HR group suggesting they may be more able to develop high level resistance. Whilst natural variations in genes related to phylogeny may mask impacts on resistance, the topoisomerase genes and other known quinolone associated resistance pathways are highly conserved across the Salmonella genus and studying a large and diverse panel of real-world isolates can be informative for mutations circulating in practice which have been selected in real conditions.

FIGURE 4

Figure 4. Quinolone resistance genes associated with different ciprofloxacin MICs amongst isolates.

We suggest an important role for substitutions within ParC and acrAB expression in the early stages of ciprofloxacin resistance development prior to the occurrence of gyrA mutations. Highly resistant isolates acquire multiple target site mutations and PMQR determinants although these are often integrated in the chromosome to minimize fitness costs and expression of acrAB-tolC is returned to wild-type levels, probably for the same reason. Whilst we are able here to identify genotypes which appear to be fit in conferring different levels of ciprofloxacin tolerance, it is still not possible to precisely assign a universal pathway in relation to the selection of resistance. Understanding the relative fitness of different genotypes in different selective conditions will help predict the evolutionary trajectory of AMR while also helping us to understand how resistance emerges in complex situations where multiple genes collaborate to confer resistance.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

H-XJ conceived and designed the experiments. M-XC, J-FZ, Y-HS, X-LL, and R-SL performed the experiments. M-XC, J-FZ, MW, and H-XJ analyzed the data. M-XC, Y-HS, and H-XJ contributed reagents, materials, and analyses tools. M-XC, LY, MW, and H-XJ wrote and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported in part by the National Natural Science Foundation of China (31772792) and the National Key Research Program of China (grant 2016YFD0501300). MW was supported by the Biotechnology and Biological Sciences Research Council (BBSRC) Institute Strategic Programme Microbes in the Food Chain BB/R012504/1 and its constituent project BBS/E/F/000PR10349.

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.663731/full#supplementary-material

References

Barry, K. E., Wailan, A. M., Sheppard, A. E., Crook, D., Vegesana, K., Stoesser, N., et al. (2019). Don’t overlook the little guy: an evaluation of the frequency of small plasmids co-conjugating with larger carbapenemase gene containing plasmids. Plasmid 103, 1–8. doi: 10.1016/j.plasmid.2019.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Baucheron, S., Imberechts, H., Chaslus-Dancla, E., and Cloeckaert, A. (2002). The AcrB multidrug transporter plays a major role in high-level fluoroquinolone resistance in Salmonella enterica serovar typhimurium phage type DT204. Microb. Drug Resist. 8, 281–289. doi: 10.1089/10766290260469543

PubMed Abstract | CrossRef Full Text | Google Scholar

Baucheron, S., Tyler, S., Boyd, D., Mulvey, M. R., Chaslus-Dancla, E., and Cloeckaert, A. (2004). AcrAB-TolC directs efflux-mediated multidrug resistance in Salmonella enterica serovar typhimurium DT104. Antimicrob. Agents Chemother. 48, 3729–3735. doi: 10.1128/aac.48.10.3729-3735.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Blair, J. M., La Ragione, R. M., Woodward, M. J., and Piddock, L. J. (2009). Periplasmic adaptor protein AcrA has a distinct role in the antibiotic resistance and virulence of Salmonella enterica serovar Typhimurium. J. Antimicrob. Chemother. 64, 965–972. doi: 10.1093/jac/dkp311

PubMed Abstract | CrossRef Full Text | Google Scholar

Blair, J. M., Richmond, G. E., and Piddock, L. J. (2014). Multidrug efflux pumps in Gram-negative bacteria and their role in antibiotic resistance. Future Microbiol. 9, 1165–1177. doi: 10.2217/fmb.14.66

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, K., Dong, N., Chan, E. W., and Chen, S. (2019). Transmission of ciprofloxacin resistance in Salmonella mediated by a novel type of conjugative helper plasmids. Emerg. Microbes Infect. 8, 857–865. doi: 10.1080/22221751.2019.1626197

PubMed Abstract | CrossRef Full Text | Google Scholar

CLSI (2019). Performance Standards for Antimicrobial Susceptibility Testing: Twenty-nine Informational Supplement M100-S29. Wayne, PA: Clinical and Laboratory Standards Institute.

Google Scholar

Eaves, D. J., Randall, L., Gray, D. T., Buckley, A., Woodward, M. J., White, A. P., et al. (2004). Prevalence of mutations within the quinolone resistance-determining region of gyrA, gyrB, parC, and parE and association with antibiotic resistance in quinolone-resistant Salmonella enterica. Antimicrob. Agents Chemother. 48, 4012–4015. doi: 10.1128/aac.48.10.4012-4015.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabrega, A., Balleste-Delpierre, C., and Vila, J. (2016). Differential impact of ramRA mutations on both ramA transcription and decreased antimicrobial susceptibility in Salmonella Typhimurium. J. Antimicrob. Chemother. 71, 617–624. doi: 10.1093/jac/dkv410

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabrega, A., Soto, S. M., Balleste-Delpierre, C., Fernandez-Orth, D., Jimenez De Anta, M. T., and Vila, J. (2014). Impact of quinolone-resistance acquisition on biofilm production and fitness in Salmonella enterica. J. Antimicrob. Chemother. 69, 1815–1824. doi: 10.1093/jac/dku078

PubMed Abstract | CrossRef Full Text | Google Scholar

Giraud, E., Cloeckaert, A., Kerboeuf, D., and Chaslus-Dancla, E. (2000). Evidence for active efflux as the primary mechanism of resistance to ciprofloxacin in Salmonella enterica serovar typhimurium. Antimicrob. Agents Chemother. 44, 1223–1228. doi: 10.1128/aac.44.5.1223-1228.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Gunell, M., Webber, M. A., Kotilainen, P., Lilly, A. J., Caddick, J. M., Jalava, J., et al. (2009). Mechanisms of resistance in nontyphoidal Salmonella enterica strains exhibiting a nonclassical quinolone resistance phenotype. Antimicrob. Agents Chemother. 53, 3832–3836. doi: 10.1128/aac.00121-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, L. H., Jensen, L. B., Sorensen, H. I., and Sorensen, S. J. (2007). Substrate specificity of the OqxAB multidrug resistance pump in Escherichia coli and selected enteric bacteria. J. Antimicrob. Chemother. 60, 145–147. doi: 10.1093/jac/dkm167

PubMed Abstract | CrossRef Full Text | Google Scholar

He, T., Wang, Y., Qian, M., and Wu, C. (2015). Mequindox resistance and in vitro efficacy in animal-derived Escherichia coli strains. Vet. Microbiol. 177, 341–346. doi: 10.1016/j.vetmic.2015.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Heisig, P. (1996). Genetic evidence for a role of parC mutations in development of high-level fluoroquinolone resistance in Escherichia coli. Antimicrob. Agents Chemother. 40, 879–885. doi: 10.1128/aac.40.4.879

PubMed Abstract | CrossRef Full Text | Google Scholar

Hohmann, E. L. (2001). Nontyphoidal salmonellosis. Clin. Infect. Dis. 32, 263–269.

Google Scholar

Hooper, D. C., and Jacoby, G. A. (2015). Mechanisms of drug resistance: quinolone resistance. Ann. N. Y. Acad. Sci. 1354, 12–31. doi: 10.1111/nyas.12830

PubMed Abstract | CrossRef Full Text | Google Scholar

Huseby, D. L., Pietsch, F., Brandis, G., Garoff, L., Tegehall, A., and Hughes, D. (2017). Mutation Supply and Relative Fitness Shape the Genotypes of Ciprofloxacin-Resistant Escherichia coli. Mol. Biol. Evol. 34, 1029–1039.

Google Scholar

Johnson, J. R., Johnston, B., Kuskowski, M. A., Sokurenko, E. V., and Tchesnokova, V. (2015). Intensity and mechanisms of fluoroquinolone resistance within the H30 and H30Rx subclones of Escherichia coli sequence type 131 compared with other fluoroquinolone-resistant E. coli. Antimicrob. Agents Chemother. 59, 4471–4480. doi: 10.1128/aac.00673-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Kronvall, G. (2010). Normalized resistance interpretation as a tool for establishing epidemiological MIC susceptibility breakpoints. J. Clin. Microbiol. 48, 4445–4452. doi: 10.1128/jcm.01101-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuang, D., Zhang, J., Xu, X., Shi, W., Chen, S., Yang, X., et al. (2018). Emerging high-level ciprofloxacin resistance and molecular basis of resistance in Salmonella enterica from humans, food and animals. Int. J. Food Microbiol. 280, 1–9. doi: 10.1016/j.ijfoodmicro.2018.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, D., Chen, K., Wai-Chi Chan, E., and Chen, S. (2015). Increasing prevalence of ciprofloxacin-resistant food-borne Salmonella strains harboring multiple PMQR elements but not target gene mutations. Sci. Rep. 5:14754.

Google Scholar

Ling, J. M., Chan, E. W., Lam, A. W., and Cheng, A. F. (2003). Mutations in topoisomerase genes of fluoroquinolone-resistant Salmonellae in Hong Kong. Antimicrob. Agents Chemother. 47, 3567–3573. doi: 10.1128/aac.47.11.3567-3573.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Loftie-Eaton, W., and Rawlings, D. E. (2012). Diversity, biology and evolution of IncQ-family plasmids. Plasmid 67, 15–34. doi: 10.1016/j.plasmid.2011.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Machuca, J., Briales, A., Labrador, G., Diaz-De-Alba, P., Lopez-Rojas, R., Docobo-Perez, F., et al. (2014). Interplay between plasmid-mediated and chromosomal-mediated fluoroquinolone resistance and bacterial fitness in Escherichia coli. J. Antimicrob. Chemother. 69, 3203–3215. doi: 10.1093/jac/dku308

PubMed Abstract | CrossRef Full Text | Google Scholar

Majowicz, S. E., Musto, J., Scallan, E., Angulo, F. J., Kirk, M., O’brien, S. J., et al. (2010). The global burden of nontyphoidal Salmonella gastroenteritis. Clin. Infect. Dis. 50, 882–889.

Google Scholar

Marcusson, L. L., Frimodt-Moller, N., and Hughes, D. (2009). Interplay in the selection of fluoroquinolone resistance and bacterial fitness. PLoS Pathog. 5:e1000541. doi: 10.1371/journal.ppat.1000541

PubMed Abstract | CrossRef Full Text | Google Scholar

Norman, A., Hansen, L. H., She, Q., and Sorensen, S. J. (2008). Nucleotide sequence of pOLA52: a conjugative IncX1 plasmid from Escherichia coli which enables biofilm formation and multidrug efflux. Plasmid 60, 59–74. doi: 10.1016/j.plasmid.2008.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Redgrave, L. S., Sutton, S. B., Webber, M. A., and Piddock, L. J. (2014). Fluoroquinolone resistance: mechanisms, impact on bacteria, and role in evolutionary success. Trends Microbiol. 22, 438–445. doi: 10.1016/j.tim.2014.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ricci, V., Blair, J. M., and Piddock, L. J. (2014). RamA, which controls expression of the MDR efflux pump AcrAB-TolC, is regulated by the Lon protease. J. Antimicrob. Chemother. 69, 643–650. doi: 10.1093/jac/dkt432

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez-Martinez, J. M., Machuca, J., Cano, M. E., Calvo, J., Martinez-Martinez, L., and Pascual, A. (2016). Plasmid-mediated quinolone resistance: two decades on. Drug Resist. Update 29, 13–29. doi: 10.1016/j.drup.2016.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, T., Yokota, S., Uchida, I., Okubo, T., Usui, M., Kusumoto, M., et al. (2013). Fluoroquinolone resistance mechanisms in an Escherichia coli isolate, HUE1, without quinolone resistance-determining region mutations. Front. Microbiol. 4:125. doi: 10.3389/fmicb.2013.00125

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, R., Swick, M. C., Ledesma, K. R., Yang, Z., Hu, M., Zechiedrich, L., et al. (2012). Temporal interplay between efflux pumps and target mutations in development of antibiotic resistance in Escherichia coli. Antimicrob. Agents Chemother. 56, 1680–1685. doi: 10.1128/aac.05693-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Toth, A., Kocsis, B., Damjanova, I., Kristof, K., Janvari, L., Paszti, J., et al. (2014). Fitness cost associated with resistance to fluoroquinolones is diverse across clones of Klebsiella pneumoniae and may select for CTX-M-15 type extended-spectrum beta-lactamase. Eur. J. Clin. Microbiol. Infect. Dis. 33, 837–843. doi: 10.1007/s10096-013-2022-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Vidovic, S., An, R., and Rendahl, A. (2019). Molecular and physiological characterization of fluoroquinolone-highly resistant Salmonella enteritidis strains. Front. Microbiol. 10:729. doi: 10.3389/fmicb.2019.00729

PubMed Abstract | CrossRef Full Text | Google Scholar

Vinue, L., Sater, M. R. A., Herriott, I., Huntley, M. H., Jacoby, G. A., and Hooper, D. C. (2019). Multiple copies of qnrA1 on an IncA/C2 plasmid explain enhanced quinolone resistance in an Escherichia coli mutant. Antimicrob. Agents Chemother. 63:e00718-19.

Google Scholar

Wang, J., Zhi, C. P., Chen, X. J., Guo, Z. W., Liu, W. L., Luo, J., et al. (2017). Characterization of oqxAB in Escherichia coli isolates from animals, retail meat, and human patients in Guangzhou, China. Front. Microbiol. 8:1982. doi: 10.3389/fmicb.2017.01982

PubMed Abstract | CrossRef Full Text | Google Scholar

WHO (2017). WHO list of Critically Important Antimicrobials for Human Medicine (WHO CIA list). Genev: WHO.

Google Scholar

Yang, L., Li, W., Jiang, G. Z., Zhang, W. H., Ding, H. Z., Liu, Y. H., et al. (2017). Characterization of a P1-like bacteriophage carrying CTX-M-27 in Salmonella spp. resistant to third generation cephalosporins isolated from pork in China. Sci. Rep. 7:40710.

Google Scholar

Zayed, A. A., Essam, T. M., Hashem, A. G., and El-Tayeb, O. M. (2015). ‘Supermutators’ found amongst highly levofloxacin-resistant E. coli isolates: a rapid protocol for the detection of mutation sites. Emerg. Microbes Infect. 4:e4.

Google Scholar

Zhang, C. Z., Ding, X. M., Lin, X. L., Sun, R. Y., Lu, Y. W., Cai, R. M., et al. (2019). The emergence of chromosomally located bla CTX-M-55 in Salmonella from foodborne animals in China. Front. Microbiol. 10:1268. doi: 10.3389/fmicb.2019.01268

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C. Z., Ren, S. Q., Chang, M. X., Chen, P. X., Ding, H. Z., and Jiang, H. X. (2017). Resistance mechanisms and fitness of Salmonella Typhimurium and Salmonella Enteritidis mutants evolved under selection with ciprofloxacin in vitro. Sci. Rep. 7:9113.

Google Scholar

Zhang, W. H., Lin, X. Y., Xu, L., Gu, X. X., Yang, L., Li, W., et al. (2016). CTX-M-27 Producing Salmonella enterica Serotypes typhimurium and Indiana are prevalent among food-producing animals in China. Front. Microbiol. 7:436. doi: 10.3389/fmicb.2016.00436

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W.-H., Zhang, C.-Z., Liu, Z.-J., Gu, X.-X., Li, W., Yang, L., et al. (2017). In vitrodevelopment of ciprofloxacin resistance of Salmonella enterica Serovars typhimurium, Enteritidis, and Indiana isolates from food animals. Microb. Drug Resist. 23, 687–694. doi: 10.1089/mdr.2016.0119

PubMed Abstract | CrossRef Full Text | Google Scholar