ABSTRACT
Stable isotope fractionation was studied during the degradation of m-xylene, o-xylene, m-cresol, and p-cresol with two pure cultures of sulfate-reducing bacteria. Degradation of all four compounds is initiated by a fumarate addition reaction by a glycyl radical enzyme, analogous to the well-studied benzylsuccinate synthase reaction in toluene degradation. The extent of stable carbon isotope fractionation caused by these radical-type reactions was between enrichment factors (ε) of −1.5 and −3.9
, which is in the same order of magnitude as data provided before for anaerobic toluene degradation. Based on our results, an analysis of isotope fractionation should be applicable for the evaluation of in situ bioremediation of all contaminants degraded by glycyl radical enzyme mechanisms that are smaller than 14 carbon atoms. In order to compare carbon isotope fractionations upon the degradation of various substrates whose numbers of carbon atoms differ, intrinsic ε (εintrinsic) were calculated. A comparison of εintrinsic at the single carbon atoms of the molecule where the benzylsuccinate synthase reaction took place with compound-specific ε elucidated that both varied on average to the same extent. Despite variations during the degradation of different substrates, the range of ε found for glycyl radical reactions was reasonably narrow to propose that rough estimates of biodegradation in situ might be given by using an average ε if no fractionation factor is available for single compounds.
Many biochemical reactions are known to cause the fractionation of stable isotopes. Molecules consisting of lighter isotopes are utilized preferentially, and consequently heavier molecules are enriched in the residual substrate pool. Well-known examples are autotrophic CO2 fixation by plants (22) and bacterial methanogenesis (14). Within the last few years, several studies investigated stable isotope fractionation during bacterial degradation of contaminants such as toluene (17), benzene (15), and chlorinated hydrocarbons (12) that were associated with a significant enrichment of 13C in the residual substrate fraction. These findings opened the door to assessing contaminant degradation at polluted sites qualitatively, and under certain conditions even quantitatively, and of outlining the future development of the contamination and its potential impact on the environment and the drinking water supply. This concept has been applied successfully in several field experiments so far (12, 16, 24, 25).
Basic features of microbial stable isotope fractionation have been examined for anaerobic and aerobic bacterial cultures with toluene as a model substrate (18). It was demonstrated that the first enzyme reaction of toluene degradation is the rate-limiting step and that this reaction is also the key process leading to isotope fractionation. Isotope effects due to the transport of toluene to and into the cells were negligible. Furthermore, it was shown that various aerobic and anaerobic degradation reactions led to characteristic degrees of fractionation (19). However, it is not possible to deduce the type of the underlying reaction mechanism from the extent of fractionation, because fractionation factors (α) can vary significantly between identical reactions. Studies dealing with the stable isotope fractionation of particular reaction mechanisms have provided information about the span of isotope fractionation effects caused by these mechanisms. Based on theoretical considerations and calculations, it was proposed that the rate-limiting step in the benzylsuccinate synthase reaction is the addition of fumarate to the benzyl radical (10).
Stable isotope fractionation during degradation reactions initiated by glycyl radical enzymes is of particular interest because this reaction mechanism seems to be predominant in the anaerobic degradation of contaminants such as aromatic hydrocarbons, alkyl phenols, and alkanes. Regardless of the electron acceptor employed, anaerobic bacterial toluene degradation proceeds via benzylsuccinate formation in all cases investigated so far (8, 26). The enzyme mechanism was first described for anaerobic toluene degradation by denitrifying bacteria, where the enzyme benzylsuccinate synthase catalyzes the addition of a fumarate molecule to the methyl group to form benzylsuccinate (3). Fumarate addition and subsequent activation and β oxidation convert the former methyl group to a carbonylic function that acts as an entry port for single electrons in the subsequent ring reduction (4). The glycyl radical formed in benzylsuccinate synthase and related enzymes derives from 5′-deoxyadenosyl that is a cleavage product of S-adenosylmethionine (7). The enzyme reaction withdraws one hydrogen atom from a highly conserved cysteine residue of the enzyme. This thiyl radical then takes one hydrogen atom from the methyl group of the aromatic substrate, which later is transferred from the cysteine residue to the Cβ of the succinyl side chain of the produced benzylsuccinate (7, 13) (Fig. 1).
Reaction mechanism proposed for benzylsuccinate synthase modified according to the work of Frey (7) for m-xylene, o-xylene, m-cresol, and p-cresol, with R1, R2, and R3 equal to H, CH3, and OH, respectively. E, polypeptide chain of the enzyme.
We investigated stable isotope fractionation in degradation pathways initiated by a glycyl radical mechanism with two bacterial strains. Desulfobacterium cetonicum, a toluene-degrading sulfate reducer, was shown to degrade m-cresol and p-cresol via fumarate addition (20, 21). Enzyme activities in cell extracts showed that the cresols are probably not converted by the toluene-degrading benzylsuccinate synthase, but whether the two different cresol isomers are attacked by the same enzyme remains undetermined.
Degradation of m-xylene and o-xylene by the sulfate-reducing strain OX39 also proceeds via fumarate addition, which was confirmed by the identification of methylbenzylsuccinate derivatives in culture supernatants (B. Morasch, unpublished data). Induction experiments with strain OX39 and with m-xylene, o-xylene, and toluene showed that every substrate needed a specific enzyme for degradation (B. Morasch, unpublished data).
The objective of the present study was to systematically investigate stable carbon isotope fractionation during radical enzyme reactions of aromatic compounds which are prominent groundwater contaminants. The study should elucidate whether this type of reaction produces consistent isotope fractionation during degradation. For better comparison of isotope fractionations of the various substrate molecules, data for the intrinsic isotope fractionation at the molecular site of the reactions are provided.
MATERIALS AND METHODS
Cultivation of bacteria.The sulfate-reducing strain OX39 was isolated from soil from a site contaminated with BTEX (benzene, toluene, ethylbenzene, and xylene) and polycyclic aromatic hydrocarbons near Stuttgart, Germany (B. Morasch, unpublished data). D. cetonicum (DSM 7267) was taken from the lab culture collection.
Bacteria were cultivated at 30°C in bicarbonate-buffered freshwater (strain OX39) or brackish (D. cetonicum) mineral medium, pH 7.4, with sulfate (10 mM) as the electron acceptor. The medium was prepared under an atmosphere of N2-CO2 (80:20) and reduced with Na2S (1 mM) (27). A sterile, anoxic FeCl2 solution was added to the medium of strain OX39 to a final concentration of 3 mM.
Strains were grown in 120-ml serum bottles half filled with mineral medium and tightly sealed with Viton rubber stoppers (Maag Technic, Dübendorf, Switzerland). Cultures for isotope fractionation experiments were inoculated with 6 ml of precultures. m-Xylene and o-xylene were injected directly into the culture bottles through the rubber stoppers with microsyringes. m-Cresol and p-cresol were added from aqueous stock solutions (100 mM). Metabolic activity was monitored by observing sulfide production (6).
Hydrocarbon analysis.Xylene concentrations were determined by high-performance liquid chromatography (Bischoff Chromatography, Leonberg, Germany) with a C18 reversed-phase column (Prontosil, 200 by 3 mm, 3-μm film thickness; Bischoff) at 30°C and by UV detection at 210 nm with a mix of acetonitrile (Chromasol V super gradient grade; Fluka, Buchs, Switzerland) and demineralized water (70:30 [vol/vol]) as eluents. Cresol concentrations were determined by using a mix of acetonitrile and ammonium phosphate buffer (100 mM, pH 3.5) (50:50 [vol/vol]) as the eluent. Culture samples were diluted 1:5 with ethanol (99.9% gradient grade) and centrifuged (20,000 × g, 5 min) to remove precipitates before analysis.
Isotope analysis. 13C/12C stable carbon isotope ratios were determined by isotope ratio-monitoring gas chromatography-mass spectrometry. The system consisted of a gas chromatograph (HP-5890; Hewlett-Packard Co., San Diego, Calif.) connected via a combustion unit (gas chromatograph-combustion interface; Finnigan, Bremen, Germany) with an isotope mass spectrometer (MAT 252; Finnigan). Samples were measured as described previously (18).
Calculations.The 13C/12C isotope ratios of the substrate were calculated as relative δ13C values (per mille) according to equation 1 below, where Rsample is the 13C/12C isotope ratio of the sample and Rstd is the isotope ratio of the international Pee Dee Belemnite standard (11). Kinetic isotope α (αC) were calculated by using equation 2, which is derived from the Rayleigh equation for closed systems (11, 23). Ct/C0 is the fraction of the substrate remaining in the sample at time t. If ln(Rt/R0) is plotted over ln(Ct/C0), for the time intervals (t), the slope of the linear regression curve gives the αC as (αC − 1). Enrichment factors (ε) are a convenient expression of stable isotope fractionation that can be retrieved directly from αC using equation 3 (5).
\mathrm|<||<|\delta|>|^|<|13|>||>|C\mathrm|<| []|<|=|>||>|\left(\frac|<|R\mathrm|<|_|<|sample|>||>||>||<|R\mathrm|<|_|<|std|>||>||>|\mathrm|<||<|-|>|1|>|\right)\mathrm|<||<|\times|>|1,000|>|(1)$$mathtex$$\[\mathrm{ln}\left(\frac{R_{t}}{R_{\mathrm{0}}}\right)\mathrm{{=}ln}\left(\frac{C_{t}}{C_{\mathrm{0}}}\right)\mathrm{{\times}({\alpha}}C\mathrm{{-}1)}\]$$mathtex$$(2)$$mathtex$$\[\mathrm{{\varepsilon}{=}1,000{\times}({\alpha}}C\mathrm{{-}1)}\]$$mathtex$$(3)
Derivation of αintrinsic.The general definition of the isotope fractionation factor α is $$mathtex$$\[\mathrm{{\alpha}{=}}k_{H}/k_{L}\]$$mathtex$$(4) where kH and kL are the reaction rate constants of the heavier and lighter isotopomeres, respectively. The intrinsic isotope fractionation factor α (αintrinsic) is defined as $$mathtex$$\[\mathrm{{\alpha}_{intrinsic}{=}}k_{H\mathrm{,intrinsic}}\mathrm{/}k_{L}\]$$mathtex$$(5) with kH,intrinsic being the reaction rate constant at the position of the molecule where the chemical reaction takes place. For our model compound toluene, the heavier atom may be located at the methyl group (C-7) or at any position of the aromatic ring. Therefore, the overall reaction rate constant (kH) is the sum of the individual reaction rate constants (kH,i) for the location of the heavy atom at any position (i) of the labeled compound divided by the number of atoms (n) of the molecule, as shown in general equation 6. $$mathtex$$\[k_{H}\mathrm{{=}(}k_{H\mathrm{,1}}\mathrm{{+}}k_{H\mathrm{,2}}\mathrm{{+}{\ldots}{+}}k_{H\mathrm{,}n}\mathrm{)/}n\mathrm{{=}}{\mathrm{{\sum}}}k_{H\mathrm{,}i}\mathrm{/}n\]$$mathtex$$(6) In the specific case of the carbon stable isotope fractionation of toluene [kH = (kH,1 + kH,2 +… + kH,7)/7], the chemical reaction (fumarate addition) takes place at the methyl group of toluene (the C-7 atom). Therefore, only the reaction rate kH,7 (when the 13C atom is located at the methyl group) is subject to a kinetic isotope effect (kH,intrinsic). For all other reaction rate constants, kH,i, the 13C atom is located elsewhere in the molecule. With the assumption that secondary isotope effects are negligible compared to primary isotope effects, the kH,i rate constants for 12C atoms at the methyl group and 13C atoms at nonreactive positions of the molecule are assumed to be identical to the reaction rate constant kL for nonlabeled molecules. In addition, we assume that the probability of there being two 13C carbon atoms per toluene molecule is statistically unlikely due to the low natural abundance of 13C. With the exception of kH,intrinsic, the kH,i reaction rate constants are therefore equal to the reaction rate constant kL for the 12C isotopomeres. Equation 6 can therefore be written as $$mathtex$$\[k_{H}\mathrm{{=}}\frac{\mathrm{1}}{n}\mathrm{{\times}}\left[k_{H\mathrm{,intrinsic}}\mathrm{{+}}k_{L}\mathrm{{\times}(}n\mathrm{{-}1)}\right]\]$$mathtex$$(7) Equation 7 inserted into equation 4 gives $$mathtex$$\[\mathrm{{\alpha}{=}}\frac{\mathrm{1}}{n}\mathrm{{\times}}\left[k_{H\mathrm{,intrinsic}}\mathrm{{+}}k_{L}\mathrm{{\times}(}n\mathrm{{-}1)}\right]\mathrm{/}k_{L}\mathrm{{=}}\frac{\mathrm{1}}{n}\mathrm{{\times}}\left[\frac{k_{H\mathrm{,intrinsic}}}{k_{L}}\mathrm{{+}(}n\mathrm{{-}1)}\right]\]$$mathtex$$(8) Equation 5 inserted into equation 8 results in $$mathtex$$\[\mathrm{{\alpha}{=}}\frac{\mathrm{1}}{n}\mathrm{{\times}}\left[\mathrm{{\alpha}_{intrinsic}{+}(}n\mathrm{{-}1)}\right]\]$$mathtex$$(9) The exponent of the Rayleigh equation 2 can therefore be written as $$mathtex$$\[\mathrm{{\alpha}{-}1{=}}\frac{\mathrm{1}}{n}\mathrm{{\times}[{\alpha}_{intrinsic}{+}(}n\mathrm{{-}1)]{-}1{=}}\frac{\mathrm{1}}{n}\ \mathrm{({\alpha}_{intrinsic}{-}1)}\]$$mathtex$$(10) or, according to equation 3, as $$mathtex$$\[\mathrm{{\varepsilon}{=}{\varepsilon}_{intrinsic}/}n\]$$mathtex$$(11)
RESULTS AND DISCUSSION
Stable 13C/12C isotope fractionation was determined during the growth of D. cetonicum with m-cresol and p-cresol that had initial δ13C signatures of −29.32
± 0.36
and −27.50
± 0.36
, respectively. Anaerobic m-cresol degradation caused a strong increase in δ13C in the residual substrate fraction, which is shown in a representative experiment where degradation of 68% of the initial 1,030 μM led to an isotope shift from −29.04
± 0.47
to −24.53
± 0.41
(Fig. 2A). The average isotope ε of −3.9
± 0.5
for anaerobic m-cresol degradation was determined from four replicates (Fig. 3). Anaerobic degradation of p-cresol produced a shift in δ13C from −27.11
± 0.17
to −23.44
± 0.25
during the degradation of 89% of the initial 1,052 μM (data not shown). The smaller increase in δ13C per amount of substrate degraded resulted in an average ε of −1.6
± 0.1
calculated from three replicates (Fig. 3). For comparison, the isotope fractionation produced during toluene degradation by D. cetonicum resulted in an ε of −2.2
± 0.4
(18).
13C enrichment in the residual substrate fraction during the anaerobic degradation of m-cresol (▴) by D. cetonicum (A) and o-xylene (•) by the sulfate-reducing strain OX39 (B). Changes in sulfide concentration (Δ and ○) and δ13C (⧫) were monitored over time. The carbon isotope composition is presented as an average of six individual measurements, with error bars indicating the standard deviations. The diagram shows the results of one representative experiment out of four and of three replicates for m-cresol and o-xylene. d, days.
Stable carbon isotope fractionation during the degradation of m-xylene (—) and o-xylene (---) by strain OX39 and of m-cresol (·-·-) and p-cresol (···) by D. cetonicum. Regression lines of the 13C/12C isotope data over the respective concentrations are plotted according to equation 2, with r2 values of 0.968, 0.893, 0.942, and 0.946 for m-xylene, o-xylene, m-cresol, and p-cresol, respectively.
Stable carbon isotope fractionation was determined during the degradation of m-xylene and o-xylene by strain OX39 in three independent growth experiments per substrate. The initial δ13C carbon signatures of the substrates were −27.66
± 0.44
for m-xylene and −28.19
± 0.06
for o-xylene. During the degradation of 82% of the initial m-xylene (203 μM), the 13C isotope became enriched to a level of −24.81
± 0.16
in one of three growth experiments (data not shown). Three replicate m-xylene degradation experiments yielded an average ε of −1.8
± 0.2
(Fig. 3). Degradation of an initial o-xylene concentration of 257 μM by strain OX39 resulted in a δ13C shift from −28.28
± 0.13
to −24.72
± 0.19
when 85% of the substrate was used (Fig. 2B). The average ε of three replicates was −1.5
± 0.1
(Fig. 3) and thus similar to the factor found for the degradation of the m-xylene isomer. In previous degradation experiments using a sediment column filled with contaminated aquifer material under sulfate-reducing conditions, carbon isotope fractionation by an ε of −1.1
for o-xylene, which was lower than the ε reported here, was determined (24). Furthermore, carbon isotope fractionation resulting in an ε of −0.9
± 0.1
was found for the degradation of 2-methylnaphthalene by a sulfate-reducing enrichment culture (9) which also initiates degradation via fumarate addition (2).
By far the most information on isotope fractionation by glycyl radical reactions is available concerning toluene degradation under various redox conditions. In a comparative study, carbon isotope α determined for bacterial pure cultures using Fe(III), NO3−, or SO42− as electron acceptors were found to be between an ε of −1.8
and an ε of −1.7
(17). These results were in the same order as those we obtained for p-cresol, m-xylene, and o-xylene. Another study reported isotope fractionation upon toluene degradation yielding ε of −0.8
and −0.5
by sulfate-reducing and methanogenic enrichment cultures, respectively, which were significantly lower than those measured in our experiments (1).
To summarize the data available on isotope fractionation during degradation reactions employing glycyl radical mechanisms, the lowest fractionation (ε = −0.5
) was found for toluene degradation under methanogenic conditions, and the strongest fractionation (ε = −3.9
) was found for m-cresol degradation by D. cetonicum. Most ε obtained were about −1.8
. It should be emphasized that isotope fractionation in the radical reactions of the benzylsuccinate synthase-type enzymes involves just one methyl carbon atom at the reactive site. However, determination of ε by compound-specific isotope analysis gives overall δ13C values of the entire molecule. This discrepancy can be overcome by calculating the fractionation factor site-specific isotopes αintrinsic, which refers to the atom of the target molecule where the enzymatic reaction takes place (equation 9).
The use of this equation allows the normalization of isotope fractionation of aromatic compounds with different numbers of carbon atoms in order to study specific reactions in detail. However, the comparison of intrinsic enrichment factors (εintrinsic) upon reactions of the benzylsuccinate synthase type in Table 1 shows the same extent of variation of the average of the intrinsic enrichment factors (εintrinsic = −12.3 ± 5.8
) as the average of the overall enrichment factors (ε = −1.67 ± 0.86
). The standard deviations are 51 and 47%, respectively (Table 1). Obviously, the specific properties of every single enzyme catalyzing the same type of reaction but taking other substrates influence the degree of isotope fractionation. Thus, the use of the εintrinsic does not lead to additional information.
Comparison of compound-specific carbon isotope enrichment factor ε to the εintrinsic calculated for the carbon atoms at the molecular site of the reaction of the respective substratea
Based on the observed isotope enrichment factor ε, we estimated the maximal molecular mass of the substrate that would still allow us to measure isotope fractionation accurately in the field. As shown above, the εintrinsic of enzyme reactions of the benzylsuccinate type are more or less equal. However, the larger the molecular mass of the substrate, the less overall isotope fractionation can be measured, because the isotope effect is diluted with an increasing number of carbon atoms. Taking an average inaccuracy of isotope analysis of ±0.5
at the detection limits, the absolute isotope shift needed to reliably detect isotope fractionation would be about 2
if about 90% of the substrate was degraded. This value correlates to an overall ε of −0.895
. If we take the average intrinsic ε (εintrinsic) from Table 1, the presumptions would be fulfilled for molecules such as toluene with seven carbon atoms (εoverall = −1.76
) or xylene with eight carbon atoms (εoverall = −1.54
). The calculated borderline would be at molecules of 13 carbon atoms, with ε of −0.949
. Beyond 14 carbon atoms, the overall ε (−0.881
) is probably too low to allow significant isotope shifts upon degradation to be analyzed. 2-Methylnaphthalene with 11 carbon atoms (Table 1) is close to the analytical limit, but isotope fractionation can be measured under ideal conditions. However, the variations of the intrinsic isotope effects reported above emphasize the insecurity of exact measurements in the field, and the limitation for field applications will probably be at 11 to 12 carbon atoms.
In situ, the enrichment of heavier isotopes along a groundwater flow can be taken as qualitative evidence for biodegradation. In addition, the extent of bacterial degradation can be quantified by combining the isotope signatures determined in the field with the laboratory-derived α and using the Rayleigh equation (equation 2) (16, 24). The ε for m-cresol, p-cresol, m-xylene, and o-xylene from this study were used together with assumed changes in the isotope signatures between δ13C levels of −30
and −20
to calculate the portion of substrate remaining compared to the initial concentration (Fig. 4). The percentage of biodegradation is defined according to equation 12 (24):
$$mathtex$$\[B\ \mathrm{[\%]{=}100{\times}(1{-}}C_{t}\mathrm{/}C\mathrm{_{0})}\]$$mathtex$$(12)
To give an impression of how accurate an average α would be, we used an ε of −1.6
± 0.86
from Table 1 to calculate the extent of biodegradation based on assumed isotope shifts. At small isotope shifts of 1% in the field, the extent of biodegradation would be 45.5% ± 16.15% (Fig. 4). The error becomes smaller as the extent of biodegradation increases; e.g., an isotope shift of 5
results in an average biodegradation of 95.1% ± 6.87% (Fig. 4). Thus, it might be possible to apply an average ε to calculate the biodegradation of methylated compounds by glycyl radical enzymes in contaminated aquifers if the shifts in δ13C are sufficiently high.
Calculation of the percentage of bacterial biodegradation based on the assumed amount of δ13C (per mille) isotope shifts. Curves with symbols are based on the isotope enrichment factors determined here for D. cetonicum degrading m-cresol (▴) and p-cresol (▾) and for strain OX39 degrading m-xylene (○) and o-xylene (□). The error range of biodegradation calculated for the average glycyl radical type enrichment factor is marked in grey.
ACKNOWLEDGMENTS
We are grateful to Stefan Haderlein for continuous support and to Christian Griebler for assistance with the calculation of error propagations.
This work was financially supported by the Deutsche Forschungsgemeinschaft (grant Schi 180/7) and by the Bundesministerium für Bildung und Forschung (grant 02WT0022).
FOOTNOTES
- Received 4 August 2003.
- Accepted 22 January 2004.
- Copyright © 2004 American Society for Microbiology
