Land-use change and management practices are normally enacted to manipulate environments to improve conditions that relate to production, remediation, and accommodation. However, their effect on the soil microbial community and their subsequent influence on soil function is still difficult to quantify. Recent applications of molecular techniques to soil biology, especially the use of 16S rRNA, are helping to bridge this gap. In this study, the influence of three land-use systems within a demonstration farm were evaluated with a view to further understand how these practices may impact observed soil bacterial communities. Replicate soil samples collected from the three land-use systems (grazed pine forest, cultivated crop, and grazed pasture) on a single soil type. High throughput 16S rRNA gene pyrosequencing was used to generate sequence datasets. The different land use systems showed distinction in the structure of their bacterial communities with respect to the differences detected in cluster analysis as well as diversity indices. Specific taxa, particularly Actinobacteria, Acidobacteria, and classes of Proteobacteria, showed significant shifts across the land-use strata. Families belonging to these taxa broke with notions of copio- and oligotrphy at the class level, as many of the less abundant groups of families of Actinobacteria showed a propensity for soil environments with reduced carbon/nutrient availability. Orders Actinomycetales and Solirubrobacterales showed their highest abundance in the heavily disturbed cultivated system despite the lowest soil organic carbon (SOC) values across the site. Selected soil properties ([SOC], total nitrogen [TN], soil texture, phosphodiesterase [PD], alkaline phosphatase [APA], acid phosphatase [ACP] activity, and pH) also differed significantly across land-use regimes, with SOM, PD, and pH showing variation consistent with shifts in community structure and composition. These results suggest that use of pyrosequencing along with traditional analysis of soil physiochemical properties may provide insight into the ecology of descending taxonomic groups in bacterial communities.
Citation: Shange RS, Ankumah RO, Ibekwe AM, Zabawa R, Dowd SE (2012) Distinct Soil Bacterial Communities Revealed under a Diversely Managed Agroecosystem. PLoS ONE 7(7): e40338. https://doi.org/10.1371/journal.pone.0040338
Editor: David L. Kirchman, University of Delaware, United States of America
Received: February 29, 2012; Accepted: June 6, 2012; Published: July 23, 2012
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Funding: The work for this publication was funded by a United States Department of Agriculture-McIntyre-Stennis Forest Service Grant (ALX-MSFSW 0219038; http://www.csrees.usda.gov/about/offices/legis/mcintirestennis.html), and a United States Department of Agriculture-National Institute of Food and Agriculture Grant (ALX-SWQ 0211874; http://www.csrees.usda.gov/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared the following interest. Scot E. Dowd is employed by Molecular Research LP. There are no patents, products in development or marketed products to declare. This does not alter the authors' adherence to all the PLoS ONE policies on sharing data and materials, as detailed online in the guide for authors.
Land use change and management in the conversion of forested, pastured, and cropped land have ecosystem-scale impacts such as: soil cycling of organic compounds –, biodiversity –, and soil nutrient dynamics . The sensitivity of microbial communities to changes in management as well as their importance to nutrient cycling is a reason why they have been considered as early indicators of change in the quality of the soil ecosystem . Although changes in bacterial communities have been reported in prior research , , advances in sequencing technology have provided researchers with the ability to assess bacterial diversity at lower costs, and quicker turnaround than prior 16S rRNA and sequencing methods. These advances have allowed researchers to enhance work in the much needed area of bacterial community structure at varying scales –.
Physical disturbance of the soil has been reported as being a crucial factor in determining soil biotic characteristics in agroecosystems . The loss/disturbance of a stratified soil microhabitat has been attributed to a decrease in the density of species that inhabit agroecosystems. It has been proposed that such soil biodiversity reductions may be negative because the recycling of nutrients and proper balance among organic matter, soil organisms and plant diversity are necessary components of a productive and ecologically balanced soil environment . For example, Acosta-Martinez et al.  evaluated the physical disturbance to bacterial communities with respect to tillage and reported that tillage reduced the bacterial diversity due to the interruption to physical diversity of the soil environment. Torsvik et al.  conducted a study comparing biodiversity of bacterial communities in perturbed soils versus relatively undisturbed, pristine soils, and found that human-induced pollution can lead to profound changes in microbial community structure by greatly reducing bacterial diversity. Other studies that have introduced organic amendments to soils have showed changes in microbial community structure in addition to biochemical changes , .
Three primary areas of land use in the southeastern United States include pine plantation, annual cropping, and livestock pastures. The incorporation of forestry (particularly pine plantation) into agroecosystems has become an important economic, social , and environmental  issue in ecosystem management of the southeastern USA. Grazing of pine forests to control understory growth is one of the management systems which have been proposed to further diversify agroecosystems, and manage fire hazards , . All three of these land use types have the potential to alter soil chemical and physical disturbance, and ultimately to influence changes in soil bacterial community structure and composition. The ecological impact of management practices has been demonstrated as being a consistent source of disturbance to soil ecosystems . Participating landowners and farmers are often engaged in multiple management practices within the same agroecosystem, creating a level of ecological complexity that is difficult to replicate in controlled settings. Hobbs and Huenneke  concluded that such disturbance has the potential to alter plant communities, but this process is poorly understood in the case of microbial species. It is therefore the objective of this study to utilize 16S rRNA sequencing techniques to elucidate the bacterial community structure and composition that exists in a single soil type exhibiting the three land use types discussed above. DNA was extracted from soils under each type of land use introduced above, amplified, and subjected to pyrosequencing. We therefore hypothesized that the specified land use types across the landscape will exhibit distinct soil bacterial communities.
Results and Discussion
Richness and Diversity Estimates
The maximum operational taxonomic units (OTUs) detected across the study site according to the observed clusters (sobs) at 3% dissimilarity was 2065 (Table 1) in the Cultivated sample 1, despite there being an obvious dominance of sequence recovery from the pastured samples. The maximum amount of OTUs observed is also reflected in the Chao and ACE values predicted for that sample. Though the highest number of OTUs was observed and predicted for this particular sample, the richness values were quite variable. All four of the estimators followed a distinct trend (pine forested<cultivated<grazed pasture) (Figures 1a and 1b). That trend exemplified an agroecosystem in which the lowest number of OTUs was found in soils sampled under the forest land use system and the highest under the grazed pasture soils. The indices reflected the same trend when calculated at 5% dissimilarity; however no differences were found between the land management systems for any of the diversity estimators.
Richness/diversity estimators (a) and rarefaction curves (b) are presented as calculated by MOTHUR at a level of 3% dissimilarity.
In the richness/diversity data presented in Figures 1a and 1b, the grazed pasture system further establishes this point as its richness trends higher than the pine forest for all three of the estimators and trended higher than the community under the cultivated system. Increased richness in semi-stable pastures can be expected, but the added effects of grazing has been reported by other researchers as serving to enhance microbial communities attributing to three major factors: organism and substrate diversification from fecal and urine deposition, stimulation rhizosphere activity as a result of mowing, and the mixing and dispersal of microbial communities through trampling –. This pattern is not only noticeable in richness, but is also reflected in the highest SOC, TN, PD, and APA values (Table 2). Though the richness values of the cultivated community are not as high as those of the grazed pasture, they do not show a significant difference. In prior studies ,  no significant differences were detected between agricultural and pastured soils. It may be assumed that bacterial community richness under cultivated management would differ due to habitat difference introduced by inorganic amendments and tillage, but Acosta-Martinez et al.  suggest diverse crop residues and root systems may be cause for increased richness in diversely planted cultivated soils.
The forested system is the most distinguishable land use regime according to the acidic pH (Table 2) and the reduced amount of OTUs detected in the communities. In figures 1a & 1b, the forested soil bacterial community consistently shows the lowest values for richness. The forested soil presents harsh conditions to resident organisms, specifically in the acidic pH and polyphenolic compounds found in loblolly pine litter . Soil pH has been found in numerous studies that utilize 16S sequencing techniques, as the driving abiotic factor shaping community structure in soils , , , . Though not directly measured in the study, polyphenolic compounds resulting from leaf litter have been reported as having a negative impact on soil bacterial communities ,  as well. This chemical disturbance creates a community with distinct abiotic and biotic characteristics, which impacts the growth and stability of the soil bacterial community.
Relative Abundance of Bacterial Phyla and Classes
Bacterial community compositions of the soils were examined at descending levels of biological classification to determine the effect of diversified land use on community membership. Detailed phylogenetic analyses grouped the soil associated bacterial sequences into 26 phyla (including unknown). The relative abundances of the 10 most abundant phyla are presented in Figure 2. The phylum distribution showed that Proteobacteria was the most dominant phyla with mean relative abundance values ranging from 39.9 % in cultivated soil, to 41.0 % in forested soil and to 45.5% in pastured soil. Other dominant phyla of note were Actinobacteria (representing 20.1–34.3% of the bacterial sequences in the samples) and Acidobacteria (which represented 4.4–20.6% of the bacterial sequences in each sample). Results from running GLM substantiate that the land uses are significantly different when considering the composition of the bacterial communities. Specific differences were observed between the cultivated and forested land use types for the phyla Acidobacteria (P<0.01), while Actinobacteria showed significant distinction (P<0.05) under the pine forested land use system.
Phyla included in this figure had relative abundance values consistently greater than 1%, as well as the abundant classes of the phylum Proteobacteria. Values presented are the mean percent. CLT = Cultivated; GRP = Grazed Pasture; PNF = Pine Plantation.
The remaining phyla accounted for fewer than 25% of the relative abundance observed and were designated as minor (reference to abundance only). Of these groups, Gemmatimonadetes showed significantly distinguishable populations in all three land use types (forested/cultivated, P<0.001; forested/pastured, P≤0.001; cultivated/pastured, P<0.05), while Chloroflexi (forested/cultivated, P<0.05; forested/pastured, P<0.05) relative abundance is distinct in the forested system, and Verrucomicrobia (pastured/cultivated, P<0.05; forested/cultivated, P<0.05), BRC1 (pastured/cultivated, P<0.001; forested/cultivated, P<0.001) and Nitrospirae (pastured/cultivated, P<0.001; forested/cultivated, P<0.001) showed significant distinction in the cultivated area. Another significant difference observed the pastured and cultivated land use was in the phyla Cyanobacteria (P≤0.005). Shannon-Wiener indices for evenness calculated at 3% dissimilarity showed that the cultivated system had the highest value of evenness (0.93) compared to the pastured (0.89) and forested (0.87) with significant differences (P<0.05) occurring between the cultivated and forested systems. When considering the difference of the composition of the bacterial community under the cultivated land use compared to the other management systems, it has been suggested  and supported  that the cultivation of a field, no matter how long ago, leaves an indelible imprint on the soil bacterial community. Through tillage, organic matter is incorporated throughout the plough layer of the soil and benefits unique microbial communities as the community reverts to an earlier (and more unstable) stage of ecological succession . These communities are described as having quick responses to conditions of feast or famine that would be created by growing season-fallow pattern in cultivated systems. Figure 2 exhibits this phenomenon as the minor phyla in the cultivated system accounted for almost 20% of the sequences detected, and showed a more even distribution of abundance. With the values of evenness presented above, the data provides further evidence of a more even community in the cultivated system disturbed by tillage.
At the class level, α-proteobacteria, Acidobacteria (class), and Actinobacteria (class) were the dominant bacterial classes found across the site. Of the classes of Proteobacteria, α-proteobacteria appeared to be the most dominant among them. The remaining classes of Proteobacteria (with exception to γ-proteobacteria) showed higher relative abundance compared to that of any other taxa at the class level (excluding class Actinobacteria). Based on GLM of transformed data, relative abundance was significantly higher for α-proteobacteria and Acidobacteria in the forested than the cultivated system, while β-proteobacteria and δ-proteobacteria were significantly higher in the cultivated system (N = 3; P<0.05). Seemingly, the pastured area was transitionary in its bacterial class abundance, as it differs only from the forested system only including an increase in Acidobacteria (N = 3; P<0.01) and a decrease in β-proteobacteria and δ-proteobacteria (N = 3; P<0.05). The observed shifts in taxonomic groups may suggest that microbial community composition changes in response to land management or environmental perturbation.
It has recently been reported that Actinobacteria and β-proteobacteria follow copiotrophic lifestyles , , while Acidobacteria can be classified as oligotrophic. The results of the other soil properties present the land use regimes as environments that could support both types of lifestyles. The copiotrophic environment could be considered as the pastured system which has significantly higher SOC, TN, and nutrient cycling due to the activity of alkaline phosphatase and phosphodiesterase. The oligotrophic environment could be then considered as the forested system which is lacking in SOM and enzymatic activity as compared to the pastured system. The trend with respect to the two trophic groups can be seen in Figure 2, in which the class Actinobacteria and the classes of Proteobacteria (with the exception of α-proteobacteria) trend down in the forested system with the opposite being observed for Acidobacteria. The instability of the cultivated community may also explain the lowest values for key taxa in the study. Even though α-proteobacteria was not explicitly designated as an oligotroph, in this study, as well as in other studies , , its members have shown the ability to outcompete most (if not all other) bacterial classes in acidic soils. What this data suggests is that the chemical ecology of the forested system may provide the foundation for significant differences in community membership as well as growth. Studies have also shown that shifts from forest to grassland soil  as well as cultivated to pasture  result in changes to bacterial community composition. The assertion of trophic lifestyles was stated as a general guideline, so a deeper look into descending taxonomic groups was observed for these particular classes to determine if this assertion holds at these levels.
A heatmap (with hierarchal clustering) (Figure 3) of the bacterial families found in the samples was generated with respect to the classes previously mentioned as following copiotrophic/oligotrophic lifestyles. The heatmap demonstrates that the top three families found in the study belong to the most abundant classes Actinobacteria and α-proteobacteria. The forested system appears to have the most distinction and was verified with ANOVA. The forested system differed from the other two systems significantly (P<0.05; n = 3) for the families: Acidobacteriaceae, Solibactereaceae, Solirubrobacterales (unclassified), Solirubrobactereaceae, Patullibactereaceae, Micromonosporaceae, Streptomycetaceae, Nocardioidaceae, Sphingomonodaceae, Alcaligenaceae, Burholderiaceae, Haliangiaceae, Polyangiaceae, Pseudomonadaceae. Out of the 42 Families graphed from these groups, less than half show abundance patterns that would support this assertion. As these are generally the most abundant families, prior assertions may be biased to those groups who are most identifiable using current methods. The evenness of the bacterial community under the cultivated system that was observed in Figure 2 and the reported Shannon evenness values is again exemplified by abundance spread further across the specified families than in any of the other systems. This is shown by the coverage of purple to orange colors in the columns labeled for the cultivated samples. The opposite is true for the Forested system, in which there is a large and visible gap in class abundance (exhibited by the large amount of black cells) in the columns designated as those belonging to the forested system. The samples of the Pastured and Cultivated systems seem to be more closely related in composition, as their clustering in the X-direction shows divergence at ∼1.05, while the Forested system diverges at ∼1.40 from the other two. Wardle et al. have also suggested that microbial community diversity/richness may not respond to cultivation in such a way that can distinguish it from a community with fewer disturbances .
A double cluster dendrogram that demonstrates the relative abundance of Families across the 9 samples across the three land use systems. Clustering in the Y-direction is indicative of abundance, not phylogenetic similarity. RA = Relative Abundance; CLT = Cultivated; PST = Grazed Pasture; FST = Grazed Pine Plantation.
Diversity and Abundance of Actinomycetes
To further demonstrate the differences in bacterial community composition, relative abundance was assessed at the level of Order for Actinobacteria as well. Table 3 displays the relative abundance and Shannon diversity indices of the most salient orders of Actinobacteria identified in the soils across the three land use types. The orders Solirubrobacterales, Actinomycetales, 0319-7L14, MC47, Acidomicrobiales, and Unclassified Actinobacteria (class) were identified as contributing substantially to the relative abundance of this class. Generally, both relative abundance and the Shannon index were lowest in the forested system for all Orders of Actinobacteria with the exception of Acidomicrobiales, which actually had the highest relative abundance and was second to the pastured system with respect to diversity.
These results suggest that pH shapes the Acidobacterial community structure more than the other environmental variables measured. Although the chemical environment seems to play a major role in the patterns of abundance, some of the taxa seem to respond to other factors. From the orders of Actinobacteria detected across the study site it is observed that Actinomycetes are the most abundant order (Table 3), and show a trend in their relative abundance and diversity across the site that reflects the amount of disturbance (pine forested<grazed pasture≤cultivated). Actinomycetes have been shown to have a preference for animal and human activity , and may be a possible biomarker in determining effects of land use change. Prior studies utilizing different primers from various geographic locations have suggested that Actinomycetes show higher abundance in agricultural and pasture soils, , , . Although Solirubrobacterales has not been extensively studied, recent studies have shown its members to be adaptive in their ability to colonize different ecosystems: fungal growing ant colonies , spinach phyllosphere , desert and Antarctic soil , . In this study, the order seems to favor physical disturbance as its highest relative abundance is in the cultivated system, and its highest diversity is in the pastured system. In this way, Solirubrobacterales performs similar to Actinomycetales with respect to disturbance.
Principle Coordinates Analysis
Unifrac metrics were used to assess community similarity between two or more samples according to their structure (weighted/quantitative) and membership (unweighted/qualitative). In the 2-dimensional plot visualized from the Unifrac weighted distance matrix principle coordinates analysis (3% dissimilarity), the samples of each system distinctively responded to the majority of the variation detected in the samples across two axes (Figure 4a). Axis 1 accounted for 34.6% of the variation, and Axis 2 accounted for 14.8% of variation. In Figure 4b, the same 2-dimensional plot was shown for the unweighted method which shows that in consideration of community membership, samples from the same type of land use system clustered together, although less distinctive (Axis 1 = 18.5%, Axis 2 = 14.1%).
PCoA plots are presented of the first two axes based on (a) weighted and (b) unweighted Unifrac distance matrices showing the quantitative and qualitative clustering of samples. The pastured system is represented by the red; the cultivated by the yellow; and the forested by the green.
The results from the Unifrac weighted and unweighted PCoA plots demonstrate the distinction that the bacterial communities under different land use management in both structure and composition. This is perhaps the most persuasive evidence that the three types of land use cultivate three distinctive bacterial communities. In the weighted PCoA plot (Figure 4a) the communities are clearly distinguished from one another quantitatively. There does appear to be some divergence within the samples from the cultivated system, as sample Cultivated 1, differs in its location on the y-axis. The divergence suggests that this sample exhibits different structure than the other two samples. This same phenomenon was evident when calculating the richness of the cultivated system, which explains the high error displayed in Figure 1a (as richness is considered a component of community structure). Regardless of its divergence, this community still seems distinct as clear separation is still detectable from the other land use types, primarily in the x-direction. The communities in the unweighted PCoA plot (constructed with respect to composition) show clear distinction as well (Figure 4b), as the samples from each of the land use systems cluster within the system. The pine forested and cultivated systems show the most dissimilarity with respect to structure and membership when considering the amount of space which separates them on the plot. Although the cultivated and pastured systems show the most similarity in physiochemical characteristics and richness/diversity values, they did show distinction in structure. Wardle  suggested that the difference in soil microbial communities of lightly disturbed (reduced/no till) and heavily disturbed (conventional tillage) systems is exemplified in the composition of the communities as seen in Figures 4a and 4b.
The Venn diagram presented in Figure S1 demonstrates the distribution of phylotypes at the 3% level of dissimilarity. The total number of phylotypes found in the entire agroecosystem is 7771. All three land use systems shared 2.2% of these phylotypes, while the most were shared between the pastured and cultivated systems (8.9%). The cultivated and forested systems share the least amount of phylotypes (2.5%), and the forested and pastured system share 4.5%.
Clustering results showed three discrete communities present under the agroecosystem. Though the grazed pasture and forested systems show the greatest dissimilarity in the data, it may be suggested that the grazed pasture and cultivated systems may show more similarity. The Venn diagram further establishes this point. Of the 11.2% of the phylotypes that the grazed pasture system shares, 8.9% is shared with the cultivated community, as compared to the 4.5% that it shares with the forested community. Even though these similarities are observed, the data also distinguishes the grazed pasture and cultivated systems in the Unifrac PCoA plots (Figures 4a & 4b). This distinction suggests that regardless of the phylotype sharing there are significant differences in the structure and membership of the communities. With the unique membership and structure elucidated by the three systems studied on a single soil type, the development of these communities with respect to management raises further questions as to the influence of environmental history and dispersal within the soil ecosystem.
In conclusion, components of this study further establish the distinction of soil bacterial communities with respect to differing land management strategies. As has been demonstrated in other studies using both 16S rRNA and PLFA techniques we notice distinction of forested soils with respect to community composition, structure, and phylogeny. The bacterial communities under in the cultivated and grazed pasture soils were similar in richness and composition, furthering the point that in a community with moderate disturbance, new individuals and groups could be introduced in a manner that promotes competition and diversity of the community, thus establishing a more stable community . Though community richness and composition observed under the cultivated and grazed pasture regimes were not distinguishable, metrics assessing phylogenetic differences showed distinct communities under all three land use types. Further analysis sheds more light on the behavior of soil bacterial Orders and Families under differing land management strategies. As has been noted in other studies, phyla and class taxa behave according to predetermined guidelines regarding the chemical ecology of microbial communities, but at the lower levels of taxa a majority of groups depart from these generalizations. Further study will be needed to actually elucidate and validate the divergence in ecology of the descending taxa.
Materials and Methods
Sites and Sampling for Physical and Biochemical Analyses
The study site was located on the Sundown Ranch demonstration farm located at 32° 26′ latitude and −87° 27′ longitude on 17 hectares of land in Perry County, Alabama, USA. The soil series of the site was Kipling clay loam (Fine, smectitic, thermic Vertic Paleudalfs). For the past 10 years this land has been used as a demonstration farm for a consortium of Alabama Land Grant Universities.
A preliminary geostatistical study of soil biochemical characteristics provided the initial evidence that soil biochemical and biological factors spatially vary with respect to land use type on this site. The preliminary study allowed for the reasonable designation of three major sampling areas that coincided with land management considering the spatial variability of the soil surface. The sampling methodology of choice was a stratified random sampling design that utilized sampling areas that coincided with specific land use types. Effort was made to obtain as statistically useful data set as possible that was well distributed throughout the site, while achieving a cost effective sampling method. The three sampling areas were constructed of sizes which were reflective of the amount of area each covered of the farm and overlaid by sampling grids. Sample numbers were assigned to each vertex in the grids, and the actual sample locations were randomly selected from these grid points. The amount of samples collected were proportional to the land area each land use system covered (loblolly pine forest = 8 ha; grazed pasture = 5 ha; cultivated = 3 ha).
The grazed pastured area is where pasture poultry are grazed twice yearly, and goats graze moderately throughout the year. This has been the management practice for the past 5 years. The loblolly pine plantation, which accounted for approximately half of the total area of the site, was designated as the forested system. It is characterized by a pine plantation at about 15 years of growth, and is lightly grazed by a goat herd. The crop-growing area was designated as the cultivated system, was used annually as a mixed cropping system (corn, beans, tomatoes, squash, and okra). This area has been under differential cropping management during the 10 years with the last 5 years being plastic mulch, inorganic fertilizer, and no herbicide.
Sampling was conducted on 14 October 2008 following the summer growing season. Soil samples of approximately 120-g field-moist weight were systematically collected from the upper 15 cm of soil at 45 sampling points in the landscape using a soil auger. In between sampling, the auger was sterilized with ethanol. Samples were preserved on ice during transport to the laboratory. Upon arrival soils were shortly stored at field moist conditions at 4°C. Prior to biochemical and physical assays soils were air-dried for 48 h, plant residues were removed by hand and the soil was sieved using a 2 mm mesh and mixed thoroughly thereafter. Five grams of soil were separated and kept at 4°C for molecular analysis.
Soil Enzyme Activity
Enzyme analysis was performed according to the methods of Tabatabai , with slight modification. The artificial substrate, p-nitrophenyl (1 mL, 0.05 M), and a pH buffer (pH values were 11 for Acid Phosphatase [APA], 6.5 for Acid Phosphatase [ACP], and 8 for Phosphodiesterase [PD]) were incubated in 25 mL glass flasks and capped at 37°C for 1 h with 1 g of soil. At the end of incubation, enzyme activity was stopped by addition of 4 mL of 0.5 M NaOH for phosphomonoesterases and 4 mL of 0.5 THAM-NaOH for phosphodiesterase followed by extraction with 1 mL of 0.5 M CaCl2. The mixture was then filtered (Whatman No. 2) and the extract analyzed using a Genesys 10 VIS spectrophotometer at (Thermo Fisher Scientific Inc., Waltham MA, USA) at 420 nm. Enzyme activity in filtrates was determined from a standard curve developed using p-nitrophenol standards. To account for non-enzymatic hydrolysis, values for controls were subtracted from sample readings. Toluene was not used in accordance with Bandick and Dick  and Elsgaard et al. , who showed that with incubation periods fewer than two hours, the absence of toluene was inconsequential to measured enzyme activity. All enzyme activities reported are expressed on a moisture-free basis.
Soil pH and SOM
Soil pH was analyzed using a 1∶2 soil/water ratio according to those methods described in McLean . The samples were then analyzed with an S500 pH Meter (following calibration at a pH of 4 and a pH of 10) [A. Daigger & Co., Vernon Hills, Illinois USA], measured to the nearest 0.01. For soil organic carbon (SOC) and total nitrogen (TN), air-dried soils were transferred to Auburn University Soil Testing Laboratory for analysis by the dry combustion method using an Elementar Vario Macro Combustion Analyzer (Elementar Americas, Inc., Mt. Laurel, NJ).
Soil Textural Analysis
Soil textural analysis was determined using the Buyocos hydrometer method for 40 g duplicate sample and <2 mm diameter. To enhance dispersion of the soil particles chemically, 100 mL of 0.05 g mL–1 Sodium Hexametaphosphate (HMP) solution was added to each sample and shaken (for mechanical dispersion) for 12 h . After shaking, the HMP solution and soil were transferred into a 1000 mL sedimentation cylinder. Deionized water was added to bring the total volume in the sedimentation cylinder to 1000 mL. To begin sedimentation process, the sedimentation cylinder containing the sample was agitated by shaking the cylinder back and forth for a minimum of 30 s. During the agitation, care was taken to ensure that particles were not stuck on the sedimentation cylinder. After agitation the cylinder was placed on the countertop, signifying time zero. Readings were then taken at elapsed times of 0.667, 3, 10, 30, 90, 120, and 720 min using a 152 H hydrometer (H-B Instrument Company, Collegeville, Pennsylvania, USA). Temperature of the suspension liquid was recorded simultaneously with hydrometer readings. To calibrate the hydrometer, a blank reading was taken in a solution containing 0.05 g mL–1 100 mL of Sodium HMP solution and 900 mL deionized water, but without soil.
DNA Extraction, Amplification, and Sequencing
For molecular analysis of community DNA, nine soil samples were randomly selected from the soil samples within each of the three strata. DNA was extracted from these soil samples for each of the three strata and composited following extraction and validation (providing three representative samples per land use system). As previously demonstrated in other 16S rRNA studies , , , DNA was extracted from approximately 0.25 g of soil (oven dry basis of field-moist soil) using the Power Soil Extraction Kit (MO BIO Laboratories, Soloana Beach, California) according to the included protocol. Extracted DNA (2 µL) was quantified using Nanodrop ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington DE), and run on 0.8% agarose gel with 0.5 M TBE buffer. The samples were then submitted to Research and Testing Laboratories (Lubbock, TX) for PCR optimization and pyrosequencing analysis. Bacterial tag-encoded FLX amplicon pyrosequencing PCR, massively parallel pyrosequencing and tag design were carried out according to procedure described previously by Dowd et al. , .
Samples were evaluated using Tag-encoded FLX amplicon pyrosequencing (bTEFAP), which has had prior description and utilization by Dowd et al. in characterizing bacterial populations in a variety of studies , , . All DNA samples were diluted to 20 ng/µl. A 20 ng (1 µl) aliquot of each sample DNA was used for a 25 µl PCR reaction with 5 min denature at 95°C, 30 cycles of 94°C 30 sec –52°C 40 sec –70°C 40 sec with a final extension of 70°C for 5 minutes. The 16S universal Eubacterial primers 28 F (5′- GGC GVA CGG GTG AGT AA) and 530 R (5′-CCG CNG CNG CTG GCA CS) Amplicons were mixed in equal volumes and purified using Agencourt Ampure beads (Agencourt Bioscience Corporation, MA, USA). In preparation for FLX sequencing (Roche, Nutley, NJ), the DNA fragments size and concentration were measured by using DNA chips under a Bio-Rad Experion Automated Electrophoresis Station (Bio-Rad Laboratories, Hercules, CA) and a TBS-380 Fluorometer (Promega Corporation, Madison, WI). A 9.6×106 sample of double-stranded DNA molecules/µl with an average size of 625 bp were combined with 9.6 million DNA capture beads, and then amplified by emulsion PCR. After bead recovery and bead enrichment, the bead-attached DNAs were denatured with NaOH, and sequencing primers were annealed. The 454 Titanium sequencing run was performed on a 70×75 GS PicoTiterPlate by using a Genome Sequencer FLX System (Roche, Nutley, NJ).
Bioinformatics and Statistical Analysis
Quality trimmed sequences were provided with the sequencing services by Research and Testing Laboratories (Lubbock, TX) following the process described in Acosta-Martinez et al. . As described in the aforementioned studies, each sequence was trimmed to utilize only high quality sequence information; tags were extracted from the FLX generated multi-FASTA file, while being parsed into individual sample specific files based upon the tag sequence. Tags which did not have 100% homology to the original sample tag designation were not considered. Sequences which were less than 250 bp after quality trimming were not considered. The B2C2 software , which is described and freely available from Research and Testing Laboratory (Lubbock, TX, USA), was used to deplete samples of definite chimeras. Further processing and OTU based analyses were then carried out using the MOTHUR v.1.19.4  suite of algorithms for sequence processing and diversity analysis, including commands for identifying/consolidating unique sequences, filtering, multiple sequence alignment, generating distance matrices, and clustering of sequences into OTUs. The resulting clusters were assessed at 3% and 5% dissimilarity to provide the data needed for diversity analysis. Based upon the literature we can expect that 0% dissimilarity in sequences will provide dramatic overestimation of the species present in a sample, based upon rarefaction . The resulting sequences were then evaluated using the classify.seqs algorithm (Bayesian method) in MOTHUR against a database derived from the Greengenes set using a bootstrap cutoff of 65%. The sequences contained within the curated 16S database were those considered of high quality based upon Greengenes  standards and which had complete taxonomic information within their annotations. Clusters at 3% and 5% were then utilized to generate rarefaction curves and the (diversity) indices ACE  and Chao1  as well as unweighted and weighted UniFrac for Principle Coordinate Analysis (PCoA) plots and a venn diagram.
All statistical analysis was performed using the SPSS package (SPSS Inc, v 17.0, Chicago, Illinois). The generalized linear model (GLM) was used to assess the means of soil physical, chemical and microbial properties among the systems followed by a Tukey’s HSD test for pairwise comparisons. Relative abundance data is presented as percentages/proportions, but prior to subjection to GLM, they were transformed using the arcsine function for normal distribution prior to analysis. NCSS package (NCSS, 2007, v 7.1.2, Kaysville, Utah) was used for cluster analysis through which double dendrograms were generated through use of the Manhattan distance method with no scaling, and the unweighted pair technique.
Raw sequences were submitted to the NCBI Sequence Read Archive (SRA) and can be found under the accession number SRA007616.
The authors would like to thank Drs. Conrad Bonsi, Desmond Mortley, and Berhanu Tameru for their suggestions and edits to the manuscript.
Conceived and designed the experiments: RS ROA RZ. Performed the experiments: RS SED. Analyzed the data: RS. Contributed reagents/materials/analysis tools: ROA SED. Wrote the paper: RS AMI.
- 1. Parfitt RL, Scott NA, Ross DJ, Salt GJ, Tate KR (2003) Landuse change effects on soil C and N transformations in soils of high N status: comparisons under indigenous forest, pasture and pine plantation, Biogeochemistry 66: 203–221.RL ParfittNA ScottDJ RossGJ SaltKR Tate2003Landuse change effects on soil C and N transformations in soils of high N status: comparisons under indigenous forest, pasture and pine plantation, Biogeochemistry66203221
- 2. Rillig MC, Ramsey PW, Morris S, Paul EA (2003) Glomalin, an arbuscular-mycorrhizal fungal soil protein, responds to land-use change. Plant Soil 253: 293–299.MC RilligPW RamseyS. MorrisEA Paul2003Glomalin, an arbuscular-mycorrhizal fungal soil protein, responds to land-use change.Plant Soil253293299
- 3. Houghton RA, Goodale C (2004) Effects of land-use change on the carbon balance of terrestrial ecosystems. In: DeFries R, Asner G, Houghton RA, editors. pp. 85–98.RA HoughtonC. Goodale2004Effects of land-use change on the carbon balance of terrestrial ecosystems.R. DeFriesG. AsnerRA Houghton8598editors.Ecosystems and Land Use Change.Geophysical Monograph Series.Vol. 153.Washington, DC: American Geophysical Union. Washington, DC: American Geophysical Union.
- 4. Dupouey JL, Dambrine E, Laffite JD, Moares C (2002) Irreversible impact of past land use on forest soils and biodiversity. Ecology 83: 2978–84.JL DupoueyE. DambrineJD LaffiteC. Moares2002Irreversible impact of past land use on forest soils and biodiversity.Ecology83297884
- 5. Honnay O, Piessens K, Van Landuyt W, Hermy M, Gulinck H (2003) Satellite based land use and landscape complexity indices as predictors for regional plant species diversity. Landscape Urban Plan 63: 241–250.O. HonnayK. PiessensW. Van LanduytM. HermyH. Gulinck2003Satellite based land use and landscape complexity indices as predictors for regional plant species diversity.Landscape Urban Plan63241250
- 6. Fu B, Chen L, Ma K, Zhou H, Wang J (2000) The relationships between land use and soil conditions in the hilly area of the loess plateau in northern Shaanxi, China. Catena 39: 69–78.B. FuL. ChenK. MaH. ZhouJ. Wang2000The relationships between land use and soil conditions in the hilly area of the loess plateau in northern Shaanxi, China.Catena396978
- 7. Kennedy AC, Stubbs TL (2006) Soil microbial communities as indicators of soil health. Ann Arid Zone 45: 287–308.AC KennedyTL Stubbs2006Soil microbial communities as indicators of soil health.Ann Arid Zone45287308
- 8. Ibekwe AM, Kennedy AC, Frohne PS, Papiernik SK, Yang CH, et al. (2002) Microbial diversity along a transect of agronomic zones. FEMS Microbiol Ecol 39: 183–191.AM IbekweAC KennedyPS FrohneSK PapiernikCH Yang2002Microbial diversity along a transect of agronomic zones.FEMS Microbiol Ecol39183191
- 9. Lauber CL, Strickland MS, Bradford MA, Fierer N (2008) The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol Biochem 40: 2407–2415.CL LauberMS StricklandMA BradfordN. Fierer2008The influence of soil properties on the structure of bacterial and fungal communities across land-use types.Soil Biol Biochem4024072415
- 10. Roesch LF, Fulthrope RR, Riva A, Casella G, Hadwin AKM, et al. (2007) Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J 1: 283–290.LF RoeschRR FulthropeA. RivaG. CasellaAKM Hadwin2007Pyrosequencing enumerates and contrasts soil microbial diversity.ISME J1283290
- 11. Acosta-Martinez V, Dowd S, Sun Y, Allen V (2008) Tag-encoded pyrosequencing analysis of bacterial diversity in a single soil type as affected by management and land use. Soil Biol Biochem 40: 2762–2770.V. Acosta-MartinezS. DowdY. SunV. Allen2008Tag-encoded pyrosequencing analysis of bacterial diversity in a single soil type as affected by management and land use.Soil Biol Biochem4027622770
- 12. Lauber CL, Hamady M, Knight R, Fierer N (2009) Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microbiol 75: 5111–5120.CL LauberM. HamadyR. KnightN. Fierer2009Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale.Appl Environ Microbiol7551115120
- 13. Brye KR, Slaton NA, Savin MC, Norman RJ, Miller DM (2003) Short-term effects of land leveling on soil physical properties and microbial biomass. Soil Sci Soc Am J 67: 1405–1417.KR BryeNA SlatonMC SavinRJ NormanDM Miller2003Short-term effects of land leveling on soil physical properties and microbial biomass.Soil Sci Soc Am J6714051417
- 14. Hendrix PH, Crossley Jr DA, Coleman DC (1990) Soil biota as components of sustainable agro-ecosystem. In: Edwards CA, Lal R, Madden P, Miller RH, House G, editors. pp. 637–654.PH HendrixDA Crossley JrDC Coleman1990Soil biota as components of sustainable agro-ecosystem.CA EdwardsR. LalP. MaddenRH MillerG. House637654editors.Sustainable Agricultural Systems.Iowa: Soil Water Conservation Society. Iowa: Soil Water Conservation Society.
- 15. Acosta-Martinez V, Dowd SE, Bell CW, Lascano R, Booker JD, et al. (2010a) Microbial Community Composition as Affected by Dryland Cropping Systems and Tillage in a Semiarid Sandy Soil. Diversity 2: 910–931.V. Acosta-MartinezSE DowdCW BellR. LascanoJD Booker2010aMicrobial Community Composition as Affected by Dryland Cropping Systems and Tillage in a Semiarid Sandy Soil.Diversity2910931
- 16. Torsvik V, Daae FL, Sandaa RA, Ovreas L (1998) Novel techniques for analyzing microbial diversity in natural and perturbed environments. J Biotechnol 64: 53–62.V. TorsvikFL DaaeRA SandaaL. Ovreas1998Novel techniques for analyzing microbial diversity in natural and perturbed environments.J Biotechnol645362
- 17. Parham JA, Deng SP, Da HN, Sun HY, Raun WR (2003) Long-term cattle manure application in soil. II. Effect of soil microbial populations and community structure. Biol Fert Soils 38: 209–215.JA ParhamSP DengHN DaHY SunWR Raun2003Long-term cattle manure application in soil. II. Effect of soil microbial populations and community structure.Biol Fert Soils38209215
- 18. Acosta-Martinez V, Harmel RD (2006) Soil microbial communities and enzyme activities under various poultry litter application rates. J Environ Qual 35: 1309–1318.V. Acosta-MartinezRD Harmel2006Soil microbial communities and enzyme activities under various poultry litter application rates.J Environ Qual3513091318
- 19. Jacobson MG (2002) Ecosystem management in the southeast United States: Interest of forest landowners in joint management across ownerships. Small-scale Forest Economics, Management and Policy 1: 71–92.MG Jacobson2002Ecosystem management in the southeast United States: Interest of forest landowners in joint management across ownerships.Small-scale Forest Economics, Management and Policy17192
- 20. Lugo AE (1997) The apparent paradox of reestablishing species richness on degraded lands with tree monocultures. Forest Ecol Manag 99: 9–19.AE Lugo1997The apparent paradox of reestablishing species richness on degraded lands with tree monocultures.Forest Ecol Manag99919
- 21. Tsiouvaras CN, Havlik NA, Bartolome JW (1989) Effects of Goats on Understory Vegetation and Fire Hazard Reduction in a Coastal Forest in California. Forest Sci 35: 1125–1131.CN TsiouvarasNA HavlikJW Bartolome1989Effects of Goats on Understory Vegetation and Fire Hazard Reduction in a Coastal Forest in California.Forest Sci3511251131
- 22. Hart SP (2001) Recent perspectives in using goats for vegetation management in the USA. J Dairy Sci 84(E. Suppl.). pp. E170–E176.SP Hart2001Recent perspectives in using goats for vegetation management in the USA. J Dairy Sci 84(E.Suppl.)E170E176
- 23. Hobbs RJ, Hueneke LF (1992) Disturbance, diversity, and invasion: Implications for conservation. Conserv Biol 6: 324–337.RJ HobbsLF Hueneke1992Disturbance, diversity, and invasion: Implications for conservation.Conserv Biol6324337
- 24. Kohler F, Hamelin J, Gillet F, Gobat J-M, Buttler A (2005) Soil microbial community changes in wooded mountain pastures due to simulated effects of cattle grazing. Plant Soil 278: 327–40.F. KohlerJ. HamelinF. GilletJ-M GobatA. Buttler2005Soil microbial community changes in wooded mountain pastures due to simulated effects of cattle grazing.Plant Soil27832740
- 25. Patra AK, Abbadie L, Clays-Josserand A, et al. (2005) Effects of grazing on microbial functional groups involved in soil N dynamics. Ecol Monogr 75: 65–80.AK PatraL. AbbadieA. Clays-Josserand2005Effects of grazing on microbial functional groups involved in soil N dynamics.Ecol Monogr756580
- 26. Sorensen LI, Mikola J, Kytoviita M-M, Olofsson J (2009) Trampling and spatial heterogeneity explain decomposer abundances in a sib-arctic grassland subjected to simulated reindeer grazing. Ecosystems 12: 830–842.LI SorensenJ. MikolaM-M KytoviitaJ. Olofsson2009Trampling and spatial heterogeneity explain decomposer abundances in a sib-arctic grassland subjected to simulated reindeer grazing.Ecosystems12830842
- 27. Jangid K, Williams MA, Franzluebbers AJ, Sanderlin JS, Reeves JH, et al. (2008) Relative impacts of land-use, management intensity and fertilization upon soil microbial community structure in agricultural systems. Soil Biol Biochem 40: 2843–2853.K. JangidMA WilliamsAJ FranzluebbersJS SanderlinJH Reeves2008Relative impacts of land-use, management intensity and fertilization upon soil microbial community structure in agricultural systems.Soil Biol Biochem4028432853
- 28. Acosta-Martinez V, Dowd SE, Sun Y, Webster D, Allen VG (2010b) Pyrosequencing analysis for characterization of soil bacterial populations as affected by an integrated livestock-cotton production system. Appl Soil Ecol 45: 13–25.V. Acosta-MartinezSE DowdY. SunD. WebsterVG Allen2010bPyrosequencing analysis for characterization of soil bacterial populations as affected by an integrated livestock-cotton production system.Appl Soil Ecol451325
- 29. Karchesy JJ, Hemingway RW (1980) Loblolly pine bark flavonoids. J Agr Food Chem 28: 222–228.JJ KarchesyRW Hemingway1980Loblolly pine bark flavonoids.J Agr Food Chem28222228
- 30. Fierer N, Jackson RB (2006) The diversity and biogeography of soil bacterial communities. P Natl Acad Sci USA 103: 626–631.N. FiererRB Jackson2006The diversity and biogeography of soil bacterial communities.P Natl Acad Sci USA103626631
- 31. Nacke H, Thurmer A, Wollherr A, Will C, Hodac L, et al. (2011) Pyrosequencing-based assessment of bacterial community structure along different management types in German forest and grassland soils. PLoS ONE 6(2): e17000.H. NackeA. ThurmerA. WollherrC. WillL. Hodac2011Pyrosequencing-based assessment of bacterial community structure along different management types in German forest and grassland soils.PLoS ONE6(2)e17000
- 32. Kraus TEC, Dahlgren RA, Zasoski RJ (2003) Tannins in nutrient dynamics of forest ecosystems – a review. Plant Soil 256: 41–66.TEC KrausRA DahlgrenRJ Zasoski2003Tannins in nutrient dynamics of forest ecosystems – a review.Plant Soil2564166
- 33. Selvakumar G, Saha S, Kundu S (2007) Inhibitory activity of pine needle tannin extracts on some agriculturally resourceful microbes. Indian J Microbiol 47: 267–270.G. SelvakumarS. SahaS. Kundu2007Inhibitory activity of pine needle tannin extracts on some agriculturally resourceful microbes.Indian J Microbiol47267270
- 34. Buckley DH, Schmidt TM (2001) The structure of microbial communities in soil and the lasting impact of cultivation. Microb Ecol 42: 11–21.DH BuckleyTM Schmidt2001The structure of microbial communities in soil and the lasting impact of cultivation.Microb Ecol421121
- 35. Osborne CA, Zwart AB, Broadhurst LM, Young AG, Richardson AE (2011) The influence of sampling strategies and spatial variation on the detected soil bacterial communities under three different land use types. FEMS Microbiol Ecol 78: 70–79.CA OsborneAB ZwartLM BroadhurstAG YoungAE Richardson2011The influence of sampling strategies and spatial variation on the detected soil bacterial communities under three different land use types.FEMS Microbiol Ecol787079
- 36. Neher DA (1999) Soil community composition and ecosystem processes: Comparing agricultural ecosystems with natural ecosystems. Agroforest Syst 45: 159–185.DA Neher1999Soil community composition and ecosystem processes: Comparing agricultural ecosystems with natural ecosystems.Agroforest Syst45159185
- 37. Fierer N, Bradford MA, Jackson RB (2007) Toward an ecological classification of soil bacteria. Ecology 88: 1354–1364.N. FiererMA BradfordRB Jackson2007Toward an ecological classification of soil bacteria.Ecology8813541364
- 38. Singh BK, Bardgett RD, Smith P, Reay DS (2010) Microorganisms and climate change: terrestrial feedbacks and mitigation options. Nat Rev Microbiol 8: 779–790.BK SinghRD BardgettP. SmithDS Reay2010Microorganisms and climate change: terrestrial feedbacks and mitigation options.Nat Rev Microbiol8779790
- 39. Wardle DA, Yeates GW, Watson RN, Nicholson KS (1995) The detritus food-web and the diversity of soil fauna as indicators of disturbance regimes in agro-ecosystems. Plant Soil 170: 35–43.DA WardleGW YeatesRN WatsonKS Nicholson1995The detritus food-web and the diversity of soil fauna as indicators of disturbance regimes in agro-ecosystems.Plant Soil1703543
- 40. Hill P, Kristufek V, Dijfkhuizen L, Boddy C, Kroetsch D, et al. (2011) Land use intensity controls actinobacterial community structure. Microb Ecol 61: 286–302.P. HillV. KristufekL. DijfkhuizenC. BoddyD. Kroetsch2011Land use intensity controls actinobacterial community structure.Microb Ecol61286302
- 41. Burke RA, Molina M, Cox JE, Osher LJ, Piccolo MC (2003) Stable carbon isotope ratio and composition of microbial fatty acids in tropical soils. J Environ Qual 32: 198–206.RA BurkeM. MolinaJE CoxLJ OsherMC Piccolo2003Stable carbon isotope ratio and composition of microbial fatty acids in tropical soils.J Environ Qual32198206
- 42. Waldrop MP, Balser TC, Firestone MK (2000) Linking microbial community composition to function in a tropical soil. Soil Biol Biochem 32: 1837–1846.MP WaldropTC BalserMK Firestone2000Linking microbial community composition to function in a tropical soil.Soil Biol Biochem3218371846
- 43. Ishak HD, Miller JL, Sen R, Dowd SE, Meyer E, et al. (2011) Microbiomes of ant castes implicate new microbial roles in the fungus-growing ant Trachymyrmex septentrionalis. Sci Rep 1: 204.HD IshakJL MillerR. SenSE DowdE. Meyer2011Microbiomes of ant castes implicate new microbial roles in the fungus-growing ant Trachymyrmex septentrionalis.Sci Rep1204
- 44. Lopez-Velasco G, Welbaum GE, Boyer RR, Mane SP, Ponder MA (2011) Changes in spinach phylloepiphytic bacteria communities following minimal processing and refrigerated storage described using pyrosequencing of 16S rRNA amplicons. J Appl Microbiol 110: 1203–1214.G. Lopez-VelascoGE WelbaumRR BoyerSP ManeMA Ponder2011Changes in spinach phylloepiphytic bacteria communities following minimal processing and refrigerated storage described using pyrosequencing of 16S rRNA amplicons.J Appl Microbiol11012031214
- 45. Saul-Tcherkas V, Steinberger Y (2011) Soil Microbial Diversity in the Vicinity of a Negev Desert Shrub–Reaumuria negevensis. Microb Ecol 61: 64–81.V. Saul-TcherkasY. Steinberger2011Soil Microbial Diversity in the Vicinity of a Negev Desert Shrub–Reaumuria negevensis.Microb Ecol616481
- 46. Chong CW, Convey P, Pearce DA, Tan IKP (2012) Assessment of soil bacterial communities on Alexander Island (in the maritime and continental Antarctic transitional zone. Polar Biol 35: 387–399.CW ChongP. ConveyDA PearceIKP Tan2012Assessment of soil bacterial communities on Alexander Island (in the maritime and continental Antarctic transitional zone.Polar Biol35387399
- 47. Connell JH (1978) Diversity in tropical rain forests and coral reefs. Science 199: 1302–1310.JH Connell1978Diversity in tropical rain forests and coral reefs.Science19913021310
- 48. Tabatabai MA (1994) Enzymes, In: Weaver RW, Augle S, Bottomly PJ, Bezdicek D, Smith S, Tabatabai MA, Wollum A, editors. Methods of Soil Analysis, Part 2: Microbiological and Biochemical properties. Madison: Soil Science Society of America. pp. 775–833.MA Tabatabai1994Enzymes, In: Weaver RW, Augle S, Bottomly PJ, Bezdicek D, Smith S, Tabatabai MA, Wollum A, editors. Methods of Soil Analysis, Part 2: Microbiological and Biochemical properties.Madison: Soil Science Society of America775833
- 49. Bandick AK, Dick RP (1999) Field management effects on soil enzyme activities. Soil Biol Biochem 31: 1471–1479.AK BandickRP Dick1999Field management effects on soil enzyme activities.Soil Biol Biochem3114711479
- 50. Elsgaard L, Anderson GH, Eriksen J (2002) Measurement of arylsulphatase activity in agricultural soils using a simplified assay. Soil Biol Biochem 34: 79–82.L. ElsgaardGH AndersonJ. Eriksen2002Measurement of arylsulphatase activity in agricultural soils using a simplified assay.Soil Biol Biochem347982
- 51. Mclean EO (1982) Soil pH and lime requirements. In: Page AL, Miller RH, Keeney , DR , editors. pp. 199–224.EO Mclean1982Soil pH and lime requirements.AL PageRH MillerKeeneyDR199224editors.Methods of Soil Analysis, Part 2.Chemical and Microbiological Properties.Madison: Soil Science Society of America. Madison: Soil Science Society of America.
- 52. Keller JM, Gee GW (2006) Comparison of American Society of Testing Materials and Soil Science Society of America hydrometer methods for particle-size analysis. Soil Sci Soc Am J 70: 1094–1100.JM KellerGW Gee2006Comparison of American Society of Testing Materials and Soil Science Society of America hydrometer methods for particle-size analysis.Soil Sci Soc Am J7010941100
- 53. Dowd SE, Sun Y, Secor PR, Rhoads DD, Wolcott BM, et al. (2008a) Survey of bacterial diversity in chronic wounds using pyrosequencing, DGGE, and full ribosome shotgun sequencing. BMC Microbiol 8: 43.SE DowdY. SunPR SecorDD RhoadsBM Wolcott2008aSurvey of bacterial diversity in chronic wounds using pyrosequencing, DGGE, and full ribosome shotgun sequencing.BMC Microbiol843
- 54. Dowd SE, Sun Y, Wolcott RD, Carroll JA (2008b) Bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP) for microbiome studies: bacterial diversity in the ileum of newly weaned Salmonella-infected pigs. Foodborne Pathog Dis 5: 459–472.SE DowdY. SunRD WolcottJA Carroll2008bBacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP) for microbiome studies: bacterial diversity in the ileum of newly weaned Salmonella-infected pigs.Foodborne Pathog Dis5459472
- 55. Feingold SM, Dowd SE, Gontcharova V, Liu C, Henley KE, et al. (2010) Pyrosequencing study of fecal microflora of autistic and control children. Anaerobe 16: 444–453.SM FeingoldSE DowdV. GontcharovaC. LiuKE Henley2010Pyrosequencing study of fecal microflora of autistic and control children.Anaerobe16444453
- 56. Gontcharova V, Youn E, Wolcott RD, Hollister EB, Gentry TJ, et al. (2010) Black Box Chimera Check (B2C2): a windows-based software for batch depletion of chimeras from bacterial 16S rRNA gene datasets. Open Microbiol J. 4: 47–52.V. GontcharovaE. YounRD WolcottEB HollisterTJ Gentry2010Black Box Chimera Check (B2C2): a windows-based software for batch depletion of chimeras from bacterial 16S rRNA gene datasets. Open Microbiol J.44752
- 57. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, et al. (2009) Introducing mothur: open source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75: 7537–7541.PD SchlossSL WestcottT. RyabinJR HallM. Hartmann2009Introducing mothur: open source, platform-independent, community-supported software for describing and comparing microbial communities.Appl Environ Microbiol7575377541
- 58. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, et al. (2006) Greengenes, a Chimera-Checked 16S rRNA Gene Database and Workbench Compatible with ARB. Appl Environ Microbiol 72: 5069–72.TZ DeSantisP. HugenholtzN. LarsenM. RojasEL Brodie2006Greengenes, a Chimera-Checked 16S rRNA Gene Database and Workbench Compatible with ARB.Appl Environ Microbiol72506972
- 59. Chao A, Ma MC, Yang MCK (1993) Stopping rules and estimation for recapture debugging with unequal failure rates. Biometrika 80: 193–201.A. ChaoMC MaMCK Yang1993Stopping rules and estimation for recapture debugging with unequal failure rates.Biometrika80193201
- 60. Chao A (1984) Nonparametric estimation of the number of classes in a population. Scand. J Stat 11: 265–270.A. Chao1984Nonparametric estimation of the number of classes in a population. Scand.J Stat11265270