Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 07 November 2025

Sec. Plant Abiotic Stress

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1679156

This article is part of the Research TopicExploring Cold Tolerance and Stress in PlantsView all 3 articles

Nitric oxide-mediated modulation of reproductive resilience under cold stress in chickpea

  • 1Department of Botany, Panjab University, Chandigarh, India
  • 2Crop Improvement Division, Indian Institute of Pulses Research, Kanpur, India
  • 3Department of Plant Sciences, Central University of Punjab, Bathinda, India
  • 4Department of Agriculture Biotechnology, Chaudhary Sarwan Kumar (CSK) Himachal Pradesh Agricultural University, Palampur, India
  • 5Department of Biotechnology, National Institute of Plant Genome Research, New, Delhi, India
  • 6The University of Western Australia Institute (UWA) of Agriculture, The University of Western Australia, Perth, WA, Australia
  • 7Sustainable Intensification Innovation Lab, Kansas State University, Manhattan, KS, United States

Chickpeas are particularly sensitive to cold stress during the reproductive phase, which can significantly impair pod set and yield. This study examined the role of sodium nitroprusside (SNP), a nitric oxide (NO) donor, in mitigating cold-induced reproductive damage in cold-tolerant (CT) and cold-sensitive (CS) chickpea genotypes. After 100 days of outdoor growth, plants were subjected to cold stress (15/8°C day/night; 12 h photoperiod) for 21 days in walk-in growth chambers during the reproductive stage of development. Control plants were maintained at 25/15°C day/night temperature. SNP treatment (1 mM) was applied exogenously each time, first two days prior to stress onset and then at seven-day intervals (three applications total). Cold stress significantly lowered endogenous NO levels in leaves, anthers, and ovules, particularly in CS genotypes, thereby leading to reduced pollen viability and germination. SNP treatment restored NO and improved reproductive performance, with stronger responses in the CS than the CT genotype. For instance, pollen germination increased by 57.9% in CS versus 17.6% in CT, and pollen viability increased by 28.0% and 13.1%, respectively. Enhanced anther function resulted in a 157.2% increase in pod set and 62.0% higher seed yield in CS. SNP also improved physiological traits, including a 43.9% increase in cellular viability, 18.6% in stomatal conductance, and 41.9% in chlorophyll content in CS genotypes. Cryoprotectants (proline, trehalose, and sucrose) accumulated in anthers, reinforcing cold resilience, while oxidative stress was simultaneously alleviated through reduced malondialdehyde, hydrogen peroxide, and electrolyte leakage, together with the upregulation of both enzymatic (superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APx), and glutathione reductase (GR)) and non-enzymatic (ascorbic acid (ASC) and reduced glutathione (GSH)) components. Notably, CS genotypes showed more pronounced improvements from SNP application than CT genotypes, particularly in terms of reproductive success and yield-related traits. These findings highlight the potential of NO donors, such as SNP, to enhance cold tolerance in chickpeas, with promising implications for safeguarding productivity under low-temperature stress, especially in sensitive cultivars.

1 Introduction

Chickpeas (Cicer arietinum L.) are a major legume crop cultivated worldwide, particularly in semi-arid regions, where they serve as a key source of dietary protein and support sustainable agricultural systems (Gaur et al., 2015). According to the Food and Agriculture Organization of the United Nations (FAO), 2025, the average global production of chickpeas is 16.51 million tons. India accounts for the largest chickpea area (10.47 Mha) and production (12.27 Mt), whereas Australia recorded the highest yield (3301.7 kg/ha) (FAOSTAT 2025). Other major producers include Pakistan (0.84 Mha, 0.24 Mt), Myanmar (0.32 Mha, 0.41 Mt), Türkiye (0.46 Mha, 0.58 Mt), Ethiopia (0.21 Mha, 0.45 Mt), Canada (0.13 Mha, 0.14 Mt), Mexico (0.075 Mha, 0.14 Mt), the United States of America (USA) (0.15 Mha, 0.21 Mt), and Russia (0.40 Mha, 0.53 However, chickpeas are highly sensitive to cold stress, typically defined as temperatures below 20/10°C (day/night), especially during their reproductive stage. This critical phase often coincides with low temperatures in many production regions (Clarke and Siddique, 2004; Nayyar et al., 2005; Berger et al., 2012). Exposure to cold during reproduction can lead to substantial yield losses owing to flower abortion, reduced pod set, and impaired seed development (Croser et al., 2003; Nayyar et al., 2005; Rani et al., 2020). These yield penalties are primarily attributed to impaired pollen viability, ineffective fertilization, and oxidative damage to floral tissues, which interfere with normal reproductive processes (Kiran et al., 2019; Rani et al., 2021; Kaur et al., 2022).

With climate change increasing the frequency and severity of extreme weather events, improving cold stress tolerance in chickpeas has become an urgent breeding and agronomic priority. One promising strategy involves the application of nitric oxide (NO), a versatile signaling molecule that plays a key role in modulating plant responses to abiotic stressors. Nitric oxide (NO) has been extensively studied for its ability to enhance stress tolerance in various species, including cold stress (Siddiqui et al., 2011; Puyaubert and Baudouin, 2014; Fancy et al., 2017). It acts as a regulator of physiological and molecular responses by modulating antioxidant defenses, maintaining cellular redox balance, and activating the expression of stress-responsive genes (Esim et al., 2014; Zhang et al., 2019; Wang et al., 2021). In cold-stressed tea plants, NO elevates the levels of protective metabolites, including gamma-aminobutyric acid (GABA), proline, and sugars (Wang et al., 2021). In wheat, NO application preserves the integrity of the photosynthetic apparatus and stabilizes the cell membranes under cold stress (Feng et al., 2021). The protective role of NO is multifaceted; it can directly scavenge reactive oxygen species (ROS), enhance the activity of antioxidant enzymes, and mitigate oxidative damage to cellular components, including lipids, proteins, and nucleic acids (Fan et al., 2015). Additionally, NO interacts with other signaling molecules to activate cold-responsive pathways and gene expression, thereby safeguarding cellular structures under stress (Feng et al., 2021; Wu et al., 2022).

NO also plays a vital role in plant reproduction, particularly in regulating pollen tube growth and guidance, which are essential for successful fertilization. This function is thought to be mediated through the cyclic guanosine monophosphate (cGMP) signaling pathway (Prado et al., 2004). Although the role of NO in reproductive development under normal conditions is well established, its involvement in protecting reproductive structures under cold stress remains largely unexplored. Most existing studies have focused on NO-mediated protection in vegetative tissues, leaving a significant gap in our understanding of its specific role in safeguarding fertility and yield during low-temperature stress in reproductive tissues.

To address this gap, the present study investigated the potential of exogenously applied NO to preserve reproductive function in chickpeas under cold stress. We hypothesized that endogenous NO in both leaves and reproductive tissues plays a critical role in maintaining fertility during cold stress. Using cold-tolerant and cold-sensitive chickpea genotypes, we evaluated the efficacy of sodium nitroprusside (SNP), a widely used NO donor, in mitigating cold-induced reproductive damage in plants. Specifically, we assessed the effects of SNP treatment on flower retention, pollen viability, pod set and seed yield. We also examined its impact on oxidative damage in anthers and ovules, which are key determinants of reproductive success. Our findings demonstrate that NO enhances cold resilience by reducing oxidative stress, increasing cryoprotective compounds, and improving the functionality of reproductive tissues. These results suggest that NO application could be a promising strategy for improving cold tolerance in chickpea and other sensitive legume crops, ultimately contributing to more stable yields under adverse environmental conditions. Therefore, this study aimed to (i) determine how endogenous NO levels and related enzymes (NOS and NR) respond in chickpea genotypes under cold stress, (ii) evaluate the efficacy of exogenous SNP application in modulating physiological, biochemical, and reproductive traits, and (iii) assess whether SNP can improve pod set and yield, particularly in cold-sensitive genotypes.

2 Materials and methods

2.1 Preliminary experiment: plant raising and stress imposition

2.1.1 Plant material and growth conditions

A preliminary experiment was conducted to evaluate endogenous NO levels and the activities of nitric oxide synthase (NOS) and nitrate reductase (NR) in the leaves, anthers, and ovules of cold-tolerant (CT, ICC 17258) and cold-sensitive (CS, GPF 2) chickpea genotypes under cold stress. Seeds were surface sterilized with 0.1% mercuric chloride for 2 min, rinsed twice with distilled water, inoculated with Mesorhizobium cicero, and sown in pots containing 8 kg of sandy loam soil (63.4% sand, 24.6% silt, 12% clay; pH 7.1) mixed with sand at a ratio of 3:1. The soil was amended with farmyard manure and tricalcium phosphate (10 mg kg−1) (Awasthi et al., 2014). The soil organic carbon content was 1.4–6.1 g/kg, with nutrient availability of 54, 43, and 158 kg ha−1 for N, P, and K, respectively.

2.1.2 Growth regime and cold stress treatment

Seeds were sown in the first week of November 2022. Plants were raised outdoors for 50 days under ambient conditions (26.9/16.6°C day/night; 1,300–1,500 μmol m−2 s−1 light; 60–70% RH). At reproductive onset (bud to pod formation), plants were exposed to 25/15°C, 12 h photoperiod, in a walk-in chamber. Another set of plants was exposed to cold stress (15/8°C (12 h photoperiod; 500 μmol m−2 s−1 light; 65–70% relative humidity (RH)), for 21 days. For this treatment, the temperature of the chamber was reduced by 2°C per day until the desired low temperature was reached. The control plants were maintained at 25/15°C under identical conditions. Recovery was achieved by raising the temperature stepwise (2°C per day) to 32/22°C for both control and stressed plants until they reached maturity.

2.1.3 Sampling and measurements

After 21 days of cold stress, NO, NOS, and NR activity were assayed in all three organs. Three biological replicates were analyzed for each genotype × treatment combination.

2.1.3.1 Endogenous nitric oxide levels

Fresh leaf tissue (approximately 0.5 g) was homogenized in 3 mL of ice-cold 50 mM acetic acid buffer (pH 3.6) containing 4% zinc diacetate using a mortar and a pestle. The homogenate was centrifuged at 10,000 × g for 15 min at 4°C, and the resulting supernatant was collected. The pellet was re-extracted with 1 mL of the same buffer and centrifuged again under the same conditions. The supernatants were combined, and 0.1 g of charcoal was added to the mixture. The solution was vortexed and filtered to obtain the final extract. Nitrite, a stable oxidation product of NO, was quantified by mixing 1 mL of the filtrate with 1 mL of Griess reagent and incubating the mixture at room temperature for 30 min. The absorbance was measured at 540 nm using a spectrophotometer. The nitric oxide concentration was calculated using a standard curve generated with sodium nitrite (NaNO2), following the protocol described by Zhou et al. (2005).

2.1.3.2 Nitrate reductase and nitric oxide synthase activities

Fresh plant tissue was homogenized in 5 mL of 10 mM phosphate buffer (pH 7.0) containing 4% (w/v) polyvinylpyrrolidone (PVP) and 1 mM ethylenediaminetetraacetic acid (EDTA). The homogenate was centrifuged at 12,000 × g for 15 min at 4°C, and the resulting supernatant was used as the crude enzyme extract.

For protein extraction, samples were homogenized in a buffer containing 100 mM HEPES–KOH (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid buffered with potassium hydroxide) (pH 7.5), 1 mM EDTA, 10% (v/v) glycerol, 5 mM 1,4-dithiothreitol (DTT), 0.5 mM phenylmethyl sulfonyl fluoride (PMSF), 0.1% (v/v) Triton X-100, 1% PVP, and 20 µM flavin adenine dinucleotide (FAD). After centrifugation at 12,000 × g for 20 min at 4°C, the supernatant was collected for spectrophotometric assays of NR and NOS activity at 520 nm and 340 nm, respectively, following the modified protocols of Sun et al. (2014) and Zhao et al. (2009).

To assess NR activity, 250 µL of enzyme extract was mixed with 250 µL of pre-warmed (25°C) assay buffer containing 50 mM HEPES–KOH (pH 7.5), 10 mM MgCl2 (magnesium chloride), 1 mM DTT, 2 mM KNO3 (potassium nitrate), and 200 µM NADH (nicotinamide adenine dinucleotide). The reaction was carried out at 30°C for 30 min and terminated by adding 50 µL of 0.5 M zinc acetate. Nitrite formation was measured by adding 1 mL of 1% sulfanilamide in 3 M HCl (hydrochloric acid) and 1 mL of 0.02% N-(1-naphthyl) ethylenediamine in 0.2 M HCl. For the nitrite assay, 1.5 mL of the filtrate was mixed with an equal volume of Griess reagent containing 1% sulfanilamide and 0.1% N-1-naphthyl ethylenediamine dihydrochloride dissolved in 5% phosphoric acid (H3PO4). The mixture was incubated at room temperature for 30 min to allow the colorimetric reaction to occur. Absorbance was measured at 540 nm, and the nitrite concentration was calculated using a standard curve prepared with sodium nitrite (NaNO2).

NOS activity was measured in a 1 mL reaction mixture containing 100 mM phosphate buffer (pH 7.0), 1 mM L-arginine, 2 mM MgCl2, 0.3 mM CaCl2, 4 µM tetrahydrobiopterin (BH4), 1 µM FAD, 1 µM flavin mononucleotide (FMN), 0.2 mM DTT, 0.2 mM NADPH, and 200 µL of protein extract. The rate of NADPH oxidation was monitored by recording the decrease in absorbance at 340 nm for 5 min. NOS activity was calculated using the extinction coefficient of NADPH (ϵ = 6.22 mM−1 cm−1).

2.2 Optimization of exogenous SNP concentration under low temperatures

A follow-up preliminary experiment was conducted to identify the optimal SNP concentration for alleviating cold stress. Plants of a CT (ICC 17258) and CS (GPF 2) chickpea genotype were grown outdoors for 50 d under ambient conditions, with temperatures ranging from 16.6 to 26.9°C, a light intensity of 1,300–1,500 μmol m−2 s−1, and 60–70% relative humidity. Seed treatment with Rhizobium prior to sowing and pot soil details are provided in Section 2.1. At the start of the reproductive phase, some plants were maintained in a chamber at 25/15°C (12 h photoperiod) to serve as a control group. Another group was subjected to a 21-day cold stress treatment, during which the chamber temperature was gradually decreased by 2°C per day until it reached the desired low temperatures (18/8°C, 15/8°C, and 13/7°C day/night) with a light intensity of 500 μmol m −2 s −1. SNP (0.5, 1, and 1.5 mM) was applied as a foliar spray using 20 mL per plant (with 0.05% Tween-20 as a surfactant) to 50-day-old plants. The first application was administered 2 days prior to cold stress, followed by two additional sprays at 7-day intervals. Following this, the temperature for both groups was slowly raised by 2°C per day until it reached 32/22°C, and the plants were maintained at this temperature until they matured. The pod set, defined as the percentage of flowers that developed into pods, was recorded at maturity. Each treatment was applied in triplicate for both genotypes under control and cold stress conditions, and samples for analysis were collected from these biological replicates.

2.3 In-depth experiments

Based on preliminary optimization (Section 2.2), 1 mM SNP was identified as the most effective concentration, particularly under moderate cold stress (15/8°C), and was therefore used in subsequent detailed experiments.

2.3.1 Experimental design and plant material

Four chickpea genotypes were selected based on their contrasting responses to cold stress in prior screenings: two cold-tolerant (CT: ICC 17258 and ICC 16349) and two cold-sensitive (CS: ICC 15567 and GPF 2). Seeds were obtained from ICRISAT (Hyderabad, India) and PAU (Ludhiana, India). Plants were raised as described in Section 2.1.

2.3.2 Treatments and SNP application

Plants were exposed to two temperature regimes:

● Control: 25/15°C (day/night, 12 h photoperiod), and

● Low-temperature stress (LT): 15/8°C (day/night, 12 h photoperiod).

Within each regime, the plants were divided into two groups: untreated and treated with SNP (1 mM). SNP was applied as a foliar spray with Tween 20 surfactant, at a volume of ~20 mL per plant, two days prior to cold stress induction, followed by two further applications at 7-day intervals.

The resulting treatments were as follows: Control, Control + SNP (1 mM), LT, and LT + SNP (1 mM).

2.3.3 Sampling and replication

Each treatment was performed in triplicate using independent biological sets. Leaves, anthers, and ovules were sampled after 21 d of cold stress for physiological and biochemical assays, whereas yield-related traits (flowering, pod set, and seed yield) were assessed at physiological maturity.

2.3.4 Endogenous nitric oxide

Endogenous NO levels were measured as described in Section 2.1.1.

2.3.5 Physiological and biochemical measurements

2.3.5.1 Electrolyte leakage

Membrane damage was estimated by measuring electrolyte leakage (EL) from the leaf discs. The samples were incubated in deionized water at 25°C for 12 h, and the initial conductivity (C1) was recorded. After heating at 80°C for 10 min, the final conductivity (C2) was measured. EL (%) was calculated as (C1/C2) × 100 (Kaushal et al., 2013).

2.3.5.2 Cellular viability

Cell viability was assessed using the TTC reduction assay based on mitochondrial activity. Fresh tissues were incubated with 0.5% TTC in phosphate buffer (pH 7.4) at 25°C in the dark. Formazan was extracted in ethanol and absorbance read at 530 nm. The results were expressed as absorbance g−1 FW (Steponkus and Lanphear, 1967).

2.3.5.3 Relative leaf water content

Relative leaf water content (RLWC) was used to assess the leaf water status (Barrs and Weatherley, 1962). Fresh leaf samples (500 mg) were immersed in distilled water for 2 h, blotted dry with filter paper, and weighed to determine the turgid weight (TW). The samples were then oven-dried at 110°C for 24 h to obtain the dry weight (DW). RLWC was calculated as:

RLWC (%)=(FWDW)/(TWDW)×100, where FW is the fresh weight.

2.3.5.4 Stomatal conductance

Stomatal conductance (SC) was measured on fully expanded leaves using a portable leaf porometer (Decagon Devices, Pullman, Washington, USA). Values are expressed in millimoles per square meter per second (mmol m−2 s−1).

2.3.5.5 Chlorophyll content

Chlorophyll was extracted from 500 mg of fresh leaf tissue using 80% acetone. The extract was centrifuged at 5,702 × g for 15 min, and the supernatant was collected for analysis. Absorbance was measured at 666 and 653 nm using a spectrophotometer. The total chlorophyll content was calculated using the formula described by Lichtenthaler and Wellburn (1983).

2.3.5.6 Chlorophyll fluorescence

The maximum quantum efficiency of photosystem II (Fv/Fm) was measured to assess the photochemical efficiency of the leaves using a chlorophyll fluorometer (OS1-FL; Opti-Sciences, Hudson, New Hampshire, USA).

2.3.6 Oxidative stress and antioxidants

2.3.6.1 Malondialdehyde concentration

MDA, a marker of lipid peroxidation, was quantified using TBA–MDA adducts. Tissue extracts in TCA were reacted with TBA in TCA, heated at 95°C, cooled, centrifuged, and the absorbance was recorded at 532 nm. MDA was calculated using an extinction coefficient of 155 mM−1 cm−1 (Heath and Packer, 1968).

2.3.6.2 Hydrogen peroxide

H2O2 was extracted in 80% acetone, reacted with titanium reagent and ammonia, centrifuged, and the pellet was dissolved in H2SO4. Absorbance was measured at 410 nm, and H2O2 content was calculated using an extinction coefficient of 0.28 mmol−1 cm−1 (Mukherjee and Choudhuri, 1983).

2.3.7 Antioxidant enzyme activities

● Superoxide dismutase (SOD): Activity was measured by monitoring the inhibition of NBT photoreduction in a riboflavin–methionine system at 560 nm. Results were expressed as units mg−1 protein (Dhindsa and Matowe, 1981).

● Catalase (CAT): Activity was measured by monitoring the decomposition of H2O2 at 410 nm (ϵ = 40 mM−1 cm−1) and expressed as mmol min−1 mg−1 of protein (Teranishi et al., 1974).

● Ascorbate peroxidase (APX): Activity was measured by monitoring ascorbate oxidation at 290 nm (ϵ = 2.8 mM−1 cm−1) and expressed as mmol min−1 mg−1 protein (Nakano and Asada, 1981).

● Glutathione reductase (GR): Activity was measured by following NADPH oxidation at 340 nm and expressed as mmol min−1 mg−1 protein (Mavis and Stellwagen, 1968).

2.3.8 Non-enzymatic antioxidants

● Ascorbic acid (ASC): Extracted in TCA, reacted with DNPH and thiourea, developed with H2SO4, and absorbance was measured at 530 nm. ASC content was quantified from a standard curve and expressed as mg g−1 DW (Mukherjee and Choudhuri, 1983).

● Reduced glutathione (GSH): Extracted in metaphosphoric acid, reacted with DTNB, NADPH, and GR enzyme, and absorbance was measured at 412 nm. The GSH content was determined from a standard curve and expressed as nmol g−1 DW (Griffith, 1980).

2.3.9 Osmolyte and sugar analysis

2.3.9.1 Proline

The proline content was estimated using the ninhydrin method. Fresh tissue was homogenized in 3% sulfosalicylic acid, centrifuged, and reacted with acid ninhydrin. After heating and toluene extraction, the absorbance was measured at 520 nm. The results were expressed as µmol g−1 DW (Bates et al., 1973).

2.3.9.2 Trehalose

Trehalose was quantified using the anthrone-sulfuric acid method. The samples were hydrolyzed in HCl, neutralized, and reacted with anthrone. Absorbance was measured at 620 nm, and the trehalose concentration was determined using a standard curve (Trevelyan et al., 1952).

2.3.9.3 Sucrose

The sucrose content was determined using a resorcinol reagent. Tissue extracts were hydrolyzed, reacted with resorcinol-HCl, and the absorbance was measured at 480 nm. The results are expressed as mg g−1 DW (van Handel, 1968).

2.3.10 Reproductive trait assessment

Pollen viability, germination, and ovule viability were evaluated according to established protocols.

● Pollen viability: Assessed using Alexander’s stain (Alexander, 1969). Viable pollen was stained purple, and non-viable pollen remained green.

● Pollen germination: Evaluated in vitro using a germination medium containing sucrose and boric acid. Germinated pollen was counted microscopically, as described by Brewbaker and Kwack (1963).

● Ovule viability was assessed by fluorescein diacetate (FDA) staining, with viable ovules fluorescing under UV light (Heslop-Harrison and Heslop-Harrison, 1970).

● Stigma receptivity: Determined using the hydrogen peroxide bubbling test, in which active peroxidases in receptive stigmas produce effervescence (Shivanna and Rangaswamy, 1992).

● Pod set: Calculated as the percentage of flowers that developed into mature pods at physiological maturity (Sharma and Nayyar, 2014).

2.4 Statistical analysis

A preliminary study assessed endogenous NO, NOS, and NR activities in two contrasting chickpea genotypes across three organs under control and cold stress conditions (Section 2.1). Each genotype was represented by three pots (each containing two plants) with three biological replicates per treatment (nine plants in total). Data were analyzed using a three-way ANOVA (Genotype × Treatment × Organ) in RStudio (R Core Team, 2023) to evaluate main and interactive effects. As the three-way interaction effects (Genotype × Treatment × Organ) showed inconsistent significance across traits, post hoc Tukey’s HSD test (p< 0.05) was applied within each organ to compare genotype × treatment combinations using the emmeans package (Lenth, 2023), and compact letter displays (CLD) were generated with multcomp (Hothorn et al., 2008) (Supplementary Figure 1). Summary statistics were processed using dplyr (Wickham et al., 2023), and visualizations were created with ggplot2 (Wickham, 2016). ANOVA tables, including degrees of freedom, mean squares, and significance levels, are provided in Supplementary Table 1.

The subsequent experiment (Section 2.2) evaluated the effect of different SNP concentrations at various temperatures on pod set (%). Expanded to include two genotypes, four temperature regimes, and four SNP concentrations (32 treatment combinations). Data were analyzed using three-way ANOVA (Genotype × Temperature × SNP). Post hoc comparisons were performed using Tukey’s test (emmeans), and CLD was generated with the multcomp package. Data handling with dplyr, plots with ggplot2. ANOVA tables including df, sums of squares, mean squares, and significance levels are provided in Supplementary Table 2.

Based on these results, the selected treatment (1 mM SNP) was further evaluated using a completely randomized block design (CRBD) comprising four contrasting genotypes, four treatments, and three organs (leaves, anthers, and ovules) (Section 2.3). Each genotype was tested in three independent replicates. For physiological and biochemical traits, one plant per replicate was analyzed (n = 3), while for yield traits, three plants per replicate were used (n = 9 per genotype per treatment). Data for physiological and biochemical traits were analyzed using a three-way ANOVA (Genotype × Treatment × Organ). As the three-way interaction effects (Genotype × Treatment × Organ) showed inconsistent significance across traits, Tukey’s HSD test (p< 0.05) was applied within each organ to compare genotype × treatment combinations, ensuring consistency across all graphical (Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10) and tabulated representations (Tables 1, 2, 3, 4). For traits not measured across all organs (leaf injury, reproductive, and yield), two-way ANOVA (Genotype × Treatment) was applied, followed by the Tukey test. Data handling and summary statistics were carried out using dplyr, and visualizations were generated with ggplot2. Complete ANOVA outputs are provided in Supplementary Table 3.

Figure 1
Bar chart showing endogenous nitric oxide (NO) levels (nmol g-1 DW) in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was subjected to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents NO concentration, and error bars show the standard error of the mean (n = 3). Lowercase letters above bars indicate statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 1. Endogenous NO levels (nmol/g DW) in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258; CT2: ICC 16349; CS1: ICC 15567; CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n = 3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 2
Bar chart illustrating superoxide dismutase (SOD) activity in three organs, leaves, anthers, and ovules, of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was subjected to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis shows enzyme activity in units per milligram protein, and error bars denote standard error (n = 3). Lowercase letters above bars indicate statistically significant differences among genotype × treatment combinations according to Tukey’s test (p < 0.05) within each organ.

Figure 2. Superoxide dismutase activity in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean + SE (n = 3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 3
Bar chart showing catalase (CAT) enzyme activity in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents catalase activity and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 3. Catalase activity in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n =3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 4
Bar chart showing ascorbate peroxidase (APx) enzyme activity in three organs, leaves, anthers, and ovules, of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents ascorbate activity, and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 4. Ascorbate activity in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n =3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 5
Bar chart showing glutathione reductase (GR) enzyme activity in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents glutathione reductase, and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 5. Glutathione reductase in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n = 3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 6
Bar chart showing ascorbic acid in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents ascorbic acid content and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 6. Ascorbic acid in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n = 3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p <0.05), within each organ.

Figure 7
Bar chart showing reduced glutathione in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents reduced glutathione content, and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 7. Reduced glutathione in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n = 3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 8
Bar chart showing proline content in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents proline content, and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 8. Proline content in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean + SE (n=3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 9
Bar chart showing trehalose content in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents trehalose content, and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 9. Trehalose content in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n = 3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Figure 10
Bar chart showing sucrose content in three organs—leaves, anthers, and ovules—of four chickpea genotypes: CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype was exposed to four treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars are grouped by organ and color-coded by genotype. The y-axis represents sucrose content, and error bars indicate standard error (n = 3). Lowercase letters above bars denote statistically significant differences among genotype × treatment interactions according to Tukey’s test (p < 0.05) within each organ.

Figure 10. Sucrose content in three organs (leaves, anthers, ovules) of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean ± SE (n = 3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05), within each organ.

Table 1
www.frontiersin.org

Table 1. Effect of cold stress on physiological parameters in various organs and contrasting chickpea genotypes under four treatments (control, control+SNP, LT, SNP+LT) and with percentage change.

Table 2
www.frontiersin.org

Table 2. Effect of cold stress on anti-oxidant traits across organs and contrasting Chickpea genotypes under four treatments (control, C+SNP, LT, SNP+LT) and with percentage change.

Table 3
www.frontiersin.org

Table 3. Effect of cold stress on cryoprotectants in various organs and contrasting chickpea genotypes under four treatments (control, control+SNP, LT, SNP+LT) and with percentage change.

Table 4
www.frontiersin.org

Table 4. Effect of cold stress on reproductive and yield traits in contrasting chickpea genotypes under the treatment condition (control, Control+SNP, LT, SNP+LT) and with percentage change.

Principal component analysis (PCA) was performed using the ggplot2, factoMineR, and factoextra packages to explore the relationships among traits under low temperature (LT) and SNP + LT treatments across all three tissues. The results were further validated using heatmaps generated with the pheatmap package (Lê et al., 2008; Wickham, 2016; Kolde, 2019; Kassambara and Mundt, 2020).

3 Results

3.1 Preliminary experiments in cold-stressed plants of cold-tolerant and cold-sensitive genotypes

This experiment evaluated endogenous NO levels and NOS and NR activities in the leaves, anthers, and ovules of CT and CS genotypes after 21 d of cold stress (as described in the Materials and Methods).

3.1.1 Endogenous nitric oxide levels

Under cold stress, the CT genotype exhibited significantly higher NO accumulation across all tissues than the CS genotype (Supplementary Figure 1A).

3.1.2 Nitric oxide synthase activity

Cold stress significantly increased NOS activity in the leaves and anthers of the CT genotype, whereas no significant change was observed in the ovules (Supplementary Figure 1B). Conversely, the CS genotype displayed a marked decline in NOS activity in all organs under cold stress, which correlated with the observed decline in endogenous NO levels.

3.1.3 Nitrate reductase activity

Under cold stress, the CT genotype showed a significant increase in NR activity in the leaves, anthers, and ovules (Supplementary Figure 1C). In contrast, NR activity in the CS genotype declined, particularly in anthers.

3.1.4 Testing the effect of SNP on cold-stressed cold-tolerant and cold-sensitive genotypes

A follow-up preliminary experiment assessed the effect of different SNP concentrations (0.5, 1.0, and 1.5 mM) under three low-temperature regimes (13/7°C, 15/8°C, and 18/8°C) on the pod set percentage in cold-stressed CT and CS genotypes (Supplementary Figure S2). Both CT and CS genotypes generally showed improved pod set when treated with SNP under cold conditions. However, at the most severe temperature regime (13/7°C), the CS genotypes exhibited limited responsiveness, consistent with the results shown in the figure.

3.2 In-depth studies

To further explore the role of NO in cold stress responses, both CT and CS chickpea genotypes were treated with 1 mM SNP under control and cold stress conditions. The following traits were evaluated.

3.2.1 Endogenous nitric oxide levels

Cold stress significantly reduced endogenous NO levels in all tissues examined (Figure 1), with the most pronounced decline observed in the anthers of both CT and CS genotypes. Under non-stress conditions, SNP application markedly increased NO levels in the leaves, anthers, and ovules of both CT and CS plants, with a stronger response observed in the CS genotypes. Under cold stress, SNP similarly elevated NO levels in all organs, with a substantially greater increase in CS genotypes, particularly in the anthers.

3.2.2 Relative leaf water content

Cold stress significantly decreased the RLWC in both CT and CS genotypes, with a more severe reduction in the CS genotypes (Table 5). Under control conditions, SNP slightly increased RLWC in CT genotypes, but no such improvement was observed in CS genotypes. However, under cold stress, SNP application enhanced the RLWC in both CT and CS genotypes, with a notably greater effect in the CS genotypes.

Table 5
www.frontiersin.org

Table 5. Effect of cold stress on leaf injury traits in contrasting chickpea genotypes under 4 treatment conditions (control, control+SNP, LT, SNP+LT) and with percentage change.

3.2.3 Stomatal conductance

Cold stress significantly decreased stomatal conductance, with a greater decline in the CS genotypes (Table 5). Under control conditions, SNP enhanced stomatal conductance in both the CT and CS genotypes. Under cold stress, SNP treatment resulted in a slight increase in stomatal conductance in CT genotypes and a more substantial improvement in CS genotypes.

3.2.4 Chlorophyll content

Leaf chlorophyll content declined under cold stress, with a more pronounced reduction in the CS genotypes (Table 5). Under control conditions, SNP application increased the chlorophyll content in both the CT and CS genotypes. When applied under cold stress, SNP further enhanced the chlorophyll content, particularly in the CS genotypes.

3.2.5 Chlorophyll fluorescence

Cold stress reduced chlorophyll fluorescence in the CT and CS genotypes, with a greater decline observed in the CT genotypes (Table 5). Under cold stress, SNP improved chlorophyll fluorescence in both CT and CS genotypes, particularly in the CS genotypes.

3.2.6 Electrolyte leakage

Cold stress significantly increased EL in the leaves, anthers, and ovules of both CT and CS genotypes (Table 1). Among the tissues, the leaves exhibited the highest increase in EL, particularly in the CS genotype. SNP treatment under cold stress significantly reduced EL in all tissues, with ovules showing the greatest reduction, particularly in the CS genotypes.

3.2.7 Cellular viability

Cold stress significantly decreased cellular viability in the leaves, anthers, and ovules of both CT and CS genotypes, with a more pronounced decline in the CS genotypes (Table 1). Ovules were the most adversely affected, particularly in the CS genotypes. SNP application under cold stress improved cellular viability across all tissues, with CS genotypes showing a more robust recovery. While leaves and anthers showed moderate improvements, ovules exhibited the most significant recovery in cellular viability, particularly in the CS genotypes.

3.2.8 Oxidative stress and antioxidant responses

3.2.8.1 Malondialdehyde

Cold stress significantly increased MDA levels in the leaves, anthers, and ovules of both CT and CS genotypes (Table 1). The increase was more pronounced in the CS genotypes, with leaves exhibiting the highest levels, followed by anthers and ovules. SNP treatment under cold stress significantly reduced MDA accumulation in all organs, with the greatest reduction observed in the CS genotypes. Among the organs, leaves, and anthers responded more effectively to SNP treatment than did the ovules.

3.2.8.2 Hydrogen peroxide

Cold stress markedly increased the H2O2 levels in all the organs (Table 1). Both CT and CS genotypes showed elevated levels, particularly in the anthers and ovules. Under control conditions, SNP had a contrasting effect, slightly increasing H2O2 levels in CT genotypes while decreasing them in CS genotypes. However, under cold stress, SNP consistently reduced H2O2 levels in all organs, with a more substantial reduction observed in CT, particularly in anthers and ovules.

3.2.8.3 Superoxide dismutase activity

Cold stress significantly upregulated SOD activity in both vegetative and reproductive organs, with CS genotypes exhibiting a sharper increase (Figure 2, Table 2). The leaves showed the highest SOD induction. Under cold stress, SNP further enhanced SOD activity, particularly in the CS genotypes. Among the organs, leaves, and ovules showed a more consistent and pronounced.

3.2.8.4 Catalase activity

Catalase activity increased across all organs under cold stress, with a more substantial rise observed in the CS genotypes than in the CT genotypes (Figure 3, Table 2). Anthers and ovules showed greater increases than leaves, particularly in the CS genotypes. SNP application under cold stress further enhanced CAT activity in all organs, with the most pronounced effects observed in the CS genotypes. The strongest enhancement was observed in ovules, followed by anthers and leaves.

3.2.8.5 Ascorbate peroxidase activity

Cold stress increased APX activity in all examined organs in both CT and CS genotypes, with ovules showing the strongest response, particularly in the CT genotypes (Figure 4, Table 2). However, CS genotypes showed a relatively lower induction of APX activity under cold stress conditions. SNP treatment further boosted APX activity in all organs, with the greatest enhancement observed in the anthers. Notably, CS genotypes showed a more marked benefit from the SNP application.

3.2.8.6 Glutathione reductase activity

Cold stress upregulated GR activity in all organs, with the anthers and ovules showing the greatest changes (Figure 5, Table 2), reflecting their increased oxidative stress vulnerability. While the CT genotypes maintained relatively higher baseline GR levels, the CS genotypes responded more strongly to SNP treatment. SNP application under cold stress markedly enhanced GR activity, particularly in the anthers of CS genotypes.

3.2.8.7 Ascorbic acid

Cold stress induced organ-specific changes in ASC levels (Figure 6, Table 2). Leaves showed moderate increases in the CT and CS genotypes, whereas the anthers and ovules of the CS genotypes exhibited notable declines. In contrast, CT genotypes maintained higher Asc levels across all organs. SNP treatment significantly increased Asc content, particularly in the anthers and ovules of the CS genotypes.

3.2.8.8 Reduced glutathione

Under cold stress, GSH levels showed genotype- and organ-specific responses (Figure 7, Table 2). In the CT genotypes, GSH levels increased in all organs. In contrast, CS genotypes exhibited a moderate increase in GSH levels in leaves, but significant reductions in anthers and ovules. SNP application under cold stress enhanced GSH accumulation in all organs, with the most pronounced increases observed in the anthers and ovules of the CS genotypes.

3.2.9 Cryoprotectants

3.2.9.1 Proline

Cold stress triggered an increase in proline levels in the leaves and anthers of CT genotypes (Figure 8, Table 3). In contrast, CS genotypes exhibited a notable reduction in proline content in the reproductive tissues, particularly the ovules. SNP application under cold stress conditions significantly enhanced proline accumulation across all organs, with a more pronounced effect in the CS genotypes. Among the organs assessed, ovules were the most adversely affected by cold stress in the CS genotypes, showing the steepest decline in proline levels. However, they also demonstrated the strongest recovery after SNP treatment.

3.2.9.2 Trehalose

Cold stress increased trehalose content in all organs of CT genotypes but markedly reduced trehalose content in CS genotypes, especially in the anthers and ovules (Figure 9, Table 3). SNP application effectively alleviated these reductions and promoted trehalose accumulation in the CS and CT genotypes, with a more pronounced effect observed in the CS genotype. Anthers benefited the most from SNP treatment in terms of trehalose restoration.

3.2.9.3 Sucrose

Under cold stress, the sucrose content showed genotype- and organ-specific responses (Figure 10; Table 3). In CS genotypes, ovules exhibited the most substantial decline, followed by leaves and anthers. In contrast, the CT genotypes showed moderate increases in sucrose content across all organs under cold stress. SNP application significantly mitigated sucrose depletion in CS genotypes and further enhanced sucrose accumulation in CT genotypes. Among the organs, leaves exhibited the greatest recovery and overall increase in sucrose content after SNP treatment, with anthers and ovules also showing considerable improvements, particularly in the CS genotype.

3.2.10 Reproductive function

Cold stress negatively impacted key reproductive traits in both CT and CS genotypes, with more pronounced effects observed in the CS genotypes (Table 4). Cold stress significantly decreased pollen viability, pollen germination, stigma receptivity, and ovule viability, with CT genotypes exhibiting moderate declines and CS genotypes exhibiting severe impairments.

SNP application significantly alleviated cold-induced reproductive injuries. Both CT and CS genotypes exhibited improved pollen viability and germination following SNP treatment, with greater recovery in the CS genotypes. Stigma receptivity and ovule viability were also enhanced by SNP, with more pronounced benefits in the CS genotypes.

3.2.11 Yield traits

Cold stress caused a significant decline in reproductive output in both CT and CS genotypes, with the CS genotypes showing greater sensitivity (Figure 11; Table 4). The number of pods per plant decreased markedly under cold stress, with the sharpest reduction observed in the CS genotypes. Similarly, the seed weight per plant decreased under cold stress, with a more severe decline in the CS genotype than in the CT genotype. Under control conditions, SNP application modestly improved pod number and seed weight in the CT and CS genotypes. However, under cold stress, SNP markedly mitigated the yield losses. The CS genotypes showed greater relative improvement, particularly in seed weight, which exhibited the most significant recovery compared with untreated cold-stressed plants.

Figure 11
The figure contains two bar graphs showing the effect of different treatments on chickpea genotypes. Panel (A) presents pod number, and panel (B) shows seed weight for four contrasting genotypes labeled CT1 (ICC 17258), CT2 (ICC 16349), CS1 (ICC 15567), and CS2 (GPF-2). Each genotype has four treatment bars: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Bars represent mean values with error bars indicating standard error. Different lowercase letters above the bars denote statistically significant differences among genotype × treatment combinations.

Figure 11. (A) Pod number and (B) Seed weight of four contrasting genotypes (CT1: ICC 17258, CT2: ICC 16349, CS1: ICC 15567, CS2: GPF-2) under different treatments: Control, Control + SNP, Low Temperature (LT), and SNP + LT. Data represent mean + SE (n =3). Different lowercase letters indicate significant differences among genotype *treatment interaction according to Tukey's test (p<0.05).

3.3 Principal component analysis and heat map

3.3.1 Leaves

The PCA of the 24 leaf traits under LT and SNP + LT explained 87.8% of the variation (Supplementary Figure 3A). PC1 (81.3%) separated stress-protective traits (endogenous NO and chlorophyll fluorescence) from damage indicators (MDA and H2O2). PC2 (6.5%) was mainly influenced by trehalose and SOD. Heat maps and clustering (Supplementary Figure 3B) showed that the CS genotypes under LT were grouped with high EL, MDA, and H2O2, which decreased after SNP treatment. CT genotypes clustered with higher protective traits, particularly under SNP + LT.

3.3.2 Anthers

The PCA of the 20 anther traits explained 94.2% of the variation (Supplementary Figure 4A). PC1 (80.5%) distinguished reproductive success traits (pollen viability and stigma receptivity) from stress markers (MDA and H2O2). PC2 (13.7%) reflected the antioxidant enzyme activities (CAT, APX, and GR). Heat maps (Supplementary Figure 4B) confirmed that the CS genotypes under LT clustered with oxidative stress traits, whereas SNP shifted them toward improved reproductive function. The CT genotype showed stronger antioxidant clustering under SNP + LT conditions.

3.3.3 Ovules

The PCA of the 20 ovule traits explained 93% of the variation (Supplementary Figure 5A). PC1 (82.1%) separated protective metabolites (proline and NO) from stress damage (MDA and EL). PC2 (10.9%) was associated with antioxidant enzymes (CAT and SOD). Heat maps (Supplementary Figure 5B) indicated that CS genotypes under LT clustered with high oxidative stress, which SNP reduced, while CT genotypes under SNP + LT clustered with stronger antioxidant responses.

The corresponding loading values for each trait and organ are listed in Supplementary Table 4. Across leaves, anthers, and ovules, PCA consistently showed that SNP reduced oxidative damage in CS genotypes and enhanced protective responses in CT genotypes, underscoring SNP’s role of SNP in mitigating cold-induced reproductive injury.

4 Discussion

Cold stress during reproductive development is a critical constraint in chickpeas, leading to impaired reproductive function and yield loss (Croser et al., 2003; Clarke and Siddique, 2004; Kaur et al., 2022). In this study, we demonstrated that nitric oxide (NO), supplied through sodium nitroprusside (SNP), significantly improved cold tolerance during the reproductive stage. Across multiple organs—leaves, anthers, and ovules—SNP enhanced antioxidant activity, reduced oxidative stress markers, promoted the accumulation of cryoprotectants, and improved reproductive traits. These coordinated responses contribute to higher pollen viability, ovule function, and pod set under low temperatures, particularly in cold-sensitive (CS) genotypes. Our findings align with earlier reports on the protective roles of NO under abiotic stress (Cai et al., 2015; Sami et al., 2021; Kaya et al., 2022) but extend this knowledge by establishing its importance in safeguarding reproductive resilience under cold stress in chickpeas.

4.1 Elucidating the mechanisms of NO underlying cold stress mitigation

4.1.1 Nitric oxide-mediated alleviation of cold stress injury

Cold stress reduced endogenous NO levels in leaves, anthers, and ovules, whereas SNP supplementation restored NO levels, particularly in CS genotypes. Similar enhancements of endogenous NO levels following exogenous SNP treatment have been reported in salt-stressed wheat (Kaya et al., 2022) and mustard (Sami et al., 2021). SNP-treated plants exhibited less injury under cold stress.

Electrolyte leakage (EL), a marker of membrane damage (Bajji et al., 2002), increased significantly in cold-stressed plants but was markedly reduced by SNP, consistent with observations in SNP-treated tomatoes under heat stress (Siddiqui et al., 2017) and water-stressed Cakile maritima (Jday et al., 2016). Cellular viability, reflecting mitochondrial activity (Aslam et al., 2022), improved with SNP treatment, likely due to NO-mediated protective mechanisms (Hajihashemi, 2021; Zhang et al., 2023). Similarly, SNP reduced cold-induced damage to chlorophyll content and fluorescence, which are key indicators of photosynthetic efficiency (Cai et al., 2015), consistent with the findings in mustard (Sami et al., 2021).

Cold-induced membrane damage and photosynthetic inhibition are largely due to the accumulation of ROS (Aslam et al., 2022; Song et al., 2024). Elevated EL, reduced viability, and impaired photosynthesis observed in this study confirm oxidative injury, which was effectively alleviated by SNP treatment. Cold stress also reduced relative leaf water content (RLWC), similar to earlier reports in rice (Dong et al., 2019), due to reduced stomatal conductance (Sun et al., 2022) and impaired hydraulic capacity in the roots (Vernieri et al., 2001, Phaseolus). In contrast, SNP maintained a better water status, consistent with its protective effects in sunflowers (Cechin et al., 2015) and Cakile maritima (Jday et al., 2016).

4.1.2 Nitric oxide and oxidative stress regulation

Cold stress triggered excessive ROS accumulation, as reflected by increased malondialdehyde (MDA), H2O2, and electrolyte leakage across the leaves, anthers, and ovules. These oxidative damages are consistent with earlier findings in chickpeas and other crops under chilling stress (Kumar et al., 2011; Shen et al., 2014; Song et al., 2024). SNP supplementation markedly reduced these effects, highlighting NO’s role of NO in maintaining the cellular redox balance. The alleviation of oxidative stress is strongly associated with enhanced antioxidant activity. SNP treatment increased the activities of superoxide dismutase (SOD), ascorbate peroxidase (APX), and glutathione reductase (GR), and improved the pools of ascorbic acid (ASC) and reduced glutathione (GSH). Similar NO-mediated improvements in enzymatic and non-enzymatic antioxidants have been reported in mustard (Sami et al., 2021), tomato (Zhao et al., 2011), barley (Chen et al., 2010), and wheat (Kaya et al., 2022). By lowering ROS markers and enhancing antioxidant machinery, NO effectively protects membranes, sustains the photosynthetic apparatus, and stabilizes reproductive tissues under cold stress, thereby reinforcing its role as a versatile stress mitigator.

4.1.3 Nitric oxide and cryoprotectant solutes

Cold stress induces a significant accumulation of proline, likely through the upregulation of biosynthetic enzymes and release from feedback inhibition (Raza et al., 2023). Proline functions as a multifunctional cryoprotectant, contributing to osmotic adjustment, reactive oxygen species (ROS) scavenging, and stress resilience (Raza et al., 2023). SNP-treated plants accumulated even higher proline levels under cold stress, indicating a positive role of NO in proline metabolism, consistent with the findings in Camellia sinensis under cold stress (Wang et al., 2021).

Trehalose, a non-reducing disaccharide with known osmoprotective and cryoprotective functions (Kosar et al., 2018), also increased under cold stress and was further enhanced by SNP, similar to the results in heat-stressed wheat (Iqbal et al., 2022). Sucrose levels increased significantly under cold stress, especially in CT genotypes, corroborating observations in wheat (Crespi et al., 1991) and tomato (Li et al., 2024), where sucrose accumulation is linked to the enhanced activity of sucrose-synthesizing enzymes (Crespi et al., 1991). SNP treatment further elevated sucrose levels across all organs, in agreement with prior reports on SNP-treated cucumber seedlings under cold stress (Feng et al., 2021).

Together, the increased accumulation of proline, trehalose, and sucrose under SNP treatment highlights NO as a central regulator of osmotic adjustment and cryoprotection, thereby supporting cold tolerance in chickpeas.

4.2 Organ-specific responses to cold stress and SNP

Among the three organs examined, anthers were the most susceptible to cold stress, showing greater sensitivity than leaves and ovules. This vulnerability directly threatens chickpea reproduction, as reflected by significant reductions in pollen viability, pollen germination, ovule viability, and stigma receptivity, findings consistent with earlier reports (Rani et al., 2021; Sharma and Nayyar, 2014). Cold-induced decline in NO levels in anthers, coupled with oxidative stress, likely disrupts reproductive processes, leading to poor fertilization and reduced pod set.

SNP treatment alleviated these effects, with the strongest improvements observed in anthers, particularly in the CS genotypes. SNP markedly increased endogenous NO levels, reduced EL, improved cellular viability, and enhanced cryoprotectant accumulation, all of which contribute to the protection of reproductive structures (Xie et al., 2022). Although leaves and ovules were less severely affected, they also benefited from SNP treatment, showing reduced oxidative damage. Overall, the results indicate that SNP-mediated NO supply confers protection across all major organs, with critical benefits for reproductive tissues.

4.3 Genotypic differences in response to cold stress and SNP

CS genotypes suffered more severe cold-induced damage than CT genotypes, as shown by reduced NO levels, higher EL, and lower cellular viability, particularly in anthers and ovules. Greater cold sensitivity in plants has previously been linked to weaker antioxidant defense activation, lower carbohydrate reserves (Kumar et al., 2011; Karami-Moalem et al., 2018, chickpea), and stronger inhibition of primary photochemistry (Shen et al., 2014, rice; Song et al., 2024, tobacco).

SNP treatment was particularly effective in CS genotypes, leading to higher NO accumulation and improved physiological performance in CS genotypes. In anthers, this increase was correlated with reduced oxidative damage, improved pollen viability, and better reproductive outcomes. Similar genotype-dependent benefits of SNP have been reported in cold-stressed tomatoes (Zhao et al., 2011) and Cd-stressed barley (Chen et al., 2010).

4.4 Implications for reproductive success and future directions

The most significant outcome of NO application was its impact on reproductive resilience under cold stress. SNP-treated plants maintained higher pollen viability, stigma receptivity, and ovule functionality, ultimately leading to improved pod set. These findings are consistent with those of earlier studies linking NO to reproductive development and fertilization processes (Prado et al., 2004; Zafra et al., 2010). While CT genotypes generally sustained higher pod set under stress, SNP markedly improved reproductive outcomes in CS genotypes, effectively narrowing the performance gap between the two. This demonstrates NO’s role of NO as a broad regulator of reproductive success, although its efficiency may vary with the genetic background. Taken together, our results support a model in which NO enhances cold tolerance through a multifaceted strategy involving ROS detoxification, boosting cryoprotectant levels, stabilizing water status, and safeguarding reproductive organ function. These coordinated processes converge to preserve reproductive capacity and yield stability at low temperatures.

Exogenous SNP application markedly improved cold stress resilience in chickpeas, with the strongest benefits observed in cold-sensitive (CS) genotypes. SNP restored endogenous NO levels in leaves, anthers, and ovules, reducing electrolyte leakage, sustaining cellular viability, and alleviating oxidative damage in the flowers. These effects were linked to enhanced antioxidant activity (SOD, APX, GR) and higher ASC and GSH levels, which supported efficient ROS detoxification.

SNP also promoted the accumulation of proline, trehalose, and sucrose, reinforcing osmotic adjustment and cryoprotection. Reproductive tissues, especially anthers, were most vulnerable to cold stress but showed the greatest recovery under SNP treatment, resulting in improved pollen viability, stigma receptivity, ovule function, and pod set. While CT genotypes maintained better tolerance, SNP narrowed the performance gap by significantly improving reproductive outcomes in the CS lines.

Overall, SNP enhanced cold tolerance through coordinated ROS detoxification, osmolyte accumulation, and reproductive organ protection, thereby supporting chickpea yield stability. Future studies should directly evaluate photosynthetic efficiency, water status, and the molecular regulation of NO-mediated antioxidant responses.

5 Conclusions

This study demonstrates that SNP treatment significantly enhances cold stress resilience in chickpea plants, particularly in cold-sensitive (CS) genotypes. Exogenous SNP application effectively elevated endogenous nitric oxide (NO) levels in anthers, leaves, and ovules, with the strongest effects observed in CS genotypes, which are typically more vulnerable to oxidative damage at low temperatures. Increased NO levels were associated with reduced electrolyte leakage and improved cellular viability, indicating that SNP mitigates membrane damage and oxidative stress. Additionally, SNP treatment boosted cryoprotectant accumulation in anthers, improving their structural integrity and function, which are key factors for maintaining pollen viability and germination under cold conditions. The observed improvements in chlorophyll content and stomatal conductance further suggest that SNP contributes to enhanced photosynthetic efficiency and water use under cold stress conditions. Overall, these findings highlight SNP’s potential of SNPs as a promising tool for improving cold tolerance, reproductive success, and yield stability in chickpea crops facing increasingly variable and challenging climates.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

SKa: Methodology, Writing – review & editing. DP: Formal analysis, Methodology, Software, Writing – original draft. UJ: Formal analysis, Software, Writing – original draft. SKu: Writing – review & editing. KDS: Writing – review & editing. SP: Writing – review & editing. KHMS: Writing – review & editing. PP: Writing – review & editing. HN: Conceptualization, Data curation, Formal analysis, Investigation, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This study was funded by the Department of Biotechnology (DBT), New Delhi, India. DP and SKa gratefully acknowledge the Department of Biotechnology (DBT), New Delhi, for awarding research fellowships for this project. HN acknowledges DBT, New Delhi, for providing funding support for the research.

Acknowledgments

Contribution Number 26-048-J from the Kansas Agricultural Experiment Station is acknowledged.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1679156/full#supplementary-material.

References

Alexander, M. P. (1969). Differential staining of aborted and non-aborted pollen. Stain Technol. 44, 117–122.

PubMed Abstract | Google Scholar

Aslam, M., Fakher, B., Ashraf, M. A., Cheng, Y., Wang, B., and Qin, Y. (2022). Plant low temperature stress: Signaling and response. Agronomy 12, 1–21. doi: 10.3390/agronomy12030702

Crossref Full Text | Google Scholar

Awasthi, R., Kaushal, N., Vadez, V., Turner, N. C., Berger, J., Siddique, K. H., et al. (2014). Individual and combined effects of transient drought and heat stress on carbon assimilation and seed filling in chickpea. Funct. Plant Biol. 41, 1148–1167. doi: 10.1071/FP13340

PubMed Abstract | Crossref Full Text | Google Scholar

Bajji, M., Kinet, J. M., and Lutts, S. (2002). The use of the electrolyte leakage method for assessing cell membrane stability as a water stress tolerance test in durum wheat. Plant Growth Regul. 36, 61–70. doi: 10.1023/A:1014732714549

Crossref Full Text | Google Scholar

Barrs, H. D. and Weatherley, P. E. (1962). A re-examination of the relative turgidity technique for estimating water deficits in leaves. Aust. J. Biol. Sci. 15, 413–428. doi: 10.1071/BI9620413

Crossref Full Text | Google Scholar

Bates, L. S., Waldren, R. P., and Teare, I. D. (1973). Rapid determination of free proline for water-stress studies. Plant Soil 39, 205–207. doi: 10.1007/BF00018060

Crossref Full Text | Google Scholar

Berger, J. D., Kumar, S., Nayyar, H., Street, K. A., Sandhu, J. S., Henzell, J. M., et al. (2012). Temperature-stratified screening of chickpea (Cicer arietinum L.) genetic resource collections reveals very limited reproductive chilling tolerance compared to its annual wild relatives. Field Crops Res. 126, 119–129. doi: 10.1016/j.fcr.2011.09.020

Crossref Full Text | Google Scholar

Brewbaker, J. L. and Kwack, B. H. (1963). The essential role of calcium ion in pollen germination and pollen tube growth. Am. J. Bot. 50, 859–865. doi: 10.1002/j.1537-2197.1963.tb06564.x

Crossref Full Text | Google Scholar

Cai, H., He, M., Ma, K., Huang, Y., and Wang, Y. (2015). Salicylic acid alleviates cold-induced photosynthesis inhibition and oxidative stress in Jasminum sambac. Turk. J. Biol. 39, 241–247. doi: 10.3906/biy-1406-35

Crossref Full Text | Google Scholar

Cechin, I., Cardoso, G. S., Fumis, T., de, F., and Corniani, N. (2015). Nitric oxide reduces oxidative damage induced by water stress in sunflower plants. Bragantia 74, 200–206. doi: 10.1590/1678-4499.353

Crossref Full Text | Google Scholar

Chen, F., Wang, F., Sun, H., Cai, Y., Mao, W., Zhang, G., et al. (2010). Genotype-dependent Effect of exogenous nitric oxide on Cd-induced Changes in antioxidative metabolism, ultrastructure, and photosynthetic performance in Barley Seedlings (Hordeum vulgare). J. Plant Growth Regul. 29, 394–408. doi: 10.1007/s00344-010-9151-2

Crossref Full Text | Google Scholar

Clarke, H. J. and Siddique, K. H. M. (2004). Response of chickpea genotypes to low temperature stress during reproductive development. Field Crops Res. 90, 323–334. doi: 10.1016/j.fcr.2004.04.001

Crossref Full Text | Google Scholar

Crespi, M. D., Zabaleta, E. J., Pontis, H. G., and Salerno, G. L. (1991). Sucrose synthase expression during cold acclimation in wheat. Plant Physiol. 96, 887–891. doi: 10.1104/pp.96.3.887

PubMed Abstract | Crossref Full Text | Google Scholar

Croser, J. S., Clarke, H. J., Siddique, K. H. M., and Khan, T. N. (2003). Low-temperature stress: implications for chickpea (Cicer arietinum L.) improvement. Crit. Rev. Plant Sci. 22, 185–219. doi: 10.1080/713610855

Crossref Full Text | Google Scholar

Dhindsa, R. S. and Matowe, W. (1981). Drought tolerance in two mosses: correlated with enzymatic defence against lipid peroxidation. J. Exp. Bot. 32, 79–91. doi: 10.1093/jxb/32.1.79

Crossref Full Text | Google Scholar

Dong, J., Zhao, J., Zhang, S., Yang, T., Liu, Q., Mao, X., et al. (2019). Physiological and genome-wide gene expression analyses of cold-induced leaf rolling at the seedling stage in rice (Oryza sativa L.). Crop J. 7, 431–443. doi: 10.1016/j.cj.2019.01.003

Crossref Full Text | Google Scholar

Esim, N., Atici, O., and Mutlu, S. (2014). Effects of exogenous nitric oxide in wheat seedlings under chilling stress. Toxicol. Ind. Health 30, 268–274. doi: 10.1177/0748233712457444

PubMed Abstract | Crossref Full Text | Google Scholar

Fan, J., Chen, K., Amombo, E., Hu, Z., Chen, L., and Fu, J. (2015). Physiological and molecular mechanisms of nitric oxide (NO) involved in Bermudagrass response to cold stress. PloS One 10, 1–14. doi: 10.1371/journal.pone.0132991

PubMed Abstract | Crossref Full Text | Google Scholar

Fancy, N. N., Bahlmann, A., and Loake, G. J. (2017). Nitric oxide functions in plant abiotic stress. Plant Cell Environ. 40, 462–472. doi: 10.1111/pce.12707

PubMed Abstract | Crossref Full Text | Google Scholar

Feng, Y., Fu, X., Han, L., Xu, C., Liu, C., Bi, H., et al. (2021). Nitric oxide functions as a downstream signal for melatonin-induced cold tolerance in cucumber seedlings. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.686545

PubMed Abstract | Crossref Full Text | Google Scholar

Food and Agriculture Organization of the United Nations (FAO) (2025). FAOSTAT statistical database. Available online at: http://www.fao.org/faostat/en/data (Accessed October 23, 2025).

Google Scholar

Gaur, P. M., Samineni, S., Tripathi, S., Varshney, R. K., and Gowda, C. L. (2015). Allelic relationships of flowering time genes in chickpea. Euphytica 203, 295–308. doi: 10.1007/s10681-014-1261

Crossref Full Text | Google Scholar

Griffith, O. W. (1980). Determination of glutathione and glutathione disulfide using glutathione reductase and 2-vinylpyridine. Anal. Biochem. 106, 207–212. doi: 10.1016/0003-2697(80)90139-6

PubMed Abstract | Crossref Full Text | Google Scholar

Hajihashemi, S. (2021). Characterization of exogenous nitric oxide effect on Crocus sativus response to different irrigation regimes. J. Plant Growth Regul. 40, 1510–1520.

Google Scholar

Heath, R. L. and Packer, L. (1968). Photoperoxidation in isolated chloroplasts: I. Kinetics and stoichiometry of fatty acid peroxidation. Arch. Biochem. Biophys. 125, 189–198. doi: 10.1016/0003-9861(68)90654-1

PubMed Abstract | Crossref Full Text | Google Scholar

Hothorn, T., Bretz, F., and Westfall, P. (2008). Simultaneous inference in general parametric models. Biom. J. 50, 346–363. doi: 10.1002/bimj.200810425

PubMed Abstract | Crossref Full Text | Google Scholar

Iqbal, N., Sehar, Z., Fatma, M., Umar, S., Sofo, A., and Khan, N. A. (2022). Nitric oxide and abscisic acid mediate heat stress tolerance through regulation of osmolytes and antioxidants to protect photosynthesis and growth in wheat plants. Antioxidants 11, 1–21. doi: 10.3390/antiox11020372

PubMed Abstract | Crossref Full Text | Google Scholar

Jday, A., Rejeb, K. B., Slama, I., Saadallah, K., Bordenave, M., Planchais, S., et al. (2016). Effects of exogenous nitric oxide on growth, proline accumulation and antioxidant capacity in Cakile maritima seedlings subjected to water deficit stress. Funct. Plant Biol. 43, 939–948. doi: 10.1071/FP15363

PubMed Abstract | Crossref Full Text | Google Scholar

Karami-Moalem, S., Maali-Amiri, R., and Kazemi-Shahandashti, S.-S. (2018). Effect of cold stress on oxidative damage and mitochondrial respiratory properties in chickpea. Plant Physiol. Biochem. 122, 31–39. doi: 10.1016/j.plaphy.2017.11.011

PubMed Abstract | Crossref Full Text | Google Scholar

Kassambara, A. and Mundt, F. (2020). Factoextra: extract and visualize the results of multivariate data analyses. R package version 1.0.7. Available online at: https://CRAN.R-project.org/package=factoextra (Accessed January 01, 2020).

Google Scholar

Kaur, A., Gupta, N., Sharma, S., Singh, P., and Singh, S. (2022). Physiological and biochemical characterization of chickpea genotypes for cold tolerance at reproductive stage. S. Afr. J. Bot. 150, 488–499. doi: 10.1016/j.sajb.2022.08.011

Crossref Full Text | Google Scholar

Kaushal, N., Awasthi, R., Gupta, K., Gaur, P., Siddique, K. H., and Nayyar, H. (2013). Heat-stress-induced reproductive failures in chickpea (Cicer arietinum) are associated with impaired sucrose metabolism in leaves and anthers. Funct. Plant Biol. 40, 1334–1349. doi: 10.1071/FP13082

PubMed Abstract | Crossref Full Text | Google Scholar

Kaya, C., Ugurlar, F., Farooq, S., Ashraf, M., AlYemeni, M. N., and Ahmad, P. (2022). Combined application of asparagine and thiourea improves tolerance to lead stress in wheat by modulating AsA-GSH cycle, lead detoxification, and nitrogen metabolism. Plant Physiol. Biochem. 190, 119–132. doi: 10.1016/j.plaphy.2022.08.014

PubMed Abstract | Crossref Full Text | Google Scholar

Kiran, A., Kumar, S., Nayyar, H., and Sharma, K. D. (2019). Low temperature-induced aberrations in male and female reproductive organ development cause flower abortion in chickpea. Plant Cell Environ. 42, 2075–2089. doi: 10.1111/pce.13536

PubMed Abstract | Crossref Full Text | Google Scholar

Kolde, R. (2019). Pheatmap: pretty heatmaps. R package version 1.0.12. Available online at: https://CRAN.R-project.org/package=pheatmap (Accessed January 04, 2019).

Google Scholar

Kosar, F., Akram, N. A., Ashraf, M., Sadiq, M., and Al-Qurainy, F. (2018). Trehalose-induced improvement in growth, photosynthetic characteristics, and levels of some key osmoprotectants in sunflower (Helianthus annuus L.) under drought stress. Pak. J. Bot. 50, 955–961.

Google Scholar

Kumar, S., Malik, J., Thakur, P., Kaistha, S., Sharma, K. D., Upadhyaya, H. D., et al. (2011). Growth and metabolic responses of contrasting chickpea (Cicer arietinum L.) genotypes to chilling stress at the reproductive phase. Acta Physiol. Plant 33, 779–787. doi: 10.1007/s11738-010-0602-y

Crossref Full Text | Google Scholar

Lê, S., Josse, J., and Husson, F. (2008). FactoMineR: an R package for multivariate analysis. J. Stat. Software 25, 1–18. doi: 10.18637/jss.v025.i01

Crossref Full Text | Google Scholar

Lenth, R. V. (2023). emmeans: Estimated Marginal Means, aka Least-Squares Means. R package version 1.9.5. Available online at: https://CRAN.R-project.org/package=emmeans (Accessed October 23, 2025).

Google Scholar

Li, S., Wang, Y., Liu, Y., Liu, C., Xu, W., Lu, Y., et al. (2024). Sucrose synthase gene SUS3 could enhance cold tolerance in tomato. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1324401

PubMed Abstract | Crossref Full Text | Google Scholar

Lichtenthaler, H. K. and Wellburn, A. R. (1983). Determinations of total carotenoids and chlorophylls a and b of leaf extracts in different solvents. Biochem. Soc Trans. 11, 591–592. doi: 10.1042/bst0110591

Crossref Full Text | Google Scholar

Mavis, R. D. and Stellwagen, E. (1968). Purification and subunit structure of glutathione reductase from bakers’ yeast. J. Biol. Chem. 243, 809–814. doi: 10.1016/S0021-9258(19)81737-4

PubMed Abstract | Crossref Full Text | Google Scholar

Mukherjee, S. P. and Choudhuri, M. A. (1983). Implications of water stress-induced changes in the levels of endogenous ascorbic acid and hydrogen peroxide in Vigna seedlings. Physiol. Plant 58, 166–170. doi: 10.1111/j.1399-3054.1983.tb04162.x

Crossref Full Text | Google Scholar

Nakano, Y. and Asada, K. (1981). Hydrogen peroxide is scavenged by ascorbate-specific peroxidase in spinach chloroplasts. Plant Cell Physiol. 22, 867–880. doi: 10.1093/oxfordjournals.pcp.a076232

Crossref Full Text | Google Scholar

Nayyar, H., Bains, T., and Kumar, S. (2005). Low temperature induced floral abortion in chickpea: relationship to abscisic acid and cryoprotectants in reproductive organs. Environ. Exp. Bot. 53, 39–47. doi: 10.1016/j.envexpbot.2004.02.011

Crossref Full Text | Google Scholar

Prado, A. M., Porterfield, D. M., and Feijó, J. A. (2004). Nitric oxide is involved in growth regulation and re-orientation of pollen tubes. Dev 131, 2707–2714. doi: 10.1242/dev.01153

PubMed Abstract | Crossref Full Text | Google Scholar

Puyaubert, J. and Baudouin, E. (2014). New clues for a cold case: nitric oxide response to low temperature. Plant Cell Environ. 37, 2623–2630. doi: 10.1111/pce.12329

PubMed Abstract | Crossref Full Text | Google Scholar

Rani, A., Devi, P., Jha, U. C., Sharma, K. D., Siddique, K. H. M., and Nayyar, H. (2020). Developing climate-resilient chickpea involving physiological and molecular approaches with a focus on temperature and drought stresses. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.01759

PubMed Abstract | Crossref Full Text | Google Scholar

Rani, A., Kiran, A., Sharma, K. D., Prasad, P. V. V., Jha, U. C., Siddique, K. H. M., et al. (2021). Cold tolerance during the reproductive phase in Chickpea (Cicer arietinum L.) is associated with superior cold acclimation ability involving antioxidants and cryoprotective solutes in anthers and ovules. Antioxidants 10, 1–27. doi: 10.3390/antiox10111693

PubMed Abstract | Crossref Full Text | Google Scholar

Raza, A., Charagh, S., Abbas, S., Hassan, M. U., Saeed, F., Haider, S., et al. (2023). Assessment of proline function in higher plants under extreme temperatures. Plant Biol. J. 25, 379–395. doi: 10.1111/plb.13510

PubMed Abstract | Crossref Full Text | Google Scholar

R Core Team (2023). R: A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing). Available online at: https://www.R-project.org (Accessed October 19, 2023).

Google Scholar

Sami, F., Siddiqui, H., and Hayat, S. (2021). Nitric Oxide-mediated enhancement in photosynthetic efficiency, ion uptake and carbohydrate metabolism that boosts overall photosynthetic machinery in mustard plants. J. Plant Growth Regul. 40, 1088–1110. doi: 10.1007/s00344-020-10166-5

Crossref Full Text | Google Scholar

Sharma, K. D. and Nayyar, H. (2014). Cold stress alters transcription in meiotic anthers of cold-tolerant chickpea (Cicer arietinum L.). BMC Res. Notes 7, 1–13. doi: 10.1186/1756-0500-7-717

PubMed Abstract | Crossref Full Text | Google Scholar

Shen, C., Li, D., He, R., Fang, Z., Xia, Y., Gao, J., et al. (2014). Comparative transcriptome analysis of RNA-seq data for cold-tolerant and cold-sensitive rice genotypes under cold stress. J. Plant Biol. 57, 337–348. doi: 10.1007/s12374-014-0183-1

Crossref Full Text | Google Scholar

Shivanna, K. R. and Rangaswamy, N. S. (1992). “In vitro germination methods,” in Pollen biology. Eds. Shivanna, K. R. and Rangaswamy, N. S. (Springer Berlin Heidelberg, Berlin, Heidelberg), 9–22.

Google Scholar

Siddiqui, M. H., Alamri, S. A., Al-Khaishany, M. Y. Y., Al-Qutami, M. A., Ali, H. M., and Khan, M. N. (2017). Sodium nitroprusside and indole acetic acid improve the tolerance of tomato plants to heat stress by protecting against DNA damage. J. Plant Int. 12, 177–186. doi: 10.1080/17429145.2017.1310941

Crossref Full Text | Google Scholar

Siddiqui, M. H., Al-Whaibi, M. H., and Basalah, M. O. (2011). Role of nitric oxide in tolerance of plants to abiotic stress. Protoplasma 248, 447–455. doi: 10.1007/s00709-010-0206-9

PubMed Abstract | Crossref Full Text | Google Scholar

Song, X., Wang, H., Wang, Y., Zeng, Q., and Zheng, X. (2024). Metabolomics combined with physiology and transcriptomics reveal how Nicotiana tabacum leaves respond to cold stress. Plant Physiol. Biochem. 208, 1–17. doi: 10.1016/j.plaphy.2024.108464

PubMed Abstract | Crossref Full Text | Google Scholar

Steponkus, P. L. and Lanphear, F. O. (1967). Refinement of the triphenyl tetrazolium chloride method of determining cold injury. Plant Physiol. 42, 1423–1426. doi: 10.1104/pp.42.10.1423

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, C., Lu, L., Liu, L., Liu, W., Yu, Y., Liu, X., et al. (2014). Nitrate reductase-mediated early nitric oxide burst alleviates oxidative damage induced by aluminium through enhancement of antioxidant defenses in roots of wheat (Triticum aestivum). New Phytol. 201, 1240–1250. doi: 10.1111/nph.12597

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, C., Zhu, L., Cao, L., Qi, H., Liu, H., Zhao, F., et al. (2022). PKS5 confers cold tolerance by controlling stomatal movement and regulating cold-responsive genes in Arabidopsis. Life 12, 1–13. doi: 10.3390/life12101633

PubMed Abstract | Crossref Full Text | Google Scholar

Teranishi, Y., Tanaka, A., Osumi, M., and Fukui, S. (1974). Catalase activities of hydrocarbon-utilizing Candida Yeasts. Agric. Biol. Chem. 38, 1213–1220. doi: 10.1080/00021369.1974.10861301

Crossref Full Text | Google Scholar

Trevelyan, W. E., Forrest, R. S., and Harrison, J. S. (1952). Determination of yeast carbohydrates with the anthrone reagent. Nature 170, 626–627. doi: 10.1038/170626a0

PubMed Abstract | Crossref Full Text | Google Scholar

van Handel, E. (1968). Direct microdetermination of sucrose. Analytical Biochem. 22, 280–283.

PubMed Abstract | Google Scholar

Vernieri, P., Lenzi, A., Figaro, M., Tognoni, F., and Pardossi, A. (2001). How the roots contribute to the ability of Phaseolus vulgaris L. @ to cope with chilling-induced water stress. J. Exp. Bot. 52, 2199–2206. doi: 10.1093/jexbot/52.364.2199

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Yu, Q., Li, Y., Li, J., Chen, J., Liu, Z., et al. (2021). Mechanisms of nitric oxide in the regulation of chilling stress tolerance in Camellia sinensis. Horticulturae 7, 1–13. doi: 10.3390/horticulturae7100410

Crossref Full Text | Google Scholar

Wickham, H. (2016). ggplot2: elegant graphics for data analysis (New York US: Springer-Verlag).

Google Scholar

Wickham, H., François, R., Henry, L., and Müller, K. (2023). dplyr: A grammar of data manipulation. Available online at: https://CRAN.R-project.org/package=dplyr (Accessed October 23, 2025).

Google Scholar

Wu, J., Wang, J., Hui, W., Zhao, F., Wang, P., Su, C., et al. (2022). Physiology of plant responses to water stress and related genes: a review. Forests 13, 1–16. doi: 10.3390/f13020324

Crossref Full Text | Google Scholar

Xie, D. L., Zheng, X.-L., Zhou, C.-Y., Kanwar, M. K., and Zhou, J. (2022). Functions of redox signaling in pollen development and stress response. Antioxidants 11, 1–17. doi: 10.3390/antiox11020287

PubMed Abstract | Crossref Full Text | Google Scholar

Zafra, A., Rodríguez-García, M. I., and Alché, J. D. D. (2010). Cellular localization of ROS and NO in olive reproductive tissues during flower development. BMC Plant Biol. 10, 1–14. doi: 10.1186/1471-2229-10-36

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, X., Liu, L., Chen, B., Qin, Z., Xiao, Y., Zhang, Y., et al. (2019). Progress in understanding the physiological and molecular responses of Populus to salt stress. Int. J. Mol. Sci. 20, 1–17. doi: 10.3390/ijms20061312

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, X., Ma, M., Wu, C., Huang, S., and Danish, S. (2023). Mitigation of heat stress in wheat (Triticum aestivum L.) via regulation of physiological attributes using sodium nitroprusside and gibberellic acid. BMC Plant Biol. 23, 1–13. doi: 10.1186/s12870-023-04321-9

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, M.-G., Chen, L., Zhang, L.-L., and Zhang, W.-H. (2009). Nitric reductase-dependent nitric oxide production is involved in cold acclimation and freezing tolerance in Arabidopsis. Plant Physiol. 151, 755–767. doi: 10.1104/pp.109.140996

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, R., Sheng, J., Lv, S., Zheng, Y., Zhang, J., Yu, M., et al. (2011). Nitric oxide participates in the regulation of LeCBF1 gene expression and improves cold tolerance in harvested tomato fruit. Postharvest Biol. Tech. 62, 121–126. doi: 10.1016/j.postharvbio.2011.05.013

Crossref Full Text | Google Scholar

Zhou, B., Guo, Z., Xing, J., and Huang, B. (2005). Nitric oxide is involved in abscisic acid-induced antioxidant activities in Stylosanthes guianensis. J. Exp. Bot. 56, 3223–3228. doi: 10.1093/jxb/eri319

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: nitric oxide, cold, stress, chickpea, SNP, pollen viability, antioxidant defense

Citation: Kaur S, Padhiar D, Jha UC, Kumar S, Sharma KD, Parida SK, Siddique KHM, Prasad PVV and Nayyar H (2025) Nitric oxide-mediated modulation of reproductive resilience under cold stress in chickpea. Front. Plant Sci. 16:1679156. doi: 10.3389/fpls.2025.1679156

Received: 04 August 2025; Accepted: 16 October 2025;
Published: 07 November 2025.

Edited by:

Zemin Wang, Gansu Agricultural University, China

Reviewed by:

Mehdi Hosseinifarahi, Islamic Azad University, Yasuj, Iran
Márcia Carvalho, University of Trás-os-Montes and Alto Douro, Portugal

Copyright © 2025 Kaur, Padhiar, Jha, Kumar, Sharma, Parida, Siddique, Prasad and Nayyar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Harsh Nayyar, aGFyc2huYXl5YXIuaG5AZ21haWwuY29t; Uday Chand Jha, dTk4MTE5ODFAZ21haWwuY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.