Complete genome sequence of Pseudomonas alcaliphila JAB1 (=DSM 26533), a versatile degrader of organic pollutants

In this study, following its isolation from contaminated soil, the genomic sequence of Pseudomonas alcaliphila strain JAB1 (=DSM 26533), a biphenyl-degrading bacterium, is reported and analyzed in relation to its extensive degradative capabilities. The P. alcaliphila JAB1 genome (GenBank accession no. CP016162) consists of a single 5.34 Mbp-long chromosome with a GC content of 62.5%. Gene function was assigned to 3816 of the 4908 predicted genes. The genome harbors a bph gene cluster, permitting degradation of biphenyl and many congeners of polychlorinated biphenyls (PCBs), a ben gene cluster, enabling benzoate and its derivatives to be degraded, and phe gene cluster, which permits phenol degradation. In addition, P. alcaliphila JAB1 is capable of cometabolically degrading cis-1,2-dichloroethylene (cDCE) when grown on phenol. The strain carries both catechol and protocatechuate branches of the β-ketoadipate pathway, which is used to funnel the pollutants to the central metabolism. Furthermore, we propose that clustering of MALDI-TOF MS spectra with closest phylogenetic relatives should be used when taxonomically classifying the isolated bacterium; this, together with 16S rRNA gene sequence and chemotaxonomic data analyses, enables more precise identification of the culture at the species level. Electronic supplementary material The online version of this article (10.1186/s40793-017-0306-7) contains supplementary material, which is available to authorized users.


Introduction
Over recent decades, significant quantities of potentially harmful chemicals have been released into the environment, creating countless numbers of contaminated sites. Major contaminants include halogenated and nitrated alicyclic, aliphatic, aromatic and polyaromatic compounds of industrial and agricultural origin. Many of these compounds have been found to have toxic, mutagenic and carcinogenic effects on living organisms. Removal of these xenobiotics usually involves physical and chemical processes, such as landfill, excavation and incineration, which are expensive and difficult to execute. An alternative approach, bioremediation, uses ubiquitous plant-microbe interactions to degrade xenobiotics [1]. Bacteria and fungi are natural recyclers capable of funneling toxic organic compounds to central metabolism intermediates [2], thereby creating harmless products [3] such as carbon dioxide and water. In addition to the enormous phylogenetic diversity of microorganisms, the richness of their metabolic activities promotes the degradation of pollutants and xenobiotics in different environments.

Organism information
Classification and features P. alcaliphila was described as a facultatively psychrophilic alkaliphilic species isolated from seawater off the coast of Rumoi, Hokkaido, Japan. The characteristics of this species are as follows: alkaliphile, incapable of growth at >40°C, catalase-and oxidase-positive and also capable of reducing nitrate to nitrite and of hydrolyzing casein and gelatin [21]. Further physiological features are listed in Table 1. The JAB1 cells are monotrichous rods as shown in Fig. 1.
The JAB1 strain was originally misidentified as P. pseudoalcaligenes [22]. The consensus 16S rRNA gene sequence, compiled from four 16S rRNA gene copies contained in the JAB1 genome, had 99.93% similarity to those of P. alcaliphila AL 15-21 T [21], P. chengduensis  [24]. Additional closest matches included P. toyotomiensis HT-3 T (99.86% similarity) and P. mendocina CH50 T (99.24% similarity). A phylogenetic tree indicates closest relatedness of the JAB1 strain to P. alcaliphila AL 15-21 T (Fig. 2). In addition, the JAB1 strain was unable to grow at 41°C, which is a typical feature of P. alcaliphila but not of other closely related pseudomonads [23]. Furthermore, whole-cell MALDI-TOF MS analysis, performed following the methodology described elsewhere [12], indicated that JAB1 spectra clustered with those of P. alcaliphila AL 15-21 T (Fig. 3). The results of MALDI-TOF MS profiling thus further confirmed the identity of the JAB1 strain as P. alcaliphila. Therefore, we propose that MALDI-TOF MS analysis be performed of the isolate and its closest phylogenetic relatives when taxonomically classifying the isolated bacterium. In addition to 16S rRNA gene sequence and chemotaxonomic data analysis, MALDI-TOF MS can provide additional information resulting in more precise classification of the isolate.

Extended feature descriptions: biodegradative capabilities
The JAB1 strain was isolated from legacy contaminated soil in Jablonné nad Orlicí in the Czech Republic as a biphenyl-degrading bacterium capable of extensive degradation of several congeners of PCBs when grown in the presence of biphenyl [22]. PCB-degradation capabilities were further determined by a resting cell assay. Briefly, microcosms (volume of 119 ml) bearing biphenylpreinduced P. alcaliphila JAB1 mid-log phase culture resuspended in mineral salt solution (20 ml per microcosm) were spiked with a 0.001% (w/v) commercial mixture of Creek, CA). PCBs were analyzed in ethyl acetate extracts according to the method described by Čvančarová M. et al. [25]. CBAs were analyzed in the extracts using GC-MS after methylation with diazomethane according to our previously published protocol [26]. The respective chemical standards for the analytes were obtained from Merck (Darmstadt, Germany), Supelco (Steinheim, Germany), TCI Europe (Zwijndrecht, Belgium) and AccuStandard (New Haven, USA). Metabolic activity of the JAB1 strain resulted in the depletion of various mono-, di-, tri-and tetra-chlorinated biphenyls as shown in Fig. 4. At the same time, 2-CBA, 3-CBA, 4-CBA, 2,3-diCBA, 2,4-diCBA and 2,5-diCBA formation (data not shown) was observed over the course of the coincubation period; these are common biodegradation intermediates of various PCB congeners and products of the upper PCB degradation pathway.
The degradation of cDCE was assessed during growth on phenol, a compound known to promote cDCE degradation in microbial cultures [27]. Microcosms (volume of 119 ml) bearing P. alcaliphila JAB1 culture (20 ml per microcosm) grown and repeatedly propagated on 1 mM phenol and 0.01 mM cDCE were monitored for phenol and cDCE depletion as described elsewhere [28] as well as the increase in microbial culture density measured spectrophotometrically. Additional cultures were prepared in the presence of 0.01 mM cDCE as sole carbon source and of 0.01 mM cDCE and 1 mM sodium pyruvate to exclude direct utilization of cDCE by JAB1.   [62] using secondary structure-based alignment of 16S rRNA gene sequences of the JAB1 strain (four copies of the 16S rRNA gene retrieved from the JAB1 genome) and its closest relatives [63]. The evolutionary distances were computed using the Kimura 2-parameter method [64]. All positions with less than 95% site coverage were eliminated; there were a total of 1375 positions in the final dataset. The tree construction method used was Neighbor-Joining [65]. The bootstrap test (1000 replicates) was used to test tree topology [66]; only bootstrap values >70 are indicated In addition, SPME-GC-MS analysis of culture headspaces was performed over time in order to identify possible volatile chlorinated intermediates of the cDCE degradation pathway.
Monitoring of cDCE levels in cultures grown on 1 mM phenol revealed consistent cDCE depletion after the first 14 h of cultivation, when phenol (the growthsupporting substrate) was completely or almost completely consumed (Fig. 5). After 72 h of cultivation, 100% cDCE depletion was attained exclusively in phenol-grown cultures; during the same time period, cDCE loss accounted for only 13.4% ± 8.8 pp. in sodium pyruvate-grown cultures, 9.7% ± 9.1 pp. for JAB1 cultures exposed solely to cDCE and 12.5% ± 1.2 pp. in sterile medium. Such minimal cDCE loss under conditions tested can be ascribed to non-biological processes. Proven total cDCE depletion only in phenol-grown cultures confirms our hypothesis that cDCE degradation occurs cometabolically. To date, 1,1-dichloroacetaldehyde was confirmed as a degradation intermediate by means of SPM-GC-MS in culture headspaces. Quantification of this molecule and its accumulation over time will be an area for future research together with the study of possible metabolic pathways to promote further transformation of 1,1-dichloroacetaldehyde.

Growth conditions and genomic DNA preparation
Following its isolation from soil, the culture was preserved in a mineral salt solution [12] with biphenyl as sole carbon source and, over the long term, in glycerol stocks prepared from growing liquid cultures. For the purposes of DNA isolation, the culture was grown overnight on plate count agar (Difco, UK) at 28°C. Genomic DNA was isolated using the PureLink™ Genomic DNA Mini Kit (Invitrogen, USA) according to the manufacturer's instructions.

Genome sequencing and assembly
The genomic DNA of P. alcaliphila JAB1 was used to prepare shotgun and 8 kb paired-end 454 sequencing libraries according to the Library Preparation Method Manual (Roche). These libraries were sequenced with the GS FLX instrument using GS FLX+ chemistry (Roche) at the Institute of Molecular Genetics AS CR (Prague, Czech Republic). The resulting 110,702 shotgun and 461,976 paired-end reads comprised a total of 221,447,771 bases and represented 41.5-fold genome coverage. These reads were assembled in Newbler 2.8 software (Roche), and gaps were closed by local assembly of selected reads in Staden software [29] into a single circular contig.

Genome annotation
The NCBI Prokaryotic Genome Annotation Pipeline [30] was used for automated genome annotation, and functional annotations were carried out by searching against KEGG [31], COG [32] and Pfam [33] databases. SignalP [34] and TMHMM [35] tools were used for the prediction of genes with signal peptides and transmembrane helices, respectively. Searches against NCBI-NR [36], RefSeq [37], UniProtKB/Swiss-prot [38] and TCDB [39] were carried out for functional assignment of selected genes.

Genome properties
The P. alcaliphila JAB1 genome consists of a single 5,340,293 bp-long chromosome with a GC content of 62.5% (Fig. 6). The 4908 predicted genes correspond to an 89% coding density. The genome contains 4773 CDS, 65 tRNA genes for all 20 amino acids, 4 rRNA operons, 4 ncRNA genes, 3 CRISPR repeats and 54 pseudogenes. Function was assigned to 3816 CDS, with 957 CDS annotated as hypothetical proteins only ( Table 3). The distribution of COG functional categories is shown in Table 4.

Insights from the genome sequence
The genome harbors approximately 50 oxygenase genes, many of which are relevant to the biodegradation capabilities of the strain, especially those encoding for mono-and di-oxygenases responsible for the hydroxylation and opening of the aromatic ring, which are crucial steps in aromatic xenobiotic degradation (see below for details). Several other genes found are associated with the heavy metal resistance of the strain (e.g. heavy metal efflux P-type ATPases), indicating its overall adaptation to contaminated environments; organic compound contamination is often accompanied by high concentrations of heavy metals [40,41]. The genome contains a complete set of genes for flagellar assembly and 50 additional genes associated with chemotaxis. Roughly 70 of the predicted CDS are putative intact or mutated transposase genes and phage-related genes.

Extended insights: degradation of aromatic compounds
In total, five regions harboring genes for aromatic compound degradation pathways were identified in the JAB1 genome as shown in Additional file 1: Table S1.
The first region, position range 1,474,745-1,486,552 in the annotated genome sequence, harbors genes encoding for a three-component benzoate 1,2-dioxygenase (benABC genes) followed by enzymes of a pathway enabling the complete degradation of benzoates and benzoate derivatives, yielding acetyl-CoA and pyruvate (catechol metacleavage pathway encoded by benD, fdx, and dmp cluster). The role played by the gene designated nahX encoding for putative ATP-cob(I)alamin adenosyltransferase in aromatic compound degradation is unclear [42]. No transcriptional regulator was found within or adjacent to this region, leading us to assume that genes from this region are not transcribed or their transcription is mediated in trans by regulatory proteins encoded elsewhere.
The second region, position range 1,518,873-1,540,279, harbors genes encoding for complete salicylate degradation via the catechol meta-cleavage pathway (sal/fdx genes), with salicylate 1-monooxygenase (encoded by salA) being responsible for opening the aromatic ring. Genes antAB encoding for large and small subunits of a putative terminal dioxygenase were also found in this region. Apart from genes presumed to be involved in the enzymatic degradation of aromatic compounds, ORFs encoding for transport proteins were found in this second region, namely a putative benzoate:H + symporter (encoded by gene benE), a multi-component ABC transporter (six ORFs, designated ORF1-6) and a porin protein (benK), which are most likely involved in aromatic compound uptake. BenK ortholog has been shown to be involved in benzoate uptake in Acinetobacter sp. ADP1 [43]. Moreover, BenE and BenK orthologs from P. putida have been demonstrated to act as benzoate uptake proteins when heterologously expressed in yeasts [44]. Nevertheless, detailed experimental evidence is lacking on the function of these proteins and their role in aromatics/aromatic xenobiotics degradation, which represents a potential challenge for future research [45,46]. The third region, position range 1,545,237-1,557,889, harbors bph genes encoding for enzymes of the complete upper and lower biphenyl degradation pathway. The structure of this operon is virtually identical to that of the bph operon found in the genome of the model PCB-degrading strain P. pseudoalcaligenes KF707 [47]. While pentadiene and its chlorinated derivatives can be degraded by enzymes encoded by the same operon (bphHIJ genes), no ORF encoding for benzoate terminal dioxygenase, essential for the completion of biphenyl degradation, is present in this region; nevertheless, such genes are found elsewhere in the JAB1 genome. Two other ORFs, tbuX and a gene designated ORF7, encoding for putative aromatic-transporter Fig. 6 Circular map of the P. alcaliphila JAB1 chromosome. From outside to the center: CDS on forward strand colored according to their COG categories, all CDS and RNA genes on forward strand, all CDS and RNA genes on reverse strand, CDS on reverse strand colored according to their COG categories, GC content, and GC skew. The map was generated using CGView [67] transmembrane proteins, were found in this region. Although the TbuX protein has been reported to be involved in the utilization of toluene in Ralstonia pickettii PKO1 [48], no detailed information exists on the precise role played by this group of proteins in aromatic compound transport. The predicted ORF7-encoded protein exhibits homology with members of the BphX family. Despite being commonly found in upper biphenyl degradation pathwayencoding supra-operonic clusters of various taxa, the function of these proteins in aromatic degradation is not known [49][50][51].
The fourth region, position range 4,397,470-4,411,693, harbors gene cluster pheKLMNOP encoding for multicomponent phenol 2-monoxygenase and an fdx gene encoding a chloroplast-type ferredoxin essential for electron delivery to the active center of the oxygenase. The adjacent dmp cluster encodes for enzymes responsible for the meta-cleavage pathway of catechol/catechol derivatives, yielding acetyl-CoA and pyruvate. In this study, we demonstrate that strain JAB1 is capable of utilizing phenol as a sole carbon source and cometabolite of cDCE degradation.
The fifth region, position range 4,445,675 to 4,469,934, contains benABCD genes, which are responsible for the transformation of benzoate into catechol, as well as cat-ABC and pcaDIJF genes, which encode the enzymes of the catechol branch of the β-ketoadipate pathway, yielding succinyl-CoA and acetyl-CoA [52]. This region also harbors pcaHGB genes for the protocatechuate branch of the β-ketoadipate pathway which enable hydroxybenzoate degradation; the pobA gene encoding 4-hydroxybenzoate 3-monooxygenase, responsible for the first step of 4hydroxybenzoate hydroxylation, is located upstream of this fifth region (positions 4,435,851 to 4,437,041). In addition, this region harbors three ORFs encoding for putative transport proteins, which are most likely involved in aromatic compound uptake, designated as ORF8 (encoding for outer membrane porin from the OprD family), pcaK (an ortholog of p-hydroxybenzoate and protocatechuate transport protein from the Major Facilitator Superfamily [53]) and benE (benzoate:H + symporter from the BenE family). Two other ORFs encoding for a cell division protein (zapE gene) [54] and a permease of unknown substrate specificity and function (designated ORF9) are also present in the fifth region in a pcaHG-zapE-ORF9 arrangement. As the intergenic region between the pcaG and zapE genes is only 1 bp in length, with zapE and ORF9 sequences even overlapping by 5 bp, we hypothesize that pcaHG-zapE-ORF9 represents a single transcriptional unit. The positioning of one or more genes involved in cell division as part of a cistron together with aromaticdegradation genes is rather unusual. The total is based on the total number of protein coding genes in the genome Since the regions harboring degradation-determining genes are flanked by transposase and conjugation proteincoding genes, and, at the same time, they generally exhibit analogous architecture and high similarity to corresponding orthologs from other pseudomonads, we hypothesize that their acquisition in the JAB1 genome was most likely due to multiple horizontal gene transfer events. However, further study is required to confirm this hypothesis.

Conclusions
Thorough microbiological experiments and whole genome sequence analysis lead us to conclude that P. alcaliphila JAB1 is a versatile degrader of organic pollutants. Its extensive degradative capabilities are enabled by a variety of genes determining the degradation of both biphenyl/biphenyl derivatives (mediated by the biphenyl degradation pathway encoded by the bph cluster) and monocyclic aromatics such as benzoate and its derivatives (halobenzoates, toluates), salicylic acid and phenol/phenolic compounds.

Additional file
Additional file 1: Table S1. Description of five regions harboring genes for aromatic compound degradation pathways were identified in the JAB1 genome.