Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2009, p. 308-315, Vol. 75, No. 2
0099-2240/09/$08.00+0 doi:10.1128/AEM.01739-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Applied Chemistry and Microbiology and Department of Applied Biology, University of Helsinki, Helsinki, Finland
Received 29 July 2008/ Accepted 23 October 2008
|
|
|---|
|
|
|---|
Recent studies show that in addition to boreal forest soils, group 1.1c Crenarchaeota also inhabit the soil of the Austrian alps (19, 20), although group 1.1b Crenarchaeota are more abundant. In the alpine soil, only group 1.1b Crenarchaeota were found in young soil uncovered by ice for only 4 years and lacking significant vegetation. However, 1.1c Crenarchaeota formed an integral part of the community in soil that had been free from ice cover for over 135 years and had a stable vegetation cover. An even richer diversity of 1.1c Crenarchaeota was observed in 9,500-year-old vegetated soil. For the vegetated and certainly mycorrhizospheric soils of the Rocky Mountains (Colorado), Oline et al. (22) showed that group 1.1c Crenarchaeota colonized the soil from the foot hill forest (1,800 m above sea level) to the alpine tundra (3,400 m above sea level), but 1.1b Crenarchaeota were still the most abundant archaea.
Crenarchaeota have been found to be associated with the roots of a variety of plant species (5, 8, 31). In a recent study (34), archaea were surveyed in roots of 13 different plant genera, collected from their natural habitats, including turf, prairie, agricultural, and forest sites in a temperate region. All the archaea detected belonged to the 1.1b group, and the diversity of these archaea did not vary significantly between plant genera from the same site, only between plant roots and bulk soil. There was more variability between the different geographical locations than between the plant genera.
Pinus sylvestris is one of the dominant tree species in Finnish forests. When grown in boreal forest humus, nonmycorrhizal P. sylvestris fine roots harbor high numbers of bacteria (37), but they appear to be nearly devoid of archaea (2, 3). In nearly all cases in our previous work, archaea were detected in P. sylvestris fine roots only when they were colonized by mycorrhizal fungi (3). Sliwinski and Goodman (34) were able to show archaeal group 1.1b colonization in roots of Pinus species in temperate forests but did not report the mycorrhizal status of the sampled roots.
In the present study, we surveyed the archaeal population compositions in two coniferous and two deciduous tree species developed in boreal forest humus. As the ectomycorrhizal fungus Paxillus involutus is symbiotic with all four tree species, it was used to test the effect of mycorrhizal colonization on the archaeal populations on the roots.
|
|
|---|
Ten individual mycorrhizal root tips (0.02 g [fresh weight {fw}]) were collected from each P. involutus inoculated seedling, and 10 to 15 individual fine roots (0.02 g [fw]) were collected from each noninoculated seedling. The mycorrhizas or fine roots obtained from one seedling were considered one mycorrhizal and one fine-root sample, respectively. Totals of 20 mycorrhizal and 20 fine-root samples, representing five replicates of each tree species, were collected with sterile forceps from the microcosms for DNA extraction. Five samples of 0.2 g (fw) uncolonized forest humus outside the mycorrhizosphere were also harvested from different microcosms. The samples were kept at –80°C until DNA extraction.
DNA extraction and PCR.
The mycorrhizal and fine-root samples were ground with sterilized quartz sand in 50 µl sterile double-distilled water to disrupt the plant and fungal tissues prior to DNA extraction. The humus samples did not require any preextraction treatment. DNA was extracted with a MoBio soil DNA extraction kit (MoBio Laboratories, Inc.) according to the manufacturer's instructions.
A nested PCR was performed with an Eppendorf MasterCycler gradient (Eppendorf) to detect archaeal 16S rRNA genes. To test the effects of primer combinations, PCRs were carried out with three different primer sets on all samples (Table 1). All primer sets targeted sequences at the beginning (forward primers) and middle (reverse primers) of the 16S rRNA gene. All primers used in this study were designed by others to be specific to and target most archaea. PCRs were performed with Phusion proofreading PCR polymerase (Finnzymes) according to the manufacturer's instructions. The sizes and quality levels of the PCR products were checked with 1% agarose-ethidium bromide gels according to standard protocols (29).
|
View this table: [in a new window] |
TABLE 1. Primer combinations used in PCR
|
Sequencing.
Isolated DGGE fragments were reamplified for sequencing with the same nested primer pair used for the original amplifications (A109a/A934b or A109a/Ar9r) and sequenced at the Haartman Institute, University of Helsinki, with primers Ar514r (11) and Ar344f (28). The sequences obtained with these primers overlapped in the middle of the fragment, as they were produced from the middle toward the ends of the gene fragment (as in reference 22). This method also resulted in sequences of different lengths, as the signal of the sequencing reaction subsided variably. After an initial sequencing of one representative band from each OTU group, two or more comigrating bands from different samples of the same type, obtained with the same primer combination, were chosen for sequencing to test the reproducibility of the DGGE method. The sequences were edited with the program package Vector NTI (InforMax) and compared to all archaeal and bacterial 16S rRNA gene sequences in the DDBJ/EMBL/GenBank databases with BLASTn. Most sequences obtained were 450 to 700 bp. The sequences were tested for chimeras with the Chimera Check tool in Ribosomal Database Project II (http://rdp.cme.msu.edu/index.jsp). No chimerical sequences were found, but all unspecific and bacterial sequences were subtracted from the analysis.
Phylogenetic analyses.
The edited sequences were aligned using the ClustalW alignment tool (http://www.ebi.ac.uk/) with default settings, except that the Gonnet matrix was used. The alignments were checked and manually edited with the BioEdit alignment editor (version 7.0.5; Tom Hall). Maximum-likelihood and maximum-parsimony analyses were performed with PAUP 4.0b (Sinauer Associates, Sunderland, MA), using heuristic search parameters for 10 repeats of stepwise addition of random sequence with a tree bisection-reconnection branch-swapping algorithm. A maximum of 10,000 trees was set for the maximum-parsimony analyses. Consensus bootstrap trees were calculated for both maximum-parsimony and maximum-likelihood phylograms with 1,000 bootstrap replicates. For the initial phylogenetic analyses, all sequences were used. However, some of the sequences were too short to be used for reasonable phylogenetic analyses and were left out of the final analysis.
Statistical analyses.
Statistical significance for differences between frequencies of occurrence of archaeal OTUs in different sample types was determined by Kruskal-Wallis nonparametric one-way analysis of variance and post hoc analysis using Dunn's multiple comparison test (InStat package; GraphPad Software, Inc.). The tests were performed on the combined results for OTUs obtained by all three primer sets. Nonarchaeal and unidentified OTUs were omitted from the analysis.
Nucleotide sequence accession numbers.
The GenBank accession numbers for the new sequences are AM903346 to AM903374.
|
|
|---|
![]() View larger version (120K): [in a new window] |
FIG. 1. Representative lanes of DGGE runs of all different sample types. Samples in the upper sectors are from fine roots (Fr; nonmycorrhizal), and those in the lower sectors are from mature mycorrhizas (Mm; inoculated with P. involutus). The tree species is marked above each sector, with the sector furthest to the right representing uncolonized humus (Uh) samples. A molecular marker (M) is shown, and the primer combination numbers (Table 1) are marked above each pair of lanes. The archaeal bands are marked with a number corresponding to the OTU codes. b, bacteria; *, unidentified.
|
Comparison of average appearances of archaeal OTUs.
The number of archaeal OTUs was statistically significantly higher in mycorrhizal samples than in fine-root samples of P. abies and B. pendula (two-tailed P = 0.007 and 0.003, respectively) (Fig. 2). The result was the same for P. sylvestris, but statistical analyses were not carried out, due to a complete absence of archaeal OTUs in fine roots. The same trend could be seen in the A. glutinosa samples with primer combinations 1 and 2, although the result was not statistically significant (two-tailed P = 0.21). Primer combination 3 specifically amplified archaeal 16S rRNA gene sequences from A. glutinosa fine-root samples, resulting in a high average number of OTUs compared with the levels for other fine-root samples. Nevertheless, the difference in number of archaeal OTUs between A. glutinosa fine-root and mycorrhizal samples observed when all primers were combined showed no statistical significance (P = 0.45).
![]() View larger version (24K): [in a new window] |
FIG. 2. Average number of archaeal OTUs per sample type, obtained from five replicate samples. The fine-root (Fr; nonmycorrhizal) and mycorrhizal (Mm; inoculated with P. involutus) samples of the same tree are shown in pairs. The tree species of each sample type is given below the x axis. Error bars show standard errors of means. The occurrences of OTUs in fine-root and mycorrhizal samples of the same tree species are shown above the column pairs as statistically significant (**) or not significant (ns). P. sylvestris fine roots could not be tested for statistical significance. The patterns on the columns represent the primer combinations used (no pattern, combination 1; waves, combination 2; vertical stripes, combination 3).
|
Phylogenetic analyses.
The majority of the 16S rRNA gene sequences obtained belonged to Crenarchaeota (Fig. 3). In addition, Euryarchaeota (Halobacterium) were found in two mycorrhizal samples (Fig. 4). The reproducibility of the DGGE method was tested by sequencing several OTUs with the same mobility from different samples. We found that although the OTUs showed similar mobility patterns in the DGGE run, the sequences were often quite different. According to DGGEs with PCR products obtained with primer combination 1, many of the fine-root OTUs were bacterial and were omitted from the study. However, all corresponding counterparts in the mycorrhizal samples were archaeal. Only two of the sequenced mycorrhizal OTUs were bacterial. Thus, all DGGE bands which were included in the analyses were with reasonable certainty archaeal.
![]() View larger version (46K): [in a new window] |
FIG. 3. Maximum-likelihood phylogram presenting members of group 1.1c Crenarchaeota. The sequences found in this study are shown in bold script. The first two letters in the sample codes represent the tree species (Ps, P. sylvestris; Pa, P. abies; Bp, B. pendula; Ag, A. glutinosa). The middle two letters of the mycorrhizal sample codes represent the fungus (Pi, P. involutus). The last two letters represent the sample type (Fr, fine root; Mm, mycorrhiza). The numbers (1, 2, or 3) after the OTU code indicate the primer combination used to obtain the OTU. Values for nodes with over 50% bootstrap support are shown. The tree is rooted by Euryarchaeota.
|
![]() View larger version (58K): [in a new window] |
FIG. 4. Maximum-likelihood phylogram presenting the Halobacteriales. The sequences of this study are shown in bold script, and the sequence codes are as described in the legend to Fig. 3. Values for nodes with over 50% bootstrap support are shown. The tree is rooted by Crenarchaeota.
|
All the new crenarchaeotal 16S rRNA gene sequences fell into group 1.1c and clustered with Crenarchaeota previously found in boreal forest soil habitats (Fig. 3). The sequences were divided into three clusters, the FFSC cluster, the nonmycorrhizal fine-root cluster, and the P. involutus mycorrhiza cluster. The sequences from fine-root and mycorrhizal samples never fell into the same cluster. All sequences obtained from fine-root samples fell within the nonmycorrhizal fine-root cluster, with 16S rRNA gene sequences presently found only in Finnish forest soils. The A. glutinosa fine roots showed a unique community of group 1.1c Crenarchaeota. These sequences were not seen in any other samples, and they formed a separate cluster within the nonmycorrhizal fine-root cluster in the phylogenetic analyses (Fig. 3). All but one of the mycorrhizal sequences fell into the P. involutus mycorrhiza cluster. The P. involutus mycorrhiza cluster contained 16S rRNA gene sequences from different forest, field, and alpine soil and freshwater environments, covering temperate and boreal areas. BpPiMm54 was the only sequence associated with the FFSC cluster containing 1.1c Crenarchaeota from boreal and temperate forest soils.
The sequences of OTUs 5 (from P. sylvestris), 6 (P. sylvestris), 7 (P. sylvestris), 8 (P. sylvestris), 14 (A. glutinosa), 47 (B. pendula), 53 (B. pendula), and 55 (P. abies) from the mycorrhizal samples fell into the P. involutus mycorrhiza cluster, and the sequences of OTUs 2 (P. abies and B. pendula), 6 (B. pendula), and 71 (A. glutinosa) from the fine-root samples fell into the nonmycorrhizal fine-root cluster (Fig. 1). However, these sequences were too short for reasonable phylogenetic analysis. The 16S rRNA gene sequences for the two euryarchaeotal samples AgPiMm51 and AgPiMm52 fell with the Halobacterium, together with similar sequences from boreal forest P. sylvestris mycorrhizas.
|
|
|---|
In our study, only three archaeal OTUs were obtained from soil currently uncolonized by mycorrhizal fungi. The sampled soil was not under the direct influence of the tree seedling and was not in contact with ectomycorrhizal hyphae, which distribute plant-derived substrates outside the immediate rhizosphere. We suggest that a rich vegetation cover with subsequent mycorrhizospheres plays an important role in the colonization of soil by 1.1c Crenarchaeota. This also agrees with the results obtained by Nicol et al. (20), who showed that in glacier foreland soil, no 1.1c Crenarchaeota were found until the soil had been ice free for at least 135 years, had a stable vegetation cover, and had a well-developed organic top layer. Our theory is supported by the fact that recently uncovered soil is primarily colonized by weakly to moderately mycotrophic plant species but that the mature flora mostly represent highly mycotrophic plant species (27). These plants are colonized by mycorrhizal fungi, and selection pressure on rhizosphere bacteria has been shown to be induced by these fungi (1, 15, 33). These more organic soils are also likely to contain more saprophytic fungi than the young mineral soils. Currently, there is no information about the relationship between saprophytic fungi and archaea. We also presume that mycorrhizal fungi influenced the archaeal phylotypes retrieved by Oline et al. (22), who sampled 1-cm-diameter and 10-cm-deep cores of mature, vegetated soils, most of which were dominated by ectomycorrhizal plant species, such as the tree species Pinus ponderosa, Pinus contorta, and Pseudotsuga menziesii. Since root and hyphal materials were not completely excluded from the samples, it is likely that mycorrhizas or mycorrhizal hyphae ended up in the soil samples when the soil was sieved with a 2-mm mesh. Kemnitz et al. (13) found that in temperate deciduous forests, 1.1c Crenarchaeota are also dominated by heavily mycorrhizal plants. The crenarchaeotal 1.1c 16S rRNA gene sequences found in mycorrhizal samples in our study were similar to the sequences retrieved from soils in the above-mentioned studies (13, 22). Only Sliwinsky and Goodman (34) did not find group 1.1c Crenarchaeota in soils vegetated by commonly mycorrhizal plants. Wallander et al. (38) have calculated that boreal forest humus contains a substantial amount of ectomycorrhizal hyphae, between 125 and 200 kg ha–1, in the organic top layer, equaling 1.25 to 2 g cm–2. Levels for all ectomycorrhizal mycelia (including ectomycorrhizal mantle structures) were estimated to be as high as 700 to 900 kg ha–1. This shows that ectomycorrhizal fungi constitute a huge part of the biomass in the forest soil associated with trees.
There was a great difference in archaeal community between the rhizospheres of some of the tree species, even when they were grown in soil originating from the same location. This contrasts with Sliwinski and Goodman's results (34), where significant differences were detected not in the rhizospheres of plants from the same site but only between different geographical locations. Our trees were all grown in the same humus, obtained from a dry pine forest. The differences in archaeal rhizosphere communities would certainly be even greater if the tree species had been grown in soil from a forest dominated by the tree in question, but then we would have been in the same situation as Sliwinski and Goodman (34), where the influence of the plant could not be separated from that of the location. Nicol et al. (20) studied the rhizospheres of different plant species in Austrian alpine glacier front soil. Minor differences were observed in the community structures of archaea associated with the roots of different plant species, but only between plants from soils of different ages. These results could be explained by the mycorrhizosphere effect, since few natural root systems are devoid of mycorrhizal fungi. In our present study, we found that the influence of the tree on the archaeal populations was minimized by the impact of the symbiotic ectomycorrhizal fungus colonizing the tree roots. Mycorrhizal fine roots had more even frequencies of archaeal OTUs, and these sequences were more similar to each other than to the sequences from uncolonized fine roots.
The mycorrhizal roots of most tree species harbored more-diverse populations of archaea than uncolonized fine roots or humus. Simon et al. (32) estimated that the doubling time of group 1.1b Crenarchaeota inhabiting tomato roots was 8 days. If the group 1.1c Crenarchaeota detected in this study grow equally slowly, mycorrhizas may present a more constant environment for these organisms than the rapidly growing root tips. It is also possible that the archaea are simply outcompeted by more-opportunistic bacteria on the root surfaces. This could explain the relatively high number of bacterial 16S rRNA gene sequences obtained from the fine-root samples. The crenarchaeotal 16S rRNA gene sequences obtained from fine roots formed a distinct cluster, separate from the mycorrhizal clusters. Interestingly, their closest neighbors were crenarchaeotal 16S rRNA gene sequences thus far found only in Finnish dry pine forest soil habitats (soil, mycorrhizospheres, and tree roots). In contrast to the crenarchaeotal 16S rRNA gene sequences detected in the mycorrhizal samples, none of those detected in the fine-root samples were similar to 16S rRNA gene sequences from other locations or habitat types.
There was no statistically significant difference in the occurrence of archaeal OTUs in the mycorrhizal samples derived from the four different tree species, and the sequences of the archaeal 16S rRNA gene OTUs all fell within the same cluster. This confirms that the mycorrhizal fungus in the roots attracts its own archaeal community independent of the one colonizing the rhizosphere, which is in agreement with our earlier results (3) showing that the archaeal population in the P. sylvestris mycorrhizosphere varies with the species of ectomycorrhizal fungus colonizing the roots. Euryarchaeota were detected in only one P. abies and one A. glutinosa mycorrhizal sample. These findings confirm our earlier results showing that, unlike some other ectomycorrhizal fungi, P. involutus does not particularly attract Euryarchaeota (3).
With the exception of A. glutinosa fine roots, the numbers of archaeal OTUs were low, and the archaeal 16S rRNA gene fragments detected appeared to be similar in the fine-root samples of the different trees. This is slightly surprising, as the different tree species influence the soil in their root and mycorrhizal vicinities by different means (26). The pH of the soil in a coniferous forest is low, and the typical pH of Finnish dry pine forest soils is between 4 and 5. Conifers, especially spruce, may decrease the pH of the soil to as low as 3.5, whereas birch has been shown to increase the soil pH up to 5.7 and enhance the cycling of nutrients (16). Deciduous litter also contains more easily leachable carbon, whereas coniferous litter is acidic and contains recalcitrant material, such as waxes, phenolics, and lignin (9, 24). The litter was not a factor in this study, because our experiment stretched over only one growing season. However, the short time of this experiment showed that A. glutinosa seedlings promoted both different numbers and different types of Crenarchaeota on their fine roots.
The A. glutinosa fine-root samples displayed the highest overall number of different archaeal OTUs in this study, followed by the A. glutinosa mycorrhizal samples (Fig. 2). The A. glutinosa fine-root samples were also the only ones that presented PCR products with primer combination 3 which formed a separate cluster within the nonmycorrhizal fine-root cluster. A. glutinosa is known for its ability to attract nodule-forming, nitrogen-fixing bacteria. No nodules were discovered on the roots of our test trees, but the A. glutinosa rhizospheric environment and consequently its microbial flora are probably specialized and different from the other tree rhizospheric microbial communities. This is evident also in the residing archaeal population. With the exception of the A. glutinosa fine-root results, archaea did not appear to thrive on tree fine roots but, in agreement with our earlier studies, preferred the mycorrhizas (2, 3). The primer set used had a great effect on the detection of archaea in the various samples. This is a known problem in studying natural, uncultured populations where even slight adjustments of procedures may alter the results. The lack of redundancy between the OTUs produced with the different primer combinations was conspicuous. However, the outer primers of primer combination 1 targeted a broad spectrum of microbial 16S rRNA genes and appeared to select for the Crenarchaeota in the P. involutus mycorrhiza cluster (Fig. 3), while the more archaea-specific outer primers selected for Crenarchaeota affiliating with the nonmycorrhizal fine-root cluster. It is possible that new primer combinations may reveal yet undiscovered archaeal groups in fine root, mycorrhizal, and uncolonized humus habitats.
Our conclusion is that different tree species have different influences on soil archaeal communities. These differences are diminished, however, by the ectomycorrhizal fungus colonizing the roots. The ectomycorrhizal fungus P. involutus attracts group 1.1c Crenarchaeota, and the crenarchaeotal populations in the mycorrhizas are different from those in the uncolonized fine roots of boreal forest trees.
Graeme Nicol and Leone Montonen are warmly thanked for their critical comments on the manuscript. Uwe Münster and Jari Valkonen are gratefully acknowledged for their support during this work.
Published ahead of print on 31 October 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»