1. Introduction
Stomata are the cellular pores on the leaf epidermis that allow plants to balance photosynthetic carbon dioxide (CO2) uptake with water vapor loss. Stomatal movements result from changes in turgor of stomatal cells (Jezek & Blatt, Reference Jezek and Blatt2017). Stomatal opening is induced by an increase of turgor pressure in guard cells (GCs), while a decrease of turgor pressure in GCs results in stomatal closure. To optimise gas exchange, stomata interpret and integrate a plethora of environmental cues such as light, humidity, temperature, CO2 concentration and even biotic factors like pathogens (Engineer et al., Reference Engineer, Hashimoto-Sugimoto, Negi, Israelsson-Nordström, Azoulay-Shemer, Rappel, Iba and Schroeder2016; Jezek & Blatt, Reference Jezek and Blatt2017; Kollist et al., Reference Kollist, Nuhkat and Roelfsema2014; Merilo et al., Reference Merilo, Jõesaar, Brosché and Kollist2014; Murata et al., Reference Murata, Mori and Munemasa2015; Sierla et al., Reference Sierla, Waszczak, Vahisalu and Kangasjärvi2016). In high light, for example, stomata of C3 and C4 plants open to provide sufficient CO2 for photosynthesis. In low light, on the other hand, less CO2 is required to saturate photosynthesis and, consequently, stomata close to limit water loss. Therefore, stomatal responsiveness and fast opening and closing kinetics can significantly contribute to plant water use efficiency (WUE) in changing environments (Lawson & Vialet-Chabrand, Reference Lawson and Vialet-Chabrand2019; McAusland et al., Reference McAusland, Vialet‐Chabrand, Davey, Baker, Brendel and Lawson2016). WUE represents the ratio of carbon assimilation and water loss and is a crucial trait for plant productivity and stress resilience (Leakey et al., Reference Leakey, Ferguson, Pignon, Wu, Jin, Hammer and Lobell2019; McAusland et al., Reference McAusland, Vialet‐Chabrand, Davey, Baker, Brendel and Lawson2016). Grasses, which include the cereals like rice, maize and wheat, show comparatively fast stomatal movements that likely contribute to more water-efficient gas exchange in changing environments (Franks & Farquhar, Reference Franks and Farquhar2007; Lawson & Matthews, Reference Lawson and Matthews2020; McAusland et al., Reference McAusland, Vialet‐Chabrand, Davey, Baker, Brendel and Lawson2016).
During the day plants face changing environmental conditions such as fluctuating ambient light intensity (Q out), temperature (T) and relative humidity (RH). Stomata mostly respond locally to environmental stimuli. This allows infrared gas analyser (IRGA)-based leaf gas exchange studies to be robust since leaves are placed in a chamber and exposed to controlled Q out, RH, T and CO2 concentration ([CO2]) regardless of the ambient conditions. Nevertheless, it has already been suggested that external ambient conditions might systemically affect local stomatal responses measured by IRGA systems (Devireddy et al., Reference Devireddy, Zandalinas, Gómez-Cadenas, Blumwald and Mittler2018; Devireddy et al., Reference Devireddy, Arbogast and Mittler2020;
Ehonen et al., Reference Ehonen, Holtta and Kangasjärvi2020). This might be particularly relevant for gas exchange studies that are performed in greenhouse or field settings with significant daily and seasonal environmental fluctuations. However, the putative systemic influence of the varying ambient conditions on gas exchange measurements is not generally accounted for.
Furthermore, gas exchange parameters such as carbon assimilation (A), stomatal conductance to water vapor (g sw), intrinsic water use efficiency (iWUE) and stomatal kinetics are influenced by anatomical traits such as stomatal density (SD) and stomatal length (SL) (Elliott-Kingston et al., Reference Elliott-Kingston, Haworth, Yearsley, Batke, Lawson and McElwain2016; Faralli et al., Reference Faralli, Matthews and Lawson2019; Haworth et al., Reference Haworth, Marino, Loreto and Centritto2021; Lawson & Blatt, Reference Lawson and Blatt2014). SD and SL are negatively correlated and vary in response to a variety of environmental conditions such as T, RH, [CO2] or Q out (Bertolino et al., Reference Bertolino, Caine and Gray2019; Franks et al., Reference Franks, Drake and Beerling2009; Zhang et al., Reference Zhang, Wang, Yang, Cui and Niu2021). The seasonal variation of environmental conditions might, therefore, affect the intraspecific plasticity of stomatal anatomical traits influencing gas exchange performance and, consequently, the results of long-term gas exchange phenotyping studies.
Here, we quantified stomatal conductance kinetics in 120 individuals of the grass model Brachypodium distachyon (Bd21-3) in a greenhouse over the course of 2 years. Simultaneously, we logged the environmental conditions in the greenhouse (Q out, T and RH) and time of the day (time) and quantified how these parameters affected the measured gas exchange parameters (A, g sw, iWUE and stomatal response kinetics). We additionally quantified anatomical traits of stomata (SD and SL) in three different seasons (summer, autumn and winter) and observed a significant impact of seasonal growth conditions on these traits. This allowed us to correlate how variations in SD and SL influence steady-state gas exchange, stomatal kinetics and maximum stomatal conductance (g smax). When calculating anatomical g smax based on anatomical traits, we realised that existing approaches to calculate maximum pore area for the double end-correction version of the equation by Franks and Farquhar (Reference Franks and Farquhar2001) did not sufficiently account for the graminoid morphology. Using quantitative morphometry of open and closed B. distachyon stomata we determined how to accurately calculate maximum pore area and pore depth. These adjustments allowed for an accurate prediction of physiological g smax based on stomatal anatomical traits in B. distachyon.
2. Material and methods
2.1. Plant material and growth conditions
B. distachyon Bd21-3 seeds were vernalised in water for 2 days at 4oC before being transferred to soil. Plants were grown in a greenhouse with 18 hr light:6 hr dark, average day temperature = 28oC, average night temperature = 25oC and average RH = 40%. We used 6 × 6 × 8 cm pots per plant filled with four parts soil (Einheitserde CL ED73) and one part vermiculite. The greenhouse is located at 49° 24’ 52.38’’ N and 8° 40’ 5.808’’ E at the Centre for Organismal Studies Heidelberg, Im Neuenheimer Feld 360, 69120 Heidelberg, Germany. Daily mean temperature (January 2019 to September 2021) varied between 3 and 5°C in December to February, 8–11°C in March to April, 13–21°C in May to September and 6–12°C in October to November (Deutscher Wetterdienst, https://cdc.dwd.de/portal/). Average daylight hours are 8.3–10.2 hr in December to February, 11.9–13.8 hr in March to April, 12.6–16.2 hr in May to September and 9.1–10-8 hr in October to November (Deutscher Wetterdienst, https://cdc.dwd.de/portal/).
2.2. Leaf-level gas exchange measurements
All measurements were performed on B. distachyon leaves 3 weeks after sowing using a LI-6800 Portable Photosynthesis System (Li-COR Biosciences Inc., Lincoln, NE, USA) equipped with a Multiphase Flash Fluorometer (6800-01A) chamber. The youngest, fully expanded leaf was measured using the 2 cm2 leaf chamber. Conditions in the LI-6800 chamber for light-response experiments were as follows: flow rate, 500 μmol/s; fan speed, 10,000 rpm; leaf temperature, 28oC; RH, 40%; [CO2], 400 μmol/mol; photosynthetic active radiation (PAR), 1,000–100–1,000–0 μmol PAR m-2 s-1 (20 min per light step) (Figure 1a). Light-response measurements of A and g sw were obtained for 120 wild-type Bd21-3 individuals between May 2019 and September 2021. Gas exchange measurements were automatically logged every minute. Relative g sw was calculated by normalising g sw to the highest g sw value observed to evaluate kinetics of stomatal response regardless of variation on absolute g sw, eliminating the influence of stomatal density and leaf area. Because B. distachyon leaves do not fill the 2 cm2 chamber, individual leaf area was measured for a subset of 35 individuals to accurately quantify absolute g sw and A. To obtain a mean approximation of gas exchange levels for the total 120 individuals, absolute g sw and A were corrected by using the average leaf area (0.64 cm2) from the data subset (n = 35). Intrinsic WUE (iWUE) was calculated as the A to g sw ratio (A/g sw). Ambient light intensity (Q out) was monitored during the measurements using an external LI-190R PAR Sensor (Li-COR Biosciences Inc., Lincoln, NE, USA) attached to LI-6800. Greenhouse temperature and RH were monitored during the experiments using an Onset HOBO U12-O12 4-channel data logger (Onset Computer Corporation, Bourne, MA, USA) that was placed next to the plants used for analysis. One-phase decay or one phase association non-linear regressions were obtained for the stomatal closure transitions (1,000–100 and 1,000–0 PAR) and stomatal opening transition (100–1,000 PAR), respectively, to determine half-time (T 50%) and rate constant (k). Maximum stomatal conductance (g smax) measurements were performed with the following conditions: flow rate, 500 μmol/s; fan speed, 10,000 rpm; leaf temperature, 28oC; RH, 68–70%; [CO2], 100 μmol/mol; PAR, 1,500 μmol PAR m-2 s-1. Gas exchange measurements were automatically logged every minute and g smax was calculated as the average of the last 5 min at steady-state.
2.3. Microscopy analysis of stomatal anatomical traits
Leaves of a subset of individuals (n = 4–6 per season; n = 5 in summer, n = 6 in autumn, n = 4 in winter) were collected after LI-6800 measurements and fixed in 7:1 ethanol:acetic acid. To prepare samples for imaging, leaf tissue was rinsed in water and mounted on slides in Hoyer’s solution. The abaxial side was imaged using a Leica DM5000B (Leica Microsystems, Wetzlar, Germany). For SD, 3–5 fields of view (0.290 mm2, 20× objective) per leaf were counted resulting in 60–160 stomata per individual. For SL and width (W A), stomata from 4 to 6 fields of view (0.0725 mm2, 40× objective) per leaf were measured resulting in 20–70 stomata per leaf.
2.4. Correlation analysis and statistics
Correlation analysis for gas exchange parameters independent of leaf area like iWUE (ratio between A and g sw) and half-time of opening/closing were performed for all 120 individuals. Correlation analysis for gas exchange parameters dependent on leaf area like g sw and A were only performed for the subset of 35 individuals, for which the accurate individual leaf area was determined. For the correlation analysis between steady-state gas exchange and environmental conditions, the last 5 min of steady-state gas exchange parameters (A, g sw and iWUE) at the second, third and fourth light step (100, 1,000 and 0 PAR) and their corresponding ambient conditions (Q out, RH and T) were averaged. Pearson correlation matrices represented as heatmaps were obtained for steady-state g sw, A, iWUE, Q out, RH, T and time at 100, 1,000 and 0 PAR. ROUT method was used to remove outliers (Motulsky & Brown, Reference Motulsky and Brown2006).
For the correlation analysis between stomatal kinetics and steady-state g sw, Pearson correlation matrices represented as heatmaps were obtained for steady-state initial/final g sw and half-time (n = 35). For the correlation between stomatal kinetics and environmental conditions, Pearson correlation matrices represented as heatmaps were obtained for half-time, initial/final T, RH Q out, and time of the day (n = 120).
Relevant correlations between different pairs of parameters were represented with linear or non-linear regressions. Significant (p < .05) and non-significant linear regressions were represented with solid and dashed lines, respectively. Non-linear regressions (quadratic function) were chosen for correlations with time of the day as we observe an axis of symmetry around noon, for which a non-linear model was biologically more appropriate.
For the correlation analysis between stomatal anatomical traits and growth environmental conditions, Pearson correlation matrices represented as heatmaps were obtained for SD, SL and environmental growth conditions (average T, average RH and day length). Pearson correlation matrices and linear regressions were obtained for the correlation analysis between stomatal anatomy and gas exchange parameters (steady-state g sw, A, iWUE and T 50%).
To test for significant differences between two groups we performed an unpaired t-test. One-way ANOVA followed by Tukey’s multiple comparison test was used when comparing more than two groups. p values are indicated directly in the graphs. All analyses were performed on GraphPad Prism version 9.1.0, GraphPad Software, San Diego, CA, USA, www.graphpad.com.
2.5. Stomatal morphometric analysis and anatomical gsmax calculations
To characterise fully open stomata, leaves were treated with 4 μM fusicoccin (Santa Cruz Biotechnology, Inc., Dallas, TX, USA) solution in opening-closing buffer (50 mM KCl and 10 mM MES-KOH). Collected leaves were dipped into 70% ethanol and infiltrated with fusicoccin solution. For infiltration, a needleless syringe was used to infiltrate the leaf tissue on the adaxial side until the tissue was visibly wet. Infiltrated leaves were then cut into smaller pieces (approx. 3–5 mm long) and incubated overnight in fusicoccin solution in the light. To analyse closed stomata, leaves were treated with 50 μM ABA (Merck, Darmstadt, Germany) solution in opening-closing buffer (50 mM KCl and 10 mM MES-KOH) as described for fusicoccin treatment and incubated overnight in the dark. Before imaging on the confocal microscope, leaves were stained with propidium iodide (1%) (Thermo Fisher Scientific, Waltham, MA, USA). Z-stacks of 30 open and 30 closed stomata from 3 different individuals each were taken on the Leica TCS SP8 confocal microscope (Leica Microsystems, Wetzlar, Germany). The obtained z-stacks were analysed using Fiji. For each stoma, pore length (PL), pore width (PW) at the centre of the pore, guard cell length (GCL), right and left guard cell width at the middle of the stoma (GCWC) and stoma width at the apices were measured (W A) on the z-sum projection image. To measure the exact pore area, each pore was manually traced with the polygon selection tool. To measure pore depth (l), the central pore part was selected with the rectangle selection tool and resliced starting at the top, avoiding interpolation. Pore depth was measured on the z-sum projection of the reslice.
The leaves assessed for physiological g smax were fixed with ethanol:acetic acid 7:1. GCL, W A and stomatal width at the centre (W C) were measured on light microscope pictures (20–40 stomata per leaf). GC width was calculated as half of W A or W C. SD was obtained by counting the number of stomata in five different areas per leaf (20× objective).
For the anatomical maximum stomatal conductance calculations, we used the anatomical g smax equation from Franks and Farquhar (Reference Franks and Farquhar2001) (see Figure 5h), where SD is the stomatal density (stomata mm-2,), a max is the maximum pore area (μm2), l is the pore depth (μm), d is the diffusivity of water in air (0.0000249 m2/s−1, at 25°C), v is the molar volume of air (0.024464 m3/mol−1, at 25°C) and π is the mathematical constant. Maximum pore area was either measured by hand-tracing the stomatal pore of fully open stomata or calculated as an ellipse (with major axis equal to pore length and minor to half the pore length) or a rectangle (pore width × pore length).
3. Results
3.1. B. distachyon shows fast and consistent stomatal gas exchange in response to changing light
The plants’ physiology including stomatal gas exchange dynamics are strongly influenced by the environmental conditions the plant is exposed to (Arve et al., Reference Arve, Terfa, Gislerød, Olsen and Torre2013; Durand et al., Reference Durand, Brendel, Buré and Le Thiec2020; Matthews et al., Reference Matthews, Vialet-Chabrand and Lawson2017). Closed-system IRGA allow gas exchange measurements within a chamber with tightly controlled environmental settings regardless of ambient conditions (Douthe et al., Reference Douthe, Gago, Ribas-Carbó, Núñez, Pedrol, Flexas, Sánchez-Moreiras and Reigosa2018). However, plants have previously developed in and acclimated to a specific environment. Furthermore, during measurements, most distal parts of the plant remain exposed to ambient environmental conditions that might significantly differ from the conditions in the IRGA chamber (Figure 1a). To quantify the consistency of stomatal responses and the influence of variable greenhouse conditions on gas exchange, we analysed gas exchange parameters and kinetics of 120 wild-type B. distachyon individuals (Bd21-3) over the course of 2 years in a partially environmentally controlled greenhouse. The ambient conditions in the greenhouse varied remarkably over the course of the 120 IRGA measurements (Supplementary Figure S1A–C) and the measurements covered a broad range of hours of the day (time) from 6 am to 7 pm (Supplementary Figure S1D). We obtained consistent and reproducible stomatal light-responses (R 2 = 0.66–0.76) (Figure 1b,e–g) despite variation in absolute stomatal conductance (g sw) levels (Figure 1c). Because the B. distachyon leaf is smaller than the 2 cm2 chamber used, g sw was corrected using the average leaf area from a data subset for which we measured and corrected for the actual individual leaf area (n = 35) (Supplementary Figure S1E), to obtain a mean approximation of gas exchange levels for the total 120 individuals. The data subset corrected with the actual leaf area (Figure 1c, magenta dots) nicely overlapped with the average correction of the 120 individuals (Figure 1c, black dots) and together revealed reasonable variation of absolute g sw values. Importantly, the 35 leaf area corrected measurements covered the range of environmental conditions observed in all 120 individuals (Supplementary Figure S1A–D).
The first light transition (1,000–100 PAR) resulted in a 70% decrease in g sw with a half-time of 6 min (T −50% I = 6.021 min) (Figure 1e). An increase in light intensity (100–1,000 PAR) induced an exponential increase in g sw with a half-time of less than 5 min (T +50% II = 4.493 min) until reaching similar g sw as in the previous high light step (1,000 PAR) (Figure 1f). Switching from 1,000 to 0 PAR resulted in strikingly fast stomatal closure with a half-time of only ~3 min (T −50% III = 3.361 min) and, thus, represented the quickest of the three light transition responses (Figure 1g). g sw was on average 0.29 ± 0.06 mol m−2 s−1 at high light and 0.10 ± 0.05 mol m−2 s−1 at low light (Figure 1c). We observed an average of 0.015 ± 0.013 mol m−2 s−1 of residual g sw in darkness (Figure 1c). At high light, iWUE was on average 77 ± 13 μmol CO2/mol H2O, whereas at low light iWUE was 51 ± 26 μmol CO2/mol H2O (Figure 1d). A was on average 21 ± 4 μmol m−2 s−1 at high light and 4 ± 1 μmol m−2 s−1 at low light (Supplementary Figure S1F).
Together, B. distachyon shows fast stomatal light responses typical for grasses, which were consistent over 120 measurements.
3.2. Quantitative effects of greenhouse environmental fluctuations on gas exchange in B. distachyon
To quantify how the different environmental conditions affected gas exchange, we performed correlation analysis between gas exchange parameters [stomatal conductance (g sw), carbon assimilation (A) and intrinsic water-use efficiency (iWUE)] and environmental conditions [temperature (T), ambient light intensity (Q out), RH, time of the day (time)]. Correlations were done separately for low light steady-state (100 PAR, Figure 2a,g), high light steady-state (1,000 PAR, Figure 2d,h), steady-state in darkness (0 PAR, Supplementary Figure S2E) and for opening and closing kinetics (Supplementary Figure S2H–K). Because exact leaf area was only measured for a subset of 35 individuals to calculate accurate absolute g sw and A, the correlation analysis between environmental parameters and absolute gas exchange parameters (i.e., g sw and A) was performed using the 35 individuals only. On the other hand, the 120 samples were used for correlation analysis between environmental conditions and leaf-area independent parameters like iWUE and stomatal kinetics (half-time).
iWUE was negatively correlated with T at both light conditions 100 and 1,000 PAR (Figure 2a,d). Increasing temperatures significantly associated with decreasing iWUE values (Figure 2b,e). iWUE also correlated with time in both light conditions (Figure 2a,d). A quadratic relation can be observed between time and iWUE (Figure 2c,f), particularly at low light, with the lowest iWUE reached at midday (Figure 2c). Similar correlations between iWUE and T or time were observed in the data subset (n = 35) (Figure 2g,h).
Regarding the steady-state gas exchange parameters, T showed a considerable influence on g sw at both low and high light (Figure 2g,h) and g sw significantly increased with rising ambient temperatures (Figure 2i,j). Q out, on the other hand, correlated with both A and g sw at high light (Figure 2h). Both A and g sw significantly increased with increasing Q out (Figure 2k,l). Together, this explained why iWUE is only correlated with T but not with Q out.
Lastly, time significantly correlated with g sw at all light conditions and also with A at low light (Figure 2g,h and Supplementary Figure S2A–E,G). No significant correlations were observed between ambient conditions (Q out, T or RH) and g sw at 0 PAR (Supplementary Figure S2E–G), even though an influence of T on g sw is suggested (Supplementary Figure S2F) as observed at 100 and 1,000 PAR (Figure 2i,j).
Finally, stomatal kinetics (i.e., half-time T 50%) significantly depended on the initial and/or final steady-state g sw (Supplementary Figure S2H), which in turn were affected by the environment (see above). In addition, our data also suggested an influence of diurnal rhythm (time) on stomatal closure speed (half-time T 50%) as stomata seem to close slower at noon (Figure 2m–o and Supplementary Figure S2I–K).
In summary, fluctuations in ambient conditions such as temperature and light intensity during measurements, and diurnal rhythm influenced steady-state gas exchange parameters and/or stomatal kinetics within a strictly controlled IRGA leaf chamber, which highlights the relevance of considering systemic effects on stomatal physiology experiments.
3.3. Seasonal changes in greenhouse growth conditions affect stomatal anatomical traits and gas exchange
While the artificial light intensity and light-darkness cycles are controlled, the contribution of ambient sunlight intensity and day length (DL), average temperature (T) and RH vary among seasons in our greenhouse (Supplementary Figure S3A,D–F). Stomatal anatomical traits such as SD and stomatal size are strongly influenced by environmental cues to which the plants are exposed to during development (Casson & Gray, Reference Casson and Gray2008; Liu et al., Reference Liu, He, Zhang, Li, Wang, Sack and Yu2018; Qi & Torii, Reference Qi and Torii2018; Terfa et al., Reference Terfa, Olsen and Torre2020). To test the plasticity and variability of stomatal anatomical traits in B. distachyon wild-type plants, we quantified SD and SL as a proxy for stomatal size from 15 individuals grown in different seasons – summer (May to June, n = 5), autumn (October to November, n = 6) and winter (January to February, n = 4) – and correlated these traits with growth conditions (T, RH and DL; Figure 3a and Supplementary Figure S3A). In summer, SD was ~40% higher than in winter (105.8 ± 13.1 vs. 74.2 ± 6.8 stomata per mm2) and SL reduced by ~10% (24.6 ± 0.7 vs. 26.8 ± 0.3 μm) (Supplementary Figure S3A–C). Consequently and as previously described (Franks & Beerling, Reference Franks and Beerling2009; Haworth et al., Reference Haworth, Marino, Loreto and Centritto2021; Zhang et al., Reference Zhang, Wang, Yang, Cui and Niu2021), we observed a strong negative correlation between SD and SL (Figure 3a,d) and a strong correlation between anatomy and environment (Figure 3a). SD and SL correlated in an opposite manner with the different environmental parameters (Figure 3a). An increase in T correlated with an increase in SD and a decrease in SL (Figure 3b,c). Since T and RH were negatively correlated, RH correlated in an opposite manner with SD and SL (Figure 3e,f). SD and SL were also inversely correlated with DL, with longer days associated with shorter stomata and higher SD (Supplementary Figure S3G,H). Overall, summer plants grown during longer days with higher ambient light intensity, higher T, and lower RH (Supplementary Figure S3D–F), developed higher SD and lower SL (Supplementary Figure S3B,C). In autumn and winter, plants grown during shorter days with lower ambient light intensity, lower T and higher RH (Supplementary Figure S3D–F), developed lower SD and higher SL (Supplementary Figure S3B,C). However, which environmental parameter primarily caused changes to stomatal anatomy is unclear.
Besides a seasonal variation on stomatal anatomical traits, we observed seasonal variation on gas exchange (Figure 3g–i). Our data suggested slightly, yet not significantly higher A and lower g sw in summer under high light compared to autumn and winter (Figure 3g,h). Yet, iWUE was significantly higher in summer than in autumn and winter under high light conditions (Figure 3i). Our data also suggests an increase in A and g sw from autumn to winter (Figure 3g,h). Under low light conditions (100 PAR), no differences in A occurred among seasons (Figure 3g). On the other hand, higher g sw was observed in autumn and winter (Figure 3h), resulting in lower iWUE at 100 PAR (Figure 3i). When measuring physiological maximum stomatal conductance (g smax) of autumn/winter plants, the anatomical offset between SD and SL seemed to compensate for stomatal gas exchange maximum capacity, even though a non-significant decrease in average g smax was observed in autumn/winter (Supplementary Figure S3I). In conclusion, stomatal anatomical traits of wild-type B. distachyon are surprisingly plastic and variable among seasons even in a semi-controlled growth environment, likely contributing to the seasonal variation on functional traits.
3.4. Stomatal anatomical traits influence gas exchange
Due to the seasonal variation in stomatal anatomy and functional traits, we quantified how the anatomical variation translates into changes in functional traits such as steady-state gas exchange parameters (at high light 1,000 PAR, low light 100 PAR and darkness 0 PAR) and stomatal kinetics (Supplementary Figure S4).
In terms of steady-state g sw, lower SD (and higher SL, to a lesser extent) are the anatomical traits associated with higher operational stomatal conductance (g sw) in B. distachyon (Figure 4a,b). Similarly, an increase in SD and a decrease in SL resulted in an increase of A at high light, while no effect was observed in 100 PAR (light limiting condition) (Figure 4c,d). Consequently, higher SD and lower SL result in higher iWUE (Figure 4e,f), even though SD had a stronger effect on iWUE (Figure 4f). Thus, the higher iWUE observed in summer (Figure 3i) might be primarily caused by the higher SD and lower SL observed in this season (Supplementary Figure S3B,C). The correlations between anatomical traits (SD and SL) and g sw were stronger in low light than in high light (Supplementary Figure S4 and Figure 4a,b) likely contributing to the higher seasonal variation in iWUE at low light than in high light (Figure 3i).
Regarding stomatal kinetics (T 50%), the effect of SL on stomatal closure and opening was non-significant (Figure 4g and Supplementary Figure S4). However, while the influence of SD on stomatal opening was also weak (Figure 4h and Supplementary Figure S4), surprisingly stronger correlations and significant effects were observed between SD and stomatal closure kinetics (T −50%) (Figure 4h and Supplementary Figure S4). High SD was strongly correlated with water-use efficiency as the increase of SD led to higher steady-state iWUE (Figure 4f) and faster stomatal closure (Figure 4h) contributing to higher water-use efficiency in changing environments. Overall, stomatal anatomy is strongly correlated with stomatal functioning and the seasonal variation in stomatal anatomy strongly contributed to seasonal variation in gas exchange.
3.5. Morphometric analysis of graminoid stomata to optimise anatomical gsmax predictions in B. distachyon
Finally, we wanted to mathematically describe the impact of the observed trade-off between SD and SL on maximum stomatal conductance (g smax) by calculating the theoretical anatomical g smax, which is based on the anatomical traits SD, maximum pore area (a max) and pore depth (l). While SD is assessed for any species simply by counting stomata per leaf area, formulae to calculate maximum pore area (a max, μm2) and pore depth (l, μm) were optimised for Arabidopsis-like stomatal morphologies and ellipsoid pores (Dow et al., Reference Dow, Berry and Bergmann2014; Franks & Farquhar, Reference Franks and Farquhar2001) (Figure 5a).
High-resolution confocal stacks of fusicoccin (fus)-treated open stomatal complexes in B. distachyon revealed hexagonal rather than elliptical pores (Figure 5b). We performed careful morphometric analysis of open and closed stomata to characterise GCL, PL, PW, GCWC and stomatal WA (Figure 5b). Furthermore, we manually traced and measured pore areas of fus-treated complexes to empirically determine a more appropriate way to calculate a max of graminoid stomata. Calculating a max as for Arabidopsis-like stomata (ellipse with major axis equal to pore length and minor to half the pore length, Figure 5a) caused significant pore area overestimation compared to the manually traced pores of fus-treated open complexes (Figure 5b,c). A rectangular rather than an ellipsoid approach still overestimated pore area (Figure 5b and Supplementary Figure S5A), as the manually measured pore areas approximated 90 ± 7% of the rectangular pore area (Supplementary Figure S5C). Thus, B. distachyon stomatal pore area can be accurately estimated using rectangular pore area calculations multiplied by a correction factor of 0.9 (Figure 5d).
Regarding pore depth (l), we observed that the mean pore depth measured from orthogonal resliced confocal stacks (3.37 ± 0.4 μm) approximates the mean GC width at the centre (3.23 ± 0.3 μm) but not the GC W A (5.92 ± 0.4 μm) in fully opened stomata (fus-treated) (Supplementary Figure S5G). In closed stomata (ABA-treated), on the other hand, the mean GC width at the centre was 2.48 ± 0.3 μm while the GC W A was 3.52 ± 0.6 μm (Supplementary Figure S5G). Therefore, if not exactly measured from orthogonal sections, then pore depth could be approximated as central GC width of fully open stomata or apical GC width of closed stomata.
We then calculated anatomical g smax using hand-traced a max and estimated a max using formulae for (a) ellipse, (b) rectangle pore and (c) rectangle pore multiplied by the correction factor 0.9. Pore depth was measured from orthogonal resliced confocal stacks (Figure 5b) and stomatal density was determined by counting stomata in 3–5 different fields of view using light microscopy. We could observe that anatomical gs max using hand-traced a max nicely correlated with anatomical g smax calculated for the rectangle pore multiplied by the correction factor 0.9 (Figure 5f). This was not the case when using ellipse a max (Figure 5e) or rectangular pore without the correction factor (Supplementary Figure S5B,C).
To determine a max from simple light microscopy pictures, where pores are hard to see, we calculated correction coefficients to estimate pore length and width from GC length. By calculating the ratios of the morphometrically determined GCL, PL and PW, we found that PL is 44 ± 3% of GCL and PW is 13 ± 3% of GCL (Supplementary Figure S5D–F). Thus, for calculations using light microscopy pictures, we estimated PL as 0.44*GCL and PW as 0.13*GCL (Figure 5g,h). To approximate pore depth (l), GC W A was used for closed or partially open stomata (½ of the stoma W A) (Figure 5g,h and Supplementary Figure S5G) and GC width at the centre was used for fully open stomata.
Next, we tested if our adjustments for anatomical g smax calculations could be used to reliably predict physiological g smax in B. distachyon. We performed IRGA-based measurements of g smax (physiological g smax) in four independent individuals, collected these exact leaf zones, and measured anatomical traits from segments of those by using both standard light microscopy (after fixation with ethanol:acetic acid 7:1) and confocal microscopy (after treatment with fusicoccin). No significant differences were found between physiological g smax and anatomical g smax based on measured anatomical parameters using light microscopy (LM) or confocal microscopy (CM) (Figure 5i and Supplementary Figure S5H). In summary, the optimised formula for accurate anatomical g smax estimation can be used to reliably predict physiological g smax in B. distachyon.
Finally, we calculated anatomical g smax for the summer, autumn and winter individuals whose anatomical traits (SD and SL) were shown in Supplementary Figure S3B,C. Even though a decrease in average anatomical g smax was observed in autumn and winter, this difference was non-significant (Figure 5j). These results match our observations in physiological g smax measurements between summer and autumn/winter plants (Supplementary Figure S3I). Ultimately, these findings suggest that the trade-off between SD and SL in wild-type B. distachyon might serve as a mechanism to maximise stomatal conductance capacity in different environments.
4. Discussion
Consistent and reproducible stomatal kinetics were observed for B. distachyon regardless of the variable greenhouse conditions. B. distachyon displayed the fast stomatal movements typical for grass species, which are faster than most non-grass species with kidney-shaped GC (Franks & Farquhar, Reference Franks and Farquhar2007; Grantz & Assmann, Reference Grantz and Assmann1991; McAusland et al., Reference McAusland, Vialet‐Chabrand, Davey, Baker, Brendel and Lawson2016; Merilo et al., Reference Merilo, Jõesaar, Brosché and Kollist2014). The quick stomatal movements of grasses like B. distachyon are associated with the graminoid morphology, where two lateral subsidiary cells (SCs) flank the central, dumbbell-shaped GCs (Gray et al., Reference Gray, Liu and Facette2020; Nunes et al., Reference Nunes, Zhang and Raissig2020; Stebbins & Shah, Reference Stebbins and Shah1960). Fast stomatal movements in grasses require SCs (Raissig et al., Reference Raissig, Matos, Anleu Gil, Kornfeld, Bettadapur, Abrash, Allison, Badgley, Vogel, Berry and Bergmann2017), which might function as specialised ion reservoirs (Raschke & Fellows, Reference Raschke and Fellows1971) and mechanically accommodate GC movement to accelerate both stomatal opening and closing (Franks & Farquhar, Reference Franks and Farquhar2007). In addition, the reduced volume-to-surface ratio of dumbbell-shaped GCs likely requires less exchange of water and ions to be pressurised (Franks & Farquhar, Reference Franks and Farquhar2007).
Furthermore, no major asymmetry between closure and opening speed was observed. This is consistent with the results from a comparison of stomatal kinetics between eight species with kidney-shaped GCs and seven species with dumbbell-shaped GCs, where species with dumbbell-shaped GCs displayed the quickest responses and showed the most similarity between opening and closure times (McAusland et al., Reference McAusland, Vialet‐Chabrand, Davey, Baker, Brendel and Lawson2016). Nonetheless, the fastest stomatal movement in B. distachyon was dark-induced stomatal closure (1,000–0 PAR). Faster stomatal closing than opening has been previously described for several species and suggested to be a water-conserving strategy (Lawson & Vialet-Chabrand, Reference Lawson and Vialet-Chabrand2019; Leakey et al., Reference Leakey, Ferguson, Pignon, Wu, Jin, Hammer and Lobell2019; McAusland et al., Reference McAusland, Vialet‐Chabrand, Davey, Baker, Brendel and Lawson2016). Fast stomatal movements are important to quickly adjust stomatal pores to avoid excess of water vapor loss through stomata (g sw) during suboptimal carbon assimilation (A) in low light. Intrinsic water use efficiency (iWUE, the ratio between A and g sw), varied between 50 and 80 μmol/mol which is consistent with the range of iWUE described for other C3 grass species such as wheat (25–65 μmol/mol) and rice (50–80 μmol/mol) (Giuliani et al., Reference Giuliani, Koteyeva, Voznesenskaya, Evans, Cousins and Edwards2013; Jahan et al., Reference Jahan, Amthor, Farquhar, Trethowan and Barbour2014).
Despite the consistent stomatal responsiveness in B. distachyon, we observed a significant influence of the time of the day on light-response stomatal kinetics. While diurnal variation on gas exchange has been well described for C3 species (de Dios, Reference Resco de Dios2017; Matthews et al., Reference Matthews, Vialet-Chabrand and Lawson2017; Miao et al., Reference Miao, Cai, Wu and Wang2021; Roussel et al., Reference Roussel, Dreyer, Montpied, Le-Provost, Guehl and Brendel2009; Stangl et al., Reference Stangl, Tarvainen, Wallin, Ubierna, Räntfors and Marshall2019; Vahisalu et al., Reference Vahisalu, Kollist, Wang, Nishimura, Chan, Valerio, Lamminmäki, Brosché, Moldau, Desikan, Schroeder and Kangasjärvi2008), the observed diurnal variation on stomatal responsiveness to changing light in the tightly regulated conditions of the IRGA chambers was compelling. Stomatal closure and opening speed were mainly affected by the time of the day and by steady-state g sw prior and/or after the change in light intensity. A similar dependence of g sw kinetics on light intensity transitions, on the time of the day and on steady-state g sw prior to light intensity changes has been recently described in Musa spp. (Eyland et al., Reference Eyland, van Wesemael, Lawson and Carpentier2021).
Steady-state g sw, on the other hand, was significantly influenced by the environmental conditions. Stomatal conductance was stimulated by increasing ambient temperatures. This effect has been observed and described as a leaf cooling mechanism to cope with higher temperatures (Gommers, Reference Gommers2020; Lamba et al., Reference Lamba, Hall, Räntfors, Chaudhary, Linder, Way, Uddling and Wallin2018; Sonawane et al., Reference Sonawane, Koteyeva, Johnson and Cousins2021; Urban et al., Reference Urban, Ingwers, McGuire and Teskey2017a; Reference Urban, Ingwers, McGuire and Teskey2017b). Yet, after exceeding a certain threshold, such high temperatures may lead to stress-induced stomatal closure (Faria et al., Reference Faria, García-Plazaola, Abadía, Cerasoli, Pereira and Chaves1996; Ikkonen et al., Reference Ikkonen, Shibaeva, Rosenqvist and Ottosen2015; Yamori et al., Reference Yamori, Suzuki, Noguchi, Nakai and Terashima2006; Zhou et al., Reference Zhou, Xu, Pan and Yu2015). In addition, increasing ambient light intensity also significantly impacted IRGA measurements by triggering increases in g sw and A levels. Recent studies have also reported systemic stomatal responses to light and heat in Arabidopsis (Devireddy et al., Reference Devireddy, Zandalinas, Gómez-Cadenas, Blumwald and Mittler2018; Reference Devireddy, Arbogast and Mittler2020) and in response to darkness and elevated CO2 in birch and poplar (Ehonen et al., Reference Ehonen, Holtta and Kangasjärvi2020). Consistently, our results suggest that B. distachyon stomata integrate both local environmental cues and systemic signals from distal parts of the plant. Therefore, it is important to monitor environmental conditions and consider their impact on gas exchange measurements, particularly in greenhouse or field studies with significant environmental fluctuations.
Apart from the effect of environmental conditions on gas exchange, different growth conditions had a major impact on anatomical traits. SD and SL inversely varied among seasons due to significant variation in ambient growth conditions. The seasonal trade-off between SD and SL mostly maintained maximum gas exchange capacity in the different growth environments. In addition, higher SD was associated with faster stomatal closure and the combination of higher SD and lower SL associated with improved iWUE in wild-type B. distachyon. This suggests that higher SD and lower SL, which are associated with improved stomatal responsiveness and more water use efficient gas exchange, are a morphological adaptation to summer. Higher stomatal densities in warmer environments were reported and associated as an ecophysiological significant response for leaf evaporative cooling (Carlson et al., Reference Carlson, Adams and Holsinger2016; Hill et al., Reference Hill, Guerin, Hill and Watling2015). In contrast, autumn and winter seasons feature shorter days, decreased light intensity and colder temperatures, which can negatively affect photosynthesis (Feng et al., Reference Feng, Raza, Li, Chen, Khalid, Du, Liu, Wu, Song, Yu, Zhang, Yuan, Yang and Yang2018; Yamasaki et al., Reference Yamasaki, Yamakawa, Yamane, Koike, Satoh and Katoh2002). Thus, a decrease in stomatal density to increase the leaf surface allocated to light harvesting, compensated by an increase in stomatal size to maintain maximum gas exchange capacity, might sustain optimal gas exchange and photosynthesis in winter. While we observed higher iWUE in B. distachyon wild-type plants with lower SL and higher SD, crop species (wheat, barley and rice) overexpressing an inhibitor of stomatal development (EPF1) show a reduction in SD and improved iWUE (Caine et al., Reference Caine, Yin, Sloan, Harrison, Mohammed, Fulton, Biswal, Dionora, Chater, Coe, Bandyopadhyay, Murchie, Swarup, Quick and Gray2018; Dunn et al., Reference Dunn, Hunt, Afsharinafar, Al Meselmani, Mitchell, Howells, Wallington, Fleming and Gray2019; Hughes et al., Reference Hughes, Hepworth, Dutton, Dunn, Hunt, Stephens, Waugh, Cameron and Gray2017). Thus, genetically modifying SD (and/or SL) in B. distachyon beyond the intraspecific range of variation might allow to improve iWUE. B. distachyon genotypes varying in single morphological traits may help to better understand the independent influence of SL and SD on gas exchange kinetics, capacity and water-use efficiency, and verify the correlations observed in this present study.
Some studies described a negative correlation between stomatal size and stomatal speed (e.g., (Drake et al., Reference Drake, Froend and Franks2013; Durand et al., Reference Durand, Brendel, Buré and Le Thiec2020; Kardiman & Ræbild, Reference Kardiman and Ræbild2018). In a study comparing different rice genotypes, stomatal size was negatively correlated with stomatal half-time (Zhang et al., Reference Zhang, Peng and Li2019), with larger stomata being faster. In a comparison of 15 different species (including 7 grass species), however, smaller stomata were associated with faster stomatal movements (McAusland et al., Reference McAusland, Vialet‐Chabrand, Davey, Baker, Brendel and Lawson2016). Other approaches using a broader range of species from different plant groups (by comparing 7 (Elliott-Kingston et al., Reference Elliott-Kingston, Haworth, Yearsley, Batke, Lawson and McElwain2016) and 31 (Haworth et al., Reference Haworth, Scutt, Douthe, Marino, Gomes, Loreto, Flexas and Centritto2018) different species) suggest that stomatal speed is not related to SL but rather positively correlated with SD, as we observed in B. distachyon. A major effect of SD and minor effect of SL on stomatal speed under fluctuating light has also been described in a study comparing different Arabidopsis genotypes (Sakoda et al., Reference Sakoda, Yamori, Shimada, Sugano, Hara-Nishimura and Tanaka2020). Higher SD could potentially cause proximity effects triggering stomata in close vicinity to react to local changes in a more coordinated manner. However, the effect of the variation of SL and SD on gas exchange speediness may vary among species and/or genotypes. Nonetheless, environmentally induced SD and SL variation and its impact on gas exchange must be considered during long term studies performed in greenhouse or field settings.
To facilitate correlations of anatomical stomatal traits to theoretical gas exchange maximum capacity (anatomical g smax) in B. distachyon, we adjusted established equations (Dow et al., Reference Dow, Berry and Bergmann2014; Franks & Farquhar, Reference Franks and Farquhar2001) to the graminoid stomatal morphology. For grass stomata, stomatal pores are rather hexagons than ellipses and the equation presented in this study accurately predicted anatomical g smax for B. distachyon. Therefore, it can be reliably used to predict gas exchange capacity of B. distachyon genotypes varying in anatomical traits. Differences between the anatomical and physiological g smax might reveal impaired stomatal signalling and thus, provide a tool to identify mutant phenotypes of stomatal function. In addition, anatomical g smax could help to weigh the effect of variations in single morphological traits (e.g., stomatal density, stomatal size, pore area) on g smax.
In conclusion, stomatal conductance kinetics are fast and consistent in the grass model species B. distachyon. Nevertheless and even though stomata primarily respond to the local environment (i.e., within the IRGA chamber), ambient light intensity, temperature and time of the day can have systemic effects on gas exchange influencing results from IRGA measurements. Stomatal anatomical traits are highly plastic and environmental-responsive and, furthermore, have a major impact on gas exchange. For that reason, the effect of growth conditions on stomatal anatomical traits must be considered in leaf-level gas exchange studies.
Acknowledgements
The authors would like to thank Michael Schilbach for managing the greenhouse and gardening support. Furthermore, the authors acknowledge all members of the lab for carefully reading and critically evaluating this manuscript.
Financial support
This work was supported by a DAAD Study Scholarship (to M.W.S.) and by the German Research Foundation (D.F.G.) Emmy Noether Programme Grant RA 3117/1-1 (to M.T.R.).
Conflicts of interest
The authors declare no potential conflicts of interest.
Authorship contributions
T.D.G.N. and M.T.R. designed the study and conceived the experiments. T.D.G.N. and M.W.S. performed the experiments. T.D.G.N., M.W.S., H.L. and M.T.R. analysed and interpreted the data. T.D.G.N., M.W.S., H.L. and M.T.R. wrote the manuscript.
Data availability statement
All gas exchange data and anatomical data used in this study can be found in Supplementary Dataset 1. Stomatal images are available upon request.
Supplementary Materials
To view supplementary material for this article, please visit http://doi.org/10.1017/qpb.2021.19.
Comments
Dear Olivier
With our great pleasure we are submitting our manuscript „Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model“. We are convinced that it fits very well into the scope of this journal. We are taking a highly quantitative approach to measure the effect of fluctuating environmental conditions on stomatal gas exchange and stomatal anatomical traits in the model grass Brachypodium distachyon.
Over the course of almost two years we performed leaf-level gas exchange measurements on 100 wild-type B. distachyon Bd21-3 individuals. While initially struggling with the fact that our greenhouse is not extremely well controlled and thus affecting our measurements, we decided to actually quantify the effects of fluctuating environmental conditions on our measurements, which occur within a leaf chamber head, where the environmental parameters are perfectly controlled. To evaluate both steady-state gas exchange parameters and stomatal opening and closing speed we employed a changing light regime. While the gas exchange responses were surprisingly robust, we could clearly identify an influence of the ambient environment on local gas exchange.
Furthermore, we quantified for 15 of the 100 plants how seasonal growth conditions affect stomatal anatomy. We found that both the number and size of stomata are extremely plastic within wild-type B. distachyon. The quantitative description of the anatomical parameters, ultimately allowed us to correlate how anatomical parameters affect gas exchange physiology. One of the surprises was that a higher stomatal density positively affects stomatal closure speediness.
Finally, we adapted a formula that allows the prediction of gas exchange capacity (gsmax) from anatomical parameters to the grass stomatal morphology. Originally, this equation was formulated to account for kidney-shaped, Arabidopsis-like stomata. We used detailed confocal stacks of open and closed stomata to accurately measure the stomatal pore and then empirically optimized the formula to calculate maximum pore area and pore depth from simple light microscopy images. These adjustments allowed us to accurately predict physiological gsmax from anatomical parameters.
To our knowledge, the dataset presented here is the first to systematically quantify intraspecific variation of stomatal gas exchange and anatomy in a semi-controlled greenhouse setting. The careful correlation analysis will aid us and others in the field to predict, which environmental parameters to consider that could affect leaf-level gas exchange measurements. Furthermore, the empirically optimized methods to calculate maximal stomatal pore area and pore depth of grass stomata will enable a more accurate prediction of gas exchange capacity in grasses.
We would like to thank you very much for considering our work and we are looking forward to hearing from you.
With my very best wishes,
Michael