Next Article in Journal
Itraconazole Loaded Biosurfactin Micelles with Enhanced Antifungal Activity: Fabrication, Evaluation and Molecular Simulation
Next Article in Special Issue
Scolopax rusticola Carrying Enterobacterales Harboring Antibiotic Resistance Genes
Previous Article in Journal
Time-Kill Analysis of Canine Skin Pathogens: A Comparison of Pradofloxacin and Marbofloxacin
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 12
Issue 10
10.3390/antibiotics12101549
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Analysis of Virulence Determinants and Genomic Islands of blaNDM-1-Producing Enterobacter hormaechei Clinical Isolates from Greece

by
Angeliki Mavroidi
1,
Konstantina Gartzonika
2,
Nick Spanakis
3,
Elisavet Froukala
3,
Christos Kittas
2,
Georgia Vrioni
3 and
Athanasios Tsakris
3,*
1
Department of Microbiology, General University Hospital of Patras, 26504 Patras, Greece
2
Department of Microbiology, Medical School, Ioannina University Hospital, 45110 Ioannina, Greece
3
Department of Microbiology, Medical School, University of Athens, 11527 Athens, Greece
*
Author to whom correspondence should be addressed.
Antibiotics 2023, 12(10), 1549; https://doi.org/10.3390/antibiotics12101549
Submission received: 23 September 2023 / Revised: 12 October 2023 / Accepted: 16 October 2023 / Published: 18 October 2023
(This article belongs to the Special Issue Antimicrobial Resistance Genes: Spread and Evolution)
Browse Figures
Review Reports Versions Notes

Abstract

:
Nosocomial outbreaks of multidrug-resistant (MDR) Enterobacter cloacae complex (ECC) are often reported worldwide, mostly associated with a small number of multilocus-sequence types of E. hormaechei and E. cloacae strains. In Europe, the largest clonal outbreak of blaNDM-1-producing ECC has been recently reported, involving an ST182 E. hormaechei strain in a Greek teaching hospital. In the current study, we aimed to further investigate the genetic make-up of two representative outbreak isolates. Comparative genomics of whole genome sequences (WGS) was performed, including whole genome-based taxonomic analysis and in silico prediction of virulence determinants of the bacterial cell surface, plasmids, antibiotic resistance genes and virulence factors present on genomic islands. The enterobacterial common antigen and the colanic antigen of the cell surface were identified in both isolates, being similar to the gene clusters of the E. hormaechei ATCC 49162 and E. cloacae ATCC 13047 type strains, whereas the two strains possessed different gene clusters encoding lipopolysaccharide O-antigens. Other virulence factors of the bacterial cell surface, such as flagella, fimbriae and pili, were also predicted to be encoded by gene clusters similar to those found in Enterobacter spp. and other Enterobacterales. Secretion systems and toxin–antitoxin systems, which also contribute to pathogenicity, were identified. Both isolates harboured resistance genes to multiple antimicrobial classes, including β-lactams, aminoglycosides, quinolones, chloramphenicol, trimethoprim, sulfonamides and fosfomycin; they carried blaTEM-1, blaOXA-1, blaNDM-1, and one of them also carried blaCTXM-14, blaCTXM-15 and blaLAP-2 plasmidic alleles. Our comprehensive analysis of the WGS assemblies revealed that blaNDM-1-producing outbreak isolates possess components of the bacterial cell surface as well as genomic islands, harbouring resistance genes to several antimicrobial classes and various virulence factors. Differences in the plasmids carrying β-lactamase genes between the two strains have also shown diverse modes of acquisition and an ongoing evolution of these mobile elements.

1. Introduction

Enterobacter hormaechei and E. cloacae are members of the Enterobacter cloacae complex (ECC), frequently recognised as the causative agents of nosocomial infections, including pneumonia, urinary tract and soft-tissue infections, septicaemia, and meningitis [1]. The species/subspecies assignment of the genus Enterobacter has often been inconsistent via routine identification techniques, 16S rRNA and housekeeping genes, while the reclassification of the species/subspecies via phylogenetic studies based on whole genome DNA-DNA hybridisations and sequencing is ongoing [2,3,4,5]. Genome-wide analysis of ECC type strains, including both clinical and environmental isolates, has revealed the presence of highly conserved genes associated with the survival of ECC strains, whereas genes involved in drug resistance and virulence were present in different genomic regions [6]. The surface of the bacterial cell is composed of common components, which contribute to the cell wall rigidity and the interaction of the bacteria with the environment, host immune response and host–pathogen interactions [7,8].
ECC are members of the ESKAPE pathogens, exhibiting resistance to several, if not all, antimicrobial classes, and are therefore described as multidrug-resistant (MDR) pathogens [9,10]; https://arpsp.cdc.gov/profile/antibiotic-resistance/mdr-enterobacter (last accessed on 1 August 2023). The most common antimicrobial resistance mechanisms in such microorganisms include the acquisition of mobile genetic elements (MGEs) from other Enterobacterales and the over-expression of efflux pumps [1]. The predominant cause of resistance to β-lactams is the production of β-lactamases, which inactivate the drug by hydrolysing the β-lactam ring. Based on the Ambler scheme, β-lactamases are grouped into four classes: the class B metallo-enzymes, also known as metallo-β-lactamases (MBLs); and the serine-reactive hydrolases (classes A, C and D) [11]. ECC strains are intrinsically resistant to cephalosporins due to the production of the chromosomal cephalosporinase blaACT-16, whereas the spread of acquired β-lactamases has been often documented, such as the class A extended-spectrum beta-lactamases (ESBLs) of the CTX-M and SHV families as well as carbapenemases, such as the VIM- and NDM-types (class B), KPC-2 (class A) and OXA-48 (class D) [1,7,8,12,13].
Nosocomial outbreaks of MDR ECC are associated with a small number of MLST STs, such as ST78, ST114 and ST171, compared with the overall diversity of the STs identified in Enterobacter spp. [12,13,14]. Similarly, with other Enterobacterales, there were no obvious associations of NDM producers and STs or well-defined blaNDM-producing international clones [14]. ST182 has been sporadically reported in several countries, while it was associated with outbreaks in the Czech Republic and Greece [15,16,17]. In the Czech Republic, the first blaNDM-4-producing ST182 isolate was reported in 2012 and a clonal expansion of these isolates was observed during 2016 [15,16]. The largest outbreak in Europe of blaNDM-1-producing ST182 isolates has been recently reported in a retrospective study of carbapenem-non-susceptible ECC, during a 6-year period (March 2016–March 2022) in the University Hospital of Ioannina, Greece [17].
In the present study, we aimed to further investigate the genetic make-up of blaNDM-1-producing ECC strains from the aforementioned outbreak in Greece. For this purpose, we have compared the whole genome sequences of two blaNDM-1-producing ECC isolates, a representative of ST182 strain EC-ML559 and strain EC-ML621, which co-produced the ESBL blaCTX-M-15 and belonged to a closely related pulsed-field gel electrophoresis (PFGE) subtype to the major subtype of the clonal outbreak. In silico characterisation and comparative genomics were performed for the gene clusters encoding the major virulence components of the cell surface of Gram-negative bacteria, such as the enterobacterial common antigen (ECA), colanic acid (CA), and O-antigen lipopolysaccharide (LPS), pili, flagella and fimbriae [7,8]. Additionally, we have explored the presence of genomic islands, which carried acquired plasmids, antimicrobial resistance and virulence genes, including secretion systems and toxin–antitoxin systems [18,19,20].

2. Results

2.1. Antimicrobial Susceptibility Testing, MLST, WGS and Type (Strain) Genome Server (TYGS) Analysis

The antibiotic susceptibility testing of strains EC-ML559 and EC-ML621 has revealed that both strains were multidrug-resistant, e.g., resistant to β-lactams (including carbapenems), β-lactam/β-lactamase inhibitor combinations (amoxicillin/clavulanic, ticarcillin/clavulanic acid, piperacillin/tazobactam, acid), quinolones (ciprofloxacin, ofloxacin), trimethoprim/sulfamethoxazole, tobramycin (Table S1). The characteristics of the WGS assemblies, MLST allelic profiles, predicted genomic islands and MGEs of strains EC-ML559 and EC-ML621 are presented in Table 1. Strain EC-ML621 was assigned to the novel MLST ST2143 (available at https://pubmlst.org/organisms/enterobacter-cloacae, last accessed on 1 August 2023), differing only in one nucleotide position in the MLST alleles (dnaA-554) compared with strain EC-ML559 of ST182 (dnaA-49), as shown via MLST and WGS. The type-based species and subspecies clustering of 31 Enterobacter spp. strains yielded 25 species clusters (average branch support of 92.5%) and 26 subspecies clusters (Figure S1). The provided query strains EC-ML559 and EC-ML621 belonged to the same cluster with the E. hormaechei type strains, being closely related to the E. hormaechei subsp. xiangfangensis LMG27195 type strain; thus, they were assigned as E. hormaechei.

2.2. Enterobacterial Common Antigen (ECA), Colanic Acid (CA) and Lipopolysaccharide (LPS) O-Antigen

In strains EC-ML559 and EC-ML621, gene clusters encoding proteins involved in the biosynthesis of the enterobacterial common antigen (ECA), colanic acid (CA) and lipopolysaccharide (LPS) O-antigens were identified (Figure 1). The ECA cluster, known as the wec operon from wecA to wecG, of strains EC-ML559 and EC-ML621 was more similar to E. hormaechei ATCC49162 than E. cloacae ATCC 13047. The ECA is a carbohydrate antigen built of repeating units of three amino sugars, the structure of which is conserved throughout Enterobacterales [21]. The ECA cluster encodes proteins involved in many steps of ECA synthesis, such as the addition of ECA polysaccharide chains to a lipid carrier and the transportation of ECA to the bacterial cell surface [21]. The CA gene cluster, within the wza and galF genes, includes genes for the production of the extracellular polysaccharide CA or M antigen, and synthesis of GDP-L-fucose, one of the precursors of CA, which is widely found within E. coli, as well as within other species of the Enterobacterales [22]. The CA cluster of strains EC-ML559 and EC-ML621 was highly similar to the CA cluster of the E. cloacae ATCC13047 type strain, whereas no CA cluster was identified in the E. hormaechei ATCC49162 type strain.
The (LPS) O-antigen biosynthetic loci were also located on gene clusters within the housekeeping genes galF and gndA (Figure 2), as described in E. coli and E. cloacae isolates [23,24]. Phylogenetic analysis of the O-antigen clusters of strains EC-ML559 and EC-ML621 was performed at the NCBI database using other ECC type strains, which revealed correlations of the O-antigen clusters with the subspecies level of the type strains (Figure S2). The O-antigen biosynthetic locus of strain EC-ML559 was highly similar to the E. hormaechei strain CUVET18-121 (GenBank: CP114979.1; region: 3040072–3053366; 99.99% identity) from Thailand, the E. hormaechei strain Ehh_23 (GenBank:CP126767.1; region: 2958890–2972184; 98.98% identity) and the E. hormaechei strain Ehh_24 (GenBank: CP126761.1; region: 2960240–2973534; 98.98% identity) from Switzerland, and 70% identical to the E. coli strain RHB04-C20 (GenBank CP055786; region: 1810262–1823551). All these clusters contained the rfbB, -D, -A, -C (or rmlB, -D, -A, -C) genes, which encode proteins for dTDP-rhamnose synthesis. Analysis of the WGS assemblies in the E. cloacae PubMLST database (avalaiable at https://pubmlst.org/organisms/enterobacter-cloacae, last accessed on 1 August 2023) has shown that the E. hormaechei strain CUVET18-121 belonged to ST182, whereas strains E. hormaechei Ehh_23 and E. hormaechei Ehh_24 were of a different ST (ST245). Notably, the O-antigen biosynthetic locus of strain EC-ML621 had a different genetic structure, being highly similar to the locus of an E. hormaechei strain Eh202_LBHALD (GenBank CP129250, region: 507656–523007) from Senegal, which belonged to ST182. Moreover, by using the Mobile Element Finder tool, it was revealed that the E. hormaechei strains CUVET18-121, Ehh_23, Ehh_24 and Eh202_LBHALD carried only the chromosomal blaACT-16.

2.3. Flagella, Fimbriae and Secretion Systems

Strains EC-ML559 and EC-ML621 possessed flag-3a loci, which encode peritrichous flagellar systems, as described previously in the E. cloacae ATCC 13047 and other Enterobacter spp. isolates [25,26]. The flag-3a loci all comprised three gene blocks occurring in the following order: block 1, flhEAB-cheZYBR-fliEFGHIJKLMNOPQR; block 2, fliZACDST; and block 3, flgNMABCDEFGHIJKL-flhDC-motAB-cheAW. Both strains harboured the fimA-Z gene clusters encoding type I fimbriae (or type I pili) biosynthesis proteins [27,28,29], which were highly similar (>99.9%) to the E. hormaechei subsp. xiangfangensis strain Eh8322_LBHALD (GenBank CP129636.1, region: 1872180–1878869) from Senegal, E. hormaechei strain CUVET18-121 (CP114979, region: 1131840–1138529) from Thailand and other E. hormaechei strains. Moreover, both strains also possessed a second fimbriae gene cluster, which was also identified in the aforementioned strains. Notably, both fimbriae gene clusters were absent from the E. hormaechei ATCC 49162 type strain.
BlastN comparisons showed that the T2SS, T6SS and T4P1 secretion systems of strains EC-ML559 and EC-ML621 were more similar to the respective gene clusters of the E. hormaechei ATCC 49162 type strain than the E. cloacae 13047 type strain (Figure 3). Both strains possessed T4P1 gene clusters encoding pili biosynthesis proteins [30,31], containing the hofQ-hof-hofO-pilN-pilM genes (GenBank, JASKGQ010000009.1; region 43949–45669). Strain EC-ML621 carried two more pili biosynthesis loci; the first locus (contig JASKGQ010000016.1, region 1: 77466–86705) was similar to the E. hormaechei strain 21KM1498 plasmid pKM98_p3 (GenBank CP126866.1) and E. ludwigii_strain CM-TZ4 (GenBank CP116349.1), while the second locus (JASKGQ010000022, region: 20641–30420) was located on the blaLAP2-carrying plasmid sequence, which was similar to the E. hormaechei subsp. xiangfangensis strain 090095 pLAP2_090095 (GenBank CP128418.1).

2.4. In Silico Prediction of Genomic Islands

In the identified genomic island (Figure S3), both EC-ML559 and EC-ML621 isolates harboured β-lactamase-carrying contigs highly similar to published plasmid sequences (Figure S4), including the E. cloacae pECL_A plasmid carrying blaTEM-1 [3]. Both isolates carried the blaOXA-1 gene, but it was not possible to determine in silico whether it was located on a plasmid or on the chromosome due to the small sizes of the obtained blaOXA-1-harbouring contigs and the abundance of the presence of blaOXA-1 genes among different plasmid types and bacterial species. The blaNDM-1-harbouring contigs of both isolates were highly similar to DNA fragments of Klebsiella pneumoniae subsp. pneumoniae strain KPX plasmid pKPX-1 (GeneBank: AP012055.1) and E. hormaechei subsp. xiangfangensis strain ST114 plasmid pLAU_ENM30_NDM1 (GeneBank: MN792917.1) [32,33] Of note, the genetic background of blaNDM-1 differed in strains EC-ML559 and EC-ML621 (Figures S4 and S5). Moreover, strain EC-ML621 also possessed blaCTX-M-14-, blaCTX-M-15- and blaLAP-2-harbouring plasmid sequences (Table 1, Figure S4). The blaCTX-M-14 gene was carried on a contig (JASKGQ010000024.1), which also possessed an IncL/M genetic element, showing 99.6% similarity to the repC, repB and the replication initiation protein-encoding gene repA of the K. pneumoniae plasmid pMU407.1 (GenBank: U27345.1). The blaLAP-2 genetic background was highly similar to the E. hormaechei subsp. xiangfangensis strain 090095 plasmid pLAP2_090095, GenBank Accession no. CP128418, WGS: JAHEVO010000000) and other Enterobacterales [34,35,36]. The latter strain belonged to a different MLST ST (ST175) and shared 33 out of 53 rMLST loci compared to strain EC-ML621, as shown by analysing the WGS assembly in the E. cloacae PubMLST database (avalaiable at https://pubmlst.org/organisms/enterobacter-cloacae, last accessed on 1 August 2023).
In the predicted genomic islands, besides the presence of β-lactamase genes, both strains EC-ML559 and EC-ML621 carried resistance genes to several antimicrobial classes (e.g., aminoglycosides, chloramphenicol, quinolones, fosfomycin, trimethoprim, sulfonamides, macrolides, quaternary ammonium compounds), heavy-metal resistance genes to arsenic copper/silver and tellurium and genes associated with bacterial virulence, including type II and type IV toxin–antitoxin systems, secretion systems (T1SS, T2SS, T4SS, T6SS and T4P), flagellin and fimbriae genes (Table 1). The origin of the transfer site (oriT), the relaxase gene, the gene encoding type IV coupling protein (T4CP) and the gene clusters for bacterial type IV secretion system (T4SS) identified in strains EC-ML559 and EC-ML621 are presented in Table 2.

3. Materials and Methods

3.1. Bacterial Isolates, Antimicrobial Susceptibility Testing, Whole Genome Sequencing (WGS) and MLST Typing

Susceptibility testing of strains EC-ML559 and EC-ML621 to various antimicrobial agents has been performed using the automated system MicroScan (Beckman Coulter, Brea, CA, USA), according to the manufacturer’s instructions. The categorisation of the isolates and interpretation of the minimum inhibition concentrations (MICs) of the antimicrobial agents were based on the European Committee on Antimicrobial Susceptibility Testing (EUCAST) criteria and breakpoints for Enterobacterales (available at: https://www.eucast.org/fileadmin/src/media/PDFs/EUCAST_files/Breakpoint_tables/v_13.0_Breakpoint_Tables.pdf, last accessed on 1 August 2023). The genome sequence of the representative ST182 outbreak strain EC-ML559 has been presented in a previous study [17]; (GenBank Accession no. JARUPS000000000.1). Strain EC-ML621 was cultured in LB medium and allowed to grow until the mid-logarithmic phase at 30 °C in a shaking incubator. DNA was extracted using a genomic DNA isolation kit (Promega, Madison, WI, USA). The concentration and quality of the DNA were determined via Nanodrop 2000 (Thermo Fisher Scientific, Dreieich, Germany). Multilocus-sequence typing (MLST) was performed via the amplification of the seven housekeeping genes, using primer sequences described in E. cloacae MLST database (PubMLST web site https://pubmlst.org/ecloacae, last accessed on 1 August 2023)). The isolate EC-ML621, the novel allele and ST identified were deposited in the PubMLST database under PubMLST id 1364 [37]; (https://pubmlst.org/bigsdb?page=info&db=pubmlst_ecloacae_isolates&id=1364, last accessed on 1 August 2023). Moreover, WGS of isolate EC-ML621 was performed on an Illumina platform (Eurofins Genomics Europe Sequencing GmbH, Konstanz, Germany). After sequence cleaning and normalisation, high-quality reads were assembled using SPAdes [38].

3.2. Whole-Genome-Based Taxonomic Analysis

Whole-genome-based taxonomic analysis of the ECC strains EC-ML559 and EC-ML621 was performed by using the free bioinformatics platform Type (Strain) Genome Server (TYGS) available at https://tygs.dsmz.de (last accessed on 23 July 2023) [39], while information on nomenclature, synonymy and associated taxonomic literature was provided by the List of Prokaryotic names with Standing in Nomenclature (LPSN, available at https://lpsn.dsmz.de; last accessed on 23 July 2023) [40]. Briefly, the TYGS analysis included the determination of closely related type strains, pairwise comparison of genome sequences and phylogenetic inference via the MASH algorithm. The resulting intergenomic distances were used to infer a balanced minimum evolution tree with branch support via FASTME 2.1.6.1, including SPR post processing. Branch support was inferred from 100 pseudo-bootstrap replicates each. The trees were rooted at the midpoint and visualised with PhyD3. The type-based species clustering using a 70% dDDH radius around each of the 28 type strains, and subspecies clustering was carried out using a 79% dDDH threshold.

3.3. In Silico Prediction of Genomic Islands, Mobile Elements, Antimicrobial Resistance and Virulence Genes

Genomic island prediction was performed for strains EC-ML559 and EC-ML621 by using the IslandViewer4 web application (available at: https://www.pathogenomics.sfu.ca/islandviewer; last accessed on 10 August 2023), which includes all genomic islands found on published sequenced bacterial and archaeal genomes that have been predicted using the currently most accurate GI prediction methods, including IslandPick, SIGI-HMM, IslandPath-DIMOB [41]. The E. hormaechei subsp. xiangfangensis LMG27195 type strain genome sequence was used as a reference sequence for the alignment of each genome sequence. The presence of mobile elements, antimicrobial resistance and virulence genes was predicted by using the Mobile Element Finder, Plasmid Finder. ResFinder and Virulence Finder tools available at: https://www.genomicepidemiology.org/services/ (last accessed on 10 August 2023), Center for Genomic Epidemiology, Technical University of Denmark, Denmark [42].
The Artemis Comparison Tool version 18.2.0 [43] and SnapGene version 5.1.5 (GSL Biotech LLC., San Diego, CA, USA) were used to manipulate and visualise the sequences studied. The oriTfinder tool [44] was used to explore the presence of conjugative regions of the self-transmissible MGEs: the origin of transfer site (oriT), the relaxase gene, the gene encoding type IV coupling protein (T4CP) and the gene cluster for bacterial type IV secretion system (T4SS). BlastN comparisons of the beta-lactamase-carrying contigs of the isolates with published plasmid sequences were performed by using the BLAST Ring Image Generator (BRIG) version 0.95 software (available at: https://brig.sourceforge.net/ (last accessed on 23 August 2023)) [45]. BlastN comparisons and visualisations of the gene clusters encoding ECA, CA, O-antigens and secretion systems of the ECC strains EC-ML559 and EC-ML621 were performed by using the Easyfig version 2.2.5 software [46], (available at: http://mjsull.github.io/Easyfig/ (last accessed on 23 August 2023)).

4. Discussion

The genus Enterobacter comprises species/subspecies often recognised as human, animal, and plant pathogens [1]. The difficulties in the species/subspecies assignment of the genus Enterobacter spp. and the definition of international clones among the diverse population structure of ECC strains are challenging for the diagnosis of infections caused by these pathogens and outbreak analysis [2,3,4,5]. Nosocomial infections caused mainly by MDR E. cloacae and E. hormaechei are frequently reported and are difficult to treat. A recent investigation of an outbreak in new-borns with E. cloacae complex isolates revealed, via WGS, E. bugandensis as the causative agent, and highlighted the need for better discrimination of Enterobacter species inside the ECC [47]. Moreover, the arsenal of these pathogens to acquire antimicrobial resistance determinants and virulence factors via mobile genetic elements from other Gram-negative bacteria further complicates the diagnosis and treatment of infections caused by ECC. In the present study, we have performed a comprehensive whole genome analysis of two representative E. hormaechei blaNDM-1-producing isolates (strains EC-ML559 and EC-ML621) recovered during a large outbreak in a Greek teaching hospital [17]. The isolates, which belonged to the closely related MLST STs 182 and 2143, were assigned as E. hormaechei via whole genome-based taxonomic analysis. In silico prediction and comparison of bacterial cell surface components and other systems involved in the pathogenesis of this microorganism were performed, including genomic islands, which often carry antimicrobial resistance genes and virulence factors, such as those involved in bacterial secretion, host determination and the colonisation of different strains.
Antigenic and phase variation mechanisms have been proposed to enable immune evasion during chronic or recurrent infection and/or the generation of variants with altered ability to colonise different niches in the host [8]. Variation in all major surface antigens is a hallmark of pathogens persistently colonising mucosal surfaces. Polysaccharides are important constituents of the bacterial cell surface and prominent antigens, including the enterobacterial common antigen, colanic acid or M-antigen and the lipopolysaccharide of Gram-negative bacteria [7,8]. Gene clusters involved in the biosynthesis of the ECA and CA were identified in both strains, EC-ML559 and EC-ML621, being similar to the respective gene clusters of the E. hormaechei ATCC 49162 and E. cloacae ATCC 13047 type strains. Gram-negative bacteria also produce LPS, which is composed of an O-antigen, an outer core, and an inner core all joined by covalent bonds [23,24]. The O-antigens are structurally and serologically diverse, allowing for escape from the immune system of the host. Additionally, antisera against O-antigens provide a basis for serotyping schemes in many Gram-negative bacteria and a basic tool utilised in outbreak investigations and epidemiological surveillance. The LPS O-antigens of E. coli have been well characterised [23], and, in a recent study, the molecular characterisation of E. cloacae O-antigens has been performed [24]. In the present study, the E. hormaechei strains EC-ML559 and EC-ML621 were found to possess different types of O-antigen clusters, although they belonged to closely related MLST STs, differing only in one MLST allele by one nucleotide difference. Moreover, the same O-antigen cluster of strain EC-ML559 was identified in E. hormaechei strains of different STs (ST182 and ST245), sharing similarities with a gene cluster found in an E. coli strain, suggesting the horizontal spread of these gene clusters via recombination and evolution in different strains of these species.
Among the components of the bacterial cell surface also associated with pathogenicity are filamentous organelles, such as flagella, pili and fimbriae [7,8]. These structures are involved in bacterial motility and adherence, biofilm formation and colonisation on infected hosts, acting as virulence factors. The bacterial flagellum consists of a filament core, which is preserved across bacterial species, while the outer domains exhibit high variability, and, in some cases, are even completely absent [25]. The Enterobacterales encode five distinct flagellar systems and the genes required for the assembly, maintenance and functioning of this flagellum are organised in five genomic clusters (flag-1 to flag-5) [26]. Both strains, EC-ML559 and EC-ML621, carried flag-3a peritrichous flagellar loci, which are restricted to members of three genera: Enterobacter, Erwinia, and Pantoea. All Gram-negative bacterial species that have been examined to date were found to produce one or more types of fimbriae with adhesive properties, which are involved in the specific adherence to receptors present on certain host cell types and the colonisation of the host [27,28,29]. In a previous study, at least 16 fimbrial gene clusters were identified in two enterohaemorrhagic E. coli O157:H7, and six of them were detected only in the E. coli O157 and O145 serotypes tested, while the fimB-H operon has been previously characterised in E. coli K12 strain and uropathogenic E coli [7,28,29]. Strains EC-ML559 and EC-ML621 also carried fimbrial genes, which were found in several E. hormaechei strains, which may contribute to pathogenesis. Type IV pili are also filamentous structures of the bacterial cell surface widespread among the β-, γ-, and δ-proteobacteria and the cyanobacteria, which contribute to motility and act as virulence factors [30,31]. In strain EC-ML621, two of the three pili biosynthesis gene clusters that were identified showed similarities with published plasmid sequences, suggesting that these clusters can be transferred horizontally between different strains or species.
Secretion system proteins aid in the transfer of a variety of virulence determinants across the bacterial cell membrane [18]. Because of the remarkable structural and functional similarities, secretion systems may be exploited for the development of novel antibacterial compounds. Recent studies have shown that secretion systems, including T6SS accessory and core components, T4SS, and multidrug resistance genes/efflux system genes seemed vital for the survival of bacteria in various environmental niches, such as humans and plants. The type IV secretion system (T4SS) is evolutionarily related to the conjugation system of bacteria and it was speculated that the low GC content and the presence of T4SS may be linked with horizontal gene transfer (HGT), which helps bacteria adapt to environmental changes and acquire antibiotic resistance. In a recent study of 952 K. pneumoniae strains, the presence of one conserved T2SS, one conserved T5SS, and two conserved T6SS was detected in more than 90% of the strains, whereas T1SS and T4SS were enriched in the hypervirulent and classical multidrug resistance pathotypes of K. pneumoniae, respectively [19]. In ECC, T2SS genes were absent in the E. kobei, E. hormaechei, and E. ludwigii, type strains previously studied, but there is a lack of experimental data to link the absence of T2SS with the pathogenicity of these strains. In the present study, the T1SS, T2SS, T4SS, T6SS and T4P secretion systems were identified in genomic islands present in both studied isolates, showing similarities with the E. hormaechei ATCC 49162 and E. cloacae ATCC 13047 type strains, whereas Type IV pili biosynthesis gene clusters (T4SS) were located on plasmid sequences of strain EC-ML621. Furthermore, both strains possessed toxin–antitoxin systems, which are abundant in bacterial genomes, in which they move via horizontal gene transfer. Specific toxin–antitoxin systems have been involved in the protection against phages, and may be involved in pathogenicity [20].
Finally, both isolates possessed plasmid sequences that carried antimicrobial resistance and virulence genes. In regard to β-lactamase genes, the blaOXA-1, blaTEM-1-carrying E. cloacae pECLA-type and the blaNDM-1 K. pneumoniae pKPX-1-type plasmid sequences were found in both strains, whereas strain EC-ML621 also carried blaCTX-M-14-, blaCTX-M-15- and blaLAP-2-harbouring plasmid sequences. The blaCTX-M-14- and blaCTX-M-15-harbouring plasmids are ubiquitous in Enterobacterales. The Ambler class A LAP-1 β-lactamase was first characterised in 2007 [34], and shows a narrow hydrolysis spectrum against β-lactams inhibited by clavulanic acid. The blaLAP-1 gene was located on the same plasmid with the quinolone resistance gene qnrS1, which was identified in E. cloacae isolates from France and Vietnam. LAP-2 is a variant of LAP-1, which was first identified in an E. cloacae clinical isolate from China in 2008 [35]. Enterobacter spp. strains carrying the blaLAP-1 gene have been reported in Spain, Algeria and Korea, while blaLAP-2 has been detected in species of the genera Klebsiella, Enterobacter and Salmonella from China, Tunisia, Norway and the Netherlands [36].
The spread of antimicrobial resistance and virulence genes is often mediated by MGEs, which are able to move within or between DNA molecules (e.g., insertion sequences, transposons and gene cassettes/integrons), and genetic elements that are able to transfer between bacterial cells, such as plasmids and integrative conjugative elements [48]. In a recent review analysis of NDM-positive strains recovered worldwide, blaNDM has been reported to be carried on plasmids (n = 355) with a variety of 20 replicon types in the Enterobacterales, a finding that highlights multiple acquisitions of blaNDM by various plasmids via horizontal transfer [14]. blaNDM-1 was located mostly on incompatibility group (Inc) replicon types IncA/C, ColE, IncFIA, IncFIB, IncF IncII, IncHI1, IncHI3, IncHIB, IncL/M, IncN2, IncP, IncR, IncT, IncX1, IncX3, and IncY. In K. pneumoniae, a comparative genomic analysis of carbapenemase-harbouring plasmids from a European survey has revealed three modes of acquisition and transmission of antibiotic resistance genes via plasmids, one plasmid/multiple lineages, multiple plasmids/multiple lineages, and multiple plasmids/one lineage [49]. In the present study, both studied strains belonged to closely related MLSTs and carried IncFII (pECL_A)-like carrying blaTEM-1, IncFII(pKPX1)-type-carrying blaNDM-1 and IncR plasmid sequences. Strain EC-ML559 harboured a Col440I-like carrying the quinolone resistance gene qnrB19. Strain EC-ML621 also harboured other plasmidic sequences found in other Enterobacterales, such as IncL/M (pMU407.1) carrying blaCTX-M-14, IncFII(Yp), IncFIB(pB171) and IncN [48,49]. For example, the IncL/M (pMU407.1) and IncFII replicon types have also been associated with blaOXA-48 and blaKPC, respectively. These findings highlight the potential of the dispersion of these replicon type plasmids, which often carry other antimicrobial resistance and virulence genes, among different species in Enterobacterales, including species of the genus Enterobacter.

5. Conclusions

In the present study, a comprehensive analysis of the WGS assemblies of two representative blaNDM-1-producing outbreak isolates from Greece has shown that these isolates possess components of the bacterial cell surface, such as the enterobacterial common antigen, the colanic antigen, LPS-O-antigen, flagella, fimbriae and pili, which play a role in the pathogenicity of Gram-negative bacteria. Moreover, the presence of genomic islands, harbouring resistance genes to several antimicrobial classes and virulence factors, such as secretion systems and toxin–antitoxin systems, was revealed. Differences in the plasmids carrying β-lactamase genes between the two strains revealed diverse modes of acquisition and ongoing evolution of these mobile elements. In conclusion, WGS analysis of ECC strains would facilitate the development of methods for accurate identification and molecular serotyping schemes. Furthermore, WGS and comparative genomics of clinical and environmental ECC would elucidate the evolution, epidemiology, prevention and treatment of infectious disease caused by these broad-host-range niche-associated species.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antibiotics12101549/s1, Table S1: Antimicrobial susceptibility testing of strains EC-ML559 and EC-ML621; Figure S1: Tree inferred with FastME 2.1.6.1. from GBDP distances calculated from genome sequences; Figure S2: Phylogenetic analysis of the O-antigen clusters of strains EC-ML559 and EC-ML621 compared with other Enterobacter cloacae complex type strains; Figure S3: Predicted genomic islands and virulence/resistance gene annotations (when available) are colored within the images based on the legend; Figure S4: BlastN comparisons of β-lactamase harboring sequences of strains EC-ML559 and EC-ML621 identified in silico with those published plasmid sequences; Figure S5: Genetic background of blaNDM-1-harbouring contigs of strains EC-ML559 and EC-ML621.

Author Contributions

Conceptualisation, A.M., G.V. and A.T.; methodology, A.M., K.G., N.S., E.F., C.K. and G.V.; formal analysis-investigation, A.M., K.G., N.S., E.F., C.K. and G.V.; writing—original draft preparation A.M. and A.T.; review and editing, A.M., K.G., G.V. and A.T.; supervision, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The nucleotide sequences of the contigs of isolate EC-ML621 have been deposited in the Genome collection of Enterobacter cloacae PubMLST database, the Joint Genome Institute (JGI), Berkeley, USA, database under analysis Project GOLD ID Ga059773 and the National Center for Biotechnology Information (NCBI), National Institutes of Health, USA, database DDBJ/ENA/GenBank under the accession JASKGQ000000000. The version described in this paper is version JASKGQ010000000. The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Davin-Regli, A.; Lavigne, J.P.; Pagès, J.M. Enterobacter spp.: Update on taxonomy, clinical aspects, and emerging antimicrobial resistance. Clin. Microbiol. Rev. 2019, 32, e00002-19. [Google Scholar] [CrossRef] [PubMed]
  2. Paauw, A.; Caspers, M.P.M.; Schuren, F.H.J.; Hall, M.A.L.-V.; Delétoile, A.; Montijn, R.C.; Verhoef, J.; Fluit, A.C. Genomic diversity within the Enterobacter cloacae complex. PLoS ONE 2008, 3, e3018. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, W.-Y.; Wong, C.-F.; Chung, K.M.-K.; Jiang, J.-W.; Leung, F.C.-C. Comparative genome analysis of Enterobacter cloacae. PLoS ONE 2013, 8, e74487. [Google Scholar] [CrossRef] [PubMed]
  4. Ciufo, S.; Kannan, S.; Sharma, S.; Badretdin, A.; Clark, K.; Turner, S.; Brover, S.; Schoch, C.L.; Kimchi, A.; DiCuccio, M. Using average nucleotide identity to improve taxonomic assignments in prokaryotic genomes at the NCBI. Int. J. Syst. Evol. Microbiol. 2018, 68, 2386–2392. [Google Scholar] [CrossRef] [PubMed]
  5. Federhen, S.; Rossello-Mora, R.; Klenk, H.-P.; Tindall, B.J.; Konstantinidis, K.T.; Whitman, W.B.; Brown, D.; Labeda, D.; Ussery, D.; Garrity, G.M.; et al. Meeting report: GenBank microbial genomic taxonomy workshop (12–13 May, 2015). Stand. Genom. Sci. 2016, 11, 15. [Google Scholar] [CrossRef]
  6. Mustafa, A.; Ibrahim, M.; Rasheed, M.A.; Kanwal, S.; Hussain, A.; Sami, A.; Ahmed, R.; Bo, Z. Genome-wide analysis of four Enterobacter cloacae complex type strains: Insights into virulence and niche adaptation. Sci. Rep. 2020, 10, 8150. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, J.; Ma, W.; Wang, X. Insights into the structure of Escherichia coli outer membrane as the target for engineering microbial cell factories. Microb. Cell Fact. 2021, 20, 73. [Google Scholar] [CrossRef] [PubMed]
  8. Van der Woude, M.W.; Bäumler, A.J. Phase and antigenic variation in bacteria. Clin. Microbiol. Rev. 2004, 17, 581–611. [Google Scholar] [CrossRef] [PubMed]
  9. De Oliveira, D.M.P.; Forde, B.M.; Kidd, T.J.; Harris, P.N.A.; Schembri, M.A.; Beatson, S.A.; Paterson, D.L.; Walker, M.J. Antimicrobial resistance in ESKAPE pathogens. Clin. Microbiol. Rev. 2020, 33, e00181-19. [Google Scholar] [CrossRef]
  10. Annavajhala, M.K.; Gomez-Simmonds, A.; Uhlemann, A.C. Multidrug-resistant Enterobacter cloacae complex emerging as a global, diversifying threat. Front. Microbiol. 2019, 10, 44. [Google Scholar] [CrossRef]
  11. Hall, B.G.; Barlow, M. Revised Ambler classification of β-lactamases. J. Antimicrob. Chemother. 2005, 55, 1050–1051. [Google Scholar] [CrossRef] [PubMed]
  12. Peirano, G.; Matsumura, Y.; Adams, M.D.; Bradford, P.; Motyl, M.; Chen, L.; Kreiswirth, B.N.; Pitout, J.D. Genomic epidemiology of global carbapenemase-producing Enterobacter spp., 2008–2014. Emerg. Infect. Dis. 2018, 24, 1010–1019. [Google Scholar] [CrossRef] [PubMed]
  13. Izdebski, R.; Baraniak, A.; Herda, M.; Fiett, J.; Bonten, M.J.; Carmeli, Y.; Goossens, H.; Hryniewicz, W.; Brun-Buisson, C.; Gniadkowski, M.; et al. MLST reveals potentially high-risk international clones of Enterobacter cloacae. J. Antimicrob. Chemother. 2015, 70, 48–56. [Google Scholar] [CrossRef] [PubMed]
  14. Wu, W.; Feng, Y.; Tang, G.; Qiao, F.; McNally, A.; Zong, Z. NDM Metallo-β-lactamases and their bacterial producers in health care settings. Clin. Microbiol. Rev. 2019, 32, e00115-18. [Google Scholar] [CrossRef] [PubMed]
  15. Papagiannitsis, C.C.; Studentova, V.; Chudackova, E.; Bergerova, T.; Hrabak, J.; Radej, J.; Novak, I. Identification of a New Delhi metallo-β-lactamase-4 (NDM-4)-producing Enterobacter cloacae from a Czech patient previously hospitalized in Sri Lanka. Folia Microbiol. 2013, 58, 547–549. [Google Scholar] [CrossRef] [PubMed]
  16. Paskova, V.; Medvecky, M.; Skalova, A.; Chudejova, K.; Bitar, I.; Jakubu, V.; Bergerova, T.; Zemlickova, H.; Papagiannitsis, C.C.; Hrabak, J. Characterization of NDM-encoding plasmids from Enterobacteriaceae recovered From Czech Hospitals. Front. Microbiol. 2018, 9, 1549. [Google Scholar] [CrossRef] [PubMed]
  17. Gartzonika, K.; Politi, L.; Mavroidi, A.; Tsantes, A.G.; Spanakis, N.; Priavali, E.; Vrioni, G.; Tsakris, A. High prevalence of clonally related ST182 NDM-1-producing Enterobacter cloacae complex clinical isolates in Greece. Int. J. Antimicrob. Agents 2023, 62, 106837. [Google Scholar] [CrossRef] [PubMed]
  18. Costa, T.R.D.; Felisberto-Rodrigues, C.; Meir, A.; Prevost, M.S.; Redzej, A.; Trokter, M.; Waksman, G. Secretion systems in Gram-negative bacteria: Structural and mechanistic insights. Nat. Rev. Microbiol. 2015, 13, 343–359. [Google Scholar] [CrossRef]
  19. Yang, M.; Zhou, X.; Bao, Y.; Zhang, Y.; Liu, B.; Gan, L.; Tao, W.; Tuo, J.; Gong, H. Comprehensive genomic analysis reveals extensive diversity of type I and Type IV secretion systems in Klebsiella pneumoniae. Curr. Microbiol. 2023, 80, 270. [Google Scholar] [CrossRef]
  20. Unterholzner, S.J.; Poppenberger, B.; Rozhon, W. Toxin-antitoxin systems: Biology, identification, and application. Mob. Genet. Elements 2013, 3, e26219. [Google Scholar] [CrossRef]
  21. Rai, A.K.; Mitchell, A.M. Enterobacterial Common Antigen: Synthesis and function of an enigmatic molecule. mBio 2020, 11, e01914-20. [Google Scholar] [CrossRef] [PubMed]
  22. Stevenson, G.; Andrianopoulos, K.; Hobbs, M.; Reeves, P.R. Organization of the Escherichia coli K-12 gene cluster responsible for production of the extracellular polysaccharide colanic acid. J. Bacteriol. 1996, 178, 4885–4893. [Google Scholar] [CrossRef] [PubMed]
  23. DebRoy, C.; Roberts, E.; Fratamico, P.M. Detection of O antigens in Escherichia coli. Anim. Health Res. Rev. 2011, 12, 169–185. [Google Scholar] [CrossRef] [PubMed]
  24. Li, Y.; Huang, J.; Wang, X.; Xu, C.; Han, T.; Guo, X. Genetic characterization of the O-Antigen and development of a molecular serotyping scheme for Enterobacter cloacae. Front. Microbiol. 2020, 11, 727. [Google Scholar] [CrossRef]
  25. Nedeljković, M.; Sastre, D.E.; Sundberg, E.J. Bacterial flagellar filament: A supramolecular multifunctional nanostructure. Int. J. Mol. Sci. 2021, 22, 7521. [Google Scholar] [CrossRef] [PubMed]
  26. De Maayer, P.; Pillay, T.; Coutinho, T.A. Flagella by numbers: Comparative genomic analysis of the supernumerary flagellar systems among the Enterobacterales. BMC Genom. 2020, 21, 670. [Google Scholar] [CrossRef] [PubMed]
  27. Mol, O.; Oudega, B. Molecular and structural aspects of fimbriae biosynthesis and assembly in Escherichia coli. FEMS Microbiol. Rev. 1996, 19, 25–52. [Google Scholar] [CrossRef]
  28. Low, A.S.; Holden, N.; Rosser, T.; Roe, A.J.; Constantinidou, C.; Hobman, J.L.; Smith, D.G.; Low, J.C.; Gally, D.L. Analysis of fimbrial gene clusters and their expression in enterohaemorrhagic Escherichia coli O157:H7. Environ. Microbiol. 2006, 8, 1033–1047. [Google Scholar] [CrossRef]
  29. Schwan, W.R. Regulation of fim genes in uropathogenic Escherichia coli. World J. Clin. Infect. Dis. 2011, 1, 17–25. [Google Scholar] [CrossRef]
  30. Mattick, J.S. Type IV pili and twitching motility. Annu. Rev. Microbiol. 2002, 56, 289–314. [Google Scholar] [CrossRef]
  31. Nudleman, E.; Kaiser, D. Pulling together with type IV pili. J. Mol. Microbiol. Biotechnol. 2004, 7, 52–62. [Google Scholar] [CrossRef] [PubMed]
  32. Huang, T.W.; Chen, T.L.; Chen, Y.T.; Lauderdale, T.L.; Liao, T.L.; Lee, Y.T.; Chen, C.P.; Liu, Y.M.; Lin, A.C.; Chang, Y.H.; et al. Copy number change of the NDM-1 sequence in a multidrug-resistant Klebsiella pneumoniae clinical isolate. PLoS ONE 2013, 8, e62774. [Google Scholar] [CrossRef] [PubMed]
  33. Merhi, G.; Amayri, S.; Bitar, I.; Araj, G.F.; Tokajian, S. Whole genome-based characterization of multidrug resistant Enterobacter and Klebsiella aerogenes isolates from Lebanon. Microbiol. Spectr. 2023, 1, e0291722. [Google Scholar] [CrossRef] [PubMed]
  34. Poirel, L.; Cattoir, V.; Soares, A.; Soussy, C.J.; Nordmann, P. Novel Ambler class A b-lactamase LAP-1 and its association with the plasmid-mediated quinolone resistance determinant QnrS1. Antimicrob. Agents Chemother. 2007, 51, 631–637. [Google Scholar] [CrossRef] [PubMed]
  35. Huang, Z.; Mi, Z.; Wang, C. A novel b-lactamase gene, LAP-2, co-produced by an Enterobacter cloacae clinical isolate in China. J. Hosp. Infect. 2008, 70, 95–96. [Google Scholar] [CrossRef] [PubMed]
  36. Kocsis, B.; Kocsis, E.; Fontana, R.; Cornaglia, G.; Mazzariol, A. Identification of blaLAP-2 and qnrS1 genes in the internationally successful Klebsiella pneumoniae ST147 clone. J. Med. Microbiol. 2013, 62, 269–273. [Google Scholar] [CrossRef] [PubMed]
  37. Jolley, K.A.; Bray, J.E.; Maiden, M.C.J. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018, 3, 124. [Google Scholar] [CrossRef] [PubMed]
  38. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  39. Meier-Kolthoff, J.P.; Göker, M. TYGS is an automated high-throughput platform for state-of-the-art genome-based taxonomy. Nat. Commun. 2019, 10, 2182. [Google Scholar] [CrossRef]
  40. Meier-Kolthoff, J.P.; SardàCarbasse, J.; Peinado-Olarte, R.L.; Göker, M. TYGS and LPSN: A database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic Acids Res. 2022, 50, D801–D807. [Google Scholar] [CrossRef]
  41. Bertelli, C.; Laird, M.R.; Williams, K.P.; Simon Fraser University Research Computing Group; Lau, B.Y.; Hoad, G.; Winsor, G.L.; Brinkman, F.S.L. IslandViewer 4: Expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 2017, 45, W30–W35. [Google Scholar] [CrossRef]
  42. Johansson, M.H.K.; Bortolaia, V.; Tansirichaiya, S.; Aarestrup, F.M.; Roberts, A.P.; Petersen, T.N. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: Mobile Element Finder. J. Antimicrob. Chemother. 2021, 76, 101–109. [Google Scholar] [CrossRef] [PubMed]
  43. Carver, T.J.; Rutherford, K.M.; Berriman, M.; Rajandream, M.A.; Barrell, B.G.; Parkhill, J. ACT: The Artemis Comparison Tool. Bioinformatics 2005, 21, 3422–3423. [Google Scholar] [CrossRef] [PubMed]
  44. Li, X.; Xie, Y.; Liu, M.; Tai, C.; Sun, J.; Deng, Z.; Ou, H.Y. oriTfinder: A web-based tool for the identification of origin of transfers in DNA sequences of bacterial mobile genetic elements. Nucleic Acids Res. 2018, 46, W229–W234. [Google Scholar] [CrossRef] [PubMed]
  45. Alikhan, N.F.; Petty, N.K.; Ben Zakour, N.L.; Beatson, S.A. BLAST Ring Image Generator (BRIG): Simple prokaryote genome comparisons. BMC Genom. 2011, 12, 402. [Google Scholar] [CrossRef] [PubMed]
  46. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A genome comparison visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef] [PubMed]
  47. Girlich, D.; Ouzani, S.; Emeraud, C.; Gauthier, L.; Bonnin, R.A.; Le Sache, N.; Mokhtari, M.; Langlois, I.; Begasse, C.; Arangia, N.; et al. Uncovering the novel Enterobacter cloacae complex species responsible for septic shock deaths in newborns: A cohort study. Lancet Microbe 2021, 2, e536–e544. [Google Scholar] [CrossRef]
  48. Partridge, S.R.; Kwong, S.M.; Firth, N.; Jensen, S.O. Mobile genetic elements associated with antimicrobial resistance. Clin. Microbiol. Rev. 2018, 31, e00088-17. [Google Scholar] [CrossRef]
  49. David, S.; Cohen, V.; Reuter, S.; Sheppard, A.E.; Giani, T.; Parkhill, J.; European Survey of Carbapenemase-Producing Enterobacteriaceae (EuSCAPE) Working Group; ESCMID Study Group for Epidemiological Markers (ESGEM); Rossolini, G.M.; Feil, E.J.; et al. Integrated chromosomal and plasmid sequence analyses reveal diverse modes of carbapenemase gene spread among Klebsiella pneumoniae. Proc. Natl. Acad. Sci. USA 2020, 117, 25043–25054. [Google Scholar] [CrossRef]
Antibiotics 12 01549 g001
Figure 1. BlastN comparisons of the ECA and CA gene clusters of strains EC-ML559 and EC-ML621.
Figure 1. BlastN comparisons of the ECA and CA gene clusters of strains EC-ML559 and EC-ML621.
Antibiotics 12 01549 g001
Antibiotics 12 01549 g002
Figure 2. BlastN comparisons of the LPS O-antigen biosynthetic loci of strains EC-ML559 and EC-ML621.
Figure 2. BlastN comparisons of the LPS O-antigen biosynthetic loci of strains EC-ML559 and EC-ML621.
Antibiotics 12 01549 g002
Antibiotics 12 01549 g003
Figure 3. BlastN comparisons of T2SS, T4SS, T6SS and T4P1 secretion systems of strains EC-ML559 and EC-ML621.
Figure 3. BlastN comparisons of T2SS, T4SS, T6SS and T4P1 secretion systems of strains EC-ML559 and EC-ML621.
Antibiotics 12 01549 g003
Table 1. Characteristics of the WGS assemblies and MLST allelic profiles, antibiotic resistance and virulence genes in the predicted genomic islands of strains EC-ML559 and EC-ML621.
Table 1. Characteristics of the WGS assemblies and MLST allelic profiles, antibiotic resistance and virulence genes in the predicted genomic islands of strains EC-ML559 and EC-ML621.
Characteristics #EC-ML559EC-ML621
WGS AssemblyGenbank Accession no. JARUPS000000000.1
Contigs: 74; total length: 4,958,007 bp,
%GC: 55.07		Genbank Accession no. JASKGQ000000000
Contigs: 100, total length: 5,153,338 bp,
%GC: 55.02
MLST ST
(allelic profile)		ST182
(49-20-19-44-90-24-32)		ST2143
(554-20-19-44-90-24-32)
PlasmidsIncFII(pECL_A), IncR
IncFII(pKPX1), Col440I		IncFII(pECLA), IncR, IncFII(Yp), IncL/M(pMU407), IncFIB(pB171), IncN
ARGsβ-lactams: blaACT-16, blaNDM-1, blaTEM-1, blaOXA-1
aminoglycosides: aph(6)-Id, aac(6′)-Ib3, aph(3′)-Ia, aac(6′)-Ib-cr, aph(3″)Ib
chloramphenicol: catB3
fosfomycin: fosA
quaternary ammonium compounds: qacE
quinolones: OqxA, OqxB, qnrB19
trimethoprim: dfrA14
sulfonamides: sul1, sul2
macrolides: mrx(A), mphR(A), mph(A)
rifampicin: ARR-3
tetracyclines: tet(D)		β-lactams: blaACT-16, blaNDM-1, blaTEM-1,
blaOXA-1, blaCTX-M-14, blaCTX-M-15, blaLAP-2
aminoglycosides: aph(6)-Id, aac(3)-IId, aadA2, aac(6′)-Ib-cr, aph(3″)-Ib
chloramphenicol: catB3
fosfomycin: fosA
quaternary ammonium compounds: qacE
quinolones: OqxA, OqxB, qnrS1
trimethoprim: dfrA12
sulfonamides: sul1
macrolides: mrx(A), mphR(A), mph(A)
HMRGsArsenic: arsA, arsB, arsD, arsH, arsR,
Copper: pcoA, pcoB, pcoC, pcoD, pcoE. pcoR, pcoS
Copper/Silver: silA, silB, silC, silE, silP, silR, silS,
Tellurium: terB, terC, terD, terW, terZ		Arsenic arsA, arsB, arsD, arsH, arsR,
Copper: pcoA, pcoB, pcoD, pcoC, pcoE, pcoR, pcoS,
Copper/Silver: silA, silB, silC, silE, silP, silR, silS,
Tellurium: terC, terD, terW
Mercury: merA, merC, merD, merE, merP, merR, merT
TA-systemstype II: CcdA, CcdB, HicA, HicB, ParD, RatA, RelB, RelE/ParE, VapB/VapC, Yaf; type IV: YeeUtype II: HicA, HicB, ParD, RatA, RelB, RelE/ParE,
VapB/Vap; type IV: YeeU, CbtA
T1SSHlyD, TolC, MacA, MacBHlyD, TolC, MacA, MacB
T2SSgspC, gspD, gspE, gspF, gspG, gspH, gspI, gspJ,
gspK, gspL, gspM
secA, secB, secD, secE, secF, secG, secM, secY
tatA, tatB, tatC, tatD, ftsY, yajC, yidC		gspC, gspD, gspE, gspF, gspG, gspH, gspI, gspJ, gspK, gspL, gspM
secA, secB, secD, secE, secF, secG, secM, secY,
tatA, tatB, tatC, tatD, ftsY, yajC, yidC
T6SStssH, tssG, tssF, tssE, tagJ, tagH, tssC, tssB, tssA, tagF, tssM, tssK, tssJtssH, tssG, tssF, tssE, tagJ, tagH, tssC, tssB, tssA, tagF, tssM, tssK, tssJ
T4P1/T4SSLocus 1: hofC, hofB, HofQ, HofP, HofO, HofN/PilN, HofM/PilMLocus 1: hofC, hofB, HofQ, HofP, HofO, HofN/PilN, HofM/PilM
Locus 2: pilN, pilO2, pilP, pilV, pilM
Locus 3: pilM, pilN, pilO2, pilP
FlagellafliZACDSTfliZACDST
FimbriaeP7E32_11190, P7E32_11195, P7E32_11200, P7E32_11205, P7E32_11210, P7E32_11215, P7E32_11220, P7E32_11225, P7E32_11230P7F73_07975, P7F73_07980, P7F73_07985, P7F73_07990, P7F73_07995, P7F73_08000, P7F73_08005, P7F73_08010, P7F73_08015
# Differences among the strains are denoted in bold. ARGs: antimicrobial resistance genes; HMRGs: heavy-metal resistance genes.
Table 2. Identification of origin of transfers in DNA sequences of bacterial mobile genetic elements in strains EC-ML559 and EC-ML621.
Table 2. Identification of origin of transfers in DNA sequences of bacterial mobile genetic elements in strains EC-ML559 and EC-ML621.
TypeContig GenBank Accession No.Region/Gene (s)
Strain EC-ML559
oriT regionJARUPS0100000332799–2897
RelaxaseJARUPS010000032traI
T4CPJARUPS010000032traD
T4SS—region 1JARUPS010000031traN, trbC, traU, traW, traC, traV, traB, traK, traL, traE, traA
T4SS—region 2JARUPS010000023traM, traK, traJ, traH, traF, traB, traH, traG
Strain EC-ML621
oriT regionJASKGQ01000002618926–18990
RelaxaseJASKGQ010000026traI
T4CPJASKGQ010000026traD
T4SS—region 1JASKGQ010000016tfc2, DotB,
traH, traI, traJ, traK, traL, traM, traN, traO, traP, traQ, traR,
traT, traU, traV, traW, traX, traY
T4SS—region 2JASKGQ010000022traH, traI, traJ, traK, traL, traM, traN, traO, traP, traQ, traR, traU,
traW, traX, traY, trbC, trbB, trbA, trbN
T4SS—region 3JASKGQ010000024traM, traY, traA, traL, traE, traK, traB, traV, traC, traW, traU,
trbC, traN, traF, traQ, traB, traH, traG, traD, traI, traX
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mavroidi, A.; Gartzonika, K.; Spanakis, N.; Froukala, E.; Kittas, C.; Vrioni, G.; Tsakris, A. Comprehensive Analysis of Virulence Determinants and Genomic Islands of blaNDM-1-Producing Enterobacter hormaechei Clinical Isolates from Greece. Antibiotics 2023, 12, 1549. https://doi.org/10.3390/antibiotics12101549

AMA Style

Mavroidi A, Gartzonika K, Spanakis N, Froukala E, Kittas C, Vrioni G, Tsakris A. Comprehensive Analysis of Virulence Determinants and Genomic Islands of blaNDM-1-Producing Enterobacter hormaechei Clinical Isolates from Greece. Antibiotics. 2023; 12(10):1549. https://doi.org/10.3390/antibiotics12101549

Chicago/Turabian Style

Mavroidi, Angeliki, Konstantina Gartzonika, Nick Spanakis, Elisavet Froukala, Christos Kittas, Georgia Vrioni, and Athanasios Tsakris. 2023. "Comprehensive Analysis of Virulence Determinants and Genomic Islands of blaNDM-1-Producing Enterobacter hormaechei Clinical Isolates from Greece" Antibiotics 12, no. 10: 1549. https://doi.org/10.3390/antibiotics12101549

APA Style

Mavroidi, A., Gartzonika, K., Spanakis, N., Froukala, E., Kittas, C., Vrioni, G., & Tsakris, A. (2023). Comprehensive Analysis of Virulence Determinants and Genomic Islands of blaNDM-1-Producing Enterobacter hormaechei Clinical Isolates from Greece. Antibiotics, 12(10), 1549. https://doi.org/10.3390/antibiotics12101549

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop