Proposal to classify Ralstonia solanacearum phylotype I strains as Ralstonia nicotianae sp. nov., and a genomic comparison between members of the genus Ralstonia

A Gram-negative, aerobic, rod-shaped, motile bacterium with multi-flagella, strain RST, was isolated from bacterial wilt of tobacco in Yuxi city of Yunnan province, China. The strain contains the major fatty acids of C16:0, summed feature 3 (C16:1ω7c and/or C16:1ω6c), and summed feature 8 (C18:1ω7c and/or C18:1ω6c). The polar lipid profile of strain RST consists of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, and unidentified aminophospholipid. Strain RST contains ubiquinones Q-7 and Q-8. 16S rRNA gene sequence (1,407 bp) analysis showed that strain RST is closely related to members of the genus Ralstonia and shares the highest sequence identities with R. pseudosolanacearum LMG 9673T (99.50%), R. syzygii subsp. indonesiensis LMG 27703T (99.50%), R. solanacearum LMG 2299T (99.28%), and R. syzygii subsp. celebesensis LMG 27706T (99.21%). The 16S rRNA gene sequence identities between strain RST and other members of the genus Ralstonia were below 98.00%. Genome sequencing yielded a genome size of 5.61 Mbp and a G + C content of 67.1 mol%. The genomic comparison showed average nucleotide identity (ANIb) values between strain RST and R. pseudosolanacearum LMG 9673T, R. solanacearum LMG 2299T, and R. syzygii subsp. indonesiensis UQRS 627T of 95.23, 89.43, and 91.41%, respectively, and the corresponding digital DNA–DNA hybridization (dDDH) values (yielded by formula 2) were 66.20, 44.80, and 47.50%, respectively. In addition, strains belonging to R. solanacearum phylotype I shared both ANIb and dDDH with strain RST above the species cut-off values of 96 and 70%, respectively. The ANIb and dDDH values between the genome sequences from 12 strains of R. solanacearum phylotype III (Current R. pseudosolanacearum) and those of strain RST were below the species cut-off values. Based on these data, we concluded that strains of phylotype I, including RST, represent a novel species of the genus Ralstonia, for which the name Ralstonia nicotianae sp. nov. is proposed. The type strain of Ralstonia nicotianae sp. nov. is RST (=GDMCC 1.3533T = JCM 35814T).

The pathogen responsible for plant bacterial wilt was first named Bacillus solanacearum (Smith, 1896) and then successively placed into the genera Bacterium, Phytomonas, Xanthomonas, Pseudomonas, and Burkholderia, before being transferred into the current genus, Ralstonia (Yabuuchi et al., 1995;Paudel et al., 2020). Although classified as R. solanacearum, distinct strains from diverse plants from various regions showed DNA-DNA hybridization values below the species cut-off value of 70%, thus R. solanacearum is considered a species complex (Genin and Denny, 2012). Based on genetic analysis of the 16S-23S internal transcribed spacer (ITS) region, endoglucanase 1 https://lpsn.dsmz.de/ (egl), and hrpB gene sequences, the R. solanacearum species complex was divided into four phylogenetically distinct groups, named phylotypes, which roughly reflect the strains' geographical origins: phylotype I strains are often sampled from Asia; phylotype II strains are isolated from the Americas, which are further subdivided into three subgroups (IIA, IIB, and IIC) (Poussier et al., 2000;Sharma et al., 2022); phylotype III are obtained from Africa; and phylotype IV are collected from south-east Asia (Prior and Fegan, 2005;Villa et al., 2005). Furthermore, Safni et al. (2014) classified the phylotype II strains as the species R. solanacearum, phylotypes I and III strains as the species R. pseudosolanacearum, and the phylotype IV strains as R. syzygii. The phylotype IV strains were initially proposed as Pseudomonas syzygii by Yabuuchi et al. (1995), transferred into the genus Ralstonia in 2004, and further sub-classified into three subspecies, R. syzygii subsp. syzygii, R. syzygii subsp. indonesiensis, and R. syzygii subsp. celebesensis (Safni et al., 2014). Recently, based on 16S rRNA gene analysis and genome sequence comparisons, Fluit et al. (2021) found that several well-known strains of phylotype I, including GMI1000, OE1-1, EP1, and CQPS-1, might represent an unnamed species within the genus Ralstonia.
In the present study, a novel pathogenic bacterium, RS T , belonging to Ralstonia was isolated from tobacco bacterial wilt. It was identified by comparing its physiological and chemical characteristics, its 16S rRNA gene, and its whole genome with those of closely related type strains. Consequently, the name Ralstonia nicotianae sp. nov. is proposed for phylotype I strains and the type strain is RS T (=GDMCC 1.3533 T = JCM 35814 T ).

Sampling site and bacterial culture
Strain RS T was isolated from bacterial wilt of flue-cured tobacco sampled in 2018 in Yuxi, Yunnan, China. The approximate geographical coordinates were 102°30′2.22″E, 24°14′22.5″N. The site is 1,630 m above sea level and has a central semi-humid subtropical cool winter plateau monsoon climate, with an average annual temperature of 15-16°C. We used the symptomatic leaves for pathogen isolation. First, the petiole part was surface-disinfected with 75% (v/v) ethanol for 5 min. Second, we removed tissues, including the epidermis layers, palisade, and spongy tissues, other than the leaf vein vascular bundle. Third, after removing the 2-3 cm part at the base, the remaining leaf vein vascular bundle was cut into several pieces and soaked in a sterile petri dish containing 10 ~ 20 mL of distilled water for 5-10 min. Finally, 100 μL of bacterial suspension with 5-6 times 10-fold dilution were spread onto a Tripheny tetrazolium chloride (TZC) agar plate. The dish was incubated under aerobic conditions at 28°C. After 36-48 h, typical circular colonies had formed, with broad white edges, strong fluidity, and a pink to light thin red liquid in the middle. The isolate was named RS T and streaked onto TZC agar at 28°C and purified three times (Figure 1). Strain RS T was stored in 20% (v/v) glycerol at −80°C for long-term preservation. It has been deposited at the Guangdong Microbial Culture Collection Center (GDMCC 1.3533 T ) and the Japan Collection of Microorganisms (JCM 35814 T

16S rRNA sequencing, and phylogenetic and genomic comparisons
The genomic DNA of strain RS T was extracted from a 2.0 mL bacterial culture using a bacterial DNA isolation kit (catalog number D3107; GBCBIO Technologies Inc., Guangzhou, China). The universal primer pair 27F and 1492R was used to amplify the 16S rRNA gene sequence (Clayton et al., 1995). The amplicon was ligated into the pMD19-T vector (Takara, Dalian, China) and sequenced following the Sanger method at Sangon Biotech Inc. (Shanghai, China). The sequences were assembled using Contig Express (Invitrogen, Waltham, MA, United States). The sequence similarity of the 16S rRNA gene from RS T was calculated using the EzBiocloud database 2 (Yoon et al., 2017). Closely related sequences of type species within the genus Ralstonia were collected from the GenBank database 3 (Sayers et al., 2021) and aligned using MUSCLE mode in the Unipro UGENE version 33.0 (Unipro, Novosibirsk, Russia) (Table 1). Trimmed sequences were used for sequence identity calculations employing the Clustal Omega multiple alignment package in the 2 https://www.ezbiocloud.net/ 3 www.ncbi.nlm.nih.gov/genbank EMBL-EBI database (Sievers and Higgins, 2021). Phylogenetic trees were constructed using the maximum-likelihood (ML), neighborjoining (NJ), and minimum-evolutions (ME) algorithms, with bootstrap values of 1,000 replications in MEGA version 11 software (Kumar et al., 2016).
For genome sequencing, the glycerol suspension was streaked onto an agar plate containing Casamino Acid-peptone-glucose (CPG) medium (Kelman, 1954). A single colony was inoculated into CPG broth and grown for 24 h at 28°C at 180 r/min, followed by centrifugation at 10,000 × g for 10 min and washing with distilled water twice (Lu et al., 2021). The total DNA of strain RS T was isolated using an EZNA bacterial DNA kit (Omega Bio-tek, Winooski, VT, United States). The genome of strain RS T was sequenced using the PacBio RS II platform (Menlo Park, CA, United States) and the Illumina HiSeq 4,000 platform (San Diego, CA, United States) at the Beijing Genomics Institute (BGI, Shenzhen, China). Four SMRT cell Zero-Mode Waveguide arrays for sequencing were used by the PacBio platform to generate the subreads set. PacBio subreads of length < 1 kb were removed (Utturkar et al., 2014). The program Pbdagcon was used for self-correction (Chin et al., 2013). The draft genomic unitigs, which are uncontested groups of fragments, were assembled using the Celera Assembler against a high quality corrected circular consensus sequence subread set (Berlin et al., 2015). To improve the accuracy of the genome sequences, GATK (Poplin et al., 2018) and SOAP tool packages (SOAP2, SOAPsnp, SOAPindel) (Li et al., 2008(Li et al., , 2009) were used to make single-base corrections. Finally, to trace the presence of any plasmid, the filtered Illumina reads were mapped using SOAP to the bacterial plasmid database. 4 Genetic comparisons between strain RS T and type strains of the genus Ralstonia were tested based on average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values using the webserver JSpecies WS version 3.9.6 (Richter et al., 2015) and the Genome-to-Genome Distance Calculator version 3.0 (Meier-Kolthoff et al., 2022), respectively. We conducted the phylogenomic analysis based on genome sequences of strain RS T and thirteen closely related strains in the genus Ralstonia using the Type (Strain) Genome Server (Meier-Kolthoff et al., 2022). Strain Cupriavidus necator N-1 T served as an outgroup. To confirm the systematic classification of all Ralstonia strains in the GenBank database that were closely related to RS T , ANI and dDDH values were tested between genome sequences of strain RS T and R. pseudosolanacearum LMG 9673 T and those from the strains mentioned above. For a quick calculation, ANI values were primarily tested using the FastANI calculator in the Genome Taxonomy Database (GTDB) (Jain et al., 2018;Parks et al., 2021). Strains with FastANI values between 93 and 97% were selected further to calculate BLAST-based ANI (ANIb), MUMmer-based ANI (ANIm), and dDDH values as described above.
For gene content analysis, the whole genome sequence of strain RS T and those of closely related species in the genus Ralstonia were obtained from the GenBank database (Sayers et al., 2021), genes were predicted using GeneMarkS version 4.28 (Besemer et al., 2001) and the genome comparisons were performed using OrthoVenn2, with a threshold e-value of 1e-5 and inflation of 1.5 (Xu et al., 2019).    The type three effectors (T3Es) of strain RS T were scanned using the Ralsto T3E database (Peeters et al., 2013;Sabbagh et al., 2019). To analyze the core effector repertoires of plant pathogenic Ralstonia species, we compared the T3Es in Ralstonia strains with the complete genome sequences from four phylotypes, including five strains of phylotype I, four strains of phylotype II, five strains of phylotype III, and three strains of phylotype IV. In addition, the secondary metabolites and antibiotic gene clusters were analyzed in silico using the online web server antiSMASH bacterial version 6.1.1 (Blin et al., 2021).
The protocol used to analyze the strains' physiological and chemical characteristics was as described previously (Lu et al., 2022a). Culture media, including CPG agar (Kelman, 1954), Kelman's TZC (García et al., 2019), Luria-Bertani (LB) agar (LA), modified selective medium South Africa (mSMSA) agar (Elphinstone et al., 1996), nutrient agar (NA), potato dextrose agar (PDA), and trypticase soy agar (TSA, BD Difco, Franklin Lakes, NJ, United States), were used to analyze the growth of the strains when incubated at 28°C for 5 days. Cells grown on TSA at 28°C for 48 h were subjected to morphological observation using transmission electron microscopy (HT7700; HITACHI, Tokyo, Japan). The motilities of the tested strains were detected on modified semi-solid motility media (SMM) containing 0.35% agar (Kelman and Hruschka, 1973). A Gram-staining kit (D008-1-1; Nanjing Jiancheng Bioengineering Institute, Nanjing, China) was used stain the strains, which was detected under a light microscope (Axio imager.Z2; Zeiss, Oberkochen, Germany). The oxygen requirement assay was conducted by streaking the test strain on TSA slope covered with mineral oil and incubating at 28°C for 5 days (Lu et al., 2022b).
To analyze the chemical sensitivity, enzymatic activity, carbon source, and substrate utilization of the tested strains, the Biolog GEN III Microplate system, API 20NE, and API ZYM kits (bioMérieux, Marcy-l'Étoile, France), respectively, were used according to the manufacturer's instructions. Oxidase activities of the strains were tested using 1% (w/v) tetramethyl-pphenylenediamine, with Escherichia coli K-12 strain MG 1655 as the negative control and Pseudomonas lijiangensis LJ2 T as the positive control. The production of bubbles detected catalase activity after adding a drop of 3% (v/v) H 2 O 2 (Lu et al., 2022a). Cells from the tested strains grown on TSA at 28°C for 48 h were harvested to analyze their fatty acid profile (Sasser, 1990). Polar lipids and isoprenoid quinones of three proposed type strains were extracted and analyzed as described previously (Minnikin and Abdolrahimzadeh, 1971;Collins and Jones, 1981  Maximum-likelihood phylogenetic tree based on 16S rRNA sequences of strain RS T and type species in the genus Ralstonia. The tree was generated using MEGA 11.0. Bootstrap values (≥50) are shown at branching points as percentages of 1,000 replicates. T indicates the type strain. Closed circles indicate that the corresponding nodes were also recovered in the trees constructed with neighbor-joining and maximum-parsimony. Strain RS T is indicated in bold. Cupriavidus necator N-1 T is used as an outgroup. The GenBank accession numbers are shown in parentheses. Bar, 0.005 substitutions per nucleotide position.

Results and discussion
3.1. 16S rRNA phylogeny 16S rRNA gene sequence analysis showed that strain RS T was closely related to the type species in the genus Ralstonia and shared the highest 16S rRNA sequence identity (99.50%) with those of R. pseudosolanacearum LMG 9673 T and R. syzygii subsp. indonesiensis LMG 27703 T , followed by R. solanacearum LMG 2299 T (99.28%) and R. syzygii subsp. celebesensis LMG 27706 T (99.21%). The 16S rRNA gene sequence identities between strain RS T and the other type strains of the genus Ralstonia were all below 98.00% ( Table 2). The ML phylogenetic tree based on 16S rRNA gene sequences revealed that strain RS T formed a coherent cluster with type species of R. pseudosolanacearum LMG 9673 T , R. solanacearum LMG 2299 T , and three subspecies, including R. syzygii subsp. indonesiensis LMG 27703 T , R. syzygii subsp. celebesensis LMG 27706 T , and R. syzygii subsp. syzygii LMG 10661 T (Figure 2). The NJ and MP trees also indicated that strain RS T was consistently separated from its closely related type species of R. pseudosolanacearum, R. solanacearum, and R. syzygii (Supplementary Figures 1, 2).

Genome features
Genome sequencing yielded 3,971,364 reads and 492 Mb of clean data, with a genome coverage depth of 87×. The overall G + C content of strain RS T was 67.10 mol%, which is in the range of 63.0-67.3 mol% (avg. 65.86 mol%) for Ralstonia species (Lu et al., 2022a). The assembled genome of strain RS T comprised 5,613,239 bp in two replicons: a circular chromosome of 3,618,724 bp and a megaplasmid of 1,994,515 bp. Strain RS T contained 5,105 predicted genes; 58 transfer RNA (tRNA) genes; four sets of 5S rRNA, 16S rRNA, and 23S rRNA genes; and 23 small RNA (sRNA) genes ( Figure 3).
Genome comparisons between strain RS T and type strains of the genus Ralstonia showed that R. pseudosolanacearum LMG 9673 T had the highest ANIb of 95.23% and ANIm of 96.14% in comparison with strain RS T , which was located in the species' border cut-off value of 95-96%, followed by R. syzygii subsp. syzygii LMG 10661 T (ANIb = 91.18% and ANIm = 93.00%), R. syzygii subsp. indonesiensis LMG 27703 T (ANIb = 91.41% and ANIm = 92.84%), and R. solanacearum K60-1 T (ANIb = 89.69% and ANIm = 91.64%) ( Table 2). The ANI values between the seven other type strains of the genus Ralstonia and strain RS T were below 90.00%. The dDDH values yielded by formula two were used for analysis (Meier-Kolthoff et al., 2022). Data revealed that R. pseudosolanacearum LMG 9673 T shared the highest dDDH value with strains RS T at 66.20%, which was below the 70% threshold value recommended for species description. The other 10 type strains had dDDH values with strain RS T below 66.00% ( Table 2).
The phylogenomic analysis showed that 14 Ralstonia strains were divided into two main phylogenomic branches. Strain RS T was positioned in the branch that contained three species, including R. pseudosolanacearum, R. syzygii, and R. solanacearum, which are the main causal agents of bacterial wilts in a wide range of plants. Strain RS T formed a coherent cluster with type species of R. pseudosolanacearum LMG 9673 T at a supported bootstrap value of 100% (Figure 4). The results allowed us to conclude that strain RS T represents a novel species in the genus Ralstonia. Interestingly, compared with strain R. solanacearum K60-1 T , R. syzygii strains were genetically closer to strain RS T and R. pseudosolanacearum LMG 9673 T .
To confirm their systematic classifications, all Ralstonia strains closely related to RS T whose genome sequences were available in the GenBank database were subjected to genome comparisons using the FastANI calculator in the Genome Taxonomy Database (GTDB) (Jain et al., 2018;Parks et al., 2021)  Localization of the principal genomic elements predicted in strain RS T . For each figure, from outer to inner: ring 1-genome size; ring 2-forward strand gene, colored according to the cluster of orthologous groups (COG) classification; ring 3-reverse strand gene, colored according to the cluster of orthologous groups (COG) classification; ring 4-forward strand ncRNA; ring 5-reverse strand ncRNA; ring 6-repeat elements; ring 7-G + C content; ring 8-GC-SKEW. Chromosome and megaplasmid of strain RS T are presented on the left and right, respectively. The phylogenomic tree showing the relationship between strain RS T and its closely related strains in the genus Ralstonia. The phylogenomic tree was generated using the Type (Strain) Genome Server (Meier-Kolthoff and Göker, 2019). The numbers above branches are Genome BLAST Distance Phylogeny approach (GBDP) pseudo-bootstrap support values >50% from 100 replicates. T indicates the type strain. Strain RS T is shown in bold. Cupriavidus necator N-1 T is used as an outgroup. Genomic assembly numbers are shown in parentheses. Bar, 0.02 changes per nucleotide position.  (Supplementary Tables 1, 2), supporting the view that these strains belong to the species of R. pseudosolanacearum. Based on the genome sequences of strains RS T , R. pseudosolanacearum LMG 9673 T , and their closely related strains using the Type (Strain) Genome Server (Meier-Kolthoff et al., 2022), a phylogenomic tree also supported the view that strain RS T and its closely associated strains of phylotype I formed a coherent cluster with five complete genome-sequenced R. pseudosolanacearum strains (Supplementary Figure 3). These results suggested that strain RS T isolated from bacterial wilt of tobacco and its closely related strains of phylotype I represent a novel species in the genus Ralstonia.

Pan-genome profile and accessory genes
Data revealed that 11 whole genome sequences from strains RS T and type species in the genus Ralstonia had a total of 55,346 proteincoding genes, with an average of 5031.45 genes for each species. The protein coding sequences were grouped into 7,159 clusters, containing 4,564 orthologous clusters and 2,595 single-copy gene clusters. All strains shared 2,657 core genome orthologs, which are mainly involved in translation, cell shape regulation, phenylacetate catabolic process, and DNA binding transcription factor activity. Strains from seven non-plant pathogenic type species within the genus Ralstonia shared 252 specific core orthologous clusters, and strains from four plant pathogenic bacteria, including RS T , R. pseudosolanacearum LMG 9673 T , R. solanacearum K60-1 T , and R. syzygii PSI 7 T , had 301 unique core orthologous gene clusters, which were absent in non-pathogenic Ralstonia. When we compared the four plant pathogenic bacteria, type strains formed 4,468 clusters, which included 2,304 core orthologs. In addition, strain RS T had 16 specific protein family clusters consisting of 47 proteins, which were not present in type strains of the other three plant pathogenic Ralstonia. Eight unique clusters in strain RS T were assigned to biological processes such as nucleobase-containing compounds, cellular aromatic and nitrogen compounds, macromolecules, cellular, primary, and heterocycle metabolic processes, and DNA transposition. Strain RS T shared 158, 160, and 100 orthologous clusters with its closely related type strains R. pseudosolanacearum LMG 9673 T , R. syzygii subsp. indonesiensis PSI 7 T , and R. solanacearum K60-1 T , respectively, suggesting that strain RS T is more closely related to R. pseudosolanacearum LMG 9673 T and R. syzygii subsp. indonesiensis PSI 7 T than to R. solanacearum K60-1 T ( Figure 5).

Core repertoire of T3Es
A total of 57 candidate T3Es, three pseudogenes (RipB, RipS5, and RipAT), and seven hypothetical effectors without assignment were found in strain RS T , among which 37 core T3Es were also found in five other closely related strains (Supplementary Table 3). A unique effector, RipAK, was present in strains of phylotype I, but absent in those belonging to phylotypes II, III, and IV. RipAK acts as a plant hypersensitive response suppressor by interacting with host pyruvate decarboxylases (PDCs), inhibiting their oligomerization and enzymatic activities (Sun et al., 2017;Wang et al., 2021). Two T3E candidates, RipBI and Hyp9, were conserved and unique among strains of phylotype II. The functional role of RipBI in Ralstonia is undefined; however, its homologous core effector XopX in Xanthomonas is a virulence factor and can both activate and suppress host immune responses in rice (Metz et al., 2005;Stork et al., 2015). Further study found that XopX induces immune responses by interacting with XopQ (Deb et al., 2020). RipBF is a phylotype IV-specific effector, which was only found in the tested R. syzygii strains; however, its functional roles remain to be determined. RipAP was absent from all R. syzygii strains, but was present in all tested strains of phylotypes I, II, and III. Although 35 candidate T3Es were present in all tested R. pseudosolanacearum strains, no specific effector was found among phylotype III strains, including LMG 9673 T , CMR15, RUN2279, RUN2474, and UW386.
Frontiers in Microbiology 10 frontiersin.org effector RipB acts as a major avirulence factor in Nicotiana species and suppresses the basal defense in susceptible hosts to promote pathogen infection (Nakano and Mukaihara, 2019;Cao et al., 2022). In addition, the nuclear-localized effectors RipAB and RipX are involved in plant immune response (Zheng et al., 2019;Sun et al., 2020). However, the functions of the other core effectors, including RipAJ, RipC1, RipAI, RipZ, RipAO, RipR, RipAM, and RipAN, still need to be determined.

Gene prediction involved in secondary metabolites
In silico secondary metabolite analysis predicted 12 gene clusters in the genome of strain RS T , referring to seven cluster types, including aryl polyene, furan, homoserine lactone (hserlactone), non-ribosomal peptide synthetase (NRPS), NRPS-like, ribosomally synthesized and post-translationally modified peptide product-like (RiPP-like), terpene, and type I polyketide synthase (T1PKS) (Supplementary Table 4). Among them, the taxonomically widespread aryl polyene, RiPP-like, and terpene biosynthetic gene clusters were also present in all tested Ralstonia strains (Nishida et al., 2005;Arnison et al., 2013;Cimermancic et al., 2014). Interestingly, more known secondary metabolite gene clusters were detected in the plant pathogenic Ralstonia species (10.87 clusters/strain) than in the non-pathogenic group (6.28 clusters/strain), which might suggest that more secondary metabolites remain to be uncovered in the latter group. Data revealed that the hserlactone and furan types were found in pathogenic Ralstonia strains. Both metabolite clusters shared less than 20% similarity to known gene clusters, suggesting that these gene clusters might encode novel metabolites and play significant roles in pathogenicity. The homoserine lactone gene cluster, which is usually involved in signal molecule production, participates in cell-to-cell communication processes. A recent study showed that the novel RasI/R quorum sensing system regulates phylotype I strain EP1 virulence. When the signal molecule N-(3-hydroxydodecanoyl)homoserine lactone synthetic gene rasI was knocked out, the mutant lost virulence toward eggplants (Yan et al., 2022). The metabolites of NRPS, NRPS-like, and T1PKS types were also only found in plant pathogenic Ralstonia, except for R. pickettii K-288 T , which only contained a T1PKS gene cluster. In silico metabolite analysis found that the beta-lactone-containing protease inhibitor (betalactone) and redox-cofactor types appeared to be conserved in non-pathogenic Ralstonia strains.
Siderophores are essential factors for both plant pathogens and beneficial microorganisms (Gu et al., 2020). Based on their chemical nature, siderophores are mainly classified into hydroxamate, catecholate, and carboxylate (Khan et al., 2018). Data revealed that strain RS T carried three putative siderophore biosynthetic clusters. Firstly, the chromosome of strain RS T encodes a micacocidin biosynthetic cluster (Kreutzer et al., 2011), which is conserved among phylotypes I, III, and IV strains, but absent in strains of phylotype II. However, a gene cluster accounting for siderophore yersiniabactin biosynthesis was only present in strains of phylotype II. Secondly, we detected a gene cluster that shared 18% similarity to the biosynthesis gene cluster of carboxylate siderophore staphyloferrin B on the megaplasmid of strain RS T (Bhatt and Denny, 2004). Thirdly, a specific NRPS and/or NRPS-like gene cluster was conserved among tested five strains of phylotype I, with a size range of 48.51-51.11 kb. It consisted of two core biosynthetic genes, seven additional biosynthetic genes, two transport-related genes, six regulatory genes, and 16 other genes on the megaplasmid of strain RS T . It shared no sequence similarity with known gene clusters; however, the MIBiG comparison showed that it had the highest similarity score of 0.65 with the catecholate siderophore Amphi-enterobactin 1 biosynthetic gene cluster in the marine organism Vibrio harveyi BAA-1116(McRose et al., 2018. However, the functional roles of this secondary metabolite biosynthetic cluster are uncharacterized. In addition, nine out of 10 tested strains from phylotypes I and III shared a specific ralsolamycin biosynthetic cluster (100% similarity), consistent with a previous report that R. solanacearum K60-1 T does not encode this specific gene cluster (Spraker et al., 2016). The diffusible lipopeptide ralsolamycin was found in a survey of secondary metabolites that induced chlamydospore formation of fungi when co-cultured with strain GMI1000 (Spraker et al., 2016). The biosynthesis locus of ralsolamycin comprises a non-ribosomal peptide synthetase and polyketide synthase hybrid gene (RmyA) and an NRPS gene (RmyB) (Spraker et al., 2016). Data showed that the PhcBSR quorum sensing (QS) system regulates ralsolamycin production . A recent study found that ralsolamycin induces the production of the fungal secondary metabolite bikaverin when Fusarium fujikuroi was exposed to ralsolamycin (Spraker et al., 2018), suggesting that this secondary metabolite can act as a communication chemical between bacteria and fungi.

Physiology and chemotaxonomy
In the API assays, strains of phylotype I, including RS T , GMI1000, and CQPS-1, shared more similarity to R. pseudosolanacearum LMG 9673 T and R. solanacearum LMG 2299 T for many phenotypic characteristics, e.g., all strains were positive for the reduction of potassium nitrate; hydrolysis of eculin; assimilation of D-glucose, potassium gluconate, malic acid, and trisodium citrate; and catalysis activity; and negative for indole production; the assimilation of Venn diagram of the distribution of shared candidate T3Es (core effectors) among strains of four phylotypes (https://bioinfogp.cnb. csic.es/tools/venny/index.html).
The polar lipids analysis showed that strain RS T contained phosphatidylethanolamine, diphosphatidylglycerol, phosphatidylglycerol, and an unidentified aminophospholipid (Supplementary Figure 4). Data from this study and our recent reports suggest that the major polar lipids of Ralstonia species are phosphatidylethanolamine and diphosphatidylglycerol (Lu et al., 2022a(Lu et al., , 2023. In addition, isoprenoid quinone analysis showed that strain RS T contained ubiquinone Q-7 (17.87%) and Q-8 (82.13%). Data from this study and our recent reports also indicated that Q-8 is the major respiratory quinone in Ralstonia species (Lu et al., 2022a(Lu et al., , 2023.

Conclusion
Based on the above polyphasic analysis, we concluded that the Ralstonia strains could be phylogenetically separated into two main clusters: non-pathogenic and plant-pathogenic bacteria. The current R. solanacearum phylotype I strain RS T and its close relatives represent a novel species in the genus Ralstonia, for which the name Ralstonia nicotianae sp. nov. is proposed.
The type strain, RS T (= GDMCC 1.3533 T = JCM 35814 T ), was isolated from a bacterial wilt of flue-cured tobacco sampled in 2018 from Yuxi, Yunnan, China. The genomic DNA G + C content of type strain RS T is 67.1 mol%. The GenBank accession number for the 16S rRNA gene sequence of the type strain is ON797093, and the genome sequence accession numbers are CP046674 and CP046675.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The GenBank accession numbers for the chromosome, megaplasmid and 16S rRNA gene sequences of strain RS T are CP046674.1, CP046675.1, and ON797093.1, respectively. The data are expressed as percentages of total fatty acids. TR, trace amount (less than 1.0% of the total); −, not detected. a Summed features are fatty acids that cannot be resolved reliably from another fatty acid using the chromatographic conditions chosen. The MIDI system groups these fatty acids together as one feature with a single percentage of the total. Summed feature 2 comprises C14 : 0 3-OH/C16 : 1 iso I. Summed feature 3 comprises C16 : 1 ω7c and/or C16 : 1 ω6c. Summed feature 8 comprises C18 : 1 ω7c and/or C18 : 1 ω6c.