International Journal of Genomics
Volume 2019 (2019), Article ID 5459190, 12 pages
https://doi.org/10.1155/2019/5459190

Research Article

Molecular Characterization of a Multidrug-Resistant Klebsiella pneumoniae Strain R46 Isolated from a Rabbit

Fei Wu,¹ Yuanyuan Ying,¹ Min Yin,¹ Yi Jiang,¹ Chongyang Wu,¹ Changrui Qian,¹ Qianqian Chen,¹ Kai Shen,¹ Cong Cheng,² Licheng Zhu,² Kewei Li,¹ Teng Xu,³ Qiyu Bao,¹ and Junwan Lu¹

¹School of Laboratory Medicine and Life Sciences/Institute of Biomedical Informatics, Wenzhou Medical University, Wenzhou, Zhejiang 325035, China
²College of Medicine and Health, Lishui University, Lishui 323000, China
³Institute of Translational Medicine, Baotou Central Hospital, Baotou 014040, China

Correspondence should be addressed to Qiyu Bao; ym888888ym@163.com and Junwan Lu; 304041474@qq.com

Received 12 January 2019; Revised 16 May 2019; Accepted 12 July 2019; Published 18 August 2019

Academic Editor: Wilfred van IJcken

Copyright © 2019 Fei Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

To investigate the mechanisms of multiple resistance and the horizontal transfer of resistance genes in animal pathogens, we characterized the molecular structures of the resistance gene-related sequences in a multidrug-resistant Klebsiella pneumoniae strain R46 isolated from a rabbit. Molecular cloning was performed to clone the resistance genes, and minimum inhibitory concentrations (MICs) were measured to determine the resistance characteristics of the cloned genes and related strains. A conjugation experiment was conducted to assess the transferability of the resistance plasmids. Sequencing and comparative genomic methods were used to analyze the structures of the resistance gene-related sequences. The K. pneumoniae R46 genome consisted of a chromosome and three resistance plasmids named pR46-27, pR46-42, and pR46-270, respectively. The whole genome encoded 34 antibiotic resistance genes including a newly identified chromosome-encoded florfenicol resistance gene named mdfA2. pR46-270, besides encoding 26 antibiotic resistance genes, carried four clusters of heavy metal resistance genes and several virulence-related genes or gene clusters. The plasmid-encoded resistance genes were mostly associated with mobile genetic elements. The plasmid with the most similarity to the floR gene-harboring plasmid pR46-27 was pCTXM-2271, a plasmid from Escherichia coli. The results of this work demonstrated that the plasmids with multidrug resistance genes were present in animal-derived bacteria and more florfenicol resistance genes such as mdfA2 could be present in bacterial populations. The resistance genes encoded on the plasmids may spread between the bacteria of different species or genera and cause the resistance dissemination.

1. Introduction

Klebsiella pneumoniae, a member of the Enterobacteriaceae, is pervasive in the natural environment and benignly colonizes the gastrointestinal tracts of humans and animals. It is an opportunistic pathogen capable of causing a wide range of diseases in humans and different animal species [1]. In animals, infections of the urinary and respiratory tracts and sepsis are the most common clinical manifestations. In addition, K. pneumoniae species are a well-documented cause of mastitis in cattle [2]; endometritis, cystitis, and liver abscess in horses; tracheitis and wounds in birds; cystitis, phlebitis, and otitis externa in dogs; mastitis and wounds in cows; and cystitis in cats, among others [3].

Notably, the prevalence of antibiotic resistance is increasing among Enterobacteriaceae, including K. pneumoniae strains isolated from animals [1]. Due to the extensive use of antibiotics in humans, veterinary medicine, and agricultural practice during the last few decades, the emergence of K. pneumoniae strains that harbor various resistance genes has increased considerably. K. pneumoniae has acquired increasingly high levels of antimicrobial drug resistance. Many publications have reported that K. pneumoniae isolated from animals, such as pigs, chickens, and dogs, carries various resistance genes and shows resistance to a variety of antibiotics, such as aminoglycosides (aadA1, aacC2, and aacA4) [4], colistin (mcr-1), amphenicols (floR) [5], sulfonamides (sul1 and sul2) [6], quinolones (qnrA, qnrB, qnrS, and qepA) [7], and β-lactams (CTX-M groups) [8]. Multidrug-resistant (MDR), extensively drug-resistant (XDR) [9], and pandrug-resistant (PDR) [10] isolates of K. pneumoniae have been reported worldwide.

The antibiotic resistance, especially multiple-drug resistance bacteria of different species in both human and animal pathogens, becomes an increasing global problem. The multidrug resistance in the microbial organisms, in particular the Enterobacteriaceae, is often due to the acquisition of resistance genes from a shared pool [11]. The resistance genes in this pool appeared to have been captured from the chromosomes of various species by the mobile genetic elements, such as insertion sequences (IS) [12], transposons [13], and integrons [14]. A great number of genotypes (subgenotypes) of the antibiotic resistance genes have been found encoded on the mobile genetic elements. These mobile genetic elements can translocate between different bacterial DNA molecules, such as chromosomes and plasmids, especially the conjugative plasmids, and can be transferred between different bacterial species or genera through horizontal gene transfer, causing resistance dissemination [15].

With the widespread use of antimicrobial agents and detergents without effective supervision and the greater incidence of animal-caused foodborne illnesses, animals, especially those in agricultural practice, may become potential reservoirs for the dissemination of antimicrobial resistance. In this work, to investigate the resistance characteristics of K. pneumoniae isolated from animals, we used molecular cloning and comparative genomic methods to identify resistance gene profiles and to assess the horizontal transfer of resistance genes in a K. pneumoniae strain isolated from a rabbit in a farm in Wenzhou, South China.

2. Materials and Methods

2.1. Bacterial Isolation and Identification

K. pneumoniae R46 was isolated on a MacConkey agar (MCA) plate with florfenicol (8 mg/L) from an anal swab sample of a rabbit from a farm in Wenzhou, South China. Species identification was conducted with the Vitek-60 microorganism auto-analysis system (BioMerieux Corporate, Lyon, France). Further verification was performed by using homology comparisons of the 16S ribosomal RNA gene and the whole genome sequence from NCBI nucleotide database (http://www.ncbi.nlm.nih.gov) via BLASTN software. Multilocus sequence typing (MLST) was performed by analyzing the housekeeping genes gapA, infB, mdh, pgi, phoE, rpoB, and tonB [16] via the online service MLST-1.8 (https://cge.cbs.dtu.dk/services/MLST/) [17].

2.2. Antimicrobial Susceptibility Testing

The agar dilution method was used to determine the minimum inhibitory concentrations (MICs) of the antibiotics against bacteria in accordance with the guidelines of the Clinical and Laboratory Standards Institute (CLSI document M100-S28, 2017) and European Committee on Antimicrobial Susceptibility Testing (EUCAST 2017). All tests were carried out in triplicate. Briefly, each isolate was diluted in the logarithmic growth period with 0.9% saline to achieve a 0.5 McFarland standard suspension, equivalent to an inoculum with . Suspensions of approximately 5 μL were plated on Mueller-Hinton (MH) agar plates (with a series of antibiotic concentrations between 0.0625 and 1,024 mg/L). The inoculated plates were then incubated for 18 hours at 37°C in a constant-temperature incubator. The MIC was defined as the lowest antibiotic concentration showing no colony growth. The commercially available E. coli strain ATCC 25922 served to ensure the quality of the plates, and E. coli DH5α carrying an empty cloning vector was used for comparison.

2.3. Cloning of K. pneumoniae R46 Resistance Genes

The genomic DNA of K. pneumoniae R46 was extracted using an AxyPrep Bacterial Genomic DNA Miniprep Kit (Axygen Scientific, Union City, CA, USA). To clone resistance genes, the primers of the ORFs with their upstream promotor regions were designed by using Primer Premier 5.0 and were synthesized by Shanghai Sunny Biotechnology Co. Ltd. (Shanghai, China) (Table 1). A standard protocol for PCR with PrimeSTAR HS DNA Polymerase (TaKaRa, Dalian, China) was performed under the following conditions: an initial cycle of 94°C for 1 min, followed by 35 cycles of 10 s at 94°C, 15 s at a specific annealing temperature, and then an extension step of 1 min/kb at 72°C, with a final extension step of 10 min at 72°C (Table 1). The PCR products were cloned into the pMD™19-T vector (TaKaRa, Dalian, China) (Table 2). The resulting recombinant plasmids were transformed into E. coli DH5α using the calcium chloride method, and bacterial colonies were grown on Luria-Bertani agar plates supplemented with ampicillin (100 mg/L). The recombinant plasmids were isolated and PCR-amplified to confirm the sizes of the inserted fragments. The cloned fragments from the transformants were further characterized by Sanger sequencing (ABI3730 Analyzer, USA) and compared to reference resistance genes using the BLASTN program.

Table 1: Primers used in this study.

Resistance genes
Primer
Sequence (5-3)b
Restriction endonuclease
Vector
Amplicon size (bp)
Annealing temperature (°C)

catB3a
P-catB3-F
GGGGTACCACTTACAGGAAACTTGGGGT
Kpn I
pMD™19-T
859
58

P-catB3-R
CGGGATCCTTAGACGGCAAACTCGAGCC
Bam HI

floRa
P-floR-F
CGGGATCCGAAGCAAAAGATAATCGGAT
Bam HI
pMD™19-T
1412
54

P-floR-R
CCAAGCTTTTAGACGACTGGCGACTTCT
Hind III

mdfA2a
P-mdfA2-F
CGGGATCCTCACATTGCTGAAACATAAACGG
Bam HI
pMD™19-T
1400
55

P-mdfA2-R
CCAAGCTTCTACCCCTGCTGCGAATTGC
Hind III

qnrB2a
P-qnrB2-F
CGGGATCCGATTTGACGCATAACCTCAT
Bam HI
pMD™19-T
713
55

P-qnrB2-R
CCAAGCTTCTAGCCAATAATCGCGATGC
Hind III

aThe primers for ORFs with the predicted promoter regions. bThe underlines present the restriction enzyme sites and their protective bases.

Table 2: Strains and plasmids used in this study.

Strains and plasmids
Relevant characteristic(s)
Reference or source

R46
A K. pneumoniae strain named R46 isolated from a rabbit
This study

DH5α
Escherichia coli DH5α was used as a host for the recombinant plasmid with the cloned resistance gene
Our lab collection

ATCC 25922
Escherichia coli ATCC 25922 is an FDA clinical isolate
Our lab collection

pMDTM19-T-ORFs/DH5α
DH5α carrying the recombinant plasmids pMDTM19-T-ORF (floR, mdfA2, qnrB3, and catB3)
This study

E. coli C600
E. coli C600 as recipient in the conjugation experiment, RIFr
Our lab collection

Plasmids

pMDTM19-T
Cloning vector for the PCR products of resistance genes of floR, mdfA2, qnrB3, and catB3, AMPr
This study

RIF: rifampin; AMP: ampicillin.

2.4. Conjugation Experiment

A conjugation experiment was carried out by biparental mating on sterile nitrocellulose filters, and rifampin-resistant E. coli C600 (EC600) was used as the recipient. Overnight cultures of the donor strain K. pneumoniae R46 and the recipient strain EC600 were mixed together in 3 mL of LB broth, harvested, and resuspended in 80 μL of LB. The mixture was spotted on a 1 cm² filter membrane, placed on an LB plate, and then incubated for mating at 37°C for 12-18 hours. The transconjugants were selected on two kinds of Mueller-Hinton agar plates: one containing 1,024 mg/L rifampin and 32 mg/L florfenicol and the other containing 1,024 mg/L rifampin and 32 mg/L chloramphenicol. The plasmid DNA of the transconjugant was extracted and analyzed by agarose gel electrophoresis and compared to the plasmid profile of K. pneumoniae R46. PCR analysis of the genes encoded in the plasmid was performed to further verify the plasmid in the transconjugants.

2.5. Sequencing and Annotations of the K. pneumoniae R46 Genome

Genomic DNA was extracted from K. pneumoniae R46 as mentioned above and sequenced with Pacific Biosciences sequencers at Annoroad Gene Technology Co. LTD (Beijing, China). In addition, a paired-end library with 300 bp insert sizes was constructed and sequenced from both ends with Illumina HiSeq-2500. Pacific Biosciences sequencing reads of approximately 10-20 kb in length were assembled by Canu v1.2 [18]. The Illumina reads were then mapped onto the assembled contigs to correct the primary assembly by using bwa0.7.13, samtools1.3, and GenomeAnalysisTK2.3.9 [19, 20]. Glimmer3.02 software with default parameters was used to predict potential open reading frames (ORFs) [21]. ORF annotations were determined by performing BLASTX comparisons with the NCBI nonredundant protein database. Comparisons of nucleotide sequences and amino acid sequences were performed by BLASTN and BLASTP, respectively [22]. BLASTP was applied to compare amino acid sequences with those in the Antibiotic Resistance Genes Database (ARDB). A map of the plasmid with GC content and GC skew was drawn with the online CGView Server (http://stothard.afns.ualberta.ca/cgview_server/) and local GView 1.7 with a visual interface [23]. The rRNA gene sequences were annotated by the online tool RNAmmer (http://www.cbs.dtu.dk/services/RNAmmer/) [24], and the tRNA sequences were annotated by the online tool tRNAscan-SE 2.0 (http://lowelab.ucsc.edu/tRNAscan-SE/) [25]. Promoter sites were determined by using Soft Berry BPROM software (http://linux1.softberry.com/berry.phtml?topic=bprom&group=programs&subgroup=gfindb). The plasmid sequences used in this study were downloaded from the NCBI nucleotide database (http://www.ncbi.nlm.nih.gov).

2.6. Comparative Genomic Analysis

Sequences containing resistance genes were obtained from the NCBI nucleotide database by the BLASTN program using the resistance gene sequences of K. pneumoniae R46 as the query. The resulting sequences were filtered, and only sequences containing resistance genes were retained. CD-HIT was used to cluster the retained sequences using the genome sequence of K. pneumoniae R46 as the reference with identity of 90% and coverage of 85% [26]. The sequence sharing the greatest similarity to the other sequences in each cluster was chosen as the candidate for ortholog analysis. Orthologous groups of the genes from the candidate sequences were identified using BLASTP [22]. Sequence retrieval, statistical analysis, and other bioinformatic tools used in this study were applied with Perl and Bioperl scripts (http://www.perl.org/).

2.7. Nucleotide Accession Numbers

The sequences are now available in the NCBI nucleotide database. The accession numbers for the chromosome and the plasmids are CP035777, CP035774, CP035775, and CP035776, respectively.

3. Results

3.1. Isolation and Identification of K. pneumoniae R46

K. pneumoniae R46 is a gram-negative bacillus with small, nonhemolytic, gray, and mucoid colonies after overnight culture on blood agar plates. Whole-genome comparative analysis combined with 16S ribosomal RNA gene homology analysis showed that the most closely related species were three Klebsiella pneumoniae strains: Klebsiella pneumoniae strain INF042 (CP024542), Klebsiella pneumoniae strain INF059 (CP024545), and Klebsiella pneumoniae strain KSB1_7J (CP024548) all with coverages of 99% and identities of 99%. MLST of seven housekeeping genes from the K. pneumoniae R46 genome showed optimum matching to the same housekeeping genes of the K. pneumoniae ST37 genome, with identities of 100.00%, 100.00%, 100.00%, 99.77%, 100.00%, 100.00%, and 100.00%, respectively. Finally, we classified the strain into the genus K. pneumoniae and named it K. pneumoniae R46 with MLST type ST37.

3.2. General Features of the K. pneumoniae R46 Genome

The K. pneumoniae R46 genome consists of a 5.12 Mb chromosome encoding 4,701 ORFs and three plasmids, namely, pR46-27, pR46-42, and pR46-270 (Table 3, Figures 1(a) and 1(b)). A total of 34 antibiotic resistance genes were identified in the whole genome, of which five, one, two, and twenty-six were encoded in the chromosome and in the plasmids pR46-27, pR46-42, and pR46-270, respectively (Table 4). The resistance genes were mainly related to the antibiotics streptomycin, aminoglycosides, β-lactams, and amphenicols. Moreover, four clusters of heavy metal resistance genes (operons related to silver, copper, arsenic, and mercury resistance) were identified in pR46-270 (Table 4). Interestingly, pR46-27 and pR46-42 belonged to the same incompatibility group IncX, with the X1 and X4 replicons, respectively, while pR46-270 was a multireplicon plasmid with three replicons: Q, FII(k), and FIB(k).

Table 3: General features of K. pneumoniae R46 genome.

Chromosome
pR46-27
pR46-42
pR46-270

Size (bp)
5,117,042
27,056
42,640
270,566

G+C (%)
57.69
50.08
41.57
51.96

Total opening reading frames
4,701
37
68
307

Known proteins
4,125 (87.75%)
25 (67.57%)
38 (55.88%)
251 (81.76%)

Hypothetical proteins
576 (12.25%)
12 (32.43%)
30 (44.12%)
56 (18.24%)

Protein coding sequence (%)
87.3
82.05
84
83.97

Average ORF length (bp)
950
600
526
740

rRNA operons
1(16s-23s-5s-5s)
1(16s-23s-5s)
5(5s-23s-16s)

tRNA
87

(a)

(b)

Figure 1: Circular maps of the plasmids pR46-27 (a) and pR46-270 (b). (a, b) Counting from the center toward the outside: (1) the inner most circle shows the position in kb. (2) GC skew (G-C/G+C), with a positive GC skew toward the outside and a negative GC skew toward the inside. (3) GC content, with an average of 50%, whereby a G+C content of more than 50% is shown toward the outside, otherwise, inward. (4) Genes encoded in the leading strands (outwards) or lagging strands (inwards). Genes with different functions are shown in different colors: red: drug resistance; yellow: mobile genetic elements; orange: transfer conjugation; green: virulence; blue: heavy metal resistance; gray: genes with unknown functions; light green: genes with other functions. (a) The plasmid pR46-27 was used as the reference genome and compared to the sequences of pCTXM-2771 (MF589339), pACN001-A (KC853434), and p160070-CTXM (MG288677). (b) R1: multidrug resistance region; R2: transfer and maintenance region; R3: heavy metal resistance and virulence region.

Table 4: Resistance genes encoded on the K. pneumoniae R46 genome.

Genome
Class of resistance genes
Resistance genes

Chromosome
Quinolone
oqxA, oqxB

Fosfomycin
fosA

β-Lactams
blaSHV-1

Amphenicols
mdfA2

pR46-27
Amphenicols
floR

pR46-42
β-Lactams
blaCTX-M-99, blaCTX-M-14

pR46-270
β-Lactams
blaOXA-1

Tetracycline
tetA, tetR

Macrolide
mef(B), mphA, mrx, mphR

Sulfonamide
2sul1, sul2, sul3

Quaternary ammonium compounds
qacH, qacEΔ1

Trimethoprim
dfrA12

Quinolone
qnrB2

Rifampin
arr-2

Aminoglycosides
aacA4, aac3, aadA1, aadA2, cmlA1, strA, strB, aphA1a

Amphenicols
catB3

Tunicamycin
tmrB

Silver
silA, silB, silC, silE, silF, silR, silS

Copper
copA, copB, copC, copD, pcoE, copG, pcoR, pcoS

Arsenic
arsR, arsD, arsA, arsB, arsC

Mercury
merA, merC, merP, merE, merD, merT, merR

3.3. MIC Results for K. pneumoniae R46 and the Cloned Resistance Genes

The MIC results for sixteen antibiotics showed that the wild-type K. pneumoniae R46 was resistant to eleven antibiotics (68.75%, 11/16), including florfenicol (>512 mg/L), chloramphenicol (256 mg/L), nalidixic acid (>32 mg/L), neomycin (>8 mg/L), ampicillin (512 mg/L), tetracycline (512 mg/L), tigecycline (4 mg/L), aztreonam (32 mg/L), gentamicin (32 mg/L), streptomycin (>32 mg/L), and kanamycin (512 mg/L). The cloned genes also conferred resistance to their corresponding antibiotics at specific levels (Table 5). Interestingly, the cloned mdfA2 gene, in addition to conferring resistance to chloramphenicol (32 mg/L), also conferred resistance to florfenicol (32 mg/L) with approximately 4-fold increases in the MIC levels over those of the controls DH5α (Table 5). The protein sharing the greatest amino acid sequence identity with MdfA2 was MdfA (85.37%, accession number AFH35853) among the function known proteins in ARDB.

Table 5: MICs (μg/mL) for K. pneumoniae R46, its cloned genes, and transconjugant (R46/EC600).

Strain
FFC
CHL
NOR
AMK
NAL
NEOa
AMP
TCY
TGCa
RIF
POL
ATM
GEN
CNX
STR
KAN

ATCC 25922
4
4
0.25
2
2
4
8
1
0.25
4
0.5
<0.5
1
0.5
8
2

E. coli C600
8
4
4
1
2
>8
4
0.5
2
1,024
0.5
<0.5
8
1
2
1

Transconjugant (R46/EC600)
512
128
8
1
16
>8
512
512
0.25
512
<0.0625
32
32
<0.0625
>32
512

K. pneumoniae R46
>512
256
8
4
>32
>8
512
512
4
16
1
32
32
1
>32
512

pMD™19-T-floR/DH5α
128
64
<0.0625
<0.0625
2
4
512
0.5
0.25
1
0.25
<0.5
1
0.5
<0.5
2

pMD™19-T-mdfA2/DH5α
32
32
0.25
1
4
4
>512
1
0.25
2
0.25
<0.5
1
0.5
4
2

pMD™19-T-qnrB2/DH5α
8
8
8
2
16
2
>512
1
0.25
8
0.5
<0.5
1
0.5
4
2

pMD™19-T-catB3/DH5α
8
8
<0.0625
2
8
4
>512
1
0.25
2
0.25
<0.5
1
0.5
4
2

E. coli DH5α
8
4
<0.0625
0.5
2
2
<1
0.5
0.25
2
0.25
<0.5
1
0.25
2
2

FFC: florfenicol; CHL: chloramphenicol; NOR: norfloxacin; AMK: amikacin; NAL: nalidixic; NEO: neomycin; AMP: ampicillin; TCY: tetracycline; TGC: tigecycline; RIF: rifampin; POL: polymyxin; ATM: aztreonam; GEN: gentamicin; CNX: cefminox; STR: streptomycin; KAN: kanamycin. aFor neomycin and tigecycline, we used EUCAST 2017 as the guideline.

3.4. Comparative Analysis of the Sequence Region Containing the mdfA Gene

A set of 979 sequences totaling approximately 9 kb in length with the mdfA gene in the center (with ≥79% nucleotide sequence identities and ≥88% sequence coverage) was retrieved from the NCBI nucleotide database. Of the 980 sequences (including one in K. pneumoniae R46), most were from Klebsiella (28.06%, 275/980), Escherichia (32.65%, 320/980), and Salmonella (30.71%, 301/980). The sequences were clustered into 53 clusters with a nucleotide sequence identity of 90% and coverage of 85%. The sequences sharing the greatest similarity in each cluster were chosen as candidates for further grouping according to the genes flanking mdfA. Finally, 11 groups were obtained. Most of the sequences (92.04%, 902/980) were grouped into 2 groups, group 7 (59.08%, 579/980) and group 11 (32.96%, 323/980), and the mdfA2-containing fragment of K. pneumoniae R46 was in group 7 (Table S1). Eleven representative sequences from 11 groups (one from each group) are illustrated according to their accession numbers (Figure 2, Table S1). The results of this gene neighborhood analysis revealed that a few of the genes downstream of mdfA were conserved and a sequence of approximately 4 kb encoding mdfA-ybjG-deoR-dacC was present in 93.06% (912/980) of the sequences. In contrast, a portion of the genes flanking mdfA from the same or different sequence groups differed greatly.

Figure 2: Comparative analysis of the mdfA-related regions of 11 representatives from 980 sequences. The representative sequences were derived from the following bacteria: CP004091 (C. sakazakii SP291 CP004091), CP026975 (E. cloacae complex bacterium FDAARGOS_77 CP026975), CP009451 (C. neteri SSMD04 CP009451), CP009459 (C. neteri ND14a CP009459), CP011574 (K. aerogenes CAV1320 CP011574), CP026192 (Enterobacteriaceae bacterium ENNIH2 CP026192), CP026193 (Enterobacteriaceae bacterium ENNIH1 CP026193), CP023504 (C. werkmanii FDAARGOS_364 CP023504), CP019113 (Enterobacter sp. SA187 CP019113), and CP005991 (Enterobacter sp. R4-368 CP005991). Homologous genes are shown in the same colors, whereas the white arrows indicate nonhomologous genes.

The sequence structures next to mdfA in the same genus were relatively conserved. Most of the sequences from Klebsiella (86.90%, 239/275) were clustered together in cluster 1, nearly all of the sequences from Escherichia (98.75%, 316/320) were in cluster 41, and almost all of the sequences from Salmonella (99.34%, 299/301) were in cluster 47 (Table S1). All of the sequences from Escherichia (320/320, 100%) in group 5 shared the same gene structure: apbC-metG-mdfA-ybjG-deoR. Among the 579 sequences in group 3 with the structure ybjJ-supH-mdfA-hp-ybjG, 576 were from Klebsiella (263) and Salmonella (301), representing 95.64% (263/275) of the Klebsiella sequences and 100% (301/301) of the Salmonella sequences (Table S1).

3.5. Comparative Analysis of the floR Gene-Carrying Plasmid pR46-27

The three plasmids sharing the highest nucleotide sequence similarities (, ) with pR46-27 were retrieved from the NCBI nucleotide database. Two of the plasmids, pCTXM-2271 (MF589339, coverage 66%) from Escherichia coli 2271 and pACN001-A (KC853434, coverage 64%) from Escherichia coli ACN001, harbored floR gene, while the other one, p160070-CTXM (MG288677, coverage 63%) from Klebsiella pneumoniae F160070, lacked the floR gene. pCTXM-2271 is 222 kb in length and encodes 271 ORFs, which is 195 kb and 234 ORFs more than pR46-27. It had 21 similar ORFs to pR46-27, accounting for 56.76% (21/37) of the ORFs in pR46-27 (Figure 1(a)). The pACN001-A plasmid is approximately 18.6 kb in length and shared nearly the same region with pR46-27 as pCTXM-2271. The ORFs of pR46-27 that matched those of the three plasmids included replication genes (bis and repB), a stabilization gene (staD), and partitioning genes (parA and parG). The floR-encoding region, flanked by a pair of IS91 sequences, had a transposon structure, and the virD2-floR-lysR gene cluster was conserved in three plasmids, excluding p160070-CTXM which does not harbor the floR gene (Figure 1(a)).

3.6. The Conjugative Plasmid pR46-270 Carries Multiple Resistance and Virulence Genes

Sequence annotation revealed that a conjugative system consisting of 24 tra genes, 5 trb genes, and a finO gene was encoded in a region approximately 36.5 kb in length in pR46-270 (Figure 1(b)). The conjugation experiment confirmed the transferability of pR46-270, facilitating the mobilization of pR46-27 (Table 5). All 26 resistance genes carried by pR46-270 were clustered in a region 47.5 kb in length (Figure 1(b), Figure 3). Interestingly, all three known sulfonamide genes (sul1 (orf00311 and orf00319), sul2 (orf00327), and sul3 (orf00284)) were found to be encoded in one K. pneumoniae plasmid for the first time in this study. Comparative genomic analysis demonstrated that the MDR region could be divided into six fragments, and all of the resistance genes were associated with the mobile genetic elements, including two class 1 integrons. The first fragment was a transposon mediated by Tn1721 (5,749 bp in length, between 177.11 and 171.36 kb) encoding two resistance genes, tetA and tetR, which showed >99% similarity to a fragment in pK1HV from a K. pneumoniae strain isolated from a healthy neonate in Vietnam. The second was a class 1 integron (11,770 bp in length, between 177.17 and 188.94 kb) carrying 5 resistance gene cassettes, dfrA12, aadA2, cmlA1, aadA1, and qacH, that showed >99% nucleotide sequence identity to a fragment in pCERC3 from an E. coli strain isolated from a healthy Australian adult. The third was a 13,942 bp fragment (between 189.68 and 203.63 kb) containing seven complete or truncated resistance genes (mphA, mrx, mphR, sul1, qnrB2, qacEΔ1, and arr-2) and three insertion sequences (IS, IS26, IS6100, and ISCR1). This fragment shared the highest identity (99%) with a sequence in plasmid pJIE137 from a K. pneumoniae strain isolated from a patient in Australia. The fourth fragment was a 7,634 bp sequence (between 204.11 and 211.74 kb) that could be divided into two parts. One part contained a truncated Tn6029 (containing repA and repC of the repABC operon) and three resistance genes (sul2 and strAB). The other part was a Tn4352 containing two copies of IS26 and aphA1a. Both parts were similar (>99% identity) to a region in plasmid p123 from an E. coli strain isolated from a canine. The fifth fragment, containing IS26, aac3, and tmrB, was a 1,934 bp sequence (between 211.97 and 223.91 kb) that showed >99% similarity to a fragment in plasmid p12181-KPC from a K. pneumoniae strain isolated from a patient in China. The last fragment was a 3,374 bp sequence (between 214.38 and 217.75 kb) composed of a class 1 integron with three resistance gene cassettes (aacA4, blaOXA-1, and catB3) that showed >99% sequence similarity to a fragment in plasmid pAUSMDU8141-1 from a Citrobacter farmeri strain isolated from a patient.

Figure 3: Schematic representation of the MDR region of plasmid pR46-270. Red arrows indicate resistance genes, yellow arrows indicate mobile genetic elements, and gray arrows indicate genes with other functions. Truncated genes are indicated by the “Δ” symbol.

In addition to the multiantibiotic resistance genes, four clusters (operons) of heavy metal resistance genes, including those for silver (silA, silB, silC, silE, silF, silR, and silS), copper (copA, copB, copC, copD, pcoE, copG, pcoR, and pcoS), arsenic (arsR, arsD, arsA, arsB, and arsC), and mercury (merA, merC, merP, merE, merD, merT, and merR), were identified in the plasmid. Moreover, several gene clusters or genes related to virulence were identified in the plasmid, such as two iron uptake systems (fecABCDEIR and iucABCDiutA), the genes of the lac operon (lacIYZ), and the fimbrial protein gene fimK.

4. Discussion

In this work, we analyzed the MDR strain K. pneumoniae R46 of ST37 isolated from a rabbit in south China. According to MLST typing, there are more than 3,000 ST types of K. pneumoniae identified worldwide (http://bigsdb.pasteur.fr), of which ST37 strains have mostly been reported to be MDR clones. The K. pneumoniae ST37 bacterium has largely been isolated from human clinical samples [27, 28]. Recently, ST37 was reported to have been isolated from animals, including companion animals [29] and chickens [30]. K. pneumoniae R46 carried 34 resistance genes and showed resistance to a variety of antibiotics, including β-lactams, amphenicols, aminoglycosides, and quinolones. In addition, this strain harbored three resistance plasmids, and 76.47% (26/34) of its resistance genes were encoded in a conjugative plasmid, all of the resistance genes were associated with mobile genetic elements. It demonstrated that horizontal gene transfer played an important role in the development of the multidrug resistance of K. pneumoniae R46.

K. pneumoniae R46 showed a high resistance level to florfenicol (). Of the eleven known florfenicol resistance genes, floR, floRv, floSt, optrA, pexA, fexB, fexA, cfr, cfr(B), cfr(C), and estDL136, only floR was found in the pR46-27 plasmid. Interestingly, an mdfA variant, named mdfA2, was encoded in the K. pneumoniae R46 chromosome and conferred resistance to florfenicol. MdfA has been reported to confer resistance to a variety of antibiotics, such as chloramphenicol [31], ciprofloxacin [32], and fluoroquinolone [33], while its resistance to florfenicol has never been documented.

To date, all mdfA variants have been found to be encoded in chromosomes, and no mobile genetic elements have been identified in their neighboring regions, indicating that mdfA genes are relatively conserved and might not be carried and transferred by mobile genetic elements. Our comparative genomic analysis demonstrated that the sequences next to the mdfA genes in different bacterial genera had different structures, while sequences of the same species or genera, to some extent, showed conserved structures. These observations indicated that sequence rearrangement rarely occurs in the mdfA-encoding region.

In addition to antibiotic agents used for human and animal therapy, a large number of other chemical substances with antibacterial activities, such as heavy metals and detergents, are used in human health care and agricultural practices. In addition to encoding 26 antibiotic resistance genes, the pR46-270 plasmid harbored 4 clusters of heavy metal resistance genes (including silver, copper, arsenic, and mercury, a total of 29 genes) and a number of other metal (iron) and virulence-related genes or gene clusters. Associations between heavy metal and antibiotic resistance have been reported previously [34]. Resistance to metals and to antibiotics might be subject to coselection [34]. Metal contamination represents a long-standing, widespread, and recalcitrant selective pressure with both environmental and clinical importance, and it potentially contributes to the maintenance and spread of antibiotic resistance factors.

5. Conclusion

Whole-genome sequencing of an animal-derived K. pneumoniae isolate harboring three resistance plasmids revealed that in addition to floR encoded in the pR46-27 plasmid, a chromosome-encoded mdfA variant, named mdfA2, also conferred florfenicol resistance. The genome encoded 34 drug resistance genes, of which 26 were encoded in the conjugative plasmid pR46-270. These 26 resistance genes were all associated with mobile genetic elements that shared homologous sequences in different bacterial genera of different origins. Four clusters of heavy metal (silver, copper, arsenic, and mercury) resistance genes and a number of virulence-related genes or gene clusters were also identified in the pR46-270 plasmid. The results of this work demonstrated that the plasmids with multidrug resistance genes were present in animal-derived bacteria and more florfenicol resistance genes such as mdfA2 could be present in bacterial populations. Active surveillance efforts are imperative for monitoring the prevalence of antibiotic-resistant bacteria, especially florfenicol-resistant strains, in animals around the world.

Abbreviations

BLAST:
Basic local alignment search tool.

Data Availability

The data related to this paper are deposited in the NCBI database.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

Junwan Lu, Yuanyuan Ying, Min Yin, Changrui Qian, Qianqian Chen, Kai Shen, and Chongyang Wu collected the strains and performed the experiments. Fei Wu, Yuanyuan Ying, and Cong Chen analyzed the experimental results. Fei Wu, Teng Xu, and Licheng Zhu performed the bioinformatic analysis. Fei Wu, Yuanyuan Ying, Qiyu Bao, and Junwan Lu co-led the writing of the manuscript. Junwan Lu, Kewei Li, Teng Xu, and Qiyu Bao designed the work. All authors read and approved the final manuscript. Fei Wu and Yuanyuan Ying contributed equally to this work.

Acknowledgments

The authors would like to acknowledge all of the study participants and individuals who contributed to the study. This work was funded by grants from the Natural Science Foundation of Zhejiang Province (grant numbers LQ17H190001 and LY14C060005), the National Natural Science Foundation of China (grant numbers 31500109 and 80215049), and the Science & Technology Project of Inner Mongolia Autonomous Region, China (grant number: 201802125).

Supplementary Materials

Table S1: 980 sequences containing mdfA gene clustered by CD-HIT with 90% identity and 85% coverage. (Supplementary Materials)

References

G. S. Davis and L. B. Price, “Recent research examining links among Klebsiella pneumoniae from food, food animals, and human extraintestinal infections,” Current Environmental Health Reports, vol. 3, no. 2, pp. 128–135, 2016.

S. Koovapra, S. Bandyopadhyay, G. Das et al., “Molecular signature of extended spectrum β-lactamase producing Klebsiella pneumoniae isolated from bovine milk in eastern and north-eastern India,” Infection, Genetics and Evolution, vol. 44, pp. 395–402, 2016.

S. Brisse and E. . Duijkeren, “Identification and antimicrobial susceptibility of 100 Klebsiella animal clinical isolates,” Veterinary Microbiology, vol. 105, no. 3-4, pp. 307–312, 2005.

Y. Guo, H. Zhou, L. Qin et al., “Frequency, antimicrobial resistance and genetic diversity of Klebsiella pneumoniae in food samples,” PLoS One, vol. 11, no. 4, article e0153561, 2016.

N. Kieffer, M. Aires-de-Sousa, P. Nordmann, and L. Poirel, “High rate of MCR-1-producing Escherichia coli and Klebsiella pneumoniae among pigs, Portugal,” Emerging Infectious Diseases, vol. 23, no. 12, pp. 2023–2029, 2017.

M. Saidani, L. Messadi, A. Soudani et al., “Epidemiology, antimicrobial resistance, and extended-spectrum beta-lactamase-producing Enterobacteriaceae in clinical bovine mastitis in Tunisia,” Microbial Drug Resistance, vol. 24, no. 8, pp. 1242–1248, 2018.

H. Wu, M. Wang, Y. Liu et al., “Characterization of antimicrobial resistance in Klebsiella species isolated from chicken broilers,” International Journal of Food Microbiology, vol. 232, pp. 95–102, 2016.

Y. Liu, Y. Yang, Y. Chen, and Z. Xia, “Antimicrobial resistance profiles and genotypes of extended-spectrum β-lactamase- and AmpC β-lactamase-producing Klebsiella pneumoniae isolated from dogs in Beijing, China,” Journal of Global Antimicrobial Resistance, vol. 10, pp. 219–222, 2017.

W. Bi, H. Liu, R. A. Dunstan et al., “Extensively drug-resistant Klebsiella pneumoniae causing nosocomial bloodstream infections in China: molecular investigation of antibiotic resistance determinants, informing therapy, and clinical outcomes,” Frontiers in Microbiology, vol. 8, p. 1230, 2017.

H. Guducuoglu, N. C. Gursoy, Y. Yakupogullari et al., “Hospital outbreak of a colistin-resistant, NDM-1- and OXA-48-producing Klebsiella pneumoniae: high mortality from pandrug resistance,” Microbial Drug Resistance, vol. 24, no. 7, pp. 966–972, 2018.

S. R. Partridge, “Analysis of antibiotic resistance regions in Gram-negative bacteria,” FEMS Microbiology Reviews, vol. 35, no. 5, pp. 820–855, 2011.

C. J. Harmer and R. M. Hall, “IS26-mediated formation of transposons carrying antibiotic resistance genes,” mSphere, vol. 1, no. 2, 2016.

Y. Liu, Y. Gao, X. Liu et al., “Transposon insertion sequencing reveals T4SS as the major genetic trait for conjugation transfer of multi-drug resistance pEIB202 from Edwardsiella,” BMC Microbiology, vol. 17, no. 1, p. 112, 2017.

M. Kaushik, S. Kumar, R. K. Kapoor, J. S. Virdi, and P. Gulati, “Integrons in Enterobacteriaceae: diversity, distribution and epidemiology,” International Journal of Antimicrobial Agents, vol. 51, no. 2, pp. 167–176, 2018.

D. Wibberg, I. P. Salto, F. G. Eikmeyer et al., “Complete genome sequencing of Acinetobacter baumannii strain K50 discloses the large conjugative plasmid pK50a encoding carbapenemase OXA-23 and extended-spectrum β-lactamase GES-11,” Antimicrobial Agents and Chemotherapy, vol. 62, no. 5, 2018.

L. Diancourt, V. Passet, J. Verhoef, P. A. D. Grimont, and S. Brisse, “Multilocus sequence typing of Klebsiella pneumoniae nosocomial isolates,” Journal of Clinical Microbiology, vol. 43, no. 8, pp. 4178–4182, 2005.

M. V. Larsen, S. Cosentino, S. Rasmussen et al., “Multilocus sequence typing of total-genome-sequenced bacteria,” Journal of Clinical Microbiology, vol. 50, no. 4, pp. 1355–1361, 2012.

S. Koren, B. P. Walenz, K. Berlin, J. R. Miller, N. H. Bergman, and A. M. Phillippy, “Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation,” Genome Research, vol. 27, no. 5, pp. 722–736, 2017.

H. Li, B. Handsaker, A. Wysoker et al., “The sequence alignment/map format and SAMtools,” Bioinformatics, vol. 25, no. 16, pp. 2078-2079, 2009.

G. A. Van der Auwera, M. O. Carneiro, C. Hartl et al., “From FastQ data to high-confidence variant calls: the Genome Analysis Toolkit best practices pipeline,” Current Protocols in Bioinformatics, vol. 43, pp. 11.10.1–11.10.33, 2013.

A. L. Delcher, K. A. Bratke, E. C. Powers, and S. L. Salzberg, “Identifying bacterial genes and endosymbiont DNA with Glimmer,” Bioinformatics, vol. 23, no. 6, pp. 673–679, 2007.

S. F. Altschul, W. Gish, W. Miller, E. W. Myers, and D. J. Lipman, “Basic local alignment search tool,” Journal of Molecular Biology, vol. 215, no. 3, pp. 403–410, 1990.

A. Petkau, M. Stuart-Edwards, P. Stothard, and G. Van Domselaar, “Interactive microbial genome visualization with GView,” Bioinformatics, vol. 26, no. 24, pp. 3125-3126, 2010.

K. Lagesen, P. Hallin, E. A. Rodland, H. H. Staerfeldt, T. Rognes, and D. W. Ussery, “RNAmmer: consistent and rapid annotation of ribosomal RNA genes,” Nucleic Acids Research, vol. 35, no. 9, pp. 3100–3108, 2007.

T. M. Lowe and P. P. Chan, “tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes,” Nucleic Acids Research, vol. 44, no. W1, pp. W54–W57, 2016.

L. Fu, B. Niu, Z. Zhu, S. Wu, and W. Li, “CD-HIT: accelerated for clustering the next-generation sequencing data,” Bioinformatics, vol. 28, no. 23, pp. 3150–3152, 2012.

C. G. Giske, I. Fröding, C. M. Hasan et al., “Diverse sequence types of Klebsiella pneumoniae contribute to the dissemination of blaNDM-1 in India, Sweden, and the United Kingdom,” Antimicrobial Agents and Chemotherapy, vol. 56, no. 5, pp. 2735–2738, 2012.

J. Yang, L. Ye, L. Guo et al., “A nosocomial outbreak of KPC-2-producing Klebsiella pneumoniae in a Chinese hospital: dissemination of ST11 and emergence of ST37, ST392 and ST395,” Clinical Microbiology and Infection, vol. 19, no. 11, pp. E509–E515, 2013.

J. Xia, L. X. Fang, K. Cheng et al., “Clonal spread of 16S rRNA methyltransferase-producing Klebsiella pneumoniae ST37 with high prevalence of ESBLs from companion animals in China,” Frontiers in Microbiology, vol. 08, p. 529, 2017.

X. Wang, Y. Liu, X. Qi et al., “Molecular epidemiology of colistin-resistant Enterobacteriaceae in inpatient and avian isolates from China: high prevalence of mcr-negative Klebsiella pneumoniae,” International Journal of Antimicrobial Agents, vol. 50, no. 4, pp. 536–541, 2017.

D. Wang, H. Liang, J. Chen, Y. Mou, and Y. Qi, “Structural and environmental features of novel mdfA variant and mdfA genes in recombinant regions of Escherichia coli,” Microbial Drug Resistance, vol. 20, no. 5, pp. 392–398, 2014.

M. C. Swick, S. K. Morgan-Linnell, K. M. Carlson, and L. Zechiedrich, “Expression of multidrug efflux pump genes acrAB-tolC, mdfA, and norE in Escherichia coli clinical isolates as a function of fluoroquinolone and multidrug resistance,” Antimicrobial Agents and Chemotherapy, vol. 55, no. 2, pp. 921–924, 2011.

J. Y. Kim, S. M. Jeon, H. Kim, M. S. Park, and S. H. Kim, “A contribution of MdfA to resistance to fluoroquinolones in Shigella flexneri,” Osong Public Health and Research Perspectives, vol. 2, no. 3, pp. 216-217, 2011.

C. Baker-Austin, M. S. Wright, R. Stepanauskas, and J. V. McArthur, “Co-selection of antibiotic and metal resistance,” Trends in Microbiology, vol. 14, no. 4, pp. 176–182, 2006.