Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2005, p. 1996-2000, Vol. 71, No. 4
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.4.1996-2000.2005
NASA-Ames Research Center, Moffett Field, California,1 National Aeronautics and Space Administration Headquarters, Washington, D.C.,2 Department of Microbiology, Montana State University, Bozeman, Montana3
Received 8 April 2004/ Accepted 27 October 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
In this study, we focused on the effects of artificial defoliation on the soil fungal community. The defoliation of plants can either increase carbon allocation to roots, thereby providing increased resources for microbial growth (25, 29), or decrease carbon allocation to roots (22, 42). This relationship depends on whether the plant is grazing tolerant (defoliation results in increased carbon to shoots) or grazing intolerant (defoliation results in increased carbon to roots) (22) and on the age of the plant (12). These relationships can be complicated, however, by root death in response to defoliation, which results in the release of nutrients that can increase the carbon and nitrogen resources available to soil communities and can alter soil pH through the release of organic acids and amino acids (37).
Fungi comprise one of the most important functional groups of soil microbes and are critical to nutrient cycling in terrestrial ecosystems (28, 53). Fungi vary in their carbon demands, and the loss of photosynthetic potential through defoliation can decrease their ability to support and maintain fungal hyphal networks (2). Most studies of defoliation effects on fungal communities have focused on grasses, and the effects can be either positive, e.g., increased root colonization by arbuscular mycorrhizae (12), or negative, e.g., reduced ability to maintain hyphal networks (2). These differences in effects depend on the severity of defoliation, the soil fertility, and the life history strategy of both the fungal species (2, 12, 37) and the defoliated plant.
Few studies have addressed defoliation or carbon sink effects in the context of conifers, although soil fungal biomass is proportionate to photosynthesizing biomass (36). Defoliation decreases fruiting body formation in Betula stands (33), decreases ectomycorrhizal (EM) infection of Douglas fir (32) and Pinyon pine (10, 17), and can alter the EM morphotype diversity of Scots pine (49). Artificial defoliation of only pines in the mixed pine-spruce system studied here can significantly alter EM fungal-community structure (9) and significantly affect the EM mututalisms of neighboring, undefoliated spruce (9). Thus, the defoliation of only one tree species in a mixed-species setting can have relatively far-reaching effects (9).
Even if the defoliation effects on fungi are limited to EM interactions, EM fungi often form extensive hyphal networks in the soil that function as nutrient-gathering instruments and effectively increase the surface area across which trees can absorb nutrients (27, 47). In this study, we used a combination of molecular genetic and culture-based methods to determine the effects of defoliation of only pines in a mixed pine-spruce forest in Yellowstone National Park on both genetic diversity and functional diversity in the soil fungal community.
We tested the null hypothesis that despite the presence of a second undefoliated tree species supplying root exudates, the defoliation of pines in mixed pine-spruce forest would not significantly affect the soil fungal community. Our rationale was that the presence of undefoliated spruce would buffer the community against the disturbance and that the experiment would provide insight into the role of system complexity in maintaining diversity in the face of disturbance. We used denaturing gradient gel electrophoresis (DGGE) to evaluate changes in soil fungal communities in response to defoliation in terms of the presence or absence and richness of soil fungal species that could be amplified using PCR. We also used filamentous-fungus plates (Fungilog; Biolog, Hayward, Calif.) (11, 16, 23) to determine effects of defoliation on soil fungal physiological biodiversity. These methods provide different measures of fungal diversity, a measure of the fungal DNA resident in the soil, and a measure of the culturable members of the fungal community. DGGE measures diversity of the fungal DNA resident in the soil, while Fungilog growth assays measure the potential physiological capabilities of the fungi that are culturable. This study is the first to combine these methods to study effects of a conifer system on the soil fungal community and to take the possible moderating effects due to above-ground species diversity into account.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Each paired treatment and control block was 5 by 3 m, and each pair was separated from the surrounding forest by at least 5 m to eliminate edge effect. Treatment and control blocks were situated directly adjacent to each other along a 30-m transect, and the three replicate paired blocks were all within 5 m of each other. This spatial arrangement was arrived at empirically (8, 9), was designed to decrease the likelihood that spatial variation in the fungal community would overshadow treatment effects, and has been used by us in the past to assess effects of artificial defoliation on the EM community of this stand (9).
Trees of both P. contorta and Picea engelmannii species had reseeded naturally and were 1 to 2 m in height, while P. contorta stems significantly outnumbered Picea engelmannii stems (P < 0.001). Control blocks contained 60 ± 17 (standard error) P. contorta individuals and 7 ± 2 Picea engelmannii individuals, while treatment blocks contained 68 ± 16 and 10 ± 2 individuals, respectively. All individuals of P. contorta in treatment plots were 50% defoliated by twisting at the base of every other needle whorl over the entire canopy.
For DGGE analysis, seven 15-ml soil samples (approximately 2 cm in diameter and 12 cm deep; number of cores, 21) were taken from within each block, preserved on dry ice for 3 days (until they were returned to the lab), and then freeze-dried for long-term storage. Roots were concentrated in the top few centimeters of each core, and soils for analysis were taken from the root zone. DNA was extracted from 10-mg soil samples (containing no root material) by using sonication and a detergent-based mini-prep (7). Detergent-based preps of small amounts of soil can provide resolution down to 1 culturable cell per 10 mg of soil by PCR-based assays (6). To detect fungal DNA in the soil extract, the universal fungal primer set ITS1F-ITS4 was used (13). All reactions were successful. The parameters for PCR were as follows: initial denaturation at 95°C for 1 min 35 s; 13 cycles of denaturation at 94°C, primer annealing for 55 s at 55°C, and polymerization for 45 s at 72°C; 9 additional cycles with the polymerization time extended to 2 min; 9 cycles with a 3-min extension; and a final 10-min polymerization step at 72°C.
Amplified fungal DNA from soils were subjected to DGGE to separate fungal taxa and make identifications of individual fungal species possible. PCR amplification products were gel purified on agarose gels (2% low-melting temperature gels and 1% unmodified gels), the bands were excised with a sterile razor blade, the gel fragments were transferred to individual 1.5-ml Eppendorf tubes with 200 µl of Tris-EDTA buffer, and the tubes were heated to 65°C for 10 min to melt the agarose gel matrix and free the DNA. To ensure consistency within samples, the analysis of sample extracts (n = 14) from the first set of paired blocks was repeated three times with each sample. There was no variation among samples; thus, each of the remaining samples was assayed once.
For each DGGE assay, 12.5 µl of a 1:10 dilution of amplified DNA in Tris-EDTA (10 ml of 1 M Tris [pH 8.0] plus 4 ml of 0.2 5 M EDTA [pH 8.0] in 986 ml of distilled H2O) was used as a PCR template in a reaction to add GC clamp tails. The PCR protocol was the same as before except that one primer (ITS1F) had an additional 40-bp GC tail.
A DCode DGGE universal detection system apparatus (Bio-Rad, Hercules, Calif.) was used to perform the DGGE analysis according to the manufacturer's specifications and as previously described (13), except that a 20 to 60% denaturant gradient was used. Individual bands were cut out and sequenced on an automated sequencer (ABI Prism 377; Perkin-Elmer, Boston, Mass.) at California State University, Hayward. Sequence data were run through the BLAST search programs of the National Center for Biotechnology Information (www4.ncbi.nlm.nih.gov/) to determine the taxonomic affiliation of each DGGE band by sequence homology.
To determine the functional diversity of the soil fungal community, we used Fungilog microtiter plates. For each assay, a 10-g soil sample (three from each of the treatment and control blocks; a total of nine samples for each experiment) was shaken in 100 ml of Ringer's solution (20) (three times for 1 min each). This solution was allowed to settle for 15 min, and 150-µl aliquots were pipetted onto the Fungilog microtiter plates. Each plate contained 95 wells with different carbon sources and a single control well with no carbon source. Plates were incubated at room temperature with antibiotics (50 µg each of streptomycin and ampicillin/ml) to ensure that results were due to fungal and not bacterial growth. Every 24 h, for 5 days, plates were scanned with an optical density plate reader (at 490 nm). This method produces 96 data points (95 substrates plus one water control) per replicate, per day, with each corresponding to a specific carbon substrate. Growth is measured colorimetrically when the tetrazolium dye in the medium accepts electrons from active electron transport. Plates were read before growth became too massive, that is, before overrunning the plates and becoming nonlinear.
The significance of the effects of defoliation on the presence or absence and abundance of soil fungal species was determined with a contingency table and
2 analysis for the whole community and Fisher's exact test for the species that exhibited differences following treatment. Changes in species richness, defined as the number of unique DGGE bands (or fungal species) per unit area (in this case per soil sample), were analyzed by using Student's t test. Principal-components analysis (PCA) was used to project the data onto a small number of statistically independent axes that account for the majority of the variability in the data. This projection results from identification of the eigenvectors of the covariance matrix of the data. Then the data were projected onto the eigendirections with the largest eigenvalues. The significance of the effects on physiological diversity, as indicated by total number of substrates utilized and rates of substrate utilization, was determined with Student's t test, with each day being analyzed individually.
| RESULTS |
|---|
|
|
|---|
2 analysis of the total community indicated no significant effect on overall soil fungal-community composition; the majority of fungal species were present at the same frequency in both treatment blocks and control blocks (Table 1). Most BLAST analyses indicated that the DGGE bands were
95% similar (excluding unknown bases) to those of EM fungi of conifers. The two exceptions were taxon 10, which had high sequence similarity (
95%) with a species of the genus Eupenicillium, a free-living ascomycete, and taxon 8, which, though clearly a basidiomycetous fungus, was not closely related to any fungus represented by a sequence in GenBank (82% match to any fungus in GenBank and to only 300 bases of more than 600 sequenced). This sample was sequenced from multiple amplifications and in both directions in order to rule out the presence of chimeras.
|
The defoliation of only pines in a mixed pine-spruce forest resulted in significant changes in physiological capacity as measured by Fungilog tests (P < 0.001), both in terms of carbon source utilization (with growth being indicated by tetrazolium dye) (Fig. 1A) and in the number of carbon substrates utilized (Fig. 1B). Both measures were significantly higher in treatment samples from 72 h until the end of the experiment. Repeated measures with a further sampling of equal size produced similar, statistically significant results.
|
|
| DISCUSSION |
|---|
|
|
|---|
DGGE analysis indicated that defoliation of only pines in experimental blocks of a mixed stand of pine and spruce can affect the abundance of individual soil fungal species but that species richness (number of species per core) and overall soil fungal-community composition are not significantly affected. Most of the fungi detected by DGGE have high sequence similarity to EM fungi, with a few being similar to free-living fungal taxa. None of the fungi with high sequence similarity to free-living fungi were affected by defoliation.
One fungus that was affected by defoliation, Cenococcum, was more abundant in controls, and a second, taxon 8, was more abundant in treatments and thus appeared to respond positively to the changes induced by defoliation. Cenococcum is a well-studied EM fungus that is very common in pine forests, although little is known about how this fungus reacts to defoliation. Since Cenococcum can form a large biomass and extensive rhizomorphs in forest soil (14, 30), a reduction in photosynthetic potential and concomitant reduction in the fixed carbon allocated to the roots of the EM plants could reduce the amount of this fungus present.
We also hypothesize that the presence of nondefoliated spruce helped maintain fungal-species richness in these nutrient-poor soils and may also have contributed to the measured increase in Cenococcum organisms. Our data suggest that roots of the nondefoliated spruce tree can support the EM fungi of the defoliated pines (9) and hence help to maintain EM fungal-community composition and infection levels. In addition, growth of fungal endophytes can be influenced by the presence of specific host plant carbohydrates (24). Thus, Cenococcum may respond more readily to spruce root exudates and increase in frequency when roots of spruce proliferate following pine defoliation (9).
The Fungilog experiments indicated significant changes in the physiological capabilities of the fungal communities. PCA clustered 8 of 9 control cores into a distinct group, while most of the treatment cores were in a looser assemblage that was distinct from the control group. These groupings were based primarily on differences in carbohydrate utilization. The single control core that was separate from the main control group was taken near the edge of adjacent treatment and control blocks, and the hyphae could have crossed over into adjoining blocks. Three of the nine treatment cores also clustered with the majority of the control cores. We hypothesize that these cores contained hyphae responding to exudates from spruce, or communities that were still in transition, despite the passage of 2 years.
Fungilog results also indicated increased physiological diversity in response to defoliation in the form of increases in both the range of utilizable substrates and their rate of utilization. This change occurred despite the presence of the second undefoliated tree species. Thus, quantitative and qualitative differences in root exudates of the two tree species may play a role in influencing fungal-community structure and function.
However, Fungilog results differ from those of the DGGE analysis, which indicated very little change in species composition and richness. The discrepancy between the two methods may be due to two factors. First, if a single fungal species has a wide range of carbon substrate usage capabilities, the relatively large changes in physiological diversity may result even with only minor changes in fungal-community composition. Fungi can have very different niche sizes, defined as the number of carbon substrates each species or strain can utilize. Under some conditions, this can range from 18 of 95 Biolog substrates utilized to 92 of 95 Biolog substrates utilized (34). Studies of fungal endophytes of conifers are consistent with these observations and indicate that strains specific to petioles, for example, have a much broader substrate utilization capability than do those restricted to needles (4). Thus, the increase in taxon 8 organisms might suffice to produce the Fungilog patterns that we observed. We will test this hypothesis by assessing the physiological capabilities of the fungus that increased in frequency when or if the fruiting bodies needed for culturing this fungus become available.
A second explanation for the differences between Fungilog results and DGGE results could be PCR amplification bias. This bias occurs when primers are used in multiplex reactions (50), particularly if the fungi present in the soil lack sequences that correspond exactly to the primers used. Our primers have been tested extensively against fungi across the Eumycota (15), so we consider this possibility unlikely. However, even when perfect sequence homology exists, other factors may inhibit amplification. For example, the secondary structure of the DNA flanking the priming site can cause DNA from some species to be preferentially amplified (26), and thus some species might not be detected, even if they are present in relatively high numbers. We compared the sequences obtained in this study with those obtained in our previous study of EM in these study plots (9). None of the EM fungal species detected in our previous study were present in our DGGE analysis of the same soils. Thus, the EM fungi that we detected in our previous study form little or no hyphal biomass in soils, or there was extensive turnover in the EM fungal community. We plan to test the primer bias hypothesis by repeating the Fungilog assays and then sequencing and identifying the fungi that are cultured on each substrate.
In summary, our data indicate that the defoliation of only pines in a mixed pine-spruce study site can significantly affect the relative abundance of some soil fungal species, though overall species richness and composition are largely unaffected. Furthermore, the defoliation had a significant effect on the physiological diversity of the soil fungal community. However, these results indicate that care must be taken when interpreting the data and that independent assessments of genetic diversity within Fungilog cultures of the type that we plan to perform may be useful. Given the increasing importance placed on studies of biodiversity and physiological diversity, and roles of biodiversity and physiological diversity in ecosystem resilience (e.g., references 1, 5, 36, 38, 43, 44, 51, and 52), it is clear that studies utilizing methods that measure both species composition and species physiological diversity are needed.
| ACKNOWLEDGMENTS |
|---|
We thank J. R. Blair for mushroom identifications, Bob Lindstrom and Ann Deutch of the Yellowstone Center for Resources, Ann Rodman of the Yellowstone Soil Survey, Alan Kerner for logistical support, and Chris Baysdorfer of California State University, Hayward, for DNA sequencing.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|