Applied and Environmental Microbiology, September 2004, p. 5057-5065, Vol. 70, No. 9
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.9.5057-5065.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Bioresources, ARC Seibersdorf Research GmbH, Seibersdorf,1 Institute for Forest Ecology, Federal Office and Research Centre for Forests, Vienna, Austria2
Received 3 November 2003/ Accepted 21 May 2004
|
|
|---|
-Proteobacteria (20%) and the Holophaga/Acidobacterium group (12%). Clones in libraries from oak-hornbeam and spruce-fir-beech forest soils were mainly related to the Holophaga/Acidobacterium group (28 and 35%), followed by members of the Verrucomicrobia (24%) and the
-Proteobacteria (27%), respectively. The soil bacterial communities in forests with distinct vegetational and soil chemical properties appeared to be well differentiated based on 16S rRNA gene phylogeny. In particular, the outstanding position of the Austrian pine forests, which are determined by specific soil conditions, was reflected in the bacterial community composition. |
|
|---|
A range of Austrian natural forests have been surveyed with respect to plant species composition, soil chemical parameters, microbial biomass, and microbially mediated nutrient transformation. Thereby, distinct soil chemical and microbiological conditions have been found to prevail in contrasting forest types (13, 14, 29; E. Hackl, M. Pfeffer, C. Donat, G. Bachmann, and S. Zechmeister-Boltenstern, submitted for publication). However, these findings have not yet been set in relation to the key players in soil nutrient turnover, i.e., the microorganisms inhabiting the different environments. The present study is focused on the bacterial communities of the forest soils and their interrelationships with vegetation and soil environmental conditions.
Soil bacteria are an essential component of the biotic community in natural forests and they are largely responsible for ecosystem functioning, because they participate in most nutrient transformations. Although the main diversity of life has been proven to be microbial, the vast majority of soil bacteria still remain unknown due to the fact that only a minor percentage of naturally occurring microorganisms can be cultured (23). To survey the phylogenetic diversity present in bacterial communities, 16S rRNA gene-based molecular techniques have commonly been applied, which circumvent limitations of culture-based studies. Forest soil bacterial diversity has been characterized by DNA sequence analysis of 16S rRNA genes, derived for instance from soils under Amazonian forest (3), under pinyon-juniper woodlands (7), and under lodgepole pine trees (4). While sequence analysis of clones from 16S rRNA gene clone libraries provides the most detailed, reliable information about the composition of microbial communities, this method is time-consuming and expensive. Terminal restriction fragment length polymorphism (TRF, or T-RFLP) analysis, in contrast, is a DNA-based method which allows rapid comparison of complex bacterial communities. TRF analysis has been applied for a wide range of environmental samples, including agricultural soils (19, 25), grassland (18), and forest soils (7) and biological soil crusts (24). Comparison of bacterial communities by TRF analysis has proven to provide information consistent with analysis of clone libraries (7).
In the present study, 16S rRNA gene-based community analysis has been applied to compare the diversity and composition of the bacterial communities in soils under six natural forest stands which differ in soil and vegetation characteristics, soil microbial biomass, and microbial nutrient turnover rates (13, 14). These include two oak-hornbeam and two spruce-fir-beech forests, which represent zonal vegetation types of the region, as well as two pine forests, which are found on azonal soils that are typically stony and water and nutrient limited. The aim of the study was to relate the bacterial community composition to the already-known characteristics of these forest soils.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Site characteristics (14) and TRF diversity data from six natural forest soils
|
|
View this table: [in a new window] |
TABLE 2. Soil chemical and microbiological characteristics of six natural forest soilsc
|
Soil sample collection.
The mineral soils under the six forest stands were sampled in Autumn 2002. Ten individual samples (0- to 10-cm depth) were collected from each forest stand at intervals of 5 m along transects of 50 m in length, which had been established during earlier studies (13, 14). The soil samples were taken to the laboratory in cooling boxes. They were sieved (2-mm mesh) and then kept at 20°C until analysis.
TRF community profiles of soil samples.
The individual soil samples were well homogenized, and DNA was extracted from approximately 0.5 g of each sample by using the UltraClean soil DNA isolation kit (MoBio Laboratories, Inc.). This DNA was then used as template for PCR with 16S rRNA gene primers 8f (10) labeled with 6-carboxyfluorescein (MWG) at the 5' end and 926r (22). Each 50-µl mixture contained 10 to 20 ng of extracted DNA, 1x reaction buffer (Gibco BRL), 200 µM (each) dATP, dCTP, dGTP, and dTTP, 3 mM MgCl2, a 0.2 µM concentration of the primers, and 2.5 U of Taq DNA polymerase (Gibco BRL). Amplifications were performed in a PTC-100 thermocycler (MJ Research, Inc.) with an initial denaturing step of 5 min at 95°C followed by 30 cycles of 30 s at 95°C, 1 min of annealing at 53°C, and 10 min of extension at 72°C. Three amplification products were prepared from each DNA sample and were then pooled to reduce PCR bias. Approximately 200 ng of fluorescently labeled PCR amplification products was digested with the restriction enzyme AluI (Gibco BRL) and subsequently purified by passage through Sephadex G-200 columns. Aliquots of 5 µl were mixed with 2 µl of loading buffer (diluted five times in deionized formamide; Fluka) and 0.3 µl of the DNA fragment length standard (Rox 500; PE Applied Biosystems Inc., Foster City, Calif.). Mixtures were denatured for 2 min at 92°C, and fluorescently labeled TRFs were then detected using an ABI 373A automated DNA sequencer (PE Applied Biosystems, Inc.) in the GeneScan mode. Lengths of labeled fragments were determined by comparison with the internal standard using the GeneScan 2.5 software package (PE Applied Biosystems Inc.).
Analysis of TRF profiles.
The DNA quantity analyzed for each of the 60 TRF community profiles was compared and standardized to the lowest quantity, according to the method of Dunbar et al. (8). TRFs of 39 to 500 bp in length and with heights of
50 fluorescence units were included in the analysis.
TRF data from individual soil samples were subjected to principal component analysis (PCA) to elucidate major variation patterns. Both fragment length and peak height were used as parameters for profile comparison. The scores of the first two components were subsequently used to compare differences between the forest stands in the TRF patterns with differences in soil chemical and microbiological properties by calculating Spearman's correlation coefficients (P
0.05). Average TRF data from 10 community profiles per forest site were subjected to cluster analysis using the Euclidean distance measure. A dendrogram was constructed using an unweighted pair-group method with arithmetic mean (UPGMA). All statistical procedures were performed using Statistica version 5.1 (1997 ed.).
Diversity statistics were calculated from standardized profiles of individual soil samples by using the number and height of peaks in each profile as representations of the number and relative abundance of phylotypes, as defined by Dunbar et al. (7). Phylotype richness (S) was calculated as the total number of distinct TRF sizes (between 39 and 496 bp) in a profile. The Shannon-Weiner diversity index (26) was calculated as follows: H =
(pi) (log2pi), where p is the proportion of an individual peak height relative to the sum of all peak heights. Simpson's diversity index (27) was calculated as follows: D = 1 [
(pi)2]. Evenness (26) was calculated from the Shannon-Weiner diversity function as follows: E = H/Hmax, where Hmax = log2(S). Diversity data were tested for significant differences among forest sites by using the nonparametric Mann-Whitney U-test (Statistica 5, version 97).
16S rRNA gene clone libraries.
Clone libraries were constructed for the forest sites Kolmberg (oak-hornbeam forest), Rothwald (spruce-fir-beech forest), and Stampfltal (Austrian pine forest). The DNA templates used for construction of clone libraries were the same DNA preparations from which 16S rRNA genes for TRF analysis were amplified. Twenty independent PCRs (10 subsamples in duplicate) were performed for each sampling site by using the primers 8f (10) and 1520r (20) and the conditions described above. 16S rRNA gene amplicons were pooled and purified by gel extraction using the Concert rapid gel extraction system (Gibco BRL). Amplicons were ligated into the TpCR 4-TOPO vector (Invitrogen), and ligation products were transformed into Escherichia coli DH5
-T1 (Invitrogen) according to the manufacturer's instructions. Approximately 50 colonies were picked for each study site, each suspended in 80 µl of Tris-EDTA buffer (pH 8.0), and incubated at 100°C for 10 min. The lysis product (1.5 µl) was added to PCR mixtures using the primers M13f (5'-GTAAAACGAGGCCAG-3') and M13r (5'-CAGGAAACAGCTATGAC-3') and the conditions described above to amplify cloned inserts. PCR products were purified with the QIAquick PCR purification system (QIAGEN) according to the manufacturer's instructions and used as templates in sequencing reactions.
DNA sequencing.
Partial 16S rRNA gene sequencing was performed by applying the primer 518r (21) by the dideoxy chain termination method (24a) by using an ABI 373A automated DNA sequencer and the ABI Prism Big Dye terminator cycle sequencing kit (PE Applied Biosystems Inc.). Sequences were subjected to BLAST analysis (1) with the National Center for Biotechnology Information (NCBI) database. Sequences were checked for the presence of chimeric artifacts using the Chimera Check program of the Ribosomal Database Project (5).
TRF analysis of 16S rRNA gene clones.
16S rRNA genes of the clones were amplified by PCR using the M13f and M13r primers (Invitrogen) and the conditions described above. Each reaction mixture (1.5 µl) was used as template in a second PCR containing forward primer 8f (10) labeled with 6-carboxyfluorescein (MWG) at the 5' end and 926r (22). Reaction conditions were the same as described above except that only 16 cycles of PCR were used instead of 30. Fluorescently labeled PCR amplification products from four to six clones (differing in TRF size) were mixed and collectively digested with the restriction enzyme AluI (Gibco BRL) and subsequently purified by passage through Sephadex G-200 columns. Aliquots of 1 µl were mixed with 3 µl of loading buffer (diluted five times in deionized formamide [Fluka]) and 0.3 µl of the DNA fragment length standard (Rox 500; PE Applied Biosystems Inc.). Mixtures were denatured for 2 min at 92°C, and fluorescently labeled TRFs were then detected using an ABI 373A automated DNA sequencer (PE Applied Biosystems Inc.) in the GeneScan mode. Lengths of labeled fragments were determined by comparison with the internal standard by using the GeneScan 2.5 software package (PE Applied Biosystems Inc.).
Nucleotide sequence accession number.
The nucleotide sequences determined in this study have been deposited in the NCBI database under accession numbers AY395085 to AY395222.
|
|
|---|
![]() View larger version (32K): [in a new window] |
FIG. 1. Representative TRF profiles of bacterial communities from mineral soils under three Austrian natural forests: oak-hornbeam forest (A), spruce-fir-beech forest (B), and Austrian pine forest (C). Fragments corresponding to dominant phylogenetic groups represented by 16S rRNA gene clones derived from the same soils are indicated and labeled with the respective peak heights.
|
50 FU averaged between 82 (Kolmberg oak-hornbeam forest) and 101 (Merkenstein pine forest). The term phylotype richness was applied to the number of peaks with intensities of
50 FU and among all soils showed highest values for the Merkenstein pine forest soils. The Shannon-Weiner index and the Simpson index were calculated as a measure of diversity, which showed only a few differences among the forest soils. However, values of the Simpson index were significantly higher for the Stampfltal pine forest. No significant differences were found in the evenness of distribution of TRFs among the forest sites studied (Table 1).
To elucidate major distribution patterns, a PCA of TRF data was performed (Fig. 2). In the score plot of the PCA, the soils of the Austrian pine forests were separated from the other forest soils along the first principal component (PC 1). TRF profiles of the pine forest soils contained fragments of 64, 65, 66, 69, 70, 71,74, 75, 124, 125, 136, 137, 140, 176, 177, and 217 bp in high abundances (>300 FU), which were heavily loaded on PC 1. Some of these fragment lengths were found to be represented by 16S rRNA gene clones from the Stampfltal pine forest soil, for instance, members of the high-G+C gram positives (70 bp),
-Proteobacteria (124 bp), and the Cytophaga/Flexibacter/Bacteroides group (140 bp), as illustrated in Fig. 1. Compared with the pine forest profiles, the profiles of the oak-hornbeam and spruce-fir-beech forest soils showed higher intensities in the fragments of 168, 220, 236, 237, and 244 bp, which had negative loadings along PC 1. 16S rRNA gene clones from the Kolmberg oak-hornbeam forest and the Rothwald spruce-fir-beech forest which corresponded with these fragment lengths were all affiliated with the Holophaga/Acidobacterium division (Fig. 1). The soils of both oak-hornbeam forests as well as the Rothwald spruce-fir-beech forest clustered together, whereas the Neuwald spruce-fir-beech forest soils were set apart from these soils along the second principal component (PC 2). Compared with profiles from the other forest soils, TRF profiles from the Neuwald forest had higher intensities of fragments with 86, 87, 88, 152, 160, 161, and 254 bp, which had positive loadings on PC 2. Differences between the forest stands in the TRF patterns corresponded with differences in the soils' chemical and microbiological properties, presented in Table 2. Significant correlations were found between PC 1 and soil pH (r = 0.76; P = 0.0001), soil C-to-N ratio (r = 0.79; P = 0.0000), the ratio of fungal to bacterial phospholipid fatty acids (PLFAs) (r = 0.83, P = 0.0000), and microbial biomass measured as SIR-C g of organic C1 (r = 0.56; P = 0.0079), FE-N g of organic C1 (r = 0.50; P = 0.0219), and total PLFAs g of organic C1 (r = 0.72; P = 0.0002). No significant correlations have been found between PC 2 and soil chemical and microbiological properties.
![]() View larger version (22K): [in a new window] |
FIG. 2. Score plot (A) and loading plot (B) of PCA of TRF data. In the loading plot, only fragments with heights of 300 FU are shown. , oak-hornbeam forest (Johannser Kogel); , oak-hornbeam forest (Kolmberg); , spruce-fir-beech forest (Rothwald); , spruce-fir-beech forest (Neuwald); , Austrian pine forest (Stampfltal); , Austrian pine forest (Merkenstein).
|
![]() View larger version (15K): [in a new window] |
FIG. 3. Dendrogram of TRF profiles from six natural forest soils generated from Euclidean distance values, using the UPGMA method. Cluster analysis was performed with average TRF values of 10 individual subsamples from the study sites. JE, Johannser Kogel; K, Kolmberg; R, Rothwald; N, Neuwald; St, Stampfltal; Me, Merkenstein.
|
The relative abundances of clones affiliated to various bacterial groups differed for the three forest soils. In the clone library derived from the Kolmberg oak-hornbeam forest, clones affiliated with the Holophaga/Acidobacterium division were most abundant, comprising 28% of total clones. Clones representing members of the Verrucomicrobia division were the second most abundant (24%), and those affiliated with the Cytophaga/Flexibacter/Bacteroides group were third most abundant (11%). Members of the
-Proteobacteria accounted for 9% of the sequences. The remaining clones were affiliated with high-G+C gram positives, Planctomycetes, ß-,
-, and
-Proteobacteria, and candidate divisions BD and OP10. One clone (clone L13K) could not be classified into the known bacterial divisions. In the clone library from the Rothwald spruce-fir-beech forest, the majority of clones were also classified as Holophaga/Acidobacterium (35%). Members of the
-Proteobacteria accounted for 27%. Minor parts of the library belonged to the Verrucomicrobiales (10%) and Planctomycetales (10%) groups as well as to members of the ß-,
-, and
-Proteobacteria (each 6%). The clone library from the Stampfltal Austrian pine forest differed from the other libraries in that the majority of clones was classified as high-G+C gram positives (49%). Again, sequences that affiliated with the
-Proteobacteria were highly abundant (20%). Twelve percent of the clones were represented by Holophaga/Acidobacterium, and the remaining 16S rRNA genes showed the highest similarity to the ß- and
-Proteobacteria, the Cytophaga/Flexibacter/Bacteroides group, and to the Verrucomicrobia.
In Table 3, the phylogenetic assignments are given for clones of the three libraries which corresponded to dominant TRFs in the community profiles. TRF lengths observed in profiles have been reported to differ from sequence-determined TRF lengths, depending on subtle differences in molecular weight, either from purine content or dye label (16). To allow comparison of TRF lengths of the clones with TRFs in community profiles, the true fragment lengths of the clones were therefore determined by subjecting them to TRF analysis. In addition to the observed TRF sizes of the clones, the TRF sizes predicted from sequence information are also shown in Table 3, providing information on TRF drift.
|
View this table: [in a new window] |
TABLE 3. Phylogenetic assignment of clones from mineral soil from three different natural forests corresponding to highly abundant TRFs in community profilesa
|
|
|
|---|
Analysis of 16S rRNA gene clone libraries allowed us to identify the most dominant members of the bacterial communities. In agreement with the results obtained by TRF analysis, clones contained in the library from the Stampfltal pine forest were affiliated with different bacterial groups than those in the libraries derived from the Kolmberg oak-hornbeam forest and the Rothwald spruce-fir-beech forest. Representing 49% of the total clones, the library from the pine forest was highly dominated by members of the high-G+C gram-positive bacteria. Similarly, in community fingerprints the corresponding TRFs were highly abundant. High-G+C gram positives were of minor importance in the Kolmberg oak-hornbeam forest and had no representatives in the library obtained from the Rothwald spruce-fir-beech forest. High-G+C gram positives were also most prevalent in 16S rRNA clone libraries from mineral forest soil samples in British Columbia, where lodgepole pine (Pinus contorta Douglas) is a major tree species (2). Nevertheless, clones representing this bacterial group were far less abundant in the British Columbian soil than in the Stampfltal samples. In clone libraries from organic forest soils (2) and from lodgepole pine rhizosphere soil (4) in British Columbia, however,
-Proteobacteria were dominant. In the Stampfltal library clones which affiliated with the
-Proteobacteria were second most abundant. High-G+C gram-positive bacteria were not detected or were of minor prevalence in clone libraries from forests other than pine forests (3, 4, 6, 9).
High-G+C gram-positive bacteria comprise a large variety of bacteria which have commonly been called actinomycetes because they include members which grow in a filamentous manner. The high prevalence of clones affiliated with this bacterial group in the Stampfltal clone library is in agreement with the fact that high concentrations of tuberculostearic acid (PLFA 10Me18:0) have previously been found in soils of the Stampfltal and Merkenstein pine forests (Hackl et al., submitted). This fatty acid is common in actinomycetes and has frequently been used to indicate actinomycetes in soils (17). Sequences retrieved from the Stampfltal pine forest soil were distributed among several phyla within the high-G+C gram positives, including filamentous actinomycetes. However, organisms related to Streptomycesa common soil microorganismwere not detected. Filamentous growth may facilitate nutrient acquisition, which would be of advantage in the pine forest soils, which are prone to desiccation and nutrient imbalances (14). In addition, high quantities of PFLAs indicative of fungi have been reported in these soils (Hackl et al., submitted). Various Stampfltal 16S rRNA gene sequences grouped with clones obtained from Australian desert soils (15) and fell into a deeply branching line of actinobacteria which was most related, though moderately, to the genus Rubrobacter.
Pine forest soils represent habitats for soil bacteria which are prone to nutrient imbalances, and of all forest soils studied they are most limited in water availability (13, 14; Hackl et al., submitted). Nevertheless, by using TRF profiling high numbers of highly dominant bacterial groups were detected as phylotypes in these soils. In particular, the Stampfltal pine forest appeared to host a highly diverse soil bacterial community, as indicated by diversity indices calculated from TRF community data. In contrast, the Rothwald spruce-fir-beech forest soil, which is most nutrient rich and high in water content, exhibited only average values of phylotype richness and community diversity. It has to be considered that diversity measures based on TRFs only give information on the diversity of the numerically dominant bacterial groups and that minor but potentially highly diverse subsets of the soil bacterial communities may not have been detected by the TRF analysis. However, the bacterial communities of the pine forest soils are set apart from those in the other forest soils in that they are composed of a multitude of highly abundant bacterial groups. A possible explanation for this result may be that these soils host a huge variety of bacteria which are able to cope with a wide range of environmental conditions, including various levels of nutrient availability. A major portion of these bacteria might be active only in short periods after rainfall or when nutrients are episodically released. Indeed, the soil moisture content and the C/N ratio have been identified as the factors which are mainly responsible for within-site variability in microbial activity in pine forest soils (13).
Bacterial communities in oak-hornbeam and spruce-fir-beech forest soils.
The oak-hornbeam and the spruce-fir-beech forest soils contained several highly abundant bacterial groups which were less dominant in the pine forest soils. Comparison of community TRF profiles with the representative 16S rRNA gene clones implied that these bacterial groups are composed of Holophaga/Acidobacterium members. The view that members of the Holophaga/Acidobacterium group may constitute a major part of the bacterial communities in oak-hornbeam and spruce-fir-beech forest soils is further supported by the fact that 16S rRNA gene clones affiliated with this group were the most abundant in libraries from both the Kolmberg and the Rothwald forest soils. Among forest soils, soils under pinyon pine-juniper woodlands have been reported to contain acidobacteria in highest numbers (6). Sequences belonging to the Holophaga/Acidobacterium division have been retrieved from a multitude of habitats, including forest soils (e.g., references 2, 4, and 6) and have revealed extraordinary phylogenetic diversity within this group (9, 25, 28). So far, only a few cultivated representatives of the Holophaga/Acidobacterium division exist, and therefore ecological and physiological data are very limited. However, the phylogenetic diversity and ubiquity of the acidobacteria suggest that they show high metabolic versatility.
The bacterial communities from the two oak-hornbeam forest soils (Johannser Kogel and Kolmberg) clustered closely in comparisons of forest soil TRF profiles, indicating similarity in community composition. Conformities in the profiles included high intensities in fragments of 148 bp in length. In the library from the Kolmberg forest soil, several clones were present which corresponded to 148-bp fragments and which were affiliated with the Verrucomicrobia group. Clones representing Verrucomicrobia members were second most abundant in the Kolmberg clone library, suggesting that Verrucomicrobia organisms may occupy important roles in the oak-hornbeam forest soils.
Similarities in bacterial community composition were observed within the two oak-hornbeam forests as well as within the two pine forests, which might reflect a selective influence of the dominant plant species at these sites on the soil bacteria. However, since the vegetation composition of the forests under study is supposed to be differentiated according to site conditions, its influence on the soil bacterial communities is closely connected with effects of other site conditions, such as climate, geological substrate, and soil type. In addition, it may be assumed that plant-specific effects through the input of litter may be more pronounced in the organic than in the mineral soil layer, which was studied. Thus, the community patterns observed in our study probably reflect the interplay of various factors, including plant species composition as well as soil type and other site conditions.
In contrast to the soil bacterial communities from the oak-hornbeam and pine forests, which seemed to be similar in composition within forest type, major differences were apparent when the soil bacterial communities from the Rothwald and the Neuwald spruce-fir-beech forests were compared. This may have been related to differences between the two spruce-fir-beech forest sites in soil type and geological substrate. While the Rothwald forest soil has been defined as planosol situated on dolomite and has been found to be especially nutrient rich with pH values of about 5, the Neuwald forest soil has been classified as stagnic luvisol situated on sandstone substrate with slow nutrient release and lower pH values (14). Compared with the Rothwald forest soil, considerably lower nutrient concentrations and slower nutrient turnover rates have been measured in the Neuwald forest soil, although both forests have a similar vegetation composition and similar topological characteristics (14). Results of the present study suggest that distinct bacterial communities are involved in nutrient cycling in the two soil environments.
Conclusions.
In the present study, 16S rRNA gene analysis proved appropriate for detecting compositional differences in the bacterial communities inhabiting various natural forest soils. Differences in bacterial community composition among forest soils were found to correspond with differences in other microbiological and soil chemistry characteristics. Our data suggest that the composition of the bacterial communities inhabiting natural forest soils is related to specific environmental conditions which go along with various types of natural forest vegetation. The comparative analysis of bacterial communities in soils under different types of natural forest has provided important baseline information about bacterial diversity and composition in undisturbed forest ecosystems. These data will allow us to go forward with studies in which the effects of management practices or environmental pollution on the soil bacterial communities and their functions will be evaluated.
We thank Alexandra Weilharter for her excellent assistance with the TRF and sequence analyses.
|
|
|---|
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»