Hostname: page-component-78c5997874-xbtfd Total loading time: 0 Render date: 2024-11-02T23:46:37.472Z Has data issue: false hasContentIssue false

Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model

Published online by Cambridge University Press:  09 March 2022

Tiago D. G. Nunes
Affiliation:
Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
Magdalena W. Slawinska
Affiliation:
Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
Heike Lindner
Affiliation:
Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
Michael T. Raissig*
Affiliation:
Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany Institute of Plant Sciences, University of Bern, Bern, Switzerland
*
Address correspondence to: Michael T. Raissig E-mail: [email protected]

Abstract

Stomata are cellular pores on the leaf epidermis that allow plants to regulate carbon assimilation and water loss. Stomata integrate environmental signals to regulate pore apertures and adapt gas exchange to fluctuating conditions. Here, we quantified intraspecific plasticity of stomatal gas exchange and anatomy in response to seasonal variation in Brachypodium distachyon. Over the course of 2 years, we (a) used infrared gas analysis to assess light response kinetics of 120 Bd21-3 wild-type individuals in an environmentally fluctuating greenhouse and (b) microscopically determined the seasonal variability of stomatal anatomy in a subset of these plants. We observed systemic environmental effects on gas exchange measurements and remarkable intraspecific plasticity of stomatal anatomical traits. To reliably link anatomical variation to gas exchange, we adjusted anatomical g smax calculations for grass stomatal morphology. We propose that systemic effects and variability in stomatal anatomy should be accounted for in long-term gas exchange studies.

Type
Original Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press in association with The John Innes Centre

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.

Figure 1. Leaf-level gas exchange measurements in response to changing light intensities reveal fast and consistent stomatal movements in Brachypodium distachyon Bd21-3. (a) Experimental setup for measuring leaf-level gas exchange parameters related with CO2 capture and H2O vapor loss by clamping a leaf in an infrared gas analyser (IRGA) chamber with controlled environmental conditions. Gas exchange is measured in changing light conditions (1,000–100–1,000–0 PAR) inducing stomatal closure in response to decreasing light intensity/darkness and stomatal opening in response to increasing light intensity. (b) Relative stomatal conductance (Rel g sw) during the light transitions (n = 120, normalised to highest g sw observed). (c) Absolute stomatal conductance (g sw) response to light transitions (in black, data from 120 individuals corrected by average leaf area of 0.64 cm2 and in magenta data from a subset of 35 individuals corrected by individual leaf area). (d) Intrinsic water-use efficiency (iWUE) response to light transitions (1,000–100–1,000 PAR) (n = 120, calculated as A/g sw). (e) One-phase decay exponential regression for the transition 1,000–100 PAR (n = 120). (f) One-phase association exponential regression for the transition 100–1,000 PAR (n = 120). (g) One-phase decay exponential regression for the transition 1,000–0 PAR (n = 120). R 2, half-time (T ±50%) and rate constant (K) are indicated. Error bars = SD.

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).

Figure 2. Influence of temperature, ambient light and diurnal rhythm on leaf-level gas exchange. (A) Correlation matrix between iWUE and environment (Q out, T, RH and time) of 120 measurements of wild-type B. distachyon (Bd21-3) at 100 PAR (second light step). (b) Linear regression between T and iWUE at 100 PAR (n = 116). (c) Non-linear regression between time and iWUE at 100 PAR (n = 116). (d) Correlation matrix between iWUE and environment (Q out, T, RH and time) of 120 measurements of wild-type B. distachyon (Bd21-3) at 1,000 PAR (third light step). (e) Linear regression between T and iWUE at 1,000 PAR (n = 119). (f) Non-linear regression between time and iWUE at 1,000 PAR (n = 120). (g) Correlation matrix between gas exchange parameters (A, g sw and iWUE) and environment (Q out, T, RH and time) of the 35 measurements (corrected by individual leaf area) of wild-type B. distachyon (Bd21-3) at 100 PAR (second light step). (h) Correlation matrix between gas exchange parameters (A, g sw and iWUE) and environment (Q out, T, RH and time) of the 35 measurements (corrected by individual leaf area) of wild-type B. distachyon (Bd21-3) at 1,000 PAR (third light step). (i) Linear regression between T and g sw at 100 PAR (n = 35). (j) Linear regression between T and g sw at 1,000 PAR (n = 35). (k) Linear regression between Q out and A at 1,000 PAR (n = 35). (l) Linear regression between Q out and g sw at 1,000 PAR (n = 35). (m) Non-linear regression between half-time of the transition 1,000–100 PAR (T −50%I) and time of the day (time) (n = 99). (n) Non-linear regression between half-time of the transition 100–1,000 PAR (T +50%II) and time of the day (time) (n = 111). (o) Non-linear regression between half-time of the transition 1,000–0 PAR (T −50%III) and time of the day (time) (n = 111). R 2 and Sy.x or p values are indicated.

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.

Figure 3. Effects of seasonal growth conditions on stomatal anatomical traits and on gas exchange. (a) Correlation matrix between stomatal length (SL), stomatal density (SD), average growth temperature (T), average growth relative humidity (RH) and day length (DL) (n = 15). (b) Non-linear quadratic (second order polynomial regression) relation between T and SD (n = 15). (c) Non-linear quadratic (second order polynomial regression) relation between T and SL (n = 15). (d) Linear relation between SD and SL (n = 15). (e) Non-linear quadratic (second order polynomial regression) relation between RH and SD (n = 15). (f) Non-linear quadratic (second order polynomial regression) relation between RH and SL (n = 15). (g) Seasonal variation on A at 1,000 (white) and 100 (grey) PAR (n = 4–6 per season; n = 5 in summer, n = 6 in autumn, n = 4 in winter). (h) Seasonal variation on g sw at 1,000 and 100 PAR (n = 4–6 per season; n = 5 in summer, n = 6 in autumn, n = 4 in winter). (i) Seasonal variation on iWUE at 1,000 and 100 PAR (n = 4–6 per season; n = 5 in summer, n = 6 in autumn, n = 4 in winter). R 2 and Sy.x or p values are indicated. Dashed lines in (B-F) indicate 95% confidence bands.

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).

Figure 4. Impact of stomatal anatomical traits on steady-state gas exchange and stomatal kinetics. (a) Linear regressions between stomatal length (SL) and g sw at 1,000 (white dots), 100 (grey dots) and 0 (black dots) PAR (n = 15). (b) Linear regressions between stomatal density (SD) and g sw at 1,000 (white dots), 100 (grey dots) and 0 (black dots) PAR (n = 15). (c) Linear regressions between SL and A at 1,000 (white dots) and 100 PAR (grey dots). (d) Linear regression between SD and A at 1,000 (white dots) and 100 PAR (grey dots). (e) Linear regressions between SL and iWUE at 1,000 PAR (white dots) and 100 PAR (grey dots) (n = 15). (f) Linear regressions between SD and iWUE at 1,000 PAR (white dots) and 100 PAR (grey dots) (n = 15). (g) Linear regressions between SL and half-time (T 50%) of the light transitions 1,000–100 (grey dots), 100–1,000 (white dots) and 1,000–0 (black dots) PAR (n = 15). (h) Linear regressions between SD and T 50% of the light transitions 1,000–100 (grey dots), 100–1,000 (white dots) and 1,000–0 (black dots) PAR (n = 15). R 2 and p values are indicated. Dashed lines indicate statistically non-significant linear regressions (p > 0.05).

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).

Figure 5. Morphometric analysis of graminoid B. distachyon stomata significantly improves anatomical g smax predictions. (a) Arabidopsis-like stoma and ellipse pore shape. (b) B. distachyon stomatal morphology traits measured; guard cell length (GCL), pore length (PL), pore width (PW), guard-cell width at the centre of the stomata (GCWC), stomatal width at the apex (W A) and pore depth (l). Pore area hand-traced (red) or geometrically defined as an ellipse (orange) or a rectangle (blue). (c) Linear relation between hand-traced pore area and ellipse pore area. (d) Linear relation between hand-traced pore area and rectangle multiplied by 0.9. (e) Linear relation of anatomical g smax calculated with hand-traced pore area and with ellipse pore. (f) Linear relation of anatomical g smax calculated with hand-traced pore area and with rectangle pore multiplied by 0.9. (g) Anatomical parameters measured using light microscopy; stomatal width at the apex (W A) and GCL. (h) Anatomical maximum stomatal conductance (anatomical g smax) equation as defined by Franks and Farquhar (Reference Franks and Farquhar2001) and B. distachyon adjustments to calculate a max (0.9*PL*PW), PL (0.44*GCL), PW (0.13*GCL) and l (W A*0.5). (i) Comparison between physiological g smax, anatomical g smax (light microscopy, LM) and anatomical g smax (confocal microscopy, CM). (j) Comparison of anatomical g smax calculated for summer, autumn and winter plants with stomatal anatomical traits represented in Supplementary Figure S3B,C. R 2 and p values are indicated.

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.

References

Arve, L. E., Terfa, M. T., Gislerød, H. R., Olsen, J. E., & Torre, S. (2013). High relative air humidity and continuous light reduce stomata functionality by affecting the ABA regulation in rose leaves. Plant, Cell Environment, 36, 382392. https://doi.org/10.1111/j.1365-3040.2012.02580.x CrossRefGoogle ScholarPubMed
Bertolino, L. T., Caine, R. S., & Gray, J. E. (2019). Impact of stomatal density and morphology on water-use efficiency in a changing world. Frontiers in Plant Science, 10, 225. https://doi.org/10.3389/fpls.2019.00225 CrossRefGoogle Scholar
Caine, R. S., Yin, X., Sloan, J., Harrison, E. L., Mohammed, U., Fulton, T., Biswal, A. K., Dionora, J., Chater, C. C., Coe, R. A., Bandyopadhyay, A., Murchie, E. H., Swarup, R., Quick, W. P., & Gray, J. E. (2018). Rice with reduced stomatal density conserves water and has improved drought tolerance under future climate conditions. The New Phytologist, 221, 371384. https://doi.org/10.1111/nph.15344 CrossRefGoogle ScholarPubMed
Carlson, J. E., Adams, C. A., & Holsinger, K. E. (2016). Intraspecific variation in stomatal traits, leaf traits and physiology reflects adaptation along aridity gradients in a South African shrub. Annals of Botany, 117, 195207. https://doi.org/10.1093/aob/mcv146 CrossRefGoogle Scholar
Casson, S., & Gray, J. E. (2008). Influence of environmental factors on stomatal development. The New Phytologist, 178, 923. https://doi.org/10.1111/j.1469-8137.2007.02351.x CrossRefGoogle ScholarPubMed
Resco de Dios, V. (2017). Circadian regulation and diurnal variation in gas exchange. Plant Physiology, 175, 34. https://doi.org/10.1104/pp.17.00984 CrossRefGoogle ScholarPubMed
Devireddy, A. R., Arbogast, J., & Mittler, R. (2020). Coordinated and rapid whole-plant systemic stomatal responses. The New Phytologist, 225, 2125. https://doi.org/10.1111/nph.16143 CrossRefGoogle ScholarPubMed
Devireddy, A. R., Zandalinas, S. I., Gómez-Cadenas, A., Blumwald, E., & Mittler, R. (2018). Coordinating the overall stomatal response of plants: Rapid leaf-to-leaf communication during light stress. Science Signaling, 11, eaam9514. https://doi.org/10.1126/scisignal.aam9514 CrossRefGoogle ScholarPubMed
Douthe, C., Gago, J., Ribas-Carbó, M., Núñez, R., Pedrol, N., & Flexas, J. (2018). Measuring photosynthesis and respiration with infrared gas analysers. In Sánchez-Moreiras, A. M. & Reigosa, M. J. (Eds.), Advances in plant ecophysiology techniques (pp. 5175). Springer International Publishing. https://doi.org/10.1007/978-3-319-93233-0_4 CrossRefGoogle Scholar
Dow, G. J., Berry, J. A., & Bergmann, D. C. (2014). The physiological importance of developmental mechanisms that enforce proper stomatal spacing in Arabidopsis thaliana. The New Phytologist, 201, 12051217. https://doi.org/10.1111/nph.12586 CrossRefGoogle ScholarPubMed
Drake, P. L., Froend, R. H., & Franks, P. J. (2013). Smaller, faster stomata: scaling of stomatal size, rate of response, and stomatal conductance. Journal of Experimental Botany, 64, 495505. https://doi.org/10.1093/jxb/ers347 CrossRefGoogle ScholarPubMed
Dunn, J., Hunt, L., Afsharinafar, M., Al Meselmani, M., Mitchell, A., Howells, R., Wallington, E., Fleming, A. J., & Gray, J. E. (2019). Reduced stomatal density in bread wheat leads to increased water-use efficiency. Journal of Experimental Botany, 70, 47374748. https://doi.org/10.1093/jxb/erz248 CrossRefGoogle ScholarPubMed
Durand, M., Brendel, O., Buré, C., & Le Thiec, D. (2020). Changes in irradiance and vapour pressure deficit under drought induce distinct stomatal dynamics between glasshouse and field-grown poplars. The New Phytologist, 227, 392406. https://doi.org/10.1111/nph.16525 CrossRefGoogle ScholarPubMed
Ehonen, S., Holtta, T., & Kangasjärvi, J. (2020). Systemic signaling in the regulation of stomatal conductance. Plant Physiology, 182, 18291832. https://doi.org/10.1104/pp.19.01543 CrossRefGoogle ScholarPubMed
Elliott-Kingston, C., Haworth, M., Yearsley, J. M., Batke, S. P., Lawson, T., & McElwain, J. C. (2016). Does size matter? Atmospheric CO2 may be a stronger driver of stomatal closing rate than stomatal size in taxa that diversified under low CO2. Frontiers in Plant Science, 7, 1253. https://doi.org/10.3389/fpls.2016.01253 CrossRefGoogle ScholarPubMed
Engineer, C. B., Hashimoto-Sugimoto, M., Negi, J., Israelsson-Nordström, M., Azoulay-Shemer, T., Rappel, W.-J., Iba, K., & Schroeder, J. I. (2016). CO2 sensing and CO2 regulation of stomatal conductance: Advances and open questions. Trends in Plant Science, 21, 1630. https://doi.org/10.1016/j.tplants.2015.08.014 CrossRefGoogle ScholarPubMed
Eyland, D., van Wesemael, J., Lawson, T., & Carpentier, S. (2021). The impact of slow stomatal kinetics on photosynthesis and water use efficiency under fluctuating light. Plant Physiology, 186, 9981012. https://doi.org/10.1093/plphys/kiab114 CrossRefGoogle ScholarPubMed
Faralli, M., Matthews, J., & Lawson, T. (2019). Exploiting natural variation and genetic manipulation of stomatal conductance for crop improvement. Current Opinion in Plant Biology, 49, 17. https://doi.org/10.1016/j.pbi.2019.01.003 CrossRefGoogle ScholarPubMed
Faria, T., García-Plazaola, J. I., Abadía, A., Cerasoli, S., Pereira, J. S., & Chaves, M. M. (1996). Diurnal changes in photoprotective mechanisms in leaves of cork oak (Quercus suber) during summer. Tree Physiology, 16, 115123. https://doi.org/10.1093/treephys/16.1-2.115 CrossRefGoogle ScholarPubMed
Feng, L., Raza, M. A., Li, Z., Chen, Y., Khalid, M. H. B., Du, J., Liu, W., Wu, X., Song, C., Yu, L., Zhang, Z., Yuan, S., Yang, W., & Yang, F. (2018). The influence of light intensity and leaf movement on photosynthesis characteristics and carbon balance of soybean. Frontiers in Plant Science, 9, 1952. https://doi.org/10.3389/fpls.2018.01952 CrossRefGoogle ScholarPubMed
Franks, P. J., & Beerling, D. J. (2009). Maximum leaf conductance driven by CO2 effects on stomatal size and density over geologic time. Proceedings of the National Academy of Sciences of the United States of America, 106, 1034310347. https://doi.org/10.1073/pnas.0904209106 CrossRefGoogle ScholarPubMed
Franks, P. J., Drake, P. L., & Beerling, D. J. (2009). Plasticity in maximum stomatal conductance constrained by negative correlation between stomatal size and density: An analysis using Eucalyptus globulus . Plant, Cell & Environment, 32, 17371748. https://doi.org/10.1111/j.1365-3040.2009.02031.x CrossRefGoogle ScholarPubMed
Franks, P. J., & Farquhar, G. D. (2001). The effect of exogenous abscisic acid on stomatal development, stomatal mechanics, and leaf gas exchange in Tradescantia virginiana . Plant Physiology, 125, 935942. https://doi.org/10.1104/pp.125.2.935 CrossRefGoogle ScholarPubMed
Franks, P. J., & Farquhar, G. D. (2007). The mechanical diversity of stomata and its significance in gas-exchange control. Plant Physiology, 143, 7887. https://doi.org/10.1104/pp.106.089367 CrossRefGoogle ScholarPubMed
Giuliani, R., Koteyeva, N., Voznesenskaya, E., Evans, M. A., Cousins, A. B., & Edwards, G. E. (2013). Coordination of leaf photosynthesis, transpiration, and structural traits in rice and wild relatives (Genus Oryza). Plant Physiology, 162, 16321651. https://doi.org/10.1104/pp.113.217497 CrossRefGoogle Scholar
Gommers, C. (2020). Keep cool and open up: Temperature-induced stomatal opening [Review of Keep Cool and Open Up: Temperature-Induced Stomatal Opening]. Plant Physiology, 182, 11881189. https://doi.org/10.1104/pp.20.00158 CrossRefGoogle Scholar
Grantz, D. A., & Assmann, S. M. (1991). Stomatal response to blue light: Water use efficiency in sugarcane and soybean*. Plant, Cell & Environment, 14, 683690. https://doi.org/10.1111/j.1365-3040.1991.tb01541.x CrossRefGoogle Scholar
Gray, A., Liu, L., & Facette, M. (2020). Flanking support: How subsidiary cells contribute to stomatal form and function. Frontiers in Plant Science, 11, 881. https://doi.org/10.3389/fpls.2020.00881 CrossRefGoogle ScholarPubMed
Haworth, M., Marino, G., Loreto, F., & Centritto, M. (2021). Integrating stomatal physiology and morphology: Evolution of stomatal control and development of future crops. Oecologia, 197, pages 867883. https://doi.org/10.1007/s00442-021-04857-3CrossRefGoogle ScholarPubMed
Haworth, M., Scutt, C. P., Douthe, C., Marino, G., Gomes, M. T. G., Loreto, F., Flexas, J., & Centritto, M. (2018). Allocation of the epidermis to stomata relates to stomatal physiological control: Stomatal factors involved in the evolutionary diversification of the angiosperms and development of amphistomaty. Environmental and Experimental Botany, 151, 5563. https://doi.org/10.1016/j.envexpbot.2018.04.010 CrossRefGoogle Scholar
Hill, K. E., Guerin, G. R., Hill, R. S., & Watling, J. R. (2015). Temperature influences stomatal density and maximum potential water loss through stomata of Dodonaea viscosa subsp. angustissima along a latitude gradient in southern Australia. Australian Journal of Botany, 62, 657665.CrossRefGoogle Scholar
Hughes, J., Hepworth, C., Dutton, C., Dunn, J. A., Hunt, L., Stephens, J., Waugh, R., Cameron, D. D., & Gray, J. E. (2017). Reducing stomatal density in barley improves drought tolerance without impacting on yield. Plant Physiology, 174, 776787. http://www.plantphysiol.org/lookup/doi/10.1104/pp.16.01844 CrossRefGoogle ScholarPubMed
Ikkonen, E. N., Shibaeva, T. G., Rosenqvist, E., & Ottosen, C. O. (2015). Daily temperature drop prevents inhibition of photosynthesis in tomato plants under continuous light. Photosynthetica, 53, 389394. https://doi.org/10.1007/s11099-015-0115-4 CrossRefGoogle Scholar
Jahan, E., Amthor, J. S., Farquhar, G. D., Trethowan, R., & Barbour, M. M. (2014). Variation in mesophyll conductance among Australian wheat genotypes. Functional Plant Biology: FPB, 41, 568580. https://doi.org/10.1071/FP13254 CrossRefGoogle ScholarPubMed
Jezek, M., & Blatt, M. R. (2017). The membrane transport system of the guard cell and its integration for stomatal dynamics. Plant Physiology, 174, 487519. https://doi.org/10.1104/pp.16.01949 CrossRefGoogle ScholarPubMed
Kardiman, R., & Ræbild, A. (2018). Relationship between stomatal density, size and speed of opening in Sumatran rainforest species. Tree Physiology, 38, 696705. https://doi.org/10.1093/treephys/tpx149 CrossRefGoogle ScholarPubMed
Kollist, H., Nuhkat, M., & Roelfsema, M. R. G. (2014). Closing gaps: Linking elements that control stomatal movement. The New Phytologist, 203, 4462. https://doi.org/10.1111/nph.12832 CrossRefGoogle ScholarPubMed
Lamba, S., Hall, M., Räntfors, M., Chaudhary, N., Linder, S., Way, D., Uddling, J., & Wallin, G. (2018). Physiological acclimation dampens initial effects of elevated temperature and atmospheric CO2 concentration in mature boreal Norway spruce. Plant, Cell & Environment, 41, 300313. https://doi.org/10.1111/pce.13079 CrossRefGoogle ScholarPubMed
Lawson, T., & Blatt, M. R. (2014). Stomatal size, speed, and responsiveness impact on photosynthesis and water use efficiency. Plant Physiology, 164, 15561570. https://doi.org/10.1104/pp.114.237107 CrossRefGoogle ScholarPubMed
Lawson, T., & Matthews, J. (2020). Guard cell metabolism and stomatal function. Annual Review of Plant Biology, 71, 273302. https://doi.org/10.1146/annurev-arplant-050718-100251 CrossRefGoogle ScholarPubMed
Lawson, T., & Vialet-Chabrand, S. (2019). Speedy stomata, photosynthesis and plant water use efficiency. The New Phytologist, 221, 9398. https://doi.org/10.1111/nph.15330 CrossRefGoogle ScholarPubMed
Leakey, A. D. B., Ferguson, J. N., Pignon, C. P., Wu, A., Jin, Z., Hammer, G. L., & Lobell, D. B. (2019). Water use efficiency as a constraint and target for improving the resilience and productivity of C3 and C4 crops. Annual Review of Plant Biology, 70, 781808. https://doi.org/10.1146/annurev-arplant-042817-040305 CrossRefGoogle ScholarPubMed
Liu, C., He, N., Zhang, J., Li, Y., Wang, Q., Sack, L., & Yu, G. (2018). Variation of stomatal traits from cold temperate to tropical forests and association with water use efficiency. Functional Ecology, 32, 2028. https://doi.org/10.1111/1365-2435.12973 CrossRefGoogle Scholar
Matthews, J. S. A., Vialet-Chabrand, S. R. M., & Lawson, T. (2017). Diurnal variation in gas exchange: The balance between carbon fixation and water loss. Plant Physiology, 174, 614623. https://doi.org/10.1104/pp.17.00152 CrossRefGoogle ScholarPubMed
McAusland, L., Vialet‐Chabrand, S., Davey, P., Baker, N. R., Brendel, O., & Lawson, T. (2016). Effects of kinetics of light‐induced stomatal responses on photosynthesis and water‐use efficiency. The New Phytologist, 211, 12091220. https://doi.org/10.1111/nph.14000 CrossRefGoogle ScholarPubMed
Merilo, E., Jõesaar, I., Brosché, M., & Kollist, H. (2014). To open or to close: Species-specific stomatal responses to simultaneously applied opposing environmental factors. The New Phytologist, 202, 499508. https://doi.org/10.1111/nph.12667 CrossRefGoogle ScholarPubMed
Miao, Y., Cai, Y., Wu, H., & Wang, D. (2021). Diurnal and seasonal variations in the photosynthetic characteristics and the gas exchange simulations of two rice cultivars grown at ambient and elevated CO2. Frontiers in Plant Science, 12, 651606. https://doi.org/10.3389/fpls.2021.651606 CrossRefGoogle ScholarPubMed
Motulsky, H. J., & Brown, R. E. (2006). Detecting outliers when fitting data with nonlinear regression – A new method based on robust nonlinear regression and the false discovery rate. BMC Bioinformatics, 7, 123. https://doi.org/10.1186/1471-2105-7-123 CrossRefGoogle ScholarPubMed
Murata, Y., Mori, I. C., & Munemasa, S. (2015). Diverse stomatal signaling and the signal integration mechanism. Annual Review of Plant Biology, 66, 369392. http://www.annualreviews.org/doi/10.1146/annurev-arplant-043014-114707 CrossRefGoogle ScholarPubMed
Nunes, T. D. G., Zhang, D., & Raissig, M. T. (2020). Form, development and function of grass stomata. The Plant Journal: For Cell and Molecular Biology, 101, 780799. https://doi.org/10.1111/tpj.14552 CrossRefGoogle ScholarPubMed
Qi, X., & Torii, K. U. (2018). Hormonal and environmental signals guiding stomatal development. BMC Biology, 16, 21. https://doi.org/10.1186/s12915-018-0488-5 CrossRefGoogle ScholarPubMed
Raissig, M. T., Matos, J. L., Anleu Gil, M. X., Kornfeld, A., Bettadapur, A., Abrash, E., Allison, H. R., Badgley, G., Vogel, J. P., Berry, J. A., & Bergmann, D. C. (2017). Mobile MUTE specifies subsidiary cells to build physiologically improved grass stomata. Science, 355, 12151218. https://doi.org/10.1126/science.aal3254 CrossRefGoogle ScholarPubMed
Raschke, K., & Fellows, M. P. (1971). Stomatal movement in Zea mays: Shuttle of potassium and chloride between guard cells and subsidiary cells. Planta, 101, 296316. https://doi.org/10.1007/BF00398116 CrossRefGoogle ScholarPubMed
Roussel, M., Dreyer, E., Montpied, P., Le-Provost, G., Guehl, J.-M., & Brendel, O. (2009). The diversity of 13C isotope discrimination in a Quercus robur full-sib family is associated with differences in intrinsic water use efficiency, transpiration efficiency, and stomatal conductance. Journal of Experimental Botany, 60, 24192431. https://doi.org/10.1093/jxb/erp100 CrossRefGoogle Scholar
Sakoda, K., Yamori, W., Shimada, T., Sugano, S. S., Hara-Nishimura, I., & Tanaka, Y. (2020). Higher stomatal density improves photosynthetic induction and biomass production in arabidopsis under fluctuating light. Frontiers in Plant Science, 11, 589603. https://doi.org/10.3389/fpls.2020.589603 CrossRefGoogle ScholarPubMed
Sierla, M., Waszczak, C., Vahisalu, T., & Kangasjärvi, J. (2016). Reactive oxygen species in the regulation of stomatal movements. Plant Physiology, 171, 15691580. https://doi.org/10.1104/pp.16.00328 CrossRefGoogle ScholarPubMed
Sonawane, B. V., Koteyeva, N. K., Johnson, D. M., & Cousins, A. B. (2021). Differences in leaf anatomy determines temperature response of leaf hydraulic and mesophyll CO 2 conductance in phylogenetically related C 4 and C 3 grass species. New Phytologist, 230, 18021814. https://doi.org/10.1111/nph.17287 CrossRefGoogle Scholar
Stangl, Z. R., Tarvainen, L., Wallin, G., Ubierna, N., Räntfors, M., & Marshall, J. D. (2019). Diurnal variation in mesophyll conductance and its influence on modelled water-use efficiency in a mature boreal Pinus sylvestris stand. Photosynthesis Research, 141, 5363. https://doi.org/10.1007/s11120-019-00645-6 CrossRefGoogle Scholar
Stebbins, G. L., & Shah, S. S. (1960). Developmental studies of cell differentiation in the epidermis of monocotyledons. II. Cytological features of stomatal development in the Gramineae. Developmental Biology, 2, 477500.CrossRefGoogle Scholar
Terfa, M. T., Olsen, J. E., & Torre, S. (2020). Blue light improves stomatal function and dark-induced closure of rose leaves (Rosa x hybrida) developed at high air humidity. Frontiers in Plant Science, 11, 1036. https://doi.org/10.3389/fpls.2020.01036 CrossRefGoogle Scholar
Urban, J., Ingwers, M., McGuire, M. A., & Teskey, R. O. (2017a). Stomatal conductance increases with rising temperature. Plant Signaling & Behavior, 12, e1356534. https://doi.org/10.1080/15592324.2017.1356534 CrossRefGoogle Scholar
Urban, J., Ingwers, M. W., McGuire, M. A., & Teskey, R. O. (2017b). Increase in leaf temperature opens stomata and decouples net photosynthesis from stomatal conductance in Pinus taeda and Populus deltoides x nigra . Journal of Experimental Botany, 68, 17571767. https://doi.org/10.1093/jxb/erx052 CrossRefGoogle Scholar
Vahisalu, T., Kollist, H., Wang, Y.-F., Nishimura, N., Chan, W., Valerio, G., Lamminmäki, A., Brosché, M., Moldau, H., Desikan, R., Schroeder, J. I., & Kangasjärvi, J. (2008). SLAC1 is required for plant guard cell S-type anion channel function in stomatal signalling. Nature, 452, 487491. https://doi.org/10.1038/nature06608 CrossRefGoogle ScholarPubMed
Yamasaki, T., Yamakawa, T., Yamane, Y., Koike, H., Satoh, K., & Katoh, S. (2002). Temperature acclimation of photosynthesis and related changes in photosystem II electron transport in winter wheat. Plant Physiology, 128, 10871097. https://doi.org/10.1104/pp.010919 CrossRefGoogle ScholarPubMed
Yamori, W., Suzuki, K., Noguchi, K. O., Nakai, M., & Terashima, I. (2006). Effects of Rubisco kinetics and Rubisco activation state on the temperature dependence of the photosynthetic rate in spinach leaves from contrasting growth temperatures. Plant, Cell and Environment, 29, 16591670. https://doi.org/10.1111/j.1365-3040.2006.01550.x CrossRefGoogle ScholarPubMed
Zhang, L., Wang, S., Yang, X., Cui, X., & Niu, H. (2021). An intrinsic geometric constraint on morphological stomatal traits. Frontiers in Plant Science, 12, 628. https://doi.org/10.3389/fpls.2021.658702 Google ScholarPubMed
Zhang, Q., Peng, S., & Li, Y. (2019). Increase rate of light-induced stomatal conductance is related to stomatal size in the genus Oryza. Journal of Experimental Botany, 70, 52595269. https://doi.org/10.1093/jxb/erz267 CrossRefGoogle ScholarPubMed
Zhou, H., Xu, M., Pan, H., & Yu, X. (2015). Leaf-age effects on temperature responses of photosynthesis and respiration of an alpine oak, Quercus aquifolioides, in southwestern China. Tree Physiology, 35, 12361248. https://doi.org/10.1093/treephys/tpv101 CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Leaf-level gas exchange measurements in response to changing light intensities reveal fast and consistent stomatal movements in Brachypodium distachyon Bd21-3. (a) Experimental setup for measuring leaf-level gas exchange parameters related with CO2 capture and H2O vapor loss by clamping a leaf in an infrared gas analyser (IRGA) chamber with controlled environmental conditions. Gas exchange is measured in changing light conditions (1,000–100–1,000–0 PAR) inducing stomatal closure in response to decreasing light intensity/darkness and stomatal opening in response to increasing light intensity. (b) Relative stomatal conductance (Rel gsw) during the light transitions (n = 120, normalised to highest gsw observed). (c) Absolute stomatal conductance (gsw) response to light transitions (in black, data from 120 individuals corrected by average leaf area of 0.64 cm2 and in magenta data from a subset of 35 individuals corrected by individual leaf area). (d) Intrinsic water-use efficiency (iWUE) response to light transitions (1,000–100–1,000 PAR) (n = 120, calculated as A/gsw). (e) One-phase decay exponential regression for the transition 1,000–100 PAR (n = 120). (f) One-phase association exponential regression for the transition 100–1,000 PAR (n = 120). (g) One-phase decay exponential regression for the transition 1,000–0 PAR (n = 120). R2, half-time (T±50%) and rate constant (K) are indicated. Error bars = SD.

Figure 1

Figure 2. Influence of temperature, ambient light and diurnal rhythm on leaf-level gas exchange. (A) Correlation matrix between iWUE and environment (Qout, T, RH and time) of 120 measurements of wild-type B. distachyon (Bd21-3) at 100 PAR (second light step). (b) Linear regression between T and iWUE at 100 PAR (n = 116). (c) Non-linear regression between time and iWUE at 100 PAR (n = 116). (d) Correlation matrix between iWUE and environment (Qout, T, RH and time) of 120 measurements of wild-type B. distachyon (Bd21-3) at 1,000 PAR (third light step). (e) Linear regression between T and iWUE at 1,000 PAR (n = 119). (f) Non-linear regression between time and iWUE at 1,000 PAR (n = 120). (g) Correlation matrix between gas exchange parameters (A, gsw and iWUE) and environment (Qout, T, RH and time) of the 35 measurements (corrected by individual leaf area) of wild-type B. distachyon (Bd21-3) at 100 PAR (second light step). (h) Correlation matrix between gas exchange parameters (A, gsw and iWUE) and environment (Qout, T, RH and time) of the 35 measurements (corrected by individual leaf area) of wild-type B. distachyon (Bd21-3) at 1,000 PAR (third light step). (i) Linear regression between T and gsw at 100 PAR (n = 35). (j) Linear regression between T and gsw at 1,000 PAR (n = 35). (k) Linear regression between Qout and A at 1,000 PAR (n = 35). (l) Linear regression between Qout and gsw at 1,000 PAR (n = 35). (m) Non-linear regression between half-time of the transition 1,000–100 PAR (T−50%I) and time of the day (time) (n = 99). (n) Non-linear regression between half-time of the transition 100–1,000 PAR (T+50%II) and time of the day (time) (n = 111). (o) Non-linear regression between half-time of the transition 1,000–0 PAR (T−50%III) and time of the day (time) (n = 111). R2 and Sy.x or p values are indicated.

Figure 2

Figure 3. Effects of seasonal growth conditions on stomatal anatomical traits and on gas exchange. (a) Correlation matrix between stomatal length (SL), stomatal density (SD), average growth temperature (T), average growth relative humidity (RH) and day length (DL) (n = 15). (b) Non-linear quadratic (second order polynomial regression) relation between T and SD (n = 15). (c) Non-linear quadratic (second order polynomial regression) relation between T and SL (n = 15). (d) Linear relation between SD and SL (n = 15). (e) Non-linear quadratic (second order polynomial regression) relation between RH and SD (n = 15). (f) Non-linear quadratic (second order polynomial regression) relation between RH and SL (n = 15). (g) Seasonal variation on A at 1,000 (white) and 100 (grey) PAR (n = 4–6 per season; n = 5 in summer, n = 6 in autumn, n = 4 in winter). (h) Seasonal variation on gsw at 1,000 and 100 PAR (n = 4–6 per season; n = 5 in summer, n = 6 in autumn, n = 4 in winter). (i) Seasonal variation on iWUE at 1,000 and 100 PAR (n = 4–6 per season; n = 5 in summer, n = 6 in autumn, n = 4 in winter). R2 and Sy.x or p values are indicated. Dashed lines in (B-F) indicate 95% confidence bands.

Figure 3

Figure 4. Impact of stomatal anatomical traits on steady-state gas exchange and stomatal kinetics. (a) Linear regressions between stomatal length (SL) and gsw at 1,000 (white dots), 100 (grey dots) and 0 (black dots) PAR (n = 15). (b) Linear regressions between stomatal density (SD) and gsw at 1,000 (white dots), 100 (grey dots) and 0 (black dots) PAR (n = 15). (c) Linear regressions between SL and A at 1,000 (white dots) and 100 PAR (grey dots). (d) Linear regression between SD and A at 1,000 (white dots) and 100 PAR (grey dots). (e) Linear regressions between SL and iWUE at 1,000 PAR (white dots) and 100 PAR (grey dots) (n = 15). (f) Linear regressions between SD and iWUE at 1,000 PAR (white dots) and 100 PAR (grey dots) (n = 15). (g) Linear regressions between SL and half-time (T50%) of the light transitions 1,000–100 (grey dots), 100–1,000 (white dots) and 1,000–0 (black dots) PAR (n = 15). (h) Linear regressions between SD and T50% of the light transitions 1,000–100 (grey dots), 100–1,000 (white dots) and 1,000–0 (black dots) PAR (n = 15). R2 and p values are indicated. Dashed lines indicate statistically non-significant linear regressions (p > 0.05).

Figure 4

Figure 5. Morphometric analysis of graminoid B. distachyon stomata significantly improves anatomical gsmax predictions. (a) Arabidopsis-like stoma and ellipse pore shape. (b) B. distachyon stomatal morphology traits measured; guard cell length (GCL), pore length (PL), pore width (PW), guard-cell width at the centre of the stomata (GCWC), stomatal width at the apex (WA) and pore depth (l). Pore area hand-traced (red) or geometrically defined as an ellipse (orange) or a rectangle (blue). (c) Linear relation between hand-traced pore area and ellipse pore area. (d) Linear relation between hand-traced pore area and rectangle multiplied by 0.9. (e) Linear relation of anatomical gsmax calculated with hand-traced pore area and with ellipse pore. (f) Linear relation of anatomical gsmax calculated with hand-traced pore area and with rectangle pore multiplied by 0.9. (g) Anatomical parameters measured using light microscopy; stomatal width at the apex (WA) and GCL. (h) Anatomical maximum stomatal conductance (anatomical gsmax) equation as defined by Franks and Farquhar (2001) and B. distachyon adjustments to calculate amax (0.9*PL*PW), PL (0.44*GCL), PW (0.13*GCL) and l (WA*0.5). (i) Comparison between physiological gsmax, anatomical gsmax (light microscopy, LM) and anatomical gsmax (confocal microscopy, CM). (j) Comparison of anatomical gsmax calculated for summer, autumn and winter plants with stomatal anatomical traits represented in Supplementary Figure S3B,C. R2 and p values are indicated.

Supplementary material: File

Nunes et al. supplementary material

Nunes et al. supplementary material

Download Nunes et al. supplementary material(File)
File 4.7 MB

Author comment: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R0/PR1

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

Review: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R0/PR2

Conflict of interest statement

Reviewer declares none.

Comments

Comments to Author: The research question posed by authors – how does seasonal variation in growth conditions affect gas-exchange measurements in greenhouse settings – is interesting. The question whether and how much it matters, what other leaves experience while the measured leaf is clamped into a leaf gas exchange device is important, and the presence and extent of systemic signalling affecting plant stomatal behaviour is debated and needs further research.

The manuscript presents interesting data and offers a valuable method advancement for calculation of grass anatomical gsmax. It also opens up exciting questions for further research. However, I feel that the data that forms the basis for the first two figures and conclusions based on respective regression analyses are not strong enough. My main concern is that in gas exchange analysis, true leaf area has been measured for only 15 plants of the measured 100 individuals and for others gs and A values are estimates based on leaf area of these 15 plants. The estimate may be good enough but it may not be – and currently the manuscript does not clearly show that the approximation is valid. I have some other concerns regarding the analyses shown in the first two figures that are detailed below.

The rest of the manuscript that focuses on the relationships of anatomical stomatal parameters with growth conditions and physiology, and suggests an improved approximation for calculation of anatomical gsmax, does not have such issues and in my opinion would be an interesting addition to literature in the field.

Specific comments and questions:

p7 lines 3-7 – If I understand correctly then accurate stomatal conductance measurements with determined leaf area for the part of the leaf in the IRGA chamber have only been carried out for 15 plants and for others the area has been estimated, meaning gs is also an estimate? In this case, what were the conditions during the measurement of these 15 plants (light, temperature, RH, time of day) and do they cover the range covered by all 100 plants? If not, there may be a bias in the estimate of stomatal conductance; would be good to bring out the data for this subset in the condition graphs shown in S1D. If my understanding is correct, then the numbers on S1A do not reflect real gs, as they have been calculated with the assumption of 2 cm2 chamber leaf area for all plants? If this is the case, this data gives no information on absolute values of gs and should not be shown; same applies to A shown in S1B (if the leaf area correction is valid – that is, the 15 plants conditions reflect the conditions throughout all experiments, estimated A based on leaf area correction could be shown instead).

Figure 2 – I am confused, whether the A and gs values used in model fitting are the values assuming 2 cm2 leaf area or the estimated gs using leaf area correction? As the leaf area seems to have been much smaller than the 2 cm2 (based on the 15 plants with known area), the first approach is not justified in my opinion. If the 15 plants accurately represent the conditions and size for all plants, gs and Anet estimates could be better for looking at these relationships. If real leaf area was known only for 15 plants, would be good to test if the relationships shown in Figure 2 hold or are similar for the 15 plants dataset. I realise it will be less powerful, but at the moment the data is a mix of known and estimated values, making it hard to make strong conclusions. From p8 lines 3-4 I understand that for the subset with known leaf area (15 plants), only the regressions between temperature and light with A and iWUE were significant and R2 values were larger than when these 15 plants were combined with data from plants with inaccurate leaf area – this makes me doubt that the regression models presented in Figure 2 are accurately representing the relationships in the 100 measured plants.

Another thing I am wondering about in Figure 2 and S2 is using multiple data points for one plant in analysis (conductance/A at 1000 PAR measured twice for each plant during the experiment) – I think it can be considered pseudoreplication that artificially increases n and can potentially bias the results of analysis.

Figure 3E and interpretation on p9, line 1 - The negative relationship between RH and SD is unexpected, as a body of evidence suggests that higher RH usually stimulates stomatal development (of course, the tested RHs then are usually much higher). It has been mostly studied in dicots, so maybe Brachypodium really is different. Or alternatively – the higher light and longer days of summer (S3) could stimulate stomatal development and the high SD at lower RH of summer would be caused by light instead.

Figure 4 – for nonsignificant relationships showing the trend line does not make much sense – as the slope is not significantly different from 0, the line is not helpful and can be misleading. I also got confused about the rate constant and half time – to my understanding they can be calculated from each other, making showing regressions of SL and SD with both of them redundant.

Minor comments:

S1D referred before any other panels in S1 – would be logical to make it A then.

I could not find reference to Figure 1 panel A

Review: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R0/PR3

Conflict of interest statement

Reviewer declares none.

Comments

Comments to Author: Nunes et al., characterized light response kinetics in 100 wild-type Brachypodium distachyon plants grown in semi controlled greenhouse conditions for almost two years. Authors used IRGA-based system for their measurements in which single leaf of the plant is positioned to the IRGA chamber with controlled environmental conditions. Authors also monitored time of the day and growth conditions (temperature, light intensity and humidity in the greenhouse) and quantified effects of these parameters on stomatal gas-exchange. Authors find that stomatal conductance was affected by temperature and light conditions in the greenhouse whereas stomatal kinetics was affected by time of the day. Previous studies have suggested that in some dicot species conditions experienced by distal parts of the plant may modify local stomatal response during IRGA measurements. Current work suggests that systemic signaling may modulate stomatal gas-exchange dynamics also in monocots.

Authors analyzed stomatal anatomical traits and showed that in Brachypodium stomatal anatomy shows clear seasonal plasticity. During the summer stomatal density was higher and it was accompanied by reduced stomatal size and this was correlated with faster stomatal closure and enhanced intrinsic water use efficiency. Finally, authors performed morphometric analysis by using confocal microscopy and defined pore dimensions for Brachypodium stomata. Further, authors developed accurate method to measure these parameters also from light microscope images, which are more convenient for large set of samples. By using these parameters, authors optimized anatomical g<sub>s</sub>max formula for grass-type stomata, and showed that it correlates well with their experimental data.

This impressive work describes in detail seasonal behavior of stomata in wild-type Brachypodium and reveals how sensitive stomatal physiology and especially stomatal anatomy are to subtle changes in the growth conditions. Experiments and analyses performed are thorough, conclusions are solid, and the manuscript is generally well written.

Minor comments:

1. Currently, it is indicated that n=15 for leaves used for anatomical analyses both in introduction and material and methods. However, it seems that for each season sample number varies between 4-6 and this is correctly indicated in the figure legends. This is slightly confusing. Perhaps it would be more straightforward to indicate sample number per season similarly throughout the article.

2. Page 9. “When measuring physiological maximum stomatal conductance (g<sub>s</sub>max) 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<sub>s</sub>max was observed in autumn/winter (Fig. S3I).”

I am curious to hear why authors choose to pool subset of autumn and winter samples for this particular analysis? (based on Figure S3 legend, 3 autumn samples and 1 winter sample have been pooled)

3. Page 9. “an increase in SL and a decrease in SD resulted in an increase of A at high light, while no effect was observed in 100 PAR (light limiting condition) (Fig.4C-D). “

Based on the data in Figure 4C-D the trends seems opposite to what was described in the manuscript. Please correct the sentence so that it will be in line with the data.

4. Page 13. “Recent studies have also reported systemic stomatal responses to light and heat in Arabidopsis (Devireddy et al., 2018, 2020; Ehonen et al., 2020). “

Ehonen et al., 2020 detected systemic effect in birch and poplar as a response to darkness and elevated CO2 treatments but actually failed to see similar effect in Arabidopsis. Please correct the sentence or the references so that the content is in line with referred articles.

5. Please indicate where these experiment have been performed. Perhaps it would be also useful to shortly explain typical seasonal weather in this place. This might be helpful for the reader since seasons can be dramatically different in different parts of the world.

Recommendation: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R0/PR4

Comments

Comments to Author: Dear Dr. Raissig,

two reviewers have now read your manuscript, and as you can see, both of them found the work novel and interesting. However, they also found a few major and several minor concerns that need to be addressed. Especially, the concern related to the low number of plants in the gas exchange analysis is critical: is there an evidence to support the idea that 15 plants is sufficient?

Looking forward to handle the revised version of your manuscript.

best wishes,

Ari Pekka Mähönen

Decision: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R0/PR5

Comments

No accompanying comment.

Author comment: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R1/PR6

Comments

No accompanying comment.

Review: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R1/PR7

Comments

Comments to Author: I think the manuscript have been greatly improved by revisions and most of my comments have been addressed.

Description of the sample number used in seasonal analysis (Figure 3) is still in part confusing. Number of samples used for this particular experiment is low and thus, it would be important to communicate it accurately.

Authors describe seasonal variation on leaf anatomy and gas exchange by using following sample groups: summer; n=5, autumn; n=6 and winter; n=4 (Figure 3). In Figure 3 panels G, H, and I Authors compare extent of seasonal variation in gas exchange traits between each seasonal groups. In the Figure 3 legend, sample number has been marked as 15 (S+A+W). This seems unusual way to describe n in this case (Figure 3 G-I), since here seasonal variation within each group is compared to the variation within another group (S vs A, S vs W, A vs W). Please indicate n so that it accurately describes number of samples in each group.

Number of data points for 1000 PAR conditions seems to be twice as high as for 100 PAR conditions for each groups. Please explain reason for this and also indicate whether or not these 1000 PAR measurements are paired. If relevant, consider also whether the statistical test used is correct.

Review: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R1/PR8

Comments

Comments to Author: The authors have carried out additional experiments and analyses to respond to questions raised in my previous review. In my opinion, the increased sample size and adjusted analyses strengthen the conclusions presented in the manuscript. The authors have found significant systemic effects of ambient conditions on gas exchange parameters measured in environmentally-controlled leaf chamber, which is an important point to consider in leaf gas exchange analyses carried out in fluctuating environments.

I have no more major questions, just a few minor points.

In Figure 1D legend: I think it should read 1000 - 100 - 1000 PAR in the brackets

Line 227: “high light steady-state (1000 PAR, Fig 2A)” – should be low light? It says 100 PAR in Fig 2A

Line 228: “steady-state in darkness (0 PAR, Fig S2A)” – should be reference to S2E?

Line 246: “influence of T on gsw is suggested (Fig. S2G)” – should be reference to S2F?

In paragraphs in lines 250-256 the authors discuss about the effects of circadian rhythms on gas-exchange, as time of day affected gas-exchange measurements. I think the term “diurnal” would be more appropriate than “circadian” about these effects. The same also applies to Figure 2 legend.

Recommendation: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R1/PR9

Comments

Comments to Author: Dear Michael,

both reviewers were satisfied with your response to their concerns and overall they appeared pleased with the revised manuscript. They both had, though, a few minor, additional comments that sounded reasonable to me. Please, address these issues, and send a revised version of the manuscript in due time.

Looking forward to hearing from you soon

best wishes,

Ari Pekka

Decision: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R1/PR10

Comments

No accompanying comment.

Author comment: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R2/PR11

Comments

No accompanying comment.

Recommendation: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R2/PR12

Comments

No accompanying comment.

Decision: Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model — R2/PR13

Comments

No accompanying comment.