In this study, we evaluated the EF colonization of A. villosum roots under different planting locations and growth ages for the first time. High-throughput sequencing technology was employed to study the abundance, composition and distribution of root EF species.
AMF and DSEs are the two main groups of EFs and play critical roles in the growth and development of host plants. Numerous studies have shown that AMF can promote plant growth and development under different habitats and nutritional conditions, particularly by enriching the phosphorus pool in the root system [32]. It has been confirmed that DSEs also have a similar function, but more evidence is needed to further elucidate the mechanism. Studies have shown that large differences exist in the colonization rates of AMF and DSEs among plants, which are affected by climate, moisture, latitude and stress factors [33–36]. For the triennial samples, AMF colonization rates of A. villosum ranged from 13.03 to 39.48% with a mean value of 27.02%, which was closer to the Longjing tea root AMF colonization rate of 27.54% [37]. Jin et al. reported a significant 70% increase in the AMF colonization rate of Solidago canadensis under dry habitats over time [38], which was in line with our results of increased AMF colonization rates of samples of BNP, WSP and YCP. However, the AMF colonization rates of the samples HHP decreased with the age, which may be related to the competition of ecological sites with other microorganisms and the use of excessive fertilizers and pesticides.
Compared to the DSE colonization rates of 21.25% of Artemisia deversa [39], annual samples of A. villosum from five sampling plots had lower DSE colonization rates, with a range of 0 ~ 9.38% and a mean value of 3.11%. However, the mean DSE colonization rate increased with increasing growth age, reaching 6.15%. In addition, the colonization rates of AMF and DSE in sample PEP-1 were relatively higher than those in other samples, supporting the hypothesis that AMF and DSE are mutually reinforcing [32]. This result was confirmed by the studies of Della Mónica et al. and Biotic et al., which reported positive correlations between AMF and DSE colonization in Trifolium repens and Achnatherum lettermanii, and both were closely related to the utilization and absorption of phosphorus [32, 40].
In this study, the difference in fungal taxa was nonsignificant among A. villosum samples from five sampling plots, and the dominant taxa at the phylum level were Ascomycota, Basidiomycota and Glomeromycota, which was the same as the composition of root fungi of Casuarina equisetifolia [35]. At the order level, Hypocreales, Chaetothyriales and Helotiales were observed in all samples, as well as the families Cordycipitaceae and Cladosporiaceae, which were identified in other plants with biocontrol potential, e.g., Fragaria ananassa [41] and Glycine max [42]. Hypocreales are widely marketed worldwide as fungicides and Cladosporiaceae fungi have potential value in the control of white rust [43, 44]. The antimicrobial fungi Phoma, saprophytic fungi Acremonium and Myrothecium and the biocontrol fungus Trichoderma were detected as the main genera in our study, which can contribute to the improvement of rhizosphere soil microstructure. The genus Beauveria was one of the dominant microbes in A. villosum roots, which not only does not cause plant diseases but is also able to prevent insect pests and exchange nitrogen sources in the host plant [45]. Therefor, the roles of Beauveria need to be viewed dialectically. It is worth mentioning that AMF and DSE were the main components of the A. villosum root system, directly supporting the colonization rates, in which Glomus of AMF and Exophiala, Cladosporium and Cladophialophora of DSE were identified as the main genera. According to previous studies, AMF and DSE were proven to be the main fungal community in Glycyrrhiza uralensis, Ferula sinkiangensis and Tinospora cordifolia, e.g., Glomus, Diversispora and Cladosporium [46–48]. The colonization of these diverse mycorrhizal fungi may form a complex network structure in or around the root system, which plays an important role in maintaining the growth of plants and resisting diseases.
The alpha diversity results showed that there was no significant difference in richness and diversity between the annual and triennial samples, while some differences existed due to spatial differences in the samples. Similar to the results reported by Huang et al., the diversity indices of C. equisetifolia endophytic fungal communities did not differ significantly between forest ages [35]. This phenomenon may be attributed to the fact that plant root endophytic fungal diversity is largely influenced by local biotic and abiotic influences [49], e.g., soil properties as well as human disturbances both cause changes in root microbial diversity [50, 51]. It is notable that fungal diversity was consistently highest in the HHP sample plot among the five plots, which may be the result of multiple factors. However, the abundance of AMF decreased in triennial samples of HHP, suggesting that it may be related to the change in soil nutrient composition over time and may require the involvement of other microorganisms for adjustment.
The community composition analysis results showed that the fungal community composition of A. villosum did not differ significantly by age, and the spatial location was the main cause of community differences. Random forest analysis revealed that the fungal taxa were responsible for the differences in community composition to some extent. In line with the report of Zhong et al[52], Nigrospora was taken as a biomarker or key taxon in this study. It was found that the AMF strain Diversispora was involved in the growth promotion of Chrysanthemum morifolium under salt stress [53], and the strain was significantly enriched in annual samples, implying that it may play a key role in this phase. In contrast, Aspergillus and Cladosporium (DSE) were enriched in triennial samples and showed significant differences from annual A. villosum samples, of which these strains were reported to be valuable for heavy metal resistance and may be related to their involvement in maintaining and promoting plant growth [54, 55]. Colonization surveys and molecular diversity analyses responded to the importance that AMF and DSEs were involved in the evolution of plant rhizosphere fungal communities and that there are specific community structure changes due to spatial and temporal differences.