Detection of Organohalide-Respiring Enzyme Biomarkers at a Bioaugmented TCE-Contaminated Field Site.

RNA-based biomarkers have been successfully detected at field sites undergoing in situ bioremediation, but the detection of expressed enzymes is a more direct way to prove activity for a particular biocatalytic process of interest since they provide evidence of potential in situ activity rather than simply confirming presence and abundance of genes in a given population by measurement of DNA copies using qPCR. Here we successfully applied shotgun proteomics to field samples from a trichloroethene (TCE)-contaminated industrial site in southern Ontario, Canada that had been bio-augmented with the commercially available KB-1TM microbial culture. The KB-1TM culture contains multiple strains of Dehalococcoides mccartyi (D. mccartyi) as well as an organohalide respiring Geobacter species. The relative abundances of specific enzymatic proteins were subsequently compared to corresponding qPCR-derived levels of DNA and RNA biomarkers in the same samples. Samples were obtained from two wells with high hydraulic connectivity to the KB-1TM-bioaugemented enhanced in situ bioremediation system, and two control wells that showed evidence of low levels of native organohalide respiring bacteria (OHRB), Dehalococcoides and Geobacter. Enzymes involved in organohalide respiration were detected in the metaproteomes of all four field samples, as were chaperonins of D. mccartyi, chemotaxis proteins, and ATPases. The most highly expressed RDase in the bioaugmentation culture (VcrA) was the most highly detected enzyme overall in the bioaugmented groundwater samples. In one background groundwater well, we found high expression of the Geobacter pceA RDase. The DNA and RNA biomarkers detected using qPCR-based assays were a set of orthologs of Dehalococcoides reductive dehalogenases (VcrA, TceA, BvcA, dehalogenase "DET1545"), and the Ni-Fe uptake hydrogenase, HupL. Within a sample, RNA levels for key enzymes correlated with relative protein abundance. These results indicate that laboratory observations of TCE-bioremediation biomarker protein expression are recapitulated in field environmental systems and that both RNA and protein biomarker monitoring hold promise for activity monitoring of in situ populations of OHRB.

Biomarkers' Responses to Reductive Dechlorination Rates and Oxygen Stress in Bioaugmentation Culture KB-1TM.

Using mRNA transcript levels for key functional enzymes as proxies for the organohalide respiration (OHR) rate, is a promising approach for monitoring bioremediation populations in situ at chlorinated solvent-contaminated field sites. However, to date, no correlations have been empirically derived for chlorinated solvent respiring, Dehalococcoides mccartyi (DMC) containing, bioaugmentation cultures. In the current study, genome-wide transcriptome and proteome data were first used to confirm the most highly expressed OHR-related enzymes in the bioaugmentation culture, KB-1TM, including several reductive dehalogenases (RDases) and a Ni-Fe hydrogenase, Hup. Different KB-1™ DMC strains could be resolved at the RNA and protein level through differences in the sequence of a common RDase (DET1545-like homologs) and differences in expression of their vinyl chloride-respiring RDases. The dominant strain expresses VcrA, whereas the minor strain utilizes BvcA. We then used quantitative reverse-transcriptase PCR (qRT-PCR) as a targeted approach for quantifying transcript copies in the KB-1TM consortium operated under a range of TCE respiration rates in continuously-fed, pseudo-steady-state reactors. These candidate biomarkers from KB-1TM demonstrated a variety of trends in terms of transcript abundance as a function of respiration rate over the range: 7.7 × 10-12 to 5.9 × 10-10 microelectron equivalents per cell per hour (µeeq/cell∙h). Power law trends were observed between the respiration rate and transcript abundance for the main DMC RDase (VcrA) and the hydrogenase HupL (R² = 0.83 and 0.88, respectively), but not transcripts for 16S rRNA or three other RDases examined: TceA, BvcA or the RDase DET1545 homologs in KB1TM. Overall, HupL transcripts appear to be the most robust activity biomarker across multiple DMC strains and in mixed communities including DMC co-cultures such as KB1TM. The addition of oxygen induced cell stress that caused respiration rates to decline immediately (>95% decline within one hour). Although transcript levels did decline, they did so more slowly than the respiration rate observed (transcript decay rates between 0.02 and 0.03 per hour). Data from strain-specific probes on the pangenome array strains suggest that a minor DMC strain in KB-1™ that harbors a bvcA homolog preferentially recovered following oxygen stress relative to the dominant, vcrA-containing strain.

Meta-analyses of Dehalococcoides mccartyi strain 195 transcriptomic profiles identify a respiration rate-related gene expression transition point and interoperon recruitment of a key oxidoreductase subunit.

A cDNA-microarray was designed and used to monitor the transcriptomic profile of Dehalococcoides mccartyi strain 195 (in a mixed community) respiring various chlorinated organics, including chloroethenes and 2,3-dichlorophenol. The cultures were continuously fed in order to establish steady-state respiration rates and substrate levels. The organization of array data into a clustered heat map revealed two major experimental partitions. This partitioning in the data set was further explored through principal component analysis. The first two principal components separated the experiments into those with slow (1.6±0.6 µM Cl-/h)- and fast (22.9±9.6 µM Cl-/h)-respiring cultures. Additionally, the transcripts with the highest loadings in these principal components were identified, suggesting that those transcripts were responsible for the partitioning of the experiments. By analyzing the transcriptomes (n=53) across experiments, relationships among transcripts were identified, and hypotheses about the relationships between electron transport chain members were proposed. One hypothesis, that the hydrogenases Hup and Hym and the formate dehydrogenase-like oxidoreductase (DET0186-DET0187) form a complex (as displayed by their tight clustering in the heat map analysis), was explored using a nondenaturing protein separation technique combined with proteomic sequencing. Although these proteins did not migrate as a single complex, DET0112 (an FdhB-like protein encoded in the Hup operon) was found to comigrate with DET0187 rather than with the catalytic Hup subunit DET0110. On closer inspection of the genome annotations of all Dehalococcoides strains, the DET0185-to-DET0187 operon was found to lack a key subunit, an FdhB-like protein. Therefore, on the basis of the transcriptomic, genomic, and proteomic evidence, the place of the missing subunit in the DET0185-to-DET0187 operon is likely filled by recruiting a subunit expressed from the Hup operon (DET0112).

Molecular biomarker-based biokinetic modeling of a PCE-dechlorinating and methanogenic mixed culture.

Bioremediation of chlorinated ethenes via anaerobic reductive dechlorination relies upon the activity of specific microbial populations--most notably Dehalococcoides (DHC) strains. In the lab and field Dehalococcoides grow most robustly in mixed communities which usually contain both fermenters and methanogens. Recently, researchers have been developing quantitative molecular biomarkers to aid in field site diagnostics and it is hoped that these biomarkers could aid in the modeling of anaerobic reductive dechlorination. A comprehensive biokinetic model of a community containing Dehalococcoides mccartyi (formerly D. ethenogenes) was updated to describe continuously fed reactors with specific biomass levels based on quantitative PCR (qPCR)-based population data (DNA and RNA). The model was calibrated and validated with subsets of chemical and molecular biological data from various continuous feed experiments (n = 24) with different loading rates of the electron acceptor (1.5 to 482 µeeq/L-h), types of electron acceptor (PCE, TCE, cis-DCE) and electron donor to electron acceptor ratios. The resulting model predicted the sum of dechlorination products vinyl chloride (VC) and ethene (ETH) well. However, VC alone was under-predicted and ETH was over predicted. Consequently, competitive inhibition among chlorinated ethenes was examined and then added to the model. Additionally, as 16S rRNA gene copy numbers did not provide accurate model fits in all cases, we examined whether an improved fit could be obtained if mRNA levels for key functional enzymes could be used to infer respiration rates. The resulting empirically derived mRNA "adjustment factors" were added to the model for both DHC and the main methanogen in the culture (a Methanosaeta species) to provide a more nuanced prediction of activity. Results of this study suggest that at higher feeding rates competitive inhibition is important and mRNA provides a more accurate indicator of a population's instantaneous activity than 16S rRNA gene copies alone as biomass estimates.