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, 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 a mild limp, tenderness, and swelling in the proximal 1/3 of her left Tibia. Blood analysis showed mild leukocytosis (11,400 cells/mL, range 4000-11000) and normal C-reactive protein (CRP) (2.2 mg/L, range 0-5). Magnetic resonance imaging (MRI) showed intra-medullary edema in the proximal tibia with enhancement following Gadolinium injection, a breach of the frontal and lateral bone cortex, and an abscess formation in the subdermal fat. The patient was admitted for a debridement operation, with no antimicrobial treatment before surgery. Pathological bone was noted by the orthopedic team, and multiple microbiology and pathology specimens were collected. The pathology report showed 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 bone and soft tissue samples.
Pathogen identification
The biopsy samples were plated on TSA 5% sheep blood agar, chocolate agar, and MacConkey agar plates (Hylab, Rehovot, Israel) at 37°c in aerobic conditions overnight and presented white-grey circular smooth colonies. Gram staining showed gram-negative coccobacilli (Fig. 1). Phenotypic tests (Methods in Supp. material A) showed this strain is catalase-negative, oxidase-negative, coagulase-negative, and Indole-negative. Colonies were tested using matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) on VITEK® MS (bioMérieux, Marcy-l'Étoile, France) which resulted in no identification of any known species. Antibiotic susceptibility based on gram-negative CLSI breakpoints and biochemical tests (Supp.Table 1) were performed on VITEK®2 (bioMérieux, Marcy-l'Étoile, France) and 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 colonies and the 16s rRNA gene was amplified and sequenced by Sanger sequencing. BLASTN analysis of the resulting sequence against GenBank generated good hits to 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 sized 1,854,057, a total length of 4,314,591 bps, and 55.3% GC content. Taxonomic analyses performed using TYGS (Type strain genome server [6] and GTDB-Tk Classify genomes 2.3.2 [7] could not assign a known species and found the closest species found were Pantoea superficialis (d0=59.5% & d4= 33% dDDH 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 our strain belongs to a still unknown species of Pantoea, which we have named Pantoea osteomyelitidis.
Pan-genome analysis was performed using Roary based on 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 less of the genomes).
The phylogenetic tree based on core genes (Fig. 2) showed high diversity within Pantoea, with many distinct lineages. Pantoea osteomyelitidissp.nov clustered with P. superficialis and P. mediterraneensis. However, there is 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.
Average nucleotide identity values based on BLASTN (ANIb) are shown as a heatmap in Fig.3. This analysis aligns mostly with the subclades within the phylogenetic tree generated according to core genome, as can be seen by values of over 80%.
Fig.3 also presents the isolation source of each species according to all reports made in GenBank. Overall, 12 species have been isolated from human samples, some depicted as pathogens and others recovered from healthy tissues. Of these, only 3 have been isolated from human samples only (P. osteomyelitidissp. 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).
The gene prediction analysis of Pantoea osteomyelitidissp.nov using Prokka version 1.14.6 [9] identified a total of 4,084 genes, comprising 3,294 proetin-coding sequences. We tested for the presence of antibiotic-resistance and virulence genes (methods in Supp. material A) and identified oqxB and CRP presence. OqxB codes for a multidrug efflux RND transporter permease subunit. Overexpression of OqxB can confer resistance to various antibiotics [10,11]. CRP gene codes for 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, NlpI codes for Lipoprotein NlpI [13]. OmpA gene codes for porin outer membrane protein A, known for its multiple functions including biofilm formation and a role in F plasmid cell conjugation [14]. EC55989_3335 codes for Hcp family type VI secretion system effector [13]. Our results from the Pan-genome analysis revealed the presence of OqxB, CRP, Nlpl, FliG, FliM, FliP and ompA within the core gene group, indicating these genes are evolutionary conserved within this genus (Supp. table.2).
An analysis of potential pathogenic proteins using PathogenFinder suggested a non-human pathogen with a 0.47 likelihood of this species to be pathogenic, by identifying 16 pathogenic protein families identified, compared to 17 non-pathogenic protein families. 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 so far been isolated only from insects.Although this analysis does not suggest a clear pathogenicity potential for Pantoea osteomyelitidissp.nov., other Pantoea species that are considered opportunistic human pathogens were also determined as non-human pathogens (Supp. table 3).
Patient treatment and follow-up
Four months following 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 specimens were again positive for Pantoea osteomyelitidissp.nov. She was treated with oral Ciprofloxacin 750 mg bid for four weeks. Six months after her surgery, the patient reported complete healing.