Whole genome sequencing uncovers a novel IND-16 metallo-β-lactamase from an extensively drug-resistant Chryseobacterium indologenes strain J31
- Tingting Wang1, 2,
- Xiawei Jiang3,
- Chunyan Feng4,
- Ang Li1,
- Huihui Dong1,
- Shaoqiang Wu†4Email author and
- Beiwen Zheng†1Email author
© The Author(s) 2016
Received: 10 August 2016
Accepted: 12 October 2016
Published: 21 October 2016
Chryseobacterium indologenes is an emerging opportunistic pathogen in hospital-acquired infection, which is intrinsically resistant to most antimicrobial agents against gram-negative bacteria. In the purpose of extending our understanding of the resistance mechanism of C. indologenes, we sequenced and analyzed the genome of an extensively antibiotic resistant C. indologenes strain, isolated from a Chinese prostate cancer patient. We also investigated the presence of antibiotic resistance genes, particularly metallo-β-lactamase (MBL) genes, and performed a comparative genomic analysis with other Chryseobacterium species.
16s rRNA sequencing indicated the isolate belongs to C. indologenes. We assembled a total of 1095M bp clean-filtered reads into 171 contigs by de novo assembly. The draft genome of C. indologenes J31 consisted of 5,830,795 bp with a GC content of 36.9 %. RAST analysis revealed the genome contained 5196 coding sequences (CDSs), 28 rRNAs, 81 tRNAs and 114 pseudogenes. We detected 90 antibiotic resistance genes from different drug classes in the whole genome. Notably, a novel bla IND allele bla IND-16 was identified, which shared 99 % identity with bla IND-8 and bla IND-10. By comparing strain J31 genome to the closely four related neighbors in the genus Chryseobacterium, we identified 2634 conserved genes, and 1449 unique genes.
In this study, we described the whole genome sequence of C. indologenes strain J31. Numerous resistance determinants were detected in the genome and might be responsible for the extensively antibiotic resistance of this strain. Comparative genomic analysis revealed the presence of considerable strain-specific genes which would contribute to the distinctive characteristics of strain J31. Our study provides the insight of the multidrug resistance mechanism in genus Chryseobacterium.
Chryseobacterium indologenes, of which the taxonomic characteristics were first described in 1983 , is a common gram-negative bacterium and belongs to the genus Chryseobacterium. It is widely found existing in natural and nosocomial environment and is occasionally isolated from the human gut . As an uncommon pathogen, it usually causes opportunistic infection to the susceptible populations such as the infant, the elder, the immuno-suppressed patient and the long-term inpatient . Since the first case of location outbreak emerged in Taiwan, China , C. indologenes has been documented to cause invasive infections, which usually result in bacteremia or pneumonia with high mortality [5–7].
Chryseobacterium indologenes has been proved to be intrinsically resistant to most antimicrobial agents often used to treat gram-negative bacteria [8, 9], but the mechanism of the multidrug resistance is not clear. The metallo-β-lactamase (MBL), which can cleave the β-lactam ring of antibiotics of the penicillin family, is thought to be closely associated with multidrug resistance of this bacterium . Since the first type of class B MBL gene, bla IND, was identified in C. indologenes, there are fifteen IND variants have been deposited in Genbank to date. These enzymes share 27–92 % identity with that of IND-1 at the amino acid sequence level [11–14]. The investigation of new bla IND allele would extend our understanding of the resistance mechanism of this bacterium.
The genomic information would provide more details about antibiotic resistance genes and help decipher the antibiotic resistance. However, only three whole genome sequences of different C. indologenes strains have been deposited in the NCBI genome database (https://www.ncbi.nlm.nih.gov/genome). In this study, we reported the whole genome sequence of one C. indologenes strain isolated from a prostate cancer patient in China and analyzed the multidrug resistance at the genomic level.
Strain information, antimicrobial susceptibility testing, DNA isolation and 16S rRNA sequencing
The urine sample was from a patient with prostate cancer in June 2014. The isolate was cultured aerobically on Columbia blood agar base plate at 37 °C. Susceptibility testing for ampicillin, amikacin, ciprofloxacin, levofloxacin, cefoperazone, nitrofurantoin, imipenem, iobramycin, piperacillin-tazobactam, colistin, cephalosporins and trimethoprim-sulfamethoxazole (TMP-SMZ) were determined by the disc diffusion method and interpreted concerning the CLSI guidelines. Late log phase cells were harvested, and genomic DNA was extracted using a DNeasy Blood & Tissue Kit (Qiagen, Germany) according to the manufacturer’s instruction. In the purpose of taxonomic identification, we amplified the 16S rRNA gene with a 16S rRNA universal primer set and the PCR product was sequenced as previously . A neighbor-joining phylogenetic tree was constructed based on the Tamura-Nei model using MEGA6 (http://www.megasoftware.net/).
Genome sequencing, assembly, and annotation
We conducted genome sequencing, assembly, and annotation by following previous studies [16, 17]. In brief, whole genome sequencing was performed using an Illumina HiSeq 2000 genomic sequencer with a 2 × 100 paired-end sequencing strategy (DNA libraries was in a size of 500-bp insertion). We removed all reads with adaptor contamination and with unknown nucleotides comprising more than 5 %. Then, low-quality reads with ambiguous sequence “N” were discarded. Subsequently, clean-filtered reads were de novo assembled into scaffolds using Velvet 1.2.07 . We used VelvetOptimiser for automatically optimizing optimum k-Mer value. And then PAGIT (Post-Assembly Genome Improvement Toolkit)  was used to extend the initial contigs and correct sequencing errors. Open reading frames (ORFs), tRNAs and rRNAs were identified using Glimmer version 3.0 , tRNAscan-SE  and RNAmmer , respectively. The functional genes were annotated and classified using the RAST (rapid annotation using subsystem technology) server  and the COGs (clusters of orthologous groups of proteins) databases . We predicted the plasmid replicons by the PlasmidFinder Tool (https://cge.cbs.dtu.dk/services/PlasmidFinder/). ISfinder (https://www-is.biotoul.fr/blast.php) was employed to search the IS sequences in the genome, with an e value of 1E−3.
Antibiotic resistance genes prediction and virulence factors analysis of J31
The protein-coding sequences were annotated by antibiotic resistance database (ARDB)  and antibiotic resistance gene-ANNOTation (ARG-ANNOT) with default parameters . We further verified all these putative ARGs through a BLAST search with cut-off e value of 1E−0.5. Using PCR sequencing, we verified IND-16 gene and downloaded other IND protein sequences from NCBI website. Multiple sequence alignments of the amino acid sequences were performed using Clustal Omega (http://www.ebi.ac.uk/Tools/msa/clustalo), and then the alignment results were visualized using BoxShade (http://www.ch.embnet.org/software/BOX_form.html). The virulence factors were predicted by VFDB database (http://www.mgc.ac.cn/VFs/main.htm) using BLAST with an e value threshold of 1E−5.
Comparative genomic analysis was performed by orthology identification method as we described previously [15, 17]. We downloaded the genome sequences used in the comparative analysis, from NCBI genome database. BLASTN were used for aligning the whole genomes between strain J31 and other five Chryseobacterium species and then genome alignment visualization was performed using BLAST ring image generator (BRIG) .
We isolated a single colony of strain J31 and purified the genomic DNA from a pure culture of the isolate. The strain persevered in the State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, Zhejiang University. The genome confirmed to C. indologenes by 16S rRNA sequencing. We constructed a phylogenetic tree from the 16s rRNA sequences and the position classified to C. indologenes. We assessed the potential contamination of the genomic library by other microorganisms by using a BLAST search against the non-redundant database.
Identification of strain J31
General genome properties
We performed whole genome sequencing by Illumina HiSeq 2000 system with 2 × 100 bp paired-end in length. After quality control, we assembled the filtered 1095M bp reads into contigs. The draft genome sequence of C. indologenes strain J31 revealed a total size of 5,830,795 bp and a GC content of 36.9 %. The assembled genome covered an average depth of 187.8-fold and contained 171 contigs, of which the largest one consisted of 476,013 bp and the length of N50 was 106,549 bp. Annotation of the contigs identified 5196 coding sequences (CDSs), 28 rRNAs, 81 tRNAs and 114 pseudogenes. BLASTN search between strain J31 and other five Chryseobacterium species revealed high similarity between Chryseobacterium isolates (Fig. 1b).
Antimicrobial susceptibility profiles and antibiotic resistance genes
The in vitro antimicrobial susceptibility testing demonstrates that the strain J31 is only susceptible to TMP-SMZ, but resistant to all the other tested antibiotics including ampicillin, amikacin, ciprofloxacin, levofloxacin, cefoperazone, nitrofurantoin, imipenem, iobramycin, piperacillin-tazobactam, colistin and cephalosporins (Additional file 1: Table S1). We screened the antibiotic resistance genes (ARGs) in the genome-wide scale in order to further explore the genetic basis of extensive resistance in this strain. In silico analysis revealed the presence of a considerable number of putative ARGs from different drug classes, 90 genes in total (Additional file 2: Table S2), which was more than that observed in Chryseobacterium oranimense . This isolate encoded 30 β-lactamase genes, which included amber class A, MBL, class C, and other novel β-lactamases genes. These ARGs might confer high-level resistance to cefepime, cefoperazone and piperacillin-tazobactam, which have been active in previous studies [11, 28]. In addition, we detected genes corresponding to rifampin, aminoglycosides, phenicols, sulfonamide, macrolide and trimethoprim resistance, which are consistent with the phenotypic results (Additional file 2: Table S2). Moreover, Multidrug efflux pumps including ABC-type transporter, MFS superfamily transporter and RND family transporter are also scattered among the genome. Therefore, these predicted ARGs might contribute to the multidrug resistance of stain J31 to the tested antibiotics and other non-tested antibiotics.
To further explore the resistance mechanism of J31, we used PlasmidFinder to detect the potential plasmids among the whole genome sequence. However, PlasmidFinder did not find any plasmid, and the result was also verified by conjugation experiments. We further investigated the other three assemblies of C. indologenes in NCBI genome database (https://www.ncbi.nlm.nih.gov/genome/) and no typing plasmid replicons could be found in all the contigs. For the genus Chryseobacterium, only two in thirty-eight draft or complete genomes consist of chromosomes and plasmids, which suggests Chryseobacterium, unlike the genus Enterococcus, might carry the resistant genes in the chromosome [16, 29].
We scanned each ARG flanking sequences in a range of 10-kb for IS sequences and junction associated proteins. Only one partial IS sequence, 107 bp ISR of ISBbi1 in IS1595 family, was found locating at 8958 bp upstream of one subclass B1 MBL (Additional file 3: Table S3). We found two integrases, a 304-aa-long site-specific tyrosine recombinase XerD at 2908 bp downstream of the penicillin-binding protein gene and a 422-aa-long hypothetical integrase at 4146 bp downstream of the putative ABC transporter gene, and two crossover junction endodeoxyribonucleases, a 184-aa-long RuvC protein at 5578 bp downstream of one β-lactamase gene and a 138-aa-long RuvA protein at 1748 bp downstream of another β-lactamase gene (Additional file 3: Table S3). Maybe there are some unknown proteins flanking the predicted ARGs which might function as the insertion elements or transposases for transferring antibiotic resistance.
A new IND-type MBL variant IND-16
We identified a novel IND-type MBL variant from in silico ARGs prediction, and designated it as IND-16. We confirmed this allele by PCR sequencing and searched the sequence against the GenBank database. The BLASTN search indicated the sequence was highly similar to bla IND genes of Chryseobacterium species. Multiple sequence alignments demonstrated that IND-16 protein shared 99 % identity with IND-8 and IND-10 and was conserved with other IND type MBLs (Additional file 4: Figure S1). According to the distinctive resistance of the strain, IND-16 was likely to contribute to the resistance of J31 to carbapenems.
Chryseobacterium indologenes is an emerging pathogen associated with indwelling devices and immunosuppression, the pathogenesis should be investigated. We performed a BLASTP search against VRDB database and found three virulence factors which contained Clp protease ClpC, molecular chaperone GroEL and ATP-dependent chaperone ClpB (Additional file 5: Table S4). In addition, we detected four conjugative transposon clusters, Tra gene clusters, in the draft genomes which inferred the strain J31 was with potential pathogenesis (Additional file 6: Table S5).
Comparative analysis with other C. indologenes strains
This study represents the first genomic features of a multidrug resistant C. indologenes isolate and demonstrates the comparative genomic analysis of strain J31 to other Chryseobacterium species. Notably, it reveals the presence of numerous resistance determinants that helped this strain resistant to many antimicrobial agents, which makes it potential candidates for nosocomial infections. Furthermore, our analysis advances the understanding of resistance in genus Chryseobacterium.
BZ and SW were involved in the overall experimental design. TW, XJ and CF performed microbiology and molecular biology experiments. XJ and AL generated and analyzed the sequencing data. BZ and TW participated in all discussions of data analysis and write the manuscript. All authors read and approved the final manuscript.
We thank Dr. Syed Shuja Zaidi and Dr. Kristen E. Murfin for their linguistic assistance during the preparation of this manuscript.
The authors declare that they have no competing interests.
Availability of supporting data
The 16S rRNA sequence and nucleotide sequence of the bla IND-16 gene of C. indologenes J31 have been deposited in GenBank with accession numbers of KR297230 and KT235893, respectively. Whole Genome Shotgun project of C. indologenes J31 has been deposited at DDBJ/EMBL/GenBank under the accession LAZY00000000. The version described in this paper is version LAZY01000000.
Ethics approval and consent to participate
This research was approved by the Research Ethics Committee of the First Affiliated Hospital, School of Medicine, Zhejiang University.
This work was supported by the National Key R&D Program of China (No. 2016YFD0501105) and the National Natural Science Foundation of China (Nos. 81301461 and 41406140).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Yabuuchi E, Kaneko T, Yano I, Moss CW, Miyoshi N. Sphingobacterium gen. nov., Sphingobacterium spiritivorum comb. nov., Sphingobacterium multivorum comb. nov., Sphingobacterium mizutae sp. nov., and Flavobacterium indologenes sp. nov.: glucose-nonfermenting gram-negative rods in CDC groups IIK-2 and IIb. Int J Syst Bacteriol. 1983;33:580–98.View ArticleGoogle Scholar
- Mukerji R, Kakarala R, Smith SJ, Kusz HG. Chryseobacterium indologenes: an emerging infection in the USA. BMJ Case Rep 2016; 2016.Google Scholar
- Ozcan N, Dal T, Tekin A, Kelekci S, Can S, Ezin O, et al. Is Chryseobacterium indologenes a shunt-lover bacterium? A case report and review of the literature. Infez Med. 2013;21:312–6.PubMedGoogle Scholar
- Hsueh PR, Hsiue TR, Wu JJ, Teng LJ, Ho SW, Hsieh WC, et al. Flavobacterium indologenes bacteremia: clinical and microbiological characteristics. Clin Infect Dis. 1996;23:550–5.View ArticlePubMedGoogle Scholar
- Nemli SA, Demirdal T, Ural S. A case of healthcare associated pneumonia caused by Chryseobacterium indologenes in an immunocompetent patient. Case Rep Infect Dis. 2015;2015:483923.PubMedPubMed CentralGoogle Scholar
- Deng L, Li MF, Li YH, Yang JL, Zhou X. Chryseobacterium indologenes in four patients with leukemia. Transpl Infect Dis. 2015;17:583–7.View ArticlePubMedGoogle Scholar
- Olbrich P, Rivero-Garvía M, Falcón-Neyra MD, Lepe JA, Cisneros JM, Marquez-Rivas J, et al. Chryseobacterium indologenes central nervous system infection in infancy: an emergent pathogen? Infection. 2014;42:179–83.View ArticlePubMedGoogle Scholar
- Chang Y-C, Lo H-H, Hsieh H-Y, Chang S-M. Identification, epidemiological relatedness, and biofilm formation of clinical Chryseobacterium indologenes isolates from central Taiwan. J Microbiol Immunol Infect. 2015;48:559–64.View ArticlePubMedGoogle Scholar
- Chen F-L, Wang G-C, Teng S-O, Ou T-Y, Yu F-L, Lee W-S. Clinical and epidemiological features of Chryseobacterium indologenes infections: analysis of 215 cases. J Microbiol Immunol Infect. 2013;46:425–32.View ArticlePubMedGoogle Scholar
- Bebrone C. Metallo-beta-lactamases (classification, activity, genetic organization, structure, zinc coordination) and their superfamily. Biochem Pharmacol. 2007;74:1686–701.View ArticlePubMedGoogle Scholar
- Zeba B, De Luca F, Dubus A, Delmarcelle M, Simporé J, Nacoulma OG, et al. IND-6, a highly divergent IND-type metallo-beta-lactamase from Chryseobacterium indologenes strain 597 isolated in Burkina Faso. Antimicrob Agents Chemother. 2009;53:4320–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamaguchi Y, Takashio N, Wachino J, Yamagata Y, Arakawa Y, Matsuda K, et al. Structure of metallo-beta-lactamase IND-7 from a Chryseobacterium indologenes clinical isolate at 1.65-A resolution. J Biochem. 2010;147:905–15.View ArticlePubMedGoogle Scholar
- Perilli M, Caporale B, Celenza G, Pellegrini C, Docquier JD, Mezzatesta M, et al. Identification and characterization of a new metallo-beta-lactamase, IND-5, from a clinical isolate of Chryseobacterium indologenes. Antimicrob Agents Chemother. 2007;51:2988–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Bellais S, Léotard S, Poirel L, Naas T, Nordmann P. Molecular characterization of a carbapenem-hydrolyzing beta-lactamase from Chryseobacterium (Flavobacterium) indologenes. FEMS Microbiol Lett. 1999;171:127–32.PubMedGoogle Scholar
- Zheng B, Li A, Jiang X, Hu X, Yao J, Zhao L, et al. Genome sequencing and genomic characterization of a tigecycline-resistant Klebsiella pneumoniae strain isolated from the bile samples of a cholangiocarcinoma patient. Gut Pathog. 2014;6:1.View ArticleGoogle Scholar
- Sharma P, Gupta SK, Diene SM, Rolain J-M. Whole-genome sequence of Chryseobacterium oranimense, a colistin-resistant bacterium isolated from a cystic fibrosis patient in France. Antimicrob Agents Chemother. 2015;59:1696–706.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng B, Zhang F, Chai L, Yu G, Shu F, Wang Z, et al. Permanent draft genome sequence of Geobacillus thermocatenulatus strain GS-1. Mar Genomics. 2014;18:129–31.View ArticleGoogle Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Swain MT, Tsai IJ, Assefa SA, Newbold C, Berriman M, Otto TD. A post-assembly genome-improvement toolkit (PAGIT) to obtain annotated genomes from contigs. Nat Protoc. 2012;7:1260–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Delcher AL, Bratke KA, Powers EC, Salzberg SL. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics. 2007;23:673–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Lagesen K, Hallin P, Rødland EA, Stærfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Aziz RK, Bartels D, Best AA, De Jongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:1.View ArticleGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinform. 2003;4:1.View ArticleGoogle Scholar
- Liu B, Pop M. ARDB-antibiotic resistance genes database. Nucleic Acids Res. 2009;37:D443–7.View ArticlePubMedGoogle Scholar
- Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L, et al. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob Agents Chemother. 2014;58:212–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Alikhan N-F, Petty NK, Ben Zakour NL, Beatson SA. BLAST ring image generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011;12:402.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsueh PR, Teng LJ, Ho SW, Hsieh WC, Luh KT. Clinical and microbiological characteristics of Flavobacterium indologenes infections associated with indwelling devices. J Clin Microbiol. 1996;34:1908–13.PubMedPubMed CentralGoogle Scholar
- Brkić DV, Zlopaša O, Bedenić B, Plečko V. Chrysobacterium gleum infection in patient with extreme malnutrition and hepatic lesion-case report. Signal vitae. 2015;10:50–2.Google Scholar