ORIGINAL RESEARCH article
The Response of Water Dynamics to Long-Term High Vapor Pressure Deficit Is Mediated by Anatomical Adaptations in Plants
- 1College of Horticulture, Northwest A&F University, Yangling, China
- 2College of Horticulture, Henan Agricultural University, Zhengzhou, China
Vapor pressure deficit (VPD) is the driver of water movement in plants. However, little is known about how anatomical adaptations determine the acclimation of plant water dynamics to elevated VPD, especially at the whole plant level. Here, we examined the responses of transpiration, stomatal conductance (gs), hydraulic partitioning, and anatomical traits in two tomato cultivars (Jinpeng and Zhongza) to long-term high (2.2–2.6 kPa) and low (1.1–1.5 kPa) VPD. Compared to plants growing under low VPD, no variation in gs was found for Jinpeng under high VPD conditions; however, high VPD induced an increase in whole plant hydraulic conductance (Kplant), which was responsible for the maintenance of high transpiration. In contrast, transpiration was not influenced by high VPD in Zhongza, which was primarily attributed to a coordinated decline in gs and Kplant. The changes in gs were closely related to stomatal density and size. Furthermore, high VPD altered hydraulic partitioning among the leaf, stem, and root for both cultivars via adjustments in anatomy. The increase in lumen area of vessels in veins and large roots in Jinpeng under high VPD conditions improved water transport efficiency in the leaf and root, thus resulting in a high Kplant. However, the decreased Kplant for Zhongza under high VPD was the result of a decline of water transport efficiency in the leaf that was caused by a reduction in vein density. Overall, we concluded that the tradeoff in anatomical acclimations among plant tissues results in different water relations in plants under high VPD conditions.
The process of water movement through soil–plant–atmosphere continuum (SPAC) is driven by atmospheric evaporative demand which can be expressed as vapor pressure deficit (VPD). Although the optimal VPD for most greenhouse crops is below 1.5 kPa (Shamshiri et al., 2016), high VPD (>2 kPa) is currently observed in greenhouses, especially during summer (Lu et al., 2015; Zhang et al., 2018). For plants grown under high VPD conditions, a central question is how they regulate transpiration (Carins Murphy et al., 2014; Allen et al., 2015; Grossiord et al., 2017). In a plant, water absorbed through the roots is transported to leaves through the xylem, finally lost via stomata by diffusion. Hence, the regulation of transpiration may occur at whole plant levels. However, the responses of physiological and anatomical traits that could influence transpiration remain largely unknown at whole plant levels during high VPD condition.
Most of the water loss by a plant occurs through stomatal apertures (Macková et al., 2013). Moreover, cuticular pathway is found to be important in regulating water loss when stomatal closure takes place (Fanourakis et al., 2013, 2019). Under steady state conditions and in the vapor phase, the transpiration rate (E) is defined mathematically as a function of stomatal conductance (gs) and VPD. Although the regulation of transpiration depends on the response of gs to VPD during the vapor phase, the efficiency of the hydraulic system determines the amount of liquid water lost to evaporation for any given soil water condition. Using an analogy of Ohm’s law, E can be expressed as the product of hydraulic conductance and water potential gradient in liquid flux. Many studies have proposed a hydraulic feedback loop to interpret the dynamic link between the liquid and gas phases (Buckley, 2005; Domec et al., 2009; Simonin et al., 2015). Thus, a coordination may exist between gs and whole plant hydraulic conductance (Kplant) with respect to water transport across the soil–plant–atmosphere continuum.
To deal with long-term environmental fluctuation, plants have evolved high plasticity in carbon allocation (Freschet et al., 2018). The changes in carbon investment generally trigger adjustments in anatomical traits involved in plant water dynamics, which occur at multiple places in the plant including stomatal and xylem tissues (Sperry et al., 1998; de Boer et al., 2016; Dewar et al., 2018). Previous studies demonstrates stomatal density decreased to prevent excessive water loss under high VPD in tomato (Lu et al., 2015; Du et al., 2019), rose (Fanourakis et al., 2012), and fava bean (Aliniaeifard et al., 2014). However, few stomata mean a reduction in the maximum potential carbon acquisition. Alternatively, plants resort to an efficient hydraulic system to withstand excessive evaporative demand. Despite a high carbon investment to xylem, these adjustments contribute to maintain carbon acquirement. Thus, a tradeoff between water loss and carbon acquirement would exist during the acclimation process of plants to high VPD. Additionally, the hydraulic system in plants shows a strong hydraulic segmentation (Cruiziat et al., 2002; Sperry and Love, 2015). Although the dynamics of hydraulic conductance of leaves, stems, and roots have been well-documented in response to soil water deficit (Domec et al., 2009; Torres-Ruiz et al., 2015), the adjustment of hydraulic structure to long-term high VPD has been poorly investigated. Therefore, systematic knowledge about acclimation at the whole plant level is critical to determining the responses of plants to high VPD and is necessary to understand the tradeoff between carbon investment in regulating water dynamics and carbon acquirement.
Tomato (Solanum lycopersicum L.) is one of the most important agricultural plants in the world. High VPD induces contrasting responses in plant water status among tomato cultivars (Zhang et al., 2018; Du et al., 2019). In the present study, two tomato cultivars, Jingpeng and Zhongza, were selected on the basis that they exhibit different responses to altered VPD (Du et al., 2019). The responses of water dynamics and anatomical traits were measured on these two cultivars after exposure to long-term high and low VPD. We hypothesized that (1) for plants with high water loss under high VPD, Kplant would increase with unaffected gs; (2) for plants with relatively low water loss under high VPD, Kplant and gs would synchronously decline; and (3) acclimation in terms of water dynamics is related to anatomical changes at multiple places in plants.
Materials and Methods
Plant Material and Growth Conditions
Seeds of Jinpeng and Zhongza were germinated and grown in plastic pots [15 cm × 10 cm (diameter × height); 1 plant/pot] containing a mixed peat-perlite substrate. The seedlings were kept in a walk-in growth chamber. The light in the chambers was given daily for 14 h at a photon flux density of 300 μmol m–2 s–1. The temperature was 28–30°C day/19–20°C night. Relative humidity was regulated between 64–70% day/77–82% night using an ultrasonic humidifier (KAJ-9.0B, Kawasima Appliance Co., Ltd., Changzhou, China) and dehumidifier (DH-702B, Chuanjing Electric Co., Ltd., Hangzhou, China). Consequently, the VPD was 1.1–1.5 kPa day/0.4–0.5 kPa night. After 5 weeks, 30 of the healthiest plants were divided into two random groups of 15 for each cultivar. For low VPD treatment, the plants were kept on previous humidity conditions. A high VPD was performed by setting 36–42% relative humidity during the day (VPD 2.2–2.6 kPa). Plants were grown for 30 d and kept well-watered during the entire growth period. New fully expanded leaves were used for measurements.
Transpiration and Stomatal Conductance
To determine transpiration at the canopy level (Ecanopy) during the photoperiod, pots were covered with plastic film and aluminum foil on the day before measurement. After at least 2 h of acclimation in the photoperiod, five pots were weighed 2 and 8 h after turning on the lights. Ecanopy was calculated by differences in weight divided by the total leaf area. Canopy stomatal conductance (gs-canopy) was determined according to a simplified inversion of the Penman–Monteith model (Monteith and Unsworth, 1990):
where R is the universal gas constant adjusted for water vapor (0.46 m3 kPa K–1 kg–1), Ta is air temperature (K), and ρ is the density of water (998 kg m–3).
Leaf level transpiration (Eleaf) and stomatal conductance (gs–leaf) were measured with a plant porometer (Yaxin-1301, Yaxin Liyi Technology Co., Ltd., Beijing, China). After at least 2 h of acclimation in the photoperiod, five leaves from different plants for each treatment and cultivar were used for measurements. VPD, light, and the temperature of the cuvette were kept at ambient levels.
Hydraulic conductance of the leaf (Kleaf), stem (Kstem), and root (Kroot) and Kplant were calculated according to Domec et al. (2009). Briefly,
where Ψsoil is the soil water potential measured with a PSYPRO Water Potential System (PSYPRO; Wescor, Inc., Logan, UT, United States), Ψleaf is the water potential of the leaf used for gas exchange measurement (transpiring leaf), Ψstem-up is the stem water potential in the upper crown section, and Ψstem–base is the stem water potential at the stem base. Ψstem was estimated from the leaf water potential of a non-transpiring leaf (achieved by covering with plastic film and aluminum foil the night before measurement) (Richter, 1997; Zsögön et al., 2015). The leaf adjacent to the transpiring leaf was used for Ψstem-up and the first true leaf was used for Ψstem-base. Measurements of Ψleaf, Ψstem-up, and Ψstem-base were performed on the same five plants used for gas exchange measurements with a pressure chamber (PMS, Corvallis, OR, United States). After measurement of hydraulic conductance, the plants were used for morphological observation.
Stomatal morphological characteristics were determined using the impression method as described by Xu and Zhou (2008). Briefly, epidermis was smeared with nail varnish in the mid-area between the midrib and lateral margin, avoiding midrib and secondary veins (Fanourakis et al., 2015b). Then, the thin film (approximately 5 mm × 5 mm) was peeled off from the leaf surface and mounted on a glass slide. All stomatal characteristics were measured on both the adaxial and abaxial sides of the leaf. For determining stomatal density (SD), five images per sampling area were taken at a magnification of 200× with a light microscope (BX50, Olympus, Tokyo, Japan). Stomatal area (SA) was measured on at least 20 stomata per sampling area at a magnification of 400× with ImageJ software (National Institutes of Health, Bethesda, MD, United States). SA was defined as the combined area of pore and a pair of guard cells following Savvides et al. (2012).
Vein and Stem Anatomical Traits
To evaluate leaf vein traits, leaflets were detached from the leaves used for Ψleaf measurements. Cleared leaflets were scanned to estimate the length of the midrib and secondary vein as well as the leaf area. Veins with an order higher than secondary were measured on 1-cm2 leaf pieces from the center of each leaf at a magnification of 40× according to Kono et al. (1982). Briefly, the lamina was boiled in 70% ethanol to remove pigment. After washing with distilled water, samples were transferred to boiling 85% lactic acid for 20 min and then spread out flat on a slide for observation. The vein density was calculated as the ratio of total vein length to the analyzed area.
To assess the xylem composition in veins, segments (0.3 cm in length) were cut from the petiole immediately below the lamina insertion point. This section was selected because it was connected to the midrib and water entering the leaflet would have to pass through this part. The middle of the plant was sampled to estimated xylem composition in the stem. After fixing in a mixture of formaldehyde, acetic acid, and alcohol for 24 h, the sample material was dehydrated in a graded ethanol-xylene series and infiltrated with paraffin. Then 12-μm thick sections were made using a rotary microtome, stained with safranin and fast green, and mounted on slides with a cover slip (Berlyn and Miksche, 1976). The cross-sectional area of the petiole and total number of vessels in the petiole were measured at magnifications of 40× and 100×, respectively. Vessel density was defined as the number of vessels per unit area. The lumen area of vessels in the leaf vein (Alumen-leaf) and stem (Alumen-stem), and wall thickness of vessels in the stem (Tw-stem), were estimated at a magnification of 400× with ImageJ (at least three different field-of-view regions). The lumen diameter of each vessel was calculated from its lumen area, assuming a circular shape. The maximum theoretical leaf vein axial hydraulic conductivity (Kleaf-max) was estimated as follows (North et al., 2013):
where N is total number of vessels in the petiole, di is the lumen diameter of each vessel, and η is the viscosity of water, further normalized by leaf area (Sack and Frole, 2006).
Root Morphological Characteristics
After measurement of hydraulic conductance, roots were detached from the plant and carefully washed. The cleaned samples were scanned and analyzed with WinRhizo software (WinRhizo, Regent Ltd., Canada).
Five plants per treatment were selected to measure plant biomass and total leaf area. The aboveground and underground dry biomass of plant was determined after drying to a constant weight at 80°C. The net assimilation rate (NAR) was calculated as follows:
where W1 and W2 are total dry weights of the whole plant at times T1 and T2; L1 and L2 are total leaf areas of the whole plant at times T1 and T2.
All statistical analyses were performed using SPSS 19.0 (IBM Corp., Armonk, NY, United States). One-way ANOVA was used to test differences between mean values (Duncan’s test). Two-way ANOVA was used to determine the main effects of cultivar, treatment, and their interactions.
Transpiration and Stomatal Conductance
The effects of VPD on Ecanopy and Eleaf were similar and cultivar specific (Figure 1). The Ecanopy and Eleaf of Jinpeng significantly increased under high VPD compared to low VPD. On the contrary, in Zhongza, these measures were both unaffected by VPD treatment. Under high VPD, gs-canopy and gs-leaf decreased by 29 and 35%, respectively, in Zhongza relative to low VPD whereas VPD had no influence on either gs-canopy or gs-leaf in Jinpeng.
Figure 1. Transpiration rate at the canopy (Ecanopy) (A) and leaf level (Eleaf) (B), stomatal conductance at the canopy (gs-canopy) (C) and leaf level (gs-leaf) (D) for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD. Data are means ± standard error (n = 5 plants). Different letters denote statistically significant differences (Duncan’s test, P < 0.05). Two-way ANOVA was used to estimate the effect of cultivar (C), treatment (T), and their interaction (C × T) (∗∗P < 0.01; ∗P < 0.05; NS, not significant).
No significant differences in SD or SA for either the adaxial or abaxial sides of the leaf were found between low and high VPD in Jinpeng (Table 1). For Zhongza, high VPD significantly decreased SD on both sides and SA on the adaxial side only. Moreover, the integrated values of SD and SA for the whole leaf also declined under high VPD compared to low VPD in Zhongza.
Table 1. Stomatal density (SD) and area (SA) for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD.
Plant Hydraulic Conductance
Compared to low VPD, Kplant and Kleaf in Jinpeng significantly increased but in Zhongza declined by 24 and 36%, respectively, under high VPD (Figure 2). No difference in Kstem was noted between low and high VPD conditions for either cultivar. Kroot increased by 46% in Jingpeng under high VPD but was similar under low and high VPD in Zhongza.
Figure 2. Hydraulic conductance of whole plant (Kplant) (A), leaf (Kleaf) (B), stem (Kstem) (C), and root (Kroot) (D) for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD. Data are means ± standard error (n = 5 plants). Different letters denote statistically significant differences (Duncan’s test, P < 0.05). Two-way ANOVA was used to estimate the effect of cultivar (C), treatment (T), and their interaction (C × T) (**P < 0.01; *P < 0.05; NS, not significant).
Relative contributions were analyzed to discern the role of Kleaf, Kstem, and Kroot in the observed changes in Kplant (Figure 3). The percentages of Kleaf, Kstem, and Kroot that made up Kplant were similar between low and high VPD in both Jinpeng and Zhongza. Kleaf accounted for 56–66% of the changes in Kplant, followed by Kstem (26%) and Kroot (13%).
Figure 3. Percentage contributions of leaf (Kleaf), stem (Kstem), and root (Kroot) hydraulic conductance to whole plant hydraulic conductance for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD. Different letters denote statistically significant differences (Duncan’s test, P < 0.05).
Leaf Vein Traits
Vein density was not affected by VPD treatment in Jinpeng, but declined in Zhongza under high VPD compared to low VPD (Figure 4). For both cultivars, high VPD had no influence on vessel density in leaf veins. Under high VPD, Alumen-leaf and Kleaf-max increased by 28 and 57%, respectively, in Jinpeng but decreased by 20 and 37%, respectively, in Zhongza.
Figure 4. Vein density (A), vessel density (B), lumen area of vessels in leaf vein (Alumen-leaf) (C), and maximum theoretical leaf vein axial hydraulic conductivity (Kleaf-max) (D) for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD. Data are means ± standard error (n = 5 plants). Different letters denote statistically significant differences (Duncan’s test, P < 0.05). Two-way ANOVA was used to estimate the effect of cultivar (C), treatment (T), and their interaction (C × T) (∗∗P < 0.01; ∗P < 0.05; NS, not significant).
Stem and Root Morphological Characteristics
Alumen-stem and Tw-stem was not significantly different between low and high VPD in either Jinpeng or Zhongza (Figure 5). Under high VPD, Jinpeng roots had 57, 33, and 17% higher volume, surface area, and average diameter, respectively, than under low VPD (Table 2). No differences in root volume, surface area, or average diameter were found in Zhongza between low and high VPD. Root total length was similar under low and high VPD in both cultivars.
Figure 5. Lumen area (Alumen-stem) (A) and wall thickness (Tw-stem) (B) of vessels in stem for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD. Data are means ± standard error (n = 5 plants). Different letters denote statistically significant differences (Duncan’s test, P < 0.05). Two-way ANOVA was used to estimate the effect of cultivar (C), treatment (T), and their interaction (C × T) (∗P < 0.05; NS, not significant).
Table 2. Root volume, surface area, average diameter, and total length for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD.
Aboveground dry weight, total dry weight and NAR was similar under high and low VPD in Jinpeng, but these parameters were significantly lower under high VPD than low VPD in Zhongza (Table 3). Although no statistical differences in underground dry weight were found between low and high VPD for both cultivars, it increased by 14% in Jinpeng and decreased by 7% in Zhongza under high VPD.
Table 3. Aboveground dry weight, underground dry weight, total dry weight, and net assimilation rate (NAR) for two tomato cultivars, Jinpeng and Zhongza, grown under low (1.1–1.5 kPa) and high (2.2–2.6 kPa) VPD.
In the soil–plant–atmosphere continuum, VPD is the driving force for water flow. High VPD induced a higher E in Jinpeng than low VPD, but E was not affected by VPD treatment in Zhongza, suggesting cultivar differences in water dynamic responses to long-term high VPD (Figure 1). Acclimation of water transport pathways for liquid and vapor phases to high VPD may separately or simultaneously occur in plants, leading to a new homeostatic state (Brodribb and Holbrook, 2006; Domec et al., 2009; Simonin et al., 2015; Fernandes-Silva et al., 2016). In the vapor phase, evaporation from the mesophyll surface to the substomatal cavity is enhanced with increasing VPD. The reduced gs in Zhongza under high VPD prevents excessive vapor diffusion out of the leaf while E is proportional to changes in VPD due to unchanged gs in Jinpeng. Additionally, water supply determines how much water can evaporate from the plant. For liquid flow moving through the plant, water supply depends on the Kplant for given soil moisture conditions. In Jinpeng, the increased E under high VPD would been maintained by a high Kplant (Supplementary Figure S1). The coupling between Kplant and E enables leaves to minimize variation in plant water status (Simonin et al., 2015). However, the coordination between Kplant and E is disrupted in Zhongza. The reduced Kplant suggests a limited water supply in Zhongza under high VPD. Thus, E for Zhongza growing under high VPD is kept at the same level as for low VPD conditions to maintain the balance between water supply and loss.
gs is tightly linked to plant water status (Bunce, 2006; Ripullone et al., 2007). According to the hydraulic feedback hypothesis, a hydraulic feedback loop can describe the mechanism driving the relationship between gs and Kplant (Buckley, 2005; Brodribb and Jordan, 2011; Savvides et al., 2012; Simonin et al., 2015). The absence of change in gs and the increase in Kplant demonstrated by Jinpeng growing under high VPD suggest that adequate water supply prevents leaf dehydration and therefore the stomata remain open. This acclimation would maximize carbon acquisition under high VPD, which was confirmed by the similar plant biomass and NAR under high and low VPD in Jinpeng (Table 3). In contrast, a coordinated declines in Kplant and gs were noted in Zhongza under high VPD (Supplementary Figure S2). Low water transport efficiency has been regarded as the initial cause of water stress symptoms (Fanourakis et al., 2012). Plants are prone to close stomata due to reductions in leaf turgor when subjected to water stress (Martins et al., 2016; Rodriguez-Dominguez et al., 2016). Thus, gs in Zhongza decreases under high VPD to minimize water loss in plants despite limiting carbon acquisition (Table 3). However, the coordination between gs and Kplant would also play a role in maintaining the integrity of xylem water transport and reducing the risk of hydraulic failure (Galmés et al., 2013; Liu et al., 2015; Salmon et al., 2015; Du et al., 2018).
Alternatively, the responses of gs to long-term high VPD could be driven by changes in stomatal size and density (Fanourakis et al., 2013; Lu et al., 2015). Franks and Beerling (2009) showed mathematically that gs is positively related to stomatal density and size based on the physics of diffusion through pores. Stomatal size and density decreased in Zhongza under high VPD but no active acclimation to high VPD occurred in Jinpeng. Thus, the reductions in stomatal density and size under high VPD are at least partially responsible for the decline in gs (Table 1 and Supplementary Figure S2). This adjustment in stomatal morphology helps plants avoid the risk of excessive water loss under conditions of high evaporative demand. Furthermore, morphological differences in stomata were more evident in the abaxial epidermis for Zhongza. Fanourakis et al. (2015a) also confirmed that operating gs was mostly situated on the abaxial surface. For plants, the regulation of gs is more effective by changing abaxial stomatal morphology because a broader distribution of stomata in abaxial epidermis.
Due to hydraulic segmentation in the plant water transport system, the hydraulic resistance of the whole plant is partitioned into its functional components related to leaves, stems, and roots (Cruiziat et al., 2002; Sperry and Love, 2015). Our study showed that Kleaf determines approximately 60% of the changes in Kplant (Figure 3). Moreover, a synchronized increase in Kplant and Kleaf was found for Jinpeng whereas a coordinated decrease in Kplant and Kleaf occurred for Zhongza under high VPD (Figure 2). Thus, Kplant was mostly dominated by Kleaf, which is consistent with previous studies even though leaves constitute less than 3% of the pathway for water flow through the whole plant (Sack et al., 2003; Domec et al., 2009; Tabassum et al., 2016).
During water movement in the leaf, leaf vein structural characteristics are a substantial constraint on Kleaf (Sack and Holbrook, 2006; Carins Murphy et al., 2014; Xiong et al., 2017). The decline in vein density of Zhongza under high VPD means not only a reduction in the surface area for the exchange of xylem water with surrounding tissue but an increase in distance for water movement from the xylem into mesophyll cells (Roth-Nebelsick et al., 2001; Sack and Frole, 2006), therefore resulting in a decrease in Kleaf. Additionally, Alumen-leaf also declined in Zhongza under high VPD (Figure 4 and Supplementary Figure S3). A small vessel diameter in general corresponds to a high resistance for water transport, but it can withstand very low negative pressure without generating an embolism (Pittermann and Sperry, 2003). In contrast, Jinpeng showed a unaffected leaf vein density and synchronized increases in Alumen-leaf and E under high VPD (Figure 4 and Supplementary Figure S3), suggesting that large vessel diameter plays a critical role in improving Kleaf and maintaining E.
Interestingly, high VPD treatment also resulted in an increase in Kroot and a concurrent increase in root volume, surface area, and average diameter for Jinpeng whereas no acclimations of Kroot and root morphology were observed in Zhongza (Table 2 and Supplementary Figure S4). Meanwhile, the large root was tightly related to high E. Modifications in root morphology would improve hydraulic properties at the soil-root interface, promoting water uptake and transport (Steudle, 2000; Domec et al., 2009). However, this active acclimation in root and leaf veins is likely to occur in plants that can maintain a high carbon acquirement because of a high carbon investment. Plants with closed stomata under high VPD would prefer to reduce the investment to aboveground according to the multiple limitation hypothesis (Farrior et al., 2013; Sellin et al., 2015). This may explain the different changes in plant biomass and NAR between Jinpeng and Zhongza under high VPD (Table 3). In the present study, no changes in Kstem, Alumen-stem and Tw-stem were found under high VPD in either cultivar (Figure 5). Thus, the major role in regulating Kplant under high VPD is attributed to the leaf and root, as previously reported for several tree species (Domec et al., 2009; Torres-Ruiz et al., 2015).
The present study indicates that different hydraulic regulation strategies are responsible for the discrepancies found in terms of water dynamics in the cultivars studied. Furthermore, the main results of the present study reinforce the idea that the responses of water dynamics to high VPD depend on acclimation at the whole plant level, particularly at the stomatal, leaf and root levels. High carbon investment in the hydraulic architecture of leaves and roots warrants a sufficient water supply to meet the requirement of water loss, maintaining gs. In contrast, plants are prone to triggering a simultaneous decrease in gs and Kleaf to prevent excessive water loss when environmental factors have a negative effect on the leaf hydraulic system. Therefore, high VPD treatment would impose a tradeoff between the water and carbon economy of the plant.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
QD, XJ, and JL conceived and designed the experiments. QD, XS, JZ, PB, and JD performed the experiments. QD and XJ analyzed the data and wrote this manuscript. JL contributed extensively to its finalization.
This work was financially supported by the China Agriculture Research System (CARS-23-C05), Project of Science and Technology Plan in Qinghai (2018-NK-123), and Project of Scientific and Technical Innovation Demonstration for Vegetable Industry in Qinghai (QNKP-2018-03-01).
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 Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00758/full#supplementary-material
Aliniaeifard, S., Matamoros, P. M., and van Meeteren, U. (2014). Stomatal malfunctioning under low VPD conditions: induced by alterations in stomatal morphology and leaf anatomy or in the ABA signaling? Physiol. Plant. 152, 688–699. doi: 10.1111/ppl.12216
Allen, C. D., Breshears, D. D., and McDowell, N. G. (2015). On underestimation of global vulnerability to tree mortality and forest die-off from hotter drought in the Anthropocene. Ecosphere 6, 129. doi: 10.1890/es15-00203.1
Brodribb, T. J., and Holbrook, N. M. (2006). Declining hydraulic efficiency as transpiring leaves desiccate: two types of response. Plant Cell Environ. 29, 2205–2215. doi: 10.1111/j.1365-3040.2006.01594.x
Brodribb, T. J., and Jordan, G. J. (2011). Water supply and demand remain balanced during leaf acclimation of Nothofagus cunninghamii trees. New Phytol. 192, 437–448. doi: 10.1111/j.1469-8137.2011.03795.x
Carins Murphy, M. R., Jordan, G. J., and Brodribb, T. J. (2014). Acclimation to humidity modifies the link between leaf size and the density of veins and stomata. Plant Cell Environ. 37, 124–131. doi: 10.1111/pce.12136
de Boer, H. J., Price, C. A., Wagner-Cremer, F., Dekker, S. C., Franks, P. J., and Veneklaas, E. J. (2016). Optimal allocation of leaf epidermal area for gas exchange. New Phytol. 210, 1219–1228. doi: 10.1111/nph.13929
Dewar, R., Mauranen, A., Makela, A., Holtta, T., Medlyn, B., and Vesala, T. (2018). New insights into the covariation of stomatal, mesophyll and hydraulic conductances from optimization models incorporating nonstomatal limitations to photosynthesis. New Phytol. 217, 571–585. doi: 10.1111/nph.14848
Domec, J. C., Noormets, A., King, J. S., Sun, G., McNulty, S. G., Gavazzi, M. J., et al. (2009). Decoupling the influence of leaf and root hydraulic conductances on stomatal conductance and its sensitivity to vapour pressure deficit as soil dries in a drained loblolly pine plantation. Plant Cell Environ. 32, 980–991. doi: 10.1111/j.1365-3040.2009.01981.x
Du, Q. J., Liu, T., Jiao, X. C., Song, X. M., Zhang, J. Y., and Li, J. M. (2019). Leaf anatomical adaptations have central roles in photosynthetic acclimation to humidity. J. Exp. Bot. 70, 4949–4962. doi: 10.1093/jxb/erz238
Du, Q. J., Xing, G. M., Jiao, X. C., Song, X. M., and Li, J. M. (2018). Stomatal responses to long-term high vapor pressure deficits mediated most limitation of photosynthesis in tomatoes. Acta Physiol. Plant. 40, 149. doi: 10.1007/s11738-018-2723-7
Fanourakis, D., Carvalho, S. M. P., Almeida, D. P. F., van Kooten, O., van Doorn, W. G., and Heuvelink, E. (2012). Postharvest water relations in cut rose cultivars with contrasting sensitivity to high relative air humidity during growth. Postharvest Biol. Technol. 64, 64–73. doi: 10.1016/j.postharvbio.2011.09.016
Fanourakis, D., Heuvelink, E., and Carvalho, S. M. P. (2013). A comprehensive analysis of the physiological and anatomical components involved in higher water loss rates after leaf development at high humidity. J. Plant Physiol. 170, 890–898. doi: 10.1016/j.jplph.2013.01.013
Fanourakis, D., Giday, H., Milla, R., Pieruschka, R., Kjaer, K. H., Bolger, M., et al. (2015a). Pore size regulates operating stomatal conductance, while stomatal densities drive the partitioning of conductance between leaf sides. Ann. Bot. 115, 555–565. doi: 10.1093/aob/mcu247
Fanourakis, D., Heuvelink, E., and Carvalho, S. M. P. (2015b). Spatial heterogeneity in stomatal features during leaf elongation: an analysis using Rosa hybrida. Funct. Plant Biol. 42, 737–745. doi: 10.1071/FP15008
Fanourakis, D., Hyldgaard, B., Giday, H., Aulik, I., Bouranis, D., Körner, O., et al. (2019). Stomatal anatomy and closing ability is affected by supplementary light intensity in rose (Rosa hybrida L.). Hortic. Sci. 46, 81–89. doi: 10.17221/144/2017-HORTSCI
Farrior, C. E., Tilman, D., Dybzinski, R., Reich, P. B., Levin, S. A., and Pacala, S. W. (2013). Resource limitation in a competitive context determines complex plant responses to experimental resource additions. Ecology 94, 2505–2517. doi: 10.1890/12-1548.1
Fernandes-Silva, A. A., López-Bernal, Á, Ferreira, T. C., and Villalobos, F. J. (2016). Leaf water relations and gas exchange response to water deficit of olive (cv. Cobrançosa) in field grown conditions in Portugal. Plant Soil 402, 191–209. doi: 10.1007/s11104-015-2786-9
Franks, P. J., and Beerling, D. J. (2009). Maximum leaf conductance driven by CO2 effects on stomatal size and density over geologic time. Proc. Natl. Acad. Sci. U.S.A. 106, 10343–10347. doi: 10.1073/pnas.0904209106
Freschet, G. T., Violle, C., Bourget, M. Y., Scherer-Lorenzen, M., and Fort, F. (2018). Allocation, morphology, physiology, architecture: the multiple facets of plant above- and below-ground responses to resource stress. New Phytol. 219, 1338–1352. doi: 10.1111/nph.15225
Galmés, J., Ochogavía, J. M., Gago, J., Roldán, E. J., Cifre, J., and Conesa, M. À (2013). Leaf responses to drought stress in Mediterranean accessions of Solanum lycopersicum: anatomical adaptations in relation to gas exchange parameters. Plant Cell Environ. 36, 920–935. doi: 10.1111/pce.12022
Grossiord, C., Sevanto, S., Borrego, I., Chan, A. M., Collins, A. D., Dickman, L. T., et al. (2017). Tree water dynamics in a drying and warming world. Plant Cell Environ. 40, 1861–1873. doi: 10.1111/pce.12991
Kono, Y., Nakata, K., and Tatsumi, J. (1982). Observations of cross veins of the second foliage leaf blade in the rice plant by use of a revised method for clearing leaves. Jap. J. Crop Sci. 51, 445–454. doi: 10.1626/jcs.51.445
Liu, Y. Y., Song, J., Wang, M., Li, N., Niu, C. Y., and Hao, G. Y. (2015). Coordination of xylem hydraulics and stomatal regulation in keeping the integrity of xylem water transport in shoots of two compound-leaved tree species. Tree Physiol. 35, 1333–1342. doi: 10.1093/treephys/tpv061
Lu, N., Nukaya, T., Kamimura, T., Zhang, D., Kurimoto, I., Takagaki, M., et al. (2015). Control of vapor pressure deficit (VPD) in greenhouse enhanced tomato growth and productivity during the winter season. Sci. Hortic. 197, 17–23. doi: 10.1016/j.scienta.2015.11.001
Macková, J., Vašková, M., Macek, P., Hronková, M., Schreiber, L., and Šantrůček, J. (2013). Plant response to drought stress simulated by ABA application: changes in chemical composition of cuticular waxes. Environ. Exp. Bot. 86, 70–75. doi: 10.1016/j.envexpbot.2010.06.005
Martins, S. C. V., McAdam, S. A. M., Deans, R. M., DaMatta, F. M., and Brodribb, T. J. (2016). Stomatal dynamics are limited by leaf hydraulics in ferns and conifers: results from simultaneous measurements of liquid and vapour fluxes in leaves. Plant Cell Environ. 39, 694–705. doi: 10.1111/pce.12668
North, G. B., Lynch, F. H., Maharaj, F. D. R., Phillips, C. A., and Woodside, W. T. (2013). Leaf hydraulic conductance for a tank bromeliad: axial and radial pathways for moving and conserving water. Front. Plant Sci. 4:78. doi: 10.3389/fpls.2013.00078
Ripullone, F., Guerrieri, M. R., Nole, A., Magnani, F., and Borghetti, M. (2007). Stomatal conductance and leaf water potential responses to hydraulic conductance variation in Pinus pinaster seedlings. Trees 21, 371–378. doi: 10.1007/s00468-007-0130-6
Rodriguez-Dominguez, C. M., Buckley, T. N., Egea, G., de Cires, A., Hernandez-Santana, V., Martorell, S., et al. (2016). Most stomatal closure in woody species under moderate drought can be explained by stomatal responses to leaf turgor. Plant Cell Environ. 39, 2014–2026. doi: 10.1111/pce.12774
Sack, L., Cowan, P. D., Jaikumar, N., and Holbrook, N. M. (2003). The ‘hydrology’ of leaves: co-ordination of structure and function in temperate woody species. Plant Cell Environ. 26, 1343–1356. doi: 10.1046/j.0016-8025.2003.01058.x
Salmon, Y., Torres-Ruiz, J. M., Poyatos, R., Martinez-Vilalta, J., Meir, P., Cochard, H., et al. (2015). Balancing the risks of hydraulic failure and carbon starvation: a twig scale analysis in declining Scots pine. Plant Cell Environ. 38, 2575–2588. doi: 10.1111/pce.12572
Savvides, A., Fanourakis, D., and van Ieperen, W. (2012). Co-ordination of hydraulic and stomatal conductances across light qualities in cucumber leaves. J. Exp. Bot. 63, 1135–1143. doi: 10.1093/jxb/err348
Sellin, A., Rosenvald, K., Õunapuu-Pikas, E., Tullus, A., Ostonen, I., and Lõhmus, K. (2015). Elevated air humidity affects hydraulic traits and tree size but not biomass allocation in young silver birches (Betula pendula). Fronti. Plant Sci. 6:860. doi: 10.3389/fpls.2015.00860
Shamshiri, R., Che Man, H., Zakaria, A., Van Beveren, P., Wan Ismail, W. I., and Ahmad, D. (2016). Membership function model for defining optimality of vapor pressure deficit in closed-field cultivation of tomato. Acta Hortic. 1152, 281–290. doi: 10.17660/actahortic.2017.1152.38
Simonin, K. A., Burns, E., Choat, B., Barbour, M. M., Dawson, T. E., and Franks, P. J. (2015). Increasing leaf hydraulic conductance with transpiration rate minimizes the water potential drawdown from stem to leaf. J. Exp. Bot. 66, 1303–1315. doi: 10.1093/jxb/eru481
Sperry, J. S., Adler, F. R., Campbell, G. S., and Comstock, J. P. (1998). Limitation of plant water use by rhizosphere and xylem conductance: results from a model. Plant Cell Environ. 21, 347–359. doi: 10.1046/j.1365-3040.1998.00287.x
Tabassum, M. A., Zhu, G., Hafeez, A., Wahid, M. A., Shaban, M., and Li, Y. (2016). Influence of leaf vein density and thickness on hydraulic conductance and photosynthesis in rice (Oryza sativa L.) during water stress. Sci. Rep. 6:36894. doi: 10.1038/srep36894
Torres-Ruiz, J. M., Diaz-Espejo, A., Perez-Martin, A., and Hernandez-Santana, V. (2015). Role of hydraulic and chemical signals in leaves, stems and roots in the stomatal behaviour of olive trees under water stress and recovery conditions. Tree Physiol. 35, 415–424. doi: 10.1093/treephys/tpu055
Xiong, D. L., Flexas, J., Yu, T. T., Peng, S. B., and Huang, J. L. (2017). Leaf anatomy mediates coordination of leaf hydraulic conductance and mesophyll conductance to CO2 in Oryza. New Phytol. 213, 572–583. doi: 10.1111/nph.14186
Zhang, D. L., Jiao, X. C., Du, Q. J., Song, X. M., and Li, J. M. (2018). Reducing the excessive evaporative demand improved photosynthesis capacity at low costs of irrigation via regulating water driving force and moderating plant water stress of two tomato cultivars. Agric. Water Manag. 199, 22–33. doi: 10.1016/j.agwat.2017.11.014
Zsögön, A., Alves Negrini, A. C., Peres, L. E. P., Nguyen, H. T., and Ball, M. C. (2015). A mutation that eliminates bundle sheath extensions reduces leaf hydraulic conductance, stomatal conductance and assimilation rates in tomato (Solanum lycopersicum). New Phytol. 205, 618–626. doi: 10.1111/nph.13084
Keywords: anatomical acclimations, hydraulics, stomatal conductance, transpiration, vapor pressure deficit
Citation: Du Q, Jiao X, Song X, Zhang J, Bai P, Ding J and Li J (2020) The Response of Water Dynamics to Long-Term High Vapor Pressure Deficit Is Mediated by Anatomical Adaptations in Plants. Front. Plant Sci. 11:758. doi: 10.3389/fpls.2020.00758
Received: 14 January 2020; Accepted: 12 May 2020;
Published: 05 June 2020.
Edited by:Juan Pedro Ferrio, Fundacion Agencia Aragonesa para la Investigacion y el Desarrollo, Spain
Reviewed by:Dimitrios Fanourakis, Technological Educational Institute of Crete, Greece
Gerardo Tapia, Institute of Agricultural Research, Chile
Copyright © 2020 Du, Jiao, Song, Zhang, Bai, Ding and Li. 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: Jianming Li, firstname.lastname@example.org
†These authors have contributed equally to this work