Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2017.01007 Microbiology Original Research Variation in Soil Microbial Community Structure Associated with Different Legume Species Is Greater than that Associated with Different Grass Species Zhou Yang 1 Zhu Honghui 2 * Fu Shenglei 3 Yao Qing 1 * 1College of Horticulture, South China Agricultural University, Guangdong Province Key Laboratory of Microbial Signals and Disease Control, Guangdong Engineering Research Center for Litchi, Guangdong Engineering Research Center for Grass Science Guangzhou, China 2State Key Laboratory of Applied Microbiology Southern China, Guangdong Provincial Key Laboratory of Microbial Culture Collection and Application, Guangdong Institute of Microbiology Guangzhou, China 3College of Environment and Planning, Henan University Kaifeng, China

Edited by: Brigitte Mauch-Mani, University of Neuchâtel, Switzerland

Reviewed by: Nicolas Legay, Institut National des Sciences Appliquées Centre Val de Loire, France; Munusamy Madhaiyan, Temasek Life Sciences Laboratory, Singapore

*Correspondence: Qing Yao yaoqscau@scau.edu.cn

Honghui Zhu zhuhh@gdim.cn

This article was submitted to Plant Microbe Interactions, a section of the journal Frontiers in Microbiology

 31 05 2017 2017 8 1007 04 03 2017 19 05 2017 Copyright © 2017 Zhou, Zhu, Fu and Yao. 2017 Zhou, Zhu, Fu and Yao

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Plants are the essential factors shaping soil microbial community (SMC) structure. When most studies focus on the difference in the SMC structure associated different plant species, the variation in the SMC structure associated with phylogenetically close species is less investigated. Legume (Fabaceae) and grass (Poaceae) are functionally important plant groups; however, their influences on the SMC structure are seldom compared, and the variation in the SMC structure among legume or grass species is largely unknown. In this study, we grew three legume species vs. three grass species in mesocosms, and monitored the soil chemical property, quantified the abundance of bacteria and fungi. The SMC structure was also characterized using PCR-DGGE and Miseq sequencing. Results showed that legume and grass differentially affected soil pH, dissolved organic C, total N content, and available P content, and that legume enriched fungi more greatly than grass. Both DGGE profiling and Miseq-sequencing indicated that the bacterial diversity associated with legume was higher than that associated with grass. When legume increased the abundance of Verrucomicrobia, grass decreased it, and furthermore, linear discriminant analysis identified some group-specific microbial taxa as potential biomarkers of legume or grass. These data suggest that legume and grass differentially select for the SMC. More importantly, clustering analysis based on both DGGE profiling and Miseq-sequencing demonstrated that the variation in the SMC structure associated with three legume species was greater than that associated with three grass species.

 legume species grass species soil microbial community structure variation within taxonomical group root exudates [2013]246 2015TX01N036 U1131001 3146421 Foundation for High-level Talents in Higher Education of Guangdong Guangdong Science and Technology Department10.13039/501100007162 NSFC-Guangdong Joint Project National Natural Science Foundation of China10.13039/501100001809 Introduction

Soil ecological processes are primarily driven by soil microbial community (SMC), and plant species is one of the important factors shaping the SMC (Berg and Smalla, 2009; Diouf et al., 2010; Ladygina and Hedlund, 2010; McLaren and Turkington, 2011). When most studies focus on the difference in the SMC associated with different species, few study addresses to what extent the difference is among different species. Plants within a closely related taxonomic group are likely to share a number of common traits, e.g., root biomass, amount, and availability of the rhizodeposits, root defensive strategies, which are the key factors shaping the SMC. Therefore, Berg and Smalla (2009) pointed out that plant phylogenetic positions can have an influence on the SMC. Recently, Bouffaud et al. (2014) demonstrated that the difference in the SMC correlated significantly with the phylogenetic distance of their plant hosts, opening a novel perspective regarding the interactions between plants and the associated SMC.

Legume (Fabaceae) and grass (Poaceae) are the two most functionally important plant groups on the globe (Croser et al., 2006; Liu et al., 2009). In agroecosystems, legume, and grass normally occur as food crops, vegetables, forages, cover crops, and weeds. Due to their wide distribution in the diverse ecosystems, the influences of legume and grass on soil ecological processes have been intensively investigated, particularly the influences on the biogeochemical cycling of nutrients. Dinesh et al. (2006) reported that long-term (12 years) cover cropping with four leguminous species significantly improved the N and C cycling driven by soil microbes, resulting in higher levels of total organic C, dissolved organic C and N, labile organic N, and etc. Cui et al. (2015) indicated that, as cover crops in a subtropical orchard, legume and grass differentially affected the hydrolysis of soil organic P.

Since the SMC is central to the nutrient biogeochemical cycling and soil fertility, better understanding of the difference in the SMC between legume and grass is particularly necessary. Many environmental factors are involved in shaping the SMC, among which plant species is one of the most intensively studied factors (Diouf et al., 2010; Ladygina and Hedlund, 2010; McLaren and Turkington, 2011). In this context, it is reasonable to speculate that legume and grass can shape respectively specific SMC, which is supported by previous reports. For example, the SMC of legume (Lotus corniculatus) was much different from that of grass (Holcus lanatus) as revealed by PLFA (Ladygina and Hedlund, 2010). In more detail, with library construction and cloning, Chen et al. (2014) showed that legume (Trifolium repens) greatly enriched alpha-proteobacteria abundance, but reduced beta-proteobacteria abundance, indicating the improved soil fertility. For fungal taxa, Benitez et al. (2016) demonstrated that two legume cover crops (Trifolium incarnatum and Vicia villosa) preferably enriched arbuscular mycorrhizal fungi. Similar results were also reported by Turner et al. (2013), who employed metatranscriptomics technique. Collectively, these studies suggest that legume and grass tend to select specific microbial taxa, which is probably beneficial to themselves (Bakker et al., 2013). Notwithstanding, few studies have been performed to date to strictly compare the differential effects of legume and grass on the SMC, especially with more plant species, and general understanding of the effects of plant taxonomic groups on their associated SMC is still poor.

Cover cropping in orchards is world widely practiced for the sustainable production, where either legume or grass cover crops is commonly grown (Cui et al., 2015). This special agroecosystem provides ideal platform for the comparison of legume and grass regarding their effects on the SMC. We hypothesized that legume and grass differentially affect the SMC, and that the species within legume or grass may differ in their effects on the SMC. In this study, we constructed mesocosms using the orchard soils to grow three legume species vs. three grass species. We monitored the soil chemical properties, quantified the abundance of bacteria and fungi, and characterized the SMC using DGGE and Miseq-sequencing. By comparing the difference between legume and grass and the variation within three species of legume or grass, we aimed to explore (a) what differences in the SMC shaped by legume and grass, and (b) to what extent the SMC varies within legume or grass species.

 Materials and methods Plant species and soils

To compare the influences of legume and grass on the SMC, three legume species [Stylosanthes guianensis (Aubl.) Sw., Trifolium pratense L., Medicago sativa L.] and three grass species (Paspalum natatum Flüggé, Festuca arundinacea L., Lolium perenne L.) were used. We selected these plant species, because each of them belongs to a different genera and all of them are commonly used as cover crops in orchards. The seeds were commercially obtained (Yifang Seed Lt. Company, Jiangsu, China).

The soils were collected from a subtropical orchard (E 112°54′19″, N 22°40′20″). The orchard was transitioned from forest land and managed for over 30 years (Wang et al., 2015), and the understory of the forest was dominated by Dicranopteris dichotoma with the coverage as high as almost 100% (Zhao et al., 2012). The soil at 0~20 cm were collected and sieved through a 5 mm pore-sized mesh for late use. The soil chemical properties were analyzed as follows: pH 5.42, organic matter content (SOM) 1.04%, total N content (TN) 0.87 g kg−1, total P content (TP) 0.89 g kg−1, total K content (TK) 15.90 g kg−1, available N content (AN, alkali-hydrolyzable N) 67.4 mg kg−1, available P content (AP) 161.0 mg kg−1, available K content (AK) 72.6 mg kg−1.

 Experimental design and sampling

Plastic pot (32 cm length × 18 cm width × 13 cm height) was used as container to construct a mesocosm containing 4.5 kg soils. Seven treatments responding to six plant species plus a control without plants were set up, with each treatment comprising four replicates. In each pot with plants, surface disinfected (10% H2O2 for 15 min) and pre-germinated seeds were sown, and then 50 seedlings were left over at 10 days after sowing. We chose this planting density much higher than that in orchards, in order to produce high root density in soils and thus to reinforce the influence of plants on the SMC. The pots were placed in glasshouse with natural day/night photoperiod and temperature (22~35°C), and watered to the moisture of 20% every 3 days.

Soils were sampled at 160 days after sowing. At harvest, the shoots were firstly removed, and then the top soil layer of about 1 cm were discarded to avoid the aboveground disturbance. All roots (including visible root fragments) in the soils were taken off with forceps. The leftover soils were homogenized with a 2 mm pore-sized mesh and divided into two aliquots. One aliquot was air dried for chemical analysis, the other aliquot was stored at −80°C for molecular analysis.

 Soil chemical analysis

Soil pH was determined with deionized water (2.5:1, w/v) using a glass electrode (Sartorius PB-10). Total organic C (TOC) and dissolved organic C (DOC) were analyzed by titration after wet oxidation with H2SO4 and K2Cr2O7 (Heanes, 1984). TN, TP, and TK were measured using the Kjeldahl method, the molybdenum blue colorimetric method and the flame photometric method, respectively (Kirk, 1950; Sommers and Nelson, 1972). AN, AP, and AK were quantified by alkali-hydrolyzed reduction diffusing method, colorimetric method and flame photometric method, as described previously (Xiong et al., 2008).

 DNA extraction and PCR-DGGE

Soil total DNA was extracted using PowerSoil® DNA Isolation kit (MoBio Laboratories Inc.) following the manufacturer's protocol (Cui et al., 2014). To obtain sufficient amount of DNA and diminish bias of DNA extraction, 3 sub-samples were used for DNA extraction respectively, then put the obtained DNA together to get a composite DNA sample. PCR-DGGE of bacterial 16S rRNA gene and fungal ITS gene were performed with these DNA solutions.

PCR amplifications were conducted in triplicate for each of the soil samples. Nested PCR was conducted with 27F/1492R and GC-341F/518R as the first and second primer set to amplify the V3 region of bacterial 16S rRNA gene (Cui et al., 2014). Similarly, to amplify the fungal ITS1 gene, ITS1F/ITS4 and GC-ITS1F/ITS2 were used as the first and second primer set according to Anderson et al. (2003). Then the PCR products were separated on a D-Code Universal Mutation Detection System (Bio-Rad Laboratories, Inc. USA) with the following conditions: V3 region of bacterial 16S rRNA gene, run on a 8% polyacrylamide gel for 900 min at a constant voltage of 70 V and at 60°C in a 45~70% horizontal denaturant gradient; ITS1 gene, run on a 8% polyacrylamide gel for 780 min at a constant voltage of 70 V and at 58°C in a 30~60% horizontal denaturant gradient. The gels were stained with SYBR Green I and visualized with Molecular Imager® Gel Doc™ XR+ (Bio-Rad Laboratories). Similarities between banding patterns and clusters were analyzed with Quantity One software.

 Real-time PCR quantification of 16S rRNA and its genes

To estimate the abundance of bacteria and fungi, the copy numbers of bacterial 16S rRNA genes and fungal ITS1 genes were quantified using real-time PCR with primers 341F/518R (Moore et al., 2011) and 5.8S/ITS1F (Fierer et al., 2005), respectively. All the real-time PCR assays were triplicately conducted on a CFX96 Optical Real-Time Detection System (Bio-Rad Laboratories, Inc. USA). The PCR reaction mix contained 10 μl 2 × PCR buffer (iQ™ SYBR Green Supermix, Bio-Rad), 2.5 μl of each primer (2 μmol l−1), 1 μl of template DNA, and sterile deionized water added to 20 μl. The qPCR conditions were as follows: an initial denaturation at 95°C for 5 min; 40 cycles of denaturation at 95°C for 15 s, annealing at 56°C (bacteria) or 55°C (fungi) for 30 s. The standard curves were generated using plasmid DNA containing one V3 of bacterial 16S rRNA gene or fungal ITS1 gene. Standard template dilution series from 1.54 × 108 to 1.54 × 102 copies (bacteria) or 8.73 × 108 to 8.73 × 102 copies (fungi) per assay was used. Melting curve analysis and agarose electrophoresis were performed to confirm the specificity of amplification products. Amplification efficiencies of 90~102% were obtained, with R2 > 0.99.

 Miseq-sequencing of 16S rRNA and its genes

In order to evaluate the SMC structure with higher resolution, Miseq-sequencing of the bacterial 16S rRNA gene V3–V4 region and fungal ITS genes ITS1 region were conducted. For this, only three replicate samples were selected from four replicate samples for each treatment, according to the DGGE profiling. For bacteria and fungi, the 515F/806R (Peiffer et al., 2013) and ITS1F/ITS2 (Mueller et al., 2014) primer sets were used, with each reverse primers fused with a unique 6-mer barcode for each samples. Firstly, PCR reactions were carried out in 30 μl reactions with 15 μl of Phusion® High-Fidelity PCR Master Mix (New England Biolabs); 0.2 μmol l−1 of forward and reverse primers, and about 10 ng template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and elongation at 72°C for 60 s, finally 72°C for 5 min. Then the PCR products were electrophoresed on 2% agarose gel for detection, and samples with bright main strip between 400~450 bp were chosen for further experiments. PCR products were mixed in equal density ratios. Then, the mixture PCR products were purified with GeneJET Gel Extraction Kit (Thermo Scientific). Sequencing libraries were generated using NEB Next® Ultra™ DNA Library Prep Kit for Illumina (NEB, USA) following manufacturer's recommendations and index codes were added. The library quality was assessed on the Qubit 2.0 Fluorometer (Invitrogen) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on an Illumina MiSeq platform and 300 bp paired-end reads were generated. Raw sequence reads were de-multiplexed, quality-filtered, processed and analyzed using QIIME (Caporaso et al., 2012). Sequences with ≥ 97% similarity were assigned to the same OUT by UCLUST clustering (Edgar, 2010). One representative sequence for each OUT were picked and used the RDP Classifier and Greengenes database (for bacterial 16S rRNA genes) and UNITE database (for fungal ITS1 genes) to annotate taxonomic information. Samples were rarified to 43,860 and 12,721 tags for bacteria and fungi, respectively, prior to downstream analyses. Because of low production of reads for SG3, we excluded this sample in the downstream analysis to avoid an underestimate with low reads number when randomly subsampling sequences to ensure a same sequence depth.

 Data analysis and statistics

Analysis of variance (ANOVA) was conducted using IBM SPSS statistics software version 21 (SPSS Inc., Chicago, IL) to test significant difference of soil chemical properties, SMC alpha diversity, and copies of bacterial 16S rRNA genes and fungal ITS genes (n = 4, except for sequencing data n = 3), and t-tests were also conducted for legume and grass (n = 12, expect for sequencing data n = 9). Principle component analysis (PCA) on the basis of soil chemical property was performed using IBM SPSS statistics software to reveal the influence of plants. Bacterial and fungal taxa at phylum level selected by plants were plotted in SigmaPlot 12.5. To further uncover biomarkers, linear discriminant analysis (LDA) effect size (LEfSe) was performed to find significantly abundant bacterial and fungal taxa within three groups (legume, grass and control). The factorial Kruskal-Wallis sum-rank test (α = 0.05) was used to identify taxa with significant differential abundances between categories (using the more strict all-against-all comparisons), and then LDA was performed to estimate the effect size of each differentially abundant taxon. Then significant taxa and others were used to generate taxonomic cladograms reflecting differences within groups (Segata et al., 2011). Cluster analysis of DGGE profiles and sequencing data at phylum level were analyzed and visualized using IBM SPSS statistics software and Past 3.11 to evaluate variations of soil samples with UPGMA algorithm.

 Sequence accession numbers

The 16S rRNA gene and ITS DNA sequence data were submitted to NCBI Sequence Read Archive (SRA) with accession number SRX2589894 and SRX2590037 for bacteria and fungi, respectively.

 Results Soil chemical properties as affected by plant species

Soil chemical analysis indicates that legume and grass affected the soil chemical properties in much different patterns. Although they did not differentially affect TOC, TP, TK, AN, AK, and C/N, their differential effects on pH, DOC, TN, AP was clearly observed (Table 1). Soil pH associated with legume was significantly lower than that associated with grass, with both higher than control. Similarly, DOC associated with legume was lower than that associated with grass. In contrast, TN and AP associated with legume were higher than those associated with grass.

Soil chemical property as affected by legume or grass.

Sample pH TOC (g kg−1) DOC (mg kg−1) TN (g kg−1) TP (g kg−1) TK (g kg−1) AN (mg kg−1) AP (mg kg−1) AK (mg kg−1) C/N
Legume SG 5.44 ab 10.3 a 174.1 ab 0.73 d 0.70 b 13.03 b 69.3 a 163.5 ab 15.6 a 14.1 a
TP 5.59 bc 9.3 a 158.6 ab 0.61 bc 0.57 a 10.92 a 64.0 a 171.5 b 19.3 ab 15.2 a
MS 5.72 c 9.3 a 171.0 ab 0.59 b 0.58 a 11.48 ab 71.9 a 166.5 ab 33.0 c 15.9 a
average 5.58** 9.6 NS 167.9* 0.65* 0.62 NS 11.81 NS 68.4 NS 167.2*** 22.6 NS 15.0 NS
Grass PN 5.75 c 9.6 a 199.0 b 0.65 c 0.70 b 12.90 ab 71.2 a 158.5 a 13.3 a 14.9 a
FA 5.75 c 8.7 a 174.1 ab 0.58 b 0.56 a 11.30 ab 69.1 a 156.8 a 16.5 ab 15.0 a
LP 5.78 c 9.9 a 186.6 ab 0.51 a 0.56 a 11.35 ab 66.1 a 158.3 a 23.1 b 19.3 b
average 5.76** 9.4 NS 186.6* 0.58* 0.61 NS 11.85 NS 68.8 NS 157.8*** 17.6 NS 16.4 NS
Control 5.24 a 8.8 a 146.1 a 0.63 bc 0.62 ab 12.20 ab 59.6 a 166.5 ab 68.0 d 14.0 a
F-value of ANOVA 12.581 2.580 2.995 25.542 5.925 3.572 1.106 5.043 176.027 6.586
P-value 0.000 0.049 0.028 0.000 0.001 0.013 0.392 0.002 0.000 0.001

Data followed with lower letters in each column are separated by one-way ANOVA (Turkey HSD multiple range test, P = 0.05). PN, Paspalum natatum; FA, Festuca arundinacea; LP, Lolium perenne; SG, Stylosanthes guianensis; TP, Trifolium pratense; MS, Medicago sativa; Control, no plants. NS and

indicate no significant difference and significant difference at P = 0.05, 0.01, 0.001, respectively, between the averages of legume and grass.

On the basis of soil chemical property, we performed PCA to evaluate the similarity of effect by three legume species or three grass species against control. Figure 1 shows that the soil chemical property associated with legume or grass was different from that without plants, while the soil chemical property associated with six different plant species did not separate from each other well.

Principle component analysis on the basis of soil chemical property. PN, Paspalum natatum; FA, Festuca arundinacea; LP, Lolium perenne; SG, Stylosanthes guianensis; TP, Trifolium pratense; MS, Medicago sativa; CK, Control without plants.

 Soil microbial community structure as affected by plant species

According to DGGE profiles, there was no significant difference in the fungal community structure between legume and grass (Table 2); however, Shannon-Weaver diversity index (H'), species evenness (J), and species richness (S) of the bacterial community associated with legume were significantly higher than those associated with grass. As a whole, the parameters of community structure associated with legume or grass were comparable to those of the control, although some individual species showed significant influences (Table 2).

Alpha diversity of the soil microbial community associated with legume and grass on the basis of DGGE profiling.

Sample H' J S
BACTERIAL COMMUNITY
Legume SG 2.99 c 0.91 c 27.25 ab
TP 2.95 c 0.88 bc 28.29 ab
MS 2.86 bc 0.85 abc 29.33 b
average 2.93*** 0.88*** 28.29***
Grass PN 2.78 abc 0.85 abc 26.53 ab
FA 2.59 a 0.80 a 25.80 a
LP 2.68 ab 0.82 ab 25.80 a
average 2.68*** 0.82*** 26.04***
Control 2.88 bc 0.87 bc 27.25 ab
F-value of ANOVA 7.496 6.559 4.167
P-value 0.000 0.001 0.007
FUNGAL COMMUNITY
Legume SG 2.14 a 0.73 a 18.75 a
TP 2.09 a 0.72 a 17.99 a
MS 2.37 b 0.78 a 21.07 a
average 2.20 NS 0.74 NS 19.27 NS
Grass PN 2.17 ab 0.75 a 18.01 a
FA 2.23 ab 0.77 a 18.01 a
LP 2.08 a 0.72 a 18.25 a
average 2.16 NS 0.75 NS 18.09 NS
Control 2.18 ab 0.75 a 18.50 a
F-value of ANOVA 5.234 2.914 2.303
P-value 0.002 0.032 0.073

Data followed with lower letters in each column are separated by one-way ANOVA (Turkey HSD multiple range test, P = 0.05). H', Shannon-Weaver diversity index; J, species evenness; S, species richness. PN, Paspalum natatum; FA, Festuca arundinacea; LP, Lolium perenne; SG, Stylosanthes guianensis; TP, Trifolium pratense; MS, Medicago sativa; Control, no plants. NS and

indicate no significant difference and significant difference at P = 0.001, respectively, between the averages of legume and grass.

According to the Miseq-sequencing, the alpha diversity (Shannon index) of bacterial community associated with legume was significantly higher than that associated with grass; however, the alpha diversity (OTU, Chao1, and Shannon index) of fungal community associated with legume was significantly lower than that associated with grass (Table 3).

Alpha diversity of the soil microbial community associated with legume and grass on the basis of Miseq-sequencing.

Sample OTU Chao1 Shannon index
BACTERIAL COMMUNITY
Legume SG 1763 ab 3084.4 a 9.19 abc
TP 1847 b 2670.6 a 9.31 c
MS 1812 ab 2600.3 a 9.25 bc
average 1807 NS 2785.1 NS 9.25***
Grass PN 1685 a 2481.8 a 9.01 a
FA 1772 ab 2569.9 a 9.14 abc
LP 1810 ab 2543.4 a 9.13 ab
average 1756 NS 2531.7 NS 9.09***
Control 1781 ab 2606.8 a 9.08 ab
F-value of ANOVA 3.056 0.674 7.647
P-value 0.040 0.673 0.001
FUNGAL COMMUNITY
Legume SG 351 b 453.0 ab 2.77 a
TP 189 a 260.2 a 2.01 a
MS 439 b 486.9 b 4.24 b
average 326* 400.0* 3.01***
Grass PN 468 b 557.6 b 4.33 b
FA 357 b 446.5 ab 4.20 b
LP 463 b 563.5 b 4.62 b
average 429* 522.5* 4.38***
Control 378 b 489.8 b 4.09 b
F-value of ANOVA 15.668 5.151 26.176
P-value 0.000 0.005 0.000

Data followed with lower letters in each column are separated by one-way ANOVA (Turkey HSD multiple range test, P = 0.05). PN, Paspalum natatum; FA, Festuca arundinacea; LP, Lolium perenne; SG, Stylosanthes guianensis; TP, Trifolium pratense; MS, Medicago sativa; Control, no plants. NS and

indicate no significant difference and significant difference at P = 0.5, 0.01, 0.001, respectively, between the averages of legume and grass.

 The abundance of bacteria and fungi as affected by plant species

To investigate the effects of legume and grass on microbial abundance, we performed qRT-PCR targeting the 16S rRNA V4 gene and ITS gene. Both legume and grass significantly increased the bacterial and fungal abundance in soils (Table 4). When no significant difference in bacterial abundance was observed between legume and grass, the fungal abundance associated with legume was significantly higher than that associated with grass (Table 4).

The abundance of bacterial 16S rRNA V4 region and fungal ITS genes in soils as affected by legume and grass.

Sample Bacteria (× 109 copies/g soil) Fungi (× 108 copies/g soil) Bacteria/fungi
Legume SG 7.66 d 7.96 d 10.4 ab
TP 4.40 b 7.50 cd 5.7 a
MS 6.64 bcd 3.44 ab 20.0 b
average 6.23 NS 6.30** 12.1 NS
Grass PN 6.96 cd 4.74 bc 14.9 ab
FA 5.46 bcd 2.81 ab 20.3 b
LP 5.26 bc 3.41 ab 16.9 b
average 5.89 NS 3.65** 17.4 NS
Control 1.26 a 1.22 a 10.7 ab
F-value of ANOVA 16.887 13.723 5.948
P-value 0.000 0.000 0.001

Data followed with lower letters in each column are separated by one-way ANOVA (Turkey HSD multiple range test's, P = 0.05). PN, Paspalum natatum; FA, Festuca arundinacea; LP, Lolium perenne; SG, Stylosanthes guianensis; TP, Trifolium pratense; MS, Medicago sativa; Control, no plants. NS and

indicate no significant difference and significant difference at P = 0.01, respectively, between the averages of legume and grass.

On the basis of Miseq-sequencing, we compared the microbial abundance between control and each plant species. Figure 2A indicates that, at phylum level, all six plant species increased the abundance of Bacteroidetes and Nitrospirae. All three grass species decreased the abundance of Verrucomicrobia; however, two (SG and TP) out of three legume species increased it (Figure 2A), indicating the difference between legume and grass. For fungi, it is clear that all six plant species increased the abundance of Asomycota, but decreased the abundance of Zygomycota (Figure 2B).

The comparison of soil microbial taxonomic profile (relative abundance, %) at phylum level between control (X-axis) and each plant species (Y-axis) on the basis of the Miseq-sequencing. (A) Bacteria; (B) Fungi.

We performed linear discriminant analysis (LDA) to reveal the microbial biomarkers of legume and grass. For bacteria, Figure 3A shows that legume enriched Actinobacteria and Nitrospirae phyla, in contrast to Bacteroidetes phylum enriched by grass. In more details, legume enriched Actinobacteria, Gammaproteobacteria, Betaproteobacteria, and Nitrospira classes, and grass enriched Sphingobacteria class; while Spartobacteria and Alphaproteobacteria classes were enriched in soils without plants (Figure 3A). For fungi, Figure 3B shows that legume, grass and control enriched Ascomycota, Basidiomycota and Zygomycota phyla, respectively. In more detail, legume enriched Sordariomycetes class, and grass enriched Dothideomycetes, Glomeromycetes, Agaricomycetes; while Eurotiomycetes was enriched in soils without plants (Figure 3B). In general, these microbial taxa can be the candidate biomarkers of legume or grass.

Cladogram indicating the phylogenetic distribution of microbial lineages as affected by legume and grass. Each circle's diameter is proportional to the given taxon's relative abundance; circles represent phylogenetic levels from phylum to genus inside out. Red indicates control without plants, green indicates grass, blue indicates legume, and yellow indicates non-significant. (A) Bacteria; (B) Fungi.

 The variation in microbial community structure among plant species

DGGE profiles were subjected to clustering analysis to reveal the similarity of the SMC structure associated with three species of legume or grass. It is obvious that both the bacterial and the fungal community structure associated with three legume species separated from each other, and in contrast, neither the bacterial nor the fungal community structure associated with three grass species separated (Figure 4). These data suggest that the SMC associated with three legume species were much different from each other, while the SMC associated with three grass species were similar to a certain degree. Clustering analysis was also performed on the basis of Miseq-sequencing (at phylum level) to validate the result from DGGE profiling, and the similar result was obtained (Figure 5).

Cluster analysis of the soil microbial community associated with legumes and grasses on the basis of DGGE profiling. (A) Bacteria of legumes; (B) Fungi of legumes; (C) Bacteria of grasses; (D) Fungi of grasses. PN, Paspalum natatum; FA, Festuca arundinacea; LP, Lolium perenne; SG, Stylosanthes guianensis; TP, Trifolium pratense; MS, Medicago sativa; CK, control without plants.

Cluster analysis of the soil microbial community associated with legumes and grasses on the basis of Miseq-sequencing. (A) Bacteria of legumes; (B) Fungi of legumes; (C) Bacteria of grasses; (D) Fungi of grasses. PN, Paspalum natatum; FA, Festuca arundinacea; LP, Lolium perenne; SG, Stylosanthes guianensis; TP, Trifolium pratense; MS, Medicago sativa; CK, control without plants.

 Discussion

Legume and grass, as important functional groups in most plant ecosystems, are much different from each other. Despite those distinguishable differences in physiological and biochemical features (Isobe et al., 2001; Haling et al., 2016; Kidd et al., 2016; Schaller et al., 2016), their differential impacts on edaphic property were explored in this study. To represent the legume group and the grass group well, three species (each from different genera) were selected for each group. We demonstrated that legume and grass not only differentially affected the soil chemical property, but also differentially shaped the SMC structure. Moreover, we found that the variation in the SMC structure associated different legume species was greater than that associated with different grass species.

Firstly, in this study, both legume and grass increased soil pH, and the effect of grass was greater than that of legume. The beneficial effect of plants in alleviating acidity of acidic soils has been reported previously (Cui et al., 2015). Legume is less effective mainly due to the release of proton during the biological N-fixation by symbionts (Haynes, 1983). The differential uptake of different nutrient forms can also modify the soil pH (Thion et al., 2016). In addition to the lower soil pH, the DOC content was also lower in the legume soils than in the grass soils. DOC is easily decomposed by microbes and thus the measured DOC is closely related to microbial activity. According to a survey of 29 sites throughout UK, van den Berg et al. (2012) indicated that DOC was positively correlated with C/N, because decomposition is limited at high at C/N unfavorable to microbial activity. In contrast, legume increased TN and AP in this study. It is readily acknowledged that TN was increased primarily because of the biological N2-fixation (BNF) associated with legume (Li et al., 2015; Bichel et al., 2016; Wu et al., 2017); however, TN is not always increased by legume because the measured TN is the outcome of the balance between N input (BNF) and root N uptake. The unchanged and even decreased TN by legume was also reported previously (Sainju et al., 2017). AP is another parameter increased by legume compared with grass in this study. In this study, we used acidic soils where FePO4 and Al PO4 are the predominant inorganic P. Therefore, the release of P in the rhizosphere of legume due to the chelation of Fe and Al by organic acids can be expected because the secretion of organic acids from legume roots is normally in a larger amount than that from grass roots (Hocking, 2001; Weisskopf et al., 2008; Maltais-Landry, 2015). Moreover, the increased mineralization of organic P by legume is also possible, because phosphatase activity is higher in the rhizosphere of legume than in the rhizosphere of grass regardless of organic P addition (Li et al., 2004). Plant nutritional strategy is an important factor involved in shaping the SMC structure (Thion et al., 2016; Guyonnet et al., 2017), however, the plant biomass and nutritional traits of different species were not monitored in this study. The nutritional strategy of these species and its influence on the SMC deserves further investigation in the future.

With Miseq-sequencing and DGGE profiling, we demonstrated that legume and grass differentially shape the SMC, including community structure and composition. The differential SMC as affected by legume and grass has been reported previously. Turner et al. (2013) indicated that both the bacterial and the fungal community structure in rhizosphere were much different between pea and wheat, as revealed by comparative metatranscriptomics. With PLFA technique, Ladygina and Hedlund (2010) showed the similar result, namely an obvious difference between Lotus corniculatus (legume) and Holcus lanatus (grass). In this study, DGGE profiling indicated that the bacterial community structure (Shannon-Weaver diversity index, species evenness, species richness) associated with legume was much different from that associated with grass, and Miseq-sequencing showed that the fungal community diversity (OTU, Chao1, Shannon index) was also much different between them. The driving force underlying the differential SMC between legume and grass can be, at least, two aspects: different root exudates and different influences on soil chemical property. The different root exudates between legume and grass have been demonstrated previously. Isobe and Tsuboki (1998) and Isobe et al. (2001) compared five legume species and four grass species, and found that legume secreted more amino acids, sugar, and flavonoids than grass. These organic compounds of low molecular weight can play critical regulation role in the SMC (de Nobili et al., 2001; Michalet et al., 2013; Bakker et al., 2015; Szoboszlay et al., 2016). Moreover, Weisskopf et al. (2008) pointed out that the secretion of organic acids by lupin was also greater than that by wheat, and this difference explained a significant proportion (15%) of the difference in their associated SMC. On the other hand, in addition to the amount of exudates, the exudate profile (or composition) of legume and grass is probably much different from each other, and consequently, can be critical in regulating the SMC. However, how the root exudate profile or how a specific organic compound regulates the SMC is far from understanding. Therefore, it is necessary to put these work forward in the future.

When considering the specific microbial taxa in this study, we found that legume enriched more fungi, which was reported previously (Turner et al., 2013). At the phylum level, it seems that grass decreased but legume increased the abundance of Verrucomicrobia. The phylum Verrucomicrobia is recently defined by phylogenetic analysis of 16S rRNA gene sequences, with very few isolates cultured so far (Hedlund, 2010). Its members are a globally distributed, abundant, and active group of soil and water bacteria (Sangwan et al., 2005; Bergmann et al., 2011). Although its function profile is largely unknown due to its unculturability, its possible roles in degrading contaminants and as indicator of soil fertility have been reported (Cardman et al., 2014; Navarrete et al., 2015). At class level, our results indicate that legume enriched Actinobacteria, Gammaproteobacteria, Betaproteobacteria, Nitrospira, and Sordariomycetes, in contrast to Sphingobacteria, Dothideomycetes, Glomeromycetes, and Agaricomycetes enriched by grass. Actinobacteria play diverse roles in their associations with plants (Barka et al., 2016); therefore, the increased abundance of Actinobacteria by legume probably contributes to more comprehensive influence of legume on soil environments over grass. Nitrospira contains one family (Nitrospiraceae), and the largest genera Nitrospira can perform complete nitrification (Daims et al., 2015).

Interestingly, our data show that the variation in the SMC structure associated with different legume species is greater than that associated with different grass species, which is supported by the conclusion of Berg and Smalla (2009), namely, monocotyledonous grass species show highly similar rhizosphere community in diverse studies. To our knowledge, this is the first report to directly reveal the difference between legume and grass in this aspect. As aforementioned, plant can shape the SMC mainly via two mechanisms, e.g., root exudates and soil chemical property. Many work indicated that soil chemical property is the critical factor driving the SMC (Muhammad et al., 2014). For instance, soil pH was demonstrated to be the predictor of soil bacterial community structure at the continental scale (Lauber et al., 2009). Both water quantity and quality of the irrigation regulated the SMC in grapefruit orchards (Bastida et al., 2017). In drained Mediterranean peaty soils, clay, bulk density, and exchangeable calcium drove ammonia-oxidizing community, while silt, ammonium, and exchangeable potassium drove nirK-type denitrifier community (Ciccolini et al., 2016). In this study, three legume species also differentially affected some soil chemical parameters, such as TN, TP, TK, AK, however, PCA indicated that the soil chemical property as a whole did not differ much among them. It is likely that other factors in addition to the variation in the soil chemical property also contributed to the significant segregation of the SMC associated with three legume species. It is the case for three grass species. Consequently, root exudates are probably the primary cause responsible for the difference between legume and grass, that is, the root exudate profile differs much among three legume species, while its does not among three grass species. The root exudate profiles were not characterized in this study; however, it is well acknowledged that flavonoids in the root exudates of legume differ much among different species (Gomaa et al., 2015). It is also the case for the amino acids and organic acids in root exudates. For example, citrate was detectable only in the root exudates of white lupin, but not of field pea and faba bean (Nuruzzaman et al., 2006), and the change patterns of carboxylate in rhizosphere were different among even three cultivars of chickpea (Wouterlood et al., 2004). It is acknowledged that the distribution patterns of secondary metabolites (including non-protein amino acids, organic acids, phenolics) in legume species do not agree with the phylogeny of the species producing them, probably due to convergent evolution, a contribution of endophytic fungi, or the horizontal gene transfer from bacteria (see review by Wink, 2013, and literature therein). Szoboszlay et al. (2016) indicated that different flavonoids (7,4′-dihydroxyflavone and naringenin) differentially regulated the microbial community. How the amino acid profile and organic acid profile affect the microbial community is largely unknown, and more work in the future is expected. Moreover, legume roots normally release much more organic compounds than grass roots, and thus, it is unquestionable that a large amount of root exudates with species-specific profile shapes the unique SMC for each legume species, but not for grass species.

 Author contributions

Designed experiments: YZ, HZ, and QY. Performed experiments: YZ. Analyzed results: YZ, SF, and QY. Wrote the manuscript: YZ, HZ, and QY.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

This study was financially supported by the Foundation for High-level Talents in Higher Education of Guangdong ([2013]246), NSFC Project (3146421), Guangdong Science and Technology Innovation Leadership Foundation (2015TX01N036), NSFC-Guangdong Joint Project (U1131001).

 References Anderson I. C. Campbell C. D. Prosser J. I. (2003). Diversity of fungi in organic soils under a moorland-Scots pine (Pinus sylvestris L.) gradient. Environ. Microbiol. 5, 1121–1132. 10.1046/j.1462-2920.2003.00522.x14641592 Bakker M. G. Bradeen J. M. Kinkel L. L. (2013). Effects of plant host species and plant community richness on streptomycete community structure. FEMS Microbiol. Ecol. 83, 596–606. 10.1111/1574-6941.1201723013423 Bakker M. G. Chaparro J. M. Manter D. K. Vivanco J. M. (2015). Impacts of bulk soil microbial community structure on rhizosphere microbiomes of Zea mays. Plant Soil 392, 115–126. 10.1007/s11104-015-2446-0 Barka E. A. Vatsa P. Sanchez L. Gaveau-Vaillant N. Jacquard C. Klenk H. P. . (2016). Taxonomy, physiology, and natural products of actinobacteria. Microbiol. Mol. Biol. Rev. 80, 1–43. 10.1128/MMBR.00019-1526609051 Bastida F. Torres I. F. Romero-Trigueros C. Baldrian P. Větrovský T. Bayona J. M. . (2017). Combined effects of reduced irrigation and water quality on the soil microbial community of a citrus orchard under semi-arid conditions. Soil Biol. Biochem. 104, 226–237. 10.1016/j.soilbio.2016.10.024 Benitez M. S. Taheri W. I. Lehman R. M. (2016). Selection of fungi by candidate cover crops. Appl. Soil Ecol. 103, 72–82. 10.1016/j.apsoil.2016.03.016 Berg G. Smalla K. (2009). Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 68, 1–13. 10.1111/j.1574-6941.2009.00654.x19243436 Bergmann G. T. Bates S. T. Eilers K. G. Lauber C. L. Caporaso J. G. Walters W. A. . (2011). The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol. Biochem. 43, 1450–1455. 10.1016/j.soilbio.2011.03.01222267877 Bichel A. Oelbermann M. Echarte L. (2016). Impact of residue addition on soil nitrogen dynamics in intercrop and sole crop agroecosystems. Geoderma. 10.1016/j.geoderma.2016.11.023. [Epub ahead of print]. Bouffaud M. L. Poirier M. A. Muller D. Moënne-Loccoz Y. (2014). Root microbiome relates to plant host evolution in maize and other Poaceae. Environ. Microbiol. 16, 2804–2814. 10.1111/1462-2920.1244224588973 Caporaso J. G. Lauber C. L. Walters W. A. Berg-Lyons D. Huntley J. Fierer N. . (2012). Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. 10.1038/ismej.2012.822402401 Cardman Z. Arnosti C. Durbin A. Ziervogel K. Cox C. Steen A. . (2014). Verrucomicrobia are candidates for polysaccharide-degrading bacterioplankton in an Arctic fjord of Svalbard. Appl. Environ. Microbiol. 80, 3749–3756. 10.1128/AEM.00899-1424727271 Chen Y. Wen X. Sun Y. Zhang J. Wu W. Liao Y. (2014). Mulching practices altered soil bacterial community structure and improved orchard productivity and apple quality after five growing seasons. Sci. Hortic. 172, 248–257. 10.1016/j.scienta.2014.04.010 Ciccolini V. Bonari E. Ercoli L. Pellegrino E. (2016). Phylogenetic and multivariate analyses to determine the effect of agricultural land-use intensification and soil physico-chemical properties on N-cycling microbial communities in drained Mediterranean peaty soils. Biol. Fertil. Soils 52, 811–824. 10.1007/s00374-016-1121-9 Croser J. S. Lülsdorf M. M. Davies P. A. Clarke H. J. Bayliss K. L. Mallikarjuna N. . (2006). Toward doubled haploid production in the Fabaceae: progress, constraints, and opportunities. Crit. Rev. Plant Sci. 25, 139–157. 10.1080/07352680600563850 Cui H. Wang C. Gu Z. Zhu H. Fu S. Yao Q. (2014). Evaluation of soil storage methods for soil microbial community using genetic and metabolic fingerprintings. Eur. J. Soil Biol. 63, 55–63. 10.1016/j.ejsobi.2014.05.006 Cui H. Zhou Y. Gu Z. Zhu H. Fu S. Yao Q. (2015). The combined effects of cover crops and symbiotic microbes on phosphatase gene and organic phosphorus hydrolysis in subtropical orchard soils. Soil Biol. Biochem. 82, 119–126. 10.1016/j.soilbio.2015.01.003 Daims H. Lebedeva E. V. Pjevac P. Han P. Herbold C. Albertsen M. . (2015). Complete nitrification by Nitrospira bacteria. Nature 528, 504–509. 10.1038/nature1646126610024 de Nobili M. Contin M. Mondini C. Brookes P. C. (2001). Soil microbial biomass is triggered into activity by trace amounts of substrate. Soil Biol. Biochem. 33, 1163–1170. 10.1016/S0038-0717(01)00020-7 Dinesh R. Suryanarayana M. A. Ghoshal Chaudhuri S. Sheeja T. E. Shiva K. N. (2006). Long-term effects of leguminous cover crops on biochemical and biological properties in the organic and mineral layers of soils of a coconut plantation. Eur. J. Soil Biol. 42, 147–157. 10.1016/j.ejsobi.2005.12.004 Diouf M. Baudoin E. Dieng L. Assigbetsé K. Brauman A. (2010). Legume and gramineous crop residues stimulate distinct soil bacterial populations during early decomposition stages. Can. J. Soil Sci. 90, 289–293. 10.4141/CJSS09034 Edgar R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. 10.1093/bioinformatics/btq46120709691 Fierer N. Jackson J. A. Vilgalys R. Jackson R. B. (2005). Assessment of soil microbial community structure by use of taxon-specific quantitative PCR assays. Appl. Environ. Microbiol. 71, 4117–4120. 10.1128/AEM.71.7.4117-4120.200516000830 Gomaa N. H. Hassan M. O. Fahmy G. M. González L. Hammouda O. Atteya A. M. (2015). Flavonoid profiling and nodulation of some legumes in response to the allelopathic stress of Sonchus oleraceus L. Acta Bot. Bras. 29, 553–560. 10.1590/0102-33062015abb0153 Guyonnet J. P. Vautrin F. Meiffren G. Labois C. Cantarel A. A. Michalet S. . (2017). The effects of plant nutritional strategy on soil microbial denitrification activity through rhizosphere primary metabolites. FEMS Microbiol. Ecol. 93:fix022. 10.1093/femsec/fix02228334144 Haling R. E. Yang Z. Shadwell N. Culvenor R. A. Stefanski A. Ryan M. H. . (2016). Growth and root dry matter allocation by pasture legumes and a grass with contrasting external critical phosphorus requirements. Plant Soil 407, 67–79. 10.1007/s11104-016-2808-2 Haynes R. J. (1983). Soil acidification induced by leguminous crops. Grass Forage Sci. 38, 1–11. 10.1111/j.1365-2494.1983.tb01614.x Heanes D. L. (1984). Determination of total organic-C in soils by an improved chromic acid digestion and spectrophotometric procedure. Commun. Soil Sci. Plan. 15, 1191–1213. 10.1080/00103628409367551 Hedlund B. P. (2010). Phylum XXIII. Verrucomicrobia phyl. nov., in Bergey's Manual® of Systematic Bacteriology, eds Krieg N. R. Staley J. T. Brown D. R. Hedlund B. P. Paster B. J. Ward N. L. . (New York: Williams & Wilkins), 795–841. Hocking P. J. (2001). Organic acids exuded from roots in phosphorus uptake and aluminum tolerance of plants in acid soils. Adv. Agron. 74, 63–97. 10.1016/S0065-2113(01)74031-X Isobe K. Tateishi A. Nomura K. Inoue H. Tsuboki Y. (2001). Flavonoids in the extract and exudate of the roots of leguminous crops. Plant Prod. Sci. 4, 278–279. 10.1626/pps.4.278 Isobe K. Tsuboki Y. (1998). Relationship between the amount of root exudate and the infection rate of arbuscular mycorrhizal fungi in gramineous and leguminous crops. Plant Prod. Sci. 1, 37–38. 10.1626/pps.1.37 Kidd D. R. Ryan M. H. Haling R. E. Lambers H. Sandral G. A. Yang Z. . (2016). Rhizosphere carboxylates and morphological root traits in pasture legumes and grasses. Plant Soil 402, 77–89. 10.1007/s11104-015-2770-4 Kirk P. L. (1950). Kjeldahl method for total nitrogen. Anal. Chem. 22, 354–358. 10.1021/ac60038a038 Ladygina N. Hedlund K. (2010). Plant species influence microbial diversity and carbon allocation in the rhizosphere. Soil Biol. Biochem. 42, 162–168. 10.1016/j.soilbio.2009.10.009 Lauber C. L. 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. 10.1128/AEM.00335-0919502440 Li Q. Song Y. Li G. Yu P. Wang P. Zhou D. (2015). Grass-legume mixtures impact soil N, species recruitment, and productivity in temperate steppe grassland. Plant Soil 394, 271–285. 10.1007/s11104-015-2525-2 Li S. Li L. Zhang F. Tang C. (2004). Acid phosphatase role in chickpea/maize intercropping. Ann. Bot. 94, 297–303. 10.1007/s11104-015-2525-215238349 Liu Q. Ge X. Chen W. Columbus J. T. (2009). Grass (Poaceae) richness patterns across China's nature reserves. Plant Ecol. 201, 531–551. 10.1007/s11258-008-9483-1 Maltais-Landry G. (2015). Legumes have a greater effect on rhizosphere properties (pH, organic acids and enzyme activity) but a smaller impact on soil P compared to other cover crops. Plant Soil 394, 139–154. 10.1007/s11104-015-2518-1 McLaren J. R. Turkington R. (2011). Plant identity influences decomposition through more than one mechanism. PLoS ONE 6:e23702. 10.1371/journal.pone.002370221858210 Michalet S. Rohr J. Warshan D. Bardon C. Roggy J. C. Domenach A. M. . (2013). Phytochemical analysis or mature tree root exudates in situ and their role in shaping soil microbial communities in relation to tree N-acquisition strategy. Plant Physiol. Biochem. 72, 169–177. 10.1016/j.plaphy.2013.05.00323727287 Moore T. A. Xing Y. Lazenby B. Lynch M. D. J. Schiff S. Robertson W. D. . (2011). Prevalence of anaerobic ammonium-oxidizing bacteria in contaminated groundwater. Environ. Sci. Technol. 45, 7217–7225. 10.1021/es201243t21786759 Mueller R. C. Paula F. S. Mirza B. S. Rodrigues J. L. M. Nüsslein K. Bohannan B. J. M. (2014). Links between plant and fungal communities across a deforestation chronosequence in the Amazon rainforest. ISME J. 8, 1548–1550. 10.1038/ismej.2013.253[10.1038/ismej.2013.253]24451208 Muhammad N. Dai Z. Xiao K. Meng J. Brookes P. C. Liu X. . (2014). Changes in microbial community structure due to biochars generated from different feedstocks and their relationships with soil chemical properties. Geoderma 226, 270–278. 10.1016/j.geoderma.2014.01.023 Navarrete A. A. Soares T. Rossetto R. van Veen J. A. Tsai S. M. Kuramae E. E. (2015). Verrucomicrobial community structure and abundance as indicators for changes in chemical factors linked to soil fertility. Antonie van Leeuwenhoek 108, 741–752. 10.1007/s10482-015-0530-326184407 Nuruzzaman M. Lambers H. Bolland M. D. Veneklaas E. J. (2006). Distribution of carboxylates and acid phosphatase and depletion of different phosphorus fractions in the rhizosphere of a cereal and three grain legumes. Plant Soil 281, 109–120. 10.1007/s11104-005-3936-2 Peiffer J. A. Spor A. Koren O. Jin Z. Tringe S. G. Dangl J. L. . (2013). Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl Acad. Sci. U.S.A. 110, 6548–6553. 10.1073/pnas.130283711023576752 Sainju U. M. Singh H. P. Singh B. P. (2017). Cover crops for enriching soil carbon and nitrogen under bioenergy sorghum, in Global Soil Security, eds Field D. J. Morgan C. L. S. McBratney A. B. (Springer International Publishing), 181–192. Sangwan P. Kovac S. Davis K. E. Sait M. Janssen P. H. (2005). Detection and cultivation of soil Verrucomicrobia. Appl. Environ. Microbiol. 71, 8402–8410. 10.1128/AEM.71.12.8402-8410.200516332828 Schaller J. Roscher C. Hillebrand H. Weigelt A. Oelmann Y. Wilcke W. . (2016). Plant diversity and functional groups affect Si and Ca pools in aboveground biomass of grassland systems. Oecologia 182, 277–286. 10.1007/s00442-016-3647-927164912 Segata N. Izard J. Waldron L. Gevers D. Miropolsky L. Garrett W. S. . (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. 10.1186/gb-2011-12-6-r6021702898 Sommers L. E. Nelson D. W. (1972). Determination of total phosphorus in soils: a rapid perchloric acid digestion procedure. Soil Sci. Soc. Am. J. 36, 902–904. 10.2136/sssaj1972.03615995003600060020x Szoboszlay M. White-Monsant A. Moe L. A. (2016). The effect of root exudate 7, 4′-dihydroxyflavone and naringenin on soil bacterial community structure. PLoS ONE 11:e0146555. 10.1371/journal.pone.014655526752410 Thion C. E. Poirel J. D. Cornulier T. De Vries F. T. Bardgett R. D. Prosser J. I. (2016). Plant nitrogen-use strategy as a driver of rhizosphere archaeal and bacterial ammonia oxidiser abundance. FEMS Microbiol. Ecol. 92:fiw091. 10.1093/femsec/fiw09127130939 Turner T. R. Ramakrishnan K. Walshaw J. Heavens D. Alston M. Swarbreck D. . (2013). Comparative metatranscriptomics reveals kingdom level changes in the rhizosphere microbiome of plants. ISME J. 7, 2248–2258. 10.1038/ismej.2013.11923864127 van den Berg L. J. L. Shotbolt L. Ashmore M. R. (2012). Dissolved organic carbon (DOC) concentrations in UK soils and the influence of soil, vegetation type and seasonality. Sci. Total Environ. 427, 269–276. 10.1016/j.scitotenv.2012.03.06922542300 Wang C. Gu Z. Cui H. Zhu H. Fu S. Yao Q. (2015). Differences in arbuscular mycorrhizal fungal community composition in soils of three land use types in subtropical hilly area of Southern China. PLoS ONE 10:e0130983. 10.1371/journal.pone.013098326107285 Weisskopf L. Bayon R. C. L. Kohler F. Page V. Jossi M. Gobat J. M. . (2008). Spatio-temporal dynamics of bacterial communities associated with two plant species differing in organic acid secretion: a one-year microcosm study on lupin and wheat. Soil Biol. Biochem. 40, 1772–1780. 10.1016/j.soilbio.2008.02.018 Wink M. (2013). Evolution of secondary metabolites in legumes (Fabaceae). S. Afr. J. Bot. 89, 164–175. 10.1016/j.sajb.2013.06.006 Wouterlood M. Cawthray G. R. Turner S. Lambers H. Veneklaas E. J. (2004). Rhizosphere carboxylate concentrations of chickpea are affected by genotype and soil type. Plant Soil 261, 1–10. 10.1023/B:PLSO.0000035568.28893.f6 Wu G. Liu Y. Tian F. Shi Z. (2017). Legumes functional group promotes soil organic carbon and nitrogen storage by increasing plant diversity. Land Degrad. Dev. 28, 1336–1344. 10.1002/ldr.2570 Xiong Y. Xia H. Li Z. Cai X. Fu S. (2008). Impacts of litter and understory removal on soil properties in a subtropical Acacia mangium plantation in China. Plant Soil 304, 179–188. 10.1007/s11104-007-9536-6 Zhao J. Wan S. Li Z. Shao Y. Xu G. Liu Z. . (2012). Dicranopteris-dominated understory as major driver of intensive forest ecosystem in humid subtropical and tropical region. Soil Biol. Biochem. 49, 78–87. 10.1016/j.soilbio.2012.02.020