Reporting two novel Kluyvera species, Kluyvera huaxiensis and Kluyvera chengduensis, isolated from human sputa

Introduction: Two Gram-stain negative, facultatively anaerobic, rod-shaped, motile, and non-spore forming bacterial strains, 142053^T and 142359^T, were isolated from sputum samples obtained from different patients.

Methods: Growth of strain 142053^T and 142359^T occurred at 8–42 °C (optimum 28–37 °C), pH 5.0–9.0 (optimum 6.0–8.0) and NaCl concentrations of 0–5% (w/v) (optimum 0–4%). Comparison of 16S rRNA gene sequences showed that strain 142053^T and 142359^T exhibited high similarity to Kluyvera intermedia NBRC 102594^T (99.48%) and Kluyvera georgiana ATCC 51603^T (98.68%), respectively.

Results: Phylogenomic tree of species from family Enterobacteriaceae revealed that strain 142053^T and 142359^T formed distinct branches within the Kluyvera clade. The digital DNA–DNA hybridization between 142053^T, 142359^T, and all type strains of Kluyvera species were ≤39.5% and ≤54.7%, respectively, and the corresponding average nucleotide identity were ≤90.56% and ≤93.85%, both below the species delineation threshold. The two strains differ phenotypically from other Kluyvera species by their lipase and arginine dihydrolase utilization and from each other by ᴅ-arabitol utilization. The predominant fatty acids identified in both strains were C_16:0, summed feature 3 (C_16:1ω7c/C_16:1ω6c), summed feature 8 (C_18:1ω7c) and C_17:0 cyclo.

Discussion: Based on the phenotypic, chemotaxonomic, phylogenetic and genomic evidence, strains 142053^T and 142359^T represent two novel species of the genus Kluyvera, for which the names Kluyvera huaxiensis sp. nov. and Kluyvera chengduensis sp. nov. are proposed, with type strains 142053^T (= GDMCC 1.4845^T = JCM 37413^T) and 142359^T (= GDMCC 1.5094^T = JCM 37737^T), respectively.

Introduction

Kluyvera is a genus of Gram-negative, rod-shaped bacteria in the family Enterobacteriaceae (Farmer et al., 1981). Currently, the genus Kluyvera contains five validly published species-Kluyvera ascorbata, Kluyvera cryocrescens, Kluyvera georgiana, Kluyvera intermedia, and Kluyvera sichuanensis and five unnamed genomospecies (1, 2, 3, 4/6, and 5) (Farmer et al., 1981; Müller et al., 1996; Pavan et al., 2005; Liu et al., 2020; Rodriguez and Gutkind, 2024). Kluyvera typically inhabits sewage, hospital sinks, and human sputum, urine, or the gastrointestinal tract (Farmer et al., 1981; Sarria et al., 2001; Li et al., 2019; Alsuheili et al., 2024; Endo et al., 2024). Currently, there is an increasing number of reports on human infections caused by this bacterium, including oral infections, bloodstream infections, enteritis, biliary tract infections, abdominal abscesses, and urinary tract infections. Clinically, these infections pose challenges such as therapeutic difficulty, prolonged treatment duration, and adverse outcomes, including patient death (Mutoh et al., 2019; Zou et al., 2019; Akiki et al., 2023; Inoue et al., 2023; Ulloa-Clavijo et al., 2023). The most frequently reported Kluyvera species associated with human infections are K. ascorbata and K. cryocrescens (Wong, 1987; Lin et al., 2002; Yoshino et al., 2016; Zou et al., 2019), although clinical cases involving K. intermedia and K. georgiana have also been reported in recent years (Inoue et al., 2023; Endo et al., 2024; Xu et al., 2025). Kluyvera possesses a large and significantly flexible gene pool (Yu et al., 2023). Notably, the extended-spectrum β-lactamase gene bla_CTX-M and the fosfomycin resistance gene fosA originated from this genus (Humeniuk et al., 2002; Rodríguez et al., 2004; Ito et al., 2018; Rodriguez et al., 2018). Alarmingly, several plasmid-borne carbapenemase genes, including bla_KPC (Tiwari et al., 2022), bla_OXA-48 (Hernández-García et al., 2018), and bla_NDM (Loiodice et al., 2023), and colistin resistance genes mcr (Strepis et al., 2021) have also been identified in Kluyvera species. Taken together, Kluyvera represents an emerging clinical and epidemiological concern due to its dual role as an opportunistic pathogen and as a reservoir of antimicrobial resistance genes. In this study, we conducted a taxonomic analysis on two clinical strains, 142053^T and 142359^T, which were found to represent two novel species within the genus Kluyvera.

Materials and methods

Strains and the study

Strain 142053^T and 142359^T were isolated from the sputum samples of two different inpatients at West China Hospital in Chengdu (30.6800°N, 104.0811°E), China, in April and August 2024, respectively. The two sputum samples were collected as part of routine care for managing suspected pulmonary infection in the two patients. The sputum samples were vortexed with 0.9% normal saline at a 1:1 ratio for 10 s. The mucoid portion was then picked up with a sterile cotton swab and inoculated onto 5% sheep blood agar (Huankai; Guangzhou, China) plates using the quadrant streak method. The plates were incubated at 37 °C for 48 h, and the colonies of the two isolates were round, smooth, convex, and white. The two isolates were stored long-term at −80 °C in nutrient broth (Huankai) containing 20% glycerol. Preliminary species identification was conducted using matrix-assisted laser desorption ionization time of flight mass spectrum (MALDI-TOF; Bruker; Billerica, MA) with the MBT Compass Library 2023 according to the manufacturer’s instructions. The study has been approved by the Ethical Committee of West China Hospital with waiver of informed consent as this study was to characterize bacterial strains.

Analysis of the 16S rRNA gene sequence

Genomic DNA of strain 142053^T and 142359^T was extracted using a bacterial genomic DNA extraction kit (Tiangen; Beijing, China). The 16S rRNA gene was amplified using universal primers 27F and 1492R (Moreno et al., 2002), and sequenced bidirectionally using Sanger sequencing (Sangon Biotech; Shanghai, China) for preliminary taxonomic identification. The assembled 16S rRNA gene sequences were compared with the NCBI 16S rRNA database using BLAST, applying a 98% identity threshold to identify closely related genera. Along with those retrieved from type strains of closely related genera, 16S rRNA gene sequences of 142053^T and 142359^T were aligned using Clustal Omega v1.2.4 (Madeira et al., 2019).

Whole genome sequencing and analysis

Genomic DNA was extracted using the same method as described above. Whole genome sequencing of strain 142053^T and 142359^T was performed by NovaSeq 6,000 system (Illumina; San Diego, CA, USA). Reads were de novo assembled into contigs using SPAdes v4.0.0 (Bankevich et al., 2012) in isolate modes. The draft genomes of strain 142053^T and 142359^T were compared with those of type strains of Kluyvera species to assess genomic relatedness using average nucleotide identity (ANI) and in silico DNA–DNA hybridization (isDDH). ANI and isDDH values were calculated using FastANI v1.34 (Jain et al., 2018), OrthoANI (Yoon et al., 2017) and Genome-to-Genome Distance Calculator (formula 2) (Meier-Kolthoff et al., 2013b) with the default settings, respectively. POCP was calculated as described previously (Qin et al., 2014). Sequences were submitted to the online TYGS platform¹ to assess their identification at the genus or species level (Meier-Kolthoff and Göker, 2019).

Genome sequences of type strains of all Enterobacteriaceae species were retrieved from GenBank. These sequences were annotated using Prokka v1.14.5 (Seemann, 2014). Orthologues were identified via PIRATE v1.0.5 (Bayliss et al., 2019) to represent the core genome. A maximum-likelihood (ML) phylogenetic tree was inferred from core-genome using IQ-tree v2.1.2 (Minh et al., 2020) with the GTR + GAMMA model and 1,000 bootstrap replicates. Additionally, neighbor-joining (NJ) and maximum parsimony (MP) trees were constructed with 1,000 bootstrap replicates (Kumar et al., 2018) to assess the robustness of the phylogeny. All genomic data analyzed in this study are available in Supplementary 1. Antimicrobial resistance determinants were predicted using ABRicate v1.0.0² against ResFinder databases. Plasmid replicon types were determined using the PlasmidFinder tool available from the Center for Genomic Epidemiology.³

Physiology and chemotaxonomy

Cell motility was examined by observing the bacterial growth and diffusion on the deep semi-solid nutrient medium of 0.3% (w/v) agar (Hopebio; Qingdao, China). Anaerobic growth was examined by streaking the bacterial cultures on brain heart infusion agar plates and placed in the GasPakTM EZ anaerobic bag (BD; Franklin Lakes, NJ, USA) at 37 °C for 3 days. After incubation in nutrient broth at 37 °C for 3 days, the flagella of strain 142053^T and 142359^T were observed using an H-7650 transmission electron microscope (Hitachi; Tokyo, Japan). Bacterial growth was examined in 5 mL aliquots of nutrient broth. These aliquots were dispensed into tubes with an inner diameter of 16 mm and incubated at temperatures of 4, 8, 18, 28, 32, 37, 42, 45, 48, and 50 °C for 7 days. Meanwhile, salt and pH tolerances were measured using nutrient broth at 37 °C for 7 days at different NaCl concentrations (0, 0.5, 1, 2, 3, 4, 5, 7.5, 10, and 15%, w/v) and at a pH unit of 4.0–12.0 (in increments of 1.0 unit), respectively (Hu et al., 2017). Catalase activity was tested by examining the production of bubbles after addition of 3% (v/v) hydrogen peroxide solution and oxidase activity was measured by using 1% tetramethyl-p-phenylenediamine dihydrochloride solution.

DNase activity was detected using DNase agar (Solarbio; Beijing, China). After incubation at 30 °C for 3 days, 1 M HCl was added for the detection. Malonate, phenylalanine deaminase, and KCN experiments were carried out using biochemical identification tubes (Huankai). Commercially available API 20E, API 50CH, and API ZYM kits (bioMérieux; Marcy l’Etoile, France) were employed to determine biochemical characteristics and enzyme activities. The assays were conducted in accordance with the manufacturer’s instructions, with Escherichia coli strain ATCC 25922 and Pseudomonas aeruginosa strain ATCC 27853 serving as controls. All experiments were performed in triplicate.

Fatty acid analysis

Cellular fatty acid analysis was conducted by the Guangdong Institute of Microbiology (Guangdong, China). In brief, bacteria were inoculated onto nutrient agar plates and incubated at 37 °C for 24 h. Fatty acid methyl esters were extracted and then analyzed via gas chromatography following the guidelines of the Sherlock Microbial Identification System (MIDI Inc.; Newark, DE, USA), as detailed in previous reports (Sasser, 1990). The peaks were automatically integrated, and the proportions of fatty acids were calculated using the MIDI identification database RTSBA6 (version 6.00; MIDI Inc.).

Antimicrobial susceptibility testing

Minimum inhibitory concentrations (MICs) of amikacin, ampicillin, ampicillin-sulbactam, aztreonam, ceftriaxone, ceftazidime, cefepime, cefotaxime, cefuroxime, chloramphenicol, ciprofloxacin, colistin, imipenem, meropenem, piperacillin-tazobactam, sulfamethoxazole-trimethoprim, and tigecycline against these two isolates was determined using the broth microdilution method in accordance with the guidelines of the Clinical and Laboratory Standards Institute (CLSI) (CLSI, 2025). Cation Adjusted Mueller-Hinton Broth (Oxoid, Basingstoke, UK) was used, and antimicrobial powders (purity ≥98%) were purchased from Meilunbio (Dalian, China) or Macklin (Shanghai, China). The breakpoints for colistin and tigecycline were interpreted based on the criteria of the European Committee on Antimicrobial Susceptibility Testing (EUCAST);⁴ otherwise, the breakpoints defined by the CLSI were applied.

Results

Two novel Kluyvera species were identified

The 16S rRNA gene sequence of strain 142053^T showed highest similarity to that of K. intermedia NBRC 102594^T (99.48%) and Citrobacter freundii ATCC 8090^T (99.22%), while that of strain 142359^T was most similar to K. georgiana ATCC 51603^T (98.68%) and Klebsiella electrica DSM 102253^T (98.11%). A maximum-likelihood tree was constructed based on 16S rRNA gene sequences, incorporating the sequences of the genera Citrobacter, Klebsiella, and Pseudocitrobacter, which are closely related to Kluyvera (Liu et al., 2020). Both strains clustered within the Kluyvera clade. However, they each form relatively long branches distinct from other Kluyvera species, suggesting potential taxonomic novelty (Supplementary Figure S1). Given the limited resolution of 16S rRNA gene-based analyses for differentiating species within the Enterobacteriaceae (Church et al., 2020), whole-genome sequencing was performed to confirm their phylogenomic placement and assess species-level distinctiveness.

Sequencing of strain 142053^T yielded 11,580,316 paired-end reads, representing 1.74 Gb of data, with a Q30 score of 93.5% and an average quality score of 26.8. These reads were then assembled into 80 contigs (N50 = 183,280 bp). The draft genome of strain 142053^T is 5,003,793 bp in size, with a G + C content of 53.49 mol%. Strain 142359^T produced 8,187,706 reads, representing 1.23 Gb of data, with a Q30 score of 95% and an average quality score of 29.2, which were assembled into 52 contigs (N50 = 341,410 bp). The draft genome of strain 142359^T is 4,907,633 bp in size, with a G + C content of 54.48 mol%. The POCP values for strains 142053^T and 142359^T against the reference strains of all Kluyvera species ranged from 75.82 to 87.03% and from 76.49 to 86.15%, respectively, supporting their assignment to the genus Kluyvera. The ANI values between strain 142053^T and the type strains of all other Kluyvera species ranged from 85.80 to 90.56%, while strain 142359^T showed ANI values between 85.45 and 93.85% relative to these type strains. All values were below the ≥95–96% ANI threshold recommended for species delineation (Meier-Kolthoff et al., 2013b) (Table 1). The OrthoANI values were also below the 95% (Supplementary Figure S2). Correspondingly, the isDDH values between strain 142053^T and other Kluyvera type strains ranged from 28.8 to 39.5%, and those for strain 142359^T ranged from 27.8 to 54.7%, well below the ≥ 70.0% threshold for species delineation (Richter and Rossello-Mora, 2009; Meier-Kolthoff et al., 2013a). Notably, strain 142053^T exhibited high genomic relatedness to the unnamed Kluyvera genomospecies 5, with isDDH and ANI values of 87.4 and 98.39%, respectively (Table 1). These results indicate that K. genomospecies 5 and K. huaxiensis represent the same species. Therefore, K. genomospecies 5 should be formally designated as K. huaxiensis.

Table 1

Table 1. Average nucleotide identity (ANI), in silico DNA–DNA hybridization (isDDH) and POCP values between strain 142053^T, 142359^T and the reference strains of other Kluyvera species.

A maximum-likelihood phylogenetic tree constructed using 265 core genes from the family Enterobacteriaceae (Figure 1) placed both strain 142053^T and 142359^T within a clade comprising members of genus Kluyvera. For independent validation, NJ, MP and TYGS-based phylogenomic trees were reconstructed (Supplementary Figures S3–S5), and the results were consistent with the ML phylogeny. This Kluyvera clade was well separated from other related genera, including the most closely related genus Pseudocitrobacter, further supporting that strain 142053^T and 142359^T belong to the genus Kluyvera. The ML phylogeny inferred from 1,668 core genes within the genus Kluyvera positioned strains 142053^T and 142359^T in two distinct clades, supporting their classification as novel species, as both were clearly separated from all currently recognized Kluyvera species (Figure 2). Taken together, these genomic and phylogenomic analyses provide robust evidence that strains 142053^T and 142359^T represent two novel species within the genus Kluyvera.

Figure 1

Figure 1. Phylogenomic tree of strain 142053^T, 142359^T, and other genera in the family Enterobacteriaceae. The core genome tree was constructed via the maximum-likelihood method based on 1,000 resamplings. Bar, 0.5 substitutions per nucleotide position.

Figure 2

Figure 2. Phylogenomic tree of strain 142053^T, 142359^T, and other species in the genera Kluyvera. The core genome tree was constructed via the maximum-likelihood method based on 1,000 resamplings. Bar, 0.5 substitutions per nucleotide position.

The physiological and biochemical characteristics of strain 142053^T and 142359^T are listed in Table 2. Briefly, both strains are positive for indole production, nitrate reduction, citrate utilization, and the activities of lysine decarboxylase, arginine dihydrolase, ornithine decarboxylase, β-galactosidase (ONPG), and catalase. These strains can assimilate aesculin, amygdalin, ᴅ-cellobiose, ᴅ-glucose, ᴅ-galactose, ᴅ-lactose, ᴅ-maltose, ᴅ-melibiose, ᴅ-sorbitol, ᴅ-trehalose, ᴅ-xylose, ʟ-arabinose, ʟ-rhamnose, mannitol, malonate, glycerol, methyl α-ᴅ-glucopyranoside, and salicin. Conversely, they are negative for acetoin production (Voges-Proskauer), DNase activity, H₂S production, oxidase activity, phenylalanine deaminase activity, and the activities of gelatinase and urease. Moreover, neither strain can utilize adonitol, dulcitol, erythritol, or inositol. In comparison with other species in genus Kluyvera, strains 142053^T and 142359^T display unique physiological characteristics. Both show weakly positive lipase activity and are positive for arginine dihydrolase. Additionally, the two strains can be distinguished from each another by ᴅ-arabitol utilization, with 142359^T but not 142053^T can metabolize ᴅ-arabitol. These distinct traits serve as markers for differentiating strains 142053^T and 142359^T from other Kluyvera species. The cellular fatty acid profiles of two novel strains were analyzed and compared with that of K. sichuanensis GDMCC 1.1872^T (Supplementary Table S1). The fatty acid composition of other Kluyvera species has not yet been reported. The major cellular fatty acids of the two strains are C_16:0, summed feature 3 (C_16:1ω7c/C_16:1ω6c), summed feature 8 (C_18:1ω7c), and C_17:0 cyclo. Transmission electron micrographs of strains 142053^T and 142359^T are shown in Supplementary Figsures S6, S7, respectively.

Table 2

Table 2. Biochemical characteristics of strain 142053^T, 142359^T, and other species within the genus Kluyvera.

The two strains were resistant to ampicillin but susceptible to ampicillin-sulbactam, aztreonam, ceftriaxone, ceftazidime, cefepime, cefotaxime, ciprofloxacin, imipenem, meropenem, piperacillin-tazobactam, amikacin, colistin, and tigecycline. Strain 142053^T was susceptible to cefuroxime but resistant to sulfamethoxazole-trimethoprim, while strain 142359^T showed the opposite pattern (Supplementary Table S2). Like other Kluyvera species, both strains possess an intrinsic chromosomal bla_CTX-M gene. Strain 142053^T encodes CTX-M-3, which is closely related to CTX-M-37 (differing by three amino acids) from Kluyvera genomospecies 5 (now to be considered K. huaxiensis) (Rodríguez et al., 2021), whereas strain 142359^T harbors the novel variant CTX-M-283 (GenBank accession no. PV575350). In addition, strain 142053^T also carries fosfomycin resistance gene fosA3, sulfonamide-resistance gene sul1, trimethoprim-resistance gene dfrA21, and a plasmid with an IncFIB (pHCM2) replicon, which likely accounts for the last two resistance markers.

Description of Kluyvera huaxiensis sp. nov.

Kluyvera huaxiensis (hua. xi. en’sis N. L. fem. Adj. huaxiensis, referring to West China [“Huaxi” in Chinese], where the type strain was recovered).

The strain 142053^T is Gram-negative, facultatively anaerobic, motile, and non-spore-forming rods, measuring 0.6–0.8 μm in width and 1.0–2.5 μm in length. After 12-h incubation on nutrient agar, colonies are round, smooth, convex, and white. Growth occurs at 8–42 °C (optimum, 28–37 °C), pH 5.0–9.0 (optimum, 6.0–8.0), and NaCl concentrations from 0 to 5% (w/v) (optimum, 0–4%). Strain 142053^T is positive for indole production, nitrate reduction, citrate utilization, and enzymatic activities including lysine decarboxylase, arginine dihydrolase, ornithine decarboxylase, β-galactosidase (ONPG), and catalase. The strain assimilates aesculin, amygdalin, ᴅ-cellobiose, ᴅ-glucose, ᴅ-galactose, ᴅ-lactose, ᴅ-maltose, ᴅ-melibiose, ᴅ-sorbitol, ᴅ-trehalose, ᴅ-xylose, ʟ-arabinose, ʟ-rhamnose, mannitol, malonate, glycerol, methyl α-ᴅ-glucopyranoside, and salicin. It is negative for acetoin production (Voges-Proskauer test), DNase activity, H₂S production, oxidase activity, phenylalanine deaminase activity, gelatinase, and urease. No utilization is observed for adonitol, dulcitol, erythritol, or inositol. Notably, the strain exhibits weak lipase activity and is positive for arginine dihydrolase, while ᴅ-arabitol metabolism is not detected. The major cellular fatty acids are C_16:0 (26.5%), summed feature 3 (C_16:1ω7c/C_16:1ω6c, 24.4%), summed feature 8 (C_18:1ω7c, 18.0%) and C_17:0 cyclo (9.5%). The genomic DNA G + C content is 53.49 mol%.

The type strain is 142053^T (= GDMCC 1.4845^T = JCM 37413^T), isolated from a clinical sputum sample at West China Hospital, Chengdu, China. The GenBank accession number for the 16S rRNA gene sequence is PV388010, and the draft genome sequence has been deposited under JBEFLV000000000 at DDBJ/EMBL/GenBank.

Description of Kluyvera chengduensis sp. nov.

Kluyvera chengduensis (cheng. du. en’sis N. L. fem. Adj. chengduensis, referring to Chengdu, China, where the type strain was recovered).

Cells of strain 142359^T are Gram-negative, facultatively anaerobic, motile, non-spore-forming rods, measuring 0.6–0.8 μm in width and 1.0–1.5 μm in length. After 12-h incubation on nutrient agar, colonies are round, smooth, convex, and white. Growth occurs at 8–42 °C (optimum, 28–37 °C), pH 5.0–9.0 (optimum, 6.0–8.0), and NaCl concentrations 0–5% (w/v) (optimum, 0–4%). Strain 142359^T is positive for indole production, nitrate reduction, citrate utilization, and the activities of lysine decarboxylase, arginine dihydrolase, ornithine decarboxylase, β-galactosidase (ONPG), and catalase. The strain assimilates aesculin, amygdalin, ᴅ-cellobiose, ᴅ-glucose, ᴅ-galactose, ᴅ-lactose, ᴅ-maltose, ᴅ-melibiose, ᴅ-sorbitol, ᴅ-trehalose, ᴅ-xylose, ʟ-arabinose, ʟ-rhamnose, mannitol, malonate, glycerol, methyl α-ᴅ-glucopyranoside, and salicin. It is negative for acetoin production (Voges–Proskauer), DNase activity, H₂S production, oxidase activity, phenylalanine deaminase, gelatinase, and urease. The strain does not utilize adonitol, dulcitol, erythritol, or inositol. Additionally, it displays weakly positive lipase activity and is positive for arginine dihydrolase. Strain 142359^T can further metabolize ᴅ-arabitol. The major cellular fatty acids are C_16:0 (24.9%), summed feature 8 (C_18:1ω7c, 22.7%), summed feature 3 (C_16:1ω7c/C_16:1ω6c, 20.0%), and C_17:0 cyclo (9.7%). The genomic DNA G + C content is 54.48 mol%.

The type strain is 142359^T (= GDMCC 1.5094^T = JCM 37737^T), isolated from a clinical sputum sample at West China Hospital, Chengdu, China. The GenBank accession number for the 16S rRNA gene sequence is PV388012, and the draft genome sequence has been deposited under JBIQOI000000000 at DDBJ/EMBL/GenBank.

Discussion

Strains 142053^T and 142359^T were preliminarily assigned to the Kluyvera species by MALDI-TOF. A maximum-likelihood tree of 265 Enterobacteriaceae core genes placed strains 142053^T and 142359^T within the Kluyvera clade, clearly separated from related genera. Phylogenomic analysis based on 1,668 Kluyvera core genes further resolved strain 142053^T and strain 142359^T into distinct branches, supporting their classification as two different species. Although the genus Kluyvera currently comprises only five validly published species, inclusion of five previously proposed but unvalidated genomospecies (1, 2, 3, 4/6, and 5) was necessary to provide a comprehensive phylogenomic framework for the genus. These genomospecies have been consistently recognized in previous taxonomic and genomic studies as distinct Kluyvera lineages, and therefore they were retained in the phylogenetic analysis (Liu et al., 2020; Rodríguez et al., 2021; Rodriguez and Gutkind, 2024). The analysis revealed that K. genomospecies 5 and K. huaxiensis clustered together and exhibited high isDDH and ANI values (Table 1; Figure 2), confirming that they represent the same species. Therefore, the entity previously designated as K. genomospecies 5 should be reclassified as K. huaxiensis. Notably, the genomic sequence named Kluyvera sp. STS39-E (accession no. GCA 040932435) likely represents a novel genomospecies within the genus Kluyvera. As for clinical significance, strain 142053^T was isolated from a patient with bilateral interstitial pulmonary fibrosis who did not exhibit any clinical manifestation of pneumonia, and hence it was considered as a colonizer. In contrast, strain 142359^T was recovered from a patient with pneumonia, which improved following treatment with piperacillin-tazobactam. However, as the same sputum sample also yielded Klebsiella pneumoniae, the clinical significance of strain 142359^T remains uncertain, and it may represent either a pathogen in a mixed infection or a colonizer in a patient with K. pneumoniae-associated pneumonia.

There are limitations of this study. First, only a single isolate was obtained for each of the two proposed species. Although both species are represented by at least two independent genomes, the phenotypic characterization based on a single isolate may not fully capture intraspecies variability. Second, phenotypic features of other Kluyvera species were derived from published studies, which may vary due to difference in experimental conditions. Third, polar lipids, quinones and whole-cell sugar compositions were not analyzed. As such, phenotypic traits used to distinguish the two proposed species from other Kluyvera species should be interpreted with caution.

Nevertheless, the core-genome phylogenetic analysis, together with chemotaxonomic data, clearly places strain 142053^T and 142359^T within the genus Kluyvera. The ANI and isDDH analyses conclusively demonstrate that these strains represent two novel species, for which the names Kluyvera huaxiensis sp. nov. and Kluyvera chengduensis sp. nov. are proposed.

Data availability statement

The 16S rRNA gene sequences for strains 142053^T are deposited in GenBank under accession no. PV388010. Additionally, strain 142359^T has 16S rRNA gene accession no. PV388012. Genome assemblies for these strains are recorded under the accession nos. JBEFLV000000000 and JBIQOI000000000, respectively.

Ethics statement

The study has been approved by the Ethical Committee of West China Hospital with waiver of informed consent as this study was to characterize bacterial strains.

Author contributions

LL: Methodology, Formal analysis, Writing – original draft, Investigation. YF: Investigation, Writing – review & editing, Data curation, Formal analysis. HW: Investigation, Writing – review & editing. ZZ: Resources, Writing – review & editing, Data curation, Project administration, Supervision, Funding acquisition, Conceptualization.

Funding

The author(s) declared that financial support was received for this work and/or its publication.The work was supported by grants from the National Key Research and Development Program of China (grant nos. 2023YFC2308800, 2022YFC2303900, and 2024YFE0106200).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1683098/full#supplementary-material

Footnotes

¹^https://tygs.dsmz.de/user_requests/new

²^https://github.com/tseemann/abricate

³^https://genomicepidemiology.org

⁴^http://www.eucast.org/

References

Akiki, M., Ali, R., Jamil, A., Slim, J., and Miller, R. (2023). Kluyvera ascorbata: an unusual cause of septic shock in a patient with urothelial cancer. Cureus 15:e51057. doi: 10.7759/cureus.51057,

PubMed Abstract | Crossref Full Text | Google Scholar

Alsuheili, A. Z., Alsubaie, M. A., Alnajjar, A. A., and Al-Dajani, N. (2024). Kluyvera Cryocrescens urinary tract infection in a neonate: a case report and literature review. Pediatr. Infect. Dis. J. 43, e461–e463. doi: 10.1097/inf.0000000000004507,

PubMed Abstract | Crossref Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021,

PubMed Abstract | Crossref Full Text | Google Scholar

Bayliss, S. C., Thorpe, H. A., Coyle, N. M., Sheppard, S. K., and Feil, E. J. (2019). PIRATE: a fast and scalable pangenomics toolbox for clustering diverged orthologues in bacteria. Gigascience 8, 1–9. doi: 10.1093/gigascience/giz119,

PubMed Abstract | Crossref Full Text | Google Scholar

Church, D. L., Cerutti, L., Gürtler, A., Griener, T., Zelazny, A., and Emler, S. (2020). Performance and application of 16S rRNA gene cycle sequencing for routine identification of Bacteria in the clinical microbiology laboratory. Clin. Microbiol. Rev. 33:e00053–19. doi: 10.1128/cmr.00053-19,

PubMed Abstract | Crossref Full Text | Google Scholar

CLSI (Ed.) (2025). Performance standards for antimicrobial susceptibility testing. 35th Edn. Wayne: CLSI supplement M100. Clinical and Laboratory Standards Institute.

Google Scholar

Endo, G., Kanai, S., Nishio, H., Hashimoto, H., Higurashi, Y., Nomura, Y., et al. (2024). Kluyvera Georgiana bacteremia due to acute cholangitis: a report of the first known case and a literature review. Intern. Med. 63, 2689–2693. doi: 10.2169/internalmedicine.3036-23,

PubMed Abstract | Crossref Full Text | Google Scholar

Farmer, J. J., Fanning, G. R., Huntley-Carter, G. P., Holmes, B., and Brenner, D. J. (1981). Kluyvera, a new (redefined) genus in the family Enterobacteriaceae: identification of Kluyvera ascorbata sp. nov. and Kluyvera cryocrescens sp. nov. in clinical specimens. J. Clin. Microbiol. 13, 919–933. doi: 10.1128/jcm.13.5.919-933.1981,

PubMed Abstract | Crossref Full Text | Google Scholar

Hernández-García, M., León-Sampedro, R., Pérez-Viso, B., Morosini, M. I., López-Fresneña, N., Díaz-Agero, C., et al. (2018). First report of an OXA-48- and CTX-M-213-producing Kluyvera species clone recovered from patients admitted in a University Hospital in Madrid, Spain. Antimicrob. Agents Chemother. 62:e01238–18. doi: 10.1128/aac.01238-18,

PubMed Abstract | Crossref Full Text | Google Scholar

Hu, Y., Feng, Y., Zhang, X., and Zong, Z. (2017). Acinetobacter defluvii sp. nov., recovered from hospital sewage. Int. J. Syst. Evol. Microbiol. 67, 1709–1713. doi: 10.1099/ijsem.0.001847,

PubMed Abstract | Crossref Full Text | Google Scholar

Humeniuk, C., Arlet, G., Gautier, V., Grimont, P., Labia, R., and Philippon, A. (2002). Beta-lactamases of Kluyvera ascorbata, probable progenitors of some plasmid-encoded CTX-M types. Antimicrob. Agents Chemother. 46, 3045–3049. doi: 10.1128/aac.46.9.3045-3049.2002,

PubMed Abstract | Crossref Full Text | Google Scholar

Inoue, T., Hadano, Y., Koga, H., Tamehiro, K., Sagara, S., Tokuda, Y., et al. (2023). Kluyvera intermedia bacteremia with septic shock: a case report. IDCases 32:e01765. doi: 10.1016/j.idcr.2023.e01765,

PubMed Abstract | Crossref Full Text | Google Scholar

Ito, R., Pacey, M. P., Mettus, R. T., Sluis-Cremer, N., and Doi, Y. (2018). Origin of the plasmid-mediated fosfomycin resistance gene fosA3. J. Antimicrob. Chemother. 73, 373–376. doi: 10.1093/jac/dkx389,

PubMed Abstract | Crossref Full Text | Google Scholar

Jain, C., Rodriguez, R. L., Phillippy, A. M., Konstantinidis, K. T., and Aluru, S. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9:5114. doi: 10.1038/s41467-018-07641-9,

PubMed Abstract | Crossref Full Text | Google Scholar

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Luo, L., Xiao, Z., Wang, G., Li, C., Zhang, Z., et al. (2019). Characterization of a carbapenem-resistant Kluyvera Cryocrescens isolate carrying Bla(ndm-1) from hospital sewage. Antibiotics (Basel, Switzerland) 8:149. doi: 10.3390/antibiotics8030149,

PubMed Abstract | Crossref Full Text | Google Scholar

Lin, J. C., Chen, C. H., Siu, L. K., and Chang, F. Y. (2002). Nosocomial outbreak of Kluyvera cryocrescens bacteremia. Infect. Control Hosp. Epidemiol. 23, 62–64. doi: 10.1086/503456,

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, L., Feng, Y., Wei, L., Qiao, F., and Zong, Z. (2020). Precise species identification and taxonomy update for the genus Kluyvera with reporting Kluyvera sichuanensis sp. nov. Front. Microbiol. 11:579306. doi: 10.3389/fmicb.2020.579306,

PubMed Abstract | Crossref Full Text | Google Scholar

Loiodice, M., Ribeiro, M., Peixe, L., and Novais, Â. (2023). Emergence of NDM-1 and KPC-3 carbapenemases in Kluyvera Cryocrescens: investigating genetic heterogeneity and acquisition routes of Bla(NDM-1) in Enterobacterales species in Portugal. J. Glob. Antimicrob. Resist. 34, 195–198. doi: 10.1016/j.jgar.2023.07.008,

PubMed Abstract | Crossref Full Text | Google Scholar

Madeira, F., Park, Y. M., Lee, J., Buso, N., Gur, T., Madhusoodanan, N., et al. (2019). The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 47, W636–w641. doi: 10.1093/nar/gkz268,

PubMed Abstract | Crossref Full Text | Google Scholar

Meier-Kolthoff, J. P., Auch, A. F., Klenk, H. P., and Goker, M. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform. 14:60. doi: 10.1186/1471-2105-14-60,

PubMed Abstract | Crossref Full Text | Google Scholar

Meier-Kolthoff, J. P., and Göker, M. (2019). TYGS is an automated high-throughput platform for state-of-the-art genome-based taxonomy. Nat. Commun. 10:2182. doi: 10.1038/s41467-019-10210-3,

PubMed Abstract | Crossref Full Text | Google Scholar

Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., von Haeseler, A., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015,

PubMed Abstract | Crossref Full Text | Google Scholar

Moreno, C., Romero, J., and Espejo, R. T. (2002). Polymorphism in repeated 16S rRNA genes is a common property of type strains and environmental isolates of the genus Vibrio. Microbiology (Reading) 148, 1233–1239. doi: 10.1099/00221287-148-4-1233,

PubMed Abstract | Crossref Full Text | Google Scholar

Müller, H. E., Brenner, D. J., Fanning, G. R., Grimont, P. A., and Kampfer, P. (1996). Emended description of Buttiauxella Agrestis with recognition of six new species of Buttiauxella and two new species of Kluyvera: Buttiauxella Ferragutiae sp. nov., Buttiauxella Gaviniae sp. nov., Buttiauxella Brennerae sp. nov., Buttiauxella Izardii sp. n. Int J Syst Bacteriol. 46, 50–63. doi: 10.1099/00207713-46-1-50

Crossref Full Text | Google Scholar

Mutoh, Y., Kobe, T., Hirano, T., Ichihara, T., Takenaka, H., Niinomi, T., et al. (2019). The first case of third-generation cephalosporins resistant Kluyvera Ascorbata biliary tract infection in Japan: a case report and review of the literature. Idcases 15:e00498. doi: 10.1016/j.idcr.2019.e00498,

PubMed Abstract | Crossref Full Text | Google Scholar

Pavan, M. E., Franco, R. J., Rodriguez, J. M., Gadaleta, P., Abbott, S. L., Janda, J. M., et al. (2005). Phylogenetic relationships of the genus Kluyvera: transfer of Enterobacter Intermedius Izard et al. 1980 to the genus Kluyvera as Kluyvera intermedia comb. nov. and reclassification of Kluyvera Cochleae as a later synonym of K. intermedia. Int. J. Syst. Evol. Microbiol. 55, 437–442. doi: 10.1099/ijs.0.63071-0,

PubMed Abstract | Crossref Full Text | Google Scholar

Qin, Q. L., Xie, B. B., Zhang, X. Y., Chen, X. L., Zhou, B. C., Zhou, J., et al. (2014). A proposed genus boundary for the prokaryotes based on genomic insights. J. Bacteriol. 196, 2210–2215. doi: 10.1128/jb.01688-14,

PubMed Abstract | Crossref Full Text | Google Scholar

Richter, M., and Rossello-Mora, R. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA 106:19126–19131. doi: 10.1073/pnas.0906412106,

PubMed Abstract | Crossref Full Text | Google Scholar

Rodriguez, M. M., Ghiglione, B., Power, P., Naas, T., and Gutkind, G. (2018). Proposing Kluyvera georgiana as the origin of the plasmid-mediated resistance gene fosA4. Antimicrob. Agents Chemother. 62:e00710–18. doi: 10.1128/aac.00710-18,

PubMed Abstract | Crossref Full Text | Google Scholar

Rodriguez, M. M., and Gutkind, G. (2024). Re-updating the taxonomy of Kluyvera genus for a better understanding of CTX-M β-lactamase origin. Microbiol. Spectr. 12:e0405423. doi: 10.1128/spectrum.04054-23,

PubMed Abstract | Crossref Full Text | Google Scholar

Rodríguez, M. M., Power, P., Naas, T., and Gutkind, G. (2021). Redefining the origin and evolution of chromosomally encoded Bla(CTX-M/KLU) in the context of a revised taxonomy of genus Kluyvera. Antimicrob. Agents Chemother. 65:e0242420. doi: 10.1128/aac.02424-20,

PubMed Abstract | Crossref Full Text | Google Scholar

Rodríguez, M. M., Power, P., Radice, M., Vay, C., Famiglietti, A., Galleni, M., et al. (2004). Chromosome-encoded CTX-M-3 from Kluyvera ascorbata: a possible origin of plasmid-borne CTX-M-1-derived cefotaximases. Antimicrob. Agents Chemother. 48, 4895–4897. doi: 10.1128/aac.48.12.4895-4897.2004,

PubMed Abstract | Crossref Full Text | Google Scholar

Sarria, J. C., Vidal, A. M., and Kimbrough, R. C. 3rd (2001). Infections caused by Kluyvera species in humans. Clin. Infect. Dis. 33, E69–E74. doi: 10.1086/322686,

PubMed Abstract | Crossref Full Text | Google Scholar

Sasser, M. (1990). "Identification of bacteria by gas chromatography of cellular fatty acids". MIDI technical note 101. Newark, DE: MIDI inc.

Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153,

PubMed Abstract | Crossref Full Text | Google Scholar

Strepis, N., Voor In 't Holt, A. F., Vos, M. C., Zandijk, W. H. A., Heikema, A. P., Hays, J. P., et al. (2021). Genetic analysis of mcr-1-carrying plasmids from gram-negative bacteria in a Dutch tertiary care hospital: evidence for intrapatient and interspecies transmission events. Front. Microbiol. 12:727435. doi: 10.3389/fmicb.2021.727435

Crossref Full Text | Google Scholar

Tiwari, A., Paakkanen, J., Österblad, M., Kirveskari, J., Hendriksen, R. S., and Heikinheimo, A. (2022). Wastewater surveillance detected Carbapenemase enzymes in clinically relevant gram-negative Bacteria in Helsinki, Finland; 2011-2012. Front. Microbiol. 13:887888. doi: 10.3389/fmicb.2022.887888,

PubMed Abstract | Crossref Full Text | Google Scholar

Ulloa-Clavijo, C., Suárez-Laurés, A., Viejo De la Cuadra, G., Galván, L., Martínez-Suárez, C., and Sánchez-Álvarez, E. (2023). Kluyvera ascorbata sepsis in a patient on hemodialysis. Nefrologia (Engl Ed) 43, 792–794. doi: 10.1016/j.nefroe.2021.09.013,

PubMed Abstract | Crossref Full Text | Google Scholar

Wong, V. K. (1987). Broviac catheter infection with Kluyvera cryocrescens: a case report. J. Clin. Microbiol. 25, 1115–1116. doi: 10.1128/jcm.25.6.1115-1116.1987,

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, L., Zhu, X., and Liu, Y. (2025). Community-acquired pneumonia caused by Kluyvera intermedia: a rare case report. Respir. Med. Case. Rep. 55:102194. doi: 10.1016/j.rmcr.2025.102194,

PubMed Abstract | Crossref Full Text | Google Scholar

Yoon, S. H., Ha, S. M., Lim, J., Kwon, S., and Chun, J. (2017). A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie Van Leeuwenhoek 110, 1281–1286. doi: 10.1007/s10482-017-0844-4,

PubMed Abstract | Crossref Full Text | Google Scholar

Yoshino, Y., Nakazawa, S., Otani, S., Sekizuka, E., and Ota, Y. (2016). Nosocomial bacteremia due to Kluyvera cryocrescens: case report and literature review. IDCases 4, 24–26. doi: 10.1016/j.idcr.2016.02.007,

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, K., Huang, Z., Lan, R., Morris, J. G. Jr., Xiao, Y., Fu, S., et al. (2023). Genomic characteristion of opportunistic pathogen Kluyvera reveals a novel CTX-M subgroup. Microorganisms 11, 1–13. doi: 10.3390/microorganisms11122836,

PubMed Abstract | Crossref Full Text | Google Scholar

Zou, S. H., Zhu, L. Y., Li, Y., and Zhang, F. G. (2019). A case of a persistent postoperative infection caused by multidrug-resistant Kluyvera ascorbata in the oral and maxillofacial region. Case Rep. Infect. Dis. 2019, 1–4. doi: 10.1155/2019/2180567,

PubMed Abstract | Crossref Full Text | Google Scholar