Complete-genome sequencing and comparative genomic characterization of blaNDM-5 carrying Citrobacter freundii isolates from a patient with multiple infections

Background The emergence and wide spread of carbapenemase-producing Enterobacteriaceae (CPE) poses a growing threat to global public health. However, clinically derived carbapenemase-producing Citrobacter causing multiple infections has rarely been investigated. Here we first report the isolation and comparative genomics of two blaNDM-5 carrying Citrobacter freundii (C. freundii) isolates from a patient with bloodstream and urinary tract infections. Results Antimicrobial susceptibility testing showed that both blaNDM-5 carrying C. freundii isolates were multidrug-resistant. Positive modified carbapenem inactivation method (mCIM) and EDTA-carbapenem inactivation method (eCIM) results suggested metallo-carbapenemase production. PCR and sequencing confirmed that both metallo-carbapenemase producers were blaNDM-5 positive. Genotyping and comparative genomics analyses revealed that both isolates exhibited a high level of genetic similarity. Plasmid analysis confirmed that the blaNDM-5 resistance gene is located on IncX3 plasmid with a length of 46,161 bp, and could successfully be transferred to the recipient Escherichia coli EC600 strain. A conserved structure sequence (ISAba125-IS5-blaNDM-5-trpF-IS26-umuD-ISKox3) was found in the upstream and downstream of the blaNDM-5 gene. Conclusions The data presented in this study showed that the conjugative blaNDM-5 plasmid possesses a certain ability to horizontal transfer. The dissemination of NDM-5-producing C. freundii isolates should be of close concern in future clinical surveillance. To our knowledge, this is the first study to characterize C. freundii strains carrying the blaNDM-5 gene from one single patient with multiple infections. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-023-09579-9.

Citrobacter freundii, belonging to the genus Citrobacter of the family Enterobacteriaceae, is rarely the causative pathogen of infections but can cause a wide spectrum of opportunistic nosocomial infections including respiratory tract, urinary tract, and bloodstream [5].Additionally, it can lead to neonatal meningitis and brain abscesses, which are associated with high mortality rates [6].The emergence of multidrug resistant C. freundii strains, particularly those that produce carbapenemase enzymes, has posed challenges for infection treatment and has become an increasing global public health concern.This is especially true for immunocompromised patients, who heavily rely on antibiotics [7].
In this study, we presented the results of completegenome sequencing and comparative genomic characterization of two NDM-5 producing C. freundii strains isolated from a patient with concurrent bloodstream and urinary tract infections.

Antimicrobial resistance profiles of both C. freundii isolates
Antimicrobial susceptibility testing of both C. freundii isolates revealed high MIC values for different drugs and different susceptibility-resistance levels depending on the drug tested.However, both isolates exhibited similar phenotypic antibiotic susceptibility, demonstrating a multi-drug resistant (MDR) characteristic.Both C. freundii isolates exhibited resistance to amoxicillin/ clavulanic acid, piperacillin/tazobactam, ceftazidime, ceftriaxone, cefepime, cefotaxime, ciprofloxacin, levofloxacin, trimethoprim/sulfamethoxazole, and gentamicin as shown in Table 1.In contrast, isolates remained susceptible to amikacin, aztreonam, fosfomycin and tigecycline (Table 1).Intermediate resistance was exhibited when isolates were cultured in the presence of imipenem, meropenem and polymixin B (Table 1).

Phenotype and genotype detection revealed the mechanism of carbapenems non-susceptibility
Positive modified carbapenem inactivation method (mCIM) and EDTA-carbapenem inactivation method (eCIM) results suggested metallo-carbapenemase production of both C. freundii isolates.Then, PCR and sequencing were used to study the molecular determinant, and the results demonstrated that both metallocarbapenemase producers were bla NDM-5 positive, thus confirming their intermediate resistance phenotype towards imipenem and meropenem.Furthermore, in silico analysis showed that both bla NDM-5 carrying isolates carried additional resistance genes conferring resistance to β-lactams (bla TEM-1B , bla OXA-1 , bla CMY-48 , bla DHA-1 ), this genetic profile correlated with their phenotypic resistance to ceftazidime, ceftriaxone, cefepime, and cefotaxime.Moreover, they carried resistance genes for aminoglycosides (aac(6')-Ib-cr, aac(3)-IId, aadA1), which correlated with phenotypic resistance to gentamicin.The presence of the qnrB4 gene indicated resistance to quinolones, which was in line with their resistance to ciprofloxacin and levofloxacin.Furthermore, the isolates carried resistance genes for trimethoprim/sulfamethoxazole (dfrA1, sul1, sul2), which were associated with their resistance to trimethoprim/sulfamethoxazole.The presence of amphenicol resistance genes (catA2, catB3), a tetracycline resistance gene tet(D), and a macrolides resistance gene mph(A) correlated with their resistance to ciprofloxacin, tetracycline, and azithromycin, respectively.Although the presence of the rifamycin resistance gene ARR-3 contributed to both isolates having MIC values exceeding 128 mg/L for rifamycin, the absence of a defined breakpoint prevented the determination of the resistance phenotype with certainty.

High degree of genetic similarity of both NDM-5-producing C. freundii isolates
An average nucleotide identity blast (ANIb) analysis, which measures the nucleotide-level genomic similarity between the coding regions of two genomes, was performed.The analysis revealed that the two isolates exhibited a remarkable similarity of over 99% (Supplementary Fig. 1).Additionally, when core single nucleotide polymorphisms (SNPs) were calculated, only five base differences were detected between the two isolates.MLST typing revealed that both NDM-5-producing C. freundii isolates belonged to the sequence typing ST 22. Singlenucleotide polymorphism (SNP)-based phylogenetic tree for both isolates and other 78 NDM-producing strains indicated that DY2007 and DY2010 showed high degree of similarity, and were closely related to various strains from different countries, including Myanmar, USA, France, Australia, Malaysia, Germany, China, and Singapore (Fig. 1, Table S1).

Characterization of both bla NDM-5 carrying plasmids
S1-PFGE and Southern blotting revealed that both C. freundii isolates contained a ~ 50 kb plasmid harbouring the bla NDM-5 gene (Fig. 2).Plasmid replicons analysis indicated that both pNDM-5 were IncX3 type with a length of 46,161 bp.The results of the conjugation assay showed the bla NDM-5 gene of both isolates could successfully be transferred to the recipient Escherichia coli EC600 strain.The results were further confirmed by PCR using bla NDM-5 specific primers and sequencing.

Whole-genome sequencing analysis of both bla NDM-5 carrying C. freundii isolates
The genomic characteristics of both bla NDM-5 carrying C. freundii isolates were shown in Table S2.DY2007 contained a circular chromosome and two plasmids with a genome size of 5,253,532 bp, while DY2010 consisted of a circular chromosome and three plasmids with a genome size of 5,260,876 bp.The average GC content of both genomes was 51.6%.The complete sequence of bla NDM-5 carrying pNDM-5 was covered using combinatorial PCR and standard Sanger sequencing to accomplish sequence integrality of contigs.It was found that both isolates harbored an identical pNDM-5 plasmid of 46,161 bp length, with a GC content of 46.65% and 65 predicted coding sequences (Fig. 3A).The plasmid carried multiple coding genes, including IncX plasmid conjugal transfer associated genes, antibiotic resistance genes, stability related genes, functional protein coding genes, mobile element associated genes, and other hypothetical genes.A search of the nr/nt database revealed a 100% identity to Escherichia coli strain WCHEC020031 plasmid pNDM5_020031 (GenBank accession number: CP033399.1)at 100% coverage, 100% identity to Klebsiella pneumoniae strain 19110124 plasmid p19110124-3 (GenBank accession number: CP064177.1)at 100% coverage, and 99.98% identity to Escherichia coli strain L53 plasmid pL53-4 (GenBank accession number: CP034737.1)at 99% coverage.Furthermore, a conserved structure sequence (ISAba125-IS5-bla NDM-5 -trpF-IS26-umuD-ISKox3) was found in the upstream and downstream of bla NDM-5 (Fig. 3B).[5][6][7]10].In this study, we reported the complete-genome sequencing and comparative genomic characterization of two bla NDM-5 carrying C. freundii isolates from an inpatient with concurrent bloodstream and urinary tract infection.To the best of our knowledge, this is the first study to characterize C. freundii strains carrying the bla NDM-5 gene from one single patient with multiple infections.Antibiotic susceptibility testing indicated the multidrug-resistant characteristic of both isolates, with resistant to classical β-lactams, aminoglycosides, macrolides, quinolones, and sulfonamides, and intermediate to polypeptide antibiotics and carbapenems, the resistant phenotype of which was consistent with resistant genotype (Table 1).
Phylogenetic analysis determined the high homology of both ST 22 isolates, which further highlighted the potential threat of systemic infections caused by this multidrug-resistant (MDR) strain.The multidrug-resistant bla NDM-5 carrying C. freundii isolates have been reported previously from South Korea and Nigeria [8,9].
Compared to vertical transmission, horizontal transfer seems to be more threatening for antimicrobial resistance.In this study, two isolates both successfully transferred bla NDM-5 and carbapenem non-susceptible phenotype to the recipient strain E. coli EC600, which confirmed the horizontal gene transfer characteristic of NDM-5.In recent years, IncX3 plasmids harboring bla NDM variants have increasingly been characterized worldwide.It has also been proved that IncX3 plasmids plays an important role in the dissemination of bla NDM-5 Fig. 1 The core-genome phylogenetic tree generated by kSNP.Blue color of the right panel indicates positive antibiotic resistance genes of the corresponding strains gene in Enterobacteriaceae [12].Tian et al. observed that horizontal gene transfer (HGT) of bla NDM-5 among distinct Enterobacteriaceae species was mostly mediated by IncX3 plasmids [16].A previous study found that bla NDM-5 might spread among humans and the environment via IncX3 plasmids in an intensive vegetable farming area in eastern China [19].Consistently, bla NDM-5 was located on an IncX3 plasmid (~ 50 kb) in this study.However, bla NDM-5 was additionally detected on IncF and IncI1 plasmids [20,21], highlighting the significant compatibility and transmission hazard associated with bla NDM-5 .
In order to further understand the evolutionary relationship of bla NDM-5 harboring C. freundii isolates, we downloaded the nucleotide sequences of 78 NDM-producing strains from NCBI and analyzed their homology with both isolates in this study.The results revealed that 17 strains from around the world were closely related to both isolates, suggesting parallel evolution of these isolates.In plasmids, genes are typically associated with mobile genetic elements such as transposons (Tn) and insertion sequences (IS).According to the wholegenome sequencing analysis of the pNDM-5 plasmids, the comparison of genetic context flanking bla NDM-5 in both isolates was mostly identical to published plasmids, i.e., ISAba125-IS5-bla NDM-5 -trpF-IS26-umuD-ISKox3.The bla NDM-5 genetic structure is widespread in Enterobacteriaceae for bla NDM horizontal transfer and has been reported in bla NDM-5 and bla NDM-9 transmission [22].Characterization of the bla NDM-5 genetic contents revealed that it was flanked by multi-insertional sequences.Of these, ISAba125 was conservative in bla NDM-5 -positive isolates.It is consistent with the discovery that ISAba125 (intact or truncated) upstream of bla NDM is common in bla NDM genetic settings [23], indicating its' important role in bla NDM transmission.NDM, a highly prevalent plasmid-borne metallo-β-lactamase, has been identified in various species of Enterobacteriaceae worldwide.Its frequent co-occurrence with ISAba125 suggests a potential origin from Acinetobacter spp., a bacterium where this association is commonly observed [24].Tn125, a composite transposon based on ISAba125, has been reported as one of the genetic elements implicated in the dissemination of bla NDM .However, in Enterobacteriaceae, Tn125 exhibits interruptions or truncations, leading to diverse genetic contexts for bla NDM [3].
Whole-genome sequencing (WGS) plays a pivotal role in clinical cases involving systemic infections.It offers a comprehensive perspective on the complete genome of the pathogen, allowing for meticulous analysis of genetic variations, resistance mechanisms, and the identification of potential virulence factors.WGS enhances our understanding of the pathogenesis of systemic infections, aids in tracing transmission routes, and assists in formulating suitable treatment strategies [25][26][27].This study bears   with Escherichia coli strain WCHEC020031 plasmid pNDM5_020031 (GenBank accession number: CP033399.1),Klebsiella pneumoniae strain 19,110,124 plasmid p19110124-3 (GenBank accession number: CP064177.1),and Escherichia coli strain L53 plasmid pL53-4 (GenBank accession number: CP034737.1).The figure was plotted using BRIG, and the bla NDM-5 gene was highlighted in red.B Genetic environment of bla NDM-5 on pNDM-5 and related plasmids.Open reading frames were indicated as arrows.Shared areas with highly similar sequences were drawn by lake green.Conjugal transfer associated genes were shown by brown; bla NDM-5 gene were indicated by red arrows; functional protein coding genes were colored by green; other antibiotic resistance genes were colored by blue

Conclusions
In this study, we sequenced and characterized the comparative genomes of two bla NDM-5 carrying C. freundii isolates from an inpatient with multiple infections, and found different antimicrobial resistant genes in a transferable IncX3-type plasmid.The dissemination of this MDR isolate should be of close concern in future clinical surveillance.

Case presentation and bacterial isolates
Two C. freundii isolates were collected from a 62 years old female bladder cancer inpatient with multiple infections of the Affiliated Dongyang Hospital of Wenzhou Medical University (Wenzhou, China) in 2020.One isolate (DY2007) was first recovered from bloodstream on January 21 and the other one (DY2010) was subsequently obtained from urinary tract on February 12.The patient presented with lower back pain, accompanied by fever and chills, and was admitted on January 20, 2020.She received intravenous administration of 2g of ceftriaxonesulbactam every 8 h for antimicrobial treatment until February 1.Additionally, on January 23, based on the strain's antimicrobial susceptibility testing results, amikacin injection was added to the treatment regimen at a dosage of 0.2g every 12 h.This treatment continued until February 1, when the patient's symptoms improved, and she requested discharge.Subsequently, on February 9, the patient was readmitted due to lower back pain and received the aforementioned amikacin treatment until February 17, when the patient made a complete recovery and was subsequently discharged.The strains were isolated using sheep blood agar cultured overnight at 37 • C and were initially identified using MALDI-TOF MS (Bio-Mérieux, France).The isolates were stored in 30% glycerol at -80 • C until further analysis.

Antimicrobial susceptibility testing
The minimum inhibitory concentrations (MICs) of 17 antibiotics, including amoxicillin/clavulanic acid, piperacillin/tazobactam, ceftazidime, ceftriaxone, cefepime, cefotaxime, ciprofloxacin, levofloxacin, imipenem, meropenem, trimethoprim/sulfamethoxazole, amikacin, gentamicin, aztreonam, fosfomycin, tigecycline, and polymixin B, were determined by the VITEK 2 system with AST-GN13 card and the agar dilution method.The breakpoint of tigecycline was interpreted according to the recommendations of the Food and Drug Administration (FDA) [28], and the breakpoint of other antibiotics were interpreted according to the Clinical and Laboratory Standards Institute (CLSI) 2021 guidelines [29].

Carbapenemases detection and molecular mechanisms identification
The mCIM and eCIM method were used to determine carbapenemase production according to the CLSI guidelines.Carbapenemase genes (bla KPC , bla NDM , bla IMP , bla VIM , and bla OXA-48 ) were identified by PCR amplification as in our previous publication [30].Positive amplification products were then sequenced for verification and subtype typing.

Plasmid characterization and conjugation assay
Plasmid sizes of the strains were determined using the S1 nuclease pulsed-field gel electrophoresis (S1-PFGE) method.The location of the bla NDM-5 gene was investigated by Southern blotting with a specific digoxigenin-labelled bla NDM-5 probe using the DIGHigh Prime DNA Labeling and Detection Starter Kit II (Roche Diagnostics, Germany) [10].Replicon types of plasmid incompatibility (Inc) groups were identified by multiplex PCR as previously described [10].The plasmid conjugation and transformation methods were performed to verify the transferability of the NDMbearing plasmid with E. coli 600 as a recipient strain.The transconjugants were then screened on BHI agar plates supplemented with 2 mg/L meropenem and were identified by MALDI-TOF MS.NDM-bearing recipient strain was confirmed by PCR and sequencing.

Whole genome sequencing
Genomic DNA of both NDM-5 producing C. freundii isolates was extracted using the QIAmp DNA Mini Kit (Qiagen, Germany).A Qubit Fluorometer (Thermo scientific, USA) was then used to determine the concentration and purity of DNA.Sequencing libraries were prepared using the Illumina Nextera XT Kit and sequenced using Illumina HiSeq 4000-PE150 platform (Illumina, USA).Raw sequencing data of both isolates were assembled using SOAP de novo software [34] and were deposited in GenBank under the following accession numbers: JAJDSQ000000000, and CP086287-CP086290, respectively.Antimicrobial resistance genes and plasmid replicon types were matched to the Center for Genomic Epidemiology (http:// www.genom icepi demio logy.org/) Resfinder and Plasmid finder databases.The gaps were covered using combinatorial PCR to accomplish sequence integrality of contigs.The RAST server (http:// rast.nmpdr.org/) was used to annotate the bacterial genomes, and the ISFinder database (https:// www-is.bioto ul.fr/) was used to identify IS elements and transposons.Multiple plasmid alignment was conducted and plotted between the bla NDM-5 -harboring plasmid (named pNDM-5) and the reference plasmid using the BLAST Ring Image Generator (BRIG) [35].Easyfig 2.2.3 was used to analyze the genetic environment surrounding the bla NDM-5 resistance gene [36].

Fig. 2
Fig.2bla NDM-5 gene location analysis.A S1-PFGE of both bla NDM-5 carrying C. freundii isolates DY2007 and DY2010.Salmonella enterica serotype H9812 was used as molecular marker.B Corresponding Southern blotting analysis using bla NDM-5 -specific probe.A and B was cropped from different gels.Full-length blots/gels are presented in Supplementary Fig.2

2
Fig.2bla NDM-5 gene location analysis.A S1-PFGE of both bla NDM-5 carrying C. freundii isolates DY2007 and DY2010.Salmonella enterica serotype H9812 was used as molecular marker.B Corresponding Southern blotting analysis using bla NDM-5 -specific probe.A and B was cropped from different gels.Full-length blots/gels are presented in Supplementary Fig.2

Fig. 3
Fig.3Genomic analyses of pNDM-5 plasmid.A Comparison of the pNDM-5 plasmid sequence identified in isolates DY2007 and DY2010 with Escherichia coli strain WCHEC020031 plasmid pNDM5_020031 (GenBank accession number: CP033399.1),Klebsiella pneumoniae strain 19,110,124 plasmid p19110124-3 (GenBank accession number: CP064177.1),and Escherichia coli strain L53 plasmid pL53-4 (GenBank accession number: CP034737.1).The figure was plotted using BRIG, and the bla NDM-5 gene was highlighted in red.B Genetic environment of bla NDM-5 on pNDM-5 and related plasmids.Open reading frames were indicated as arrows.Shared areas with highly similar sequences were drawn by lake green.Conjugal transfer associated genes were shown by brown; bla NDM-5 gene were indicated by red arrows; functional protein coding genes were colored by green; other antibiotic resistance genes were colored by blue