Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2003, p. 1800-1809, Vol. 69, No. 3
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.3.1800-1809.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Biological Sciences, University of Essex, Wivenhoe Park, Colchester CO4 3SQ, England, United Kingdom
Received 5 August 2002/ Accepted 25 November 2002
|
|
|---|
|
|
|---|
Management systems, such as tillage regimens, have also been investigated. Curci et al. (11) discovered that enzyme activity was higher in the uppermost 20 cm of soil in plots tilled by shallow ploughing and scarification than in those tilled by deep ploughing. Boddington and Dodd (2) investigated the effect of soil disturbance on the spore density, species richness, and extraradical mycelium lengths of arbuscular mycorrhizal fungi. All factors were reduced in disturbed soil in comparison with undisturbed soil, further indicating that microbial communities are influenced by management practices. Although organic farms have been shown to have slightly higher levels of organic matter and carbon than neighboring conventional farms (14, 37, 53), only limited research has investigated the structures and compositions of microbial communities following a switch to organic farming practices.
Using nucleic acid methodology, clear differences have been observed between communities from different management regimens. A ranking of the complexity of community DNA regimes from grasslands by DNA-DNA hybridization resulted in different complexities, in the order unimproved > semiimproved > improved grassland soil communities (9, 10). The overall levels of management and inorganic nutrient amendment designated the status of agricultural grassland pasture sites, improved being the most managed and unimproved being the least managed. This relationship held true for soils from three separate sites (9, 10). McCaig et al. (39) separated improved and unimproved grassland soils on the basis of their 16S ribosomal DNA (rDNA) clone diversity, with improved soil being less diverse (concurring with the findings of Clegg et al. [10] for the same sites).
Significant changes in soil microbial structure with vegetation change have been detected by use of molecular methods. Nusslein and Tiedje (45) showed that a change in the vegetative cover of a Hawaiian soil from forest to pasture led to a significant change in the composition of the soil bacterial community upon analysis of G+C distributions and corresponding clone group abundances. Peak G+C contents of the forest soil differed from those of the pasture soil, and none of the dominant phylotypes found in the forest soil were detected in the pasture soil. However, these Hawaiian soils are young and relatively geographically isolated (44) and therefore may be more susceptible to change than well-established soils. Indeed, Buckley and Schmidt (5), using 16S rRNA-targeted oligonucleotide probes to quantify the abundance of rRNA from major phylogenetic lineages, observed remarkably similar microbial community structures in soils that shared similar long-term histories of agricultural management practices despite differences in above-ground community compositions and different recent land-use management practices.
In the present research, we investigated soil microbial communities at a number of geographically distant agricultural sites (separated by over 65 km) comprising two distinct soil types at farms located in eastern England. By analyzing the physical, chemical, and biological characteristics of the soils, we sought to identify causal relationships between land-use management practices and soil health. A polyphasic approach was used. In addition to characterizing the soils by traditional analyses, the total and active members of the bacterial communities were described by molecular analysis. Finally, to incorporate metabolic activity, the physiological profiles of the bacterial communities were determined.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Field site information July 2000
|
Soil, physical, chemical, and biological analyses.
Labile phosphates, nitrates, and sulfates in the soil were measured by ion chromatography with a Dionex Corporation (Sunnyvale, Calif.) chromatograph. The organic matter concentration was measured by loss on ignition with a muffle furnace (1), and total C and N concentrations were measured with a model 2400 CHN analyzer (Perkin-Elmer, Buckinghamshire, United Kingdom). Water-holding capacity was determined by using the 0-bar method of saturation followed by free draining (7). Total soil bacterial counts were measured by using acridine orange staining (26), and viable counts were determined after plating on 0.1x tryptone soy agar (55).
DNA and RNA extraction.
DNA and RNA were coextracted (as described by Steffan et al. [59]) from 0.5 g of soil in the presence of 0.5 ml of phosphate buffer (100 mM, pH 7.5) and 0.5 ml of Tris-equilibrated phenol. Mechanical lysis was performed by using 0.5 g of acid-washed, UV-treated glass beads (150 to 212 µm in diameter; Sigma/Aldrich). The soil-bead mixture was vortexed for 30 s by using a Whirlimixer (Fisons Scientific) and placed on ice to prevent enzymatic degradation. The mixture was vortexed two more times for 30 s each time, stored on ice, and centrifuged (Eppendorf) for 10 min at 12,000 x g to separate the aqueous layer from the soil-phenol layer. Sequential phenol-chloroform purification was performed by adding 1 volume of phenol and 1 volume of chloroform-isoamyl alcohol (24:1) to the DNA solution, mixing the solution by shaking for 1 min, and centrifuging the mixture for 5 min at 12,000 x g. It was at this stage that crude DNA was quantified by agarose gel eletrophoresis to enable equivalent concentrations of DNA to be added as templates to the subsequent PCR. DNA mass ladders (Bioline, London, United Kingdom) were run alongside the environmental DNA. Following analysis with Phoretix 1D Advanced Analysis software (Nonlinear Dynamics Ltd., Newcastle, United Kingdom), calibration was produced from the known standards and the concentrations of the unknowns were calculated from the calibration. Further purification of DNA and RNA for PCR was achieved after agarose gel electrophoresis as follows: the DNA and RNA bands were excised from the gels and purified by using Hybaid recovery DNA purification kit II and a Bio 101 RNaid kit (BIO101, Vista, Calif.), respectively, in accordance with the manufacturers' instructions. DNA and RNA were eluted in nuclease-free sterile water and stored at -20°C.
DGGE analysis.
PCR of RNA and DNA was performed with 16S rDNA universal bacterial denaturing gradient gel electrophoresis (DGGE) primers (synthesized by Invitrogen Custom Primers, Paisley, United Kingdom) 2 and 3 (43) throughout to amplify the V3 hypervariable region of 16S rRNA genes (65). Additionally, these primers exclude regions of secondary structure which are located downstream and which are known to inhibit reverse transcription (RT) (67); these primers are thus also suitable for use in RT-PCR amplifications. Prior to RT-PCR, RNA was treated with DNase RQ1 (Promega, Southampton, United Kingdom) in accordance with the manufacturer's instructions. RT then was performed, and a cDNA copy was synthesized from the rRNA by using Superscript II (Gibco BRL, Paisley, United Kingdom) before PCR with primers 2 and 3. To avoid the possibility of contaminating DNA in the RNA extracts used for RT-PCR, control PCRs containing DNase-treated nucleic acid extracts that were not previously subjected to RT were performed. PCR mixtures for the above-described primer set consisted of 5 µl of PCR buffer (Bioline), 2 mM MgCl2, 250 µM each deoxynucleoside triphosphate, 400 nM each primer, 1 U of Biopro DNA polymerase (Bioline), 1 µl of template DNA, and sterile nuclease-free water to a final volume of 50 µl. PCR conditions were 5 min at 95°C; 30 cycles of 30 s at 94°C, 30 s at 55°C, and 1 min at 72°C; and a final extension for 10 min at 72°C.
The PCR products were examined by DGGE. The gels were formed as stated by Muyzer et al. (43), with the exception that an acrylamide gradient of 8 to 12% was used in addition to the 40 to 60% urea-formamide gradient. The gels were poured by using the gradient delivery system supplied with the Bio-Rad DCode system and were run for 5 h at 180 V and a constant temperature of 58°C. The gels were fixed overnight with 10% (vol/vol) ethanol- 0.5% (vol/vol) glacial acetic acid. Staining was performed with a 0.1% (wt/vol) silver nitrate solution, and developing was performed with a solution containing 0.01% (wt/vol) sodium borohydride and 0.4% (vol/vol) formaldehyde. The gels were fixed with 0.75% (wt/vol) sodium carbonate and preserved with a solution containing 25% (vol/vol) ethanol and 10% (vol/vol) glycerol. The gels were scanned by using an Epson GT9600 scanner and analyzed by using the Phoretix 1D Advanced Analysis package.
Terminal restriction fragment length polymorphism (T-RFLP) analysis.
PCR was performed with the fluorescence-labeled oligonucleotides 63F (5'-CAGGCCTAACACATGCAAGTC-3') (38, 48) and 518R (5'-CGTATTACCGCGGCTGCTCG-3') (34); these were labeled at the 5' end with the phosphoramidite dyes FAM and HEX (Applied Biosystems, Foster City, Calif.), respectively. These primers generate an amplicon including the V1, V2, and V3 hypervariable regions (65). PCR mixtures consisted of the ingredients described above. PCR conditions were 2 min at 94°C; 30 cycles of 1 min at 94°C, 1 min at 55°C, and 2 min at 72°C; and a final extension for 10 min at 72°C. Following the removal of unincorporated nucleotides with QIAquick columns (Qiagen), the products were digested at 37°C for 3 h with restriction enzyme AluI or CfoI. These enzymes, together with primer 63F, were previously demonstrated to generate larger numbers of terminal restriction fragments (T-RFs) than other enzyme combinations; primer 518R provides greater discrimination than primer 1389R, which was previously shown to be less discriminatory than primers located in the 5' end of the gene (48). Restriction digests were then mixed with 2 µl of deionized formamide and 0.5 µl of a ROX-labeled Genescan 500-bp internal size standard (Applied Biosystems), denatured for 5 min by boiling, and immediately transferred to ice. Triplicate samples were loaded onto an ABI310 automated genetic analyzer (Applied Biosystems) for capillary electrophoresis, and the resulting profiles were analyzed with Genescan software. Fragments with a peak height of less than 50 were excluded from subsequent analyses.
BIOLOG analysis.
A 100-ml soil suspension was generated from an original 10-g aliquot (dry-weight equivalent) in 0.25x Ringer's solution (BDH, Poole, United Kingdom) (24), thus achieving a 10-1 dilution. This solution was vigorously shaken for 10 min on an armed shaker (Stuart Scientific) to dislodge bacteria before dilution to 10-2. Low-speed centrifugation (1,500 x g for 10 min) was used to remove soil particles. A 150-µl aliquot of the bacterial suspension was inoculated into each well of a BIOLOG, Inc. (Hayward, Calif.), Eco-Plate. The plate was read with a Dias microplate reader (Dynex Technologies) over a 7-day period at a wavelength of 600 nm. Readings at day 0 were subtracted from subsequent readings to eliminate background color generated from the substrates and the bacterial suspension. In addition, readings generated from the control wells (no substrate; tetrazolium dye only) were also subtracted. The BIOLOG data were analyzed as described by Garland (22). The average well color development (AWCD) value was calculated for each sample at each time point by dividing the sum of the optical density data by 31 (number of substrates). Further analysis was performed on sample data when the AWCD value was approximately 1. The data were normalized by dividing the absorbance values by the AWCD values to reduce bias between samples with different inoculum densities.
General statistical analysis.
The Shannon index (H') (56) (H' = -
pi x ln pi) was calculated from the bacterial community profiles, where pi is the proportion of members that a particular species contributes to the total in the sample. This value was obtained by normalization of the volume data derived from the DGGE gels by using the Phoretix 1D Advanced Analysis package for each phylotype (band) within each community and representation as a proportion of the total volume data for that community. This software eliminates the background and automatically detects peaks when noise levels and minimum peak thresholds are set. Equitability (J) (50) (J = H'/H'max) was then calculated for a bacterial community, where H'max is the possible maximum diversity that could be obtained from that community, thus giving a measure of community heterogeneity unbiased by sample size. This principle has been used in other studies of microbial ecology (30, 40). UPGMA (unweighted pair-group method with arithmetic mean) dendrograms were constructed by using the Dice-Sorensen similarity index (50a) with either Phoretix or MultiVariate Statistical Package (MVSP) version 3.12h (GeoMem, Blairgowrie, United Kingdom) (see Fig. 3). Indices were similarly generated from the T-RFLP and substrate utilization data. Statistical significance was determined by analysis of triplicates with analysis of variance and Tukey tests (Microsoft Excel version 7a); a P value of <0.05 was required to establish a significant difference between samples. Alternatively, when duplicate samples were analyzed, all data were cited. Principal-component analysis (PCA) was performed by using Minitab for Windows version 13.20.
![]() View larger version (15K): [in a new window] |
FIG. 3. UPGMA dendrogram constructed from the similarity matching data (Dice-Sorensen index) produced from the DGGE profiles of 16S rDNA amplified from soil and generated by using MVSP version 3.12h. The scale bar represents percent similarity. Duplicate samples were analyzed.
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. Physical, chemical, and biological characteristics of soil from sample sitesa
|
![]() View larger version (16K): [in a new window] |
FIG. 1. Total carbon content (obtained with a CHN analyzer) plotted against organic matter content (determined by loss on ignition [LOI]) (1). Symbols: , EC2; , EC3; +, RMA; , RMB; , RMC; , RCA; , RCC; and x, RCD. Groupings were supported by ordering of principal components derived from T-RFLP profiles (see the text and Fig. 7).
|
![]() View larger version (16K): [in a new window] |
FIG. 2. PCA ordination of data (axes 1 and 2) generated from the substrate utilization profiles produced from the optical densities measured from the BIOLOG Eco-Plate. Symbols are as defined in the legend to Fig. 1. Groupings were supported by soil physicochemical data.
|
|
View this table: [in a new window] |
TABLE 3. Measures of heterogeneity for the substrate utilization and nucleic acid profiles for the soil communities at the eight sample sitesa
|
![]() View larger version (10K): [in a new window] |
FIG. 4. PCA ordering of data (axes 1 and 2) generated from the DGGE profiles of amplified bacterial 16S rDNA. Symbols are as defined in the legend to Fig. 1. Groupings were supported by cluster analysis of DGGE profiles (Fig. 3).
|
![]() View larger version (15K): [in a new window] |
FIG. 5. UPGMA dendrogram constructed from the similarity matching data (Dice-Sorensen index) produced from the DGGE profiles of 16S rRNA amplified from soil and generated by using MVSP version 3.12h. The scale bar represents percent similarity. Duplicate samples were analyzed.
|
![]() View larger version (18K): [in a new window] |
FIG. 6. PCA ordering of data (axes 1 and 2) generated from all T-RFLP fragments (5' and 3' T-RF frequencies generated after AluI and CfoI digestion) of 16S rDNA amplified from soil. The mean and standard error number of T-RFs analyzed per sample was 172.4 ± 4.0; the mean and standard error number of T-RFs per primer-enzyme combination was 43.1 ± 1.2. The eight soil samples analyzed in triplicate yielded 24 T-RFLPs. Symbols are as defined in the legend to Fig. 1.
|
![]() View larger version (21K): [in a new window] |
FIG. 7. UPGMA dendrogram constructed from similarity matching data (Dice-Sorensen index) for all T-RFLP fragments (5' and 3' T-RF frequencies generated after AluI and CfoI digestion) of 16S rDNA amplified from soil and generated by using MVSP version 3.12h. The mean and standard deviation number of T-RFs analyzed was as stated in the legend to Fig. 6. The scale bar represents percent similarity.
|
|
|
|---|
Land-use management has been shown to have an effect on both the physical and the chemical compositions of the soil and on the structure of the microbial community (2, 14, 37, 39, 53). In the present study, however, the organically managed fields did not show significantly greater levels of organic matter than the fields under GAP management. Similarly, elevated bacterial numbers were not found in organically managed fields. Fields with high levels of organic matter did, however, show elevated nitrate concentrations.
The substrate utilization data separated the EC fields from all other fields as the most dissimilar, as did profiles generated from the total (rDNA) and active (rRNA) bacterial communities by DGGE and by T-RFLP analysis of products amplified from rDNA. The EC soils were distinguished from the other soils on the basis of soil type and physicochemical properties, in particular, by a significant increase in the inorganic carbon content due to the chalky nature of the soils and by a water-holding capacity significantly lower than those of all other soils in the study.
Despite variations in many chemical and biological soil characteristics, five fields of the typical stagnogley soil type on two geographically distinct farms (RMB, RMC, RCA, RCC, and RCD) exhibited almost identical rDNA-derived DGGE profiles. Likewise, these soil communities formed a similar cluster after analysis of rDNA T-RFLP profiles. Therefore, it is likely that the total bacterial community compositions have been determined primarily by the underlying soil chemistry and structure rather than by the different management practices or cropping regimens at these sites.
Spatial homogeneity in soil over distances of several hundred meters was also revealed in Dutch meadows by 16S rRNA and rDNA fingerprinting (temperature gradient gel electrophoresis [TGGE]) by Felske and Akkermans (17). Almost identical bacterial profiles were generated between samples, suggesting that despite heterogeneity in soil properties, the dominant members of the microbial community appeared to be ubiquitous in this area (as determined with universal bacterial primers). Buckley and Schmidt (5) also highlighted the fact that in the published literature on the average abundance of 16S rDNA clones from the alpha, beta, and gamma Proteobacteria and Actinobacteria, originating from diverse soils from three different continents, samples showed remarkable similarities (3, 4, 33, 35, 36, 39). This finding suggests that certain characteristics of soil environments can lead to overall similarities. Spatial and temporal variations in soil chemical, physical, and biological properties, such as water content, nutrient and mineral availability, pH, and soil structure, contribute to immense heterogeneity and complexity of niches available within the soil environment (41, 64). This variability results in the formation of many discrete microhabitats, thereby conferring suitable habitats on a much smaller scale; habitat specificity is therefore not prohibitive. In addition, microorganisms have the ability to persist despite unfavorable conditions for growth, thereby permitting almost a universal presence of species, provided that initial dispersal was permitted (21). These factors may explain the lack of differentiation in the total bacterial community structures on two geographically distinct farms (>65 km apart), observed in the present study.
The T-RFLP data derived from DNA were similar to the data obtained from DGGE, forming three discrete clusters (Fig. 3 and 7). T-RFLP analysis provided slightly greater discrimination than DGGE, probably due to the inclusion of the V1, V2, and V3 hypervariable regions (65), whereas primers 2 and 3 covered only the V3 region. This factor, in addition to the combination of results from two restriction enzymes and both forward and reverse primers, may also explain the greater variability between replicates seen with T-RFLP analysis than with DGGE. However, the PCA ordination of the T-RFs produced slightly different groupings than the cluster dendrograms constructed from similarity matching data (Fig. 6 and 7). Samples RMA, RCC, and RCD were grouped together by PCA due to the more even distribution of phylotypes in their communities than in those of other samples, as determined by the equitability index (Table 3). These groupings were similar to those produced by the plotting of total carbon content against loss-on-ignition data (Fig. 1), perhaps suggesting a relationship between organic and inorganic carbon ratios and distribution within bacterial communities. The fact that ordination takes into account not only the presence and absence of species but also the distribution in the community highlights the need for comprehensive and combinational analysis methods. When Buckley and Schmidt (5) examined five fields, on the same site, that had been cultivated long term but with five different historically cultivated treatments, they found that the microbial community structures (T-RFLP analysis of 16S rDNA) were remarkably similar among these fields. Similarly, comparing improved and unimproved grassland soils, McCaig et al. (40) did not observe differences between the diversities of individual unimproved and improved grassland communities on the basis of 16S rDNA-derived DGGE profiles. However, DGGE banding profiles obtained from triplicate samples, pooled prior to analysis, indicated that there was less evenness in improved soils, suggesting that selection for specific bacterial groups had occurred.
As it has been well established that most of the bacterial soil community is inactive, a description of community structure on the basis of rDNA diversity provides a historical perspective rather than an immediate description of the present soil bacterial community (18, 64). Many bacterial species are known to vary their ribosome numbers in accordance with their cellular activities (66); thus, studying a bacterial community by rRNA-based approaches permits a description of the active community. The rRNA-derived DGGE profiles showed greater discrimination than either of the rDNA-based methods. Treatment differences were highlighted in the different relative abundances of sequences between EC2 (organic) and EC3,and between RMB (manure treated) and RMC. The five fields (RMB, RMC, RCA, RCC, and RCD) previously found to form one cluster on the basis of rDNA profiles (Fig. 3 and 7) were further separated by farm location into two clusters (Fig. 5). Fields RCC and RCD (typical stagnogley soil type with leguminous cropping) were also separated from field RCA (typical stagnogley soil type with nonleguminous cropping). Therefore, the rRNA diversity analysis would seem to be more indicative of short-term changes in land-use management.
In a study of soil bacterial community profiles in five Dutch grassland meadows, Felske et al. (19), using multiple competitive RT-PCR and DGGE analyses, discovered an increase in the total number of rRNAs after these fields were taken out of agricultural production. However, although the multiple competitive RT-PCR analysis indicated activity shifts for the predominant soil bacteria during the subsequent vegetational succession, quantitative dot blot hybridization failed to detect differences at a higher taxonomic level; although the vegetation clearly changed, the general composition of the bacterial community remained remarkably stable. Similarly, Buckley and Schmidt (5) failed to detect significant differences in the microbial group rRNA abundances in five fields despite differences in chemical inputs, tillage, plant community composition, and productivity. Therefore, it would appear that agricultural soil habitats are immensely stable environments where competition among microbial organisms is not severe, allowing species to persist throughout fluctuations in above-ground vegetation and in soil chemical properties.
Despite exhibiting the highest substrate utilization diversities, the EC fields showed the lowest relative diversities on the basis of DGGE analysis of rDNA (Table 3). Smalla et al. (58) discovered, by comparing the TGGE profiles of the original inoculum with those in the BIOLOG wells following incubation, that fast-growing bacteria that adapted to high substrate concentrations were numerically dominant in the BIOLOG wells. Therefore, the comparatively high diversities of sequences generated from rRNA-derived DGGE profiles of the EC soils may more accurately reflect the functional community. The high substrate utilization diversities in these soils may suggest an overlap between the metabolically active members of the soil community, as assessed from rRNA, and those active in the BIOLOG wells.
There seemed to be a strong link between the elevated nitrate concentrations within these soils and the significant increases in organic matter in these soils compared to their nonleguminous counterparts. This finding may indicate the presence of nitrogen-fixing bacteria and, in conjunction with the higher organic matter concentrations, may advocate the use of leguminous crop planting for improving soil health. The benefit to agriculture of nitrogen fixation from nodulated legumes has long been established (31, 57, 62, 68). In long-term field experiments on loamy sand and sandy loam, leguminous cropping stimulated microbial activity in the rhizosphere more than cereals, maize, or oil flax. Moreover, normal GAP farming may also decrease nitrogen fixation activity. Hoflich et al. (28) demonstrated that the leghemoglobin content of pea nodules (an indicator of nitrogen fixation activity) was reduced by the application of high levels of nitrogen in crop rotation. Herbicides widely applied for the suppression of weed growth in white clover seed crops have also been shown to be toxic to the growth of the nodule-forming bacterium Rhizobium trifolii and to the nitrogen-fixing symbiosis of these organisms (8).
In conclusion, the polyphasic approaches used to analyze bacterial community structure were generally in broad agreement, with samples being separated primarily by soil type. The community within the EC fields was discriminated by all methods, including BIOLOG. The community within field RMA was likewise discriminated by rDNA and rRNA analyses. Cluster analysis of the T-RFLP data and the rDNA-derived DGGE data grouped the remaining fields on the basis of soil type. Thus, we hypothesize that soil type is the overriding factor in community determinations at these sites, followed by the planting of, specifically, leguminous crops rather than other management practices. The bacterial profiles generated by RT-PCR of rRNA further separated the soil by geographic location and also revealed small treatment differences, providing the greatest discrimination between sites and emphasizing the importance of studying both total and active communities. Soil has been shown to have an immense capacity for diversity and therefore a large buffering capacity before the results of management practices will affect the dominant members of the community. However, the longer-term impacts of management practices may be much more significant, and their effects upon the fungal community were not examined in this study; therefore further attention is warranted on this subject.
We are grateful for the financial support of this research by Unilever and Syngenta.
|
|
|---|
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»