Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2006, p. 5342-5348, Vol. 72, No. 8
0099-2240/06/$08.00+0 doi:10.1128/AEM.00400-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Institute of Groundwater Ecology, GSFNational Research Center for Environment and Health, Neuherberg, Germany,1 Department of Bioremediation, UFZCenter for Environmental Research, Leipzig-Halle, Germany,2 Department of Biogeochemistry, Max Planck Institute for Terrestrial Microbiology, Marburg, Germany3
Received 17 February 2006/ Accepted 8 June 2006
|
|
|---|
|
|
|---|
New developments in stable-isotope approaches, particularly nucleic acid-based stable isotope probing (SIP) (for a recent review, see reference 4) have facilitated the targeted detection and identification of microbes metabolizing defined, mostly low-molecular-weight carbon sources in soil (6, 14, 17, 19). However, carbon sources in a soil environment usually are mixtures of different biomacromolecules, typical of biomass itself. The fate of biomacromolecule mixtures may differ considerably from that of individual components, because (i) a substrate mixture may result in competition, inhibition, and priming effects that change microbial degradation kinetics and (ii) some of the substrates may not be easily accessible to microbes due to chemical or physical protection mechanisms within the soil structure. Here, we followed the fate of fully 13C-labeled Escherichia coli cells as a proxy for labile biomass introduced into the soil (e.g., by manure or feces). The sequestration of this biomass-derived carbon pool was followed into the soil microbial food web, and microbial populations specifically active in its initial mineralization and assimilation were identified by SIP of rRNA.
|
|
|---|
13C-labeled biomass was added to microcosms at
1.3 x 108 cells per g dry weight (gdw) of soil as determined by microscopic cell counting and most-probable-number analysis in Luria-Bertani medium to check viability (7). The soil material used, the Ah horizon of a haplic phaeozem, was sampled at the long-term agricultural research site "Ewiger Roggenbau" (Halle, Germany). A detailed description of soil characteristics, incubation setup, experimental conditions, analytical procedures, and sampling is provided elsewhere (7). The microcosm bioreactor system was modified from that described in reference 31 and comprised circulating humidified airflow, CO2 traps, and temperature control (20°C). Pressure loss caused by CO2 removal was compensated for by addition of oxygen from a flexible gas bag acting as a self-dosing system. In order to obtain representative soil samples and to average microscale heterogeneities of the soil, multiple independent vertical soil cores were sampled at different locations in each reactor at every sampling time. The cores were pooled, homogenized, and subsampled for analysis. Soil samples were taken after 1 h and after 1, 2, 4, and 8 weeks of incubation. Unlabeled E. coli cells were used for a 12C control series. Samples for molecular analyses were frozen immediately (20°C) and stored until analysis. The reactor content was remixed after sampling to ensure homogeneous distribution of the remaining material. Trapped CO2 was analyzed for quantity and isotopic composition as published (18).
RNA extraction and centrifugation.
RNA was extracted from the soil samples using a previously described lysis protocol (13). Five hundred nanograms of RiboGreen-quantified rRNA extracts was density resolved by equilibrium density gradient centrifugation in cesium trifluoroacetate (CsTFA) (16). After centrifugation, the density of collected gradient fractions was determined refractometrically, and rRNA was precipitated for subsequent community analyses (13). Control gradients were conducted with rRNA from unamended soil, from soil incubated with unlabeled E. coli cells, and from unlabeled and 13C-labeled E. coli cells to calibrate the centrifugation system for the range of expected buoyant densities.
Quantitative and qualitative community analyses.
Bacterial and fungal rRNA templates precipitated and reeluted from gradient fractions were quantified by quantitative reverse transcription-PCR using assays previously described (13, 15) and standardized to relative ng µl1 units using a dilution series of pure-culture rRNA extracts of E. coli and Fusarium oxysporum. Terminal restriction fragment length polymorphism (T-RFLP) fingerprinting of bacterial communities was conducted as published elsewhere (14, 15). Additional fingerprints were generated as follows. Digested amplicons (
50 ng in 10 µl) were desalted by using DyeEx spin columns (QIAGEN). Desalted digests (1 µl) were mixed with 13 µl of Hi-Di formamide (Applied Biosystems) containing a 400-fold dilution of a 6-carboxy-X-rhodamine-labeled MapMarker 1000 ladder (Bio-Ventures, Murfreesboro, Tenn.), denatured (3 min at 95°C), cooled on ice, and size separated on a 3730 DNA analyzer (Applied Biosystems). Electrophoresis was performed with POP-7 polymer in a 50-cm capillary array under the following conditions: 10-s injection time, 2-kV injection voltage, 7-kV run voltage, 66°C run temperature, and 63-min analysis time. Electropherograms were analyzed using the GeneMapper 3.5 software package (Applied Biosystems). Cloning and sequencing of bacterial and fungal communities in resolved gradient fractions were done as described previously (15). Cloned reverse transcription-PCR amplicons were fully sequenced, resulting in
900 bp for bacteria and
600 bp for fungi. Phylogenetic trees were reconstructed using the ARB software package (11) and maximum-likelihood algorithms as described previously (12).
Nucleotide sequence accession numbers.
Sequence data were deposited with GenBank under accession numbers DQ643629 to DQ643792.
|
|
|---|
1.3 x 108 cells per gdw final concentration. The amount of carbon added corresponded to
36 µg 13C per gdw of soil, to
0.2% of the indigenous soil organic carbon (
19 mg Corg per gdw), to
26% of the soil bacterial biomass, and to
18% of the natural 13C content. Within 2 and 4 weeks of incubation, more than 30 and 40% of the added carbon was oxidized to 13CO2, respectively. Over the entire incubation period of 32 weeks, around 55% of the label was recovered as 13CO2, indicating that one-half of the added biomass was not oxidized but rather was sequestered into the indigenous soil microbiota and nonliving soil organic matter (7).
Identification of 13C-assimilating indigenous bacteria.
In order to identify distinct members of the soil microbiota that had been predominantly active in incorporating the added 13C label, 13C-enriched (or "heavy") rRNA molecules were separated from unlabeled bulk rRNA by density gradient centrifugation. RNA for gradients was extracted after 1 h and after 1, 2, and 8 weeks of incubation. At 1 h, highly 13C-labeled bacterial rRNA was present (Fig. 1) at a buoyant density (BD) matching that of fully labeled E. coli rRNA (Fig. 1, dashed box). In the controls, such "heavy" rRNA was not found for bacteria. Just 1 h after addition, the "heavy" rRNA detected was likely to originate directly from intact or lysed E. coli cells, as the amendment had certainly not been completely degraded already at this early time point. T-RFLP fingerprinting of bacterial populations in gradient fractions ranging from "light" to "heavy" confirmed that, 1 h after amendment, the 497-bp terminal restriction fragment (T-RF) characteristic of E. coli rRNA dominated the "heavy" fractions (Fig. 2). In the unlabeled E. coli control gradient for this time point, fingerprints of both "light" rRNA and unlabeled rRNA tailing into the "heavy" fractions were highly alike, and the 497-bp T-RF was not allocated to specific fractions.
![]() View larger version (26K): [in a new window] |
FIG. 1. Quantitative profiles of RNA-SIP centrifugation gradients. Shown are the distributions of bacterial (solid circles) and fungal (open circles) small-subunit (SSU) rRNA in density gradient fractions of rRNA extracted from soil microcosms incubated with unlabeled (12C, top row) and 13C-labeled E. coli biomass (bottom row) at successive times after amendment as specified. Domain-specific template distribution within gradient fractions was quantified by quantitative reverse-transcription-PCR (13). Dashed box, banding of E. coli pure-culture rRNA in density centrifugation gradients, showing the quantitative distribution of fully 13C-labeled (solid squares) and unlabeled (open squares) E. coli SSU rRNA in separate CsTFA density gradients.
|
![]() View larger version (25K): [in a new window] |
FIG. 2. T-RFLP fingerprints of density-resolved bacterial communities. Profiles were generated from representative "light" and "heavy" rRNA gradient fractions of the labeled (13C) and unlabeled (12C) soil incubations as specified. CsTFA BDs (g ml1) of fractions are given in brackets. T-RFs in "heavy" rRNA (after 1 week) characteristic for "heavy" fractions are specified, along with their respective fragment lengths, and were subsequently identified by cloning and sequencing.
|
After 1 week of incubation, rRNA from a "heavy" gradient fraction (BD of 1.815 g ml1) was cloned in order to identify the putative key players in biomass degradation represented by the "heavy" T-RFs. A total number of 99 clones were randomly selected, sequenced, and phylogenetically analyzed. As indicated already by the T-RFLP fingerprint of the corresponding gradient fraction (Fig. 2), the composition of phylotypes in the clone library was rather diverse (Table 1). However, a number of clones (n = 23) within distinct phylogenetic lineages were identified that were specifically represented within the T-RFs characteristic for "heavy" rRNA (Fig. 3; Table 1): a cluster of sequences of uncultivated
-Proteobacteria within the Xanthomonadaceae, most closely related to Lysobacter spp. (179- and 233-bp T-RFs), diverse sequence types within the Myxococcales (445- and 457-bp T-RFs), and some uncultivated members of the phylum Bacteroidetes (451- and 486-bp T-RFs). Furthermore, an abundant 96-bp T-RF was observed in "heavy" rRNA after 1 week of incubation, but only one Bacteroidetes clone was found to display such a fragment. Since bias by undersampling in clone libraries must be taken into account, it is also possible that this peak represents Bacteroidetes members specifically allocated to the "heavy" rRNA. Remarkably, all three lineages identified by T-RFLP to be selectively enriched in "heavy" rRNA (Lysobacter spp., Myxococcales, and Bacteroidetes) are well known to comprise unicellular gliding bacterial micropredators (22-24).
|
View this table: [in a new window] |
TABLE 1. Composition of the "heavy" bacterial clone librarya
|
![]() View larger version (29K): [in a new window] |
FIG. 3. Phylogenetic affiliation of representative bacterial 16S rRNA clones found to be specific for "heavy" rRNA fractions. Clones are marked in boldface, and their predicted T-RFs given in parentheses. Clones were generated from "heavy" rRNA (1.815 g ml1) after 1 week of degradation of added 13C-E. coli biomass in soil. The scale bar represents 10% sequence divergence. GenBank accession numbers of reference sequences are given.
|
Over time, differences between the fingerprints of "light" and "heavy" rRNA became less pronounced (Fig. 2). This indicates that the 13C label was no longer specifically retained by certain bacterial populations but rather was spreading, via trophic networks, within the entire soil microbiota. Highly labeled bacterial rRNA remained detectable even after 8 weeks of incubation (Fig. 1). Of the
45% of added labeled carbon remaining in the soil after this time (7), around half was still present within the soil microbial community itself and persisted within the living fraction of the soil carbon pool. The other half of the remaining carbon appears to be present in nonliving soil organic matter, as suggested by phospholipid fatty acid analysis (7).
Dynamics of fungal rRNA within gradient fractions.
In contrast to the active bacteria, fungal populations in the soil appeared to be much more dynamic in biomass turnover. Here, the addition of the E. coli cells caused a dramatic increase of fungal rRNA detectable in relation to bacterial templates within 1 h (Fig. 1). However, intriguingly, most fungal rRNA was intrinsically distributed within "heavier" gradient fractions than bacterial templates in both 13C-labeled and unlabeled E. coli incubation series. Distinctions between the two treatments were not significant, with most fungal templates after 1 h banding at
1.81 g/ml in the 13C gradient and slightly below in the 12C gradient. The amounts of fungal rRNA measured were largest after 1 week of incubation but decreased subsequently, indicating a decline of the stimulatory effect by labile biomass addition.
In order to preliminarily identify the detected fungi, three small clone libraries were generated from different fractions after 1 h and after 2 weeks (Table 2). All libraries were dominated by Microascaceae closely related to Petriella and Pseudallescheria spp. (Fig. 4). A further frequent group of clones was related to Chaetomium spp., but significant differences between density-resolved fungal communities were not evident from these libraries.
|
View this table: [in a new window] |
TABLE 2. Composition of fungal clone librariesa
|
![]() View larger version (23K): [in a new window] |
FIG. 4. Phylogenetic affiliation of fungal 18S rRNA clones from gradient fractions. Representative clones within the Sordariomycetes retrieved from "heavy" (Y04Fu and I04Fu clones) and "light" (I08Fu clones) fractions are shown in boldface. The scale bar represents 10% sequence divergence. GenBank accession numbers of reference sequences are given.
|
|
|
|---|
The experiment was conducted using only one soil type and one bacterial biomass substrate. Also, due to the laborious nature of rRNA-SIP, replicate centrifugation gradients or data from replicate incubations are not presented. Therefore, and without testing other soils and settings, quantitative estimates of the role of the identified microbes in terrestrial C fluxes cannot be deduced at present. The amount of 13C label amended to the soil (36 µg biomass 13C per g soil) corresponded to only 2.8 µmol of 13C per g soil and thus represents a small amount of label in comparison to other rRNA-SIP studies (32). However, the amount of biomass added certainly will cause alterations in the activity and composition of intrinsic microbiota (e.g., priming effects [8]), as shown here in both unlabeled and 13C-labeled time series, for bacteria and fungi. Nevertheless, the amount of labeled biomass added represented only 0.2% of the indigenous soil organic carbon and therefore can still be considered "close to in situ", especially regarding the rather patchy applications of labile bacterial biomass introduced into agricultural soil environments by, e.g., manure or feces.
Despite this small substrate pulse, diverse members of the Myxococcales, the Xanthomonadaceae, and the Bacteroidetes were identified to have incorporated the added carbon with high specificity after 1 week of incubation. All of these indigenous soil bacteria are more or less closely related to organisms capable of gliding motility. Many members of these lineages have been described as micropredators, capable of actively lysing their living bacterial prey (22-24). It has been suggested that they play a key role in the turnover of biomass carbon and other labile organic material in soil ecosystems (21), but direct and cultivation-independent evidence has not been available yet. Evidently, it is not possible to infer a functional trait as complex as gliding motility by mere 16S rRNA sequence similarity. However, this functional trait would provide the detected bacteria with preferential access to the added cellular substrate and thus a unifying competitive advantage, despite being phylogenetically diverse. The data suggest that these microbes were selectively incorporating 13C from the added cellular organic material, regardless of whether actively lysing the added 13C-E. coli cells or degrading dead biomass.
In contrast to bacteria, the 13C label could not be clearly tracked into fungal components of the soil microbial food web. Intriguing questions are why the detected fungal rRNA was mostly distributed in fractions "heavier" than those for bacterial rRNA also in the 12C control incubations and why the addition of 13C-E. coli did not cause a more pronounced shift in BD for fungal templates. Earlier studies report the location of unlabeled Fusarium sp. rRNA at densities similar to unlabeled bacterial rRNA (13), and shifts in BD for indigenous soil fungal communities upon 13C labeling have been observed (15). The contrasting observations of our present report illustrate the urgent need for a more exhaustive investigation of the distribution of unlabeled and 13C-labeled eukaryotic rRNA templates than can be expected for CsTFA centrifugation gradients.
Fraction-specific differences in the qualitative assemblage of fungal rRNA templates were not evident. It is unlikely that rRNA-based fungal fingerprinting techniques (for a recent summary, see reference 1) would have revealed clear distinctions. Rather, we hypothesize that the total indigenous fungal community was stimulated from the decomposing E. coli biomass, and the Microascaceae and Chaetomiaceae were most abundant. However, they did not incorporate the added label into rRNA to an extent clearly differentiable from the controls. This clearly prevents any predatory interpretation of these fungi, albeit this has been described for others (27) and cannot be totally excluded. Information on the ecological role of the Microascaceae in natural habitats is scarce (2), and such fungi are mainly known as opportunistic pathogens (26). But generally, fungi are as likely as bacteria to utilize decomposing biomass. Here, the most surprising finding was the apparent small time span required by the Petriella relatives to form increased amounts of rRNA after addition of E. coli biomass, thus highlighting a very dynamic role of these fungi in the degradation of available microbial biomass in the soil.
A lack of specific quantification and fingerprinting assays has prevented us from including a further group of microeukaryotes in the analyses: the protozoa. In addition to bacteria and fungi, protozoa would certainly have contributed to the degradation of the added E. coli biomass, but unfortunately detailed information cannot be provided for this experiment. Current efforts focus on the development of suitable assays that make distinct protozoan lineages accessible for (RNA-) SIP studies, and future work will show whether 13C label stemming from added bacterial biomass can be traced also to these eukaryotic lineages. Already in an earlier study, terrestrial cercomonads were shown to have incorporated carbon stemming from 13C-labeled methylotrophs (15); however, their detection was not targeted.
In conclusion, with the help of rRNA-SIP and a complex 13C substrate, we provide insight into the "black box" of a soil microbial food web and direct evidence for the identity of specific key players therein. Putatively motile bacterial micropredators were shown to play a pivotal role in the sequestration and turnover of the added labile biomass. Bacterial gliding motility may be a key trait controlling ecological niche positioning of these microbes, as it grants preferential access to attractive nondiffusible substrates. In consequence, the detected "heavy" bacterial lineages could be considered an ecologically unified functional clade. Further investigations regarding other soil types and water availabilities may show whether this holds true as a general concept for soil bacterial biomass turnover. Also, experiments with more-recalcitrant carbon sources may show whether RNA-SIP is applicable to the turnover of more-stable carbon pools, such as humified soil organic matter.
We thank Bianca Pommerenke (MPI Marburg) for dedicated technical assistance and an anonymous reviewer for valuable suggestions.
|
|
|---|
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»