Two novel Raoultella species associated with bleeding cankers of broadleaf hosts, Raoultella scottia sp. nov. and Raoultella lignicola sp. nov.

Abstract

Seventeen Gram-negative, facultatively anaerobic bacterial strains were isolated from bleeding cankers of various broadleaf hosts and oak rhizosphere soil in Great Britain. The strains were tentatively identified as belonging to the genus Raoultella based on 16S rRNA gene sequencing. Multilocus sequence analysis (MLSA), based on four protein-encoding genes (fusA, leuS, pyrG, and rpoB), separated the strains into three clusters within the Raoultella genus clade. The majority of strains clustered with the type strain of Raoultella terrigena, with the remaining strains divided into two clusters with no known type strain. Whole genome sequencing comparisons confirmed these two clusters of strains as belonging to two novel Raoultella species which can be differentiated phenotypically from their current closest phylogenetic relatives. Therefore, two novel species are proposed: Raoultella scottia sp. nov. (type strain = BAC 10a-01-01^T = LMG 33072^T = CCUG 77096^T) and Raoultella lignicola sp. nov. (type strain = TW_WC1a.1^T = LMG 33073^T = CCUG 77094^T).

1 Introduction

Raoultella was proposed as a novel genus in 2001 to contain mostly environmental “Klebsiella-like organisms” previously classified as Klebsiella terrigena, Klebsiella ornithinolytica and Klebsiella planticola (Drancourt et al., 2001). In 2014 a fourth species, Raoultella electrica isolated from biofilms of a microbial fuel cell, was proposed (Kimura et al., 2014) to bring the total of validly named species to four. There is an ongoing debate and some controversy around the phylogeny of Raoultella in relation to Klebsiella, with some taxonomists proposing that Raoultella should be reunified with Klebsiella due to the position of the genus clade within Klebsiella in a phylogenomic tree based on core protein sequences (Ma et al., 2021). The GTDB database also lists the Raoultella species under the genus Klebsiella (Parks et al., 2022), whereas another study suggests that Klebsiella is polyphyletic and could be divided into two clades, with Raoultella retained as a separate genus (Liu et al., 2020). Raoultella are ubiquitous in the environment and isolated from plants, water and soil although there are increasing reports of these species as clinical pathogens, especially R. planticola and R. ornithinolytica (Sękowska, 2017; Appel et al., 2021). Additionally, Raoultella species have been shown to successfully degrade fatty compounds, crude oil and biopolymers, and could play a potential role in protecting the natural environment through the degradation of toxic substances (Sękowska, 2019). For example, a novel Raoultella strain (isolated from heavy metal-polluted soil) demonstrated plant growth-promoting characteristics in experiments with pak choi, in addition to its cadmium-binding properties through ion exchange and chelation (Xu et al., 2019). Likewise, an encapsulated R. planticola strain applied to cotton seedlings under salt stress was shown to promote germination, biomass and photosynthetic pigment production (Wu et al., 2014). While these examples indicate the positive effects of Raoultella species on plant growth, they are predominantly associated with epiphytic lifestyles. Positive effects on plant protection by endophytic Raoultella species have been identified, such as the nematicidal activity of R. ornithinolytica against the pine wilt nematode (Shanmugam et al., 2018), although these are less frequently reported in literature.

In the present study, 17 strains tentatively identified as Raoultella species based on partial 16S rRNA gene sequencing, were isolated from bleeding cankers of various broadleaf hosts and oak rhizosphere soil. The bleeds were sampled as part of a survey of bleeding cankers and wetwood of broadleaf hosts in Great Britain.¹ Further molecular characterisation of the strains demonstrated that eight of these strains could not be assigned to an existing Raoultella species and following a polyphasic approach, were identified as belonging to two novel species for which the names Raoultella scottia sp. nov. and Raoultella lignicola sp. nov. are proposed.

2 Materials and methods

2.1 Isolation of bacteria

Bacterial were isolated from exudate, necrotic wood tissue or the rhizosphere soil. Swab samples were collected from exudate of the cankers, or from small bark panels (10 cm × 10 cm) of the bleeding lime trees. Swabs were rehydrated in sterile 1/4 Ringers, spread on Luria-Bertani (LB) agar (Oxoid) and eosin methylene blue (EMB) agar (Merck) and were incubated at 35°C under anaerobic conditions for 4 days. Wood chips were removed from the advancing lesion front of the lime panels with a sterile scalpel, plated on LB agar and incubated at 25°C for 7 days. Single colonies were obtained through re-streaking on LB agar and incubation at 25°C. Bacterial strains were isolated from the rhizosphere soil as previously described (Maddock et al., 2022). A list of strains used in the study along with their isolation source and location is presented in Supplementary Table 1.

2.2 Genotypic characterisation

DNA was extracted from all strains with an alkali lysis method (Niemann et al., 1997) followed by storage at −20°C. In order to assign a tentative identification, partial 16S rRNA gene sequencing was performed on all strains using the primer *pD (16F536), and the almost complete 16S rRNA gene was sequenced for the proposed type strains of the novel species using standard primers and PCR conditions (Coenye et al., 1999). Partial 16S rRNA gene sequences were compared against the EzBioCloud database to obtain a tentative identity (Yoon et al., 2017). Multilocus sequence analysis (MLSA) based on partial sequencing of four protein-encoding genes (fusA, leuS, pyrG and rpoB) was performed on the strains tentatively identified as Raoultella species using the published primers and conditions (Delétoile et al., 2009). Consensus sequences for both the 16S rRNA and protein-encoding genes were generated, aligned and trimmed in BioEdit v7.2.5 (Hall, 1999). The resulting lengths were: 16S rRNA—1,389 bp, fusA—633 bp, leuS—642 bp, pyrG—306 bp and rpoB—501 bp. Phylogenetic trees were constructed with the maximum likelihood method and 1,000 bootstrap replicates in PhyML 3.0 (Guindon et al., 2010) following Smart Model Selection (Lefort et al., 2017).

2.3 Genomic characterisation

Whole genome sequencing was performed on five representative strains (selected based on clusters observed from the MLSA results) from different hosts and locations: BAC 10a-01-01^T, Txe 2.2, WB_B2P2.3, TW_WC1a.1^T and BAC 2a-02-02. Sequencing was carried out by MicrobesNG (Birmingham, UK) using the Illumina HiSeq platform following DNA extraction by enzymatic cell lysis and purification with SPRI (Solid Phase Reversible Immobilization) beads. Adapters were trimmed from reads using Trimmomatic 0.30 with a sliding window quality cut-off of Q15 (Bolger et al., 2014) and de novo assembly was performed using SPAdes version 3.11.1 (Bankevich et al., 2012). The resulting contigs were annotated in the prokaryotic genome annotation pipeline (PGAP) (Tatusova et al., 2016) for GenBank submission, as well as in Prokka (Seemann, 2014).

Pairwise phylogenomic comparisons were carried out between the proposed novel species and existing Raoultella species, using genome blast distance phylogeny (GBDP) on the Type (Strain) Genome Server (TYGS) (Meier-Kolthoff and Göker, 2019). Accurate intergenomic distances were inferred under the algorithm “trimming” and distance formula d₅ (Meier-Kolthoff et al., 2013), with 100 distance replicates each. The resulting intergenomic distances were used to infer a balanced minimum evolution tree with branch support via FASTME 2.1.6.1 including SPR postprocessing (Lefort et al., 2015) and 100 pseudo-bootstrap replicates. The tree was rooted at the midpoint (Farris, 1972). DNA-DNA similarity comparisons were conducted with both average nucleotide identity (ANI), using FastANI (Jain et al., 2018), and in silico DNA-DNA hybridization (isDDH) with the Genome-to-Genome Distance Calculator (GGDC) (Meier-Kolthoff et al., 2022). Average amino acid identity (AAI) and pairwise percentage of conserved proteins (POCP) were calculated using the AAI-matrix calculator from the Kostas Lab (Rodriguez-R and Konstantinidis, 2016) and the script pocp.rb (Hoelzer, 2020) which follows the approach by Qin et al. (2014), respectively.

2.4 Genome features

Amino acid sequences from the whole genomes of BAC 10a-01-01^T, Txe 2.2, WB_B2P2.3 and TW_WC1a.1^T were examined to assess the pathogenic potential of their proteome. DIAMOND version v2.1.6 (Buchfink et al., 2021) was used to query the genomes against the Virulence Factor Database (VFDB) (Liu B. et al., 2022), accessed 13th October 2023. A query cut-off of 97% and a percentage identity equal to or greater than 50 were used to ensure high sequence alignment identification (Doonan et al., 2019). To assess the microbe-plant interaction potential of the two novel species, the amino acid annotations were queried against the “plant bacterial only interaction factors” (PIFAR-Pred) and “plant growth-promoting traits” (PGPT-Pred) databases using BlastP + HMMER Aligner/Mapper through the PLant-associated BActeria web resource (PLaBAse) (Martínez-García et al., 2016; Patz et al., 2021; Samad et al., 2022). Following open reading frame identification with orfipy (Singh and Wurtele, 2021), detection of potential type III secretion system effectors (T3SS) was performed using the Effectidor web server (Wagner et al., 2022).

2.5 Phenotypic and chemotaxonomic characterisation

Phenotypic and chemotaxonomic characterisation was carried out on all existing Raoultella species type strains and BAC 10a-01-01^T, Txe 2.2, WB_B2P2.3 and TW_WC1a.1^T from the present study. Bacterial strains were grown on LB agar for 24 h at 30°C to examine colony morphology, while the temperature range was determined following growth on LB agar at 4, 10, 25, 30, 37 and 41°C for 3 d. Salt and pH tolerance were tested in saline-free nutrient broth supplemented with NaCl (1–11%) in increments of 1% w/v, and in LB broth with the pH adjusted from 5 to 9 in increments of 1 as previously described (Brady et al., 2022). Catalase and oxidase activity were determined by bubble production in 3% v/v H₂O₂ and by staining with Kovács reagent (1% tetra-methyl-p-phenylenediamine dihydrochloride), respectively. Cell morphology, size and motility were determined using a light microscope and the CellSens software v 1.11 (Olympus Life Science, Tokyo, Japan) following growth in LB broth at 25°C for ∼3 h. The flagella arrangement was visualized by transmission electron microscopy (FEI Tecnai 12 120kV BioTwin Spirit TEM) following negative staining as previously published (Brady et al., 2022).

API 20E, API 50 CHB/E (bioMérieux) and GEN III GN/GP microplate (Biolog) assays were performed on the selected strains as per the manufacturer’s instructions. API 20E and API 50 CHB/E galleries were incubated at 37°C for 24 h and 30°C for 48 h, respectively, while GEN III microplates were incubated at 30°C for 18 h. Fatty acid methyl ester (FAME) profiles were determined by Fera Science Ltd. (York, UK) using the protocol based on the Sherlock Microbial Identification System Version 6.4 (MIDI Inc.). Strains were cultivated on TSA at 28°C for 24 h prior to fatty acid extraction and the results obtained were compared against the RTSBA6 6.21 library.

3 Results and discussion

3.1 Bacterial isolation

Strains were isolated over two years (2020–2021) from bleeding cankers of Liriodendron tulipifera (tulip tree), Tilia x europaea (common lime), Tilia x moltkei (Von Moltke’s lime), Fagus sylvatica (common beech) and Platanus x acerifolia (London plane) (Figure 1) across several counties of Great Britain. External symptoms included dark brown to black bleeding exudate, staining the outer bark while the underlying inner bark was necrotic and stained rust to brown in color. Additionally, several strains were isolated from the rhizosphere soil of both healthy and diseased oak suffering from Acute Oak Decline in Surrey, Great Britain.

FIGURE 1

3.2 Genotypic characterisation

Pairwise similarity values for the partial 16S rRNA gene sequences of the strains demonstrated highest similarity to R. terrigena (99.3–100%), and to Klebsiella, Kluyvera, Citrobacter and Lelliottia species (∼99.0%). In the maximum likelihood tree based on concatenated partial sequences of fusA, leuS, pyrG and rpoB genes, the strains tentatively identified as Raoultella species formed three clusters (Figure 2). A large proportion of strains, including those isolated from London plane, oak, beech, lime and rhizosphere soil, clustered with the type strain of Raoultella terrigena with 100% bootstrap support. The remaining strains were divided into two clusters, which were situated between the R. terrigena cluster and the other existing Raoultella species. The first contained seven strains from tulip and lime trees and clustered with 100% bootstrap support, basal to the R. terrigena cluster. The second contained a single strain isolated from lime (TW_WC1a.1^T) on its own branch, although without significant bootstrap support. Neither of these two clusters contained any reference strains, suggesting they belong to two novel species. The 16S rRNA gene tree (Supplementary Figure 1) demonstrated a polyphyletic topology for the genus Raoultella, with R. terrigena clustering between strains of one of the proposed novel species, and the remaining Raoultella species clustering at a lower level with the proposed single strain species, although there is little to no bootstrap support for these clusters. It should be noted that the 16S rRNA gene Sanger sequence for BAC 10a-01-01^T differa from that produced from the whole genome sequence by several base pairs. However, Raoultella species are known to have approximately eight copies of the 16S rRNA gene (Stoddard et al., 2015) and it has been previously demonstrated that sequence differences can exist between copies, especially within the Proteobacteria (Ibal et al., 2019) which may account for these nucleotide differences.

FIGURE 2

3.3 Genomic characterisation

Genome sequences were submitted to GenBank under the BioProject number PRJNA955701. Assembly accession numbers and genome metrics are listed in Supplementary Table 2. Prokka annotation of the sequences resulted in genomes of 4.39–5.93 Mbp, with a G + C content of 56.0 to 57.5 mol %. The phylogenomic tree presents a clear phylogeny for the genus Raoultella with current species and proposed novel species all contained in a well-supported clade (Figure 3). BAC 10a-01-01^T, WB_B2P2.3 and Txe 2.2 formed a single distinct cluster, with R. terrigena indicated as the current closest phylogenetic neighbor, while TW_WC1a.1^T clustered with the remaining Raoultella species at a lower level but also with 100% support. Results from the phylogenomic analyses agree with those of Liu et al. (2020), with species of the genus Klebsiella divided into two clades and Raoultella species constituting their own discrete clade. Furthermore, the topology of the phylogenomic tree is congruent with a previous phylogeny based on non-paralogous ribosomal MLST loci (Bray et al., 2022) where the genus Raoultella is shown to separate Klebsiella into two clades, the “pneumoniae-africana-variicola” clade and the “oxytoca-indica-grimontii-michiganensis” clade.

FIGURE 3

The ANI and isDDH values (Table 1) confirmed that BAC 10a-01-01^T, WB_B2P2.3 and Txe 2.2 belong to a single taxon with values exceeding 98.9 and 92.8%, respectively, which are well above the suggested species boundary cut-off values (Goris et al., 2007; Meier-Kolthoff et al., 2013). Additionally, both the ANI and isDDH values were higher between R. terrigena and the first proposed novel species reflecting their close phylogenetic relationship. TW_WC1a.1^T, comprising the second proposed novel species, displayed ANI values of 85.6–87.6% to the first proposed novel species and the existing Raoultella species.

Intra-genus AAI values of > 90% were observed between the proposed novel species and existing Raoultella species (Supplementary Table 3). Four representative Klebsiella species were included in the AAI and POCP analyses. The AAI values between Raoultella species, including the proposed novel species, and Klebsiella species were in the range of 87.4–89.9%. Several previous studies have used AAI to delineate genera with suggested boundary cut-off values, however, these do seem to differ between bacterial families (Walter et al., 2017; Wirth and Whitman, 2018; Nicholson et al., 2020; Zheng et al., 2020) and there is no single boundary that can be applied across all families and genera. Based on the values observed in the present study, the boundary cut-off value for the Enterobacteriaceae could be substantially higher than the 74–79% proposed for other families (Nicholson et al., 2020; Soutar and Stavrinides, 2022). The POCP values were less informative in terms of genus delineation as they ranged from 65 to 89.9% within Raoultella, 40.5 to 79.8% between Raoultella and Klebsiella species, and 46.7 to 84.0% within Klebsiella.

3.4 Genome features

Queries of the amino acid sequences against the VFDB revealed that several virulence genes are present in both novel species. The majority of these are associated with adherence, effector delivery systems, immune modulation and nutritional/metabolic factors. Genes for the type II, IV and VI secretion systems were found in the first novel Raoultella species (BAC 10a-01-01^T, Txe 2.2 and WB_B2P2.3) while only genes for the type II secretion system were found in the second novel Raoultella species (TW_WC1a.1^T). Results from the genome features queries are listed in Supplementary Table 4. The lack of type III secretion system genes for both species was confirmed by the results from Effectidor. The absence of a type III secretion system, and the relatively low numbers of virulence genes for biofilm formation, invasion and response regulation coupled with the low percentage of pectate cell wall degrading enzymes and lipopolysaccharides suggest that the two novel species do not have high pathogenic potential. Furthermore, the identified gene categories have a limited association to plant hosts. For example, many of the genes associated with immune modulation such as lpx genes that significantly contribute to lipopolysaccharide biosynthesis (Moosavian et al., 2020) are homologous to those experimentally identified in human pathogens such as Haemophilus influenzae. This could potentially be seen as a limitation of the VFDB, however, based on the presence of important plant pathogens such as P. syringae in the dataset it would seem more likely that the novel species analyzed in this work have limited ability to interact with their host. Likewise, the Hcp, clpV and dotU type VI secretion proteins identified here show the highest homology to those experimentally verified in Klebsiella pneumoniae. However, this is likely due to evolutionary lineage and does not limit the importance of these type VI secretion system genes which are known to play crucial roles in colonization, growth, competition, virulence, biofilm formation and environmental iron acquisition in plant pathogens such as Acidovorax citrulli and Xanthomonas phaseoli pv. manihotis (Montenegro Benavides et al., 2021; Fei et al., 2022). While the novel species identified in this work may not resemble plant pathogens, they still appear to have homologous genes required for colonization of the plant and/or lesion.

Relatively low percentages of plant growth-promoting traits also suggest that the two novel species do not play a significant role in improving plant nutrition and, as they were isolated from the lesion microbiome, this scenario is unlikely.

While the identified virulence genes limit the pathogenic potential of the two novel species, they do, however, highlight an interesting comparison. R. terrigena has previously been isolated from both the gut of ghost moth larvae (Thitarodes xiaojinensis), and from root and seed tissue of alpine bistort (Polygonum viviparum), which acts as the host plant during larvae maturation. R. terrigena was shown to degrade quercetin produced by the plant, enabling the growth of another moth larvae species which is normally inhibited by this secondary metabolite (Liu G. et al., 2022). With many of the virulence features identified in the two novel Raoultella species being related to immune modulation and nutritional/metabolic factors, the potential for a similar role for these bacteria and herbivorous insects that feed on weakened host plants should not be ignored. This model could be of importance when considering the correlation of Agrilus biguttatus and oak, a plant-pathogen-pest interaction in which a healthy host is able to limit larval development but a weakened host suffering from Acute Oak Decline is not (Brown et al., 2017).

3.5 Phenotypic and chemotaxonomic characterisation

Following incubation on LB agar at 30°C for 24 h, colonies were round, smooth with entire edges, 3–4 mm in diameter, slightly convex, translucent, non-pigmented and cream in color. Strains exhibited weak growth at 10°C but grew readily at 25–37°C, and only R. planticola and R. ornithinolytica grew well at 41°C. Strains from all Raoultella species grew in the pH range of 5–9 and could tolerate salt in supplemented broth up to 8% and weakly at 9–10%. All species were positive for catalase and negative for oxidase activity. Cells from both novel species are short rods (1.0–1.3 μm × 1.6–2.6 μm) and are found singularly, in pairs or rarely in chains, and are non-motile but possess fimbriae (Supplementary Figure 2).

The two proposed novel species could be differentiated from each other, and existing members of Raoultella, by their reactions to acetoin production, ornithine decarboxylase and urease, as well as the utilization of several carbon sources. Distinguishing characteristics for the genus are listed in Table 2, and a list of carbon sources utilized by all current members of Raoultella is presented in Supplementary Table 5. The fatty acid composition for the two proposed novel species were similar to those of the existing Raoultella species, with C_16:0, C_18:1 ω7c and summed feature 3 (C_16:1 ω7c and/or C_16:1 ω6c) comprising the major fatty acids. Complete fatty acid profiles are presented in Table 3.

4 Conclusion

Despite the role of several Raoultella species in causing human infections including urinary tract infections, sepsis and bacteraemia (Kaya et al., 2015; Seng et al., 2016; Mihu et al., 2023) there is a clear association of this genus with the natural environment with R. planticola and R. terrigena strains often isolated from botanical, aquatic and soil environments (Grimont and Grimont, 2015). The isolation of two novel Raoultella species associated with bleeding cankers of broadleaf hosts reinforces this association. Due to the low number of identified virulence factors in both novel species, and their sporadic isolation from bleeding cankers in various tree species, it is unlikely that they play a role in the symptom development. Instead, it is likely that they are part of the tree hosts’ microbiome or wider natural environment.

The future of the genus Raoultella is currently unclear as taxonomists fail to agree on whether the genus Klebsiella is polyphyletic, or monophyletic with the inclusion of Raoultella in the Klebsiella genus clade. There is little to differentiate these two genera based on phenotype or ecology, and their phylogeny varies depending on whether it is based on core genome sequences, core protein sequences or rMLST loci (Liu et al., 2020; Ma et al., 2021; Bray et al., 2022). As additional novel species are classified as belonging to either Raoultella or Klebsiella, the true taxonomy of these genera may become more apparent. The identification of two novel Raoultella species in the present study based on phylogenetic, genomic and phenotypic data, provides support for the genus to remain independent from Klebsiella for the time being. The names Raoultella scottia sp. nov. (type strain = BAC 10a-01-01^T = LMG 33072^T = CCUG 77096^T) and Raoultella lignicola sp. nov. (type strain = TW_WC1a.1^T = LMG 33073^T = CCUG 77094^T) are proposed for the strains isolated from bleeding cankers of broadleaf hosts. The description of Raoultella lignicola sp. nov. is currently based on a single strain and, as additional strains are isolated, the biochemical profile may change.

4.1 Emended description of the genus Raoultella (Drancourt et al., 2001)

Raoultella (Ra.oul.tel’la. M.L. dim. suff. tella; N.L. fem. dim. n. Raoultella, named after the French bacteriologist Didier Raoult, Université de la Méditerranée, Marseille, France)

The description is based on data from Drancourt et al. (2001), Kimura et al. (2014) and this study.

Gram-negative short capsulated rods with rounded ends (0.9–1.3 μm × 1.6–2.6 μm), facultatively anaerobic, oxidase-negative and catalase positive. Cells occur singularly, in pairs and sometimes in chains, non-motile but possess fimbriae. On LB agar, colonies are cream-colored, round, smooth with entire edges, 3–4 mm in diameter, slightly convex, and translucent but form large mucoid colonies on McConkey agar. Strains can grow at temperatures between 10 and 37°C, growth at 41°C is variable depending on the species. Strains can grow within the pH range 5–9 and in the presence of 0–8% (w/v) NaCl. Positive for β-galactosidase activity, lysine decarboxylase and citrate utilization, but negative for H₂S production and tryptophan deaminase and gelatinase activity. Arginine dihydrolase, ornithine decarboxylase and urease activity are variable along with acetoin and indole production. Nitrate is reduced to nitrite. Utilize glucose, mannose, fructose, galactose, D-glucose-6-phosphate and D-fructose-6-phosphate.

Major fatty acids include C_16:0, C_18:1 ω7c and summed feature 3 (C_16:1 ω7c and/or C_16:1 ω6c).

Frequently isolated from environmental sources including water, soil, plants and woody tissues of broadleaf trees. Several species are considered opportunistic pathogens and cause pneumonia and bacteraemia in oncologic and immune-compromised patients.

The type species is Raoultella planticola.

4.2 Description of Raoultella scottia sp. nov.

Raoultella scottia (scot’ti.a. N.L. fem. n. scottia, named after Lewis Scott for his contribution and support of Acute Oak Decline research.)

The description is as given for the genus with the following additional characteristics. Negative for arginine dihydrolase, urease activity and indole production, but positive for ornithine decarboxylase and acetoin production. Acid is produced from glycerol, L-arabinose, D-ribose, D-xylose, D-adonitol, D-galactose, D-glucose, D-fructose, D-mannose, L-sorbose, inositol, D-mannitol, D-sorbitol, methyl-αD-glucopyranoside, N-acetylglucosamine, arbutin, esculin ferric citrate, salicin, D-cellobiose, D-maltose, D-lactose, D-melibiose, D-saccharose, D-trehalose, D-raffinose, gentiobiose, L-fucose, D-arabitol, potassium gluconate and potassium 5-ketogluconate; while production from dulcitol and D-tagatose is variable.

The type strain is BAC 10a-01-01^T = LMG 33072^T = CCUG 77096^T, and was isolated from a bleeding canker of Liriodendron tulipifera (tulip tree) in Great Britain.

Accession numbers for the whole genome and 16S rRNA gene sequences are JARXNH000000000 and OQ813756, respectively.

4.3 Description of Raoultella lignicola sp. nov.

Raoultella lignicola (li.gni.co.la. L. masc./fem. n. suff. -cola, dweller; N.L. masc./fem. n. lignicola, wood inhabitant.)

The description is as given for the genus with the following additional characteristics. Negative for arginine dihydrolase, ornithine decarboxylase and urease activity, and indole and acetoin production. Acid is produced from glycerol, L-arabinose, D-ribose, D-xylose, D-adonitol, D-galactose, D-glucose, D-fructose, D-mannose, L-sorbose, inositol, D-mannitol, D-sorbitol, methyl-αD-glucopyranoside, N-acetylglucosamine, arbutin, esculin ferric citrate, salicin, D-cellobiose, D-maltose, D-lactose, D-melibiose, D-saccharose, D-trehalose, D-raffinose, gentiobiose, L-fucose, D-arabitol, potassium gluconate and potassium 5-ketogluconate.

The type strain is TW_WC1a.1^T = LMG 33073^T = CCUG 77094^T, and was isolated from a bleeding canker of Tilia x europaea (common lime) in Great Britain.

Accession numbers for the whole genome and 16S rRNA gene sequences are JARXNK000000000 and OQ813758, respectively.

References

• 1

  AppelT. M.Quijano-MartínezN.De La CadenaE.MojicaM. F.VillegasM. V. (2021). Microbiological and clinical aspects of Raoultella spp.Front. Public Health9:686789. 10.3389/fpubh.2021.686789

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S.et al (2012). SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing.J. Comput. Biol. J. Comput. Mol. Cell Biol.19455–477. 10.1089/cmb.2012.0021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data.Bioinformatics302114–2120. 10.1093/bioinformatics/btu170

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BradyC.AsselinJ. A.BeerS.BrurbergM. B.CramptonB.VenterS.et al (2022). Rahnella perminowiae sp. nov., Rahnella bonaserana sp. nov., Rahnella rivi sp. nov. and Rahnella ecdela sp. nov., isolated from diverse environmental sources, and emended description of the genus Rahnella.Int. J. Syst. Evol. Microbiol.72:005190. 10.1099/ijsem.0.005190

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BrayJ. E.CorreiaA.VargaM.JolleyK. A.MaidenM. C. J.RodriguesC. M. C. (2022). Ribosomal MLST nucleotide identity (rMLST-NI), a rapid bacterial species identification method: Application to Klebsiella and Raoultella genomic species validation.Microb. Genom.8:849. 10.1099/mgen.0.000849

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BrownN.JegerM.KirkS.WilliamsD.XuX.PautassoM.et al (2017). Acute oak decline and Agrilus biguttatus: The co-occurrence of stem bleeding and D-shaped emergence holes in Great Britain.Forests8:87. 10.3390/f8030087

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BuchfinkB.ReuterK.DrostH. G. (2021). Sensitive protein alignments at tree-of-life scale using DIAMOND.Nat. Methods18366–368. 10.1038/s41592-021-01101-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  CoenyeT.FalsenE.VancanneytM.HosteB.GovanJ. R. W.KerstersK.et al (1999). Classification of Alcaligenes faecalis-like isolates from the environment and human clinical samples as Ralstonia gilardii sp. nov.Int. J. Syst. Bacteriol.49405–413. 10.1099/00207713-49-2-405

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  DelétoileA.DecréD.CourantS.PassetV.AudoJ.GrimontP.et al (2009). Phylogeny and identification of Pantoea species and typing of Pantoea agglomerans strains by multilocus gene sequencing.J. Clin. Microbiol.47300–310. 10.1128/JCM.01916-08

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  DoonanJ.DenmanS.PachebatJ. A.McDonaldJ. E. (2019). Genomic analysis of bacteria in the Acute Oak decline pathobiome.Microb. Genom.5:15. 10.1099/mgen.0.000240

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  DrancourtM.BolletC.CartaA.RousselierP. (2001). Phylogenetic analyses of Klebsiella species delineate Klebsiella and Raoultella gen. nov., with description of Raoultella ornithinolytica comb. nov., Raoultella terrigena comb. nov. and Raoultella planticola comb. nov.Int. J. Syst. Evol. Microbiol.51925–932. 10.1099/00207713-51-3-925

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  FarrisJ. S. (1972). Estimating phylogenetic trees from distance matrices.Am. Nat.106645–668.

  • Google Scholar
  • View reference in article

• 13

  FeiN.JiW.YangL.YuC.QiaoP.YanJ.et al (2022). Hcp of the Type VI Secretion System (T6SS) in Acidovorax citrulli group II strain Aac5 Has a dual role as a core structural protein and an effector protein in colonization, growth ability, competition, biofilm formation, and ferric iron absorption.Int. J. Mol. Sci.23:9632. 10.3390/ijms23179632

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  GorisJ.KonstantinidisK. T.KlappenbachJ. A.CoenyeT.VandammeP.TiedjeJ. M. (2007). DNA-DNA hybridization values and their relationship to whole-genome sequence similarities.Int. J. Syst. Evol. Microbiol.5781–91. 10.1099/ijs.0.64483-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  GrimontP. A. D.GrimontF. (2015). “Klebsiella,” in Bergey’s manual of systematics of Archaea and bacteria, edsTrujilloM.DedyshS.DeVosP.HedlundB.KämpferP.RaineyF. A.et al (New York, NY: American Cancer Society), 10.1002/9781118960608.gbm01150

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  GuindonS.DufayardJ.-F. F.LefortV.AnisimovaM.HordijkW.GascuelO. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0.Syst. Biol.59307–321. 10.1093/sysbio/syq010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HallT. A. (1999). BioEdit: A user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT.Nucleic Acids Symp. Ser.4195–98.

  • Google Scholar
  • View reference in article

• 18

  HoelzerM. (2020). Calculation of percentage of conserved proteins. Available online at: https://github.com/hoelzer/pocp(accessed March 19, 2024).

  • Google Scholar
  • View reference in article

• 19

  IbalJ. C.PhamH. Q.ParkC. E.ShinJ. H. (2019). Information about variations in multiple copies of bacterial 16S rRNA genes may aid in species identification.PLoS One14:e0212090. 10.1371/journal.pone.0212090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  JainC.Rodriguez-RL. M.PhillippyA. M.KonstantinidisK. T.AluruS. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries.Nat. Commun.9:5114. 10.1038/s41467-018-07641-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  KayaS.BayramoğluG.SönmezM.KöksalI. (2015). Raoultella ornithinolytica causing fatal sepsis.Brazil. J. Infect. Dis.19230–231. 10.1016/j.bjid.2014.12.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  KimuraZ. I.ChungK. M.ItohH.HiraishiA.OkabeS. (2014). Raoultella electrica sp. nov., isolated from anodic biofilms of a glucose-fed microbial fuel cell.Int. J. Syst. Evol. Microbiol.641384–1388. 10.1099/ijs.0.058826-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  LefortV.DesperR.GascuelO. (2015). FastME 2.0: A comprehensive, accurate, and fast distance-based phylogeny inference program.Mol. Biol. Evol.322798–2800. 10.1093/molbev/msv150

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  LefortV.LonguevilleJ.-E.GascuelO. (2017). SMS: Smart model selection in PhyML.Mol. Biol. Evol.342422–2424. 10.1093/molbev/msx149

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  LiuB.ZhengD.ZhouS.ChenL.YangJ. (2022). VFDB 2022: A general classification scheme for bacterial virulence factors.Nucleic Acids Res.50D912–D917. 10.1093/nar/gkab1107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  LiuG.CaoL.HanR. (2022). Plant quercetin degradation by gut bacterium Raoultella terrigena of ghost moth Thitarodes xiaojinensis.Front. Microbiol.13:1079550. 10.3389/fmicb.2022.1079550

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  LiuL.FengY.WeiL.QiaoF.ZongZ. (2020). Precise species identification and taxonomy update for the Genus Kluyvera with reporting Kluyvera sichuanensis sp. nov.Front. Microbiol.11:579306. 10.3389/fmicb.2020.579306

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  MaY.WuX.LiS.TangL.ChenM.AnQ. (2021). Proposal for reunification of the genus Raoultella with the genus Klebsiella and reclassification of Raoultella electrica as Klebsiella electrica comb. nov.Res. Microbiol.172:103851. 10.1016/j.resmic.2021.103851

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  MaddockD.KileH.DenmanS.ArnoldD.BradyC. (2022). Description of three novel species of Scandinavium: Scandinavium hiltneri sp. nov., Scandinavium manionii sp. nov. and Scandinavium tedordense sp. nov., isolated from the oak rhizosphere and bleeding cankers of broadleaf hosts.Front. Microbiol.13:1011653. 10.3389/fmicb.2022.1011653

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  Martínez-GarcíaP. M.López-SolanillaE.RamosC.Rodríguez-PalenzuelaP. (2016). Prediction of bacterial associations with plants using a supervised machine-learning approach.Environ. Microbiol.184847–4861. 10.1111/1462-2920.13389

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  Meier-KolthoffJ. P.GökerM. (2019). TYGS is an automated high-throughput platform for state-of-the-art genome-based taxonomy.Nat. Commun.10:2182. 10.1038/s41467-019-10210-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  Meier-KolthoffJ. P.AuchA. F.KlenkH.-P.GökerM. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions.BMC Bioinform.14:60. 10.1186/1471-2105-14-60

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  Meier-KolthoffJ. P.CarbasseJ. S.Peinado-OlarteR. L.GökerM. (2022). TYGS and LPSN: A database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes.Nucleic Acids Res.50D801–D807. 10.1093/nar/gkab902

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  MihuA. G.SusanM. M.StrautiC. N.MotM. D.MuresanuH. D.BaltaC.et al (2023). First case of Raoultella planticola urinary tract infection reported in western Romania.Medicina59:506. 10.3390/medicina59030506

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  Montenegro BenavidesN. A.AlvarezB. A.Arrieta-OrtizM. L.Rodriguez-RL. M.BoteroD.TabimaJ. F.et al (2021). The type VI secretion system of Xanthomonas phaseoli pv. manihotis is involved in virulence and in vitro motility.BMC Microbiol.21:14. 10.1186/s12866-020-02066-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  MoosavianM.EmamN.PletzerD.SavariM. (2020). Rough-type and loss of the LPS due to lpx genes deletions are associated with colistin resistance in multidrug-resistant clinical Escherichia coli isolates not harbouring mcr genes.PLoS One15:e0233518. 10.1371/journal.pone.0233518

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  NicholsonA. C.GulvikC. A.WhitneyA. M.HumrighouseB. W.BellM. E.HolmesB.et al (2020). Division of the genus Chryseobacterium: Observation of discontinuities in amino acid identity values, a possible consequence of major extinction events, guides transfer of nine species to the genus epilithonimonas, eleven species to the genus Kaistella, and three species to the genus Halpernia gen. nov., with description of Kaistella daneshvariae sp. nov. and Epilithonimonas vandammei sp. nov. derived from clinical specimens.Int. J. Syst. Evol. Microbiol.704432–4450. 10.1099/ijsem.0.003935

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  NiemannS.PühlerA.TichyH. V.SimonR.SelbitschkaW. (1997). Evaluation of the resolving power of three different DNA fingerprinting methods to discriminate among isolates of a natural Rhizobium meliloti population.J. Appl. Microbiol.82477–484. 10.1046/j.1365-2672.1997.00141.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  ParksD. H.ChuvochinaM.RinkeC.MussigA. J.ChaumeilP.-A.HugenholtzP. (2022). GTDB: An ongoing census of bacterial and archaeal diversity through a phylogenetically consistent, rank normalized and complete genome-based taxonomy.Nucleic Acids Res.50D785–D794. 10.1093/nar/gkab776

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  PatzS.GautamA.BeckerM.RuppelS.Rodríguez-PalenzuelaP.HusonD. (2021). PLaBAse: A comprehensive web resource for analyzing the plant growth-promoting potential of plant-associated bacteria.bioRxiv[Preprint]. bioRxiv 2021.12.13.472471.

  • Google Scholar

• 41

  QinQ. L.XieB.Bin, ZhangX. Y.ChenX. L.ZhouB. C.et al (2014). A proposed genus boundary for the prokaryotes based on genomic insights.J. Bacteriol.1962210–2215. 10.1128/JB.01688-14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  Rodriguez-RL. M.KonstantinidisK. T. (2016). The enveomics collection: A toolbox for specialized analyses of microbial genomes and metagenomes.PeerJ4:e1900. 10.7287/peerj.preprints.1900v1

  • CrossRef
  • Google Scholar

• 43

  SamadA.NemanjaK.SaschaP.UlrikeL.BoykeB.CathrinS.et al (2022). Two new Rhizobiales species isolated from root nodules of common sainfoin (Onobrychis viciifolia) show different plant colonization strategies.Microbiol. Spectr.10e1099–e1022. 10.1128/spectrum.01099-22

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  SeemannT. (2014). Prokka: Rapid prokaryotic genome annotation.Bioinformatics302068–2069. 10.1093/bioinformatics/btu153

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  SękowskaA. (2017). Raoultella spp.–clinical significance, infections and susceptibility to antibiotics.Folia Microbiol.62221–227. 10.1007/s12223-016-0490-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  SękowskaA. (2019). The many faces of Raoultella spp.Postepy Hig. Med. Dosw.73713–720. 10.5604/01.3001.0013.6377

  • CrossRef
  • Google Scholar

• 47

  SengP.BoushabB. M.RomainF.GourietF.BruderN.MartinC.et al (2016). Emerging role of Raoultella ornithinolytica in human infections: A series of cases and review of the literature.Int. J. Infect. Dis.4565–71. 10.1016/j.ijid.2016.02.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  ShanmugamG.DubeyA.PonpandianL. N.RimS. O.SeoS. T.BaeH.et al (2018). Genomic insights into nematicidal activity of a bacterial endophyte, Raoultella ornithinolytica MG against pine wilt nematode.Plant Pathol. J.34250–255. 10.5423/PPJ.OA.10.2017.0217

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  SinghU.WurteleE. S. (2021). orfipy: A fast and flexible tool for extracting ORFs.Bioinformatics373019–3020. 10.1093/bioinformatics/btab090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  SoutarC. D.StavrinidesJ. (2022). Phylogenomic analysis of the Erwiniaceae supports reclassification of Kalamiella piersonii to Pantoea piersonii comb. nov. and Erwinia gerundensis to the new genus Duffyella gen. nov. as Duffyella gerundensis comb. nov.Mol. Genet. Genomics297213–225. 10.1007/s00438-021-01829-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  StoddardS. F.SmithB. J.HeinR.RollerB. R. K.SchmidtT. M. (2015). rrnDB: Improved tools for interpreting rRNA gene abundance in bacteria and archaea and a new foundation for future development.Nucleic Acids Res.43D593–D598. 10.1093/nar/gku1201

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  TatusovaT.DiCuccioM.BadretdinA.ChetverninV.NawrockiE. P.ZaslavskyL.et al (2016). NCBI prokaryotic genome annotation pipeline.Nucleic Acids Res.446614–6624. 10.1093/nar/gkw569

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  WagnerN.AvramO.Gold-BinshtokD.ZerahB.TeperD.PupkoT. (2022). Effectidor: An automated machine-learning-based web server for the prediction of type-III secretion system effectors.Bioinformatics382341–2343.

  • Google Scholar

• 54

  WalterJ. M.CoutinhoF. H.DutilhB. E.SwingsJ.ThompsonF. L.ThompsonC. C. (2017). Ecogenomics and taxonomy of Cyanobacteria phylum.Front. Microbiol.8:2132. 10.3389/fmicb.2017.02132

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  WirthJ. S.WhitmanW. B. (2018). Phylogenomic analyses of a clade within the roseobacter group suggest taxonomic reassignments of species of the genera Aestuariivita, Citreicella, Loktanella, Nautella, Pelagibaca, Ruegeria, Thalassobius, Thiobacimonas and Tropicibacter, and the proposal.Int. J. Syst. Evol. Microbiol.682393–2411. 10.1099/ijsem.0.002833

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  WuZ.GuoL.ZhaoY.LiC. (2014). Effect of free and encapsulated Raoultella Planticola Rs-2 on cotton growth promotion under salt stress.J. Plant Nutr.371187–1201. 10.1080/01904167.2014.881865

  • CrossRef
  • Google Scholar

• 57

  XuS.XingY.LiuS.HuangQ.ChenW. (2019). Role of novel bacterial Raoultella sp. strain X13 in plant growth promotion and cadmium bioremediation in soil.Appl. Microbiol. Biotechnol.1033887–3897. 10.1007/s00253-019-09700-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  YoonS. H.HaS. M.KwonS.LimJ.KimY.SeoH.et al (2017). Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies.Int. J. Syst. Evol. Microbiol.671613–1617. 10.1099/ijsem.0.001755

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  ZhengJ.WittouckS.SalvettiE.FranzC. M. A. P.HarrisH. M. B.MattarelliP.et al (2020). A taxonomic note on the genus Lactobacillus: Description of 23 novel genera, emended description of the genus Lactobacillus beijerinck 1901, and union of Lactobacillaceae and Leuconostocaceae.Int. J. Syst. Evol. Microbiol.702782–2858. 10.1099/ijsem.0.004107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar