Phylogenomic Analyses of the Genus Pseudomonas Lead to the Rearrangement of Several Species and the Definition of New Genera

Simple Summary Pseudomonas represents a very important bacterial genus that inhabits many environments and plays either prejudicial or beneficial roles for higher hosts. However, there are many Pseudomonas species which are too divergent to the rest of the genus. This may interfere in the correct development of biological and ecological studies in which Pseudomonas are involved. Thus, we aimed to study the correct taxonomic placement of Pseudomonas species. Based on the study of their genomes and some evolutionary-based methodologies, we suggest the description of three new genera (Denitrificimonas, Parapseudomonas and Neopseudomonas) and many reclassifications of species previously included in Pseudomonas. Abstract Pseudomonas is a large and diverse genus broadly distributed in nature. Its species play relevant roles in the biology of earth and living beings. Because of its ubiquity, the number of new species is continuously increasing although its taxonomic organization remains quite difficult to unravel. Nowadays the use of genomics is routinely employed for the analysis of bacterial systematics. In this work, we aimed to investigate the classification of species of the genus Pseudomonas on the basis of the analyses of the type strains whose genomes are currently available. Based on these analyses, we propose the creation of three new genera (Denitrificimonas gen nov. comb. nov., Neopseudomonas gen nov. comb. nov. and Parapseudomonas gen nov. comb. nov) to encompass several species currently included within the genus Pseudomonas and the reclassification of several species of this genus in already described taxa.


Introduction
Pseudomonas is one of the most diverse and adaptable prokaryotic genera and their metabolic versatility allows their members to survive in many different environments [1,2]. Members of the genus Pseudomonas have been identified in human and animal related sources, plants, soil, water environments, psychrophilic environments, and other environmental niches and hosts [3][4][5]. Also, some species of this genus are known to play relevant roles in their hosts, such as P. aeruginosa, which causes human lung infections [6], species belonging to the P. fluorescens lineage, which are able to promote plants growth [7,8], or even some diverse species suggested to interfere with insects' biology [9,10], amongst many other cases. The multiplicity of environments where Pseudomonas grow and diversify has led to the broad evolution of its members, making it one of the most diverse bacterial genera. Pseudomonas' diversity, together with the fact that these bacteria are tremendously versatile and easy to grow under laboratory conditions, has led to the continuous discovery of new species of this genus [4]. However, the large number of species belonging to Pseudomonas and related taxa has made their taxonomic classification challenging, involving constant reshaping and reclassifications [11,12].
Moreover, the correct assignment of bacterial species is extremely important for making correct assumptions when carrying out microbiological studies, especially those with ecological relevance [13]. Assigning bacterial functions to wrongly classified taxa could cause future research to be based on false statements [14]; thus, the correct taxonomic allocation of members of this widely spread genus is highly desirable.
During the last two decades, the systematic classification of Pseudomonas has been based primarily on 16S rRNA gene sequence comparisons complemented by the construction of phylogenetic trees based on rpoB, rpoD, and gyrB gene sequences [11,15]. Recently, the use of genomic tools has facilitated taxonomic analysis of the genus Pseudomonas and the description of many novel Pseudomonas species, i.e., [12,16,17]. In these types of studies, the average nucleotide identity (ANI) values and digital DNA-DNA hybridization (dDDH) have both been employed for validating and confirming the taxonomic relatedness among similar and different species, which is always based on the thresholds suggested by several authors [18,19].
Currently, the use of genomics has become so commonplace that genome sequences are usually required for describing and validating novel taxa [19]. As a result, the use of genomics in taxonomy has become known as 'phylogenomics' and refers to the computation of phylogenetic trees based on genome-scale approaches. In addition to the overall genome related index (OGRI), which mainly comprises ANI and dDDH indices, several tools and algorithms have recently been developed for the purpose of carrying out phylogenomics [19,20], and have been used to organize and classify a diverse number of taxa through genome-based strategies [21][22][23][24]. In this sense, ANIb values over 95-96% and dDDH values over 70% indicate that a certain pair of bacterial strains represents the same species [19,20]. Also, some authors suggested ANIb < 75-76% for genera delimitation [25,26]. Despite this, there are no well-defined thresholds for this purpose. In addition, average amino acid identities (AAI) are becoming to be calculated in taxonomy for novel genera descriptions; as an example, Ma et al. [27] suggested AAI > 86% as the threshold for genera delimitation in the family Enterobacteriaceae. However, there is no consensus on AAI values for the delimitation of species or upper taxonomic levels. Thus, phylogenetics are crucial in these cases, since these help to decipher the evolutionary divergences among different clades.
In addition, Hesse et al. [28] provided the genome sequence of many type strains of the Pseudomonas genus and used them to conduct evolutionary studies. In fact, last year, Lalucat et al. [12] used many of these genomes for constructing the phylogeny of Pseudomonas, organizing clades into the following diverse groups: P. anguilliseptica, P. straminea, P. putida, P. syringae, P. lutea, P. asplenii, P. fluorescens, P. pertucinogena, P. aeruginosa, P. resinovorans, P. linyingensis, P. oleovorans, P. stutzeri, and P. oryzihabitans. Also, in a recent study including 494 Pseudomonas genomes, some non-type strains for which genomes have been published are suggested to be mis-classified and that should represent different species of Pseudomonas or even different genera [29]. However, although many species are distantly related to the core of the genus, there are no published reports suggesting their reclassification in other genera.
Thus, the present study applies phylogenomic approaches to revise the taxonomic organization of the genus Pseudomonas through the analysis of public genomes from type strains of species currently included into the genus Pseudomonas and other genera of the family Pseudomonadaceae. Furthermore, we analyzed the 16S rRNA and housekeeping genes sequences to confirm these rearrangements. The obtained results are in agreement with those of the genome-based phylogeny and OGRI analysis. Altogether, our analyses support the reclassification of several species into new taxa.
Based on the results of this study, we propose the creation of three novel genera to encompass several species currently included into the genus Pseudomonas and the reclassification of some Pseudomonas species into the genera Chryseomonas, Stenotrophomonas, and Xanthomonas.

Genome Sequence Data and Annotation
The List of Prokaryotic Names with Standing in Nomenclature (LPSN; https://lpsn. dsmz.de/, accessed on 31 August 2020), was used to search for validated type strains of the family Pseudomonadaceae, and, when available, their genomic sequences were downloaded from the NCBI database (https://www.ncbi.nlm.nih.gov/, accessed on 31 August 2020). In addition, and owing to their high similarity with some Pseudomonas strains, genomes of the genera Stenotrophomonas and Xanthomonas were also downloaded for taxonomic purposes. Similarly, the genomes of the type strains of other genera were downloaded for further use as phylogenetic tree outgroups.
All genome sequences were uploaded in batch mode to RAST (v2.0) [30] and annotated using the default settings.
Basic statistics for each genomic sequence were obtained using QUAST (v5.0.2) [31] and the graphics were performed in R using the ggplot2 package [32].

Phylogenetics and Phylogenomics
Trees based on the 16S rRNA gene and MLSA based on different housekeeping genes were created using MEGA (v7) [33] as previously described [16,34]. The 16S rRNA and housekeeping genes sequences from all the analyzed type strains were retrieved from RAST annotations when available, or from the NCBI database when the complete sequence of a gene was unavailable in the genome.
In order to assess the similarity of the nucleotide sequences of the 16S rRNA genes of Pseudomonas type strains, their sequences were uploaded to EzBioCloud [35], and percentage of similarity data for each strain was recorded.
Phylogenomic trees were constructed using the UBCG (v3) tool (default settings) [36] which created alignments (with MAFFT) and trees based on 92 housekeeping genes. Then, these trees were visualized and edited in The Interactive Tree Of Life (iTOL) tool (v5) [37].

Comparative Genomics and a Genome Wide Association Study (GWAS)
For the purpose of carrying out comparative analyses, genomes were annotated using prokka (v1.14.6) [44] and their annotations in General Feature Format (GFF) were used to perform pan-genome calculations and comparisons. To do this, we ran PPanGGOLiN (v1.1.96) using the default settings, except for the MMSeq2 identity threshold employed for clustering, which was set at 0.7 [45]. Then, PPanGGOLiN scripts were used to create a protein presence/absence matrix table in the csv format. Pan-genome plots were created with the roary_plots.py script (https://sanger-pathogens.github.io/Roary/, accessed on 1 February 2021). GWAS analyses between proteins of different sets of genomes belonging to different phylogenetic clades were performed using Scoary (v1.6.16) [46]; the input consisted in the gene presence/absence matrix, traits tables (0-1 binary code) that identified genomes with clades. The output tables were examined to identify proteins exclusively present or absent in certain phylogenetic clades. The COG family of these proteins were assigned using eggNOG-mapper (v5) [47].

Phylogenomics in Pseudomonadaceae
To obtain a more global picture of the taxonomic organization of the genus Pseudomonas, genomes of type strains from the family Pseudomonadaceae were downloaded, whose genomic data are summarized in Table S1.
The phylogenomic tree built on the basis on 92 housekeeping genes from the genomes of all type strains of the genus Pseudomonas and related genera is shown in Figure 1. This tree distinguishes some Pseudomonas lineages and clades that are differentiated from the big clade that clusters most Pseudomonas type strains, including that of P. aeruginosa, the type species of the genus ( Figure 1). Some of these more divergent type strains have a different genome size and G + C mol% content ( Figure 2). OGRIs values shows that they have lower ANIb values with those type strains included in the main Pseudomonas clade than some thresholds that have been suggested for genera delimitation (75-76%) [25,26] (Table S2). On the other hand, some pairs of type species display higher ANIb or dDDH values than those suggested for species differentiation, which are 95-96% in the case of ANIb and 70% in the case of dDDH.
Additionally, we found that the encoded proteins in the genomes of the most divergent clades or lineages differ substantially from those of the remaining Pseudomonas ( Figure 3). Out of the protein clusters generated in the pangenomic analysis, the P. pertucinogena clade had 316 protein families which were unique to this lineage, as they were absent in any other Pseudomonas species (Figure 3 and Figure S1). Likewise, P. caeni DSM 24390 T , P. hussainii JMC 19513 T , and the P. luteola clade had 1462, 1784, and 708 unique proteins, respectively. To know what metabolic categories differentiate these divergent clades from the remaining Pseudomonas, these proteins were assigned to Clusters of Orthologous Groups of proteins (COGs). The majority of them were either hypothetical proteins or proteins included within the categories "Cell wall/membrane/envelope biogenesis", "Amino acid transport and metabolism", "Signal transduction mechanisms" or "Transcription" ( Figure S1). On the contrary, as depicted in Figure 3, several protein families abundant in Pseudomonas belonging to the main clade were absent in the P. pertucinogena and P. luteola clades.  Phylogenomic tree of the Pseudomonadaceae type strains based on 92 concatenated housekeeping genes. The tree was built with UBCG, which uses MAFFT to create the multi-gene alignment and FastTree for computing the tree. Red bars represent GC (mol%) content for each genome. Blue circles represent bootstrap values, which indicate the number of individual phylogenetic trees of each of the 92 genes that support each branch (Gene Support Index). Bootstrap value of 92 means that the branch is supported by all UBCG phylogenies.

Figure 2.
Size and G + C (mol%) content of the genomes of the Pseudomoans type strains. Boxplots represents these genomic features of the genus (orange) and those that remain when exclude the genomes of the type strains of P. hussaini, P. caeni, P. geniculata, P. cissicola, and those of the P. luteola and P. pertucinogena clades (purple).

Figure 2.
Size and G + C (mol%) content of the genomes of the Pseudomoans type strains. Boxplots represents these genomic features of the genus (orange) and those that remain when exclude the genomes of the type strains of P. hussaini, P. caeni, P. geniculata, P. cissicola, and those of the P. luteola and P. pertucinogena clades (purple).

Taxonomic Status of Phylogenetically Distant Lineages and Clades
Considering all above-mentioned analyses, a case-by-case summary of the situation of those phylogenetically distant lineages and clades to the main Pseudomonas group is presented below. Based on that, we resolve a more appropriate taxonomic classification for each of them.

Taxonomic Status of Phylogenetically Distant Lineages and Clades
Considering all above-mentioned analyses, a case-by-case summary of the situation of those phylogenetically distant lineages and clades to the main Pseudomonas group is presented below. Based on that, we resolve a more appropriate taxonomic classification for each of them.
3.2.1. Pseudomonas geniculata P. geniculata ATCC 19374 T forms a lineage in the Pseudomonadaceae phylogenomic tree highly divergent from other Pseudomonas strains and those of the related genera. Thus, we compared its 16S rRNA gene sequence with those of the type strains held in the NCBI database, finding that it is more closely related to Stenotrophomonas species than to Pseudomonas ones. This was already pointed out previously: Anzai et al. [11] and Ramos et al. [48], who placed it within the genus Stenotrophomonas, although this reclassification has never been validly published.
Hence, we performed diverse approaches to settle the taxonomic placement of this type species. The 16S rRNA-based phylogeny placed P. geniculata ATCC 19374 T within the genus Stenotrophomonas, together with S. maltophilia NCTC 10257 T and S. pavanii LMG 25348 T (Figure 4c), sharing with them a 99.6 and 99.8% sequence similarity, respectively. These results are in concordance with the phylogeny based on concatenated 16S rRNA and gyrB genes ( Figure 4b) and the UBCG tree ( Figure 4a). The ANIb and dDDH values shared among P. geniculata ATCC 19374 T and Stenotrophomonas strains were below the threshold values used for species delimitation [18,19] (Tables S3 and S4). Concretely, values of ANIb and dDDH between P. geniculata ATCC 19374 T and S. maltophilia NCTC 10257 T are 92.37% and 48.5%, respectively, and those between P. geniculata ATCC 19374 T and S. pavanii LMG 25348 T are 42.3% and 90.75%. Altogether, the results of OGRI, phylogenetic, and phylogenomic analyses allow us to propose the reclassification of P. geniculata as Stenotrophomonas geniculata comb. nov. shared among P. geniculata ATCC 19374 T and Stenotrophomonas strains were below the threshold values used for species delimitation [18,19] (Tables S3 and S4). Concretely, values of ANIb and dDDH between P. geniculata ATCC 19374 T and S. maltophilia NCTC 10257 T are 92.37% and 48.5%, respectively, and those between P. geniculata ATCC 19374 T and S. pavanii LMG 25348 T are 42.3% and 90.75%. Altogether, the results of OGRI, phylogenetic, and phylogenomic analyses allow us to propose the reclassification of P. geniculata as Stenotrophomonas geniculata comb. nov. Hu et al. [49] already reported many phenotypic and chemotaxonomic similarities between P. cissicola CCUG 18839 T and Xanthomonas species. This was also already pointed out by Anzai et al. [11]. Later, Parkinson et al. [50] performed phylogenetic approaches which located P. cissicola in the Xanthomonas citri subsp. citri clade. Hu et al. [49] already reported many phenotypic and chemotaxonomic similarities between P. cissicola CCUG 18839 T and Xanthomonas species. This was also already pointed out by Anzai et al. [11]. Later, Parkinson et al. [50] performed phylogenetic approaches which located P. cissicola in the Xanthomonas citri subsp. citri clade.
The phylogenetic trees of this work were built UBCG ( Figure 4a) and those built with the sequences of both the 16S rRNA gene and the concatenated sequences of the housekeeping genes gyrB, rpoD, dnaK, and fyuA ( Figure 5) clustered P. cissicola CCUG 18839 T together with Xanthomonas citri LMG9322 T , with X. campestris, X, axonopodis, X. euvesicatoria, and X. perforans also closely related. P. cissicola CCUG 18839 T and Xanthomonas citri LMG9322 T share ANIb and dDDH values of 98.3% and 86.7%, respectively (Tables S4 and S5), which are above the thresholds for species differentiation [18,19]. These values are lower than 94.0% and 55.8%, respectively between P. cissicola CCUG 18839 T and any of the remaining Xanthomonas species (Tables S4 and S5). Therefore, the high similarity between the type strains of P. cissicola and Xanthomonas citri supports the reclassification of the species P. cissicola into the genus Xanthomonas as a later synonym of Xanthomonas citri. ria, and X. perforans also closely related. P. cissicola CCUG 18839 T and Xanthomonas citri LMG9322 T share ANIb and dDDH values of 98.3% and 86.7%, respectively (Tables S4 and  S5), which are above the thresholds for species differentiation [18,19]. These values are lower than 94.0% and 55.8%, respectively between P. cissicola CCUG 18839 T and any of the remaining Xanthomonas species (Tables S4 and S5). Therefore, the high similarity between the type strains of P. cissicola and Xanthomonas citri supports the reclassification of the species P. cissicola into the genus Xanthomonas as a later synonym of Xanthomonas citri.

Pseudomonas caeni
P. caeni DSM 24390 T was located on a separate branch from the remaining Pseudomonas species in all phylogenetic trees constructed in this study (Figures 1 and 6). The analysis of the ANIb values (lower than 73.4% in all cases) and AAI values (lower than 70.9%) also showed large differences between the genomes of P. caeni DSM 24390 T and those of the remaining members of the family Pseudomonadaceae (Tables S2 and S6 and Figure 7). Indeed, genome size and G + C mol% content of P. caeni DSM 24390 T are extremely divergent from other type strains of the genus Pseudomonas (Figure 2), as it has been recently reported [28]. Therefore, given the results obtained in the present work P. caeni should be transferred to a new genus for which we propose the name Denitrificimonas gen. nov. because of its denitrifying potential, being Denitrificimonas caeni gen. nov. comb. nov. its type species.

Pseudomonas caeni
P. caeni DSM 24390 T was located on a separate branch from the remaining Pseudomonas species in all phylogenetic trees constructed in this study (Figures 1 and 6). The analysis of the ANIb values (lower than 73.4% in all cases) and AAI values (lower than 70.9%) also showed large differences between the genomes of P. caeni DSM 24390 T and those of the remaining members of the family Pseudomonadaceae (Tables S2 and S6 and Figure 7). Indeed, genome size and G + C mol% content of P. caeni DSM 24390 T are extremely divergent from other type strains of the genus Pseudomonas (Figure 2), as it has been recently reported [28]. Therefore, given the results obtained in the present work P. caeni should be transferred to a new genus for which we propose the name Denitrificimonas gen. nov. because of its denitrifying potential, being Denitrificimonas caeni gen. nov. comb. nov. its type species.

P. hussainii
P. hussainii JMC 19513 T is well differentiated from its closest related taxa in all the constructed phylogenetic and phylogenomic trees of the family Pseudomonadaceae ( Figures  1 and 6). The ANIb and AAI values between P. hussaini JMC 19513 T and other Pseudomonas species are below 75.3% and 70%, respectively (Tables S2 and S6 and Figure 7). Therefore, we propose the reclassification of the species P. hussainii in the novel genus Parapseudomonas gen. nov. as Parapseudomonas hussainii gen. nov. comb. nov.

Pseudomonas luteola Clade
This clade encompasses the species P. luteola, P. asuensis, P. zeshuii, and P. duriflava. The phylogenetic analyses based on the 16S rRNA genes, gyrB, rpoB, and rpoD housekeeping genes and 92 housekeeping genes selected using the UBCG tool (Figures 1 and 6) clearly separate the type species of this clade from all the remaining genera included in the family Pseudomonadaceae. The ANIb and AAI values between each one of the type strains of this clade and those of the other Pseudomonadaceae type strains were lower than

Conclusions
This work highlights the relevance of the use of genomics in prokaryotic taxonomy. The phylogenetic and phylogenomic analyses performed in this study provide an overview of the genus Pseudomonas showing misclassified species or lineages. Phylogenetics and phylogenomics, together with OGRI values, outperform the resolution of 16S rRNA and housekeeping genes analyses to clarify the taxonomic organization of a broad genus such as Pseudomonas. These analyses allowed the reclassification of some previous Pseudomonas species into newly defined genera and into other already defined ones.
The description is as given for Denitificimonas caeni, which is the type species. The genus has been separated from Pseudomonas based on the physiology and phylogenetic analyses of genome and 16S rRNA gene sequences.   (Figures 1 and 6). The ANIb and AAI values between P. hussaini JMC 19513 T and other Pseudomonas species are below 75.3% and 70%, respectively (Tables S2 and S6 and Figure 7). Therefore, we propose the reclassification of the species P. hussainii in the novel genus Parapseudomonas gen. nov. as Parapseudomonas hussainii gen. nov. comb. nov.

Pseudomonas luteola Clade
This clade encompasses the species P. luteola, P. asuensis, P. zeshuii, and P. duriflava. The phylogenetic analyses based on the 16S rRNA genes, gyrB, rpoB, and rpoD housekeeping genes and 92 housekeeping genes selected using the UBCG tool (Figures 1 and 6) clearly separate the type species of this clade from all the remaining genera included in the family Pseudomonadaceae. The ANIb and AAI values between each one of the type strains of this clade and those of the other Pseudomonadaceae type strains were lower than 75% and 73% (Tables S2 and S6 and Figure 7) and the similarity values of 16S rRNA gene sequences were equal to or lower than 97%.
At the same time, the type species of this group (Figures 1 and 6) showed among them ANIb values ranging from 76% to 95%, except in the case of P. zeshuii KACC 15471 T and P. luteola NBRC 103146 T , which showed an ANIb value of 97.9% and a dDDH value of 92.0% (Tables S2 and S3). Therefore, as previously suggested by Lalucat et al. [12], P. zeshuii is a later synonym of P. luteola.
The species P. luteola was transferred to the genus Chryseomonas in year 1987 based on its DNA relatedness with Chryseomonas polytricha [51] and transferred again to the genus Pseudomonas in 1997 based on the 16S rRNA gene sequence homology with this genus [52]. However, all analyses performed in this study showed that the clade of P. luteola corresponds to a different genus from Pseudomonas. Thus, the initial reclassification of the species into the genus Chryseomonas, as C. luteola, should be maintained. Based on this, we propose to transfer to the genus Chryseomonas the species P. asuensis and P. duriflava as Chryseomonas asuensis comb. nov. and Chryseomonas duriflava comb. nov.

Pseudomonas pertucinogena Clade
The clade of P. pertucinogena currently includes 18 Pseudomonas species (see Figures 1 and 6 and Table S2). The phylogenetic trees based on genome sequences, 16S rRNA gene sequences and the sequences of concatenated housekeeping genes gyrB, rpoB, and rpoD located all species of the P. pertucinogena clade in a divergent taxonomic group within the family Pseudomonadaceae (Figure 6). The 16S rRNA gene similarities, ANIb and AAI values between any of the type strains of the clade and those of the remaining type strains of the family are lower than 95%, 75%, and 67%, respectively (Tables S2 and S6 and Figure 7), which suggests that the clade is substantially divergent from the Pseudomonas genus. This was also reported by Lalucat et al. [12] who suggested that the cluster of P. pertucinogena should be considered a different genus of the Pseudomonadaceae family based on phylogenetic approaches. All species from this group showed ANIb values among them lower than 95%, except in the case of P. abyssi and P. gallaeciensis which showed an ANIb value of 97,6 and a dDDH value of 80.1% (Table S3). These results agree with those of Lalucat et al. [12], who suggested the synonymy of the names P. abyssi and P. gallaeciensis. Since the name P. abyssi has priority, P. gallaeciensis is a later synonym of P. abyssi. Considering the results obtained in this work, we propose the creation of the novel genus Neopseudomonas gen. nov. for the remaining 17 species of the P. pertucinogena clade, with 17 new combinations.

Conclusions
This work highlights the relevance of the use of genomics in prokaryotic taxonomy. The phylogenetic and phylogenomic analyses performed in this study provide an overview of the genus Pseudomonas showing misclassified species or lineages. Phylogenetics and phylogenomics, together with OGRI values, outperform the resolution of 16S rRNA and housekeeping genes analyses to clarify the taxonomic organization of a broad genus such as Pseudomonas. These analyses allowed the reclassification of some previous Pseudomonas species into newly defined genera and into other already defined ones.
The description is as given for Denitificimonas caeni, which is the type species. The genus has been separated from Pseudomonas based on the physiology and phylogenetic analyses of genome and 16S rRNA gene sequences. The description is as given for Pseudomonas caeni [54] with the following modification. The genomic G + C content of the type strain is 48.3% and its genomic size is approximately 3.02 Mbp. The type strain is HY-14 = DSM 24390 = KCTC 22292 = CCTCC AB208156.
The description is as given for Parapseudomonas hussainii, which is the type species. The genus has been separated from Pseudomonas based on phylogenetic analyses of genome and 16S rRNA gene sequences. The description is as given for Chryseomonas asuensis [56] with the following modification. The genomic G + C is 53.6% and its genomic size is approximately 5. Gram negative, motile, and non-spore-forming bacteria. Cells are rod-shaped. Aerobic. Catalase and oxidase positive. The major fatty acids are those from Summed feature 8 (C 18:1 ω6c and/or C 18:1 ω7c) and the main respiratory quinone is Q9. The G + C content as calculated from genome sequences is approximately 60% and its genome size, 4.0 Mbp. The type species is Neopseudomonas pertucinogena comb. nov.  The description is as given for Pseudomonas aestusnigri [59] with the following modification. The genomic G + C is 60.4% and its genomic size is about 3.83 Mbp. The type strain is VGXO14 = CCUG 64165 = CECT 8317.  The description is as given for Pseudomonas salegens [71] with the following modification. The genomic G + C is 57.7% and its genomic size is approximately 3.80 Mbp. The type strain is GBPy5 =CECT 8338 = IBRC-M 10762. The description is as given for Pseudomonas salina [72] with the following modification. The genomic G + C is 57.5% and its genomic size is approximately 4.26 bp. The type strain is XCD-X85 = CGMCC 1.12482 = JCM 19469. The description is as given for Pseudomonas xinjiangensis [73] with the following modification. The genomic G + C is 60.7% and its genomic size is about 3.54 Mbp. The type strain is S3-3 = CCTCC AB 207151 = DSM 23391 = NRRL B-51270.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/biology10080782/s1, Figure S1: Heatmap of COG categories of unique proteins of the divergent Pseudomonas type strains, Table S1: Genomic features of the type strains used in this study, Table S2: ANIb values shared among Pseudomonas type strains that are proposed as novel genera and representative type strains belonging to different Pseudomonas and Pseudomonadaceae clades, Table S3: ANIb values shared among different Stenotrophomonas type strains and P. geniculata ATCC 19374 T . Table S4: Digital DNA-DNA Hybridization (dDDH) values of different type strains that are studied for its taxonomic relocation. Table S5: ANIb values shared among different Xanthomonas type strains and P. cissicola CCUG 18839 T . Table S6: AAIb values shared among Pseudomonas type strains that are proposed as novel genera and representative type strains belonging to different Pseudomonas and Pseudomonadaceae clades.