Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2000, p. 2732-2742, Vol. 66, No. 7
Max-Planck-Institut für terrestrische
Mikrobiologie, D-35043 Marburg, Germany
Received 1 February 2000/Accepted 11 April 2000
The population dynamics of Archaea after flooding of an
Italian rice field soil were studied over 17 days. Anoxically incubated rice field soil slurries exhibited a typical sequence of reduction processes characterized by reduction of nitrate, Fe3+, and
sulfate prior to the initiation of methane production. Archaeal population dynamics were followed using a dual approach involving molecular sequence retrieval and fingerprinting of small-subunit (SSU)
rRNA genes. We retrieved archaeal sequences from four clone libraries
(30 each) constructed for different time points (days 0, 1, 8, and 17)
after flooding of the soil. The clones could be assigned to known
methanogens (i.e., Methanosarcinaceae,
Methanosaetaceae, Methanomicrobiaceae, and
Methanobacteriaceae) and to novel euryarchaeotal (rice
clusters I, II, and III) and crenarchaeotal (rice clusters IV and VI)
lineages previously detected in anoxic rice field soil and on rice
roots (R. Grosskopf, S. Stubner, and W. Liesack, Appl. Environ.
Microbiol. 64:4983-4989, 1998). During the initiation of
methanogenesis (days 0 to 17), we detected significant changes in the
frequency of individual clones, especially of those affiliated with the
Methanosaetaceae and Methanobacteriaceae.
However, these findings could not be confirmed by terminal restriction
fragment length polymorphism (T-RFLP) analysis of SSU rDNA amplicons.
Most likely, the fluctuations in sequence composition of clone
libraries resulted from cloning bias. Clonal SSU rRNA gene sequences
were used to define operational taxonomic units (OTUs) for T-RFLP
analysis, which were distinguished by group-specific TaqI
restriction sites. Sequence analysis showed a high degree of
conservation of TaqI restriction sites within the different
archaeal lineages present in Italian rice field soil. Direct T-RFLP
analysis of archaeal populations in rice field soil slurries revealed
the presence of all archaeal lineages detected by cloning with a
predominance of terminal restriction fragments characteristic of rice
cluster I (389 bp), Methanosaetaceae (280 bp), and
Methanosarcinaceae/rice cluster VI (182 bp). In general,
the relative gene frequency of most detected OTUs remained rather
constant over time during the first 17 days after flooding of the soil.
Most minor OTUs (e.g., Methanomicrobiaceae and rice cluster
III) and Methanosaetaceae did not change in relative
frequency. Rice cluster I (37 to 30%) and to a lesser extent rice
cluster IV as well as Methanobacteriaceae decreased over
time. Only the relative abundance of Methanosarcinaceae (182 bp) increased, roughly doubling from 15 to 29% of total archaeal gene frequency within the first 11 days, which was positively correlated to the dynamics of acetate and formate concentrations. Our
results indicate that a functionally dynamic ecosystem, a rice field
soil after flooding, was linked to a relatively stable archaeal
community structure.
In most anaerobic ecosystems,
methanogenic Archaea significantly affect carbon cycling,
since fermentable substrates are degraded completely to CO2
and CH4 via the anaerobic food chain, with methanogenesis as the predominant terminal respiratory process (41).
Important habitats for methane-forming Archaea are wetlands,
including rice field soils. The latter are estimated to contribute
about 25% to the budget of global methane emissions (36)
and therefore have a significant impact on global warming
(25).
In periodically submerged wetlands such as rice field soils,
CH4 production is initiated shortly after flooding, only
after inorganic electron acceptors such as oxygen, nitrate,
manganese(IV), iron(III), and sulfate have been reduced sequentially.
Although biogeochemical processes after the flooding of rice field soil have been studied in detail, the factors controlling the initiation of
methanogenesis are not fully understood (10). The sequence of reduction processes is best described by the thermodynamic theory,
which predicts preferential reduction of available electron acceptors
with the most positive redox potential (35, 50). For a
variety of soils, this sequence of reduction processes has been
demonstrated after flooding (1, 32, 33). However, in most
soils, initiation of methane production deviates from the predicted
concept of sequential reduction. Methane formation, albeit at trace
amounts, could be observed directly after flooding of rice field soils
despite the presence of oxidants such as nitrate, iron(III), and
sulfate (40, 48, 49). Similarly, methane formation was
observed at redox potentials of 0 to 70 mV after flooding of oxic soils
(33) and in pure cultures of Methanosarcina barkeri (14), although it is generally assumed that
methanogens become active only at redox potentials below Understanding the function of ecosystems requires comprehension of both
biogeochemical processes and microbial community structures (10,
44). With respect to the archaeal populations involved in the
initiation of methanogenesis after the flooding of soils, the central
questions are which factors control the diversity and abundance of
methanogens and how, if at all, populations alter in response to
changing environmental conditions.
Few data are available about the population dynamics of methanogens
after the flooding of soils. Using most-probable-number counts, it has
been shown for geographically diverse rice field soils that the total
population size of methanogens remains rather constant even after
periods of drainage and aeration (30) and when monitored
over 2 years (3). However, the methanogenic activity may
undergo a shift from hydrogenotrophy to acetotrophy during the
initiation of methane production, as suggested by methyl fluoride
inhibition experiments (40). Currently, we can only speculate whether this shift in activity is also reflected in methanogenic population dynamics.
The known limitations of cultivation-dependent approaches
(26) have hampered the retrieval of numerically abundant and
functionally relevant microorganisms from complex environmental
communities (2). The slow-growing or difficult-to-cultivate
methanogenic Archaea especially may require extremely long
incubation times (>1 year) to obtain most-probable-number cell counts
(16). The application of molecular methods (i.e., the
analysis of small-subunit [SSU] rRNA genes as universal markers
[17, 47]) for community structure analysis facilitates
detection and identification of microorganisms that are difficult to
cultivate or may be regarded as unculturable from our current knowledge
of the growth requirements of these microorganisms. With these methods
at hand, the study of population structures becomes feasible for
methanogenic ecosystems.
Molecular surveys of SSU rDNA sequences in anoxic rice field soil
(9, 16) and on rice roots (15, 24) revealed a
high archaeal diversity, including known methanogenic families
(Methanosarcinaceae, Methanosaetaceae,
Methanomicrobiaceae, and
Methanobacteriaceae) as well as phylogenetically novel
eury- and crenarchaeotal lineages, termed rice clusters I to VI (RC-I
to RC-VI). The phylogenetic positions of only two of these clusters
(RC-I and -II) and their presence in methanogenic enrichment cultures
(24) suggest an affiliation with methanogenic taxa.
The ecophysiology of members of all other clusters (RC-III to RC-VI) is
still unknown.
In this investigation, we studied the dynamics of archaeal populations
after the flooding of rice field soil to establish a link between
structure (which species) and function (metabolic activities) of the
archaeal populations involved in the initiation of methanogenesis. Soil
slurries were incubated anoxically over 17 days and sampled for
chemical and molecular analyses. We followed changes in the archaeal
community structure by analyzing SSU ribosomal DNA (rDNA) clone
libraries derived from soil DNA extracts. In parallel, we determined
relative gene frequencies as analyzed by terminal restriction fragment
length polymorphism (T-RFLP) (27) analysis of archaeal SSU
rDNA PCR fragments (9). Our results indicate a surprisingly
stable archaeal community structure as analyzed by T-RFLP, considering
the changes associated with the flooding of the soil. Only
Methanosarcina-like populations were found to increase
significantly in relative gene frequency.
Soil samples and slurry experiments.
Soil was sampled in
February 1997 from a drained rice field of the Italian Rice Research
Institute near Vercelli (Po River valley, Italy). Soil parameters were
described previously (8). The soil was air dried and stored
at room temperature. Preparation of the soil and sieving (mesh size, 2 mm) were done as previously described (8). Soil slurries (3 g of dry soil and 3 g of sterile double-distilled water) were
prepared in triplicate for each time point in 12.5-ml serum vials which
were sealed with butyl rubber septa. Anoxic incubations of soil
slurries for each time point were started at daily intervals by the
addition of anoxic water with a syringe. The headspaces of vials were
flushed with N2 for 10 min, and vials were statically
incubated at 25°C. At the end of the experiment, the vials were
analyzed for gases (CH4, H2, and
CO2). Slurry samples were taken, frozen, and stored
( Chemical analyses.
Pore water samples were analyzed by ion
chromatography for the determination of nitrate and sulfate
(5). Acetate and formate were analyzed by high-pressure
liquid chromatography (23). Fe2+ was analyzed in
slurry samples using the ferrozin reaction method (1). Gas
samples (CH4, H2, and CO2) were
taken from the headspaces of the flasks and measured by gas
chromatography (40). Low H2 concentrations were
determined by using a reducing gas detector (RGD2; Trace Analytical,
Menlo Park, Calif.) (40).
Extraction of soil DNA and PCR amplification.
Total DNA was
extracted from one soil sample for each time point using a direct lysis
technique for cell disruption modified from Moré et al.
(31) as described previously (19). Briefly, cells
in soil samples (300 µl of slurry) were lysed by bead-beating (45 s
at 6.5 m s
0099-2240/00/$04.00+0
Copyright © 2000, American Society for Microbiology. All rights reserved.
Archaeal Population Dynamics during Sequential
Reduction Processes in Rice Field Soil
![]()
ABSTRACT
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
![]()
INTRODUCTION
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
150 mV
(50). Apparently, methanogenic Archaea in soils
are activated much faster than was previously thought (40).
![]()
MATERIALS AND METHODS
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
20°C) for analyses of biogeochemical parameters (Fe2+,
NO3
, NO2
,
SO42
, and fatty acids) and molecular analyses.
1) in the presence of sodium dodecyl
sulfate solution. DNA was purified from the supernatant with ammonium
acetate, isopropanol, and ethanol precipitation steps. Finally, the DNA
extract was further purified using a silica matrix-based purification
protocol (Prep-a-gene; Bio-Rad, Munich, Germany). Aliquots of DNA
extracts were analyzed by standard gel electrophoresis to verify
extraction. DNA concentrations were quantified fluorometrically
(PicoGreen dsDNA quantitation kit; Molecular Probes, Leiden, The Netherlands).
Cloning and sequencing. Archaeal SSU rDNA amplicons were cloned in Escherichia coli JM109 using the pGEM-T Easy Vector System (Promega, Mannheim, Germany) according to the manufacturer's instructions. Randomly selected clones were further analyzed and sequenced as described by Rotthauwe et al. (38).
Sequence data analysis and phylogenetic placement. Raw sequence data were assembled and checked with the Lasergene software package (DNASTAR, Madison, Wis.). The phylogenetic analysis (i.e., alignment and treeing) was performed by using the ARB software package (version 2.5b; O. Strunk and W. Ludwig, Technische Universität München, Munich, Germany, http://www.biol.chemie.tu-muenchen.de/pub/ARB/) as described previously in detail (15). Briefly, new sequences were fitted into an alignment of SSU rDNA sequences containing 758 partial (>700 nucleotides) or complete archaeal sequences retrieved from public databases (6) using the automated alignment tools of the ARB software package. When necessary, alignments were corrected manually. Chimera formation has been reported to occur during PCR amplification of mixed SSU rDNA (22, 46). The terminal 400 sequence positions of the 5' and 3' ends of the archaeal SSU rDNA sequences were used for fractional treeing to identify possible chimeras based on mismatches of the secondary SSU rRNA structure (28).
For treeing, almost full-length SSU rDNA sequences (>1,400 bases) were selected to construct an archaeal base frequency filter (50 to 100% similarity), which was subsequently used to generate an initial neighbor-joining tree using selected eury- and crenarchaeotal SSU rDNA reference sequences (>1,400 bases). Sequences from rice field soil clone libraries (635 to 770 nucleotides in length) were then added to this tree using the ARB parsimony tool, which allows the addition of short sequences to phylogenetic trees without changing global tree topologies (29). In addition, the tree topology was evaluated using fastDNAml as implemented in ARB. Public databases were screened for close relatives of rice field soil SSU rDNA clones as described previously (15). Rarefaction analysis of SSU rDNA clone libraries was performed using the Analytical Rarefaction software (version 1.2; S. M. Holland, University of Georgia, Athens, http://www.uga.edu/~strata/Software.html).T-RFLP analysis. T-RFLP analysis was performed as described by Chin et al. (9) with minor modifications. PCR assays (100 µl) were used for the amplification of archaeal SSU rDNA genes following the above PCR protocol. From 400 to 800 pg of template DNA was used. PCR products were labeled with 6-carboxyfluorescein (FAM) at the 5' end using the Ar912r primer. Alternatively, the 5' FAM-labeled Ar109f primer was used. Prior to digestion, PCR product concentrations were determined photometrically. DNA (50 ng), 3 U of TaqI (Promega), 1 µl of the appropriate 10× incubation buffer (Promega), and 1 µg of bovine serum albumin were combined in a total volume of 10 µl and digested for 2 h at 65°C. Fluorescently labeled terminal restriction fragments (T-RFs) were size separated on an ABI 373A automated sequencer (PE Applied Biosystems). T-RFLP electropherograms were analyzed by peak area integration of the different T-RFs (GeneScan 2.1 software; PE Applied Biosystems). The percent fluorescence intensity represented by single T-RFs was calculated relative to the total fluorescence intensity of all T-RFs to obtain a measure of relative SSU rRNA gene frequency.
Nucleotide sequence accession numbers. The sequences of SSU rDNA rice field clones of the four clone libraries (AS 00, AS 01, AS 08, and AS 17) were deposited in the GenBank, EMBL, and DDBJ databases under accession numbers AF225574 to AF225693.
| |
RESULTS |
|---|
|
|
|---|
Sequential reduction in rice field soil slurries.
Reduction
processes in the rice field soil started immediately after the onset of
anoxia. Nitrate was present at a low concentration (116 µM) initially
and fell below the detection limit within 24 h (Fig.
1A). Fe2+ increased rapidly
after the start of incubation and approached a maximum amount of 80 µM Fe2+ (per gram dry weight [dw] of soil) after 7 days. Sulfate in the pore water increased from 390 to 490 µM during
the first day, probably due to HCO3
or
Fe3+ reduction-mediated desorption from ferric iron
minerals (49), but was rapidly reduced subsequently and
could not be detected after 4 days. After only 1 day of incubation,
methane formation was detectable, albeit at low methane mixing ratios
(40). Within the first 4 days, methane increased
exponentially, as revealed by logarithmic plotting (Fig. 1B), although
electron acceptors other than CO2 were still present (Fig.
1A). After 6 days of incubation, methane was produced at a constant
rate of 0.5 µmol day
1 (g dw)
1.
|
Diversity of Archaea in rice field soil slurries.
We constructed four clone libraries (days 0, 1, 8, and 17) of archaeal
SSU rDNA fragments, which were PCR amplified from soil DNA, to assess
the archaeal diversity present, identify dominant populations, and
follow population dynamics during sequential reduction processes. More
than 30 clones of each library were randomly selected and analyzed by
comparative SSU rDNA sequencing. Within a total of 125 clones analyzed,
5 chimeric sequences were detected and excluded from further analyses.
Phylogenetic placement of the nonchimeric clones revealed that they all
fell within known eury- and crenarchaeotal lineages (Fig. 2 and
3), i.e.,
major methanogenic groups such as the Methanosarcinaceae,
Methanosaetaceae, Methanomicrobiaceae, and
Methanobacteriaceae, as well as the as yet uncultured
Archaea assigned to RC-I, -II, -III, -IV, and -VI (9,
15, 16). Most clone sequences were closely related to archaeal
rice field clones (>97% sequence similarity) reported earlier, but
some phylogenetically more distant sequences were also detected. For
example, the Methanosaeta-like clone AS 08-32 differed by
more than 5% in sequence from its closest relative, S15-24
(9). Clone AS 08-10, grouping with the
Methanomicrobiaceae, differed 6% from its closest relative
S15-30, and clone AS 01-23 (Methanobacteriaceae) exhibited
6.5% sequence difference from ARR21 (15), an SSU rDNA clone
from rice roots.
|
|
|
Phylotype richness.
The analysis of total community diversity
by cloning approaches may be biased by undersampling (43).
Therefore, we determined the diversity coverage of our archaeal SSU
rDNA clone libraries by rarefaction analysis (18). This
calculation allows the determination of the number of clones necessary
to cover the expected number of species
E(Sn) in a sample of n
individuals selected at random from a community containing N
individuals and S species. Since SSU rDNA clones in this
study originated from the same soil, we combined all clones sampled at
different time points with the rationale of estimating the total
archaeal diversity in the soil independent of community structure
shifts over time. The expected number of different SSU rDNA sequences
was plotted against the number of clones retrieved (Fig.
5). We defined sequences with >99%
sequence similarity (<8 different nucleotides per ~800 bp of SSU
rDNA fragment) as belonging to the same species. However, even SSU rDNA
sequences differing up to 3% may represent a single species
(42). Therefore, rarefaction calculations were performed at
both the 99% and 97% sequence similarity levels. Species rarefaction calculations for 120 clones reached saturation at neither 99% (56 different species) nor 97% (25 different species) sequence similarity levels, indicating undetected archaeal diversity.
The SSU rDNA sequences detected in our clone libraries fell into nine major clone groups (Fig. 4). By plotting the expected number of groups against the number of SSU rDNA clones retrieved, rarefaction calculations were shown to approach saturation at 120 clones, indicating sufficient diversity coverage at this level of resolution. According to these calculations, a random sample of 30 clones may
contain seven to eight of the different major lineages, and in fact we
detected six, seven, seven, and eight groups in the four clone
libraries.
|
T-RFLP analysis. It is well known that analysis of the microbial community structure by cloning may be biased (43, 45). Thus, the significant changes in archaeal community structure as detected by clone frequencies during sequential reduction processes may have been subject to cloning biases as well. Therefore, we analyzed PCR-amplified SSU rDNA fragments directly by T-RFLP to monitor archaeal community structure changes.
Analysis of community structures by T-RFLP results in a pattern of T-RFs that may be unrelated to the phylogeny of the SSU rDNA sequences studied, i.e., sequences with different group affiliations can be represented in the same operational taxonomic unit (OTU) (27). Therefore, we analyzed 267 archaeal clones obtained from Vercelli rice field samples in this study and in independent previous studies (9, 15, 16, 24) for group-specific TaqI restriction sites by sequence data comparison. The combined data set revealed a high degree of conservation of different group-specific restriction sites within the abundant and also within some of the less abundant clone groups (Table 1). For example, 95% of 77 Methanosarcina-like clones exhibited a 182-bp T-RF, 97% of 35 Methanosaeta-like clones showed a common T-RF of 280 bp, and 89% of 62 RC-I clones exhibited a 389-bp T-RF. Unfortunately, each of the three most frequent clone groups exhibited a common T-RF with one of the less frequent groups. Methanosarcina-like and RC-VI clones mainly shared a T-RF of 182 bp (9), Methanosaeta-like and RC-V clones exhibited a common 280-bp T-RF, and the majority of the RC-I and RC-II clones contributed to the 389-bp OTU. However, of the 120 clones analyzed in this study, only 1 (AS 00-12) fell within RC-II and none fell within RC-V. Thus, the 280-bp and 389-bp OTUs predominantly represented Methanosaeta-like and RC-I SSU rDNA sequences, respectively, and only to a negligible extent RC-II and RC-V sequences. In contrast, although Methanosarcina-like clones were clearly most abundant (40 of 120 clones), small numbers of the crenarchaeotal RC-VI clones were present in all clone libraries (n = 6), and therefore this population was likely to be represented in the 182-bp OTU. A typical archaeal community fingerprint obtained by T-RFLP analysis of SSU rDNA amplicons of slurry DNA is shown in Fig. 6A. All OTUs predicted by sequence data analysis (Table 1) could be detected in the T-RFLP electropherogram and assigned to the corresponding archaeal lineages.
|
|
Population dynamics of Archaea.
Direct T-RFLP analysis
was performed with PCR amplicons from soil DNA extracts for all time
points of our slurry experiment. By integrating the fluorescence
intensity of the different OTUs, relative archaeal SSU rRNA gene
frequencies were quantified (Fig. 7). The
major archaeal groups detected via T-RFLP analysis were identical to
the most abundant lineages found in clone libraries, e.g.,
Methanosarcinaceae, RC-I, and Methanosaetaceae.
But while the clone libraries displayed large shifts in relative
composition over time (see above, Fig. 4), gene frequencies as
determined by T-RFLP were rather constant, indicating a potential bias
inherent in the procedure of cloning. Most prominent was the 389-bp OTU assigned to RC-I, but its frequency remained relatively constant over
time, varying only between 30 and 35% of total archaeal gene frequency. The 182-bp OTU representing Methanosarcina-like
and RC-VI populations was the only OTU which increased significantly in
intensity during sequential reduction processes. A steady increase in
the 182-bp OTU from 15% on day 0 to 29% of total archaeal gene frequency on day 13 was accompanied by an apparent relative decrease in
especially the 88-bp, 389-bp, and 754-bp OTUs.
|
| |
DISCUSSION |
|---|
|
|
|---|
After flooding of rice field soils, methane production is typically initiated shortly after the reduction of electron acceptors such as nitrate, Fe3+, and sulfate. The sequence of these biogeochemical processes (metabolic activities) and the initiation of methane production observed in this study were in good agreement with previous reports (1, 39, 48, 49). Here, we focused on establishing a link between the observed changing biogeochemical processes (which ecosystem function) upon soil flooding to changes in archaeal populations, especially the methanogens (which species).
The different molecular approaches used resulted in pronounced differences in community composition determined during the initiation of methanogenesis. Analysis of SSU rDNA clone libraries (Fig. 4) suggested significant shifts in archaeal community compositions, whereas T-RFLP fingerprints exhibited a rather constant population structure when followed closely at daily intervals (Fig. 7). For example, after only 1 day of incubation, Methanosaeta-like and Methanobactericeae clones increased 10- and 5-fold in abundance, respectively. Subsequently, the same populations either decreased fourfold or disappeared completely after 7 days (Fig. 4).
Several lines of evidence suggest that most likely the cloning step was biased, resulting in large apparent fluctuations in the population structure. First, Methanosaeta spp. are characterized by slow growth (20). In fact, enumeration of acetotrophic methanogens in rice field soil required more than 40 weeks of incubation for the detection of these fastidious microorganisms (16). Although the in situ growth rates of these Archaea in soils are largely unknown, it is unlikely that the Methanosaeta-like populations increased 10-fold within 1 day in soil slurry experiments. Moreover, high concentrations of methanogenic substrates, such as acetate, present during the first 8 days of incubation (Fig. 1B) are known to be more favorable for development of the faster-growing Methanosarcina spp., which are adapted to higher acetate concentrations than are Methanosaeta spp. (20).
The strongest argument for a cloning bias, however, results from the comparison of the two methods employed in this study, cloning and T-RFLP analysis of SSU rDNA PCR products. We observed large shifts in community structure only with the cloning approach. Likewise, large fluctuations in archaeal SSU rDNA sequence abundances in clone libraries were observed in a recent study of archaeal communities before and after a temperature shift in methanogenic rice field soil (9). Again, T-RFLP analysis revealed less pronounced differences in community structure than were suggested by cloning analysis. On the other hand, in a study of a bioreactor sample, a combination of cloning and subsequent analysis of clones by ARDRA performed in three replications suggested that cloning appeared to be rather reproducible (12).
Several factors possibly influencing the abundance of sequences in clone libraries have been pinpointed, especially the PCR and the cloning procedure itself (for a review, see reference 45).
We cannot exclude an underlying PCR bias possibly affecting both methods, cloning and T-RFLP, because both methods rely virtually on the same PCR protocol (with the exception that primer Ar912r was fluorescently labeled for T-RFLP). Since the observed large fluctuations in community composition over time were restricted to cloning analysis, we were prompted to conclude that the cloning step in particular was subject to bias.
Possible bias related only to the cloning procedure may be attributed to random error resulting from undersampling of diversity (43), toxicity of vector inserts to the transformed host, or the choice of cloning kit. Rainey et al. (unfortunately, no original data were provided in this short communication [37]) reported differences in SSU rDNA clone library composition depending on the cloning system used, i.e., blunt-end versus sticky-end cloning. In our study we used TA cloning for separating SSU rDNA PCR products amplified from environmental samples. Although widely used (34), the influence of bias related to this cloning protocol has not yet been analyzed with respect to the sequence distribution in SSU rDNA clone libraries.
Despite the principal caveat that clone abundance in libraries does not necessarily reflect population abundance in environmental samples (17), our SSU rDNA clone libraries contained valuable information about the sequence diversity of Archaea in rice field soil at different time points (i.e., days 0, 1, 8, and 17) after flooding. We detected most of the archaeal groups (Methanosarcinaceae, Methanosaetaceae, Methanomicrobiaceae, Methanobacteriaceae, RC-I to -IV, and RC-VI) detected previously in studies of rice field soils from Vercelli focusing either on the diversity (15, 16, 24) or on community structure changes after a temperature shift (9). We failed to detect sequences related to the novel RC-V previously shown to be present in bulk soil of 90-day-old microcosms (15). However, it is still unknown whether RC-V as well as the other as yet uncultured Archaea assigned to RC-I to -VI (15, 16) are methanogenic. Their relevance during the initiation of methanogenesis in rice field soil remains to be elucidated. It is striking that although the soil was sampled in different years and certainly not at exactly the same site near the Vercelli Rice Research Institute, the majority of our SSU rDNA clones were closely related or even identical to previously reported archaeal sequences. Only a few clones were more distantly related, with up to 7% sequence differences from their closest previously reported relatives.
As indicated by rarefaction analysis, our clone libraries (n = 120) did not fully cover the phylotype richness of Archaea in Vercelli rice field soil (Fig. 5). However, rarefaction analysis focusing on the level of lineages detected revealed that it was unlikely to detect sequences representing additional or novel archaeal lineages (Fig. 5). Moreover, the high similarities of our SSU rDNA sequences to those from previous studies on Vercelli rice field soil also suggested that we would not detect novel phylotypes by analyzing additional clones.
By using a large data set comprising 267 clones from our study and those obtained with a similar PCR assay in the previous studies mentioned above, we were able to confirm previously detected restriction sites for the major phylogenetic lineages present in Vercelli rice field soil (Table 1) (9). Moreover, we could show a high degree of conservation of group-specific T-RFs for most of the lineages present (Table 1). For example, 95% of frequently found clones related to Methanosarcina spp. had a T-RF of 182 bp and 89% of all RC-I clones had a T-RF of 389 bp. On the other hand, clones related to Methanosarcina spp. shared the 182-bp T-RF with RC-VI clones. These RC-VI clones were also a frequent group and thus had to be considered for the analysis of relative gene frequency of this OTU.
We emphasize that the group specificity of restriction sites described here is valid only in the context of sequences retrieved from the same environment. Other sequences not detected in Vercelli rice field soil by PCR or cloning and belonging to other taxonomic groups may share common restriction sites with sequences retrieved from rice field soil samples (e.g., Methanospirillum sp. and RC-I) due to the conserved nature of the SSU rRNA gene. Knowing these limitations but equipped with the detected degree of conservation of TaqI restriction sites, we were able to expand the utility of the T-RFLP analysis considerably.
All predicted OTUs were actually found in direct T-RFLP analysis of slurry samples. For analyzing archaeal population dynamics, we extended the T-RFLP analysis from a qualitative assessment of OTUs to a quantitative assessment of relative gene frequencies by integrating peak areas of individual OTUs. T-RFLP analysis was highly reproducible, as indicated by the relatively constant gene frequencies of most OTUs measured at daily intervals over 17 days (Fig. 7). PCR replications (n = 3) were also highly reproducible, as indicated by a low average standard deviation (SD) (<1%) for all T-RFs. A similar approach was used successfully by Suzuki et al. (43) to determine relative gene frequencies of fluorescently labeled SSU rDNA PCR fragments of different lengths.
In general, the relative gene frequency of most detected OTUs (Fig. 7) remained rather constant over time during the first 17 days after flooding of the soil. Most minor populations (i.e., Methanomicrobiaceae and RC-III) and Methanosaetaceae apparently did not change at all in frequency. Notably RC-I (37 to 30%) and to a lesser extent RC-IV as well as Methanobacteriaceae gene frequencies decreased over time. The remarkable stability of the archaeal populations after the flooding of the soil is in contrast to the changes in environmental parameters such as redox potential (40), availability of electron acceptors (Fig. 1A), and availability of electron donors such as acetate and hydrogen (Fig. 1B). Denitrification intermediates were probably too low in concentration in our slurry experiments to affect methanogenesis by nitrogen oxide inhibition (21), as indicated by the low initial nitrate concentration (<116 µM).
The most striking point was the significant increase in relative gene frequency of the 182-bp OTU comprising Methanosarcinaceae and RC-VI populations from 15 to 29% over 13 days (Fig. 7). By using a novel T-RFLP assay (i.e., with the fluorescently labeled primer Ar109f) (Fig. 6B), which allowed the differentiation of these two groups, we could demonstrate that in fact only the relative gene frequencies of the Methanosarcinaceae increased. Theoretically, we have to consider that an increase in relative gene frequency of the Methanosarcinaceae could be linked to a decrease in the absolute gene frequencies of the other populations, caused, for instance, by predatory protozoa (11) and phages (4). However, it is more likely that growth of Methanosarcina-like populations occurred because of the energetically permissive concentrations of acetate and hydrogen during the initiation of methanogenesis (Fig. 1B). We detected a similar temporal pattern of community structure change in a study analyzing the effect of soil aggregate size on the initiation of methanogenesis (37a). By T-RFLP analysis, Chin et al. (9) also detected an increase in the OTU comprising Methanosarcinaceae and RC-VI upon anoxic incubation of methanogenic rice field soil and an even stronger increase in the RC-I-related OTU. The different incubation temperature (30°C) used in their study as well as other factors (i.e., soil sample, acetate accumulation to higher values, and lower methane production rate) may be responsible for differences in community structure changes from what was found in our study.
After the complete reduction of sulfate (day 4), acetate concentrations accumulated to a second maximum of about 600 µM around day 6, whereas hydrogen concentrations remained below 6 Pa. Apparently, the acetotrophic activity of methanogens was the limiting factor during the initiation of methanogenesis and led to the observed accumulation of acetate. This is in agreement with results from methyl fluoride inhibition experiments in anoxic rice field soil, which suggested that hydrogenotrophic methanogenesis dominates the initial phase of methane production (40). In addition, Chidthaisong and Conrad (7) also reported that acetate utilization by methanogens was slower than acetate production from glucose in rice field soil during the reduction phase (i.e., in the presence of sulfate and Fe3+). The second increase in acetate concentration was well correlated with the steady increase in the Methanosarcinaceae OTU (Fig. 7). From day 6 on, acetate concentrations decreased consistent with acetotrophic methanogenesis, e.g., by Methanosarcina spp.
Another important observation with respect to a link between metabolic activity and population dynamics was that the Methanosarcinaceae still increased (day 13; Fig. 7) when pore water concentrations of acetate (Fig. 1B) were well below the known threshold (>200 to 300 µM) for growth of Methanosarcina spp. Methanosaetaceae populations apparently could not benefit from their lower threshold for acetate and grow, as indicated by constant gene frequencies of the 280-bp T-RF. Most likely, the versatile Methanosarcina spp. switched to other electron donors, such as hydrogen or methanol (not measured). Interestingly, formate was present in a concentration of up to 120 µM between days 9 and 13 (Fig. 1B), and Methanosarcinaceae did not increase further after formate was depleted below the detection limit. Although formate cannot be utilized by the known Methanosarcina spp., formate is probably in equilibrium with H2/CO2 as catalyzed by fermentative bacteria.
Typically, only small amounts of methane are formed, probably from H2/CO2, during the initiation of methanogenesis (40, 48). The small amounts of methane (~0.56 µmol; Fig. 1B) that we detected during the first 3 days after flooding of the soil are contradicted by a significant increase in gene frequency (~4%) of Methanosarcina-like populations. The small amounts of methane formed were probably not sufficient to explain the observed increase in Methanosarcina-like populations by de novo growth. Further investigations are necessary to find out which methanogens actually become active during the important phase of initiation of methanogenesis. In a forthcoming study, we will therefore address this question by targeting directly the rRNA of Archaea. The rRNA content of cells is generally accepted as an indicator of the metabolic activity of microbial populations.
In summary, we observed that the population structure of Archaea remained remarkably constant over time after flooding of rice field soil. Only Methanosarcina-like populations increased significantly in gene frequency relative to the other populations. Probably due to a high organic matter content of the soil utilized, fermentation processes were vigorous. This was indicated by the initially high levels of acetate and hydrogen (day 1, Fig. 1). These conditions were obviously favorable for Methanosarcina-like populations. Even the switch from hydrogenotrophic to acetotrophic methanogenesis was apparently not accompanied by a population shift but instead by a shift in activity of the Methanosarcina-like populations, indicating that the same population could manage different ecosystem functions. We conclude that a functionally extremely dynamic ecosystem, a flooded rice field soil, was linked to a relatively stable archaeal community structure. This is in agreement with the temporal stability of the population structure observed in other environments such as cyanobacterial hot spring mats (13) as well as oxic (19) and anoxic rice field soils (9, 37a). On the other hand, ecosystems such as an apparently functionally stable methanogenic bioreactor fed with glucose may exhibit an extremely dynamic community structure (12), indicating that a general principle for microbial community dynamics relative to ecosystem function is not yet feasible.
| |
ACKNOWLEDGMENTS |
|---|
This study was supported by a grant from the Deutsche Forschungsgemeinschaft (DFG) within the SFB395 and the Max-Planck-Society.
We thank Bianca Wagner for excellent technical assistance, Peter F. Dunfield for English language suggestions, and Ralf Conrad for continuing support. We are indebted to Steven M. Holland for providing the Analytical Rarefaction 1.2 software, available at http://www.uga.edu/~strata/Software.html.
| |
FOOTNOTES |
|---|
* Corresponding author. Mailing address: Max-Planck-Institut für terrestrische Mikrobiologie, Karl-von-Frisch Strasse, D-35043 Marburg, Germany. Phone: 49-6421-178 830. Fax: 49-6421-178 809. E-mail: friedric{at}mailer.uni-marburg.de.
| |
REFERENCES |
|---|
|
|
|---|
| 1. | Achtnich, C., F. Bak, and R. Conrad. 1995. Competition for electron donors among nitrate reducers, ferric iron reducers, sulfate reducers, and methanogens in anoxic paddy soil. Biol. Fertil. Soils 19:65-72. |
| 2. |
Amann, R. I.,
W. Ludwig, and K. H. Schleifer.
1995.
Phylogenetic identification and in situ detection of individual microbial cells without cultivation.
Microbiol. Rev.
59:143-169 |
| 3. | Asakawa, S., and K. Hayano. 1995. Populations of methanogenic bacteria in paddy field soil under double cropping conditions (rice-wheat). Biol. Fertil. Soils 20:113-117. |
| 4. |
Ashelford, K. E.,
M. J. Day,
M. J. Bailey,
A. K. Lilley, and J. C. Fry.
1999.
In situ population dynamics of bacterial viruses in a terrestrial environment.
Appl. Environ. Microbiol.
65:169-174 |
| 5. | Bak, F., G. Scheff, and K. H. Jansen. 1991. A rapid and sensitive ion chromatographic technique for the determination of sulfate and sulfate reduction rates in fresh-water lake-sediments. FEMS Microbiol. Ecol. 85:23-30[CrossRef]. |
| 6. |
Benson, D. A.,
I. Karsch-Mizrachi,
D. J. Lipman,
J. Ostell,
B. A. Rapp, and D. L. Wheeler.
2000.
GenBank.
Nucleic Acids Res.
28:15-18 |
| 7. |
Chidthaisong, A.,
B. Rosenstock, and R. Conrad.
1999.
Measurement of monosaccharides and conversion of glucose to acetate in anoxic rice field soil.
Appl. Environ. Microbiol.
65:2350-2355 |
| 8. | Chin, K. J., and R. Conrad. 1995. Intermediary metabolism in methanogenic paddy soil and the influence of temperature. FEMS Microbiol. Ecol. 18:85-102[CrossRef]. |
| 9. |
Chin, K. J.,
T. Lukow, and R. Conrad.
1999.
Effect of temperature on structure and function of the methanogenic archaeal community in an anoxic rice field soil.
Appl. Environ. Microbiol.
65:2341-2349 |
| 10. |
Conrad, R.
1996.
Soil microorganisms as controllers of atmospheric trace gases (H2, CO, CH4, OCS, N2O, and NO).
Microbiol. Rev.
60:609-640 |
| 11. | Couteaux, M. M., and J. F. Darbyshire. 1998. Functional diversity amongst soil protozoa. Appl. Soil Ecol. 10:229-237[CrossRef]. |
| 12. |
Fernandez, A.,
S. Y. Huang,
S. Seston,
J. Xing,
R. Hickey,
C. Criddle, and J. M. Tiedje.
1999.
How stable is stable? Function versus community composition.
Appl. Environ. Microbiol.
65:3697-3704 |
| 13. | Ferris, M. J., and D. M. Ward. 1997. Seasonal distributions of dominant 16S rRNA-defined populations in a hot spring microbial mat examined by denaturing gradient gel electrophoresis. Appl. Environ. Microbiol. 63:1375-1381[Abstract]. |
| 14. | Fetzer, S., and R. Conrad. 1993. Effect of redox potential on methanogenesis by Methanosarcina barkeri. Arch. Microbiol. 160:108-113[CrossRef]. |
| 15. |
Grosskopf, R.,
S. Stubner, and W. Liesack.
1998.
Novel euryarchaeotal lineages detected on rice roots and in the anoxic bulk soil of flooded rice microcosms.
Appl. Environ. Microbiol.
64:4983-4989 |
| 16. |
Grosskopf, R.,
P. H. Janssen, and W. Liesack.
1998.
Diversity and structure of the methanogenic community in anoxic rice paddy soil microcosms as examined by cultivation and direct 16S rRNA gene sequence retrieval.
Appl. Environ. Microbiol.
64:960-969 |
| 17. | Head, I. M., J. R. Saunders, and R. W. Pickup. 1998. Microbial evolution, diversity, and ecology: a decade of ribosomal RNA analysis of uncultivated microorganisms. Microb. Ecol. 35:1-21[CrossRef][Medline]. |
| 18. | Heck, K. L., Jr., G. Van Belle, and D. Simberloff. 1975. Explicit calculation of the rarefaction diversity measurement and the determination of sufficient sample size. Ecology 56:1459-1461[CrossRef]. |
| 19. |
Henckel, T.,
M. Friedrich, and R. Conrad.
1999.
Molecular analyses of the methane-oxidizing microbial community in rice field soil by targeting the genes of the 16S rRNA, particulate methane monooxygenase, and methanol dehydrogenase.
Appl. Environ. Microbiol.
65:1980-1990 |
| 20. |
Jetten, M. S. M.,
A. J. M. Stams, and A. J. B. Zehnder.
1992.
Methanogenesis from acetate a comparison of the acetate metabolism in Methanothrix soehngenii and Methanosarcina spp.
FEMS Microbiol. Rev.
88:181-197.
|
| 21. | Klüber, H. D., and R. Conrad. 1998. Inhibitory effects of nitrate, nitrite, NO and N2O on methanogenesis by Methanosarcina barkeri and Methanobacterium bryantii. FEMS Microbiol. Ecol. 25:331-339[CrossRef]. |
| 22. |
Kopczynski, E. D.,
M. M. Bateson, and D. M. Ward.
1994.
Recognition of chimeric small-subunit ribosomal DNAs composed of genes from uncultivated microorganisms.
Appl. Environ. Microbiol.
60:746-748 |
| 23. | Krumbock, M., and R. Conrad. 1991. Metabolism of position-labelled glucose in anoxic methanogenic paddy soil and lake sediment. FEMS Microbiol. Ecol. 85:247-256[CrossRef]. |
| 24. | Lehmann-Richter, S., R. Grosskopf, W. Liesack, P. Frenzel, and R. Conrad. 1999. Methanogenic Archaea and CO2-dependent methanogenesis on washed rice roots. Environ. Microbiol. 1:159-166[CrossRef][Medline]. |
| 25. | Lelieveld, J., P. J. Crutzen, and F. J. Dentener. 1998. Changing concentration, lifetime and climate forcing of atmospheric methane. Tellus Ser. B Chem. Phys. Meteorol. 50:128-150. |
| 26. | Liesack, W., P. H. Janssen, F. A. Rainey, N. L. Ward-Rainey, and E. Stackebrandt. 1997. Microbial diversity in soil: the need for a combined approach using molecular and cultivation techniques, p. 375-439. In J. D. van Elsas, J. T. Trevors, and E. M. H. Wellington (ed.), Modern soil microbiology. Marcel Dekker, New York, N.Y. |
| 27. | Liu, W. T., T. L. Marsh, H. Cheng, and L. J. Forney. 1997. Characterization of microbial diversity by determining terminal restriction fragment length polymorphisms of genes encoding 16S rRNA. Appl. Environ. Microbiol. 63:4516-4522[Abstract]. |
| 28. | Ludwig, W., S. H. Bauer, M. Bauer, I. Held, G. Kirchhof, R. Schulze, I. Huber, S. Spring, A. Hartmann, and K. H. Schleifer. 1997. Detection and in situ identification of representatives of a widely distributed new bacterial phylum. FEMS Microbiol. Lett. 153:181-190[CrossRef][Medline]. |
| 29. | Ludwig, W., O. Strunk, S. Klugbauer, N. Klugbauer, M. Weizenegger, J. Neumaier, M. Bachleitner, and K. H. Schleifer. 1998. Bacterial phylogeny based on comparative sequence analysis. Electrophoresis 19:554-568[CrossRef][Medline]. |
| 30. | Mayer, H. P., and R. Conrad. 1990. Factors influencing the population of methanogenic bacteria and the initiation of methane production upon flooding of paddy soil. FEMS Microbiol. Ecol. 73:103-112[CrossRef]. |
| 31. |
Moré, M. I.,
J. B. Herrick,
M. C. Silva,
W. C. Ghiorse, and E. L. Madsen.
1994.
Quantitative cell lysis of indigenous microorganisms and rapid extraction of microbial DNA from sediment.
Appl. Environ. Microbiol.
60:1572-1580 |
| 32. | Patrick, W. H., Jr., and C. N. Reddy. 1978. Chemical changes in rice soils, p. 361-379. International Rice Research Institute, Los Banos, Phillipines. |
| 33. | Peters, V., and R. Conrad. 1996. Sequential reduction processes and initiation of CH4 production upon flooding of oxic upland soils. Soil Biol. Biochem. 28:371-382[CrossRef]. |
| 34. | Pitulle, C., and N. R. Pace. 1999. T-cloning vector for plasmid-based 16S rDNA analysis. Biotechniques 26:222-224[Medline]. |
| 35. | Ponnamperuma, F. N. 1972. The chemistry of submerged soils. Adv. Agron. 24:29-96. |
| 36. | Prinn, R. G. (ed.). 1994. Global atmospheric-biospheric chemistry, p. 1-18. In Global atmospheric-biospheric chemistry. Plenum, New York, N.Y. |
| 37. | Rainey, F. A., N. Ward, I. Sly, and E. Stackebrandt. 1994. Dependence on the taxon composition of clone libraries for PCR amplified, naturally occurring 16S rDNA, on the primer pair and the cloning system used. Experientia 50:796-797. |
| 37a. | Ramakrishnan, B., T. Lueders, R. Conrad, and M. Friedrich. Effect of soil aggregate size on methanogenesis and archaeal community structure in anoxic rice field soil. FEMS Microbiol. Ecol., in press. |
| 38. | Rotthauwe, J. H., K. P. Witzel, and W. Liesack. 1997. The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. Appl. Environ. Microbiol. 63:4704-4712[Abstract]. |
| 39. | Roy, R., and R. Conrad. 1999. Effect of methanogenic precursors (acetate, hydrogen, propionate) on the suppression of methane production by nitrate in anoxic rice field soil. FEMS Microbiol. Ecol. 28:49-61. |
| 40. | Roy, R., H. D. Klüber, and R. Conrad. 1997. Early initiation of methane production in anoxic rice soil despite the presence of oxidants. FEMS Microbiol. Ecol. 24:311-320[CrossRef]. |
| 41. | Schink, B. 1992. Syntrophism among prokaryotes, p. 276-299. In A. Balows, H. G. Trüper, M. Dworkin, W. Harder, and K.-H. Schleifer (ed.), The prokaryotes. Springer, New York, N.Y. |
| 42. |
Stackebrandt, E., and B. M. Goebel.
1994.
Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology.
Int. J. Syst. Bacteriol.
44:846-849 |
| 43. |
Suzuki, M. T.,
M. S. Rappe, and S. J. Giovannoni.
1998.
Kinetic bias in estimates of coastal picoplankton community structure obtained by measurements of small-subunit rRNA gene PCR amplicon length heterogeneity.
Appl. Environ. Microbiol.
64:4522-4529 |
| 44. | Tiedje, J. M., S. Asuming-Brempong, K. Nusslein, T. L. Marsh, and S. J. Flynn. 1999. Opening the black box of soil microbial diversity. Appl. Soil Ecol. 13:109-122[CrossRef]. |
| 45. |
von Wintzingerode, F.,
U. B. Goebel, and E. Stackebrandt.
1997.
Determination of microbial diversity in environmental samples pitfalls of PCR-based rRNA analysis.
FEMS Microbiol. Rev.
21:213-229[CrossRef][Medline].
|
| 46. | Wang, G. C. Y., and Y. Wang. 1997. Frequency of formation of chimeric molecules is a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes. Appl. Environ. Microbiol. 63:4645-4650[Abstract]. |
| 47. | Ward, D. M., M. M. Bateson, R. Weller, and A. L. Ruff-Roberts. 1992. Ribosomal rRNA analysis of microorganisms as they occur in nature. Adv. Microb. Ecol. 12:219-286. |
| 48. | Yao, H., and R. Conrad. 1999. Thermodynamics of methane production in different rice paddy soils from China, the Philippines and Italy. Soil Biol. Biochem. 31:463-473[CrossRef]. |
| 49. | Yao, H., R. Conrad, R. Wassmann, and H. U. Neue. 1999. Effect of soil characteristics on sequential reduction and methane production in sixteen rice paddy soils from China, the Philippines, and Italy. Biogeochemistry 47:269-295[CrossRef]. |
| 50. | Zehnder, A. J. B., and W. Stumm. 1988. Geochemistry and biogeochemistry of anaerobic habitats, p. 1-38. In A. J. B. Zehnder (ed.), Biology of anaerobic microorganisms. Wiley, New York, N.Y. |
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»