Soil Rhizosphere Microbial Communities and Enzyme Activities under Organic Farming in Alabama

Evaluation of the soil rhizosphere has been limited by the lack of robust assessments that can explore the vast complex structure and diversity of soil microbial communities. Our objective was to combine fatty acid methyl ester (FAME) and pyrosequencing techniques to evaluate soil microbial community structure and diversity. In addition, we evaluated biogeochemical functionality of the microbial communities via enzymatic activities of nutrient cycling. Samples were taken from a silt loam at 0–10 and 10–20 cm in an organic farm under lettuce (Lactuca sativa), potato (Solanum tuberosum), onion (Allium cepa L), broccoli (Brassica oleracea var. botrytis) and Tall fescue pasture grass (Festuca arundinacea). Several FAMEs (a15:0, i15:0, i15:1, i16:0, a17:0, i17:0, 10Me17:0, cy17:0, 16:1ω5c and 18:1ω9c) varied among the crop rhizospheres. FAME profiles of the soil microbial community under pasture showed a higher fungal:bacterial ratio compared to the soil under lettuce, potato, onion, and broccoli. Soil under potato showed higher sum of fungal FAME indicators compared to broccoli, onion and lettuce. Microbial biomass C and enzyme activities associated with pasture and potato were higher than the other rhizospheres. The lowest soil microbial biomass C and enzyme activities were found under onion. Pyrosequencing revealed significant differences regarding the maximum operational taxonomic units (OTU) at 3% dissimilarity level (roughly corresponding to the bacterial species level) at 0–10 cm (581.7–770.0) compared to 10–20 cm 309 (563.3–727.7) soil depths. The lowest OTUs detected at 0–10 cm were under broccoli (581.7); whereas the lowest OTUs found at 10–20 cm were under potato (563.3). The predominant phyla (85%) in this soil at both depths were Bacteroidetes (i.e., Flavobacteria, Sphingobacteria), and Proteobacteria. Flavobacteriaceae and Xanthomonadaceae were predominant under broccoli. Rhizobiaceae, Hyphomicrobiaceae, and Acidobacteriaceae were more abundant under pasture compared to the cultivated soils under broccoli, potato, onion and lettuce. This study found significant differences in microbial community structure and diversity, and enzyme activities of nutrient cycling in this organic farming system under different rhizospheres, which can have implications in soil health and metabolic functioning, and the yield and nutritional value of each crop.


Introduction
Organic farming obtained recognition since the end of the 20th Century due to increased public awareness of the negative impacts of intensive agricultural practices (i.e., fertilization, tillage, and pesticide application) on the environment and food safety [1][2][3][4][5][6].Organic farming practices attempt to preserve environmental biodiversity by crop diversification and avoiding chemical products using organic based composting techniques for fertilizer (i.e., manure) [7,8].The pioneers of the organic farming movement in 1940 emphasized as its main focus to maintain or improve soil health, the capacity of soil to function as a vital living system within land-use boundaries [9], which is dependent on soil decomposition processes and nutrient cycling [6][7][8] driven by the soil microbial communities.Microbial communities are key drivers in ecosystem functioning and sustainability because they are the main source of necessary enzymes capable of decomposing plant-derived compounds [7].Microbial communities associated to the rhizosphere, defined as the volume of soil adjacent to and influenced by plant root, are important to plant health and soil fertility [10].The rhizosphere of most crops can exert an influence on the soil microbial communities within the first 0.1 m of soil [11].Root exudates from different plants can stimulate the growth of unique bacterial and fungal populations in the vicinity of roots [12].Thus, determining the size and composition of soil microbial communities associated with different rhizospheres under organic farming practices is essential in order to evaluate above-and below-ground agroecosystem health and functioning.
Studies evaluating the microbial communities associated with different rhizospheres should include a combination of different techniques to explore the composition at different levels including molecular techniques [13][14][15][16] and fatty acid profiling such as Phospholipid fatty acid analysis (PLFA) or Fatty acid methyl ester (FAME) methods [17][18][19].Fatty acid profiling can provide information of the overall soil microbial community structure by describing the abundance of major microbial groups with indicator fatty acids for fungal populations (16:1ω5c, 18:3ω6c, 18:1ω9c and 18:2ω6c) and bacterial populations including Gram-positive bacteria (i.e., i15:0, a15:0, i17:0, a17:0), Gram-negative bacteria (i.e., cy17:0, cy19:0, i13:0 3OH and i17:0 3OH) and actinomycetes (10Me16:0, 10Me17:0 and 10Me18:0).Along with fatty acid profiling, more detailed exploration of microbial community composition, to the species level, has been possible with molecular fingerprinting techniques.Recent advances in sequencing technology, such as the development of pyrosequencing, has been described as a promising technique to better characterize microbial diversity in natural environments compared with traditional molecular biological techniques involving cloning [20][21][22].Pyrosequencing coupled with the 454 Genome Sequencer FLX System (Roche, Nutley, New Jersey) provides a fast and massive sequencing approach (i.e., 400-600 MB output per run in a day), which can provide several unique sequences or operational taxonomic units (OTUs; roughly corresponding to the species level) in soils.
In Alabama, organic farming is a newly emerging practice.For example, the first organic farm (Tune Farm, Falkville, AL) was certified in 2003.Thus, information is needed on the microbial communities associated to organic rhizospheres to help producers in the selection of sustainable crop rotations and management practices in this region.We evaluated soil microbial communities under an undisturbed pasture of Tall fescue grass (Festuca arundinacea) and from lettuce (Lactuca sativa), potato (Solanum tuberosum), onion (Allium cepa L), and broccoli (Brassica oleracea var.botrytis), which were previously under potato for the prior two years.The microbial community composition was evaluated using the ester linked (EL)-FAME and pyrosequencing analyses under the different rhizospheres.We also explored the relationship of the FAME indicators of fungal and bacterial populations, and bacterial phyla distribution from pyrosequencing with other soil properties such as microbial biomass C and N, total C and N, and soil biogeochemical functionality assessed using enzyme activities of nutrient cycling.By coupling robust information regarding the vast diversity of the soil microbial community via FAME profiling and pyrosequencing with an assessment of biogeochemical functionality via enzymatic activities of nutrient cycling, we aim to increase the understanding of microbial communities associated to fundamental biogeochemical processes under different rhizospheres.

Description of Experimental Site and Management
The Tune Organic Farm is located in Morgan County, Falkville, Alabama [−86°54'20" (−86.905572)west longitude; 34°20'9" (34.335921) north latitude; elevation].The land at this site is characterized by 0 to 4% slopes.The soil is characterized as silt loam (12.8% clay, 71.8% silt and 15.4% sand; fine-loamy, mixed, semiactive), and it is classified as a mesic Typic Fragiudult.This soil has an average pH of 6.5 and 2.25 total C kg −1 soil.The farm is located in a region with a mean annual precipitation of 1406 mm (55 in) with most of the precipitation occurring from January through May and October through December.The total area evaluated comprises 4.4 ha (11 acres) under irrigation with county water with spigot located on the soil surface and covered with mulch between 25.4 mm (1 inch) and 58.8 mm (2.31 inches).Prior to receiving certification as an organic production system in 2003, Tune Organic Farm had practiced uncertified organic production system since 1988.
Prior to this study, this farm only underwent organic potato and pasture (predominantly Festuca arundinacea grasses) production (2006)(2007)(2008).The soil was conventionally tilled (raised beds) in October 2006, followed by hand weed control.The farm utilized alternative cropping rotation management practices using winter crops as recommended by USDA, such as rye (Secale cereal L.), hairy vetch (Vicia hirsuta) or buckwheat (Fagopyrum esculentum).In 2007, the farm underwent a minimum biannual minimum tillage (0-15 cm) occurring every November after the winter crop harvest; disking and leveling with a field cultivator in April before summer crop planting.Mulch tillage tilling to depths of 25 and 51 mm (~1 and 2 inches) was incorporated using a rotary field cultivator before planting.Potato (Yukon Gold cv.), lettuce (Pirat Butterhead cv.), broccoli (Green Sprouting Calabrese cv.) and onion (Mustang Yellow cv.) were planted in early March.

Soil Sampling
Soil surface samples were taken within the first week of May 2009 under potato, lettuce, broccoli, onion and an undisturbed pasture.Except for pasture, this sampling corresponds to 9 weeks after planting.For each crop, three sampling locations were randomly selected.A composite soil sample was taken at each sampling location within the rhizosphere only (samples were not taken within furrows) by combining three soil cores (2.5 cm dia.) and then split into 0-10 and 10-20 cm depths.

Microbial Community Structure and Diversity According to FAME Analysis
The ester linked (EL)-FAME method was performed as described by Schutter and Dick [23].The following 4 steps were performed: (1) saponification and methylation of ester-linked fatty acids by incubation of 3 g of soil in 15 mL of 0.2 M KOH in methanol at 37 °C for 1h.During that time, the samples are vortexed every 10 min, and an addition of 3 mL of 1.0 M acetic acid to neutralize the pH of the mixture at the end of incubation; (2) Partitioning of the FAMEs into an organic phase by adding 10 mL of hexane followed by centrifugation at 480 × g for 10 min; (3) Transfer of the hexane layer to a clean glass test tube and evaporation of the hexane under a stream of N 2 ; and (4) Dissolve FAMEs in 0.5 mL of 1:1 hexane:methyl-tert butyl ether and transferred to a GC vial for analysis.Using this protocol, the extracted FAMEs were analyzed in a 6890 GC Series II (Hewlett Packard, Wilmington, DE, USA) equipped with a flame ionization detector and a fused silica capillary column (25 m × 0.2 mm) using H 2 (ultra high purity) as the carrier gas.The temperature program was ramped from 170 °C to 250 °C at 5 °C min −1 as previously done by Acosta-Martinez et al. [21].Fatty acids were identified and quantified by comparison of retention times and peak areas to components of MIDI standards.FAME concentrations (nmol g −1 soil) were calculated by comparing peak areas to an analytical standard (19:0, Sigma Chemical Co., St. Louis, MO, USA) calibration curve.The FAMEs are described by the number of C atoms, followed by a colon, the number of double bonds and then by the position of the first double bond from the methyl (ω) end of the molecule.Cis isomers are indicated by c, and branched fatty acids are indicated by the prefixes i and a for iso and anteiso, respectively.Other notations are Me for methyl, OH for hydroxyl and cy for cyclopropane.

Bacterial Community Structure and Diversity According to Pyrosequencing Analysis
DNA was extracted from approximately 0.5 g of soil (oven dry basis of field-moist soil) using the Fast DNA Spin Kit for soil (QBIOgene, Carlsbad, CA, USA) following the manufacturer's instructions.The DNA extracted (1 µL) was quantified using Nanodrop ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA).The integrity of the DNA extracted from the soils was confirmed by running DNA extracts on 0.8% agarose gel with 0.5 × TBE buffer (45 mM Tris-borate, 1 mM EDTA, pH 8.0).The 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing technique was used as described by Dowd et al. [24] and Acosta-Martínez et al. [21,22] for characterization of soil bacterial community composition.We used the new bacterial tag-encoded FLX-Titanium amplicon pyrosequencing (bTEFAP) approach, which is based upon the same original principles of the method described by Acosta-Martínez et al. [21,22], but instead utilizing the Titanium reagent procedures [25] and a one-step PCR rather than a two-step labeling reaction, a mixture of Hot Start and HotStar high fidelity taq polymerases, and amplicons originating from the 27F region numbered in relation to E. coli rRNA.Quality trimmed sequences (~Q20) averaging 410 bp with range of 250-520 bp obtained from the FLX sequencing run were processed using a custom scripted bioinformatics pipeline (Research and Testing Laboratory, Lubbock, TX).Sequences which were less than 250 bp after quality trimming were not considered.Samples were depleted of definite chimeras using B2C2 software, which is described and freely available from Research and Testing Laboratory (Lubbock, TX, USA; www.researchandtesting.com/B2C2.html)[26].The resulting sequences averaging 4216 sequences per sample were then evaluated using BLASTn [27] against a monthly updated custom 16S database derived from GenBank [28].Phylogenetic assignments were based upon NCBI taxonomic designations.

Additional Soil Analyses
Microbial biomass carbon (MBC) and nitrogen (MBN) were determined on a 15 g oven-dry equivalents field-moist soil sample (<2 mm) in 0.5 M K 2 SO 4 (1:5 w/v) after chloroform fumigation [29,30].Total organic C and N concentrations in the extracts from the fumigated (24 h) and non-fumigated (control) soil were quantified by a CN analyzer (Shimadzu Model TOC-V/ CPH -TN, Japan).The non-fumigated control values were subtracted from the fumigated values.The MBC and MBN were calculated using a k EC factor of 0.45 [31] and k EN factor of 0.54 [32], respectively.Each sample had duplicate analyses and results are expressed on a moisture-free basis.Moisture was determined after drying at 105 °C for 48 h.
A subsample of the soil was air-dried (<2 mm) for selected analyses.For example, soil pH was measured using a 1:2.5 air-dried soil to water solution ratio using an Orion conductivity meter (model-160).Enzyme activities important for C (β-glucosidase), C and N (β-glucosaminidase), P (acid phosphatase) and S (arylsulfatase) cycling were evaluated using 1 g of air-dried soil with their appropriate substrate and incubated for 1 h (37 °C) at their optimal pH as described previously [33,34].Soil texture was assessed using the pipette method [35].Total C, N, and S were determined in subsamples of air-dried soil, sieved to <180 µm, by automated dry combustion method using a vario Max CNS analyzer.

Statistical Analysis
To compare FAME indicators under different rhizospheres in this organic production system on a multivariate scale, we used principal components analysis (PCA) using PC-ORD software (version 5, MjM Software Design).In the PCA, square root transformations were performed on all FAMEs to create a more normally distributed data set and to reduce the coefficient of variation among FAMEs.The square root transformation is similar to log transformations, which is commonly used in ecological studies [36].The PCAs were performed using cross-products matrix with variance/covariance centered and calculating scores for FAMEs by weighted averaging.In addition, a non-metric multidimensional scaling (NMS) plot, non-parametric method, was used to provide graphical ordination of the FAME indicators, microbial biomass and enzyme activities under different rhizospheres with two synthetic axes.Data were transformed by the different rhizospheres using the 'general relativization' procedure in PC-ORD [36].'Relativization' in multivariate analysis removes the potentially strong influence of absolute abundance on community data.Indicator species' analysis and correlations with the main matrix are presented when appropriate to facilitate the identification of variables that are important descriptors of differences among soil microbial community (FAME) structure, microbial biomass and biogeochemical cycling [37,38].
A double dendogram was performed using comparative functions and multivariate hierarchical clustering methods of NCSS 2007 (NCSS, Kaysville, Utah) based upon the top 27 most abundant bacterial classes including the weighted-pair group clustering method and the Manhattan distance method with no scaling.For bacterial diversity and richness analysis, multiple sequence alignment was performed using MUSCLE (with parameter -maxiters 1, -diags1 and -sv) [39].Based on the alignment, a distance matrix was constructed using DNAdist from PHYLIP version 3.6 with default parameters from Felsenstein [40,41].These pairwise distances served as input to DOTUR [42] for clustering the sequences into OTUs of defined sequence similarity that ranged from 0% to 20% dissimilarity.At 3% dissimilarity, we expect that we can obtain accurate estimation of the majority of species present in a given sample.The clusters based upon dissimilarity of 3%, served as OTUs for generating predictive rarefaction models and for making calculations with the richness (diversity) indexes Ace and Chao1 [43] in DOTUR.The heat map colors indicated a variation on the relative percentage of bacteria in this soil that ranged from <0.00% in black up to 85% in red.Manhattan distance for the relative percentage data is calculated between rows j and k, where: and P is the percentage of similarity between samples:

Selected Soil Properties
At 0-10 cm, total C ranged in this soil from 2.03 to 3.00 g kg −1 soil while total N ranged from 0.17 to 0.28 g kg −1 soil (Table 1).Soil pH (~6.80) was significantly higher under broccoli, onion and lettuce compared to pasture and potato (P < 0.05).Soil associated with the onion rhizosphere exhibited lower levels of MBC and MBN compared to the other rhizospheres.Enzyme activities varied across the different rhizospheres, but trends differed depending on the enzyme.For example, β-glucosidase activity was higher under potato (264 mg PN kg −1 soil h −1 ) compared to the other rhizospheres (121-159 mg PN kg −1 soil h −1 ), whereas β-glucosaminidase activity was highest in soil under pasture, broccoli and potato (54.50-64 mg PN kg −1 soil h −1 ) compared to onion and lettuce (28.10-35.30mg PN kg −1 soil h −1 ).Acid phosphatase activity was higher in soil under broccoli (534 mg PN kg −1 soil h −1 ), followed by pasture and potato (average: 443 mg PN kg −1 soil h −1 ), and the lowest activities were found under onion and lettuce (average: 362.5 mg PN kg −1 soil h −1 ).Arylsulfatase activity was more similar among the rhizospheres (144-227 mg PN kg −1 soil h −1 ) but was lowest under onion.In summary, soil properties (overall) exhibited higher levels at 0-10 cm compared to 10-20 cm in this organic production system, except for soil pH.Values in parenthesis () are the standard error of the mean (n = 3).

Microbial Community Structure using FAME Analyses
Using PCA to compare the FAME profiles within the rhizospheres at 0-10 cm, we found significant overall differences in the microbial community structure among the rhizospheres within the first PC axis, which accounted for the greatest variation among these samples (90.7%) (Figure 1).A significant correlation was found between the Eigenvector scores of soil under potato, broccoli, and pasture, which were significantly different (P < 0.007) from soil under onion.The scores obtained for the FAME profiles under broccoli were not significantly different from those obtained under lettuce, and the same applies for the scores obtained under lettuce compared to onion.The trends observed in the PCA were due to distinct differences in the soil fungal FAME indicators among the rhizospheres (Table 2).This trend was found for the fungal FAMEs: pasture > potato > broccoli > lettuce > onion.Soil bacterial FAME indicators, Gram+, Gram-, and actinomycetes, were higher under potato, followed by pasture, broccoli, and lettuce, and the lowest under onion.The only FAME indicator for Gram-populations detected in this soil was cy17:0, which ranged across all rhizospheres from 36.93 to 70.53 nmols g −1 soil.Pasture and broccoli were not significantly different from one another for any bacterial FAME indicator.Soil under potato showed significantly greater abundance of actinomycetes (32.8 nmols g −1 soil) indicators compared to onion and lettuce (14.87 and 20.00 nmols g −1 soil, respectively).Overall, our results indicate that there are more significant differences in the fungal community (sum ranged 345-854) than bacterial community (sum ranged 417-779) among these rhizospheres.Values in parenthesis () are the standard error of the mean (n = 3).
Our non-metric multidimensional scaling (NMDS) plot revealed correlations among the FAME indicators for bacterial and fungal populations, microbial biomass and enzyme activities under the rhizospheres evaluated (Figure 2).The first axis explained 89% of the variability among this data set.The NMDS plot displayed soil fungal FAMEs (i.e., F13) positively correlated to soil under pasture while β-glucosidase activity (E2) and FAME indicators for Gram− populations (i.e., F7) were most correlated to soil under onion.FAME indicators for actinomycetes (i.e., F8), some bacterial FAMEs (F5 and F8) and β-glucosidase activity (E2) were correlated to soil under lettuce.Gram+ bacterial FAMEs were correlated to soil under potato.Soil under broccoli grouped near the centroid of the plot, and demonstrated little to no correlation with all soil properties relative to the other rhizospheres.As a result of these trends it seems obvious that soil fungi demonstrate either no relationship or negative relationships with all the other rhizospheres except pasture.Furthermore, lettuce and onion demonstrate similar soil properties, as both are characterized as having distinctly different soil microbial properties from pasture or potato rhizospheres.

Figure 2.
Non-metric multidimensional scaling (NMDS) plot of the relationship between soil microbial community composition (according to 13 FAMEs), enzyme activities and microbial biomass under different rhizospheres in an organic production system.The figure shows FAMEs for bacterial populations according to F1-5 for Gm+ (Gram+), F6-7 for Gm− (Gram−) and F8-10 for actinomycetes (Actino) while the fungal populations are represented by F11-13.Fungal:bacterial ratios are shown as F:B.The enzyme activities are represented (Enz) with E labels such as E1 for arylsulfatase, E2 for β-glucosidase, E3 for β-glucosaminidase and E4 acid phosphatase.Microbial biomass is shown as MB (MC for microbial biomass C or MN for microbial biomass N).
The bacterial diversity indexes, ACE and Chao1, at 3% dissimilarity were higher at 0-10 cm than at 10-20 cm under pasture, potato, and lettuce (Table 4).In addition, these indexes tended to be higher under potato than the other rhizospheres at 0-10 cm.At the lower depth of 10-20 cm, soil under onion showed the highest levels of these indexes among these rhizospheres.Table 3. Top predominant bacterial families under different rhizospheres in an organic production system (0-10 cm).In bold are the bacteria lower in pasture and those with (*) represent higher bacterial abundance in the rhizosphere specified.

Microbial Community Composition According to the FAME Analyses, and Linkage to Enzyme Activity Trends
Our findings demonstrate distinct contrasts in both soil microbial community characteristics and enzyme activities across the rhizospheres evaluated.The first distinct difference found in the microbial community structure was the highest fungal:bacterial ratio under pasture compared to the other rhizospheres.These results are consistent with other studies reporting higher fungal populations, including mycorrhiza, under non-disturbed soils compared to cultivated soils [21,44].Soil under potato also showed higher fungal populations compared to the rhizospheres of lettuce, broccoli, and onion, which may be explained by the lower soil pH (close to 6) under potato.Blagodatskaya and Anderson [45] reported that fungi are found in higher abundance in acidic soils than alkaline soils, and that soil pH influenced the percent fungal-to-bacterial ratio to a greater degree than substrate in their study.This farm had been under potato production for the two years prior to establishing different crops, which may have caused a long-lasting accumulation of certain substrates or fungal spores involved in arbuscular mycorrhizae (AM) associations and/or decreased the soil pH enough to encourage higher fungal populations under potato compared to the other crops.The findings of higher fungal populations under pasture and potato compared to lettuce, onion, and broccoli can have an impact on soil processes.Previous studies have reported that soil microbial communities with higher fungal: bacterial ratios have been associated to potential benefits in C sequestration as fungal populations tend to have higher C assimilation efficiencies than bacterial populations and store higher amounts of the C they metabolize [46].
The different rhizospheres showed differences in the group of enzyme activities evaluated, which may reveal variations in C, N, P, and S biogeochemical cycling among these rhizospheres.The fact that the potato rhizosphere demonstrated the highest fungal FAMEs, but also highest activities of β-glucosidase and β-glucosaminidase are of ecological significance to C and N cycling, respectively [46,47].However, soil under broccoli demonstrated the highest acid phosphatase activity, followed by pasture, which have important implications of higher P transformation and turnover under broccoli and pasture.Higher levels of acid phosphatase activity have been shown in other pasture soils by Speir and Ross [48] but no study highlighted higher levels of this enzyme activity under broccoli compared to other vegetables.Various studies in respect to relationships between phosphatase activity and mineralization of organic P in the rhizosphere may be attributed to the complex interaction of many factors including plant species, soil properties (form and level of P, pH, soil moisture), and time (different stages of plant growth) [49,50].
Another distinctive trend in this organically managed system was the lowest levels of FAME indicators for Gram+ and Gram− bacterial populations and fungal populations detected under onion compared to the other rhizospheres.Total FAMEs and microbial biomass (MBC and MBN), indicators for microbial population size, were lower under onion.As microbial communities are the main source of enzymes in soil, the lower microbial biomass under onion may explain the lower enzyme activities also found [14].These findings suggested lower biogeochemical cycling potential under onion compared to broccoli, lettuce, potato, and pasture.

Soil Bacterial Communities According to Pyrosequencing and Linkage to Enzyme Activity Trends
Assessing soil bacterial diversity using pyrosequencing promises to expand our understanding of the distribution of major bacterial groups in soil, which contribute to essential soil processes including soil organic matter incorporation, aggregation and nutrient cycling [21,51].However, we are only beginning to explore how plants influence the make-up of their own rhizosphere microbial communities [52], and how it may be linked to soil functioning.Pyrosequencing allowed us to explore up to 6,000 sequences and revealed between 494-868 unique operational taxonomic units (OTU; roughly corresponding to the species level) at the 3% dissimilarity level for this soil under organic production.This technique allowed us to explore almost 7 times higher resolution (number of sequences) than previously possible [20,53].In agreement with other soil microbial ecology research efforts, this study detected higher overall abundances of the major bacterial phyla in soil: Acidobacteria, Actinobacteria, Proteobacteria, Verrucomicrobia, Bacteroidetes, Firmicutes, Chloroflexi, Planctomycetes and Gemmatimonadetes [54,55].However, our study provided data demonstrating the varying distribution among the major bacterial phyla across soils with different rhizospheric characteristics.For example, other studies have shown a high predominance of Acidobacteria among other bacterial phyla [56].In this soil, the three most predominant bacteria have been linked to soil organic matter dynamics, which showed this trend: Bacteroidetes (Flavobacteria and Sphingobacteria) > Proteobacteria > Actinobacteria.Previous studies have emphasized that Flavobacteriaceae are recognized for their role in soil metabolic functioning [57], and it is possible that they could provide a strong contribution to the mineralization of primary-produced organic matter in this organic soil [58].The higher percentage of Proteobacteria in soils is well known and it has been reported in several other studies that this bacteria encompasses an enormous level of morphological, physiological and metabolic diversity that could play a significant role in global C, N, and S cycling [21,59].Actinobacteria, also found predominant in this organic soil, has been recognized as important in organic matter dynamics [22].
Previous studies have indicated that root exudates and plant species strongly determine the bacterial composition in the rhizosphere, producing plant genotype-specific community structures in the same soil [60][61][62][63].In addition to the fact that pasture demonstrated the highest fungal:bacterial ratios compared to the other rhizospheres, we found distinct differences in bacterial communities in pasture compared to the other rhizospheres.Rhizobiaceae, Hyphomicrobiaceae, Bradyrhizobiaceae, Acetobacteriaceae, Conexibacteriaceae, and Burkholderiacea were up to 2 times higher under pasture compared to the other rhizospheres.These Proteobacteria (Bradyrhizobiaceae, Hyphomicrobiaceae, and Acetobacteriaceae) represent key bacteria involved in soil C and N cycling [21,59].The predominance of many bacteria under pasture, including Acidobacteriaceae, may reflect the combined effects of very different root systems and their exudates (i.e., type and amounts), and differences in soil properties (i.e., lower soil pH, less disturbance) compared to the other rhizospheres.
Soil under potato revealed higher abundance of bacteria belonging to the β-Proteobacteria (class) (i.e., Rhodocyclaceae, Comamonadaceae, and Oxalobacteraceae), which has been found to have a significant role in C cycling.Other studies have linked β-Proteobacteria to the degradation of halogenated compounds [64], biological phosphate removal [65], denitrification and reduction of (per) chlorates [66,67], and N fixation [68].The overall findings of a different microbial community structure (i.e., higher fungal FAMEs and certain Proteobacteria) under potato may explain the higher β-glucosidase and β-glucosaminidase activities (involved in C and N cycling, respectively) compared to lettuce, onion and broccoli.Previously, Manter et al. [69] showed highly diverse and cultivar specific bacterial endophyte community in potato roots using pyrosequencing; their study included the potato variety studied here.The study by Manter et al. [69] reported the presence of Rheinheimera, Dydobacter, Devosia, Pedobacter, and Pseudoxanthomonas as endophyte communities of potato and emphasized these were not previously identified as endophytes for this crop.However, we did not find any of these bacteria in the soil evaluated under potato, which could be due to several factors including the soil, climate, and environmental interactions.
High abundances of bacteria in the Streptomycetaceae family were associated with onion rhizosphere; which also demonstrated the lowest soil microbial biomass and enzyme activities compared to the other rhizospheres.This seems to be an interesting finding because bacteria in the Streptomycetaceae family are known to produce secondary metabolites which can present antagonistic activities against phytopathogens and could provide plant protection [70] and influence the microbial populations in the competitive rhizosphere surroundings [71,72].For example, studies have reported evidence on Streptomycetaceae members' potential biocontrol towards infection by pathogenic fungi populations of Fusarium oxysporum, the causal agent of onion wilt [73].In addition, onions produce polyphenols, a secondary metabolite, in the daily food intake [74,75] and a rich source of flavonoids [76], which have shown anti-fungal activity in vitro [77].This could explain the low presence of fungal FAME indicators in the soil rhizosphere under onion.Our hypothesis of onion's negative influence on the surrounding microbial populations may also be supported by the increases in bacterial diversity indexes with depth.Bacterial diversity indexes under onion were within the ranges of the other rhizospheres at 0-10 cm, but the indexes became higher under onion compared to the other rhizospheres at 10-20 cm; possibly due to less influence of the onion roots at the lower soil depth.
In this organic farming system under similar management, variations in soil moisture under the different rhizospheres could also explain the differences observed in bacterial distribution.Soil moisture was correlated with higher levels of soil total N, total C, total S, MBC, and enzyme activity supporting previous research findings [78,79].Previous studies have also reported the influence of soil moisture on enzyme activities indirectly through increasing microbial growth and substrate availability [80].This pattern was followed by total C, MBN, and acid phosphatase, similar to findings obtained by Burke et al. [81] in different systems.Singh et al. [82] reported that soil moisture stimulates microorganism activities that enabled organic matter decomposition.In addition, the differences in the root systems and exudates must play another dominant effect in the microbial communities, exemplified by the lowest microbial communities and enzyme activities under onion compared to the other rhizospheres under same management.

Conclusions
Our study demonstrated distinct differences among the rhizospheres evaluated here by combining a high-throughput technique as pyrosequencing with FAME profiling.Distinct differences were also found on the metabolic functional capacity of the microbial community under the different rhizospheres according to enzyme activities of nutrient cycling.Soil under potato and pasture showed the highest fungal:bacterial ratios compared to soil under lettuce, onion, and broccoli.In contrast to pasture and potato, onion showed the lowest microbial biomass, sum of fungal and bacterial FAME indicators, and enzyme activities.Onion rhizosphere also showed a high abundance of the Streptomycetaceae family, which are known to produce secondary metabolites that can exert antagonistic activities against phytopathogens.The possibility that this bacterial family may have reduced important members of this soil microbial community, as reflected in the enzyme activities and microbial biomass, should be explored.Our study approach may serve as the basis for future studies trying to identify microbial communities important in soil processes and health, plant productivity and the nutritional value of crops.

Figure 1 .
Figure 1.Principal component analysis (PCA) of soil microbial community composition according to 10 fatty acid methyl ester (FAME) indicators for fungal and bacterial populations under different rhizospheres in an organic production system.

Figure 3 .
Figure 3. Double hierarchal dendogram to evaluate bacteria distribution in this soil as affected by the different rhizospheres in an organic production system.

Table 1 .
Selected soil properties studied under different rhizospheres in an organic production system.

Table 2 .
Soil microbial community composition according to FAME indicators under different rhizospheres in an organic production system (0-10 cm).

Table 4 .
Soil bacterial diversity indexes studied under different rhizospheres in an organic production system.