Plant growth at the silking stage
At the silking stage, root dry weight was significantly affected (P ≤ 0.01) by G, Cd, Si and G ⋅ Cd (Fig. 1a). Under Cd0, Si1 significantly declined root dry weight by 15.2% in Zhongke11; under Si0 and Si1, Cd1 stress significantly increased root dry weight by 25.9% and 25.8 in Shengrui999. Aboveground dry weight was significant affected (P ≤ 0.05) by G and Cd (Fig. 1b), under Si0, Cd stress significantly increased aboveground dry weight by 7.9% in Zhongke11. Plant height was significant affected (P ≤ 0.01) by G, Cd, Si, G ⋅ Cd and G ⋅ Cd ⋅ Si (Fig. 1c), under Cd0, Si1 significantly increased plant height by 9.4% in Zhongke11. In Shengrui999, under Si0, Cd1 significantly reduced plant height by 13.5%; Under Cd stress, Si1 significantly increased plant height by 9.4%. Root length was significant affected (P ≤ 0.05) by Si (Fig. 1e), under Cd0, Si1 significantly declined root length by 35.3% in Zhongke11. Root surface area was significant affected (P ≤ 0.05) by G and Cd (Fig. 1f), under Cd0, Si1 significantly declined root surface area by 31.0% in Zhongke11. Root volume was significant affected (P ≤ 0.01) by G and Cd (Fig. 1g), under Si1, Cd1 significantly increased root volume by 61.1% in Zhongke11. Average root diameter was significant affected (P ≤ 0.05) by Cd and G ⋅ Cd (Fig. 1h), under Si0, Cd stress significantly increased average root diameter by 23.0% in Zhongke11.
Plant growth at the maturity stage
At maturity, root dry weight was significant affected (P ≤ 0.01) by genotype (G) and G ⋅ Cd (Table 2). Under Si0, Cd stress significantly declined root dry weight by 20.8% in Zhongke11; but Cd stress significantly increased root dry weight by 38.9% and 57.9% under Si0 and Si1 in Shengrui999, respectively. Aboveground dry weight was significant affected (P ≤ 0.01) by genotype (G) and Si (Table 2). Under Cd stress, Si application resulted in a significant increase the aboveground biomass by 11.9% in Shengrui999. Grain yield was significant affected (P ≤ 0.05) by Cd and Si (Table 2). Under nil Cd, the application of Si resulted in a significant increase the grain yield by 10.5% in Zhongke11; Under Si1, Cd stress significantly declined grain yield by 9.82% in Zhongke11. Under Cd stress, Si application significant increased grain yield by 14.4% in Shengrui999.
Table 2
Dry weight and ear characteristics of maize genotypes (Zhongke11 and Shengrui999) at the maturity stage in response to silicon (Si0, 0 mg kg−1; Si1, 200 mg kg−1) and cadmium (Cd0, 0 mg kg−1; Cd1, 20 mg kg−1) applications.
Genotype | Treatment | Dry weight | Grain Yield (g plant−1) | 100 seed weight (g) | Ear length (cm) | Ear thickness (cm) | Ear rows | Bare top length (mm) |
Root (g plant−1) | Aboveground (g plant−1) |
Zhongke11 | Cd0Si0 | 28.5 a | 205 ab | 99.7 bc | 24.3 ab | 18.4 | 47.6 | 14.7 ab | 16.3 bc |
| Cd1Si0 | 22.6 bc | 200 bc | 96.9 bc | 23.0 c | 17.5 | 47 | 14.7 ab | 25.8 a |
| Cd0Si1 | 25.7 ab | 218 a | 110 a | 24.0 bc | 17.1 | 48.9 | 15.3 a | 7.21 e |
| Cd1Si1 | 22.8 bc | 207 ab | 99.3 bc | 24.7 ab | 18.3 | 47.5 | 13.3 abc | 11.0 de |
Shengrui999 | Cd0Si0 | 13.5 d | 188 cd | 96.8 bc | 24.5 ab | 19.2 | 46.4 | 12.7 bc | 17.4 bc |
| Cd1Si0 | 18.7 c | 185 d | 90.0 c | 23.1 c | 18.4 | 45.6 | 12.0 c | 18.7 b |
| Cd0Si1 | 12.7 d | 200 bc | 107 ab | 25.5 a | 17.1 | 48 | 13.3 abc | 14.4 bcd |
| Cd1Si1 | 20.0 c | 207 ab | 103 ab | 24.3 ab | 18.6 | 48.4 | 13.3 abc | 12.9 cd |
ANOVA | G | ** | ** | ns | ns | ns | ns | ** | ns |
| Cd | ns | ns | * | ** | ns | ns | ns | ** |
| Si | ns | ** | ** | ** | ns | ns | ns | ** |
| G ⋅ Cd | ** | ns | ns | ns | ns | ns | ns | ** |
| G ⋅ Si | ns | ns | ns | ns | ns | ns | ns | ** |
| Cd ⋅ Si | ns | ns | ns | * | ns | ns | ns | * |
| G ⋅ Cd ⋅ Si | ns | ns | ns | ns | ns | ns | ns | ns |
For each parameter across genotypes, mean data (± SE, n=3) with different letters indicate significant difference (P ≤ 0.05). ANOVA results for the main factors (genotype, G; silicon, Si; cadmium Cd) and their interactions (G ⋅ Cd, G ⋅ Si, Cd ⋅ Si and G ⋅ Cd ⋅ Si) are given for each parameter (*, P ≤ 0.05; **, P ≤ 0.01; ns, non-significant). |
Cd, Si and Cd ⋅ Si was significant affected (P ≤ 0.05) 100-seed weight (Table 2). Under Si0, Cd stress significant reduced 100-seed weight by 5.6% and 5.8% in Zhongke11 and Shengrui999, respectively. Under Cd stress, Si application significant increased 100-seed weight by 7.7% and 5.3% in Zhongke11 and Shengrui999, respectively. Bare top length was significant affected (P ≤ 0.05) by Cd, Si, G ⋅ Cd, G ⋅ Si and Cd ⋅ Si (Table 2). In Zhongke11, under Si0, Cd stress significant increased bare top length by 58.1%; Under Cd0, Si application significant reduced bare top length by 55.8%; Under Cd1, Si application significant reduced bare top length by 57.3%. In Shengrui999, under Cd stress, Si application significant reduced bare top length by 30.6% (Table 2) (Fig. S1).
Cd concentration and accumulation
At the silking stage, root Cd concentration was significant affected (P ≤ 0.05) by G, Cd, G ⋅ Cd, G ⋅ Si and G ⋅ Cd ⋅ Si (Fig. 2a); Stem Cd concentration was significant affected (P ≤ 0.05) by G, Si and Cd ⋅ Si (Fig. 2b); Leaf Cd concentration was significant affected (P ≤ 0.05) by G, Cd, Si, G ⋅ Cd, and Cd ⋅ Si (Fig. 2c). Under Cd1, Si significantly reduced stem Cd concentration by 18.3%, but significantly increased leaf Cd concentration by 110% in Zhongke11. Under Cd1, Si significantly reduced root and stem Cd concentration by 15.5% and 19.0%, respectively, but significantly increased leaf Cd concentration by 15.5% in Shengrui999.
At the maturity stage, root Cd concentration was significant affected (P ≤ 0.05) by Cd, Si, G ⋅ Si and G ⋅ Cd ⋅ Si (Fig. 2e); Stem Cd concentration was significant affected (P ≤ 0.05) by G, Cd, Si, G ⋅ Cd and Cd ⋅ Si (Fig. 2f); Leaf Cd concentration was significant affected (P ≤ 0.01) by G, Cd, Si, G ⋅ Cd, G ⋅ Si, Cd ⋅ Si and G ⋅ Cd ⋅ Si (Fig. 2g); Grain Cd concentration was significant affected (P ≤ 0.01) by Cd, Si, and Cd ⋅ Si (Fig. 2h). Under Cd1, Si significantly reduced leaf and grain Cd concentration by 31.6% and 14.4%, respectively, but significantly increased root Cd concentration by 57.1% in Zhongke11. Under Cd1, Si significantly reduced root, stem, leaf and grain Cd concentration by 23.6%, 17.2%, 9.1% and 21.4%, respectively, in Shengrui999. Cd concentration in seeds ranged from 0.05−0.06 µg g−1 when grown in Cd amended soil, which is under the safety threshold for human health in maize (0.1 mg kg−1) (Codex Alimentarius Commission 2014).
At the silking stage, root Cd accumulation was significant affected (P ≤ 0.05) by G, Cd, Si, G ⋅ Cd, G ⋅ Si, Cd ⋅ Si and G ⋅ Cd ⋅ Si (Fig. S2, Table. S1); Stem Cd content was significant affected (P ≤ 0.05) by Cd, Si and Cd ⋅ Si; Leaf Cd content was significant affected (P ≤ 0.01) by G, Cd, Si, G ⋅ Cd, and Cd ⋅ Si. Under Cd1, Si1 significantly reduced stem Cd content by 24.9%, but significantly increased root and leaf Cd content by 37.9% and 108.5%, respectively, in Zhongke11. Under Cd1, Si1 significantly reduced root and stem Cd content by 11.3% and 15.8%, respectively, but significantly increased leaf Cd content by 20.8% in Shengrui999.
At the maturity stage, root Cd accumulation was significant affected (P ≤ 0.01) by G, Cd, Si, G ⋅ Cd, G ⋅ Si, Cd ⋅ Si and G ⋅ Cd ⋅ Si (Fig. S2, Table S1); Stem Cd content was significant affected (P ≤ 0.05) by Cd; Leaf Cd content was significant affected (P ≤ 0.01) by G, Cd, Si, G ⋅ Cd, G ⋅ Si, Cd ⋅ Si, and G ⋅ Cd ⋅ Si; Grain Cd content was significant affected (P ≤ 0.05) by Cd, Si, and Cd ⋅ Si; Under Cd1, Si significantly reduced leaf and grain Cd content by 30.5% and 12.2%, respectively, but significantly increased root Cd content by 88.0% in Zhongke11. Under Cd1, Si significantly reduced root and grain Cd content by 18.6% and 10.2%, respectively, in Shengrui999.
Cd bioconcentration and translocation factors
At the silking stage, Si supply increased the Cd bioconcentration factor (BCF) in leaf in Zhongke11 by 200% (Fig. 3c), but reduced the BCF in root (Fig. 3a), stem (Fig. 3b) and leaf (Fig. 3c) by 35.2%, 37.1% and 18.2%, respectively, in Shengrui999. At the maturity stage, Si application increased root (Fig. 3e) Cd BCF by 69.5%, but reduced leaf (Fig. 3g) Cd BCF by 26.2% in Zhongke11; Si application reduced root (Fig. 3e), stem (Fig. 3f) by 29.8%, leaf (Fig. 3g) and grain (Fig. 3h) Cd BCF 35.1%, 29.8%, 22.8% and 33.3%, respectively, in Shengrui999. Roots had significantly higher Cd BCF than other organs at the silking and maturity stages.
At the silking stage, under Cd stress, Si supply increased the leaf Cd translocation factor (TF) in Zhongke11 and Shengrui999 by 85.2% and 50.0%, respectively (Fig. 4c), but reduced stem Cd TF by 29.4% in Zhongke11 (Fig. 4b), and reduced root Cd TF by 35.2% in Shengrui999 (Fig. 4a). At the maturity stage, Si application increased root Cd TF by 69.5%, reduced stem and leaf Cd TF by 39.3% and 55.9%, respectively, in Zhongke11 (Fig. 4f, g); and reduced root Cd TF by 35.1%, increased leaf Cd TF by 18.9% in Shengrui999 (Fig. 4e, g). Roots had significantly higher Cd TF than other organs at the silking and maturity stages (Fig. 4).
Soil Cd concentration and pH
At the silking and maturity stages, soil Cd concentration was significant affected (P ≤ 0.01) by Cd, G ⋅ Si and G ⋅ Cd ⋅ Si (Fig. 5a, c); Si supply reduced soil Cd concentration in Zhongke11, but lowered in Shengrui999. At the both growth stages, pH was significant affected (P ≤ 0.05) by Si and G ⋅ Cd ⋅ Si (Fig. 5b, d). In Shengrui999, under Cd stress, Si supply increased soil pH at the both growth stages, but soil without Si had lower pH.
Plant Cd tolerance, daily intake and health risk assessment
At the silking stage, under Cd stress, Zhongke11 increased Cd tolerance by 6.0% under Si0 and no effect under Si1, and Shengrui999 increased Cd tolerance by 6.0% (Si0) and 5.6% (Si1) (Fig. 6a). At maturity, under Cd stress, Zhongke11 reduced Cd tolerance index by 5.1% (Si0) and 6.0% (Si1); Shengrui999 increased Cd tolerance index by 1.2% (Si0) and 7.1% (Si1) (Fig. 6c). Shengrui999 showed higher tolerance to Cd stress than Zhongke11 at the both silking and maturity stages (Fig. 6a, c).
The application of Si significantly decreased both daily intake of Cd (DIM) and health risk index (HRI) under Cd stress. The daily intake of Cd by adults ranged from 2.34E−05 (Cd1Si1) to 2.98E−05 (Cd1Si0) in Shengrui999 (Fig. 6b). A similar trend was observed for Cd health risk index. Under Cd stress, the application of Si decreased Cd health risk index by 14.4% (Zhongke11) and 21.4% (Shengrui999) (Fig. 6d).
Correlations between root and ear parameters and Cd accumulation
At the silking stage, the correlation among root parameters (root length, root surface area, root volume and average root diameter), plant organs biomass, Cd concentration, bioaccumulation factor, translocation factor was analyzed (Table S2). Pearson’s correlation analysis showed that root parameters (root length, root surface area, and root volume) had significant negative correlation with pH (P ≤ 0.05). Root volume had significant positive correlation with root, stem and soil Cd concentration, root and stem bioaccumulation and translocation factor (P ≤ 0.05). Average root diameter had significant positive correlation with root, stem and soil Cd concentration, root bioaccumulation and translocation factor (P ≤ 0.05).
At the maturity stage, the correlation among grain yield, ear parameters (100-seed weight, ear length, ear thickness and bare top length), plant organs Cd concentration, bioaccumulation factor and translocation factor were analyzed (Table S3). Pearson’s correlation analysis showed that grain yield had significant negative correlation with root, stem, leaf and grain Cd concentration; root, stem, leaf and grain bioaccumulation factor, root and grain translocation factor, DIM and HRI (P ≤ 0.05); 100-seed weight had significant negative correlation with stem, leaf, grain and soil Cd concentration, stem, leaf and grain bioaccumulation and translocation factor, DIM and HRI (P ≤ 0.05). 100-seed weight and ear thickness had significant positive correlation with pH. Bare top length had significant negative correlation with pH (P ≤ 0.05).
Principal component analysis of growth and physiological traits
At the silking stage, the PCA identified four principal components (PCs) (Table S4). PC1 and PC2 accounted for 54.9% and 16.7% of the variation, respectively. PC1 separated the effects of Cd treatment, and PC2 separated the effects of genotype treatment (Fig. S3a). The Cd concentration (root, stem, leaf and soil), content (root, stem, leaf and soil), bioaccumulation factor (root, stem and leaf), and translocation factor (root and stem) were the key factors in PC1. Root biomass, plant height, total root surface area, total root volume, pH, and leaf translocation factor were the key factors in PC2 (Table S4).
At the maturity stage, PCA identified five PCs (Table S5). PC1 and PC2 accounted for 63.2% and 11.0% of the variation, respectively (Fig. S3b). PC1 separated the effects of Cd treatment. The Cd concentration (stem, grain and soil), content (stem and grain), translocation factor (stem, leaf and grain), daily intake of metal and Cd health risk index were the key factors in PC1. Stem biomass, grain yield and ear rows were the key factors in PC2 (Table S5).