Modulation of Antioxidant Defense System Is Associated with Combined Drought and Heat Stress Tolerance in Citrus

Drought and high temperatures are two major abiotic stress factors that often occur simultaneously in nature, affecting negatively crop performance and yield. Moreover, these environmental challenges induce oxidative stress in plants through the production of reactive oxygen species (ROS). Carrizo citrange and Cleopatra mandarin are two citrus genotypes with contrasting ability to cope with the combination of drought and heat stress. In this work, a direct relationship between an increased antioxidant activity and stress tolerance is reported. According to our results, the ability of Carrizo plants to efficiently coordinate superoxide dismutase (SOD), ascorbate peroxidase (APX), catalase (CAT), and glutathione reductase (GR) activities involved in ROS detoxification along with the maintenance of a favorable GSH/GSSG ratio could be related to their relative tolerance to this stress combination. On the other hand, the increment of SOD activity and the inefficient GR activation along with the lack of CAT and APX activities in Cleopatra plants in response to the combination of drought and heat stress, could contribute to an increased oxidative stress and the higher sensibility of this citrus genotype to this stress combination.

Introduction

Environmental stresses cause large economic losses in agriculture every year, constraining crop yield and production. Owing to the consequences of the climate change, different combinations of abiotic stress conditions are severely impacting on plants in the natural field (Mittler, 2006; Suzuki et al., 2014; Zandalinas et al., 2017a). Although research on plants is traditionally based on the study of the responses to single abiotic factors, further effort has been made over the last years to analyze plant responses to different combined stresses, either abiotic or biotic (Suzuki et al., 2014; Zandalinas et al., 2017a). Particularly, drought and heat are considered one of the most frequent abiotic stress combinations that drastically affect global agricultural systems [International Panel of Climate Change (IPCC, 2014)].

Reactive oxygen species are normally produced as a result of aerobic metabolism. However, metabolic imbalances produced by changes in environmental conditions promote the over-accumulation of ROS (Suzuki et al., 2012). In general, abiotic stresses that limit CO₂ availability due to stomatal closure enhance the accumulation of ROS. Interestingly, while ROS, such as H₂O₂, are considered important signal transduction molecules (Baxter et al., 2014; Mittler, 2016), they are also toxic, causing extensive cellular damage and inhibition of photosynthesis (Choudhury et al., 2016). To prevent damage, ROS accumulation is mitigated by the antioxidant machinery including ROS-scavenging enzymes and increased levels of antioxidants such as AsA and GSH, components of the so-called Halliwell-Asada cycle (Mittler et al., 2004). One of the key enzymes of the antioxidant defense system is the SOD, which constitutes the first level of defense against superoxide radicals. SOD-catalyzed ${\mathrm{O}}_{\mathit{2}}^{\mathrm{•-}}$ dismutation renders H₂O₂ as a reaction product, which in turn is removed by APX and CAT activities (Mittler et al., 2004). APX reduces H₂O₂ using AsA as the electron donor and the balance between GSH and GSSG is critical for maintaining a favorable redox status for the detoxification of H₂O₂. In addition, GR, the rate-limiting enzyme of AsA–GSH cycle, keeps the GSH/GSSG ratio favorable for AsA reduction (Foyer and Noctor, 2005).

Several studies have reported that the ability of plants to balance ROS production and scavenging is associated to a higher tolerance to different environmental stresses (Hernandez et al., 2000; Lin et al., 2004; Arbona et al., 2008; Martinez et al., 2016). The accumulation of high amount of ROS-response transcripts in plants subjected to different combinations of stress factors, reflects the relevance of ROS as an important component of acclimation pathways during combined stresses (reviewed in Suzuki et al., 2014). For example, it has been suggested the key role of cytosolic APX1 for the acclimation of plants to a combination of drought and heat (Koussevitzky et al., 2008). In that work, Arabidopsis mutants deficient in this enzyme (apx1), were found to be highly sensitive to this stress combination. Furthermore, ROS–ABA interactions are very important for plant acclimation to stress combination. In this way, previous reports have shown that mutants impaired in the function of the ABA and ROS-regulated protein phosphates 2C (PP2Cs) (abi-1) were sensitive to the combined impact of drought and heat, as well as salinity and heat (Suzuki et al., 2016; Zandalinas et al., 2016a). Furthermore, several studies have reported that the expression of different ROS-scavenging enzymes and the accumulation of different antioxidants exhibit a unique mechanism of response during stress combination that is different than that found in response to each individual stress (Rizhsky et al., 2002, 2004; Srivastava et al., 2012; Prasch and Sonnewald, 2013; Rivero et al., 2013; Pandey et al., 2015; Jin et al., 2016).

We recently demonstrated the different ability of two citrus genotypes, Carrizo citrange and Cleopatra mandarin, to tolerate drought and heat applied alone or in combination. Therefore, physiological responses in terms of gas exchange parameters and chlorophyll fluorescence, evidenced the higher susceptibility of Cleopatra mandarin to combined drought and heat conditions (Zandalinas et al., 2016b). Moreover, metabolite profiling of leaves of both citrus genotypes in response to combined drought and heat revealed that the accumulation of secondary metabolites with antioxidant function is associated to sensitivity to this stress combination (Zandalinas et al., 2017b). Therefore, the higher sensitivity of Cleopatra plants required a higher accumulation of protective metabolites oriented to mitigate the damaging effects of stress, including flavonols, flavones, and limonoids (Zandalinas et al., 2017b). However, the role of the antioxidant defense involving ROS-scavenging enzymes in the tolerance of citrus plants to combined drought and heat is currently unknown. Previous reports have associated the ability to modulate the antioxidant system with the tolerance of citrus plants to waterlogging (Arbona et al., 2008; Hossain et al., 2009), salinity (Arbona et al., 2003), or WS (Wu et al., 2006). In general, these investigations concluded that coordinated antioxidant activity associated to increased activities of SOD and CAT, along with a modulation of the AsA–GSH cycle, allowed citrus plants to reduce stress-induced oxidative damage.

The aim of the present work was to determine the importance of the modulation of the antioxidant system in citrus tolerance to the combination of drought and high temperatures. To achieve this, oxidative metabolism and related antioxidants were studied in two citrus genotypes (Carrizo citrange and Cleopatra mandarin) with different ability to cope with this combined stresses (Zandalinas et al., 2016b).

Materials and Methods

Plant Material and Growth Conditions

Carrizo citrange (Poncirus trifoliata L. Raf. × Citrus sinensis L. Osb.) and Cleopatra mandarin (Citrus reshni Hort. Ex Tan.) plants were purchased from a commercial nursery (Beniplant S.L., Penyíscola, Spain). One-year-old seedlings of both citrus genotypes were grown in plastic pots filled with perlite and watered three times a week with a half-strength Hoagland solution under greenhouse conditions, with natural photoperiod and day and night temperature averaging 25.0 ± 3.0°C and 18.0 ± 3.0°C, respectively. Then, plants were maintained for 2 weeks in growth chambers to acclimate to a 16-h photoperiod at 25°C and relative moisture at approximately 80%. Temperature and relative moisture were recorded regularly with a portable USB datalogger (OM-EL-WIN-USB, Omega, NJ, United States).

Stress Treatments and Experimental Designs

A 24-h experiment of combined drought and heat was performed with both types of plants (Figure 1). High temperatures (40°C) were firstly imposed for 7 days to a group of plants, maintaining another group at 25°C as control. After imposing the temperature treatment, severe WS conditions were applied by transplanting a group of plants grown at 25 or at 40°C to dry perlite. Therefore, four experimental groups for each citrus genotype were established: well-watered plants at 25°C (CT) and at 40°C (HS) and plants subjected to WS at 25°C (WS) and at 40°C (WS+HS). Leaves with an intermediate position in the canopy were harvested and immediately submerged in liquid N₂.

FIGURE 1

FIGURE 1. Experimental design used to subject Carrizo and Cleopatra plants to drought (WS), heat stress (HS), and a combination of drought and heat stress (WS+HS) with details of period times for each stress treatment. Percentages of affected leaves in Carrizo (CC) and Cleopatra (CM) subjected to WS, HS, and WS+HS are also indicated.

Proline Concentration

Proline analysis was performed as described by Bates et al. (1973) with some modifications. Briefly, 50 mg of ground leaf tissue was extracted in 5 ml of 3% sulfosalicylic acid (Panreac, Barcelona, Spain) by sonication for 30 min. After centrifuging at 4000 × g for 20 min at 4°C, 1 ml of the supernatant was mixed with 1 ml of glacial acetic acid and ninhydrin reagent (Panreac) in a 1:1 (v:v) ratio. The reaction mixture was incubated in a water bath at 100°C for 1 h and subsequently centrifuged at 2000 × g for 5 min at 4°C. Finally, absorbance was read at 520 nm. A standard curve was assayed with pure proline (Sigma-Aldrich, St. Louis, MO, United States).

Leaf Water Status

Relative water content of citrus leaves was calculated using adjacent leaves, which were weighed to obtain a leaf M_f. Leaves were allowed to rehydrate overnight in an opaque beaker filled with water. Therefore, they were reweighed to obtain M_t. Finally, leaves were dried at 80°C for 48 h to obtain M_d. RWC was calculated as [(M_f - M_d) × (M_t - M_d)^-1] × 100 according to Morgan (1984).

Malondialdehyde Concentration

Malondialdehyde content was measured following the procedure of Hodges et al. (1999) with modifications. Ground leaf tissue (0.2 g) were extracted in 2 mL 80% ethanol by sonication for 30 min and, after that, centrifuged at 12000 × g for 10 min. Different aliquots of the supernatant were mixed either with 20% trichloroacetic acid or with a mixture of 20% trichloroacetic acid and 0.5% thiobarbituric acid. Both mixtures were incubated in a water bath at 90°C for 1 h. After cooling samples in ice, homogenates were centrifuged at 2000 × g for 10 min at 4°C. Lastly, the absorbance at 440, 534, and 600 nm of supernatants was read. The MDA concentration in the extracts was calculated as follows:

(1) [(Abs 532_+TBA) - (Abs 600_+TBA) - (Abs 532_-TBA - Abs 600_-TBA)] = A.

(2) [(Abs 440_+TBA - Abs 600_+TBA) × 0.0571] = B

(3) MDA equivalents (nmol ml^-1) = (A - B/157,000) × 10⁶

MDA concentration was expressed as nmol MDA per gram of fresh weight.

Gene Expression

The specific primers used for the amplification of each gene are included in Supplementary Table S1. qRT-PCR analyses were performed in a StepOne Real-Time PCR system (Applied Biosystems, Foster City, CA, United States). The reaction mixture contained 1 μL of cDNA, 5 μL of SYBRGreen (Applied Biosystems) and 1 μM of each gene-specific primer pair in a final volume of 10 μL. The thermal profile used to analyze the relative gene expression consisted of 10 min at 95°C for pre-incubation, followed by 40 cycles of 10 s at 95°C for denaturation, 10 s at 60°C for annealing and 20 s at 72°C for extension. Amplicon specificity of the PCR reaction was evaluated by the presence of a single peak in the dissociation curve after the amplification steps. The expression levels of all genes was normalized against the expression of two endogenous control genes (tubulin and actin) based on previous housekeeping selection for citrus tissues (Mafra et al., 2012) and the relative expression were calculated by using REST (Pfaffl et al., 2002). For all genes studied, the reference sample was the expression value obtained at the non-stressed samples and set at zero.

Antioxidant Enzyme Activities

About 100 mg of frozen ground leaf tissue were extracted in 2 mL of phosphate buffer in a ball mill (MillMix20, Domel, Železniki, Slovenija). After centrifugation 14000 × g at 4°C for 10 min, supernatant was recovered. Different buffers were used for enzyme extractions as follows: for APX, 50 mM phosphate buffer (pH 7.8) supplemented with 1 mM sodium ascorbate and 1 mM EDTA; for SOD, 50 mM phosphate buffer (pH 6.8) with 1.33 mM diethyl-diamino-pentaacetic acid; finally, CAT and GR were extracted in 50 mM phosphate buffer (pH 6.8 and pH 7.5, respectively). The APX activity was assayed following the depletion in absorbance at 290 nm due to AsA consumption. The SOD activity was determined following the ${\mathrm{O}}_{\mathit{2}}^{\mathrm{•-}}$-induced reduction of nitroblue tetrazolium using the xanthine–xanthine oxidase system. CAT was determined using the hydrogen peroxide-dependent reduction of titanium chloride. The GR activity was studied following the increase in absorbance at 412 nm during 2 min as result of the production of the adduct DTNB-GSH after GSSG reduction. The reaction was initiated by adding a suitable aliquot of enzyme extract and the increment in absorbance was recorded during 3 min at 265 nm. Soluble protein content was determined according to Bradford (1976) using BSA as a standard. Enzyme activity was expressed as U mg^-1 protein. Further details on enzyme assays are provided in Hossain et al. (2009).

Ascorbate and Glutathione Levels

Procedures for AsA and GSH determinations are described in Hossain et al. (2009). In short, AsA assay is derived from the reduction of Fe³⁺ to Fe²⁺ in acidic solution by AsA. Fe²⁺ forms a red chelate with bipyridyl that absorbs at 525 nm. DHA was calculated by subtracting AsA from total AsA. The DTNB-GSSG reductase recycling process was used to calculate both total (GSH+GSSG) and GSSG levels.

Statistical Analyses

Data are means of three independent determinations and were subjected to analysis of variance (ANOVA) using a two-way ANOVA with the interaction genotype × stress followed by Tukey post hoc test (P < 0.05) when a significant difference was detected.

Results

Leaf Damage

As shown in Figure 1, 24-h of drought applied individually induced visible leaf damage in both citrus genotypes (30 and 28% of Carrizo and Cleopatra leaves, respectively, were injured). Carrizo plants subjected to continuous HS (40°C) were slightly affected, showing only 5 and 7% of total leaves damaged at 6 days and at the end of the experiment, respectively. On the contrary, after 4 days of heat regime, 26% of Cleopatra leaves were damaged, reaching about 60% at the end of the experiment. Plants subjected to a combination of WS+HS showed the highest percentage of leaf damage in both citrus genotypes. Hence, 39 and 88% of leaves were affected by the combined stresses in Carrizo and Cleopatra, respectively (Figure 1).

Water Status

Leaf RWC of Carrizo and Cleopatra plants subjected to drought, HS and a combination of WS+HS was measured (Table 1). WS+HS conditions similarly decreased leaf RWC in Carrizo and Cleopatra: in Carrizo plants subjected to WS and HS, RWC reached 60 and 75% respect to control values, respectively. In Cleopatra plants, RWC decreased to 60 and 69% (with respect to controls) in response to WS and HS, respectively. Interestingly, stress combination had an additive impact on this parameter, showing the greatest decrease (43 and 39% with respect to control values in Carrizo and Cleopatra, respectively; Table 1).

TABLE 1

TABLE 1. Relative water content (RWC) of Carrizo and Cleopatra leaves subjected to drought (WS), heat (HS), and their combination (WS+HS).

Proline Concentration

Endogenous proline levels were examined in leaves and roots of both citrus genotypes in response to individual and combined stresses (Figure 2). In general, basal proline content of both Cleopatra leaves and roots almost doubled the levels observed in Carrizo. Furthermore, proline concentration in Carrizo leaves significantly increased respect to control values in response to individual stresses. In addition, stress combination induced the highest proline concentration in this genotype. Proline content only increased in response to WS and WS+HS in Cleopatra leaves (Figure 2A). On the other hand, significant increments of proline levels were observed in Carrizo roots subjected to WS (2.2-fold) and especially to WS+HS (3.1-fold), whereas HS did not impact on proline build-up. Finally, proline levels increased similarly in Cleopatra roots in response to WS and WS+HS (about two-fold) and HS caused a reduction of its levels below control values (Figure 2B).

FIGURE 2

FIGURE 2. Proline accumulation in leaves (A) and roots (B) of Carrizo and Cleopatra plants subjected to drought (WS), heat (HS), and a combination of drought and heat stress (WS+HS). Different letters denote statistical significance at p ≤ 0.05. G, genotypes; S, stress treatment; G×S, interaction genotype × stress treatment. ^∗P < 0.05; ^∗∗P < 0.01; ^∗∗∗P < 0.001; ns, no statistical differences.

MDA Concentration

Oxidative damage in terms of MDA concentration in response to drought, HS and the combination of WS+HS was studied in leaves and roots of both citrus genotypes (Figure 3). MDA accumulated in Carrizo leaves in response to WS and more prominently in response to WS+HS. On the contrary, Cleopatra leaves increased MDA content in response to HS and especially during WS+HS (Figure 3A). MDA accumulation pattern in roots was different between both citrus genotypes. Whereas WS induced MDA accumulation only in Cleopatra, HS slightly increased its accumulation in both citrus plants. Strikingly, stress combination resulted in a minor MDA accumulation in Carrizo roots whereas in Cleopatra roots, it resulted in a strong MDA accumulation (Figure 3B).

FIGURE 3

FIGURE 3. Malondialdehyde (MDA) accumulation in leaves (A) and roots (B) of Carrizo and Cleopatra plants subjected to drought (WS), heat (HS), and a combination of drought and heat stress (WS+HS). Different letters denote statistical significance at p ≤ 0.05. G, genotypes; S, stress treatment; G×S, interaction genotype × stress treatment. ^∗P < 0.05; ^∗∗P < 0.01; ^∗∗∗P < 0.001; ns, no statistical differences.

Antioxidant Enzymatic Activity

Under all conditions (control or stress), the SOD activity was significantly higher in Carrizo (five-fold) compared to Cleopatra plants. Imposition of individual and combined stresses had no significant impact on SOD activity in Carrizo leaves, whereas Cleopatra plants showed a two-fold and three-fold increment of this enzymatic activity in response to individual and combined stresses, respectively (Figure 4A). Furthermore, the relative expression of the gene encoding SOD-CuZn in Carrizo was up-regulated under individual stress conditions. In Cleopatra leaves, an accumulation of SOD-CuZn transcripts was observed in response to HS and WS+HS treatments. In addition, SOD-Fe transcripts slightly accumulated in response to HS in Carrizo and in response to WS and WS+HS in Cleopatra (Figure 4B).

FIGURE 4

FIGURE 4. Effects of drought (WS), heat (HS), and a combination of drought and heat stress (WS+HS) on SOD activity (A) and transcript expression (B) in leaves of Carrizo and Cleopatra plants. Different letters denote statistical significance at p ≤ 0.05. G, genotypes; S, stress treatment; G×S, interaction genotype × stress treatment. ^∗P < 0.05; ^∗∗P < 0.01; ^∗∗∗P < 0.001; ns, no statistical differences. Scale for gene expression is log₂ of the mean values after normalization against control plants.

Similar to SOD, CAT activity was more than three-fold higher in Carrizo than in Cleopatra in all conditions studied. In response to individual drought and HS, CAT activity did not change with respect to control values in leaves of both citrus genotypes. Interestingly, under stress combination, CAT activity increased in Carrizo and decreased in Cleopatra compared to control levels (Figure 5A). Analysis of the relative accumulation of CAT transcripts in Carrizo revealed a remarkable up-regulation under individual and especially under combined stress conditions. Contrarily, CAT was down-regulated in Cleopatra leaves, particularly under WS and WS+HS (Figure 5B).

FIGURE 5

FIGURE 5. Effects of drought (WS), heat (HS), and a combination of drought and heat stress (WS+HS) on CAT activity (A) and transcript expression (B) in leaves of Carrizo and Cleopatra plants. Different letters denote statistical significance at p ≤ 0.05. G, genotypes; S, stress treatment; G×S, interaction genotype × stress treatment. ^∗P < 0.05; ^∗∗P < 0.01; ^∗∗∗P < 0.001; ns, no statistical differences. Scale for gene expression is log₂ of the mean values after normalization against control plants.

Ascorbate peroxidase activity significantly increased in response to HS and the combination of WS+HS with respect to control conditions in Carrizo leaves, whereas in Cleopatra a significant increment in APX activity was observed only in response to HS (Figure 6A). Moreover, the relative expression of cytosolic APX was up-regulated under HS and especially under WS and WS+HS in Carrizo, whereas only HS and WS+HS induced the accumulation of APX transcripts in Cleopatra (Figure 6B).

FIGURE 6

FIGURE 6. Effects of drought (WS), heat (HS), and a combination of drought and heat stress (WS+HS) on APX activity (A) and cytosolic APX transcript expression (B) in leaves of Carrizo and Cleopatra plants. Different letters denote statistical significance at p ≤ 0.05. G, genotypes; S, stress treatment; G×S, interaction genotype × stress treatment. ^∗P < 0.05; ^∗∗P < 0.01; ^∗∗∗P < 0.001; ns, no statistical differences. Scale for gene expression is log₂ of the mean values after normalization against control plants.

In Carrizo plants, WS significantly increased the GR activity whereas neither HS nor WS+HS had effect on it. In contrast, in Cleopatra plants, WS and WS+HS increased GR activity and HS did not change this enzymatic activity respect to control levels (Figure 7A). Nevertheless, GR transcript number increased under all stress conditions studied in both genotypes, mainly in Carrizo leaves under stress combination (Figure 7B).

FIGURE 7

FIGURE 7. Effects of drought (WS), heat (HS), and a combination of drought and heat stress (WS+HS) on GR activity (A) and transcript expression (B) in leaves of Carrizo and Cleopatra plants. Different letters denote statistical significance at p ≤ 0.05. G, genotypes; S, stress treatment; G×S, interaction genotype × stress treatment. ^∗P < 0.05; ^∗∗P < 0.01; ^∗∗∗P < 0.001; ns, no statistical differences. Scale for gene expression is log₂ of the mean values after normalization against control plants.

AsA and GSH Pool

Under combined stress, tAsA and AsA levels increased in Carrizo and Cleopatra leaves with respect to control values (Table 2). Moreover, Cleopatra showed a higher tAsA content than Carrizo under combined stress conditions. However, DHA content only increased in Cleopatra leaves in response to stress combination. In addition, no significant alteration in leaf redox AsA/DHA ratio was observed within each citrus genotype (Table 2). Additionally, in response to stress combination, Carrizo and Cleopatra leaves accumulated significant higher levels of tGSH, GSH, and GSSG respect to control values (Table 3). Furthermore, HS induced an accumulation of tGSH, GSH, and GSSG compared to control conditions only in Carrizo leaves. GSH/GSSG ratio increased in Cleopatra leaves upon imposition of WS with respect to control values and higher values in this ratio were found in CT and WS conditions respect to Carrizo values (Table 3).

TABLE 2

TABLE 2. Ascorbate (AsA), total ascorbate (tASA), and dehydroascorbate (DHA) content in Carrizo and Cleopatra leaves subjected to drought (WS), heat (HS), and their combination (WS+HS).

TABLE 3

TABLE 3. Total glutathione (tGSH), reduced glutathione (GSH), and oxidized glutathione (GSSG) content in Carrizo and Cleopatra leaves subjected to drought (WS), heat (HS), and their combination (WS+HS).

Discussion

Abiotic stresses including high temperatures, drought or different combinations of environmental challenges, induce metabolic imbalances that can cause an oxidative stress in plant cells. This effect results in the generation and accumulation of ROS, promoting oxidation of cellular components, hindering metabolic activities and affecting organelle integrity (Suzuki et al., 2012). In citrus plants, it has been proposed that environmental cues such as waterlogging, Cu toxicity, salinity or drought induce oxidative damage (Arbona et al., 2003, 2008; Wu et al., 2006; Hossain et al., 2009; Hippler et al., 2016), highlighting the importance of modulating the antioxidant system efficiently to cope with these abiotic stresses. In the present work, the antioxidant machinery of two citrus genotypes, Carrizo citrange and Cleopatra mandarin, with contrasting ability to tolerate the combination of drought and heat (Zandalinas et al., 2016b) was investigated to correlate differences in the modulation of the antioxidant system with tolerance to this stress combination. In this sense, Cleopatra constitutes a genotype more sensitive than Carrizo to drought combined with heat according to data presented in Figure 1 and also reported in Zandalinas et al. (2016b). Therefore, the percentage of damaged leaves in response to heat or a combination of drought and heat was significantly higher in Cleopatra than in Carrizo (Figure 1), demonstrating the higher ability of Carrizo to deal with stresses that involve high temperatures. Moreover, a positive correlation between stress sensitivity and proline accumulation was observed. Hence, Cleopatra accumulated higher amounts of proline in leaves and roots in response to the combination of drought and heat compared to Carrizo. Additionally, both citrus genotypes accumulated higher levels of this metabolite in response to combined stress factors, a more damaging situation with respect to individual stresses. This result matches other works in which a higher stress pressure exerts a major proline accumulation (Claussen, 2005; Kaur and Asthir, 2015) due to its protective roles, including maintenance of redox balance and radical scavenging, maintenance of protein structure and contribution to reduce cell membrane damage (Shao et al., 2008; Szabados and Savouré, 2010).

In addition, oxidative damage (estimated by MDA accumulation) was also higher in leaves and roots of Cleopatra in response to combined stresses (Figure 3), suggesting that the extent of oxidative damage is directly linked to susceptibility of citrus plants to the combination of WS+HS. The increased SOD and CAT activities of Carrizo in both basal and stress conditions compared to Cleopatra (Figures 4A, 5A) could be related to an active and efficient antioxidant response that might be involved in maintaining a lower MDA content (and oxidative stress, Gill and Tuteja, 2010) especially under the combination of drought and heat, and therefore helping citrus plants to cope with the combined stresses. On the contrary, whereas SOD activity of Cleopatra leaves increased in response to stress imposition, CAT activity did not change in response to individual stresses and even decreased below control levels under WS+HS conditions. In addition to CAT, APX removes H₂O₂ and this reaction has been previously reported to be a crucial process for the tolerance of plants to combined drought and heat (Koussevitzky et al., 2008). In our work, APX activity was significantly induced by HS and WS+HS in Carrizo leaves (Figure 6A), suggesting an efficient H₂O₂ scavenging ability under these stress conditions. However, in Cleopatra, only HS significantly induced an increased APX activity (Figure 6A). However, under high temperatures, this enzyme activity could be insufficient to scavenge the excess of H₂O₂ when CAT activity is not activated (Figure 5A), rendering an increased oxidative damage. Additionally, APX dismutase H₂O₂ using AsA as the electron donor (Foyer and Noctor, 2011). Both citrus genotypes showed increases in leaf AsA and tAsA contents in response to WS+HS, suggesting that the accumulation of AsA could be related with a stronger stress pressure. In addition, Cleopatra showed higher tAsA and DHA levels as well as a lower AsA/DHA ratio with respect to Carrizo during this stress combination (Table 2), which are according to the lower APX activity observed during WS+HS in this citrus genotype (Figure 6A).

Accurate modulation of the GSH cycle is involved in maintaining a favorable GSH/GSSG ratio required for cellular redox regulation. In this way, GR activity could effectively recycle GSH at the expense of NADPH (Foyer and Noctor, 2011). The pattern observed for GR activity in Figure 7 indicates that Carrizo plants, despite increasing tGSH, GSH, and GSSG levels in response to HS and WS+HS (Table 3), preserved the GR activity as well as the GSH/GSSG ratio around control values, probably as a result of a lower incidence of the oxidative damage. In contrast, in Cleopatra leaves, the reduction in GSH/GSSG ratio with respect to control values, especially under WS+HS (Table 3), suggests an impairment of GSH recycling. This result points to a better ROS non-enzymatic detoxification system and to an efficient GSH recycling in Carrizo plants compared to Cleopatra, with no apparent NADPH limitation. In this sense, it has been previously reported that maintenance of a more GSSG status could be a consequence of an enhanced ROS accumulation (Foyer and Noctor, 2011). Our results are in accordance with this statement since MDA specially accumulated in Cleopatra in response to WS+HS (Figure 3). Furthermore, the activation of GR activity observed in Cleopatra under this stress combination (Figure 7A) might be insufficient to keep a proper GSH/GSSH ratio, leading to a lower ability for ROS detoxification (Arbona et al., 2008) and, as a result, to a higher sensitivity of this citrus genotype to WS+HS. These results also demonstrated the previous hypothesis, suggesting that a deficient antioxidant system in Cleopatra plants under the combination of drought and heat would lead to an enhanced activation of secondary metabolites with antioxidant properties including flavonols, flavones and limonoids to supplement the antioxidant deficiency and mitigate the damaging effects of stress (Zandalinas et al., 2017b). However, all these metabolic strategies, including proline accumulation (Figure 2), do not seem to be effective as Cleopatra mandarin suffered important damage under WS+HS conditions.

According to our data, the combination of drought and heat negatively impacted both citrus genotypes (Figure 1) but the effective activation of the antioxidant machinery was associated to the ability to tolerate this stress combination. As a result, in Cleopatra plants in response to WS+HS, the increment of SOD activity (Figure 4A) along with the decline in CAT activity (Figure 5A) and the lack of APX activity increase (Figure 6A), compared to control values could be partially responsible of its increased oxidative damage and sensitivity to the combination of drought and high temperatures. In contrast, the ability of Carrizo plants to efficiently activate antioxidant enzymes involved in ROS detoxification along with preserving a favorable GSH/GSSG ratio would be partially related to genotype tolerance to combined stresses. This work provides physiological basis for directing future genetic programs to improve the antioxidant system of Cleopatra mandarin, a genotype that has been very useful as a rootstock for plants cultivated under conditions of water scarcity. However, its future use can be seriously compromised in a scenario of climatic change due to the high sensibility to combined conditions of heat and drought.

Author Contributions

SZ and DB performed the research. AG-C and VA supervised the project and provided funding. SZ, DB, and AG wrote the manuscript and prepared figures. SZ, DB, VA, and AG-C revised the final version. All authors have read and approved the final version of the manuscript.

Funding

This work was supported by the Ministerio de Economía (MINECO) and Universitat Jaume I through grants n° AGL2016-76574-R and UJI-B2016-23/UJI-B2016-24, respectively. SZ was supported by a predoctoral fellowship from Universitat Jaume I.

Conflict of Interest Statement

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.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00953/full#supplementary-material

Abbreviations

ABA, abscisic acid; APX, ascorbate peroxidase; AsA, ascorbate; CAT, catalase; CT, control; GR, glutathione reductase; GSH, glutathione; GSSG, oxidized glutathione; HS, heat stress; M_d, dry mass; MDA, malondialdehyde; M_f, fresh mass; M_t, turgid mass; PP2Cs, protein phosphatases 2C; REST, Relative Expression Software Tool; ROS, reactive oxygen species; RWC, relative water content; SOD, superoxide dismutase; WS, water stress; WS+HS, water stress and heat stress combination.

References

Arbona, V., Flors, V., Jacas, J., García-Agustín, P., and Gómez-Cadenas, A. (2003). Enzymatic and non-enzymatic antioxidant responses of Carrizo citrange, a salt-sensitive citrus rootstock, to different levels of salinity. Plant Cell Physiol. 44, 388–394. doi: 10.1093/pcp/pcg059

PubMed Abstract | CrossRef Full Text | Google Scholar

Arbona, V., Hossain, Z., López-Climent, M. F., Pérez-Clemente, R. M., and Gómez-Cadenas, A. (2008). Antioxidant enzymatic activity is linked to waterlogging stress tolerance in citrus. Physiol. Plant. 132, 452–466. doi: 10.1111/j.1399-3054.2007.01029.x

PubMed Abstract | 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Baxter, A., Mittler, R., and Suzuki, N. (2014). ROS as key players in plant stress signalling. J. Exp. Bot. 65, 1229–1240. doi: 10.1093/jxb/ert375

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254.

Google Scholar

Choudhury, F. K., Rivero, R. M., Blumwald, E., and Mittler, R. (2016). Reactive oxygen species, abiotic stress and stress combination. Plant J. doi: 10.1111/tpj.13299 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Claussen, W. (2005). Proline as a measure of stress in tomato plants. Plant Sci. 168, 241–248. doi: 10.1016/j.plantsci.2004.07.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Foyer, C. H., and Noctor, G. (2005). Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell 17, 1866–1875. doi: 10.1105/tpc.105.033589

PubMed Abstract | CrossRef Full Text | Google Scholar

Foyer, C. H., and Noctor, G. (2011). Ascorbate and glutathione: the heart of the redox hub. Plant Physiol. 155, 2–18. doi: 10.1104/pp.110.167569

PubMed Abstract | CrossRef Full Text | Google Scholar

Gill, S. S., and Tuteja, N. (2010). Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48, 909–930. doi: 10.1016/j.plaphy.2010.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez, J. A., Jimenez, A., Mullineaux, P., and Sevilia, F. (2000). Tolerance of pea (Pisum sativum L.) to long-term salt stress is associated with induction of antioxidant defences. Plant Cell Environ. 23, 853–862. doi: 10.1046/j.1365-3040.2000.00602.x

CrossRef Full Text | Google Scholar

Hippler, F. W. R., Cipriano, D. O., Boaretto, R. M., Quaggio, J. A., Gaziola, S. A., Azevedo, R. A., et al. (2016). Citrus rootstocks regulate the nutritional status and antioxidant system of trees under copper stress. Environ. Exp. Bot. 130, 42–52. doi: 10.1016/j.envexpbot.2016.05.007

CrossRef Full Text | Google Scholar

Hodges, D. M., DeLong, J. M., Forney, C. F., and Prange, R. K. (1999). Improving the thiobarbituric acid-reactive-substances assay for estimating lipid peroxidation in plant tissues containing anthocyanin and other interfering compounds. Planta 207, 604–611. doi: 10.1007/s004250050524

PubMed Abstract | CrossRef Full Text | Google Scholar

Hossain, Z., López-Climent, M. F., Arbona, V., Pérez-Clemente, R. M., and Gómez-Cadenas, A. (2009). Modulation of the antioxidant system in citrus under waterlogging and subsequent drainage. J. Plant Physiol. 166, 1391–1404. doi: 10.1016/j.jplph.2009.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

IPCC (2014). “Climate change 2014: synthesis report” in Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Core Writing Team], eds R. K. Pachauri and L. A. Meyer (Geneva: IPCC), 151.

Google Scholar

Jin, R., Wang, Y., Liu, R., Gou, J., and Chan, Z. (2016). Physiological and metabolic changes of purslane (Portulaca oleracea L.) in response to drought, heat, and combined stresses. Front. Plant Sci. 6:1123. doi: 10.3389/fpls.2015.01123

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaur, G., and Asthir, B. (2015). Proline: a key player in plant abiotic stress tolerance. Biol. Plant. 59, 609–619. doi: 10.1007/s10535-015-0549-3

CrossRef Full Text | Google Scholar

Koussevitzky, S., Suzuki, N., Huntington, S., Armijo, L., Sha, W., Cortes, D., et al. (2008). Ascorbate peroxidase 1 plays a key role in the response of Arabidopsis thaliana to stress combination. J. Biol. Chem. 283, 34197–34203. doi: 10.1074/jbc.M806337200

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, K.-H. R., Weng, C.-C., Lo, H.-F., and Chen, J.-T. (2004). Study of the root antioxidative system of tomatoes and eggplants under waterlogged conditions. Plant Sci. 167, 355–365. doi: 10.1016/j.plantsci.2004.04.004

CrossRef Full Text | Google Scholar

Mafra, V., Kubo, K. S., Alves-Ferreira, M., Ribeiro-Alves, M., Stuart, R. M., Boava, L. P., et al. (2012). Reference genes for accurate transcript normalization in citrus genotypes under different experimental conditions. PLoS ONE 7:e31263. doi: 10.1371/journal.pone.0031263

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez, V., Mestre, T. C., Rubio, F., Girones-vilaplana, A., Moreno, D. A., Mittler, R., et al. (2016). Accumulation of flavonols over hydroxycinnamic acids favors oxidative damage protection under abiotic stress. Front. Plant Sci. 7:838. doi: 10.3389/fpls.2016.00838

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittler, R. (2006). Abiotic stress, the field environment and stress combination. Trends Plant Sci. 11, 15–19. doi: 10.1016/j.tplants.2005.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittler, R. (2016). ROS are good. Trends Plant Sci. 22, 11–19. doi: 10.1016/j.tplants.2016.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittler, R., Vanderauwera, S., Gollery, M., and Van Breusegem, F. (2004). Reactive oxygen gene network of plants. Trends Plant Sci. 9, 490–498. doi: 10.1016/j.tplants.2004.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgan, J. A. (1984). Interaction of water supply and N in wheat. Plant Physiol. 76, 112–117.

PubMed Abstract | Google Scholar

Pandey, P., Ramegowda, V., and Senthil-Kumar, M. (2015). Shared and unique responses of plants to multiple individual stresses and stress combinations: physiological and molecular mechanisms. Front. Plant Sci. 6:723. doi: 10.3389/fpls.2015.00723

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W., Horgan, G. W., and Dempfle, L. (2002). Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 30, e36.

PubMed Abstract | Google Scholar

Prasch, C. M., and Sonnewald, U. (2013). Simultaneous application of heat, drought, and virus to Arabidopsis plants reveals significant shifts in signaling networks. Plant Physiol. 162, 1849–1866. doi: 10.1104/pp.113.221044

PubMed Abstract | CrossRef Full Text | Google Scholar

Rivero, R. M., Mestre, T. C., Mittler, R., Rubio, F., Garcia-Sanchez, F., and Martinez, V. (2013). The combined effect of salinity and heat reveals a specific physiological, biochemical and molecular response in tomato plants. Plant Cell Environ. 37, 1059–1073. doi: 10.1111/pce.12199

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizhsky, L., Liang, H., and Mittler, R. (2002). The combined effect of drought stress and heat shock on gene expression in tobacco. Plant Physiol. 130, 1143–1151. doi: 10.1104/pp.006858.then

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizhsky, L., Shulaev, V., and Mittler, R. (2004). Measuring programmed cell death in plants. Methods Mol. Biol. 282, 179–189.

Google Scholar

Shao, H., Chu, L., Shao, M., Jaleel, C. A., and Mi, H. (2008). Higher plant antioxidants and redox signaling under environmental stresses. C. R. Biol. 331, 433–441. doi: 10.1016/j.crvi.2008.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Srivastava, G., Kumar, S., Dubey, G., Mishra, V., and Prasad, S. M. (2012). Nickel and ultraviolet-B stresses induce differential growth and photosynthetic responses in Pisum sativum L. seedlings. Biol. Trace Elem. Res. 149, 86–96. doi: 10.1007/s12011-012-9406-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, N., Basil, E., Hamilton, J. S., Inupakutika, M. A., Madhuri, A., Zandalinas, S. I., et al. (2016). ABA is required for plant acclimation to a combination of salt and heat stress. PLoS ONE 11:e0147625. doi: 10.1371/journal.pone.0147625

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, N., Koussevitzky, S., Mittler, R., and Miller, G. (2012). ROS and redox signalling in the response of plants to abiotic stress. Plant Cell Environ. 35, 259–270. doi: 10.1111/j.1365-3040.2011.02336.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, N., Rivero, R. M., Shulaev, V., Blumwald, E., and Mittler, R. (2014). Abiotic and biotic stress combinations. New Phytol. 203, 32–43. doi: 10.1111/nph.12797

PubMed Abstract | CrossRef Full Text | Google Scholar

Szabados, L., and Savouré, A. (2010). Proline: a multifunctional amino acid. Trends Plant Sci. 15, 89–97. doi: 10.1016/j.tplants.2009.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Q. S., Xia, R. X., and Zou, Y. N. (2006). Reactive oxygen metabolism in mycorrhizal and non-mycorrhizal citrus (Poncirus trifoliata) seedlings subjected to water stress. J. Plant Physiol. 163, 1101–1110. doi: 10.1016/j.jplph.2005.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zandalinas, S. I., Balfagón, D., Arbona, V., Gómez-Cadenas, A., Inupakutika, M. A., et al. (2016a). ABA is required for the accumulation of APX1 and MBF1c during a combination of water deficit and heat stress. J. Exp. Bot. 67:erw299. doi: 10.1093/jxb/erw299

PubMed Abstract | CrossRef Full Text | Google Scholar

Zandalinas, S. I., Rivero, R. M., Martínez, V., Gómez-Cadenas, A., and Arbona, V. (2016b). Tolerance of citrus plants to the combination of high temperatures and drought is associated to the increase in transpiration modulated by a reduction in abscisic acid levels. BMC Plant Biol. 16, 105. doi: 10.1186/s12870-016-0791-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zandalinas, S. I., Mittler, R., Balfagón, D., Arbona, V., and Gómez-Cadenas, A. (2017a). Plant adaptations to the combination of drought and high temperatures. Physiol. Plant. doi: 10.1111/ppl.12540 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Zandalinas, S. I., Sales, C., Beltran, J., Gómez-Cadenas, A., and Arbona, V. (2017b). Activation of secondary metabolism in citrus plants is associated to sensitivity to combined drought and high temperatures. Front. Plant Sci. 7:1954. doi: 10.3389/fpls.2016.01954

PubMed Abstract | CrossRef Full Text | Google Scholar