The [beta]-Lactamase Gene Profile and a Plasmid-Carrying Multiple Heavy Metal Resistance Genes of Enterobacter cloacae.
Bacteria of the Enterobacter cloacae complex (ECC), which comprises six species, namely, E. cloacae, E. asburiae, E. hormaechei, E. kobei, E. ludwigii, and E. nimipressuralis , are widely distributed in nature. As pathogens, ECC species are highly adapted to the environment and are able to contaminate hospital medical devices. Currently, E. cloacae and E. hormaechei are most frequently isolated from human clinical specimens, and E. cloacae is among the Enterobacter sp. that have most commonly caused nosocomial infections in the last decade . Furthermore, E. cloacae has assumed clinical importance and has emerged as a major human pathogen; it accounts for up to 5% of hospital-acquired bacteremia cases, 5% of nosocomial pneumonia cases, 4% of nosocomial urinary tract infections, and 10% of postsurgical peritonitis cases .
Owing to the low-level but inducible expression of a chromosomal ampC gene encoding the AmpC [beta]-lactamase, E. cloacae is intrinsically resistant to ampicillin, amoxicillin-clavulanate, and first-generation cephalosporins . Generally, the resistance of E. cloacae to third-generation cephalosporins is caused by its overproduction of the AmpC [beta]-lactamases when the production of this cephalosporinase is inducible in the presence of strong [beta]-lactam antibiotics (cefoxitin and imipenem); thus, treatment with third-generation cephalosporins may promote the development of AmpC-overproducing mutants. AmpC-producing organisms become resistant to almost all [beta]-lactam antibiotics, with the exception of cefepime, cefpirome, and carbapenems. Most chromosomal ampC genes are inducible in the presence of certain agents such as cefoxitin and imipenem. Inducible AmpC expression is regulated by AmpR in the presence of two other gene products, namely, AmpD and AmpG. The regulation of AmpC production has been historically understood to require three proteins: AmpG, a plasma membrane-bound permease; AmpD, a cytosolic peptidoglycan-recycling amidase; and AmpR, the transcriptional regulator of AmpC. Derepression has been associated previously with structural defects within the ampD gene or with decreased ampD expression. Derepression represents the inability of AmpR to keep AmpC expression at constitutively low wild-type levels . As a result, the AmpC enzyme confers resistance to third-generation cephalosporins and is not inhibited by common [beta]-lactamase inhibitors. However, fourth-generation cephalosporins still retain activity against most Enterobacteriaceae strains.
In addition to therapeutic antibiotic agents, 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. Recently, concerns have been raised regarding coselection for antibiotic resistance among bacteria exposed to disinfectants and heavy metals (particularly copper, zinc, and mercury) used in some livestock species as growth promoters and therapeutic agents . Enterobacteriaceae (including E. coli, K. pneumoniae, and E. cloacae) are highly adept at acquiring resistance genes to all disinfectants, heavy metals, and antibiotics through horizontal gene transfer between different bacteria within the environment; such genes include extended-spectrum beta-lactamases (ESBLs), copper and arsenic resistance systems (the pco and ars operons), and enzymes that hydrolyze cephalosporins (AmpC enzymes) [7,8]. Many gram-negative organisms (such as E. coli, E. cloacae, and K. pneumoniae) encode broad-substrate efflux pumps [6, 7, 9], and a variety of multidrug pumps that have activity against disinfectants are similarly encoded by some gram-positive organisms, including Staphylococcus aureus [9, 10]. Besides the efflux pumps, other mechanisms such as chemisorption facilitating cadmium to bind to the bacterial surface also played a role in the heavy metal (cadmium) resistance . Alternatively, in both gram-positive and gram-negative bacteria, mechanisms of acquired resistance to disinfectants may be associated with efflux pump-encoding genes introduced on mobile genetic elements or, in gram-negative bacteria, with mutations causing the constitutive overexpression of efflux pumps. Compared with antibiotics and disinfectants, heavy metals (copper, arsenic, zinc, and mercury) are very persistent in the environment and may accumulate in soil, water, and sediments from agricultural practices as well as from other sources such as aquacultural and industrial effluents . Like the mechanisms of resistance to disinfectants, efflux pumps can expel heavy metal ions; such pumps include elements of the czc system, which encodes a pump for zinc, cobalt, and cadmium, and pcoA, which is an element of a copper extrusion system (W. ).
In this study, we identified [beta]-lactamase genes in 212 clinical
E. cloacae isolates and sequenced a [bla.sub.MIR] gene-carrying strain. Molecular analyses were performed to analyze the function of the resistance genes, and a comparative genomics analysis was conducted to elucidate the potential horizontal gene transfer patterns of genes related to both antibiotic and heavy metal resistance between bacteria of different species or genera. Our analysis revealed the distinct structure of a large plasmid carrying multiple clusters of heavy metal resistance genes that, to our knowledge, have not been described previously in E. cloacae.
2. Materials and Methods
2.1. Bacterial Strain Collection, Genomic DNA Extraction, and High-Throughput Sequencing. A total of 212 nonduplicate clinically significant E. cloacae strains were isolated from the First Affiliated Hospital of Wenzhou Medical University (Zhejiang, China) between 2008 and 2012. These isolates were identified by a VITEK 60 microbial autoanalyzer (bio-Merieux, Lyon, France). Bacteria and plasmids used in this study are listed in Table 1. Among the isolates, 31, 36, 43, 32, and 70 strains were isolated in the years 2008, 2009, 2010, 2011, and 2012, respectively. All strains were resistant to a minimum of one or two antibiotics. For pooled genomic DNA sequencing, each clinical strain was incubated in 5 mL of Luria-Bertani (LB) broth at 37[degrees]C for approximately 16 h to obtain a concentration equivalent to an optimum optical density ([OD.sub.600] = 15 [+ or -] 0.2). The cultures were pooled, and genomic DNA was extracted from 100 mL of the mixed bacteria using an AxyPrep Bacterial Genomic DNA Miniprep Kit (Axygen Scientific, Union City, CA, USA). The pooled genomic DNA was sequenced with a HiSeq 2500 DNA sequencer at Annoroad Gene Technology Co. Ltd. (Beijing, China). Reads derived from the HiSeq 2500 sequencing were initially assembled de novo with SOAPdenovo software to obtain contigs of the genome sequences. Glimmer software (http://ccb.jhu.edu/software/glimmer) was used to predict protein-coding genes with potential open reading frames (ORF) > 150 bp in length. BLASTX (https://blast.ncbi.nlm .nih.gov) was used to annotate the predicted protein-coding genes against the nonredundant protein database with an e-value threshold of 1E-5.
2.2. Collection of Reference Sequences for Resistance Genes and Mapping of Sequencing Reads to the Reference Genes. The nucleotide sequences of all [beta]-lactamase genes, including those encoding for Ambler class A, B, C, and D 3-lactamases, were obtained from the GenBank nucleotide database, and the high-throughput sequencing reads were mapped to the reference sequences of the [beta]-lactamase genes as previously described . The relative abundance (sequencing depth) of a certain gene was calculated as the cumulative nucleotide length of the mapped short reads on the gene divided by the gene size.
2.3. Screening of [beta]-Lactamase Gene-Positive Strains and Cloning and Phylogenetic Analysis of [bla.sub.MIR] Genes. The E. cloacae strains were screened by PCR amplification for the presence of [beta]-lactamase genes as previously described [15-18]. The primers used for the cloning of complete ORFs contained a pair of flanking restriction endonuclease adapters (BamHI for the forward primers and HindlU for the reverse primers) and were designed using the Primer Premier 5.0 software package; the primers are shown in Table 2. Genomic DNA was extracted from each of the 212 clinical E. cloacae isolates using the AxyPrep Bacterial Genomic DNA Miniprep kit (Axygen Scientific, Union City, CA, USA) and was used as the template for PCR amplification. Positive amplification products were verified by sequencing with an ABI 3730 automated sequencer (Shanghai Sunny Biotechnology Co. Ltd., Shanghai, China), and the sequencing results were compared with Basic Local Alignment Search Tool (BLAST) algorithms (https://blast .ncbi.nlm.nih.gov/blast.cgi). The amplicons of the bla ORFs were digested with the corresponding restriction endonucleases and ligated into the pET-28a vector with a T4 DNA ligase cloning kit (Takara Bio Inc., Dalian, China). The recombinant plasmid was transformed into competent E. coli BL21 cells using the calcium chloride method and cultured on LB agar plates supplemented with ampicillin (200 [micro]g/mL), and the cloned ORFs were confirmed by sequencing. For the phylogenetic analysis, all [bla.sub.MIR] gene sequences were collected from the NCBI nucleotide database using [bla.sub.MIR] as the key search term. A phylogenetic tree of the MIR amino acid sequences from both the database and this work was reconstructed by the maximum likelihood method, and the resulting trees were analyzed with bootstrap values of 100 replicates using MEGA 6.0 (https://www.megasoftware.net/).
2.4. Antimicrobial Susceptibility Testing. Antimicrobial susceptibility testing was conducted for all tested antibiotics by the agar dilution method, and the minimum inhibitory concentrations (MICs) were interpreted based on the Clinical and Laboratory Standards Institute (CLSI) breakpoint criteria (CLSI, 2018) (https://clsi.org/standards/products/ packages/m02-m07-m100-package/). Strain ATCC 25922 was used as the quality control strain. The 14 antibiotics (or antibiotics combination) used in this work were cephamycins (cefoxitin and cefminox), semisynthetic broad-spectrum penicillins (ampicillin and piperacillin), a first-generation cephalosporin (cefazolin), third-generation cephalosporins (ceftazidime, cefoperazone, cefotaxime, and ceftriaxone), a fourth-generation cephalosporin (cefoselis), a monobactam (aztreonam), an aminoglycoside (kanamycin), and combinations of antibiotics with [beta]-lactamase inhibitors (piperacillin/ tazobactam and imipenem/cilastatin sodium hydrate) (Table 3).
2.5. Whole Genome Sequencing (WGS) of Y546 and Comparative Genomics Analysis. The E. cloacae Y546 genomic DNA was extracted with the AxyPrep Bacterial Genomic DNA Miniprep Kit (Axygen Scientific, Union City, CA, USA) and sequenced with Illumina HiSeq 2500 and Pacific Biosciences sequencers at Annoroad Gene Technology Co. Ltd. (Beijing, China). Reads derived from HiSeq 2500 sequencing were initially assembled de novo with SOAPdenovo software to obtain genome sequence contigs. Reads of approximately 10-20 kb in length from the Pacific Biosciences sequencing were mapped onto the primary assembly for contig scaffolding. Gaps were filled either by remapping the short reads from the HiSeq 2500 sequencing or by sequencing the gap PCR product. Potential ORFs were predicted and annotated using Glimmer3 (http://www.cbcb .umd.edu/software/glimmer) and BASys , respectively. GC view software was used to construct the basic genomic features. Annotations were revised using UniProt (http://www. uniprot.org/) and BLAST (https://blast.ncbi.nlm.nih.gov/ blast.cgi). Plasmid replicons and plasmid incompatibility groups were predicted using Plasmid Finder (https://cge.cbs. dtu.dk//services/PlasmidFinder/). Furthermore, the multilocus sequence typing (MLST) database for E. cloacae (https://pubmLst.org/ecloacae/) was used to determine the sequence type of E. cloacae Y546.
The plasmid and chromosomal genomic sequences used in this study for the whole genome comparative analysis were downloaded from the NCBI database (http://www.ncbi.nlm.nih.gov). The plasmid and chromosome sequences were selected based on the comparison of the whole genome sequence (pY546) against the sequences of plasmids and chromosomes available in the NCBI database; a cutoff value (max score) of approximately 8700 was defined. For the comparative analysis of the heavy metal gene cluster regions on the pY546 plasmid, the sequences containing corresponding gene clusters with sequence identities of [greater than or equal to] 80% with respect to those encoded on pY546 were obtained from the NCBI
nucleotide database by BLASTn algorithms. Multiple sequence alignments were performed using MAFFT . Comparisons of the nucleotide sequences were made using BLASTn. Insertion sequences were predicted using ISfinder . Orthologous groups of genes from plasmids or chromosomes were identified using BLASTp and Inparanoid . Additional bioinformatics software was written using Python (https://www.python.org/) and Biopython .
2.6. Pulsed-Field Gel Electrophoresis (PFGE). The clonal relatedness of MIR-producing E. cloacae isolates was assessed by XbaI (Takara Bio Inc., China) PFGE. Briefly, genomic DNA fragments were resolved on a 1% agarose (SeaKem Gold Aga-rose, Lonza) gel at 14[degrees]C, and electrophoresis was conducted at 6 V/cm using a CHEF PFGE instrument (Bio-Rad, USA) at a pulse time gradient of 2.25-55.5 s and a total run time of 18 h. Salmonella enterica serovar H9812 was used as a control. Cluster analysis was performed using an unweighted pair-group method with arithmetic (UPGMA) means. Isolates were allocated into genetic similarity clusters using a similarity cutoff value of 80% .
2.7. Nucleotide Sequence Accession Numbers. The sequences of the chromosome and the plasmid of Y546 and the [bla.sub.MIR] genes in this work have been submitted to NCBI GenBank with accession numbers of CP032916 (Y546), CP032915 (pY546), MK033024 ([bla.sub.MIR-CG34]), MK033023 ([bla.sub.MIR-Y490]), MK033025 ([bla.sub.MIR-CG76]), MK033026 ([bla.sub.MIR-Y546]), MK033021 ([bla.sub.MIR-CG85]), and MK033022 ([bla.sub.MIR-Y482]), respectively.
3.1. Mapping and Screening Results for the [beta]-Lactamase Genes in the Sequenced Bacteria. A total of 75 [beta]-lactamase gene sequences were collected from the database (Supplementary Table S1). The pooled genomic DNA sequences of the 212 isolated strains generated approximately 34.1 gigabases. All reads ranged from 100 to 110 nucleotides in length. The mapping of the sequencing reads onto the reference sequences yielded the identification of resistance genes; in addition, the quantity of mapped reads on a specific reference could suggest the relative abundance of the reads in the sequenced samples. This analysis showed that the samples contained a total of 12 hits related to [beta]-lactamase resistance genes. The most abundant gene was [bla.sub.TEM], which had a sequence depth of 466.15 (Supplementary Table S2). The other genotypes with a greater abundance were [bla.sub.SHV], [bla.sub.DHA], and [bla.sub.CTX-M] (especially the [bla.sub.CTX-M-9] and [bla.sub.CTX-M-1] subtypes), while the genotypes with a lower abundance were [bla.sub.Z], [bla.sub.VEB], [bla.sub.KLUC], [bla.sub.MIR], [bla.sub.SFO], [bla.sub.AZECL], [bla.sub.OXA], and [bla.sub.ACT]. The screening results for the [bla.sub.VEB], [bla.sub.Z], [bla.sub.AZECL], and [bla.sub.MIR] genotypes revealed that among the 212 strains, only 0.47% (1/212; Y412), 0.94% (2/212; CG3 and CG4), 0.94% (2/212; Y411 and CG90), and 2.83% (6/212; CG34, Y490, CG76, Y546, CG85, and Y482) carried [bla.sub.VEB], [bla.sub.Z], [bla.sub.AZECL], and [bla.sub.MIR], respectively.
3.2. Cloning and Functional Detection of the Resistance Genes. Fourteen antimicrobial agents were used to detect the MIC levels of the [bla.sub.MIR]-positive wild-type strains and the recombinant strains expressing the cloned [bla.sub.MIR] genes (pET28a-[bla.sub.MIR]/BL21). The MICs of the 14 antimicrobial agents against these strains are shown in Table 3. The MICs for all 6 [bla.sub.MIR]-positive wild-type strains (CG34, Y490, CG76, Y546, CG85, and Y482) demonstrated that they were resistant to 4 commonly used broad-spectrum beta-lactam antibiotics, including ampicillin, a first-generation cephalosporin (cefazolin), and cephamycins (cefmenoxime and cefoxitin), and E. cloacae CG76 displayed higher resistance levels than the other strains to all antibiotics tested. Like the host wild-type strains, the 6 recombinant strains expressing the cloned [bla.sub.MIR] genes (pET28a-[bla.sub.MIR]/BL21) were resistant to ampicillin, cefazolin, cefmenoxime, and cefoxitin. In addition, the recombinants with two other extended-spectrum [beta]-lactamase (ESBL) genes, namely, [bla.sub.SHV-12] and [bla.sub.CTX-M-9a], encoded on the Y546 chromosome displayed higher hydrolytic activity against these four [beta]-lactam antibiotics.
3.3. Clonal Relatedness of the [bla.sub.MIR]-Positive Strains and Genotypes and the Phylogenetic Tree Analysis of the [bla.sub.MIR] Genes. All 6 [bla.sub.MIR]-positive strains (Y490, Y482 CG34, CG76, CG85, and Y546) had distinct XbaI PFGE patterns (Figure 1), indicating that the prevalence of [bla.sub.MIR]-positive isolates was caused by disseminated gene transfer. Sequencing results showed that the [bla.sub.MIR] ORFs from strains CG85 and Y482 belonged to the [bla.sub.MIR-17] genotype, while the ORFs from strains CG34, CG76, Y490, and Y546 matched [bla.sub.MIR-5], [bla.sub.MIR-21], [bla.sub.MIR-3], and [bla.sub.MIR-20], respectively. These ORFs had 99% amino acid similarity to their respective reference sequences. To further analyze the evolutionary relationship of the 6 [bla.sub.MIR] genes identified in this work with other [bla.sub.MIR] genes, we performed a multiple sequence alignment on a total of 30 [bla.sub.MIR] variants including the 6 [bla.sub.MIR] genes identified in this study. The multiple-sequence alignment identified the Pro380Leu variant in [bla.sub.MIR-Y546], [bla.sub.MIR-Y482], and [bla.sub.MIR-Y490] and the Ala381Gln variant in [bla.sub.MIR-CG76], [bla.sub.MIR-CG85], [bla.sub.MIR-Y546], [bla.sub.MIR-Y482], and [bla.sub.MIR-Y490]. The Asn206His variant was identified only in [bla.sub.MIR-CG34] (Figure 2). The phylogenetic analysis (Figure 3) showed that with the exception of 2 sequences from CG85 and Y482 that had the same amino acid sequence and were located in the same branch, the 6 MIR proteins were located in unique branches.
3.4. General Features of the Y546 Genome. The genome of E. cloacae strain Y546 consists of a chromosome and a plasmid (pY546); the general features of the Y546 genome are shown in Table 4. The chromosome is 4.78 Mb in length, harbors 4312 ORFs, and has an average GC content of 56.02%. In addition to an ampC gene [bla.sub.MIR], the chromosome also encodes two other extended-spectrum [beta]-lactamase (ESBL) genes, namely, [bla.sub.CTX-M-9] and [bla.sub.SHV-12]. MLST determined that E. cloacae strain Y546 contains the leuS-90, rpoB-20, gyrB-127, dnaA-120, fusA-25, and rplB-12 alleles and belongs to the sequence type ST466. The pY546 plasmid, an IncHI1B plasmid, is 208,740 bp in length and encodes 232 ORFs (Supplementary Table S3), of which 56.89% (132/232) encode proteins with known functions, and it contains a number of accessory modules identified at different sites within the backbone. pY546 contained an incomplete copper resistance operon (pcoBCDE/cusRS, ORF91-96) and several clusters of genes related to resistance to other metals, including arsenic (arsABCDR, ORF100-105), tetrathionate (ttrABCDRS, ORF142-146) and tellurite (terCDEF, ORF153-156). This plasmid also encodes numerous metallic ion metabolism and transfer proteins, such as a potassium transporter (Kef, ORF109), a fluoride ion transporter (CrcB, ORF113), divalent cation transporters (ORF111 and ORF112), and voltage-gated chloride channel proteins (ORF120 and ORF121). In addition, TonB-dependent receptors (TBDRs, ORF12-14), which are involved in the uptake of essential nutrients, are identified in pY546.
3.5. Comparative Genomic Analysis of the pY546 Plasmid Genome. A total of six sequences having greater than 40% nucleotide sequence identity to the pY546 sequence were retrieved from GenBank. Four of these were plasmid sequences, namely, pKPN-332 of K. pneumoniae strain KPNIH39 (49%, CP014763.1), pKPN-3967 of K. pneumoniae strain KPNIH49 (47%, CP026186.1), plasmid unnamed 1 of K. pneumoniae strain KSB2_1B (44%, CP024507.1), and pKPN-262 of K. pneumoniae subsp. pneumoniae KPNIH27 (44%, CP007734.1). The other two were chromosome sequences of E. coli S43 (41%, CP010237.1) and E. coli MEM (40%, CP012378.1) (Figure 4). Comparative genomics analysis showed that pY546 is approximately 80-100 kb smaller than any of the three named plasmids (pKPN-332, pKPN-3967, and pKPN-262). The sequences of 102 genes (43.4%, 102/235) on pY546 showed high similarity (>90%) with those on each of the three named plasmids. These plasmids shared a conserved backbone sequence with pY546; the backbone included the replication initiation gene (repA), stable maintenance genes (parAB or sopAB), DNA mismatch repair system genes (mutS), and so on. On the other hand, all these plasmids possessed their own variable regions, which mainly included heavy metal (cop, ars, and ter) resistance gene clusters and hypothetical genes. The tetrathionate resistance genes (ttrABCRS) encoded on pY546 were not present in the three named plasmids. The named plasmids also had multiple copies of common mobile elements, such as transposons and insertion elements (IS). In the two E. coli (S43 and MEM) chromosome sequences, the pco gene cluster (pcoABCDRSE) was detected in MEM but not in S43. Nevertheless, arsBCD, which is a part of the arsenic resistance determinants, was present in both S43 and MEM. The tetrathionate resistance genes were unique to pY546, while the chromosome of S43 carried the ter operon, but the MEM chromosome did not (Figure 4). The ter operon terC-DEF was detected and in the same orientation in the three plasmids pKPN-332, pKPN-3967, and pKPN-262.
3.6. Comparative Analysis of Copper and Arsenic Resistance Gene Regions on the pY546 Plasmid. Comparative analysis of an 8.7 kb fragment of pY546 encoding both copper (pco) and arsenic (ars) operons showed that the 5 sequence fragments with the highest similarity to that of pY546 were 4 fragments from the pKO_JKo3_1 plasmid of Klebsiella oxytoca JKo3 (100%, AP014952.1), the pKPN1705-1 plasmid of Klebsiella quasivariicola KPN1705 (100%, CP022824.1), the CSK29544_3p plasmid of Cronobacter sakazakii ATCC 29544 (100%, CP011050.1), and the pKPN-262 plasmid of Klebsiella pneumoniae subsp. pneumoniae KPNIH27 (82%, CP007734.1) and 1 fragment from the chromosome of Escherichia coli MEM (CP012378.1, 82%) (Figure 5). The fragment sharing the highest sequence identity with that of pY546 (from E. cloacae) was in the CSK29544_3p plasmid of Enterobacter sakazakii. All sequences except for pY546 contained the complete copper (pco) operon structure. pY546, however, contained an incomplete copper (pco) operon with a truncated pcoB ((h)pcopB) gene and without pcoA gene. Four (pY546, CSK29544_3p, pKO_Jko3_1, and pKPN1705-1) sequences contained the complete ars operon gene clusters. Moreover, the latter two plasmids (pKO_Jko3_1 and pKPN1705-1) contained an additional functional gene, namely, arsH, which encoded an organoarsenical oxidase (NADPH-dependent FMN reductase)  and conferred resistance to trivalent forms of organoarsenic compounds. Two of the plasmids (pKO_Jko3_1 and pKPN1705-1) were also identified to contain two copies of arsA and arsD with inverted orientations. The CSK29544_3p plasmid contained the same gene arrangement and content as pY546, but pY546 and CSK29544_3p contained fewer genes than the other two plasmids (pKO_Jko3_1 and pKPN1705-1). The pY546 and CSK29544_3p plasmids lacked arsH and contained only one copy each of arsAD, which was oriented oppositely in these two plasmids. On the other hand, the arsBCRH gene cluster was identified in pKPN-262, while the MEM chromosome contained arsBCR but lacked arsH.
The production of [beta]-lactamases is the predominant [beta]-lactam resistance mechanism in gram-negative bacteria. The Ambler molecular classification categorizes these [beta]-lactamases into four enzyme classes, namely, A, B, C, and D. Class A, C, and D enzymes all possess an active site serine, whereas class B [beta]-lactamases are metalloenzymes with a [Zn.sup.2+] ion(s) in the active site . In this study, a total of 12 [beta]-lactamase-encoding genes, including 7 class A [beta]-lactamase genes ([bla.sub.SHV], [bla.sub.CTX-M], [bla.sub.Z], [bla.sub.VEB], [bla.sub.KLUC], [bla.sub.SFO], and [bla.sub.TEM]), 4 class C [beta]-lactamase genes ([bla.sub.MIR], [bla.sub.DHA], [bla.sub.ACT], and [bla.sub.AZECL-29]), and 1 class D [beta]-lactamase gene ([bla.sub.OXA]) were identified in 212 E. cloacae isolates from a teaching hospital in South China. Over the past years, a variety of metalloenzymes (NDM- and IMP-type) have been found in E. cloacae and have contributed to infectious outbreaks in China  and Japan . However, we did not detect any genes encoding class B metalloenzymes in these 212 E. cloacae isolates.
The class A enzymes are regarded as extended-spectrum [beta]-lactamases (ESBLs) that can hydrolyze extended-spectrum cephalosporins; they are inhibited by clavulanic acid and are spreading widely among Enterobacteriaceae. The CTX-M enzymes are replacing the SHV and TEM enzymes as the most prevalent type of ESBL in Enterobacteriaceae [29, 30]. Additional clinically relevant types of ESBLs include the VEB, PER, GES, TLA, IBC, SFO-1, BES-1, and BEL-1 types. The SFO-1 [beta]-lactamase was first reported in 1988 in a clinical E. cloacae isolate in Japan, and it confers resistance to third-generation cephalosporins (Matsumoto & Inoue, 1999); VEB was reported in China during an outbreak of infection caused by E. cloacae . No document has yet reported the identification of the [bla.sub.Z] gene that encodes class A enzymes in E. cloacae, and the [bla.sub.Z] gene has been found only once in Staphylococcus aureus . In this work, however, we isolated two E. cloacae strains (CG3 and CG4) that carried the [bla.sub.Z] gene.
In addition to the [bla.sub.MIR] and [bla.sub.DHA] genes identified in this work, genes encoding AmpC enzymes belonging to Ambler class C and Bush-Jacoby group 1 include [bla.sub.CMY], [bla.sub.FOX], [bla.sub.LAT], [bla.sub.ACT], [bla.sub.MOX], [bla.sub.ACC], and [bla.sub.BIL] and their derivatives (H. ). The emergence of AmpC-producing Enterobacter spp. has been observed globally in health care-associated settings and in the community . However, unlike most of the AmpC genes, [bla.sub.MIR] has been found only in some strains of several Enterobacter spp., mainly in strains of E. cloacae. A total of 24 [bla.sub.MIR] nucleotide sequences (between the [bla.sub.MIR-1] and [bla.sub.MIR-21] subtypes) are available in the NCBI nucleotide database; these sequences mainly came from the ECC, such as E. cloacae and E. aerogenes, as well as strains of K. pneumoniae and E. coli. In this work, 6 MIRs belonged 5 subtypes, including [bla.sub.MIR-3], [bla.sub.MIR-5], [bla.sub.MIR-17], [bla.sub.MIR-21] and [bla.sub.MIR-20]. Although the [bla.sub.MIR] gene has been primarily identified on bacterial chromosomes, it is also encoded on plasmids. The AmpC [beta]-lactamase [bla.sub.MIR-1] was first described in K. pneumoniae plasmids ; [bla.sub.MIR-1] confers resistance to penicillins and broad-spectrum cephalosporins, including cefoxitin and ceftibuten, but not to cefepime, cefpirome, meropenem, or imipenem. The resistance features of [bla.sub.MIR] genes in human and animal isolates were different from those of some plasmid-encoded AmpC-type [beta]-lactamase genes, such as [bla.sub.DHA] and [bla.sub.CMY], and have been reported worldwide to hydrolyze third-generation cephalosporins [35, 36]. The [bla.sub.MIR] gene identified in this work was encoded on the chromosome and showed high sequence identity with other homologous [bla.sub.MIR] genes found in other Enterobacteriaceae. Like the other previously reported MIRs, they showed resistance to ampicillin, cefazolin, cefmenoxime, and cefoxitin, but sensitive to fourth-generation cephalosporins (cefoselis) and monobactam (aztreonam).
To adapt to environmental changes, bacteria often harbor genes conferring resistance to toxic metal compounds; these genes include those encoding copper and arsenic ion transportation systems . Copper sulfate is a common feed supplement for pigs, chickens, and calves worldwide. The copper-binding operon system (PcoBCDE and CusR/S), which is known to transport copper-derived compounds out of the bacterial cell to balance the concentration of copper salts, was elaborated on the pRJ1004 plasmid of E. coli isolates from piggeries in which the animals were provided food supplemented with copper sulfate . Despite the identity of arsenic and arsenite compounds as high-toxicity compounds that are neither used in agriculture nor found in either the community or the hospital sector, the presence of arsenic resistance determinants on Enterobacteriaceae plasmids, especially the IncH-type plasmids, has been described before . Three prototypes of ars operons, including the three-, four-, and five-gene arsenic resistance determinants, namely, arsABC, arsABCD, and arsABCDR, respectively, have been well documented, although novel resistance mechanisms have also been described . The arsABCDR operon is related to resistance to arsenic-derived compounds, including arsine, arsenic, arsenite, and arsenate. Heavy metal resistance genes or gene clusters have been widely identified in different genera of both gram-positive and gram-negative bacteria and are encoded on both chromosomes and plasmids. In this work, on the plasmid pY546, we found four clusters of genes conferring resistance to heavy metals, such as arsenic (arsABCDR), tetrathionate (ttrABCDRS), and tellurite (terCDEF) as well as an incomplete copper resistance operon (pcoBCDE/cusRS). We must expect that bacteria have adapted to heavy metals with an increasing frequency.
In this work, through high-throughput sequencing, we identified twelve [beta]-lactamase genotypes in 212 clinical E. cloacae isolates; of these, [bla.sub.Z] has not yet been reported in E. cloacae. Furthermore, whole genome analysis of the [bla.sub.MIR]-carrying E. cloacae strain Y546 demonstrated that the strain harbored a large plasmid carrying a variety of gene clusters and genes, such as heavy metal resistance gene clusters (e.g., thepco, ars, ter, and ttr operons), conferring resistance to antimicrobials. Comparative genomics analysis showed that the sequences sharing the highest similarity to pY546 were plasmids from K. pneumoniae strains (44-49% similarity) and the chromosome of E. coli (40-41% similarity) and that the sequence fragments with the highest similarity to heavy metal resistance gene clusters on pY546 were from other plasmids and other chromosome sequences. The colocalization of antibiotic resistance genes and heavy metal resistance genes in the genomes of clinical pathogens, which may facilitate the persistence, coselection, and dissemination of these genes between different bacterial species or genera, is alarming and needs further surveillance.
The data used to support the findings of this study are included within the article.
Highlights. Twelve [beta]-lactamase genotypes, including [bla.sub.SHV], [bla.sub.TEM], [bla.sub.DHA], [bla.sub.CTX-M], [bla.sub.Z], [bla.sub.VEB], [bla.sub.KLUC], [bla.sub.MIR], [bla.sub.SFO], [bla.sub.AZECL-29], [bla.sub.OXA], and [bla.sub.ACT], were identified in 212 E. cloacae genomes. [bla.sub.Z] was found in E. cloacae for the first time. [bla.sub.MIR] was identified in six E. cloacae strains that were not clonally related. The gene clusters related to resistance to heavy metals (such as copper, arsenic, and tellurite) were identified to be encoded on the pY546 plasmid. The sequences with the highest similarity to pY546 were on plasmids from K. pneumoniae strains and on E. coli chromosomes; in addition, the sequence fragments with the highest similarity to the heavy metal resistance gene clusters on pY546 were from other sources.
Conflicts of Interest
The authors declare that there are no conflicts of interest in this work.
Chongyang Wu and Chaoqin Lin contributed equally to this work.
The work was funded by grants from the National Natural Science Foundation of China (31500109, 81501808, and 80215049) and the Science and Technology Foundation of Wenzhou City (Y20170205).
Table S1: the reference [beta]-lactamase gene sequences collected from GenBank. Table S2: the mapping result of [beta]-lactamase resistance genes in the E. cloacae pooled genomic sequences. Table S3: annotation result of the pY546 genome. (Supplementary Materials)
 M. L. Mezzatesta, F. Gona, and S. Stefani, "Enterobacter cloacae complex: clinical impact and emerging antibiotic resistance," Future Microbiology, vol. 7, no. 7, pp. 887-902, 2012.
 N. Petrosillo, M. Vranic-Ladavac, C. Feudi et al., "Spread of Enterobacter cloacae carrying bla NDM-1, bla CTX-M-15, bla SHV-12 and plasmid-mediated quinolone resistance genes in a surgical intensive care unit in Croatia," Journal of Global Antimicrobial Resistance, vol. 4, pp. 44-48, 2016.
 S. Wang, S.-Z. Xiao, F.-F. Gu et al., "Antimicrobial susceptibility and molecular epidemiology of clinical Enterobacter cloacae bloodstream isolates in Shanghai, China," Plos One, vol. 12, no. 12, article e0189713, 2017.
 K. Zhou, W. Yu, X. Cao et al., "Characterization of the population structure, drug resistance mechanisms and plasmids of the community-associated Enterobacter cloacae complex in China," The Journal of Antimicrobial Chemotherapy, vol. 73, no. 1, pp. 66-76, 2018.
 D. Balasubramanian, H. Kumari, and K. Mathee, "Pseudomonas aeruginosa AmpR: an acute-chronic switch regulator," Pathogens and Disease, vol. 73, no. 2, pp. 1-14, 2014.
 A. D. Wales and R. H. Davies, "Co-selection of resistance to antibiotics, biocides and heavy metals, and its relevance to foodborne pathogens," Antibiotics, vol. 4, no. 4, pp. 567-604, 2015.
 D. Deus, C. Krischek, Y. Pfeifer et al., "Comparative analysis of the susceptibility to biocides and heavy metals of extended-spectrum [beta]-lactamase-producing Escherichia coli isolates of human and avian origin, Germany," Diagnostic Microbiology and Infectious Disease, vol. 88, no. 1, pp. 88-92, 2017.
 C. J. Slipski, G. G. Zhanel, and D. C. Bay, "Biocide selective TolC-independent efflux pumps in Enterobacteriaceae," The Journal of Membrane Biology, vol. 251, no. 1, pp. 15-33, 2018.
 Y. Zhai, Z. He, Y. Kang et al., "Complete nucleotide sequence of pH11, an IncHI2 plasmid conferring multi-antibiotic resistance and multi-heavy metal resistance genes in a clinical Klebsiella pneumoniae isolate," Plasmid, vol. 86, pp. 26-31, 2016.
 M. A. Argudin, B. Lauzat, B. Kraushaar et al., "Heavy metal and disinfectant resistance genes among livestock-associated methicillin-resistant Staphylococcus aureus isolates," Veterinary Microbiology, vol. 191, pp. 88-95, 2016.
 S. Z. Abbas, M. Rafatullah, K. Hossain, N. Ismail, H. A. Tajarudin, and H. P. S. Abdul Khalil, "A review on mechanism and future perspectives of cadmium-resistant bacteria," International journal of Environmental Science and Technology, vol. 15, no. 1, pp. 243-262, 2018.
 A. Waseem, J. Arshad, F. Iqbal, A. Sajjad, Z. Mehmood, and G. Murtaza, "Pollution status of Pakistan: a retrospective review on heavy metal contamination of water, soil, and vegetables," BioMed Research International, vol. 2014, Article ID 813206, 29 pages, 2014.
 W. Deng, Y. Quan, S. Yang et al., "Antibiotic resistance in Salmonella from retail foods of animal origin and its association with disinfectant and heavy metal resistance," Microbial Drug Resistance, vol. 24, no. 6, pp. 782-791, 2018.
 T. Xu, J. Ying, X. Yao et al., "Identification and characterization of two novel bla(KLUC) resistance genes through large-scale resistance plasmids sequencing," PLoS One, vol. 7, no. 10, article e47197, 2012.
 J. P. R. Furlan and E. G. Stehling, "Detection of [beta]-lactamase encoding genes in feces, soil and water from a Brazilian pig farm," Environmental Monitoring and Assessment, vol. 190, no. 2, p. 76, 2018.
 L. Haubert, I. S. Kroning, M. A. Iglesias, and W. P. da Silva, "First report of the Staphylococcus aureus isolate from subclinical bovine mastitis in the South of Brazil harboring resistance gene dfrG and transposon family Tn916-1545," Microbial Pathogenesis, vol. 113, pp. 242-247, 2017.
 S. D. Lahiri, M. R. Johnstone, P. L. Ross, R. E. McLaughlin, N. B. Olivier, and R. A. Alm, "Avibactam and class C [beta]-lactamases: mechanism of inhibition, conservation of the binding pocket, and implications for resistance," Antimicrobial Agents and Chemotherapy, vol. 58, no. 10, pp. 5704-5713, 2014.
 L. X. Zhu, Z. W. Zhang, D. Liang et al., "Multiplex asymmetric PCR-based oligonucleotide microarray for detection of drug resistance genes containing single mutations in Enterobacteriaceae," Antimicrobial Agents and Chemotherapy, vol. 51, no. 10, pp. 3707-3713, 2007.
 G. H. Van Domselaar, P. Stothard, S. Shrivastava et al., "BASys: a web server for automated bacterial genome annotation," Nucleic Acids Research, vol. 33, Supplement 2, pp. W455-W459, 2005.
 K. Katoh and D. M. Standley, "MAFFT multiple sequence alignment software version 7: improvements in performance and usability," Molecular Biology and Evolution, vol. 30, no. 4, pp. 772-780, 2013.
 P. Siguier, J. Perochon, L. Lestrade, J. Mahillon, and M. Chandler, "ISfinder: the reference centre for bacterial insertion sequences," Nucleic Acids Research, vol. 34, no. 90001, pp. D32-D36, 2006.
 M. Remm, C. E. V. Storm, and E. L. L. Sonnhammer, "Automatic clustering of orthologs and in-paralogs from pairwise species comparisons," Journal of Molecular Biology, vol. 314, no. 5, pp. 1041-1052, 2001.
 P. J. A. Cock, T. Antao, J. T. Chang et al., "Biopython: freely available Python tools for computational molecular biology and bioinformatics," Bioinformatics, vol. 25, no. 11, pp. 1422-1423, 2009.
 S. Roussel, B. Felix, N. Vingadassalon et al., "Staphylococcus aureus strains associated with food poisoning outbreaks in France: comparison of different molecular typing methods, including MLVA," Frontiers in Microbiology, vol. 6, article 882, 2015.
 T. Tsubouchi and Y. Kaneko, "Draft genome sequence of the arsenic-resistant bacterium Brevundimonas denitrificans [TAR-002.sup.T]," Genome Announcements, vol. 5, no. 47, 2017.
 K. Bush and G. A. Jacoby, "Updated functional classification of beta-lactamases," Antimicrobial Agents and Chemotherapy, vol. 54, no. 3, pp. 969-976, 2010.
 C. Liu, S. Qin, H. Xu et al., "New Delhi metallo-[beta]-lactamase 1(NDM-1), the dominant carbapenemase detected in carbapenem-resistant Enterobacter cloacae from Henan Province, China," PLoS One, vol. 10, no. 8, article e0135044, 2015.
 K. Hayakawa, T. Miyoshi-Akiyama, T. Kirikae et al., "Molecular and epidemiological characterization of IMP-type metallo-[beta]-lactamase-producing Enterobacter cloacae in a large tertiary care hospital in Japan," Antimicrobial Agents and Chemotherapy, vol. 58, no. 6, pp. 3441-3450, 2014.
 E. Juhasz, L. Janvari, A. Toth, I. Damjanova, A. Nobilis, and K. Kristof, "Emergence of VIM-4- and SHV-12-producing Enterobacter cloacae in a neonatal intensive care unit," International Journal of Medical Microbiology, vol. 302, no. 6, pp. 257-260, 2012.
 T. Naas, G. Cuzon, A. L. Robinson et al., "Neonatal infections with multidrug-resistant ESBL-producing E. cloacae and K. pneumoniae in neonatal units of two different hospitals in Antananarivo, Madagascar," BMC Infectious Diseases, vol. 16, no. 1,p. 275,2016.
 X. Jiang, Y. Ni, Y. Jiang et al., "Outbreak of infection caused by Enterobacter cloacae producing the novel VEB-3 beta-lactamase in China," Journal of Clinical Microbiology, vol. 43, no. 2, pp. 826-831, 2005.
 H. Deng, H.-B. Si, S.-Y. Zeng et al., "Prevalence of extended-spectrum cephalosporin-resistant Escherichia coli in a farrowing farm: ST1121 clone harboring IncHI2 plasmid contributes to the dissemination of blaCMY-2," Frontiers in Microbiology, vol. 6, 2015.
 F. I. Mohd Khari, R. Karunakaran, R. Rosli, and S. Tee Tay, "Genotypic and phenotypic detection of AmpC [beta]-lactamases in Enterobacter spp. isolated from a teaching hospital in Malaysia," PLoS One, vol. 11, no. 3, article e0150643, 2016.
 G. A. Papanicolaou, A. A. Medeiros, and G. A. Jacoby, "Novel plasmid-mediated beta-lactamase (MIR-1) conferring resistance to oxyimino- and alpha-methoxy beta-lactams in clinical isolates of Klebsiella pneumoniae," Antimicrobial Agents and Chemotherapy, vol. 34, no. 11, pp. 2200-2209, 1990.
 Y. Maeyama, Y. Taniguchi, W. Hayashi et al., "Prevalence of ESBL/AmpC genes and specific clones among the third-generation cephalosporin-resistant Enterobacteriaceae from canine and feline clinical specimens in Japan," Veterinary Microbiology, vol. 216, pp. 183-189, 2018.
 Y. Xie, L. Tian, G. Li et al., "Emergence of the third-generation cephalosporin-resistant hypervirulent Klebsiella pneumoniae due to the acquisition of a self-transferable [bla.sub.DHA-1]-carrying plasmid by an ST23 strain," Virulence, vol. 9, no. 1, pp. 838-844, 2018.
 T. J. Lawton, G. E. Kenney, J. D. Hurley, and A. C. Rosenzweig, "The CopC family: structural and bioinformatic insights into a diverse group of periplasmic copper binding proteins," Biochemistry, vol. 55, no. 15, pp. 2278-2290, 2016.
 M. W. Gilmour, N. R. Thomson, M. Sanders, J. Parkhill, and D. E. Taylor, "The complete nucleotide sequence of the resistance plasmid R478: defining the backbone components of incompatibility group H conjugative plasmids through comparative genomics," Plasmid, vol. 52, no. 3, pp. 182-202, 2004.
 N. S. Chauhan, S. Nain, and R. Sharma, "Identification of arsenic resistance genes from marine sediment metagenome," Indian Journal of Microbiology, vol. 57, no. 3, pp. 299-306, 2017.
Chongyang Wu, (1,2) Chaoqin Lin, (1,2) Xinyi Zhu, (1,2) Hongmao Liu, (1,2) Wangxiao Zhou, (2) Junwan Lu, (3) Licheng Zhu, (3) Qiyu Bao (iD), (2) Cong Cheng (iD), (3) and Yunliang Hu (iD) (1,2)
(1) The Second Affiliated Hospital and Yuying Children's Hospital of Wenzhou Medical University, Wenzhou, Zhejiang 325000, China
(2) School of Laboratory Medicine and Life Sciences/Institute of Biomedical Informatics, Wenzhou Medical University, Wenzhou 325035, China
(3) College of Medicine and Health, Lishui University, Lishui 323000, China
Correspondence should be addressed to Cong Cheng; firstname.lastname@example.org and Yunliang Hu; email@example.com
Received 13 July 2018; Accepted 17 October 2018; Published 20 December 2018
Academic Editor: Joao Paulo Gomes
Caption: Figure 1: Pulsed-field gel electrophoresis (PFGE) analysis of the 6 [bla.sub.MIR]-positive E. cloacae isolates. PFGE result showed that all 6 [bla.sub.MIR]-positive isolates had a totally distinct PFGE pattern.
Caption: Figure 3: Phylogenetic tree of 30 MIR amino acid sequences. The "(a)" symbols indicate the MIRs of this study.
Caption: Figure 4: Complete sequence of the pY546 plasmid and comparative genomic analysis of the pY546 plasmid sequence with other sequences. The circles (from innermost to outermost) represent (i) the scale in kb; (ii) the cumulative GC skew; (iii) the GC content; (iv) the annotated coding sequences with selected genes indicated according to the gene function: heavy metal resistance genes in red arrows, transposase genes, IS elements in bottle-green arrows, and hypothetical proteins in dark gray arrows; (v) circles (from inside to outside) representing three homologous plasmids (pKPN-332, CP014763.1; pKPN-3967, CP026186.1; and pKPN-262, CP007734.1) and two chromosome fragments (S43, CP010237.1, and MEM, CP012378.1), respectively.
Caption: Figure 5: Comparative analysis of the copper and arsenic resistance gene clusters of pY546 and other sequences from different bacteria. The homologous gene clusters among plasmids pY546, CSK29544_3P (CP011050.1), pKO_JKo3_1 (AP014952.1), pKPN1705-1 (CP022824.1), pKPN-262 (CP007734.1), and MEM (CP012378.1) with the copper resistance gene clusters in bottle green, and the arsenic resistance gene cluster in orange. Annotated coding sequences are displayed as arrows. Coding sequences are colored based on their assigned gene functions.
Table 1: Bacteria and plasmids used in this study. Relevant Reference Strain characteristic(s) or source BL21 E. coli BL21 was used as  a host for the cloned [bla.sub.MIR] gene ATCC25922 E. coli ATCC25922 is FDA  clinical isolate CG34 [bla.sub.MIR]-producing This study isolate of E. cloacae CG34 CG76 [bla.sub.MIR]-producing This study isolate of E. cloacae CG76 CG85 [bla.sub.MIR]-producing This study isolate of E. cloacae CG85 Y546 [bla.sub.MIR]-producing This study isolate of E. cloacae Y546 Y482 [bla.sub.MIR]-producing This study isolate of E. cloacae Y482 Y490 [bla.sub.MIR]-producing This study isolate of E. cloacae Y490 pET28a/BL21 BL21 carrying expression This study vector pET28a, km (r) pET28a- BL21 carrying recombinant This study [bla.sub.MIR-CG34]/BL21 plasmid pET28a-[bla.sub.MIR-CG34] pET28a- BL21 carrying recombinant This study [bla.sub.MIR-Y490]/BL21 plasmid pET28a-[bla.sub.MIR-Y490] pET28a- BL21 carrying recombinant This study [bla.sub.MIR-CG76]/BL21 plasmid pET28a-[bla.sub.MIR-CG76] pET28a- BL21 carrying recombinant This study [bla.sub.MIR-Y546]/BL21 plasmid pET28a-[bla.sub.MIR-Y546] pET28a- BL21 carrying recombinant This study [bla.sub.MIR-CG85]/BL21 plasmid pET28a-[bla.sub.MIR-CG85] pET28a- BL21 carrying recombinant This study [bla.sub.MIR-Y482]/BL21 plasmid pET28a-[bla.sub.MIR-Y482] pET28a- BL21 carrying recombinant This study [bla.sub.CTX-M-9a-Y546]/ plasmid pET28a- BL21 [bla.sub.CTX-M_9a_Y546] pET28a- BL21 carrying recombinant This study [bla.sub.SHV-12-Y546]/ plasmid pET28a- BL21 [bla.sub.SHV-12-Y546] km: kanamycin; (r): resistant. Table 2: Primers used in the study for the detection of [beta]- lactamase-encoding genes. Primer Target(s) Sequence (5'-3') veb-sf [bla.sub.VEB] GATTGCTTTAGCCGTTTTGTC veb-sr ATCGGTTACTTCCTGTTGTTGTTTC z-sf [bla.sub.Z] ACAGTTCACATGCCAAAGAGT z-sr CTTACCGAAAGCAGCAGGTG mir-sf [bla.sub.MIR] GCCGCACCGATGTCCGAAAAA mir-sr GGTTTAAAGACCCGCGTCGTCATGG azecl-29-sf [bla.sub.AZECL-29] GTCTTTACGCTAACACCAGCATCGG azecl-29-sr TCAGCATTTCCCAGCCCAATC veb-ff [bla.sub.VEB] CGGGATCC#ATGAAAATCGTAAAAAGGATATTAT veb-fr CCCAAGCT#TTATTTATTCAAATAGTAATTCCACG z-ff [bla.sub.Z] CGGGATCC#ATGAAAAAGTTAATATTTTTAATTG z-fr CCCAAGCT#TTAAAATTCCTTCATTACACTCTTG mir-ff [bla.sub.MIR] CGGGATCC#ATGATGACAAAATCCCTAAGCTGTG mir-fr CCCAAGCT#TTACTGCAGCGCGTCGAGGATACGG azecl-29-ff [bla.sub.AZECL-29] CGGAATTC#ATGATGAAAAAAAACCTAAGCTGTG azecl-29-fr CCCAAGCT#TTACTGCAGCGCGTCGAGGATACG ctx-m-9a-ff [bla.sub.CTX-M-9a] CGGGATCC#ATGGTGACAAAGAGAGTGCAACGGA ctx-m-9a-fr CCAAGCTT#TTACAGCCCTTCGGCGATGATTCTC shv-12-ff [bla.sub.SHV-12] CGGGATCC#GTGGTTATGCGTTATATTCGCCTGT shv-12-fr CCAAGCTT#TTAGCGTTGCCAGTGCTCGATCAGC Annealing Amplicon temperature size Primer ([degrees]C) (bp) Reference veb-sf 50 452  veb-sr z-sf 50 479  z-sr mir-sf 50 545  mir-sr azecl-29-sf 50 381  azecl-29-sr veb-ff 50 918 This study veb-fr z-ff 50 864 This study z-fr mir-ff 66 1164 This study mir-fr azecl-29-ff 60 1164 This study azecl-29-fr ctx-m-9a-ff 60 876 This study ctx-m-9a-fr shv-12-ff 60 867 This study shv-12-fr Underlined sequences denote restriction endonuclease sites. sf: forward screening primer; sr: reverse screening primer; ff: forward full-length primer; fr: reverse full-length primer. Note: Sequences denote restriction endonuclease sites indicated with #. Table 3: The MIC values of antibacterial drugs for the strains ([micro]g/mL). Strains AMP KAN CAZ CTX ATCC25922 4 2 0.25 0.0625 BL21 0.5 4 0.0625 0.0156 pET-28a/BL21 * 0.5 32 0.0625 0.0156 CG34 256 0.5 0.125 0.125 Y490 256 2 0.5 0.25 CG76 >1024 8 128 128 Y546 512 2 2 1 CG85 >1024 4 1 1 Y482 64 2 0.25 0.125 pET-28a-[bla.sub.MIR] (CG34)/BL21 * 128 64 4 4 pET-28a-[bla.sub.MIR] (Y490)/BL21 * 128 64 4 4 pET-28a-[bla.sub.MIR] (CG76)/BL21 * 128 64 4 4 pET-28a-[bla.sub.MIR] (Y546)/BL21 * 128 64 4 4 pET-28a-[bla.sub.MIR] (CG85)/BL21 * 128 64 4 4 pET-28a-[bla.sub.MIR] (Y482)/BL21 * 128 64 4 4 pET-28a-[bla.sub.CTX-M-9a]/BL21 * 8 64 0.25 4 pET-28a-[bla.sub.SHV-12]/BL21 * 8 64 2 4 Strains CPZ CMN CEF CFZ ATCC25922 0.125 1 0.125 2 BL21 0.0156 1 0.0313 2 pET-28a/BL21 * 0.0156 1 0.0313 2 CG34 0.0625 1024 0.0625 512 Y490 1 >1024 0.125 512 CG76 64 1024 16 >1024 Y546 2 >1024 0.25 1024 CG85 1 >1204 0.25 1024 Y482 0.25 1024 256 512 pET-28a-[bla.sub.MIR] (CG34)/BL21 * 1 32 0.125 256 pET-28a-[bla.sub.MIR] (Y490)/BL21 * 0.5 32 0.0625 256 pET-28a-[bla.sub.MIR] (CG76)/BL21 * 1 32 0.0313 256 pET-28a-[bla.sub.MIR] (Y546)/BL21 * 1 32 0.0625 256 pET-28a-[bla.sub.MIR] (CG85)/BL21 * 1 32 0.125 256 pET-28a-[bla.sub.MIR] (Y482)/BL21 * 1 32 0.125 256 pET-28a-[bla.sub.CTX-M-9a]/BL21 * 1 4 0.03 64 pET-28a-[bla.sub.SHV-12]/BL21 * 1 4 0.25 64 Strains FOX CRO ATM PRL ATCC25922 4 0.0313 0.125 2 BL21 2 0.0156 0.0156 0.5 pET-28a/BL21 * 2 0.0156 0.0156 0.5 CG34 1024 0.0313 0.25 2 Y490 1024 0.5 0.0313 2 CG76 1024 512 32 512 Y546 1024 2 1 8 CG85 >1024 1 0.25 4 Y482 0.0625 0 0.125 2 pET-28a-[bla.sub.MIR] (CG34)/BL21 * 32 4 1 8 pET-28a-[bla.sub.MIR] (Y490)/BL21 * 32 4 1 8 pET-28a-[bla.sub.MIR] (CG76)/BL21 * 32 4 1 8 pET-28a-[bla.sub.MIR] (Y546)/BL21 * 32 4 1 8 pET-28a-[bla.sub.MIR] (CG85)/BL21 * 32 4 1 8 pET-28a-[bla.sub.MIR] (Y482)/BL21 * 32 4 1 8 pET-28a-[bla.sub.CTX-M-9a]/BL21 * 4 0.03 0.03 2 pET-28a-[bla.sub.SHV-12]/BL21 * 4 0.03 0.03 32 Strains PTZ IMC ATCC25922 4 0.25 BL21 0.5 0.5 pET-28a/BL21 * 0.5 0.5 CG34 2 4 Y490 4 2 CG76 32 8 Y546 8 4 CG85 4 4 Y482 2 0.5 pET-28a-[bla.sub.MIR] (CG34)/BL21 * 8 1 pET-28a-[bla.sub.MIR] (Y490)/BL21 * 4 1 pET-28a-[bla.sub.MIR] (CG76)/BL21 * 4 1 pET-28a-[bla.sub.MIR] (Y546)/BL21 * 8 1 pET-28a-[bla.sub.MIR] (CG85)/BL21 * 8 1 pET-28a-[bla.sub.MIR] (Y482)/BL21 * 8 1 pET-28a-[bla.sub.CTX-M-9a]/BL21 * 1 1 pET-28a-[bla.sub.SHV-12]/BL21 * 2 2 * IPTG was added to a final concentration of 1 mmol/L to standard Mueller-Hinton (M-H) agar plates. All MIC determinations were performed by agar dilution assays. FOX: cefoxitin; CAZ: ceftazidime; CTX: cefotaxime; AMP: ampicillin; KAN: kanamycin; CPZ: cefoperazone; CMN: cefminox; CEF: cefoselis; CFZ: cefazolin; CRO: ceftriaxone; ATM: aztreonam; PRL: piperacillin; PTZ: piperacillin/tazobactam; IMC: imipenem/cilastatin sodium hydrate. Table 4: General features of the Y546 genome. Chromosome Plasmid Size (bp) 4,787,919 208,741 GC content (%) 56.02 52.63 ORFs 4312 234 Known proteins 3476 (80.6%) 137 (58.5%) Hypothetical proteins 836 (19.4%) 97 (41.5%) Protein coding (%) 88.02 82.80 Average ORF length (bp) 977 738 Average protein length (aa) 324 245 tRNAs 85 0 rRNA operons (16S-23S-5S) *7 0 16S-23S-5S-5S Figure 2: Comparison of the MIR amino acid sequences from 6 E. cloacae isolates. * 10 ** 20 30 CG34 1 MMTKSLSCAL LLSVASSAFA APMSEKQLAE CG75 1 MMTKSLSCAL LLSVTSSAFA APMSEKQLAE CG85 1 MMTKSLSCAL LLSVASSAFA APMSEKQLAE Y482 1 MMTKSLSCAL LLSVASSAFA APMSEKQLAE Y490 1 MMTKSLSCTL LLSVASSAFA APMSEKQLAE Y546 1 MMTKSLSCAL LLSVTSSAFA APMSEKQLAE 40 50 60 CG34 1 VVERTVTPLM NAQAIPGMAV AVIYQGQPHY CG75 1 VVERTVTPLM NAQAIPGMAV AVIYQGQPHY CG85 1 VVERTVTPLM NAQAIPGMAV AVIYQGQPHY Y482 1 VVERTVTPLM NAQAIPGMAV AVIYQGQPHY Y490 1 VVERTVTPLM NAQAIPGMAV AVIYQGQPHY Y546 1 VVERTVTPLM NAQAIPGMAV AVIYQGQPHY 70 80 CG34 1 FTFGKADVAA NKPVTPQTLF CG75 1 FTFGKADVAA NKPVTPQTLF CG85 1 FTFGKADVAA NKPVTPQTLF Y482 1 FTFGKADVAA NKPVTPQTLF Y490 1 FTFGKADVAA NKPVTPQTLF Y546 1 FTFGKADVAA NKPVTPQTLF 90 100 * 110 CG34 81 ELGSISKTFT GVLGGDAIAR GEIALGDPVA CG75 81 ELGSISKTFT GVLGGDAIAR GETALGDPVA CG85 81 ELGSISKTFT GVLGGDAIAR GEIALGDPVA Y482 81 ELGSISKTFT GVLGGDAIAR GEIALGDPVA Y490 81 ELGSISKTFT GVLGGDAIAR GEIALGDPVA Y546 81 ELGSISKTFT GVLGGDAIAR GEIALGDPVA 120 130 140 CG34 81 KYWPELTGKQ WQGIRMLDLA TYTAGGLPLQ CG75 81 KYWPELTGKQ WQGIRMLDLA TYTAGGLPLQ CG85 81 KYWPELTGKQ WQGIRMLDLA TYTAGGLPLQ Y482 81 KYWPELTGKQ WQGIRMLDLA TYTAGGLPLQ Y490 81 KYWPELTGKQ WQGIRMLDLA TYTAGGLPLQ Y546 81 KYWPELTGKQ WQGIRMLDLA TYTAGGLPLQ 150 160 CG34 81 VPDEVTDTAS LLRFYQNWQP CG75 81 VPDEVTDTAS LLRFYQNWQP CG85 81 VPDEVTDTAS LLRFYQNWQP Y482 81 VPDEVTDTAS LLRFYQNWQP Y490 81 VPDEVTDTAS LLRFYQNWQP Y546 81 VPDEVTDTAS LLRFYQNWQP 170 180 190 CG34 161 QWKPGTTRLY ANASIGLFGA LAVKPSGMSY CG75 161 QWKPGTTRLY ANASIGLFGA LAVKPSGMSY CG85 161 QWKPGTTRLY ANASIGLFGA LAVKPSGMSY Y482 161 QWKPGTTRLY ANASIGLFGA LAVKPSGMSY Y490 161 QWKPGTTRLY ANASIGLFGA LAVKPSGMSY Y546 161 QWKPGTTRLY ANASIGLFGA LAVKPSGMSY 200 * 210 220 CG34 161 EQAMTTRVFK PLKLDNTWIN VPKAEEAHYA CG75 161 EQAMTTRVFK PLKLDHTWIN VPKAEEAHYA CG85 161 EQAMTTRVFK PLKLDHTWIN VPKAEEAHYA Y482 161 EQAMTTRVFK PLKLDHTWIN VPKAEEAHYA Y490 161 EQAMTTRVFK PLKLDHTWIN VPKAEEAHYA Y546 161 EQAMTTRVFK PLKLDHTWIN VPKAEEAHYA 230 240 CG34 161 WGYREGKAVH VSPGMLDAEA CG75 161 WGYREGKAVH VSPGMLDAEA CG85 161 WGYREGKAVH VSPGMLDAEA Y482 161 WGYREGKAVH VSPGMLDAEA Y490 161 WGYREGKAVH VSPGMLDAEA Y546 161 WGYREGKAVH VSPGMLDAEA 250 * 260 * 270 CG34 241 YGVKTNVKDM ASWLIANMKP DSLQAPSLKQ CG75 241 YGVKTNVKDM ASWVIANMKP DSLQAPSLKQ CG85 241 YGVKTNVKDM ASWVIANMKP DSLQAPSLKQ Y482 241 YGVKTNVKDM ASWVIANMKP DSLQAPSLKQ Y490 241 YGVKTNVKDM ASWLIANMKP DSLHAPSLKQ Y546 241 YGVKTNVKDM ASWLIANMKP DSLQAPSLKQ 280 290 * CG34 241 GIALAQSRYW RVGAMYQGLG WEMLNWPVDA CG75 241 GIALAQSRYW RVGAMYQGLG WEMLNWPVDA CG85 241 GIALAQSRYW RVGAMYQGLG WEMLNWPVDV Y482 241 GIALAQSRYW RVGAMYQGLG WEMLNWPVDV Y490 241 GIALAQSRYW RVGAMYQGLG WEMLNWPVDA Y546 241 GIALAQSRYW RVGAMYQGLG WEMLNWPVDA 310 320 CG34 241 KTVVGGSDNK VALAPLPVAE CG75 241 KTVVGGSDNK VALAPLPVAE CG85 241 KTVVGGSDNK VALAPLPVAE Y482 241 KTVVGGSDNK VALAPLPVAE Y490 241 KTVVGGSDNK VALAPLPVAE Y546 241 KTVVGGSDNK VALAPLPVAE 330 340 350 CG34 321 VNPPAPPVKA SWVHKTGSTG GFGSYVAFIP CG75 321 VNPPAPPVKA SWVHKTGSTG GFGSYVAFIP CG85 321 VNPPAPPVKA SWVHKTGSTG GFGSYVAFIP Y482 321 VNPPAPPVKA SWVHKTGSTG GFGSYVAFIP Y490 321 VNPPAPPVKA SWVHKTGSTG GFGSYVAFIP Y546 321 VNPPAPPVKA SWVHKTGSTG GFGSYVAFIP 360 370 380 CG34 321 EKQLGIVMLA NKSYPNPARV EAAYRILDAL CG75 321 EKQLGIVMLA NKSYPNPARV EAAYRILDAL CG85 321 EKQLGIVMLA NKSYPNPARV EAAYRILDAL Y482 321 EKQLGIVMLA NKSYPNPARV EAAYRILDAL Y490 321 EKQLGIVMLA NKSYPNPARV EAAYRILDAL Y546 321 EKQLGIVMLA NKSYPNPARV EAAYRILDAL CG34 321 Q CG75 321 Q CG85 321 Q Y482 321 Q Y490 321 Q Y546 321 Q
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Wu, Chongyang; Lin, Chaoqin; Zhu, Xinyi; Liu, Hongmao; Zhou, Wangxiao; Lu, Junwan; Zhu, Licheng; Bao|
|Publication:||International Journal of Genomics|
|Date:||Jan 1, 2018|
|Previous Article:||Genomic and Postgenomic Approaches to Understand Environmental Microorganisms.|
|Next Article:||Common DNA Variants Accurately Rank an Individual of Extreme Height.|