Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2008, p. 6405-6416, Vol. 74, No. 20
0099-2240/08/$08.00+0 doi:10.1128/AEM.00197-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Department of Civil and Environmental Engineering, 4105 Seamans Center, The University of Iowa, Iowa City, Iowa 52242,1 Department of Energy Joint Genome Institute, 2800 Mitchell Drive, Walnut Creek, California 94598,2 Los Alamos National Laboratory, Los Alamos, New Mexico 87545,3 Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada T6G 2P5,4 School of Molecular and Microbial Biosciences, Building G08, University of Sydney, Sydney, New South Wales 2006, Australia5
Received 22 January 2008/ Accepted 12 August 2008
|
|
|---|
|
|
|---|
Polaromonas strains can biodegrade important groundwater contaminants, including petroleum hydrocarbons (6, 24, 38) and chlorinated ethenes (6, 10, 38). Petroleum hydrocarbons include aliphatic, aromatic, and polycyclic compounds, all of which are naturally occurring but problematic due to their toxicity and persistence (3). Chlorinated ethenes are essentially xenobiotic and are synthesized for the manufacture of plastics (vinyl chloride) and for use as solvents (tetrachloroethene, trichloroethene, and cis-dichloroethene [cDCE]). Cleanup of sites contaminated with chlorinated ethenes and/or petroleum hydrocarbons by using physical or chemical methods is difficult and expensive, so there has been intense interest in finding bacteria capable of degrading these pollutants for use in bioremediation.
Polaromonas sp. strain JS666 was isolated from a granular activated-carbon filter treating chloroethene-contaminated groundwater and is capable of growth-coupled cDCE oxidation and cometabolic degradation of diverse chlorinated aliphatic pollutants (10). As a result, strain JS666 is potentially useful for bioaugmentation of chlorinated-ethene-contaminated sites because growth of JS666 on cDCE in the subsurface would essentially convert the pollutant into a catalyst (i.e., new cells), leading to theoretically increasing degradation rates over time. Transfer of the biodegradative activity of JS666 to in situ conditions requires a better understanding of the physiology and biochemistry of the organism, especially with respect to the enzymes involved in organochlorine degradation, the ability to degrade or resist other pollutants, the range of nutrient sources used, and the types of environmental conditions preferred.
The metabolic and phylogenetic novelty of Polaromonas sp. strain JS666 and its potential for use in bioremediation provided the impetus for complete genome sequencing. Here, we describe a bioinformatic and phenotypic analysis that reveals diverse metabolic capabilities of JS666 that may be useful for bioremediation and biocatalysis. We also describe the roles of mobile genetic elements in shaping the genome and metabolic abilities of the bacterium.
|
|
|---|
3-kb inserts), pMCL200 (
7-kb inserts), or pCC1Fos (Epicentre, Madison, WI) (
35-kb inserts). Ligation reactions were used to transform Escherichia coli DH10B, and the resulting colonies containing plasmids were picked into 384-well plates containing LB and glycerol. Template DNA for sequencing was produced by amplifying plasmid subclones using rolling-circle amplification (Templiphi; GE Healthsciences, Piscataway, NJ) as described previously (14) or Sprintprep (Agencourt Biosciences, Beverly, MA) magnetic-bead DNA purification. Subclone inserts were sequenced from both ends using universal primers and Energy Transfer (GE Healthsciences, Piscataway, NJ) or Big Dye (ABI, Foster City, CA) terminator chemistry.
Sequence assembly and gap closure.
A total of 78,720 sequencing reads from the three libraries (3 kb, 33,792; 8 kb, 39,936; 40 kb, 4,992) were attempted and assembled with PHRAP (18) to produce a draft assembly containing 63 contigs. An additional 970 finishing reads were sequenced to resolve misassemblies, close gaps, and polish low-quality regions. The finished genome consists of one chromosome (5.2 Mb) and two plasmids (360 kb and 338 kb) containing no gaps and an error rate of less than 1 in 50,000 bases. Dupfinisher (C. S. Han and P. Chain, presented at the International Conference on Bioinformatics and Computational Biology, June 26-29, 2006) and transposon bombing (Epicentre Biotechnologies, Madison, WI) were used to resolve misassemblies.
Sequence analysis and automated annotation.
Gene modeling was performed by using the Critica (4), Glimmer (12), and Generation (http://compbio.ornl.gov/public/tools) modeling packages. The results were combined, a BLASTP search of the translations versus GenBank's nonredundant (NR) database was conducted, and the best NR match was used to select the preferred gene model. If there was not a BLAST match, the longest gene model was retained. Overlapping gene models were manually examined, and genes with BLAST matches were retained. The revised gene-protein set was searched against the KEGG GENES, Pfam, PROSITE, PRINTS, PRoDom, and COGs databases, in addition to the BLASTP versus NR databases. Genes were functionally categorized using the KEGG and COGs hierarchies. Manual revisions were made to the machine annotation as appropriate. Chromosome and plasmid maps were generated using the program CGview (40). The IslandPath interface (20) was used to search for genomic islands based on dinucleotide bias and unusual G+C content (default settings, ±3.48% G+C), presence of an integrase gene, and proximity to a tRNA gene. Inteins were identified by the presence of a homing endonuclease domain and flanking conserved regions (A, B, F, and G) in the N- and C-terminal splicing domains; these analyses were performed by Francine Perler of New England Biolabs (34).
Growth experiments.
Either minimal salts medium (MSM) (10) or R2A medium (Oxoid) was used for culture growth. Fed-batch growth substrate tests were done in MSM, using test compounds (
97% purity) as sole sources of carbon and energy. Polaromonas sp. strain JS666 cells grown to mid-exponential phase on succinate, cDCE, or ethanol were inoculated (optical density at 600 nm [OD600] = 0.01 to 0.07) into 40 ml of MSM in 160-ml serum bottles with air headspace. Each compound was tested over a range of starting concentrations from 0.2 to 10 mM, with the concentration of volatile compounds estimated as if they were entirely dissolved in the aqueous phase. Bottles were sealed with Teflon-coated butyl rubber septa and aluminum crimp caps (Wheaton) and incubated at 20 to 22°C with shaking at 190 to 250 rpm. At least two independent growth experiments were conducted for each substrate, including a control with no carbon source and uninoculated controls for compounds yielding colored breakdown products (catechol, protocatechuate, ferulate, and gentisate). Chemolithotrophic thiosulfate and carbon monoxide oxidation cultures were also prepared as described above with the following exceptions: sodium thiosulfate (0 to 40 mM) or CO (10% of the headspace) was added as the energy source, and sodium bicarbonate (10 mM) was added as a carbon source. Duplicate sets of growth experiments were conducted for both lithotrophic and mixotrophic growth; in the latter case, ethanol (5 to 10 mM) was provided as a limited source of organic carbon.
Anaerobic nitrate reduction cultures were prepared as described above, except 160 ml of MSM was used, leaving no headspace in the serum bottles, and the medium was purged with nitrogen gas to create anaerobic conditions. Ethanol (5 mM) was provided as a carbon source, and samples were removed while the bottles were being purged with nitrogen gas.
Auxanography plates were prepared with 1.5% Noble agar in MSM (made with high-performance liquid chromatography grade water). Exponentially growing succinate cultures were diluted in MSM to an OD600 of 0.1, and 100 µl was spread on each plate. Solid-phase test compounds (5 mg) were deposited on the edges of the plates to allow a concentration gradient to form in the agar, and the starting locations of the compounds were marked. The plates were sealed with Parafilm, inverted, and incubated at 20 to 22°C for 2 to 8 weeks. All compounds were tested with duplicate plates in two independent experiments, as were the no-carbon and no-inoculum controls.
PFGE.
Agarose plugs were prepared and pulsed-field gel electrophoresis (PFGE) was performed as described previously (28), except that succinate-grown JS666 cells were suspended in PFGE grade agarose at an OD600 of 3.
Analytical methods.
Growth was indirectly estimated in liquid cultures by using the OD600. Following auxanography, cells were scraped from the plates and then resuspended in MSM and washed, and total protein was analyzed by a previously described UV protein assay (10). Heptane, hexane, octane, 1,2-dichloroethane, and cDCE were analyzed in headspace samples (100 µl) by gas chromatography under isothermal conditions with flame ionization detection and a J&W Scientific GSQ (30-m by 0.53-mm inside diameter) column. Nitrate and nitrite concentrations in anaerobic cultures were measured on a Dionex DX600 ion chromatograph with a Dionex IonPac AS18 hydroxide-selective anion-exchange column.
Nucleotide sequence accession numbers.
The finished, annotated sequence was deposited in GenBank under accession numbers CP000316 to CP000318 and is also available at http://genome.jgi-psf.org/finished_microbes/poljs/poljs.info.html.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Features of the Polaromonas sp. strain JS666 genomea
|
![]() View larger version (42K): [in a new window] |
FIG. 1. The chromosome of Polaromonas sp. strain JS666. Approximate locations of important gene clusters are indicated outside of the first circle, which depicts coding regions on the plus strand colored by functional category: translation, ribosomal structure, and biogenesis (yellow); transcription (magenta); DNA replication, recombination, and repair (red); cellular processes (sky blue); signal transduction (lavender); energy production and conversion (green); carbon and carbohydrate metabolism (blue); amino acid transport and metabolism (orange); nucleotide transport and metabolism (light red); coenzyme metabolism (pink); lipid metabolism (cyan); secondary-metabolite biosynthesis, transport, and catabolism (dark green); general function prediction only (brown); function unknown (dark gray); and no COG (light gray). Moving toward the center, the second circle depicts coding regions on the minus strand (with the same color scheme as the plus strand). The third circle shows G+C content (deviation from average), and the fourth circle G+C skew in purple and olive. Genomic islands (GI1 to GI8) and the Mu phage are indicated by black bars between the third and fourth circles. The scale (in kbp) is indicated on the innermost circle.
|
![]() View larger version (26K): [in a new window] |
FIG. 2. (A) The 360-kb plasmid pPol360 of Polaromonas sp. strain JS666. (B) The 338-kb plasmid pPol338 of Polaromonas sp. strain JS666. Approximate locations of important metabolic gene clusters, plasmid replication, and conjugation genes are indicated outside of the first circle, which depicts coding regions on the plus strand colored by functional category: plasmid transfer (orange); plasmid replication and maintenance (yellow); transport (light green); DNA replication, recombination, and repair (red); metabolism (cyan); biosynthesis (dark blue); regulators (purple); cell envelope and membrane biogenesis (pink); general function prediction (brown); conserved hypothetical (dark gray); and no COG (light gray). Moving toward the center, the second circle depicts coding regions on the minus strand (with the same color scheme as the plus strand). The third circle shows G+C content (deviation from average), and the fourth circle G+C skew in purple and olive. The scale (in kbp) is indicated on the innermost circle. The base pair numbering for both plasmids begins near the putative origins of replication (i.e., repB [Bpro_5513] in pPol338 and trfA [Bpro_5045] in pPol360). (C) PFGE of genomic DNA (gDNA) from Polaromonas sp. strain JS666 depicting the presence of two extrachromosomal genetic elements. L, 50-kb lambda ladder.
|
![]() View larger version (20K): [in a new window] |
FIG. 3. Functional distribution of COGs among replicons in the JS666 genome sequence. COG designations: J, translation and ribosomal structure/biogenesis; K, transcription; L, DNA replication, recombination, and repair; C, energy production/conversion; E, amino acid transport/metabolism; G, carbohydrate transport/metabolism; H, coenzyme metabolism; I, lipid metabolism; P, inorganic-ion transport/metabolism; Q, secondary metabolism; M, cell envelope biogenesis outer membrane; N, cell motility and secretion; O, posttranslational modification and protein turnover; T, signal transduction mechanisms; U, intracellular trafficking and secretion; R, general function prediction only; and S, function unknown.
|
-like initiator (Bpro_5045; 28% amino acid identity to trfA of RP4/RK2) associated with 11 direct repeats of a 22-bp iteron (MRNCGTGTTNGWCCGCTAATAG), parA and parB genes, and an exonuclease.
![]() View larger version (34K): [in a new window] |
FIG. 4. (A) Comparison of the putative plasmid replication regions of pPol338 and pPol360. (B) Alignment of the putative plasmid conjugation gene regions of pPl338 and pPol360 with the conjugation gene regions from the hydrogen catabolic plasmid HG1, IncP- and Ti plasmids, and the integrative conjugative element R391.
|
and Ti plasmids, but the closest database matches (25 to 30% amino acid identity) are to gene products from the hydrogen catabolic plasmid pHG1 (Trb), integrative conjugative elements (R391, SXT; Tra), and IncH plasmids (R27 and R478; Trh). The pil genes of the JS666 plasmids (pilLNOPQRST1MT2V) are similar in arrangement to pilus biosynthesis genes of Pseudomonas aeruginosa elements (pKLC102 and PAPI-1) and the IncI1 plasmid R64. Both pPol338 and pPol360 carry multiple putative stabilization systems, including Burkholderia-like parAB gene pairs (Bpro_5043-Bpro_5044 and Bpro_5514-Bpro_5515) (15), solitary parB genes (Bpro_5053 and Bpro_5341), resolvases (Bpro_5040, Bpro_5516, Bpro_5529, and Bpro_5548), and DNA methylases (Bpro_5027 and Bpro_5389). pPol338 also carries a stbB antitoxin homolog (Bpro_5318) and a higAB-like killer-antidote gene pair (Bpro_5380-Bpro_5381). Other putative "backbone" genes present on both pPol338 and pPol360 include primase (Bpro_5218 and Bpro_5493), helicase (Bpro_5077 and Bpro_5490), exonuclease VIII (Bpro_5046 and Bpro_5358), and ribonucleotide reductase (Bpro_5207 and Bpro_5497) genes.
Other mobile genetic elements.
Diverse transposases are present in the JS666 genome, consisting of 55 genes of 40 different types (see Table S1A in the supplemental material). They include representatives of 14 insertion sequence (IS) families, a Mu-like transposase, and four Tn3-like transposases. Associations between transposases and genes of interest for biodegradation and biocatalysis were noted in many cases (see Table S2 in the supplemental material and below). There are eight putative genomic islands (see Table S1B in the supplemental material) and a 39-kb Mu family phage (Fig. 1; see Table S1C in the supplemental material). Several of the genomic islands contain genes of interest to biodegradation and biocatalysis, including aromatic-compound utilization (GI1), haloacid dehalogenase (GI2), and cytochrome P450 (GI5). Two types of putative retrons (type II introns) are present in the JS666 genome. The first type is present in four chromosomal copies (Bpro_949, Bpro_4601, Bpro_4617, and Bpro_4899) and contains a reverse transcriptase of 564 to 569 amino acids (aa) with a C-terminal endonuclease domain, while the second type has two copies on pPol338 (Bpro_5394 and Bpro_5396) associated with an IS605 family transposase (Bpro_5398) and a reverse transcriptase of 456 to 461 aa. JS666 carries two inteins, one (351 aa) in the DnaB helicase (Bpro_3033) and the other (380 aa) in a ribonucleoside reductase (Bpro_4588). Both retrons and inteins have potential uses in biotechnology, for example, as a source of "homing" endonucleases (39). Three copies of the 132- to 169-bp miniature inverted-repeat transposable element (MITE) "Nezha" are present in the JS666 genome (50). This MITE is related to cyanobacterial MITEs and is not found in other bacterial genomes, suggesting a recent horizontal gene transfer event from cyanobacteria to JS666 (50).
Central metabolism and C1 and C2 metabolism.
The JS666 genome contains complete glycolytic and Entner-Doudoroff pathways and a complete tricarboxylic acid cycle. The pentose phosphate pathway appears to be missing 6-phosphogluconolactonase. Glyoxalate pathway genes encoding isocitrate lyase (Bpro_2101) and malate synthase (Bpro_4517) are present, as expected for an organism isolated on a C2 compound (cDCE). Oxalate and glycolate catabolic genes were also noted (see Table S2 in the supplemental material). Two ribulose bisphosphate carboxylase homologs (Bpro_0032 and Bpro_0093) are in the genome; however, a phosphoribulokinase homolog was not found, suggesting that JS666 does not fix carbon via the Calvin cycle. Alternative pathways for carbon fixation, such as the C4-dicarboxylic acid cycle, the 3-hydroxypropionate pathway, the reverse tricarboxylic acid cycle, and the acetyl-coenzyme A (CoA) pathway, are either missing or incomplete, suggesting that JS666 is not capable of autotrophic CO2 fixation.
The presence of a monofunctional (CoxMSL-like) carbon monoxide dehydrogenase (Bpro_0576 to Bpro_0578) and the lack of a complete CO2 fixation pathway suggest that JS666 uses carbon monoxide (CO) as an energy source in the presence of an organic carbon source (i.e., mixotrophy). However, experiments with 10% CO in the headspace and 0 to 10 mM ethanol as a carbon source failed to show growth of JS666 on CO as a carbon and/or energy source. Genes encoding a putative tungsten-containing formate dehydrogenase are present, suggesting growth on formate as a carbon and/or energy source (see Table S2 in the supplemental material) (this was not tested experimentally).
Alternative electron donors and acceptors.
The presence of a gene cluster similar to the E. coli K-12 nar operon (narGHJI) (Bpro_4592 to Bpro_4599) (see Table S2 and Fig. S1 in the supplemental material) raises the possibility that JS666 uses nitrate as an alternative electron acceptor for growth. Incubating JS666 under anaerobic conditions with 5 mM ethanol and 5 to 35 mM sodium nitrate resulted in nitrite production (0.5 to 0.8 mM) but failed to show the use of nitrate as an electron acceptor for growth.
The coincidence of a putative hydrogenase gene cluster related to E. coli membrane-bound respiratory hydrogenases 3 and 4 (Bpro_0916 to Bpro_0921) (see Table S2 in the supplemental material) and several possible formate dehydrogenase genes (Bpro_0937 to Bpro_0942) suggests that JS666 possesses a fermentative formate hydrogenlyase system that converts formate into molecular hydrogen and carbon dioxide. However, this possibility was not tested experimentally.
The presence of potential sulfur oxidation genes (Bpro_2625 to Bpro_2631), encoding proteins with 25 to 50% identity to those from mixotrophic sulfur-oxidizing Paracoccus pantotrophus (16) could indicate the ability to grow chemolithotrophically by oxidizing thiosulfate, sulfite, hydrogen sulfide, and/or elemental sulfur in the presence of an organic carbon source (see Table S2 and Fig. S2 in the supplemental material). However, growth tests with thiosulfate (0 to 40 mM) as an electron donor and bicarbonate (10 mM) or ethanol (0 to 5 mM) as a carbon source were negative (Table 2).
|
View this table: [in a new window] |
TABLE 2. Growth substrate range of Polaromonas sp. strain JS666
|
Metabolism of haloalkanes and haloacids.
Microbial hydrolytic dehalogenase enzymes are important in the metabolism of diverse pollutants, such as chlorinated alkanes and aromatics (23). JS666 possesses putative haloalkane biodegradation genes, the products of which are 48 to 75% identical to those in the dichloroethane-degrading Xanthobacter autotrophicus strain GJ10 (Fig. 1, 2, and 5; see Table S2 in the supplemental material) (22, 46). They include homologs of haloalkane dehalogenase DhlA (Bpro_0521, Bpro_0547, and Bpro_2447), haloalkanoate dehalogenase DhlB (Bpro_0530 and Bpro_5186), the regulator DhlR (Bpro_0533 and Bpro_5183), and the uptake protein DhlC (Bpro_0531 and Bpro_5185). The presence of an apparently complete set of dichloroethane catabolic genes in JS666 is surprising, since previous tests (10) showed that cDCE-grown cells attacked dichloroethane but did not grow on the substrate. Further attempts to test the growth of JS666 on dichloroethane in this study were also negative (Table 2). Nearly identical dhlRBC clusters are located on the chromosome (Bpro_0530 to Bpro_0533) and on pPol360 (Bpro_5183 to Bpro_5186) within a duplicated 9.9-kb region that shares 97% nucleotide identity. The duplication could have been mediated by nearby transposases (Bpro_518 and Bpro_5194-Bpro_5195). Duplication of haloalkanoate dehalogenase genes may have been significant for the evolution of cDCE biodegradation in JS666.
|
View larger version (9K): [in a new window] |
FIG. 5. Organization of haloalkane and haloacid degradation genes of JS666 and X. autotrophicus strain GJ10. dhlA, haloalkane dehalogenase; dhlR, sigma-54-dependent transcriptional regulator; hyp, hypothetical protein; dhlC, putative uptake protein; dhlB, haloalkane dehalogenase; aldB, chloroacetaldehyde dehydrogenase; PQQ deH, pyrroloquinoline quinone alcohol dehydrogenase. Bpro gene numbers are shown above the JS666 genes, GenBank nucleotide accession numbers are shown above the GJ10 genes, and percentage amino acid identities to JS666 homologs are shown under the GJ10 genes.
|
![]() View larger version (24K): [in a new window] |
FIG. 6. Phylogenetic tree depicting the relationship between known dehalogenases and dehalogenase homologs in the JS666 genome. An alignment of 360 amino acids (including gaps) was made with ClustalX and converted into a neighbor-joining tree, which was visualized with TreeView using the haloalkane dehalogenase DhlA as the outgroup. The filled circles at nodes indicate bootstrap values of >95%, while the open circles at nodes indicate bootstrap values of 75% to 95%. Bootstrap values lower than 75% are not shown. The bar represents 10% sequence difference.
|
The n-alkane pathway is most likely initiated by the Bpro_5301 product, a probable cytochrome P450 with 76% amino acid identity to a Mycobacterium hydroxylase that attacks C5 to C10 alkanes (43). Downstream gene products (Bpro_5303 to Bpro_5305) are 51 to 53% similar to aldehyde dehydrogenase AlkH, alcohol dehydrogenase AlkJ, and acyl-CoA synthetase AlkK in Pseudomonas putida GPo1 (Fig. 7A). A second alcohol dehydrogenase gene of the pyrroloquinoline quinone type (Bpro_5034) is also present, suggesting the ability to metabolize both primary and secondary alcohols, as observed in Pseudomonas butanovora (48). JS666 did not use ethane, propane, or hexane for growth, but growth on heptane and octane was observed (Table 2).
![]() View larger version (27K): [in a new window] |
FIG. 7. (A) Organization of alicyclic acid biodegradation genes of JS666 and R. palustris CGA009. badH, 2-hydroxycyclohexanecarboxyl-CoA dehydrogenase; badI, 2-ketocyclohexanecarboxyl-CoA hydrolase; aliB, cyclohexanecarboxyl-CoA dehydrogenase; badK, enoyl-CoA hydratase; aliA, cyclohexanecarboxylate-CoA ligase. JS666 gene numbers (Bpro) are shown above the JS666 genes, CGA009 gene numbers (RPA) are shown above the genes, and percentage amino acid identities to JS666 homologs are shown under the CGA009 genes. (B) Organization of n-alkane biodegradation genes of JS666, P. putida GPo1, and Mycobacterium sp. strain HXN-1500. alkK, acyl-CoA synthetase; alkJ, alcohol dehydrogenase; alkH, aldehyde dehydrogenase; PQQ deH, pyrroloquinoline quinone alcohol dehydrogenase; aphG, cytochrome P450 alkane hydroxylase; aphH, ferredoxin reductase; aphI, ferredoxin. JS666 gene numbers (Bpro) are shown above the JS666 genes, GenBank nucleotide accession numbers are shown above the GPo1 and HXN-1500 genes, and percentage amino acid identities to JS666 homologs are shown under the GPo1 and HXN-1500 genes. (C) Organization of cyclic alcohol biodegradation genes of JS666 and Acinetobacter sp. strain SE19. R1, sigma-54-dependent regulator; chnC, 6-hexanolactone hydrolase; chnB, cyclohexanone monooxygenase; chnA, cyclohexanol dehydrogenase; chnD, 6-hydroxyhexanoic acid dehydrogenase; chnE, 6-oxo-hexanoate dehydrogenase; chnR, araC-like regulator; chnZ, hypothetical protein; chnY, pilin invertase; chnX, hypothetical protein. Bpro gene numbers are shown above the JS666 genes, GenBank nucleotide accession numbers are shown above the SE19 genes, and percentage amino acid identities to JS666 homologs are shown under the SE19 genes.
|
The oxidation of cyclohexanol to adipic acid is an industrially important reaction for nylon synthesis and is catalyzed by a Baeyer-Villiger monooxygenase capable of performing regio- and stereoselective oxygenations (45). Plasmid pPol338 encodes a Baeyer-Villiger monooxygenase (Bpro_5565) as part of a putative cyclohexanol biodegradation gene cluster (Fig. 2 and 7C; see Table S2 in the supplemental material) similar to that of Acinetobacter sp. strain SE19 (9). Growth experiments confirmed that strain JS666 uses cyclohexanol as a carbon and energy source (Table 2).
Aromatic metabolism.
Diverse natural and xenobiotic aromatic compounds are oxidized via the β-ketoadipate and gentisate pathways (19), and the corresponding genes are found in many bacterial genomes. Phenylpropenoid compounds, such as ferulate, are used by plants to make lignin, which is the second most abundant polymer on Earth (after cellulose). In JS666, multiple clusters of central β-ketoadipate pathway genes are present (Bpro_626 to Bpro_628, Bpro_1228 to Bpro_1331, Bpro_1855 to Bpro_1857, Bpro_1901 to Bpro_1906, and Bpro_5171 to Bpro_5173), in addition to gene clusters encoding related peripheral catabolic pathways, including ferulate (Bpro_3111 to Bpro_3118 and Bpro_3385 to Bpro_3387), benzoate (Bpro_1617 to Bpro_1633) (see Fig. S3 in the supplemental material), 4-nitrobenzoate (Bpro_682 and Bpro_5125), phthalate (Bpro_1908 to Bpro_1912), and hydroxyquinol (Bpro_2063-Bpro_2064) (Fig. 1; see Table S2 in the supplemental material).
Ferulate.
Strain JS666 possesses putative ferulate catabolic genes similar to those described in Pseudomonas sp. strains HR199 (32) and KT2440 (Fig. 1 and 8), but the ferulate operon of JS666 lacks an obvious vanillate dehydrogenase (vdh) gene and has four additional genes between aat and acd. The JS666 ferulate operon is flanked by two regulatory genes, one of which resembles FerR of KT2440. While a ferulate pathway that proceeds through vanillin has been proposed for strain HR199 (32) and also seems plausible for KT2440 (25), the JS666 ferulate pathway more likely proceeds via vanillyl-CoA due to the lack of vdh and the presence of a thioesterase gene (see Fig.S4 in the supplemental material). Growth experiments confirmed that JS666 could grow on ferulate (Table 2).
|
View larger version (12K): [in a new window] |
FIG. 8. Organization of ferulate and vanillate biodegradation genes of JS666 and P. putida KT2440. R1, GntR-type regulator; ech, enoyl-CoA hydratase; fcs, feruloyl-CoA synthetase; aat, β-ketothiolase; deH, dehydrogenase; th, thioesterase; hypo, hypothetical protein; acd, acyl-CoA dehydrogenase; R2 and ferR, marR-type regulators; vanA, vanillate demethylase oxygenase subunit; vanB, vanillate demethylase reductase subunit; vanR, gntR-like regulator; vdh, vanillate dehydrogenase. JS666 gene numbers (Bpro) are shown above the JS666 genes, KT2440 gene numbers (PP) are shown above the KT2440 genes, and percentage amino acid identities to JS666 homologs are shown under the KT2440 genes.
|
![]() View larger version (20K): [in a new window] |
FIG. 9. Organization of protocatechuate biodegradation genes of JS666, Cupriavidus metallidurans CH34, and A. baylyi ADP1. pcaQ, LysR family regulator; pcaH, protocatechuate 3,4-dioxygenase beta subunit; pcaG, protocatechuate 3,4-dioxygenase alpha subunit; pcaR, pobR family regulator; pcaI, β-ketoadipate:succinyl-CoA transferase alpha subunit; pcaJ, β-ketoadipate:succinyl-CoA transferase beta subunit; pcaF, β-ketoadipyl-CoA thiolase; pcaB, β-carboxy-cis,cis-muconate lactonizing enzyme; pcaC, -carboxymuconolactone decarboxylase; pcaD, β-ketoadipate enol-lactone hydrolase. JS666 gene numbers (Bpro) are indicated above the JS666 genes, CH34 gene numbers (Rmet) are shown above the CH34 genes, ADP1 gene numbers (ACIAD) are shown above the ADP1 genes, and percentage amino acid identities to JS666 homologs are shown under the CH34 and ADP1 genes.
|
Benzoate.
A cluster of 16 genes encodes the putative aerobic benzoate biodegradation pathway in strain JS666 (Fig. 1; see Table S1 in the supplemental material). The cluster bears similarity to the aerobic benzoate degradation genes in Azoarcus evansii strain KB 740 (see Fig. S3 in the supplemental material), as well as the anaerobic genes of R. palustris strain CGA009. Growth on benzoate was not observed in liquid media (Table 2) but was observed via auxanography (photo not shown). Total protein on benzoate-fed auxanography plates ranged from 0.92 to 0.96 mg per plate, while total protein on control plates (with no exogenous carbon source added) ranged from 0.44 to 0.70 mg per plate.
Phthalate.
A putative phthalate biodegradation gene cluster in strain JS666 consists of genes encoding a phthalate 4,5-dioxygenase (Bpro_1909 and Bpro_1910), 4,5-dihydro-4,5-dihydroxyphthalate dehydrogenase (Bpro_1908), and 4,5-dihydroxyphlatate decarboxylase (Bpro_1912). The encoded enzymes show 45 to 79% identity to those encoded on the P. putida PHT plasmid (30). The expected product of 4,5-dihydroxyphlatate decarboxylase is protocatechuate, which is expected to be further metabolized in strain JS666 (see Fig. S4 in the supplemental material). However, growth on o-phthalate was not observed (Table 2). The "phthalate" genes of JS666 could be involved in utilization of iso (meta)-phthalate or tere (para)-phthalate isomers, rather than ortho-phthalate or, alternatively, they may encode phthalate ester degradation, which requires the ester form for transport or induction.
Hydroxyquinol.
Bpro_2064 and Bpro_2063 encode a putative hydroxyquinol 1,2-dioxygenase and 3-hydroxy-cis,cis-muconate reductase, respectively (Fig. 1; see also Table S1 in the supplemental material). These encoded enzymes are expected to convert hydroxyquinol into β-ketoadipate, which is expected to be further metabolized (see Fig. S4 in the supplemental material). The possibility that JS666 uses hydroxyquinol as a growth substrate was not tested in liquid media. Auxanography experiments revealed evidence of growth on hydroxyquinol (see Fig. S6 in the supplemental material). Total protein on hydroxyquinol-fed plates ranged from 0.51 to 0.53 mg per plate, while total protein on control plates (with no exogenous carbon source added) ranged from 0.44 to 0.70 mg per plate. Therefore, growth of JS666 on hydroxyquinol remains unresolved.
Another major pathway for aromatic metabolism is the gentisate pathway, which is used when ring cleavage occurs between carboxyl- and hydroxyl-substituted carbon atoms (19). Gentisate biodegradation genes were noted in the genome (Bpro_981-Bpro_982, Bpro_991, and Bpro_2125-Bpro_2126), and growth on gentisate was confirmed (Table 2). Aromatic compounds that appear to be degraded via the gentisate pathway in strain JS666 include salicylate (Bpro_980 and Bpro_983 to Bpro_985) and 3-hydroxybenzoate (Bpro_2128 and Bpro_3591).
Salicylate and 3-hydroxybenzoate.
The chromosomally located salicylate pathway in strain JS666 (Fig. 1) bears similarity to the pathway described in naphthalene-degrading Ralstonia sp. strain U2 (51) (Fig. 10). Growth experiments in liquid media confirmed that JS666 used salicylate as a growth substrate (Table 2). The putative JS666 salicylate 5-hydroxylase subunits (Bpro_0983 and Bpro_0984) are 79% (large subunit) and 59% (small subunit) identical to the U2 salicylate 5-hydroxylase. One major difference between the two gene clusters is that a naphthalene dioxygenase is cotranscribed with the salicylate degradation genes in strain U2, while in strain JS666, there is no evidence for a naphthalene dioxygenase, but a 2-component oxygenase similar to vanillate demethylase is encoded immediately downstream (Fig. 10). One plausible hypothesis (not tested experimentally) is that the JS666 oxygenase enables oxidative demethylation of o-anisic acid (2-methoxyphenol) to yield salicylate. The observation of growth on gentisate, along with the discovery of additional putative salicylate/3-hydroxybenzoate hydroxylases in the JS666 genome (Bpro_2128 and Bpro_3591), raised the possibility that JS666 could grow on 3-hydroxybenzoate (see Fig. S4 in the supplemental material). This was confirmed experimentally (Table 2).
|
View larger version (11K): [in a new window] |
FIG. 10. Organization of salicylate/gentisate biodegradation genes in JS666 and Ralstonia sp. strain U2. R1, LysR family regulator; nagAa, ferredoxin reductase; nagI, gentisate 1,2-dioxygenase; nagK, fumarylpyruvate hydrolase; nagG, salicylate 5-hydroxylase large oxygenase component; H, salicylate 5-hydroxylase small oxygenase component (nagH); Ab, ferredoxin nagAb; nagL, maleylpyruvate isomerase; R2, TetR regulator; oxyg, oxygenase; ABC, ABC transporter; reduc, reductase; nagR, LysR family regulator; nagAc, napththalene dioxygenase large subunit; nagAd, napththalene dioxygenase small subunit. JS666 gene numbers (Bpro) are shown above the JS666 genes, GenBank nucleotide accession numbers are shown above the U2 genes, and percentage amino acid identities to JS666 homologs are shown under the U2 genes.
|
Opine metabolism.
The biosynthesis and catabolism of opines (unusual amino acid derivatives) is central to the interaction of the pathogen Agrobacterium with host plants, and opine catabolism is also found in other plant-associated bacteria. A large cluster of genes likely to encode the uptake and catabolism of opines (octopine, nopaline, and mannopine) is found on the JS666 chromosome (Fig. 1; see also Table S2 in the supplemental material). The opine genes Bpro_4331 to Bpro_4362 were located in a putative genomic island (GI8) (see Table S1B in the supplemental material) that also contains genes for amino acid metabolism and transport. The data suggest that JS666 could occupy a plant-associated niche, which would have relevance for potential applications in combined microbial/plant bioremediation.
|
|
|---|
Catabolic pathways were assigned only in cases where whole clusters of catabolic genes were present (not single genes) and where good BLAST matches for the majority of such genes were obtained. In several cases, growth substrates predicted by genomic analyses were confirmed experimentally (Table 2). An auxanography approach confirmed that JS666 uses catechol and benzoate as growth substrates, even though liquid growth experiments were negative. Despite extensive washing and dilutions, we observed apparent pinpoint colonies on control plates (i.e., JS666 plated on MSM-Noble agar with no exogenous carbon source). This suggests that JS666 grows on trace constituents in the solid media. However, differential growth on plates amended with a substrate in comparison to the no-carbon control plates was confirmed by protein analysis.
Our inability to observe growth on other predicted substrates could indicate that the actual substrates are substituted variants or that the conditions used for testing were not optimal (e.g., toxic substrates that support growth only in a narrow concentration range). We have evidence for this in the case of salicylate, which initially tested negative at 2 mM but was positive at 0.6 mM. In cases where we were unable to demonstrate growth on the predicted compounds experimentally, likely catabolic pathways were still presented to provide a basis for future work, which could, for example, aim to test close structural analogs of the compounds we have experimentally tested (e.g., chlorobenzoate instead of benzoate) or substrates that were unavailable due to cost/legal restrictions (e.g., phenylacetate). In addition, identification of putative catabolic genes and pathways is important to provide targets for heterologous expression and protein purification because such strategies do not depend on whole pathways being functional, only on the activities of single or few gene products.
The replication and transfer genes of the JS666 plasmids are not closely related to those of experimentally studied plasmids, but some similarities are apparent. The presence of an IncP
-like replication determinant (trfA) linked to haloacid dehalogenase genes in pPol360 highlights the major role played by IncP plasmids in bacterial adaptation to pollutants (7, 13, 42). The similarity of pPol338 replication genes to those of Burkholderia small chromosomes indicates that gene transfer between Polaromonas and Burkholderia is likely, which is of interest due to the importance of members of the genus Burkholderia as plant and human pathogens and agents of bioremediation (8). The presence of extensive and well-conserved plasmid transfer and mobilization functions on both pPol338 and pPol360 is notable, suggesting both elements are independently capable of conjugative transfer. Establishing the host range and transfer properties of the JS666 plasmids is significant for bioremediation due to their likely roles in moving catabolic genes among bacteria in situ. The biology of the JS666 plasmids is also of interest for the development of new cloning vectors for betaproteobacteria.
Transposable elements play a major role in biodegradation pathway evolution and are often found associated with catabolic plasmids (29, 41, 49). While we did not attempt to define the inverted repeats associated with ISs and transposons, abundant and diverse transposable elements are clearly present in the JS666 genome, especially on the two plasmids. The duplication of dehalogenase genes on pPol360 and the chromosome suggests a recent transposase-mediated event, which may have been involved in the evolution of halogenated hydrocarbon metabolism (e.g., cDCE). A large region of pPol360 (76 kb) containing the β-ketoadipate, hydroxylaminobenzoate, and haloalkanoate metabolic genes is flanked by identical IS3 family transposases (Bpro_5122 and Bpro_5204) and could represent a compound (type I) catabolic transposon. Transposase genes are found throughout the alkane/cycloalkane/cyclic-alcohol gene region of pPol338, suggesting that ISs have also contributed to the evolution of aliphatic hydrocarbon catabolic abilities in JS666.
Based on our genome analysis, we propose that the acquisition of plasmids pPol338 and pPol360 played a role in the evolution of the cDCE biodegradation pathway in JS666. These large plasmids have persisted in JS666 despite an extended enrichment and isolation process (10), providing circumstantial evidence that they contribute to growth on cDCE. The presence of an alkane monooxygenase gene on pPol338 is consistent with the prior observation that cDCE-grown JS666 cells produce epoxyethane from ethene (10). A plausible cDCE assimilation pathway could arise from the intersection of oxygenase and dehydrogenase genes from alkane pathways (pPol338) and hydrolytic dehalogenase genes from haloalkanoate pathways (pPol360 and chromosomal). We have observed that the cDCE-degrading phenotype is unstable in media lacking cDCE (T. E. Mattes and N. V. Coleman, unpublished data), which would be consistent with plasmid- or transposon-borne catabolic genes, but this observation could alternatively reflect poor regulatory control over a newly evolved metabolic pathway. Another possibility is that gene amplification processes (36) contribute to the unstable phenotype.
The JS666 genome encodes many enzymes that are potentially useful for biocatalysis, including the cytochrome P450 monooxygenase (Bpro_5299 to Bpro_5301) (44), the ferulate/vanillate system (Bpro_3111 to Bpro_3118 and Bpro_3385 to Bpro_3387) (17, 35), and the cyclohexanone monooxygenase-like gene (Bpro_5565) (2). Metabolic engineering of the alkane, ferulate, and cyclohexanol biodegradation pathways of JS666 could theoretically be used to produce the fragrances limonene and vanillin and the nylon monomer adipic acid, respectively. Other JS666 genes of strong biotechnological interest not discussed in detail here include a polyketide synthase region on pPol338, which is similar to cyanobacterial toxin gene clusters; several P450 oxygenases (Bpro_1581, Bpro_4920, and Bpro_2807); a halogenase (Bpro_3311); an amidase (Bpro_571); and a nitrilase (Bpro_3376).
We thank Andy Holmes and Erich Lanka for helpful discussions.
Published ahead of print on 22 August 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»