Assessment of anaerobic toluene biodegradation activity by bssA transcript/gene ratios

Genetic redundancy in the naphthalene-degradation pathway of Cycloclasticus pugetii strain PS-1 enables response to varying substrate concentrations

Polycyclic aromatic hydrocarbon (PAH) contamination in marine environments range from low-diffusive inputs to high loads. The influence of PAH concentration on the expression of functional genes [e.g. those encoding ring-hydroxylating dioxygenases (RHDs)] has been overlooked in PAH biodegradation studies. However, understanding marker-gene expression under different PAH loads can help to monitor and predict bioremediation efficiency. Here, we followed the expression (via RNA sequencing) of Cycloclasticus pugetii strain PS-1 in cell suspension experiments under different naphthalene (100 and 30 mg L-1) concentrations. We identified genes encoding previously uncharacterized RHD subunits, termed rhdPS1α and rhdPS1β, that were highly transcribed in response to naphthalene-degradation activity. Additionally, we identified six RHD subunit-encoding genes that responded to naphthalene exposure. By contrast, four RHD subunit genes were PAH-independently expressed and three other RHD subunit genes responded to naphthalene starvation. Cycloclasticus spp. could, therefore, use genetic redundancy in key PAH-degradation genes to react to varying PAH loads. This genetic redundancy may restrict the monitoring of environmental hydrocarbon-degradation activity using single-gene expression. For Cycloclasticus pugetii strain PS-1, however, the newly identified rhdPS1α and rhdPS1β genes might be potential target genes to monitor its environmental naphthalene-degradation activity.

Substrate-independent expression of key functional genes in Cycloclasticus pugetii strain PS-1 limits their use as markers for PAH biodegradation

Microbial degradation of petroleum hydrocarbons is a crucial process for the clean-up of oil-contaminated environments. Cycloclasticus spp. are well-known polycyclic aromatic hydrocarbon (PAH) degraders that possess PAH-degradation marker genes including rhd3α, rhd2α, and pahE. However, it remains unknown if the expression of these genes can serve as an indicator for active PAH degradation. Here, we determined transcript-to-gene (TtG) ratios with (reverse transcription) qPCR in cultures of Cycloclasticus pugetii strain PS-1 grown with naphthalene, phenanthrene, a mixture of these PAHs, or alternate substrates (i.e., no PAHs). Mean TtG ratios of 1.99 × 10^-2, 1.80 × 10^-3, and 3.20 × 10^-3 for rhd3α, rhd2α, and pahE, respectively, were measured in the presence or absence of PAHs. The TtG values suggested that marker-gene expression is independent of PAH degradation. Measurement of TtG ratios in Arctic seawater microcosms amended with water-accommodated crude oil fractions, and incubated under in situ temperature conditions (i.e., 1.5°C), only detected Cycloclasticus spp. rhd2α genes and transcripts (mean TtG ratio of 4.15 × 10^-1). The other marker genes-rhd3α and pahE-were not detected, suggesting that not all Cycloclasticus spp. carry these genes and a broader yet-to-be-identified repertoire of PAH-degradation genes exists. The results indicate that the expression of PAH marker genes may not correlate with PAH-degradation activity, and transcription data should be interpreted cautiously.

Does It Pay Off to Explicitly Link Functional Gene Expression to Denitrification Rates in Reaction Models?

Environmental omics and molecular-biological data have been proposed to yield improved quantitative predictions of biogeochemical processes. The abundances of functional genes and transcripts relate to the number of cells and activity of microorganisms. However, whether molecular-biological data can be quantitatively linked to reaction rates remains an open question. We present an enzyme-based denitrification model that simulates concentrations of transcription factors, functional-gene transcripts, enzymes, and solutes. We calibrated the model using experimental data from a well-controlled batch experiment with the denitrifier Paracoccous denitrificans. The model accurately predicts denitrification rates and measured transcript dynamics. The relationship between simulated transcript concentrations and reaction rates exhibits strong non-linearity and hysteresis related to the faster dynamics of gene transcription and substrate consumption, relative to enzyme production and decay. Hence, assuming a unique relationship between transcript-to-gene ratios and reaction rates, as frequently suggested, may be an erroneous simplification. Comparing model results of our enzyme-based model to those of a classical Monod-type model reveals that both formulations perform equally well with respect to nitrogen species, indicating only a low benefit of integrating molecular-biological data for estimating denitrification rates. Nonetheless, the enzyme-based model is a valuable tool to improve our mechanistic understanding of the relationship between biomolecular quantities and reaction rates. Furthermore, our results highlight that both enzyme kinetics (i.e., substrate limitation and inhibition) and gene expression or enzyme dynamics are important controls on denitrification rates.

Combined Use of Diagnostic Fumarate Addition Metabolites and Genes Provides Evidence for Anaerobic Hydrocarbon Biodegradation in Contaminated Groundwater

The widespread use of hydrocarbon-based fuels has led to the contamination of many natural environments due to accidental spills or leaks. While anaerobic microorganisms indigenous to many fuel-contaminated groundwater sites can play a role in site remediation (e.g., monitored natural attenuation, MNA) via hydrocarbon biodegradation, multiple lines of evidence in support of such bioremediation are required. In this study, we investigated two fuel-contaminated groundwater sites for their potential to be managed by MNA. Microbial community composition, biogeochemical indicators, fumarate addition metabolites, and genes diagnostic of both alkane and alkyl-monoaromatic hydrocarbon activation were assessed. Fumarate addition metabolites and catabolic genes were detected for both classes of hydrocarbon biodegradation at both sites, providing strong evidence for in situ anaerobic hydrocarbon biodegradation. However, relevant metabolites and genes did not consistently co-occur within all groundwater samples. Using newly designed mixtures of quantitative polymerase chain reaction (qPCR) primers to target diverse assA and bssA genes, we measured assA gene abundances ranging from 10⁵–10⁸ copies/L, and bssA gene abundances ranging from 10⁵–10¹⁰ copies/L at the sites. Overall, this study demonstrates the value of investigating fuel-contaminated sites using both metabolites and genes diagnostic of anaerobic hydrocarbon biodegradation for different classes of hydrocarbons to help assess field sites for management by MNA.

Benzylsuccinate Synthase is Post-Transcriptionally Regulated in the Toluene-Degrading Denitrifier Magnetospirillum sp. Strain 15-1

The facultative denitrifying alphaproteobacterium Magnetospirillum sp. strain 15-1 had been isolated from the hypoxic rhizosphere of a constructed wetland model fed with toluene. This bacterium can catabolize toluene anaerobically but not aerobically. Here, we used strain 15-1 to investigate regulation of expression of the highly oxygen-sensitive glycyl radical enzyme benzylsuccinate synthase, which catalyzes the first step in anaerobic toluene degradation. In cells growing aerobically with benzoate, the addition of toluene resulted in a ~20-fold increased transcription of bssA, encoding for the catalytically active subunit of the enzyme. Under anoxic conditions, bssA mRNA copy numbers were up to 129-fold higher in cells growing with toluene as compared to cells growing with benzoate. Proteomics showed that abundance of benzylsuccinate synthase increased in cells growing anaerobically with toluene. In contrast, peptides of this enzyme were never detected in oxic conditions. These findings show that synthesis of benzylsuccinate synthase was under stringent post-transcriptional control in the presence of oxygen, which is a novel level of regulation for glycyl radical enzymes.

Toluene biodegradation in the vadose zone of a poplar phytoremediation system identified using metagenomics and toluene-specific stable carbon isotope analysis

Biodegradation is an important mechanism of action of phytoremediation systems, but performance evaluation is challenging. We applied metagenomic molecular approaches and compound-specific stable carbon isotope analysis to assess biodegradation of toluene in the vadose zone at an urban pilot field system where hybrid poplars were planted to remediate legacy impacts to an underlying shallow fractured bedrock aquifer. Carbon isotope ratios were compared spatio-temporally between toluene dissolved in groundwater and in the vapor phase. Enrichment of ¹³C from toluene in the vapor phase compared to groundwater provided evidence for biodegradation in the vadose zone. Total bacterial abundance (16S rRNA) and abundance and expression of degradation genes were determined in rhizosphere soil (DNA and RNA) and roots (DNA) using quantitative PCR. Relative abundances of degraders in the rhizosphere were on average higher at greater depths, except for enrichment of PHE-encoding communities that more strongly followed patterns of toluene concentrations detected. Quantification of RMO and PHE gene transcripts supported observations of active aerobic toluene degradation. Finally, spatially-variable numbers of toluene degraders were detected in poplar roots. We present multiple lines of evidence for biodegradation in the vadose zone at this site, contributing to our understanding of mechanisms of action of the phytoremediation system.

Transcriptional Regulation of the Peripheral Pathway for the Anaerobic Catabolism of Toluene and m-Xylene in Azoarcus sp. CIB

Alkylbenzenes, such as toluene and m-xylene, are an important class of contaminant hydrocarbons that are widespread and tend to accumulate in subsurface anoxic environments. The peripheral pathway for the anaerobic oxidation of toluene in bacteria consists of an initial activation catalyzed by a benzylsuccinate synthase (encoded by bss genes), and a subsequent modified β-oxidation of benzylsuccinate to benzoyl-CoA and succinyl-CoA (encoded by bbs genes). We have shown here that the bss and bbs genes, which are located within an integrative and conjugative element, are essential for anaerobic degradation of toluene but also for m-xylene oxidation in the denitrifying beta-proteobacterium Azoarcus sp. CIB. New insights into the transcriptional organization and regulation of a complete gene cluster for anaerobic catabolism of toluene/m-xylene in a single bacterial strain are presented. The bss and bbs genes are transcriptionally coupled into two large convergent catabolic operons driven by the PbssD and PbbsA promoters, respectively, whose expression is inducible when cells grow anaerobically in toluene or m-xylene. An adjacent tdiSR operon driven by the PtdiS promoter encodes a putative two-component regulatory system. TdiR behaves as a transcriptional activator of the PbssD, PbbsA, and PtdiS promoters, being benzylsuccinate/(3-methyl)benzylsuccinate, rather than toluene/m-xylene, the inducers that may trigger the TdiS-mediated activation of TdiR. In addition to the TdiSR-based specific control, the expression of the bss and bbs genes in Azoarcus sp. CIB is under an overimposed regulation that depends on certain environmental factors, such as the presence/absence of oxygen or the availability of preferred carbon sources (catabolite repression). This work paves the way for future strategies toward the reliable assessment of microbial activity in toluene/m-xylene contaminated environments.

The sequence capture by hybridization: a new approach for revealing the potential of mono-aromatic hydrocarbons bioattenuation in a deep oligotrophic aquifer

The formation water of a deep aquifer (853 m of depth) used for geological storage of natural gas was sampled to assess the mono-aromatic hydrocarbons attenuation potential of the indigenous microbiota. The study of bacterial diversity suggests that Firmicutes and, in particular, sulphate-reducing bacteria (Peptococcaceae) predominate in this microbial community. The capacity of the microbial community to biodegrade toluene and m- and p-xylenes was demonstrated using a culture-based approach after several hundred days of incubation. In order to reveal the potential for biodegradation of these compounds within a shorter time frame, an innovative approach named the solution hybrid selection method, which combines sequence capture by hybridization and next-generation sequencing, was applied to the same original water sample. The bssA and bssA-like genes were investigated as they are considered good biomarkers for the potential of toluene and xylene biodegradation. Unlike a PCR approach which failed to detect these genes directly from formation water, this innovative strategy demonstrated the presence of the bssA and bssA-like genes in this oligotrophic ecosystem, probably harboured by Peptococcaceae. The sequence capture by hybridization shows significant potential to reveal the presence of genes of functional interest which have low-level representation in the biosphere.

The ecology of anaerobic degraders of BTEX hydrocarbons in aquifers

The degradation of benzene, toluene, ethylbenzene and xylene (BTEX) contaminants in groundwater relies largely on anaerobic processes. While the physiology and biochemistry of selected relevant microbes have been intensively studied, research has now started to take the generated knowledge back to the field, in order to trace the populations truly responsible for the anaerobic degradation of BTEX hydrocarbons in situ and to unravel their ecology in contaminated aquifers. Here, recent advances in our knowledge of the identity, diversity and ecology of microbes involved in these important ecosystem services are discussed. At several sites, distinct lineages within the Desulfobulbaceae, the Rhodocyclaceae and the Gram-positive Peptococcaceae have been shown to dominate the degradation of different BTEX hydrocarbons. Especially for the functional guild of anaerobic toluene degraders, specific molecular detection systems have been developed, allowing researchers to trace their diversity and distribution in contaminated aquifers. Their populations appear enriched in hot spots of biodegradation in situ. ¹³C-labelling experiments have revealed unexpected pathways of carbon sharing and obligate syntrophic interactions to be relevant in degradation. Together with feedback mechanisms between abiotic and biotic habitat components, this promotes an enhanced ecological perspective of the anaerobic degradation of BTEX hydrocarbons, as well as its incorporation into updated concepts for site monitoring and bioremediation.

Aerobic Toluene Degraders in the Rhizosphere of a Constructed Wetland Model Show Diurnal Polyhydroxyalkanoate Metabolism

Constructed wetlands (CWs) are successfully applied for the treatment of waters contaminated with aromatic compounds. In these systems, plants provide oxygen and root exudates to the rhizosphere and thereby stimulate microbial degradation processes. Root exudation of oxygen and organic compounds depends on photosynthetic activity and thus may show day-night fluctuations. While diurnal changes in CW effluent composition have been observed, information on respective fluctuations of bacterial activity are scarce. We investigated microbial processes in a CW model system treating toluene-contaminated water which showed diurnal oscillations of oxygen concentrations using metaproteomics. Quantitative real-time PCR was applied to assess diurnal expression patterns of genes involved in aerobic and anaerobic toluene degradation. We observed stable aerobic toluene turnover by Burkholderiales during the day and night. Polyhydroxyalkanoate synthesis was upregulated in these bacteria during the day, suggesting that they additionally feed on organic root exudates while reutilizing the stored carbon compounds during the night via the glyoxylate cycle. Although mRNA copies encoding the anaerobic enzyme benzylsuccinate synthase (bssA) were relatively abundant and increased slightly at night, the corresponding protein could not be detected in the CW model system. Our study provides insights into diurnal patterns of microbial processes occurring in the rhizosphere of an aquatic ecosystem.

IMPORTANCE

Constructed wetlands are a well-established and cost-efficient option for the bioremediation of contaminated waters. While it is commonly accepted knowledge that the function of CWs is determined by the interplay of plants and microorganisms, the detailed molecular processes are considered a black box. Here, we used a well-characterized CW model system treating toluene-contaminated water to investigate the microbial processes influenced by diurnal plant root exudation. Our results indicated stable aerobic toluene degradation by members of the Burkholderiales during the day and night. Polyhydroxyalkanoate synthesis in these bacteria was higher during the day, suggesting that they additionally fed on organic root exudates and reutilized the stored carbon compounds during the night. Our study illuminates microbial processes occurring in the rhizosphere of an aquatic ecosystem.

Modeling of the Reaction Mechanism of Enzymatic Radical C-C Coupling by Benzylsuccinate Synthase

Molecular modeling techniques and density functional theory calculations were performed to study the mechanism of enzymatic radical C–C coupling catalyzed by benzylsuccinate synthase (BSS). BSS has been identified as a glycyl radical enzyme that catalyzes the enantiospecific fumarate addition to toluene initiating its anaerobic metabolism in the denitrifying bacterium Thauera aromatica, and this reaction represents the general mechanism of toluene degradation in all known anaerobic degraders. In this work docking calculations, classical molecular dynamics (MD) simulations, and DFT+D2 cluster modeling was employed to address the following questions: (i) What mechanistic details of the BSS reaction yield the most probable molecular model? (ii) What is the molecular basis of enantiospecificity of BSS? (iii) Is the proposed mechanism consistent with experimental observations, such as an inversion of the stereochemistry of the benzylic protons, syn addition of toluene to fumarate, exclusive production of (R)-benzylsuccinate as a product and a kinetic isotope effect (KIE) ranging between 2 and 4? The quantum mechanics (QM) modeling confirms that the previously proposed hypothetical mechanism is the most probable among several variants considered, although C–H activation and not C–C coupling turns out to be the rate limiting step. The enantiospecificity of the enzyme seems to be enforced by a thermodynamic preference for binding of fumarate in the pro(R) orientation and reverse preference of benzyl radical attack on fumarate in pro(S) pathway which results with prohibitively high energy barrier of the radical quenching. Finally, the proposed mechanism agrees with most of the experimental observations, although the calculated intrinsic KIE from the model (6.5) is still higher than the experimentally observed values (4.0) which suggests that both C–H activation and radical quenching may jointly be involved in the kinetic control of the reaction.

Functional Gene Markers for Fumarate-Adding and Dearomatizing Key Enzymes in Anaerobic Aromatic Hydrocarbon Degradation in Terrestrial Environments

Anaerobic degradation is a key process in many environments either naturally or anthropogenically exposed to petroleum hydrocarbons. Considerable advances into the biochemistry and physiology of selected anaerobic degraders have been achieved over the last decades, especially for the degradation of aromatic hydrocarbons. However, researchers have only recently begun to explore the ecology of complex anaerobic hydrocarbon degrader communities directly in their natural habitats, as well as in complex laboratory systems using tools of molecular biology. These approaches have mainly been facilitated by the establishment of a suite of targeted marker gene assays, allowing for rapid and directed insights into the diversity as well as the identity of intrinsic degrader populations and degradation potentials established at hydrocarbon-impacted sites. These are based on genes encoding either peripheral or central key enzymes in aromatic compound breakdown, such as fumarate-adding benzylsuccinate synthases or dearomatizing aryl-coenzyme A reductases, or on aromatic ring-cleaving hydrolases. Here, we review recent advances in this field, explain the different detection methodologies applied, and discuss how the detection of site-specific catabolic gene markers has improved the understanding of processes at contaminated sites. Functional marker gene-based strategies may be vital for the development of a more elaborate population-based assessment and prediction of aromatic degradation potentials in hydrocarbon-impacted environments.

Identification of biomarker genes to predict biodegradation of 1,4-dioxane

Bacterial multicomponent monooxygenase gene targets in Pseudonocardia dioxanivorans CB1190 were evaluated for their use as biomarkers to identify the potential for 1,4-dioxane biodegradation in pure cultures and environmental samples. Our studies using laboratory pure cultures and industrial activated sludge samples suggest that the presence of genes associated with dioxane monooxygenase, propane monooxygenase, alcohol dehydrogenase, and aldehyde dehydrogenase are promising indicators of 1,4-dioxane biotransformation; however, gene abundance was insufficient to predict actual biodegradation. A time course gene expression analysis of dioxane and propane monooxygenases in Pseudonocardia dioxanivorans CB1190 and mixed communities in wastewater samples revealed important associations with the rates of 1,4-dioxane removal. In addition, transcripts of alcohol dehydrogenase and aldehyde dehydrogenase genes were upregulated during biodegradation, although only the aldehyde dehydrogenase was significantly correlated with 1,4-dioxane concentrations. Expression of the propane monooxygenase demonstrated a time-dependent relationship with 1,4-dioxane biodegradation in P. dioxanivorans CB1190, with increased expression occurring after over 50% of the 1,4-dioxane had been removed. While the fraction of P. dioxanivorans CB1190-like bacteria among the total bacterial population significantly increased with decrease in 1,4-dioxane concentrations in wastewater treatment samples undergoing active biodegradation, the abundance and expression of monooxygenase-based biomarkers were better predictors of 1,4-dioxane degradation than taxonomic 16S rRNA genes. This study illustrates that specific bacterial monooxygenase and dehydrogenase gene targets together can serve as effective biomarkers for 1,4-dioxane biodegradation in the environment.