Genomic characterization of a polyvalent hydrocarbonoclastic bacterium Pseudomonas sp. strain BUN14

Bioremediation offers a viable alternative for the reduction of contaminants from the environment, particularly petroleum and its recalcitrant derivatives. In this study, the ability of a strain of Pseudomonas BUN14 to degrade crude oil, pristane and dioxin compounds, and to produce biosurfactants, was investigated. BUN14 is a halotolerant strain isolated from polluted sediment recovered from the refinery harbor on the Bizerte coast, north Tunisia and capable of producing surfactants. The strain BUN14 was assembled into 22 contigs of 4,898,053 bp with a mean GC content of 62.4%. Whole genome phylogeny and comparative genome analyses showed that strain BUN14 could be affiliated with two validly described Pseudomonas Type Strains, P. kunmingensis DSM 25974T and P. chloritidismutans AW-1T. The current study, however, revealed that the two Type Strains are probably conspecific and, given the priority of the latter, we proposed that P. kunmingensis DSM 25974 is a heteronym of P. chloritidismutans AW-1T. Using GC-FID analysis, we determined that BUN14 was able to use a range of hydrocarbons (crude oil, pristane, dibenzofuran, dibenzothiophene, naphthalene) as a sole carbon source. Genome analysis of BUN14 revealed the presence of a large repertoire of proteins (154) related to xenobiotic biodegradation and metabolism. Thus, 44 proteins were linked to the pathways for complete degradation of benzoate and naphthalene. The annotation of conserved functional domains led to the detection of putative genes encoding enzymes of the rhamnolipid biosynthesis pathway. Overall, the polyvalent hydrocarbon degradation capacity of BUN14 makes it a promising candidate for application in the bioremediation of polluted saline environments.

www.nature.com/scientificreports/ l3-ketoadipate pathway (ortho-cleavage pathway) is the key route for the catabolism of a wide variety of aromatic compounds through the protocatechuate branch (pea genes) and the catechol branch (cat genes). The phylogeny of the genus Pseudomonas is complex, with some 219 validly published and correctly named species 11 . However, the current phylogeny is subject to considerable debate 12 and many of the named species may be considered as conspecific. For example, Pseudomonas kunmingensis HL22-2T and P. chloritidismutans AW-1T which are closely related to P. stutzeri may not deserve separate species status 13,14 .
To date, the genome sequences of eighteen strains affiliated to P. kunmingensis and one affiliated to P. chloritidismutans strain are available via the NCBI database 15 . While the genes and pathways associated with n-alkane degradation have been reported for P. chloritidismutans 16,17 , similar information for Pseudomonas kunmingensis is not available.
Pseudomonas sp. BUN14, isolated from hydrocarbon-polluted sediments from the Bizerte coast refinery harbour, north Tunisia, exhibits the capacity for hydrocarbon degradation. In this study, we combined genomic analysis and degradation experiments to demonstrate that Pseudomonas strain BUN14 has potential for application in the bioremediation of hydrocarbon contaminated marine environments.

Results and discussion
Isolation, phylogenetic assignment, and characterization of hydrocarbonoclastic bacterial strain BUN14. Strain BUN14 was isolated on ONR7a mineral medium supplemented with 1% crude oil as the sole carbon source. The comparison of the strain BUN14 16S rRNA gene sequence with those of validly described strains in the EzBioCloud 18,19 revealed that strain BUN14 showed highest similarity (99.23%) to that of P. kunmingensis HL22-2 T Phylogenetic analysis based on the 16S rRNA gene of Pseudomonas species placed P. kunmingensis HL22-2 T as the closest relative of BUN14 (Fig. S1). Growth of strain BUN14 on plates containing TSA (Trypticase soy agar) medium showed that the bacterium formed circular, opaque yellow colonies. Cells were rod-shaped and approximately 0.6 ± 0.1 μm in diameter and 1.8-2.0 μm in length (Fig. S2).
Optimization of surfactant activities. When grown on mineral medium ONR7a supplemented with crude oil and vegetable oil as sole carbon sources, BUN14 strain produced biosurfactants, as indicated by the oil spreading, drop collapse and emulsification tests (Fig. S2). Use of the Cetyl Tri Ammonium Bromide (CTAB)-Methylene blue agar method indicated that BUN14 produced anionic biosurfactant (Fig. S2). Optimal conditions for growth and biosurfactant production were assessed using response surface methodology (RSM) experiments. The analysis of variance for the fitted mathematical models shows that the regression sum of squares was statistically significant (P < 0.01) (Table S5). Response surface plots, showing optimal biosurfactant production conditions, are shown in Fig. 1. Statistical analyses of central composite design (CCD) experiments demonstrated that the percentage of oil substrate, NaCl concentration, inoculum size and incubation time all affected the biosurfactant production. Based on desirability function, the optimal BUN14 growth (OD 600nm , 0.54), biosurfactant production (OD 625nm , 2.58) and activities (emulsion index E24, 40.09% and oil-displacement ODA, 20.40 cm 2 ), were obtained after approx. 7 days of cultivation crude oil, 2.50% NaCl concentration inoculum size source (Fig. S3). Accordingly, as for previous studies 20, 21 substrate, NaCl concentration, inoculum size and incubation time significantly affected biosurfactant production (0.01 < p < 0.05).
The basic structure of the isolated biosurfactant was evaluated by Fourier Transform InfraRed (FT-IR) spectrometry and compared to a reference Pseudomonas aeruginosa (Fig. 2). The peaks appearing at 3278.56 (Fig. 2a) and 3270.42 cm −1 (Fig. 2b) denoted the presence of -OH stretching (free hydroxyl groups of rhamnose rings) of hydroxyl group. The adsorption peaks at 2924.64 (Fig. 2b) and 3000.0 cm −1 (Fig. 2a) indicated the presence of terminal methyl group of aliphatic stretching bands (CH2, CH3). The absorption peaks at 1097.35 (Fig. 2a) and 1049.1 cm −1 (Fig. 2b) confirmed the presence of C-O-C vibrations (rhamnose ring). The area between 1495.92 and 1150.02 cm −1 (Fig. 2a) represented C-H and OH deformation vibrations typical for carbohydrates, as in the rhamnose units of the biosurfactant. The apparent similarity of the main functional groups determined by FTIR spectrometry of the commercial rhamnolipid (R90) from Pseudomonas aeruginosa (AGAE Technologies, Corvallis, OR, USA) and the isolated biomolecule from strain BUN14 (Fig. 2), suggested that the BUN14 biosurfactant product was composed on rhamnose rings with long hydrocarbon chains.
Members of the genus Pseudomonas are known for biosynthesis of rhamnolipids with biosurfactant properties 20,22 . Hydrocarbon degradation. Strain BUN14 could grow on and utilize various hydrocarbons, including pyrene, naphthalene, phenanthrene, carbazole, dibenzofuran, dibenzothiohene, biphenyl, pristane, fluoranthene, crude oil, octadecane and tetradecane as sole carbon and energy sources (Fig. S4a-c). These findings are in agreement with previous reports demonstrating that many Pseudomonas species are capable of utilizing hydrocarbons as carbon and energy sources [23][24][25] . The kinetics of hydrocarbon biodegradation in liquid media was determined using GC-FID analysis. The total degradation of petroleum TERHCs by strain BUN14 was estimated at 90% after 21 days with rapid degradation of almost all alkanes C 12 -C 36 (Fig. S4d,e). Of the more recalcitrant aliphatic compounds, 22%, 40%, 30% and 14% of pristane, naphthalene, DBT and DBF, respectively, were degraded by strain BUN14 (Fig. S4d) over a period of 21 days. The ability of various bacteria to degrade naphthalene, including Sphingomonas, Pseudomonas and Acidovorax, has been reported 24, 26-29 . General features of the draft genome of strain BUN14. The (Table S1). Evaluation of the predicted proteins using BUSCO showed that 99.8% of gammaproteobacteria conserved single copy orthologs were present in the strain BUN14 draft genome with 0 duplication (Table 1). Functional annotation of the predicted proteome showed that ~ 87% of the proteins could be assigned to EggNOG orthologous genes (OGs). The overall distribution revealed that the category metabolism was overrepresented, comprising ~ 35% of the OGs compared to cellular processes and signalling, and information storage and processing, which comprise ~ 23% and 17%, respectively (Table S1).   S8). Two important species discriminatory genomic metrics, ANI and iDDH, provided further support for the clustering of these strains. Interestingly, the two type strains (DSM 25974T and AW-1), strains BUN14 and CCUG 36651 share ANI and iDDH values above 96 and 70%, respectively, suggesting that they belong to the same species [30][31][32][33] . The ANI range of 86.74-86.88% between the four strains and P. stutzeri ATCC 17588 T and the distinct clustering of the latter in the presented phylogeny provide strong evidence for the separation of the above www.nature.com/scientificreports/ four from the P. stutzeri group. Mehboob et al. 17 reported similar ANI values to the current data for AW-1 and strains P. stutzeri except P. stutzeri CCUG 29243 which share ANI value of 97% with AW-1. However, our work included only the type strain of P. stutzeri. P. chloritidismutans AW-1. The latter was first described as a novel species in 2002 based on DDH, physiological and biochemical data 14 and later proposed as conspecific with P. stutzeri, specifically P. stutzeri genomovar 3, based on multigenic phylogeny and by disputing the validity of the phenotypic features that separate the two taxa 34 . Subsequent studies 17 and this work clearly demonstrate that P. chloritidismutans AW-1 and P. stutzeri ATCC 17588 T are distinct species 34 . P. kunmingensis DSM 25974 T35 was proposed as a distinct species based on a combination of 16S rRNA gene analysis, phenotypic characterisation and DDH values 13 . However, the analysis omitted P. chloritidismutans AW-1 from the comparison. For instance, the low 16S rRNA gene similarity value reported between DSM 25974 T and closely related taxa was only possible without AW-1. The full-length 16S rRNA genes of AW-1 and DSM 25974 T (data not shown), extracted from their respective genomes, share a similarity of 99.74%. Moreover, the phylogenetic and the genomic metrics presented here indicate that DSM 25974 T and AW-1 share ANI and DDH values that are below the threshold for species separation [36][37][38] . Considering the priority of P. chloritidismutans AW-1 T and the multiple lines of phylogenetic and genomic evidence, we suggest that P. kunmingensis DSM 25974 T is a heterotypic synonym of P. chloritidismutans AW-1 T . Similarly, the strains BUN14 and CCUG 36651, and perhaps CCUG 29243 17 , may be more appropriately affiliated with P. chloritidismutans.

Genomic determinant of hydrocarbon degradation in strain BUN14 and comparisons with closely related Pseudomonas spp.
To identify the proteins potentially implicated in hydrocarbon degradation, the proteome of strain BUN14 was annotated using BlastKOALA 37 . Of the 4551 proteins, 2564 (56.34%) could be assigned to 2062 KEGG Orthologues (KOs) including 154 linked to xenobiotic biodegradation and metabolism (Fig. S9). Of the 154 proteins, 44 have been associated with the pathways for complete degradation of benzoate and naphthalene (Table S2). Gene product prediction using FGENESB revealed that 39 of these proteins are encoded in four genomic loci, comprising five putative operons, while the genes of five proteins appear to be transcribed as independent transcription units. The strain BUN14 genome contained six genes (dmpP, dmpO, dmpN, dmpM, dmpL and dmpK) encoding enzymes that catalyse the conversion of benzene to catechol (KEGG module: M00548; Fig. 4). The genome also contained the genes benA-xylX, benB-xylY, benC-xylZ and benD-xylL, which encode the enzymes for degradation of benzoate to catechol via catechol/methylbenzoate (KEGG module: M00551) (Fig. 4). Two distinct genomic loci contained bphH, bphI, bphJ, dmpB, dmpC, dmpD, dmpH, mhpE, mhpF and praC, encoding enzymes for the complete pathway for degradation of catechol (KEGG module: M00569, meta-cleavage) to yield Acetyl-CoA (Fig. 4), and the catB, catC, catA, and pcaD genes that encode enzymes for conversion of catechol (KEGG module: M00568; ortho-cleavage) to 3-oxoadipate (Fig. 4). The complete degradation pathway of benzoates and  Although partial degradation of dioxin was detected in culture experiments and the analysis of genome showed only some of genes of the dioxin degradation pathway; praC, dmpH bphI, bphH, bphJ, mhpE and mhpF (Table S2), associated with the conversion of 2-hydroxymuconate to acetyl-CoA. However, dxnA/dbfA, dbfB and dxnB which encode enzymes involved in the conversion of dibenzo-p-dioxin to 2-hydroxymuconate have not been identified. Genomic analysis of Pseudomonas putida strain B6-2 revealed genes for the complete biphenyl degradation pathway (represented by the gene cluster bphA-D and pbhH-K) 40 . Like strain BUN14, the genome sequence of Sphingomonas. wittichii contains only putative genes that code for enzymes in the initial part of the dioxin degradative pathway 41 . www.nature.com/scientificreports/ The strain BUN14 proteome included the nine proteins (nahD, nahE, nahC, nahF, nahB, nahAd, nahAc, nahAb and nahAa) that catalyse the complete pathway for the degradation of naphthalene (KEGG module: M00534) to salicylate (Fig. 5). Comparisons of the strain BUN14 genome with those of 15 related Pseudomonas species, including its closest relatives, P. kunmingensis DSM 25974 T , CCUG 36651 and P. chloritidismutans AW-1, showed that only BUN14 and P. balearica DSM 6083T harboured the complete pathway for naphthalene degradation (Fig. 6). Apart from DSM 18231T, which has an orthologue of nahD, all the compared strains lacked the orthologues of nahD and nahE. However, DSM 18231T lacked other genes associated with the naphthalene degradation pathway, and the other Pseudomonas genomes lacked between two and all nine genes of the complete pathway. The catechol degradation pathway has been shown to be the major catabolic route by which many bacteria biodegrade naphthalene 42 and a repertoire of genes linked to the compete naphthalene degradation pathway has been reported in many taxa, including Pseudomonas, Paraburkholderia, Alcaligenes and Rhodococcus 23,42-45 . The BUN14 genome encoded a membrane bound acyl-CoA desaturase (WP_104098302.1; FADS-like; alkB like), rubredoxin-2 (WP_003283319.1) and rubredoxin-NAD(+) reductase (WP_104098084.1), which are likely to determine its ability to degrade n-alkane. Orthologues of these three proteins were also found in the proteomes of nine of the compared genomes including the closest relatives of BUN14, P. kunmingensis DSM 25974 T , CCUG 36651 and P. chloritidismutans AW-1 (Fig. S10). Unlike P. chloritidismutans, no alkane 1-monooxygenase (AlkB; WP_023446487.1) genes were identified in the strain BUN14 genome. Most AlkB proteins have been shown to carry membrane fatty acid desaturase (FADSs), the actions of which depend on rubredoxin and rubredoxin reductase, encoded on a distinct genomic region 10 . Both rubredoxin and rubredoxin reductase homologues were identified in the BUN14 genome, suggesting that either an alternative, and unknown, enzyme catalyses the initial activation of alkanes, or that the sequence homology of the BUN14 putative AlkB gene was too low to be detected by Blast analysis. www.nature.com/scientificreports/ Given the identification of putative rhamnolipids in the extracted biosurfactant fraction of BUN14 cultures, the proteome was scanned for proteins potentially involved in the synthesis of this group of compounds. Blastp searches against 3-(3-hydroxydecanoyloxy) decanoate synthase RhlA (Q51559), rhamnosyl transferase RhlB (Q51560) and Rhamnosyltransferase 2 RhlC (Q9I4K5) of P. aeruginosa did not yield orthologous of these proteins. The BUN14 proteome was subsequently annotated for glycosyltransferase (GT1 and GT2) signatures using dbCAN2 46 and for conserved domains using the NCBI conserved domains database search tool 47 . This combined strategy revealed four proteins (Table S3) that showed similar domain architecture to P. aeruginosa RhlA, RhlB and RhlC. WP_104098192.1 (296 aa) which shared only 13.6% similarity with Q51559 (295 aa) of P. aeruginosa, harboured conserved domains similar to Q51559 and was identified as a specific hit for the alpha/beta hydrolase superfamily.Only one BUN14 protein, WP_104098902.1 (310 aa), was predicted to be a member of the glycosyltransferase family GT1. This protein shared 22.1% similarity with Q51560 (325 aa) of P. aeruginosa and the domain architecture of both proteins was similar. Seven BUN14 proteins were predicted to be members of the glycosyltransferase family GT2. Of these, two hypothetical proteins, WP_042926909.1 (291 aa) and WP_003300529.1 (357 aa) were predicted to harbour functional domains similar to those of Q9I4K5 (325 aa) of P. aeruginosa. However, WP_042926909.1 and WP_003300529.1 share only 20.5 and 22% similarity, respectively, with Q9I4K5. Therefore, WP_104098192.1 and WP_104098902.1 could potentially carry the molecular functions of the RhlA and RhlB proteins, respectively, while WP_042926909.1 and WP_003300529.1 are potential RhlC functional homologues. Indeed, the genetic characterization of rhamnolipid biosynthesis pathways is hindered by high levels of sequence diversity in the key genes, as seen in the low levels or absence of homology in the rhlA-C genes/proteins in P. aeruginosa, Burkholderia and related organisms 48 .

Conclusion
This study has demonstrated the versatility of Pseudomonas strain BUN14 in the degradation of a wide range of polyaromatic and aliphatic hydrocarbons. The genome of strain BUN14 was fully sequenced to identify the genetic determinants of its functional capacity, including the ability to utilize various hydrocarbons and to produce biosurfactants. The complete degradation pathways of some hydrocarbons such as benzoate and naphthalene were identified in the BUN14 genome, while the absence of some key enzyme genes suggested that the genome may harbor genes encoding novel functionalities, or functional homologues with very low sequence homology. On the basis of whole genome phylogeny and other genomic metrics (ANI and iDDH), the current study demonstrates that P. chloritidismutans AW-1 T is the closest relative of strain BUN14. Furthermore, P. chloritidismutans AW-1 T and P. kunmingensis DSM 25974 are conspecific and given that the former has been validly described prior to P. kunmingensis DSM 259, we suggest that DSM 25974 strain should be considered as a heteronym of P. chloritidismutans AW-1 T .

Materials and methods
Strain isolation and carbon source utilization. Strain BUN14 was isolated by enrichment culture in mineral salts medium (ONR7a) supplemented with 1% crude oil as a sole carbon source, after incubation for 21 days at 30 °C with shaking (150 × g) 5  1% of crude oil was used as only energy and carbon source 45,49 . To determine the phylogenetic affiliation of the strain, the 16S rRNA gene was amplified, sequenced, and queried using the 16S rRNA identification module in EzBioCloud 18,19 . The 16S rRNA sequences of the closest type strain relatives of BUN14 obtained from the Ezbiocloud database were aligned using Mafft 50 . Gaps were removed using the default setting in Gblocks v.0.91b 51 . A phylogenetic tree was constructed by the neighbour joining method and the tree topology was evaluated by performing bootstrap analysis of 1000 data sets using MEGA6.0 52 .
The cellular morphology of BUN14 was visualized using a JCM-5700 Scanning Electron Microscope, resolution 0.6 nm, specimen size 5 mm ∅ × 0.6 mm high, with a Gatan Digital Micrograph imaging system and SE & BS detectors 53 . BUN14 cultures were centrifuged (10 min 14,000g) and bacterial biomass was fixed in 2.5% gluteraldehyde in 0.075 M K-phosphate buffer (pH 7.4) for 2 h at room temperature. The preparation of the Scanning Electron Microscope (SEM) was performed according to Stanton et al. 53 . To test the ability of strain BUN14 to utilize hydrocarbon substrates as the sole carbon and energy source, the organism was inoculated onto solid ONR7a agar media containing individual hydrocarbons (pristane, phenanthrene, pyrene, naphthalene, fluoranthene, dibenzothiophene, dibenzofuran, squalene, and carbazole: 50 mM final concentration) and incubated for 7 days at 30 ± 1 °C. Growth of bacterial colonies was considered as a positive indication of the capacity to use the substrate as the sole carbon and energy source.
Optimization of surfactant activities. Strain BUN14 was screened for biosurfactant production using the Cetyl Trimethyl Ammonium Bromide agar plate assay (CTAB), the drop collapse test, emulsion index (E 24 ) and the oil-displacement test (ODA) 5,54,55 . For the optimization of surfactant activities, response surface methodology (RSM) using central composite design (CCD) was performed (Table S4). A CCD composed of 29 experiments have been planned (Table S5). The effects of various independent variables (substrate (waste frying oil) concentration (X1), NaCl concentration (X2), inoculum size (X3) and incubation time (X4)) on the biosurfactant production yield (Response Y) were evaluated at three levels (Table S6). The specific codes for each independent variable and range are given in Table S6. The relationship between these variables and the biosurfactant production yield was defined by the following model: www.nature.com/scientificreports/ where Y are the response (emulsion index (E 24 ) and oil displacement activity (ODA); Xj is the variables of studied factors and b 0 , b j , b jk , and b jj : model coefficients.
The five replicates at the center point were carried out in order to estimate the pure error variance. The significance of the fitted model was tested by the means of the analysis of variance (ANOVA). The relationship between the response and the experimental variables was illustrated graphically by plotting the response surfaces. The NemrodW software was used for experimental design and statistical analysis 56 . The optimum conditions for maximum biosurfactant production was determined using the desirability functions 57 .
The extraction of the biosurfactant was performed using liquid-liquid extraction 20 . The supernatants of BUN14 cultures were collected after centrifugation (12,000 rpm/20 min at 4 °C) adjusted to pH 2.0 with 6 N HCl and left overnight at 4 °C. The precipitates surfactant was collected by centrifugation at 12,000 rpm for 30 min at 4 °C. For additional purification, the crude biosurfactant was extracted at three successive washes with a mixture of the chloroform-methanol (2:1, v/v) and concentrated using rotary evaporation at 40 °C 20 . Functional group evaluation of the extracted biosurfactant were performed using Fourier Transform Infra-Red spectroscopy (Perkin Elmer FTIR model 2000) as described previously 20,58 . Infrared absorption spectra were obtained over the range of 400-4000 cm −1 with a resolution of 4 cm.
Hydrocarbon degradation. For hydrocarbon degradation analysis, strain BUN14 was cultured in ONR7a liquid mineral medium supplemented with 5 g L −1 Na-acetate for 72 h in shaking (150 × g) at 30 ± 1 °C. Bacterial cells in logarithmic phase were collected by centrifugation (10 min, 14,000g), washed twice in phosphate buffered saline (PBS 1×; 140 mM NaCl, 2.7 mM KCl, 4.3 mM Na 2 HPO 4 ·7H 2 O and 1.5 mM KH 2 PO 4 ), resuspended and inoculated (~ 10 6 cells ml −1 measured by the DAPI count method) into 50 ml ONR7a liquid mineral medium supplemented with sterile Arabian Light Crude Oil (1%, v/v), pristane (1% v/v) and 50 ppm (final concentration) of naphthalene, DBT and DBF (Sigma Aldrich, Milano-Italy). Cultures containing the same amount of hydrocarbon but without inoculation were used as abiotic controls. Cultures were incubated at 30 ± 1 °C for 21 days with shaking. Total Extracted and Resolved Hydrocarbons and their derivatives (TERHCs) were extracted from cultures using dichloromethane (Sigma-Aldrich, Milan; 10% v/v) following the 3550C EPA (Environmental Protection Agency) procedure as previously reported 59,60 . Degradation rates were quantified by a Master GC DANI Instruments GC-FID (Development ANalytical Instruments DANI Instruments S.p.A., Milan, Italy), equipped with SSL injector and FID detection. The extent of biodegradation was expressed as the percentage of hydrocarbon degraded compared to the abiotic control 60 .
BUN14 genome sequencing, assembly, annotation and analysis. DNA extraction was performed on mid-log phase cells by sodium dodecyl sulfate (SDS)-proteinase K treatment with an additional equal volume of chloroform/isoamyl alcohol (24:1 v/v). Purified genomic DNA was sequenced on an Illumina MiSeq platform (MRDNA, Clearwater, Tx, USA). The 11,335,376 paired reads were filtered according to read quality, and reads below a mean quality, score of 23 were removed using prinseq-lite software. The reads were assembled using SPAdes 61 . The genome of strain BUN14 was structurally annotated using PROKKA 62 and FGENESB (http:// www. softb erry. com/). The CGView Server 63 was used for circular representation of multiple genomes. The draft genome of strain BUN14 was used as the reference genome and was compared with genomes of P. kunmingensis DSM 25974, P. kunmingensis CCUG 36651 and P.chloritidismutans AW-1. Functional annotation was accomplished using a combination of RAST 63 , BlastKOALA 37 , KEGG 64,65 and EggNOG 66 . Orthologous relationships among the predicted proteomes of strain BUN14 and close relatives were determined using orthofinder 67 . Single copy orthologues were aligned using mafft 50 , trimmed using Gblocks v.0.91b 51 and maximum likelihood phylogeny was constructed using iq-tree 68 . Orthologous average nucleotide identity (OrthoANI) and in silico DNA-DNA hybridization (iDDH) were determined using OAT 69 and GGDC 38 , respectively. Comparisons of genomic regions among related Pseudomonas species were accomplished using SimpleSynteny 70 and a heat map was generated using ClustVis 71 .
Nucleotide sequence accession number. The genome sequence has been deposited in the Bioproject and Biosample Genomes online database under PRJNA420855 and SAMN08122818 accession numbers, respectively.
Whole Genome Shotgun sequence data has been deposited at DDBJ/ENA/GenBank under the accession number PISM00000000.