- Research
- Open access
- Published:
Providencia pseudovermicola sp. nov.: redefining Providencia vermicola and unveiling multidrug-resistant strains from diabetic foot ulcers in Egypt
BMC Microbiology volume 25, Article number: 238 (2025)
Abstract
Background
Providencia species are concerning due to their intrinsic resistance to colistin and tigecycline, complicating the treatment of multidrug-resistant (MDR) infections.
Methods
In the current study, two MDR isolates, DFU6 and DFU52T, were recovered from infected diabetic foot ulcers in Egypt in 2024. Following their initial identification as Providencia stuartii using VITEK® 2 and MALDI-TOF-MS, the isolates were subjected to whole-genome sequencing via DNBseq.
Results
While the 16S rRNA gene showed 100% similarity to that of Providencia vermicola, phylogenomic analysis against the type strains in the TYGS database, including P. vermicola DSM 17385T confirmed that these isolates represent a distinct species within the genus, further supported by overall genome-relatedness indices (ORGIs). This discrepancy prompted us to revise the taxonomy of all published genomes of P. vermicola strains (n = 59) which revealed misidentification of at least 56 strains that are unrelated to the type strain of this species. DFU6 and DFU52T carried novel sequence types (ST29 and ST41, submitted to PubMLST) and harbored multiple resistance genes. Both strains contained the qnrD1 gene on a small, non-mobilizable plasmid. DFU52T possessed a conjugative plasmid encoding blaCMY−6, blaNDM−1, rmtC, aac(6’)-Ib10, sul1, aph(3’)-Ia, and qacEΔ1. DFU6 carried an ISEcp1-associated blaCTX−M−14, along with aadA, dfrA1, lnuF in a class 2 integron, and armA, msrE, and mphE on a resistance plasmid. Both isolates also featured a pathogenicity island (PAI) integrated into the pheV gene with fimbriae-encoding genes.
Conclusion
Following our reassessment of the taxonomic classification of all P. vermicola strains with published genomes, we propose reclassifying certain strains, including DFU6 and DFU52T, into distinct species for which we propose the name Providencia pseudovermicola sp. nov. We recommend DFU52T (= CCASU-2024-72) as the type strain for the novel species. We also shed light on the public health threat of this novel species as a human pathogen that harbours carbapenem and aminoglycoside resistance genes on mobile genetic elements.
Introduction
Providencia spp. are non-lactose fermenting, Gram-negative rods belonging to the family Morganellaceae within the order Enterobacterales [1]. Members of the genus Providencia can thrive in diverse environments [2] and have been isolated from various organisms [3].
Providencia is a rapidly evolving genus, with at least 18 species identified to date as per the List of Prokaryotic names with Standing in Nomenclature (LPSN), available at: https://lpsn.dsmz.de/genus/providencia. Although they are primarily opportunistic pathogens that predominantly affect immunocompromised and hospitalized patients [4], Proviencia spp. pose a significant public health concern. This is due to their ability to cause severe infections in vulnerable patients, combined with their intrinsic resistance to several last-resort antibiotics, such as colistin and tigecycline [5]. Moreover, growing evidence points to the emergence of extended-spectrum ß-lactamase (ESBL)- and carbapenemase-producer strains, further limiting treatment options [6,7,8,9,10,11]. The genus Providencia has been associated with various types of infections, most commonly urinary tract infections, gastroenteritis, bloodstream infections, and wound infections [12]. Providencia infections have been associated with high mortality rates [13]. Among all known Providencia species, Providencia rettgeri, Providencia stuartii, and Providencia alcalifaciens are the most common clinical MDR species and are known to be the primary cause of several nosocomial outbreaks [7, 12, 14, 15].
Remarkably, six new Providencia species have been discovered in the past five years, driven by advancements in whole genome sequencing (WGS) [16,17,18,19,20]. However, these developments have also raised concerns about the taxonomy of some Providencia species and the validity of species assignments for certain type strains. Dong, et al. [20] have recently redefined the species Providencia thailandensis as P. stuartii based on their average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values. A few years earlier, Andolfo, et al. [21] published the genome sequence of the Providencia vermicola type strain DSM 17385T. Using multilocus sequence typing with different marker sets, they revealed that none of the previously published P. vermicola genomes accurately represent this taxonomic species. The authors hence called for a reevaluation of the taxonomy of some P. vermicola strains with published genomes.
As the first report in Egypt, we present the genomes of two multidrug-resistant (MDR) Providencia sp. isolates recovered from infected diabetic foot ulcers (DFUs) in diabetic patients in Egypt. Genome-based phylogenetic analysis revealed that these isolates cluster with P. vermicola genomes previously reported by Andolfo, et al. [21] to be unrelated to the type strain of this species. After reassessing the taxonomic classification of all P. vermicola strains with published genomes, we propose the reclassification of certain P. vermicola strains, including ours, into a distinct species. We suggest the name Providencia pseudovermicola sp. nov. to reflect the previous misclassification and present one of our strains, DFU52T (= CCASU-2024-72), as the type strain for this novel Providencia species.
Materials and methods
Collection and preliminary identification of bacterial isolates
DFU6 and DFU52T were obtained from the microbiology laboratory of the National Institute of Diabetes and Endocrinology (NIDE). DFU6 was collected in February 2024 from a 52-year-old female patient, and DFU52T was obtained in March 2024 from a 43-year-old female patient. The isolates were received by the microbiology laboratory for microbiological identification and antimicrobial susceptibility testing as part of the routine clinical care of the patients admitted to the institute with infected DFUs. After identification using conventional microbiological techniques, the isolates were identified by Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) (BioMérieux; Marcy l’Etoile, France) and the VITEK® 2 automated system (bioMérieux, Lyon, France), using Gram-negative identification cards (GN ID card).
Antimicrobial susceptibility testing
The susceptibility of the two strains to a large panel of antimicrobial agents was evaluated using the Kirby-Bauer disc diffusion method. The following antimicrobial discs (Oxoid, UK) were utilized: amoxicillin/clavulanic acid (20/10 μg), piperacillin/tazobactam (100/10 μg), cefoxitin (30 μg), ceftriaxone (30 μg), cefotaxime (30 μg), cefepime (30 μg), aztreonam (30 μg), meropenem (10 μg), gentamicin (10 μg), amikacin (30 μg), ciprofloxacin (5 μg), tetracycline (30 μg), tigecycline (15 μg), chloramphenicol (30 μg), nitrofurantoin (300 μg), fosfomycin (200 μg), and trimethoprim/sulfamethoxazole (1.25/23.75 μg). Additionally, the broth microdilution assay was employed to determine the minimum inhibitory concentrations (MICs) of colistin (Sigma-Aldrich, St Louis, MO, USA) across a concentration range of 128–0.125 μg/ml. All susceptibility tests were conducted and interpreted according to the Clinical and Laboratory Standards Institute (CLSI) guidelines [22] for all antimicrobial agents, except for tigecycline, for which the susceptibility breakpoints recommended by EUCAST v14.0 for Enterobacterales were applied [23]. Escherichia coli ATCC 25922 and Pseudomonas aeruginosa ATCC 27853 served as quality control strains. The susceptibility profiles of the isolates were confirmed using the VITEK® 2 automated system for antimicrobial susceptibility testing, employing AST-GN73 cards following the manufacturer’s instructions.
Whole genome sequencing-based identification
Whole genome sequencing, assembly, and annotation
DNA was extracted following the manufacturer’s instructions using QIAamp DNA Kits (Qiagen, Hilden, Germany). WGS and library construction were carried out by BGI Tech Solutions Hong Kong Co., Ltd., China, using DNBseq™ sequencing technology. Quality control and trimming of the raw reads were completed in the pre-assembly stage using BGI’s SOAPnuke software [24]. The filtered reads were de novo assembled using the genome assembly tool provided by the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) [25]. The assemblers Unicycler version 0.4.8 [26] and PlasmidSPAdes were used for the genome and plasmids assembly, respectively. The resulting draft genome was submitted to the NCBI Prokaryotic Genome Annotation Pipeline for gene annotation. Gene sequences were assigned to Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology (KO) groups using GhostKOALA, which maps genes to KEGG pathways to identify their biological roles [27].
16S rRNA-based identification
The full-length 16S rRNA gene sequence (1,454 bp) was obtained from the WGS of strain DFU52T, while only a partial sequence (913 bp) was available for strain DFU6. Both sequences were submitted to the National Center for Biotechnology Information (NCBI) database (www.ncbi.nlm.nih.gov), with accession numbers PQ592516.1 and PQ592533.1, respectively. Pairwise alignment of the two sequences was performed using the Nucleotide Basic Local Alignment Search Tool (BLASTn) [28]. The gene sequences were also analyzed using the 16S rRNA-based identification service on the EzBioCloud server [10] to determine closely related type strains. Further alignment was carried out with entries from the NCBI using the Nucleotide Basic Local Alignment Search Tool (BLASTn). This search targeted the NCBI 16S ribosomal RNA Database for Bacteria and Archaea type strains, a curated resource that includes representative 16S rRNA sequences from bacterial and archaeal type strains [11]. To further ensure comprehensive results, a BLASTn search was also conducted against the NCBI nucleotide collection (nr/nt) database to find the closest species whose genes were not present in the 16S rRNA database.
Genome-based taxonomy
To begin the genome-based taxonomy analysis, the Type Strain Genome Server (TYGS) provided by DSMZ (https://tygs.dsmz.de/) was used (accessed on November 10, 2024) to identify the closest related type strain genomes. Using the MASH algorithm [29], a rapid estimation of intergenomic similarity was obtained by comparing the draft genome sequence with all type strains available in the TYGS database. The most similar type strains were selected to construct a phylogenomic tree. For the phylogenomic analysis, pairwise comparisons among the best-matching genomes were performed using the Genome BLAST Distance Phylogeny (GBDP) method [30], with precise intergenomic distances calculated using the ‘trimming’ algorithm and the d5 distance formula [30]. A balanced minimum evolution tree was then built using FASTME 2.1.6.1 with SPR post-processing [31]. Branch support was evaluated through 100 pseudo-bootstrap replicates. The resulting tree was rooted at the midpoint and visualized with PhyD3 [32]. Species and subspecies clustering is done using dDDH cut-off values of 70% and 79%, respectively. In the meantime, the genome-based taxonomy of the strains was investigated by tetra correlation search (TCS) against JSpeciesWS continuously updated reference database, accessed on November 9, 2024 [33]. After confirming that our strains belong to the genus Providencia, their taxonomic placement relative to the type strains of all known Providencia species was verified by calculating the Overall Genome Relatedness Indices (OGRIs), as previously recommended [34]. The calculated OGRIs included the Average Nucleotide Identity (ANI) using the alignment algorithms MUMmer (ANIm) and BLAST+ (ANIb), along with the correlation index of the Tetra-nucleotide signatures (Tetra). All calculations were performed using JSpeciesWS tools [33]. Additionally, digital DNA–DNA hybridization (dDDH) values were determined using the Genome-to-Genome Distance Calculator (GGDC) 3.0, accessible at https://ggdc.dsmz.de/ggdc.php [35].
Taxonomic reassessment of the species P. vermicola
As of our analysis date (November 10, 2024), only 59 P. vermicola genomes were identified in the NCBI genome database, all of which were downloaded for a genome-based taxonomic reassessment (Supplementary Table 1). Our analysis began with constructing a preliminary phylogenetic tree for the P. vermicola genomes, including our strains DFU6 and DFU52T, using the reference sequence alignment-based phylogeny builder (REALPHY) [36]. This free online tool enables phylogenetic inference from WGS data. All P. vermicola genomes were mapped to the reference genome P. vermicola DSM17385T (GenBank accession: JAGSPI010000019.1) via Bowtie2, with multiple sequence alignments generated from these mappings. From these alignments, a SNP-based phylogenetic tree was inferred using PhyML, employing Klebsiella pneumoniae subsp. pneumoniae MGH 78,578 (GenBank accession: NC_009648.1) as an outgroup. The tree was visualized using the interactive tree of life (iTOL) online tool v6.7 (https://itol.embl.de/itol.cgi). Genomes with branch lengths greater than zero (33 genomes) were then uploaded to TYGS (which accepts a maximum of 50 user genomes), along with our strains DFU6 and DFU52T, for species and subspecies clustering.
Multilocus sequence typing (MLST)
The multilocus sequence types (MLST) of the isolates were determined using the PubMLST database for Providencia spp. This database employs the MLST scheme outlined by Juneja and Lazzaro [37], which is based on five housekeeping genes: fusA, gyrB, ileS, lepA, and leuS.
Resistome, virulome, and mobilome analysis
The web portal Resistance Gene Identifier (RGI) version 6.0.3, hosted by the Comprehensive Antibiotic Resistance Database (CARD) [38], was utilized to identify genes that exhibit at least 80% coverage and 95% identity with known resistance genes listed in the CARD database version 3.3.0. Virulence genes were identified using VFanalyzer, a pipeline for screening known and potential VFs in bacterial genomes. The tool is provided by the Virulence Factors of Pathogenic Bacteria (VFDB) database [39]. The annotated contigs were also examined using SnapGene software version 5.1.3.1 from Insightful Science (http://www.snapgene.com) to explore the genetic environment surrounding various resistance and virulence genes.
Mobile genetic elements were identified using MobileElementFinder tool v1.0.3 [40] hosted by the Center for Genomic Epidemiology. VRprofile web server [23] was also used for the same purpose together with predicting the mobility of the identified plasmid sequences.
The draft genomes were mapped against the genome of P. vermicola strain P13 (GenBank accession: CP097327.1) using IslandViewer4 webtool (http://www.pathogenomics.sfu.ca/islandviewer/) [41] to identify the genomic islands. The integrative and conjugative elements (ICEs) were identified using ICEfinder, a web-based tool hosted by the ICEberg database that includes comprehensive information about bacterial ICEs [42]. Prophage regions were identified using PHASTEST (PHAge Search Tool with Enhanced Sequence Translation) webtool [43].
Phenotypic analysis of DFU52T, the type strain of P. pseudovermicola sp. nov.
Gram staining was performed to confirm the Gram reaction and examine the morphology of the isolate. The isolate was cultured under aerobic and anaerobic conditions on various media, including MacConkey agar, Luria Bertani (LB) agar, Trypticase Soy agar (TSA), and Müller-Hinton agar (MHA). The growth was evaluated in Trypticase Soy broth across a range of temperatures (4 °C to 45 °C), pH levels (4–10), and salinity concentrations (0–10%). Sporulation ability was assessed using thermal shock at 80 °C for 30 min. Enzymatic activities, such as catalase and oxidase, were tested following previously described methods [44]. The complete biochemical profile of DFU52T was determined using the VITEK® 2 system.
Results
Preliminary identification and antimicrobial susceptibility profiles
Both DFU6 and DFU52T were identified as P. stuartii using the VITEK® 2 and MALDI-TOF-MS identification system.
Both DFU6 and DFU52T were classified as MDR, exhibiting non-susceptibility to at least one antimicrobial agent in three or more antimicrobial classes, following the MDR phenotype definition provided before [45]. The complete antimicrobial susceptibility profiles of the isolates are presented in Table 1.
Genome assembly and annotation metrics
The assembly features of the genomes of DFU6 and DFU52T revealed a genome size of 4,416,300 bp and 4,423,558 bp, respectively, with G + C content values of 40.78% for DFU6 and 41.16% for DFU52T. DFU6 draft genome contained 50 contigs with an N50 of 735,893 bp and an L50 of 2, while DFU52T had 72 contigs with an N50 of 172,528 bp and an L50 of 9. Regarding genome quality, the completeness was 99.9% for DFU6 and 100% for DFU52T, with contamination levels of 0.4% and 0.7%, respectively. Annotation features showed a total of 4,024 genes in DFU6 and 4,084 genes in DFU52T, including 3,881 and 3,937 protein-coding CDSs and 79 and 78 RNA genes, respectively. DFU6 and DFU52T both carried type I-F CRISPR-associated protein-coding genes, with 4 and 3 CRISPR arrays, respectively. The functional annotations of the draft genomes are shown in Supplementary Table 2.
16S rRNA-based identification
The pairwise alignment of the partial sequence of the 16S rRNA gene of DFU6 showed 100% identity to that of DFU52T. Upon uploading the full 16S rRNA sequence of DFU52T to the Ezbiocloud 16S-based ID tool, the highest similarity score (99.66%) was obtained with P. manganoxydans LLDRA6. Alignment of the full sequence against the NCBI rRNA type strains/16S ribosomal RNA database revealed the highest match with P. huaxiensis strain WCHPr000369 (NR_174258.1) with 98% coverage and 99.08% identity, followed closely by P. vermicola strain OP1 (NR_042415.1) with 98% coverage and 99.21% identity.
nBLAST analysis of the full 16S rRNA sequence against the nucleotide (nr/nt) database showed 100% identity and coverage with seven Providencia strains, four of which were identified to species level as P. vermicola. Other strains identified as P. vermicola showed slightly lower identity percentages (99.29–99.87%). This discrepancy in 16S rRNA gene-based identification indicated a need for further analysis using whole-genome data.
Genome-based phylogeny and taxonomy
The phylogenomic tree presented in Fig. 1, which includes DFU6 and DFU52T alongside type strains from the TYGS database, demonstrates that DFU6 and DFU52T form a distinct cluster at the species and subspecies levels. Notably, they do not cluster with any type strains from other species in the database, including the P. vermicola type strain DMS17385T.
Phylogenomic tree constructed using GBDP distances derived from the genomes of TYGS type strains most closely related to DFU6 and DFU52T. Clusters representing species and subspecies are defined by dDDH thresholds of 70% and 79%, respectively. Strains from the current study are highlighted in bold red within the tree
All genomes found in the JSpeciesWS genome database showed a Z-Score of less than 0.999 set by the server as the cut-off value for species assignment. The highest Z-score was found for P. stuartii 50655837 with a value of 0.98781. The genomes of the type strains of all currently known Providencia species were included in all-versus all OGRIs analysis as shown in Fig. 2. All OGRIs of DFU6 and DFU52T to all type strains including P. vermicola type strain DMS17385T, were below the species delineation thresholds.
Taxonomic reassessment of the species P. vermicola
The significant discrepancies between the 16S rRNA analysis and genome-based taxonomy, regarding the similarity to P. vermicola, prompted additional analyses to confirm the taxonomic positioning of P. vermicola genomes available in the NCBI database, whose 16S rRNA genes shared 100% identity with those of our strains.
The preliminary SNP-based phylogenetic tree, illustrated in Supplementary Fig. 1, was constructed using our strains in conjunction with all available P. vermicola genomes (n = 59). This analysis facilitated the selection of 33 P. vermicola genomes with branch lengths greater than zero for TYGS species and subspecies clustering, alongside type strains deposited in the TYGS database. As illustrated in Fig. 3, the P. vermicola genomes were categorized into three distinct species clusters. Notably, the majority of P. vermicola genomes clustered with our strains, which are distinctly separated from the cluster containing the P. vermicola type strain DMS 17385T. Our analysis confirms that at least 56 strains were misidentified as P. vermicola and that their dDDH values compared to P. vermicola DSM 17385T are far below the species delineation cut-off values. Based on our analysis, only fly-401 and fly-1209 correctly belong to the species P. vermicola. Additionally, the strains fly-1010, fly-1073, and strain S1-B1-56 belong to a novel species that has yet to be identified.
Multilocus sequence types of DFU6 and DFU52T
DFU6 carried novel alleles of gyrB (allele 61), ileS (allele 73), lepA (allele 60), and leuS (allele 61), resulting in the unique allelic profile [fusA9 - gyrB61 - ileS73 - lepA60 - leuS61], which was submitted to the MLST database and designated as ST29. Similarly, DFU52T carried the novel alleles ileS (allele 74) and leuS (allele 60), forming the distinct allelic profile [fusA9 - gyrB31 - ileS74 - lepA9 - leuS60], which was submitted to the MLST database and designated as ST41.
Phylogenomic tree constructed using GBDP distances derived from DFU6, DFU52T, and other P. vermicola genomes retrieved from the NCBI database alongside most closely related TYGS type strains. Clusters representing species and subspecies are defined by dDDH thresholds of 70% and 79%, respectively. Strains from the current study are highlighted in bold red within the tree
Resistance genes and associated MGEs
The genomes of DFU6 and DFU52T contained a diverse array of resistance genes, with some predicted to be located on the chromosomes and others associated with MGEs. Genes encoding efflux pumps from all six currently recognized efflux pump classes were identified in both genomes. Additional genes with the potential to confer resistance to various classes of antimicrobial agents are detailed in Table 2.
DFU52T harbored a 140,281 bp IncA/C2 plasmid, designated pDFU52T_1, as depicted in Fig. 4. This plasmid was predicted to be conjugative and carried seven resistance genes conferring resistance to β-lactams, including carbapenems, as well as aminoglycosides, sulfonamides, and disinfectants. Many of the resistance genes were associated with MGEs, enabling potential independent mobility. The aminoglycoside resistance gene aph(3’)-Ia was flanked by IS26 within the composite transposon Tn4352. Additionally, the genes aac(6’)-Ib10, qacEdelta1, and sul1 were located within a class I integron, while blaCMY−6 was linked to the insertion sequence ISEcp1.
Both DFU6 and DFU52T harboured the qnrD1 gene on a small, 2,122 bp non-mobilizable plasmid (Fig. 4). This plasmid exhibited high similarity to other small plasmids found in bacterial strains from the genera Providencia, Salmonella, Proteus, Morganella, and Moellerella.
Circular maps illustrating the resistance plasmids carried by DFU52T. The unlabelled arrows indicate open reading frames that encode hypothetical proteins. This figure was generated using the Proksee online tool [46]
In DFU6, the resistance genes aadA, dfrA1, and lnuF were associated with a class 2 integron and were carried on a 84,533 contig that showed 99.99% similarity to P. rettgeri YPR31 plasmid pYPR31 (GenBank accession: CP053897.1) and Proteus mirabilis HBRJC7FQ plasmid pHBRJC7 (GenBank accession: MK630213.1). Moreover, the genes armA, msrE, and mphE were identified on 10,530 bp contig that showed 99.99% similarity to resistance plasmids carried by Enterobacteriaceae species. The full sequence of the resistance plasmids carried by DFU6 could not be identified. A genomic island (GI) was identified at the 3’ end of the mnmE gene, which encodes a tRNA modification GTPase, in DFU6. A partial sequence of the GI, spanning 21,926 base pairs, was identified (Supplementary Fig. 2). This sequence exhibited the highest similarity to Proteus mirabilis genomic island 2 (PGI2) (GenBank accession: MK847916.1), with 99% coverage and 99.46% identity. Only one resistance gene (dfrA15) was associated with the class I integron carried by the genomic island.
Interestingly, DFU6 carried an ISEcp1-mediated insertion of blaCTX−M−14 with a 5 bp target site duplication (TCCTT), as shown in Fig. 5. The genes were inserted into a fimbrial protein-coding gene with no similarity in the NCBI nucleotide database.
Gene maps depicting the ISEcp1-mediated blaCTX−M−14 insertion in DFU6 (B), compared to P. stuartii strain CMC-4104 chromosome (GenBank Accession: CP095443.1 - region: 1234044..1241727) (A) and P. vermicola P8538 chromosome (GenBank accession: CP048796.1– region: 1167541..1176163) (C). Arrows represent open reading frames and are labelled by their predicted protein products. Grey panels between maps correspond to the similarity percent between the chromosomal regions. The figure was created by Easyfig version 2.2.5
While most resistance genes were predicted to be associated with MGEs, the genes aac(2’)-Ia, ampC, and catA were identified on the chromosomes of DFU6 and DFU52T and showed no association with MGEs.
Virulence genes and associated MGEs
The functional annotation of the draft genomes revealed that both DFU6 and DFU52T carried at least 74 genes associated with secretion systems and 50 genes encoding bacterial motility proteins. DFU6 encoded 371 transporter-related genes, slightly more than the 364 identified in DFU52T. Additionally, DFU6 harbored 29 genes linked to prokaryotic defense systems, compared to 30 in DFU52T, as detailed in Supplementary Table 2.
BLAST analysis of the draft genomes against the VFDB was performed to identify genes associated with experimentally confirmed virulence factors. This allowed the identification of virulence factors involved in adherence, motility, secretion systems, metal uptake, immune evasion, and stress survival. The genomes of DFU6 and DFU52T encoded genes for at least five types of fimbriae and pili, including type I fimbriae, type 3 fimbriae, type IV pili, E. coli common pilus (ECP), and P fimbriae. Additionally, the genomes harbored genes for five secretion systems: type I, type III and type VI secretion systems, as well as the Sec and Tat secretion pathways. All virulence genes identified in the genomes of DFU6 and DFU52T are listed in Supplementary Table 3.
Additional virulence genes identified from the annotations include those encoding the HlyD family type I secretion (T1SS) periplasmic adaptor subunit, the type I secretion system permease/ATPase, and the type IV secretion system. Both genomes also contained a hemolysin family protein with 59% amino acid sequence identity to the HlyA hemolysin protein, as well as a ShlB/FhaC/HecB family hemolysin secretion/activation protein. Additionally, the quorum sensing QseBC two-component system was encoded by both genomes.
An ICE was identified in the draft genomes of DFU6 and DFU52T, showing 99.99% identity to PAIs found in P. stuartii and other members of the order Enterobacteriales, including Escherichia coli, Klebsiella pneumoniae, Morganella morganii, and Proteus mirabilis. The island was integrated into the phenylalanine tRNA gene (pheV), with total sizes of 92,644 bp in DFU6 and 92,668 bp in DFU52T, sharing a sequence identity of 99.92%. As depicted in Fig. 6, the islands contained fimbriae-coding loci with predicted protein products showing 58% identity to F17 fimbriae produced by Escherichia coli 111KH86. They also harbored genes encoding a TonB-dependent receptor, the type IV pilus biosynthesis protein PilL, and the trimeric autotransporter adhesin of Proteus (TaaP). Additionally, the islands carried biosynthetic gene clusters for type I polyketide synthase (T1PKS) and non-ribosomal peptide synthase (NRPS) production. Genes coding ArdC family antirestriction proteins were also found within the island.
Genetic map illustrating the organization of the PAI carried by DFU6 and DFU52T. Red arrows indicate genes involved in the conjugation and integration of the island, black arrows denote virulence genes, orange arrows represent insertion sequences, blue arrows highlight genes encoding ArdC family antirestriction proteins, purple arrows correspond to core biosynthetic genes, green arrows signify additional biosynthetic genes, light blue arrows indicate transport-related genes, and white arrows represent all other genes
Prophage regions
PHASTEST analysis of the DFU6 draft genome identified two intact prophage regions, with genes most commonly matching Entero_phiT5282H (NC_049429) and Cronob_phiES15 (NC_018454) (Fig. 7). In contrast, the DFU52T draft genome contained three questionable prophage regions, with genes most commonly matching Burkho_BcepB1A (NC_005886), Cronob_phiES15 (NC_018454), and Escher_HK75 (NC_016160). A detailed description of the prophage regions identified in DFU6 and DFU52T is provided in Supplementary Table 4.
Phenotypic features of DFU52T, the type strain of P. pseudovermicola sp. nov
DFU52T are Gram-negative, non-spore-forming rods that exhibited robust growth under both aerobic and anaerobic conditions on all tested culture media. The isolates demonstrated the ability to grow across the full range of tested temperatures, pH levels, and salinity concentrations. DFU52T tested negative for catalase and oxidase activities. Biochemical features of DFU52T compared to the type strains of other Providencia species are shown in Table 3. The full biochemical profile of DFU52T as defined by VITEK® 2 system is shown in Supplementary Table 5.
Discussion
DFUs are among the most prevalent and severe complications associated with diabetes mellitus, often resulting from peripheral nerve damage and arterial disease [47]. These ulcers frequently become infected, with over half showing signs of infection at the time of presentation, leading to significant morbidity and mortality [48]. The treatment of infected DFUs is challenging due to several factors, including patients’ compromised immune status and poor peripheral circulation. In addition, the polymicrobial nature of the wounds, biofilm formation by infecting pathogens, and the widespread antimicrobial resistance among causative organisms further complicates DFU management [49]. Recent reports from around the world [50,51,52,53] including Egypt [54, 55], have documented a rising prevalence of Providencia species isolated from DFUs. These infections exacerbate the challenges in managing DFUs due to their intrinsic resistance to several classes of antimicrobial agents [56]. Several key mechanisms mediate intrinsic resistance in Providencia species: (a) the inducible expression of AmpC β-lactamases and AAC(2’)-Ia provides resistance to most penicillins, first-generation cephalosporins, and aminoglycosides [57, 58]; (b) tetracycline resistance is driven by the constitutive activity of multidrug efflux pumps [59]; (c) Polymyxin resistance is associated with cell envelope modifications that either inhibit colistin binding to lipid targets or involve lipid A alterations that reduce its binding affinity [60].
In the current study, two MDR Providencia sp. isolates were recovered from infected DFUs in diabetic patients from Egypt. A notable discrepancy was observed between phenotypic identification methods and 16S rRNA-based identification. This prompted us to conduct a more in-depth genome-based taxonomic analysis. The 16S rRNA gene analysis revealed 100% identity to P. vermicola strains; however, in the GBDP tree generated by TYGS, our strains did not cluster with the type strain of this species. As previously recommended [34, 61], the OGRIs for our strains were calculated in comparison to the type strains of all known Providencia species, and all values were well below the established thresholds for species delineation: 95% for ANI, 70% for dDDH, and 0.999 for Tetra [62]. While the reliability of 16S rRNA-based identification has been questioned in several studies [20, 63, 64], the classification error of our isolates was attributed to the misclassification of most P. vermicola strains with genomes deposited in the NCBI database, leading to confusion in the taxonomy of this species.
In 2021, Andolfo, et al. [21] published the complete genome of the P. vermicola type strain DSM 17385T. MLST analysis in the same study using different marker sets revealed that none of the P. vermicola genomes published at the time, including the strains P8538, LLDRA6, and G1, truly belonged to this taxonomic species. By the time of our study, the number of genomes in the NCBI database assigned to P. vermicola had grown to 59, with 57 genomes submitted after the publication by Andolfo, et al. [21]. Of these, only two strains were correctly classified as P. vermicola, while three strains appear to represent a novel species that requires further characterization. The remaining strains, whose genomes clustered with ours, belong to a distinct species, for which we propose the name P. pseudovermicola. We attribute this misclassification to the earlier publication of the genome for Providencia strain P8538 [65] before that of the P. vermicola type strain DSM 17385T [21]. To address this taxonomic confusion, we designate DFU52T (= CCASU-2024-72) as the type strain of the novel species P. pseudovermicola, accompanied by its full biochemical characterization and genome analysis.
DFU6 and DFU52T were found to carry acquired resistance genes for multiple classes of antimicrobial agents most of them were associated with mobile genetic elements such as plasmids, integrons, genomic islands, and insertion sequences. Given that Providencia species are inherently resistant to many antimicrobial classes, including last-line agents such as colistin and tigecycline [66], the acquisition of additional resistance genes may contribute to the emergence of extensive drug resistance or pandrug resistance. DFU52T harbored a conjugative plasmid containing six resistance genes, including blaCMY6, blaNDM−1 and rmtC, which confer resistance to cephalosporins, carbapenems, and all members of the aminoglycoside family. Resistance plasmids of Providencia species carrying carbapenemase-coding genes have been widely reported [9, 11, 67, 68], underscoring the high capacity of this genus to acquire and disseminate resistance genes. Interestingly, DFU6 carried a unique ISEcp1-mediated insertion of blaCTX−M−14. ISEcp1 was reported by Poirel, et al. [69] to mobilize the adjacent blaCTX−M genes by a transpositional mechanism. This was reported in many Gram-negative species [70, 71] and was confirmed here in DFU6 through the identification of the target site duplication, a hallmark of transposition [72]. Both isolates carried qnrD gene on a small non-mobilizable plasmid commonly found in the members of the family Morganellaceae, including strains from several Providencia species [21, 73,74,75].
The phenotypic antimicrobial susceptibility profiles of DFU52T and DFU6 aligned with their genotypic characteristics. DFU52T demonstrated susceptibility exclusively to cefepime, aztreonam, and trimethoprim/sulfamethoxazole. Cefepime is a widely recommended treatment for Enterobacteriales infections [76]. The susceptibility of the NDM-producing DFU52T to aztreonam can be attributed to the ability of aztreonam to evade Metallobetalactamase-mediated hydrolysis [77]. Its sensitivity to trimethoprim/sulfamethoxazole supports recent efforts to reintroduce older antimicrobial agents into clinical use [78, 79]. Although DFU52T carried the sulfonamide resistance gene sul1, no trimethoprim resistance genes were detected. Conversely, DFU6, which harbored dfrA1 and dfrA15 trimethoprim resistance genes, was resistant to trimethoprim/sulfamethoxazole. Additionally, the genome of DFU6 did not carry any carbapenemases-coding genes, making it sensitive to meropenem.
While the virulence factors of most Providencia species remain largely unexplored, key contributors to the pathogenicity of Providencia species that have been studied include cellular adherence, fimbriae production, secretion systems, and urease activity [3, 80]. Most research in this area has focused on Providencia stuartii. While most virulence factors have been studied in the context of UTIs [3], the most common Providencia infection, some are also relevant to wound infections and DFUs. Genes coding for Type 3 (MR/K) fimbriae, Type I fimbriae, Type 4 pilli, E. coli common pilus (ECP), and P fimbriae were carried by the genomes of DFU6 and DFU52T. These structures have been previously reported to play crucial roles in host cell adhesion and biofilm formation [81,82,83]. Biofilm formation is a critical step in the pathophysiology of DFUs. It significantly contributes to disease progression, lesion chronicity, antibiotic resistance emergence, and wound healing challenges [84]. Motility also plays a crucial role in skin infections, enabling movement across the skin surface and facilitating relocation into wounds [85]. At least five motility-related loci were found in the genomes of DFU6 and DFU52T. Genes encoding six types of secretion systems were identified, including T1SS, T3SS, T6SS, the Sec secretion pathway, and the Tat secretion pathway. Secretion systems transport substrates including adhesins, enzymes, siderophores, and toxins across bacterial membranes and play an important role in colonizing various eukaryotic hosts [80, 86, 87]. Although the role of the urease enzyme in UTIs caused by Providencia species is well-documented [3], its contribution to wound infections remains poorly understood. Urease catalyzes the conversion of urea into ammonia, leading to a significant increase in the pH of the surrounding tissue. This can make wounds more susceptible to opportunistic infections, promoting polymicrobial infections, a common feature of DFUs [88]. Urease also contributes to pathologic effects by causing cytotoxicity to host cells, stimulating inflammatory responses, and enhancing biofilm colonization and persistence [89, 90]. The TCS QseBC quorum-sensing regulatory system, encoded by the genomes of DFU6 and DFU52T, is a global regulator of biofilm growth, bacterial motility, and virulence in E. coli and K. pneumoniae [91].
DFU6 and DFU52T carried a PAI with 99.95% sequence identity to the ICE known as ICEPm1. This PAI is commonly associated with P. mirabilis, P. stuartii, and M. morganii [92]. The self-transmissibility of ICEPm1 has been previously confirmed [93], and the island has shown sequence matches in other species, including E. coli and K. pneumoniae. It carries genes encoding iron acquisition siderophores and Taap a surface adhesin and autoagglutinin that contribute to the virulence of pathogenic Gram-negative bacteria by facilitating adherence, biofilm formation, invasion, survival within eukaryotic cells, serum resistance, and cytotoxicity [94, 95].
Despite possessing CRISPR systems, both DFU6 and DFU52T harbored prophage genomes. Together with resistance and virulence-associated MGEs, prophages play a crucial role in enhancing the evolutionary potential of the genus Providencia.
As we propose DFU52T as the type strain of the novel species P. pseudovermicola, we provide a comprehensive biochemical profile to facilitate its identification and differentiation from other Providencia species. The biochemical analysis reveals that all Providencia species, including DFU52T, are negative for H2S production and predominantly lack the enzymes ß-galactosidase, ornithine decarboxylase, and lysine decarboxylase, suggesting that these traits may be characteristic of the genus. Glucose fermentation is a consistent metabolic feature across all species, underscoring its essential role in the core metabolic pathways of the genus. However, the ability to ferment other sugars, such as D-mannitol, sucrose, D-xylose, and cellobiose, varies among species. Additionally, variability is observed in traits such as indole production, citrate utilization, and urease activity, which may correspond to species-specific ecological niches or adaptive strategies. These differences provide valuable markers for distinguishing P. pseudovermicola DFU52T from other members of the genus.
Conclusion
This study reports the identification of a novel Providencia species, P. pseudovermicola sp. nov., isolated from infected DFUs in Egypt. Phylogenomic analyses confirmed P. pseudovermicola as a distinct species within the genus Providencia, resolving misidentifications of numerous strains previously classified as P. vermicola. We propose DFU52T (= CCASU-2024-72) as the type strain for this novel species. The detection of MGEs carrying carbapenem and aminoglycoside resistance genes, coupled with intrinsic resistance to colistin and tigecycline, highlights the significant public health threat posed by the novel species. These findings emphasize the need for enhanced identification and monitoring of Providencia species in clinical settings. Moreover, the application of genomic tools to refine species classification and trace the dissemination of antimicrobial resistance genes is essential. Future research should focus on elucidating the pathogenic mechanisms and epidemiological distribution of P. pseudovermicola, along with its clinical impact on human infections.
Description of Providencia pseudovermicola sp. nov.
pseudovermicola. N.L. fem. adj. pseudovermicola, referring to the incorrect classification as Providencia vermicola, with “pseudo-” meaning false or incorrect and “vermicola” referring to the original misidentification.
The novel species P. pseudovermicola sp. nov. exhibited characteristics typical of the genus Providencia, including fermentative metabolism, facultative anaerobic growth, and Gram-negative, rod-shaped cell morphology. It lacked spore formation and tested positive for catalase activity and negative oxidase activity. The strain demonstrated the ability to grow across a wide range of temperatures (4–45 °C), salinity levels (0–10% NaCl), and pH values (4–10). It can grow well on different culture media including MacConkey agar, LB agar, TSA, and MHA.
On TSA, it formed small, round colonies with a smooth, creamy appearance. Acid production was observed from D-glucose and D-mannitol, but not from D-sorbitol, D-cellobiose, D-xylose, or sucrose. The strain tested positive for citrate utilization and urease activity, but negative for ß-galactosidase, ornithine decarboxylase, lysine decarboxylase, indole production, and H2S production. These biochemical traits distinguish P. pseudovermicola from closely related species within the genus.
The type strain, DFU52T (= CCASU-2024-72), was isolated from an infected DFU in Egypt and has been deposited in the Culture Collection Ain Shams University.
Data availability
The complete sequence of the 16S rRNA gene from DFU52T has been submitted to the NCBI GenBank database with the accession number PQ592516. In addition, a partial sequence of the same gene associated with DFU6 has been submitted under the accession number PQ592533. Both draft genomes are accessible in the BioProject PRJNA1181229, with the respective accession numbers JBIZWE000000000 for DFU6 and JBIZWF000000000 for DFU52T.
Abbreviations
- ANI:
-
Average nucleotide identity COGs: Clusters of orthologous group
- ANIb:
-
Average nucleotide identity calculated by the alignment algorithms BLAST +
- ANIm:
-
Average nucleotide identity calculated by MUMmer
- dDDH:
-
The digital DNA–DNA hybridization
- DFU:
-
Diabetic foot ulcer
- GBDP:
-
Genome BLAST distance phylogeny
- GI:
-
Genomic island
- ICE:
-
Integrative and conjugative element
- MDR:
-
Multidrug-resistant
- MLST:
-
Multilocus sequence typing
- OGRIs:
-
Overall genome-relatedness indices
- PAI:
-
Pathogenicity island
- Tetra:
-
Tetra-nucleotide signatures
- TYGS:
-
Type strain genome server
- WGS:
-
Whole genome sequencing
References
Adeolu M, Alnajar S, Naushad S, R SG. Genome-based phylogeny and taxonomy of the ‘enterobacteriales’: proposal for enterobacterales Ord. Nov. Divided into the Fam.lies enterobacteriaceae, Erwiniaceae Fam. Nov., pectobacteriaceae Fam. Nov., Yersiniaceae Fam. Nov., hafniaceae Fam. Nov., Morganellaceae Fam. Nov., and Budviciaceae Fam. Nov. Int J Syst Evol Microbiol. 2016;66(12):5575–99. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijsem.0.001485
O’Hara CM, Brenner FW, Miller JM. Classification, identification, and clinical significance of proteus, Providencia, and Morganella. Clin Microbiol Rev. 2000;13(4):534–46. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/CMR.13.4.534
Manos J, Belas R. The Genera Proteus, Providencia, and Morganella. In: The Prokaryotes: a Handbook on the Biology of Bacteria Volume 6: Proteobacteria: Gamma Subclass. Edited by Dworkin M, Falkow S, Rosenberg E, Schleifer K-H, Stackebrandt E. New York, NY: Springer New York; 2006: 245–269.
Galac MR, Lazzaro BP. Comparative pathology of bacteria in the genus Providencia to a natural host, drosophila melanogaster. Microbes Infect. 2011;13(7):673–83. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.micinf.2011.02.005
Yuan C, Wei Y, Zhang S, Cheng J, Cheng X, Qian C, et al. Comparative genomic analysis reveals genetic mechanisms of the variety of pathogenicity, antibiotic resistance, and environmental adaptation of Providencia genus. Front Microbiol. 2020;11:572642. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2020.572642
Barrios H, Garza-Ramos U, Reyna-Flores F, Sanchez-Perez A, Rojas-Moreno T, Garza-Gonzalez E, et al. Isolation of carbapenem-resistant NDM-1-positive Providencia rettgeri in Mexico. J Antimicrob Chemother. 2013;68(8):1934–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jac/dkt124
Liu M, Yi N, Wang X, Wang R. Analysis of resistance genes of carbapenem-resistant Providencia rettgeri using whole genome sequencing. BMC Microbiol. 2023;23(1):283. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12866-023-03032-3
Abdallah M, Alhababi R, Alqudah N, Aldyyat B, Alharthy A. First report of carbapenem-resistant Providencia stuartii in Saudi Arabia. New Microbes New Infect. 2018;26:107–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.nmni.2018.09.007
Iwata S, Tada T, Hishinuma T, Tohya M, Oshiro S, Kuwahara-Arai K, et al. Emergence of carbapenem-resistant Providencia rettgeri and Providencia stuartii producing IMP-type metallo-beta-lactamase in Japan. Antimicrob Agents Chemother. 2020;64(11). https://doiorg.publicaciones.saludcastillayleon.es/10.1128/AAC.00382-20
Tshisevhe VS, Lekalakala MR, Tshuma N, van Janse S, Mbelle N. Outbreak of carbapenem-resistant Providencia rettgeri in a tertiary hospital. S Afr Med J. 2016;107(1):31–3. https://doiorg.publicaciones.saludcastillayleon.es/10.7196/SAMJ.2016.v107.i1.12002
Bocanegra-Ibarias P, Duran-Bedolla J, Silva-Sanchez J, Garza-Ramos U, Sanchez-Perez A, Garza-Gonzales E, et al. Identification of Providencia spp. Clinical isolates co-producing carbapenemases IMP-27, OXA-24, and OXA-58 in Mexico. Diagn Microbiol Infect Dis. 2024;109(1):116246. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.diagmicrobio.2024.116246
Malviya M, Kale-Pradhan P, Coyle M, Giuliano C, Johnson LB. Clinical and drug resistance characteristics of Providencia infections. Microorganisms. 2024;12(10). https://doiorg.publicaciones.saludcastillayleon.es/10.3390/microorganisms12102085
Edwards F, Edwards KH, Tabah A, Harris PNA, Laupland KB. Differentiation of Providencia species bloodstream infections: a population-based analysis. Diagn Microbiol Infect Dis. 2025;111(1):116602. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.diagmicrobio.2024.116602
Mnif B, Ktari S, Chaari A, Medhioub F, Rhimi F, Bouaziz M, et al. Nosocomial dissemination of Providencia stuartii isolates carrying Bla OXA-48, Bla PER-1, Bla CMY-4 and qnrA6 in a Tunisian hospital. J Antimicrob Chemother. 2013;68(2):329–32. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jac/dks386
Shah MM, Odoyo E, Ichinose Y. Epidemiology and pathogenesis of Providencia alcalifaciens infections. Am J Trop Med Hyg. 2019;101(2):290–3. https://doiorg.publicaciones.saludcastillayleon.es/10.4269/ajtmh.18-0376
Yang W, Chen J, Yang F, Ji P, Shen S, Yin D, et al. Identification of a novel Providencia species showing multi-drug-resistant in three patients with hospital-acquired infection. Int J Antimicrob Agents. 2024;64(2):107211. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ijantimicag.2024.107211
Dong X, Xiang Y, Yang P, Wang S, Yan W, Yuan Y, et al. Novel Providencia xianensis Sp. nov.: a multidrug-resistant Sp.cies identified in clinical infections. Eur J Clin Microbiol Infect Dis. 2024;43(7):1461–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10096-024-04821-y
Li Z, Liao F, Ding Z, Chen S, Li D. Providencia manganoxydans Sp. nov., a Mn(II)-oxidizing bacterium isolated from heavy metal contaminated soils in Hunan Province, China. Int J Syst Evol Microbiol. 2022;72(8). https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijsem.0.005474
Zhou K, Liang J, Dong X, Zhang P, Feng C, Shi W, et al. Identification and characterization of a novel chromosomal aminoglycoside 2’-N-Acetyltransferase, AAC(2’)-If, from an isolate of a novel Providencia species, Providencia wenzhouensis R33. Front Microbiol. 2021;12:711037. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2021.711037
Dong X, Jia H, Yu Y, Xiang Y, Zhang Y. Genomic revisitation and reclassification of the genus providencia. mSphere. 2024;9(3):e0073123. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/msphere.00731-23
Andolfo G, Schuster C, Gharsa HB, Ruocco M, Leclerque A. Genomic analysis of the nomenclatural type strain of the nematode-associated entomopathogenic bacterium Providencia vermicola. BMC Genomics. 2021;22(1):708. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12864-021-08027-w
CLSI. Performance standards for antimicrobial susceptibility testing; 30th edition. CLSI supplement M100. Wayne, PA: Clinical and Laboratory Standards Institute.; 2020.
Breakpoint tables for. interpretation of MICs and zone diameters [http://www.eucast.org]
Chen Y, Chen Y, Shi C, Huang Z, Zhang Y, Li S, et al. SOAPnuke: a mapreduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience. 2018;7(1):1–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/gigascience/gix120
Olson RD, Assaf R, Brettin T, Conrad N, Cucinell C, Davis JJ, et al. Introducing the bacterial and viral bioinformatics resource center (BV-BRC): a resource combining PATRIC, IRD and vipr. Nucleic Acids Res. 2023;51(D1):D678–89. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkac1003
Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13(6):e1005595. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pcbi.1005595
Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428(4):726–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jmb.2015.11.006
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S0022-2836(05)80360-2
Ondov BD, Treangen TJ, Melsted P, Mallonee AB, Bergman NH, Koren S, et al. Mash: fast genome and metagenome distance estimation using MinHash. Genome Biol. 2016;17(1):132. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13059-016-0997-x
Meier-Kolthoff JP, Auch AF, Klenk HP, Goker M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics. 2013;14:60. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/1471-2105-14-60
Lefort V, Desper R, Gascuel O. FastME 2.0: a comprehensive, accurate, and fast distance-based phylogeny inference program. Mol Biol Evol. 2015;32(10):2798–800. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/molbev/msv150
Kreft L, Botzki A, Coppens F, Vandepoele K, Van Bel M. PhyD3: a phylogenetic tree viewer with extended phyloxml support for functional genomics data visualization. Bioinformatics. 2017;33(18):2946–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/bioinformatics/btx324
Richter M, Rossello-Mora R, Oliver Glockner F, Peplies J. JSpeciesWS: a web server for prokaryotic species circumscription based on pairwise genome comparison. Bioinformatics. 2016;32(6):929–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/bioinformatics/btv681
Riesco R, Trujillo ME. Update on the proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. Int J Syst Evol Microbiol. 2024;74(3). https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijsem.0.006300
Meier-Kolthoff JP, Carbasse JS, Peinado-Olarte RL, Goker M. TYGS and LPSN: a database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic Acids Res. 2022;50(D1):D801–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkab902
Bertels F, Silander OK, Pachkov M, Rainey PB, van Nimwegen E. Automated reconstruction of whole-genome phylogenies from short-sequence reads. Mol Biol Evol. 2014;31(5):1077–88. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/molbev/msu088
Juneja P, Lazzaro BP. Providencia sneebia Sp. Nov. And Providencia Burhodogranariea Sp. Nov., isolated from wild drosophila melanogaster. Int J Syst Evol Microbiol. 2009;59(Pt 5):1108–11. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijs.0.000117-0
Alcock BP, Huynh W, Chalil R, Smith KW, Raphenya AR, Wlodarski MA, et al. CARD 2023: expanded curation, support for machine learning, and resistome prediction at the comprehensive antibiotic resistance database. Nucleic Acids Res. 2023;51D1:D690–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkac920
Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y, et al. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2005;33(Database issue):D325–328. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gki008
Johansson MHK, Bortolaia V, Tansirichaiya S, Aarestrup FM, Roberts AP, Petersen TN. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: mobileelementfinder. J Antimicrob Chemother. 2021;76(1):101–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jac/dkaa390
Bertelli C, Laird MR, Williams KP, Simon Fraser University Research, Computing G, Lau BY, Hoad G, et al. IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 2017;45(W1):W30–5. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkx343
Bi D, Xu Z, Harrison EM, Tai C, Wei Y, He X, et al. ICEberg: a web-based resource for integrative and conjugative elements found in bacteria. Nucleic Acids Res. 2012;40(Database issue):D621–626. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkr846
Wishart DS, Han S, Saha S, Oler E, Peters H, Grant JR, et al. PHASTEST: faster than PHASTER, better than PHAST. Nucleic Acids Res. 2023;51(W1):W443–50. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkad382
Liu Y, Du J, Lai Q, Zeng R, Ye D, Xu J, et al. Proposal of nine novel species of the Bacillus cereus group. Int J Syst Evol Microbiol. 2017;67(8):2499–508. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijsem.0.001821
Magiorakos AP, Srinivasan A, Carey RB, Carmeli Y, Falagas ME, Giske CG, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect. 2012;18(3):268–81. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1469-0691.2011.03570.x
Grant JR, Enns E, Marinier E, Mandal A, Herman EK, Chen CY, et al. Proksee: in-depth characterization and visualization of bacterial genomes. Nucleic Acids Res. 2023;51(W1):W484–92. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkad326
Singh N, Armstrong DG, Lipsky BA. Preventing foot ulcers in patients with diabetes. JAMA. 2005;293(2):217–28. https://doiorg.publicaciones.saludcastillayleon.es/10.1001/jama.293.2.217
Ndosi M, Wright-Hughes A, Brown S, Backhouse M, Lipsky BA, Bhogal M, et al. Prognosis of the infected diabetic foot ulcer: a 12-month prospective observational study. Diabet Med. 2018;35(1):78–88. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/dme.13537
Yang L, Rong GC, Wu QN. Diabetic foot ulcer: challenges and future. World J Diabetes. 2022;13(12):1014–34. https://doiorg.publicaciones.saludcastillayleon.es/10.4239/wjd.v13.i12.1014
Jneid J, Benamar S, Pagnier I, Levy P-Y, Lavigne J-P, La Scola B. Draft genome sequence of Providencia heimbachae, isolated from a diabetic foot ulcer. Genome Announcements. 2016;4(2). https://doiorg.publicaciones.saludcastillayleon.es/10.1128/genomea.00276-00216
Veve MP, Mercuro NJ, Sangiovanni RJ, Santarossa M, Patel N. Prevalence and predictors of Pseudomonas aeruginosa among hospitalized patients with diabetic foot infections open forum. Infect Dis. 2022;9(7). https://doiorg.publicaciones.saludcastillayleon.es/10.1093/ofid/ofac297
Jain SK, Barman R. Bacteriological profile of diabetic foot ulcer with special reference to drug-resistant strains in a tertiary care center in North-East India. Indian J Endocrinol Metabol. 2017;21(5):688–94. https://doiorg.publicaciones.saludcastillayleon.es/10.4103/ijem.IJEM_546_16
Saseedharan S, Sahu M, Chaddha R, Pathrose E, Bal A, Bhalekar P, et al. Epidemiology of diabetic foot infections in a reference tertiary hospital in India. Brazilian J Microbiol. 2018;49(2):401–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.bjm.2017.09.003
Saseedharan S, Sahu M, Chaddha R, Pathrose E, Bal A, Bhalekar P, et al. Epidemiology of diabetic foot infections in a reference tertiary hospital in India. Braz J Microbiol. 2018;49(2):401–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.bjm.2017.09.003
Abu-El-Azayem AK, Nashaat N, Dwedar RA, Fekry KM, Bassyouni RH, Hegab AS. Microbiological profile of diabetic foot infections %J. Microbes Infect Dis. 2024;5(4):1530–40. https://doiorg.publicaciones.saludcastillayleon.es/10.21608/mid.2024.295059.1977
Guan J, Bao C, Wang P, Jing Y, Wang L, Li X, et al. Genetic characterization of four groups of chromosome-borne accessory genetic elements carrying drug resistance genes in Providencia. Infect Drug Resist. 2022;15:2253–70. https://doiorg.publicaciones.saludcastillayleon.es/10.2147/IDR.S354934
Stock I, Wiedemann B. Natural antibiotic susceptibility of Providencia Stuartii, P. rettgeri, P. alcalifaciens and P. rustigianii strains. J Med Microbiol. 1998;47(7):629–42. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/00222615-47-7-629
Swiatlo E, Kocka FE. Inducible expression of an aminoglycoside-acetylating enzyme in Providencia stuartii. J Antimicrob Chemother. 1987;19(1):27–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jac/19.1.27
Yaghoubi S, Zekiy AO, Krutova M, Gholami M, Kouhsari E, Sholeh M, et al. Tigecycline antibacterial activity, clinical effectiveness, and mechanisms and epidemiology of resistance: narrative review. Eur J Clin Microbiol Infect Dis. 2022;41(7):1003–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10096-020-04121-1
Samonis G, Korbila IP, Maraki S, Michailidou I, Vardakas KZ, Kofteridis D, et al. Trends of isolation of intrinsically resistant to colistin Enterobacteriaceae and association with colistin use in a tertiary hospital. Eur J Clin Microbiol Infect Dis. 2014;33(9):1505–10. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10096-014-2097-8
Chun J, Rainey FA. Integrating genomics into the taxonomy and systematics of the bacteria and archaea. Int J Syst Evol Microbiol. 2014;64(Pt 2):316–24. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijs.0.054171-0
Chun J, Oren A, Ventosa A, Christensen H, Arahal DR, da Costa MS, et al. Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. Int J Syst Evol Microbiol. 2018;68(1):461–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijsem.0.002516
Mulet M, Lalucat J, Garcia-Valdes E. DNA sequence-based analysis of the pseudomonas species. Environ Microbiol. 2010;12(6):1513–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1462-2920.2010.02181.x
Hamed SM, Amer MA. Pseudocitrobacter cyperus, a novel bacterial species recovered from Cyperus alternifolius in Egypt. BMC Microbiol. 2025;25(1):20. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12866-024-03710-w
Lupande-Mwenebitu D, Khedher MB, Khabthani S, Rym L, Phoba MF, Nabti LZ, et al. First genome description of Providencia vermicola isolate bearing NDM-1 from blood culture. Microorganisms. 2021;9(8). https://doiorg.publicaciones.saludcastillayleon.es/10.3390/microorganisms9081751
Abdallah M, Balshi A. First literature review of carbapenem-resistant Providencia. New Microbes New Infect. 2018;25:16–23. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.nmni.2018.05.009
Oikonomou O, Liakopoulos A, Phee LM, Betts J, Mevius D, Wareham DW. Providencia stuartii isolates from Greece: co-carriage of cephalosporin (blaSHV-5, blaVEB-1), carbapenem (blaVIM-1), and aminoglycoside (rmtB) resistance determinants by a multidrug-resistant outbreak clone. Microb Drug Resist. 2016;22(5):379–86. https://doiorg.publicaciones.saludcastillayleon.es/10.1089/mdr.2015.0215
Sabtcheva S, Stoikov I, Ivanov IN, Donchev D, Lesseva M, Georgieva S, et al. Genomic characterization of carbapenemase-producing enterobacter hormaechei, serratia marcescens, citrobacter freundii, Providencia Stuartii, and Morganella Morganii clinical isolates from Bulgaria. Antibiot (Basel). 2024;13(5). https://doiorg.publicaciones.saludcastillayleon.es/10.3390/antibiotics13050455
Poirel L, Lartigue MF, Decousser JW, Nordmann P. ISEcp1B-mediated transposition of blaCTX-M in Escherichia coli. Antimicrob Agents Chemother. 2005;49(1):447–50. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/AAC.49.1.447-450.2005
Shawa M, Furuta Y, Mulenga G, Mubanga M, Mulenga E, Zorigt T, et al. Novel chromosomal insertions of ISEcp1-bla(CTX-M-15) and diverse antimicrobial resistance genes in Zambian clinical isolates of Enterobacter cloacae and Escherichia coli. Antimicrob Resist Infect Control. 2021;10(1):79. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13756-021-00941-8
Hamed SM, Mohamed HO, Ashour HM, Fahmy LI. Comparative genomic analysis of strong biofilm-forming Klebsiella pneumoniae isolates uncovers novel ISEcp1-mediated chromosomal integration of a full plasmid-like sequence. Infect Dis (Lond). 2024;56(2):91–109. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/23744235.2023.2272624
Munoz-Lopez M, Garcia-Perez JL. DNA transposons: nature and applications in genomics. Curr Genomics. 2010;11(2):115–28. https://doiorg.publicaciones.saludcastillayleon.es/10.2174/138920210790886871
Guillard T, Cambau E, Neuwirth C, Nenninger T, Mbadi A, Brasme L, et al. Description of a 2,683-base-pair plasmid containing QnrD in two Providencia rettgeri isolates. Antimicrob Agents Chemother. 2012;56(1):565–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/AAC.00081-11
Kraychete GB, Campana EH, Picao RC, Bonelli RR. qnrD-harboring plasmids in Providencia spp. Recovered from food and environmental Brazilian sources. Sci Total Environ. 2019;646:1290–2. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2018.07.378
Mazzariol A, Kocsis B, Koncan R, Kocsis E, Lanzafame P, Cornaglia G. Description and plasmid characterization of < em > qnrd determinants in < em > proteus mirabilis and < em > morganella Morganii. Clin Microbiol Infect. 2012;18(3):E46–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1469-0691.2011.03728.x
Fouad A, Gill CM, Simner PJ, Nicolau DP, Asempa TE. Cefepime in vivo activity against carbapenem-resistant enterobacterales that test as cefepime susceptible or susceptible-dose dependent in vitro: implications for clinical microbiology laboratory and clinicians. J Antimicrob Chemother. 2023;78(9):2242–53. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jac/dkad229
Li X, Zhang J, Wang J, Long W, Liang X, Yang Y, et al. Activities of Aztreonam in combination with several novel beta-lactam-beta-lactamase inhibitor combinations against carbapenem-resistant Klebsiella pneumoniae strains coproducing KPC and NDM. Front Microbiol. 2024;15:1210313. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2024.1210313
Mostafa SH, Saleh SE, Hamed SM, Aboshanab KM. Febrile illness of bacterial etiology in a public fever hospital in Egypt: high burden of multidrug resistance and WHO priority gram negative pathogens. Germs. 2022;12(1):75–85. https://doiorg.publicaciones.saludcastillayleon.es/10.18683/germs.2022.1308
Theuretzbacher U, Van Bambeke F, Canton R, Giske CG, Mouton JW, Nation RL, et al. Reviving old antibiotics. J Antimicrob Chemother. 2015;70(8):2177–81. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jac/dkv157
Pais SV, Kim E, Wagner S. Virulence-associated type III secretion systems in Gram-negative bacteria. Microbiol (Reading). 2023;169(6). https://doiorg.publicaciones.saludcastillayleon.es/10.1099/mic.0.001328
Blumer C, Kleefeld A, Lehnen D, Heintz M, Dobrindt U, Nagy G, et al. Regulation of type 1 fimbriae synthesis and biofilm formation by the transcriptional regulator lrha of Escherichia coli. Microbiol (Reading). 2005;151(Pt 10):3287–98. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/mic.0.28098-0
Murphy CN, Mortensen MS, Krogfelt KA, Clegg S. Role of Klebsiella pneumoniae type 1 and type 3 fimbriae in colonizing silicone tubes implanted into the bladders of mice as a model of catheter-associated urinary tract infections. Infect Immun. 2013;81(8):3009–17. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/IAI.00348-13
Berry JL, Pelicic V. Exceptionally widespread nanomachines composed of type IV Pilins: the prokaryotic Swiss army knives. FEMS Microbiol Rev. 2015;39(1):134–54. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/femsre/fuu001
Pouget C, Dunyach-Remy C, Pantel A, Schuldiner S, Sotto A, Lavigne JP. Biofilms in diabetic foot ulcers: significance and clinical relevance. Microorganisms. 2020;8(10). https://doiorg.publicaciones.saludcastillayleon.es/10.3390/microorganisms8101580
Zegadlo K, Gieron M, Zarnowiec P, Durlik-Popinska K, Krecisz B, Kaca W, et al. Bacterial motility and its role in skin and wound infections. Int J Mol Sci. 2023;24(2). https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms24021707
Green ER, Mecsas J. Bacterial secretion systems: an overview. Microbiol Spectr. 2016;4(1). https://www.doi.org/10.1128/microbiolspec.VMBF-0012-2015
Klein JA, Predeus AV, Greissl AR, Clark-Herrera MM, Cruz E, Cundiff JA, et al. Pathogenic diversification of the gut commensal Providencia alcalifaciens via acquisition of a second type III secretion system. Infect Immun. 2024;92(10):e0031424. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/iai.00314-24
Rippke F, Berardesca E, Weber TM. pH and microbial infections. Curr Probl Dermatol. 2018;54:87–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1159/000489522
Li X, Lu N, Brady HR, Packman AI. Biomineralization strongly modulates the formation of proteus mirabilis and Pseudomonas aeruginosa dual-species biofilms. FEMS Microbiol Ecol. 2016;92(12). https://doiorg.publicaciones.saludcastillayleon.es/10.1093/femsec/fiw189
Konieczna I, Zarnowiec P, Kwinkowski M, Kolesinska B, Fraczyk J, Kaminski Z, et al. Bacterial urease and its role in long-lasting human diseases. Curr Protein Pept Sci. 2012;13(8):789–806. https://doiorg.publicaciones.saludcastillayleon.es/10.2174/138920312804871094
Lv J, Zhu J, Wang T, Xie X, Wang T, Zhu Z, et al. The role of the two-component QseBC signaling system in biofilm formation and virulence of hypervirulent Klebsiella pneumoniae ATCC43816. Front Microbiol. 2022;13:817494. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2022.817494
Flannery EL, Mody L, Mobley HL. Identification of a modular pathogenicity island that is widespread among urease-producing uropathogens and shares features with a diverse group of mobile elements. Infect Immun. 2009;77(11):4887–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/IAI.00705-09
Flannery EL, Antczak SM, Mobley HL. Self-transmissibility of the integrative and conjugative element ICEPm1 between clinical isolates requires a functional integrase, relaxase, and type IV secretion system. J Bacteriol. 2011;193(16):4104–12. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/JB.05119-11
Mil-Homens D, Fialho AM. Trimeric autotransporter adhesins in members of the Burkholderia cepacia complex: a multifunctional family of proteins implicated in virulence. Front Cell Infect Microbiol. 2011;1:13. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fcimb.2011.00013
Alamuri P, Lower M, Hiss JA, Himpsl SD, Schneider G, Mobley HL. Adhesion, invasion, and agglutination mediated by two trimeric autotransporters in the human uropathogen proteus mirabilis. Infect Immun. 2010;78(11):4882–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/IAI.00718-10
Acknowledgements
Not applicable.
Funding
Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).
The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.
Author information
Authors and Affiliations
Contributions
SMH, MAA, MMD, RM, AA, AIG, and INH: Conceptualization, Methodology, Investigation, Resources; SMH and MAA: Writing– Original Draft, MMD, RM, AA, AIG, and INH: Writing– Review & Editing.
Corresponding author
Ethics declarations
Ethical approval
The study was conducted in accordance with the Declaration of Helsinki and ethical approval was obtained from the Ethics Committee of the Faculty of Pharmacy, MSA University (Approval Number: M6/HEC6/S2024PD), and the National Institute of Diabetes and Endocrinology (Approval Number: IDE 00315). The need for patients’ consent to participate was waived by the Institutional Review Boards of MSA university and NIDE as the isolates were taken as part of the routine clinical care of patients and no personal data were published.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hamed, S.M., Darwish, M.M., Monir, R. et al. Providencia pseudovermicola sp. nov.: redefining Providencia vermicola and unveiling multidrug-resistant strains from diabetic foot ulcers in Egypt. BMC Microbiol 25, 238 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12866-025-03927-3
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12866-025-03927-3