Original Research ARTICLE
Heterogeneity of Stomatal Pore Area Is Suppressed by Ambient Aerosol in the Homobaric Species, Vicia faba
- 1Department of Botany and Plant Sciences, Kearney Agricultural Center, University of California, Riverside, Riverside, CA, United States
- 2Institute of Crop Science and Resource Conservation, University of Bonn, Bonn, Germany
Stomatal pore area is heterogeneous across leaf surfaces. This has been considered as “patchy stomatal conductance,” and may have substantial implications for photosynthetic efficiency. Aerosols have always been important elements of plant environments, but their effects on stomatal control of plant water relations, and stomatal heterogeneity specifically, have not been considered. Here we evaluate the spatial coordination of pore area in the glabrous and homobaric leaves of Vicia faba grown under two aerosol treatments and measured at four levels of VPD. We construct a large dataset (n > 88,000 discrete comparisons) of paired pore areas and distances between the pores. Plants were grown in ambient urban air and in filtered air (FA) to determine the effect of ambient aerosol on stomatal properties. Pore area exhibited spatial organization, as well as considerable variability among closely co-located pores. The difference between pore areas was positively correlated with the distance between the pores, in both aerosol treatments and at all VPDs. However, aerosol deposition reduced both the magnitude of variability between pores and the rate at which this variability increased with pore separation distance. These data support previous conclusions that deposition of hygroscopic aerosol may create a thin aqueous film across the leaf surface that connects neighboring stomata to each other and to the leaf interior. Aerosol impacts on stomatal heterogeneity and gas exchange are not adequately considered in current assessments of stomatal control.
Heterogeneous stomatal opening across the surface of individual leaves may result in stomatal “patchiness,” random variation, or spatially coherent trends in pore area across the surface. Stomatal heterogeneity has been observed in many species (Downton et al., 1988, 1990; Terashima et al., 1988; Sharkey and Seemann, 1989; Terashima, 1992; Pospisilova and Santrucek, 1994; Haefner et al., 1997; Weyers and Lawson, 1997; Beyschlag and Eckstein, 1998; Mott and Buckley, 1998, 2000). Non-uniform pore areas and resulting uneven conductance for CO2 and water vapor complicate standard measurements and calculations of gas exchange parameters. Heterogeneity observed at all stages of leaf elongation in Rosa shows that this is not a transient developmental feature (Fanourakis et al., 2015). Under conditions of high boundary layer conductance, stomatal heterogeneity may be detrimental to gas exchange efficiency, but under conditions of low wind or large leaves, which reduce boundary layer conductance, and under conditions of low overall stomatal conductance, heterogeneity may improve photosynthetic efficiency (Buckley et al., 1999). This reflects the effects of evaporative cooling of the leaf on fluxes and gradients of both water and CO2, and the transport efficiency of those pores that remain widely open. The significance and ubiquity of stomatal heterogeneity and its relationship with environmental conditions (Cheeseman, 1991; Gunasekera and Berkowitz, 1992; Cardon et al., 1994) as well as the role of aerosol deposition (Burkhardt and Grantz, 2017) in generating such heterogeneity have not been adequately considered.
Areas of coordinated stomatal behavior may range from a few mm to a few cm in extent (Terashima, 1992; Mott and Buckley, 1998, 2000; Mott and Peak, 2007) and are often bounded by veins, particularly in heterobaric leaves (Siebke and Weis, 1995; Haefner et al., 1997; Weyers and Lawson, 1997; Mott and Buckley, 2000). Patterns of heterogeneity can be transient, with patches of coherent stomatal response migrating within and even between areoles of heterobaric leaves (Kamakura et al., 2012), even though gas diffusion is restricted to areoles defined by vasculature with bundle sheath extensions (Terashima, 1992; Mott and Buckley, 2000).
Spatially coherent stomatal behavior is also observed in homobaric leaves (Kappen et al., 1987; Loreto and Sharkey, 1990; Mott and Parkhurst, 1991; Terashima, 1992; Eckstein et al., 1996; Mott and Buckley, 2000; Kamakura and Furukawa, 2008). This heterogeneity is less patchy and more characterized by trends across the leaf surface, from leaf base to tip (Nardini et al., 2008), margin to midrib (Weyers and Lawson, 1997), and more generally across the lamina (Terashima et al., 1988). Patches observed as 14CO2 fixation in homobaric V. faba (Terashima et al., 1988) were larger than in many heterobaric species and could be considered gradients across the leaf, with greater heterogeneity within the patches than in heterobaric species (Spence, 1987; Terashima et al., 1988). Stomatal aperture in homobaric Commelina communis (Smith et al., 1989) exhibited gradients as high as 3 μm per mm and 20 μm from leaf edge to midrib. This behavior may be observed even in epidermal peels exposed to uniform physical and chemical conditions. The distribution of pore opening often approximates a normal distribution (e.g., V. faba; Laisk et al., 1980). This variance within and between leaves obscures spatial patterns across individual leaves (Weyers and Lawson, 1997), so that heterogeneity in homobaric leaves remains poorly characterized.
Uneven levels of conductance across the surfaces of individual leaves have been associated with discrepancies in estimation of intercellular CO2 (Ci). Values calculated from measured fluxes of CO2 and water vapor exceeded directly measured values by up to 15% (Buckley et al., 1999; Boyer, 2015a,b; Hanson et al., 2016). This suggested a pathway for water efflux that was unavailable to CO2 influx, potentially peristomatal transpiration from around the stomatal pores (Maier-Maercker, 1983; Grantz, 1990), or cuticular transpiration from across the epidermis between the pores (Kerstiens, 1996; Boyer et al., 1997; Hanson et al., 2016). Alternatively, stomatal heterogeneity could distort the calculation of Ci (Terashima et al., 1988) as measurements of steady state gas exchange reflect the leaf-wise average of highly dynamic heterogeneous individual pores and patches of pores (Siebke and Weis, 1995).
Increased evaporative demand (leaf to air vapor pressure difference; VPD) reduces stomatal conductance and increases stomatal heterogeneity (Sharkey and Seemann, 1989; Downton et al., 1990; Haefner et al., 1997). Experimentally imposed local changes in humidity at the leaf surface altered the apertures of stomata within the affected area, and those of stomata up to 0.4 cm away, outside of the directly affected area. The linkage was apparently through epidermal turgor, although this mechanism may operate most effectively over short distances (Haefner et al., 1997; Mott et al., 1997, 1999; Mott and Franks, 2001). Variability in mesophyll biochemistry (Osmond et al., 1999) and xylem water relations may become more effective at larger scales of leaf or branch (Buckley and Mott, 2000).
Aerosol deposition (Burkhardt and Grantz, 2017) may affect both the asymmetric flux pathways for water and CO2 (Grantz et al., 2018) and stomatal heterogeneity. Deposition of hygroscopic, particularly chaotropic, aerosol (Tsigaridis et al., 2006; Pringle et al., 2010; Burkhardt and Grantz, 2017) reduces surface tension and results in development of thin liquid films on leaf surfaces (Eiden et al., 1994; Dutcher et al., 2010; Burkhardt and Hunsche, 2013; Burkhardt and Grantz, 2017; Fernandez et al., 2017). The liquid films penetrate into stomatal pores (Eichert et al., 1998, 2008; Basi et al., 2014; Kaiser, 2014), providing a liquid phase linkage between the saturated leaf apoplast and the dry atmosphere. This pathway is not under diffusional (i.e., stomatal) control and thus water loss from this pathway through evaporation at the leaf surface increases as VPD increases, even in species such as V. faba that exhibit strong closing response to increasing VPD (Pariyar et al., 2013; Grantz et al., 2018). The presence of liquid films can be visualized by electrical conductance measurements (Burkhardt and Eiden, 1994; Burkhardt et al., 1999; Burkhardt and Hunsche, 2013) and by electron micrography (Burkhardt and Grantz, 2017; Burkhardt et al., 2018; Grantz et al., 2018). Aerosol deposition reduces stomatal apertures (Burkhardt et al., 2001a; Burkhardt, 2010; Burkhardt and Grantz, 2017; Grantz et al., 2018) while increasing both water flux (Burkhardt et al., 2001a; Grantz et al., 2018) and minimum (cuticular) leaf conductance (Burkhardt and Pariyar, 2014, 2016; Grantz et al., 2018).
We hypothesized that aerosol-induced surface moisture may link individual pores across the leaf surface by providing a more homogeneous hydraulic and vapor pressure environment, thereby reducing stomatal heterogeneity despite potential desiccating effects of liquid phase water loss and the previously documented reduction in pore area (Grantz et al., 2018). Here we characterize the distribution of stomatal pore areas, a subject of previous consideration (Laisk et al., 1980; Spence, 1987; Gorton et al., 1989; Terashima, 1992; Weyers and Lawson, 1997) and the role of aerosol deposition on this heterogeneity, which has not previously been considered. We analyze a previously available dataset of 3600 direct microscopic pore area measurements (Grantz et al., 2018) to create a novel database of 88,200 discrete pore to pore comparisons of distances between the pores (d) and the differences between their pore areas (ΔA). We evaluate the distribution of pore areas, the local and larger scale heterogeneity of stomatal opening, and the effect of four levels of VPD and two levels of ambient aerosol on these characteristics.
Materials and Methods
Plants of Vicia faba (L.) were grown from seed as described previously (Grantz et al., 2018) in plastic pots in greenhouses at the University of Bonn, Germany. Plants received complete nutrient solution (Ferty 3, Planta Düngemittel GmbH, Hohenstauf, Germany) weekly and irrigation as needed.
Plants were randomly assigned either to a greenhouse ventilated with ambient air (AA) or an adjacent greenhouse ventilated with filtered air (FA). Filtration removed nearly all particles (Pariyar et al., 2013; Grantz et al., 2018). The aerosol was typical of ambient particulate matter in central Europe (about 35% ionic; see Grantz et al., 2018 and references therein). Aerosol concentration in the AA greenhouse (cloud chamber condensation nuclei counter; TSI 3783; TSI, Shoreview, MN, United States) was 6–7 × 109 particles m–3. This was reduced in the FA greenhouse by 99% to 5–10 × 106 m–3; confirming previous measurements (e.g., Pariyar et al., 2013). Other environmental parameters including temperature, relative humidity and concentrations of ozone were similar in the greenhouses (see Supplementary Figures S1, S2 in Grantz et al., 2018). Plants were exposed to natural daylength and sunlight (approximately 70% of ambient irradiance).
Measurements were obtained 5 weeks after planting on one leaflet of leaf 5 or 6 (youngest fully expanded; mean leaflet area about 32 cm2) when plants had been exposed for 3–4 weeks in the greenhouses. Leaflets remained attached to intact plants during all measurements.
Measurement of Pore Area at Controlled VPD
Plants were transported from Bonn to Kiel, Germany. Pore area measurements (see Grantz et al., 2018 for further details) were obtained while leaf to air VPD was controlled in a flow-through gas exchange system with an integrated inverted video microscope (Kaiser and Kappen, 1997, 2000; Kaiser, 2009). The gas exchange cuvette was held at constant temperature (25°C ± 0.1°C) and irradiance (PPFD = 450 ± 25 μmol m–2 s–1). Dew point was held at 23.15, 19.0, 14.0, and 5.0°C (±0.05°C), yielding VPD of 0.33, 0.97, 1.6, and 2.3 kPa and RH of 90%, 68%, 50%, and 27%, at abaxial leaf temperature (25°C; Type K thermocouple, 0.075 mm). Air circulation in the cuvette (1 m s–1) yielded laminar boundary layer conductance of 1300 mmol m–2 s–1 (Kaiser, 2009).
Nine plants were used for each of the AA and FA treatments (Experiment 1, Grantz et al., 2018). Measurements on plants subjected to the two aerosol treatments were alternated on successive days. Pore area was measured on 50 stomata per leaflet, on the abaxial surface, within the 1 cm2 area of a predefined grid. The sample area was located away from major veins and leaf margins. The 50 stomata represented about 0.2% of the stomata in the 1 cm2 viewing area. The location (pore center) and focus depth (narrowest part of the pore center) of each pore was stored electronically (Kaiser, 2009) to allow rapid re-imaging at each VPD.
Data were stored with plant, leaf and pore identifiers, pore area, and pore coordinates in the Cartesian framework defined across the microscope stage (Kaiser, 2009). The area of each pore was determined as an idealized ellipse determined by the directly measured length and width of each pore, taken as the major and minor axes (Grantz et al., 2018). The difference in pore area (ΔA) between pairs of stomata was evaluated as a function of their separation distance (d) across the leaf surface.
Values of d were calculated from the coordinates of the center of each pore using the Pythagorean Theorem executed in a custom Python program (v. 3.6; Python Software Foundation1). Unique pair comparisons [C(50,2) = (50)(50-1)/2 = 1225 per leaflet] were evaluated. With 9 leaflets/aerosol treatment, 2 treatments, and 50 pores/leaflet, there were 900 unique pores and with 4 levels of VPD there were 88,200 unique comparisons of ΔA vs. d. The maximum value of d was about 1.4 cm, representing the diagonal of the 1 cm2 square examined on each leaf.
Regressions of ΔA on d were calculated in Sigma Plot, v. 13.0 (Systat Software, San Jose, CA, United States; Regression Wizard). Normality [Kolmogorov–Smirnov (K–S) Test] and equality of variance (Brown–Forsythe Test) were generally not satisfied (P < 0.050). The sole exception was FA at low VPD, which was normally distributed. Mean ΔA (Table 1) and regression coefficients (Table 2) were evaluated by Two-Way ANOVA (Sigma Plot) with mean separation using the Holm-Sidak Multiple Range Test.
Table 2. Linear regression coefficients of the relationships between the difference in pore area (ΔA) and the distance between the pores (d).
The imaging system facilitated collection of a large number of clear photomicrographs (Figure 1). At all levels of stomatal opening, both within each VPD and across the levels of VPD, the pores were well characterized as ellipses, allowing the area (A) of each pore to be calculated from the directly measured length and width. Leaves of V. faba exhibited uniformly kidney-shaped guard cells (Figure 1), lacking specialized stomatal subsidiary cells and the bundle sheath extensions that play a role in stomatal patchiness in heterobaric species.
Figure 1. Typical stomatal pores whose areas (A) were measured. Representative example at low VPD (0.33 kPa). Micrograph dimensions 300 × 300 μm.
The distances (d) between pores were distributed approximately normally (Figure 2). Neither the shape of the distribution nor the magnitude of the distance scale differed between leaves exposed during leaf development to ambient aerosol (AA) or to FA (cf. Figures 2A,B). These characteristics were also stable across the levels of VPD (not shown), despite the potential for leaf shrinkage due to reduced leaf water content (not measured) at elevated VPD.
Figure 2. Frequency distribution of pore separations (d) among the 50 stomata measured on each of 9 leaves (n = 11,025 per panel). Distance classes: 0 ≤ d < 1500 to 10500 ≤ d < 12000, centered at 750–11250. For d ≥ 12000 (n = 20 for AA and 55 for FA) frequencies were pooled with 10500 ≤ d < 12000. (A) Ambient air; (B) filtered air.
The large data set available for this study allowed a robust characterization of the distributions of pore area (A) in the two aerosol treatments over a range of VPDs (Figure 3). At low VPD (0.33 kPa) the distributions of A in both AA (Figure 3A) and FA (Figure 3B) treatments were approximately bell shaped, although only the FA treatment (Figure 3B) formally satisfied the Kolmogorov–Smirnov test for normality (K–S statistic of 0.661). All other distributions exhibited K–S < 0.001.
Figure 3. Frequency distribution of pore areas (A) at different levels of evaporative demand (VPD), in leaves exposed to ambient (gray bars; A,C,E,G) or filtered (hatched bars; B,D,F,H) air. Area classes: 0 ≤ A < 50 to 400 ≤ A < 450, centered at 25–425. The smallest visible bars represent 1 pore.
The distributions of A were similar in AA and FA leaves at all levels of VPD (cf. Figures 3A–H). VPD did not differ during plant growth between the AA and FA treatments but was varied experimentally during measurement of pore areas. Skewing toward the origin increased with VPD (Figure 3) so that the quasi-bell shaped distribution was no longer evident at 1.6 kPa (Figures 3E,F) or 2.3 kPa (Figures 3G,H) in either aerosol treatment.
In most cases, increasing VPD reduced pore area, i.e., with greater deviation below the 1:1 line (Figure 4). In FA at 0.97 kPa the pattern was different. The less open pores closed in response to increased VPD, but pores initially exhibiting A > 190 μm2 (Figure 4A; red circles) did not close below the initial A. Individual stomata behaved consistently over the range of VPDs. Those with small pore areas at 0.33 kPa exhibited small areas at higher VPD. Many that were only slightly open at low VPD closed completely at higher VPD. With increasing VPD, the number of such stomata and their initial pore areas increased (horizontal part in lower left end of each curve; Figure 4).
Figure 4. The relationship between stomatal pore area at elevated VPD [(A) 0.97 kPa; (B) 1.6 kPa; (C) 2.3 kPa] and pore area measured at low VPD (0.33 kPa). The AA treatment is shown in black, and the FA treatment in red. The solid line represents the 1:1 line.
Leaves exposed to ambient aerosol also exhibited greater skewing toward the origin. This reduction in magnitude was reflected in reduced mean ΔA between pores in AA relative to FA observed at all levels of VPD (Table 1). In general, the AA treatment closed more substantially than the FA as VPD was increased, particularly at greater initial pore area at low VPD (Figure 4). Among stomata with areas below 180 μm2 at 0.33 kPa, there was little effect of aerosol on the relationships between pore area at low VPD and at higher VPD. This point of divergence between AA and FA was similar at all VPD (Figures 4A–C). Above this value of initial opening, pore areas in FA and AA diverged, with AA smaller than FA. This was observed from 0.33 kPa to 0.97 kPa (Figure 4A), even though the response of FA to VPD was minimal among the more open pores. The sensitivity to aerosol increased with pore area observed at low VPD, but did not change substantially with increasing VPD.
We had hypothesized that stomatal heterogeneity, both ΔA between closely co-located pores and the increase of ΔA with d between more distant pores, would be reduced by aerosol deposition. The suggested mechanism, involving liquid films on the leaf surface, is illustrated in Figure 5 showing an open stoma at the left and a closed stoma at the right. Regardless of pore area, water evaporates from the surface of the leaf (blue arrows), potentially drawing both from the surface water (Figure 5; blue lines) and from soil, tissue, and apoplastic storage (continuity indicated by green lines). Conventional transpiration from the apoplast (green arrows) occurs through the open pore (left side) but this diffusive transport is almost completely blocked by stomatal closure (right side).
Figure 5. The hypothesized role of thin liquid films in reducing microenvironmental variability across the leaf surface. Blue lines represent liquid formed by condensation to deposited hygroscopic aerosol on the leaf surface, including that which penetrates the stomatal pores; blue arrows represent the evaporative flux of water derived from this liquid. Green lines represent conventional apoplastic liquid water derived from the soil-root-xylem pathway passing through and around mesophyll cells (m); green arrows represent the much larger evaporative flux of this water when stomata are open (left side) and its elimination by stomatal closure (right side). Evaporation occurs from aerosol-condensed water with both open (left side) and closed (right side) stomata, as liquid continuity between leaf interior and exterior is maintained largely independently of stomatal pore area determined by guard cell (gc) turgor. The leaf boundary layer (BL) is uniformly humidified across the leaf surface and the osmotic environment of the cuticle (c) and underlying epidermis (e) is made more homogeneous by the presence of a liquid film. All fluxes (blue and green arrows) eventually pass through the BL and enter the mixed free atmosphere (TA), representing water loss from the leaf.
In contrast, a continuous liquid film developing from deliquescence of deposited hygroscopic aerosol may spread across large areas of the leaf and penetrate into the stomatal pores. This liquid path through the pore reduces stomatal control of water loss. This may degrade epidermal water status by enhancing water loss, while also reducing the variability in the hydraulic and humidity environments across the leaf surface (Figure 5). While reduced epidermal water content might increase stomatal heterogeneity, the uniform environment under the liquid film might lead to more uniform epidermal water status and reduce heterogeneity.
We observed a reduction in heterogeneity in the ambient aerosol treatment. The intercept of ΔA on d was significantly lower in AA than FA at all VPD (Table 2A), reflecting reduced heterogeneity among closely co-located pores (extrapolated to zero separation). Exposure to aerosol also significantly reduced both mean ΔA (Table 1 and Figure 6) and median ΔA (Figure 7) at all levels of VPD. The considerable local variability is illustrated by the closely positioned pores in Figure 1, exhibiting ΔA of about 25% despite the guard cells bordering on the same epidermal cell.
Figure 6. Relationships between the paired difference between pore area (ΔA) and the distance separating the pores (d) in ambient air (AA; gray bars; A,C,E,G) and filtered air (FA; hatched bars, B,D,F,H) at different levels of evaporative demand (VPD), presented as mean ± s.e. Distance classes as in Figure 2.
Figure 7. Relationships between the paired difference between pore area (ΔA) and the distance separating the pores (d) in ambient air (AA; gray bars; A,C,E,G) and filtered air (FA; hatched bars, B,D,F,H) at different levels of evaporative demand (VPD), presented as median. Distance classes as in Figure 2.
Although there were fewer pores at small and large separations (Figure 2), the distributions of both mean ΔA (Figure 6) and median ΔA (Figure 7) revealed a consistent increase with d over the entire range. This was observed in both AA and FA leaves, but the increase in ΔA with d was reduced by aerosol (Figures 6, 7). The slope of ΔA on d was significantly reduced in AA relative to FA at VPD of 0.33 kPa and 1.6 kPa (Table 2B; c.f. Figures 6G,H, 7G,H). Slopes did not differ at 0.97 kPa or 2.3 kPa. This reflects a general reduction in the increase in heterogeneity with distance within the 1 cm2 sample area caused by aerosol exposure. The median analysis is presented, along with the mean, because the data often varied from normality.
The impact of VPD was less consistent. In AA, mean ΔA (Table 1) and the intercept of ΔA at d = 0 (Table 2A) increased with initial increase of VPD from 0.33 kPa through 0.97 to 1.6 kPa but declined at 2.3 kPa. In contrast, in FA, the mean and intercept of ΔA increased with VPD from 0.33 to 0.97 kPa, but began to decline at 1.6 kPa and declined further at 2.3 kPa. The pattern of response of ΔA over all pores (mean) and among closely spaced pores (intercept) was similar in the aerosol treatments, but the decline began at lower VPD in the absence of aerosol. In contrast, the slope of ΔA on d, reflecting larger scale heterogeneity, did not differ consistently between levels of VPD (Table 2B).
Heterogeneity of stomatal pore area was composed of two components, variability among pore areas that was not related to the distance between the pores (Table 2A; left-most bars, Figures 6, 7), and variability which trended systematically with d (Table 2B; distributions, Figures 6, 7). The first type of variability was expressed both as local heterogeneity between closely co-located pores, and as a component of the magnitude of ΔA between more distant pores. Large values of ΔA (>35 μm2) were observed among neighboring pores in both AA and FA. In all cases, local variability in ΔA among neighboring pores was greater than the increase in ΔA with d across the observation area of 1.4 cm (Figures 6, 7).
This study had two objectives, to characterize the distribution of individual pore areas over a range of VPD levels, and to elucidate the role of ambient aerosol in stomatal regulation of plant water relations. The first objective is advanced through analysis of an unusually large number of non-destructively imaged pore areas, acquired previously (Grantz et al., 2018). Previous experiments have approached this objective by imaging many fewer pores, either the same pores repeated over time in epidermal peels (Gorton et al., 1989) or different pores at different times in leaf impressions (Smith et al., 1989). We advance the second objective with a novel analysis of aerosol impacts on stomatal heterogeneity. We showed previously that transpiration per unit stomatal aperture and minimum leaf conductance both increased following aerosol deposition (Burkhardt et al., 2001a; Grantz et al., 2018). No previous analyses have considered whether aerosol deposition might also affect the uniformity of stomatal opening.
Pore Area Distribution
Stomatal opening is characterized by a strong element of randomness that often results in a quasi-normal (bell-shaped) distribution of pore areas (Gorton et al., 1989; Terashima, 1992; Weyers and Lawson, 1997). Our observations confirm earlier reports, including in V. faba (Laisk et al., 1980; Kappen et al., 1987; Spence, 1987), of large within-leaf variability with quasi-normally distributed pore areas, particularly at VPD = 0.33 kPa and 0.97 kPa. There was large heterogeneity even among nearly contiguous pores, as observed previously in V. faba and other species. For example, adjacent stomata of homobaric tobacco (Nicotiana tabacum) ranged from fully closed to fully open (Pospisilova and Santrucek, 1994).
At higher VPD (1.6 kPa and 2.3 kPa), the distributions skewed toward the origin, with a larger proportion of closed pores and reduced mean and median values, as described previously (Grantz, 1990; Grantz et al., 2018). These highly skewed distributions may also belong to theoretically normal distributions that extend into the imaginary territory of negative pore area (Laisk et al., 1980). The maintenance of the bell-shaped distribution even as mean and median pore area declined, reflects synchronous, parallel, and potentially coordinated responses, of similar magnitude by pores of different initial areas and in different locations across the leaf (Saxe, 1979; Spence et al., 1983; Kappen et al., 1987; Spence, 1987). This coherent behavior resembles an emergent property (Mott and Buckley, 2000; Mott and Peak, 2007) but the mechanism of such coordination remains unknown. Elevated VPD increases transpiration and degrades epidermal water status, suggesting metabolic and hydropassive responses of guard cells to epidermal water relations (Mott and Franks, 2001; Buckley, 2005, 2019).
The deposition of hygroscopic aerosol on the leaf surface may act similarly. The presence of thin liquid films lining stomatal pores provides a non-diffusive liquid pathway for water loss (Burkhardt et al., 2001b; Grantz et al., 2018), potentially increasing transpiration and degrading epidermal water status. Hygroscopic, particularly chaotropic, aerosol enhances the formation of such films by deliquescence and facilitates their spread across the leaf surface and into stomatal pores by reducing the surface tension of the liquid on the leaf surface (Monteith, 1957; Eiden et al., 1994; Dutcher et al., 2010; Burkhardt and Hunsche, 2013; Burkhardt and Grantz, 2017; Fernandez et al., 2017). Electrical conductance measurements (Burkhardt and Eiden, 1994; Burkhardt et al., 1999; Burkhardt and Hunsche, 2013) and electron micrography (Grantz et al., 2018) demonstrate these films and their penetration into stomatal pores (Eichert et al., 1998, 2008; Basi et al., 2014; Kaiser, 2014). The resulting hydraulic linkage of apoplast to the leaf boundary layer is associated with a reduction of mean pore area and skewing of areas toward the origin (Burkhardt et al., 2001a, 2012; Burkhardt, 2010; Pariyar et al., 2013; Grantz et al., 2018).
The current study demonstrates an apparent 3-way synergy between elevated VPD, initial stomatal opening at low VPD, and deposition of aerosol, consistent with a cumulative effect on epidermal water relations. At pore areas above about 180 μm2 (at VPD = 0.33 kPa) the dynamics of AA and FA differ as VPD increased. The role of initial opening may reside in the conditions required to sustain liquid water on the leaf surface. This depends on the humidity of the BL which is substantially affected by transpiration. The role of aerosol deposition is to induce water accumulation on the leaf surface from unsaturated air, which becomes more significant at elevated VPD. In the absence of aerosol (FA), condensation on the leaf surface requires approximately 100% RH, but the presence of hygroscopic aerosol (AA) reduces this requirement substantially (e.g., to 75% RH for deliquescence of NaCl aerosol; Burkhardt and Grantz, 2017).
Pore Area Heterogeneity
Vicia faba is homobaric, without bundle sheath extensions or specialized stomatal subsidiary cells. This removes within-areole coordination of gaseous and water potential environments and suggests weaker coordination among stomata at local scale and potentially greater coordination at larger scale than in heterobaric species. In the present study, variability among closely co-located pores was large. This baseline level of ΔA was a dominant component of mean and median ΔA at all levels of pore separation and at all levels of VPD. At larger scale, there was a consistent increase of ΔA with d at all VPD and both aerosol treatments. However, the increase in ΔA with distance was less than the local variability observed between closely spaced stomata. This reflects poor coordination between individual pores that was degraded further with greater separation within the 1.4 cm scale of our observations.
Degraded epidermal water status is associated with increased stomatal heterogeneity (Sharkey and Seemann, 1989; Downton et al., 1990; Haefner et al., 1997). However, in the present study, VPD did not consistently increase the magnitude of ΔA among neighboring pores, increasing only at moderate levels then declining.
The increased water loss per unit stomatal opening observed in the AA treatment (Grantz et al., 2018), suggested that aerosol could increase stomatal heterogeneity. However, this was not observed. Aerosol decreased heterogeneity among neighboring pores at all levels of VPD and reduced the rate of increase of ΔA with d at all VPD.
The liquid film at the leaf surface in the presence of aerosol may link stomata over potentially large areas of the leaf surface by making the near-surface micro-environment more homogeneous. This unifying effect may be more significant than the effect of the liquid film on increased water loss and potential impacts on epidermal water status.
Theories of stomatal optimization (Cowan and Farquhar, 1977) and of cavitation avoidance (Sperry et al., 2017), implicitly consider stomata as a population, whether at the leaf, branch or larger scale (Mott and Peak, 2007). Under conditions of high boundary layer conductance, heterogeneity may be detrimental, but in low wind or with large leaves it may improve gas exchange efficiency. This may explain similar growth of plants grown in AA or FA conditions (unpublished observations) in the slowly ventilated greenhouse, but degraded water use efficiency of AA leaves in the rapidly stirred cuvette of a gas exchange system (Pariyar et al., 2013). By reducing heterogeneity, aerosol deposition could reduce the errors in calculated values of Ci, potentially a factor that distinguishes gas exchange measurements made under controlled environment conditions from those made in more aerosol-rich field conditions.
Stomatal heterogeneity and patchy conductance remain enigmatic (Gunasekera and Berkowitz, 1992). Stomatal patchiness represents stronger stomatal coordination at local scale and poorer coordination at larger scale than predicted by true randomness (Sharkey and Seemann, 1989; Downton et al., 1990; Haefner et al., 1997; Lawson et al., 1998; Mott and Buckley, 2000; Mott and Peak, 2007). While heterogeneity among individual pore areas has been considered as noise (e.g., Cheeseman, 1991) it has more recently been seen as informative regarding the physiological bases of stomatal opening mechanics (Laisk et al., 1980), responses to environment (Mott et al., 1997; Mott and Buckley, 2000) and of the signal processing required for attainment of quasi-optimal stomatal behavior (Cheeseman, 1991; Cardon et al., 1994; Siebke and Weis, 1995; Kaiser and Kappen, 1997; Mott and Peak, 2007). Environmental drivers of coordination among individual members of populations of stomatal pores, including light and VPD, have been considered critical components of the stomatal regulatory system. Aerosol deposition to leaves has not been considered as often, but may play an important role.
We demonstrate, using a very large data set, that stomatal pore areas may be described as quasi-normally distributed, even as mean and median values change with VPD and aerosol deposition. We show that VPD and aerosol deposition both reduce mean pore area and the differences among pores, both at local and greater scale within individual leaves of homobaric V. faba. Synergy of these two ubiquitous environmental factors with the previously unexplored factor of initial stomatal opening, suggests linkage through impacts on epidermal water relations. The aerosol exposure in this experiment was typical of regional European particulate pollution (Putaud et al., 2010), and the study material a short-lived herbaceous annual. These effects may be greater in more hirsute and longer-lived leaves, and in more heavily polluted environments (Kuki et al., 2008; Yang et al., 2015). Aerosol impacts on leaf water relations and gas exchange, mediated by the degree of stomatal opening and its heterogeneity, may be more significant than commonly realized. Further research focused on the control elements regulating stomatal pore area (e.g., Buckley, 2019) and on their interactions with aerosol deposition (e.g., Burkhardt and Grantz, 2017) may be increasingly relevant to characterization of local and global budgets of water and carbon.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
DG and JB conceived of the experiments. DG conducted the gas exchange and pore dimension experiments and analyzed the data. MK developed the software for spatial analysis. All authors contributed to interpretation of the results and wrote the manuscript.
The authors acknowledge funding from the University of California, Riverside Research Allocation Process, the Deutsche Forschungsgemeinschaft (DFG; BU 1099/7-2), and from USDA NIFA Project No. CA-R-BPS-5005-H.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank H. Kaiser for interesting discussions leading up to this analysis as well as considerable assistance with the measurements.
Basi, S., Burkhardt, J., Noga, G., and Hunsche, M. (2014). Hygroscopic salts support the stomatal penetration of glyphosate and influence its biological efficacy. Weed Biol. Manag. 14, 186–197. doi: 10.1111/wbm.12046
Buckley, T. N., and Mott, K. A. (2000). Stomatal responses to non-local changes in PFD: evidence for long distance hydraulic interactions. Plant Cell Environ. 23, 301–309. doi: 10.1046/j.1365-3040.2000.00552.x
Burkhardt, J., Basi, S., Pariyar, S., and Hunsche, M. (2012). Stomatal penetration by aqueous solutions – an update involving leaf surface particles. New Phytol. 196, 774–787. doi: 10.1111/j.1469-8137.2012.04307.x
Burkhardt, J., Kaiser, H., Goldbach, H., and Kappen, L. (1999). Measurements of electrical leaf surface conductance reveal recondensation of transpired water vapour on leaf surfaces. Plant Cell Environ. 22, 189–196. doi: 10.1046/j.1365-3040.1999.00387.x
Burkhardt, J., and Pariyar, S. (2014). Particulate pollutants are capable to ‘degrade’ epicuticular waxes and to decrease the drought tolerance of Scots pine (Pinus sylvestris L.). Environ. Pollut. 184, 659–667. doi: 10.1016/j.envpol.2013.04.041
Burkhardt, J., Zinsmeister, D., Grantz, D. A., Vidic, S., Sutton, M. A., Hunsche, M., et al. (2018). Camouflaged as ‘degraded wax’, hygroscopic aerosols contribute to leaf desiccation, tree mortality, and forest decline. Environ. Res. Lett. 13:085001. doi: 10.1088/1748-9326/aad346
Cardon, Z. G., Mott, K. A., and Berry, J. A. (1994). Dynamics of patchy stomatal movements, and their contribution to steady-state and oscillating stomatal conductance calculated with gas-exchange techniques. Plant Cell Environ. 17, 995–1007. doi: 10.1111/j.1365-3040.1994.tb02033.x
Cheeseman, J. M. (1991). PATCHY: simulating and visualizing the effects of stomatal patchiness on photosynthetic CO2 exchange studies. Plant Cell Environ. 14, 593–599. doi: 10.1111/j.1365-3040.1991.tb01530.x
Downton, W. J. S., Loveys, B. R., and Grant, W. J. R. (1988). Non-uniform stomatal closure induced by water stress causes putative non-stomatal inhibition of photosynthesis. New Phytol. 110, 503–509. doi: 10.1111/j.1469-8137.1988.tb00289.x
Eichert, T., Kurtz, A., Steiner, U., and Goldbach, H. E. (2008). Size exclusion limits and lateral heterogeneity of the stomatal foliar uptake pathway for aqueous solutes and water suspended nanoparticles. Physiol. Plant. 134, 151–160. doi: 10.1111/j.1399-3054.2008.01135.x
Eiden, R., Burkhardt, J., and Burkhardt, O. (1994). Atmospheric aerosol particles and their role in the formation of dew on the surface of plant leaves. J. Aerosol Sci. 25, 367–376. doi: 10.1016/0021-8502(94)90087-6
Fernandez, V., Bahamonde, H. A., Peguero-Pina, J. J., Gil-Pelegrin, E., Sancho-Knapik, D., Gil, L., et al. (2017). Physico-chemical properties of plant cuticles and their functional and ecological significance. J. Exp. Bot. 19, 5293–5306. doi: 10.1093/jxb/erx302
Grantz, D. A., Zinsmeister, D., and Burkhardt, J. (2018). Ambient aerosol increases minimum leaf conductance and alters the aperture–flux relationship as stomata respond to vapor pressure deficit (VPD). New Phytol. 219, 275–286. doi: 10.1111/nph.15102
Hanson, D. T., Stutz, S. S., and Boyer, J. S. (2016). Why small fluxes matter: the case and approaches for improving measurements of photosynthesis and (photo)respiration. J. Exp. Bot. 67, 3027–3039. doi: 10.1093/jxb/erw139
Kaiser, H. (2009). The relation between stomatal aperture and gas exchange under consideration of pore geometry and diffusional resistance in the mesophyll. Plant Cell Environ. 32, 1091–1098. doi: 10.1111/j.1365-3040.2009.01990.x
Kaiser, H., and Kappen, L. (1997). In situ observations of stomatal movements in different light-dark regimes: the influence of endogenous rhythmicity and long-term adjustments. J. Exp. Bot. 48, 1583–1589. doi: 10.1093/jexbot/48.313.1583
Kaiser, H., and Kappen, L. (2000). In situ observation of stomatal movements and gas exchange of Aegopodium podagraria L. In the understorey. J. Exp. Bot. 51, 1741–1749. doi: 10.1093/jexbot/51.351.1741
Kamakura, M., Kosugi, Y., Muramatsu, K., and Muraoka, H. (2012). Simulations and observations of patchy stomatal behavior in leaves of Quercus crispula, a cool-temperate deciduous broad-leaved tree species. J. Plant Res. 125, 339–349. doi: 10.1007/s10265-011-0460-8
Kuki, K. N., Oliva, M. A., Pereir, A. E. G., Costa, A. C., and Cambraia, J. (2008). Effects of simulated deposition of acid mist and iron ore particulate matter on photosynthesis and the generation of oxidative stress in Schinus terebinthifolius Radii and Sophora tomentosa L. Sci. Total Environ. 403, 207–214. doi: 10.1016/j.scitotenv.2008.05.004
Mott, K. A., and Franks, P. J. (2001). The role of epidermal turgor in stomatal interactions following a local perturbation in humidity. Plant Cell Environ. 24, 657–662. doi: 10.1046/j.0016-8025.2001.00705.x
Mott, K. A., Shope, J. C., and Buckley, T. N. (1999). Effects of humidity on light-induced stomatal opening: evidence for hydraulic coupling among stomata. J. Exp. Bot. 50, 1207–1213. doi: 10.1093/jxb/50.336.1207
Nardini, A., Gortan, E., Ramani, M., and Salleo, S. (2008). Heterogeneity of gas exchange rates over the leaf surface in tobacco: an effect of hydraulic architecture? Plant Cell Environ. 31, 804–812. doi: 10.1111/j.1365-3040.2008.01798.x
Osmond, C. B., Kramer, D., and Luttge, U. (1999). Reversible, water stress-induced non-uniform chlorophyll fluorescence quenching in wilting leaves of Potentilla reptans may not be due to patch stomatal responses. Plant Biol. 1, 618–624. doi: 10.1111/j.1438-8677.1999.tb00272.x
Pariyar, S., Eichert, T., Goldbach, H. E., Hunsche, M., and Burkhardt, J. (2013). The exclusion of ambient aerosols changes the water relations of sunflower (Helianthus annuus) and bean (Vicia faba) plants. Environ. Exp. Bot. 88, 43–52. doi: 10.1016/j.envexpbot.2011.12.031
Pringle, K. J., Tost, H., Pozzer, A., Poschl, U., and Lelieveld, J. (2010). Global distribution of the effective aerosol hygroscopicity parameter for CCN activation. Atmos. Chem. Phys. 10, 5241–5255. doi: 10.5194/acp-10-5241-2010
Putaud, J.-P., Van Dingenen, R., Alastuey, A., Bauer, H., Birmili, W., Cyrys, J., et al. (2010). A European aerosol phenomenology 3. Physical and chemical characteristics of particulate matter from 60 rural, urban, and kerbside sites across Europe. Atmos. Environ. 44, 1308–1320. doi: 10.1016/j.atmosenv.2009.12.011
Saxe, H. (1979). A structural and functional study of the coordinated reactions of individual Commelina communis L. stomata (Commelinaceae). Am. J. Bot. 66, 1044–1052. doi: 10.1002/j.1537-2197.1979.tb06320.x
Sharkey, D. T., and Seemann, J. R. (1989). Mild water stress effects on carbon-reduction-cycle intermediates, ribulose bisphosphate carboxylase activity, and spatial homogeneity of photosynthesis in intact leaves. Plant Physiol. 89, 1060–1065. doi: 10.1104/pp.89.4.1060
Spence, R. D. (1987). The problem of variability in stomatal responses, particularly aperture variance, to environmental and experimental conditions. New Phytol. 107, 303–315. doi: 10.1111/j.1469-8137.1987.tb00182.x
Spence, R. D., Sharpe, P. J. H., Powell, R. D., and Rogers, C. A. (1983). Epidermal and guard cell interactions on stomatal aperture in epidermal strips and intact leaves. Ann. Bot. 52, 1–12. doi: 10.1093/oxfordjournals.aob.a086538
Sperry, J. S., Venturas, M. D., Anderegg, W. R. L., Mencuccini, M., Mackay, D. S., Wang, Y., et al. (2017). Predicting stomatal responses to the environment from the optimization of photosynthetic gain and hydraulic cost. Plant Cell Environ. 40, 816–830. doi: 10.1111/pce.12852
Terashima, I., Wong, S.-C., Osmond, C. B., and Farquhar, G. D. (1988). Characterization of non-uniform photosynthesis induced by abscisic acid in leaves having different mesophyll anatomies. Plant Cell Physiol. 29, 385–394.
Tsigaridis, K., Krol, M., Dentener, F. J., Balkanski, Y., Lathiere, J., Metzger, S., et al. (2006). Change in global aerosol composition since preindustrial times. Atmos. Chem. Phys. 6, 5143–5162. doi: 10.5194/acp-6-5143-2006
Keywords: aerosol, climate change, gas exchange, humidity, particulate matter, patchy stomata, vapor pressure deficit
Citation: Grantz DA, Karr M and Burkhardt J (2020) Heterogeneity of Stomatal Pore Area Is Suppressed by Ambient Aerosol in the Homobaric Species, Vicia faba. Front. Plant Sci. 11:897. doi: 10.3389/fpls.2020.00897
Received: 13 December 2019; Accepted: 02 June 2020;
Published: 25 June 2020.
Edited by:Graham Dow, ETH Zürich, Switzerland
Reviewed by:Dimitrios Fanourakis, Technological Educational Institute of Crete, Greece
Andreas M. Savvides, Agricultural Research Institute, Cyprus
Peter Lehmann, ETH Zürich, Switzerland
Copyright © 2020 Grantz, Karr and Burkhardt. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David A. Grantz, email@example.com