- 1Institute of Microbiology, Jiangxi Academy of Sciences, Nanchang, China
- 2Nanchang Institute of Technology, Nanchang, China
Soil amendments, including various types of fertilizers, are often used to control the uptake of heavy metals such as cadmium in cropping fields. The influence of these amendments on other members of the agroecosystem, such as arbuscular mycorrhizal fungi (AMF), remains less well investigated. Here, we established an experiment with the application of woody peat organic fertilizer and phosphate rock powder to examine its effects on AMF communities in two cadmium-contaminated vegetable crop fields (cucumber and pepper). We found that the application of phosphate rock powder enhanced soil phosphorus content, while the application of woody peat organic fertilizer enhanced soil nitrogen content, but neither influenced AMF abundance. We also found little influence of either amendment on measures of AMF diversity, except in one case where the Shannon index of diversity was lower in pepper fields amended with phosphate rock powder. We did, however, find significant shifts in the community composition and relative abundances of AMF taxa in the two vegetable fields, primarily as a result of shifts in the soil pH and nitrogen content.
1 Introduction
As a result of ongoing anthropogenic activities, agricultural soils are being continually degraded (Durdu et al., 2023). Notably, heavy metals are accumulating in agricultural soils globally because of increased uses of pesticides and fertilizers, as well as other industrial actions (Rashid et al., 2023). This heavy metal contamination can negatively influence crop production through a number of pathways, and can also influence human health. For example, cadmium (Cd) is a non-essential element for plant growth, but owing to its high chemical activity in soil, can be easily absorbed by plants and accumulate within plant tissues, ultimately posing risks to human health through consumption (McLaughlin et al., 2021). Levels of Cd have increased in agricultural soils in recent decades (Six and Smolders, 2014; Shi et al., 2019), and now reducing the Cd content and activity in soils, as well as decreasing its absorption and accumulation by plants, has become a priority for safe production of edible agricultural products (Mubeen et al., 2023).
For lightly and moderately polluted soils, a number of soil amendments, including organic material and phosphate, have been commonly applied to reduce the bioavailability of heavy metals in the soil through the mechanisms of adsorption, complexation and precipitation (Gao J. et al., 2023; Wang F. et al., 2022; Cui et al., 2023). While there has been work to understand the impact of these activities on farmland soil, such as the positive effects of reducing plant heavy metal content and crop yield (Feng and Cheng, 2023), less attention has been paid to the influence of these amendments on environmental and ecological risks, such as the impact of these measures on soil microbial community structure and their functions (Holland et al., 2018). For example, how these amendments might influence crop interactions with soil microbes, such as arbuscular mycorrhizal fungi (AMF) that form symbiotic relationships with over 70% of terrestrial plants (Brundrett and Tedersoo, 2018), obtaining carbon from plants in exchange for helping plants to absorb mineral nutrients and water (Smith and Smith, 2011).
In agriculture, AMF can promote crop yields while reducing the application of chemical fertilizers, controlling disease occurrence, and enhancing the host plant’s adaptability to drought, salinity and heavy metal pollution (Hao et al., 2022). AMF can also reduce the uptake and transport of Cd in crops through various mechanisms, including regulating plant and root growth, binding heavy metals with root and hyphal exudates, increasing antioxidant enzyme activities, interacting with soil beneficial microorganisms, and immobilizing heavy metals via extraradical hyphae as well as their improvement on soil aggregation (Wang, 2017; Wang H. et al., 2022; Moreno Jiménez et al., 2024). Studies on the application of AMF inoculant, combined with other remediation strategies, to mitigate Cd contamination in soils has shown promising results for various crops (Cakmak et al., 2023; Zulfiqar et al., 2023; Liu et al., 2024).
Given that AMF are widespread soil fungi in agricultural ecosystems, native AMF communities may also significantly influence the Cd remediation efficacy of soil amendments. Therefore, evaluating the interactions between soil amendments and native AMF communities is essential when developing strategies to mitigate Cd contamination. Findings from studies examining the effects of soil environmental changes on AMF under conditions that do not involve or explicitly consider Cd contamination remain inconclusive and highly context-dependent (van der Heyde et al., 2017; Mohamed et al., 2023). For example, some studies report that AMF abundance are enhanced by organic amendments or reduced by excessive phosphate, while others find more nuanced or contrasting outcomes (Treseder, 2004; Park et al., 2024). The diversity of outcomes highlights the complexity of AMF interactions with soil environments and the necessity for more targeted research that considers the specific conditions under which AMF are studied. This includes exploring how soil amendments can influence the structure and function of AMF communities in the context of Cd contamination.
Here, we used an experimental approach to explicitly investigate the effects of soil amendments on the AMF community in Cd contaminated vegetable fields. We examined changes in both soil chemistry and the AMF community, and their associations, to understand the primary factors influencing AMF community. We used two types of vegetable crops, pepper (Capsicum annuum) and cucumber (Cucumis sativus), and two types of soil amendments (phosphate rock powder and woody peat organic fertilizer). In all, our research will provide insights into the responses, protection and potential application of AMF communities in farmland soils.
2 Materials and methods
2.1 Site
Our study took place in Jiangxi Province, China (30°31′N, 115°56′E), characterized by a subtropical continental monsoon climate. The annual average temperature is 17.5°C and the average annual precipitation is 1,600 mm. The soil type in the area is typical red soil, with Cd levels ranging from 0.61 to 1.23 mg kg−1. According to the national risk control standard for soil contamination of agricultural land, the risk screening values for Cd contamination in vegetable fields are 0.3 mg kg−1 for soils with a pH below 7.5 and 0.6 mg kg−1 for soils with a pH above 7.5 (Chinese National Standardization Administration, 2018). Therefore this site is considered to be light to moderately Cd contaminated. The site of our field experiment is cultivated with various vegetables, including peppers, cucumbers, Chinese cabbage, cauliflower and tomatoes.
2.2 Experimental design
In June 2022, we established 13 plots of 1 m × 20 m in each of two adjacent land parcels (26 plots in total), one where peppers were cultivated and the other where cucumbers were cultivated. At each site, we established a completely randomized design with three treatments: (1) the addition of phosphate rock powder (P, with five replicates), with the particle size of 75 μm, consisting of 13.9% P2O5 and 0.033 mg kg−1 of Cd; (2) the addition of woody peat organic fertilizer (O, with five replicates), containing ≥45% organic matter, ≥10% humic acid, and ≥ 5% N + P2O5 + K2O, with a Cd content of 0.068 mg kg−1; (3) a control (C, no amendment was added, with three replicates). Both soil amendments were applied at a rate of 0.3 kg·m−2, following the application rates commonly used in local Cd-contaminated soil remediation practices.
After crops were harvested in October, we selected five random locations within each plot to collect soil cores. These locations were systematically distributed along the longer edges of each plot to ensure representative sampling. We determined a sampling depth of 20 cm based on the average tillage depth commonly employed in this vegetable fields. To collect an adequate volume of root segments for assessing mycorrhizal colonization, we used a large-diameter soil corer (10 cm). We thoroughly mixed the collected soil to form a single sample, with a total volume of ~3 kg. We picked at least 30 small root segments from each soil sample and preserved in a 50% v/v alcohol solution for subsequent analysis of mycorrhizal colonization. Then, we removed two 10 g soil samples which were frozen at −20°C. We used one of these samples for DNA sequencing and the other for the analysis of hyphal density. Finally, we used the remaining soil for chemical analysis after drying at room temperature.
By analyzing the Cd content in cucumber and pepper fruits (unpublished data), we observed that Cd concentrations in fruits from amended fields were consistently below the national food safety threshold of 0.05 mg kg−1 fresh weight. Additionally, while the application of soil amendments resulted in a Cd reduction of 7 to 18% compared to the control fields, this decrease, though notable, did not reach statistical significance (p > 0.05). These results highlight the potential, albeit limited, efficacy of soil amendments in mitigating Cd accumulation in edible plant tissues.
2.3 Soil chemistry
We analyzed soil organic matter using the potassium dichromate titration method; total nitrogen and available (alkali-hydrolyzable) nitrogen using a Kjeldahl nitrogen analyzer; total phosphorus and Cd by inductively coupled plasma optical emission spectrometry; available phosphorus extracted with hydrochloric acid-ammonium fluoride and determined by the molybdenum antimony colorimetric method; and pH by the potentiometric method with a water-to-soil quality ratio of 2.5:1.
2.4 AMF abundance
We estimated AMF abundance including both mycorrhizal colonization rate and hyphal length density using the methods in Brundrett et al. (1994). We counted AMF intraradical structures using the intersecting lines method (McGonigle et al., 1990) with 200 fields of view examined per sample and calculated mycorrhizal colonization rate (%) as the ratio of the number of fields of view containing AMF structures to the total number of fields. We counted the extraradical hyphal length using the grid intersection method (Tennant, 1975) with 100 fields of view observed per filter membrane. We calculated the total hyphal length on each filter using the formula:
where h represent the hyphal length (mm), Σc is the total number of intersections between hyphae and grid lines, Σs was the total number of fields observed (100 in this study), a was the side length of each small square in a 10-part grid (0.05 mm in this study), and g was the side length of the grid (0.5 mm in this study). Finally, the hyphal length density was converted to the density per unit dry weight of soil based on the soil moisture content.
2.5 DNA extraction and sequencing
We used high-throughput sequencing to analyze the diversity and compositional shifts of AMF communities. To do so, we extracted total DNA from 300 mg of soil samples, from which we assessed the concentration and purity. We used a nested PCR to amplify the 18S rDNA small subunit (SSU) regions using the primer pairs AML1/AML2 (Lee et al., 2008) for the first amplification and AMV4.5NF/AMDGR for the second amplification (Sato et al., 2005). We recovered the PCR products, purified and eluted them, and used electrophoresis for detection, which we then quantified. DNA extraction and sequencing were performed by Shanghai Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China), utilizing the Illumina MiSeq platform (Illumina, United States) for sequencing, with a minimum sequencing depth of 45,000 reads per sample.
2.6 Sequencing data processing
We quality-controlled the raw sequencing reads using Trimmomatic (Bolger et al., 2014) and assembled using FLASH software (Magoč and Salzberg, 2011). We removed chimeras using UCHIME software (Edgar et al., 2011). We clustered sequences subjected to quality control into Operational Taxonomic Units (OTUs) at a 95% similarity threshold using UPARSE (Edgar, 2013), and classified and annotated species using the RDP classifier,1 followed by comparison and identification with the MaarjAM database (Öpik et al., 2010).
2.7 Statistical analysis
We used one-way ANOVA to compare differences in soil chemical composition, AMF colonization rate and hyphal length density, and the OTU numbers and Shannon diversity index of AMF communities among different treatments for each of the two vegetable fields. To examine the relationship between AMF abundance and soil factors, we used Pearson correlation analysis. To compare the AMF community composition among treatments, we ordinated composition using non-metric multidimensional scaling (NMDS) based on Bray–Curtis distances using the vegan R package (Oksanen et al., 2017). We then used the ‘adonis’ function to conduct a permutational multivariate analysis of variance (PERMANOVA) to compare AMF community structure among the treatments. We employed a Variance Inflation Factor (VIF) analysis to test for collinearity among variables when analyzing the factors influencing community composition using the ‘vif’ function in the car package (Fox and Weisberg, 2019). Next, we used the ‘envfit’ function to fit the selected soil factors to the NMDS ordination. We used the Kruskal-Wallis H test to analyze the differences in relative abundance of the top 20 most abundant OTUs across all treatments in each vegetable field in response to the application of soil amendments compared to the control. Differences among groups were assessed using Scheffe post hoc tests. AMF families typically exhibit conserved functional traits and responses to environmental changes (Chagnon et al., 2013), which makes family-level analysis relevant for understanding broad ecological patterns. By integrating both family and OTU-level analyses, we aimed to achieve a more nuanced and comprehensive perspective on how soil amendments impact the AMF community. We used the pheatmap package (Kolde, 2019) to calculate the Pearson correlation coefficients between soil chemistry and the AMF families, as well as differences among the top 20 most abundant OTUs. We conducted all statistical analyses in R 4.1.0 (R Core Team, 2022).
3 Results
3.1 Soil chemistry
We found that the addition of woody peat organic fertilizer and phosphate rock powder influenced several soil chemical variables compared to control treatments in both vegetable fields (Table 1). Specifically, both total and available nitrogen content increased in both vegetable fields with the application of organic fertilizer, while available phosphorus content increased in the pepper field with the application of phosphate rock powder. Soil pH was inherently lower in cucumber fields than that in pepper fields. The addition of organic fertilizer led to lower pH in pepper, but not cucumber fields. Finally, we found no differences in other soil chemical variables among the three treatments or between the two types of vegetable fields.
3.2 AMF abundance
There was no difference in mycorrhizal colonization rates in both pepper and cucumber following the addition of soil amendments compared to the control treatment (Figure 1A) (F2,10 = 0.9, p = 0.44; F2,10 = 1.5, p = 0.28). Likewise, we found no difference in the hyphal length densities between the treatments (Figure 1B) (F2,10 = 2.9, p = 0.10; F2,10 = 0.4, p = 0.70).
Figure 1. Effects of applying nutrient amendments on AMF colonization rate (A) and hyphal length density (B). Bars represent means ± SE, with different letters indicating significant differences (p < 0.05). C, control; P, phosphate rock powder; O, organic fertilizer.
We found a significant negative correlation between mycorrhizal colonization rate and soil total nitrogen content (r = −0.48, p = 0.014), available nitrogen (r = −0.44, p = 0.024), as well as the soil nitrogen-to-phosphorus ratio (r = −0.50, p = 0.010). No significant correlation was found between the AMF abundance and the other soil variables (Table 2).
Table 2. Correlation between mycorrhizal colonization rate and hyphal length density and soil chemical variables across the three treatments and two vegetable fields.
3.3 AMF community
We found a total of 1,567,232 optimized amplicon sequences from the 26 plots. After subsampling to a minimum sample sequence of 45,829, we assigned a total of 1,191,554 sequences to 193 taxa belonging to 7 families and 10 genera within the Glomeromycota phylum. Of these, 125 OTUs were identified in pepper fields and 149 OTUs were identified in cucumber fields. Across both vegetable crops, members of the Glomeraceae were the most dominant (76.3% of all sequences), followed by Paraglomeraceae (14.4%), Archaeosporaceae (7.4%) and Gigasporaceae (1.5%). These four families together comprised 99.5% of all sequences, while members of the Claroideoglomeraceae, Acaulosporaceae, and Ambisporaceae comprised less than 0.5%.
The effects of soil amendments on the α-diversity of AMF communities varied with the type of amendments and vegetable species (Figure 2). We found that neither soil amendment influenced the number of AMF taxa in the soils of both vegetable fields (Figure 2A). For the Shannon index, we found a significant decrease in the pepper field after the addition of phosphate rock powder, but no such effect in the cucumber field, or with the application of organic fertilizer in either field (Figure 2B).
Figure 2. Effects of applying nutrient amendments on AMF OTU richness (A) and Shannon diversity index (B). Bars represent means ± SE, with different letters indicating significant differences (p < 0.05). C, control; P, phosphate rock powder; O, organic fertilizer.
At the family level, we found that the application of soil amendments led to a decrease in the relative abundance of the most dominant Glomeraceae in both vegetable fields (Figure 3). In the pepper field, Archaeosporaceae, which were absent in the control treatment increased to more than 35% in response to the addition of both amendments. Likewise, the relative abundance of Gigasporaceae, which was rare in the control treatment increased with the addition of organize fertilizer in both vegetable fields. In the cucumber field, the relative abundance of Paraglomeraceae increased substantially from the control treatment to the addition of both amendments.
Figure 3. Effects of applying amendments on the relative abundance of AMF families. C, control treatment; P, phosphate rock powder; O, organic fertilizer.
Among the top 20 most abundant OTUs in the pepper field, four decreased in relative abundance after the addition of phosphate rock powder, while two decreased and one increased after the addition of organic fertilizer (Figure 4). In the cucumber field, one taxon decreased in response to phosphate rock powder addition, and one taxon increased and one decreased in response to organic fertilizer (Figure 4).
Figure 4. Differences in relative abundances of the dominant OTUs in each of the two vegetable fields applying soil amendments compared to the control. P, phosphate rock powder; O, organic fertilizer. *p < 0.05, **p < 0.01.
We found significant differences in the composition of the AMF community among the different vegetable fields and soil amendment treatments (p < 0.01). These differences were predominantly associated with variation in soil nitrogen and the nitrogen-phosphorus ratio (Figure 5) (p < 0.001). Other soil factors significantly associated with changes in AMF community structure included soil Cd content (p < 0.01) and pH (p < 0.05). In contrast, variations in soil phosphorus had a relatively minor contribution to these differences.
Figure 5. NMDS ordination showing changes in AMF community composition. Points represent each replicate of each treatment, while ellipses describe 95% confidence intervals for a given treatment. Arrows indicate the direction and degree of correlations between soil chemical characteristics and the NMD axes. C, control treatment; P, phosphate rock powder; O, organic fertilizer. TN, total nitrogen; AN, available nitrogen; TP, total phosphorus; AP, available phosphorus; N: P, total nitrogen-to-phosphorus ratio; Cd, total cadmium.
At the family level, the difference in relative abundance of the most abundant taxa, Glomeraceae, was significantly correlated with soil total nitrogen. The relative abundance of Gigasporaceae was significantly affected by several factors including Cd, total nitrogen, available nitrogen, available phosphorus content, and the nitrogen-phosphorus ratio. The Archaeosporaceae was influenced by soil pH and total phosphorus content (Figure 6). At the OTU level, of the top 20 most abundant OTUs, 12 showed significant correlations with soil factors (Figure 6). While most (11/12) of these OTUs was associated with only one soil factor, one was associated with four soil factors. Of the soil factors, available nitrogen was most strongly associated with variation in the most OTUs (Figure 6).
Figure 6. Correlation heat map between AMF families and dominant OTUs and soil chemistry. The colors indicate the correlation coefficient (r). TN, total nitrogen; AN, available nitrogen; TP, total phosphorus; AP, available phosphorus; N:P, total nitrogen-to-phosphorus ratio; Cd, total cadmium. *p < 0.05, **p < 0.01.
4 Discussion
Overall, it was not surprising that the application of organic fertilizer and phosphate rock powder as soil amendments led to significant increases in soil nitrogen and phosphorus levels (Table 1). The increase in soil nutrients could allow plants to rely more on their own root system for nutrient acquisition, reducing their dependence on AMF (Treseder, 2004; Singh et al., 2022) because the carbon costs of AMF associations may outweigh the nutritional benefits they provide (Johnson, 2010). However, despite the significant negative correlation between soil nitrogen and AMF abundance, we found no influence of the application of soil amendments on mycorrhizal colonization rates and hyphal length density (Figure 1). This is likely because the ambient levels of nitrogen and phosphorus are high in the soil, suggesting that nutrient acquisition may not be a primary role of AMF in symbiosis with plants. Indeed, the positive impacts of AMF on soil aggregation, water uptake and transport, and disease resistance can also play an important role in the plant-AMF symbiosis, but be less influenced by nutrient availability (Delavaux et al., 2017). Moreover, the relatively high ambient total phosphorus and organic carbon in the soil may have mitigated their responses to the common amounts of amendment addition (7% of total phosphorus and 2% of total organic carbon). This could partly explain the lack of significant changes in AMF abundance.
Despite the lack of response of AMF abundance to soil amendments, we did find lower Shannon diversity (but not richness) of AMF communities after applying phosphate rock powder. The adverse effects of nutrient additions on AMF diversity are consistent with results in many other ecosystems, including grasslands (Santos et al., 2006), alpine meadows (Liu et al., 2012), tropical rainforests (Camenzind et al., 2014), and other agroecosystems (Lin et al., 2012; Gosling et al., 2013). This is also consistent with the NMDS ordination results which indicated that applying phosphate rock powder led to a phylogenetically more dispersed AMF community (Figure 5). This community is more likely to be structured by competition, permitting taxa with various traits to respond to nutrient variations in different ways (Webb et al., 2002). However, the lack of response of AMF diversity to the nitrogen inputs from organic fertilizer may have been because the differential responses of AMF taxa to nitrogen and phosphorus addition (Camenzind et al., 2014).
Even when there was no influence of amendments on abundance or richness of AMF, we did consistently find shifts in the AMF community structure among treatments (Figure 5). Correspondingly, at both the family and OTU levels, there were notable changes in the relative abundance of dominant groups after the application of soil amendments, leading to increases of some taxa and decreases of others (Figures 3, 4). Different AMF taxa vary in their allocation to the intraradical structures that obtain carbon from plant roots and extra-radical hyphae that are responsible for obtaining nutrients from the soil (Maherali and Klironomos, 2012). As a result, taxa with a higher proportion of intraradical allocation (such as the Glomeraceae) tend to be more competitive in carbon resource acquisition (Chagnon et al., 2013), and thus tend to be more dominant. However, while we found that Glomeraceae dominated in the control treatment, their dominance diminished with soil amendments, while the abundance of Gigasporaceae, which has a greater advantage in acquiring soil nutrient resources, increased (Figure 3). This suggests that processes other than the availability of soil nutrients more likely regulates AMF community structure. Furthermore, it is possible that the effects of soil amendments on certain AMF groups may become more pronounced over an extended period. Due to the constraints imposed by the vegetable cropping system, the duration of our study was limited to 4 months, which represents one of the key limitations of this research.
Although we cannot full disentangle mechanisms leading the observed shifts in AMF community composition with the application of organic fertilizer, the declines in soil pH, observed specifically in pepper fields, may play an important role. Soil pH regulates both nutrient availability (Hartemink and Barrow, 2023) and AMF growth and reproduction (Wang et al., 1993; Van Aarle et al., 2002). As a result, AMF taxa with variation in nutrient and growth traits (Hart and Reader, 2002; Koch et al., 2017) can shift in their relative abundances under different levels of soil pH (Davison et al., 2021; Bodenhausen et al., 2023), such that soil pH may be a key driver of changes in AMF community composition (Figures 5, 6). A lack of pH change under organic fertilizer addition in cucumber fields suggests that crop-specific factors such as root exudates, nutrient uptake patterns, or differential microbial activity might mediate the soil’s response to amendments (Vives-Peris et al., 2020; Baggs et al., 2023; Shu et al., 2023). This discrepancy highlights the complexity of plant–soil-microbe interactions and underscores the importance of considering crop-specific responses when evaluating the effects of soil amendments.
At our study, not only was the ambient nitrogen and phosphorus content in the soil high, but so was the content of Cd (average of 0.7 mg/kg), leading to high risk for Cd pollution (Chinese National Standardization Administration, 2018). Although the response of AMF to variation in Cd content is mixed (Weissenhorn et al., 1993; Krishnamoorthy et al., 2015), we found no evidence for an influence of Cd on AMF abundance. However, we did find an influence of Cd on the composition and structure of AMF communities (Figures 5, 6). While we cannot explicitly disentangle the mechanisms underlying this effect, it is likely related to the differential responses of different AMF taxa to Cd (Gai et al., 2018; Chen et al., 2023). For example, the growth of Gigasporaceae family was markedly inhibited by Cd exposure (Figure 6). Existing evidence suggests that AMF communities in low P or high nitrogen-to-phosphorus ratio environments tend to favor Gigasporaceae species likely due to their competitive life history traits (Johnson et al., 2003; Chagnon et al., 2013). This is consistent with the significant increase in their abundance following the addition of organic matter in this study, whereas their response to phosphate addition was negligible (Figure 3). These findings suggest that the response of Gigasporaceae to varying resource inputs may be largely unaffected by the inhibitory effects of Cd, indicating a potential resilience of this family to heavy metal stress in specific environmental contexts. While the two soil amendments did not significantly influence the total Cd concentration in the soil, they have the potential to mitigate plant Cd uptake by decreasing Cd bioavailability (Wang F. et al., 2022). However, the absence of data on variations in available Cd content constrains our ability to comprehensively assess the extent to which these amendments may have impacted AMF abundance and community composition through the modulation of Cd bioavailability.
The compositional shifts of AMF communities we observed indicate that the use of soil amendments can have certain ecological impacts and/or risks (van der Heijden et al., 1998; Rillig and Mummey, 2006; Hawkins et al., 2023). Specifically, the increased abundance of Gigasporaceae in response to organic matter addition may enhance nutrient use efficiency in crops, improve soil aggregation, and boost carbon sequestration due to their significant investment in extraradical hyphae (Hart and Reader, 2002; Chagnon et al., 2013). In contrast, a decline in Glomeraceae could increase the risk of excessive Cd accumulation in agricultural products, as they generally demonstrate higher remediation efficiency than Gigasporaceae (Ferrol et al., 2009). Given the significant gaps in our understanding of the functional diversity of AMF communities, it remains challenging to predict the specific consequences of these changes (Wang and Rengel, 2024; Martin and van der Heijden, 2024) and how these changes in AMF communities will influence production in heavy metal-contaminated farmland. Nevertheless, the application of AMF inoculants on Cd contaminated farmland soils is considered an effective strategy to reduce plant Cd uptake and enhance crop yields (Liu et al., 2018; Gao Y. et al., 2023). Therefore, the practical application of our study lies in our ability to select representative taxa with differential responses to soil amendments, then explore their roles in reducing plant Cd uptake, increasing plant yields and improving soil aggregation, and to ultimately identify fungal strains with potential for widespread application. However, we did find that AMF responses to soil amendments differed between the two different kinds of crops in our study, underscoring the need for considering crop-dependent strategies for AMF inoculation in agricultural soils, where targeted AMF inoculation could be tested for its efficacy in reducing Cd uptake and ensuring food safety.
5 Conclusion
We examined the response of AMF communities to the application of soil amendments intended to allow safer utilization of crops in Cd contaminated vegetable fields. We found that in two types of Cd-contaminated vegetable fields (peppers and cucumbers), the AMF community composition and structure changed with the addition of soil amendments, but no change in the overall abundance. Interestingly, this response was primarily associated with changes in soil nitrogen and the nitrogen-to-phosphorus ratio, rather than phosphorus itself, reflecting the important role of the absolute levels of nutrients in regulating AMF communities. Additionally, we preliminarily identified certain AMF taxa that were more responsive to soil amendments. These findings provide valuable insights into exploring and identifying more efficient AMF strains for Cd remediation. However, further research is needed to elucidate the mechanisms underlying these responses and determine whether the selected taxa play specific roles in regulating plant heavy metal uptake. Long-term field trials across a wider range of crops and soil conditions are crucial to ensure the scalability and reliability of these strategies for broader agricultural applications.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found: https://www.ncbi.nlm.nih.gov/, PRJNA1131945.
Author contributions
LL: Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. YS: Data curation, Formal analysis, Investigation, Writing – review & editing. WX: Investigation, Visualization, Writing – review & editing. XW: Validation, Writing – review & editing. ZX: Investigation, Methodology, Resources, Writing – review & editing. YL: Formal analysis, Writing – review & editing. XS: Conceptualization, Data curation, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Natural Science Foundation of Jiangxi Province (20224BAB205009 and 20232BAB213081), the Key Research and Development Program of Jiangxi Academy of Sciences (2022YSBG22009, 2022YSBG10005, and 2022YSBG21009) and the National Natural Science Foundation of China (32060253).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Baggs, E. M., Cairns, J. E., Mhlanga, B., Petroli, C. D., Chamberlin, J., Karwat, H., et al. (2023). Exploiting crop genotype-specific root-soil interactions to enhance agronomic efficiency. Front. Soil Sci. 3:1125604. doi: 10.3389/fsoil.2023.1125604
Bodenhausen, N., Hess, J., Valzano, A., Deslandes-Hérold, G., Waelchli, J., Furrer, R., et al. (2023). Predicting soil fungal communities from chemical and physical properties. J. Sustain. Agric. Environ. 2, 225–237. doi: 10.1002/sae2.12055
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Brundrett, M., Melville, L., and Peterson, L. (1994). Practical methods in mycorrhiza research. Sidney, British Columbia: Mycologue Publications.
Brundrett, M. C., and Tedersoo, L. (2018). Evolutionary history of mycorrhizal symbioses and global host plant diversity. New Phytol. 220, 1108–1115. doi: 10.1111/nph.14976
Cakmak, I., Lambers, H., Grant, C. A., and Zhao, F. (2023). Arbuscular mycorrhizal fungi: key players in avoiding cadmium accumulation in food crops. Plant Soil 484, 13–32. doi: 10.1007/s11104-022-05802-w
Camenzind, T., Hempel, S., Homeier, J., Horn, S., Velescu, A., Wilcke, W., et al. (2014). Nitrogen and phosphorus additions impact arbuscular mycorrhizal abundance and molecular diversity in a tropical montane forest. Glob. Chang. Biol. 20, 3646–3659. doi: 10.1111/gcb.12618
Chagnon, P.-L., Bradley, R. L., Maherali, H., and Klironomos, J. N. (2013). A trait-based framework to understand life history of mycorrhizal fungi. Trends Plant Sci. 18, 484–491. doi: 10.1016/j.tplants.2013.05.001
Chen, Y., Zhou, X., Wang, Z., Su, X., Liu, F., Tian, X., et al. (2023). Cd contamination determined assembly processes and network stability of AM fungal communities in an urban green space ecosystem. Sci. Total Environ. 899:166372. doi: 10.1016/j.scitotenv.2023.166372
Chinese National Standardization Administration (2018). GB 15618–2018: Soil environmental quality – risk control standard for soil contamination of agricultural land. Beijing: China Environment Publishing Group.
Cui, W., Liu, Y., Li, W., Pei, L., Xu, S., Sun, Y., et al. (2023). Remediation agents drive bacterial community in a cd-contaminated soil. Toxics 11:53. doi: 10.3390/toxics11010053
Davison, J., Moora, M., Semchenko, M., Adenan, S. B., Ahmed, T., Akhmetzhanova, A. A., et al. (2021). Temperature and pH define the realised niche space of arbuscular mycorrhizal fungi. New Phytol. 231, 763–776. doi: 10.1111/nph.17240
Delavaux, C. S., Smith-Ramesh, L. M., and Kuebbing, S. E. (2017). Beyond nutrients: a meta-analysis of the diverse effects of arbuscular mycorrhizal fungi on plants and soils. Ecology 98, 2111–2119. doi: 10.1002/ecy.1892
Durdu, B., Gurbuz, F., Koçyiğit, H., and Gurbuz, M. (2023). Urbanization-driven soil degradation; ecological risks and human health implications. Environ. Monit. Assess. 195:1002. doi: 10.1007/s10661-023-11595-x
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Feng, H., and Cheng, J. (2023). Whole-process risk management of soil amendments for remediation of heavy metals in agricultural soil-a review. Int. J. Environ. Res. Public Health 20:1869. doi: 10.3390/ijerph20031869
Ferrol, N., González-Guerrero, M., Valderas, A., Benabdellah, K., and Azcón-Aguilar, C. (2009). Survival strategies of arbuscular mycorrhizal fungi in cu-polluted environments. Phytochem. Rev. 8, 551–559. doi: 10.1007/s11101-009-9133-9
Gai, J. P., Fan, J. Q., Zhang, S. B., Mi, N. N., Christie, P., Li, X. L., et al. (2018). Direct effects of soil cadmium on the growth and activity of arbuscular mycorrhizal fungi. Rhizosphere 7, 43–48. doi: 10.1016/j.rhisph.2018.07.002
Gao, Y., An, T., Kuang, Q., Wu, Y., Liu, S., Liang, L., et al. (2023). The role of arbuscular mycorrhizal fungi in the alleviation of cadmium stress in cereals: a multilevel meta-analysis. Sci. Total Environ. 902:166091. doi: 10.1016/j.scitotenv.2023.166091
Gao, J., Han, H., Gao, C., Wang, Y., Dong, B., and Xu, Z. (2023). Organic amendments for in situ immobilization of heavy metals in soil: a review. Chemosphere 335:139088. doi: 10.1016/j.chemosphere.2023.139088
Gosling, P., Mead, A., Proctor, M., Hammond, J. P., and Bending, G. D. (2013). Contrasting arbuscular mycorrhizal communities colonizing different host plants show a similar response to a soil phosphorus concentration gradient. New Phytol. 198, 546–556. doi: 10.1111/nph.12169
Hao, Z., Xie, W., and Chen, B. (2022). Application of arbuscular mycorrhizal fungi in agriculture: research progress and challenges. Sci. Technol. Rev. 40, 87–98. doi: 10.3981/j.issn.1000-7857.2022.03.008
Hart, M. M., and Reader, R. J. (2002). Taxonomic basis for variation in the colonization strategy of arbuscular mycorrhizal fungi. New Phytol. 153, 335–344. doi: 10.1046/j.0028-646X.2001.00312.x
Hartemink, A. E., and Barrow, N. J. (2023). Soil pH – nutrient relationships: the diagram. Plant Soil 486, 209–215. doi: 10.1007/s11104-022-05861-z
Hawkins, H.-J., Cargill, R. I. M., Van Nuland, M. E., Hagen, S. C., Field, K. J., Sheldrake, M., et al. (2023). Mycorrhizal mycelium as a global carbon pool. Curr. Biol. 33, R560–R573. doi: 10.1016/j.cub.2023.02.027
Holland, J. E., Bennett, A. E., Newton, A. C., White, P. J., McKenzie, B. M., George, T. S., et al. (2018). Liming impacts on soils, crops and biodiversity in the UK: a review. Sci. Total Environ. 610-611, 316–332. doi: 10.1016/j.scitotenv.2017.08.020
Johnson, N. C. (2010). Resource stoichiometry elucidates the structure and function of arbuscular mycorrhizas across scales. New Phytol. 185, 631–647. doi: 10.1111/j.1469-8137.2009.03110.x
Johnson, N. C., Rowland, D. L., Corkidi, L., Egerton-Warburton, L. M., and Allen, E. B. (2003). Nitrogen enrichment alters mycorrhizal allocation at five Mesic to semiarid grasslands. Ecology 84, 1895–1908. doi: 10.1890/0012-9658(2003)084[1895:NEAMAA]2.0.CO;2
Koch, A. M., Antunes, P. M., Maherali, H., Hart, M. M., and Klironomos, J. N. (2017). Evolutionary asymmetry in the arbuscular mycorrhizal symbiosis: conservatism in fungal morphology does not predict host plant growth. New Phytol. 214, 1330–1337. doi: 10.1111/nph.14465
Krishnamoorthy, R., Kim, C.-G., Subramanian, P., Kim, K.-Y., Selvakumar, G., and Sa, T.-M. (2015). Arbuscular mycorrhizal fungi community structure, abundance and species richness changes in soil by different levels of heavy metal and metalloid concentration. PLoS One 10:e0128784. doi: 10.1371/journal.pone.0128784
Lee, J., Lee, S., and Young, J. P. W. (2008). Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. FEMS Microbiol. Ecol. 65, 339–349. doi: 10.1111/j.1574-6941.2008.00531.x
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
Liu, L., Li, J., Yue, F., Yan, X., Wang, F., Bloszies, S., et al. (2018). Effects of arbuscular mycorrhizal inoculation and biochar amendment on maize growth, cadmium uptake and soil cadmium speciation in cd-contaminated soil. Chemosphere 194, 495–503. doi: 10.1016/j.chemosphere.2017.12.025
Liu, Y., Shi, G., Mao, L., Cheng, G., Jiang, S., Ma, X., et al. (2012). Direct and indirect influences of 8 yr of nitrogen and phosphorus fertilization on Glomeromycota in an alpine meadow ecosystem. New Phytol. 194, 523–535. doi: 10.1111/j.1469-8137.2012.04050.x
Liu, H., Wang, H., Nie, Z., Tao, Z., Peng, H., Shi, H., et al. (2024). Combined application of arbuscular mycorrhizal fungi and selenium fertilizer increased wheat biomass under cadmium stress and shapes rhizosphere soil microbial communities. BMC Plant Biol. 24:359. doi: 10.1186/s12870-024-05032-5
Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Maherali, H., and Klironomos, J. N. (2012). Phylogenetic and trait-based assembly of arbuscular mycorrhizal fungal communities. PLoS One 7:e36695. doi: 10.1371/journal.pone.0036695
Martin, F. M., and van der Heijden, M. G. A. (2024). The mycorrhizal symbiosis: research frontiers in genomics, ecology, and agricultural application. New Phytol. 242, 1486–1506. doi: 10.1111/nph.19541
McGonigle, T., Miller, M., Evans, D., Fairchild, G., and Swan, J. (1990). A new method which gives an objective measure of colonization of roots by vesicular—arbuscular mycorrhizal fungi. New Phytol. 115, 495–501. doi: 10.1111/j.1469-8137.1990.tb00476.x
McLaughlin, M. J., Smolders, E., Zhao, F., Grant, C., and Montalvo, D. (2021). Managing cadmium in agricultural systems. In Advances in Agronomy. ed. D. L. Sparks Cambridge, MA : Academic Press, 166, 1–129.
Mohamed, O., Said, E. K., Miloud, S., Abdellatif, H., El Hassan, A., and Rachid, B. (2023). Effect of agricultural management practices on diversity, abundance, and infectivity of arbuscular mycorrhizal fungi: a review. Symbiosis 91, 33–44. doi: 10.1007/s13199-023-00948-2
Moreno Jiménez, E., Ferrol, N., Corradi, N., Peñalosa, J. M., and Rillig, M. C. (2024). The potential of arbuscular mycorrhizal fungi to enhance metallic micronutrient uptake and mitigate food contamination in agriculture: prospects and challenges. New Phytol. 242, 1441–1447. doi: 10.1111/nph.19269
Mubeen, S., Ni, W., He, C., and Yang, Z. (2023). Agricultural strategies to reduce cadmium accumulation in crops for food safety. Agriculture 13:471. doi: 10.3390/agriculture13020471
Oksanen, J., Blanchet, F., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2017). Vegan: community ecology package. R package v.2.4-3.
Öpik, M., Vanatoa, A., Vanatoa, E., Moora, M., Davison, J., Kalwij, J., et al. (2010). The online database Maarj AM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota). New Phytol. 188, 223–241. doi: 10.1111/j.1469-8137.2010.03334.x
Park, S. H., Kang, B. R., Kim, J., Lee, Y., Nam, H. S., and Lee, T. K. (2024). Enhanced soil fertility and carbon dynamics in organic farming systems: the role of arbuscular mycorrhizal fungal abundance. J. Fungi 10:598. doi: 10.3390/jof10090598
R Core Team (2022). R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.
Rashid, A., Schutte, B. J., Ulery, A., Deyholos, M. K., Sanogo, S., Lehnhoff, E. A., et al. (2023). Heavy metal contamination in agricultural soil: environmental pollutants affecting crop health. Agronomy 13:1521. doi: 10.3390/agronomy13061521
Rillig, M. C., and Mummey, D. L. (2006). Mycorrhizas and soil structure. New Phytol. 171, 41–53. doi: 10.1111/j.1469-8137.2006.01750.x
Santos, J. C., Finlay, R. D., and Tehler, A. (2006). Molecular analysis of arbuscular mycorrhizal fungi colonising a semi-natural grassland along a fertilisation gradient. New Phytol. 172, 159–168. doi: 10.1111/j.1469-8137.2006.01799.x
Sato, K., Suyama, Y., Saito, M., and Sugawara, K. (2005). A new primer for discrimination of arbuscular mycorrhizal fungi with polymerase chain reaction-denature gradient gel electrophoresis. Grassl. Sci. 51, 179–181. doi: 10.1111/j.1744-697X.2005.00023.x
Shi, T., Zhang, Y., Gong, Y., Ma, J., Wei, H., Wu, X., et al. (2019). Status of cadmium accumulation in agricultural soils across China (1975–2016): from temporal and spatial variations to risk assessment. Chemosphere 230, 136–143. doi: 10.1016/j.chemosphere.2019.04.208
Shu, X., Liu, W., Huang, H., Ye, Q., Zhu, S., Peng, Z., et al. (2023). Meta-analysis of organic fertilization effects on soil bacterial diversity and community composition in agroecosystems. Plan. Theory 12:3801. doi: 10.3390/plants12223801
Singh, S. K., Wu, X., Shao, C., and Zhang, H. (2022). Microbial enhancement of plant nutrient acquisition. Stress Biol. 2:3. doi: 10.1007/s44154-021-00027-w
Six, L., and Smolders, E. (2014). Future trends in soil cadmium concentration under current cadmium fluxes to European agricultural soils. Sci. Total Environ. 485-486, 319–328. doi: 10.1016/j.scitotenv.2014.03.109
Smith, S. E., and Smith, F. A. (2011). Roles of arbuscular mycorrhizas in plant nutrition and growth: new paradigms from cellular to ecosystem scales. Annu. Rev. Plant Biol. 62, 227–250. doi: 10.1146/annurev-arplant-042110-103846
Tennant, D. (1975). A test of a modified line intersect method of estimating root length. J. Ecol. 63, 995–1001. doi: 10.2307/2258617
Treseder, K. K. (2004). A meta-analysis of mycorrhizal responses to nitrogen, phosphorus, and atmospheric CO2 in field studies. New Phytol. 164, 347–355. doi: 10.1111/j.1469-8137.2004.01159.x
Van Aarle, I. M., Olsson, P. A., and Söderström, B. (2002). Arbuscular mycorrhizal fungi respond to the substrate pH of their extraradical mycelium by altered growth and root colonization. New Phytol. 155, 173–182. doi: 10.1046/j.1469-8137.2002.00439.x
van der Heijden, M. G., Klironomos, J. N., Ursic, M., and Moutoglis, P. (1998). Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396, 69–72. doi: 10.1038/23932
van der Heyde, M., Ohsowski, B., Abbott, L. K., and Hart, M. (2017). Arbuscular mycorrhizal fungus responses to disturbance are context-dependent. Mycorrhiza 27, 431–440. doi: 10.1007/s00572-016-0759-3
Vives-Peris, V., de Ollas, C., Gómez-Cadenas, A., and Pérez-Clemente, R. M. (2020). Root exudates: from plant to rhizosphere and beyond. Plant Cell Rep. 39, 3–17. doi: 10.1007/s00299-019-02447-5
Wang, F. (2017). Occurrence of arbuscular mycorrhizal fungi in mining-impacted sites and their contribution to ecological restoration: mechanisms and applications. Crit. Rev. Environ. Sci. Technol. 47, 1901–1957. doi: 10.1080/10643389.2017.1400853
Wang, F., Cheng, P., Zhang, S., Zhang, S., and Sun, Y. (2022). Contribution of arbuscular mycorrhizal fungi and soil amendments to remediation of a heavy metal-contaminated soil using sweet sorghum. Pedosphere 32, 844–855. doi: 10.1016/j.pedsph.2022.06.011
Wang, F., and Rengel, Z. (2024). Disentangling the contributions of arbuscular mycorrhizal fungi to soil multifunctionality. Pedosphere 34, 269–278. doi: 10.1016/j.pedsph.2023.12.015
Wang, G. M., Stribley, D. P., Tinker, P. B., and Walker, C. (1993). Effects of pH on arbuscular mycorrhiza I. Field observations on the long-term liming experiments at Rothamsted and Woburn. New Phytol. 124, 465–472. doi: 10.1111/j.1469-8137.1993.tb03837.x
Wang, H., Zhao, X., Zhang, J., Lu, C., and Feng, F. (2022). Arbuscular mycorrhizal fungus regulates cadmium accumulation, migration, transport, and tolerance in Medicago sativa. J. Hazard. Mater. 435:129077. doi: 10.1016/j.jhazmat.2022.129077
Webb, C. O., Ackerly, D. D., McPeek, M. A., and Donoghue, M. J. (2002). Phylogenies and community ecology. Annu. Rev. Ecol. Syst. 33, 475–505. doi: 10.1146/annurev.ecolsys.33.010802.150448
Weissenhorn, I., Leyval, C., and Berthelin, J. (1993). Cd-tolerant arbuscular mycorrhizal (AM) fungi from heavy-metal polluted soils. Plant Soil 157, 247–256. doi: 10.1007/BF00011053
Keywords: arbuscular mycorrhizal fungi, soil amendments, heavy metals, cadmium, vegetable fields
Citation: Li L, Shi Y, Xia W, Wang X, Xin Z, Liao Y and Sun X (2024) Soil amendments altered arbuscular mycorrhizal fungal communities in cadmium-contaminated vegetable fields. Front. Microbiol. 15:1470137. doi: 10.3389/fmicb.2024.1470137
Edited by:
Xian-Chun Zeng, China University of Geosciences Wuhan, ChinaReviewed by:
Fayuan Wang, Qingdao University of Science and Technology, ChinaMariana Erasmus, University of the Free State, South Africa
Shalini Dhiman, Guru Nanak Dev University, India
Jinglong Li, The University of Queensland, Australia
Copyright © 2024 Li, Shi, Xia, Wang, Xin, Liao and Sun. 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) and the copyright owner(s) 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: Zaijun Xin, emFpanVueGluQDEyNi5jb20=; Xiaoyan Sun, eGlhb3lhbl9zdW4wNUAxNjMuY29t