Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2002, p. 2484-2494, Vol. 68, No. 5
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.5.2484-2494.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Max-Planck-Institut für Terrestrische Mikrobiologie, D-35043 Marburg, Germany
Received 30 November 2001/ Accepted 26 February 2002
|
|
|---|
0.4 Pa), indicating the presence of a competitive population of hydrogenotrophic sulfate-reducing bacteria (SRB). This was paralleled by a suppressed activity of the hydrogenotrophic RC-I methanogens as indicated by the lowest SSU rRNA quantities detected in all experiments. Full inhibition of methanogenesis only became apparent when acetate was depleted to nonpermissive thresholds (<5 µM) after 10 days. Apparently, a competitive, acetotrophic population of SRB was not present initially, and hence, acetotrophic methanosarcinal populations were less suppressed than under ferrihydrite amendment. In conclusion, although methane production was inhibited effectively under both mitigation regimens, different methanogenic populations were either suppressed or stimulated, which demonstrates that functionally similar disturbances of an ecosystem may result in distinct responses of the populations involved. |
|
|---|
The mitigative potential of electron acceptor supplementation in the form of amorphous Fe(III) oxides (i.e., ferrihydrite [1, 5, 21]), sulfate-containing fertilizers (i.e., ammonium sulfate [9, 25, 50]), or other sulfate-containing amendments (i.e., gypsum or phosphogypsum [12, 26]) has been demonstrated in field and laboratory studies and reductions in CH4 emission rates of >70% have been reported (13, 21). The outcompetition of methanogens for common substrates like H2 or acetate by ferric iron-reducing bacteria (FRB) or sulfate-reducing bacteria (SRB) is well understood on a thermodynamic and process-based level. Conversely, little is known about the effects of alternative electron acceptor supplementation on the community structure of methanogens. How rapidly can inhibitory effects be observed, and are they reflected in the structure and dynamics of the methanogenic population or in rRNA expression levels? This knowledge is important for the better understanding of mitigation mechanisms and hence also for possible implementations for rice agriculture.
The aim of this study was to monitor the effects of different mitigation strategies on the methanogenic community in rice field soil using tools of molecular microbial ecology. Slurries of rice field soil were amended with ferrihydrite and gypsum and, together with an unamended control series, were incubated anoxically for 30 days. Small subunit (SSU) ribosomal DNA (rDNA) and rRNA dynamics of major methanogenic lineages, including the novel methanogens of rice cluster I (RC-I) (17, 31), were followed by terminal restriction fragment length polymorphism (T-RFLP) profiling. Furthermore, we determined SSU rRNA expression levels by PCR-independent quantitative rRNA hybridization. Our results revealed pronounced differences in methanogenic community structure and activity which were well correlated to biogeochemical measurements and dependent on soil treatment.
|
|
|---|
Extraction of nucleic acids.
DNA from soil samples was extracted as described previously (18) and purified using polyvinyl-polypyrrolidone spin columns (19). Total RNA was extracted using the protocol of Lüdemann et al. (30), with minor modifications. Six-hundred-microliter slurry aliquots were placed into 2-ml FastRNA tubes with lysing matrix B (Qbiogene, Heidelberg, Germany) containing 700 µl of precooled TPM buffer (50 mM Tris-HCl [pH 7.5], 1.7% [wt/vol] polyvinylpyrrolidone, 10 mM MgCl2) (15). The mixture was shaken for 40 s at maximum speed in a FastPrep 120 bead beater (Qbiogene) and cell debris and soil particles were sedimented by centrifugation (4 min; 20,000 x g; 4°C). The supernatant was saved, while the pellet was reextracted with 700 µl of a phenol-based lysis buffer (50 mM Tris-HCl [pH 7.5], 10 mM EDTA, 1% [wt/vol] sodium dodecyl sulfate [SDS], 6% [vol/vol] water-saturated phenol) by a second round of bead beating. After centrifugation, the combined supernatants were extracted with 500 µl each of water-saturated phenol, phenol-chloroform-isoamyl alcohol (25:24:1 [vol/vol/vol]), and finally chloroform-isoamyl alcohol (24:1 [vol/vol]) using 2-ml Phase Lock Gel Heavy cups (Eppendorf, Hamburg, Germany) for the last two steps to facilitate the clean separation of aqueous and organic phases and thereby maximize RNA yield. Total nucleic acids were precipitated from the aqueous phase with 3 volumes of ethanol and 0.1 volume of 3 M sodium acetate by freezing (-80°C overnight) and centrifugation (1 h; 20,000 x g; 4°C). The precipitate was washed once with 250 µl of cooled 70% ethanol, dried, and resuspended in 50 µl of EB buffer (10 mM Tris-HCl [pH 8.5]) (Qiagen, Hilden, Germany). RNA extracts were further purified by gel filtration with Sephadex G-75 spin columns (36) to remove residual humic acids. Extracted nucleic acids were checked by standard gel electrophoresis on a 1% agarose gel by using DNA SmartLadder (Eurogentec, Seraing, Belgium) or Escherichia coli 16S and 23S rRNA (Roche Diagnostics, Mannheim, Germany) as a reference and by ethidium bromide staining. Furthermore, nucleic acid extracts were quantified by standard UV photometry (BioPhotometer; Eppendorf).
PCR and RT-PCR amplification.
PCR amplification of archaeal SSU rRNA genes for T-RFLP was done with slurry DNA extracts as described before (32) using the primer combination of Ar109f (5'-ACK GCT CAG TAA CAC GT-3') (16) and 5' 6-carboxyfluorescein-labeled Ar912rt (5'-GTG CTC CCC CGC CAA TTC CTT TA-3') (MWG Biotech, Ebersberg, Germany). This modification of the original Ar915r primer (54) carries a 3-bp elongated 3' end and was designed to minimize nonspecific primer binding during reverse transcription (RT), especially as the selectivity of the original oligonucleotide has recently been shown to be questionable (53). The elongation of the primer has little effect on the absolute numbers of archaeal target sequences with zero to four mismatches in our ARB database (data not shown) but drastically reduces the numbers of nontarget bacterial matches with three and four mismatches (from n = 186 to 8 and 4,213 to 195, respectively). Consequently, it allows a more efficient discrimination between archaeal target and bacterial nontarget templates. This is especially important during RT, since RNA-DNA duplexes are more stable than DNA-DNA structures (48) and a low-stringency temperature of 48°C had to be used. The specificity of archaeal rRNA-targeted T-RFLP fingerprints was verified by comparison of rRNA- and rDNA-derived profiles generated with both the original Ar915r and the modified Ar912rt oligonucleotides. While only defined terminal restriction fragments (T-RFs) were detected for rDNA profiles with both primers (profiles were practically identical), the new primer eliminated putative nontarget amplicons (e.g., undefined T-RFs) below the detection threshold for rRNA profiles (data not shown).
For RT, coextracted DNA was removed from the slurry RNA extracts by DNase digestion of 1:10-diluted extracts. We added 5 µl of soil RNA to 5 µl of 10x DNase digestion buffer (200 mM Tris-HCl [pH 7.5], 100 mM MgCl2, 20 mM CaCl2), 10 U of RNase-free DNase I (Roche Diagnostics), and diethyl pyrocarbonate-treated H2O to a total volume of 50 µl. Incubation was performed for 1 h at 37°C, and heating to 80°C for 10 min stopped the reaction. Complete removal of DNA from the RNA extracts was verified by PCR (see above). RT-PCR was performed using the Access one-tube RT-PCR system (Promega, Hilden, Germany). The reaction mixture contained, in a total volume of 50 µl, 1x avian myeloblastosis virus-Tfl reaction buffer, 1 mM MgSO4, 200 µM deoxynucleoside triphosphates, 0.5 µM (each) primer (see above), 5 U of avian myeloblastosis virus reverse transcriptase, and 5 U of Tfl DNA polymerase. One microliter of DNase-treated, 1:10-diluted soil RNA extract was added as template. RT was carried out for 45 min at 48°C, followed immediately by 28 amplification cycles (30 s, 94°C; 45 s, 52°C; 90 s, 68°C) and a terminal extension step (5 min, 68°C). After amplification, PCR aliquots (5 µl) were visualized by standard agarose gel electrophoresis. Amplicons were purified with the MinElute PCR purification kit (Qiagen), re-eluted in 25 µl of EB buffer, and quantified by UV photometry.
T-RFLP analysis.
Since rDNA-targeted T-RFLP analysis has been shown earlier to be highly reproducible with replicate DNA extractions (33, 37), DNA was extracted only once from each time point and analyzed by triplicate T-RFLP profiling. T-RFLP analyses of rRNA were from triplicate RNA extractions for each time point. Restriction digests (TaqI; Promega) and T-RFLP analysis were performed as previously described (31). The digest (1.25 µl) was mixed with 0.95 µl of formamide loading dye (Amersham Pharmacia Biotech, Freiburg, Germany) and 0.3 µl of GeneScan-1000 (ROX) size standard (Applied Biosystems, Weiterstadt, Germany), denatured for 3 min at 95°C, and placed immediately on ice. T-RFLP analysis was performed on an ABI Prism 377 DNA sequencer (Applied Biosystems) in GeneScan mode. The fragment mixture (1.9 µl) was loaded into the wells of a 36-well shark tooth comb on a 36-cm 5% (wt/vol) polyacrylamide gel containing 6 M urea and 1x TBE buffer (89 mM Tris-borate, 2 mM EDTA).
Statistical analysis of T-RFLP data.
For statistics, T-RFLP data were evaluated as recently described (33, 44) using SYSTAT 10 software (SPSS Inc., Chicago, Ill.). Principal component analysis (PCA) was performed on T-RFLP profiles (14) with samples (soil treatment, time point, and targeted nucleic acids) as rows and the relative fluorescence intensities of all T-RFs as columns. A covariance data matrix was extracted with pairwise deletion and varimax factor rotation. Data reduction provided a two-factorial ordering of the variance of T-RFLP profiles, which was plotted on a map.
Quantitative rRNA hybridization.
Nonradioactive slot blot quantification of environmental RNA extracts was performed using the digoxigenin (DIG) system (Roche Diagnostics) as recommended by the manufacturer and described by Manz et al. (34), with minor modifications. In order to analyze all 18 different rRNA extracts (three treatments, six time points) together with an appropriate standard dilution series on one membrane for each probe, the RNAs originally extracted in triplicate from each sample were blotted in duplicate using a Hoefer PR 648 slot blot manifold (48 slots; Amersham Pharmacia Biotech). All hybridization, washing, and detection steps were done in 15-cm hybridization tubes (Ochs, Göttingen, Germany) in an OV4 hybridization oven (Whatman Biometra, Göttingen, Germany). After prehybridization, 50 pmol of DIG-labeled probe (MWG Biotech) was hybridized to the membrane-bound rRNA overnight at 40°C in 3 ml of hybridization buffer (5x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate], 5x Denhardt's solution, 1% blocking reagent [Roche], 0.1% N-lauroylsarcosine, 0.1% SDS). Washing steps were carried out for 30 min in 1x SSC-1% SDS (pH 7) twice at 40°C and once at the probe-specific stringent washing temperatures (45, 54, 57). DIG-labeled probes were then detected with Fab fragments of anti-DIG antibodies conjugated to alkaline phosphatase (Roche Diagnostics) and the ECF fluorescence substrate (Amersham Pharmacia Biotech) on a Storm 860 gel and blot imaging system (Molecular Dynamics, Sunnyvale, Calif.) (55). Blot images were quantitatively evaluated using ImageQuant 5.0 (Molecular Dynamics).
Oligonucleotide probes and reference RNAs.
The following oligonucleotide probes were used in this study: S-*-Univ-1392-a-A-15 (38), S-D-Arch-0915-a-A-20 (54), S-F-Mbac-0310-a-A-22, S-F-Msae-0825-a-A-23, S-G-Msar-0821-a-A-24 (45), and S-*-RCI-0645-a-A-23 (57). E. coli 16S and 23S rRNAs (Roche) were used as universal standards. Methanobacterium bryantii DSM 863T, Methanosaeta concilii DSM 3671T, and Methanosarcina barkeri DSM 800T were obtained from the Deutsche Sammlung von Mikroorganismen und Zellkulturen (DSMZ; Braunschweig, Germany) and grown in standard media as recommended by the DSMZ. RNA was extracted from harvested cell pellets following the above protocol for use as archaeal and family-specific standards. For the quantification of RC-I Archaea, 16S rRNA of a nearly full-length RC-I clone (AS04-16; accession number AJ308972) was in vitro transcribed (42) as previously described (57).
|
|
|---|
120 µM (Fig. 1D), while amendment of gypsum drastically elevated sulfate concentrations to up to 6.5 mM (Fig. 1G). High sulfate reduction rates were observed only after day 8 and were 0.5 µmol day-1 g-1 on average (Table 1). After day 26, the amended sulfate was depleted. Methane formation started directly after flooding (47) in all three experiments and increased to similar levels of
140 Pa on day 3 (Fig. 2). In the unamended soil, methane mixing ratios continued to increase and production rates were a constant 0.5 µmol day-1 g-1 on average after day 10. Under ferrihydrite amendment, methane production was strongly suppressed after day 3 (Fig. 1D and 2), but mixing ratios still increased with low constant production rates of 80 nmol day-1 g-1 thereafter. In the gypsum-amended soil, an inhibition of methane production became apparent only after day 6 but was absolute between days 11 and 22 (Fig. 1G and 2; Table 1). Only after the end of sulfate reduction (day 26) did methane mixing ratios start to increase again in low amounts. Total methane emissions during the experiment were reduced by 85% under ferrihydrite amendment and by 69% under gypsum amendment.
![]() View larger version (33K): [in a new window] |
FIG. 1. Time course of reduction processes in unamended (A, B, C), ferrihydrite-amended (D, E, F), and gypsum-amended (G, H, I) rice field soil slurries. (A, D, G) Concentrations of SO42- and Fe2+ and CH4 partial pressure. (B, E, H) H2 partial pressure and concentration of acetate. (C, F, I) Concentrations of propionate, isovalerate, and caproate (means ± standard deviations [SD]; n = 3).
|
|
View this table: [in a new window] |
TABLE 1. Rates of CH4 production, Fe3+ reduction, and SO42- reduction in unamended and ferrihydrite- and gypsum-amended rice field soil slurries within 30 days of anoxic incubation; inhibition time points of different processes; and contribution of different respiratory processes to total reducing equivalents formed
|
![]() View larger version (22K): [in a new window] |
FIG. 2. Time course of CH4 partial pressures (logarithmic scale) in unamended and electron acceptor-amended rice field soil slurries (means ± SD; n = 3).
|
92 Pa on day 1 in all parallel experiments (Fig. 1). It was depleted to
1 Pa on day 4 in all three experiments, but remained below 1 Pa (0.4 Pa on average after day 5) only under gypsum amendment (Fig. 1H). In the other two series, H2 partial pressures increased again to a second transient maximum of
10 Pa around day 9 and decreased to low levels of
0.5 Pa towards the end of incubation (Fig. 1B and E). Maximum acetate concentrations (
2 mM) also accumulated on day 1 after flooding (Fig. 1). Acetate consumption was most efficient in the ferrihydrite series, and concentrations rapidly decreased to a constant
60 µM from day 4 onward (Fig. 1E). In the other two series, depletion was much slower and reached concentrations below 100 µM only after 8 to 10 days. But while acetate slowly decreased from 60 to 25 µM after day 11 in the unamended slurries (Fig. 1B), it was depleted below the detection limit (
5 µM) between days 12 and 26 in the gypsum series (Fig. 1H). Other volatile fatty acids, like propionate, isovalerate, and caproate, transiently accumulated to substantial concentrations in the unamended experiment (Fig. 1C). Under ferrihydrite supplementation, depletion of all three compounds was accomplished on day 4 (Fig. 1F), while in the gypsum series, only the consumption of propionate was similarly efficient (Fig. 1I).
T-RFLP fingerprinting of archaeal SSU rDNA and rRNA dynamics.
Nucleic acids of all three flooding experiments were extracted on days 1, 5, 10, 16, 22, and 30 after flooding and analyzed by T-RFLP. Relative population dynamics within the archaeal communities of the three flooding experiments as revealed by combined rDNA- and rRNA-targeted T-RFLP analysis are shown in Fig. 3. Structure and compositional shifts of the natural archaeal community (Fig. 3A) were similar to those found recently (32, 43). Most prominent were the T-RFs characteristic for RC-I (392 bp) and Methanosarcina spp. (185 bp), both contributing
30% relative amplicon frequency (RAF) on average. While the methanosarcinal T-RF more than doubled from 17 to 39% between day 1 and day 10, the RC-I T-RF decreased from 37 to 25% concomitantly. Also abundant were T-RFs of 283 bp (Methanosaetaceae,
11%) and 91 bp (Methanobacteriaceae,
9%). The T-RFLP profiling of SSU rRNA revealed even more pronounced dynamics of Methanosarcina spp., since their T-RF increased to 47% RAF on day 10 (Fig. 3B). The methanosaetal T-RF was always more abundant in rRNA profiles and, in contrast to stable rDNA abundances, doubled its RAF from 13% on day 5 to 26% on day 30. The RC-I T-RF remained stable (
21%) but was always less abundant in rRNA profiles. rRNA-based RAFs of the Methanobacteriaceae were only initially comparable to rDNA ratios but decreased to
3% from day 10 onward.
![]() View larger version (64K): [in a new window] |
FIG. 3. Archaeal population dynamics in unamended (A, B), ferrihydrite-amended (C, D), and gypsum-amended (E, F) rice field soil slurries as determined by SSU rDNA (A, C, E) and SSU rRNA (B, D, F) targeted T-RFLP analysis. Shown is the percentage of individual T-RFs relative to the total integrated fluorescence (mean - SD; n = 3 T-RFLP replicates for rDNA and 3 extraction replicates for rRNA). Only T-RFs representing major methanogenic lineages are shown, and minor archaeal T-RFs described previously (6, 32, 44) have been combined as "Diverse." MB, Methanobacteriaceae; MS, Methanosarcinaceae; MX, Methanosaetaceae; RC-I, rice cluster I; RAF, relative amplicon frequency.
|
27% already after 22 days. Differences were also detected for the RC-I SSU rRNA RAF, which increased from 19 to 27% within the first 5 days after flooding and remained constant (
23%) thereafter. The addition of gypsum, on the other hand, resulted in archaeal rDNA dynamics that were quite similar to those of unamended soil. Again, a strong but less marked increase of Methanosarcina species amplicon ratios was observed with time, while all other T-RFs remained relatively stable (Fig. 3E). Differences to the unamended soil became evident only by rRNA profiling (Fig. 3F). Methanosarcinal frequencies reached a maximum after a strong initial increase already on day 5, decreased thereafter, but increased again slightly on day 30. The frequencies of RC-I rRNA, however, were only
17% on average and were always reduced compared to both unamended and ferrihydrite-amended rice field soil slurries.
Statistical analysis.
PCA is a preferred method for reduction of complex data sets and has been frequently used to facilitate the comparison of microbial communities with data derived from amplified gene restriction patterns (14, 41). PCA of T-RFLP profiles was done to visualize patterns of time-dependent community shifts within the three differently treated soils by generalizing variations in archaeal SSU rDNA and rRNA pools. To enhance visual clarity, PCA scores of averaged triplicate T-RFLP results were plotted. The total variance of T-RFLP fingerprints was reduced to two hypothetical principal components, which explained 62 and 28% of the variance, respectively. A factorial map (Fig. 4) shows the position, in ordination space, of archaeal SSU rDNA and rRNA fingerprints of the three soils and consecutive time points. PCA allowed for a clear separation of rDNA- and rRNA-derived profiles in ordination space. Furthermore, data points from day 1 after flooding were grouped for both rDNA and rRNA and thus indicate similar communities with almost identical T-RFLP patterns. Sequential rDNA dynamics within the natural and gypsum-amended slurries were also comparable, e.g., shifts in principal component ordination ran in parallel and end points were again closely grouped. However, PCA clearly showed that archaeal rDNA dynamics under ferrihydrite amendment were distinct and strongly reduced. For rRNA samples, successive shifts of principal component scores were more pronounced than for rDNA fingerprints (i.e., rRNA pools were more dynamic). Here, the unamended and gypsum-amended time series again displayed similar dynamics but were better discriminated than for rDNA profiles. Both were again clearly differentiated from the archaeal rRNA fingerprints of the ferrihydrite slurries.
![]() View larger version (24K): [in a new window] |
FIG. 4. Principal component ordination of the overall variance between averaged triplicate T-RFLP results generated with archaeal SSU rDNA and rRNA amplicons of natural and ferrihydrite- and gypsum-amended rice field soils. The abundance of all original T-RFs detected was included in this analysis. Successive time points (days) for each treatment and nucleic acid are connected by lines and indicated by numbers to illustrate time-dependent population shifts.
|
6.6 µg g-1 to
3.7 µg g-1 within the 30 days of incubation in all three flooding experiments (Fig. 5A). Total amounts of Archaea-specific rRNA (Fig. 5B) remained relatively constant (
1.45 µg g-1) in unamended rice field soil over time, while they decreased in both ferrihydrite- and gypsum-amended slurries to
1 µg g-1 (day 30). Relative to the total extractable rRNA (universal probe), the percentage of archaeal rRNA was
22% on day 1 in all slurries and 36, 30, and 27% on day 30 in the unamended and ferrihydrite- and gypsum-amended slurries, respectively.
![]() View larger version (26K): [in a new window] |
FIG. 5. Methanogenic population dynamics in unamended and electron acceptor-amended rice field soil slurries as determined by quantitative slot blot hybridization of environmental rRNA extracts (means ± SD; n = 2 extractions per time point). (A) Total rRNA; (B) archaeal rRNA; (C) Methanosarcina spp.; (D) Methanosaetaceae; (E) Methanobacteriaceae; (F) rice cluster I. Abbreviations of probes are from reference 3.
|
126 ng g-1 on average), with the exceptions of slightly elevated levels (168 ng g-1) on day 5 under ferrihydrite amendment (Fig. 5E). Finally, we followed the rRNA dynamics of the yet uncultured methanogens of RC-I (31), using an oligonucleotide probe designed to specifically target this novel lineage (57). rRNA amounts of RC-I in the unamended soil were very stable, although they increased slowly from 361 to 409 ng g-1 during the experiment (Fig. 5F). Dynamics were similar in the gypsum-amended slurries; however, absolute levels were constantly
100 ng g-1 below those in natural soil. Surprisingly, under ferrihydrite addition, RC-I rRNA levels increased drastically from 280 to 462 ng g-1 between day 1 and day 5 after flooding and decreased again to
330 ng g-1 towards the end of the experiment. The strongest transient rRNA dynamics of the novel RC-I methanogens were thus detected under Fe(III) supplementation. |
|
|---|
The population dynamics of methanogens were followed by PCR-T-RFLP analysis of archaeal SSU rDNA and rRNA as well as by quantitative rRNA hybridization probing. All major trends detected by rRNA-based T-RFLP analysis were also found by the PCR-independent quantitative hybridization results. By analyzing rDNA and rRNA in parallel, we were able to compare population dynamics (SSU rDNA frequency) and activity shifts of community members as indicated by changes in rRNA content. We found that rRNA expression patterns may reveal population activity (4, 35) much more sensitively than is possible by following the relative SSU rDNA frequency. For example, the increase in relative rRNA abundance of Methanosarcina spp. (Fig. 3 and 5) in unamended soil was much more pronounced than the increase detectable for SSU rDNA (Fig. 3), which suggests both activation and growth of this population. The caveat of measuring the relative gene abundance may be that a strong increase of one population (i.e., Methanosarcina spp.) could mask a small absolute increase of the gene copy number in another population (32). In the case of Methanosaeta spp., we were unable to detect an increase of rDNA ratio by PCR-T-RFLP (Fig. 3), but both RT-PCR-T-RFLP (Fig. 3) and hybridization probing (Fig. 5) revealed a doubling of rRNA amounts. The most parsimonious explanation appears to be that Methanosaeta spp. were strongly activated but did not grow significantly during the first 30 days after flooding. In fact, these findings may represent an important principle for the population dynamics of methanogens in rice field soil: only the fast growing and metabolically versatile Methanosarcina spp. are capable of growing under the transiently high substrate concentrations within the first days after flooding (32). In laboratory culture, exponentially growing bacteria exhibit a tight coupling between growth rate and ribosome synthesis, whereas DNA synthesis rates are regulated more gradually (35). In natural environments, under conditions of low nutrient availability or starvation, it is largely unknown how the synthesis of ribosomes (and rRNA) is actually linked to growth (i.e., cell division, DNA replication) of especially slow growing and uncultivated microorganisms. We hypothesize that Methanosaeta spp. in rice field soil initially increase their rRNA content only, while their relative SSU rDNA frequency, and thus most likely their population size, remains rather constant. Similarly, the activation of Methanosaeta spp. in a bioreactor was initially detected only by a significant increase in rRNA abundance (11).
Addition of 1.5% ferrihydrite extended the phase of iron reduction to over 30 days compared to
8 days in unamended soil. This was accompanied by a rapid but incomplete inhibition of methanogenesis after day 3 and by significant effects on the population dynamics of methanogens, indicating that a potent population of FRB successfully competed with methanogens and was indigenously present in the soil. However, the mechanism of inhibition of methanogenesis by ferrihydrite addition in rice field soil is not totally clear. Previous studies have shown that ferric iron reducers can outcompete methanogens for H2 in rice field soil (1) by lowering the H2 partial pressure to thresholds (<0.4 Pa) that are thermodynamically not permissive for methanogenesis (29). In principle, competition for acetate should be similar to competition for H2 (5), but the experimental data available are contradictory with respect to the underlying mechanism of inhibition. FRB were shown to either maintain a lower threshold for acetate in Texas rice field soil than methanogens (52) or have a similar threshold as found in Italian rice field soil (2).
In our experiments, neither the H2 partial pressures nor the acetate thresholds were lower in ferrihydrite-amended slurries than in unamended soil, albeit methanogenesis was strongly (but not completely) inhibited (Table 1). The second transient accumulation of H2 (
10 Pa on day 9), when iron was reduced at a high rate, indicates that H2 was a less competitive substrate than fatty acids for FRB. In contrast, acetate and other short- chain fatty acids were consumed more rapidly in ferrihydrite-amended slurries than in controls, but threshold concentrations for acetate were similar and did not decrease below 60 µM. This concentration is still permissive for acetotrophic methanogenesis by Methanosaeta spp. (22, 62). The lack of full inhibition of methanogenesis may be attributable to a decreasing bioavailability of Fe(III) as indicated by the progressive lowering of the iron reduction rate from day 7 onward (Fig. 3D). Only
200 µmol of Fe3+ g-1 of an available
240 µmol g-1 (
100 µmol of indigenous Fe3+ g-1 plus 140 µmol of added Fe3+ g-1) was reduced after 30 days. It cannot be ruled out that the synthetic ferrihydrite utilized for amendment was less bioavailable for FRB due to a higher degree of crystallinity than the indigenous ferrihydrite (51). The low acetate concentrations observed corroborate both the increase of methanosaetal rRNA (Fig. 3D and 5D), which was strongest under ferrihydrite amendment, and the suppressed dynamics of Methanosarcina spp. (Fig. 3D), as indicated by rDNA and rRNA analysis. Most likely, the acetotrophic Methanosaeta spp. contributed substantially to the low amounts of methane produced in spite of ferrihydrite-mediated inhibition, whereas Methanosarcina spp. were apparently outcompeted for both hydrogen and acetate by FRB after day 4. The minor initial increase of Methanosarcina species ratios was likely to be related to excess substrate availability until day 4 (Fig. 1).
Unexpectedly, RC-I methanogens were significantly stimulated temporarily (up to day 5) only under ferrihydrite amendment as indicated by rRNA analysis (Fig. 3D and 5F), as were members of the Methanobacteriaceae, albeit to a lesser extent (Fig. 5E). Enrichment cultures of the novel RC-I methanogens have been shown to grow on H2-CO2 (31), and the members of the Methanobacteriaceae are also hydrogenotrophic (62). The transient stimulation of these hydrogenotrophic methanogens was paralleled by an increase in the total archaeal rRNA content between days 1 and 5, which was strongest under ferrihydrite amendment. In the same samples, no surplus increase of rRNA of the acetotrophic Methanosarcina and Methanosaeta spp. was observed. Furthermore, in the ferrihydrite-amended slurries, initial sulfate concentrations were much lower than in the control, probably due to sulfate anion adsorption to ferric iron minerals (39). Hence, the activity of SRB, major competitors for H2 (8), was also reduced. Thus, under ferrihydrite amendment, it is likely that a much larger fraction of the initially produced H2 was available for hydrogenotrophic methanogens than in unamended soil. However, the strong transient increase, especially of RC-I rRNA, most likely cannot be explained considering the energy available from the small amount of methane actually produced within the first 5 days under ferrihydrite amendment. Rather, we speculate that under these conditions, energy conservation of hydrogenotrophic methanogens might have been coupled to Fe(III) reduction, a process which has been shown to occur in hyperthermophilic methanogens such as Methanopyrus kandleri and Methanococcus thermolithotrophicus (56). If this was true, based on the molecular fingerprinting of methanogenic populations, our results might indicate a previously unknown effect of ferric iron fertilization: besides being inhibited due to substrate competition for methanogenesis, methanogens might even be stimulated to reduce Fe(III) instead of CO2.
The addition of 0.15% gypsum extended the phase of sulfate reduction to over 26 days after flooding. However, the immediate metabolic potential of the intrinsic SRB population in Italian rice field soil was not sufficient to inhibit methanogens by substrate competition, as indicated by the delayed inhibition of methanogenesis from day 6 on only. SRB required more than twice as much time than FRB to induce similar inhibitory effects on the methanogens. In contrast to this longer lag phase, SRB were capable of completely inhibiting methanogenesis after day 10, most likely after a sufficient population size was established, until sulfate was finally depleted. Interestingly, the activation of SRB and the gradual inhibition of methanogens were paralleled by two phases of substrate depletion. H2 was depleted to mixing ratios as low as 0.3 Pa already after day 4 under gypsum amendment (Fig. 1H) and a second transient accumulation of up to 10 Pa as in the ferrihydrite-amended and unamended slurries was not observed. Thus, a competitive population of H2-utilizing SRB was indigenously present in Italian rice field soil and, upon addition of sulfate, capable of outcompeting methanogens for H2 directly, which has been observed earlier (1, 2). Conversely, acetate-utilizing SRB were capable of completely outcompeting acetotrophic methanogens only after
11 days of incubation, when acetate became depleted below detection (<5 µM). This was also the phase of constant sulfate reduction rates and complete inhibition of methanogenesis (Fig. 1G). Most likely, acetate-utilizing SRB were initially present at lower numbers than H2 utilizers or as spores, i.e., Desulfotomaculum spp. (59), which became active only after germination.
T-RFLP, PCA, and quantitative rRNA hybridization showed that methanogenic population dynamics under gypsum amendment were not as different from those of the unamended soil as from those under ferrihydrite addition. We relate this to the delayed inhibition of namely the acetotrophic methanogens, resulting in "intermediate" methanosarcinal rRNA dynamics and nearly unchanged rDNA dynamics. Interestingly, RC-I rRNA levels as revealed by T-RFLP and quantitative rRNA analysis were significantly lower under gypsum amendment, which is basically a reversion of the effect observed on day 5 under ferrihydrite addition and therefore underscores the potential of H2-utilizing SRB for outcompeting hydrogenotrophic methanogens.
In this study, we induced major functional shifts within the rice field soil microbial community by supplementing alternative electron acceptors in the form of ferrihydrite and gypsum and thereby promoting respiratory processes other than methanogenesis. The enhanced activity of FRB and SRB resulted in an almost complete inhibition of methanogenesis under conditions of limiting substrate and nonlimiting electron acceptor availability. Steady-state methane production rates in unamended slurries after day 10 (0.5 µmol day-1 g-1) corresponded to sulfate reduction rates in gypsum-amended slurries after day 8 (0.5 µmol day-1 g-1) and to iron reduction rates in ferrihydrite-amended slurries after day 7 (4.8 µmol day-1 g-1). Considering the electron uptake potential of eight electrons per CO2 and SO42- and one electron per Fe3+, only the amounts of sulfate reduced perfectly matched the quantity of methane, which was not produced under inhibition. The balance of ferrihydrite reduction showed that FRB oxidized electron donors other than acetate and H2, since more ferrihydrite was reduced than expected (Table 1); this may be another reason for the lower efficiency of inhibition of methanogenesis under ferrihydrite amendment. The total amount of gypsum added was only 1/10 that of the ferrihydrite amendment, but still the mitigative effects were comparable (69 and 85% methane reduction, respectively). More research is needed to fully evaluate the long-term effects of alternative electron acceptor supplementation and their mitigative potential for rice agriculture.
We thank Bianca Wagner for expert technical assistance and Ralf Conrad for continuing support.
|
|
|---|
-subunit (mcrA) genes in rice field soil and enrichment cultures reveal the methanogenic phenotype of a novel archaeal lineage. Environ. Microbiol. 3:194-204.[CrossRef][Medline]
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»