
==== Front
Plant Physiol
Plant Physiol
plphys
Plant Physiology
0032-0889
1532-2548
Oxford University Press US

38833587
10.1093/plphys/kiae320
kiae320
Research Article
Signaling and Response
AcademicSubjects/SCI01270
AcademicSubjects/SCI01280
AcademicSubjects/SCI02286
AcademicSubjects/SCI02287
AcademicSubjects/SCI02288
Stomatal CO2 responses at sub- and above-ambient CO2 levels employ different pathways in Arabidopsis
https://orcid.org/0000-0001-9521-9045
Koolmeister Kaspar Institute of Technology, University of Tartu, Nooruse 1, Tartu 50411, Estonia
Institute of Bioengineering, University of Tartu, Nooruse 1, Tartu 50411, Estonia

https://orcid.org/0000-0003-1767-6413
Merilo Ebe Institute of Technology, University of Tartu, Nooruse 1, Tartu 50411, Estonia

https://orcid.org/0000-0002-6392-859X
Hõrak Hanna Institute of Technology, University of Tartu, Nooruse 1, Tartu 50411, Estonia

https://orcid.org/0000-0002-6895-3583
Kollist Hannes Institute of Bioengineering, University of Tartu, Nooruse 1, Tartu 50411, Estonia
Université Paris-Saclay, CNRS, INRAE, Université Evry, Université Paris Cité, Institute of Plant Sciences Paris-Saclay (IPS2), 91190 Gif sur Yvette, France

Author for correspondence: hanna.horak@ut.ee (H.H.)
Kaspar Koolmeister and Hanna Hõrak contributed equally.

H.H. designed the study, H.H., K.K., and E.M. performed experiments, H.H., K.K., E.M., and H.K. analyzed data, H.H., K.K., and H.K. wrote the manuscript, all authors commented, edited and approved the final manuscript.

The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (https://academic.oup.com/plphys/pages/General-Instructions) is Hanna Hõrak.

Conflict of interest statement. None declared.

9 2024
04 6 2024
04 6 2024
196 1 608620
29 8 2023
12 5 2024
21 6 2024
© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Stomatal pores that control plant CO2 uptake and water loss affect global carbon and water cycles. In the era of increasing atmospheric CO2 levels and vapor pressure deficit (VPD), it is essential to understand how these stimuli affect stomatal behavior. Whether stomatal responses to sub-ambient and above-ambient CO2 levels are governed by the same regulators and depend on VPD remains unknown. We studied stomatal conductance responses in Arabidopsis (Arabidopsis thaliana) stomatal signaling mutants under conditions where CO2 levels were either increased from sub-ambient to ambient (400 ppm) or from ambient to above-ambient levels under normal or elevated VPD. We found that guard cell signaling components involved in CO2-induced stomatal closure have different roles in the sub-ambient and above-ambient CO2 levels. The CO2-specific regulators prominently affected sub-ambient CO2 responses, whereas the lack of guard cell slow-type anion channel SLOW ANION CHANNEL-ASSOCIATED 1 (SLAC1) more strongly affected the speed of above-ambient CO2-induced stomatal closure. Elevated VPD caused lower stomatal conductance in all studied genotypes and CO2 transitions, as well as faster CO2-responsiveness in some studied genotypes and CO2 transitions. Our results highlight the importance of experimental setups in interpreting stomatal CO2-responsiveness, as stomatal movements under different CO2 concentration ranges are controlled by distinct mechanisms. Elevated CO2 and VPD responses may also interact. Hence, multi-factor treatments are needed to understand how plants integrate different environmental signals and translate them into stomatal responses.

Arabidopsis stomatal closure in response to CO2 concentration elevation relies on different molecular mechanisms at sub- and above-ambient CO2 levels.

Estonian Research Council 10.13039/501100002301 PSG404 Institute of Technology 10.13039/100008961 PRG719 PRG1620 PRG433 PRG2004 Estonian Research Infrastructure Roadmap project TAIM European Regional Development Fund 10.13039/501100008530 Center of Excellence in Molecular Cell Engineering
==== Body
pmcIntroduction

Atmospheric CO2 concentration has nearly doubled within the past 150 years. As a result, the global average temperature has increased by 1.5 °C and relative air humidity decreased across many vegetated areas (Yuan et al. 2019). Increased air temperature and decreased air humidity lead to rising values of vapor pressure deficit (VPD), the difference between actual and saturated air vapor pressures, and this increases transpiration from plants. In order to survive such conditions, plants need to adjust their water management by regulating stomatal conductance. Each stoma is composed of two guard cells and an opening formed between them; the size of the opening is controlled by increasing or decreasing guard cell turgor pressure. Guard cells adjust their turgor pressure in response to various abiotic stimuli to balance water loss and CO2 uptake for photosynthesis. Understanding how stomata respond to changing CO2 concentrations and increasing VPD is needed for breeding climate-ready crops.

Plant stomata close in response to elevated CO2 concentration and open in response to decreased CO2 concentration. Ultimately, elevated CO2 activates the S-type anion channel SLOW ANION CHANNEL-ASSOCIATED 1 (SLAC1), which causes rapid stomatal closure (Negi et al. 2008; Vahisalu et al. 2008). CO2 can enter guard cells through the PLASMA MEMBRANE INTRINSIC PROTEIN 2 (PIP2) plasma membrane channel or by transmembrane diffusion (Katsuhara and Hanba 2008; Wang et al. 2016). Carbonic anhydrases BETA CARBONIC ANHYDRASE 1 (CA1) and BETA CARBONIC ANHYDRASE 4 (CA4) accelerate the conversion of intracellular CO2 to bicarbonate (HCO3−), which can act as a second messenger (Hu et al. 2010). In guard cells, CO2/HCO3− promotes interaction between the protein kinases MITOGEN-ACTIVATED PROTEIN KINASE 4 (MPK4)/MITOGEN-ACTIVATED PROTEIN KINASE 12 (MPK12), and the Raf-like kinase HIGH LEAF TEMPERATURE 1 (HT1), leading to HT1 inhibition, which is an essential step in the regulation of stomatal responses to CO2 (Hashimoto et al. 2006; Hashimoto-Sugimoto et al. 2016; Hõrak et al. 2016; Jakobson et al. 2016; Takahashi et al. 2022; Yeh et al. 2023). HT1 phosphorylates the CONVERGENCE OF BLUE LIGHT AND CO2 (CBC1)/CONVERGENCE OF BLUE LIGHT AND CO2 2 (CBC2) Raf-like kinases that function downstream of HT1 and these Raf kinases can inhibit the S-type anion channel activation via a currently unknown mechanism (Hõrak et al. 2016; Hiyama et al. 2017; Hayashi et al. 2020). Stomata in HT1-deficient plants do not respond to CO2 concentration changes while carbonic anhydrase, MPK12 and SLAC1 mutants exhibit impaired stomatal CO2 responses (Hashimoto et al. 2006; Vahisalu et al. 2008; Hu et al. 2010; Hashimoto-Sugimoto et al. 2016; Hõrak et al. 2016; Jakobson et al. 2016).

Elevated VPD increases transpiration and reduces epidermal turgor that due to mechanical interactions between guard and epidermal cells in angiosperms leads to faster light-induced stomatal opening (Mott et al. 1999; Pichaco et al. 2024). To prevent wilting, stomata close under elevated VPD. Abscisic acid (ABA) is a drought-induced plant stress hormone and an important stomatal regulator (Cutler et al. 2010). VPD has a direct effect on ABA concentration as increased ABA levels in angiosperms were observed 20 min after increasing VPD from 0.7 to 1.5 kPa (McAdam and Brodribb 2015), promoting the conclusion that ABA may be involved in elevated VPD-induced stomatal closure in angiosperms (McAdam et al. 2016). The protein kinase OPEN STOMATA 1 (OST1) and the leucine-rich receptor-like pseudokinase GUARD CELL HYDROGEN PEROXIDE-RESISTANT 1 (GHR1) are activated in the presence of ABA and trigger anion efflux through the major guard cell slow-type anion channel SLAC1 (Brandt et al. 2012; Hua et al. 2012; Sierla et al. 2018). OST1 and GHR1 are both involved in elevated VPD- and CO2-induced stomatal closure response (Xue et al. 2011; Merilo et al. 2018; Sierla et al. 2018; Hsu et al. 2021; Jalakas et al. 2021b).

Understanding CO2-induced plant stomatal closure responses is essential for future plant breeding. Due to changing climate conditions, it is also important to understand if and how stomatal CO2 regulation is affected by elevated VPD levels. Previous work in grasses suggests that elevated VPD levels reduce both stomatal conductance and stomatal sensitivity to CO2 concentration changes (Morison and Gifford 1983) but the interactions of CO2 and humidity responses in dicots remain poorly understood.

Under light, CO2 concentration inside the leaf is usually below ambient levels and this causes stomatal opening. Stomata close when CO2 concentration inside the leaf increases to ambient levels, and an additional rise in CO2 concentration to above-ambient levels causes further stomatal closure (Brodribb et al. 2009; Hõrak et al. 2017). Thus, stomatal closure response exists within sub-ambient as well as in above-ambient CO2 levels; however, it is not clear whether these responses are controlled by the same regulators. In some studies, CO2-induced stomatal closure is defined as the response to an increase in CO2 concentration from ambient to above-ambient levels (Franks and Britton-Harper 2016; Hõrak et al. 2017), while some studies use a [CO2] change from sub-ambient to above-ambient levels (Azoulay-Shemer et al. 2015; Chater et al. 2015). Data from previous studies comparing CO2 responses in ferns and angiosperms suggest that stomatal responses to CO2 are different, when changing CO2 levels in the sub-ambient or above-ambient ranges (Brodribb et al. 2009; Hõrak et al. 2017). To address the underlying mechanisms of CO2-induced stomatal closure at different CO2 transitions under normal and elevated VPD conditions, we studied plants deficient either in guard cell anion channel SLAC1 and its activation (slac1-3, ghr1-3, and ost1-3) or in the CO2-specific stomatal signaling branch regulated by MPK12 and HT1 kinases (mpk12-4, ht1-2, and ht1-8D) and carbonic anhydrases CA1 and CA4 (ca1ca4). Our results show that different stomatal regulators have a different degree of importance in CO2-induced stomatal closure in sub-ambient and above-ambient CO2 levels and are also differently affected by VPD.

Results

Stomatal closure kinetics are different between sub-ambient to ambient and ambient to above-ambient [CO2] transitions

We analyzed stomatal responses to CO2 in the sub-ambient and above-ambient concentration ranges in the model plant Arabidopsis (Arabidopsis thaliana) to clarify whether these responses are controlled by the same or by different regulators. Four different CO2 transition sequences were used (Fig. 1, A, C, E, and G); in two setups, we applied high VPD (2.3 kPa) as an additional factor before and throughout CO2 treatments (Fig. 1, A and E). Experiments were started with ambient (400 ppm) CO2 concentration and each consecutive CO2 treatment lasted for 2 h. This approach allowed us to investigate plant stomatal response to different CO2 transitions—stomatal opening in response to CO2 transition from 400 to 100 ppm (from here on referred to as 400–100; Fig. 1, A, C, and G), from 800 to 400 ppm (800–400; Fig. 1C), or from 800 to 100 ppm (800–100; Fig. 1E), and stomatal closure in response to CO2 transition from 100 to 400 ppm (100–400; Fig. 1A), from 400 to 800 ppm (400–800; Fig. 1, A, C, and E), or from 100 to 800 ppm (100–800; Fig. 1G). We analyzed the amplitude and speed of stomatal responses, whereas the speed was defined here as stomatal 75% response time (see Methods and Supplementary Fig. S1).

Figure 1. Kinetics of CO2-induced stomatal responses in wild-type Arabidopsis with regular and elevated vapor pressure deficit (VPD). Col-0 wild-type Arabidopsis stomatal response to sequential changes in CO2 concentration under regular (A, C, E, G) and elevated (A, E) VPD conditions. Mean stomatal conductance ± standard error of the mean (SEM) is shown. B, D, F, H) Boxplot of 75% response time (min) of stomatal response to CO2 concentration changes. Boxes represent 25% to 75% quartiles and median as the horizontal lines, whiskers indicate the smallest and largest values, and points show individual plant values. B, F) Statistically significantly different groups are indicated with different letters (Two-way ANOVA with Tukey post hoc test, P < 0.05). D, H) Statistically significantly different groups are indicated with different letters (One-way ANOVA with Tukey post hoc test, P < 0.05). Sample size was 5 in (A, B, C, D, E, F) and 14 in (G, H). Start of the first treatment was between 11:30 and 12:30.

Stomatal closure responses had clearly different 75% response times in Col-0 wild-type plants at different CO2 transitions: the 400–800 stomatal closure was faster than 100–400 closure under normal VPD (Fig. 1, A and B). The rapid 400–800 response was consistent throughout all experimental setups (Fig. 1, A to F). The observed different kinetics of these CO2-induced stomatal closure responses suggest that they could be regulated by different components. Stomatal opening responses to sub-ambient CO2 levels, 400–100 (Fig. 1, A to D, G, and H) and 800–100 (Fig. 1, E and F) had relatively slow 75% response times, whereas the 800–400 response, opening from above-ambient to ambient CO2 levels, had lower response range and faster 75% response time (Fig. 1, C and D). This suggests that stomatal opening to sub-ambient CO2 concentrations is a much stronger, albeit slower, response than stomatal opening during the above-ambient to ambient [CO2] change.

High VPD but not the order of CO2 transitions affects stomatal CO2-response kinetics

To study whether high VPD affects CO2 responses, we conducted the 400–100–400–800 and the 400–800–100 CO2 transitions under conditions where plants were first acclimatized to increased VPD (2.3 kPa) for ∼3 h and then subjected to changes in CO2 levels under the elevated VPD conditions. Under such conditions, both the 400–100 and 800–100 stomatal opening responses were significantly faster than under normal VPD (Fig. 1, A, B, E, and F). Stomatal closure in sub-ambient to ambient CO2 concentration range was also enhanced under high VPD conditions, resulting in a shorter stomatal 75% response time during the 100–400 ppm [CO2] transition compared with normal VPD (Fig. 1, A and B), whereas we were not able to detect an effect of high VPD on the 400–800 ppm [CO2] transition 75% response time (Fig. 1, A, B, E, and F). These results suggest that during high VPD stress, plant stomata could be primed for faster movements in the sub-ambient to ambient CO2 concentration range.

In our first experiments, the 400–100 opening stimulus was applied before the 400–800 closure. To test whether the exposure to sub-ambient CO2 levels had an effect on stomatal CO2-responsiveness, we also applied the CO2 transitions in reverse order (400–800–400–100, Fig. 1C). Stomatal 75% response times during the 400–100 transition appeared slightly shorter when it was the last transition (Fig. 1, C and D) than when it was the first (Fig. 1, A, B, G, and H), whereas the 400–800 transition response speed was unaffected by previous treatments (Fig. 1, A, B, and C to F). Thus, there were no major effects of the order of CO2 concentration transitions on stomatal CO2-response kinetics.

The CO2-specific pathway components MPK12 and carbonic anhydrases are more involved in sub-ambient to ambient than in above-ambient CO2-induced stomatal closure

To better understand the role of the CO2-signalling module comprising MPK12, HT1 and carbonic anhydrases CA1 and CA4 in stomatal CO2 responses at different CO2 levels, we analyzed the mpk12–4, ht1-2, ht1-8D, and ca1a4 in a similar experimental set-up as described for wild-type plants in Fig. 1A. All the CO2-signalling mutants showed very little stomatal closure in magnitude compared with wild-type plants, whereas 75% response time was significantly reduced only in the ca1ca4 mutant in response to the 100–400 transition (Fig. 2, A, C, and E). The 75% response times for the ht1-2 and ht1-8D mutants in this CO2 transition are not informative due to hardly any stomatal response in these mutants (Fig. 2, A, C, and E). The response of mpk12-4 to the 100–400 transition was as fast as in wild-type, but significantly smaller in magnitude (Fig. 2, A, C, and E). In response to the 400–800 transition, the HT1 mutants had a very weak stomatal response with small magnitude and slow speed, whereas stomatal closure in the mpk12-4 and ca1ca4 plants was slower than in WT, but larger in magnitude (Fig. 2, A, C, and E). Very weak responses of the HT1 mutants suggest that HT1 is necessary to initiate stomatal responses to CO2 concentration changes in both ambient and sub-ambient levels. MPK12 and carbonic anhydrase mutants display stronger stomatal responses in the above-ambient than sub-ambient CO2 range (Fig. 2, A, C, and E), suggesting that respective signaling components have a more important role in sub-ambient compared with the above-ambient CO2 concentrations.

Figure 2. CO2 pathway mutants retain response patterns under regular and elevated vapor pressure deficit (VPD) conditions. A and B) Stomatal response to CO2 concentration changes from 400 to 100 parts per million (ppm), 100 to 400 ppm, and 400 to 800 ppm under regular (A) and elevated (B) VPD conditions, mean stomatal conductance ± SEM is shown. C, D) Boxplot of stomatal conductance (gs) change (mmol m−2 s−1) in response to CO2 concentration changes from 400 to 100 ppm, 100 to 400 ppm, and 400 to 800 ppm, respectively. E, F) Boxplot of 75% response time (min) of stomatal response to CO2 concentration changes from 400 to 100 ppm, 100 to 400 ppm, and 400 to 800 ppm, respectively. C to F) Boxes represent 25% to 75% quartiles and median as the horizontal lines, whiskers indicate the smallest and largest values, and points show individual plant values. Statistically significantly different groups are marked with different letters (One-way ANOVA with Tukey post hoc test, P < 0.05). A to F) Sample size was 5 for all plant lines. VPD during experiments was 0.9 kPa in (A, C, E) and 2.3 kPa in (B, D, F). Start of the first treatment was between 11:30 and 12:30. Col-0 data are the same as used in Fig. 1, and experiments with Col-0 and mutant lines were done together.

Elevated VPD led to lower stomatal conductance in all studied mutants (Fig. 2, A and B), and thus smaller magnitudes of CO2-induced changes in stomatal conductance (Fig. 2, C and D). Elevated VPD had some effects on the patterns of CO2-responsiveness in the CO2-signalling mutants (Fig. 2, A to F). The ca1ca4 and mpk12-4 stomatal response times were similar to WT in the 400–800 transition under elevated VPD (Fig. 2F), and the mpk12-4 mutant had a significantly slower than wild-type response to the 100–400 transition only under elevated VPD. Thus, the relatively larger degree of importance of MPK12 and carbonic anhydrases in regulating stomatal responses in the sub-ambient CO2 ranges was more pronounced under elevated VPD conditions.

Stomatal opening response to 400–100 transition in all of the studied CO2-signaling mutants was lower in magnitude and slower in response time, irrespective of VPD, whereas the difference in response time between wild-type and mutants was larger under elevated VPD (Fig. 2, A to F). Thus, stomatal opening response to sub-ambient [CO2] involves all of the CO2-signaling pathway components represented in this study, including the carbonic anhydrases CA1 and CA4, and the HT1 and MPK12 kinases.

SLAC1 and GHR1 are more important for above-ambient CO2-induced stomatal closure

The SLAC1 anion channel is a major component in the activation of stomatal closure response. Thus, we examined CO2 responses across different CO2 concentration ranges in plants deficient in the SLAC1 anion channel (slac1-3) and in OST1 or GHR1: proteins involved in SLAC1 activation (ost1-3, ghr1-3; Fig. 3). Stomatal 75% response time of slac1-3 plants was longer in both, 100–400 and 400–800 CO2 transitions (Fig. 3, A and E) but the latter was more affected as the 75% response time difference compared with wild-type plants was notably larger in the 400–800 transition. Yet, the magnitude of stomatal closure between wild-type Col-0 and slac1-3 was only different in the 100–400 transition (Fig. 3C), while their 400–800 response magnitude was similar. Under elevated VPD, slac1-3 stomatal response amplitude was similar to wild-type both during the 100–400 and 400–800 transitions, although stomatal response times of slac1-3 were still significantly longer than wild-type on both transitions (Fig. 3, D and F). Together, these results show that SLAC1 is important in stomatal closure in both sub-ambient and above-ambient [CO2] ranges, but response speed tends to be more severely impacted in the 400–800 transition, suggesting a more prominent role for SLAC1 in ensuring fast above-ambient CO2-induced stomatal closure.

Figure 3. CO2-response patterns in anion channel activation mutants are affected by vapor pressure deficit (VPD) conditions. A and B) Stomatal response to CO2 concentration changes from 400 to 100 parts per million (ppm), 100 to 400 ppm, and 400 to 800 ppm in regular (A) and elevated (B) VPD conditions, mean stomatal conductance ± SEM is shown. C, D) Boxplot of stomatal conductance (gs) change (mmol m−2 s−1) in response to CO2 concentration changes from 400 to 100 ppm, 100 to 400 ppm, and 400 to 800 ppm, respectively. E, F) Boxplot of 75% response time (min) of stomatal response to CO2 concentration changes from 400 to 100 ppm, 100 to 400 ppm, and 400 to 800 ppm, respectively. C to F) Boxes represent 25% to 75% quartiles and median as the horizontal lines, whiskers indicate the smallest and largest values, and points show individual plant values. Statistically significantly different groups are marked with different letters (One-way ANOVA with Tukey post hoc test, P < 0.05). A to F) Sample size was 5 for all plant lines. VPD during experiments was 0.9 kPa in (A, C, E) and 2.3 kPa in (B, D, F). Start of the first treatment was between 11:30 and 12:30. Col-0 data are the same as used in Fig. 1, and experiments with Col-0 and mutant lines were done together.

Plants deficient in SLAC1 activation via OST1 or GHR1 also showed impaired CO2-induced stomatal closure both in the 100–400 and in the 400–800 transitions (Fig. 3, A, C, and E). The ost1-3 mutant had long stomatal response times in both the 100–400 and 400–800 transitions, whereas the ghr1-3 response was similar to wild type in the 100–400 transition, but slower and very small in magnitude during the 400–800 transition (Fig. 3, A, C and E). Thus, both SLAC1-activating proteins are involved in CO2 responses in all tested concentration ranges, but GHR1, like SLAC1, appears to contribute more toward above-ambient CO2-induced stomatal closure.

As in other mutants (Fig. 2), elevated VPD lowered stomatal conductance and stomatal response magnitude, but tended to increase stomatal opening speed (Fig. 3). The ost1-3 75% response times and magnitude remained similar in both VPD conditions, in line with its VPD-insensitivity (Fig. 3, C to F, Merilo et al. 2013). In ghr1-3, stomatal response to 400–800 was either absent or extremely weak under elevated VPD (Fig. 3, B and D).

Sub-ambient CO2-induced stomatal opening response magnitude was similar to wild-type in all of the studied anion channel activation mutants, irrespective of VPD (Fig. 3, C and D). Stomatal 75% response time under regular VPD was similar to wild-type in all but ost1-3 (Fig. 3E), whereas under elevated VPD conditions, both ghr1-3 and ost1-3 had longer stomatal opening 75% response times (Fig. 3F). Therefore, regulation of SLAC1 is less important for sub-ambient CO2-induced stomatal opening than the guard cell CO2-specific signaling pathway. Only under elevated VPD, the 400–100 transition response was slower compared to the wild-type in ghr1-3 (Fig. 3F), indicating an interaction between CO2 and VPD signaling in stomatal opening responses.

Shifting CO2 levels from 100 to 800 ppm masks the differences between stomatal behavior in 100–400 and 400–800 ppm [CO2] transitions

To further study how stomatal movements differ depending on CO2 concentrations, we did additional gas-exchange measurements, during which CO2 concentration was changed directly from sub-ambient (100 ppm) to above-ambient (800 ppm). Previously, in Figs. 2 and 3 we saw different stomatal response characteristics for 100–400 and 400–800 CO2 transitions for different mutants that we failed to detect during the 100–800 CO2 transition (Fig. 4). For example, mpk12-4 stomatal response was small in amplitude in the 100–400 transition, but much larger in the 400–800 transition (Fig. 2, A and C), whereas in the 100–800 transition the mpk12-4 stomatal response amplitude was similar to wild type (Fig. 4A). Similarly, in the ca1ca4 plants 100–400 stomatal response was smaller and 400–800 response larger than in wild-type in amplitude (Fig. 2C), yet the 100–800 response magnitude was similar to wild-type (Fig. 4C). The ghr1-3 mutant had slower response to the 400–800 transition (Fig. 3E), whereas its 75% response time was similar to wild-type in the 100–800 transition (Fig. 4F). These results demonstrate that differences in stomatal responses between plant lines can remain undiscovered depending on CO2 concentrations that are used for experiments.

Figure 4. CO2 transition from 100 to 800 parts per million (ppm) masks different responses present in 100 to 400 and 400 to 800 ppm CO2 transitions. A and B) Stomatal response to CO2 concentration changes from 400 to 100 ppm and 100 to 800 ppm, mean stomatal conductance ± SEM is shown. C, D) Boxplot of stomatal conductance (gs) change (mmol m−2 s−1) in response to CO2 concentration changes from 400 to 100 ppm and 100 to 800 ppm, respectively. E, F) Boxplot of 75% response time (min) of stomatal response to CO2 concentration changes from 400 to 100 ppm and 100 to 800 ppm, respectively. C to F) Boxes represent 25% to 75% quartiles and median as the horizontal lines, whiskers indicate the smallest and largest values, and points show individual plant values. Statistically significantly different groups are marked with different letters (One-way ANOVA with Tukey post hoc test, P < 0.05). A to F) Sample size was 14 for Col-0; 6 for ht1–2; 7 for ht1-8D; 8 for ghr1-3, mpk12-4 and ost1-3; 9 for slac1-3 and ca1ca4. VPD during experiments was 0.9 kPa in (A to F). Start of the first treatment was between 11:30 and 12:30. Col-0 data are the same as used in Fig. 1, and experiments with Col-0 and mutant lines were done together.

Similar to CO2-induced stomatal closure experiments, [CO2] ranges were important also in stomatal opening assays (Figs. 2–3 vs Fig. 5). For example, response to the 400–100 transition in mpk12-4 was slower and smaller in magnitude compared with wild-type plants (Fig. 2, C and E), but in the 800–100 experiments, mpk12-4 had normal response amplitude and wild-type-like 75% response time (Fig. 5, C and E). The slac1-3 plants had similar to wild-type 400–100 stomatal opening speed and magnitude (Fig. 3, C and E), but were slower in the 800–100 response (Fig. 5F). We also observed a significantly slower opening response to the 800–400 transition in slac1-3 (Supplementary Fig. S2F), which indicates that the reduced stomatal opening speed of the 800–100 response in this mutant is caused by slower opening in the 800–400 range. Thus, although wild-type plants have no discernible differences between the 400–100 and the 800–100 stomatal opening responses (Fig. 1, A, B, E, and F), the molecular mechanisms are at least partly different for these CO2 transitions.

Figure 5. Stomatal opening in 800 to 100 parts per million (ppm) CO2 transition masks different responses present in sub-ambient and above-ambient CO2 concentrations. A and B) Stomatal response to CO2 concentration changes from 400 to 800 ppm and 800 to 100 ppm, mean stomatal conductance ± SEM is shown. C, D) Boxplot of stomatal conductance (gs) change (mmol m−2 s−1) in response to CO2 concentration changes from 400 to 800 ppm and 800 to 100 ppm, respectively. E, F) Boxplot of 75% response time (min) of stomatal response to CO2 concentration changes from 400 to 800 ppm and 800 to 100 ppm, respectively. C to F) Boxes represent 25% to 75% quartiles and median as the horizontal lines, whiskers indicate the smallest and largest values, and points show individual plant values. Statistically significantly different groups are marked with different letters (One-way ANOVA with Tukey post hoc test, P < 0.05). A to F) Sample size was 5 for all plant lines. VPD during experiments was 0.9 kPa in (A to F). Start of the first treatment was between 11:30 and 12:30. Col-0 data are the same as used in Fig. 1, experiments with Col-0 and mutant lines were done together.

We combined information from previously analyzed CO2 transitions in the mutants used in our study with our results (Table 1). Our findings mostly confirm previous results, where available, with the exception of some differences in response amplitude in the slac1-3, ost1-3, and ca1ca4 mutants that can be explained by longer treatment duration in this work that allowed the slower stomatal responses of these mutants to reach amplitudes similar to wild-type plants. We also found a faster stomatal response to sub-ambient to above-ambient CO2 transition in the ht1-2, potentially explained by different parameters used to describe response speed in different studies. Our experiments add information on the ambient to sub-ambient, above-ambient to ambient, sub-ambient to above-ambient, above-ambient to sub-ambient, and sub-ambient to ambient CO2 transitions that have not been systematically addressed before in all the mutants analyzed here.

Table 1. Stomatal CO2 responses for the studied mutants from the previously published studies (Hashimoto et al. 2006 [1]; Vahisalu et al. 2008 [2]; Hu et al. 2010 [3]; Xue et al. 2011 [4]; Laanemets et al. 2013 [5]; Merilo et al. 2013 [6]; Hu et al. 2015 [7]; Matrosova et al. 2015 [8]; Hashimoto-Sugimoto et al. 2016 [9]; Hõrak et al. 2016 [10]; Jakobson et al. 2016 [11]; Sierla et al. 2018 [12]; Takahashi et al. 2022 [13]; Yeh et al. 2023 [14]) and from this study

Plant line	Ambient to sub-ambient	Sub-ambient to ambient	Ambient to above-ambient	Sub-ambient to above-ambient	Above-ambient to ambient	Above-ambient to sub-ambient	
Amplitude	Speed	Amplitude	Speed	Amplitude	Speed	Amplitude	Speed	Amplitude	Speed	Amplitude	Speed	
ht1-2	+[1,10]/+	+[1,10]/+	na/+	na/−	+[1,8–10]/+	+[1,8–10]/+	+[4]/+	+[4]/−	na/+	na/−	+[8,9]/+	+[8,9]/+	
ht1-8D	+[10]/+	+[10]/+	na/+	na/+	+[10,13]/+	+[10,13]/+	na/+	na/+	na/+	na/+	+[13]/+	+[13]/+	
mpk12-4	+[14]/+	+[14]/+	na/+	na/−	+[11,14]/+	+[11,14]/+	na/−	na/+	+[14]/+	+[14]/+	−[13]/−	−[13]/−	
ca1ca4	+[1,7,8]/+	+[3,7,8]/+	na/+	na/+	+[3,7]/+	+[3,7]/+	+[3,7]/−	+[3,7]/+	na/+	na/+	+[3,7]/+	+[3,7]/+	
slac1-3	−[5]/−	−[5]/−	na/+	na/+	+[2,6]/−	+[2,6]/+	na/+	na/+	+[2]/−	+[2]/+	na/−	na/+	
ghr1-3	na/−	na/−	na/+	na/−	+[12]/+	+[12]/+	na/+	na/−	na/+	na/−	na/−	na/−	
ost1-3	+[4]/−	+[4]/+	na/+	na/+	+[4,6]/−	+[4,6]/+	+[4]/+	+[4]/+	+[6]/+	+[6]/+	+[4]/−	+[4]/+	
Published results/our results, “+” means different from wild type and “−” means not different from wild type.

Discussion

Here, we show that stomatal closure in response to an increase in CO2 concentration, which occurs both at sub-ambient and at above-ambient CO2 concentration ranges, is regulated by both the CO2-specific and SLAC1-related pathways under all CO2 ranges, but these pathways have a different degree of importance under different CO2 concentration ranges. Stomatal closure in response to the 400–800 CO2 concentration range is faster (Fig. 1B), while the 100–400 response has greater amplitude (Figs. 1A and 2C). These differences may be explained by higher guard cell volume and larger stomatal apertures under lower CO2 levels, leading to slower responses with larger amplitudes, similar to slower responses of larger stomata (Drake et al. 2013; Kübarsepp et al. 2020). However, the relative contribution of the components involved in stomatal CO2-signalling is also different in the 400–800 and 100–400 stomatal closure responses (Figs. 2 to 3). Increasing CO2 abruptly from 100 to 800 masked the presence of two processes with different kinetics (Figs. 2 to 4), indicating the necessity to analyze stomatal closure responses to elevated CO2 levels separately in the ambient to above-ambient and sub-ambient to ambient CO2 concentration ranges.

The HT1 kinase is required for plant stomatal CO2 signaling; plants with impaired HT1 function are nearly insensitive to all CO2 concentration changes (Hashimoto et al. 2006; Hashimoto-Sugimoto et al. 2016; Hõrak et al. 2016). HT1 together with MPK12 or MPK4 forms a primary CO2 sensing complex, where CO2/bicarbonate triggers the interaction of MPKs with HT1, and this leads to inhibition of the HT1 kinase activity (Takahashi et al. 2022). In our experiments, mpk12-4 plants had disrupted stomatal response amplitude to the 100–400 CO2 transition, while the 400–800 amplitude was unaffected (Fig. 2, A and C). In Tõldsepp et al (2018)mpk12 mpk4GC double-mutants, where MPK4 expression is suppressed only in guard cells, had no 400–800 CO2 response, yet mpk4GC single-mutants responded to the 400–800 CO2 transition similar to wild-type plants. Takahashi et al. (2022) showed MPK12 mutants to have the 400–800 response similar to the current study, although neither of these studies tested the 100–400 response. Thus, it seems that the MPK12/MPK4–HT1 complex largely loses functionality if HT1 is impaired (Fig. 2), same happens if both MPK12 and MPK4 are missing from guard cells (Tõldsepp et al. 2018), but losing only MPK12 preferentially affects the 100–400 CO2 transition (Fig. 2). Therefore, while in the 400–800 transition, the lack of MPK12 is likely compensated by MPK4, and MPK4 cannot effectively replace the function of MPK12 at sub-ambient CO2 levels. This could mean that while MPK12 and MPK4 both can form a CO2/bicarbonate sensing complex with HT1, their affinity for CO2/bicarbonate may be different.

CA1 and CA4 also affected stomatal responsiveness more in the 100–400 CO2 range, with slow and shallow stomatal response in the ca1ca4 mutant (Fig. 2). The strong stomatal closure in above-ambient 400–800 CO2 transition in this mutant (Fig. 2, Supplementary Fig. S2, A, C and E) might be related to increased autonomous CO2 conversion to HCO3− in the elevated CO2 environment due to shifting of the reaction balance toward bicarbonate production under above-ambient CO2 levels or by the involvement of other carbonic anhydrases, as recently demonstrated (Sun et al. 2022).

SLAC1 is important for stomatal closure responses to both elevated CO2 transitions, but compared with wild type, the response speed of slac1–3 in the 400–800 transition was more affected than in the 100–400 transition (Fig. 3, A and E). In slac1–3 plants, there was small in magnitude and slow stomatal closure and opening in the 400–800–400 CO2 transitions, yet stomatal opening in response to the 400–100 transition was strong and similar in speed to wild-type plants (Supplementary Fig. S2, B, D and F), further supporting the major role of SLAC1 in stomatal responses under ambient to above-ambient CO2 concentration changes. In the 100–400 CO2 transition, stomatal responsiveness may be partly compensated by other ion channels that are functional in slac1–3. In addition to the S-type anion channels like SLAC1, stomatal closure is also affected by R-type anion channels, such as QUICK-ACTIVATING ANION CHANNEL 1 (QUAC1, Meyer et al. 2010; Imes et al. 2013). Jalakas et al. (2021a) demonstrated that the quac1-1 slac1-3 double mutant and quac1-1 slac1-3 slah1-3 triple-mutant had no response to the 400–800 CO2 transition, while individually slac1-3 had a weak response and quac1-1 stomatal response was similar to wild-type plants. SLAC1 HOMOLOGUE 3 (SLAH3) is another S-type anion channel contributing to stomatal closure (Zhang et al. 2016) and could potentially compensate for the lack of SLAC1 in the 100–400 CO2 response, although it does not affect the 400–800 stomatal CO2 response (Jalakas et al. 2021a). Future studies should address the 100–400 CO2 response in mutants deficient in major S- and R-type anion channels to better understand their potential role in the 100–400 CO2 response.

Stomatal opening triggered by decreased CO2 levels, from 400 to 100 and from 800 to 100, had larger magnitude than the 800–400 response (Fig. 1, A, C, E, and G). Responses to 100 ppm final CO2 concentration had a consistently large amplitude and slow response rate in different experimental setups and were not affected by the order of CO2 treatments. Thus, stomatal opening in response to sub-ambient CO2 concentrations is a very prominent response. This is in line with Merilo et al. (2014), where two stimuli with an opposing effect on stomata were simultaneously applied, such as darkness and low CO2, or low CO2 and elevated VPD. In such combinations, Arabidopsis stomata always opened in response to sub-ambient CO2 levels, further indicating that reduction of CO2 is a strong and prevailing signal. The sub-ambient CO2-induced stomatal opening, similar to the light-induced opening, likely involves H+ ATPase activation (Inoue and Kinoshita 2017), changes in sugar and starch metabolism (Flütsch and Santelia 2021), and suppression of stomatal closure via e.g. inhibiting anion channel activation (Marten et al. 2007). The combined activation of all these processes may explain the slow stomatal opening rate in response to sub-ambient CO2 levels.

Elevated VPD triggers ABA biosynthesis (McAdam et al. 2016), which can potentially increase stomatal responsiveness to CO2 due to interactions of CO2 and ABA signaling pathways (Raschke 1975; Merilo et al. 2013; Chater et al. 2015). In previous experiments, elevated VPD has been shown to either increase stomatal responsiveness to elevated CO2, potentially via enhanced ABA levels (Bunce 1998), or decrease it, potentially due to reduced stomatal apertures under elevated VPD (Morison and Gifford 1983; Talbott et al. 2003). Our results showed an enhanced stomatal response speed in the 400–100, 100–400, and 800–100 CO2 concentration transitions but not in the already faster 400–800 transition (Fig. 1, A, B, E, and F). Elevated VPD has been shown to accelerate stomatal opening in light in angiosperms due to reduced back-pressure of epidermal cells on guard cells (Mott et al. 1999; Pichaco et al. 2024), our results of faster stomatal opening in response to sub-ambient [CO2] under elevated VPD are in line with this (Fig. 1, B and E). Additionally, elevated VPD also accelerated stomatal closure responses to elevated CO2 levels in the sub-ambient to ambient concentration ranges in wild-type Arabidopsis (Fig. 1B). Faster stomatal closure may be caused by increased ABA levels under elevated VPD conditions (McAdam and Brodribb 2015) or may be explained by the smaller steady-state stomatal conductance caused by smaller stomatal apertures under elevated VPD that can adjust faster in response to environmental changes.

Steady-states of stomatal conductance of all plant lines were decreased by elevated VPD (Figs. 2, A and B, 3, A and B), confirming that VPD is an important factor for steady-state stomatal conductance (Grossiord et al. 2020; López et al. 2021). Mutants with disrupted stomatal ABA response, ghr1-3, ost1-3, and slac1-3, also had lower steady-state stomatal conductances under elevated VPD (Fig. 3, A and B). Their high VPD-induced decrease in steady-state stomatal conductance could be caused by ABA-independent active processes (e.g. ABA-independent OST1 activation, whose contribution to stomatal closure under high VPD increases in time according to Jalakas et al. (2021b)), or hydropassive stomatal closure. These data indicate that elevated ABA levels alone are not sufficient to explain the decrease of stomatal conductance under elevated VPD conditions, as discussed before (Wang et al. 2017; Merilo et al. 2018; Yaaran et al. 2019). CO2-induced stomatal responses under elevated VPD were sometimes faster and mostly had a smaller amplitude (Figs. 2, B, D and F, 3, B, D and F). However, the response to the 400–800 CO2 transition disappeared completely in the ghr1-3 plants under elevated VPD conditions (Fig. 3, B and D). GHR1 contributes to stomatal closure in response to both elevated VPD (Hsu et al. 2021) and CO2 (Hõrak et al. 2016; Sierla et al. 2018). Thus, CO2 and VPD responses may interact in ghr1-3: if the relatively small elevated VPD-induced stomatal closure already occurred in the ghr1-3 plants subjected to elevated VPD (Sierla et al. 2018; Hsu et al. 2021), no further response to CO2 elevation was triggered.

Here we have focused on guard cell-specific stomatal CO2-signaling components and their different contribution to stomatal closure responses in the sub-ambient to ambient and ambient to above-ambient CO2 concentration transitions. In addition to these components, it is likely that signal mediators outside guard cells also contribute to different stomatal CO2-response characteristics under different CO2 levels. While guard cells in isolated epidermis can respond to elevated CO2 levels by a decrease in stomatal aperture (Webb et al. 1996; Chater et al. 2015), signals from mesophyll are needed for strong CO2-induced stomatal closure (Mott et al. 2008; Fujita et al. 2013). Mesophyll processes, e.g. photosynthesis and sugar metabolism, are known to impact stomatal behavior (Lawson and Matthews 2020), and their contribution to different stomatal CO2-response patterns under different CO2 concentration ranges merits further study.

Elevated temperatures caused by climate change increase the evaporative demand of the atmosphere manifested as higher VPD levels, which increases transpiration and triggers stomatal closure to avoid wilting. Together, changes in atmospheric VPD and CO2 levels are among the greatest agricultural challenges of the future, yet how these stimuli together affect plant stomatal behavior is poorly understood. Here, we show that while elevated VPD negatively affects steady-state stomatal conductances, it has little effect on stomatal CO2-responsiveness. Nevertheless, in some genetic backgrounds, we found an interaction between CO2 and VPD treatments, indicating that the simultaneous effects of these factors on stomatal behavior merit further study. We also show that stomatal response to elevated CO2 has different kinetics under sub-ambient and above-ambient CO2 concentration ranges and its known regulators contribute to a different degree under different CO2 concentration transitions. Thus, to better understand stomatal responses to CO2 it is necessary to carefully consider CO2 levels and experimental setups.

Materials and methods

Plant lines and growth conditions

Arabidopsis (Arabidopsis thaliana) accession Col-0 and the following mutants in the same genetic background were used for experiments: slac1-3 (Vahisalu et al. 2008), ost1-3 (Yoshida et al. 2002), ghr1-3 (Sierla et al. 2018), ht1-2 (Hashimoto et al. 2006), ht1-8D (Hõrak et al. 2016), mpk12-4 (Jakobson et al. 2016), and ca1ca4 (Hu et al. 2010). Plants were grown in 4:2:3 v/v peat:vermiculite:water mixture at 12/12 photoperiod with 150 µmol m−2 s−1 light in controlled-environment growth cabinets (AR-66LX; Percival Scientific; MCA1600, Snijders Scientific) at 70% relative humidity and day-time temperature of 23 °C (VPD 0.84 kPa) and nighttime temperature 18 °C (VPD 0.62 kPa). Plant age at experiment time was ∼25 days.

Gas-exchange measurements

Measurements of stomatal conductance were carried out with a temperature-controlled custom-built gas-exchange device (Kollist et al. 2007; Hõrak et al. 2017). Plants were inserted into measurement cuvettes and allowed to acclimate for 1 to 2 h at normal air humidity (VPD, 0.9 kPa) or at lower air humidity (VPD, 2.3 kPa), 24 °C air temperature and 400 ppm CO2. Experiments with various CO2 transitions were carried out and in some cases both under normal and at high VPD conditions, as shown in Fig. 1. The first CO2 treatment was applied approximately at noon (11:30 to 12:30), stomatal conductance was always followed for 2 h under each treatment (Fig. 1).

Data analysis

Magnitude of stomatal closure response was calculated as the absolute difference in stomatal conductance between last time point before treatment and at the end of 2 h of treatment at a given CO2 level (displayed as “Change in gs”). To describe the kinetic characteristics of stomatal responses, we identified the difference between stomatal conductance at the last time point before treatment and at the time point when maximal change in stomatal conductance had occurred for any given response as the total 100% response and calculated the time when 75% of the total stomatal response was achieved (see also illustration in Supplementary Fig. S1). The 75% response time describes the overall response speed across the 2-h treatment time but not the initial kinetics of the stomatal responses. One-way ANOVA with Tukey post hoc test was used for statistical analyses as indicated in the figure legends, and P < 0.05 was considered statistically significant. Statistical analyses were carried out using Past 4.0 (Hammer et al. 2001) and Statistica 7.1 (Stat. Soft. Inc).

Accession numbers

AGI accession numbers for genes studied in this article are AT1G12480 (SLAC1), AT3G01500 (CA1), AT1G70410 (CA4), AT2G46070 (MPK12), AT1G62400 (HT1), AT4G33950 (OST1), and AT4G20940 (GHR1).

Supplementary data

The following materials are available in the online version of this article.

Supplementary Figure S1 . Scheme explaining the calculation process for 75% stomatal response time using Col-0 and ghr1-3 ambient to above-ambient CO2 stomatal responses as an example.

Supplementary Figure S2 . Plant stomatal responses to 400–800–400–100 parts per million (ppm) CO2 concentration changes.

Supplementary Material

kiae320_Supplementary_Data

Funding

This work was supported by the Estonian Research Council grant PSG404 to H.H., Basic funding from the Institute of Technology, PRG719, and PRG1620 to E.M. and PRG433 and PRG2004 to H.K.; and Estonian Research Infrastructure Roadmap project TAIM; and European Regional Development Fund via Center of Excellence in Molecular Cell Engineering.

Data availability

The data underlying this article will be shared on reasonable request to the corresponding author.

Dive Curated Terms

The following phenotypic, genotypic, and functional terms are of significance to the work described in this paper:

MPK12 Gramene: AT2G46070

MPK12 Araport: AT2G46070

OST1 Gramene: AT4G33950

OST1 Araport: AT4G33950

SLAC1 Gramene: AT1G12480

SLAC1 Araport: AT1G12480

MPK4 Gramene: AT4G01370

MPK4 Araport: AT4G01370

stomatal opening Planteome: TO:0020097

CBC2 Gramene: AT5G50000

CBC2 Araport: AT5G50000

SLAH3 Gramene: AT5G24030

SLAH3 Araport: AT5G24030

ABA CHEBI: CHEBI:2365

GHR1 Gramene: AT4G20940

GHR1 Araport: AT4G20940

CBC1 Gramene: AT3G01490

CBC1 Araport: AT3G01490

HT1 Gramene: AT1G62400

HT1 Araport: AT1G62400
==== Refs
References

Azoulay-Shemer T , PalomaresA, BagheriA, Israelsson-NordstromM, EngineerCB, BargmannBOR, StephanAB, SchroederJI. Guard cell photosynthesis is critical for stomatal turgor production, yet does not directly mediate CO2- and ABA-induced stomatal closing. Plant J. 2015:83 (4 ):567–581. 10.1111/tpj.12916 26096271
Brandt B , BrodskyDE, XueS, NegiJ, IbaK, KangasjärviJ, GhassemianM, StephanAB, HuH, SchroederJI. Reconstitution of abscisic acid activation of SLAC1 anion channel by CPK6 and OST1 kinases and branched ABI1 PP2C phosphatase action. Proc Natl Acad Sci U S A. 2012:109 (26 ):10593–10598. 10.1073/pnas.1116590109 22689970
Brodribb TJ , McAdamSAM, JordanGJ, FeildTS. Evolution of stomatal responsiveness to CO2 and optimization of water-use efficiency among land plants. New Phytol. 2009:183 (3 ):839–847. 10.1111/j.1469-8137.2009.02844.x 19402882
Bunce JA . Effects of humidity on short-term responses of stomatal conductance to an increase in carbon dioxide concentration. Plant Cell Environ. 1998:21 (1 ):115–120. 10.1046/j.1365-3040.1998.00253.x
Chater C , PengK, MovahediM, DunnJA, WalkerHJ, LiangY-K, McLachlanDH, CassonS, IsnerJC, WilsonI, et al Elevated CO2-induced responses in stomata require ABA and ABA signaling. Curr Biol. 2015:25 (20 ):2709–2716. 10.1016/j.cub.2015.09.013 26455301
Cutler SR , RodriguezPL, FinkelsteinRR, AbramsSR. Abscisic acid: emergence of a core signaling network. Annu Rev Plant Biol. 2010:61 (1 ):651–679. 10.1146/annurev-arplant-042809-112122 20192755
Drake PL , FroendRH, FranksPJ. Smaller, faster stomata: scaling of stomatal size, rate of response, and stomatal conductance. J Exp Bot. 2013:64 (2 ):495–505. 10.1093/jxb/ers347 23264516
Flütsch S , SanteliaD. Mesophyll-derived sugars are positive regulators of light-driven stomatal opening. New Phytol. 2021:230 (5 ):1754–1760. 10.1111/nph.17322 33666260
Franks PJ , Britton-HarperZJ. No evidence of general CO2 insensitivity in ferns: one stomatal control mechanism for all land plants? New Phytol. 2016:211 (3 ):819–827. 10.1111/nph.14020 27214852
Fujita T , NoguchiK, TerashimaI. Apoplastic mesophyll signals induce rapid stomatal responses to CO2 in Commelina communis. New Phytol. 2013:199 (2 ):395–406. 10.1111/nph.12261 23560389
Grossiord C , BuckleyTN, CernusakLA, NovickKA, PoulterB, SiegwolfRTW, SperryJS, McDowellNG. Plant responses to rising vapor pressure deficit. New Phytol. 2020:226 (6 ):1550–1566. 10.1111/nph.16485 32064613
Hammer Ø , HarperDAT, RyanPD. PAST: paleontological statistics software package for education and data analysis. Palaeontol Electron. 2001:4 (1 ):9pp. http://palaeo-electronica.org/2001_1/past/issue1_01.htm
Hashimoto-Sugimoto M , NegiJ, MondaK, HigakiT, IsogaiY, NakanoT, HasezawaS, IbaK. Dominant and recessive mutations in the Raf-like kinase HT1 gene completely disrupt stomatal responses to CO2 in Arabidopsis. J Exp Bot. 2016:67 (11 ):3251–3261. 10.1093/jxb/erw134 27034327
Hashimoto M , NegiJ, YoungJ, IsraelssonM, SchroederJI, IbaK. Arabidopsis HT1 kinase controls stomatal movements in response to CO2. Nat Cell Biol. 2006:8 (4 ):391–397. 10.1038/ncb1387 16518390
Hayashi M , SugimotoH, TakahashiH, SekiM, ShinozakiK, SawasakiT, KinoshitaT, InoueS. Raf-like kinases CBC1 and CBC2 negatively regulate stomatal opening by negatively regulating plasma membrane H+-ATPase phosphorylation in Arabidopsis. Photochem Photobiol Sci. 2020:19 (1 ):88–98. 10.1039/c9pp00329k 31904040
Hiyama A , TakemiyaA, MunemasaS, OkumaE, SugiyamaN, TadaY, MurataY, ShimazakiK. Blue light and CO2 signals converge to regulate light-induced stomatal opening. Nat Commun. 2017:8 (1 ):1284. 10.1038/s41467-017-01237-5 29101334
Hõrak H , KollistH, MeriloE. Fern stomatal responses to ABA and CO2 depend on species and growth conditions. Plant Physiol. 2017:174 (2 ):672–679. 10.1104/pp.17.00120 28351911
Hõrak H , SierlaM, TõldseppK, WangC, WangY-S, NuhkatM, ValkE, PechterP, MeriloE, SalojärviJ, et al A dominant mutation in the HT1 kinase uncovers roles of MAP kinases and GHR1 in CO2-induced stomatal closure. Plant Cell. 2016:28 (10 ):2493–2509. 10.1105/tpc.16.00131 27694184
Hsu P-K , TakahashiY, MeriloE, CostaA, ZhangL, KernigK, LeeKH, SchroederJI. Raf-like kinases and receptor-like (pseudo)kinase GHR1 are required for stomatal vapor pressure difference response. Proc Natl Acad Sci U S A. 2021:118 (47 ):e2107280118. 10.1073/pnas.2107280118
Hu H , Boisson-DernierA, Israelsson-NordströmM, BöhmerM, XueS, RiesA, GodoskiJ, KuhnJM, SchroederJI. Carbonic anhydrases are upstream regulators of CO2-controlled stomatal movements in guard cells. Nat Cell Biol. 2010:12 (1 ):87–93. 10.1038/ncb2009 20010812
Hu H , RappelW-J, OcchipintiR, RiesA, BöhmerM, YouL, XiaoC, EngineerCB, BoronWF, SchroederJI. Distinct cellular locations of carbonic anhydrases mediate carbon dioxide control of stomatal movements. Plant Physiol. 2015:169 (2 ):1168–1178. 10.1104/pp.15.00646 26243620
Hua D , WangC, HeJ, LiaoH, DuanY, ZhuZ, GuoY, ChenZ, GongZ. A plasma membrane receptor kinase, GHR1, mediates abscisic acid- and hydrogen peroxide-regulated stomatal movement in Arabidopsis. Plant Cell Online. 2012:24 (6 ):2546–2561. 10.1105/tpc.112.100107
Imes D , MummP, BöhmJ, Al-RasheidKAS, MartenI, GeigerD, HedrichR. Open stomata 1 (OST1) kinase controls R–type anion channel QUAC1 in Arabidopsis guard cells. Plant J. 2013:74 (3 ):372–382. 10.1111/tpj.12133 23452338
Inoue S , KinoshitaT. Blue light regulation of stomatal opening and the plasma membrane H+-ATPase. Plant Physiol. 2017:174 (2 ):531–538. 10.1104/pp.17.00166 28465463
Jakobson L , VaahteraL, TõldseppK, NuhkatM, WangC, WangY-S, HõrakH, ValkE, PechterP, SindarovskaY, et al Natural variation in Arabidopsis Cvi-0 accession reveals an important role of MPK12 in guard cell CO2 signaling. PLoS Biol. 2016:14 (12 ):e2000322. 10.1371/journal.pbio.2000322 27923039
Jalakas P , NuhkatM, VahisaluT, MeriloE, BroschéM, KollistH. Combined action of guard cell plasma membrane rapid- and slow-type anion channels in stomatal regulation. Plant Physiol. 2021a:187 (4 ):2126–2133. 10.1093/plphys/kiab202 34009364
Jalakas P , TakahashiY, WaadtR, SchroederJI, MeriloE. Molecular mechanisms of stomatal closure in response to rising vapour pressure deficit. New Phytol. 2021b:232 (2 ):468–475. 10.1111/nph.17592 34197630
Katsuhara M , HanbaYT. Barley plasma membrane intrinsic proteins (PIP aquaporins) as water and CO2 transporters. Pflüg Arch—Eur J Physiol. 2008:456 (4 ):687–691. 10.1007/s00424-007-0434-9
Kollist T , MoldauH, RasulovB, OjaV, RämmaH, HüveK, JaspersP, KangasjärviJ, KollistH. A novel device detects a rapid ozone-induced transient stomatal closure in intact Arabidopsis and its absence in abi2 mutant. Physiol Plant. 2007:129 (4 ):796–803. 10.1111/j.1399-3054.2006.00851.x
Kübarsepp L , LaanistoL, NiinemetsÜ, TaltsE, TosensT. Are stomata in ferns and allies sluggish? Stomatal responses to CO2, humidity and light and their scaling with size and density. New Phytol. 2020:225 (1 ):183–195. 10.1111/nph.16159 31479517
Laanemets K , WangY-F, LindgrenO, WuJ, NishimuraN, LeeS, CaddellD, MeriloE, BroscheM, KilkK, et al Mutations in the Slac1 anion channel slow stomatal opening and severely reduce K+ uptake channel activity via enhanced cytosolic [Ca2+] and increased Ca2+ sensitivity of K+ uptake channels. New Phytol. 2013:197 (1 ):88–98. 10.1111/nph.12008 23126621
Lawson T , MatthewsJ. Guard cell metabolism and stomatal function. Annu Rev Plant Biol. 2020:71 (1) :273–302. 10.1146/annurev-arplant-050718-100251 32155341
López J , WayDA, SadokW. Systemic effects of rising atmospheric vapor pressure deficit on plant physiology and productivity. Glob Change Biol. 2021:27 (9 ):1704–1720. 10.1111/gcb.15548
Marten H , HedrichR, RoelfsemaMRG. Blue light inhibits guard cell plasma membrane anion channels in a phototropin-dependent manner. Plant J. 2007:50 (1 ):29–39. 10.1111/j.1365-313X.2006.03026.x 17319842
Matrosova A , BogireddiH, Mateo-PeñasA, Hashimoto-SugimotoM, IbaK, SchroederJI, Israelsson-NordströmM. The Ht1 protein kinase is essential for red light-induced stomatal opening and genetically interacts with Ost1 in red light and Co2 -induced stomatal movement responses. New Phytol. 2015:208 (4 ):1126–1137. 10.1111/nph.13566 26192339
McAdam SAM , BrodribbTJ. The evolution of mechanisms driving the stomatal response to vapor pressure deficit. Plant Physiol. 2015:167 (3 ):833–843. 10.1104/pp.114.252940 25637454
McAdam SAM , SussmilchFC, BrodribbTJ. Stomatal responses to vapour pressure deficit are regulated by high speed gene expression in angiosperms. Plant Cell Environ. 2016:39 (3 ):485–491. 10.1111/pce.12633 26353082
Merilo E , JõesaarI, BroschéM, KollistH. To open or to close: species-specific stomatal responses to simultaneously applied opposing environmental factors. New Phytol. 2014:202 (2 ):499–508. 10.1111/nph.12667 24392838
Merilo E , LaanemetsK, HuH, XueS, JakobsonL, TulvaI, Gonzalez-GuzmanM, RodriguezPL, SchroederJI, BroschèM, et al PYR/RCAR receptors contribute to ozone-, reduced air humidity-, darkness-, and CO2-induced stomatal regulation. Plant Physiol. 2013:162 (3 ):1652–1668. 10.1104/pp.113.220608 23703845
Merilo E , YarmolinskyD, JalakasP, ParikH, TulvaI, RasulovB, KilkK, KollistH. Stomatal VPD response: there is more to the story than ABA. Plant Physiol. 2018:176 (1 ):851–864. 10.1104/pp.17.00912 28986421
Meyer S , MummP, ImesD, EndlerA, WederB, Al-RasheidKAS, GeigerD, MartenI, MartinoiaE, HedrichR. AtALMT12 represents an R-type anion channel required for stomatal movement in Arabidopsis guard cells. Plant J. 2010:63 (6 ):1054–1062. 10.1111/j.1365-313X.2010.04302.x 20626656
Morison JIL , GiffordRM. Stomatal sensitivity to carbon dioxide and humidity: a comparison of two C3 and two C4 grass Species. Plant Physiol. 1983:71 (4 ):789–796. 10.1104/pp.71.4.789 16662909
Mott KA , ShopeJC, BuckleyTN. Effects of humidity on light-induced stomatal opening: evidence for hydraulic coupling among stomata. J Exp Bot. 1999:50 (336 ):1207–1213. 10.1093/jxb/50.336.1207
Mott KA , SibbernsenED, ShopeJC. The role of the mesophyll in stomatal responses to light and CO2. Plant Cell Environ. 2008:31 (9 ):1299–1306. 10.1111/j.1365-3040.2008.01845.x 18541006
Negi J , MatsudaO, NagasawaT, ObaY, TakahashiH, Kawai-YamadaM, UchimiyaH, HashimotoM, IbaK. CO2 regulator SLAC1 and its homologues are essential for anion homeostasis in plant cells. Nature. 2008:452 (7186 ):483–486. 10.1038/nature06720 18305482
Pichaco J , ManandharA, McAdamSAM. Mechanical advantage makes stomatal opening speed a function of evaporative demand. Plant Physiol. 2024:195 (1 ):370–377. 10.1093/plphys/kiae023 38217870
Raschke K . Simultaneous requirement of carbon dioxide and abscisic acid for stomatal closing in Xanthium strumarium L. Planta. 1975:125 (3 ):243–259. 10.1007/BF00385601 24435438
Sierla M , HõrakH, OvermyerK, WaszczakC, YarmolinskyD, MaierhoferT, VainonenJP, SalojärviJ, DenessioukK, LaanemetsK, et al The receptor-like pseudokinase GHR1 is required for stomatal closure. Plant Cell. 2018:30 (11 ):2813–2837. 10.1105/tpc.18.00441 30361234
Sun P , IsnerJ-C, Coupel-LedruA, ZhangQ, PridgeonAJ, HeY, MenguerPK, MillerAJ, SandersD, McgrathSP, et al Countering elevated CO2 induced Fe and Zn reduction in Arabidopsis seeds. New Phytol. 2022:235 (5 ):1796–1806. 10.1111/nph.18290 35637611
Takahashi Y , BosmansKC, HsuP-K, PaulK, SeitzC, YehC-Y, WangY-S, YarmolinskyD, SierlaM, VahisaluT, et al Stomatal CO2/bicarbonate sensor consists of two interacting protein kinases, Raf-like HT1 and non-kinase-activity activity requiring MPK12/MPK4. Sci Adv. 2022:8 (49 ):eabq6161. 10.1126/sciadv.abq6161
Talbott LD , RahvehE, ZeigerE. Relative humidity is a key factor in the acclimation of the stomatal response to CO2. J Exp Bot. 2003:54 (390 ):2141–2147. 10.1093/jxb/erg215 12867546
Tõldsepp K , ZhangJ, TakahashiY, SindarovskaY, HõrakH, CeciliatoPHO, KoolmeisterK, WangY-S, VaahteraL, JakobsonL, et al Mitogen-activated protein kinases MPK4 and MPK12 are key components mediating CO2-induced stomatal movements. Plant J. 2018:96 (5 ):1018–1035. 10.1111/tpj.14087 30203878
Vahisalu T , KollistH, WangY-F, NishimuraN, ChanW-Y, ValerioG, LamminmäkiA, BroschéM, MoldauH, DesikanR, et al SLAC1 is required for plant guard cell S-type anion channel function in stomatal signalling. Nature. 2008:452 (7186 ):487–491. 10.1038/nature06608 18305484
Wang Y , HillsA, Vialet-ChabrandS, PapanatsiouM, GriffithsH, RogersS, LawsonT, LewVL, BlattMR. Unexpected connections between humidity and Ion transport discovered using a model to bridge guard cell-to-leaf scales. Plant Cell. 2017:29 (11 ):2921–2939. 10.1105/tpc.17.00694 29093213
Wang C , HuH, QinX, ZeiseB, XuD, RappelW-J, BoronWF, SchroederJI. Reconstitution of CO2 regulation of SLAC1 anion channel and function of CO2-permeable PIP2; 1 aquaporin as CARBONIC ANHYDRASE4 interactor. Plant Cell. 2016:28 (2 ):568–582. 10.1105/tpc.15.00637 26764375
Webb AAR , McAinshMR, MansfieldTA, HetheringtonAM. Carbon dioxide induces increases in guard cell cytosolic free calcium. Plant J. 1996:9 (3 ):297–304. 10.1046/j.1365-313X.1996.09030297.x
Xue S , HuH, RiesA, MeriloE, KollistH, SchroederJI. Central functions of bicarbonate in S-type anion channel activation and OST1 protein kinase in CO2 signal transduction in guard cell. EMBO J. 2011:30 (8 ):1645–1658. 10.1038/emboj.2011.68 21423149
Yaaran A , NeginB, MoshelionM. Role of guard-cell ABA in determining steady-state stomatal aperture and prompt vapor-pressure-deficit response. Plant Sci. 2019:281 :31–40. 10.1016/j.plantsci.2018.12.027 30824059
Yeh C-Y , WangY-S, TakahashiY, KuuskK, PaulK, ArjusT, YadlosO, SchroederJI, IlvesI, Garcia-SosaAT, et al MPK12 in stomatal CO2 signaling: function beyond its kinase activity. New Phytol. 2023:239 (1 ):146–158. 10.1111/nph.18913 36978283
Yoshida R , HoboT, IchimuraK, MizoguchiT, TakahashiF, AronsoJ, EckerJR, ShinozakiK. ABA-Activated SnRK2 protein kinase is required for dehydration stress signaling in Arabidopsis. Plant Cell Physiol. 2002:43 (12 ):1473–1483. 10.1093/pcp/pcf188 12514244
Yuan W , ZhengY, PiaoS, CiaisP, LombardozziD, WangY, RyuY, ChenG, DongW, HuZ, et al Increased atmospheric vapor pressure deficit reduces global vegetation growth. Sci Adv. 2019:5 (8 ):eaax1396. 10.1126/sciadv.aax1396 31453338
Zhang A , RenH-M, TanY-Q, QiG-N, YaoF-Y, WuG-L, YangL-W, HussainJ, SunS-J, WangY-F. S-Type anion channels SLAC1 and SLAH3 function as essential negative regulators of inward K+ channels and stomatal opening in Arabidopsis. Plant Cell. 2016:28 (4 ):949–965. 10.1105/tpc.15.01050 27002025
