All local, national, and international guidelines and legislation were adhered to for the use of plants in this study. The reagents and the equipment used are listed in the Table of Materials.
Experimental layout and material
A pot experiment was arranged at the rooftop of the agricultural engineering department building of Khwaja Fareed University of Engineering and Information Technology, Pakistan (28.42° N, 70.30° E) to investigate the role of vermiculite and Azospirillum brasilense to mitigate cadmium (Cd) toxicity in wheat crop during the winter season of 2023–2024. A completely randomized design (CRD) was used for this experiment, with three replications. An approved wheat variety of “Dilkash-20” was obtained from the Adoptive Research Farm, Rahim Yar Khan, which is suitable for the local region. A pot size 1 x 1 x 1 feet (30.48 cm × 30.48 cm × 30.48 cm) was used, having 12 kg of soil. Seed sterilization was performed before the experiment in the agronomy lab with 70% ethanol. Seeds were planted on 10th November, 2023. 8 wheat seeds were placed in each pot, and 4 plants were maintained at the tillering stage. Vermiculite was obtained from the local market and applied during the addition of the soil at a rate of 150 g/kg of soil. Azospirillum brasilense was obtained from the Government College University Lahore (GCUL), Lahore, Pakistan. Azospirillum brasilense was applied with the seed inoculation in wheat crop by the standard procedure as reported Zaheer et al.17. Experiment was consist of 9 treatments (T0 = Control (No soil amendment, No Cd toxicity), T1 = Cd toxicity 100 mg/kg of soil, T2 = Cd 100 mg/kg of soil + Vermiculite, T3 = Cd 100 mg/kg of soil + A. brasilense, T4 = Cd 100 mg/kg of soil + Vermiculite + A. brasilense, T5 = Cd toxicity 200 mg/kg of soil, T6 = Cd 200 mg/kg of soil + Vermiculite, T7 = Cd 200 mg/kg of soil + A. brasilense, T8 = Cd 200 mg/kg of soil + Vermiculite + A. brasilense). The average monthly temperature in °C and humidity in % at the experimental site are shown in Table 1. Before the start of the experiment, soil analysis was done, and soil properties are shown in Table 2. Irrigation was applied as per requirement at tillering, booting, anthesis, and grain filling stage, and fertilizers were also applied uniformly as per recommendations. Irrigation was applied manually with the use of a water jar and maintained 70% soil moisture in each pot, having no drought stress. Fertilizers were applied as per recommendations of the agriculture department (120-80-60 kg NPK ha-1). Fertilizers were calculated as per the size of each pot. Full doses of P and K and half a dose of N were applied at sowing time, and the remaining half dose of N was applied at the tillering stage. Destructive sampling for root and shoot length measurements (30 days after sowing) was performed on separate pots designated for early-stage harvest; these pots were not used for final yield measurements, so early destructive sampling did not affect later yield data. A schematic timeline of the experiment is as follows: Day 0 (soil preparation, Cd spiking, vermiculite mixing, seed inoculation, sowing), Day 10–30 (root sampling for qPCR), Day 30 (root and shoot length), Day 85 (anthesis: physiological measurements), and Day 120 (harvest: yield parameters). Plant samples were collected at three time points. At 10, 20, and 30 days after sowing (DAS), roots were collected from separate pots. At 30 DAS, separate pots were destructively harvested for root and shoot length. At 85 DAS (anthesis), physiological parameters were measured non-destructively. At 120 DAS, final yield parameters were recorded. Pots used for early destructive sampling (10–30 DAS) were not used for yield measurements.
Cadmium was applied as cadmium chloride (CdCl₂•2.5H₂O, analytical grade, ≥98% purity). The required amounts for target concentrations of 100 mg Cd kg⁻1 soil and 200 mg Cd kg⁻1 soil were calculated based on the dry weight of soil (12 kg per pot). CdCl₂ was dissolved in 500 mL of deionized water per pot and mixed thoroughly with the soil using a mechanical mixer for 15 min to ensure homogeneous distribution. After mixing, the spiked soil was placed in plastic-lined pots and allowed to equilibrate for 14 days at room temperature (25 ± 2 °C) with periodic mixing every 3 days to stabilize Cd speciation and allow aging. Soil moisture was maintained at 70% of field capacity during the equilibration period. After equilibration, representative soil samples (10 g per pot) were collected, air-dried, digested with HNO₃/HClO₄ (4:1, v/v), and analyzed for total Cd concentration using an atomic absorption spectrophotometer to verify actual Cd levels. Measured Cd concentrations were 98.4 ± 3.2 mg kg⁻1 for the 100 mg kg⁻1 treatment and 197.6 ± 5.1 mg kg⁻1 for the 200 mg kg⁻1 treatment, confirming that target levels were achieved. Physicochemical properties and heavy metal content of Vermiculite used in this study are shown in Table 3.
| Month | Temperature (°C) | Humidity (%) |
| November, 2023 | 28 | 56 |
| December, 2023 | 18 | 62 |
| January, 2024 | 17 | 63 |
| February, 2024 | 18 | 32 |
| March, 2024 | 21 | 25 |
| April, 2024 | 24 | 20 |
Table 1: Average temperature and humidity of the experimental site during 2023–2024.
| Parameters | Values (2023-24) |
| Organic matter | 0.56% |
| pH | 7.42 |
| EC | 227 µS/cm |
| T.S.S. | 0.54% |
| Available-P | 5.21 ppm |
| Available-K | 119 ppm |
| Saturation percentage | 34 |
| Soil separates: Sand, Silt, Clay (%) | 36%, 40 %, 24 % |
Table 2: Soil properties of the soil used for the experiment before the sowing of the crop.
| Property | Value | Metal | Concentration (mg kg⁻¹) |
| Grade | Horticultural grade | Cadmium (Cd) | < 0.01 (below detection limit) |
| Particle size | 2-4 mm | Lead (Pb) | 2.1 |
| pH (1:5 in H₂O) | 7.2 | Copper (Cu) | 1.8 |
| Cation exchange capacity (CEC) | 85.4 cmolc kg⁻¹ | Zinc (Zn) | 5.3 |
| Water-holding capacity (WHC) | 450% (w/w) | - | - |
| Bulk density | 0.12 g cm-3 | - | - |
| Application rate | 150 g kg⁻¹ soil (15% w/w) | - | - |
| Amount per pot (12 kg soil) | 1.8 kg | - | - |
Table 3: Physicochemical properties and heavy metal content of Vermiculite used in this study.
Studied traits
All plant growth-related parameters, such as Plant height, spike length, number of spikelets per spike, grains per spike, 1000-grain weight, and grain yield per plant, were recorded manually with the use of a meter rod and a weight balance by using a standard procedure. Biochemical analysis (determination of antioxidant enzyme activities) and physiological analysis (photosynthetic rate, stomatal conductance, chlorophyll content, electrolyte leakage (EL), crop growth rate (CGR), and leaf area index (LAI)) were taken at the anthesis stage after 85 days of sowing. Root and shoot length were noticed after 30 days of wheat plant growth, when plants were removed to maintain 4 plants in the pot after 20 days of sowing. An infrared gas analyzer (Cl-340) and chlorophyll meter (CL-1) were used to observe the photosynthetic rate, stomatal conductance, and chlorophyll contents, respectively. Electrolyte leakage (EL) was measured with the equation reported by Mussell and Staples18 having standard procedure reported by Zaheer et al.14. EL was calculated as: EL (%) = (Initial EC / Final EC) × 100, where Initial EC is the electrical conductivity measured after 24 h of incubation at 25 °C, and Final EC is the conductivity measured after autoclaving the samples at 121 °C for 20 min.
Dissolved organic nitrogen (DOC), dissolved organic carbon (DOC), and soil organic carbon (SOC) were measured with the standard procedure19. LAI and CGR were noticed as per the procedure20. N, P, and K uptake by the wheat plant was noticed, as reported in the procedure described by the previous report21.
Statistical analysis
All the data for this research were collected from three replications, compiled, and averaged. Data was analyzed statistically with the use of Statistic 8.1 and Microsoft Excel 2010. Data were analyzed using one-way ANOVA followed by Fisher's Least Significant Difference (LSD) test at p ≤ 0.05 for mean separation. Different letters showing with the data/results show the statistical difference between each treatment.