- 1Instituto Ciencias Agroalimentarias Animales y Ambientales (ICA3), Campus Colchagua, Universidad de O′Higgins, San Fernando, Chile
- 2Centro Tecnológico de la Planta Forestal, Instituto Forestal, Sede Biobío, San Pedro de la Paz, Chile
- 3United States Department of Agriculture Forest Service, Rocky Mountain Research Station, Moscow, ID, United States
- 4Laboratorio de Análisis y Modelamiento de Geoinformación, Departamento de Manejo de Bosques y Medio Ambiente, Facultad de Ciencias Forestales, Universidad de Concepción, Concepción, Chile
The search for drought tolerant species or cultivars is important to address water scarcity caused by climate change in Mediterranean regions. The anisohydric–isohydric behavior concept has been widely used to describe stomatal regulation during drought, simply in terms of variation of minimal water potential (Ψmin) in relation to pre-dawn water potential (Ψpd). However, its simplicity has sometimes failed to deliver consistent results in describing a complex behavior that results from the coordination of several plant functional traits. While Prunus dulcis (almond) is known as a drought tolerant species, little information is available regarding consistent metrics to discriminate among cultivars or the mechanisms underlying drought tolerance in almond. Here we show a sequence of plant stomatal, hydraulic, and wilting responses to drought in almonds, and the main differences between anisohydric and isohydric cultivars. In a pot desiccation experiment we observed that stomatal closure in P. dulcis is not driven by loss in turgor or onset of xylem cavitation, but instead, occurs early in response to decreasing Ψmin that could be related to the protection of the integrity of the hydraulic system, independently of cultivar. Also, we report that anisohydric cultivars of P. dulcis are characterized by maximum stomatal conductance, lower water potentials for stomatal closure and turgor loss, and lower vulnerability to xylem cavitation, which are traits that correlated with metrics to discriminate anisohydric and isohydric behavior. Our results demonstrate that P. dulcis presents a strategy to avoid cavitation by closing stomata during the early stages of drought. Future research should also focus on below-ground hydraulic traits, which could trigger stomatal closure in almond.
Introduction
The effects of climate change, such as an increase in temperature and altered patterns of frequency, distribution, and intensity of precipitation, have been extensively characterized (Cavin et al., 2013). Areas with Mediterranean climate, such as central Chile, are especially vulnerable to climate change because more frequent and extreme drought events are expected within the next decades (Schär et al., 2004; Bambach et al., 2013; del Pozo et al., 2019). Central Chile will experience a 20% decline in precipitation by 2050 coupled with an average annual increase of 3–4°C in temperature (Hannah et al., 2013), putting at risk the feasibility of cultivation of several agronomic species. Thus, drought leading to water stress in plants is the major factor reducing agricultural productivity (Chaves et al., 2003), which is one of the main economic activities of Mediterranean central Chile. To face drought in the long term, new suitable species and cultivars with drought tolerance characteristics are being identified and considered as options for the agricultural sector. One new, drought tolerant species being considered for central Chile is almond (Prunus dulcis L.) (Torrecillas et al., 1996), historically cultivated in other Mediterranean regions due to its capability to withstand water stress (Marsal et al., 1997; García-Tejero et al., 2018). Morphological and structural leaf adaptations that deliver protection against excessive water loss, such as reduced leaf area and stomatal density, thicker leaf cell walls, and increased cuticle thickness have been observed in almond cultivars (Camposeo et al., 2011; Oliveira et al., 2018). Physiologically, a gradual decrease in photosynthesis and stomatal conductance with increasing water stress is a characteristic response of drought-adapted plants and it has been frequently observed in almonds (Romero et al., 2004).
At a physiological level, drought is one of the most studied stresses in plants. The primary effect of drought is the reduction in stomatal conductance to decrease water losses by transpiration and preventing drops in leaf water potential (ΨL), which in turn produces a reduction in carbon assimilation (Flexas et al., 2004, 2013), negatively affecting plant growth and crop production. Thus, stomatal regulation is key to control water loss and ΨL during water stress, but different species and cultivars have evolved other adaptive responses to limit dramatic decreases in ΨL, leading to the isohydric and anisohydric behavior concept (Tardieu and Simonneau, 1998). In brief, isohydric species or cultivars maintain higher ΨL due to stomatal closure when the soil dries, whereas anisohydric species or cultivars sustain higher stomatal conductance allowing a drop in ΨL with the progression of water deficit (Tardieu et al., 1996). These behaviors do not represent a simple dichotomy but rather a continuum of strategies in stomatal control. This classification in the behavior of stomata during water stress has led to the search for anisohydric species or cultivars, which could hypothetically sustain higher carbon assimilation rates during drought, and it has been extensively studied in agronomically relevant species such as Vitis vinifera (Schultz, 2003; Soar et al., 2006; Prieto et al., 2010). On the contrary, in almond, the identification of drought-tolerant genotypes has only been focused on physiological changes in response to specific drought treatments, or limited to a narrow range in soil water potentials (Isaakidis et al., 2004; Rouhi et al., 2007; Yadollahi et al., 2011). This does not allow characterization of the stringency of stomatal control with respect to plant water status and the proper identification of anisohydric and isohydric cultivars.
Despite this very simplified classification of stomatal regulation during drought, it has been observed that the anisohydric or isohydric behavior of cultivars, and therefore their classification, depends largely on the soil water content and its hydraulic properties. Tramontini et al. (2014) showed that water relations of V. vinifera (Syrah, near-anisohydric, and Cabernet Sauvignon, near-isohydric) cultivars varied significantly depending on soil hydraulic conductivity, and Hochberg et al. (2018) and Domec and Johnson (2012) indicated that the shift between anisohydric and isohydric behavior depends on soil water potential where the plants water relation parameters are measured in relation to soil properties. Regardless of such inconsistencies in the definition of the anisohydric and isohydric behavior, some methodologies have been described to isolate the effect of environmental variation. The hydroscape, developed by Meinzer et al. (2016), describes the stringency of stomatal regulation of ΨL during the drying of soil. This concept incorporates ranges of pre-dawn water potential (Ψpd) and mid-day minimum water potential (Ψmin) over which stomata are effective in controlling the drop of ΨL during soil desiccation until Ψpd = Ψmin and describes a landscape where the plant is able to sustain photosynthesis until full stomatal closure. Thus, larger hydroscapes are indicative of a higher degree of anisohydry. This methodology has been previously assessed in almond rendering consistent results with other traits related to stomatal stringency during water stress (Álvarez-Maldini et al., 2021).
Hydroscapes have also been correlated with other plant functional traits linked to anisohydric behavior. The water potential at turgor loss point (ΨTLP) is a strong predictor of an overall drought tolerance across species and biomes (Bartlett et al., 2012); a more negative ΨTLP is indicative of anisohydry and higher drought tolerance, a result initially reported by Meinzer et al. (2016) in different species. Likewise, through a meta-analysis, Fu and Meinzer (2018) found a strong correlation between the larger hydroscape areas and the lower ΨTLP values in species across different biomes. Also, Li et al. (2019) reported anisohydric Eucalyptus species with larger hydroscape areas were strongly linked to lower values of ΨTLP, water potential at stomatal closure (Pgs90), maximum stomatal conductance (gsmax), and leaf and stem hydraulic vulnerability to embolism (PL50 and Px50). However, the hydroscape methodology, which has been mainly used to describe stomatal behavior of species distributed across various ecosystems, has yet to be tested broadly across a range of species and cultivars of agricultural interest and in relationship to other hydraulic and functional plant traits. Despite that Meinzer et al. (2016) and Li et al. (2019) have consistently shown a relationship between hydroscape area and traits such as ΨTLP in species from various environments, this link has yet to be proven among cultivars within a species. Given their reduced genetic variability, compared to different species adapted to distinct environments, changes in stomatal behavior and its relationship with other plant traits could be less clear. But ΨTLP and other functional traits have shown to be strong indicators of anisohydric behavior in V. vinifera (Tombesi et al., 2014; Albuquerque et al., 2020; Knipfer et al., 2020), which is in agreement with the multiple mechanisms that are integrated in the hydroscape that regulate plant water status and stomatal behavior.
Despite that this methodology has shown promise to properly characterize anisohydric and isohydric behavior, which in turn delivers information regarding cultivar selection, it is labor-intensive research to examine the regulation of water potential during prolonged periods of soil desiccation. Thus, correlating results from hydroscape analysis with informative functional traits such as ΨTLP, gsmax, or water potential at 50% loss of hydraulic conductivity (P50) would be less time-consuming and a valuable tool to help in the selection of cultivars adapted to drought prone environments. Also, identifying mechanisms and traits underlying drought resistance is essential to predict the future impacts of drought on productive species cultivated in the Mediterranean climate (Martin-StPaul et al., 2017). Considering the current need to assess almond cultivars for cultivation in central Chile, our objective was to describe anisohydric and isohydric behavior in almond cultivars. We used a pot desiccation experiment to explore dynamics of water potential, stomatal conductance, and hydroscape area toward better understanding the relationship of hydroscape area and plant functional traits.
Materials and methods
Plant material and growth conditions
The experiment was conducted between October of 2021 and February of 2022 at the experimental station in the Instituto de Ciencias Agroalimentarias Animales y Ambientales (ICA3) of the Universidad de O′Higgins, Chile (latitude: –34.61°, longitude: –70.99°, elevation: 352 m). The plant material corresponded to three cultivars of almond, Avijor (Ferragnès × Tuono), Isabelona (Blanquerna × Bella d′Aurons), and Soleta (Blanquerna × Bella d′Aurons), that were grafted onto rootpac-20 rootstocks (Prunus besseyi × Prunus cerasifera) 1 year earlier. Plants were donated by Agromillora Sur nursery (Río Claro, Maule region, Chile, latitude: –35.19°, longitude: –71.25°) in October of 2021 and transplanted into 7 L plastic pots (20 cm width × 20 cm length × 25 cm height) using a 1:1 v:v mixture of peat and perlite as substrate. An equal amount of substrate (by weight) was added to each pot (the oven-dry weight of substrate is used in Equation 1, below). The split-plot experiment had two irrigation treatments (WW and PD; whole plots) and three cultivars (Soleta, Isabelona, and Avijor; split plots). We had four blocks (each with an irrigation line) for each whole plot. Block corresponds to random effect, while cultivar and irrigation treatment are fixed effects. Each block had three experimental units (cultivars) and each experimental unit was composed by four pots with an individual plant. Thus, we employed 96 pots [2 irrigation treatments × 4 blocks × 3 cultivars × 4 pots (experimental units) of each cultivar]. Plants were irrigated daily by drip irrigation until they reached container capacity and maintained in this condition for acclimation for 1 month (November 2021) before onset of the irrigation treatments.
Irrigation treatments
After the 1-month acclimation period, the two irrigation treatments commenced. The first irrigation treatment, well-watered (WW), maintained plants at container capacity for the remainder of the experiment. For the second irrigation treatment, pot desiccation (PD), irrigation ceased allowing for a progressive decrease in substrate water content due to plant transpiration (Álvarez-Maldini et al., 2021). To ensure that water losses corresponded to plant transpiration only and not to evaporation from the substrate, the surface of each pot was covered with a plastic film, allowing only the plant stem to emerge. Pot weights were subsequently recorded at midday for the remainder of the experiment. The gravimetric substrate water content (GWC, %) was calculated according to Equation 1:
Where, Pot is the weight of the pot (plus substrate) at each measurement time, Potdry is the weight of the pot plus the initial, oven-dried weight of the substrate, and Potwet is the weight of the pot plus the substrate at container capacity.
Midday and pre-down water potential and pressure-volume curves
Throughout application of the irrigation treatments, pre-dawn water potential (Ψpd) and minimum mid-day water potential (Ψmin) were measured at least eight times for each cultivar since the beginning of the PD treatment, procuring water potential measurements along a wide range of GWC. For leaf water potential measurements, we randomly selected two plants per experimental unit and irrigation treatment (thus, 48 plants total; 2 plants × 2 irrigation treatments × 3 cultivars × 4 replicates). One fully developed leaf from the upper third of the canopy was excised with a razor blade and ΨL was measured with a Scholander pressure chamber model 1505D-EXP (PMS Instruments, Albany, OR, United States) (Scholander et al., 1965). The Ψpd was measured between 06:00 and 07:00 a.m., before dawn, and Ψmin was measured between 13:00 and 14:00 p.m. local time, and both parameters were measured until Ψpd = Ψmin was reached, which was indicative of the limits of stomatal control over plant water status.
Pressure volume (PV) curves were generated from plants of the WW treatment, for which one plant per experimental unit was randomly selected (4 plants total; 1 irrigation treatment × 3 cultivars × 1 plant of each cultivar). Leaf samples were collected before sunrise considering fully developed leaves from the upper third of the plant canopy and cut with a razor blade at the petiole and immediately stored in a sealed bag with damp paper and transported to the laboratory. Leaf PV curves were measured according to the protocol described by Tyree and Hammel (1972). Thus, leaves were allowed to slowly dehydrate under laboratory conditions, and ΨL and fresh mass were measured periodically with an analytical balance (Boeco model A0021E, Hamburg, Germany). The water potential at turgor loss point (ΨTLP) was identified as the inflection point of the 1/ΨL vs. relative water content (RWC) curve. Mean modulus of elasticity (ε) was estimated as the slope of turgor potential (Ψp) vs. RWC in the phase from full turgor to turgor loss point. Capacitance at full turgor (CFT) was calculated from the slope of the linear portion from RWC and ΨL before the ΨTLP, normalized to saturated water content at leaf area.
Gas exchange measurements
On the same days and plants that measurements of Ψpd and Ψmin were recorded, gas exchange was measured at midday with a CIRAS-3 portable photosynthesis system equipped with a CFM-3 chlorophyll fluorescence module (PP Systems, Amesbury, MA, United States). The CO2 concentration in the leaf cuvette was adjusted to 400 ppm, the leaf temperature was maintained at 25 ± 1°C, and the PAR was set to 1,500 μmol photons m–2 s–1. Leaves were acclimated to cuvette conditions for at least 5 min before each measurement. Then, net photosynthesis (AN), stomatal conductance (gs), transpiration rate (E), and sub-stomatal CO2 concentration (Ci) were measured. The instantaneous water use efficiency (iWUE) was calculated as the ratio between AN and E. The maximum photosynthetic and stomatal conductance rates (Amax and gsmax, respectively) was considered as AN and gs measured in plants before onset of the irrigation treatments.
Stem vulnerability to xylem cavitation
Stem vulnerability curves were constructed using the air injection method with a double ended pressure sleeve connected to a Scholander pressure chamber (model described above). We used three plants of each cultivar following the protocol of Ennajeh et al. (2010). We procured the use of branches longer than 40 cm in length to accommodate maximum vessel length to avoid an open-vessel artifact (Cochard et al., 2013). Before the stem was inserted into the pressure-sleeve, the middle portion of bark was removed with a razor blade to facilitate air entry to xylem conduits. Leaves and side branches were removed, and cuts were sealed with glue and parafilm. A flexible tube connected to a solution tank filled with distilled water was connected to the basal end of the stem. The tank was installed 60 cm in height. Then, the stems were flushed at 0.06 MPa for 45 min to remove air bubbles and the maximum conductivity of the stem was measured (Kmax). Flow was measured gravimetrically by collecting the water from the distal end in a pre-weighed, 2-ml Eppendorf tube filled with cotton wool. Flow measurements were recorded for 2-min intervals.
After measuring Kmax, the pressure in the chamber was progressively increased, at 5-min intervals, to 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 5.0, 6.0, and 7.0 MPa. Stem hydraulic conductance (Kh) was measured after slowly releasing the air and allowing the stem to equilibrate for 3-min, then flow measurements were recorded during a 2-min interval. The percentage loss of conductivity (PLC) was calculated using Equation 2 as follows:
The relationship between stem hydraulic conductance and PLC was fitted with a Weibull curve using the fitplc package (Duursma and Choat, 2017) in the R software (version 4.2.0, R Core Team, 2022). The water potential corresponding to a 50% loss in conductivity (P50) and bootstrapped confidence intervals (CI) were calculated.
Metrics to determine anisohydric and isohydric behavior
Considering leaf water relations and gas exchange measurements described above, different metrics to assess anisohydric and isohydric behavior were calculated.
First, the slope of the relationship between Ψpd and Ψmin (σ) was calculating using linear regression according to the methodology described by Martínez-Vilalta et al. (2014).
Second, the hydroscape area (hereafter, hydroscape) was measured following methodology described in Meinzer et al. (2016) and Álvarez-Maldini et al. (2021). In brief, the hydroscape is the area comprising the Ψpd vs. Ψmin regression line and a 1:1 line, which is calculated according to Equation 3:
Where a is the intercept of the Ψpd vs. Ψmin regression line, representing the most negative Ψmin when Ψpd = 0, and b is the intersection of Ψpd vs. Ψmin and the 1:1 line, corresponding to the water potential at Ψpd = Ψmin which indicates the limit of stomatal control to prevent further decrease in leaf water potential.
Third, the water potential at stomatal closure (Ψgs90) was calculated according to Li et al. (2019). Briefly, the stomatal conductance (gs, described in section “Gas exchange measurements”) was plotted against Ψmin and fitted with a weighted polynomial regression to obtain Ψgs90 using the fitplc package from R software (version 4.2.0) (Duursma and Choat, 2017).
Statistical analysis
To assess the relationship between Ψpd and Ψmin, a first-order kinetic model was fitted using PROC REG procedure (SAS Institute Inc., Cary, NC, United States) (Álvarez-Maldini et al., 2021).
Cultivar effects on Amax, gsmax, Emax, iWUEmax, GWC, πo, ΨTLP, RWCTLP, ε, and CFT were assessed using one-way analyses of variance using generalized mixed models using procedure PROC GLIMMIX (SAS Institute Inc., Cary, NC, United States) with selection of distribution considering the Akaike Information Criteria. Differences among means were determined using a Tukey (HSD) test for multiple comparisons.
To evaluate the effect of cultivar and irrigation treatment on AN, E, gs, iWUE, Ψmin and Ψpd an ANOVA analysis was performed according to our experimental design as described above. Differences among means were determined by using a Tukey (HSD) test for multiple comparisons.
To assess the relationship between gs and Ψmin during PD, a first-order kinetic model was fitted using the PROC NLIN procedure (SAS Institute Inc., Cary, NC, United States) with the Gauss-Newton method through a derivative-free algorithm. The cultivar effect was evaluated using the extra sums of squares principle (Bergerud, 1996).
Estimation and inference of regression with piecewise linear model was used to perform the analysis of the water potential (WP) curve. The model generates a linear estimate and calculates the boundaries between linear phases and corresponding slope values according to the methodology described by Muggeo (2003) and adapted by Muggeo et al. (2014) to generate models by different segments with continuous intersection point and with unequal slopes. In our analysis we generated three segments to intersection points or boundaries (Θ1 and Θ2). All the parameters for the model were fitted using a least-squares solver implemented in R software (R Core Team) by “segmented” package v1.6-0 (Muggeo, 2022).
All visualizations were made using SigmaPlot 14 (Systat Software Inc., San Jose, CA, United States).
Results
Initial plant-water relations and gas exchange among cultivars
At the end of the acclimation period and before onset of the irrigation treatments, significant differences were observed in maximum values in the gas exchange parameters of gsmax, Emax, and iWUEmax among cultivars, while no differences were observed in Amax (Table 1). Soleta presented a significantly higher gsmax, followed by Isabelona, and lastly Avijor; a similar trend was observed in Emax, with Soleta and Isabelona reaching significantly higher values, followed by Avijor. Consistently as consequence of higher stomatal conductance, Soleta and Isabelona displayed significantly lower iWUEmax, and Avijor with significantly higher iWUEmax (Table 1).
Table 1. Mean (± standard error) values of leaf photosynthetic traits of Prunus dulcis cultivars measured before the imposition of the pot desiccation treatment.
At the onset of the irrigation treatments the pot water content of Avijor, Isabelona, and Soleta was similar (102.3 ± 1.4%, 99.8 ± 4.0%, and 100.2 ± 2.1%, respectively).
Hydraulic traits of cultivars
In relation to the pressure-volume curve analysis, cultivar was significant for ΨTLP (Table 2; Soleta < Avijor = Isabelona) but not for SWC, πo, RWCTLP, CFT, or ε traits. We observed a trend (p = 0.0791) for higher πo in Soleta, followed by Isabelona, and then Avijor.
Table 2. Values of key metrics to describe anisohydric and isohydric behavior in Prunus dulcis cultivars and of hydraulic traits derived from pressure volume curves (mean ± standard deviation), stem vulnerability curves (mean, CI in brackets) and leaf water potential causing 90% of stomatal closure (Ψgs90) for different almond cultivars.
Water potential at 50% loss of conductivity (P50) showed different values among cultivars. Soleta and Isabelona had similar P50 values (–3.73 MPa and –3.80 MPa, respectively) that were lower than Avijor (–2.97 MPa) (Table 2).
Different behavior of cultivars during pot desiccation
The PD experiment lasted 50 days between December of 2021 until January of 2022. Water potential (Ψmin and Ψpd) and gas exchange parameters were measured 1, 9, 15, 20, 23, 27, 35, 44, and 50 days after the beginning of the PD treatment. During the PD treatment imposition, GWC content decreased steadily, reaching pot water contents of 35.9 ± 3.7%, 31.9 ± 1.9%, and 30.7 ± 4.8% in Avijor, Isabelona, and Soleta, respectively. Pot water content in the WW treatment was 89.6 ± 5.8% in Avijor, 90.0 ± 1.3% in Isabelona, and 86.6 ± 4.8% in Soleta.
Among cultivars, all leaf water relations and gas exchange parameters were affected by the PD treatment. At conclusion of the PD treatment, AN was significantly affected by cultivar (p < 0.0001), and iWUE was only affected by the PD treatment (p = 0.0277) with significantly lower values in stressed plants. Minimum water potential was independently affected by the irrigation treatment and cultivar, while both factors significantly interacted to affect E, gs and Ψpd (Table 3). At the end of the PD treatment and for all cultivars, Ψmin decreased significantly compared with the WW treatment, reaching minimum values in Soleta followed by Avijor and by Isabelona. In all the parameters affected by the irrigation × cultivar interaction (E, gs and Ψpd), the PD treatment significantly decreased values compared to the WW treatment. In regard to gs, at the end of the PD treatment, Isabelona reached the lowest value, with higher gs in Soleta. Transpiration rate (E) was higher in Soleta in the WW treatment compared with Avijor and Isabelona, but no differences among cultivars were observed in the PD treatment (Table 3). The lowest values of Ψpd at the end of the PD treatment for Avijor, Isabelona, and Soleta were –2.93 ± 0.29, –2.39 ± 0.29, and –3.38 ± 0.12 MPa, respectively (Table 3).
Table 3. Mean (± standard deviation) values of leaf photosynthetic traits and water potential at conclusion of the pot desiccation (PD treatment), source of variation, and p-values of Prunus dulcis cultivars subjected to PD and well-watered (WW) irrigation treatments.
Although, as previously described, the PD treatment affected leaf water potential and gas exchange parameters at the end of the experiment, the desiccation treatment influenced the behavior of stomatal conductance differently, with decreasing Ψmin of each cultivar during the first stages of the desiccation treatment (Figure 1). Thus, different models were fitted for Avijor, Isabelona, and Soleta. A faster decrease in gs was observed in Avijor at higher Ψmin values. On the contrary, Soleta sustained higher stomatal conductance for a longer period (at more negative Ψmin) during the progression of PD, while Isabelona displayed an intermediate behavior between Soleta and Avijor (Figure 1). This pattern was similar for water potential at stomatal closure (Ψgs90) (Table 2), with Soleta and Isabelona displaying lower values of Ψgs90 (–2.14 MPa and –2.25 MPa, respectively) and Avijor displaying a higher value (–1.72 MPa), although the large confidence intervals of Avijor overlapped with the Ψgs90 of the other two cultivars.
Figure 1. Effect of midday minimum water potential (Ψmin) on stomatal conductance (gs) during the pot desiccation treatment of Prunus dulcis cultivars. Symbols indicate measured data and lines correspond modeled data for each cultivar. Fitted models and corresponding R2-value for all cultivars are shown in the panel.
Regarding the metrics of stringency of stomatal control, the slope of the linear regression fitted to the relationship of Ψpd vs. Ψmin (σ) revealed that σ only ranged between 0.811 MPa MPa–1 in Isabelona and 0.922 MPa MPa–1 in Avijor (Table 2), with none of the cultivars displaying a perfect isohydric or anisohydric behavior (slopes 0 or 1, respectively).
Despite not observing perfect anisohydric or isohydric behavior, the hydroscape area ranged from 8.74 MPa2 in Avijor to 9.223 MPa2 in Soleta (Figure 2 and Table 2), with Isabelona displaying an intermediate value (Table 2). Thus, cultivar rankings based on hydroscape area from the more isohydric to the more anisohydric is: Isabelona, Avijor, and Soleta.
Figure 2. Trajectories of pre-dawn water potential (Ψpd) vs. minimum midday water potential (Ψmin) during the pot desiccation treatment for the Prunus dulcis cultivars Avijor (A), Isabelona (B), and Soleta (C). Solid black line represents the 1:1 line, dashed red line is the fitted linear regression of measured data. Different symbols in colors represent different individual plants measured. Hydroscape is represented by the area bounded by the regression and the 1:1 line.
As previously demonstrated by Knipfer et al. (2020), the relationship between Ψpd and Ψmin during drought is not linear, and is constituted by three phases of plant dehydration with break points between each phase that are indicative of the thresholds of stomatal closure and leaf turgor loss. Accordingly, we were able to perform the analysis of the water potential curve (WP curve) during pot desiccation for each cultivar and identify three distinct phases along the WP curve (Figure 3 and Table 4), although it displayed a different shape than the one described by Knipfer et al. (2020). According to the slope of the three phases (Table 4), a steeper decrease in water potential was observed during the first stages of the PD treatment, which continued until –0.803 MPa in Avijor, –0.334 MPa in Isabelona, and –0.997 MPa in Soleta, corresponding to the boundary between the first and second phase (Θ1). Following Θ1, reductions in water potential were less pronounced in the second and third phases, observed through lower slope values (β2 and β3) (Table 4). The boundary between the second and the third phase (Θ2) was higher in Avijor (–1.553 MPa), Isabelona (–2.941 MPa), and Soleta (–3.513 MPa) (Figure 3).
Figure 3. Adjustment of the piecewise linear regression (PLR) model on the relationship between pre-dawn water potential (Ψpd) and minimum mid-day water potential (Ψmin) on the Prunus dulcis cultivars Avijor (A), Isabelona (B), and Soleta (C). Black circles indicate data from each individual plant, blue line is the fitted model. Vertical solid lines are boundaries between phases I and II (Θ1) and between phases II and III (Θ2), corresponding SE values are displayed in gray lines. Model output parameters are summarized in Table 4.
Table 4. Summary of output parameters from the piecewise linear regression (PLR) model as described by Knipfer et al. (2020) used for the analysis between Ψpd and Ψmin for Prunus dulcis cultivars.
Discussion
The anisohydric and isohydric concept has been used for several decades as a framework to describe different behaviors of plant water-relations during drought. However, this concept is not exempt from criticism due to common lack in consistency, with species or cultivars shifting from isohydric to anisohydric behavior, and vice versa (Hochberg et al., 2018). Isohydric and anisohydric behavior are simply described in terms of changes in Ψmin in response to Ψpd, although it represents a whole-plant hydraulic strategy resulting from the coordination and trade-offs among different plant functional traits (Meinzer et al., 2017; Fu et al., 2019; Henry et al., 2019). Thus, using plant functional traits to predict species distribution, dynamics, and responses to environmental change has been important in physiological ecology (Adier et al., 2014), but is now rapidly gaining usefulness in the search for cultivars of interest in the agricultural industry. To fully comprehend the behavior in stomatal regulation during water stress of three almond cultivars, we assessed several metrics of anisohydric and isohydric behavior and leaf and stem functional traits for a more comprehensive framework to characterize plant response to drought (Skelton et al., 2015).
Initially, we assessed the slope of the Ψmin vs. Ψpd (σ) as previously described by Martínez-Vilalta et al. (2014), but none of the cultivars displayed neither perfect isohydric nor anisohydric behavior (0 or 1 values, respectively). Although σ showed a tendency of Avijor to have a more isohydric behavior in this study, it has failed to correlate with other metrics of stomatal sensitivity in 44 species (Martínez-Vilalta and Garcia-Forner, 2017) and it has proven to be inconsistent to rank species/cultivars along a continuum of isohydric to anisohydric depending of the values of Ψpd (Poni et al., 2007; Domec and Johnson, 2012; Hochberg et al., 2018). We also used the hydroscape metric to assess the anisohydric and isohydric behavior of the three almond cultivars. According to Meinzer et al. (2016), the hydroscapes describe a water potential landscape over which plants operate before drought-induced stomatal closure, with anisohydry characterized by larger areas. Our results regarding hydroscape values are similar to the more anisohydric species reported by Meinzer et al. (2016) with 8.4 MPa2 in Quercus garryana, and higher to the values reported by Li et al. (2019) in Eucalyptus species. Thus, according to this metric the rank of our cultivars from more isohydric to more anisohydric is Avijor, Isabelona, and Soleta. Despite the above-mentioned issues of σ, both metrics (hydroscapes and σ) ranked the cultivars equally. These consistent results could be related to the wide range in Ψpd we measured. A wide range of hydroscape areas have been reported among different species (Meinzer et al., 2016; Li et al., 2019), while in this study a lower range in area variation among cultivars (8.13 MPa2 vs. 9.23 MPa2) was observed. Regardless of this lower variation, cultivars also displayed differential behavior observed in stomatal responses to decreasing water availability (Figure 1). While Avijor had a faster decrease in stomatal conductance with decreasing Ψmin, Soleta sustained higher stomatal conductance at lower values of Ψmin; this points toward an isohydric behavior in Avijor and an anisohydric behavior in Soleta. This agrees with Tombesi et al. (2014) who reported stomatal closure at higher leaf and stem water potentials in a near-isohydric V. vinifera cultivar (Montepulciano) vs. the anisohydric cultivar (Sangiovese).
Despite reports that species with higher sensitivity to stomatal closure during leaf dehydration also have higher values of gsmax (Skelton et al., 2015; Henry et al., 2019), our results showed that the more anisohydric cultivar Soleta, which sustained higher gs at lower Ψmin, also presented higher gsmax, Emax, and lower iWUE (Table 1) under conditions of high water availability. This concurs with Meinzer et al. (2017), who reported that anisohydric species have faster kinetics of stomatal opening and activation of photosynthesis with greater stomatal conductance, photosynthetic capacity, and lower iWUE, characteristic of species from dry regions that aim to maximize the utilization of unpredictable and rare precipitation events (Li et al., 2016, 2018). This trait has also been observed between anisohydric and isohydric Vitis cultivars (Dayer et al., 2020). These results indicate that anisohydric cultivars such as Soleta are characterized by physiological traits linked to drought tolerant species that maximize carbon gain during periods of high water availability.
As expected, the imposition of the PD experiment caused an increase in iWUE in all cultivars at the beginning of the experiment (Supplementary Figure 1). However, during the first 30 days of the PD treatment, the more anisohydric cultivars Soleta and Isabelona sustained higher iWUE than the isohydric-behaving Avijor (Supplementary Figure 1). Higher iWUE during PD in Soleta and Isabelona are explained mainly by maintenance of higher photosynthesis rates than Avijor during PD. Although it has been previously reported than isohydric cultivars are expected to present higher iWUE (Schultz, 1996; Poni et al., 2007), our results agree with Pou et al. (2012) who indicates that anisohydric Vitis cultivars show higher WUE due to maintenance of photosynthesis stability.
The water potential at stomatal closure (Ψgs90) is also a trait that can accurately quantify anisohydric and isohydric behavior (Skelton et al., 2015; Hochberg et al., 2018), where a more negative Ψgs90 indicates that plants can sustain carbon assimilation over a wider range of leaf water potentials, which is in direct relationship with larger hydroscape areas. In agreement, we observed that the more anisohydric cultivars Soleta and Isabelona had lower Ψgs90, in contrast with the more isohydric cultivar Avijor (Table 2); a similar trend was also observed in ΨTLP values among cultivars. The water potential at turgor loss point has been widely used as an indicator of drought tolerance, as leaves with lower ΨTLP maintain open stomata, hydraulic conductance, and are able to grow under dryer conditions (Bartlett et al., 2012; Scoffoni et al., 2012; Trifilò et al., 2015). Recently, lower ΨTLP values have been positively correlated with larger hydroscapes areas, and thus anisohydry (Meinzer et al., 2016; Li et al., 2019), which supports our results regarding the contrasting anisohydric behavior of Soleta and the isohydric behavior of Avijor. Our values of ΨTLP and Ψgs90 were also similar to the ones reported by Hernandez-Santana et al. (2016) in almond (–2.26 and –2.14 MPa, respectively), which are in the range of values displayed by drought-tolerant species.
Although no significant differences were observed in π0 among cultivars (p = 0.0791) (Table 2), the trend in values with higher π0 in Soleta and lower in Avijor, and the absolute lack of differences in ε, support the notion that contrast observed in ΨTLP among cultivars is mainly driven by osmotic adjustment and not by modifications of mechanical properties of the cell wall. This is consistent with a previous meta-analysis (Bartlett et al., 2012) indicating that π0 is the main driver of ΨTLP across species, and that ε does not play a role in determining drought-tolerance. Regardless of the relationship between ΨTLP and Ψgs90, our results suggest that ΨTLP is not the main factor triggering stomatal closure because, independently of the cultivar, complete stomatal closure occurred before turgor loss (Ψgs90 > ΨTLP). This indicates that stomatal closure is uncoupled from bulk leaf water status. Similar results were reported in anisohydric and isohydric Vitis cultivars (Tombesi et al., 2014) and several woody species adapted to sites with contrasting resource availability (Henry et al., 2019).
Consistent with the previously described traits, the P50 displayed that same variation among almond cultivars, with more negative values in Soleta and Isabelona than in Avijor (Table 2), indicating a higher vulnerability to cavitation in Avijor. Information regarding the relationship between vulnerability to cavitation and anisohydric and isohydric behavior is scarce for cultivars, but Alsina et al. (2007) reported large variability in PCL50 across Vitis cultivars (–1 to –3 MPa). Reports in different species have found a correlation between a higher xylem vulnerability to cavitation with isohydric behavior, such as in Acer pseudoplatanus and Corylus avellana (Li et al., 2016). It also has been proposed that xylem cavitation is the hydraulic signal that triggers stomatal closure in several species (Salleo et al., 2000; Nardini et al., 2001; Zufferey et al., 2011). On the contrary, the differences between Ψgs90 and ΨTLP in our results showed that stomatal closure occurred before the 50% loss of xylem hydraulic conductivity in the three cultivars (Ψgs90 > P50). Recent evidence also supports the notion that stomatal closure precedes xylem cavitation at community and global scales (Klein, 2014; Bartlett et al., 2016; Martin-StPaul et al., 2017) and even among cultivars (Carminati and Javaux, 2020; Dayer et al., 2020; Gambetta et al., 2020). These results indicate that the hydraulic safety margin (HSM), defined as the remainder of Ψgs90 minus P50, is generally positive and thus constitutes a strategy to avoid lethal embolism.
As another metric to describe the dynamics of leaf water potential during pot desiccation, we performed an analysis of the water potential curve using a piecewise linear regression according to Knipfer et al. (2018). Here, authors describe the triphasic nature of the WP curve and the prediction of the first boundary (Θ1) between phases I and II that would correspond to the threshold at which leaf gas exchange is substantially decreased, and the second boundary (Θ2) matching the ΨTLP. We were able to model the three phases of the WP curve and predict its corresponding boundaries and slopes (Table 4), but the shape of our curve diverged from the one described by Knipfer et al. (2020). Specially, a steeper decrease in water potential was observed in phase II in our cultivars while Knipfer et al. (2020) showed a minor reduction in Ψmin with the continuous decrease in Ψpd. A faster dry-down in our experiment could explain the steeper reduction in Ψmin during the PD treatment, indicating that factors that affect the rates of water loss and transpiration also influence dynamics of water potential. Also, the almond experiment described by Knipfer et al. (2020) was performed in 56-L potted plants, which can affect the pace of the dry-down, and thus, the dynamics of water potential, explaining the steeper decline in Ψmin during phase II of our experiment. Differences in cultivar/rootstock combination could also affect the shape of the WP curve. Our predicted boundaries, Θ1 and Θ2, did not match stomatal closure nor ΨTLP, respectively, as indicated by Knipfer et al. (2020). It is interesting to note, however, the smaller distance between Θ1 and Θ2 boundaries in the isohydric Avijor (Figure 3A), compared to the larger distance described between the two boundaries in the more anisohydric cultivars Isabelona and Soleta (Figures 3B,C). Along the three phases of the WP curve, a consistent reduction in the slope of each phase was observed (Table 4), thus the faster reduction in Ψmin occurred during the initial stages of the PD treatment until Θ1, which could have triggered the onset of stomatal closure.
Considering all the traits described above, we observed responses to the PD treatment that suggest almond has an avoidance strategy rather than a tolerance strategy for hydraulic failure as described by Hochberg et al. (2017). First, the steep decrease in Ψmin in response to decrease in Ψpd (Figure 1) eventually led to stomatal closure (Ψgs90) at water potentials higher than the loss of cell turgor (ΨTLP) and preceded development of cavitation in xylem conduits (P50) (Figure 4). This concurs with recent research stating that drought-resistant plants close their stomata at water potentials higher than that at which substantial embolism occurs (Martin-StPaul et al., 2017). Second, we observed leaf shedding in all cultivars, especially Avijor, which further supports an avoidance strategy toward protection of perennial tissues. These responses did, however, prevent us from observing a direct measurement of Ψmin = Ψpd. Thus, for hydroscape calculation, Ψmin = Ψpd was estimated by extrapolation of water potential values at the end of the PD experiment, rendering unreliable Ψmin = Ψpd values well beyond limits of stomatal control. Then, because the WP curves did not follow a strictly linear fit (Table 4 and Figure 3), we estimated Ψmin = Ψpd as the intersection between a 1:1 line and the modeled piecewise regression, which yielded Ψmin = Ψpd values for Avijor (–6.15 MPa), Isabelona (–4.33 MPa), and Soleta (–4.83 MPa) that were closer to the observed water potential values for turgor loss and loss of hydraulic conductivity.
Figure 4. Relationships among hydroscape area and traits linked to anisohydric and isohydric behavior in the Prunus dulcis cultivars Avijor Isabelona, and Soleta. Ψgs90: leaf water potential causing 90% of stomatal closure, ΨTLP: water potential at turgor loss point, and P50: water potential at 50% loss of hydraulic conductivity. Values of maximum transpiration rate (Emax) are shown for each cultivar.
Despite the above-mentioned issues when calculating hydroscapes, a clear relationship was observed between hydroscape area and traits related to anisohydric behavior. Thus, Soleta cultivar with a near-anisohydric behavior displayed larger hydroscape area with lower values of Ψgs90, ΨTLP, and P50, contrary to the behavior displayed by Avijor (Figure 4). Thus, considering the above-mentioned traits for Soleta and Isabelona with an increased iWUE during the first stages of drought indicate, both cultivars are suitable candidates to be established under Mediterranean climate conditions of central Chile, especially in areas with restricted access to irrigation water.
Furthermore, our results regarding Soleta cultivar are in agreement with Álvarez et al. (2020) who indicated that this cultivar displayed physiological responses linked to drought tolerance when growing with its original root system (self-rooted). This result was, however, diminished when plants were grafted onto the dwarfing rootstock Rootpac-20. These results correlate with Ben Yahmed et al. (2016) that described a poor adaptation of this rootstock to Mediterranean or rain-fed conditions. Thus, considering these reports, it is important to highlight the influence of the rootstock on drought tolerance, and that the anisohydric behavior of Soleta and Isabelona cultivars may have been attenuated by the Rootpac-20 rootstock. Future work should focus on the effects of rootstock on anisohydric and isohydric behavior and consider the use of self-rooted plants for cultivation under Mediterranean conditions.
Conclusion
Our results from a pot desiccation treatment revealed that almond cultivars followed a similar response to drought: decreasing Ψmin led to stomatal closure to avoid hydraulic disfunction. The different thresholds of water potentials indicate that neither turgor loss nor xylem cavitation are the hydraulic signal to induce stomatal closure, but rather that loss in soil hydraulic conductivity could be the main trigger of stomatal closure. Thus, future research regarding almond and the differences in anisohydric and isohydric behavior among its cultivars should account for below- and above-ground hydraulic traits.
Cultivars did, however, display differences in anisohydric and isohydric behavior. Anisohydric cultivars, Soleta and Isabelona, displayed larger hydroscape areas. Moreover, they presented different boundaries along the WP curve, higher gsmax and sustained higher stomatal conductance at lower Ψmin, lower ΨTLP and vulnerability to xylem cavitation (P50), and higher iWUE than the isohydric cultivar Avijor. Because ΨTLP, P50, and gsmax are easily measured in well-watered plants, their assessment could be readily used to discriminate anisohydric and isohydric plants under consideration for agronomic development in Mediterranean climates without the laborious need for controlled drought experiments. Selecting plants more anisohydric in nature, such as Soleta in this research, with their ability to withstand cavitation better than their isohydric cohorts, may be preferred given expected changes in Mediterranean climates.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
CÁ-M, MP, and MA conceived and designed the experiment. CÁ-M, DE, and FA performed physiological measurements. MA and SS: contributed to data analysis. CÁ-M and RD wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Agencia Nacional de Investigación y Desarrollo (ANID) Proyecto Fondecyt de Iniciación (no. 11200807).
Acknowledgments
We thank to Agromillora Sur S.A, especially to Natalia Torres and Cristobal Calvo, for kindly providing the plant material to conduct this investigation.
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.
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.2022.974050/full#supplementary-material
References
Adier, P. B., Salguero-Gómez, R., Compagnoni, A., Hsu, J. S., Ray-Mukherjee, J., Mbeau-Ache, C., et al. (2014). Functional traits explain variation in plant lifehistory strategies. Proc. Natl. Acad. Sci. U.S.A. 111, 740–745. doi: 10.1073/PNAS.1315179111/SUPPL_FILE/PNAS.201315179SI.PDF
Albuquerque, C., Scoffoni, C., Brodersen, C. R., Buckley, T. N., Sack, L., and McElrone, A. J. (2020). Coordinated decline of leaf hydraulic and stomatal conductances under drought is not linked to leaf xylem embolism for different grapevine cultivars. J. Exp. Bot. 71, 7286–7300. doi: 10.1093/jxb/eraa392
Alsina, M., De Herralde, F., Aranda, X., Savé, R., and Biel, C. (2007). Water relations and vulnerability to embolism are not related: Experiments with eight grapevine cultivars. Vitis 46, 1–6.
Álvarez, S., Martín, H., Barajas, E., Rubio, J. A., and Vivaldi, G. A. (2020). Rootstock effects on water relations of young almond trees (cv. Soleta) when subjected to water stress and rehydration. Water 12:3319. doi: 10.3390/W12123319
Álvarez-Maldini, C., Acevedo, M., and Pinto, M. (2021). Hydroscapes: A useful metric for distinguishing iso-/anisohydric behavior in almond cultivars. Plants 10, 1–11. doi: 10.3390/plants10061249
Bambach, N., Meza, F. J., Gilabert, H., and Miranda, M. (2013). Impacts of climate change on the distribution of species and communities in the Chilean Mediterranean ecosystem. Reg. Environ. Chang. 13, 1245–1257. doi: 10.1007/S10113-013-0425-7/TABLES/2
Bartlett, M. K., Klein, T., Jansen, S., Choat, B., and Sack, L. (2016). The correlations and sequence of plant stomatal, hydraulic, and wilting responses to drought. Proc. Natl. Acad. Sci. U.S.A. 113, 13098–13103. doi: 10.1073/PNAS.1604088113/SUPPL_FILE/PNAS.1604088113.SD03.XLSX
Bartlett, M. K., Scoffoni, C., and Sack, L. (2012). The determinants of leaf turgor loss point and prediction of drought tolerance of species and biomes: A global meta-analysis. Ecol. Lett. 15, 393–405. doi: 10.1111/j.1461-0248.2012.01751.x
Ben Yahmed, J., Ghrab, M., and Ben Mimoun, M. (2016). Eco-physiological evaluation of different scion-rootstock combinations of almond grown in Mediterranean conditions. Fruits 71, 185–193. doi: 10.1051/FRUITS/2016003
Bergerud, W. A. (1996). Introduction to logistic regression models: With worked forestry examples: biometrics information handbook no.7. Working Paper. 157. Victoria, BC: British Columbia Ministry of Forests.
Camposeo, S., Palasciano, M., Vivaldi, G. A., and Godini, A. (2011). Effect of increasing climatic water deficit on some leaf and stomatal parameters of wild and cultivated almonds under Mediterranean conditions. Sci. Hortic. 127, 234–241. doi: 10.1016/j.scienta.2010.09.022
Carminati, A., and Javaux, M. (2020). Soil rather than xylem vulnerability controls stomatal response to drought. Trends Plant Sci. 25, 868–880. doi: 10.1016/j.tplants.2020.04.003
Cavin, L., Mountford, E. P., Peterken, G. F., and Jump, A. S. (2013). Extreme drought alters competitive dominance within and between tree species in a mixed forest stand. Funct. Ecol. 27, 1424–1435. doi: 10.1111/1365-2435.12126
Chaves, M. M., Maroco, J. P., and Pereira, J. S. (2003). Understanding plant responses to drought - from genes to the whole plant. Funct. Plant Biol. 30, 239–264. doi: 10.1071/FP02076
Cochard, H., Badel, E., Herbette, S., Delzon, S., Choat, B., and Jansen, S. (2013). Methods for measuring plant vulnerability to cavitation: A critical review. J. Exp. Bot. 64, 4779–4791. doi: 10.1093/jxb/ert193
Dayer, S., Herrera, J. C., Dai, Z., Burlett, R., Lamarque, L. J., Delzon, S., et al. (2020). The sequence and thresholds of leaf hydraulic traits underlying grapevine varietal differences in drought tolerance. J. Exp. Bot. 71, 4333–4344. doi: 10.1093/jxb/eraa186
del Pozo, A., Brunel-Saldias, N., Engler, A., Ortega-Farias, S., Acevedo-Opazo, C., Lobos, G. A., et al. (2019). Climate change impacts and adaptation strategies of agriculture in Mediterranean-Climate Regions (MCRs). Sustainability 11:2769. doi: 10.3390/SU11102769
Domec, J. C., and Johnson, D. M. (2012). Does homeostasis or disturbance of homeostasis in minimum leaf water potential explain the isohydric versus anisohydric behavior of Vitis vinifera L. cultivars? Tree Physiol. 32, 245–248.
Duursma, R. A., and Choat, B. (2017). fitplc - an R package to fit hydraulic vulnerability curves. J. Plant Hydraul. 4:e002. doi: 10.20870/jph.2017.e002
Ennajeh, M., Vadel, A. M., Cochard, H., and Khemira, H. (2010). Comparative impacts of water stress on the leaf anatomy of a drought-resistant and a drought-sensitive olive cultivar. J. Hortic. Sci. Biotechnol. 85, 289–294. doi: 10.1080/14620316.2010.11512670
Flexas, J., Bota, J., Cifre, J., Escalona, J. M., Galmés, J., Gulías, J., et al. (2004). Understanding down-regulation of photosynthesis under water stress: Future prospects and searching for physiological tools for irrigation management. Ann. Appl. Biol. 144, 273–283. doi: 10.1111/j.1744-7348.2004.tb00343.x
Flexas, J., Scoffoni, C., Gago, J., and Sack, L. (2013). Leaf mesophyll conductance and leaf hydraulic conductance: An introduction to their measurement and coordination. J. Exp. Bot. 64, 3965–3981. doi: 10.1093/jxb/ert319
Fu, X., and Meinzer, F. C. (2018). Metrics and proxies for stringency of regulation of plant water status (iso/anisohydry): A global data set reveals coordination and trade-offs among water transport traits. Tree Physiol. 39, 122–134. doi: 10.1093/treephys/tpy087
Fu, X., Meinzer, F. C., Woodruff, D. R., Liu, Y. Y., Smith, D. D., McCulloh, K. A., et al. (2019). Coordination and trade-offs between leaf and stem hydraulic traits and stomatal regulation along a spectrum of isohydry to anisohydry. Plant Cell Environ. 42, 2245–2258. doi: 10.1111/pce.13543
Gambetta, G. A., Herrera, J. C., Dayer, S., Feng, Q., Hochberg, U., and Castellarin, S. D. (2020). The physiology of drought stress in grapevine: Towards an integrative definition of drought tolerance. J. Exp. Bot. 71, 4658–4676. doi: 10.1093/jxb/eraa245
García-Tejero, I. F., Gutiérrez Gordillo, S., Souza, L., Cuadros-Tavira, S., and Durán Zuazo, V. H. (2018). Fostering sustainable water use in almond (Prunus dulcis Mill.) orchards in a semiarid Mediterranean environment. Arch. Agron. Soil Sci. 65, 164–181. doi: 10.1080/03650340.2018.1492113
Hannah, L., Roehrdanz, P. R., Ikegami, M., Shepard, A. V., Shaw, M. R., Tabor, G., et al. (2013). Climate change, wine, and conservation. Proc. Natl. Acad. Sci. U.S.A. 110, 6907–6912. doi: 10.1073/PNAS.1210127110/SUPPL_FILE/PNAS.201210127SI.PDF
Henry, C., John, G. P., Pan, R., Bartlett, M. K., Fletcher, L. R., Scoffoni, C., et al. (2019). A stomatal safety-efficiency trade-off constrains responses to leaf dehydration. Nat. Commun. 10:3398. doi: 10.1038/s41467-019-11006-1
Hernandez-Santana, V., Rodriguez-Dominguez, C. M., Fernández, J. E., and Diaz-Espejo, A. (2016). Role of leaf hydraulic conductance in the regulation of stomatal conductance in almond and olive in response to water stress. Tree Physiol. 36, 725–735. doi: 10.1093/treephys/tpv146
Hochberg, U., Rockwell, F. E., Holbrook, N. M., and Cochard, H. (2018). Iso/Anisohydry: A plant–environment interaction rather than a simple hydraulic trait. Trends Plant Sci. 23, 112–120. doi: 10.1016/j.tplants.2017.11.002
Hochberg, U., Windt, C., Ponomarenko, A., Zhang, Y.-J., Gersony, J., Rockwell, F., et al. (2017). Stomatal closure, basal leaf embolism, and shedding protect the hydraulic integrity of grape stems. Plant physiol. 174, 764–775. doi: 10.1104/pp.16.01816
Isaakidis, A., Sotiropoulos, T., Almaliotis, D., Therios, I., and Stylianidis, D. (2004). Response to severe water stress of the almond (Prunus amygdalus) ‘ferragnès’ grafted on eight rootstocks. N. Z. J. Crop Hortic. Sci. 32, 355–362. doi: 10.1080/01140671.2004.9514316
Klein, T. (2014). The variability of stomatal sensitivity to leaf water potential across tree species indicates a continuum between isohydric and anisohydric behaviours. Funct. Ecol. 28, 1313–1320. doi: 10.1111/1365-2435.12289
Knipfer, T., Bambach, N., Hernandez, M. I., Bartlett, M. K., Sinclair, G., Duong, F., et al. (2020). Predicting stomatal closure and turgor loss in woody plants using predawn and midday water potential. Plant Physiol. 184, 881–894. doi: 10.1104/pp.20.00500
Knipfer, T., Barrios-Masias, F. H., Cuneo, I. F., Bouda, M., Albuquerque, C. P., Brodersen, C. R., et al. (2018). Variations in xylem embolism susceptibility under drought between intact saplings of three walnut species. Tree Physiol. 38, 1180–1192. doi: 10.1093/treephys/tpy049
Li, S., Feifel, M., Karimi, Z., Schuldt, B., Choat, B., and Jansen, S. (2016). Leaf gas exchange performance and the lethal water potential of five European species during drought. Tree Physiol. 36, 179–192. doi: 10.1093/TREEPHYS/TPV117
Li, X., Blackman, C. J., Choat, B., Duursma, R. A., Rymer, P. D., Medlyn, B. E., et al. (2018). Tree hydraulic traits are coordinated and strongly linked to climate-of-origin across a rainfall gradient. Plant Cell Environ. 41, 646–660. doi: 10.1111/pce.13129
Li, X., Blackman, C. J., Peters, J. M. R., Choat, B., Rymer, P. D., Medlyn, B. E., et al. (2019). More than iso/anisohydry: Hydroscapes integrate plant water use and drought tolerance traits in 10 eucalypt species from contrasting climates. Funct. Ecol. 33, 1035–1049. doi: 10.1111/1365-2435.13320
Marsal, J., Girona, J., and Mata, M. (1997). Leaf water relation parameters in almond compared to hazelnut trees during a deficit irrigation period. J. Am. Soc. Hortic. Sci. 122, 582–587. doi: 10.21273/JASHS.122.4.582
Martínez-Vilalta, J., and Garcia-Forner, N. (2017). Water potential regulation, stomatal behaviour and hydraulic transport under drought: deconstructing the iso/anisohydric concept. Plant Cell Environ. 40, 962–976. doi: 10.1111/pce.12846
Martínez-Vilalta, J., Poyatos, R., Aguadé, D., Retana, J., and Mencuccini, M. (2014). A new look at water transport regulation in plants. New Phytol. 204, 105–115. doi: 10.1111/NPH.12912
Martin-StPaul, N., Delzon, S., and Cochard, H. (2017). Plant resistance to drought depends on timely stomatal closure. Ecol. Lett. 20, 1437–1447. doi: 10.1111/ele.12851
Meinzer, F. C., Smith, D. D., Woodruff, D. R., Marias, D. E., McCulloh, K. A., Howard, A. R., et al. (2017). Stomatal kinetics and photosynthetic gas exchange along a continuum of isohydric to anisohydric regulation of plant water status. Plant Cell Environ. 40, 1618–1628. doi: 10.1111/pce.12970
Meinzer, F. C., Woodruff, D. R., Marias, D. E., Smith, D. D., McCulloh, K. A., Howard, A. R., et al. (2016). Mapping ‘hydroscapes’ along the iso- to anisohydric continuum of stomatal regulation of plant water status. Ecol. Lett. 19, 1343–1352. doi: 10.1111/ele.12670
Muggeo, V. M. R. (2003). Estimating regression models with unknown break-points. Stat. Med. 22, 3055–3071. doi: 10.1002/SIM.1545
Muggeo, V. M. R. (2022). “segmented”: Regression models with break-points/change-points (with Possibly Random Effects) estimation. R Package version 1.6-0. Available online at: https://cran.r-project.org/web/packages/segmented/segmented.pdf
Muggeo, V. M. R., Atkins, D. C., Gallop, R. J., and Dimidjian, S. (2014). Segmented mixed models with random changepoints: a maximum likelihood approach with application to treatment for depression study. Stat. Modell. 14, 293–313. doi: 10.1177/1471082X13504721
Nardini, A., Tyree, M. T., and Salleo, S. (2001). Xylem cavitation in the leaf of Prunus laurocerasus and its impact on leaf hydraulics 1. Plant Physiol. 125, 1700–1709. doi: 10.1104/pp.125.4.1700
Oliveira, I., Meyer, A., Afonso, S., and Gonçalves, B. (2018). Compared leaf anatomy and water relations of commercial and traditional Prunus dulcis (Mill.) cultivars under rain-fed conditions. Sci. Hortic. 229, 226–232. doi: 10.1016/j.scienta.2017.11.015
Poni, S., Bernizzoni, F., and Civardi, S. (2007). Response of “Sangiovese” grapevines to partial root-zone drying: Gas-exchange, growth and grape composition. Sci. Hortic. 114, 96–103. doi: 10.1016/J.SCIENTA.2007.06.003
Pou, A., Medrano, H., Tomàs, M., Martorell, S., Ribas-Carbó, M., and Flexas, J. (2012). Anisohydric behaviour in grapevines results in better performance under moderate water stress and recovery than isohydric behaviour. Plant Soil 359, 335–349. doi: 10.1007/s11104-012-1206-7
Prieto, J. A., Lebon, É, and Ojeda, H. (2010). Stomatal behavior of different grapevine cultivars in response to soil water status and air water vapor pressure deficit. OENO One 44, 9–20. doi: 10.20870/oeno-one.2010.44.1.1459
R Core Team (2022). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org
Romero, P., Navarro, J. M., García, F., and Botía Ordaz, P. (2004). Effects of regulated deficit irrigation during the pre-harvest period on gas exchange, leaf development and crop yield of mature almond trees. Tree Physiol. 24, 303–312. doi: 10.1093/TREEPHYS/24.3.303
Rouhi, V., Samson, R., Lemeur, R., and Van Damme, P. (2007). Photosynthetic gas exchange characteristics in three different almond species during drought stress and subsequent recovery. Environ. Exp. Bot. 59, 117–129. doi: 10.1016/j.envexpbot.2005.10.001
Salleo, S., Nardini, A., Pitt, F., and Lo Gullo, M. A. (2000). Xylem cavitation and hydraulic control of stomatal conductance in Laurel (Laurus nobilis L.). Plant Cell Environ. 23, 71–79. doi: 10.1046/J.1365-3040.2000.00516.X
Schär, C., Vidale, P. L., Lüthi, D., Frei, C., Häberli, C., Liniger, M. A., et al. (2004). The role of increasing temperature variability in European summer heatwaves. Nature 427, 332–336. doi: 10.1038/nature02300
Scholander, P. F., Hammel, H. T., Bradstreet, E. D., and Hemmingsen, E. A. (1965). Sap pressure in vascular plants. Science 148, 339–346. doi: 10.1126/SCIENCE.148.3668.339
Schultz, H. R. (1996). Water relations and photosynthetic responses of two grapevine cultivars of different geographical origin during water stress. Acta Hortic. 427, 251–266. doi: 10.17660/ACTAHORTIC.1996.427.30
Schultz, H. R. (2003). Differences in hydraulic architecture account for near-isohydric and anisohydric behaviour of two field-grown Vitis vinifera L. cultivars during drought. Plant Cell Environ. 26, 1393–1405.
Scoffoni, C., McKown, A. D., Rawls, M., and Sack, L. (2012). Dynamics of leaf hydraulic conductance with water status: Quantification and analysis of species differences under steady state. J. Exp. Bot. 63, 643–658. doi: 10.1093/jxb/err270
Skelton, R. P., West, A. G., and Dawson, T. E. (2015). Predicting plant vulnerability to drought in biodiverse regions using functional traits. Proc. Natl. Acad. Sci. U.S.A. 112, 5744–5749. doi: 10.1073/pnas.1503376112/-/DCSupplemental.www.pnas.org/cgi//10.1073/pnas.1503376112
Soar, C. J., Speirs, J., Maffei, S. M., Penrose, A. B., McCarthy, M. G., and Loveys, B. R. (2006). Grape vine varieties Shiraz and Grenache differ in their stomatal response to VPD: Apparent links with ABA physiology and gene expression in leaf tissue. Aust. J. Grape Wine Res. 12, 2–12. doi: 10.1111/j.1755-0238.2006.tb00038.x
Tardieu, F., Lafarge, T., and Simonneau, T. (1996). Stomatal control by fed or endogenous xylem ABA in sunflower: Interpretation of correlations between leaf water potential and stomatal conductance in anisohydric species. Plant Cell Environ. 19, 75–84. doi: 10.1111/j.1365-3040.1996.tb00228.x
Tardieu, F., and Simonneau, T. (1998). Variability among species of stomatal control under fluctuating soil water status and evaporative demand: Modelling isohydric and anisohydric behaviours. J. Exp. Bot. 49, 419–432. doi: 10.1093/jxb/49.special_issue.419
Tombesi, S., Nardini, A., Farinelli, D., and Palliotti, A. (2014). Relationships between stomatal behavior, xylem vulnerability to cavitation and leaf water relations in two cultivars of Vitis vinifera. Physiol. Plant. 152, 453–464. doi: 10.1111/ppl.12180
Torrecillas, A., Alarcón, J. J., Domingo, R., Planes, J., and Sánchez-Blanco, M. J. (1996). Strategies for drought resistance in leaves of two almond cultivars. Plant Sci. 118, 135–143. doi: 10.1016/0168-9452(96)04434-2
Tramontini, S., Döring, J., Vitali, M., Ferrandino, A., Stoll, M., and Lovisolo, C. (2014). Soil water-holding capacity mediates hydraulic and hormonal signals of near-isohydric and near-anisohydric Vitis cultivars in potted grapevines. Funct. Plant Biol. 1119–1128. doi: 10.1071/FP13263
Trifilò, P., Nardini, A., Gullo, M. A. L., Barbera, P. M., Savi, T., and Raimondo, F. (2015). Diurnal changes in embolism rate in nine dry forest trees: Relationships with species-specific xylem vulnerability, hydraulic strategy and wood traits. Tree Physiol. 35, 694–705. doi: 10.1093/treephys/tpv049
Tyree, M. T., and Hammel, H. T. (1972). The measurement of the turgor pressure and the water relations of plants by the pressure-bomb technique. J. Exp. Bot. 23, 267–282. doi: 10.1093/JXB/23.1.267
Yadollahi, A., Arzani, K., Ebadi, A., Wirthensohn, M., and Karimi, S. (2011). The response of different almond genotypes to moderate and severe water stress in order to screen for drought tolerance. Sci. Hortic. 129, 403–413. doi: 10.1016/j.scienta.2011.04.007
Keywords: drought, functional traits, xylem vulnerability to cavitation, hydroscapes, leaf water potential, stomatal conductance
Citation: Álvarez-Maldini C, Acevedo M, Estay D, Aros F, Dumroese RK, Sandoval S and Pinto M (2022) Examining physiological, water relations, and hydraulic vulnerability traits to determine anisohydric and isohydric behavior in almond (Prunus dulcis) cultivars: Implications for selecting agronomic cultivars under changing climate. Front. Plant Sci. 13:974050. doi: 10.3389/fpls.2022.974050
Received: 20 June 2022; Accepted: 04 August 2022;
Published: 25 August 2022.
Edited by:
Raul Antonio Sperotto, Universidade do Vale do Taquari - Univates, BrazilReviewed by:
Duncan Smith, University of Wisconsin-Madison, United StatesValmor João Bianchi, Federal University of Pelotas, Brazil
Tom Gradziel, University of California, Davis, United States
Copyright © 2022 Álvarez-Maldini, Acevedo, Estay, Aros, Dumroese, Sandoval and Pinto. 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: Carolina Álvarez-Maldini, carolina.alvarez@uoh.cl