Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2008, p. 2728-2739, Vol. 74, No. 9
0099-2240/08/$08.00+0 doi:10.1128/AEM.02199-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Department of Civil and Environmental Engineering, University of California, Berkeley, California 94720-1710,1 North Wind Inc., 1425 Higham St., Idaho Falls, Idaho 83402,2 CDM, 1331 17th Street, Suite 1200, Denver, Colorado 80202,3 Malcolm Pirnie, Inc., 2000 Powell Street, Suite 1180, Emeryville, California 94608,4 Earth Sciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, California 947205
Received 26 September 2007/ Accepted 22 February 2008
|
|
|---|
|
|
|---|
Complete reductive dechlorination of chlorinated ethenes to the innocuous end product ethene has been demonstrated only by organisms in the genus "Dehalococcoides" (33). To date, there are six described strains of Dehalococcoides, each of which has different dechlorination capabilities (1, 8, 9, 25, 27, 34). Of the six, the genomes of strain 195 (32), strain CBDB1 (16), and strain BAV1 (http://www.jgi.doe.gov/) have already been sequenced and annotated. From the annotations and other studies, insights have been gained into the physiology and biochemical pathways of Dehalococcoides spp., including the metabolically diverse reductive dehalogenase (RDase) genes that code for the proteins that are responsible for the respiration of chlorinated ethenes. For example, the tceA gene present in strains 195 and FL2 codes for a protein that catalyzes the sequential metabolic transformation of TCE to cis-DCE (cDCE) and VC and the cometabolic generation of ethene from VC (9, 21, 22, 32), while the gene product of vcrA, found in strains VS and GT, reduces DCE isomers and VC (27, 34) and the product of bvcA found in strain BAV1 is implicated in metabolic transformation of VC to ethene; however, the protein has not yet been fully characterized (8, 15).
In addition to the three aforementioned genes, whose functions are known, multiple nonidentical putative RDase genes have also been found in the different isolates (12, 15, 16, 32) and in a Dehalococcoides-containing enrichment culture (37). From the genome annotations, 17 intact putative RDase genes are present in strain 195 (32), while strain CBDB1, an organism that metabolizes chlorobenzene and polychlorinated dibenzo-p-dioxins, possesses 32 putative RDase genes (16). For strains BAV1 and FL2, 11 and 14 putative RDase genes, respectively, have been identified (12, 15) (http://www.jgi.doe.gov/). However, the substrate specificity, function, and regulation of these putative RDase genes are mostly unknown at this time.
Monitoring of the presence of Dehalococcoides strains at solvent-contaminated field sites has typically focused on the 16S rRNA gene, using such methods as quantitative PCR (qPCR), end-point PCR, and terminal restriction fragment length polymorphism (6, 10, 18, 23, 29), while the RDase genes are less frequently monitored (3, 31). However, analysis of the 16S rRNA gene sequences of different Dehalococcoides isolates, environmental samples, and enrichment cultures has revealed that despite differences in chlorinated-ethene metabolic capabilities, the Dehalococcoides 16S rRNA gene sequences are highly conserved (5, 9, 10, 16, 34). Consequently, metabolic function cannot be inferred from Dehalococcoides phylogeny, and the detection of 16S rRNA genes alone at field sites could lead to potential misinterpretations if the Dehalococcoides spp. present do not actually have the specific metabolic capabilities of interest. Therefore, genes that are specific to functions of interest can serve as useful biomarkers whose quantification can enable effective monitoring of different Dehalococcoides activities.
Monitoring the presence or absence of target genes may also be inadequate in some cases to predict microbial activity, since the presence of a gene does not guarantee its expression or its activity (19). Further, it has been demonstrated previously for strain 195 that growth and dechlorination activity can decouple (25), resulting in cell and gene concentrations that remain static while dechlorination continues. Consequently, quantifying expression of functionally important genes may give a more robust surrogate indicator of microbial activity. Previous findings have indicated that RDase genes are expressed in response to the presence of different chlorinated organics (7, 13, 14, 17, 30, 37), and RDase gene expression has been correlated to the physiological ability of Dehalococcoides spp. to dechlorinate beyond cDCE (17). Previous experiments have also shown that the half-life of RDase transcripts in the absence of chlorinated ethenes is between 4 and 6 h (17). The short-lived nature of the RDase mRNA suggests that it could serve as a sensitive biomarker to indicate the activity of Dehalococcoides spp.
In this study, qPCR and reverse transcription (RT-qPCR) were applied to quantify the DNA and RNA of targeted RDase genes (tceA, bvcA, and vcrA) and 16S rRNA genes of Dehalococcoides spp. at a TCE-dense nonaqueous-phase liquid (DNAPL) site during different phases of a bioremediation application. Biostimulation with whey to promote the activity of indigenous microorganisms, bioaugmentation with a Dehalococcoides-containing enrichment culture, and varying of the concentration of whey were strategies implemented at the site to enhance dechlorination. Nucleic acid data collected from the groundwater were used to diagnose the in situ physiology of Dehalococcoides spp. and to evaluate the progress of sequential treatments.
|
|
|---|
Periodic whey powder injections as described in Table 1 were performed as part of the treatment strategy within the two treatment plots. During injection, groundwater was extracted from the extraction wells and recirculated through a whey injection system, where whey powder was fed in-line at the desired feed rate, mixed using an in-line actuator, and reinjected into the injection wells to the final concentrations described in Table 1. This procedure typically took 3 h, and the recirculation system was only used during injections as a source of water for the injections. Ambient groundwater flow was the primary mechanism of whey distribution following injections. A groundwater sample was collected from each monitoring well on each treatment plot for molecular analysis prior to biostimulation with whey amendment as a baseline (April 2005), 1 month following whey amendment (July 2005), 1 month following bioaugmentation (August 2005), 4 months following bioaugmentation and modification of whey amendment (November 2005), 7 months following bioaugmentation (February 2006), and finally 9 months following bioaugmentation and 2 months after cessation of whey amendments (April 2006).
|
View this table: [in a new window] |
TABLE 1. Dates and concentration of whey injected at treatment plots 1 and 2
|
Groundwater samples were also collected to measure geochemical composition, organic carbon, redox parameters (data not shown), and chlorinated ethenes. Groundwater samples were collected using dedicated Teflon tubing in each port of the continuous-multichannel-tubing sampling systems and pumping using a peristaltic pump. Purging was performed by pumping at approximately 100 ml/min and until water quality indicator parameters (temperature, pH, conductivity, dissolved oxygen, and oxidation-reduction potential) stabilized based on measurements with a water quality instrument (Quanta G by Hydrolab) and a flow-through cell. For the volatile organic carbon results presented here, samples were sent to a fixed lab (Severn Trent Laboratories) for analysis according to standard EPA method SW846 8260B.
DNA and RNA isolation.
Groundwater samples intended for genomic DNA (gDNA) isolation were collected from monitoring wells (MW) into autoclaved 1-liter bottles and were shipped overnight on ice to the laboratory at the University of California at Berkeley. Samples were stored at 4°C upon arrival and processed within 48 h. Between 100 ml and 1,000 ml of groundwater was filtered through a 0.2-µm, surfactant-free cellulose acetate filter (Fisher, Houston, TX) that was subsequently vortexed vigorously for 5 min in 2 ml of the sample groundwater in order to dislodge filtered particles. The resultant liquid was centrifuged at 21,000 x g at 4°C for 10 min, the supernatant was discarded, and the cell pellet was stored at –80°C until extraction. gDNA was isolated from frozen cell pellets using the UltraClean microbial DNA isolation kit (Mo Bio Laboratories, Carlsbad, CA) according to the manufacturer's instructions. A negative control of 0.2-µm-filtered and autoclaved Milli-Q water was extracted in parallel with each set of samples. The gDNA concentration was quantified using the PicoGreen double-stranded-DNA quantitation kit (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions.
Groundwater samples intended for RNA isolation were filtered on site at the time of sampling. Between 275 ml and 8 liters of groundwater was pressured through a sterile 142-mm-diameter, 0.22-µm hydrophilic polyvinylidene fluoride filter (Millipore, Billerica, MA) via a tripod collection system. The filtration process typically took less than 5 min. The filters were then placed into sterile 50-ml conical tubes and immediately frozen on dry ice. Samples were shipped overnight in dry ice to the University of California at Berkeley and upon arrival stored at –80°C until extraction. Total RNA was isolated from the frozen filters using the RNA PowerSoil total RNA isolation kit (Mo Bio Laboratories, Carlsbad, CA) according to the manufacturer's instructions. Prior to the first step of the extraction procedure, the filters in the 50-ml conical tubes were cut into smaller pieces with scissors while the tubes were submerged in liquid nitrogen. After RNA isolation, contaminating DNA was removed by applying two successive DNase treatments via the DNA-free kit (Ambion, Austin, TX) according to the manufacturer's recommendation. The purified total RNA concentration was quantified with the RiboGreen RNA quantification kit (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. In order to deal with the limitations inherent in imperfect DNA and RNA extraction efficiencies in this study, a single method was applied throughout the study (by the same operator for consistency) and results were compared over time for relative interpretation rather than absolute quantification.
Quantification of genes and gene expression.
To quantify the three identified RDase genes (tceA, vcrA, and bvcA) and the 16S rRNA gene of Dehalococcoides spp., qPCR was applied to the gDNA as described previously (11). Briefly, each 25-µl reaction mixture contained 3 µl of sample or 10-fold serially diluted standard, 1x TaqMan Universal PCR master mix (Applied Biosystems, Foster City, CA), 0.7 µM of the appropriate forward and reverse primers (Table 2), and 0.2 µM of TaqMan probe (Applied Biosystems, Foster City, CA) (Table 2). Standard curves for qPCR were constructed using a single linearized plasmid containing all three RDase genes and the 16S rRNA gene fragments as described by Holmes et al. (11). All samples and standards were analyzed in triplicate, and gene concentrations were reported as copies of gene per liter of groundwater. Selected samples were routinely diluted 10- to 100-fold to test for the presence of PCR inhibitors in the groundwater, and a linear response was always obtained from the dilution series (see Fig. S2 in the supplemental material).
|
View this table: [in a new window] |
TABLE 2. Primer and probe sequences used in this studya
|
RNA standards for the tceA and vcrA transcripts were in vitro transcribed as described previously (13, 17). Standards for the bvcA and 16S rRNA transcripts were synthesized in vitro in this study using the MEGAscript T7 kit (Ambion, Austin, TX) according to the manufacturer's instructions. Input to the RNA synthesis kit was a PCR product that contained the qPCR primers and probe target region and the T7 RNA polymerase promoter sequence. The PCR product for in vitro transcription was synthesized as follows. The T7 sequence was incorporated upstream of the coding sequence via a forward primer that was identical in sequence to the qPCR forward primer except with the T7 sequence added to the 5' end (Table 2). The reverse primer (Table 2) was downstream of the qPCR target region to generate an amplicon of at least 1,000 bp to enhance the efficiency of in vitro transcription. The T7 sequences incorporating the forward primer and the reverse primer were used to amplify gDNA from strain 195 or strain BAV1 in a PCR according to the following protocol: an initial denaturation step at 94°C (12 min), followed by 35 cycles of 94°C (60 s), 50°C (45 s), and 72°C (120 s), with a final extension at 72°C for 12 min. The 30-µl PCR contained 1x PCR buffer II, 2.5 mM MgCl2, 0.13 mg/ml bovine serum albumin, 1 mM of deoxynucleoside triphosphate (dNTP) mix, 0.05 µM of each primer, and 1.5 U of AmpliTaq Gold DNA polymerase (Applied Biosystems, Foster City, CA).
Following the PCR and the subsequent in vitro transcription, the synthesized RNA standard was purified using the MEGAclear kit (Ambion, Austin, TX) according to the manufacturer's instructions. Contaminating DNA within the synthesized transcripts was removed via the DNA-free kit (Ambion, Austin, TX). Complete DNA removal was confirmed by triplicate RT-qPCR reactions where the reverse transcriptase had been omitted and no fluorescence signal was detected after 40 qPCR cycles. The RNA standard concentration was quantified using the RiboGreen RNA quantification kit (Invitrogen, Carlsbad, CA), and the theoretical copy number was calculated according to the size of the input PCR amplicon and an average molecular weight of 340 Da per RNA nucleotide.
RDase gene cDNA clone library.
In order to identify additional highly expressed RDase genes in the field samples, degenerate primers (15) were used to detect expression of putative RDase genes by means of a clone library. Total RNA from sample MW 1D4 (February 2006) was reverse transcribed using the SuperScript III First-Strand synthesis system (Invitrogen, Carlsbad, CA) and PCR amplified. In a reverse-transcription reaction of 10 µl, 2 µl of sample RNA was first incubated with 50 ng of random hexamers and 1 mM of dNTP mix at 65°C for 5 min and later on ice for 1 min. A 10-µl cDNA synthesis mix containing 2x RT buffer, 10 mM MgCl2, 0.02 M dithiothreitol, 40 U RNaseOUT, and 200 U of reverse transcriptase was then combined with the previous 10-µl sample and incubated at 25°C for 20 min, 50°C for 60 min, and 65°C for 5 min. Parallel reactions with no reverse transcriptase and with no template were prepared as negative controls. Following reverse transcription, 2 µl of the cDNA product was used in a 25-µl PCR that contained 1x Ex Taq buffer, 0.8 mM dNTP mix, 0.5 µM (each) of the degenerate primers (15), and 0.625 U Ex Taq hot-start DNA polymerase (Takara Mirus Bio, Madison, WI). The same thermal cycling conditions as above were applied except that 40 cycles were used and the annealing temperature was 48°C. For sample MW 1D4, the expected PCR product of
1,700 bp was detected on an Agilent 2100 bioanalyzer (Agilent, Palo Alto, CA) using DNA 7500 chips according to the manufacturer's protocol, and the negative controls showed no band.
The clone library was constructed with the TOPO TA cloning kit (with the pCR2.1-TOPO vector) (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. A total of 75 plasmid inserts were later digested with the restriction enzymes HhaI or MspI (Promega, Madison, WI) at 34°C for 3 h, and the resulting fragments were analyzed on an Agilent 2100 bioanalyzer (Agilent, Palo Alto, CA) using DNA 7500 chips. Distinct clones were selected for sequencing at Elim Biopharmaceuticals, Inc. (Hayward, CA). Primers for sequencing of internal regions of a clone were designed using the Primer Express v2.0 software program (Applied Biosystems), and the sequenced reads were assembled using the Cap3 Assembly software program (http://pbil.univ-lyon1.fr/cap3.php). After alignment and visual checking, the near-full-length RDase sequences were compared to sequences in GenBank using BLASTn and BLASTx (http://www.ncbi.nlm.nih.gov/).
Quantification of the FtL-RDase-1638 gene.
Two primer and probe sets (Table 2) were used to target the front and rear regions of the FtL-RDase-1638 gene retrieved from the cDNA clone library. Two distinct regions were targeted in order to obtain an accurate placement of this gene, since a nucleotide BLAST analysis against the genome of strain BAV1 using the IMG 2.3 software at JGI (http://www.jgi.doe.gov/) indicated a possible truncation of an FtL-RDase-1638-like gene by an integrase in that isolate. Standards for the qPCR analysis of the FtL-RDase-1638I and FtL-RDase-1638II gene segments were synthesized by IDT (Coralville, IA) as two separate 80-bp double-stranded DNAs that contained the primer and probe target sequence. The theoretical copy number was calculated using the length of the duplex and an average molecular weight of 660 Da per DNA nucleotide.
Nucleotide sequence accession numbers.
Sequences from the RDase gene cDNA clone library have been submitted to the GenBank database under accession no. EU137841 to EU137844.
|
|
|---|
![]() View larger version (32K): [in a new window] |
FIG. 1. Molar concentrations of chlorinated ethenes and ethene in groundwater samples over the 1-year monitoring period for treatment plots 1 and 2. Concentrations are reported as averages for eight sampling locations within each treatment plot, and error bars represent analytical error. The date and concentration of whey injected into the respective treatment plots are indicated on the top axis, and manipulations implemented at the site are indicated by arrows on the graphs. For clarity, only one of the two sampling events performed within 2 weeks of each other during March 2005 is plotted on the graph.
|
Quantifying RDase and 16S rRNA genes.
qPCR was used to monitor the concentrations of 16S rRNA, tceA, bvcA, and vcrA genes of Dehalococcoides spp. in treatment plots 1 and 2 over the 1-year monitoring period spanning biostimulation and bioaugmentation activities (Fig. 2). In order to focus on the temporal trends of the data, the respective gene concentrations were averaged spatially across the eight sampling wells on each treatment plot. At the baseline sampling event (April 2005), the Dehalococcoides 16S rRNA gene was detected at 1.2 x 104 ± 0.20 x 104 copies/liter and 1.1 x 105 ± 0.13 x 105 copies/liter in plots 1 and 2, respectively (Fig. 2). Redox conditions at the site during the baseline sampling event indicated predominantly aerobic conditions, with plot 2 exhibiting slightly more reducing conditions than plot 1 (data not shown). Genes associated with VC degradation, bvcA and vcrA, were present in both treatment plots 1 and 2, while the tceA gene was detectable only in plot 2 initially.
![]() View larger version (47K): [in a new window] |
FIG. 2. Dynamics of the Dehalococcoides (Dhc) 16S rRNA gene and RDase gene concentrations over the 1-year monitoring period for treatment plots 1 and 2. Manipulations implemented at the site are indicated by arrows on the graphs. Data at each time point are averages for samples from the eight monitoring wells in the respective plot, and each error bar represents one standard deviation for the qPCR method.
|
Following the July 2005 sampling, bioaugmentation was implemented at the site with an enrichment culture that had a Dehalococcoides 16S rRNA gene concentration of 4.3 x 107 ± 0.72 x 107 copy/ml, and it contained the tceA (5.1 x 107 ± 0.28 x 107 copy/ml) and vcrA (3.6 x 107 ± 0.12 x 107 copy/ml) genes but no bvcA gene. For 3 months after bioaugmentation (August to November 2005), while all monitored genes increased to some extent, the most significant change was the order-of-magnitude increase in the vcrA gene in both plots (Fig. 2).
By February 2006, 3 months after the whey injection strategy was altered, the 16S rRNA gene concentrations in plot 1 increased by 2 orders of magnitude while those in plot 2 increased by 19-fold (Fig. 2). In parallel with the significant increase in the 16S rRNA gene was the increase in the three RDase genes, with the bvcA gene concentration continuing to dominate in numbers. During the 1-year biostimulation-and-bioaugmentation period, the 16S rRNA gene concentrations increased by more than 3 orders of magnitude in plots 1 and 2 (Fig. 2).
Accounting for the dominant Dehalococcoides cells.
The 16S rRNA gene sequences of distinct Dehalococcoides strains are highly conserved, making it impossible to differentiate metabolic functions using this gene. Primers and probes that target the 16S rRNA gene in this study were designed to target 16S rRNA sequences of all known Dehalococcoides strains. The correlation between concentrations of 16S rRNA gene copies and the sum total of the tceA, vcrA, and bvcA genes is plotted in Fig. 3 in an attempt to identify the dominant Dehalococcoides cells within the populations. As illustrated in the regression analysis in Fig. 3, there was good linear correlation (r2 = 0.99) between the 16S rRNA gene and the sum total of the three RDase gene concentrations, and the two variables were within a factor of 1.4 ± 0.064 (slope of regression line ± 95% confidence interval) of each other over concentrations that span several orders of magnitude. With the exception of a few samples analyzed during the July 2005 sampling in plot 2, the data were distributed randomly around the 1:1 ratio and the differences were within the analytical reproducibility and uncertainty of the qPCR technique. The close tracking of copy numbers between functional and ribosomal genes was observed both during the early phases of treatment, when bvcA was the only dominant RDase gene, and during the later phases of treatment, when both the vcrA and bvcA genes were dominant.
![]() View larger version (16K): [in a new window] |
FIG. 3. The sums of RDase gene concentrations are compared against the Dehalococcoides (Dhc) 16S rRNA gene concentrations on a log-log scale. Each data point represents a sample collected at a monitoring well from treatment plots 1 and 2 during each sampling event. Each error bar represents one standard deviation for the qPCR method. The dashed line represents the hypothetical 1:1 correlation between the two variables.
|
![]() View larger version (47K): [in a new window] |
FIG. 4. (A) Expression profile of the Dehalococcoides (Dhc) 16S rRNA gene and the RDase genes during the February 2006 and April 2006 sampling events. Data were calculated from triplicate RT-qPCRs, and each error bar represents one standard deviation. Labels on the x axis represent monitoring well designations. (B) Gene concentrations at the corresponding monitoring wells. Data were calculated from triplicate qPCRs, and each error bar represents one standard deviation.
|
|
View this table: [in a new window] |
TABLE 3. RDase transcripts recovered from a cDNA-clone library
|
![]() View larger version (40K): [in a new window] |
FIG. 5. Concentrations of the FtL-RDase-1638 gene are compared against Dehalococcoides (Dhc) 16S rRNA gene and RDase gene concentrations. The data represent two sampling events for monitoring well 2B4. FtL-RDase-1638I and FtL-RDase-1638II represent regions from bp 381 to 449 and from bp 1056 to 1129, respectively, on the sequenced clone. Data were calculated from triplicate qPCRs, and each error bar represents one standard deviation.
|
|
|
|---|
The total number of Dehalococcoides 16S rRNA gene copies present in the field samples was compared to the sum of identified RDase gene copies (tceA, bvcA, and vcrA) in order to determine whether the dominant Dehalococcoides cells present at the site could be identified by quantification of the three identified RDase genes. The accuracy of the qPCR comparison was facilitated by the application of a four-gene plasmid standard that has been shown to minimize analytical variability among individual genes within single samples (11). Because of the highly conserved nature of the 16S rRNA gene among Dehalococcoides strains and its single copy in the genome, it can be assumed that the measured 16S rRNA gene concentration represents the total Dehalococcoides concentration within a sample. Furthermore, it was assumed that each of the three identified RDase genes was present as a single copy in the genome, since this has been consistently observed in previously characterized Dehalococcoides isolates (9, 15, 27, 32; http://www.jgi.doe.gov/). Hence, significant discrepancies between the sum total of the three RDase genes and the 16S rRNA gene concentrations would then suggest the presence of Dehalococcoides cells not detectable by bvcA-, vcrA-, and tceA-targeted primers/probes. In this study, the close agreement observed between the 16S rRNA gene concentration and the RDase gene sum across several orders of magnitude of detection (Fig. 3) suggests that the majority of Dehalococcoides cells at this field site carry a tceA, vcrA, or bvcA gene and that the number of Dehalococcoides cells with none of these RDase genes is relatively insignificant. Furthermore, the three RDase genes were determined to be present in separate Dehalococcoides cells, since independent variation in the ratio of gene concentration over time indicates that colocation of any of the genes is unlikely (Fig. 2). While a complete accounting of Dehalococcoides cells by quantifying known RDase genes was possible at this site and for other enrichment cultures (11, 31), this might not always be the case, as has been demonstrated with other environmental samples (31).
The simultaneous tracking of three Dehalococcoides RDases at this site also provided insights into their potential interactions under field conditions. In general, Dehalococcoides populations also compete against other TCE-respiring organisms outside of the Dehalococcoides genus that are only able to carry out incomplete TCE respiration to cDCE (33). In fact, a previous modeling study demonstrated that growth of Dehalococcoides populations could be restricted by other PCE- or TCE-to-cDCE dehalorespirers under certain growth conditions (2), with the competitive advantage of Dehalococcoides populations lying in their ability to metabolize cDCE and VC.
Dynamically, while the entire Dehalococcoides population showed a significant increase over time, the extent of the increase was unique for each targeted functional gene (Fig. 2). For both treatment plots, the tceA-carrying cells showed about a 2-orders-of-magnitude increase throughout the study but never represented more than 0.2 to 4.5% of the total Dehalococcoides population even though significant reduction of TCE was measured (Fig. 2), suggesting that these cells were outcompeted by other TCE-respiring organisms and other Dehalococcoides strains. In contrast, the vcrA-carrying cells in plot 1 increased by more than 4 orders of magnitude over the 1-year monitoring period, while the bvcA-carrying cells exhibited a more than 3-orders-of-magnitude increase (Fig. 2). The ratio of bvcA-to-vcrA cells changed from the initial value of 4:1 to 1:1 toward the end of the monitoring period in plot 1. For plot 2, both cells increased by more than 3 orders of magnitude and the 2:1 ratio that was initially observed between the bvcA- and vcrA-carrying cells was maintained toward the end. The tceA-, bvcA-, and vcrA-carrying cells occupied a unique niche in this subsurface environment as the only currently known organisms that can reduce chlorinated ethenes beyond DCE (33), and their significant overall increases were reflected in the production of VC and ethene during the latter phases of treatment.
In addition to biostimulation, bioaugmentation with a Dehalococcoides-containing culture was also implemented at this site. The rationale for implementing bioaugmentation was to ensure that the field site contained the appropriate metabolic capability for complete reductive dechlorination to ethene. In hindsight, this strategy might not have been necessary, since the three characterized genes that are vital for converting chlorinated ethenes to ethene were indigenous to the site. The addition of the tceA- and vcrA-containing cells from the bioaugmentation culture created a redundancy in metabolic functions, and it was not possible to differentiate the contribution from the two sources to the overall gene concentration in order to specifically evaluate the performance of the bioaugmentation culture in the field environment. In contrast, because the bvcA gene was absent in the bioaugmentation culture, the measured increase in that gene concentration could be solely attributed to biostimulation of the indigenous cells.
Although net increases in 16S rRNA gene concentrations are indicative of increased Dehalococcoides cell numbers, the presence of a gene does not always guarantee expression (19). Further, changes in the gene concentration from net cell growth require considerable amounts of time for Dehalococcoides cells, as evidenced by doubling times reported under optimized laboratory conditions of 19 h for strain 195 (25) and 2.2 days for strain BAV1 (8). The doubling time in the less-controlled field conditions might be considerably longer. Laboratory experiments have also indicated that the half-life of Dehalococcoides RDase transcripts is on the order of 4 to 6 h in the absence of chlorinated ethenes (17) and the up-regulation response time after chlorinated ethenes exposure is approximately 3 h (14). This relatively short response time is ideal for diagnosing cell physiology.
Although a significant increase in the overall gene concentration was observed at the field site during the 1-year monitoring period, there were two periods during which significant growth was not detected (between July and November 2005 and between February and April 2006). However, the transcript data measured for the RDase genes and 16S rRNA gene during one of those periods verified that Dehalococcoides cells were indeed physiologically active even though no significant net increase in gene concentration was observed (Fig. 4). The quantity of the 16S rRNA transcripts was comparable to that for actively dechlorinating D. ethenogenes strain 195 cells grown in the laboratory (P. K. H. Lee and L. Alvarez-Cohen, unpublished data). Furthermore, the expression data indicated variable levels of activity among the Dehalococcoides populations. While the concentration of the bvcA gene was the highest in groundwater samples, its expression level was unexpectedly lower than those of both the vcrA and tceA genes when the transcripts were normalized to the respective gene copy number. Transcript-per-gene data normalize for the effects of cell concentration on the expression profile. The bvcA expression data suggested that not all bvcA genes were actually in metabolically active cells. On the other hand, the higher vcrA transcript-per-gene levels indicated that the vcrA-carrying cells were relatively more active. The expression of the bvcA and vcrA transcripts was consistent with the chlorinated ethene data, which showed significant increases in VC and ethene during the February and April 2006 samplings. Interestingly, the tceA transcripts were not uniformly detected in all samples that contained the tceA gene. The lack of expression suggests that not all tceA-carrying cells were metabolically active and they might not have played a major role in the TCE-to-VC reduction. These data could further explain why the tceA strain lagged behind the other Dehalococcoides cells in overall cell numbers.
In our examination of expressed RDase genes at the site detected using the cDNA clone library, it was surprising to recover a highly expressed gene (FtL-RDase-1638) that shared a high percent identity (99%) with a putative RDase gene in strain CBDB1 (Table 3). The close accounting between Dehalococcoides 16S rRNA gene copies and the three targeted RDase genes suggested that it was unlikely that strain CBDB1, an isolate that does not carry the tceA, vcrA, or bvcA gene targeted here, represented a significant fraction of the Dehalococcoides population (Fig. 3). Quantitative analysis of the copy numbers of the FtL-RDase-1638 gene using two distinct regions of the sequence suggested that this gene is likely present in the genomes of bvcA-carrying cells at this site (Fig. 5), with the possibility that it is also present in other low-abundance Dehalococcoides cells. The discrepancy of less than a factor of 3.5 between the absolute quantities of the FtL-RDase-1638 gene fragments is probably due to the application of two separately quantified standards in qPCR analysis (11). Clone libraries of a Dehalococcoides-containing enrichment, KB-1, also indicated the presence and expression of a gene (KB1RdhAB5) similar to FtL-RDase-1638 when KB-1 was exposed to TCE, cDCE, or VC (37). While the genomic location of this gene was not determined, bvcA-carrying cells are present in KB-1. Interestingly, although an FtL-RDase-1638-like gene was not experimentally identified as 1 of the 11 putative RDase genes in strain BAV1 (15), a nucleotide-nucleotide BLAST analysis against the genome of strain BAV1 identified a 645-bp fragment that has 95% similarity to the amino terminus-coding region of the FtL-RDase-1638 gene (http://www.jgi.doe.gov/). According to the annotation of strain BAV1, this 645-bp fragment is followed downstream by an integrase and a transposase (see URL above), indicating that this RDase gene has been truncated in this strain. Based on the data in this study together with the KB-1 and strain BAV1 results, the FtL-RDase-1638 gene is likely present in the bvcA-carrying genome.
In summary, this study demonstrates that quantifying genes and transcripts of the 16S rRNA gene and three characterized RDase genes of Dehalococcoides during different stages of bioremediation can provide important and complementary information about the dynamics and activity of these organisms. Since Dehalococcoides spp. are the only organisms known to convert TCE to innocuous ethene, having techniques that can monitor and optimize their activity can substantially promote the success of bioremediation projects. The type of data collected in this study can also be useful to regulatory agencies to evaluate the effectiveness of in situ biological processes.
This research was supported by Environmental Security Technology Certification Program (ESTCP) projects ER-0218 and ER-0318, National Science Foundation grant no. BES-01-0504244, and NIEHS Superfund Basic Research Project ES04705-19.
Published ahead of print on 7 March 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»