Effects of various biochar treatments on lettuce growth
The effects of cadmium (Cd) stress, in conjunction with raw biochar (BC) and cystamine-modified biochar (Cys-BC) treatments, on lettuce growth are presented in Table 1. Cd exposure significantly inhibits lettuce growth, particularly at the highest concentration of 5 mg kg− 1. This inhibition is evident through reductions in plant height, leaf length, and fresh weight. Under high Cd stress (T0 − 5), fresh weight decreases by 51.77% compared to the control group (CK), which is not exposed to Cd. However, treatments raw BC (C1.2−b) and Cys-BC (T1.2−b) mitigate these adverse effects, resulting in increases in plant fresh weight of 18.99% and 43.07%, respectively, compared to the T0 − 5 treatment.
Table 1
Effects of different biochar treatments on lettuce growth parameters.
Treatment | Plant height (cm) | Stem diameter (mm) | Leaf number | Leaf width (cm) | Leaf length (cm) | Fresh weight (g·plant− 1) |
CK | 23.80 ± 0.24a | 9.85 ± 0.24A | 13.40 ± 0.80a | 9.80 ± 0.22A | 12.7 ± 0.28a | 143.28 ± 1.04A |
T0 − 1 | 16.81 ± 0.47e | 5.58 ± 0.18C | 11.20 ± 0.75b | 6.72 ± 0.56C | 9.52 ± 0.54d | 88.28 ± 0.84F |
T0 − 3 | 16.09 ± 0.19ef | 5.34 ± 0.04C | 10.80 ± 0.40b | 6.22 ± 0.10C | 8.31 ± 0.29e | 79.71 ± 1.19H |
T0 − 5 | 15.78 ± 0.10f | 5.09 ± 0.17C | 10.40 ± 0.49b | 5.85 ± 0.3C | 7.07 ± 0.24f | 69.93 ± 1.41I |
C1.2−1 | 22.36 ± 0.34cd | 9.06 ± 0.27B | 12.80 ± 0.4a | 8.88 ± 0.35B | 11.85 ± 0.65abc | 108.61 ± 1.12D |
C1.2−3 | 22.77 ± 0.57bc | 9.07 ± 0.34B | 12.60 ± 0.49a | 8.84 ± 0.51B | 11.35 ± 0.48bc | 97.67 ± 0.71E |
C1.2−5 | 21.65 ± 0.50d | 9.06 ± 0.56B | 12.80 ± 0.4a | 8.53 ± 0.27B | 10.97 ± 0.33c | 83.21 ± 1.61G |
T1.2−1 | 23.48 ± 0.42ab | 9.55 ± 0.49AB | 13.40 ± 0.49a | 9.08 ± 0.61AB | 11.92 ± 0.47abc | 123.19 ± 0.96B |
T1.2−3 | 22.58 ± 0.42c | 9.28 ± 0.24AB | 13.00 ± 0.63a | 9.31 ± 0.52AB | 12.23 ± 0.28ab | 111.64 ± 1.5C |
T1.2−5 | 22.60 ± 0.31c | 8.96 ± 0.2B | 13.00 ± 0.63 | 9.24 ± 0.55AB | 11.25 ± 0.86bc | 100.05 ± 1.18E |
CK group represents the control group without Cd or biochar added to the soil. The T/Ca−b treatments refer to the use of cystamine-modified (T) and raw (C) biochar, with "a" representing the biochar addition level (0 ~ 1.2%), and "b" representing the concentration of added Cd (1, 3, 5 mg kg− 1) in the soil. |
Effects of various biochar treatments on root morphology
Scanning images of the lettuce root system under various biochar treatments (Fig. S1) reveal significant reductions in lateral root development under different levels of Cd stress. Table 2 presents a quantitative analysis of root morphology, including total root length, surface area, volume, and tip count. In the T0 − 5 treatment, total root length decreases by 48.72% compared to the control (CK) during early growth. However, both raw BC and Cys-BC treatments result in notable improvements. The Cys-BC treatments significantly increase root surface area and volume by 1.06 and 2.32 times, respectively, compared to the T0 − 5 treatment. With prolonged biochar application, all root parameters further improve, particularly in the Cys-BC treatments, where total root length, surface area, and volume are 1.10, 1.09, and 1.07 times higher, respectively, compared to the T0 − 5 treatments.
Table 2
Quantitative analysis of root morphology in lecture under different biochar treatments.
Days | Treatment | Root length (cm) | Root surface area (cm2) | Total root volume (cm3) | Root diameter (mm) | Tips count |
14 | CK | 123.28 ± 5.47g | 48.42 ± 2.46fg | 1.541 ± 0.03de | 1.035 ± 0.11abc | 376.33 ± 28.99a |
T0 − 5 | 63.32 ± 3.21h | 25.07 ± 1.35h | 0.787 ± 0.13e | 0.858 ± 0.146bc | 159.53 ± 23.71h |
C1.2−5 | 115.37 ± 7.91g | 35.06 ± 1.84gh | 0.848 ± 0.13e | 0.967 ± 0.163abc | 275.32 ± 12.53g |
T1.2−5 | 138.29 ± 4.09fg | 45.72 ± 3.54fgh | 1.203 ± 0.44de | 1.052 ± 0.127abc | 398.67 ± 58.62ef |
21 | CK | 186.01 ± 8.45de | 62.10 ± 4.29ef | 1.650 ± 0.66de | 1.083 ± 0.147abc | 504.53 ± 34.25d |
T1.2−5 | 167.73 ± 10.10ef | 44.17 ± 2.78fgh | 0.925 ± 0.46e | 0.838 ± 0.164c | 467.52 ± 42.53e |
C1.2−5 | 202.26 ± 13.54de | 73.49 ± 5.98de | 2.104 ± 0.88cd | 1.145 ± 0.131abc | 427.53 ± 62.53e |
T1.2−5 | 214.19 ± 17.93d | 90.95 ± 4.35cd | 3.074 ± 0.50bc | 1.352 ± 0.154a | 566.32 ± 32.53cd |
28 | CK | 290.54 ± 15.56c | 112.11 ± 10.87bc | 3.442 ± 0.57b | 1.228 ± 0.112abc | 725.62 ± 42.53ab |
T0 − 5 | 222.25 ± 10.35d | 91.525 ± 8.32cd | 2.999 ± 0.49bc | 1.311 ± 0.224ab | 532.43 ± 62.53d |
C1.2−5 | 361.00 ± 19.47b | 123.61 ± 14.77b | 3.368 ± 0.44b | 1.090 ± 0.151abc | 640.65 ± 52.53b |
T1.2−5 | 467.37 ± 25.63a | 191.04 ± 12.60a | 6.214 ± 1.45a | 1.301 ± 0.197a | 828.41 ± 62.53a |
Effects of various biochar treatments on photosynthesis
Table 3 illustrates a negative correlation between increasing Cd concentrations and chlorophyll content in lettuce leaves. Chlorophyll-a levels decrease by 17.53–46.75%, while chlorophyll-b levels decline by 19.67–54.10%. These reductions reflect the plant's adaptive response to Cd stress, as it modulates chlorophyll composition to mitigate the toxic effects of Cd. The significant decline in chlorophyll-b indicates a strategic adaptation aimed at maintaining photosynthetic efficiency under Cd stress. Both raw and Cys-BC treatments contribute to restoring or maintaining the standard chlorophyll a/b ratio during Cd-induced stress. The C1.2−3 treatment increases chlorophyll a and b levels by 13.86% and 16.67%, respectively, compared to the T0 − 3 control. More pronounced improvements are observed in the T1.2−5 treatment, with chlorophyll-a increasing by 51.22% and chlorophyll-b by 67.86%, relative to the T0 − 5 control.
Table 3
Effects of different biochar treatments on chlorophyll pigment content in lettuce.
Treatments | Chlorophyll a content (mg g− 1 FW) | Chlorophyll b content (mg g− 1 FW) | Chl a/ Chl b |
CK | 1.54 ± 0.14a | 0.61 ± 0.12A | 2.52 ± 0.11bc |
T0 − 1 | 1.27 ± 0.08g | 0.49 ± 0.02EF | 2.59 ± 0.01d |
T0 − 3 | 1.01 ± 0.07h | 0.36 ± 0.11FG | 2.81 ± 0.03e |
T0 − 5 | 0.82 ± 0.04i | 0.28 ± 0.04G | 2.93 ± 0.07e |
C1.2−1 | 1.35 ± 0.08c | 0.53 ± 0.05CD | 2.55 ± 0.10bc |
C1.2−3 | 1.15 ± 0.04e | 0.42 ± 0.02EF | 2.74 ± 0.09cd |
C1.2−5 | 0.93 ± 0.02f | 0.32 ± 0.01G | 2.91 ± 0.14a |
T1.2−1 | 1.51 ± 0.12b | 0.60 ± 0.01AB | 2.52 ± 0.04bc |
T1.2−3 | 1.47 ± 0.10c | 0.57 ± 0.02BC | 2.58 ± 0.07bc |
T1.2−5 | 1.24 ± 0.08d | 0.47 ± 0.01DE | 2.64 ± 0.12ab |
The effects of different biochars on the photosynthetic parameters of the lecture are presented in the supplementary materials. Fig. S2 and S3 illustrate the effectiveness of biochar treatments in alleviating the negative impacts of Cd stress on lettuce photosynthesis, with Cys-BC exhibiting superior performance.
Effect of various biochar treatments on oxidative stress markers
Malondialdehyde (MDA)
Figure 1a shows that MDA levels significantly increase under Cd stress (T0 − 5). In the absence of biochar, MDA concentrations rise to 2.17 and 3.02 times in above- and below-ground parts, compared to the control (CK). However, the introduction of biochar treatments-specifically, raw biochar (C1.2−b) and Cys-BC (T1.2−b)-markedly reduces MDA levels. In below-ground parts, MDA levels are reduced by 31.25% and 55.28%, respectively, and in above-ground parts by 25.24% and 43.57%, demonstrating the potential of Cys-BC in mitigating oxidative stress.
Hydrogen peroxide (H₂O₂)
Figure 1b indicates that H₂O₂ levels are 2.18 and 2.33 times higher under Cd stress compared to CK in above- and below-ground parts, respectively. The T1.2−b treatment reduces H₂O₂ levels by 46.45% in above-ground parts and 56.62% in below-ground parts. This highlights the efficacy of Cys-BC in reducing oxidative damage induced by Cd stress.
Glutathione (GSH)
Figure 1c shows a decrease in GSH content under Cd stress, with levels 2.38-fold and 2.02-fold lower than CK in above- and below-ground parts, respectively. T1.2−b treatment increases GSH levels by 1.30, 2.63, and 1.74 times in below-ground parts, compared to CK, T0 − b, and C1.2−b treatments, respectively. This increase suggests that Cys-BC plays a vital role in enhancing the plant’s antioxidant defense system under Cd stress.
Cysteine (Cys)
Figure 1d details that Cd stress induces a substantial increase in cysteine content, with levels in below-ground parts rising 3.08-fold compared to CK. The T1.2−b treatment further amplifies Cys content, with increases of 4.71, 1.09, and 1.32 times in above-ground parts and 13.86, 1.42, and 2.31 times in below-ground parts, compared to CK, T0 − b, and C1.2−b treatments. This indicates that Cys-BC contributes to boosting cysteine production, a key player in plant defense against Cd toxicity.
Effects of various biochar treatments on antioxidant enzyme activity
Superoxide Dismutase (SOD)
As illustrated in Fig. 2a, Cd stress (T0 − 5) results in a significant reduction in SOD activity, with decreases of 55.17% in above-ground parts and 61.77% in below-ground parts, compared to control (CK). In contrast, the Cys-BC treatment (T1.2−5) leads to substantial increases in SOD activity, with an increase of 73.83% in above-ground parts and 99.39% in below-ground parts. These findings highlight the beneficial effect of Cys-BC in enhancing the enzymatic defense against oxidative stress.
Peroxidase (POD)
Figure 2b illustrates a reduction in POD activity of 47.77% and 47.20% in the above-ground and below-ground parts under Cd stress (T0 − 5). The application of raw biochar treatment (C1.2−5) results in a slight improvement in POD activity, with increases of 35.60% in the above-ground and 23.79% in the below-ground parts. Notably, Cys-BC (T1.2−b) demonstrates the most significant enhancement, with POD activity increased by 64.37% in the above-ground and 70.42% in the below-ground parts, highlighting its role in enhancing the plant’s antioxidant capacity.
Catalase (CAT) Activity: Under Cd stress, CAT activity decreases by 56.15% in above-ground tissues and 47.50% in below-ground tissues, as seen in Fig. 2c. Raw biochar (C1.2−5) enhances CAT activity, resulting in increases of 58.56% in above-ground tissues and 12.12% in below-ground tissues. However, Cys-BC (T1.2−5) demonstrate the most significant, increasing CAT activity by 86.30% in above-ground tissues and 71.52% in below-ground tissues. This underscores its effectiveness in mitigating oxidative damage through enhanced CAT activity.
Glutathione reductase (GR)
Figure 2d illustrates that Cd stress has a more pronounced effect on GR activity in below-ground tissues, resulting in reductions of 13.49% in above-ground parts and 45.31% in below-ground parts compared to CK. Raw biochar (C1.2−5) moderately enhances GR activity, increasing it by 31.62% in above-ground parts and 15.13% in below-ground parts. The Cys-BC treatment (T1.2−5) shows the most substantial improvement, with GR activity rising by 67.18% in above-ground parts and 50.39% in below-ground parts. This further supports the capacity of Cys-BC to strengthen the plant's enzymatic defense mechanisms.
Effects of various biochar treatments on cadmium transfer rate in lettuce
Comparing Fig. 3a and 3b reveals that Cd accumulates predominantly in the roots of lettuce, exhibiting a higher transfer rate from soil to roots than from roots to shoots. The Cys-BC treatment is particularly effective in reducing Cd transfer rates, as evidenced by the lower slopes and intercepts in Fig. 3c and 3d compared to those observed with. The root-to-shoot Cd translocation slope ratio under the Cys-BC treatment is 3.03, significantly higher than the 2.63 ratio observed with raw biochar. This indicates that Cys-BC effectively limits Cd movement to the edible parts of the plant, thereby enhancing food safety.
Effects of various biochar treatments on soil properties
Figure 4a and 4b illustrate how the addition of raw BC and Cys-BC alters soil Cd fractions. Both treatments lead to a significant reduction in the proportions of acid-extractable, oxidizable, and reducible Cd forms, while markedly increasing the residual Cd fraction. The residual Cd content rises from 15% ± 0.8% in the non-biochar treatments to 44% ± 0.5% with raw BC treatments, and further increases to 67% ± 1% with Cys-BC treatments. These findings indicate that Cys-BC is more effective than raw BC in converting active Cd into stable forms, achieving a 50% reduction in acid-extractable and reducible Cd compared to raw BC.
Figure 4c and 4d illustrate the effects on soil pH and cation exchange capacity (CEC). Both raw BC and Cys-BC significantly enhance soil pH. Raw BC raises soil pH by 0.25 to 0.79 units, while Cys-BC increases it by 0.10 to 0.32 units. Furthermore, both treatments improve soil CEC. The application of raw BC results in an increase in CEC of 3.69 and 5.23 mol kg− 1 at dosages of 1.2% and 1.8%, respectively. Cys-BC demonstrates an even greater effect, with increases of 5.79 and 6.37 mol kg− 1 at the same dosages.
Table 4 presents the linear regression models that correlate soil pH and CEC with DTPA-extractable Cd. In treatments with raw BC, soil pH has a substantial influence on DTPA-extractable Cd (P = 0.001–0.003), as indicated by a high correlation coefficient (R² = 0.995–0.999). In contrast, in treatments with Cys-BC treatments, the influence of soil pH on available Cd is less significant (P = 0.073–0.135), with a lower correlation (R²=0.747–0.859). While soil CEC correlates with DTPA-extractable Cd in both treatments, the correlation is stronger in Cys-BC treatments (R² = 0.9865–0.9972, P = 0.0014–0.0068). These results suggest that raw BC primarily reduces Cd availability by increasing soil pH, whereas Cys-BC exerts its effect mainly by enhancing soil CEC. The improved CEC in Cys-BC treatments may be attributed to the structural properties of Cys-BC, particularly the abundance of functional groups such as amino, thiol, aldehyde, and hydroxyl groups on its surface, which enhance its cation exchange capacity.
Table 4
Linear regression models correlating soil pH and CEC with DTPA-extractable Cd.
Soil | Treatments | Cd levels (mg kg− 1) | Fitting equation | Correlation (R2) | P value |
pH | Raw BC | 1 | Y = -0.102*X + 0.991 | 0.999 | 0.001** |
3 | Y = -0.316*X + 3.093 | 0.997 | 0.001** |
5 | Y = -0.499*X + 5.052 | 0.995 | 0.003** |
Cys-BC | 1 | Y = -0.359*X + 2.814 | 0.859 | 0.073* |
3 | Y = -1.379*X + 10.640 | 0.851 | 0.077* |
5 | Y = -1.808*X + 14.300 | 0.747 | 0.136 |
CEC | Raw BC | 1 | Y = -0.021*X + 0.4777 | 0.964 | 0.018* |
3 | Y = -0.066*X + 1.497 | 0.946 | 0.027* |
5 | Y = -0.103*X + 2.533 | 0.942 | 0.030* |
Cys-BC | 1 | Y = -0.035*X + 0.6342 | 0.987 | 0.007** |
3 | Y = -0.137*X + 2.285 | 0.994 | 0.003** |
5 | Y = -0.191*X + 3.505 | 0.997 | 0.001** |
Effects of various biochar treatments on microbial diversity
Table 5 demonstrates that the goods coverage indices for all samples exceed 98%, confirming sufficient sequencing depth. The Shannon diversity index decreases by 8.83% for bacteria and 15.81% for fungi in the T0 − 5 treatment compared to the control (CK) group. Similarly, the Chao1 indices indicate reductions of 8.71% for bacteria and 13.07% for fungi in T0 − 5 relative to CK. The ACE indices also decline by 8.69% for bacteria and 13.14% for fungi. In contrast, the application of Cys-BC significantly mitigates the adverse of cadmium on soil microbial diversity, resulting in increases in the Shannon, Chao1, and ACE indices of 12.81%, 14.48%, and 17.15%, respectively.
Table 5
Alpha indices statistics for bacteria and fungi in soil sample groups
Soil microbe | Treatments | Number of OTUs | Diversity index | Richness estimator | Goods coverage | Phylogenetic diversity |
Shannon | Simpson | Chao1 | ACE |
Bacteria | CK | 2432 | 9.390 | 0.996 | 3078.653 | 3072.826 | 0.982 | 165.653 |
T0 − 5 | 2268 | 8.561 | 0.995 | 2842.903 | 2805.170 | 0.981 | 162.321 |
C1.2−5 | 2607 | 9.566 | 0.996 | 3121.403 | 3096.870 | 0.981 | 175.993 |
T1.2−5 | 2686 | 9.658 | 0.997 | 3254.430 | 3286.196 | 0.982 | 182.432 |
Fungi | CK | 987 | 7.222 | 0.988 | 1134.655 | 1073.938 | 0.997 | 306.808 |
T0 − 5 | 968 | 6.080 | 0.998 | 986.443 | 932.886 | 0.994 | 312.511 |
C1.2−5 | 1022 | 7.346 | 0.982 | 1109.757 | 1121.158 | 0.997 | 311.814 |
T1.2−5 | 1148 | 7.719 | 0.969 | 1276.26 | 1289.536 | 0.997 | 354.215 |
Figure 5a illustrates significant variability in bacterial species abundance between raw BC and CK. The predominant bacterial classes identified include Gammaproteobacteria, Deltaproteobacteria, Planctomycetacia, and Bacilli. Both Gammaproteobacteria and Deltaproteobacteria are instrumental in mitigating Cd toxicity through the degradation of harmful substances and the facilitation of heavy metal transformation (Nawaz et al. 2024). Planctomycetacia plays a vital role in the cycling of nutrients, particularly nitrogen and carbon. Bacilli play a crucial role in plant resistance by producing antioxidants and bioactive compounds.
In contrast, Fig. 5b highlights the significant variability in bacterial species abundance between Cys-BC and CK. The primary bacterial classes observed include Deltaproteobacteria, Nitrospira, and Chloroflexia. Deltaproteobacteria play a crucial role in sulfur and organic matter metabolism, potentially reducing the bioavailability and toxicity of heavy metals such as cadmium. Nitrospira are essential to the nitrogen cycle, particularly in the process of nitrite oxidation, which helps maintain soil nutrient balance. Chloroflexia significantly contributes to photosynthesis and the degradation of organic matter, thereby promoting soil health by facilitating the breakdown of organic materials.
Figure 5c and 5d demonstrate significant variability in fungal species abundance across different treatments. The key fungal classes identified include Mucoromycota-Mucoromycete, Basidiomycota, Tremellomycetes, and Leotiomycetes. Mucoromycota play a crucial role in organic matter decomposition and nutrient cycling, with Mucoromycete potentially forming symbiotic relationships with plants to enhance nutrient uptake and improve tolerance to heavy metals. Basidiomycota includes several important decomposers capable of breaking down complex organic matter, such as lignin and cellulose, thereby enhancing soil quality under Cd stress through enhanced biodegradation (Kumar et al. 2023). Tremellomycetes are often associated with saprophytic and symbiotic relationships that facilitate organic matter decomposition and nutrient release, while they may also degrade organic pollutants by producing extracellular enzymes in response to Cd stress. Finally, Leotiomycetes are primarily involved in humus formation and the maintenance of soil health, supporting plant survival in challenging environments.