3.1 Plant growth and biomass production
Plant leaf area increased with increasing CA concentration and decreased with increasing TA and MA concentrations (Fig. 1). Under Cd stress (CK), root, shoot, and plant biomass decreased by 64.66%, 24.86%, and 39.19%, respectively. After the application of CA at different concentrations, root, shoot, and plant biomass increased significantly compared with CK, namely by 153.96%, 38.39%, and 64.29%, respectively, with the addition of 15 mmol/kg CA. However, the addition of TA and MA at low concentrations significantly facilitated plant growth. Root, shoot, and plant biomass increased by 107.74%, 14.27%, and 31.86%, respectively, when 1 mmol/kg TA was added, and increased by 57.86%, 38.26%, and 45.03%, respectively, after the addition of 1 mmol/kg MA compared with CK. Low concentrations of CA, TA, and MA could also stimulate plant stem growth, which increased by 29.81%, 28.45%, and 49.81%, respectively, under the addition of 1 mmol/kg CA, TA, and MA.
3.2 Chlorophyll fluorescence parameters and fast light response curves
The addition of LMWOAs with higher concentration could significantly improve the plant chlorophyll content (Fig. S1). The F0 and Fm values under Cd stress increased by 74.63% and 6.33%, respectively, whereas the Fv/Fm and Fv/F0 values decreased by 31.53% and 63.84%, respectively (Fig. 2a). The addition of CA could reduce the values of F0 and Fm and elevate those of Fv/Fm and Fv/F0. At a dose of 15 mmol/kg CA, the values of F0 and Fm decreased by 31.62% and 39.68%, whereas the Fv/Fm and Fv/F0 values were 1.36 and 2.18 times those of CK. The values of F0 and Fm increased continuously, and those of Fv/Fm and Fv/F0 decreased gradually with the application of TA and MA from low to high concentrations. Compared with CK, the F0 and Fm values at 15 mmol/kg TA increased by 13.68% and 65.78%, whereas the Fv/Fm and Fv/F0 values decreased by 13.04% and 23.14%, respectively. Meanwhile, the values of F0 and Fm increased by 25.64% and 55.57%, and that of Fv/Fm and Fv/F0 decreased by 6.52% and 12.31% with the addition of 15 mmol/kg MA. According to the Fv/Fm image (Fig. 2a), CKK was dark blue without Cd stress, and CK showed a light green color under Cd stress. With the application of CA, TA, and MA from low to high concentrations, the blue color deepened in the CA treatment, whereas the green color deepened in the treatments with TA and MA.
The fast light response curves were also measured, and the values of ETR, Y(II), and Y(NPQ) decreased by 39.61%, 35.60% and 30.70%, respectively, at the light intensity of 231 µmol m− 2·s− 1 under Cd stress (Fig. 2b, 2c, 2d, 2e). The addition of CA increased the values of ETR, Y(II), and (NPQ) by 53.50%, 51.42%, and 38.70%, respectively, at 15 mmol/kg CA. However, the application of TA and MA at low concentrations increased the values of ETR, Y(II), and Y(NPQ) by 4.35%, 9.29%, and 1.57%, respectively, when 1 mmol/kg TA was added and by 17.96%, 18.81%, and 1.59%, respectively, for 1 mmol/kg MA compared with CK. With increasing CA concentrations, the Y(NO) value decreased gradually by 41.06% at 15 mmol/kg CA. However, the value of Y(NO) increased with the increase of TA and MA concentrations by 21.64% and 15.91%, respectively, with the addition of 15 mmol/kg TA and MA.
3.3 Cd uptake by the plants
Metal accumulation was evaluated to reveal the effects of LMWOAs on Cd uptake in Bidens pilosa L. The Cd accumulated in roots and shoots increased gradually with the application of CA, TA, and MA from low to high concentrations. The Cd concentrations in roots and shoots were 75.22 and 109.49 mg/kg in the treatment of CK (Fig. 3a, 3b). The Cd concentrations in roots were 109.65, 152.07, and 108.07 mg/kg, and that in shoots were 186.38, 161.94, and 155.30 mg/kg, respectively, when 15 mmol/kg CA, TA, and MA were added. The transport, absorption, and enrichment of Cd by Bidens pilosa L. varied greatly with the application of LMWOAs. The TF value increased continuously with increasing CA concentrations, with an increase by 16.45% with the addition of 15 mmol/kg CA (Table 1). However, the TF value increased by 12.61% and 5.73%, respectively, with the addition of 1 mmol/kg TA and MA and decreased by 27.03% and 1.55%, respectively, after the application of 15 mmol/kg TA and MA. The PR value increased gradually with increasing CA concentrations and was 169.21% higher in the treatment with the addition of 15 mmol/kg CA. Higher PR values were found in the plants treated with 5 mmol/kg TA and MA, increasing by 71.28% and 63.11%, respectively. The BCF value increased continuously with increasing CA, TA, and MA concentrations, with 77.98%, 78.33%, and 64.49% higher levels, respectively, with the addition of 15 mmol/kg.
Table 1
Effects of exogenous LMWOAs addition on the bioaccumulation, transportation and extraction of heavy metal
| TF | RE | BCF |
CK | 1.460 ± 0.111a | 0.100 ± 0.001a | 2.808 ± 0.022a |
CA1 | 1.555 ± 0.136acd | 0.125 ± 0.006b | 2.749 ± 0.091a |
CA5 | 1.606 ± 0.103bcd | 0.183 ± 0.014c | 4.039 ± 0.081b |
CA15 | 1.700 ± 0.005b | 0.268 ± 0.006d | 4.998 ± 0.143c |
TA1 | 1.644 ± 0.086bc | 0.142 ± 0.004e | 3.219 ± 0.015d |
TA5 | 1.481 ± 0.003ad | 0.171 ± 0.001f | 4.353 ± 0.115e |
TA15 | 1.065 ± 0.014e | 0.168 ± 0.007fg | 5.007 ± 0.283c |
MA1 | 1.544 ± 0.031acd | 0.126 ± 0.006b | 2.655 ± 0.080a |
MA5 | 1.530 ± 0.071acd | 0.163 ± 0.001g | 3.997 ± 0.172b |
MA15 | 1.437 ± 0.011a | 0.141 ± 0.001e | 4.619 ± 0.165f |
3.4 Speciation distribution of Cd in the soil
The acid-soluble, reducible, oxidizable, and residual states of Cd accounted for 38.60%, 12.73%, 12.89%, and 35.79%, respectively, in CK (Fig. 4). The addition of LMWOAs significantly affected the speciation distribution of Cd in the soil and enhanced the proportions of the acid-soluble state. The proportions of acid-soluble Cd in the soil increased by 10.26%, 6.78%, and 6.58%, respectively, whereas those of the reducible state decreased by 3.00%, 2.97%, and 0.95%. The proportions of the residual state decreased by 8.00%, 4.07%, and 5.86%, respectively, when 15 mmol/kg CA, TA, and MA were added. However, no significant difference was found in the proportions of oxidizable Cd with the application of LMWOAs.
3.5 Soil bacterial community structure
To explore the influences of the three LMWOAs on the rhizosphere bacterial composition under Cd stress, bacterial 16Sr DNA was sequenced in the experiments treated with 15 mmol/kg CA, TA, and MA. The Chao, Ace, and Shannon index values of the rhizosphere bacterial communities decreased by 9.79%, 8.39%, and 3.16%, respectively, under Cd stress (Table 2). The application of 15 mmol/kg CA, TA, and MA reduced the diversity and abundance of rhizosphere bacteria, and the inhibitory effect of CA was the weakest. The values of the Chao, Ace, and Shannon indices decreased by 6.39%, 5.38%, and 4.28%, respectively with CA addition, and decreased by 21.85%, 21.59%, and 11.28%, respectively, with TA application. When MA was added, they decreased by 11.68%, 10.57%, and 5.15%, respectively. The percentages of Proteobacteria, Bacteroidoia, and Gemmatimonadota of rhizosphere bacteria decreased by 10.16%, 2.78%, and 2.21%, whereas those of Acidobacteriota and Chloroflexi increased by 8.96% and 4.84%, respectively, under Cd stress (Fig. 5a). With the application of LMWOAs, the proportions of Proteobacteria and Bacteroidota increased significantly, whereas those of Acidobacteriota, Chloroflexi, and Gemmatimonadota decreased compared to CK.
Table 2
Effects of exogenous LMWOAs addition on OTU abundance and diversity index of bacteria communities in Cd contaminated soil
| chao | ace | shannon | simpson | coverage |
CKK | 2413.21 ± 9.92b | 2397.64 ± 3.95b | 5.70 ± 0.09a | 0.013 ± 0.002a | 0.988 ± 0.002a |
CK | 2176.94 ± 125.95a | 2196.52 ± 80.93a | 5.52 ± 0.27a | 0.016 ± 0.005a | 0.986 ± 0.003a |
CA | 2037.79 ± 82.64c | 2078.40 ± 54.09c | 5.28 ± 0.04b | 0.018 ± 0.001a | 0.988 ± 0.001a |
TA | 1706.72 ± 0.54e | 1716.66 ± 13.37e | 4.90 ± 0.22c | 0.032 ± 0.008b | 0.989 ± 0.001a |
MA | 1922.53 ± 48.46d | 1964.35 ± 34.29d | 5.23 ± 0.14b | 0.017 ± 0.002a | 0.987 ± 0.001a |
The proportions of the genera Ensifer, Ramlibacter, Devosia, and Pseudoxanthomonas in the rhizosphere soil were reduced due to Cd stress, whereas they were enhanced by the application of LMWOAs (Fig. 5b). Interestingly, the abundances of Actinophytocola, Azotobacter, Pseudomonas, and Ramlibacter, which were extremely low in the CKK and CK treatments, were high in the treatments with LMWOA addition. The three LMWOAs had different impacts on the bacterial genera in the rhizosphere. The addition of CA could significantly increase the abundances of Azotobacter and Pseudomonas, whereas the abundance levels of Actinophytocola and Ensifer were improved in the TA treatments.
3.6 Redundancy analysis
The CA, TA, and MA levels in Bidens Pilosa L. varied with the addition of exogenous LMWOAs, and the contents in the plant leaves were determined (Fig. S2). The antioxidant levels (CAT, SOD, GSH), osmotic adjustment substance levels (Tsp, Pro), and oxidative damage index levels (MDA, H2O2, O2−) of Bidens Pilosa L. were also measured after exogenous LMWOA addition (Figs. S3, S4, S5). The correlation between LMWOAs in plants and the related parameters was described by redundancy analysis (Fig. 6a). In the RDA diagram, axis 1 explained 58.10% of the variation of the physiological index, and axis 2 explained 16.48% of the variation. In the photosynthetic system, CA and MA addition was significantly negatively correlated with Fm and F0 and positively correlated with F0/Fm and Fv/F0. In the antioxidant system, CA and MA addition was significantly positively correlated with SOD and GSH and negatively correlated with MDA, H2O2, and O2−. In the osmotic regulation system, CA and MA additions were significantly positively correlated with total soluble protein and proline levels.
The contents of CA, MA, and TA in plant roots and rhizosphere soil enzyme activity (sucrase activity, urease activity, protease activity, phosphatase activity, catalase activity) were also determined (Fig. S6, S7), and there were correlations between LMWOAs in plant roots and soil enzyme activities and bacterial abundance (Fig. 6b). In the RDA diagram, axis 1 explained 59.25% of the variation in the physiological index, and axis 2 explained 13.83% of this variation. The addition of CA and MA was positively correlated with protease, urease, and phosphatase levels, whereas TA addition was positively correlated with sucrose levels. The application of CA, TA, and MA could reduce catalase activity in the soil. Based on Pearson´s correlation heat map analysis, LMWOAs and soil enzymes had similar effects on bacterial phyla and genera (Fig. 6c, 6d), and it can be concluded that CA, MA, and protease have similar effects on bacterial communities. The effect of TA on the bacterial community was similar to those of urease, phosphatase, and sucrase.