Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2006, p. 3504-3514, Vol. 72, No. 5
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.5.3504-3514.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Environmental Microbiology,1 Division of Structural Biology, GBFGerman Research Centre for Biotechnology, Mascheroder Weg 1, D-38124 Braunschweig, Germany2
Received 9 November 2005/ Accepted 13 March 2006
|
|
|---|
subunits were targeted, since they are responsible for the substrate specificities of these enzymes. SSCP functional genotype fingerprinting revealed a substantial diversity of oxygenase genes in three differently BTEX-contaminated soil samples, and sequence analysis indicated that in both the soil DNA and the bacterial isolates, genes for oxygenases related to the isopropylbenzene (cumene) dioxygenase branch of the toluene/biphenyl oxygenase subfamily were predominant among the detectable genotypes. The peptide sequences of the two most abundant
subunit sequence types differed by only five amino acids (residues 258, 286, 288, 289, and 321 according to numbering in cumene dioxygenase
subunit CumA1 of Pseudomonas fluorescens IP01). However, a strong correlation between sequence type and substrate utilization pattern was observed in isolates harboring these genes. Two of these residues were located at positions contributing, according to the resolved crystal structure of cumene dioxygenase from Pseudomonas fluorescens IP01, to the inner surface of the substrate-binding pocket. Isolates containing an
subunit with isoleucine and leucine at positions 288 and 321, respectively, were capable of degrading benzene and toluene, whereas isolates containing two methionine substitutions were found to be incapable of degrading toluene, indicating that the more bulky methionine residues significantly narrowed the available space within the substrate-binding pocket. |
|
|---|
The aerobic degradation of aromatic compounds is frequently initiated by Rieske nonheme iron oxygenases, which catalyze the incorporation of two oxygen atoms into the aromatic ring to form arene cis-diols (30), followed by a dehydrogenation reaction catalyzed by a cis-dihydrodiol dehydrogenase to give catechol or substituted catechols which serve as substrates for oxygenolytic aromatic ring cleavage. Rieske nonheme iron oxygenases are multicomponent enzyme complexes composed of a terminal oxygenase component (iron-sulfur protein [ISP]) and different electron transport proteins (a ferredoxin and a reductase or a combined ferredoxin-NADH-reductase) (13). The catalytic iron-sulfur proteins are heteromultimers, comprising a large (
) and a small (ß) subunit, with the former containing a Rieske-type [2Fe-2S] cluster, a mononuclear nonheme iron oxygen activation center, and a substrate-binding site (13, 27) and being responsible for substrate specificity (30). Comparison of the amino acid sequences of the terminal oxygenase
subunits (ISP
) revealed that they form a family of diverse but evolutionarily related sequences, and four distinct major lineages have been identified (30). Although none of the enzymes is completely specific, a broad correlation between the grouping in toluene/biphenyl, naphthalene, benzoate, or phthalate subfamilies and the native substrates oxidized by the subfamily members can be observed. However, various enzyme-engineering studies on biphenyl, benzene, chlorobenzene, and naphthalene dioxygenases also showed that only slight sequence differences in the amino acid sequences of ISP
can be associated with dramatic changes in substrate specificity or regiospecificity (8, 24, 50, 69). Information on Rieske nonheme iron oxygenases has been obtained mainly from cultured strains. However, various culture-independent studies suggest that distinct and potentially numerically dominant sequence types of as-yet-unrecognized functional genes exist in natural and anthropogenically influenced environments and that the sequences from cultured strains are not likely to represent the microbial functional gene diversity in the environment (63, 64, 68).
Members of the toluene/biphenyl oxygenase subfamily have been reported to be involved in the degradation of benzene, chlorobenzenes, toluene, ethylbenzene, isopropylbenzene, and biphenyl, and phylogenetic analysis of their ISP
peptide sequences (47) indicated the presence of several distinct sequence clusters, such as benzene/toluene oxygenases, isopropylbenzene oxygenases, and biphenyl oxygenases. The amino acid residues reported to be responsible for the substrate specificity and regioselectivity of members of the subfamily are located within a central core of the
subunit (20, 27, 69). ISP
DNA sequence alignments had identified highly conserved regions for use as primer binding sites in PCR, including conserved regions flanking the central core, thus facilitating molecular studies on the diversity of these catabolic genes in the environment (38, 64, 68). Detection and differentiation of toluene/biphenyl oxygenases by using PCR-based genetic profiling techniques targeting regions responsible for substrate specificity could provide hints on the overall metabolic network acting cooperatively for in situ biodegradation. Monitoring the structure and diversity of predominant catabolic gene polymorphisms can give indications on functional adaptations, as optimized catabolic enzyme variants are likely to be selected for functioning under the existing environmental conditions.
PCR-single-strand conformation polymorphism (SSCP) analysis has been demonstrated to be a powerful technique for elucidating the structure of complex microbial communities as well as the diversity of catabolic genes (37, 52-54). Based on sequence-dependent folding of single-stranded DNA (ssDNA) into secondary-structure conformations, which exhibit different electrophoretic mobilities in nondenaturing polyacrylamide gels, DNA molecules of the same length but with different nucleotide sequences can be separated (54). Here we describe the application of the PCR-SSCP technique for the simultaneous detection and discrimination of ISP
gene segments retrieved from pure culture strains exhibiting different aromatic hydrocarbon-degradative capabilities and from DNA extracted from soil samples differently contaminated with benzene, toluene, ethylbenzene, and xylenes (BTEX). The translation products of ISP
gene segments recovered from SSCP fingerprints were analyzed on the basis of structural information of related dioxygenase enzymes in order to evaluate the functional significance of structurally important amino acid residues.
|
|
|---|
gene segment differing from that reported previously (GenBank accession no. AJ010057) was detected (see Results). The nearly complete 16S rRNA gene sequence (1,455 nucleotides) obtained showed 100% sequence identity with 16S rRNA genes from Pandoraea spp. but differed from the sequence previously reported for strain JB1 (X92188) by seven nucleotides. Additionally, 36 bacterial isolates obtained from BTEX-contaminated soil samples (37), the majority of which were capable of mineralizing benzene, were analyzed for the presence of ISP
genes of the toluene/biphenyl subfamily. |
View this table: [in a new window] |
TABLE 1. Reference strains used in this study and results of the PCR targeting ISP genes of the toluene/biphenyl oxygenase subfamily
|
Chemicals.
Chemicals were purchased from Aldrich Chemie (Steinheim, Germany), Fluka AG (Buchs, Switzerland), or Merck AG (Darmstadt, Germany). Benzene-cis-dihydrodiol (cis-3,5-cyclohexadiene-1,2-diol) was obtained from Sigma-Aldrich Chemie (Deisenhofen, Germany).
Oxygen uptake measurements with whole cells.
Transformation of benzene, toluene, phenol, and benzene-cis-dihydrodiol was determined polarographically by measuring O2 consumption with a Clark-type oxygen electrode (Bachhofer, Reutlingen, Germany). Benzene-grown cells were harvested during late-exponential-phase growth, washed, and resuspended to a turbidity (A546) of 30 in 50 mM phosphate buffer (pH 7.5). Of this suspension, 0.1 ml was added to 2.9 ml phosphate buffer (50 mM, pH 7.5) saturated with air. After 5 min of constant endogenous oxygen uptake, the reaction was started by adding the substrate (100 mM aqueous stock solutions or solutions in dimethyl sulfoxide in the cases of benzene and toluene) to a final concentration of 0.2 mM. Uptake rates were determined at 25°C and corrected for endogenous O2 uptake. In the case of benzene-cis-dihydrodiol transformation, metabolism was confirmed by high-performance liquid chromatography analysis. Samples were taken during oxygen uptake measurements, and 10 µl of cell-free supernatants was analyzed with a Shimadzu high-performance liquid chromatography system (LC-10AD liquid chromatograph, DGU-3A degasser, SPD-M10A diode array detector, and FCV-10AL solvent mixer) equipped with an SC125/LiChrospher column (Bischoff, Leonberg, Germany) by use of an aqueous solvent system (flow rate, 1 ml/min) containing 0.01% (vol/vol) H3PO4 (87%) and 30% (vol/vol) methanol.
DNA extraction.
DNA from eight reference strains (Burkholderia xenovorans LB400, Cupriavidus necator H850, Pseudomonas sp. strain Cam-1, Pseudomonas pseudoalcaligenes KF707, Pandoraea sp. strain JB1*, Pseudomonas putida F1, Ralstonia sp. strain PS12, and Rhodococcus globerulus P6) was extracted using the QIAGEN Genomic DNA Tip 100/G system (QIAGEN, Hilden, Germany). DNA concentrations in the extracts were quantified using a PicoGreen double-stranded DNA (dsDNA) quantitation kit (Molecular Probes, Leiden, The Netherlands) and adjusted to a final concentration of 1 ng/µl with Tris-HCl buffer (10 mM, pH 8.0). To extract DNA from other reference strains and from the BTEX isolates, single colonies were boiled and centrifuged (39) and the supernatants were used either 10-fold or 100-fold diluted in Tris-HCl buffer (10 mM, pH 8.0) as template DNA in PCR. Total DNA from soil samples (wet weight, 500 mg) was extracted using a FastDNA spin kit for soil (Bio101 Systems, Q-BIOgene, Heidelberg, Germany) and DNA concentrations in extracts were determined as described above. The soil extracts, containing approximately 2 ng/µl DNA, were used either directly or 10-fold diluted in Tris-HCl buffer (10 mM, pH 8.0) as template DNA in PCR.
Evaluation of primer binding sites suitable for amplification of ISP
gene segments.
DNA and amino acid sequences of the
subunits of Rieske nonheme iron oxygenases (ISP
) were retrieved from GenBank (http://www.ncbi.nlm.nih.gov/) and aligned with ClustalW (http://www.ebi.ac.uk/clustalw/). A subset of 31 ISP
genes (GenBank accession numbers AF006691, AJ293587, U53507, D37828, AF049851, AJ010057, U47637, AB086835, M86348, U95054, M83673, AY027651, E04215, D88020, X80041, U27591, D32142, U24277, AF148496, AJ006307, J04996, U15298, U78099, D17319, U51165, AF452376, AB121977, AY554272, AJ536756, AY831463, and AB193045), comprising the majority of currently known members of the toluene/biphenyl subfamily (30), was examined for conserved regions for use as PCR primer binding sites. Amino acid residues which were reported to have an influence on substrate specificity of oxygenases of the toluene/biphenyl subfamily are encoded by ISP
gene segments flanked by conserved sites at approximately positions 670 and 1170 of the bphA1 gene from B. xenovorans LB400 (8, 26, 46, 60, 69). Two degenerate primer sets targeting these conserved regions were used to amplify 525- to 531-bp segments of ISP
genes. The primer set BPDOXF/BPDOXR was designed by Yeates et al. (68); however, DNA sequence alignments revealed a critical Y-A mismatch between reverse primer BPDOXR and the ISP
genes bphA1 from R. globerulus P6 and Rhodococcus erythropolis TA421 and bpdC1 from Rhodococcus sp. M5 (GenBank accession numbers X80041, D88020, and U27591, respectively) at position 3 from the primer 3' end. To increase their specificity for amplification of ISP
genes of the toluene/biphenyl subfamily, these primers were modified. The modified forward and reverse primers (bphAf668-3 [5'-GT TCC GTG TAA CTG GAA RTW YGC-3'] and bphAr1153-2 [5'-CCA GTT CTC GCC RTC RTC YTG HTC-3'], respectively) comprised a pool of degenerate primers in which a minimum of 14 bases at the 3' end corresponded to the aligned toluene/biphenyl ISP
gene sequences. For PCR-SSCP analysis, ISP
gene segments were amplified with the modified primer set by using a 5'-end-phosphorylated reverse primer (54). The forward primer, bphAf371B (5'-GGC TTT CAC CTG CAS YTA YCA YGG-3') (targeting positions 348 to 371 in bphA1 of LB400) was used together with bphAr1153-2 to amplify 799- to 826-bp ISP
segments, which comprised the primer binding site for bphAf668-3.
PCR amplification of ISP
gene segments.
PCR amplification was performed in a final volume of 50 µl containing 1x PCR buffer (QIAGEN, Hilden, Germany) supplemented with 1.5 mM MgCl2, 200 µM of each dNTP (Bioline, Luckenwalde, Germany), 0.5 µM of each primer (MWG-Biotech, Ebersberg, Germany), 2.5 U of HotStarTaq polymerase (QIAGEN), and 2 µl of template DNA. The PCR program comprised an initial denaturation for 15 min at 95°C followed by 32 cycles of 40 s at 94°C, 40 s at 58°C, and 60 s at 72°C and a final elongation for 10 min at 72°C. For amplification of soil DNA, the cycle number was increased to 35 cycles (sample 1Y or 3Y) or 40 cycles (sample 1X). The PCR conditions for amplification with BPDOXF and BPDOXR were as described previously (68). PCR products (2 µl of the PCR mixture) were separated by agarose gel electrophoresis (1.5% agarose, 1x Tris-acetate-EDTA buffer; 1 h at 90 V) and visualized by ethidium bromide staining. PCR yield estimates were obtained by quantitative image analysis of digitized agarose gel images using Scion Image Beta version 4.0.2 (Scion Corporation, Frederick, MD).
PCR-SSCP analysis of ISP
gene segments.
PCR products were purified with a MinElute PCR purification kit (QIAGEN) as recommended by the manufacturer and eluted with 20 µl Tris-HCl buffer (10 mM, pH 8.0). The phosphorylated strands were removed by digestion with 12.5 U lambda exonuclease (New England Biolabs, Beverly, MA) in a final volume of 25 µl for 1 h at 37°C (54). The digestion was stopped by the addition of 5 µl sodium acetate (3 M, pH 5.5), and the ssDNA was purified with a MinElute PCR purification kit and eluted with 20 µl sterile water. The ssDNA concentration was assessed by comparing the fluorescence intensities of ssDNA sample bands to the fluorescence intensities of dsDNA standard bands (Low DNA Mass Ladder; Invitrogen) after agarose gel electrophoresis (1.5% agarose, 1x Tris-acetate-EDTA buffer; 1 h at 90 V) and ethidium bromide staining. Prior to the SSCP procedure, samples were diluted with sterile water to yield an ssDNA concentration equivalent to 1.7 ng/µl dsDNA per sample. PCR-negative controls were used without further dilution. Aliquots of 17 µl were evaporated for 15 min at 30°C in a vacuum concentrator to yield a final volume of approximately 4 µl. The samples were mixed with 4 µl SSCP loading buffer (95% [vol/vol] formamide, 10 mM NaOH, 0.025% [wt/vol] bromophenol blue, and 0.025% [wt/vol] xylene cyanol), incubated for 2 min at 95°C in a water bath, and placed on ice water for 5 min. The influence of gel temperature and MDE gel solution (Cambrex Bioscience, Rockland, ME) concentration on the resolution performance of the SSCP method was tested with ssDNA ISP
segments generated from six reference strains (LB400, H850, Cam-1, KF707, JB1*, F1, PS12, and P6) in individual runs using the following conditions (gel temperature/MDE gel solution concentration): (i) 20°C/0.6x, (ii) 24°C/0.6x, (iii) 28°C/0.6x, (iv) 28°C/0.7x, and (v) 30°C/0.7x. Electrophoresis under optimized conditions was performed on 21- by 21-cm 0.7x MDE gels in 1x Tris-borate-EDTA as the running buffer at 400 V and 10 mA for 20 h at 30°C in a Pharmacia Macrophor electrophoresis unit 2010-001 connected to a circulating water bath (Lauda Ecoline RE104 thermostat). The gels were silver stained according to the procedure of Bassam et al. (5) and dried at room temperature. ISP
SSCP fingerprints from pure cultures were reproduced in at least two independent SSCP analyses. Silver-stained gels were scanned on a flatbed scanner (Epson Expression 1600 PRO) to create 8-bit TIFF images, which were then analyzed with GelCompar II software (Applied Maths, Kortrijk, Belgium).
Isolation and reamplification of single-strand products from SSCP profiles.
Selected single-strand products were excised from silver-stained polyacrylamide gels and transferred to microtubes containing 50 µl of elution buffer (10 mM Tris-HCl [pH 8.8], 50 mM KCl, 1.5 mM MgCl2, 0.1% Triton X-100). The tubes were incubated for 20 min at 95°C, and 3 µl of eluted ssDNA was subjected to PCR reamplification using the same primers and conditions used to generate the original dsDNA segments.
Cloning of PCR-amplified ISP
segments.
PCR-reamplified products from selected SSCP bands and ISP
gene segments amplified with the primer set bphAf371B/bphAr1153-2 were purified with a QIAquick PCR purification kit (QIAGEN), ligated into pGEM-T-vector system (Promega), and transformed into Escherichia coli JM109 competent cells (Promega) as specified by the manufacturers. Plasmid inserts were amplified by PCR with vector-specific primers M13 forward and M13 reverse (Promega), annealing in the flanking cloning region.
DNA sequencing and sequence analyses.
PCR products were purified with a QIAquick PCR purification kit and sequenced using an ABI PRISM BigDye Terminator v1.1 Ready Reaction cycle sequencing kit (Applied Biosystems, Darmstadt, Germany). Primers used for sequencing reactions were the same as those used in the original PCR, and sequencing was performed on an ABI PRISM 3100 genetic analyzer (Applied Biosystems). Raw sequence data from both strands were assembled with Sequencher software version 4.0.5 (Gene Codes Corporation, Ann Arbor, MI), and consensus sequences were translated in frame (omitting the primer sequences) with BioEdit software version 5.0.9 (http://www.microbiology.ncsu.edu/resources/index.html). For cloned PCR products, only consensus sequences deduced from at least three identical or nearly identical clones (difference of up to four bases) were included in the sequence analysis in order to reduce potential sequence diversity biases due to PCR and/or sequencing errors. DNA similarity searches were performed with the BLASTN algorithm and public nucleotide databases (http://www.ncbi.nlm.nih.gov/). DNA and protein alignments were conducted with ClustalW. Phylogenetic analyses were done with MEGA software version 3.0 (http://www.megasoftware.net/), by using the neighbor-joining method with the p-distance model and pairwise deletion of gaps/missing data. A consensus tree was inferred from a total of 1,000 bootstrap trees generated for each data set.
AFDRA and ARDRA.
The presence of a catechol 2,3-dioxygenase (C23O) gene was determined by using PCR amplification with primers C23O-ORF-F and C23O-ORF-R, and C23O types were identified by amplified functional DNA restriction (AFDRA) analysis as reported previously (36). Amplified ribosomal DNA restriction analysis (ARDRA) was performed by using protocols described previously (37).
Nucleotide sequence accession numbers.
The ISP
sequences and the 16S rRNA gene sequence of Pandoraea sp. strain JB1* are available under GenBank accession numbers DQ166964 to DQ167022 and DQ336933 to DQ336947.
|
|
|---|
genes. Using this primer set, single products of the expected size were amplified from 11 reference strains known to possess an ISP
gene of the toluene/biphenyl oxygenase subfamily (Table 1). No products were amplified from reference strains containing
subunits of the naphthalene (P. putida G7 [58] and P. stutzeri AN10 [11]) or benzoate/toluate (P. putida mt-2 [67] and MT53 [40]) oxygenase subfamilies and other reference strains obviously devoid of genes encoding Rieske nonheme iron oxygenase
subunits of the toluene/biphenyl subfamily (Table 1). However, compared to the other positively tested target genes, amplification of bphA1 from R. globerulus P6 resulted in a considerably lower product yield, as judged from quantitative image analysis of ethidium bromide-stained agarose gels (data not shown). To increase the specificity of the primers for amplification of ISP
genes of the toluene/biphenyl subfamily, the primers were modified to match with an updated DNA alignment comprising 31 reference sequences. PCR performed with the modified primers resulted in similar product yields for all reference strains containing target genes.
SSCP analysis of ISP
gene segments.
A PCR-SSCP method was established for the detection and discrimination of ISP
genes of the toluene/biphenyl oxygenase subfamily using
subunit gene segments amplified from reference strains with the primer set bphAf668-3/bphAr1153-2. Five different combinations of gel temperature and MDE gel solution concentration were tested to determine the optimum conditions for differentiation of the approximately 530-base single-strand ISP
segments, and best resolution and selectivity of SSCP profiles were achieved using an MDE gel concentration of 0.7x and a gel temperature of 30°C, allowing the discrimination of the predominant bands of 11 reference ssDNA products (Fig. 1). Single-strand ISP
gene segments of identical mobilities were obtained from B. xenovorans LB400 and C. necator H850, which was expected, as both strains were reported to contain ISP
genes identical in sequence in the fragment analyzed here (38). The developed PCR-SSCP method was able to separate even highly homologous ISP
gene segments, as indicated by the different migration patterns of the single strands originating from Pseudomonas sp. strain Cam-1 and P. pseudoalcaligenes KF707, which differed by only one nucleotide. Sequences retrieved from predominant bands were identical to those previously reported, except that of the Pandoraea sp. strain JB1* bphA1 sequence, which differed in seven nucleotides from the sequence reported for Burkholderia sp. strain JB1 (AJ010057). The ISP
gene fragment recovered from the PCR-SSCP profile of the biphenyl degrader Pseudomonas sp. strain IC (14) was found to be identical to bphA1 from Pseudomonas sp. strain B4 (U95054). SSCP analysis of the majority of reference strains resulted in the identification of one ISP
gene only; however, a second ISP
gene could be observed with C. necator H850 (Fig. 1, lane 3, lower band), the deduced amino acid sequence of which exhibited the highest levels of similarity to putative ISP
proteins from Bordetella bronchiseptica strain RB50 (BX640439) and Variovorax paradoxus strain E4C (AF209470) (84% to 86% amino acid sequence identity).
![]() View larger version (16K): [in a new window] |
FIG. 1. PCR-SSCP analysis of ISP gene segments obtained from reference strains. The following strains were analyzed: P. putida F1 (todC1, lane 1); B. xenovorans LB400 (bphA1, lane 2); C. necator H850 (bphA1, lane 3); Pseudomonas sp. strain IC (bphA1, lane 4); R. globerulus P6 (bphA1, lane 5); P. pseudoalcaligenes KF707 (bphA1, lane 6); P. aeruginosa JI104 (bnzA, lane 7); Ralstonia sp. strain PS12 (tecA1, lane 8); Pseudomonas sp. strain Cam-1 (bphA1, lane 9); Pandoraea sp. strain JB1* (bphA1, lane 10); and P. fluorescens IP01 (cumA1, lane 11). Lane 12 shows the PCR-negative control. Molecular weight marker III (M) (250 ng per lane; Roche) was used to normalize the band patterns.
|
genes in isolates from BTEX-contaminated soil samples.
genes of the toluene/biphenyl oxygenase subfamily using primers bphAf668-3 and bphAr1153-2, and products of the expected size could be amplified with DNA extracted from 33 out of the 36 bacterial isolates (Table 2). ISP
gene segments could not be amplified from the isolates IA1YICDA and IA1YICDB, which were capable of mineralizing toluene, as well as m- and p-xylene, but not benzene (37). This phenotype is similar to that of strains harboring the well-studied toluene degradation pathway encoded by TOL plasmids (12, 67), where degradation is initiated by oxidation of the methyl group. Additionally, isolate 3YC1b failed to give a PCR product. Its capability to grow on benzene, however, suggests that it either harbors an ISP
gene exhibiting low gene sequence homology with the primers used here or possibly initiates the degradation by means of two successive monooxygenation reactions. |
View this table: [in a new window] |
TABLE 2. Substrate utilization profiles and genotypic characteristics of BTEX isolatesa
|
gene segments revealed the presence of only one predominant single-strand product per isolate and allowed five different electrophoretic mobilities to be distinguished (Fig. 2A, band positions A, B, C, D, and E). The nucleotide sequence of the most frequently occurring band, designated ISP
type A, was, as verified by sequencing, identical in 17 out of the 33 positively tested isolates, and the derived translation product showed highest identity to the ISP
peptide (IpbAa) of P. putida RE204 isopropylbenzene dioxygenase (22) (Fig. 3). The second predominant ISP
type, referred to as ISP
type C, was observed in 13 isolates. All but one of these isolates contained an ISP
gene sequence that differed from sequence type A by only six nucleotides, corresponding to a difference of five amino acids in the deduced peptide sequence. The ISP
sequence retrieved from isolate 1XB1, designated ISP
type C**, differed from sequence type C by a single nucleotide, resulting in one amino acid difference. Three further ISP
sequence types, detectable in one isolate each, were found to be either highly homologous to ISP
types A and C (sequence type E) or more closely related to the ISP
sequence cluster encoding biphenyl dioxygenases of gram-positive strains and (chloro-)benzene/toluene dioxygenases (ISP
types B and D) (Fig. 3). ISP
PCR-SSCP analysis performed on artificial mixtures of genomic DNA extracted from three strains harboring different ISP
types showed that each component was amplified by PCR, and product yields reflected the abundance of each ISP
gene in the template mixture (see Fig. S1 in the supplemental material).
![]() View larger version (44K): [in a new window] |
FIG. 2. SSCP profiles of ISP gene segments obtained from BTEX-degrading isolates (A) and BTEX-contaminated soil samples (B). Bands selected for reamplification (by PCR), cloning, and sequencing are indicated by letters. Molecular weight marker III (Roche) (250 ng per lane) was used to normalize the band patterns. ISP types were named according to their band positions as types A through E for the isolate profiles and types A through P for the soil profiles. NTC, PCR nontemplate control.
|
![]() View larger version (27K): [in a new window] |
FIG. 3. Phylogenetic analysis of deduced peptide sequences from ISP gene segments recovered from SSCP profiles (comprising 158 to 161 amino acids, corresponding to positions 224 to 384 of BphA1 from B. xenovorans LB400). Characterized ISP peptides of the toluene/biphenyl oxygenase subfamily (IpbAa, Pseudomonas putida RE204 [AF006691]; EbdAa, Pseudomonas putida 01G3 [AJ293587]; IpbA1, Pseudomonas sp. strain JR1 [U53507]; CumA1, Pseudomonas fluorescens IP01 [D37828]; EdoA1, Pseudomonas fluorescens CA-4 [AF049851]; BphA1, Burkholderia sp. strain JB1 [AJ010057]; BphA1, Comamonas testosteroni B-356 [U47637]; BphA1, Comamonas testosteroni TK102 [AB086835]; BphA1, Pseudomonas sp. strain KKS102 [D17319]; BphA1, Burkholderia xenovorans LB400 [M86348]; BphA1, Pseudomonas sp. strain B4 [U95054]; BphA1, Pseudomonas pseudoalcaligenes KF707 [M83673]; BphA, Pseudomonas sp. strain Cam-1 [AY027651]; BnzA, Pseudomonas aeruginosa JI104 [E04215]; BphA1, Rhodococcus erythropolis TA421 [D88020]; BphA1, Rhodococcus globerulus P6 [X80041]; BpdC1, Rhodococcus sp. strain M5 [U27591]; BphA1, Rhodococcus sp. strain RHA1 [D32142]; IpbA1, Rhodococcus erythropolis BD2 [U24277]; BedC1, Pseudomonas putida ML2 [AF148496]; McbAa, Ralstonia sp. strain JS705 [AJ006307]; TodC1, Pseudomonas putida F1 [J04996]; TcbAa, Pseudomonas sp. strain P51 [U15298]; TecA1, Ralstonia sp. strain PS12 [U78099]; ISP , Variovorax paradoxus E4C [AF209469]; HcaE, Bordetella bronchiseptica RB50 [NC_002927]; ISP , Cupriavidus necator H850 [2nd SSCP band]) and ISP peptides from bacterial isolates (bold italics) and soil samples (bold) are shown. The ISP peptides obtained from soil (soil-) were designated according to the procedure of sequence recovery, i.e., direct sequencing of reamplified SSCP products (d) or sequencing of cloned PCR-amplified products (c), followed by band position (A through P) and sample origin (1X, 1Y, or 3Y). In cases where more than one consensus sequence was retrieved by cloning, the different sequence types are indicated by lowercase letters (a, b, or c) after the band position designation. The ISP peptide sequences from isolates (iso-) were deduced from directly sequenced reamplification products. The tree was constructed by using the neighbor-joining method in MEGA with the p-distance model and pairwise deletion of gaps/missing data. The scale bar indicates the percent divergence. The names of ISP peptide clusters were assigned according to the native substrates used for isolation of their hosts and characterization of their activity. pos., positive; neg., negative.
|
type and the substrate utilization profile observed with the respective isolates (Table 2), as all isolates harboring ISP
type A could grow on benzene and toluene as a sole source of carbon and energy, while all but 1 of the 13 isolates containing ISP
type C were able to grow on benzene but not on toluene. To provide evidence that degradation of benzene by the isolates is initiated by dioxygenation rather than two successive monooxygenations, oxygen uptake rates by two representative isolates containing ISP
type A and two isolates containing ISP
type C as well as that by the isolate containing ISP
type B were assessed (Table 3). All isolates exhibited a significant oxygen uptake with benzene cis-dihydrodiol, indicating that this compound is dehydrogenated and subsequently subject to ring cleavage, whereas oxygen uptake with phenol was absent. While all tested isolates showed significant oxygen uptake with benzene, which had been used as growth substrate, only the two isolates harboring ISP
type A were active with toluene. |
View this table: [in a new window] |
TABLE 3. Oxygen uptake rates of selected benzene-degrading isolates
|
genes upstream of the region used for PCR-SSCP, the isolates were subjected to a second PCR analysis using the primer set bphAf371B/bphAr1153-2. PCR amplification of ISP
gene segments was observed with the same 33 isolates described above (data not shown), and sequencing of the cloned products from seven strains harboring ISP
type A, four strains harboring ISP
type C, and the one strain harboring ISP
type C** (Table 2) revealed that they were identical in the 300-bp DNA stretch upstream of the sequence derived from SSCP analysis, indicating that sequence differences were confined to a region of the carboxy-terminal portion of the
subunits that was previously identified to be of crucial importance for determining substrate specificity (8, 24, 26, 60, 69).
ISP
gene diversity in BTEX-contaminated soil samples.
The developed PCR-SSCP technique was used to assess the diversity and distribution of Rieske nonheme iron oxygenases of the toluene/biphenyl subfamily in differently contaminated soil samples from a BTEX-contaminated site. Comparison of the SSCP profiles revealed significant differences in the catabolic gene structure between the three soil samples (Fig. 2B). Bands that were detectable against a low surrounding background could be sequenced directly from reamplified ssDNA, resulting in unambiguous identification of the single-strand products. The majority of bands with high relative intensities and identical mobilities, like bands A and C, were found to contain identical single-stranded DNA species. Several of the closely migrating single-strand products represented identical ssDNA sequence types (bands G1/G2, H1/H2, K1/K2, and L1/L2), suggesting that the same single strand formed two conformers migrating with slightly different mobilities. In the case of faint bands, direct sequencing of reamplified single-strand products did not result in conclusive sequences. To reveal the identity of the ISP
types contained in these bands, the reamplification products were cloned and sequenced. For each band, one (or two distinct) consensus sequence(s) could be inferred from at least three clones containing nearly or completely identical sequence types.
Phylogenetic analysis of the peptide sequences deduced from directly sequenced PCR-amplified SSCP products revealed the majority of the soil peptides to be related to ISP
sequences of ethylbenzene/isopropylbenzene dioxygenases (Fig. 3), with the nucleotide sequences retrieved from bands A, C, and E being identical to the corresponding ISP
segments retrieved from isolates. The ISP
gene sequences retrieved from bands G1/G2 and L1/L2, respectively, differed by nine nucleotides; however, none of them caused a difference in the deduced peptide sequences, which were most closely related to TodC1 from the toluene-degrading bacterium Thauera sp. strain DNT-1 (57). The ISP
sequence types soil-dK, soil-dP, and soil-dD, which were, among the soil samples tested, detectable only in the highly contaminated soil samples 1X and 1Y, grouped together in a divergent branch of the toluene/biphenyl oxygenase subfamily and exhibited less than 75% identity to previously described ISP
peptide sequences. The peptide sequence inferred from an intense band, which was detected only in the SSCP profile of soil 3Y, was found to be identical to benzene dioxygenase ISP
from P. putida ML2 (62). ISP
sequences closely related to P. pseudoalcaligenes KF707 (61) and Pseudomonas sp. strain Cam-1 (45) biphenyl dioxygenase ISP
peptides were inferred from bands H1/H2 and J.
|
|
|---|
gene homologues amplified from bacterial isolates and environmental samples. Several primer sets suitable for amplification of Rieske nonheme iron oxygenase ISP
gene segments have been reported previously (4, 34, 38, 64, 68). Since it has been noted that amino acid residues important for the substrate specificity of Rieske nonheme iron oxygenases of the toluene/biphenyl subfamily are concentrated in the C-terminal portion of the
subunits (26, 51, 69), primer sets enabling amplification of this region are of special interest. SSCP analysis has been shown to be a sensitive method for the discrimination of gene homologues that is able to distinguish sequence differences as small as one base substitution in several hundred bases (33, 37, 54). However, it has also been demonstrated that the sensitivity of PCR-SSCP in the detection of single base substitutions decreased with segments larger than 300 bases (33, 55). Thus, the primer binding sites previously selected by Yeates et al. (68) were chosen, as these produced the shortest possible gene segment encoding the major portion of amino acids crucial for substrate specificity. Adjusting the reported primer sequences to match the toluene/biphenyl dioxygenase subfamily members together with careful optimization of the PCR-SSCP conditions allowed detection and differentiation of all reference ISP
gene segments down to single base differences.
PCR-SSCP analysis of ISP
gene segments obtained from the three differently contaminated soil samples showed that there was a high degree of catabolic genotype diversity. Sequence analysis indicated that the detectable polymorphisms were dominated by ISP
sequence types related to the ethylbenzene/isopropylbenzene oxygenase branch rather than by sequences related to the expected toluene/benzene oxygenase branch (Fig. 3). In particular, ISP
types A and C not only were observed in all three soil samples, suggesting that the environmental conditions were selective for strains possessing these ISP
types, but also were predominant in benzene- and/or toluene-degrading bacterial isolates previously isolated from these soil samples (Fig. 2A and B).
This finding was unexpected as, to our knowledge, Rieske nonheme iron oxygenases closely related to the ethylbenzene/isopropylbenzene oxygenase branch were previously not observed to be prevalent in BTEX-contaminated environments. Previous investigations, however, have concentrated on the presence and abundance of ISP
genes more closely related to the toluene/benzene and biphenyl dioxygenase branches by using DNA:DNA hybridization (59) or PCR-based assays (15, 44). The use of specific DNA probes and oligonucleotide primers may have failed to detect genes exhibiting lower levels of homology with the targeted genes. More-detailed studies on the diversity of Rieske nonheme iron oxygenases in pristine soil and soil with a history of aromatic hydrocarbon pollution (68) as well as in soil and groundwater at an aromatic hydrocarbon-contaminated site (63, 64) have been performed by analyzing clone libraries of PCR products which were amplified with degenerate primers targeting the same primer binding sites as used here (68) or primer binding sites flanking the ISP
gene region which encodes the Rieske-type [2Fe-2S] cluster (63, 64). In these studies, amplification of sequence types showing no clear affiliation to previously described ISP
proteins (68) or of a significant portion of gene segments which did not code for Rieske nonheme iron oxygenase ISP
proteins (63) was reported. In contrast to these findings, the vast majority of the peptide sequences retrieved from the BTEX-contaminated soil samples analyzed here showed a clear membership to the toluene/biphenyl oxygenase subfamily, with only a small fraction of the clones (less than 10%) prepared from faint SSCP bands exhibiting no clear pattern of relationship to previously described Rieske nonheme iron oxygenase ISP
proteins. A comparison of the ISP
peptide sequences obtained in this study with those reported by Taylor et al. (63, 64) suggests that the structures of dioxygenases putatively involved in pollutant degradation are clearly different in the organic waste landfill soil and groundwater samples and the industrial site analyzed here. Common to both environments seems to be the low abundance of genes closely related to the toluene/benzene dioxygenase branch of the toluene/biphenyl subfamily, suggesting that the previously characterized dioxygenases may not play significant roles in the sites under investigation.
As indicated above, the ISP
sequence types A and C were not only abundant in all three contaminated soil samples analyzed in this study but also predominant among the benzene-degrading isolates recovered from these soil samples. Overall, ISP
gene segments of the toluene/biphenyl subfamily were amplifiable from 33 out of 34 isolates that were capable of growth on benzene, and biochemical analysis supported the notion that, at least in the five representative strains analyzed, benzene metabolism is initiated by dioxygenation, yielding benzene dihydrodiol as an intermediate, rather than monooxygenation, which would result in phenol as intermediate (Table 3). In a recent study, the presence of alternative modes of benzene/toluene activation has been suggested for a subset of the isolates analyzed here, whereas genes encoding Rieske nonheme iron oxygenases were not detected (34). Evidently, the primers designed by Hendrickx et al. (34) for detection of toluene/biphenyl ISP
genes failed to amplify the genes reported here. Moreover, there was no correlation between the growth phenotype and the catabolic genotypes detected by Hendrickx et al. (34) (Table 3). As benzene lacks a side chain, this substrate cannot be transformed by xylene monooxygenases (genotype xylM). Furthermore, neither the strains harboring a gene coding for a phenol hydroxylase
subunit (genotype tbmD) nor strain 1YB2, which contained a gene fragment probably encoding an
subunit of aromatic monooxygenases (genotype tmoA), showed significant oxygen uptake with phenol. Therefore, degradation of benzene via phenol as the intermediate seems rather unlikely.
The fact that all isolates harboring genes similar to those encoding isopropylbenzene dioxygenases were capable of mineralizing benzene but only a subset was able to grow on toluene was on the first view astonishing, as all previously described isolates harboring ethylbenzene/isopropylbenzene dioxygenase genes (P. putida RE204, P. putida 01G3, Pseudomonas sp. strain JR1, Pseudomonas fluorescens IP01, P. fluorescens CA-4, and Thauera sp. strain DNT-1) were reported to mineralize toluene (1, 16-19, 22, 57). However, there was a strong correlation between the growth phenotype and the ISP
type detected (Table 2), with only isolates containing ISP
type C, B, or D unable to grow on toluene, suggesting that the differences in substrate specificity are related to structural differences in the
subunits of the terminal oxygenase components. Alignments of ISP
peptides retrieved from BTEX-contaminated soil and isolates to those of related proteins (see Fig. S2 in the supplemental material) revealed that the residues His234 and His240 (numbered according to P. fluorescens IP01 cumene dioxygenase [IP01-CumDO]), which are involved in the coordination of the mononuclear iron in the active site of IP01-CumDO and RHA1-BphDO (20, 27), as well as Asp231, which was postulated to play the role of an electron transfer bridge (20, 27), were completely conserved among the ISP
peptides retrieved from bacterial isolates and environmental samples. Among the 14 residues proposed to define the inner surface of the substrate-binding pocket in IP01-CumDO (20), 8 were found to be strictly conserved among the ISP
sequences retrieved in this study (Gln227, Phe228, Met232, Ala235, His323, Leu333, and Phe378, in addition to the above-mentioned His240). The isopropylbenzene dioxygenase-like ISP
sequences derived from the soil samples and isolates in this study differed from the IP01-CumDO sequence at up to three residues defining the inner surface of the substrate-binding pocket (see Fig. S2 in the supplemental material). Strikingly, isolates harboring an A321L ISP
variant (ISP
types A and E) were capable of mineralizing both toluene and benzene, whereas those harboring I288M and A321M variants (ISP
type C and C**) were capable of mineralizing benzene only. Modeling performed on the basis of the IP01-CumDO crystal structure (pdb-id 1WQL [20]), using "O" (35) for manual replacement of side chains and choice of the best-fitting rotamer, revealed that replacing Ile288 and Ala321 in IP01-CumDO with bulky methionine significantly narrowed the available space in the substrate-binding pocket (data not shown). This gives a possible explanation for the different phenotypes and the negligible oxygen uptake of isolates 1YB2 and 3YdBTEX1b with toluene as the substrate. Also, ISP
types B and D as well as K1/K2 and P contained bulky residues (methionine and phenylalanine, respectively) in positions 288 and 321, and the two single isolates containing ISP
type B and D were also capable of mineralizing benzene only. Taken together, these results indicate a significant contribution of the amino acid side chains in positions 288 and 321 to the shape and size of the substrate-binding pocket. Bagnéris et al. (3) also had suggested that less bulky amino acids in position 321 should facilitate access to toluene, which, as a substrate, is larger than benzene.
In the well-characterized isopropylbenzene-degradative pathways of P. putida RE204 and P. fluorescens IP01 (1, 22), the genes coding for the oxygenase subunits are organized in operons together with genes encoding a cis-dihydrodiol dehydrogenase and an extradiol dioxygenase (C23O) of the subfamily I.3.A (23). In the isolates 1YB2 and 1YdBTEX2, which had been analyzed in more detail and obviously do not express any 3-isopropylcatechol 2,3-dioxygenase (subfamily I.3.A) during growth on benzene (37), benzene mineralization seems to involve an isopropylbenzene dioxygenase-like Rieske nonheme oxygenase for initial activation and a subfamily I.2.A C23O (36, 37) for catechol cleavage. As shown in Table 2, in addition to strain 1YB2, 12 further independently selected isolates (affiliated to three distinct Pseudomonas lineages) harbored the same combination of a type C (or C**) ISP
gene and a type A C23O gene. While the isolates harboring ISP
type C contained the C23O type A exclusively, the strains harboring ISP
type A were found to contain different C23O subfamily I.2.A genes or failed to yield a PCR product with primers targeting C23O subfamily I.2.A genes; the latter result indicated that these strains obviously contained a C23O of other subfamilies. The high degree of sequence homology between ISP
type A and C on the one hand and C23O type A2 and A1 on the other suggests a close evolutionary relationship between the two gene variants. The flexible combination of ISP
type A genes with different types of C23O genes indicates that the pathways were assembled in a modular fashion, in which the C23O genes were retrieved from different sources, possibly driven by the selective conditions of the environment, and the occurrence of identical or nearly identical genes and pathway module combinations in different hosts indicates that both the ethylbenzene/isopropylbenzene dioxygenase and the catechol 2,3-dioxygenase are localized on mobile genetic elements and have been spread through horizontal gene transfer.
C23Os of subfamily I.2.A have previously not been reported to be involved in the degradation of benzene or toluene via a dioxygenolytic route. They have, however, been commonly described as involved in phenol (56), toluate (12, 67), and salicylate (28) degradation. The respective catabolic operons comprising the C23O I.2.A gene have a rather stable structure and comprise genes coding for both the hydrolytic branch and the oxalocrotonate branch of the meta-cleavage pathway (66). Recruitment of both branches of the meta-cleavage pathway can be assumed to be advantageous to the presence of the hydrolytic branch only (as present in strains like RE204), particularly for benzene degradation (32). Both branches were expressed in isolates 1YB2 and 1YdBTEX2 (37).
In summary, bacteria from the examined BTEX-contaminated soil samples seem to possess a genotypic flexibility, which is important for their adaptation and evolution while facing challenging and continuously changing conditions in these ecosystems. Screening methods targeting genes encoding key enzymes in pollutant degradation as molecular markers, such as PCR-SSCP analyses, can be utilized to determine both the biodegradative potential and predominant catabolic pathways in an environment. Further studies on catabolic gene diversity and genetic flexibility will certainly be important elements in our understanding of natural bioremediation processes.
This research was supported by project BIOTOOL (GOCE-003998) of the European Commission and PI 312/2-1 of the Deutsche Forschungsgemeinschaft.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
subunit. J. Bacteriol. 182:5495-5504.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»