Abstract

Soil bacteria exhibit short-term variations in community structure, providing an indication of anthropogenic disturbances. In this study, microbial biomass carbon (MBC), potentially mineralizable nitrogen (PMN), community level physiological profiling (CLPP), and culture-dependent DGGE (CD DGGE) fingerprinting of the 16S rRNA gene were used to compare microbial communities in organic farm and pasture soils subjected to differing agronomic treatments. Correlation analyses revealed significant relationships between MBC, PMN, and data derived from microbial community analyses. All measures separated soil types but varied in their ability to distinguish among treatments within a soil type. Overall, MBC, PMN, and CLPP were most responsive to compost and manure amendments, while CD DGGE resolved differences in legume cropping and inorganic fertilization. The results support the hypothesis that culturable soil bacteria are a responsive fraction of the total microbial community, sensitive to agronomic perturbations and amenable to further studies aimed at linking community structure with soil functions.

1. Introduction

Microorganisms play essential roles in organic matter decomposition, nutrient cycling, and plant productivity [1, 2]. Parameters that integrate diverse microbial populations into a single measure, such as microbial biomass carbon (MBC) and potentially mineralizable nitrogen (PMN), historically have proven to be useful and are widely employed measures of soil quality [3, 4]. Microbial biomass C encompasses a small labile fraction of total soil organic carbon that responds actively to changes in soil fertility, supports soil aggregation, and can be related to environmental factors such as climate, soil moisture, texture, and organic matter quality [5]. Potentially mineralizable nitrogen provides an index of a soil’s nitrogen-supplying capacity and has been positively correlated with other chemical and physical indicators of soil quality [4].

Various measures of functional and structural diversity in microbial communities have been proposed as appropriate indicators of changing soil quality [3, 6]. Community-level physiological profiling (CLPP) measures soil functional diversity by characterizing the relative utilization of a suite of carbon substrates. Community-level physiological profiling is a culture-based enrichment method that primarily characterizes and selects for fast-growing organisms that may be distinct from dominant bacteria in soil inocula [7, 8]; therefore, the ecological significance of CLPP data sometimes is questioned [9, 10]. However, numerous studies have shown that CLPP is appropriately responsive to a wide variety of agronomic perturbations [1118].

Several methodologies have been developed to characterize microbial structural diversity that rely on nucleic acid extraction followed by gene amplification and community fingerprinting [19, 20]. Fingerprinting techniques cannot resolve total species richness and only profile a small fraction of the total bacterial diversity depending on the resolution of each individual method [21, 22]. A current paradigm states that culture-independent molecular methods provide a better representation of soil bacterial communities, since theoretically they include novel, numerically abundant, but difficult-to-culture groups such as the Acidobacteria and Verrucomicrobia [2326]. However, an alternate hypothesis suggests that fingerprints of readily culturable bacteria may be more responsive to anthropogenic perturbations since they often exhibit rapid growth and produce substantial biomass [27]. Recently, we investigated a plate wash method to construct culture-dependent denaturing gradient gel electrophoresis (CD DGGE) profiles of the 16S rRNA gene [28]. We reported that CD DGGE profiled distinct fractions of the soil bacterial community that were not captured using culture-independent DGGE. In addition, the CD DGGE gels showed clear banding patterns of sufficient complexity and variability to investigate differences in microbial communities. The objectives of the present paper were to: (1) examine the responsiveness of CD DGGE fingerprints from pasture and organic farm soils that have been exposed to differing fertility regimes and (2) explore potential relationships between the CD DGGE profiles and MBC, PMN, and CLPP.

2. Materials and Methods

2.1. Description of Field Sites

The West Virginia University (WVU) Organic Research Farm is located in Monongalia County, West Virginia, United States of America ( N, W). Soils on the farm are Dormont and Guernsey silt loams (fine loamy mixed, superactive, mesic, Oxyaquic Hapludalfs) [29]. The market garden has a completely randomized design with four plots (4.9 × 7.6 m) per treatment [30]. Treatments included fertility input (compost and green manure) and annual crop rotation among four plant families (Fabaceae, Solanaceae, Cucurbitaceae, and Asteraceae). High-input plots received composted dairy manure (composted 1 : 1 cow manure and leaf litter applied at 22.4 Mg based on wet weight), and green manure comprised of approximately equal parts rye (Secale cereale) and vetch (Vicia villosa) applied at approximately 10 Mg and tilled into the soil 1–3 days prior to planting. Low-input plots received only green manure applied as described above. Treatments were applied to plots in late April and early May. Details of compost analyses, soil fertility, and crop yield in these plots have previously been reported by Childers [30].

The West Virginia University Reedsville Experiment Farm is located in Preston County, West Virginia, United States of America ( N, W). Soils on the farm are mapped in the Gilpin (fine-loamy, mixed, semiactive, mesic Typic Hapludults) and Wharton (fine-loamy, mixed, mesic Aquic Hapludults) soil series [29]. Over the past decade, this experimental pasture has been limed to maintain a pH circa 6.5 and fertilized with two levels of inorganic fertilizer inputs, that is, high fertility (P 134 kg and K 390 kg ) or low fertility (P 67 kg and K 195 kg ). The predominant pasture grasses were orchard grass (Dactylis glomerata), Kentucky bluegrass (Poa pratensis), and meadow fescue (Festuca elatior) interspersed with white clover (Trifolium repens) and red clover (Trifolium pratense). The percentage of legumes averaged 36% and 10% in high- and low-fertility plots, respectively [31].

2.2. Soil Sampling

On the WVU Organic Research Farm composite (12 cores), surface soil samples (15 cm depth × 2.5 cm diameter) were collected in June and August from each of high- and low-input plot that had been planted with beans (Phaseolus vulgaris), peas (Pisum sativum), and tomatoes (Lycopersicon esculentum) and peppers (Capsicum annuum). Ten soil cores (16 cm depth × 2 cm diameter) were collected from each plot at ~1.5 meter intervals along two transects at the midpoint between crop rows. Soil samples were bulked, placed in polyethylene bags, transported on ice to the laboratory, sieved (<2 mm), and refrigerated at 4°C.

On the Reedsville Experimental Farm, fenced experimental subplots were established in June on adjacent high- and low-fertility pasture sites and amended with a manure treatment. Each plot was subdivided into quadrants (1 m × 1 m), which were randomly assigned to receive a manure patch or remain as an untreated control. Manure patches (25 cm diameter, 5 cm height) were fabricated from fresh manure collected from the same fields. Five replicate quadrants from each treatment were sampled at days 0, 21, 40, 63, and 123 after site establishment. On day 0, only the non-manured (control) soils were sampled. Four to five of the samples were used for analysis, depending on the type of analysis. Twelve soil core samples (15 cm depth × 2.5 cm diameter) were collected from each quadrant. Core samples were bulked, placed in polyethylene bags, transported on ice to the laboratory, sieved (<2 mm), and refrigerated at 4°C.

2.3. Microbial Biomass Carbon and Potentially Mineralizable Nitrogen

Microbial biomass carbon was extracted from moist soil (80% water-filled porosity) according to the method of Islam and Weil [32]. Total carbon in soil extracts was measured using a Dohrmann DC-190 TOC analyzer (Rosemount Analytical Inc., Santa Clara, Calif, USA). PMN was determined according to the methods described in Keeney [33]. Ammonia in each sample was measured using a Lachat flow injection analyzer (Hach Company, Loveland, Colo, USA) using the phenolate method.

2.4. Community Level Physiological Profiling

The ability of soil microbial communities to utilize a variety of carbon sources was assessed using community-level physiological profiling originally developed by Garland and Mills [34]. Sterile waring blenders were used to homogenize (1 min, max speed, 3x with intermittent cooling on ice) 20 g of soil in 180 mL Winogradsky Salts Solution (WSS) [35]. Homogenate was diluted (10−3) and used to inoculate (100 μL well−1) BIOLOG ECO Plate MicroPlates (Biolog Inc., Hayward, Calif, USA). Each plate contained triplicates of 31 individual carbon substrates, plus 3 control wells lacking a carbon source. The microplates were incubated at 25°C for five days. Utilization of the carbon source was monitored by measuring absorbance at 590 nm. Readings were taken every 24 hours using an automated plate reader (SpectraMAX 340 pc, Molecular Devices, Sunnyvale, Calif, USA). Absorbance values from the late-log phase of the growth curves (63 hours) were used to calculate diversity indices and principal component as described in Section 2.6.

2.5. Culture-Dependent Denaturing Gradient Gel Electrophoresis

A detailed description of this procedure is presented in Edenborn and Sexstone [28]. Briefly, twenty grams (dry weight) of composite soil samples from each plot were placed in sterile Waring blenders and shaken for one minute at high speed (3x with intermittent cooling on ice) in 180 mL of sterile WSS. Blended suspensions were diluted (10−2) in sterile WSS, and 100 μL was plated onto each of four replicate agar plates (final dilution = 10−3) and incubated aerobically for two weeks at 25°C. Cultured cells were washed sequentially from four replicate R2A plates using WSS (4 mL) and a sterile disposable inoculating loop. The resulting cell suspensions were vortexed and frozen (−20°C). DNA was extracted from aliquots of each thawed cell suspension (1.8 mL) using the MoBio Microbial DNA extraction kit (MoBio Labs, Carlsbad, Calif, USA). The variable V3 region of the 16S rRNA gene from the domain Eubacteria was amplified using the PRBA338F primer with PRUN518R primer as previously described [28]. Composite PCR products (~150 ng) were loaded into a polyacrylamide gel (8%) with a 40% (16.8 g urea; 16 mL formamide/100 mL) to 60% (25.2 g urea; 24 mL formamide/100 mL) denaturing gradient and electrophoresed for 14–16 hours at 50 V and 60°C using the DCode Universal Mutation Detection System (BioRad). The gel was visualized by silver staining according to the protocol of Caetano-Anolles and Gresshoff [36]. The developed gel was scanned using Hewlett Packard ScanJet 7400c (Hewlett-Packard Co., Palo Alto, Calif, USA) at a resolution of 600 dpi and saved in uncompressed TIFF format for further analysis.

In pasture soils, temporal changes in CD DGGE fingerprints were determined following application of manure pats to pasture soil treatment plots. Prior to application of manure (day 0), five replicate soil samples were collected from high- and low-fertility plots and the culturable bacterial communities characterized as the manure pats decomposed. Additional samples were collected on days 21, 40, and 63.

2.6. Data Analysis

Differences between MBC organic farm (OF) and pastures (PAS) soils were tested using an analysis of variance (ANOVA) in a 2 × 3 factorial design (factor A: high and low input (OF) or high-and low-fertility (PAS); factor B: time-crop combination (OF) or time-manure amendment combination (PAS)). PMN in organic farm soils was compared using a 2 × 2 Factorial (Factor A: high-and low-input; Factor B: time-crop combination). For pasture soils, the same factorial design used for MBC was used for PMN. Tukey’s HSD test was used for post-ANOVA pairwise comparisons.

Total well color development (TWCD) for BIOLOG ECO plates was calculated as the sum of absorbance values determined after 63 hours of incubation. Average well color development (AWCD) was calculated from the mean for those 31 absorbance values. Substrate diversity (H) was calculated as , where is the proportion of total microbial activity on a particular carbon source. Substrate richness (S) was calculated from the total number of wells with absorbance values greater than 0.25. Substrate evenness ( ) was calculated as (S). Data were normalized by dividing each substrate absorbance value by the AWCD [34]. Robust principal component analysis (PCA) was done on a covariance matrix of normalized data using alpha = 0.

DGGE gels were analyzed using Quantity One gel analysis software (BioRad, Hercules, Calif, USA). Analysis of the bands was done by setting background subtraction at 15 using a rolling disk method and generating intensity profiles to evaluate and compare peak size, shape, and location. Bands were picked manually for a limit of intensity greater than 0.05. Similarity matrices were calculated with Sorensen’s index based on the presence/absence of peaks. Analysis of similarity (ANOSIM) was used to determine significant differences between DGGE banding profiles. In ANOSIM, an R value greater than 0 indicates that objects (DGGE profiles) are more dissimilar between groups than within groups, and a P value indicates the level of significance [37]. Enhanced nonmetric multidimensional scaling (NMDS) was used to ordinate DGGE data [37, 38]. The ability of NMDS to provide a good representation of similarity data in two-dimensional space was evaluated based on Kruskal’s stress values. Stress values <0.15 indicate that the ordination plots presented were useful representations of the CD DGGE data [37].

Multiple and simple correlations analysis was used to assess associations between principal components, NMDS axis scores, MBC, and PMN data. Microbial biomass carbon and PMN were chosen as correlates because they are considered to be broad indicators of soil quality, are relatively easy to measure and interpret, and are linked soil functions such as nutrient cycling, decomposition of organic amendments, physical stabilization of aggregates, soil productivity, and N supplying potential [39].

Analysis of variance (ANOVA) and correlation analyses were performed using SAS (SAS Institute Inc., Cary, NC). Principle component analysis (PCA), ANOSIM and NMDS were performed using the R package Vegan [40].

3. Results

3.1. Microbial Biomass Carbon and Potentially Mineralizable Nitrogen

In the present study, levels of MBC were found to be ten times higher in pasture soils compared with the organic farm (Table 1). In both systems, amendment with manure or compost at different times resulted in significant differences in MBC. Microbial biomass carbon levels in organic farm soils were significantly affected by input ( , ) and the time-crop combination ( , ) with no interactions. The MBC in high-input plots (June) varied significantly for all treatments except high-input plots planted with legumes (August) (Table 1). In pastures, MBC was significantly affected by manure amendment and time ( , ), but not fertilizer amendment ( , ). The interaction was not significant ( , ). MBC was significantly higher in September than in May, and the greatest amount of MBC was found in plots that had been amended with manure (Table 1). Clear statistical differences between inorganic fertilizer applications were not demonstrated; however, lower levels of MBC generally were observed in low- compared with high-fertility plots, both in May and September.

Pasture soils on average contained four times more PMN than organic farm soils (Table 1). In organic farm soils, a significant difference in PMN was found due to input ( , ), but neither time-crop combination ( , ) nor the interaction of these two factors was significant ( , ). When individual means were compared, no differences were found in PMN between different organic farm soils (Table 1). In pastures, both fertilizer input ( , ) and time-manure amendment ( , ) were found to have a significant effect on PMN without any interaction ( , ). Potentially mineralizable nitrogen was significantly higher in the plots with manure compared to fertilizer input (Table 1).

3.2. Community Level Physiological Profiling

Shannon diversity indices ranged between 3.26–3.35 and 3.27–3.29; richness values ranged between 28–30 and 28–29; evenness between 0.975–0.982, and 0.975–0.977 for the organic farm and pasture soils, respectively. There were no significant differences ( ) in functional diversity indices based on soil type or agronomic treatment.

Ordination of the CLPP data using principal component analysis clearly distinguished low-input organic farm soils from other soils. We first performed PCA on all 31 substrates to identify the most strongly loaded, noncollinear substrates. Ten of these substrates were then used for a subsequent PCA [41]. These substrates were L-arginine, D-cellobiose, cyclodextrin, i-erythritol, 2-hydroybenzoic acid, D-lactose, L-phenylalanine, phenylethylamine, L-threonine, and D-xylose. Ordination of these data is presented in Figure 1 where the distinct cluster formed by low-input organic farm soils can be visualized. In this analysis, 75% of the variance was accounted for by the first three principal components. The first principal component was most strongly correlated with the utilization of L-phenylalanine (−0.94), D-lactose (−0.68), and L-threonine (−0.61). Principal component 2 was most strongly correlated with 2-hydroxybenzoic acid (0.87), and principal component 3 was most strongly correlated with phenylethylamine (0.95) and L-arginine (0.69).

3.3. Culture-Dependent Denaturing Gradient Gel Electrophoresis

CD DGGE profiles of 16S rDNA showed clear differences between culturable bacterial communities in organic farm and pasture soils (Figure 2(a)). Culturable bacterial communities in the two soil types were highly distinct (ANOSIMS; , ). Ordination of the data using NMDS demonstrated that legume cropping (OF) and fertilizer (PAS) had more impact on the structure of culturable soil bacterial communities than organic matter amendments (compost or manure) (Figure 2(b)).

At each sampling time, PCR products from three replicate plots per treatment were combined to create a composite profile. Significant differences were found between the high- and low-fertility plots on day 0 (ANOSIMS; , ). Nonmetric multidimensional scaling analysis of these data is presented in Figure 3(a). There were significant differences in bacterial communities over time (ANOSIMS; , ), but no significant differences were observed due to treatment ( , ). Ordination of the data also revealed distinct clustering based on sampling time (Figure 3(b)).

In June and August, soils collected from three replicate plots per treatment at the organic farm were used for CD DGGE analysis. There were no significant differences in CD DGGE profiles from different treatments in June (ANOSIMS; , ); however, significant differences were observed in August ( , ). Bacterial communities from tomato and pepper plots (high and low fertility) collected in both June and August clustered together (Figure 4(a)). In contrast, the communities from soils that had been planted with legumes clustered separately from those from soils planted with tomatoes and peppers and from each other (Figure 4(b)). Kruskal’s stress values for all analyses were <0.15.

3.4. Correlations between MBC, PMN, CLPP, and CD DGGE

Multiple correlation analysis was used to explore potential relationships between the principal component scores from the PCA analysis of CLPP data (Figure 1) and other measured variables. Significant correlations were observed between the CLPP principal components, MBC ( , ) and PMN ( , ). When organic farm and pasture sample were analyzed separately, the correlations improved ( , ) for MBC and PMN from both systems. Multiple correlations between the CLPP principal components (Figure 1) and the axis scores from NMDS of CD DGGE data from the organic farm samples in August 2004 (Figure 4(b)) also were highly significant (dimension 1: , ; dimension 2: , ). Multiple correlation analysis was not performed on pasture data because samples for CD DGGE and CLPP were collected on different sampling dates. Simple correlations were performed for both agronomic systems using axis scores from NMDS of CD DGGE data and the MBC and PMN data. In both soil systems, significant correlations were found between the DGGE axis scores and MBC (Figure 5) but not between DGGE axis scores and PMN (data not shown).

4. Discussion

Ideal soil quality indicators are those that reflect ecosystem processes and integrate physical, chemical, and biological properties of soil. They are also sensitive to environmental disturbances and are practical and economical for a wide-variety of users [42]. Microbial community structure is of potential use as an indicator of soil quality [3]; however, the linkage to ecosystem processes and soil status is not always evident in ways that can be easily integrated into predictive models [43]. Exploring the composition of microbial communities in soil often depends on computation and inspection of derived parameters from multivariate models and the construction of ordination plots. This approach often does not provide easily interpretable indices for the quantification and comparison of soil microbial diversity. Therefore, we attempted to correlate derived parameters from multivariate analyses with more straightforward measures of soil change and determined if profiles of culturable soil bacterial communities provide information about short-term changes in soil quality that are not characterized using conventional microbial soil quality indicators.

The organic farm and the perennial pasture represent contrasting research plots that have received regular soil amendments and plant management over the past decade. Long-term changes in soil properties suggest an overall improvement in the soil quality at both sites [30, 31, 44]; however, short-term effects of specific agronomic treatments have been more difficult to resolve. In this study, we explored the responsiveness of culturable soil bacterial communities to soil management and their potential relationship to more commonly employed measures of soil perturbation (MBC, PMN, and CLPP). Culturable bacteria were sensitive and early indicators of change in metal-contaminated environments [27, 45].

We tested the hypothesis that fast-growing culturable bacterial communities, as reflected by CD DGGE profiles, are sensitive to short-term changes in dynamic agronomic systems and that these responses can be linked to other microbial measures of soil status. Obvious differences were observed between CD DGGE fingerprints of 16S rRNA genes from the two soil types. However, CD DGGE profiles were not uniformly sensitive to within-site treatment differences. Rather, we found that CD DGGE measured differences in soils that were amended with inorganic fertilizers and/or cropped with legumes, while MBC, PMN, and CLPP were more responsive to additions of manure and compost to soils. It is quite possible that the fertilizer effect we observed in pasture soils also involved a legume effect. For the past decade, experimental pastures at Reedsville, WVa, have been subjected to equal grazing pressure and animal waste deposition but have received two different levels of phosphorus, potassium, and lime amendments [44]. The higher phosphorous treatment has been shown to favor the persistence of pasture legumes, such as red and white clover, in these soils [46]. Girvan et al. [17] hypothesized that soil type and legume cropping were the overriding factors controlling the structure of bacterial communities in agricultural soils. Additional studies support the observation that legumes and nitrogen availability can influence bacterial community [47, 48]. Conventional fertilization with inorganic phosphates also has been shown to increase the diversity and activity of free-living nitrogen-fixing bacteria compared with organic amendments [49].

We found positive correlations between CLPP principal component scores and MBC and NMDS axis scores and MBC in both systems. Previous studies also have identified correlations between the structure and function of culturable microbial communities and physicochemical and vegetative parameters [45], CLPP data, gross N flux, and soil respiration [5052] and organic matter content, pH, and sodium [53]. Although correlation does not imply causation, correlative relationships suggest that these variables respond in a similar manner to environmental disturbances [50] and provide the foundation to explore potential mechanisms that may link structure and function. For example, CLPP employs microbial growth on potentially important carbon substrates, and we observed significant correlations between the CD DGGE and CLPP. Exploration of the CLPP data using PCA analysis suggested that bacteria using aromatic substrates such as L-phenylalanine and 2-hydroxybenzoic may represent useful physiological groups to further explore using CD DGGE. Langenheder et al. [54] have suggested that bacterial community structure may be related to narrow functions such as degradation of semilabile carbon compounds rather than to broad functions such as respiration and biomass production. If this is accurate, then characterizing culturable bacterial communities with specific physiological capabilities or narrow niche preferences may improve our ability to distinguish between soils that have been impacted by anthropogenic perturbations. Recently, Ritz [55] has argued that culturable communities have no utility in contemporary microbial ecology and cannot be related to soil fertility. Our results do not support this conclusion. Rather, we agree with Nichols [56] that microbial cultivation is a re-emerging frontier that can provide information not directly observable from culture-independent and sequencing efforts alone. Although bacterial growth on different types of solid media is inherently selective, the ability to partition complex soil bacterial communities into definable physiological subsets using a variety of isolation media and incubation conditions should be viewed as a potential strength of the CD DGGE technique, since selected physiological groups may be most sensitive and responsive to environmental change. It is insufficient to merely catalog genes to understand how environmental factors govern shifts in community structure. Natural selection acts on phenotype, and physiological characteristics are often linked to more than one gene and are regulated by multiple external factors. Classifying microorganisms into meaningful functional groups and characterizing the physiological and ecological redundancies within these groups potentially can make community data more manageable and relevant [43].

5. Conclusions

The organic farm and pasture represent systems with contrasting soil types, land use, and fertility regimes. In this study, we used MBC, PMN, CLPP, and CD DGGE of the 16S rRNA gene to compare microbial communities as a function of agronomic amendments. There were significant correlations between MBC and PMN and variables derived from analysis of the functional and genetic diversity of soil communities; however, each measurement technique exhibited a different response to the fertility amendments. Overall, MBC, PMN, and CLPP were most responsive to additions of manure and compost to soils; however, only CD DGGE revealed clear differences in soils that were amended with inorganic fertilizers and/or cropped with legumes. Therefore, different measures of microbial community structure do not provide redundant information in soil quality analysis, and culturable bacterial communities can provide useful indicators of select agronomic perturbations.

The minority status of culturable bacteria does not preclude the possibility that these communities are useful to advance our understanding of how microbial communities are linked to environmental processes. Culture-based techniques such CD DGGE and CLPP have the potential to provide critical insight into the mechanisms that link structural and functional diversity, to explore ways to assess physiological and ecological redundancies within functional groups, and to advance our understanding of how microbial diversity contributes to soil quality.

Acknowledgments

The authors thank Dr. George Seidel for assistance with statistical analyses. This work was supported by grants from USDA and by the West Virginia University Davis College of Agriculture, Forestry, and Consumer Sciences through funding appropriated under the Hatch Act.