Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 45134.
Published online 2017 March 22. doi:  10.1038/srep45134
PMCID: PMC5361092

The influence of environmental factors on communities of arbuscular mycorrhizal fungi associated with Chenopodium ambrosioides revealed by MiSeq sequencing investigation


Arbuscular mycorrhizal fungi (AMF) affect multiple ecosystem functions and processes, the assemblages of which vary across ecosystems. However, the influences of environmental factors on AMF communities which may shape these communities are still largely unknown. In this study, AMF communities from roots and rhizosphere soils of Chenopodium ambrosioides in different natural soils were investigated. The root habitat showed significantly smaller numbers of OTUs and lower community richness compared to the rhizosphere soil habitat. Most OTUs in the root habitat were shared by the soil habitat from the same sampling site, indicating that rhizosphere soils represent a pool of AMF species, a fraction of which is recruited by plants. Most of the AMF in root habitats were Glomeraceae, suggesting recruitment preferences of AMF by plants. The relative contributions of environmental factors to explain variations in AMF community composition and phylogenetic structure were assessed. The results revealed soil properties predominantly explained the variation, followed by geographic and climate parameters which explained a small fraction independently, while the host plant showed few explanations. Overall, our results indicated that soil and root habitats as well as soil characters, especially pH, nitrogen and micronutrients (Zn and Cu) affected AMF communities significantly.

The arbuscular mycorrhizal fungi (AMF, Glomeromycota) can form mutualistic symbioses with most terrestrial plants1. AM mutualism is an ancient plant–microbe interaction, which evolved very early in the evolution of terrestrial plants and is considered to have been crucial for the successful colonization of land by plants2,3,4. This mutualistic symbiosis enhances plant nutrient acquisition from soil and enhances tolerance to diverse biotic and abiotic stresses5,6. For example, mycorrhizal plants exhibit greater uptake of water and mineral nutrients than non-mycorrhizal plants1; AMF protect host plants against heavy metals (HMs) and enhance their HM tolerance7,8.

Environmental parameters and host plant species could affect AMF communities in terms of their composition and distribution9,10,11,12,13. Cheng et al. reported that a higher soil phosphorus (P) level is associated with lower mycorrhizal root colonization rates and lower AMF diversity14. Various AMF communities colonize the same host plant depending on the soil P concentration, according to Gosling et al.15. Thus mineral nutrients are important environmental factors. Furthermore, soil pH has been suggested to have a marked effect on AMF communities in agroecosystems and crops10. It has also been revealed that soil texture rather than P may influence AMF community structure16. Several studies have focused on the influence of HMs on AMF communities17,18. However, most of these studies were confined to local geographical scales and a limited set of soil characters. Thus, it would be helpful to consider a greater number of parameters of environment from a wider geographical area to evaluate the effects of environmental factors on AMF communities.

Many case studies have relied on microscopic identification of AMF spores from the field19,20,21. However, some AMF species sporulate infrequently or not at all in the field, and so might not be identified22,23. Moreover, spore density in soils does not necessarily reflect the AMF population colonizing the plant roots22,24. Recently, molecular analyses of DNA extracted from soils or roots have been carried out to evaluate AMF communities25,26. High-throughput sequencing technologies, such as 454 and Miseq amplicon pyrosequencing, facilitate efficient characterization of AMF communities by sequencing specific amplicons27,28. These technologies enable more detailed analyses of root or soil AMF communities.

Here, we investigated AMF communities in the rhizosphere soils and the roots of Chenopodium ambrosioides, a common plant species with a broad distribution in China, to identify the factors that affect AMF communities among different areas and habitats. To accomplish this, we collected root and rhizosphere soil samples from different geographical regions. Using high-throughput sequencing of a fragment of the small subunit (SSU) rRNA gene, we evaluated the composition of AMF communities in both root and rhizosphere soil habitats and investigated the relationships between the compositions of AMF communities and soil physical, chemical, and biological characteristics.


Soil chemical characteristics

The majority of the physicochemical characteristics differed significantly among soil types (Table 1). All soil samples were neutral (pH 7.04–7.62), with the exception of GP (pH 5.62). The XC sample had the highest NO3-N level. The TL and WH samples had considerably higher OM and total N levels but lower NH4+-N level than the other samples. Although lower total P, K, and Zn levels were detected in SS, TL, and WH compared with DN, GP, and XC, the available P, K and Zn levels were similar among all of them. Both the available and total Mn concentrations in SS, TL, and WH were significantly lower than those of DN, GP, and XC. The location, soil physical properties, plant parameters and climate characters of all sampling sites are shown in Supplementary Table S1.

Table 1
Chemical properties of sampled soils.

Overall sequencing information and taxonomic richness

A total of 1,465,650 paired-end sequences were obtained from the 36 libraries using the AMV4.5NF/AMDGR primer set. These sequences were overlapped to obtain high-quality tag sequences, and the dominant length distribution was 200−250 bp (>93%), with an average length of 218 bp. The details of the tag sequences obtained from each of the 36 samples are provided in Supplementary Table S2.

The taxonomic distributions of the obtained sequences in rhizosphere soil and root samples are shown in Supplementary Fig. S1. Fungal sequences other than Glomeromycota were detected. However, the majority of fungal sequences amplified with AMV4.5NF/AMDGR belonged to Glomeromycota (705,323 sequences, corresponding to 60.1% of the total). The second and third largest groups were Basidiomycota (15.2%) and Chytridiomycota (15.1%), respectively. The lowest numbers of Glomeromycota sequences were obtained for root and soil samples of XC (R-XC and S-XC), while the percentage of Glomeromycota sequences obtained from other samples was >40%. The Glomeromycota sequences were extracted for further analyses. Specifically, pyrosequencing of XC3 and WH3 failed because few Glomeromycota sequences were obtained, and so these two samples were excluded from further analyses.

AMF community richness and diversity

All of the rarefaction curves tended to reach saturation (Supplementary Fig. S2), revealing that the data volume of sequenced reads was sufficient to detect the majority of sequence types. Marked variation in total OTU number among the samples was suggested by the rarefaction curve. The samples from rhizosphere soil habitats had a larger number of OTUs (20.5–55.5) compared to those from root habitats (4.0–19.7) (Table 2). Consistently, the richness indices, Chao1 and ACE, showed that the samples from root habitats had the lowest number of AMF OTUs, while the samples from rhizosphere soil habitats had the highest number, with significant differences between habitats. However, no significant differences of Shannon diversity were detected between the root and rhizosphere soil habitats.

Table 2
Richness estimators and diversity indices at a 97% identity threshold.

AMF community composition

Significant variation in AMF community composition at the genus level was detected among both root and soil samples (Fig. 1). Sequences that could be classified were affiliated with nine AMF genera, while those that could not be classified were assigned as others. Glomus and Rhizophagus were the most abundant genera in all root samples, but their relative levels differed. A greater number of genera were present in samples from rhizosphere soil habitats; the most abundant genera were Glomus, Funneliformis, and Claroideoglomus.

Figure 1
Abundance percentages of AMF genera for all samples.

A PCoA based on the OTU composition is shown in Fig. 2a. Of the variation in the AMF communities, 26.6% and 13.3% could be explained by the first and second principal components, respectively. The samples from root and soil habitats at the same site clustered together (Fig. 2a). This was confirmed by the hierarchical cluster analysis, which showed clusters of samples from the same site (Fig. 2b).

Figure 2
Principal coordinates analysis (a) and hierarchical clustering (b) with Bray–Curtis distances of AMF communities. R, root samples; S, soil samples.

Indicator species analyses were used to identify specific OTUs associated with the rhizosphere soil or root habitat. Sixteen OTUs were found in rhizosphere soil samples. Among them, 11 OTUs were specific for the rhizosphere soil habitat (Table 3). However, no OTUs associated with the root habitat were detected.

Table 3
OTUs associated with the rhizosphere soil habitat.

In total, 210 OTUs were found in the AMF communities based on 97% species identity (Fig. 3). The phylogenetic placement of OTUs in different habitats was determined using the phylogenetic tree (Fig. 3). In general, the OTUs in rhizosphere soils were evenly scattered in the phylogenetic tree. However, most of AMF detected in root habitats were Glomeraceae while few AMF belonging to Acaulosporacea, Diversisporaceae or Gigasporaceae were found. Claroideoglomeraceae, Paraglomeraceae were identified in the roots of DN, GP and TL. Although the assemblages of AMF communities were markedly different, most of the OTUs in root samples were shared by the corresponding soil samples from the same site (Fig. 4).

Figure 3
Phylogenetic tree of OTUs in different habitats.
Figure 4
Venn diagram depicting OTUs that are shared or unique to root and soil samples.

Relationship between AMF community structure and environmental factors

Relationships of AMF richness (indicated by the ACE index) and phylogenetic diversity (revealed by the Faith’s index) with plant, soil and climate variables were analyzed (Supplementary Table S3). We found that the soil NH4+-N level correlated negatively with the AMF richness and phylogenetic diversity (Kendall test, P < 0.05). RDA, the multivariate analysis based on constrained ordination, was applied to analyze the influence of soil properties on the AMF community composition in C. ambrosioides rhizospheres (Fig. 5). The first two RDA components explained 31.7% of the total variation. The results showed that total N (P = 0.007), pH (P = 0.029), Zn (P = 0.033), and Cu (P = 0.047) explained the AMF community in soil habitats most and had a significant correlation between each variable and the ordination scores (Fig. 5).

Figure 5
Redundancy analysis (RDA) of AMF communities and soil properties.

To further analyze the relative importance of soil characters, plant parameters, geographic distance, and climate properties in predicting AMF community composition and phylogenetic structure, a stepwise distance-based redundancy analysis (db-RDA) was carried out. In contrast to traditional RDA, an advantage of db-RDA is that it enables the inclusion of environmental factors and the testing of their interaction using non-Euclidean distance matrices. Totally, 71.6% and 60.8% of the variations of the soil AMF community composition and phylogenetic structure were explained by the whole set of the selected variables respectively (Fig. 6). A very large fraction of variations of AMF community composition (62.8%) and phylogenetic structure (41.0%) could be explained by soil properties, followed by geographic and climate variables (31.5% for community composition and 16.5% for phylogenetic structure respectively). The plant variables explained a relative small proportion of the variations in both community composition and phylogenetic structure.

Figure 6
The contributions of environmental factors to the variation of AMF communities based on a stepwise distance-based redundancy analysis (db-RDA).

Permutational multivariate analysis of variance (PerMANOVA) was also carried out to examine the relative importance of each single environmental factor to the AMF community composition and phylogenetic structure (Supplementary Table S4). pH showed weak associations with the AMF community composition and phylogenetic structure, while the soil available Cu content revealed significant associations with them. Besides, the PCNM vector (principal components of neighbor matrices, representing the spatial factors) and MAP (mean annual precipitation) showed significant associations to the AMF community composition. For the AMF community phylogenetic structure, soil total Cu, Zn, and Cd contents as well as MAT (mean annual temperature) were significantly correlated.


It is commonly accepted that more than 80% of terrestrial plants are colonized by AMF, with which they form associations29,30. The extraradical hyphal networks produced by AMF can alter plant physiology, enhance plant nutrient absorption and translocation, and increase plant resistance to HMs31,32,33. Root habitats showed a significantly smaller number of AMF OTUs compared to rhizosphere soil habitats, which was confirmed by the lower Chao1 and ACE richness indices of root samples. Similar results have previously been reported18,34. Besides, the PCoA based on the OTU composition showed that the samples from root and soil habitats at the same site clustered together, which was consistent with the hierarchical cluster analysis results, indicating similarities in the AMF community composition among root and soil samples from the same site and variation among samples from different sites. The indicator species analyses revealed 16 OTUs associated with rhizosphere soils while no OTUs with the root habitat were detected. Besides, most OTUs were shared by roots and rhizosphere soils from the same site, while few OTUs were shared among sites. These results indicated that AMF in rhizosphere soils could be considered to represent a pool of species, a fraction of which is recruited by plants35,36.

Xu et al. detected a strong influence of the plant community on AMF communities in soil, indicating preferences in plant-AMF associations37. These preferences have been observed on both local and global scale systems36,38,39, which might result from that host plants exhibited preferential allocation of photosynthates to more beneficial AMF partners40. Consistently, Saks et al. showed that root-colonizing AMF represent a phylogenetically clustered subset of AMF available in soil41. The conclusion is also proved in this study by the result that most of the AMF in root habitats were Glomeraceae. These results indicated that there might be niche preferences among AMF affecting AMF community composition34,42. Different to the study of AMF communities by Xu et al. which covered a wide range of vegetation types37, all soil samples were collected from rhizosphere of a single plant with similar size and age to reduce the plant factors affecting the AMF community in this study. Similar sampling strategy was also used for the study of AMF communities in semiarid Mediterranean soils34. We did not detect significant influences of the single plant on AMF communities in soils, indicating that the sampling strategy successfully reduce the biotic factors affecting the AMF distribution. Furthermore, by comparing influences of a single plant on AMF communities in roots and rhizosphere soils, we suggest that the effect of AMF communities in soils by plant communities might result from the different composition of plant communities in which each single plant shows the specific recruitment preference of AMF partners.

Glomus and Rhizophagus were the most abundant genera in root samples. Yang et al. reported that almost all of the sequences found in the roots of Elsholtzia splendens were Glomus species43. Long et al. revealed that although Ambispora, Kuklospora, and Glomus dominated in the rhizosphere soils of Phytolacca americana, only Glomus was detected in the roots44. Our results are in agreement with these previous studies. Glomus has been found to be dominant in various habitats, such as HM-polluted soils45, grassland soils46, and agricultural soils47. The Rhizophagus group is also a general AMF group that has been found in diverse host species and environments48,49,50. However, almost all of the fungi detected in R-DN and R-XC (root samples from DN and XC) were Glomus, the relative levels of which were much higher than those in soil habitats. On the contrary, the relative level of Glomus in R-GP (root samples from GP) was markedly lower than that in S-GP (soil samples from GP), while Rhizophagus was significantly more abundant in R-GP than S-GP. The Rhizophagus level was also higher in the root samples from TL and WH compared to the corresponding soil samples, while no obvious recruitment was detected for R-SS. These results indicated that the AMF recruitment preferences by C. ambrosioides differed among sampling sites. Several explanations for the differences in the AMF present in the roots and rhizosphere soils of the same plant have been proposed, including different AMF life strategies, differential sporulation dynamics, and seasonal changes in the AMF community50,51,52.

It has been shown that geographical distance influences AMF community at large spatial scales, especially in global-scale studies37. However, at local or landscape scales, soil abiotic factors are the key driver in shaping AMF community composition10,11. In our study, the geographic and climate parameters could independently explain a small fraction (<10%) of variances of AMF community composition and phylogenetic structure, while a large fraction of variances (>32%) could be explained independently by soil variables. These results revealed the complexity of factors regulating AMF communities, and the distribution patterns of AMF communities could not be completely explained by soil heterogeneity.

Despite high levels of Mn in soils of DN, GP and XC, no significant differences in the species diversity and richness of AMF communities were detected between these three samples and others. However, the richness indices of the R-GP and R-XC were markedly lower than those of the other samples, revealing that Mn affects the AMF community to some extent. These results are similar to Wei et al., which reported that root colonization and AMF diversity were negatively correlated with soil Mn concentration17. These results suggest also that other environmental factors may affect the AMF community and the effects of Mn on the AMF community may be confounded by those of other factors. Indeed, various factors influence the AMF community41. For example, Yang et al. showed that HM contamination is not the only soil parameter that affects the AMF community, and soil pH, Pb, Zn, Cd, and OM levels also have a great influence18; Wei et al. concluded that changes in AMF diversity and colonization are not solely attributed to soil Mn concentration, while soil properties, especially N concentration, were also closely related to it53.

In this study, pH was found to be a significant factor influencing AMF communities. Similarly, some studies concluded that pH is the key environmental factor shaping AMF communities54. Soil acidity is one of the most important drivers of microbial communities, particularly for AMF11,55. Soil pH can directly change the physiological status of indigenous AMF, alter their ecological niches, and it could also indirectly influences the AMF community by regulating soil nutrient bioavailability, and impacting the mobility and sorption of metals. Several studies have reported strong effects of soil Zn and Cu levels on AMF abundance and diversity in soils18,34,56. These nutrients play important roles in plant metabolic processes, and their uptake is influenced by AMF57,58. It has been shown that total N is related to changes in the composition of soil microbial communities59. Avio et al. showed N-fertilization was the main factor shaping AMF communities60 while Van Diepen et al. revealed N addition significantly altered the AMF community structure61. The lack of nitrogen could cause inhibition of sporulation62, while high N availability can change nutritional processes in AMF and alter the abundance of AMF phylotypes63. Besides, the N need of a plant could facilitate colonization by AMF. These features may explain the reasons why total N, pH, Zn, and Cu affect AMF community in C. ambrosioides rhizosphere.


The AMF communities of the roots and rhizosphere soils of C. ambrosioides in six areas were investigated. Both habitats (root or rhizosphere soil) and soil properties are the key environmental factors affecting AMF communities with C. ambrosioides. Total N, pH, Zn, and Cu levels play significant roles in triggering AMF populations. Overall, C. ambrosioides rhizospheric AMF communities and their influencing factors were illustrated, which contributed a better understanding of the AMF community shaping and related ecological factors.


Root and soil sampling

Samples were collected at six sites in Yunnan, Anhui, and Guangxi Provinces, China. Two sites in Yunnan province were located in Dounan Town (DN, 23°36′ N, 103°41′ E), Yanshan County, and Xincheng Town (XC, 23°38′ N, 102°27′ E), Shiping County. Three sites were located in Anhui Province, including Tongling (TL, 30°49′ N, 117°51′ E), Wuhu (WH, 31°19′ N, 118°26′ E) and Susong (SS, 30°9′ N, 116°10′ E) Counties. Another site was located in Guiping County (GP, 23°31′ N, 110°16′ E), Guanxi province. A single plant, C. ambrosioides, was investigated to reduce the number of biotic factors affecting the AMF community. All samples were collected in August 2015. Three healthy, similarly sized plants were randomly selected and collected at each of the six sites. Root systems were carefully excavated. Using clean tweezers and a brush, the soils bound to the surface of roots were carefully removed. The removed soils were defined as rhizosphere soil samples for DNA extraction. Then the roots were wrapped in tissue paper and stored in sterile Ziploc bags containing silica gel. Topsoil samples (0–15 cm depth) were collected from beneath each plant in at least three directions for chemical property analysis. To remove aboveground plant materials, roots, and stones, the mixed and homogenized soil samples were passed through a 2 mm sieve. After packing in sterile Ziploc bags, these soil samples were transported to the laboratory and stored at 4 °C to analyze physicochemical properties or at −80 °C for DNA extraction.

Soil physicochemical properties

Soils were air-dried at room temperature, ground into a powder, and passed through a 0.15 mm nylon sieve. The soil pH was determined using a glass electrode with soil suspended in 0.01 M CaCl2 (soil: solution ratio, 1:5). The nitrate N (NO3-N), ammonia N (NH4+-N), and soil organic matter (OM) levels in the soil samples were determined according to our previous study59. Available K in soil was extracted with ammonium acetate64, and then determined by flame photometry (AAS novAA 400, Analytik Jena AG, Jena, Germany). Using sodium bicarbonate, the available P in soil was extracted and measured using the molybdenum blue method according to Watanabe and Olsen65. The total N level in soil was determined using the micro-Kjeldahl method66, and measured using a continuous flow analyzer (AA3, Bran + Luebbe, Hamburg, Germany). Soil samples for analyzing the levels of total P, total K, and HM (Mn, Zn, Cd, Cu, and Pb) were dried at 105 °C for 6 h, passed through a 0.15 mm nylon sieve, and then digested in a mixture of HNO3 and HClO4 (4:1, v/v)67. The total concentrations of K, P, Mn, Cu, and Zn were determined by inductively coupled plasma optical emission spectrometry (ICP-OES; 710series, Agilent Technologies, Palo Alto, CA, USA), and inductively coupled plasma mass spectrometry (ICP-MS; NexION 300X, Perkin Elmer, Norwalk, CT, USA) was used to measure the Cd and Pb concentrations. To analyze available HM, 20 mL modified Morgan’s solution (1 M ammonium acetate; pH 4.8) was added to 4 g soil in 100 mL Erlenmeyer flasks. The mixture was shaken on a rotary shaker at 150 rpm for 15 min, and then clear solutions were obtained by filtering through filter paper. The extracts were analyzed in terms of HM concentrations by ICP-MS59. The chemical properties of three soil sample replicates were analyzed independently.

DNA extraction and PCR amplification

Soil DNA was extracted from 0.5 g soil using a FastDNA SPIN Kit following the manufacturer’s instructions (MP Biomedicals, Santa Ana, CA, USA). Root DNA was extracted from 0.1 mg roots using the Plant DNA Extraction Kit (Tiangen Biotech, China) according to the manufacturer’s instructions. In total, 36 samples comprising 18 soil samples and 18 root samples were subjected to DNA extraction. After dissolving in 50 μL TE buffer and quantifying by spectrophotometry, the extracted DNA was stored at −20 °C for further analysis. A conserved AMF-specific primer set, AMV4.5NF (5′-AAGCTCGTAGTTGAATTTCG-3′)/AMDGR (5′-CCCAACTATCCCTATTAATCAT-3′), was used to amplify the partial AMF 18S rRNA gene fragment68. A 6 bp error-correcting barcode was included in the reverse primer to characterize the samples55. The polymerase chain reaction (PCR) was carried out in a 20 μL mixture, including 10 ng DNA template, 0.4 μL FastPfu Polymerase (TransGen Biotech, Beijing, China), 0.8 μL each 5 μM primer, 2 μL 2.5 mM dNTPs, 4 μL 5× FastPfu Buffer, and 0.2 μL bovine serum albumin (BSA; Takara Biotechnology, Dalian, China). The PCR conditions were as follows: 95 °C for 3 min; 28 cycles of denaturation at 95 °C for 30 s, primer annealing at 55 °C for 30 s, and extension at 72 °C for 45 s, followed by a final extension for 10 min at 72 °C.

Illumina MiSeq sequencing

To reduce potential early-round PCR errors, three independent PCR products for each sample were combined to construct a PCR amplicon library. The PCR products were subjected to agarose gel electrophoresis and purified using an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) according to the manufacturer’s protocol. Then the amplified DNA was quantified using QuantiFluor™-ST (Promega, Madison, WI, USA). The purified amplicons were paired-end (PE) sequenced on an Illumina MiSeq platform (Shanghai Biozeron Biotechnology Co., Ltd, Shanghai, China) according to standard protocols. In total, 36 sequencing libraries were constructed and independently sequenced.

Processing of sequencing data

The raw Illumina sequencing data were quality-filtered using the Trimmomatic software69. The reads were truncated at any site that received an average quality score <20 over a 50 bp sliding window, and the truncated reads shorter than 50 bp were discarded. Then, PE reads were assembled according to their overlap sequence with a minimum overlap length of 10 bp, discarding reads that could not be assembled. Sequences that contained more than one ambiguous character or two nucleotide mismatches in the primers, and those with a mismatch ratio within the overlap region of more than 0.2 were removed. The clean sequences were analyzed using the QIIME pipeline70. Chimeric sequences were identified and removed using UCHIME71. Operational taxonomic unit (OTU) grouping was performed using UPARSE at 97% similarity72, and then the representative sequences obtained for each OTU were assigned to taxonomic data using the RDP classifier at a 70% threshold73. The rarefaction curves, indices of ACE and Chao1, and Shannon diversity were analyzed using the Mothur software74. The Bray–Curtis distances were calculated using the QIIME pipeline70, and a principal coordinate analysis (PCoA) was performed using R ( based on the Bray–Curtis similarities. A Venn diagram based on unique and shared OTUs was produced using R to characterize the differences and similarities among the AMF communities. The indicator species analyses were performed to test whether there were specific OTUs associated with the rhizosphere soil habitat or root habitat, and the indicator value index was used to measure the associations75. The analyses were performed using the indcspecies package implemented in R with a permutation test (999 permutations)76. To examine the relationships between AMF communities and soil properties, redundancy analysis (RDA) was applied. This analysis was conducted with CANOCO for Windows77 with 999 permutations of Monte Carlo permutation tests.

Phylogenetic analysis

A neighbor-joining tree containing the type sequences of all OTUs was constructed using MEGA v6.0 with 1000 replicates78. All sequences were aligned, concatenated, and manually adjusted using Geneious Pro v4.8.3 ( The best-fit model for the datasets was selected using jModelTest v279.

Statistical analyses

PCNM vectors were calculated using the ‘pcnm’ function of the ‘vegan’ package with the R language. Prior to the calculation, GPS coordinates were converted to UTM coordinates in kilometers. PCNM vectors were used as explanatory spatial variables for db-RDA80. The weighted Unifrac distance matrix was used as a measure to determine the AMF phylogenetic structure, which was calculated using ‘unifracs’ function of ‘GUniFrac’ package with the R language. db-RDA81, which is known as constrained analysis of principal coordinates, was used to investigate the variations of AMF communities that were attributable to environmental factors. By a stepwise db-RDA, the contributions of environmental factors to the variation of soil AMF communities were summarized using Bray-Curtis dissimilarity and Unifrac distance matrices separately. The forward-selection for the environmental factors in three groups (soil, geographic and climate, and plant variables) was performed independently with the adjR2thresh stopping criterion82, and then the contribution of each of the groups or the combined groups was determined. The amount of variances explained by the individual and combined groups was tested using Monte Carlo permutation tests (999 permutations). The stepwise db-RDA was performed using ‘capscale’ function of ‘vegan’ package with the R language. PerMANOVA for the relationship between AMF community composition and each environmental variable was carried out separately based on Bray-Curtis dissimilarity distance and weighted Unifrac distance matrices using ‘adonis’ function of ‘vegan’ package with the R language. Tukey’s and Kendall test was used for multiple comparisons (P < 0.05), which were performed using SPSS version 19.0 (SPSS Inc., Chicago, IL, USA).

Additional Information

How to cite this article: Xu, X. et al. The influence of environmental factors on communities of arbuscular mycorrhizal fungi associated with Chenopodium ambrosioides revealed by MiSeq sequencing investigation. Sci. Rep. 7, 45134; doi: 10.1038/srep45134 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:


The research was funded by the National Key Research and Development Program of China (2016YFD0800800), and the National Natural Science Foundation of China (31501689 and 41571307).


The authors declare no competing financial interests.

Author Contributions C.C., J.J. and Z.S. designed the experiment; X.X., Z.Z., Z.S. and Y.C. performed the laboratory experiment and data analysis; X.X. and C.C. wrote the manuscript; Z.S. revised the manuscript.


  • Van der Heijden M. G. A. et al. . Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396, 69–72 (1998).
  • Simon L., Bousquet J., Lévesque R. C. & Lalonde M. Origin and diversification of endopmycorrhizal fungi and coincidence with vascular plants. Nature 363, 67–69 (1993).
  • Remy W., Taylor T. N., Hass H. & Kerp H. Four hundred-million-year-old vesicular arbuscular mycorrhizae. Proc. Natl. Acad. Sci. USA 91, 11841–11843 (1994). [PubMed]
  • Redecker D., Kodner R. & Graham L. E. Glomalean fungi from the Ordovician. Science 289, 1920–1921 (2000). [PubMed]
  • Dalpe Y. Mycorrhizae: a potential tool for plant protection but not a panacea. Phytoprotection 86, 53–59 (2005).
  • Smith S. E., Facelli E., Pope S. & Smith F. A. Plant performance in stressful environments: interpreting new and established knowledge of the roles of arbuscular mycorrhizas. Plant Soil 326, 3–20 (2010).
  • Doubková P., Suda J. & Sudová R. The symbiosis with arbuscular mycorrhizal fungi contributes to plant tolerance to serpentine edaphic stress. Soil Biol. Biochem. 44, 56–64 (2012).
  • Lehmann A., Veresoglou S. D., Leifheit E. F. & Rillig M. C. Arbuscular mycorrhizal influence on zinc nutrition in crop plants-A meta-analysis. Soil Biol. Biochem. 69, 123–131 (2014).
  • Oehl F. et al. . Soil type and land use intensity determine the composition of arbuscular mycorrhizal fungal communities. Soil Biol. Biochem. 42, 724–738 (2010).
  • Hazard C. et al. . The role of local environment and geographical distance in determining community composition of arbuscular mycorrhizal fungi at the landscape scale. ISME J. 7, 498–508 (2013). [PMC free article] [PubMed]
  • Jansa J., Erb A., Oberholzer H. R., Smilauer P. & Egli S. Soil and geography are more important determinants of indigenous arbuscular mycorrhizal communities than management practices in Swiss agricultural soils. Mol. Ecol. 23, 2118–2135 (2014). [PubMed]
  • Torrecillas E. et al. . Modularity reveals the tendency of arbuscular mycorrhizal fungi to interact differently with generalist and specialist plant species in gypsum soils. Appl. Environ. Microb. 80, 5457–5760 (2014). [PMC free article] [PubMed]
  • Davison J. et al. . Global assesment of arbuscular mycorrhizal fungus diversity reveals very low endemism. Science 394, 970–973 (2015). [PubMed]
  • Cheng Y., Ishimoto K., Kuriyama Y., Osaki M. & Ezawa T. Ninety-year-, but not single, application of phosphorus fertilizer has a major impact on arbuscular mycorrhizal fungal communities. Plant Soil 365, 397–407 (2013).
  • Gosling P., Mead A., Proctor M., Hammond J. P. & Bending G. D. Contrasting arbuscular mycorrhizal communities colonizing different host plants show a similar response to a soil phosphorus concentration gradient. New Phytol. 198, 546–556 (2013). [PMC free article] [PubMed]
  • Moebius-Clune D. J., Moebius-Clune B. N., van Es H. M. & Pawlowska T. E. Arbuscular mycorrhizal fungi associated with a single agronomic plant host across the landscape: community differentiation along a soil textural gradient. Soil Biol. Biochem. 64, 191–199 (2013).
  • Wei Y. et al. . Molecular diversity of arbuscular mycorrhizal fungi associated with an Mn hyperaccumulator—Phytolacca americana, in Mn mining area. Appl. Soil Ecol. 82, 11–17 (2014).
  • Yang Y. et al. . Community structure of arbuscular mycorrhizal fungi associated with Robinia pseudoacacia in uncontaminated and heavy metal contaminated soils. Soil Biol. Biochem. 86, 146–158 (2015).
  • Gai J. P., Christie P., Feng G. & Li X. L. Twenty years of research on biodiversity and distribution of arbuscular mycorrhizal fungi in China: a review. Mycorrhiza 16, 229–239 (2006). [PubMed]
  • Wetzel K., Silva G., Matczinski U., Oehl F. & Fester T. Superior differentiation of arbuscular mycorrhizal fungal communities from till and no-till plots by morphological spore identification when compared to T-RFLP. Soil Biol. Biochem. 72, 88–96 (2014).
  • de Oliveira Freitas R. et al. . Arbuscular mycorrhizal fungal communities along a pedohydrological gradient in a Central Amazonian terra firme forest. Mycorrhiza 24, 21–32 (2014). [PubMed]
  • Clapp J. P., Young J. P. W., Merryweather J. W. & Fitter A. H. Diversity of fungal symbionts in arbuscular mycorrhizas from a natural community. New Phytol. 130, 259–265 (1995).
  • Sanders I. R. Plant and arbuscular mycorrhizal fungal diversity-are we looking at the relevant levels of diversity and are we using the right techniques? New Phytol. 164, 415–418 (2004).
  • Hempel S., Renker C. & Buscot F. Differences in the species composition of arbuscular mycorrhizal fungi in spores, root and soil communities in grassland ecosystem. Environ. Microbiol. 9, 1930–1938 (2007). [PubMed]
  • Chen Y. et al. . Six-year fertilization modifies the biodiversity of arbuscular mycorrhizal fungi in a temperate steppe in Inner Mongolia. Soil Biol. Biochem. 69, 371–381 (2014).
  • Krishnamoorthy R., Kim K., Kim C. & Sa T. Changes of arbuscular mycorrhizal traits and community structure with respect to soil salinity in a coastal reclamation land. Soil Biol. Biochem. 72, 1–10 (2014).
  • Lekberg Y. et al. . 454-Sequencing reveals stochastic local reassembly and high disturbance tolerance within arbuscular mycorrhizal fungal communities. J. Ecol. 100, 1–10 (2012).
  • Öpik M. et al. . Global sampling of plant roots expands the described molecular diversity of arbuscular mycorrhizal fungi. Mycorrhiza 23, 411–430 (2013). [PubMed]
  • Barea J. M. & Jeffries P. Arbuscular mycorrhizas in sustainable soil plant systems In Mycorrhiza Structure, Function, Molecular Biology and Biotechnology (eds Hock B. & Varma A.) 521–560 (Springer, 1995).
  • Brundrett M. C. Coevolution of roots and mycorrhizas of land plants. New Phytol. 154, 275–304 (2002).
  • Miransari M., Bahrami H. A., Rejali F. & Malakouti M. J. Using arbuscular mycorrhiza to reduce the stressful effects of soil compaction on wheat (Triticum aestivum L.) growth. Soil Biol. Biochem. 40, 1197–1206 (2008).
  • Punaminiya P. et al. . Symbiotic role of Glomus mosseae in phytoextraction of lead in vetiver grass [Chrysopogon zizanioides (L.)]. J. Hazard. Mater. 177, 465–474 (2010). [PubMed]
  • Fellbaum C. R. et al. . Fungal nutrient allocation in common mycorrhizal networks is regulated by the carbon source strength of individual host plants. New Phytol. 203, 646–656 (2014). [PubMed]
  • del Mar Alguacil M., Torres M. P., Montesinos-Navarro A. & Roldán A. Soil characteristics driving arbuscular mycorrhizal fungal communities in semiarid Mediterranean soils. Appl. Environ. Microbiol. 82, 3348–3356 (2016). [PMC free article] [PubMed]
  • Johnson N. C., Rowland D. L., Corkidi L., Egerton-Warburton L. M. & Allen E. B. Nitrogen enrichment alters mycorrhizal allocation at five mesic to semiarid grasslands. Ecology 84, 1895–190 (2003).
  • Davison J., Öpik M., Daniell T. J., Moora M. & Zobel M. Arbuscular mycorrhizal fungal communities in plant roots are not random assemblages. FEMS Microb. Ecol. 78, 103–115 (2011). [PubMed]
  • Xu T. et al. . Plant community, geographic distance and abiotic factors play different roles in predicting AMF biogeography at the regional scale in northern China. Environ. Microbiol. Rep. 8, 1048–1057 (2016). [PubMed]
  • Zobel M. & Öpik M. Plant and arbuscular mycorrhizal fungal (AMF) communities – which drives which? J. Veg. Sci. 25, 1133–1140 (2014).
  • Kivlin S. N., Hawkes C. V. & Treseder K. K. Global diversity and distribution of arbuscular mycorrhizal fungi. Soil Biol. Biochem. 43, 2294–2303 (2011).
  • Bever J. D., Richardson S. C., Lawrence B. M., Holmes J. & Watson M. Preferential allocation to beneficial symbiont with spatial structure maintains mycorrhizal mutualism. Ecol. Lett. 12, 13–21 (2009). [PubMed]
  • Saks Ü. et al. . Root-colonizing and soil-borne communities of arbuscular mycorrhizal fungi in a temperate forest understory. Botany 92, 277–285 (2014).
  • Varela-Cervero S. et al. . The composition of arbuscular mycorrhizal fungal communities differs among the roots, spores and extraradical mycelia associated with five Mediterranean plant species. Environ. Microbiol. 17, 2882–2895 (2015). [PubMed]
  • Yang R., Zan S., Tang J., Chen X. & Zhang Q. Variation in community structure of arbuscular mycorrhizal fungi associated with a Cu tolerant plant-Elsholtzia splendens. Appl. Soil Ecol. 44, 191–197 (2010).
  • Long L. et al. . Molecular community analysis of arbuscular mycorrhizal fungi associated with five selected plant species from heavy metal polluted soils. Eur. J. Soil Biol. 46, 288–294 (2010).
  • Zarei M. et al. . Molecular diversity of arbuscular mycorrhizal fungi in relation to soil chemical properties and heavy metal contamination. Environ. Pollut. 158, 2757–2765 (2010). [PubMed]
  • Birgander J., Rousk J. & Olsson P. A. Comparison of fertility and seasonal effects on grassland microbial communities. Soil Biol. Biochem. 76, 80–89 (2014).
  • Dai M. et al. . 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 (2014).
  • Öpik M. et al. . Divergent arbuscular mycorrhizal fungal communities colonize roots of Pulsatilla spp. in boreal Scots pine forest and grassland soils. New Phytol. 160, 581–593 (2003).
  • Helgason T., Merryweather J. W., Young J. P. W. & Fitter A. H. Specificity and resilience in the arbuscular mycorrhizal fungi of a natural woodland community. J. Ecol. 95, 623–630 (2007).
  • Oehl F. et al. . Distinct sporulation dynamics of arbuscular mycorrhizal fungal communities from different agroecosystems in long-term microcosms. Agr. Ecosyst. Environ. 134, 257- 268 (2009).
  • Denison R. F. & Kiers E. T. Life histories of symbiotic rhizobia and mycorrhizal fungi. Curr. Biol. 21, 775–785 (2011). [PubMed]
  • López-García Á., Azcón-Aguilar C. & Barea J. M. The interactions between plant life form and fungal traits of arbuscular mycorrhizal fungi determine the symbiotic community. Oecologia 176, 1075–1086 (2014). [PubMed]
  • Wei Y. et al. . Diversity of Arbuscular Mycorrhizal Fungi Associated with a Sb Accumulator Plant, Ramie (Boehmeria nivea), in an Active Sb Mining. J. Microbiol. Biotechno. 25, 1205–1215 (2015). [PubMed]
  • Bainard L. D. et al. . 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 (2015).
  • Da Silva I. R. et al. . Diversity of arbuscular mycorrhizal fungi along an environmental gradient in the Brazilian semiarid. Appl. Soil Ecol. 84, 166–175 (2014).
  • Ban Y., Zhouying X., Zhang H., Chen H. & Tang M. Soil chemestry, traslocation of heavy metals, and mycorrhizal fungi associated with six plant species growing on lead-zinc mine tailings. Ann. Microbiol. 65, 503–515 (2015).
  • Cavagnaro T. R. Arbuscular mycorrhizas and their role in plant zinc nutrition in Mycorrhizal fungi: use in sustainable agriculture and land restoration (eds Solaiman Z. M., Abbott L. K. & Varma A.) 189–200 (Springer, 2014).
  • Lehmann A. & Rillig M. C. Arbuscular mycorrhizal contribution to copper, manganese and ironnutrient concentrations in crops - a meta-analysis. Soil Biol. Biochem. 81, 147–158 (2015).
  • Chen C. et al. . Microbial communities of an arable soil treated for 8 years with organic and inorganic fertilizers. Biol. Fertil. Soils 52, 455–467 (2016).
  • Avio L. et al. . Impact of nitrogen fertilization and soil tillage on arbuscular mycorrhizal fungal communities in a Mediterranean agroecosystem. Soil Biol. Biochem. 67, 285–294 (2013).
  • Van Diepen L. T., Lilleskov E. A. & Pregitzer K. S. Simulated nitrogen deposition affects community structure of arbuscular mycorrhizal fungi in northern hardwood forests. Mol. Ecol. 20, 799–811 (2011). [PubMed]
  • Treseder K. K. & Allen M. F. Direct nitrogen and phosphorus limitation of arbuscular mycorrhizal fungi: a model and field test. New Phytol. 155, 507–515 (2002).
  • Porras-Alfaro A., Herrera J., Natvig D. O. & Sinsabaugh R. L. Effect of long-term nitrogen fertilization on mycorrhizal fungi associated with a dominant grass in a semiarid grassland. Plant Soil 296, 65–75 (2007).
  • Merwin H. D. & Peech M. Exchangeability of soil potassium in sand, silt and clay fractions as influenced by the nature of the complementary cation. Soil Sci. Soc. Am. J. 15, 125–128 (1950).
  • Watanabe F. S. & Olsen S. R. Test of an ascorbic acid method for determining phosphorus in water and NaHCO3 extracts from soil. Soil Sci. Soc. Am. J. 29, 677–678 (1965).
  • Yeomans J. C. & Bremner J. M. A rapid and precise method for routine determination of organic carbon in soil. Commun. Soil Sci. Plant 19, 1467–1476 (1988).
  • Bansal S. & Kapoor K. K. Vermicomposting of crop residues and cattle dung with Eisenia foetida. Bioresour. Technol. 73, 95–98 (2000).
  • Lumini E., Orgiazzi A., Borriello R., Bonfante P. & Bianciotto V. Disclosing arbuscular mycorrhizal fungal biodiversity in soil through a land-use gradient using a pyrosequencing approach. Environ. Microbiol. 12, 2165–2179 (2010). [PubMed]
  • Bolger A. M., Lohse M. & Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014). [PMC free article] [PubMed]
  • Caporaso J. G. et al. . QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336 (2010). [PMC free article] [PubMed]
  • Edgar R. C., Haas B. J., Clemente J. C., Quince C. & Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200 (2011). [PMC free article] [PubMed]
  • Edgar R. C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998 (2013). [PubMed]
  • Caporaso J. G. et al. . Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA 1, 4516–4522 (2011). [PubMed]
  • Schloss P. D. et al. . Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microb. 75, 7537–7541 (2009). [PMC free article] [PubMed]
  • Dufrêne M. & Legendre P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol. Monograph 67, 345–366 (1997).
  • De Cáceres M. & Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology 90, 3566–3574 (2009). [PubMed]
  • Etten E. V. Multivariate analysis of ecological data using CANOCO. Austral. Eco. 30, 486–487 (2005).
  • Tamura K., Stecher G., Peterson D., Filipski A. & Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729 (2013). [PMC free article] [PubMed]
  • Darriba D., Taboada G. L., Doallo R. & Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat. Methods 9, 772–772 (2012). [PMC free article] [PubMed]
  • Dray S., Legendre P. & Peres-Neto P. R. Spatial modelling: a comprehensive framework for principal coordinate analysis of neighbour matrices (PCNM). Ecol. Model. 196, 483–493 (2006).
  • Legendre P. & Anderson M. J. Distance-based redundancy analysis: Testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 69, 1–24 (1999).
  • Blanchet F. G., Legendre P. & Borcard D. Forward selection of explanatory variables. Ecology 89, 2623–2632 (2008). [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group