3.1. Ambient and elevated CO2 and O3
The ambient CO2 was 399 to 400 ppm during the crop growing period of 2020–2021 and 400–401 ppm during 2021–2022. The elevated CO2 was maintained at 590 to 595 ppm, and 592 to 596 ppm during the crop-growing period of 2020–2021 and 2021–2022 respectively. Similarly, the recorded ambient ozone level was 21–24 ppb, and 22–25 ppb during the crop growing period of 2020–2021 and 2021–2022 respectively. The maintained elevated O3 levels were 61–63 ppb and 60–62 ppb during the crop growing period of 2020–2021 and 2021–2022.
Changes in canopy temperature (Tc) and gas exchange parameters
The canopy temperature (Tc) for the three mustard varieties (PDZM 31, PB, and PM 30) remained relatively consistent across the two seasons (2020–2021 and 2021–2022), irrespective of the elevated conditions. This suggests that Tc may not be as sensitive to changes in O3 and CO2 levels as other physiological parameters.
However, the photosynthetic rate (Pn) of these varieties exhibited significant variability under different treatments as shown in Fig. 2a. In the first season (2020–2021), ambient conditions yielded the highest Pn, particularly in the PM 30 variety, followed by PDZM 31 and PB. This outcome indicates that ambient conditions might be more conducive to photosynthesis in mustard. In contrast, elevated O3 (eO3) conditions led to a substantial reduction in Pn, with the most pronounced effect seen in PDZM 31, demonstrating its sensitivity to ozone stress. Interestingly, elevated CO2 (eCO2) conditions enhanced Pn across all varieties, with PB showing the most significant increase. This enhancement under eCO2 conditions aligns with the expected response of C3 plants to increased CO2 levels, suggesting a potential mitigating effect on the detrimental impacts of ozone. The combined treatment effects of elevated CO2 and O3 (eCO2*eO3) generally resulted in a decreased photosynthetic rate, with the most notable effect observed in PM 30. This indicates that the positive effect of CO2 enrichment on photosynthesis can be offset by the simultaneous increase in ozone levels. The second season (2021–2022) followed a similar trend, with ambient conditions favoring higher Pn. The negative impact of eO3 was evident again, particularly in the PB variety, which showed the most significant decline.
Stomatal conductance (gs) was also significantly impacted under elevated O3 and CO2 conditions and their interaction, as depicted in Fig. 2b. The highest gs was observed under ambient conditions, particularly in PM 30. Elevated O3 caused a marked decline in gs across all varieties, with PB showing the greatest decrease. This decline under eO3 conditions highlights the stress imposed by ozone on the stomatal functioning of mustard plants.
Transpiration rate (E) in the Fig. 2c, mirrored the trends observed in gs, with reductions under ambient conditions and increases under elevated O3. This relationship between E and gs is expected, as stomatal conductance directly influences transpiration. Elevated CO2 conditions resulted in the lowest percentage loss of transpiration rate, especially in the PM 30 variety, suggesting that increased CO2 may lead to more efficient water use in mustard plants.
Water use efficiency (WUE) varied significantly across treatments and seasons, as shown in Fig. 2d. Elevated O3 conditions generally resulted in a loss of WUE, particularly in PDZM 31 and PB, while eCO2 conditions improved WUE, notably in PB and PM 30. However, PDZM 31 showed a negative response to eCO2 in the first season, suggesting variety-specific responses to elevated CO2. The interaction of eCO2 and eO3 often led to a loss in WUE, with the least impact observed in PM 30. This indicates that while CO2 enrichment can improve WUE, its combination with elevated O3 levels may negate this benefit.
These results collectively demonstrate that elevated CO2 and O3 levels have complex and varied impacts on mustard physiology, influenced by both environmental conditions and plant variety. The differential responses of PDZM 31, PB, and PM 30 to these treatments provide valuable insights into the adaptive mechanisms of mustard plants to changing atmospheric conditions.
3.2 Changes in free proline, antioxidant enzymes and MDA activity
In the first season, the exposure of mustard varieties to elevated conditions notably influenced the concentration of free proline compared to ambient conditions. Free proline, often accumulated in plants under stress as a protective measure, saw the most significant reduction in PDZM 31, followed by PM 30 and PB. This decrease under elevated conditions suggests that these varieties might be employing other stress-response mechanisms over proline accumulation. Intriguingly, under elevated O3 (eO3) conditions, all varieties showed an increase in proline concentration, indicating an O3-induced stress response. PB, in particular, displayed a lesser increase, suggesting a possible inherent resilience or different stress response mechanism in this variety.
Across both seasons, antioxidant enzyme activities demonstrated interesting patterns. Catalase (CAT), a key enzyme in detoxifying reactive oxygen species (ROS), showed negligible effects under elevated conditions in all varieties. However, under eO3, there was a marked increase in CAT activity, especially in PDZM 31. This could be indicative of enhanced oxidative stress in plants due to ozone exposure, prompting a stronger antioxidant response. Elevated CO2 (eCO2) conditions, on the other hand, led to a decreased CAT activity compared to eO3, suggesting that eCO2 might be less stressful or that it modulates the oxidative stress differently. The interaction of eCO2 and eO3 showed an increase in CAT activity from ambient levels but a decline compared to eO3 alone, indicating a complex interplay of these gases on oxidative stress responses.
Other antioxidant enzymes like peroxidase (POX), ascorbate peroxidase (APX), and superoxide dismutase (SOD) also exhibited significant increases under elevated conditions, further underscoring the elevated oxidative stress experienced by the plants. PM 30 demonstrated a broader range of antioxidant protection, followed by PDZM 31 and PB. This variation among varieties suggests different capacities or strategies to counteract oxidative stress.
MDA content, an indicator of lipid peroxidation and cellular damage, significantly increased under eO3 conditions across all varieties, reflecting oxidative damage. The similar trends of MDA content under eCO2 and eO3*eCO2 conditions suggest that these conditions either exacerbate or do not mitigate the oxidative stress induced by eO3.
The analysis of these biochemical parameters across two seasons and under various elevated conditions (as depicted in Figs. 3 and 4) highlights the differential responses of mustard varieties to environmental stress. It underscores the complex interplay of external conditions on plant physiology and the varying coping mechanisms among different varieties. This understanding is crucial for predicting plant behavior under changing climatic conditions and for breeding more resilient crop varieties.
3.3 Changes in Growth Parameters
The study's exploration of growth parameters, including Leaf Area Index (LAI) and Relative Growth Rate (RGR), alongside yield parameters like crop index and seed yield, revealed significant variances under different treatment conditions (eO3, eCO2, eO3*eCO2, and control). These findings, detailed in Table 1, provide a nuanced understanding of how elevated environmental conditions impact mustard crop growth and productivity.
Table 1
Growth and seed yield characteristics of mustard varieties grown in elevated ozone, elevated carbon dioxide and their interaction
2020–2021
|
Treatment
|
Stages
|
eO3
|
eCO2
|
eO3*eCO2
|
Ambient
|
SEm
|
CD (0.05)
|
Parameter
|
Genotype
|
PDZM 31
|
PB
|
PM 30
|
PDZM 31
|
PB
|
PM 30
|
PDZM 31
|
PB
|
PM 30
|
PDZM 31
|
PB
|
PM 30
|
LAI
|
45 DAS
|
0.58f
|
1.26e
|
1.33e
|
2.20a
|
2.00b
|
2.10ab
|
1.48d
|
1.58d
|
1.52d
|
1.81c
|
1.86c
|
1.80c
|
0.042
|
0.129
|
60 DAS
|
1.63i
|
2.58fg
|
2.46fg
|
3.04cd
|
3.64ab
|
3.94a
|
1.97h
|
2.91de
|
2.70ef
|
2.35g
|
3.35bc
|
3.35bc
|
0.102
|
0.328
|
75 DAS
|
1.69g
|
2.44e
|
2.16f
|
2.85bc
|
3.05a
|
2.96ab
|
1.83g
|
2.63d
|
2.34e
|
2.35e
|
2.78cd
|
2.64d
|
0.047
|
0.154
|
90 DAS
|
0.79f
|
0.93e
|
0.74f
|
1.48bc
|
1.59a
|
1.52ab
|
0.96e
|
1.02e
|
0.93e
|
1.32d
|
1.38cd
|
1.36d
|
0.033
|
0.107
|
RGR
(g/m2/day)
|
45–60 DAS
|
0.05e
|
0.05de
|
0.05de
|
0.05cd
|
0.06ab
|
0.06a
|
0.05de
|
0.05cd
|
0.06abc
|
0.05de
|
0.06bc
|
0.06ab
|
0.001
|
0.004
|
60–75 DAS
|
0.04bcd
|
0.04d
|
0.04cd
|
0.05a
|
0.04abcd
|
0.04ab
|
0.04abc
|
0.04bcd
|
0.04bcd
|
0.05ab
|
0.04abcd
|
0.04bcd
|
0.002
|
0.007
|
75–90 DAS
|
0.02ef
|
0.02f
|
0.03def
|
0.03abcd
|
0.04ab
|
0.04a
|
0.02ef
|
0.03def
|
0.03bcde
|
0.03cde
|
0.04abc
|
0.04abc
|
0.003
|
0.008
|
Crop index
|
Harvest
|
30.56e
|
35.89de
|
38.87bcd
|
38.84bcd
|
45.89ab
|
46.45a
|
34.83de
|
37.69cde
|
42.17abcd
|
36.61de
|
41.85abcd
|
44.96abc
|
2.307
|
7.524
|
Seed yield (g/m2)
|
Harvest
|
153.00f
|
202.50e
|
247.50cd
|
240.00d
|
344.00a
|
326.75b
|
172.25f
|
249.50cd
|
262.25c
|
203.75e
|
294.75b
|
295.75b
|
6.728
|
19.38
|
2021–2022
|
Treatment
|
Stages
|
eO3
|
eCO2
|
eO3*eCO2
|
Ambient
|
SEm
|
CD (0.05)
|
Parameter
|
Genotype
|
PDZM 31
|
PB
|
PM 30
|
PDZM 31
|
PB
|
PM 30
|
PDZM 31
|
PB
|
PM 30
|
PDZM 31
|
PB
|
PM 30
|
LAI
|
45 DAS
|
0.41h
|
1.20g
|
1.23fg
|
2.34a
|
1.95bc
|
2.07b
|
1.43e
|
1.41ef
|
1.48e
|
1.48e
|
1.71d
|
1.82cd
|
0.062
|
0191
|
60 DAS
|
1.78f
|
2.48de
|
2.38de
|
3.01c
|
3.60a
|
3.60a
|
1.93f
|
3.19bc
|
2.67d
|
2.35e
|
3.19bc
|
3.45ab
|
0.088
|
0.301
|
75 DAS
|
1.61h
|
2.17f
|
2.17f
|
2.76bc
|
2.81b
|
2.93a
|
1.79g
|
2.47d
|
2.29e
|
2.32e
|
2.65c
|
2.66c
|
0.039
|
0.120
|
90 DAS
|
0.69f
|
0.98de
|
0.89e
|
1.32c
|
1.48ab
|
1.57a
|
0.96de
|
1.05d
|
1.01d
|
1.32c
|
1.48ab
|
1.41bc
|
0.036
|
0.123
|
RGR
(g/m2/day)
|
45–60 DAS
|
0.05ef
|
0.05f
|
0.05cd
|
0.05cd
|
0.06bcd
|
0.06a
|
0.05ef
|
0.05ef
|
0.06bc
|
0.05de
|
0.05cd
|
0.06ab
|
0.001
|
0.004
|
60–75 DAS
|
0.03c
|
0.03c
|
0.03d
|
0.04ab
|
0.04a
|
0.04bc
|
0.03c
|
0.04bc
|
0.03cd
|
0.04bc
|
0.04ab
|
0.03c
|
0.001
|
0.004
|
75–90 DAS
|
0.02bcd
|
0.02bcd
|
0.02d
|
0.03abc
|
0.03a
|
0.03ab
|
0.02bcd
|
0.03abc
|
0.02cd
|
0.03abc
|
0.03ab
|
0.03bcd
|
0.002
|
0.006
|
Crop index
|
Harvest
|
36.82d
|
38.46cd
|
38.86bcd
|
48.51a
|
47.76a
|
46.53ab
|
38.63cd
|
38.28cd
|
41.78abcd
|
44.66abc
|
44.24abcd
|
45.10abc
|
2.238
|
7.767
|
Seed yield (g/m2)
|
Harvest
|
148.00h
|
215.00f
|
247.75de
|
241.25e
|
383.25a
|
324.00b
|
174.00g
|
252.25de
|
261.75d
|
212.00f
|
286.75c
|
303.00c
|
6.803
|
19.64
|
eO3- elevated Ozone, eCO2- elevated Carbon dioxide, eO3*eCO2- elevated Ozone & elevated Carbon dioxide interaction
PDZM 31- Pusa Double zero mustard 31, PB- Pusa Bold, PM 30- Pusa Mustard 30
45 DAS- Vegetative stage, 60 DAS- Flowering stage, 75 DAS- Peak flowering and seed filling stage and 90 DAS- Physiological maturity stage
LAI- Leaf area index, RGR- Relative growth rate
|
Exposure to elevated ozone (eO3) and its combination with elevated CO2 (eO3*eCO2) notably reduced LAI and RGR compared to plants exposed only to elevated CO2 (eCO2) and control conditions. This outcome indicates that ozone stress negatively impacts leaf area development and overall plant growth, a trend commonly observed in plants subjected to oxidative stress caused by elevated O3 levels. The impact of ozone on LAI and RGR can be attributed to its interference with photosynthesis and stomatal conductance, leading to reduced carbon assimilation and growth.
The fluctuating pattern of LAI across different developmental stages (45, 60, 75, and 90 DAS) under all treatments points to a dynamic response of mustard plants to changing environmental conditions. The initial decrease in LAI at 45 DAS, followed by an increase at 60 DAS and subsequent decrease at 75 and 90 DAS, suggests that the plant's growth responses are stage-specific. PM 30 showed resilience at the initial growth stage (45 DAS), while Pusa Bold fared better at later stages under eO3 stress.
Interestingly, eCO2 conditions facilitated an increase in LAI, particularly noticeable in PDZM 31 and PM 30 at the earlier growth stages (45 and 60 DAS) and in PB at later stages (75 and 90 DAS). This increase under eCO2 conditions aligns with the expected response of C3 plants, which typically exhibit enhanced growth and leaf expansion due to increased availability of CO2 for photosynthesis.
The interaction treatment (eO3*eCO2) revealed an improvement in LAI, particularly in PB, suggesting a mitigating effect of elevated CO2 on the detrimental impact of ozone. This finding is significant as it indicates that increased CO2 levels might help buffer the negative effects of ozone on plant growth.
RGR, a crucial indicator of overall plant growth efficiency, was highest under eCO2 conditions, especially in PM 30, which was comparable with PB. This suggests that elevated CO2 positively impacts mustard's growth rate, likely due to enhanced photosynthetic activity and carbon assimilation. Conversely, a significant decline in RGR under eO3 conditions across both seasons reaffirms the inhibitory impact of ozone on plant growth.
3.4 Yield Parameters
The crop index and seed yield vary significantly between treatments, varieties and two seasons. In both seasons the crop index, the elevated O3 decreased by 16-17.5%, 13–14% and 13.6% in PDZM 31, PB and PM 30 varieties, respectively, at elevated CO2, the crops responded to increase (6.1–8.6%, 7.9–9.7% and 3.2–3.3% in PDZM 31, PB and PM 30) as compared to ambient condition. The seed yield decreased by 14.3–22%, 12.5–13% and 16.7% in the varieties respectively, at elevated O3 levels as compared to ambient, but increased significantly at elevated CO2 (Table 1). However, the ozone and carbon dioxide interactions in crop index and seed yield showed a reduced percentage decrease in interaction treatments in different varieties of mustard.
Three-way ANOVA (Table SI 2) was carried out to obtain the significant levels of canopy temperature (Tc) and gas exchange parameters (Pn, gs, E and WUE). In terms of the year, a significant difference was observed in photosynthetic rate (Pn) and stomatal conductance (gs) with respect to p ≤ 0.05 significant levels, the treatment effects of Pn, gs, E and WUE were highly significant to the levels p ≤ 0.001 followed by Tc which was slightly significant to the level p ≤ 0.05. The varietal effects concerning photosynthetic rate (Pn) were highly significant to the levels p ≤ 0.001. The treatment*varietal effects under two-way ANOVA were significantly concerning in WUE p ≤ 0.01. With respect to gas exchange parameters and canopy temperature, there was no significant difference observed.
A significant variation was observed from three-way ANOVA for antioxidant enzymes between treatment, years and varieties in Table SI 3. A highly significant difference with respect to year was observed in Catalase and Superoxide dismutase to the levels p ≤ 0.001. In terms of treatment and variety, the antioxidant enzyme activity namely proline, catalase, ascorbate peroxidase, peroxidase, superoxide dismutase and lipid peroxidation were highly significant towards the significant levels p ≤ 0.001. The two-way interaction of Treatment and variety was highly significant in considering APX, POX and SOD to the significant levels of p ≤ 0.001 whereas the effects concerning year*varieties were significant towards SOD to the significant levels p ≤ 0.01 and least significant towards CAT to the levels p ≤ 0.05. The three-way ANOVA was significant to superoxide dismutase (SOD) to P ≤ 0.01 significant levels.
Likewise, under different treatments, varieties of mustard with two seasons (2020–2021 & 2021–2022), (Table SI 4) signified that the LAI at 45 DAS, 75 DAS, RGR at 60–75 and 75–90 DAS showed a very high significance of p ≤ 0.001 and crop index with high significance (p ≤ 0.01) for the two seasons. The treatment effect is very highly significant in LAI, RGR, crop index and seed yield. The varietal effects were highly significant under LAI, RGR at 45–60 DAS and in seed yield (p ≤ 0.001), and high significance was observed in crop index (p ≤ 0.01) and RGR at 60–75 DAS was p ≤ 0.05 significant. A very high performance (p ≤ 0.001) of treatment and varietal interaction was found in LAI at 45, 75 DAS and in seed yield. The significant levels under the variety and year interaction were observed in LAI at 75 DAS (p ≤ 0.001), RGR (p ≤ 0.01) and crop index (p ≤ 0.05)
Differences in canopy temperature, gas exchange parameters, free proline, antioxidant enzyme activity, MDA, growth and seed yield among the different treatments (elevated levels of O3, CO2 and O3*CO2), among different varieties (PDZM 31, PB and PM 30) were analysed by LDA (Linear Discriminant Analysis). Furthermore, LDA exhibited significant effects on mustard varieties (Fig. 5) and treatments (Fig. 6) with the measured parameters. From the results, it was evident that the elevated O3*CO2 interaction altered the effects caused by elevated O3 with respect to canopy temperature, gas exchange parameters, free proline, antioxidant enzyme activity, MDA, growth and seed yield.