Antibiotic susceptibility testing and resistance gene screening of Providencia strains
A total of 45 nonredundant strains of five Providencia species were analyzed, and the strains included 29 Providencia rustigianii isolates, 5 Providencia alcalifaciens isolates, 5 Providencia stuartii isolates, 4 Providencia rettgeri isolates, and 2 Providencia vermicola isolates. The bacteria from 22 freshwater fish were all P. rustigianii isolates, while the 2 P. vermicola isolates were obtained from the animal specimens only. Thirteen isolates from patients were species of P. alcalifaciens, P. rettgeri and P. stuartii (Table S1).
The MIC results for 15 antibiotics indicated that the 45 isolates had higher resistance rates to levofloxacin (77.8%, 35/45) and tetracycline (53.3%, 24/45) and lower resistance rates to chloramphenicol (35.6%, 16/45), imipenem 14 (31.1%, 14/45), norfloxacin (24.4%, 11/45) and florfenicol (17.8%, 8/45), while more than 90% of the isolates were sensitive to the other 9 antibiotics, especially meropenem with all the strains sensitive to it. The clinical isolates had a wider resistance spectrum and showed resistance to 13 antibiotics (ranging from 7.7% to 100%). The animal isolates showed resistance to 8 antibiotics (ranging from 10.0% to 80.0%), while the freshwater fish bacteria only showed resistance to 5 antibiotics (ranging from 4.5% to 50.0%). Notably, although the clinical strains showed a wider resistance spectrum, the animal strains showed higher resistance rates than the clinical strains in 5 of 9 antibiotics, which they both showed resistance to, especially to florfenicol (the animal strains with a resistance rate of 50%, while the human clinical isolates was only 23.1%) (Table 3, Table S2).
Detection of the florfenicol resistance genes demonstrated that of the 7 genes (floR, fexA, fexB, optrA, estDL136,cfr and pexA) screened, only the floR gene yielded positive results in both human (3/13, 23.1%) and animal bacteria (5/10, 50%). The floR-positive strains had high MIC levels for florfenicol (> 32 mg/L) and chloramphenicol, while the floR-negative isolates showed notably lower MIC levels (with all of them exhibiting MIC levels < 8 mg/L) (Table S2).
General features of the P. rettgeri R39 genome
To analyze the molecular mechanism governing the florfenicol resistance of the bacteria from animal sources, the whole genome of P. rettgeri isolate R39, which exhibited the highest florfenicol MIC level (Table S2), was sequenced. The R39 genome was determined to consist of a circular chromosome and 4 circular plasmids. The chromosome was 4,502,622 bp in size with 40.15% GC content encoding 4,029 ORFs. The four plasmids were named pR39-33 (32,936 bp, 32 ORFs), pR39-8 (8,234 bp, 8 ORFs), pR39-2.8 (2,888 bp, 6 ORFs) and pR39-2.7 (2,705 bp, 3 ORFs), respectively (Table 4, Fig. 1). The whole genome was observed to encode 9 resistance genes, of which 4 (aadA1, sat2, catB2 and drfA1) were located on the chromosome, while the other 5 floR genes were located on the plasmids, with 4 floR genes on pR39-33 and one floR gene on pR39-8 (Fig. 1a, Fig. 1b). According to the drug susceptibility results, the drug resistance gene profile was consistent with the drug resistance spectrum of the bacterium (Table S2).
pR39-33, a plasmid composed of four copies of the pR39-8 sequence
PFGE analysis showed that R39 contained one plasmid of approximately 33 kb in length (Fig. 2a), while in the initial sequencing assembly, we only obtained the circular chromosome and three circular plasmids pR39-8, pR39-2.8, pR39-2.7 and a contig measuring approximately 25 kb in length. Annotation of the 25-kb contig sequence demonstrated that it contained three copies and more sequences of pR39-8 in the form of a tandem repeat. To determine if the largest plasmid was composed of four copies of pR39-8 sequences, the restriction physical map of the largest plasmid was analyzed. The physical map of pR39-8 demonstrated that it contained EcoRV (at position 4,381) and BamHI (at position 7,937) restriction sites 3,556 bp apart. When the plasmid was digested with both BamHI and EcoRV, two fragments of 3,556 bp and 4,678 bp were produced. We digested the plasmids extracted from R39 with BamHI or EcoRV alone, or with the combination of BamHI and EcoRV, and the results were obtained as predicted. When digested with BamHI or EcoRV alone, a single band of approximately 8 kb was obtained. When digested with BamHI + EcoRV, two bands of approximately 3.5 kb and 4.5 kb were observed (Fig. 2b). The digestion results of the recovered pR39-33 with BamHI or EcoRV alone or with both BamHI and EcoRV also yielded the expected results (Fig. 2c). This result confirmed that the structure of pR39-33 consisted of a tandem repeat of 4 copies of pP39-8 sequences.
The floR genes encoded on the plasmids
The plasmids (pR39-33 and pR39-8) cured R39 derivative (R39ΔpR39-33ΔpR39-8) exhibited notably decreased florfenicol MIC values. The florfenicol MIC decreased more than 128-fold from > 512 mg/L (the wild-type strain R39) to 4 mg/L (the plasmids cured R39, R39ΔpR39-33ΔpR39-8) (Table 5). The recombinant floR ORF together with its promoter region (pUC18–pro-floR/DH5α) conferred a much higher MIC level to florfenicol (256 mg/L) than the recipient E. coli DH5α (8 mg/L) (Table 5).
Comparative genomic analysis of the floR gene encoding plasmids
As mentioned above, R39 carried 2 resistant plasmids (pR39-33 and pR39-8), and pR39-33 is composed of a tandem repeat of four copies of the pR39-8 sequence. pR39-8 encoded 7 complete ORFs and a truncated insertion sequence (ΔIS91). These seven genes included a floR gene, 2 replication-related genes (repA and repB), 2 mobilization-related genes (mobA and mobC) and 2 hypothetical protein genes. The floR gene was related to a transposon-like structure (ΔIS91-hp-floR).
The five plasmids with the highest similarity (> 90% identity and > 83% coverage) to pR39-8 were retrieved from the NCBI nucleotide database (https://www.ncbi.nlm.nih.gov/nucleotide/). These plasmids have five or more genes that are the same as those encoded on pR39-8. The six plasmids described in this study were from bacteria of four different genera. Specifically, both plasmid pRCADGH-1 and plasmid pCCK381 were from Pasteurella multocida, the plasmids pM3446F and unamed1 were from Actinobacillus pleuropneumoniae, pMh1405 is a Mannheimia haemolytica plasmid, and pR39-8 described in this work was from a strain of P. rettgeri. Of the five plasmid sequences obtained from the database, pRCADGH-1 exhibited the same sequence as pR39-8, but the plasmids were from two different genera, that is, Pasteurella multocida and Providencia rettgeri, respectively (Fig. 3, Table S3).