Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2002, p. 326-334, Vol. 68, No. 1
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.1.326-334.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
David S. Treves,2 L.-Y. Wu,1,2 Terry L. Marsh,2 Robert V. ONeill,1 Anthony V. Palumbo,1 and James M. Tiedje2
Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831,1 Center for Microbial Ecology, Michigan State University, East Lansing, Michigan 488242
Received 11 June 2001/ Accepted 9 October 2001
|
|
|---|
|
|
|---|
To identify possible determinants that control microbial community structure in soil, we characterized microbial communities from 29 soil samples using a small-subunit (SSU) rRNA-based cloning approach. This allowed us to evaluate more than 9,000 restriction fragment length polymorphism (RFLP) profiles and to examine whether diversity profiles were correlated with particular environmental conditions. Our findings show that surface soils from low-carbon sites exhibited a uniform diversity pattern quite distinct from the pattern of saturated subsurface soils, which exhibited the more common competitive diversity pattern found in most biological communities. In contrast, microbial communities from high-carbon sites displayed the uniform diversity pattern regardless of depth or the water content of the sample. We hypothesized that spatial isolation can explain the community structure differences at the low-carbon sites but that high soil carbon contents can overcome such differences, as we observed at the high-carbon site.
|
|
|---|
The Delaware site is on Dover Air Force Base in an area where the vegetation is grass and cultivated shrubs (near well IR-07D; 39°7.374'N, 75°29.047'W). This site was contaminated with low levels of chloroethenes (17). The soil is primarily coarse sands with some lenses enriched in clay and is classified as a mesic Hapludult. The zones analyzed were all sands. The water table at the time of sampling was 5.5 m below the surface. The number of CFU, was 1.1 x 107 CFU g1 in the surface soil, but the concentration dropped sharply with depth and was <2.0 x 102 CFU g1 in samples obtained at depths of more than 6 m. The oxygen contents of water samples obtained at the Dover site are typically 2 to 3 mg/liter (7).
The Kellogg Biological Station surface samples were taken from noncultivated poplar plots in the long-term ecological research study site in southwest Michigan. The soil at this site is fine to coarse loamy sand and is classified as mesic Typic Hapludalf.
The Cannelton site is a wetland in northern Michigan on the shore of the St. Marys River. The soil is sandy and is classified as a frigid Entic Hapaquod. A tannery operated nearby from 1900 to 1958 and contaminated the area with high levels of chromium. The water table was 1.8, 1.1, and 1.4 m below the surface for the J, H, and O sample series, respectively. On the basis of color samples H17-2 and -3 showed some evidence of reducing conditions. The mean viable plate counts on R2A medium for Cannelton sites J, H, and O were 3.5 x 107 and 1.4 x 106 CFU g1 for the surface and saturated subsurface samples, respectively. Chromium concentrations in the Cannelton soils were determined with an inductively coupled plasma-mass spectrometer. Briefly, soils were extracted with nitric acid in pressurized vessels heated in a microwave (10). The extracted liquid was diluted, and metal concentrations were determined by inductively coupled plasma-mass spectrometry with a hexapole collision cell at the ICP-HEX-MS Laboratory at Michigan State University. Calibration standards for analyses were prepared by using distilled, deionized water and stock standards (J. T. Baker Analyzed). Standards and sample dilutions were prepared under identical solution conditions. All chemicals and reagents used were analytical metal grade or better.
Soil coring was done by using steam-cleaned Lexan core liners and a sonic drilling technique at the Virginia and Delaware sites. At the Cannelton site coring was done with a split-spoon sampler driven by hand, and greater depths were reached by hand auguring before sampling. Kellogg Biological Station samples were collected from surface soils only with a standard soil corer. Soil from the center of each core was removed with a sterile spatula, collected in sterile Whirl-pak bags, and either immediately frozen (Abbotts Pit and Dover samples) or dispensed into sterile plastic bottles, transported on ice, and then frozen within 12 to 48 h (Cannelton and Kellogg Biological Station samples).
DNA extraction.
DNA was extracted from 5 or 10 g of soil by grinding, freeze-thawing, and using sodium dodecyl sulfate for cell lysis (37). This method effectively lyses various types of bacteria, such as gram-positive bacteria and methanogens. The crude DNA was purified by the minicolumn purification method (37), except that the DNA was eluted twice from the resin column with 50 µl of hot (80°C) water.
SSU rRNA gene amplification and cloning.
The oligonucleotides used to amplify bacterial SSU rRNA genes were the fD1-rP1 primer set described by Weisburg et al. (35) with modifications to the linker sequences (39). All stock preparations used for PCR amplification were made and all procedures were performed with the precautions suggested by Kwok and Higuchi (14). One microliter of purified DNA (1 to 20 ng) was used as the template in a 20-µl reaction mixture. The conditions used to amplify SSU rRNA genes were those described previously (39).
The concentrations of PCR-amplified SSU rRNA gene products were estimated by comparing the band intensities on agarose gels to the band intensities for known concentrations of standard lambda DNA. The amplified PCR products were directly ligated to the PCR II vector obtained from Invitrogen (San Diego, Calif.). Ligation and transformation were carried out according to the manufacturers instructions. The ratios of inserts to vectors were 0.5:1 to 1:1. Two microliters of a ligation reaction mixture was transformed by the heat pulse method into 50 µl of Escherichia coli Top 10F' competent cells (Invitrogen). The transformed cells were then incubated in 450 µl of SOC medium at 37°C for 30 min, after which the entire SOC-cell mixture was plated onto Luria-Bertani agar plates with ampicillin.
SSU rDNA RFLP analysis.
All white colonies were picked and screened for SSU ribosomal DNA (rDNA) inserts, which were amplified (20 µl) with the primers specific for the polylinker of vector pCR II (38). One-third of the remaining 18 µl of PCR amplified products was digested with 0.1 U of MspI plus 0.1 U of RsaI or with 0.1 U of HhaI plus 0.1 U of HaeIII (Gibco BRL Life Technologies, Gaithersburg, Md.) overnight at 37°C. The resulting RFLP products were separated by gel electrophoresis in 3.5% Metaphor agarose (FMC Bioproducts, Rockland, Maine) in 1x Tris-borate-EDTA by using 7 V/cm for 4 h. The gel was stained with 0.5 µg of ethidium bromide per ml and visualized by UV excitation. RFLP patterns were compared by using the Molecular Analyst program (Bio-Rad, Hercules, Calif.). Clones with unique patterns were considered operational taxonomic units (OTUs), and their dominance patterns were evaluated as abundance profiles.
Diversity indices.
The reciprocal of Simpsons index (1/D) and the log series index were chosen to characterize the microbial communities in our soil samples because they have good to moderate discriminating ability and are used widely in ecological studies (18). The use of 1/D instead of the original formulation of Simpsons index ensures that an increase in the reciprocal index reflects an increase in diversity (18).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Site characteristics and OTUs for 29 soil samples from four geographically distinct sites
|
![]() View larger version (32K): [in a new window] |
FIG. 1. Representative microbial community diversity patterns based on OTU abundance in surface, vadose, and saturated soils from low- and high-carbon sites. The values in boxes are community diversity values based on the reciprocal of Simpsons index.
|
To provide a second estimate of community diversity, we applied the log series index (
) to our samples. This index is most sensitive to species richness (18) but has been recommended as a reasonable diversity measure even when species abundance does not follow the log series distribution (28). Visual inspection of the profiles suggested that samples with 1/D index values below about 50 and
values below 100 show typical dominance profiles (Fig. 2A). One low-carbon vadose zone sample (DB-1) showed a dominance profile, while the other two vadose zone samples did not, although their index values were lower than those of the surface soils (Fig. 2A).
![]() View larger version (17K): [in a new window] |
FIG. 2. Plots of the reciprocal of Simpsons index (1/D) versus the log series index for the microbial communities at different soil depths from low-carbon (A) and high-carbon (B) sites. The values in the box are index values that indicate dominance in community structure.
|
Microbial community structure at high-carbon sites.
In contrast to the communities at the low-carbon sites, the microbial communities at the high-carbon sites showed a uniform diversity pattern, regardless of depth (Fig. 1D to F). Not surprisingly, the 1/D values for these samples were similar to the low-carbon surface sample values and did not decrease with increasing depth of the sample (Fig. 1). Also, a comparison of the reciprocal of Simpsons index versus
showed that no high-carbon samples approached index values reflecting dominance profiles (Fig. 2B).
Besides having high total organic carbon contents, the Cannelton samples differed markedly in their chromium(III) contents, which was the dominant form of chromium at this site. This provided an opportunity to examine the effects of this contaminant on microbial community structure and diversity. The two samples with low chromium levels (M24-1 and O22-4) showed greater diversity than the samples with higher chromium levels based on the reciprocal of Simpsons index (Fig. 3). Beyond this, however, there was no clear relationship between microbial community diversity and chromium content. It is noteworthy that soils with very high chromium concentrations (20 to 30% Cr) and with severely stunted vegetation still had a reasonably high microbial diversity index.
![]() View larger version (15K): [in a new window] |
FIG. 3. Comparison of microbial community diversity (based on Simpsons index) and chromium levels at the high-carbon Cannelton site.
|
5 ng). Hence, the difference in patterns should be attributable to differences in rDNA diversity in the templates. To further evaluate this possibility, we estimated the number of chimeras by selecting, at random, 30 to 40 clones from each of the Abbotts Pit and Dover libraries and sequenced both the 5' and 3' ends. We then constructed independent phylogenetic trees and compared the branching orders to identify discrepancies which were indicative of different sources for the 3' and 5' sequences. The proportions of chimeric sequences were 5.2 and 7.4% for Abbotts Pit saturated subsurface and surface soil samples respectively, and 5.5 and 10.5% for Dover saturated subsurface and surface soil samples, respectively. These results suggest that the chimeric structure, at least for dissimilar genotypes, was not a severe problem in our results. Probably the second most significant potential artifact that could increase apparent diversity is the cloning of heteroduplexes. We found that the proportions of heteroduplexes were 3.6 and 4.5% for Abbotts Pit and Dover samples, respectively. Other potential errors that could artificially increase diversity, such as rrn heterogeneity (8) in the same organism and PCR fidelity, have been shown to occur less frequently (22, 23). In summary, while we believe that the libraries used contained some artificial OTUs, we think that the extents of the errors should not exceed 10 to 20% and that they should be similar in libraries yielding the two types of profiles. Since the differences in the indices for the dominance and uniform profiles are much larger (i.e., 2 to 3 orders of magnitude), we believe that the distinction between the two patterns is real and reflects a difference in underlying community structure. Indeed, the median value for the reciprocal of Simpsons index for the low-carbon surface samples was 5,490, compared to a value of just 15 for the low-carbon saturated samples.
How different are different OTUs?
We also used the 5' and 3' end sequence information to evaluate the differences between OTUs that varied by only one or two bands when they were cut with four tetrameric restriction enzymes. This should be the minimum level of phylotype difference detected by this method. OTUs that differed by one or two bands had SSU rRNA gene sequence similarities of 97.0% ± 2.4% at the 5' end and 97.6% ± 2.1% at the 3' end. Differences in SSU rDNA sequences of 2.5% or more correspond to <70% DNA-DNA hybridization (25), a key criterion indicating different bacterial species (34). This interpretation is supported by in silico digestion of SSU rDNA database sequences (20). Hence, in this study a unique OTU suggested a different species for nonartifactual clones. A more detailed phylogenetic analysis of clones from Abbotts Pit and Dover surface soils also supported this interpretation (J. Zhou, H. Huang, B. Xia, D. S. Treves, L. J. Hauser, R. J. Mural, A. V. Palumbo, and J. M. Tiedje, unpublished data).
Mechanisms producing the noncompetitive diversity pattern in low-carbon soils.
We propose that the dominance structure seen in the low-carbon saturated subsurface samples is a common result of competitive interactions. Some species are eliminated by competitive exclusion, and a few species, better adapted to the conditions, dominate. The lack of dominance observed in the surface soils is highly unusual in ecological communities. This extreme pattern implies that the microbial community must be experiencing an almost complete lack of competition. Many genotypes are maintained, and yet no genotype appears to be superior.
Since we hypothesize that the dominant diversity pattern observed in the low-carbon saturated subsurface soils is shaped by competitive interactions within the communities, we refer to such a diversity pattern as a competitive diversity pattern. This pattern has also been observed in the microbial communities in aquatic environments, such as sea vent biofilms (19), and is very common in plant and animal communities. In contrast, we refer to the uniform diversity pattern observed in the surface soils from the low-carbon sites and from all of the high-carbon sites as a noncompetitive diversity pattern. This pattern has also been observed in other surface soil samples, such as arctic tundra soil (38), Hawaiian soils (21), and Wisconsin agricultural soil (3).
The following four mechanisms may produce a noncompetitive diversity pattern: (i) superabundant resources (if more than enough resources are available, there is no competition and communities evolve toward high diversity), (ii) resource heterogeneity (if resources are available in many different forms, populations can avoid competition by specialization), (iii) spatial isolation (if the habitat is subdivided into many separate pockets of resources, populations can avoid competition by physical isolation), and (iv) nonequilibrium conditions (more populations can be maintained under fluctuating environmental conditions).
We examined the first three of these potential mechanisms with a simple model that includes the resource uptake by, transfer between, and loss by the microbes (Mi) adhering to a single soil particle (i) and the nutrient resource (Ri) available in the adjacent soil water (Fig. 4). A simple representation of each particle system is given by:
![]() |
![]() |
![]() View larger version (12K): [in a new window] |
FIG. 4. Conceptual model of two soil particles showing the relationships between microbial biomass (M) and nutrient (R) exchanges. Although a small sample of soil would contain a large number of similar particle systems, this minimal representation suffices for our purposes. cR1 and dR2, nutrient fluxes between soil particles through water film, which also permits migration of microorganisms between particles (albeit slow); I, nutrient input from rainfall; eR1 and eR2 nutrient fluxes that become unavailable for microbial uptake, either because of gravitational loss or because of immobilization in compounds that the organisms cannot metabolize.
|
A major difference between the surface and the saturated subsurface is water content, suggesting that free water that interconnects soil particles may be a key factor in determining diversity patterns. Genotypes with higher growth rates could produce larger populations since the limiting resource remains connected through diffusion. The bacteria reduce resources (Ri) by utilization, but diffusion from neighboring sites replaces this loss and bacteria (Mi) continue to grow. The populations that extract nutrients most efficiently dominate. With a continuous water film, to achieve the noncompetitive pattern, all uptake and loss coefficients (i.e., ai and bi) must be equal (data not shown). However, if we make a more realistic assumption and change the relative value of the uptake rate (ai), the population levels change in direct proportion (Fig. 5A). Hence, the populations cannot remain the same, and the noncompetitive pattern cannot be achieved with unrestricted solute diffusion.
![]() View larger version (19K): [in a new window] |
FIG. 5. Effects of spatial isolation on population dynamics. (A) No spatial isolation with unequal uptake and loss coefficients. When I = 5, a1 = 0.2, a2 = b1 = b2 = d = e = 0.1, c = 0.2, K = 10, and R1 = R2 = 25, M1 approaches 1.428 while M2 remains 0.714 at equilibrium. Thus, two populations cannot be equally abundant, and hence noncompetitive patterns cannot be achieved. (B) Spatial isolation. The results are similar to those described above, except that when I = c = d = 0.0, M1 = 0.129 and M2 = 0.112 after 100 time intervals; that is, initially, M1 has a slightly faster growth rate. However, M1 = M2 = 0.021 after 500 time intervals. The initial growth advantage of M1 is temporary since it uses up its isolated supply of nutrients more rapidly. Thus, the sizes of the isolated populations become nearly identical, as seen in the noncompetitive pattern, even if the populations differ in terms of their growth dynamics.
|
Mechanisms producing the noncompetitive diversity pattern in high-carbon soils.
Since the saturated subsurface samples from the high-carbon Cannelton site show a noncompetitive diversity profile, it is difficult to invoke the spatial isolation hypothesis as the sole explanation for the observations. Instead, it appears that the high carbon content and accompanying resource heterogeneity are a force driving microbial community structure at this site. The soils at this site are marshy, organic soils with very high carbon contents (up to 10 to 20% organic C). Thus, the Cannelton soil environment differs dramatically in terms of resource diversity and amount from the low-carbon subsurface soils, which contain less than 0.1% organic carbon.
It is thought that fluctuations in environmental conditions, which can also lead to more evenness in populations, are not a prominent factor in soils, especially below the surface. Even temperature, which does vary annually at the surface, does not vary at the depths of the deeper samples obtained in this study. While environmental variation cannot be ruled out as a factor influencing soil microbial diversity, it does not seem to be as significant in soil as the other two factors identified in this study.
Implications of the causes of high soil microbial diversity.
There are several important implications of the noncompetitive profile and the spatial isolation hypothesis. First, spatial isolation helps explain the incredibly high microbial diversity seen in soil on a small scale (e.g., 1 to 5 g). The high diversity seen in this study is consistent with the high diversity found in DNA reannealing studies which revealed at least 5 x 103 genotypes g of soil1 (30). Second, spatial isolation would enhance the probability of successful colonization by any alien or introduced strain but make it difficult for that strain or any indigenous strain to achieve dominance. This has important practical implications since it would be difficult for any invader to achieve dominance unless spatial isolation can be eliminated. This might be important in bioremediation of hazardous waste sites, in sequestration of C and N, in biocontrol of root diseases, and in enhancing N2 fixation by some forest or crop plants. Third, spatial isolation would favor development of microbial genotypes endemic to an area (5), which would further enhance microbial diversity on a global scale.
While spatial isolation may explain the diversity patterns that we observed at the low-carbon sites, our data from the Cannelton site suggests that large amounts of heterogeneous organic carbon also play a significant role in structuring microbial communities. If this is true, then an invading species can also become established in this environment easily; however, whereas we predict that high degrees of spatial isolation could be destroyed by the addition of water, this treatment would have little effect on a community with such rich and diverse carbon resources. The absence of a reduction in diversity in the soils with extremely high chromium contents is surprising and suggests that this carbon-rich environment is able to counteract the expected decrease in microbial diversity.
While our study supports the conclusion that very high microbial diversity is maintained in soil and the conclusion that there is an apparent lack of competition in low-carbon surface soil communities and high-carbon soil communities, further testing of the spatial isolation and resource heterogeneity hypotheses is needed. Understanding the forces that shape soil microbial communities should help in practical management of soil communities and in the wise use and maintenance of microbial diversity.
This research was funded by the Natural and Accelerated Bioremediation Research program, Biological and Environmental Research, U.S. Department of Energy. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the Department of Energy under contract DE-AC05-96OR22464.
Present address: Institute of Environmental Science, Zhongshan University, Guangzhou, China, 510275. ![]()
|
|
|---|
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»