Seeds of isogenic tomato cv. AC and its ABA-deficient tomato mutant (flacca) (Solanum lycopersicum) were provided by the Lancaster Environment Centre (Lancaster University, United Kingdom) and grown in climate-controlled greenhouses at the Faculty of Science, University of Copenhagen, Taastrup, Denmark. Due to the impairment in the oxidation of ABA-aldehyde to ABA, flacca had lower endogenous ABA concentrations compared to AC (Sagi et al., 2002). At the 4-leaf stage, the seedlings were transplanted to 4 L pots filled with 2.2 kg of peat material (Plugg-och Såjord-Dry matter ca.110 kg m^–3, organic matter >95%, pH 5.5–6.5 and EC 1.5–2.5 mS cm^–1). In total, 4 weeks after transplanting, perlite was used to cover the soil surface to minimize soil evaporation, and fertilizers as (NH₄)₂SO₄ (2.6 g) and H₂KPO₄ (1.5 g) per pot were added together with irrigation water to each pot to avoid nutrient deficiency.

From sowing, the plants were grown in two greenhouse cells with CO₂ concentrations of 400 ppm (ambient CO₂, a[CO₂]) and 800 ppm (elevated CO₂, e[CO₂]), respectively. The [CO₂] in the cells was sustained by pure CO₂ emission from a bottle tank and distributed evenly by the internal ventilation system. The [CO₂] in the cells was monitored every 6 s by a CO₂ Transmitter Series GMT220 (Vaisala Group, Helsinki, Finland). The average daily [CO₂] in each cell during treatments are shown in Figure 1. The climate conditions in the two glasshouse cells were set at: 20/16 ± 2°C day/night air temperature, 60 ± 2% relative humidity (RH), 16 h photoperiod and >300 μmol m^–2 s^–1 photosynthetically active radiation (PAR) supplied by sunlight plus LED lamps (Philips GreenPower LED toplighting, Denmark). The average temperature, RH, vapor pressure deficiency (VPD), and daily [CO₂] in the cells during the experiment are shown in Supplementary Figure 1.

The sowing and growth of AC and flacca were conducted consecutively, and all treatments for the two genotypes were carried out separately. In each cell and genotype, 40 plants were established and constantly irrigated to 90% of the pot-holding capacity. Exogenous ABA application started 20 days after transplanting for three successive days. On each day, half of the plants were sprayed with 20 μM ABA (Sigma) twice a day on the whole plant at 10:00 and 15:00 h, respectively; the other half was sprayed with deionized water as non-ABA-treated control. All the solutions contained 0.1% (v/v) Tween-20 and 0.1% (v/v) ethanol.

In total, 1 day after finishing exogenous ABA application, four plants of each treatment (eight treatments in total) were harvested as initial control, then progressive soil drying was commenced. During soil drying, for each genotype, half of the plants were well watered to 95% of the pot’s water-holding capacity, and the remaining half was subjected to progressive soil drying by withholding irrigation from pots until all the transpirable soil water was depleted. Soil water content in the pot was expressed as the fraction of transpirable soil water (FTSW) (see below). In addition to the first harvest as initial control when FTSW was ca. 0.95, there were two destructive harvests during the soil drying period: the second was conducted when FTSW was ca. 0.50, and the last harvest was conducted at the end of soil drying when FTSW below 0.1. A total of four biological replicates (four individual plants) for each treatment were harvested. Due to the different water consumption rates between the two genotypes and the two [CO₂] levels, AC grown at a[CO₂] and e[CO₂] spent 10 and 9 days to reach the end of drought stress, respectively; flacca grown at both [CO₂] environments spent 13 days. Therefore, plants among different treatments were harvested on different days but at the same three FTSW points. The experimental treatments for AC and flacca were successively conducted, this might have caused the slightly different climatic conditions (see Supplementary Figure 1) between the two genotypes.

Soil water status was measured daily by weighing the pots with an Analytical Balance (Sartorius Model QA35EDE-S) at 15:30 h and expressed as FTSW. The daily value of FTSW was estimated as the ratio between the amounts of transpirable soil water remaining in the pots and the total transpirable soil water (TTSW). TTSW was defined as the difference of pot weight between full water holding capacity and when all the transpirable soil water was depleted, and calculated as:

where WT_f is the pot weight at full water holding capacity (ca. 3.6 kg)., and WT_e is the pot weight at the end of soil drying (ca. 1.8 kg).

Compared to our previous study (Wei et al., 2020), here the definition of FTSW was different which resulted in the longer duration for flacca plants to deplete the TTSW in relation to AC plants. The reason behind this was that exogenous ABA application would affect the stomatal conductance (g_s) for both genotypes, thus using a common TTSW value would allow us to compare the g_s response among the treatments on the same base. Therefore, for all treatments, the FTSW was calculated as follows:

where WT_n is the pot weight on a given date, and WT_e is pot weight at the end of soil drying.

A total of 1 day after exogenous ABA application (before soil drying), stomatal pore aperture (SA, μm²) was determined on four replicates (10 stomata for one replicates) following the method described by Yan et al. (2012). Stomata were observed under a LEITZ DMRD microscope camera system (Leica Microscope and System GmbH, D 35530, Wetzlar, Germany) equipped with a digital camera. Stomatal pore aperture length (W_a) and pore aperture width (L_a) were measured by ImageJ software [Version 1.51k, Wayne Rasband, National Institutes of Health, United States, Java 1.6.0–24 (64 bit)]. Then SA was calculated as: (π × W_a × L_a)/4.

During progressive soil drying, g_s (mol m^–2 s^–1) were measured on the last youngest upper canopy fully expanded leaves between 9:00 and 12:00 h with a portable photosynthetic system (LiCor-6400XT, LI-Cor, Lincoln, NE, United States). Measurements were performed on one leaf per plant and four biological replicates for each treatment at 22°C cuvette temperature, 1500 μmol m^–2 s^–1 PAR, and [CO₂] of 400 ppm for a[CO₂] and 800 ppm for e[CO₂] growth environments, respectively.

At each harvest, after gas exchange measurements, the same leaf was excised for the determination of midday leaf water potential (Ψ_leaf, MPa) with a scholander-type pressure chamber (Soil Moisture Equipment Corp., Santa Barbara, CA, United States) following the method described by Liu et al. (2006). Then the excised leaves were cut into two parts, frozen in liquid nitrogen, and stored at −80°C for later determination of osmotic potential (Ψ_π, MPa) and leaf ABA concentration ([ABA]_leaf, μg g^–1 FW). Ψ_π was measured at 20°C with a psychrometer (C-52 sample chambers, Wescor Inc., Logan, UT, United States) connected to a microvoltmeter (HR-33T, Wescor Inc., Logan, UT, United States). Turgor pressure (Ψ_p, MPa) was then calculated as follows: Ψ_leaf–Ψ_π. Leaf dry weight (LDW) (g) and stem dry weight (SDW) (g) were determined at each harvest.

At each harvest, xylem saps were collected with a scholander-type pressure chamber (AGRSCI, KVL, Denmark) according to Liu et al. (2006), then stored at −80°C for determination of xylem ABA concentrations ([ABA]_xylem, pmol ml^–1). [ABA]_leaf and [ABA]_xylem were determined by Enzyme-linked immunosorbent assay following the protocol of Asch (2000).

The responses of g_s to soil drying was described by a linear-plateau model (Faralli et al., 2019):

where g_{s ini} denotes initial g_s; C denotes the FTSW threshold at which y started to diverge from g_{s ini}; and A was the slope of the linear equation.

For ABA-treated flacca plants, before FTSW declining to the threshold, there was a significant increasing trend of g_s during progressive soil drying. Therefore, when FTSW > C, the linear-plateau model was modified as follows:

where g_{s max} denotes maximum g_s when FTSW declined to the threshold (C), and D was the slope of the linear equation before the threshold. g_{s max} of non-ABA-treated flacca was obtained from the average of the individual replicated values when FTSW reached C. In addition, for all flacca plants, g_{s end} indicated the g_s at the end of drought treatment.

The parameters y, C, A, and D were estimated by PROC NLIN of PC SAS 9.4 (SAS Institute Inc., Cary, NC, United States, 2002–2012) and the coefficient of determination (r²) was calculated. Statistical comparison of each parameter obtained from the linear-plateau regression between treatments was performed by t-test using MedCalc statistical software 19.0.7.

Data were statistically analyzed using Microsoft Excel, SAS 9.4 (SAS Institute Inc., Cary, NC, United States, 2002–2012), SPSS 22.0 software (IBM SPSS Software, Armonk, NY, United States), and GraphPad Prism 9 software. Two-way analyses of variance (ANOVA) were performed to analyze the effects of [CO₂] and exogenous ABA ([ABA]) on plant dry weight, SA, water relation variables, and ABA concentrations in AC and flacca. To compare the decreasing trends of FTSW during soil drying between control treatment and exogenous ABA application, the statistical differences were analyzed by analysis of covariance (ANCOVA, days after withholding water as a covariate). One-way ANOVA (Tukey’s test) and Student’s t-test were conducted to determine the significant differences between treatments. Principle component analysis (PCA) of g_s, water relations, and leaf and xylem ABA concentrations in AC and flacca were performed in R version 4.0.0 (R Core Team, 2020). The relationship between g_s and [ABA]_leaf, g_s, and Ψ_p were evaluated by linear regression. r² of the regression lines were calculated and statistical differences on the slopes and intercepts of regression lines between treatments were performed by ANCOVA.

In AC, exogenous ABA application did not have a significant influence on LDW and SDW, while e[CO₂] increased both LDW and SDW. In flacca, the ABA-treated plants had higher LDW and SDW when FTSW = 0.50 and 0.05 and did not respond to e[CO₂] except for an increased SDW at the last harvest (Table 1).

After withholding irrigation, the FTSW in the pots of drought-stressed plants decreased steadily (Figure 1). In AC, only for plants grown at e[CO₂], exogenous ABA slowed the soil water depletion rate during progressive soil drying. In flacca, at both [CO₂] levels, exogenous ABA had a significant influence on the rates of soil water depletion. However, from 0 to 10 days, FTSW declined slower in the ABA-treated plants compared with the non-ABA-treated plants, then it declined faster on the last 3 days.

In both AC and flacca, e[CO₂] decreased stomatal pore aperture (SA), though it was more pronounced in AC (Figure 2). One day after ABA application (before imposing soil drying), exogenous ABA decreased SA in both genotypes, and interactions of [CO₂] and [ABA] were observed, indicating that the effect of e[CO₂] on SA was eliminated by exogenous ABA application, and the depressive effect of exogenous ABA was more pronounced in flacca than in AC.

In AC, during progressive soil drying, plants grown at e[CO₂] had significantly lower initial stomatal conductance (g_{s ini}) compared with their a[CO₂]-grown counterparts, and g_s started to decline when FTSW dropped to an FTSW threshold (C). The ABA-treated plants had constantly lower g_{s ini} than the non-ABA-treated pants, and C was also slightly lowered by exogenous ABA application (Figures 3A,B and Table 2).

In flacca, 1 day after exogenous ABA application, depression on g_{s ini} ranged from 50 to 58% at a[CO₂] and e[CO₂] environment, respectively. During soil drying, there were clear increasing trends of g_s in the ABA-treated plants before FTSW reaching to the threshold, and when FTSW decreased to C, only the ABA-treated plants grown at a[CO₂] still had lower maximum g_s (g_{s max}) than the non-ABA-treated plants. Moreover, the FTSW thresholds were slightly advanced by exogenous ABA application. At the end of soil drying, exogenous ABA treatment still lowered g_s (g_{s end}) of plants grown at both [CO₂] levels, but all flacca plants possessed significantly higher g_{s end} than AC. In addition, e[CO₂] decreased g_s in non-ABA-treated flacca and slightly advanced C (Figures 3C,D and Table 2).

A total of 1 day after exogenous ABA application, before imposing soil drying (i.e., FTSW = 0.95), the ABA-treated AC plants had higher leaf water potential (Ψ_leaf) and osmotic potential (Ψ_π), whereas the ABA-treated flacca plants possessed higher Ψ_leaf and turgor pressure (Ψ_p). In addition, slight interactions of [CO₂] and [ABA] were observed in Ψ_leaf and Ψ_p of AC when FTSW = 0.95, showing that the effects of exogenous ABA became less significant when plants were grown at e[CO₂]. At the end of soil drying (i.e., FTSW = 0.05), all these plant water relation variables significantly decreased by drought stress, and exogenous ABA had no influence on water relations of flacca, whereas it still slightly increased Ψ_leaf in AC. Moreover, e[CO₂] improved Ψ_leaf, Ψ_π, and Ψ_p of AC plants compared to a[CO₂], though these effects being diminished along with soil drying (Figure 4 and Table 3).

In total, 1 day after exogenous ABA application, the ABA-treated AC plants possessed higher leaf ABA concentration [ABA]_leaf than the non-ABA-treated plants, while the ABA-treated flacca had both higher [ABA]_leaf and xylem sap ABA concentration ([ABA]_xylem) compared with the non-ABA-treated controls. During soil drying, when FTSW = 0.50, there were interactions of [ABA] and [CO₂] on [ABA]_leaf of both genotypes, which was relatively higher in the ABA-treated plans grown under e[CO₂] than under a[CO₂]. Moreover, at the end of soil drying (i.e., FTSW = 0.05), significantly higher [ABA]_xylem was observed in the ABA-treated flacca than the non-ABA-treated flacca plants, accompanied by an interaction of [CO₂] and [ABA]. Drought stress increased both [ABA]_leaf and [ABA]_xylem at the end of soil drying though being less significant in flacca (Figure 5 and Table 3).

The relationships of the FTSW threshold at which g_s started to decline (C) and [ABA]_leaf at the second harvest (FTSW value was close to C) in AC and flacca are shown in Figure 6. The results showed that only in flacca, C was linearly correlated with [ABA]_leaf, whereas C of AC did not respond to increasing [ABA]_leaf. Moreover, no obvious relationship between C and [ABA]_xylem was observed in either genotype.

Principle component analysis plot of g_s, water relations, leaf, and xylem ABA concentrations of well-watered AC and flacca is depicted in Figure 7. PC1 and PC2 axes explained 70.3% of cluster formation, with 49.1% attributed to PC1 and 21.2% to PC2. Overall, AC under all treatments and ABA-treated flacca were clustered toward the right side of the PCA plot, whereas all non-ABA-treated flacca were clustered to the left in the same direction as the g_s vector, indicating non-ABA-treated flacca had higher g_s than other plants. The clustering of the [AC, 800 ppm, +ABA] indicated a covariation between Ψ_leaf and [ABA]_leaf under this treatment, which was negatively associated with g_s. [ABA, flacca] and [AC, 800 ppm, CON] treatments were grouped along with PC2, showing a close correlation with [ABA]_xylem and Ψ_p, and no correlation with g_s. Furthermore, [AC, 400 ppm] was opposed by Ψ_p and [ABA]_xylem, being clustered to the bottom of the plot in the same direction as Ψ_π vector.

In AC, there were linear relationships between g_s and [ABA]_leaf, g_s and Ψ_π, where g_s declined linearly with increasing [ABA]_leaf and decreasing Ψ_π (Figure 8). Exogenous ABA application decreased the slope of the regression line of g_s and [ABA]_leaf, while e[CO₂] had no influence on it. Regarding the regression line of g_s and Ψ_π, the slopes differed between ABA-treated and non-ABA-treated plants, whereas e[CO₂] only decreased the slope in non-ABA-treated plants. In flacca, there was no significant relationship between g_s and [ABA]_leaf under either treatment, while g_s decreased linearly with declining Ψ_π. In addition, the intercepts of these regression lines differed between [ABA] treatments due to the initial low levels of g_s.

It is well known that e[CO₂] decreases g_s while stimulating carbon assimilation, thus improving plant drought tolerance (Li S. et al., 2020). In tomato plants, previous studies have shown that the plant response to e[CO₂] was attenuated or absent in ABA-deficient mutant flacca, affirming the obligatory role of ABA in CO₂ signaling pathway (Fang et al., 2019; Wei et al., 2020). Consistent with this, here LDW and SDW of AC plants were improved by e[CO₂] during progressive soil drying, but these stimulations were less significant in flacca. Upon exogenous ABA application, retarded plant growth of flacca was rescued under both watering conditions, while growth of AC plants did not benefit from exogenous ABA treatment (Table 1), which was consistent with previous study on tomato and confirmed the important role of ABA in maintaining the shoot development (Sharp et al., 2000; Aroca et al., 2008).

It has been reported that application of exogenous ABA inhibited water loss and reduced the rate of soil drying due to stomatal closure (Hossain et al., 2015; He et al., 2019). In the present study, in both ABA-treated AC and flacca, slower soil water depletion during progressive soil drying was recorded, which coincided with the decreases of SA and g_s in relation to the non-ABA-treated controls (Figures 1–3). However, in the drought-stressed AC plants, only those grown under e[CO₂] possessed slower rate of soil water depletion upon exogenous ABA application compared with the non ABA-treated plants, indicating a significant interaction of e[CO₂] and exogenous ABA on plant water consumption when soil water was limited. Although in AC e[CO₂] accelerated water consumption due to stimulated plant growth, exogenous ABA application had no influence on plant growth (Figure 1 and Table 1). Therefore, the interactive effect of e[CO₂] and exogenous ABA on water consumption rate might be ascribed to the altered stomatal behavior, as both factors could induce stomatal closure.

It is worth noting that upon exogenous ABA treatment the depressive effect of e[CO₂] on SA was eliminated in both genotypes (Figure 2). Considering that e[CO₂] had exerted inhibiting effects on stomatal aperture, it could be assumed that exogenous ABA application overrode the effect of e[CO₂] on stomatal movement. However, the common view is that ABA can amplify the effects of [CO₂] on stomatal behavior (Engineer et al., 2016), and e[CO₂]-grown soybean showed a stronger response to endogenous ABA (Gray et al., 2016), though ABA and [CO₂] signal transductions also could be independent (Hsu et al., 2018). As our previous studies demonstrated that severe drought stress overrode the e[CO₂] effects, here exogenous ABA application exhibited the same influences as drought stress. Furthermore, decreases in SA and g_s were more pronounced in flacca than AC (Figures 2, 3). Namely, the stomata of flacca exhibited greater sensitivity to exogenous ABA than AC due to ABA deficit. Interestingly, after the onset of soil drying, before the occurrence of stomatal closure, there were significant increasing trends of g_s in the ABA-treated flacca while the ABA-treated AC possessed a constantly lower g_s before the FTSW thresholds (Figure 3), indicating that the effects of exogenous ABA application lasted for a longer period in AC than in flacca. Early studies have found that the abnormal stomatal behavior in flacca was associated with the high level of ethylene resulting from an ABA deficit (Tal et al., 1979), and other chemical signals could also counteract the ABA signaling, as is the case with cytokinin (Prerostova et al., 2018), which might accelerate the degradation process of exogenous ABA in flacca plants. In addition, flacca is known to be impaired in the oxidation of ABA aldehyde to ABA (Sagi et al., 2002), vanishing effects of exogenous ABA might be ascribed to the redistribution of ABA to other plant tissues. These results suggest the necessity to explore the regulation of ABA homeostasis in the two genotypes in future studies. Despite the significantly lower g_{s ini} in ABA-treated flacca, the g_{s end} of flacca was still higher than AC under all treatments (Figure 3 and Table 2), which indicated that flacca could not fully close stomata under severe drought stress after exogenous ABA application. In another ABA-deficient mutant sitiens, a decrease in the leaf gas exchange rate by exogenous ABA (100 μM) was significant after 52 days of mild drought stress, but it still possessed a higher level of transpiration rate than wild type (Aroca et al., 2008). Therefore, exogenous ABA priming before soil drying could not sufficiently induce stomatal closure under severe drought stress.

Previous studies on tomato revealed that e[CO₂] retarded stomatal closure during soil drying due to reduced stomatal sensitivity to ABA (Liu et al., 2019; Wei et al., 2020). However, here we did not find e[CO₂] delayed stomatal closure in AC, the reason behind this discrepancy is unknown, which might be ascribed to the varied experimental conditions between the different studies. By contrast, in flacca plants, exogenous ABA application sensitized the stomatal response to drought though only significant under a[CO₂] (Figures 3C,D and Table 2), suggesting that the stomatal closure in ABA-deficient mutant could be advanced by exogenous ABA treatment, which might contribute to water-saving under drought stress. In addition, e[CO₂] decreased g_s in flacca despite the slight influence on SA, and these effects were absent in ABA-treated flacca. Although the retarded stomatal response to e[CO₂] in flacca was attributed to ABA deficit (Wei et al., 2020; Li et al., 2021), exogenous ABA application could not recover the response due to its strongly induced stomatal closure.

Previous studies have suggested that exogenous ABA application is beneficial for plants to improve drought tolerance as a result of induced stomatal closure, enhanced water relations, and a higher OA (Du et al., 2013; He et al., 2019). In the present study, 1 day after exogenous ABA application, higher Ψ_leaf, Ψ_π, and Ψ_p were observed compared to the non-ABA-treated plants (Figure 4), confirming the positive effects of exogenous ABA on plant water relations. In an early study with exogenous ABA application (10 μM) on flacca, the authors found that the impaired shoot growth could be restored but no improved Ψ_leaf was observed (Sharp et al., 2000). As stimulation on expansive growth by ABA is primarily related to hydraulic control (Tardieu et al., 2015), here the improved water status in ABA-treated flacca might be responsible for the restored growth (Table 1). Moreover, the positive effects of exogenous ABA on plant water relations were more pronounced in flacca than in AC (Figure 4), coinciding with the stomatal responses in the two genotypes. In wild-type tomato and ABA-deficient mutant sitiens, Aroca et al. (2008) found that the genotypes showed different gene regulatory patterns in response to exogenous ABA, including ABA biosynthesis-related and aquaporin-encoding genes, which might be associated with the high sensitivity to ABA in flacca plants. Surprisingly, at the end of soil drying, the improved water relations, by application of exogenous ABA, vanished in both genotypes, including Ψ_p (Figure 4 and Table 3). Thereby, the decreased g_s in ABA-treated flacca at the end of soil drying could not be explained by modulated water relations. In addition, OA is also considered as a vital trait for maintaining leaf turgor, and increased OA accompanied by improved leaf water relations have also been reported in exogenous ABA-treated soybean (He et al., 2019). However, here we did not find the influence of exogenous ABA on OA in both genotypes during the whole period of soil drying. Therefore, more detailed mechanisms related to the ABA-modulated plant water relations merit future studies.

The positive effects of e[CO₂] on plant growth and physiology are reported to be more pronounced under drought stress which, however, depends on the severity of drought stress (Leakey et al., 2006; Gray et al., 2016). In the present study, the well-watered AC plants grown at e[CO₂] had an improved Ψ_leaf, Ψ_π, and Ψ_p in relation to those grown under a[CO₂], but these effects were eliminated during soil drying (Figure 4 and Table 3), in line with previous findings (Wei et al., 2020). As e[CO₂]-grown plant already possessed improved water relations, effects of exogenous ABA on AC became less significant at e[CO₂]. In Figure 8, e[CO₂] rendered a less sensitivity of stomata to Ψ_π, but this effect was absent in ABA-treated AC. Collectively, these results implied an offsetting effect of [CO₂] and exogenous ABA on plant water relations as well as stomatal behavior.

In the present study, 1 day after exogenous ABA application, significant increases of [ABA]_leaf in AC and flacca were observed, and the increase was more pronounced in flacca than in AC, though it was gradually diminished during soil drying (Figure 5 and Table 3). As oxidative degradation of ABA could occur rapidly (Zhang et al., 1995), in tomato plants, when ABA-treated plants were exposed to long-term drought (52 days), ABA contents in plants did not differ between the ABA-treated and non-ABA-treated plants despite the accumulated exogenous ABA effects on water relations still existing (Aroca et al., 2008). By contrast, here, we found that the vanishing effect on water relations was accompanied by decreased accumulation of [ABA]_leaf in ABA-treated flacca. More interestingly, the interactions of [CO₂] and [ABA] on [ABA]_leaf of AC and flacca indicated that the ABA degradation process might be disturbed under e[CO₂], which could explain the reduced ABA sensitivity in the drought-stressed tomato plants gown under e[CO₂] (Yan et al., 2017; Wei et al., 2020). In ABA-treated flacca, an increase in [ABA]_xylem was observed especially at the end of soil drying, which was accompanied by an interaction of [CO₂] and exogenous ABA (Figure 5 and Table 3), implying that e[CO₂] exerted an effect on ABA homeostasis. It is well documented that ABA metabolism, including biosynthesis and degradation, ABA recirculation, and exudation processes are all involved in the regulation of ABA homeostasis (Hartung et al., 2005). Moreover, redistribution of ABA within plant organs has often been reported, and translocation of ABA from leaves to roots can be intensified when roots are experiencing drought (Hartung et al., 2005; Ikegami et al., 2009; Ernst et al., 2010). Taking together, except for the ABA degradation pathway, decrease of [ABA]_leaf in the ABA-treated flacca could also be ascribed to the redistribution of ABA between shoot and root considering its inability of ABA synthesis in roots.

It is well known that tomato plants are likely to exhibit isohydric behaviors (Moshelion et al., 2015), while ABA-deficient mutant flacca and sitiens could be transpired more for longer periods during soil drying than AC (Aroca et al., 2008; Wei et al., 2020), representing a near-anisohydric characteristic. Furthermore, rapid ABA biosynthesis might facilitate isohydric behavior (Moshelion et al., 2015). In the present study, before the occurrence of stomatal closure (i.e., when FTSW = 0.50), ABA-treated flacca grown at e[CO₂] possessed relatively high [ABA]_leaf (Figure 5), which was found to be closely correlated with the FTSW threshold at which g_s started to decrease (C) (Figure 6). These results raise the possibility that anisohydric behaviors in the ABA-deficient mutant could be shifted toward near-isohydric behaviors through the combined effects of e[CO₂] and exogenous ABA, representing a restored stomatal drought response. The change in the water management strategies has been reported to be activated by different soil water statuses (Zhang et al., 2012), also the overexpression of a key tonoplast aquaporin (Sade et al., 2009), which could be linked to ABA.

Previous studies have shown that shoot-to-root ABA transport played a role in the regulation of water flux and induction of stress-resistant genes in roots (Ikegami et al., 2009). In the near-anisohydric grapevine cultivar Syrah, the authors found that it showed relatively high catabolism of ABA in xylem sap, which might lower its hypersensitivity to water stress (Dayer et al., 2020). Therefore, it could be assumed that reloading of ABA into xylem sap in flacca plants might stimulate some specific metabolic process, causing feedback on stomatal movements, which might explain the low g_s in ABA-treated plants at the end of soil drying but without [ABA]_leaf accumulation. For example, early studies on flacca have shown that low ABA content stimulated ethylene production, which could be restored to normal levels with exogenous ABA (Tal et al., 1979; Sharp et al., 2000). As ethylene production is often increased by drought stress, it might be associated with the restored drought response in ABA-treated flacca as noticed in the present study. In addition, exogenous ABA-induced other metabolic adaption has been reported widely, including stimulated resistant protein patterns (Zhou et al., 2014), enhanced energy storages, and activities of antioxidant enzymes (Latif, 2014; Gai et al., 2020), and these possibilities merit further studies.

Under well-watered conditions, the ABA-treated flacca plants possessed the same leaf characteristics as AC in the global PCA plots (Figure 7). However, despite the restored leaf gas exchange rates, water relation characteristics, and ABA concentrations in flacca (Figures 3–5), the PCA plot indicated that, under well-watered conditions, exogenous ABA could still not induce a correlation between g_s and [ABA]_leaf in flacca, but showed an interaction with Ψ_p and [ABA]_xylem (Figure 7). Furthermore, during the whole period of soil drying, g_s was still not correlated with [ABA]_leaf but only responded to Ψ_p (Figure 8). Therefore, a question arises as to how exogenous ABA application affected flacca’s stomatal movements considering the no involvement of [ABA]_leaf. Pantin et al. (2013) demonstrated that ABA can induce stomatal closure through an indirect hydraulic effect on water permeability within leaf vascular tissues. In the exogenous ABA-fed near-anisohydric grapevine cultivar, the existence of indirect ABA effects on g_s has also been reported, which was associated with the ABA catabolism (Dayer et al., 2020). Therefore, it could be assumed that the exogenous ABA decreased g_s in flacca plants through an indirect hydraulic effect. Taking together, future research should focus on exploring the correlation of ABA metabolic process and ABA indirect hydraulic effects on stomatal behavior.