Raf-like kinases and receptor-like (pseudo)kinase GHR1 are required for stomatal vapor pressure difference response

Significance With the continuing increase in global temperatures, plants transpire more water due to the increasing vapor pressure deficit. Stomatal pores in plants close rapidly in response to the rising vapor pressure deficit to counteract water loss. We demonstrate that mutations in the stomatal CO2 signaling pathway do not impair the response to an increase in vapor pressure difference (VPD). Osmotic stress causes cytoplasmic Ca2+ transients in guard cells. Nevertheless, we show that diverse investigated higher-order calcium-signaling mutants do not affect the VPD response. We reveal that B3 family Raf-like protein kinases and a plasma membrane receptor-like protein GHR1 function in the elusive leaf-to-air VPD-mediated stomatal closure pathway. Notably, ghr1 mutant alleles disrupt the classical “wrong-way” stomatal VPD response.

relative air humidity j vapor pressure deficit j abscisic acid j Raf-like MAP kinase kinase kinases j GUARD CELL HYDROGEN PEROXIDE-RESISTANT1 S tomata are located in the epidermis of aerial tissues of plants and consist of a pair of guard cells that bend apart under increased turgor to open a gas-exchange pore. By adjusting their pore apertures in response to abiotic stresses and stimuli, including water status, temperature, and carbon dioxide (CO 2 ) concentrations, stomata regulate transpiration and CO 2 influx into leaves. By influencing the rate of transpiration and hence leaf water status, changes in relative air humidity trigger the regulation of stomatal apertures (1,2). Lowering the relative humidity in the air surrounding leaves increases the leaf-toair vapor pressure difference (VPD). This increase in VPD, in turn, results in an initially enhanced transpiration rate, but stomata then respond by partially closing, thereby reducing water loss of plants. Two mechanisms can contribute to high-VPDtriggered stomatal closure: 1. A hydropassive mechanism suggests that the passive equilibration of guard-cell water potential with VPD-induced reduction in external water potential reduces guard-cell turgor and thereby the stomatal pore aperture (3)(4)(5)(6)(7)(8). 2. A hydroactive mechanism suggests that the stomatal closing response to increased VPD requires a metabolically and iontransport-driven/guard-cell signal-transduction-driven reduction in guard-cell solute concentration (2,5,7). This involves signal transduction in guard cells, leading to guard-cell ionchannel and transporter regulation, guard-cell ion efflux, and a reduction in guard-cell solutes, resulting in turgor reduction and stomatal closing (2,5).
These two mechanisms do not exclude one another, and both likely contribute to stomatal closing, depending on the species and conditions (2,5), with the hydroactive mechanism considered to be of greater importance in angiosperm species with an epidermal mechanical advantage over guard cells, notably in grasses (9)(10)(11)(12). However, the molecular mechanisms of VPD sensing and VPD-induced signal transduction during the hydroactive response remain to a large degree unknown.

Significance
With the continuing increase in global temperatures, plants transpire more water due to the increasing vapor pressure deficit. Stomatal pores in plants close rapidly in response to the rising vapor pressure deficit to counteract water loss. We demonstrate that mutations in the stomatal CO 2 signaling pathway do not impair the response to an increase in vapor pressure difference (VPD). Osmotic stress causes cytoplasmic Ca 2+ transients in guard cells. Nevertheless, we show that diverse investigated higher-order calcium-signaling mutants do not affect the VPD response. We reveal that B3 family Raf-like protein kinases and a plasma membrane receptor-like protein GHR1 function in the elusive leaf-to-air VPD-mediated stomatal closure pathway. Notably, ghr1 mutant alleles disrupt the classical "wrong-way" stomatal VPD response.
Slow (S-type) anion channels (13) and rapid (R-type) anion channels (14) play important roles in mediating and regulating guard-cell ion efflux and stomatal closure in response to upstream signaling events. These anion channels mediate anion efflux and cause plasma-membrane depolarization in guard cells (13,14). The S-type anion channel SLAC1 and the R-type anion channel ALMT12/QUAC1 function in stomatal closure (15)(16)(17)(18). Mutants in the SLAC1 gene show impairment in stomatal closure in response to several environmental stimuli, including abscisic acid (ABA), CO 2 , ozone, and low relative humidity (17,18).
The plant hormone ABA plays a central role in regulating stomatal closure (19)(20)(21). ABA-mediated stomatal regulation evolved 400 million years ago in moss species (22). ABA triggers stomatal closure through direct binding to PYR/PYL/RCAR receptors, which, in turn, leads to inhibition of clade A type 2C protein phosphatases (PP2Cs) (23)(24)(25). Inhibition of PP2Cs contributes to the activation of OST1/SnRK2.6 protein kinase (26,27), and then activated OST1/SnRK2.6 in guard cells stimulates S-type anion channels encoded by SLAC1 (28)(29)(30)(31)(32). The role of ABA in high-VPD-mediated hydroactive stomatal closure has been investigated in several ABA biosynthesis mutants and ABA signaling mutants in Arabidopsis. The resulting models have been, to a degree, debated. Early studies on ABA-deficient (aba1-1) and ABAinsensitive (abi1-1 and abi2-1) mutants suggested that ABA and ABA signaling are not necessary for high VPD-induced stomatal closure (33). This is supported by the investigation of transgenic lines expressing a guard-cell-targeted dominant abi1-1 protein phosphatase, which showed a wild-type (WT)-like stomatal response to a rapid increase in VPD, including in time-resolved gas-exchange analyses (34). Overexpression of the Arabidopsis abi1-1 mutant allele in gray poplar, Populus x canescens (Aiton) Sm., exhibited a defect in steady-state stomatal conductance in response to a VPD increase (35). In another study, an ABA biosynthesis mutant (aba2-13) and an OST1 protein kinase mutant (ost1-4) were identified by genetic screening for mutants defective in low-humidity-mediated steady-state leaf temperature regulation and thus were impaired in low-humidity-mediated stomatal closure (36). It has been proposed that high VPD/low humidity rapidly triggers ABA biosynthesis to induce stomatal closure (37,38). However, whole-plant gas-exchange analyses showed that stomata close in response to high-VPD stimulation in several ABAdeficient mutants defective in key steps of ABA biosynthesis, including the most severe nced3/nced5 double mutant (5,39,40), indicating that a direct increase in the ABA concentration of guard cells may not be essential for initiating high VPD-induced stomatal closure. Furthermore, a step change to low humidity (high VPD) caused a similar percentage degree of reduction in the guard-cell K + content in WT and the ABA biosynthesis aba3-1 mutant, indicating a functional high-VPD response in the ABAdeficient mutant (37). Moreover, an ABA receptor, pyr1/pyl1/pyl2/ pyl4/pyl5/pyl8 hextuple mutant, which disrupts ABA-induced stomatal closure (5,41), exhibited a significantly reduced initial response rate and increased half-response time of stomatal closure during high-VPD stimulation (5). These results indicate that ABA signaling may contribute to the stomatal closing rate of the high-VPD response and steady-state stomatal apertures, but rapid enhancement of ABA signaling is not essential for high VPD-induced stomatal closing to be initiated and can proceed, albeit at a slower rate, in the above mutants (5).
It has been suggested that young leaves acquire stomatal ABA sensitivity by pre-exposure to ABA or low relative humidity for 1 d (42). Analyses of ABA concentration increases in response to drought or low humidity have shown that substantial rises in the ABA concentration appear to occur after 3 or more hours of stress treatment (43), whereas high VPD-induced stomatal closing typically occurs within 15 min. Studies have shown that guard cells in nonstressed leaves have measurable and higher ABA concentrations than mesophyll cells (44)(45)(46)(47). Recent studies suggest that steady-state basal ABA signaling occurs in guard cells (46)(47)(48), and a model has been proposed that may explain many of the above findings, in which basal ABA signaling may contribute to the stomatal VPD response by amplifying the hydroactive stomatal closure pathway (20), similar to CO 2 -triggered stomatal closure (46,47). Furthermore, several studies suggest a key role for ABA in adjusting the steady-state stomatal conductance (34,35). In the present study, we pursued analyses of other guard-cell signal-transduction pathways and higher-order mutants given the limited knowledge of early VPD sensing and signaling mechanisms.
In contrast to the findings of a functional high-VPD response in ABA receptor hextuple mutants, the high-VPD response was strongly impaired in the OST1 protein kinase mutant, ost1-3 (5), consistent with earlier research (36). It has been suggested that an ABA-independent pathway is involved in OST1/SnRK2 kinase activation during water deficiency (49,50). However, apart from the requirement of the OST1/SnRK2.6 protein kinase, the molecular and genetic mechanisms that mediate VPD sensing and signal transduction in the rapid hydroactive response remain elusive.
Given the dearth of mutants and mechanisms known to affect the rapid hydroactive VPD-induced stomatal closing response, in the present study, we screened and investigated potential signaling mechanisms to determine mechanisms that mediate the hydroactive stomatal VPD response. We analyzed mutants in several guard-cell signal-transduction pathways to determine which genes, mechanisms, and pathways contribute to the VPD response and which gene combinations alone are not required for the VPD response. The time-dependent kinetics of stomatal VPD responses were investigated by intact leaf gas-exchange analyses in mutants in guard-cell signaling pathways. Differential contributions of the S-type anion channel SLAC1 and the R-type anion channel ALMT12 were found by analyzing the high-VPD response in slac1-3, almt12-1, and slac1-3/almt12-1 mutants. The involvement of osmotic/mechanical sensing and calcium signaling in the stomatal high-VPD response was investigated in the putative osmotic/mechanical senor osca1-2/1.3/2.2/2.3/3.1 quintuple and pathogen-associated molecular patterns-activated calcium channels osca1.3/1.7 double mutants, glutamate receptor-like channels glr3.2/3.3 double mutant, Ca 2+ -permeable cyclic nucleotide-gated channels cngc5/6 double mutant, cngc20 and cngc19/20crispr double mutants, calcium-dependent protein kinase cpk3/5/6/11/23 and cpk3/4/5/6/11 quintuple mutants, and calcineurin-B like proteins cbl1/4/5/8/9 quintuple and cbl2/3rf double mutants. The stomatal high-VPD responses were also analyzed in strongly CO 2insensitive mitogen-activated protein (MAP) kinase mpk12/ mpk4GC double-mutant alleles showing that early CO 2 signaling does not participate in stomatal VPD signal transduction. Interestingly, mutant alleles in the receptor-like (pseudo)kinase GHR1 greatly impaired the stomatal VPD response and impaired the classical "wrong-way" response to stomata to VPD increases. We further investigated whether the recently identified B3-family and B4-family Raf-like MAP Kinase Kinase Kinases (M3Ks) (51), which function in OST1 activation (52)(53)(54), are required for stomatal closure in response to VPD elevation and show that B3-family, but not B4-family, Raf-like M3Ks function in the stomatal VPD response, with a central role for M3Kδ5/RAF6.

Results
High VPD-Induced Stomatal Closure Is Slowed in ost1 Kinase Mutant and Guard-Cell Anion-Channel slac1 Single and slac1/ almt12 Double Mutants. The OST1/SnRK2.6 protein kinase is an important positive transducer of stomatal closing and an activator of anion channels during stomatal closure. It has been shown that ost1/snrk2.6 mutants are impaired in high VPD-induced stomatal closure (5,36,40). In controls, under the imposed experimental conditions, ost1-3 mutant leaves showed a much higher basal stomatal conductance at low VPD (∼0.8 kPa) compared to WT plants (SI Appendix, Fig. S1A). Although a transition to elevated VPD (∼2.0 kPa), caused a time-dependent reduction in stomatal conductance in both ost1-3 mutant and WT leaves, the final stomatal conductances were larger in ost1-3 mutant leaves, and the time course of the response was consistently slower in independent datasets (SI Appendix, Fig. S1A). The initial high-VPD stomatal closing responses were fitted with exponential one-phase decay functions. The average rate constant in ost1-3 mutant leaves was significantly smaller than in WT controls (SI Appendix, Fig.  S1B; ost1-3: 0.10 6 0.01 min À1 vs. WT: 0.28 6 0.06 min À1 , P < 0.03), and the average half-response times were significantly larger in ost1-3 mutant leaves than WT (SI Appendix, Fig. S1C; ost1-3: 8.9 6 0.6 min vs. WT: 5.1 6 0.4 min, P < 0.002). These results suggest that high VPD-induced stomatal closure in ost1-3 was slower than in WT, consistent with previous research (5,36,40). Similar results were also found in an independent set of experiments (SI Appendix, Fig. S1 D-F). These data show that altered stomatal response kinetics to VPD can be analyzed with our system by direct comparisons to WT plants grown in parallel and analyzed within the same datasets (Materials and Methods).
Activation of anion channels in guard cells plays an important role in stomatal closure. It has been shown that the slowtype anion channel SLAC1 is required for the high-VPD response (18). As a further test of our system, we analyzed the high-VPD response in slac1-3 mutant leaves. Similar to the findings of Vahisalu et al. (18), slac1-3 mutant leaves showed a slower response to high VPD ( Fig. 1A and SI Appendix, Figs. S2 and 3). However, unlike ost1-3 mutant leaves, the steadystate stomatal conductances at low-VPD and high-VPD conditions in slac1-3 were comparable with WT under the imposed conditions (Fig. 1A). In another independent set of experiments, slac1-3 mutant leaves exhibited slightly increased steady-state stomatal conductances compared to parallel-grown WT controls (SI Appendix, Fig. S3A).
It has been shown that the rapid-activating R-type anion channel, ALMT12/QUAC1, mediates anion efflux in guard cells during ABA-triggered, darkness-triggered, and high-CO 2 -triggered stomatal closure (15,16,55). To investigate whether ALMT12/QUAC1 is required for the high-VPD response, altm12-1 mutant leaves were analyzed. Unexpectedly, no obvious difference in the high-VPD response was found between altm12-1 mutant and WT leaves (Fig. 1B). We further analyzed the high-VPD response in slac1-3/almt12-1 double mutants. Interestingly, slac1-3/almt12-1 double-mutant leaves showed a higher steady-state stomatal conductance before and after high-VPD shifts compared to WT leaves and exhibited a slower stomatal closure in response to high-VPD stimulation, similar to slac1-3 mutant leaves ( Fig. 1C and SI Appendix, Fig. S3). Average rate constants in slac1-3 and slac1-3/almt12-1 mutant leaves were significantly smaller than in WT controls ( Fig. 1D; slac1-3: 0.13 6 0.04 min À1 vs. WT: 0.31 6 0.04 min À1 , P < 0.05 and slac1-3/almt12-1: 0.10 6 0.00 min À1 vs. WT: 0.31 6 0.04 min À1 , P < 0.009 and SI Appendix, Fig. S3 C and D; slac1-3: 0.11 6 0.01 min À1 vs. WT: 0.27 6 0.03 min À1 , P < 0.001 and slac1-3/ almt12-1: 0.10 6 0.00 min À1 vs. WT: 0.33 6 0.03 min À1 , P < 0.001), and average half-response times were larger in slac1-3 and slac1-3/almt12-1 mutant leaves than in WT controls ( These experiments provide further evidence that SLAC1 and AMLT1/QUAC1 synergistically function in controlling the steady-state stomatal conductance. Furthermore, the finding that stomatal closing occurs in slac1-3/almt12-1 mutant leaves suggests that additional anion-channel genes could function in the VPD response. Investigation of Potential Roles of Calcium-Linked Signaling Mechanisms in High-VPD-Triggered Stomatal Closure. It remains unclear how guard cells perceive and respond to increases in VPD to trigger stomatal closure. High VPD initially increases transpiration rates and leads to a rapid reduction of apoplastic water content in leaves (12). The alteration of water content and/or the initial transient increase in stomatal conductance, termed the wrong-way VPD response (56, 57) (e.g., Fig. 1), may trigger Fig. 1. slac1-3/almt12-1 double-mutant leaves have a higher stomatal conductance, whereas a delayed stomatal vapor pressure deficit response is associated with defective SLAC1 in slac1-3 single and slac1-3/almt12-1 double mutants. (A-C) Time-resolved stomatal conductance (G s ) in response to increases in VPD in Col-0 (WT), slac1-3, almt12-1, and slac1-3/ almt12-1 mutant leaves under high light (450 μmolÁm À2 Ás À1 red light combined with 50 μmolÁm À2 Ás À1 blue light). VPDs are shown on top of the time axis in each panel. WT controls shown in A-C are the same because mutants and WT were grown in parallel and analyzed within the same experimental sets. (D and E) Average rate constants and average half-response times in response to VPD elevations in A-C and other figures are calculated and presented as described in Materials and Methods. Data represent mean 6 SEM. Different letters represent significant differences (P < 0.05) analyzed by one-way ANOVA followed by Holm-Sidak multiple comparisons (Results). WT: n = 4 leaves; slac1-3: n = 3 leaves; and almt12-1 and slac1-3/almt12-1: n = 5 leaves. Each analyzed leaf response was from a separate plant for all figures. osmotic/mechanical-induced cytosolic calcium concentration changes in guard cells to trigger stomatal closure. We tested the effect of 500 mM sorbitol (as osmoticum) on the cytosolic Ca 2+ concentration on preopened stomata in leaf epidermal strips of 4-to 5-week-old Arabidopsis plants expressing the Cameleon YC3.60 sensor placed under the control of the pGC1 promoter (58). Both the addition of sorbitol (hyperosmotic stress) and removal of sorbitol (hypo-osmotic stress) caused clear F€ orster resonance energy transfer (FRET) of the Cameleon sensor, with a decrease of cyan fluorescent protein (CFP) fluorescence and an increase in cpVenus fluorescence ( Fig. 2A). When representing the data as normalized ratios (DR/R 0 ), which depict cytosolic Ca 2+ concentration changes, rapid and transient cytosolic [Ca 2+ ] increases were observed after both the hyperosmotic and hypo-osmotic stimuli (Fig. 2B). Of note, the removal of sorbitol induced a higher maximum Ca 2+ transient compared to the one observed after its addition (Fig. 2C).
To further investigate the potential roles of calciumpermeable channels in high-VPD-triggered stomatal closure, we focused on putative calcium channels with higher expression levels in guard-cell protoplasts, including glutamate-like receptors genes GLR3.2 and GLR3.3 and cyclic nucleotide-gated channel genes CNGC5, CNGC6, and CNGC20 based on transcriptome analyses (SI Appendix, Fig. S4 B-D). Stomatal conductance responses to an increase in VPD were analyzed in glr3.2/3.3, cngc5-1/6-1, and cngc20 T-DNA insertion mutants and a deletion mutant of two homologous tandem repeat genes CNGC19 and CNGC20, cngc19/20crispr (Fig. 2 F-I and SI Appendix, Fig. S5 C-F). Notably, no obvious differences in the VPD responses were observed compared to parallel-grown WT leaves in all of these mutants.
High VPD and High CO 2 Use Different Pathways to Trigger Stomatal Closure. Defects in several ABA biosynthesis and ABA signaling mutants in Arabidopsis increase steady-state stomatal conductances and thus shift VPD responses to higher stomatal conductance values (5, 37). Recent time-resolved stomatal conductance analyses have led to the model that the rapid stomatal conductance response to VPD elevation occurs in parallel to the ABA signaling pathway (5,34). A signaling pathway parallel to ABA signaling is the CO 2 signaling pathway (46). MAP kinases MPK4 and MPK12 are essential for high-CO 2 -induced stomatal closure, but not for ABA-induced stomatal closure (75,76). To investigate whether VPD and CO 2 share a common signaling pathway upstream of anion-channel activation, the stomatal response to VPD elevation was analyzed in two independent mpk12/4GC double-mutant plant lines, in which guard cell MPK4 was silenced by an artificial micro-RNA driven by the MPK12 promoter in the mpk12-4 mutant background (75,77). These mpk12/4GC double-mutant alleles are the strongest presently known specific high-CO 2 -response impaired mutants in Arabidopsis (75). These analyses showed that the mpk12/4GC double-mutant leaves exhibited a much higher steadystate stomatal conductance than WT leaves (Fig. 3 A and B), consistent with previous findings (75). Previous research has shown a slightly lower stomatal index and no clear effect on stomatal density in mpk12/4GC double mutants (75). Leaves of mpk12/4GC double mutants responded rapidly to VPD elevations with response times similar to WT (Fig. 3 A and B). Relative changes in stomatal conductance at early time points in response to VPD increase (10 min) were even larger in these mutant leaves than in WT controls ( Fig. 3 C and D). No difference in average rate constants was found between mpk12/4GC double mutants and WT (   (75). These findings that absolute high VPD-induced stomatal closing responses were stronger in mpk12/ 4GC-2 double mutants than in WT controls provide evidence that high VPD and high CO 2 trigger stomatal closure by different early signaling pathways.
As the slowing of the VPD response in these quadruplemutant alleles was significantly, but only slightly, slower, we investigated these mutants in the independent laboratory of E.M. in genotype-blinded experiments using a distinct wholeplant gas-exchange system. Higher stomatal conductances and in part increased half-response times in m3kδ1/δ5/δ6/δ7 quadruple mutants were independently found in intact wholeplant gas-exchange analyses (SI Appendix, Fig. S8). These results suggest that B3-family Raf-like M3Ks play a role in stomatal closure induced by VPD elevation. GHR1 Is Required for High VPD-Induced Stomatal Closure. GHR1 encodes a receptor-like (pseudo)kinase, which is involved in activating the S-type anion channel SLAC1 to trigger stomatal closure (77,79,80). The involvement of GHR1 in the stomatal high-VPD response was analyzed in ghr1-3 and ghr1-6 exon-insertion mutant leaves. Leaves of the strong ghr1-3 allele and the weak ghr1-6 allele mutants (80) showed higher steady-state stomatal and Raf-like m3kδ5 mutant under high light (450 μmolÁm À2 Ás À1 red light combined with 50 μmolÁm À2 Ás À1 blue light). WT controls shown in B and C are the same because mutants and WT were grown in parallel and analyzed within the same experimental sets. (D and E) Rate constants and half-response times in response to VPD elevations were calculated and presented. Data represent mean 6 SEM. n = 5 for each genotype. Different letters represent significant differences (P < 0.05) analyzed by one-way ANOVA followed by Holm-Sidak multiple comparisons. conductances in both high-VPD and low-VPD conditions compared to WT controls (Fig. 7 A and B and SI Appendix, Fig. S9 A  and B). Although high-VPD treatment triggered a degree of stomatal closure in ghr1-3 and ghr1-6 mutant leaves, ghr1 mutant leaves exhibited a clear difference in response kinetics compared to WT leaves under both high-light (500 μmolÁm À2 Ás À1 ; Fig. 7 A and B and SI Appendix, Fig. S9 A and B) and low-light (150 μmolÁm À2 Ás À1 ; SI Appendix, Fig. S9C) experimental conditions. Rate constants in the strong ghr1-3 mutant leaves were significantly smaller than in WT leaves (Fig. 7 C and D; ghr1- . Taken together, our results show that GHR1 is required for the rapid stomatal closure induced in response to VPD increases. Furthermore, in contrast to all other mutants investigated in the present study, ghr1 mutants showed no initial transient stomatal conductance increase (wrong-way response) as found under high-light conditions (Fig. 7 A and B).

Discussion
The continuing elevation in global temperature is suggested to lead to a trend of rising vapor pressure deficit, though VPD is variable depending on conditions (81). Increases in VPD affect transpiration and regulate stomatal conductance of terrestrial plants (82)(83)(84)(85). The molecular mechanisms of VPD sensing and signal transduction remain to a large degree unknown. Previous studies have focused on the roles of ABA in the VPD response (introduction). In the present study, we have investigated mutants affected in several stomatal responses. Mutants that disrupt specific known signaling mechanisms or pathways were investigated.
Analyses of the ost1-3 single mutant (SI Appendix, Fig. S1) confirmed that OST1/SnRK2.6 is a central regulator of hydroactive high VPD-induced stomatal closure (5,36,39) and provided an internal control for our analysis system. The OST1 protein kinase regulates both S-type anion channels and the R-type anion channel ALMT12/QUAC1 in guard cells or Xenopus oocytes in response to several stimuli (28,(86)(87)(88). In the present study, analyses of slac1-3 single-mutant leaves were used as an additional system control, confirming that the S-type anion channel SLAC1 contributes to high-VPD-mediated stomatal closing (Fig. 1 A, D, and E and SI Appendix, Fig. S3 A, C, and D) (18).
Notably, the single almt12-1 mutant exhibited high VPD-induced stomatal closing kinetics similar to WT plants ( Fig. 1 B, D, and E). Interestingly, however, slac1-3/almt12-1 double-mutant leaves showed a dramatic increase in the steadystate stomatal conductance and a slowing of the high-VPD response, which was kinetically similar to the slac1-3 single mutant ( Fig. 1 C, D, and E and SI Appendix, Fig. S3). These results indicate that S-type anion channels play a rate-limiting role, and R-type anion channels function together with S-type anion channels in setting the steady-state stomatal conductance Fig. 7. Stomatal VPD responses, including the initial transient increase in stomatal conductance, are impaired in ghr1 mutants. (A and B) Time-resolved stomatal conductance (G s ) in response to increases in VPD in Col-0 (WT) and ghr1-3, and ghr1-6 mutant leaves under high light (450 μmolÁm À2 Ás À1 red light combined with 50 μmolÁm À2 Ás À1 blue light). VPDs are shown on top of the time axis in each panel. (C-F) Average rate constants (C and D) and average halfresponse times (E and F) in response to VPD elevations were calculated and presented as described. Data represent mean 6 SEM. **P < 0.01 (as analyzed by Student's t test compared to WT control). WT: n = 4 leaves in A and n = 5 leaves in B; ghr1-3: n = 5 leaves; and ghr1-6: n = 5 leaves. Data represent one of two independent sets of experiments. Similar results were found in an independent set of experiments in each mutant allele (SI Appendix, Fig. S8).
under the imposed conditions. Given redundancies and nonlinearities in signal-transduction networks (89,90), our data do not unequivocally exclude a contribution of R-type anion channels to the high VPD-induced stomatal closing response. The slow, but continued, stomatal closing response in slac1/almt12 double-mutant leaves suggests that additional anion channels and mechanisms function in parallel to these channels.
In this study, higher steady-state stomatal conductances in m3kδ5 single-mutant and m3kδ1/δ5/δ6/δ7 quadruple-mutant alleles were not attributable to increases in stomatal index (SI Appendix, Fig. S10). Furthermore, the slower response to an increase in VPD in m3kδ1/δ5/δ6/δ7 quadruple mutants was not simply caused by the increased initial steady-state stomatal conductance in these quadruple mutants for several reasons. First, we found that the larger half-response time in m3kδ1/δ5/δ6/δ7 quadruple mutants was consistently observed, even at low light, with reduced initial steady-state stomatal conductance levels (SI Appendix, Fig. S7). Second, unlike m3kδ1/δ5/δ6/δ7 quadruple-mutant alleles, no significant differences in halfresponse times were found in additional controls in WT plants exposed to high-light and low-light conditions under 400 parts per million (ppm) CO 2 or 150 ppm CO 2 (SI Appendix, Fig.  S11). Third, a rapid stomatal response to increases in VPD was observed in the early CO 2 signaling mpk12/4GC mutant lines that also exhibited high initial steady-state stomatal conductance (Fig. 3). Therefore, our findings suggest that B3 M3Ks are involved in accelerating the hydroactive signaling response to an increase in VPD.
Unlike the m3kδ5 single-mutant and m3kδ1/δ5/δ6/δ7 quadruple-mutant alleles, no obvious high-VPD-triggered stomatal response phenotype was observed in m3kδ1/δ7 double, m3kδ1/ δ6/δ7 triple, and in an OK-quatdec quattuordecuple mutant (Fig.  4D). There are 80 M3K genes encoded in the Arabidopsis genome (97)(98)(99). The OK-quatdec line includes mutations in 14 Raf-like kinase genes, including several osmotic stress-activated B4 family members (53). The present findings indicate a less important role of these B4-family Raf-kinase members in the VPD response, compared to the osmotic stress response in whole seedlings. The present results further suggest that M3Kδ5/RAF6 plays an important role in VPD-induced stomatal regulation.
A previous study suggested that the receptor-like pseudokinase GHR1 mediates stomatal closure in response to ABA, CO 2 , ozone, and darkness, but may not be required for high-VPD-triggered stomatal closure (80). In the present study, the severe ghr1-3 mutant allele showed a clear impairment in high VPD-induced stomatal closure ( Fig. 7 and SI Appendix,  Fig. S9). GHR1 can directly activate the SLAC1 anion channel in Xenopus oocytes when these two proteins are artificially linked via BiFC constructs (79,80). GHR1 is required for stomatal closure in response to ABA, hydrogen peroxide, ozone, high CO 2 , and low light intensity (79,80). GHR1-mediated activation of the S-type anion channel SLAC1 is negatively regulated by the protein phosphatase 2C, ABI2, in Xenopus oocytes (77,79). As we show that early CO 2 signaling is not required for the stomatal high-VPD response (Fig. 3), the defect in the stomatal high-VPD response in ghr1 mutant leaves is likely caused by a GHR1 function downstream of early CO 2 signaling or a different signaling pathway that makes use of GHR1. It has been suggested that the receptor-like pseudokinase GHR1 serves as a scaffolding protein to recruit the calcium-dependent protein kinase CPK3 for regulating SLAC1 anion-channel activation (80). Higher-order cpk mutants analyzed here included cpk3 mutation, but did not affect the high-VPD response (SI Appendix, Fig. S6 A and B). Therefore, GHR1 may recruit additional, yet to be discovered signaling components to regulate the stomatal high-VPD response.
Classical research showed that high VPD/low humidity triggers a transient stomatal opening before stomatal closing (56,(100)(101)(102). Notably, this transient increase in stomatal conductance in response to an increase in VPD (wrong-way response) may in part result from a simplified assumption of saturated vapor pressure inside leaves during stomatal conductance recordings at high stomatal conductances (103). The wrongway response was found in WT and most of the studied mutants, including ost1-3, slac1-3, slac1-3/almt12-1, and mpk12/ 4GC mutants, with high steady-state stomatal conductances. These data indicate that the larger stomatal conductance of ghr1 leaves is not solely due to the above-simplified assumption. Interestingly, the lack of a clear wrong-way response observed in both ghr1 mutant alleles points to an additional role of GHR1. It has been proposed that the initial hydropassive loss of turgor in epidermal cells adjacent to guard cells during high-VPD/low-humidity stimulation drives the transient stomatal opening (2,104). In a previous study, a transient stomatal opening could only be observed in stomata with living subsidiary cells in Tradescantia pallida or when guard cells were surrounded by living epidermal cells in Vicia faba (3). Therefore, the lack of high-VPD-triggered transient stomatal opening in ghr1 mutants suggests a function of GHR1 in the mechanical interaction between guard cells and neighboring pavement cells.
It remains to be elucidated whether mechanical and hyperosmolarity-sensing mechanisms are involved in high VPD-induced stomatal closure. OSCA channels have been reported to be mechanical-and osmotic-sensing calcium channels in higher plants, with structural analyses resolving possible mechanisms (59)(60)(61)(62)(63)(64). In the present study, we generated higher-order osca mutants, based on guard-cell transcriptome analyses from independent laboratories (SI Appendix, Fig. S4) (58,65). However, the stomatal responses to high VPD were not altered in osca1-2/1.3/2.2/2.3/3.1 quintuple-mutant leaves (Fig. 2D), and hyperosmotic stress was insufficient to activate the rice ortholog OsOSCA2;1 when expressed alone in Xenopus oocytes (64). It has been reported that the activation of Raflike kinases or SnRK2s by hyperosmolarity was not altered in an osca septuple mutant (53). These results suggest that high VPD-induced stomatal closure may be independent of OSCAs and, in particular, the five more guard-cell-expressed OSCAs analyzed here. A recent study has identified an additional OSCA that is not highly expressed in guard cells, OSCA1.7, but is together with OSCA1.3 required for PAMP-induced stomatal closing during pathogen signaling (66). However, we found WT-like high-VPD responses in osca1.3/1.7 double-mutant leaves. We show that guard-cell cytosolic Ca 2+ concentrations transiently increase in response to both hyperosmolality and hypo-osmolality (Fig. 2 A-C), but whether Ca 2+ signaling is required for the stomatal high-VPD response remains to be investigated. Although no obvious stomatal closure phenotype in response to VPD elevation was found in the investigated calcium-permeable channel mutants, cpk quintuple, and higher-order cbl mutants (Fig. 2 D-G and SI Appendix, Fig.  S6), we do not exclude a possible role of calcium in stomatal high-VPD signal transduction. Additional higher-order mutants in Ca 2+ signaling components may be needed to impair the stomatal high-VPD response.
In summary, the identification of B3-family Raf-like kinases and the transmembrane receptor-like GHR1 in the stomatal VPD response advances the understanding of molecular and genetic mechanisms that mediate the hydroactive stomatal VPD response. Interestingly, the initial transient stomatal opening wrong-way response to a VPD increase that precedes stomatal closing (56,100) was disrupted in ghr1 mutant alleles, indicating a function of GHR1 in guard-cell-pavement-cell communication during VPD signaling. Furthermore, our data exclude the involvement of the early CO 2 signal-transduction pathway in rapid high VPD-induced stomatal closure. Further research will be required to elucidate the upstream high VPD-sensing mechanisms, which are shown here to require B3-family Raf-like kinases and the receptor-like (pseudo) protein kinase GHR1. Understanding the detailed molecular and genetic mechanisms underlying the VPD response is relevant for models that predict plant transpiration in response to climate change and could become relevant for the future development of plants with enhanced resilience to climate change.
Time-Resolved Intact VPD-Dependent Leaf Stomatal Conductance Experiments. Stomatal conductance recordings from intact leaf 8 to leaf 12 of 5.5-to 8-wk-old plants in different batches of experiments were conducted starting 1 to 2 h after growth-chamber light onset. An LI-6800 Portable Photosynthesis System with an integrated Multiphase Flash Fluorometer (6800-01A; Li-Cor Inc.) was used to measure gas exchange in Arabidopsis leaves. Leaves were clamped and kept at 400 ppm ambient CO 2 ; 21°C heat-exchanger temperature; 0.8 kPa VPD; 450 μmolÁm À2 Ás À1 red light combined with 50 μmolÁm À2 Ás À1 blue light in most of the experiments, as done in other studies analyzing the VPD response (5); and 500 μmolÁs À1 incoming airflow rate for 1.5 to 2.5 h until stomatal conductance stabilized. Some experiments presented in SI Appendix were performed under lower light intensity at 135 μmolÁm À2 Ás À1 red light combined with 15 μmolÁm À2 Ás À1 blue light, as described in the figure legends. For stomatal responses to VPD shifts, stomatal conductances were recorded every 30 s at ∼0.8 kPa VPD (65 to 72% relative air humidity), followed by ∼2.0 kPa VPD (20 to 32% relative air humidity) in most experiments, as indicated in the figures. VPD values were computed internally by the gas-exchange system. Stomatal conductance measurements in response to VPD shifts are affected by multiple parameters and therefore should be viewed as semiquantitative relative responses. Therefore, parallel-grown WT control plants were investigated in all experiments, which enabled direct resolution of non-WT-like impaired responses in mutants. Furthermore, the following measures were taken in these experiments. Matching between sample and reference infrared gas analyzers was executed before each experiment and 50 min after shifting to high VPD to adjust the errors due to the changes in humidity. The stomatal conductances at high VPD in each experiment were corrected by the differences in average stomatal conductances between 5 min before matching and 3 to 8 min after matching. The values of stomatal conductances were plotted from 2 min to 40 min or 50 min after VPD increase as in the figures. The first 2min data points were omitted because of the inaccurate read and calculation of stomatal conductances before VPD reached more reliable values. High-VPD response rate constants were approximated by exponential one-phase decays starting from 2 min after shifting to high VPD to time points that that showed the minimum stomatal conductances in the investigated genotypes and in independent experiments using GraphPad Prism 7.0e software. Stomatal conductances at the following time points after high-VPD shifts were analyzed during curve fitting as follows: 2 to ∼10 min: WT in Fig. 5D; 2 to ∼12 min: WT in Figs. 1D, 3E, 6C, and 7 C and D and SI Appendix, Figs. S3C and S9 D and E and almt12-1 in Fig. 1D; 2 to ∼14 min: WT in Fig. 6D and SI Appendix, Figs. S1B, S3D, S7C, and S9F and m3kδ5 mutant alleles in Fig. 5D; 2 to ∼18 min: mpk12/4GC mutant alleles in Fig. 3E; 2 to ∼20 min: m3kδ1/δ5/δ6/δ7 mutant alleles in Fig. 6 C and D and SI Appendix, Fig. S7C; 2 min to the end of the experiment: WT in SI Appendix, Fig. S1E, ost1-3 in SI Appendix, Fig. S1 B and E, slac1-3 in Fig. 1D and SI Appendix, Fig. S3C;,slac1-3/almt12-1 in Fig. 1D and SI Appendix, Fig. S3D, and ghr1 mutant alleles in Fig. 7 C and D and SI Appendix, Fig. S9 D-F. The times needed to reach 50% of the recorded maximum stomatal conductance reductions starting from 2 min after shifting to high VPD were calculated as half-response times.
In-Gel Kinase Assays. Transient Raf-M3Ks expression in mesophyll cell protoplasts and the following in-gel kinase assays were performed as described (54). In brief, 30 μg of mesophyll cell protoplasts isolated from 3-to 4-wk-old Arabidopsis leaves were transfected with 20 μg of pUC18 plasmids carrying 35S:Raf-M3K-GFP:nosT using the polyethylene glycol-mediated method (107). After overnight incubation in incubation buffer ( Analyses of Guard-Cell Cytosolic Calcium Dynamics. For analyses of cytosolic Ca 2+ dynamics in guard cells, we used A. thaliana Col-0 plants expressing Yellow Cameleon YC3.60 under control of the pGC1 promoter (58). Plants were grown on soil under long-day conditions (16/8 h at a photon fluence rate of 75 μmolÁm À2 Ás À1 and a temperature of 20°C). Small pieces of mature leaves of 4-wk-old plants were glued to the cover slide by using a medical adhesive (Hollister Inc.); upper cell layers were gently removed with a razor blade, thus isolating epidermal strips. The leaf epidermal strips were preincubated in a buffer solution (5 mM KCl, 10 mM MES, and 50 μM Ca 2+ , pH 6.15, adjusted with Tris base) for 2 to 3 h in white light (photon fluence rate of 125 μmolÁm À2 Ás À1 ) to induce the stomatal opening to get turgid guard cells. The epidermal strips were then mounted on an open-top imaging chamber and perfused with the same solution (5 mM KCl, 10 mM MES, and 50 μM Ca 2+ , pH 6.15) as described (58,108). To induce the hyperosmotic stress, after 5 min, the strips were perfused with the standard solution supplemented with 500 mM sorbitol. Epidermal strips were kept in this solution for 10 min, and then the hypertonic solution was replaced again with the standard buffer solution, without sorbitol, to induce the hypo-osmotic stress.
To perform Ca 2+ imaging experiments, images of epidermises with stomatal guard cells were acquired with a Nikon TE300 inverted fluorescence microscope, with a FLUORESCENCE Module ILLUMINATOR (TE-FM) epi-fluorescence attachment (Nikon Inc.) using a Nikon 60× Plan Apo oil objective with a numerical aperture (n.a.) of 1.4. Excitation light was produced by a 75-W xenon fluorescent lamp (Osram) that was attenuated by 97% (3% light transmission) by using both 4× and 8× neutral density filters to reduce exposure of the fluorescent reporters and cells to epifluorescence excitation. Cameleon YC3.60 specific excitation and detection was achieved by using a yellow fluorescent protein/CFP filter set: 440/20 nm excitation, 485/40 nm emission for CFP, and 535/30 nm emission for cpVenus with a 455 nm dichroic long-pass filter (Chroma Technology). The filter wheel, shutter, and CoolSNAP chargecoupled device camera from Photometrics (Roper Scientific) were controlled by using Metafluor System software (MDS Inc.). Images were acquired every 3 s, on a time range of >30 min.
Fluorescence intensity was determined over regions of interest (ROIs) that corresponded to single guard cells. cpVenus and CFP emissions of the analyzed Raf-like kinases and receptor-like (pseudo)kinase GHR1 are required for stomatal vapor pressure difference response