A generally healthy 37-year-old woman first presented to our hospital complaining of five months of left leg pain exacerbated by physical activity, tenderness, and swelling, without fever. Past medical history was notable for a tibial fracture at the age of 15, which was sustained by a fall over a piece of wood in the field and accompanied by a skin laceration. This fracture was treated conservatively, and the patient did not recall if antimicrobial treatment was prescribed. After an asymptomatic period of approximately 20 years, she began feeling a dull ache in her shin that was exacerbated over two years. Physical examination revealed mild limp, tenderness, and swelling in the proximal 1/3 of her left tibia. Blood analysis revealed mild leukocytosis (11,400 cells/mL, range 4000–11000) and normal C-reactive protein (CRP) (2.2 mg/L, range 0–5) levels. Magnetic resonance imaging (MRI) showed intramedullary edema in the proximal tibia with enhancement following gadolinium injection, a breach of the frontal and lateral bone cortex, and abscess formation in the subdermal fat. The patient was admitted for debridement, with no antimicrobial treatment before surgery. Pathological bone was noted by the orthopedic team, and multiple microbiological and pathological specimens were collected. Pathology revealed fragments of bone and fibrous tissue with mostly lymphoplasmacytic inflammatory infiltrate, hemorrhage, and focally purulent granulation tissue. Few areas showed new bone formation with osteoid and surrounding osteoblasts, features compatible with chronic osteomyelitis. Bacterial cultures were positive in 4/8 of the bone and soft tissue samples.
The biopsy samples were plated on TSA 5% sheep blood agar, chocolate agar, and MacConkey agar plates (Hylab, Rehovot, Israel) at 37°C under aerobic conditions overnight, and white‒gray circular smooth colonies formed. Gram staining revealed gram-negative coccobacilli (Fig. 1). Phenotypic tests (Methods in Supp. material A) showed that this strain is catalase negative, oxidase negative, coagulase negative, and indole negative. Colonies were tested via matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) via a VITEK® MS (bioMérieux, Marcy-l'Étoile, France), which resulted in the identification of no known species. Antibiotic susceptibility tests, which are based on gram-negative CLSI breakpoints and biochemical tests (Supp. Table 1), were performed on VITEK®2 (bioMérieux, Marcy-l'Étoile, France), which indicated sensitivity to all antibiotics tested, including amoxicillin/clavulanic acid, cefuroxime, chloramphenicol, ampicillin, ceftazidime, ceftriaxone, ciprofloxacin, gentamicin, piperacillin/tazobactam, trimethoprim/sulfamethoxazole, meropenem, amikacin, and ertapenem.
DNA was extracted from the colonies, and the 16S rRNA gene was amplified and sequenced via Sanger sequencing. BLASTN analysis of the resulting sequence against GenBank generated good hits for different species of Erwiniaceae. Therefore, we could only conclude that this strain belongs to Erwiniaceae. Whole-genome sequencing reads (Methods described in Supp. material A) were assembled to generate 17 contigs (N50 = 570,650), the longest contig size was 1,854,057, the total length was 4,314,591 bps, and the GC content was 55.3%. Taxonomic analyses performed via TYGS (Type strain genome server [6] and GTDB-Tk Classify genomes 2.3.2 [7] could not assign a known species and revealed that the closest species found were Pantoea superficialis (d0 = 59.5% & d4 = 33% dDDH (DNA-DNA hybridization) values by TYGS and 87.6% average nucleotide identity (ANI) by GTDB-Tk), followed by Pantoea mediterraneensis (d0 = 46%, d4 = 26.9% dDDH and 84% ANI). dDDH values of over 70% and ANI values of over 95% indicate the same species [6, 8], suggesting that our strain belongs to a still unknown species of Pantoea, which we have named Pantoea osteomyelitidis.
Pan-genome analysis was performed via Roary on the basis of the core genome alignment of Pantoea osteomyelitidis. sp. nov. with known Pantoea species according to Crosby et al. [5] and additional recently identified species (GenBank accessions in Supp. Table 1, full methods in Supp. material A). This analysis identified a total of 47,455 genes: 1,464 core genes (found in 90% of the genomes), 4,185 shell genes (found in 15–89% of the genomes) and 41,806 cloud genes (found in 15% or fewer of the genomes).
The phylogenetic tree based on core genes (Fig. 2) revealed high diversity within Pantoea, with many distinct lineages. Pantoea osteomyelitidis sp. nov. clustered with P. superficialis and P. mediterraneensis. However, there was a divergence between the common ancestor of this clade and the common ancestor of all other Pantoea species on the tree, indicating the evolutionary distance between this clade and the other clades.
The average nucleotide identity values obtained via BLASTN (ANIb) are shown as a heatmap in Fig. 3. This analysis aligns mostly with the subclades within the phylogenetic tree generated according to the core genome, as shown by values greater than 80%.
Figure 3 also presents the isolation source of each species according to all reports in GenBank. Overall, 12 species have been isolated from human samples, some depicted as pathogens and others recovered from healthy tissues. Of these, 3 have been isolated only from human samples (P. osteomyelitidis sp. nov., P. mediterraneensis and P. mediterraneensis_A), and the rest have been isolated from additional hosts or environments (P. eucrina, P. dispersa, P. conspicua, P. brenneri, P. piersonii, P. septica, P. alvi, P. anthophila and P. agglomerans).
Gene prediction analysis of Pantoea osteomyelitidis sp. nov. via Prokka version 1.14.6 [9] identified a total of 4,084 genes, comprising 3,294 protein-coding sequences. We tested for the presence of antibiotic resistance and virulence genes (Methods in Supp. material A) and identified the presence of oqxB and CRP. OqxB codes for a multidrug efflux RND transporter permease subunit. OqxB overexpression can confer resistance to various antibiotics [10, 11]. The CRP gene encodes a regulator that represses MdtEF multidrug efflux pump expression, and mutations in this gene were found to be associated with increased resistance [12]. The following virulence factors were detected: FliG, FliM, FliP, Nlpl, OmpA, and EC55989_3335. FliG and FliM are associated with the flagellar motor switch, FliP is involved in flagellar biosynthesis, and NlpI codes for lipoprotein NlpI [13]. The OmpA gene encodes porin outer membrane protein A, which is known for its multiple functions, including biofilm formation and a role in F plasmid cell conjugation [14]. EC55989_3335 encodes an Hcp family type VI secretion system effector [13]. Our results from the pangenome analysis revealed the presence of OqxB, CRP, Nlpl, FliG, FliM, FliP and ompA within the core gene group, indicating that these genes are evolutionarily conserved within this genus (Supp. Table 2).
An analysis of potential pathogenic proteins via PathogenFinder suggested a nonhuman pathogen with a 0.47 likelihood of being pathogenic, as 16 pathogenic protein families were identified, whereas 17 nonpathogenic protein families were identified. By examining all the Pantoea genomes used in our analysis, only 3 were determined by PathogenFinder as human pathogens: P. alvi (0.576), P. mediterraneensis_A (0.535) and P. bathycoeliaeme (0.585). While P. alvi and P. mediterraneensis_A have been isolated from healthy human samples, P. bathycoeliaeme has thus far been isolated only from insects. Although this analysis does not suggest a clear pathogenicity potential for Pantoea osteomyelitidis sp. nov., other Pantoea species that are considered opportunistic human pathogens were also determined to be nonhuman pathogens (Supp. Table 3).
Patient treatment and follow-up
Four months after her open bone biopsy, and upon identification of the causative organism, the patient underwent extensive debridement with bone cement supplement using calcium sulfate resorbable beads with tobramycin. One out of 5 tissue samples was again positive for Pantoea osteomyelitidis sp. nov. She was treated with 750 mg bid oral ciprofloxacin for four weeks. Six months after her surgery, the patient reported complete healing.