Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2006, p. 1476-1486, Vol. 72, No. 2
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.2.1476-1486.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Geological Survey of Denmark and Greenland (GEUS), Department of Geochemistry,1 Royal Veterinary and Agricultural University (KVL), Department of Natural Sciences,2 Royal Veterinary and Agricultural University (KVL), Department of Ecology, Copenhagen, Denmark3
Received 22 August 2005/ Accepted 22 November 2005
|
|
|---|
|
|
|---|
The degradation of phenoxy herbicides such as 2,4-dichlorophenoxyacetic acid (2,4-D), MCPA, and related compounds has been studied intensively during the last 40 years, with the majority of the work being done on 2,4-D. This, among other studies, has revealed detailed knowledge about the degradation and mineralization kinetics (27, 38), the catabolic pathways (6, 21, 46), the genes and enzymatic systems involved in the degradation (6, 20, 24, 35), and the phylogenetic composition of isolated microbial degraders (25, 53).
Numerous kinetic models have been developed to describe the mineralization of xenobiotic compounds in the environment (for a review, see reference 17). Typically, the models are divided into different forms of degradation kinetics that either take the growth of microorganisms into account or not. Furthermore, different models that take different forms of growth into account, depending on initial cell and substrate concentrations, have been developed (7, 52).
The soil bacterium Ralstonia eutropha JMP134, originally isolated from an Australian agricultural soil sample (44), contains genes for the complete pathway of 2,4-D and MCPA degradation. These genes are the ones most extensively studied and have become a model for the study of phenoxy acid degradation (11, 55). This strain harbors plasmid pJP4, which contains all of the structural and regulatory genes needed to convert these compounds to 2-chloromaleylacetic acid, which is further metabolized by chromosomally encoded gene products to finally yield CO2 and chloride (14, 15, 22, 31, 32, 34). The first step in the degradation pathway is initiated by an
-ketoglutarate-dependent dioxygenase, encoded by the tfdA gene, which cleaves the acetate side chain to produce the corresponding phenol 4-chloro-2-methylphenol (MCP) (Fig. 1) (18, 19, 54). Further degradation steps are initiated and regulated by gene products of the genes tfdB (45), tfdCDEF (35, 45), tfdR, and tfdS (37, 39).
![]() View larger version (12K): [in a new window] |
FIG. 1. Catabolic activity of the TfdA enzyme in the degradation of MCPA. -KG, -ketoglutarate.
|
-proteobacteria containing tfdA genes related to that of R. eutropha JMP134. The second group consists of
-proteobacteria that are closely related to Bradyrhizobium spp. and contains the tfdA-like gene tfdA
(60% identical to the canonical tfdA gene). The strains in this group were first isolated from pristine environments in Hawaii, Canada, and Chile (30) and subsequently from arable soil in Japan (24). Members of the third group are
-proteobacteria belonging to a species of the genus Sphingomonas. They are also characterized by using the tfdA
gene product during the initial step of the pathway instead of the TfdA enzyme (20). As well as this, several strains in the second and the third groups are known to contain the cadA gene, which is related to the tftA gene involved in the degradation of 2,4,5-trichlorophenoxyacetic acid (25). The main objective of this paper was to describe the functional diversity of MCPA degraders in a profile of MCPA-treated Danish soil. Furthermore, PCR using primers targeting the three groups of phenoxy acid degradation genes was performed to investigate whether the genes were present and which primer set revealed the most specific amplification product. In light of this, one degenerated primer set (this work) targeting the canonical tfdA gene for real-time quantitative PCR was chosen. The relative increase in the tfdA gene concentration was monitored during a degradation experiment in microcosms. Denaturing gradient gel electrophoresis (DGGE) analyses of the tfdA amplification product were used to elucidate whether only one or several homologues of the tfdA gene were involved in the MCPA degradation. This revealed a number of clearly separated bands, which were sequenced and compared to known tfdA sequences. We further examined the mineralization kinetics of MCPA depending on the effects of the initial substrate and cell concentrations in top- and subsoil. The topsoil was suggested to contain a higher initial cell concentration than the subsoil. The initial cell concentration and the growth of MCPA degraders were determined and monitored during a degradation period by quantitative real-time PCR of the tfdA gene. Microbial growth during degradation of MCPA in natural soil has previously been suggested but has only been quantified using less accurate techniques (12, 41).
|
|
|---|
Soils.
An organic topsoil sample (10 to 25 cm) and a subsoil sample (35 to 50 cm) from a sandy profile were used for the experiments. The profile was taken from the Danish agricultural experimental station in Jyndevad, located on the outwash plain in southern Jutland. The soil has been used for organic farming since 1995, and additionally, it has not been treated with any form of phenoxy acid herbicide since at least 1988. The soils were collected and mixed manually in October 2003 and stored in the dark for 2 months at 5°C before they were passed through a 2-mm sieve. The soils were thereafter stored in the dark at 18°C until 10 days before use. Defrosting and acclimatization were performed in the dark for 10 days at 10°C (41) (for further soil characteristics, see Table 1).
|
View this table: [in a new window] |
TABLE 1. Soil characteristics
|
Subsequently, 20 g (dry weight) of the topsoil and 15 g (dry weight) of the subsoil were transferred to 250-ml Pyrex red cap bottles (airtight). The soils were spiked with 2.3 ml (for the topsoils) and 1.5 ml (for the subsoils) of the respective stock solutions. Stock solutions with MCPA were carefully mixed into the soils with a sterile glass pipette, and the microcosms were incubated at 10°C in the dark. Triplicate bottles were set up for each soil scenario for eight sampling points, and in order to prevent the development of anaerobic conditions in the microcosms, they were aerated by leaving the flasks open for 20 min in a sterile hood every 14 days.
Microcosm setup for mineralization experiment.
A mineralization experiment was set up with 14C-labeled MCPA at the same concentrations described above. Stock solutions were prepared as described above for the degradation experiment. [U-14C]MCPA (specific activity of 159.7 µCi/mg, radiochemical purity of >95%; Izotop, Budapest, Hungary) was then added to the stock solutions in trace amounts to achieve an initial radioactivity in each microcosm of 60,000 dpm. The mineralization experiment was performed essentially as described previously by Mortensen and Jacobsen (41), except that 10 g of soil was used and no soil was removed during the experiment. We have previously determined the mass balance for the experimental setup, and recoveries of 93% ± 1.5% were found (5, 9).
MCPA extraction and quantification.
On days 0, 6, 12, 22, 33, 50, 68, and 115 after incubation, triplicates of each soil-MCPA degradation series (Top2.3, Top20, and Sub2.3) were extracted by pressurized liquid extraction for the determination of MCPA amounts. The extractions were performed with an ASE 200 system (Dionex, Sunnyvale, CA). The soils from the degradation assay were packed in 33-ml extraction cells sufficient to contain 40 to 50 g soil. To assist the extractions, Ottawa sand (Fisher scientific, Loughborough, United Kingdom), sufficient to fill out the remaining volume (approximately 1:1 [wt/wt]), was mixed into the soil samples. The entire extraction process was carried out at 50°C. The solvent used for the extractions was a 50:50 (vol/vol) mixture of methanol-water (high-performance liquid chromatography-grade methanol; Fisher Scientific, Loughborough, United Kingdom). Other extraction conditions were as follows: a 40% solvent volume was used for a 10-min static time in one extraction cycle and a 60-s purge time with nitrogen at 1,500 lb/in2 pressure (1 lb/in2 = 6,894.76 Pa). The extraction method was optimized for soil samples spiked with MCPA, resulting in extraction efficiencies of 80% ± 2% for the topsoil and 92% ± 4% for the subsoil.
A soil extract of 15 to 20 ml was obtained in 50-ml glass vials. After being thoroughly mixed, an aliquot was filtered through a 0.45-µm polytetrafluoroethylene filter (Titan2; Sun-Sri) to remove particles. To avoid overloading of the liquid chromatography (LC) mass spectrometry system, the filtrates were diluted with a methanol-water (50:50) solution to obtain MCPA concentrations in the range of 50 to 150 µg/liter.
The quantification of MCPA in the soil extracts was performed using an LC tandem mass spectrometry (MS/MS) system. Twenty microliters of the samples was injected, and the compounds were separated on a 100- by 2.1-mm Xterra RP18 column (3.5-µm particle size) from Waters (Milford, MA) using a Waters Alliance model 2695 high-performance liquid chromatography system. A mobile phase composed of methanol-water-0.2% acetic acid (78:20:2) at 25°C with a flow rate of 0.18 ml/min was used. For detection and quantification of the MCPA compound, a Quattro Ultima triple quadrupole mass spectrometer from Micromass (Manchester, United Kingdom) equipped with an electrospray ionization unit was used. Electrospray ionization in the negative mode was used, while further conditions were as follows: capillary voltage was 3.20 kV, cone voltage was 80 V, and collision energy for MS/MS was 20 eV. Furthermore, nitrogen was used as the desolvation gas at a flow rate of 600 liters/h. A limit of detection of 1 µg MCPA liter1 was determined for the LC-MS/MS system.
DNA extraction.
Soils for DNA extractions were set up parallel to the degradation experiment. Red cap bottles with 20 g (dry weight) of soil amended with 0, 2.3, and 20 mg MCPA kg1 (dry weight) soil for the topsoil (Top0, Top2.3, and Top20) and 0 and 2.3 mg MCPA kg1 (dry weight) soil for the subsoil (Sub0 and Sub2.3) were set up in triplicate. The amount of soil sampled for DNA extractions was 0.5 g (wet weight) to
0.4 g (dry weight), which was accomplished using a sterile stainless steel sampler designed to collect a composite sample of 0.5 g soil from 5 to 10 distinct spots in the bottle. The same bottle was used throughout the experiment and was aerated parallel to the bottles in the degradation experiment as well as for sample collection on the same days as in the degradation experiment. Whole-community DNA was extracted from the 0.5-g wet soil samples with the FastDNA SPIN kit for soil (Bio101, Inc., Carlsbad, CA). The protocol recommended by the manufacturer was followed, except that the bead-beating step was modified and a freeze-thaw step was included. The beat-beating step was modified to four 30-s pulses at speed 4 instead of one 30-s pulse at speed 5.5 in the FastPrep FP 120 instrument (Bio101), and the freeze-thaw steps were performed for 1 h at 80°C and subsequently for 30 min at 37°C. DNA was finally eluted in 100 µl DNase/RNase-free water (Bio101, Carlsbad, CA) and stored at 80°C.
Primer testing.
All PCRs were carried out as real-time PCR assays in an iCycler iQ (Bio-Rad, Hercules, CA). In order to select the most effective primer set, five primer sets were tested (Table 2) using two different commercial PCR master mix kits. tfdA* and tfdA** were the most intensively tested primers. tfdA* yielded a 307-bp fragment, while the tfdA** yielded a 210-bp fragment. tfdA** was designed in connection with this work by BLAST alignment analysis (3) of 22 known tfdA gene sequences found with R. eutropha JMP134 as the search sequence. The GenBank accession numbers for the 22 tfdA genes were as follows: M16730, AY238497, AY238496, U43197, U43276, AF439768, AY078159, AF029344, U32188, U25717, U87394, U22499, AF181982, AF182758, AF176240, U65531, U43196, AY238495, AY238494, AY238493, AY238492, and AB074491. Primers were designed to attach to most conserved regions, and since the primer set was designed on the basis of tfdA gene sequences belonging to class I and class III, it was expected that it was highly specific for these genes. Primers were obtained from MWG Biotech (Ebersberg, Germany). The first of the two commercial SYBR green PCR master mix kits we tested was iQ SYBR green supermix (Bio-Rad, Hercules, CA) containing PCR buffer, 0.4 mM of each deoxynucleoside triphosphate, 50 U/ml of iTaq DNA polymerase, and 6 mM MgCl2. This kit was used mainly with the tfdA* primer set and through the testing of tfdA
*, tfdA
**, and cadA primer sets. The other kit we tested was the QuantiTect SYBR green PCR kit (QIAGEN, Crawley, United Kingdom) containing deoxynucleoside triphosphate mix, HotStar Taq DNA polymerase, PCR buffer, Rox, and 2.5 mM MgCl2. This kit was only used for PCR with the tfdA* and tfdA** primer sets. The 25-µl reaction mixtures contained 0.4 µM of each primers, 12.5 µl of the respective SYBR green mix, 25.5 µg bovine serum albumin (Amersham Bioscience, Buckinghamshire, United Kingdom), 2.5 µl of 1:10-diluted DNA extract, and RNase/DNase-free water to complete the 25-µl volume. The DNA extracts were diluted 1:10 in RNase/DNase-free water to reduce the influence of PCR-disturbing factors from the soil, such as humic acids.
|
View this table: [in a new window] |
TABLE 2. PCR primers
|
Further investigations of the PCR products were done by gel electrophoresis of 8-µl PCR products on a 1.5% agarose gel in 1x Tris-acetate-EDTA buffer. The gels were stained in ethidium bromide and visualized under UV light.
In case the gels revealed more than one band, a Southern blot hybridization was performed as described previously by Jacobsen (26) and Sambrook and Russell (49), except that the digestion reaction was excluded. The 32P tfdA probe for the hybridization was generated using [32P]dCTP (Amersham, United Kingdom) and the Random Primed DNA Labeling kit (Roche Applied Science, Penzberg, Germany) based on the method described previously by Feinberg and Vogelstein (16). A PCR product derived from DNA of the tfdA-containing strain B. cepacia DB01(pRO101) with the tfdA-intern primer set (Table 2) was used as a template.
Quantitative real-time PCR.
After optimization of primers, quantitative real-time PCR was performed only with the tfdA** primer set as described above on the triplicate samples from days 0, 6, 12, 22, 33, 50, 68, and 115. All quantitative real-time PCRs were set up in triplicate and included three negative control samples in each PCR setup.
Standards for the quantitative PCR were prepared using the well-characterized phenoxy acid degrader R. eutropha JMP134. After inoculation into 0.5 g of the topsoil in amounts of 8 x 106, 8 x 105, 8 x 104, 8 x 103, 8 x 102, and 8 x 101 cells g1 soil, the DNA was extracted from the soils as described above, and the extracts were diluted 10-fold to reduce the effect of humic acid disturbances. R. eutropha JMP134 was used as the positive control strain as well.
DGGE analysis.
PCR for the DGGE analysis was performed as described above with the QuantiTect mastermix kit. The tfdA** primer set was used with a GC clamp attached to the 5' end of the forward primer (Table 2). DGGE analyses were performed on real triplicates of DNA extracts from days 0, 6, 12, 33, and 68 in order to secure enough wells for one data series on one gel. The technique described previously by Muyzer et al. (42) was used for separating and visualizing the tfdA PCR products of the DNA extracts. The specific details of the method used are as follows. Gels contained 8% acrylamide and a urea-formamide gradient of 55 to 70% (where 100% denaturant contained 7 M urea and 40% formamide). Electrophoresis was carried out at 70 V for 17 h in 1x Tris-acetate-EDTA buffer using a D-code apparatus (Bio-Rad, Hercules, CA). After electrophoresis, the gels were stained with SYBR gold (Molecular Probes, Eugene, OR) for 45 min, followed by a brief rinse with MilliQ water, and then visualized under UV light and finally photographed. Bands that appeared to have a unique pattern for the soil scenarios were stabbed with sterile pipette tips, which were placed in 20 µl RNase/DNase-free water and rinsed repeatedly. These served as templates in further PCRs using the tfdA** primer set (Table 2) without the GC clamp to obtain enough DNA material for sequencing. The nucleotide sequencings were performed by MWG Biotech (Ebersberg, Germany).
Data analysis.
The amounts of extracted MCPA were corrected for the extraction efficiency, after which linear regression analyses were performed on the steep segment of the degradation curves to fit it to zero-order kinetics with the following equation:
![]() | (1) |
For the results from the mineralization experiment, the calculated cumulative mineralization values were corrected for background radioactivity. After this, nonlinear regression analyses were performed as described previously by Mortensen and Jacobsen (41).
Nucleotide sequence accession numbers.
The GenBank accession numbers for the DNA sequences reported in this study are DQ272405, DQ272406, DQ272407, DQ272408, DQ272409, DQ272410, DQ272411, DQ272412, DQ272413, DQ272414, DQ272415, and DQ272416 for A1, A2, A3, A4, A5, A6, A7, A8, A9, B1, B2, and B3, respectively.
|
|
|---|
![]() View larger version (15K): [in a new window] |
FIG. 2. Degradation of MCPA in the three soil scenarios measured by LC-MS/MS. Linear regression lines of the steep segments are shown. Solid line, Top2.3; dashed line, Top20; dotted line, Sub2.3. Error bars indicate standard errors.
|
|
View this table: [in a new window] |
TABLE 3. Parameter estimates from the zero-order kinetics modela of the steep part of the degradation curves and calculated DT50 values
|
![]() View larger version (19K): [in a new window] |
FIG. 3. Mineralization of 14C-MCPA in the three soil scenarios. (A) Curves show triplicate cumulative data from Top2.3, Top20, and Sub2.3. The dashed line shows a nonlinear fit to the exponential form of the 3/2-order model. Due to the large standard deviation, the plot of the triplicates for Sub2.3 was not fitted. (B) Single plots of the Sub2.3 replicates. The dotted line shows nonlinear fits to the linear form of 3/2-order model, and the dashed line shows a nonlinear fit to the exponential form of the 3/2-order model (plots 2 and 3). Plot 1 was not satisfactorily fitted to any model. Error bars in A indicate standard errors.
|
Primer testing.
To make the real-time PCR quantification technique as accurate as possible, the choice of primer set was critical. Therefore, a well-characterized primer set covering a broad range of known tfdA sequences was chosen (59). Unfortunately, we found that this primer set was likely to generate diverse bonds, indicating that the PCR product contained genes of more than one length (Fig. 4A). One band (band 2) had the same length, approximately 350 bp, as the positive control strain, R. eutropha JMP134 (band 3).
![]() View larger version (129K): [in a new window] |
FIG. 4. Primer optimization. (A) Standard agarose gel electrophoresis of real-time PCR amplification products using the tfdA* primer; the iQ SYBR green supermix and soil DNA from days 0, 6, 12, 33, 50, 68, and 115 of Top20 were used as a template. The standard is R. eutropha JMP134 inoculated into the soil in concentrations of 8 x 106, 8 x 105, and 8 x 104 cells/g soil, with 8 x 106 cells/g soil shown to the left. Bands 1 and 2 represent the most characteristic bands, as described in the text. Band 3 represents the strain known to contain the tfdA gene R. eutropha JMP134. The ladder is a 100-bp DNA ladder. (B) Southern hybridization analysis of the agarose gel shown in A using a tfdA probe. Similarity can be detected for bands 2 and 3.
|
Consequently, we used the tfdA** primer set for further quantitative real-time PCR assays. Primer testing of primer sets tfdA
* (51), tfdA
** (24), and cadA (25) were performed as well (Table 2). Only tfdA
** amplified a PCR product that may be possible to use for quantitative real-time PCR (data not shown).
Quantitative real-time PCR analysis.
The quantitative real-time PCR analysis indicated that growth occurred in all three soil scenarios during the degradation period (Fig. 5). In Top2.3, the tfdA concentration increased to a maximum of 3.0 x 104 genes g1 soil at day 12 (Fig. 5A), while the maximum gene concentrations in Top20 and Sub2.3 were 7.0 x 105 and 2.6 x 104 genes g1 soil, obtained at days 22 to 33 and 50 to 68, respectively (Fig. 5B and C). The number of genes is based on the assumption that R. eutropha JMP134(pJP4) grown on rich media and used for real-time PCR standard curve preparation only contains one copy of the tfdA gene. The increase in the number of tfdA genes is inversely related to the concentration of MCPA in the three soils (Fig. 5A to C). Lag phases were clearly expressed and were shortest for Top2.3 and Top20 and longest for Sub2.3. It is interesting, however, that the increase in the number of tfdA genes occurred before the exponential mineralization phases did. It is also interesting that the maximum concentrations of the tfdA genes were reached while half of the MCPA remained and that a decrease in concentrations of tfdA genes was subsequently observed in the topsoils.
![]() View larger version (20K): [in a new window] |
FIG. 5. Nonfitted mineralization and degradation curves compared with the log number of tfdA genes using the tfdA** primer set in the soil. (A) Top2.3; (B) Top20; (C) Sub2.3. Error bars indicate standard errors.
|
![]() View larger version (19K): [in a new window] |
FIG. 6. Melting profiles of real-time PCR amplification products of soil DNA from days 0, 12, and 22 in Top20 using the tfdA-II primer set and the QuantiTect SYBR green PCR kit. Furthermore, the melting profile of the PCR amplification product from the tfdA-containing strain R. eutropha JMP134 is shown. The profile displays the negative first derivative of temperature versus relative fluorescence units (RFU) [d(RFU)/dT] plotted against temperature.
|
![]() View larger version (88K): [in a new window] |
FIG. 7. Functional diversity of MCPA degrader genes using PCR/DGGE analysis with the tfdA** primer set. The gene product was separated using 55 to 70% urea gels. DNA extracts from days 0, 6, 12, 33, and 68 were analyzed in real triplicates. For the marker profile (M), DNA extracts from soil samples inoculated with the tfdA-containing strain R. eutropha JMP134 were used. A1 to A10 indicate bands selected for sequence analysis from topsoil DNA, and B1 to B3 indicate bands selected for sequence analysis from subsoil DNA. A shows the DGGE profile of DNA extracted from Top2.3, B shows the DGGE profile of DNA extracted from Top20, and C shows the DGGE profile of DNA extracted from Sub2.3. D shows the DGGE profile of DNA extracted from Top0, and E shows the DGGE profile of DNA extracted from Sub0.
|
![]() View larger version (28K): [in a new window] |
FIG. 8. Sequence alignment of the DNA obtained with the tfdA** primer set. The position number relates to GenBank accession number M16730 (R. eutropha JMP134). The class I tfdA gene sequence represented by R. eutropha JMP134 and the class III tfdA sequence represented by Burkholderia cepacia pIJB were both obtained from GenBank. A2 and A9 correspond to the sequences obtained at day 68 and day 0, respectively. These products were run on a DGGE gel and stabbed prior to sequencing as indicated in the legend of Fig. 7. A2 is highly homologous to bands A1, A3 to A7, A10, and B1 to B3, while A9 and A8 are homologous. Black shading shows highly conserved regions. Designations in brackets are GenBank database accession numbers.
|
The sequences obtained from the bands stabbed from the DGGE gels were compared to existing sequences encoding
-ketoglutarate-dioxygenase obtained from GenBank, and a phylogenetic tree was constructed (data not shown). The phylogenetic analysis showed that bands A1 to A7, A10, and B1 to B3 (Fig. 7) were located in the same cluster as the class III tfdA genes showing 96 to 99% homology to the indigenous genes in this class (Fig. 8). Bands A8 and A9, on the other hand, were located in the same cluster as the class I tfdA genes showing 98 and 99% homology, respectively, to the R. eutropha JMP134 tfdA gene (Fig. 8). Furthermore, the GC content of bands A1 to A7, A10, and B1 to B3 was approximately 63%, while it was 69% for bands A8 and A9.
|
|
|---|
-I primer set yielded a PCR product for the same soil. This indicates that the latter group of phenoxy acid degraders may be involved in the degradation of 2,4-D and probably also in the degradation of MCPA, especially in soils that are not frequently exposed to the herbicide. However, during the test of primers targeting the canonical tfdA gene, the tfdA
gene, and the cadA gene, we observed that PCR products were yielded with the canonical tfdA primer and with the tfdA
primer but not with the cadA primer set. This indicates that organisms that use the tfdA
gene product in MCPA metabolism may be present as well.
Growth among bacterial cells and increases in catabolic genes during degradation of phenoxy acids have been reported elsewhere previously (see, e.g., references 12, 28, and 36), but only Lee et al. (36) used a highly accurate quantitative real-time PCR method. They reported concentrations of 2.0 x 107 tfdA genes g1 in activated sludge sequentially amended with 300 mg 2,4-D kg1. Concerning the higher amount of carbon substrate, this corresponds to what was found in the present study. de Lipthay et al. (12), however, reported slightly lower increasing concentrations of the tfdA gene in water samples amended with approximately 20 mg 2,4-D liter1, as they observed maximum concentrations of
103 genes ml1.
Analysis of the functional diversity of the indigenous MCPA degrader community by DGGE on the tfdA genes extracted has to our knowledge not previously been published. In this study, results from two different methods both indicated a change in the community of the bacteria containing the tfdA gene during a degradation period of 68 days. First, we saw a change in the maximum of the melting profiles of the PCR products from the quantitative PCR assay of the topsoil DNA from 90°C to 87.5°C. This indicates that the tfdA PCR products from the organisms grown on MCPA had a lower GC content and therefore were different from the genes present in the soil not exposed to MCPA. This corresponds with the results from the DGGE analyses in which the PCR products from the DNA sampled during the first days of the experiment had a tendency to denaturize at a higher denaturation grade. Interestingly, the maximum of the melting profiles from the first sampling days of the experiment were similar to that of the standard strain R. eutrophus JMP134. This was seen on the DGGE gels as well, where a major part of the DNA from the first sampling days seemed to denaturize in the lower part of the gels where DNA from the R. eutropha JMP134 also denaturized. The presence of tfdA genes in the control topsoil was interesting as well. The soil used for the experiments had not been exposed to phenoxy acid herbicides for at least 15 years; however, a background population of tfdA genes is maintained even in the absence of appropriate pesticide substrates. These tfdA genes did not increase in numbers during the degradation period, which indicates that they were placed in populations that did not grow actively.
Phylogenetic analyses of nucleotide sequences stabbed from bands on the gels revealed a large degree of homology to what was seen on the melting profiles and on the DGGE gels. All the sequences stabbed from the DGGE gels belonged to the first group of phenoxy acid-degrading organisms. The sequences stabbed from the bands located at the same denaturation grade as R. eutropha strain JMP134 were closely related to this strain and other class I strains as well. In contrast, the sequences stabbed from the bands that increased during the degradation period were related to the class III tfdA genes. In general, the role of those organisms containing class I tfdA genes in this particular soil remains unclear. A likely and quite interesting explanation may be that they are organisms that can only perform one or a few steps of the degradation pathway, while the class III tfdA gene-containing strains are capable of degrading the compound completely. If the class I gene-containing organisms are only capable of degrading the first step, for instance, it may not be likely to grow as a result of the degradation process. If the class III gene-containing bacteria, on the other hand, are capable of degrading MCPA completely, they may be able to grow as a result of the degradation process. This explanation is further supported by the fact that we actually were able to detect the degradation product MCP at days 0 to 22 of the degradation experiment in Top2.3 and at days 22 and 33 in Top20. Neither this compound nor the class I tfdA gene was found in Sub2.3. Until now, the majority of the research on catabolic genes involved in phenoxyacetic acid degradation has been done with respect to 2,4-D, and in general, it has been accepted that organisms containing the class I tfdA gene were capable of degrading this particular compound completely. However, it may be suggested that minor differences in the metabolism between MCPA and 2,4-D have made the class I tfdA gene-containing bacteria in this particular soil incapable of degrading MCPA completely. Classes I to III were previously defined by Fulthorpe et al. (20) and Itoh et al. (24) based on homology and whether the tfdA gene is carried on a plasmid or not. The class III tfdA genes are 77% identical to the class I genes, whereas class II only consists of Burkholderia sp. strain RASC and a few strains that are 76% identical to the class I strains and 80% identical to the class III strains (24). As is the case for the strains containing the class I gene, the organisms containing the class III gene belong to a broad range of ß- and
-proteobacteria. Even though the diversity of the tfdA gene is not yet fully understood, the DGGE analysis indicated that up to five different tfdA genes were involved in the degradation of MCPA. However, no clear phylogenetic differences could be determined between those strains represented by the strong bands on the DGGE gels. This is probably due to comigration and other types of bias likely to be observed during DGGE analysis (43). Vallaeys et al. (60) reported that strains isolated from the same soil microcosm containing highly homologous tfdA genes were different based on their 16S rRNA genes. This may improve the hypothesis that several different species may contain the catabolic gene and be involved in the degradation and that the five bands sequenced from the DGGE gels represent tfdA genes from different strains of MCPA-degrading organisms.
The kinetics of mineralization and degradation observed in this study seem to be slightly slower, especially in the subsoil experiment, than what was reported previously by Caux et al. (8) and Mortensen and Jacobsen (41) but equal to what was reported previously by Helweg (23). This rather slow degradation/mineralization may be partly due to the high MCPA concentrations used, which, as seen in this study, results in a longer lag phase. As adaptation has been shown to increase the rate of mineralization (12), a more likely explanation may be the lack of exposure to MCPA during the last 15 years. Similar results were seen for the soil not previously exposed to phenoxy acid herbicides used by Helweg (23).
For the topsoils, the mineralization curve is shifted towards longer reaction times compared with the degradation curves, indicating metabolite formation during degradation. Even though the delay of mineralization compared to degradation was quite small, we detected the first degradation product, MCP, in Top2.3 and Top20 but not in Sub2.3. However, no proper accumulation of the metabolite was observed, as the compound was not detected after days 22 and 33 in Top2.3 and Top20, respectively. It was rather surprising that we did not detect the metabolite in the subsoil, as the delay of the mineralization was more pronounced for this scenario. It is possible that other metabolites evolve in the subsoil, as we only analyzed this particular compound. This indicates that the degradation of MCP may be rate limiting in topsoils. This is further supported by the increase in the number of tfdA genes, which occurred in advance of the exponential phase of the mineralization. Detection of MCP and other unidentified metabolites during the degradation of MCPA have been reported in recent studies as well (10, 29).
In general, the amounts of MCPA mineralized in topsoils did not exceed a maximum of 60 to 70% of the initial concentrations. The remaining 30 to 40% can be explained by the incorporation of 14C from MCPA into biomass and into soil organic matter. Mineralization potentials of this magnitude have been reported elsewhere previously (41, 61). Additionally, there is a tendency towards a slightly higher total mineralization in the subsoil experiment than in the topsoil, which can be explained by a larger degree of adsorption in the topsoil than in the subsoil. This was previously shown by Jensen et al. (29) to result in a larger bioavailability in the subsoil and thereby higher mineralization potentials, corresponding to the findings previously reported by Mortensen and Jacobsen (41) and Willems et al. (61).
Using kinetic modeling, we found that the mineralization in Top2.3 and Top20 was best described by the exponential form of the 3/2-order model, in agreement with results reported previously Mortensen and Jacobsen (41), who used the exponential form for soils from three depths amended with 1 mg MCPA kg1 soil. Similarly, Reffstrup et al. (47) fitted the mineralization of the phenoxypropionic herbicide mecoprop amended at 5 to 500 mg kg1 to the exponential form of the 3/2-order model as well, while experiments amended with 0.0005 to 0.5 mg mecoprop kg1 were fitted to the linear form of the 3/2-order model. This indicates that the initial concentration of the substrate may be important for whether growth occurs or not. However, this was not observed in the present study, which may be due to substrate concentrations that were too high. The fact that Top20 fitted better to the model than Top2.3 may be a consequence of the relatively higher degree of growth in Top20 than in Top2.3 compared to the substrate concentration.
In the subsoil, the result of the present study differs from what was found previously by Fomsgaard (17) and Mortensen and Jacobsen (41). They found that growth was more likely to occur in subsoil than in topsoil because the initial population of degraders in topsoils may be able to degrade the pesticide without significant growth.
Conclusion.
In summary, we found that growth of microbial degraders occurs during the degradation of MCPA. The growth was more pronounced during the degradation of high concentrations (20 mg kg1) than during the degradation of low concentrations (2.3 mg kg1). This led to better fitting of mineralization data of the high-concentration scenario to a model taking exponential growth into account. In addition, we found that a shift in the degrader population occurred during the degradation period. Class I tfdA genes were found to be present among the naturally occurring microbial population in the topsoil, which had not been exposed to MCPA for 15 years. The concentration of these genes did not increase during the degradation period. Contrary to this, we found that degraders containing tfdA genes belonging to class III grew up and became dominant. We found that at least five different tfdA-containing genotypes grew up during the degradation period. Thus, the results of an analysis of functional genes present in soil prior to exposure with a contaminant might not necessary reflect the actual population carrying out the degradation. Last, in soils containing an initial load of tfdA class I genes, the metabolite MCP was detected during the degradation experiment, but no proper accumulation of metabolites from the MCPA degradation was observed.
We thank M. Nicolaisen for phylogenetic analysis of tfdA genes, L. F. Nielsen for initial optimization of real-time PCR, M. Andersen for the valuable technical assistance, and A. Z. Nielsen and K. L. Demant for their useful suggestions on the manuscript.
|
|
|---|
-ketoglutarate-dependent dioxygenase. J. Bacteriol. 175:2083-2086.
-Proteobacteria. Appl. Environ. Microbiol. 68:3449-3454.
and cadA, homologous with genes encoding 2,4-dichlorophenoxyacetic acid-degrading proteins. Appl. Environ. Microbiol. 70:2110-2118.
-ketoglutarate dioxygenase from Burkholderia sp. strain RASC. Appl. Environ. Microbiol. 62:2464-2469.[Abstract]
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»