ORIGINAL RESEARCH article
Impacts of Fertilization Regimes on Arbuscular Mycorrhizal Fungal (AMF) Community Composition Were Correlated with Organic Matter Composition in Maize Rhizosphere Soil
- Jiangsu Provincial Key Lab for Organic Solid Waste Utilization, Nanjing Agricultural University, Nanjing, China
The understanding of the response of arbuscular mycorrhizal fungi (AMF) community composition to fertilization is of great significance in sustainable agriculture. However, how fertilization influences AMF diversity and composition is not well-established yet. A field experiment located in northeast China in typical black soil (Chernozem) was conducted and high-throughput sequencing approach was used to investigate the effects of different fertilizations on the variation of AMF community in the rhizosphere soil of maize crop. The results showed that AMF diversity in the maize rhizosphere was significantly altered by different fertilization regimes. As revealed by redundancy analysis, the application of organic manure was the most important factor impacting AMF community composition between samples with and without organic manure, followed by N fertilizer and P fertilizer inputs. Moreover, the organic matter composition in the rhizosphere, determined by GC–MS, was significantly altered by the organic manure amendment. Many of the chemical components displayed significant relationships with the AMF community composition according to the Mantel test, among those, 2-ethylnaphthalene explained the highest percentage (54.2%) of the variation. The relative contents of 2-ethylnaphthalene and 2, 6, 10-trimethyltetradecane had a negative correlation with Glomus relative abundance, while the relative content of 3-methylbiphenyl displayed a positive correlation with Rhizophagus. The co-occurrence patterns in treatments with and without organic manure amendment were analyzed, and more hubs were detected in the network of soils with organic manure amendment. Additionally, three operational taxonomic units (OTUs) belonging to Glomerales were identified as hubs in all treatments, indicating these OTUs likely occupied broad ecological niches and were always active for mediating AMF species interaction in the maize rhizosphere. Taken together, impacts of fertilization regimes on AMF community composition were correlated with organic matter composition in maize rhizosphere soil and the application of manure could activate more AMF species to interact with other species in the maize rhizosphere. This knowledge can be valuable in regulating the symbiotic system of plants and AMF, maintaining the health and high yields of crops and providing a primary basis for rational fertilization.
An increase in fertilizer application raises a question about how nutrient-enriched soil influences soil microbial communities. Some studies have shown that the balanced fertilization, including supplementation of nitrogen (N), phosphorus (P), and potassium (K), changed the soil microbial community and enhanced the efficient metabolism of microorganisms because of the increase in nutrient availability (Lin et al., 2012). This can lead to high soil productivity in most situations; however, excessive and inappropriate chemical fertilizer application can cause a series of environmental problems and soil degradation. In contrast to chemical fertilization, the application of organic manure to agricultural ecosystems has been used to increase soil organic matter (SOM) content and accordingly, to improve soil physical, chemical and biological qualities, which subsequently enhanced crop yield (Celik et al., 2004; Walker et al., 2004; Clemente et al., 2005). The application of organic manure can especially stimulate the microbial population and diversity (Wander et al., 1994; Zhong et al., 2010). The microorganisms are important component of soil ecosystems that characterize soil fertility (Lueders et al., 2006); thus, it is important to understand the effects of organic and inorganic fertilizer applications on soil microbial communities. For example, soil fungal diversity was found to be a good and sensitive indicator of soil fertility (He et al., 2008). The arbuscular mycorrhizal fungi (AMF) are a group of important soil organisms mediating multiple functions in agro-ecosystems.
Arbuscular mycorrhizal fungi are ancient root symbionts, and their first appearance coincides with the emergence of land plants (Bonfante and Genre, 2008). These organisms form a root symbiosis with approximately 80% of terrestrial plant species and improve nutrient and water uptake as well as pathogen resistance of their hosts in exchange for plant-assimilated carbon (Smith and Read, 2010). Therefore, it is increasingly accepted that AMF play an important role in the agro-ecosystem function. The demand of agricultural soils for certain nutrients could be reflected by changes in AMF community composition and/or diversity (Lin et al., 2012). However, how the anthropogenic disturbances drive AMF community shifts has not been examined (Rillig and Mummey, 2006; Van Der Heijden et al., 2008).
All AMF belong to the phylum Glomeromycota (Schüßler and Walker, 2010), and most are obligate symbionts and not culturable alone. The identification of AMF taxa in the environment is therefore not dependent on conventional culturing methods. Instead, molecular approaches (e.g., targeted polymerase chain reaction (PCR) followed by high-throughput amplicon sequencing) have become a widely used methodology for studying AMF communities (Gorzelak et al., 2012).
In this study, we used amplicon sequencing (Illumina MiSeq sequencing) to survey AMF community composition, diversity and interactions in response to organic and inorganic fertilizer applications in the rhizosphere soil of the maize crop. The objective of this study were to clarify how different fertilization regimes affect AMF community composition in the rhizosphere soil of maize crop and which was the most important key factor changing AMF community composition in the black soils.
Materials and Methods
Description of the Experiment
The field fertilization experiment was conducted in 2013 at the Modern Agricultural Science and Technology Demonstration Station (126°50′33″ E, 45°50′44″ N), Heilongjiang Academy of Agricultural Sciences, Harbin City, Heilongjiang Province. This region has a typical monsoon climate with an annual average temperature of 3.6°C and 486.4–543.6 mm of precipitation. The cropping regime is dominated by one maize crop per year. There were four experimental treatments: (1) NK, only mineral N and K fertilizers were applied; (2) NPK, only mineral N, P, and K fertilizer were applied; (3) NPKM, mineral N, P, and K plus organic manure was applied; and (4) M, only organic manure was applied. There were three replicates of each treatment as three independent plots with a completely randomized design. The chemical fertilizers were applied at the annual rate of 165 kg N ha-1, 60 kg P2O5 ha-1, and 75 kg K2O ha-1 as urea, superphosphate, and potassium chloride, respectively. In the NPKM treatment, the chemical fertilizers were applied at rates of 132 kg N ha-1, 60 kg P2O5 ha-1, and 75 kg K2O ha-1, respectively, plus 3,800 kg ha-1 of organic manures, in which the total N applied was equal to that in NPK treatment but it was divided with 20% organic N and 80% inorganic N. In the M treatment, only 19,020 kg ha-1 of organic manure was applied with the equal N applied to that in NPK treatment. The organic manure contained 40% moisture, 47.8% organic matter, 1.86% TN, 3.11% P2O5, and 0.85% K2O.
Rhizosphere Soil Sampling and DNA Extraction
The rhizosphere soil samples were collected from each treatment in October, 2014. One composite rhizosphere sample collected from each plot consisted of the roots of five randomly selected maize plants. The roots were shaken vigorously to separate soil that was not tightly adhering to the roots. Then, 5 g of plant roots with firmly adhering soil was resuspended in 30 ml of distilled water and treated in a vortex for 1 min at high speed. The soil in the distilled water was used as the rhizosphere soil for analysis. After vacuum freeze drying, one part was stored at 4°C for GC–MS analysis, and the other part was stored at -80°C for DNA extraction.
Three rhizosphere samples from each fertilizer treatment were used for DNA extraction. The total soil DNA was extracted from 0.25 g of freeze-dried soil using a PowerSoil DNA Isolation Kit (Mo Bio Laboratories, Inc., Carlsbad, CA, USA) according to the manufacturer’s instructions. The extracted DNA was evaluated using a 1% agarose gel.
18S rRNA Gene Amplification and Illumina Sequencing
Polymerase chain reaction amplifications were conducted with the AMV4.5NF (5′-AAGCTCGTAGTTGAATTTCG-3′)/AMDGR (5′-CCCAACTATCCCTATTAATCAT-3′) primer set that amplifies the 18S rRNA gene. The primer set was selected in terms of the reproducibility and the ability to accurately describe AMF communities (Lumini et al., 2010). The reverse primer included a 6-bp error-correcting barcode unique to each sample. Amplicon sequencing was performed using the Illumina MiSeq platforms at Genesky Biotechnologies, Inc. (Shanghai, China).
We used the Quantitative Insights Into Microbial Ecology (QIIME) toolkit (Caporaso et al., 2010) and the UPARSE pipeline (Edgar, 2013) to treat raw high-throughput sequencing data. Barcodes and the standard primer sets were excluded. In brief, sequences below the quality score of 25 and fewer than 200 bp in length were excluded (Lin et al., 2012). Illumina MiSeq sequencing data were pretreated to remove the chimeras from the datasets. After optimizing the sequences, there were 1,214,071 valid sequences. The dominant length distribution was approximately 221–240 bp.
Then, we used the UPARSE pipeline to make an OTU table. The sequences were binned into OTUs using a 97% identity threshold, and the most abundant sequence from each OTU was selected as a representative sequence for that OTU. We used the SILVA database (SSU 119) and kmer searching1 to assign taxonomic data to each representative sequence. However, many fungal sequences that differed from Glomeromycota were obtained. Overall, 165,805 Glomeromycota sequences (13.5% of the total) were obtained for the primer pair AMV4.5NF/AMDGR in this study. The other major fungal groups were Chytridiomycota, for which 338,609 sequences (27.7% of the total) were detected, followed by Basidiomycota, Zygomycota, Ascomycota, and Blastocladiomycota (240,671 sequences, 19.7%; 105,528 sequences, 8.6%; 26,405 sequences, 2.2%; 12,249 sequences, 1.0%, respectively) (Table 1). Based on 97% species similarity, 42 OTUs of Glomeromycota were obtained out of 873 OTUs. Each sample was rarefied to the identical number of Glomeromycota reads (6,622) for downstream analyses. The ACE index, Shannon index and rarefaction curves were calculated by QIIME. All sequences were deposited in the NCBI Sequence Read Archive (SRA) database (accession numbers SRP078384).
TABLE 1. The proportional distributions of the sequences from each fungal phylum detected with the primer set AMV4.5NF/AMDGR for all soil samples.
Dissolved Organic Matter Extraction from the Rhizosphere Soil for GC–MS Analysis
Ethyl acetate-extracted dissolved organic matter (DOM) was employed in this study for GC–MS analysis. Owing to the high diversity of soil DOM, the ethyl acetate-extracted DOM does not contain all the complex soil SOM, but it provides a useful and thorough overview of likely molecular classes in DOM.
We used the approach of Song et al. (2016), which was modified by Pramanik et al. (2000) and Li et al. (2012). In brief, part of each soil sample (0.5 g) and 5 ml ethyl acetate were added to a centrifuge tube. After shaking at 170 rpm at 30°C for 1 h, the suspension was filtered (0.45 μm) into a clean centrifuge tube and was concentrated to 0.5 ml at 35°C using a vacuum rotary evaporator. The concentrated samples were then stored for GC–MS analysis.
The extracts of rhizosphere soil were analyzed by GC–MS. The initial oven temperature of 60°C was held for 3 min, increased at a rate of 5°C min-1 to 150°C, further increased at a rate of 3°C min-1 to 240°C, and then held for 10 min. The mass spectrometer was operated in the electron ionization mode at 70 eV with a source temperature of 230°C. A continuous scan from 40 to 500 m/z was used. Helium was used as the carrier gas at a linear velocity of 1.0 ml min-1. The mass spectra of the extracts were compared with those in the National Institute of Standards and Technology (NIST) database (Version 2.0).
Statistical analyses were conducted by the IMB SPSS statistical software package version 20 (IBM Corporation, New York, NY, USA). Data were analyzed with one-way ANOVA, and Duncan’s least significant difference (P < 0.05) was used to compare the means for each variable. Independent-samples t-test was performed to test significance (P < 0.05) of the means of ACE index and Shannon index between the treatment group with organic amendment (+M, including NPKM and M treatments) and the treatment group without organic amendment (-M, including NK and NPK treatments).
The GC–MS data were expressed as the percentage of each ethyl acetate-extracted DOM in each sample, which was calculated by dividing the sum of all peak areas per sample into its corresponding peak area. In all tests (Mantel test, Monte Carlo permutation test, Pearson correlations, and Permutational multivariate analysis), a P-value < 0.05 was considered statistically significant. Statistical procedures were carried out with R (version 2.15.0). Redundancy analysis (RDA), Monte Carlo permutation test and Permutational multivariate analysis were performed in R with the vegan package. In addition, we employed presence and absence matrix as the environmental factors for carrying out RDA analysis and Monte Carlo permutation test.
Mantel test was performed using the mantel function in R with the vegan package to calculate the correlation between the AMF community composition and ethyl acetate-extracted DOM composition. The Pearson correlation between the relative abundances of AMF genera and the related ethyl acetate-extracted DOM from the GC–MS analysis (P < 0.05 Mantel test) was determined using R.
Network analyses were carried out to better comprehend the interrelationship and interaction within the AMF community. We calculated all possible Pearson’s rank correlation coefficients via mothur (version 1.29.2) for analyzing the networks. We filtered the correlation data with a cut-off at an absolute r-value of 0.6 ~ 0.93, and then used a significant P-value of <0.05 to filter the data from the above step to improve the veracity of the networks. The networks were visualized with the interactive platform Gephi with Fruchterman–Reingold layout (Bastian et al., 2009), and a set of measures (average clustering coefficient, average path length, and modularity) was calculated (Newman, 2006).
Effect of Fertilizations on α-Diversity of AMF Communities
The number of AMF OTUs significantly increased with the application of organic manure (M) compared to the other treatments (Supplementary Figure S1). Shannon index and ACE index were calculated based on OTU matrix of each sample in different fertilization treatments (Figure 1). The ACE index was significantly higher in the M treated soil.
FIGURE 1. The AMF ACE index and Shannon index for different treatments. Bars indicated the SDs. NK, soil treated with chemical N and K fertilizer; NPK, soil treated with chemical N, P, and K fertilizer; NPKM, soil treated with chemical fertilizer (N, P, and K) plus organic manure; and M, soil treated with only organic manure. Different letters on the error bars indicate significant differences at P < 0.05.
Effect of Fertilizations on AMF Community Composition
Variations in AMF community composition were clearly shown in different treatments. The first principal component explained 62.22% of the variation in the AMF community, which was predominantly sorted along the variable present or absent organic fertilization (Figure 2). Therefore, manure fertilization was the main reason to shift the AMF community in M and NPKM treatments compared to those in NK and NPK treatments. According to the Permutational multivariate analysis for the AMF community (Supplementary Table S1), the treatments without manure amendment (-M), including NK and NPK, were significantly (r2 = 0.41559, P = 0.003) different from the treatments with manure amendment (+M), including NPKM and M, in terms of community composition level. Hence, we divided the four treatments into two groups (-M and +M) for subsequent analyses.
FIGURE 2. Redundancy analysis (RDA) of the AMF community based on OTU matrix from different treatments. The different color points represent different treatments. The three arrows represent three environment factors: M (+/-) represents treatments with (+) and without (-) organic manure; P (+/-) represents treatments with (+) and without (-) P fertilizer; N (+/-) represents treatments with (+) and without (-) N fertilizer. According to Monte Carlo permutation test the significances are: M (+/-): r2 = 0.7695, P = 0.002; P (+/-): r2 = 0.5048, P = 0.04; N (+/-): r2 = 0.6472, P = 0.022. NK, soil treated with chemical N and K fertilizer; NPK, soil treated with chemical N, P, and K fertilizer; NPKM, soil treated with chemical fertilizer (N, P, and K) plus organic manure; and M, soil treated with only organic manure. The number (1, 2, and 3) following the treatment indicates replicates of each treatment.
Concomitantly, the application of organic manure was the key factor bringing about changes in AMF community composition in the maize rhizosphere and had the most explained variance (r2 = 0.7695, P = 0.002) with Monte Carlo permutation test (Supplementary Table S2). The second most influencing factor was N fertilizer (r2 = 0.6472, P = 0.022) followed by P fertilizer (r2 = 0.5048, P = 0.04).
Effect of Fertilizations on Ethyl Acetate-Extracted DOM Composition of Rhizosphere
According to the RDA, the first principal component explained 57.88% of the variation in ethyl acetate-extracted DOM composition. Shifts in the ethyl acetate-extracted DOM composition between the +M and -M treatments were detected, suggesting that organic manure caused the maximum changes in ethyl acetate-extracted DOM composition (Figure 3). The permutational multivariate analysis indicated that the chemical components in the -M treatments (NK and NPK) were significantly parted from that in the +M treatments (NPKM and M) (Supplementary Table S3). The difference in ethyl acetate-extracted DOM composition between treatments was very similar to the results of the AMF community composition. Hence, there was a potential relationship between ethyl acetate-extracted DOM composition and AMF community composition.
FIGURE 3. Redundancy analysis of the identified ethyl acetate-extracted DOM. The different fertilizer treatments for individual samples. The three arrows represent three environment factors: M (+/-) represents treatments with (+) and without (-) organic manure; P (+/-) represents treatments with (+) and without (-) P fertilizer; N (+/-) represents treatments with (+) and without (-) N fertilizer. According to Monte Carlo permutation test the significances are: M (+/-): r2 = 0.8608, P = 0.002; P (+/-): r2 = 0.4612, P = 0.73; N (+/-): r2 = 0.3676, P = 0.136. NK, soil treated with chemical N and K fertilizer; NPK, soil treated with chemical N, P, and K fertilizer; NPKM, soil treated with chemical fertilizer (N, P, and K) plus organic manure; and M, soil treated with only organic manure. The number (1, 2, and 3) following the treatment indicates replicates of each treatment.
Similarly to AMF community composition variations, in the Monte Carlo permutation test (Supplementary Table S4), organic manure was the most significant factor influencing the ethyl acetate-extracted DOM composition (r2 = 0.8608, P = 0.002) while N and P fertilizers did not show significant impact. These results showed that the AMF community and the ethyl acetate-extracted DOM were similar in distribution and the application of organic manure was the key factor changing the AMF community composition and the ethyl acetate-extracted DOM composition rather than the application of chemical fertilizer.
Relationship between AMF Community Composition and Ethyl Acetate-Extracted DOM Composition
Ethyl acetate-extracted DOM of rhizosphere soils significantly (r = 0.5292, P = 0.002) correlated with the AMF community determined by Mantel tests. In total, 59 organic compounds were identified by GC–MS (Supplementary Figure S2), among those, 11 organic compounds were identified as being significantly associated with the AMF community (Table 2), including 6-tert-butyl-4-methylphenol (r = 0.284, P < 0.05), pentamethyl phenyl (r = 0.402, P < 0.05), octadecanoic acid (r = 0.479, P < 0.05), 3-methylbiphenyl (r = 0.330, P < 0.05), 2-nonadecanone O-methyloxime (r = 0.390, P < 0.05), 1,4-dimethylnaphthalene (r = 0.534, P < 0.05), 2-ethylnaphthalene (r = 0.736, P < 0.01), octadecane (r = 0.344, P < 0.05), hexacosane (r = 0.358, P < 0.01), 2,6,10-trimethyltetradecane (r = 0.468, P < 0.01) and 3-(2, 6, 6 -trimethyl-1-cyclohexene-1-yl) allylalcohol (r = 0.555, P < 0.01). Most of the 11 related molecules were long-chain alkanes and benzenes. Meanwhile, octadecanoic acid and 3-(2, 6, 6-trimethyl-1-cyclohexene-1-yl) allylalcohol were not detected in the NPKM treatment, while 1, 4-dimethylnaphthalene was only detected in the NK treatment. However, 2-ethylnaphthalene was not detected in the NK treatment. In addition, neither octadecane nor hexacosane was detected in the +M (M and NPKM) treatments, but both were detected in the -M (NK and NPK) treatments. This showed that the application of different fertilizers drove the variation in the AMF community by changing the composition of the organic matter.
TABLE 2. Mantel correlations between the identified ethyl acetate-extracted DOM and the AMF community.
Pearson correlations between the identified ethyl acetate-extracted DOM and AMF genera showed that two types of organic compounds were correlated with Glomus and one organic compound was correlated with Rhizophagus (Figure 4). 2-ethylnaphthalene (r = -0.6159, P < 0.05) and 2, 6, 10-trimethyltetradecane (r = -0.7106, P < 0.001) were negatively correlated with the relative abundance of Glomus and 3-methylbiphenyl (r = 0.6039, P < 0.05) was positively correlated with the relative abundance of Rhizophagus.
Differences in AMF Co-Occurrence between Treatments with and without Organic Amendment
Based on the aforementioned results, we divided the four treatments into two treatments (+M and -M). Hence, we used the two treatments for network analysis. From Figure 5, it is apparent that the results of the network analysis were observably different between the +M treatments (NPKM and M) compared to the -M treatments (NK and NPK). In the -M treatment, the network had 36 nodes and 80 edges, and the modularity was 0.540 with four communities, whilst for the +M treatment, the network presented 35 nodes and 86 edges, and the modularity was 0.438 with five communities (Figure 5). Here, we arbitrarily defined the node having more than seven edges as the network hub, which was active in mediating interspecies interactions. The number of active hubs in the +M treatment was more than in the -M treatment, indicating that +M had more active species interacting with other species within the AMF community. All the hubs were classified as the order of Glomerales (Supplementary Table S5). In addition, more modules were found in the +M treatment. Meanwhile, three OTUs serving as hubs were shared by both networks, including OTU9, OTU31, and OTU217. It was possible that these three OTUs were active for mediating AMF species interaction and would not change their roles in the network regardless of fertilization regimes. The three OTUs belonged to different branches of the phylogenetic tree (Figure 6). OTU9 and OTU 217 belonged to Glomus and OTU31 belonged to Claroideoglomus.
FIGURE 5. Network analysis of the -M and +M treatments based on Pearson correlations. The different color nodes belong to different modules. -M: the treatments without organic manure with modularity resolution of 0.540 and 4 communities; +M: the treatments with organic manure with modularity resolution of 0.438 and five communities.
FIGURE 6. Neighbor-joining phylogenetic tree of representative sequences of OTUs that are in both of the two networks obtained in this study and referenced sequences from the NCBI database.
The majority of sequences that were detected using the primer set AMV4.5NF/AMDGR were not from Glomeromycota, even though it was reported to be acceptable in several previous studies (Lumini et al., 2010; Lin et al., 2012; Dai et al., 2014; Bainard et al., 2015). We obtained nearly 170000 AMF sequences, forming 42 OTUs in all treatments. The rarefaction curves for all four treatments reached a plateau, indicating that the sequencing depth was sufficient to analyze the results (Supplementary Figure S1). In addition, the number of OTUs obtained were similar to that found in other pyrosequencing efforts from grasslands (Lumini et al., 2010; Dumbrell et al., 2011), suggesting that maize rhizosphere soils also contain abundant AMF resources compared to other ecosystems. In contrast to other reports (Oehl et al., 2005; Hijri et al., 2006; Wang et al., 2008), high-throughput sequencing data gave us a comprehensive insight into the AMF community in the agro-ecosystem in this study.
In spite of the fact that many articles have reported that the application of organic manures had negative impacts on AMF diversity (Sainz et al., 1998; Jacquot et al., 2000; Gryndler et al., 2008), several other studies have highlighted the positive influence of organic manure application on AMF population and diversity (Douds and Reider, 2003; Borie et al., 2008). Further, Albertsen et al. (2006) reported that AMF in general seemed to thrive in soil amended with organic matter. According to our results showed by ACE index and rarefaction curves (Figure 1; Supplementary Figure S1), the richness of AMF was increased significantly in the M treated soil, but it was not significantly different between the soils in +M and -M. The organic fertilizers have been shown to promote Glomus species in agricultural soils (Gryndler et al., 2006; Vestberg et al., 2011). As the soil treated with both chemical and organic fertilizers slightly, but not significantly, increased the AMF richness and diversity compared to soils with chemical-only fertilizer (Figure 1), these different responses of AMF to organic matter and alterations in AMF in relation to dosage of organic amendments were in accordance with Alguacil et al. (2011). Though, our results showed no significant difference in the α-diversity (ACE and Shannon indexes) between the -M and +M treatments according to the t-test (Figure 1), the AMF between-habitat diversity (β-diversity) was significantly separated between them (Figure 2), indicating that the changes in α-diversity and β-diversity should not occur at the same time (Van Diepen et al., 2011; Xiong et al., 2014). Some researchers reported that the change of community structure was typically correlated with shifts in functional behavior (Van Diepen et al., 2011; Fierer et al., 2013; Griffiths and Philippot, 2013). According to our results, significant difference was observed in the community structure between -M and +M treatments, but not in the richness. Therefore, these results implied that the AMF functioning probably responds first to the application of exogenous substances and the changes in microbial richness may frequently take a longer time in the maize rhizosphere soil.
P fertilizer was the main cause of changes in the soil AMF community according to the previous studies (Lin et al., 2012). However, in this study, the organic manure showed the utmost influence on the changes in AMF community composition instead of P fertilizer (Figure 2; Supplementary Table S2). Meanwhile, we determined that the variation in ethyl acetate-extracted DOM composition of the maize rhizosphere soils was predominantly impacted by organic amendment (Figure 3; Supplementary Table S4). It should be interesting to speculate that the application of organic manure was a key factor driving changes in the composition of AMF community by altering the ethyl acetate-extracted DOM composition in the maize rhizosphere soil.
Dissolved organic matter is the main substrate and energy source for microbes. The application of organic manure most likely changed the composition of the organic matter in the maize rhizosphere first, and later, the changes in organic matter composition might affect the growth and composition of AMF. There is growing evidence that the microbial composition is driven to a large extent by deterministic selection through environmental factors (Stegen et al., 2012; Valentín-Vargas et al., 2012) and the chemical nature of soil carbon drive the structure and functioning of soil microbial communities (Ng et al., 2014). Some studies have shown that AMF can directly take advantage of simple organic matter (Govindarajulu et al., 2005; Jin et al., 2005). Therefore, as biotrophic, some species of AMF that were efficient in using organic matter acquired more resources to grow and reproduce, and thereby became the dominant species after fertilization. From another perspective, root exudates are the most important sources for rhizosphere microorganisms. Yuan et al. (2016) reported that fungi out-competed bacteria in utilizing the root exudates. It is widely reported that the different fertilization regimes could change the secretion of root exudates (Carvalhais et al., 2011; Yoneyama et al., 2013; Kumar et al., 2016), which could alter the AMF community composition as well. Based on the above reasons, the AMF community composition in the maize rhizosphere soil was expectedly changed by the shift of fertilization regime, which was in accordance with the finding that organic fertilization alters the community composition of root associated fungi in Pisum sativum (Yu et al., 2013).
Microbes may consume or produce DOM, but the prevalent DOM composition may also (and in turn) select for a certain microbial community (Osterholz et al., 2016). It is widely supported that the specific carbon compounds trigger the growth of certain bacterial strains (Gómez-Consarnau et al., 2012), but also acknowledges the “generalists” – a functional redundancy of metabolic capabilities can lead to a minor influence of DOM quality on bacteria structure (Mou et al., 2007). In our study, the Mantel test did not allow resolving any causal relationships between AMF and ethyl acetate-extracted DOM composition; rather we provided a viewpoint of the effect of ethyl acetate-extracted DOM on AMF composition in the maize rhizosphere soil. In our study, 11 types of organic compounds were significantly correlated with AMF community composition. This indicated that not all organic matters affected the AMF community composition in the maize rhizosphere and that only some AMF were sensitive in their response to organic matter composition (Figure 4). Similarly, some studies (Yuan et al., 2012; Raza et al., 2015) have reported that organic compounds such as naphthyl, alkanes, and benzenes could affect the growth and spore germination of one species of fungi (Fusarium oxysporum).
The interaction of OTUs within habitat is one kind of internal representations of changes in AMF community structure. As caused by different fertilization regimes, the results (showed in Figure 5; Table 2) were enough to exposit that changes of the AMF community composition was correlated with the ethyl acetate-extracted DOM constituent. It is well-established that some rhizosphere compounds can regulate the interaction between microbes (van Dam and Bouwmeester, 2016; Zhu et al., 2016). Song et al. (2016) reported that ethyl acetate-extracted DOM would mediate microbial interaction in rhizosphere soils. Thus, it would be reasonable to propose that the changed constituent of ethyl acetate-extracted DOM attributed from the variation of fertilization would also lead to shifts in soil AMF co-occurrence patterns, and conducting a network analysis would be worthwhile to display the difference and similarity of AMF interactions between the soils with and without organic fertilizer (+M vs. -M). In our network analysis, the more active nodes were detected in the +M treatments than in the -M treatments (Figure 5). The application of organic manure activated more species of AMF involved in the maize rhizosphere soil. Soils amended with organic manure enhanced soil fertility (Chaudhry et al., 2012) and activated a more diverse group of soil microbes as compared with soils conventionally applied with inorganic fertilizers (Ai et al., 2015). Yu et al. (2013) reported that G. mosseae and G. caledonium which belonged to Glomus were relatively more abundant with increasing amount of the organic fertilizer applied, whereas Paraglomus sp. was more abundant in treatments with low dosage of organic fertilizer. The fact that different AMF had contrasting response to organic matter might indicate certain differential life strategies of AMF. Hence, some AMF species that were not active originally became active in mediating AMF species interaction after the application of organic manure. These were the major reasons of changes in AMF community composition in the maize rhizosphere with different fertilizer regimes. Only three OTUs were same in both networks: OTU9, OTU31, and OTU217. These OTUs were likely to occupy broad ecological niches and were always active for mediating AMF species interaction in the maize rhizosphere, therefore, the changes in the environment caused by fertilizations did not affect the activity of these three OTUs. In addition, they all belonged to Glomerales (Figure 6; Supplementary Table S5), which is in agreement with other studies (Oehl et al., 2003; Hijri et al., 2006; Beauregard et al., 2013) showing that Glomerales is widespread and is the dominant order in the soil. Moreover, at a higher phylogenetic resolution, OTU9 and OTU217 belonged to the genus Glomus. Glomus can easily survive and reproduce via mycelium, mycorrhizal spores or fragments. Therefore, Glomus could be more resistant and resilient to disturbances in the ecological environment thus often to play an important role in executing ecological function (Giovannetti et al., 1999; Daniell et al., 2001), such as mediating AMF species interaction.
Fertilization regimes have impacts on AMF community composition and organic matter composition in the maize rhizosphere soil. A key factor in changing the AMF community composition and the organic matter composition was the application of organic fertilizer rather than chemical N and P fertilizers. The organic matter composition and the AMF community composition were significantly correlated in the maize rhizosphere. The organic amendment also resulted in the activation of more OTUs and more complex interspecific interactions in the maize rhizosphere soil, and this active role was commonly found to be played by Glomerales. However, to gain a better understanding of manipulating the AMF community, examinations on how and to what extent the rhizospheric DOMs impact the AMF community following the addition of soil amendments, which we predicted to impact the community in this study, should be conducted in the future.
CZ, NL, and SG proposed and organized the overall project. CZ performed the majority of the experiments. JG and MW gave assistance in lab work and laboratory analyses. CZ and NL wrote the main manuscript text. SG and QS contributed insightful discussions. All authors reviewed the manuscript.
Financial supports from the National Basic Research Program of China (2013CB127403, 2015CB150500) and China Postdoctoral Science Foundation (2016T90473) are acknowledged.
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.
We are grateful to many graduate students and staff involved in maintaining the field plots and collecting soil samples who are not listed as co-authors. We also thank Dr. Waseem Raza from Nanjing Agricultural University and professor Zed Rengel from University of Western Australia for polishing the language of this paper.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.01840/full#supplementary-material
Ai, C., Liang, G., Sun, J., Wang, X., He, P., Zhou, W., et al. (2015). Reduced dependence of rhizosphere microbiome on plant-derived carbon in 32-year long-term inorganic and organic fertilized soils. Soil Biol. Biochem. 80, 70–78. doi: 10.1016/j.soilbio.2014.09.028
Albertsen, A., Ravnskov, S., Green, H., Jensen, D. F., and Larsen, J. (2006). Interactions between the external mycelium of the mycorrhizal fungus glomus intraradices and other soil microorganisms as affected by organic matter. Soil Biol. Biochem. 38, 1008–1014. doi: 10.1016/j.soilbio.2005.08.015
Alguacil, M. M., Torrecillas, E., Caravaca, F., Fernández, D. A., Azcón, R., and Roldán, A. (2011). The application of an organic amendment modifies the arbuscular mycorrhizal fungal communities colonizing native seedlings grown in a heavy-metal-polluted soil. Soil Biol. Biochem. 43, 1498–1508. doi: 10.1016/j.soilbio.2011.03.026
Bainard, L. D., Dai, M., Gomez, E. F., Torres-Arias, Y., Bainard, J. D., Sheng, M., et al. (2015). Arbuscular mycorrhizal fungal communities are influenced by agricultural land use and not soil type among the chernozem great groups of the Canadian Prairies. Plant Soil 387, 351–362. doi: 10.1007/s11104-014-2288-1
Bastian, M., Heymann, S., and Jacomy, M. (2009). “Gephi: an open source software for exploring and manipulating networks,” in Proceedings of the International AAAI Conference on Weblogs and Social Media, Menlo Park, CA. 361–362.
Beauregard, M., Gauthier, M.-P., Hamel, C., Zhang, T., Welacky, T., Tan, C., et al. (2013). Various forms of organic and inorganic P fertilizers did not negatively affect soil-and root-inhabiting AM fungi in a maize–soybean rotation system. Mycorrhiza 23, 143–154. doi: 10.1007/s00572-012-0459-6
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Carvalhais, L. C., Dennis, P. G., Fedoseyenko, D., Hajirezaei, M. R., Borriss, R., and von Wiren, N. (2011). Root exudation of sugars, amino acids, and organic acids by maize as affected by nitrogen, phosphorus, potassium, and iron deficiency. J. Plant Nutr. Soil Sci. 174, 3–11. doi: 10.1002/jpln.201000085
Celik, I., Ortas, I., and Kilic, S. (2004). Effects of compost, mycorrhiza, manure and fertilizer on some physical properties of a Chromoxerert soil. Soil Till. Res. 78, 59–67. doi: 10.1016/j.still.2004.02.012
Chaudhry, V., Rehman, A., Mishra, A., Chauhan, P. S., and Nautiyal, C. S. (2012). Changes in bacterial community structure of agricultural land due to long-term organic and chemical amendments. Microbiol. Ecol. 64, 450–460. doi: 10.1007/s00248-012-0025-y
Clemente, R., Walker, D. J., and Bernal, M. P. (2005). Uptake of heavy metals and as by Brassica juncea grown in a contaminated soil in Aznalcollar (Spain): the effect of soil amendments. Environ. Pollut. 138, 46–58. doi: 10.1016/j.envpol.2005.02.019
Dai, M., Hamel, C., Bainard, L. D., Arnaud, M. S., Grant, C. A., Lupwayi, N. Z., et al. (2014). Negative and positive contributions of arbuscular mycorrhizal fungal taxa to wheat production and nutrient uptake efficiency in organic and conventional systems in the Canadian prairie. Soil Biol. Biochem. 74, 156–166. doi: 10.1016/j.soilbio.2014.03.016
Daniell, T., Husband, R., Fitter, A., and Young, J. (2001). Molecular diversity of arbuscular mycorrhizal fungi colonising arable crops. FEMS Microbiol. Ecol. 36, 203–209. doi: 10.1111/j.1574-6941.2001.tb00841.x
Dumbrell, A. J., Ashton, P. D., Aziz, N., Feng, G., Nelson, M., Dytham, C., et al. (2011). Distinct seasonal assemblages of arbuscular mycorrhizal fungi revealed by massively parallel pyrosequencing. New Phytol. 190, 794–804. doi: 10.1111/j.1469-8137.2010.03636.x
Fierer, N., Ladau, J., Clemente, J. C., Leff, J. W., Owens, S. M., Pollard, K. S., et al. (2013). Reconstructing the microbial diversity and function of pre-agricultural tallgrass prairie soils in the United States. Science 342, 621–624. doi: 10.1126/science.1243768
Gómez-Consarnau, L., Lindh, M. V., Gasol, J. M., and Pinhassi, J. (2012). Structuring of bacterioplankton communities by specific dissolved organic carbon compounds. Environ. Microbiol. 14, 2361–2378. doi: 10.1111/j.1462-2920.2012.02804.x
Govindarajulu, M., Pfeffer, P. E., Jin, H., Abubaker, J., Douds, D. D., Allen, J. W., et al. (2005). Nitrogen transfer in the arbuscular mycorrhizal symbiosis. Nature 435, 819–823. doi: 10.1038/nature03610
Gryndler, M., Larsen, J., Hršelová, H., Řezáèová, V., Gryndlerová, H., and Kubát, J. (2006). Organic and mineral fertilization, respectively, increase and decrease the development of external mycelium of arbuscular mycorrhizal fungi in a long-term field experiment. Mycorrhiza 16, 159–166. doi: 10.1007/s00572-005-0027-4
Gryndler, M., Sudová, R., Püschel, D., Rydlová, J., Janoušková, M., and Vosátka, M. (2008). Cultivation of high-biomass crops on coal mine spoil banks: can microbial inoculation compensate for high doses of organic matter? Bioresour. Technol. 99, 6391–6399. doi: 10.1016/j.biortech.2007.11.059
He, J. Z., Zheng, Y., Chen, C. R., He, Y. Q., and Zhang, L. M. (2008). Microbial composition and diversity of an upland red soil under long-term fertilization treatments as revealed by culture-dependent and culture-independent approaches. J. Soil Sediment. 8, 349–358. doi: 10.1007/s11368-008-0025-1
Hijri, I., Sýkorová, Z., Oehl, F., Ineichen, K., Mäder, P., Wiemken, A., et al. (2006). Communities of arbuscular mycorrhizal fungi in arable soils are not necessarily low in diversity. Mol. Ecol. 15, 2277–2289. doi: 10.1111/j.1365-294X.2006.02921.x
Jacquot, E., Van Tuinen, D., Gianinazzi, S., and Gianinazzi-Pearson, V. (2000). Monitoring species of arbuscular mycorrhizal fungi in planta and in soil by nested PCR: application to the study of the impact of sewage sludge. Plant Soil 226, 179–188. doi: 10.1023/A:1026475925703
Jin, H., Pfeffer, P., Douds, D., Piotrowski, E., Lammers, P., and Shachar-Hill, Y. (2005). The uptake, metabolism, transport and transfer of nitrogen in an arbuscular mycorrhizal symbiosis. New Phytol. 168, 687–696. doi: 10.1111/j.1469-8137.2005.01536.x
Li, Z. F., Yang, Y. Q., Xie, D. F., Zhu, L. F., Zhang, Z. G., and Lin, W. X. (2012). Identification of autotoxic compounds in fibrous roots of Rehmannia (Rehmannia glutinosa Libosch.). PLoS ONE 7:e28806. doi: 10.1371/journal.pone.0028806
Lin, X., Feng, Y., Zhang, H., Chen, R., Wang, J., Zhang, J., et al. (2012). Long-term balanced fertilization decreases arbuscular mycorrhizal fungal diversity in an arable soil in North China revealed by 454 pyrosequencing. Environ. Sci. Technol. 46, 5764–5771. doi: 10.1021/es3001695
Lueders, T., Kindler, R., Miltner, A., Friedrich, M. W., and Kaestner, M. (2006). Identification of bacterial micropredators distinctively active in a soil microbial food web. Appl. Environ. Microbiol. 72, 5342–5348. doi: 10.1128/AEM.00400-06
Lumini, E., Orgiazzi, A., Borriello, R., Bonfante, P., and Bianciotto, V. (2010). Disclosing arbuscular mycorrhizal fungal biodiversity in soil through a land-use gradient using a pyrosequencing approach. Environ. Microbiol. 12, 2165–2179.
Ng, E. L., Patti, A. F., Rose, M. T., Schefe, C. R., Wilkinson, K., Smernik, R. J., et al. (2014). Does the chemical nature of soil carbon drive the structure and functioning of soil microbial communities? Soil Biol. Biochem. 70, 54–61. doi: 10.1016/j.soilbio.2013.12.004
Oehl, F., Sieverding, E., Ineichen, K., Mäder, P., Boller, T., and Wiemken, A. (2003). Impact of land use intensity on the species diversity of arbuscular mycorrhizal fungi in agroecosystems of Central Europe. Appl. Environ. Microbiol. 69, 2816–2824. doi: 10.1128/AEM.69.5.2816-2824.2003
Oehl, F., Sieverding, E., Ineichen, K., Ris, E. A., Boller, T., and Wiemken, A. (2005). Community structure of arbuscular mycorrhizal fungi at different soil depths in extensively and intensively managed agroecosystems. New Phytol. 165, 273–283. doi: 10.1111/j.1469-8137.2004.01235.x
Osterholz, H., Singer, G., Wemheuer, B., Daniel, R., Simon, M., Niggemann, J., et al. (2016). Deciphering associations between dissolved organic molecules and bacterial communities in a pelagic marine system. ISME J. 10, 1717–1730. doi: 10.1038/ismej.2015.231
Pramanik, M. H. R., Nagai, M., Asao, T., and Matsui, Y. (2000). Effects of temperature and photoperiod on phytotoxic root exudates of cucumber (Cucumis sativus) in hydroponic culture. J. Chem. Ecol. 26, 1953–1967. doi: 10.1023/A:1005509110317
Raza, W., Yuan, J., Wu, Y. C., Rajer, F., Huang, Q. W., and Shen, Q. R. (2015). Biocontrol traits of two Paenibacillus polymyxa strains SQR-21 and WR-2 in response to fusaric acid, a phytotoxin produced by Fusarium species. Plant Pathol. 64, 1041–1052. doi: 10.1111/ppa.12354
Sainz, M., Taboada-Castro, M., and Vilarino, A. (1998). Growth, mineral nutrition and mycorrhizal colonization of red clover and cucumber plants grown in a soil amended with composted urban wastes. Plant Soil 205, 85–92. doi: 10.1023/A:1004357330318
Song, Y., Zhu, C., Raza, W., Wang, D. S., Huang, Q. W., Guo, S. W., et al. (2016). Coupling of the chemical niche and microbiome in the rhizosphere: implications from watermelon grafting. Front. Agric. Sci. Eng. 3, 249–262. doi: 10.15302/J-FASE-2016105
Stegen, J. C., Lin, X., Konopka, A. E., and Fredrickson, J. K. (2012). Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 6, 1653–1664. doi: 10.1038/ismej.2012.22
Valentín-Vargas, A., Toro-Labrador, G., and Massol-Deya, A. A. (2012). Bacterial community dynamics in full-scale activated sludge bioreactors: operational and ecological factors driving community assembly and performance. PLoS ONE 7:e42524. doi: 10.1371/journal.pone.0042524
Van Der Heijden, M. G., Bardgett, R. D., and Van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
Van Diepen, L. T., Lilleskov, E. A., and Pregitzer, K. S. (2011). Simulated nitrogen deposition affects community structure of arbuscular mycorrhizal fungi in northern hardwood forests. Mol. Ecol. 20, 799–811. doi: 10.1111/j.1365-294X.2010.04969.x
Vestberg, M., Kahiluoto, H., and Wallius, E. (2011). Arbuscular mycorrhizal fungal diversity and species dominance in a temperate soil with long-term conventional and low-input cropping systems. Mycorrhiza. 21, 351–361. doi: 10.1007/s00572-010-0346-y
Walker, D. J., Clemente, R., and Bernal, M. P. (2004). Contrasting effects of manure and compost on soil pH, heavy metal availability and growth of Chenopodium album L. in a soil contaminated by pyritic mine waste. Chemosphere 57, 215–224. doi: 10.1016/j.chemosphere.2004.05.020
Wander, M., Traina, S., Stinner, B., and Peters, S. (1994). Organic and conventional management effects on biologically active soil organic matter pools. Soil Sci. Soc. Am. J. 58, 1130–1139. doi: 10.2136/sssaj1994.03615995005800040018x
Wang, Y. Y., Vestberg, M., Walker, C., Hurme, T., Zhang, X., and Lindström, K. (2008). Diversity and infectivity of arbuscular mycorrhizal fungi in agricultural soils of the Sichuan province of mainland China. Mycorrhiza 18, 59–68. doi: 10.1007/s00572-008-0161-x
Xiong, J., Sun, H., Peng, F., Zhang, H., Xue, X., Gibbons, S. M., et al. (2014). Characterizing changes in soil bacterial community structure in response to short-term warming. FEMS Microbiol. Ecol. 89, 281–292. doi: 10.1111/1574-6941.12289
Yoneyama, K., Xie, X. N., Kisugi, T., Nomura, T., and Yoneyama, K. (2013). Nitrogen and phosphorus fertilization negatively affects strigolactone production and exudation in sorghum. Planta 238, 885–894. doi: 10.1007/s00425-013-1943-8
Yuan, J., Raza, W., Shen, Q. R., and Huang, Q. W. (2012). Antifungal activity of Bacillus amyloliquef aciens NJN-6 volatile compounds against Fusarium oxysporum f. sp. cubense. Appl. Environ. Microbiol. 78, 5942–5944. doi: 10.1128/AEM.01357-12
Zhong, W., Gu, T., Wang, W., Zhang, B., Lin, X., Huang, Q., et al. (2010). The effects of mineral fertilizer and organic manure on soil microbial community and diversity. Plant Soil 326, 511–522. doi: 10.1007/s11104-009-0099-6
Keywords: fertilization regimes, organic composition, Illumina MiSeq sequencing, AMF community composition, network analysis
Citation: Zhu C, Ling N, Guo J, Wang M, Guo S and Shen Q (2016) Impacts of Fertilization Regimes on Arbuscular Mycorrhizal Fungal (AMF) Community Composition Were Correlated with Organic Matter Composition in Maize Rhizosphere Soil. Front. Microbiol. 7:1840. doi: 10.3389/fmicb.2016.01840
Received: 14 July 2016; Accepted: 01 November 2016;
Published: 16 November 2016.
Edited by:Jeanette M. Norton, Utah State University, USA
Reviewed by:Hui Li, Institute of Applied Ecology (CAS), China
Jian Yang, China University of Geosciences, China
Copyright © 2016 Zhu, Ling, Guo, Wang, Guo and Shen. 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.
*Correspondence: Ning Ling, email@example.com