Original Research ARTICLE
Transpiration Response of Cotton to Vapor Pressure Deficit and Its Relationship With Stomatal Traits
- 1Adaptive Cropping Systems Laboratory, Beltsville Agricultural Research Center, United States Department of Agriculture-Agricultural Research Service-Northeast Area, Beltsville, MD, United States
- 2Oak Ridge Institute for Science and Education, Oak Ridge, TN, United States
Many studies have demonstrated that the cotton in warm environments is vulnerable to water-limitations thus reducing the yield. A number of plant traits have been recommended to ameliorate the effects of water deficits on plant growth and yield. Limitation on maximum transpiration rate (TR) under high vapor pressure deficit (VPD), usually occurs during midday, is often considered as a water conservation trait. The genotypes with this trait are desirable in high VPD environments where water deficits commonly develop in the later part of the growing season. Our objective of the study was to find the genotypic variation for the trait limited TR under high VPD and also to study leaf temperature, water potential, photosynthesis, and stomatal conductance responses. Also, our objective was also to study the structural changes in the stomatal traits when exposed to long term high VPD conditions and involvement in such responses. In the present study, 17 cotton genotypes were studied for their (TR) response to various VPD environments under well irrigated conditions. Out of 17, eight genotypes limited TR after approximately 2 kPa VPD and rest of them increased their TR with increased VPD. Five selected genotypes with different TR response to increasing VPD were further studied for gas exchange and stomatal properties. All genotypes, irrespective of exhibiting limited TR at high VPD, reduced stomatal conductance, photosynthesis and water potential at high VPD of 3.3 kPa. The genotypes with limited TR modified their stomatal traits mostly on the adaxial surface with frequent and small stomata under high VPD. The genotypes with limited TR also exhibited an increase in epidermal cell expansion and stomatal index at contrasting VPD gradients to effectively balance the liquid and vapor phase conductance to limit TR at high VPD.
Drought is the major limiting factor for cotton production, affecting growth, productivity, and fiber quality (Parida et al., 2007). Nearly one-third of the arable land around the globe suffers from water limitation, which is estimated to double by 2050 (Vicente-Serrano et al., 2012). Therefore, it is essential to produce or identify drought tolerant varieties to improve cotton yield in water-limited as well as well-watered environments. Nevertheless, only a few plant traits have actually been considered to improve drought resistance despite the many physiological traits involved in plant responses to drought. Several adaptive strategies have evolved in plants to cope with drought stress especially in environments where water deficit occurs frequently. One of the strategies is to limit transpiration rate (TR) under high vapor pressure deficit (VPD) environments. VPD, is the difference between the amount of moisture in the air and the moisture air can hold when it is saturated. It combines the effect of temperature and RH, has an effect on transpiration of the plants. High VPD, which usually occurs usually in the midday to end of the day, influences the water balance of the plant thus affecting photosynthesis and growth. Limiting transpiration in this situation will help the plant to conserve water for use later in the crop growing seasons when drought develops. Limited TR in response to high VPD has been observed in several crop species due to genotypic variability (Fletcher et al., 2007; Devi et al., 2010; Gholipoor et al., 2010; Kholová et al., 2010; Schoppach and Sadok, 2013) and found to increase crop yield by 75% in water-limited areas based on simulations in crop models (Sinclair et al., 2005, 2010; Sermons et al., 2012; Drake et al., 2013; Shekoofa et al., 2016).
High VPD under water limited conditions worsen the stress effects on plants either by increasing the TR (anisohydric) or by reducing carbon uptake (isohydric). Even though high VPD in general increases the diffusion process, plants regulate transpiration through stomatal closure or by means of other physiological responses (Aasamaa and Sõber, 2011). Stomata regulate water status of the plant conservatively by regulating the transpiration so that it matches to the atmosphere and the proportions of the soil–plant hydraulic system (Attia et al., 2015). Plants reduce stomatal conductance in response to high VPD between the leaf and atmosphere to maintain water potential as a feedforward mechanism (Bunce, 1997). By reducing stomatal conductance to water vapor, plants minimize water loss and maintain hydration of plant cells as VPD increases. There are many studies showing that the high VPD reduces stomatal conductance thereby affecting photosynthesis and growth (Wong, 1993; Bunce, 1997; Ottosen et al., 2002; Bunce, 2003; Ben-Asher et al., 2013). On the contrary, some studies have proposed that this high VPD impact on stomatal conductance may not influence photosynthesis (Carins Murphy et al., 2014; Shibuya et al., 2017).
Irrespective of the fact that the conductance of water flux depends on the atmospheric conditions surrounding the leaf, it is the leaf anatomy that controls the conductance of liquid water and water vapor (Carins Murphy et al., 2014). Stomata are gateway for gas exchange between the leaf and the surrounding environment. The equilibrium of water use and carbon assimilation in leaves basically relies upon water transport through the plant and the resultant coordination with the stomatal system (Kaiser, 2009; Brodribb and Jordan, 2011). The differences in stomatal anatomy and structure influence transpiration and photosynthetic activities demonstrating adaptations of the plants to various environmental variables (Larcher, 2003). Many studies have shown the modification of stomatal density or anatomical changes in stomata in response to various environmental factors such as CO2 (Driscoll et al., 2006), drought (Galmés et al., 2007; Xu and Zhou, 2008), light (Martins et al., 2014) and evaporative demand (Aliniaeifard and van Meeteren, 2014). Plasticity of stomatal, epidermal and vein density in response to variation in VPD has been observed in some studies (Torre et al., 2003; Carins Murphy et al., 2014). Irrespective of several studies on stomatal physiology, the complex mechanisms involved in the adjustment of stomatal traits such as size and aperture which regulate stomatal conductance are still poorly understood. In cotton, VPD is considered as a major promoting factor for TR (Isoda and Wang, 2002). Differences in TRs of cotton were observed in humid and arid conditions with increasing transpiration in arid environments (Isoda and Wang, 2002; Wang et al., 2004; Duursma et al., 2013; Conaty et al., 2014). Both TR and stomatal conductance showed strong response to varying VPD and photosynthesis was relatively insensitive (Duursma et al., 2013). However, no studies have investigated cotton genotype differences in TR response to various VPD environments. Hence, our objective in this study was to identify the genotypes limiting TR under high VPD and identify genotypic variation for this trait in cotton, which can be implemented in the development of drought tolerant varieties. Also, to understand the physiological responses involved to control TR, especially stomatal characters. In this study, therefore, we report observations on the response of transpiration of 17 cotton genotypes to differing VPD levels. Five genotypes with differences in their transpiration response to different VPD levels were selected out of 17 for a detailed study. The selected genotypes with differential VPD-TR response were subjected to high and low VPD levels to determine differences in the response patterns of the stomatal system. Our hypothesis is that the genotypes with differences in their TR response to VPD will show differences in their stomatal properties. The cultivars that limit TR at high VPD may modify stomatal and epidermal features when acclimated to high VPD in order to manage their TR. Conversely, the same difference would not be prominent in genotypes whose TR response to increasing VPD was linear and not restricted at higher VPD.
Materials and Methods
Total of three experiments was conducted to study the TR response to VPD and stomatal traits of cotton. Out of the three experiments, two experiments were to study the TR response of 17 cotton genotypes to VPD. The third experiment was conducted to study the stomatal traits of the selected genotype, after acclimatization of the genotypes for 10 days to low and high VPD environments. Gas exchange and water potential parameters were also observed in the same experiment. The experimental details from 1 to 3 were detailed in Table 1.
TABLE 1. Average temperature (°C) and RH (%) obtained in growth chambers and the average VPD calculated from temperature and RH data in all three experiments.
Screening of Genotypes for TR Response to VPD (Experiments 1 and 2)
Seventeen cotton genotypes (Table 1) were selected for this investigation and seeds of these genotypes were obtained from USDA-ARS, Lubbock, Texas. Before the measurement of VPD response, the plants were grown in pots in a greenhouse at the Beltsville Agriculture Research Center (BARC), ARS-USDA, and Beltsville, MD, United States. The pots (21 cm diameter and depth) were filled with Vermiculite (Miracle-Gro Lawn Products, Inc., Marysville, OH, United States) up to brim and planted with two seeds per pot. After about a week, each pot was thinned to one plant. The plants were grown under well-watered conditions, and the air temperature in the greenhouse was regulated at 27°C/21°C (day/night). Humidity in the greenhouse was not controlled, and the values measured were in between 60 and 70% in all three experiments during the daytime. Humidity in the greenhouse during night time was between 70 and 80%. The calculated VPD values in the greenhouse ranged from 1.00 to 1.45 kPa during daytime and in the night around 0.6 kPa. Grow lights were used from 6 to 9 AM as a supplement to natural day light received by the plants to bring the total light received by the plants to 1000 μmol m-2 s-1. Plants were also watered with full strength Hoagland nutrient solution every alternate day. Approximately after 30 days (6–7 true leaf stage), the plants were moved to growth chambers where different temperatures and humidity levels were set to maintain desired VPD levels as listed in Supplementary Table S1. The photosynthetic photon flux density (PPFD, 16 h d-1) of 800 μmol m-2 s-1 was maintained in each chamber. The humidity levels are maintained by using different air sources. To maintain low humidity powerex 40 hp oil-less rotary scroll air compressors (Powerex, Harrison, OH, United States) were used to dry the air. For the high humidity full cone misting nozzles (McMaster-Carr, Princeton, NJ, United States) were used and these nozzles provided small atomized spray drop lets for the humidification of the chamber to create desired high humidity levels. The maintained temperatures and humidity levels in the growth chambers in each experiment and the corresponding VPDs were listed in Supplementary Table S1. The plants were acclimated to growth chamber 1 day before transpiration measurements began. In all the growth chambers the temperature was set to the 27°C/21°C (day/night) as maintained in the greenhouse. The transpiration response was measured gravimetrically in all 17 genotypes in experiments 1 and 2 with six replicates in each experiment. In the evening of the day before the measurements began all plants were watered to saturation to avoid any water limitation stress.
In the morning of the day of measurements, the soil surface was covered with aluminum foil to avoid soil evaporation during the measurements of plant water loss. To avoid the effect of the growth chamber the plants were shuffled between the growth chambers before the measurements. The plants subjected to acclimation for an hour to each VPD and then the initial weight of the plant was recorded by weighing. The plants were exposed for an hour to each VPD level and then reweighed to measure the final weight. The difference between initial and final was calculated as transpiration of the plant. The same procedure was continued for all VPD levels and the plants were harvested. The leaf area of each plant was measured using leaf area meter (LI- 3100, Li-Cor, Lincoln, NE, United States). The TR of each plant was calculated as a ratio of transpiration per leaf area. Due to the limited space available in the growth chamber, 17 cotton genotypes were divided into three sets. First and second set with six genotypes in each set and third set including the rest of the five genotypes. The experiments were run with six replicates per genotype. In each set, a total of six chambers were used, each chamber including one replicate from each genotype.
Investigation of Stomatal Traits in Response to VPD (Experiment 3)
Leaf Gas Exchange
Leaf gas exchange parameters were measured only for the selected genotypes in the third experiment. A total of five genotypes were selected from experiments 1 and 2 based on their TR response to VPD. In the third experiment, the plants were grown in the greenhouse as mentioned above in experiments 1 and 2 until 23 days after emergence (vegetative stages 4–5 true leaf stage) and then moved to growth chambers where the experiment was conducted. Six VPD treatments 0.9, 1.3, 1.7, 1.9, 2.7, and 3.3 kPa were maintained in six chambers and a total of three replicates per genotype per treatment was included. The temperatures and humidity levels maintained in each chamber were listed in Supplementary Table S1. The photosynthetic photon flux density (PPFD, 16 h d-1) of 800 μmol m-2 s-1 was maintained in each chamber. The measurements were taken after 10 days of transferring the plants to the growth chamber. The plants were not rotated among the chambers as constant VPD was maintained in each chamber where there might be an effect of the growth chamber.
The gas exchange parameters were measured on a fully developed youngest leaf (7–8 true leaf stage), which had developed (after unfolding) after transferring the plants to growth chamber of particular VPD in each genotype. Net photosynthetic rate per unit leaf area (A), stomatal conductance (gs), and TR (E) were measured using a 0.26 chamber connected to a portable photosynthesis system (LI-6400, Li-Cor, Inc., Lincoln, NE, United States). The temperature of Li-Cor leaf chamber were closely matched to the growth chamber and CO2 level was set at 400 ppm. The flow rate of Li-Cor was adjusted between 300 and 500 μmol s-1 to obtain the close RH levels as in growth chamber (see Supplementary Table S1).
Leaf Water Potential (Ψ) and Leaf to Air Temperature Differences (LT-T°C)
Leaf water potential was measured at midday (11.30–12.30 h) using the youngest, fully expanded leaves (those for which gas exchange was measured) to different VPDs after taking the gas exchange and leaf temperature measurements. A Wescor psypro meter (Wescor, Pullman, WA, United States), eight channel water potential data logger was used for leaf water potential measurement, using three replicates for each treatment per genotype. The leaf temperature was estimated by using handy Infrared meter Flir TG165 (Flir systems, Inc., Nashua, NH, United States) by pointing to the whole leaf. The difference between leaves to air temperature was calculated as the difference between leaf temperature measured and the air temperature maintained in the growth chamber.
Stomatal traits were measured in the five genotypes in the plants exposed to low VPD of 0.9 kPa and high VPD of 3.3 kPa. The impressions were collected from the leaves selected for which gas exchange was measured. The impressions for stomatal traits were taken after measuring gas exchange, leaf temperature and leaf water potential. The impression approach was used to determine leaf stomatal properties. The adaxial and abaxial epidermis of the leaf was cleaned using a degreased cotton ball and then carefully layered with nail varnish in the mid-area between the central vein and the leaf edge, for approximately 1–2 min (Xu and Zhou, 2008). The thin film (approximately 10 mm × 10 mm) was peeled off from the leaf surface using clear tape and then immediately stick on to the microscopic slide. Numbers of stomata and epidermal cells for each film strip were counted using Nikon Eclipse E600 (Nikon Inc., Melville, NY, United States) equipped with a Nikon Digital Camera DXM1200. The area of the stomata and the epidermal cell was measured using Nikon NIS elements D. The length and width of the stomatal aperture were also measured using the analytical parameters of the NIS elements D. The stomatal density (SD) and epidermal cell density (ED) were expressed as the number of stomata per unit leaf area. The stomatal ratio was calculated as the ratio of stomata from adaxial leaf surface to the abaxial leaf surface. The leaf stomatal index was estimated using the formula [s/(e + s)] × 100, where s is stomatal density and e is epidermal cell density (Xu and Zhou, 2008).
Transpiration data of each genotype from the two experiments exposed to each VPD level were averaged (n = 12) and plotted against corresponding VPDs obtained. A two-segment linear regression (Prism 6.0, GraphPad Software Inc., San Diego, CA, United States) for TR vs. VPD was performed for all the genotypes. If the slopes were not significantly different (p < 0.05) using segmental linear regression, a simple linear regression was applied to all the data. The genotypes expressing the two-segment response, the regression analysis yielded the VPD breakpoint (X0) between the two linear segments as well as the slope of each segment. The standard error for the VPD breakpoint and the two slopes was an output of the regression analysis. The gaseous traits (n = 3) measured at different VPDs were also averaged for each genotype to perform two-segmental linear regression. If the regression was not obtained then a simple linear regression was used.
Water potential, LT-T and all stomatal traits were analyzed using ANOVA and the difference between the means of the genotypes for all the parameters was performed using Tukey’s Kramer test.
Genotypic Variation in TR Response to VPD
Seventeen genotypes of cotton were screened for their TR response to VPD in two repeated sets of experiments. In both the experiments, the distribution of temperature and relative humidity during the experiments was in the range of 23–32°C and 35–75%, respectively. The range of VPD treatments including all genotypes was between 0.74 ± 0.01 and 3.02 ± 0.03 kPa. The variation in VPD was mainly due to variation in the humidity levels established in the chambers as a result of differing air sources (see section “Materials and Methods”). The desired levels of VPD were obtained in each experiment (Table 1).
A clear distinction among genotypes in the response of TR to VPD was observed (Table 2). Some genotypes were well characterized by the two-segmental analysis with a break point (BP) (X0), and other genotypes exhibited a linear increase in TR for the range of VPD tested (Figure 1 and Table 2). Out of the 17 genotypes tested, eight were found to display a BP in TR increase as VPD was increased (Table 2). An example of two segmental linear regression was provided in the Figures 1A–C. The R2 for all data of each of these eight genotypes based on the two-segment regression ranged from 0.86 to 0.98. The value for the BP ± SE ranged from 1.58 ± 0.24 to 2.29 ± 0.15 kPa, with an average of 1.98 ± 0.09 kPa for the eight genotypes. The secondary slope of these eight genotypes above the BP ranged from -11.7 to + 14.7 mg H2O m-2 s-1 kPa-1. The distribution of slopes among eight genotypes were six positive and two negatives implying limited TR after BP (Table 2).
TABLE 2. Outputs from analysis of the two-segment linear regression and linear regression models for 17 cotton genotypes tested in experiments 1 and 2.
FIGURE 1. Transpiration rate (TR) of cotton genotype to vapor pressure deficit (VPD). TR (mg H2O m-2 s-1) response to increasing VPD (kPa) of the six cotton genotypes. (A) OL220 (B) DP555 BG RR and (C) LKT 57 represents genotypes with two segmental regression. The genotypes increased TR until VPD breakpoint (BP) X0 (kPa) and then reduced TR. The values of X0 indicated in each panel from (A–C). Panels from (D–F) represents genotypes 06-46-153P, CS 50 and Siokra L23, respectively, with linear TR response to increasing VPD.
The VPDvsTR data of nine genotypes did not fit the two-segment model, and hence, a single linear regression model was used (Figures 1D–F and Table 2). The linear regression generally fit these data well with the R2 ranging from 0.76 to 0.92. The genotypes transpired with an increase in VPD in the range of 24.6 to 36.1 mg H2O m-2 s-1 kPa-1. These slopes were similar to the slope 1 of the genotypes with BP which were in the range of 23.1–39.8 mg H2O m-2 s-1 kPa-1 (Table 2).
Gas Exchange, Water Potential and Leaf to Air Temperature Differences
Gas exchange, water potential, and leaf temperature were measured in the experiment 3 which was conducted with five selected genotypes from experiments 1 and 2 with and without limited TR. In the third experiment with VPD held between 0.9 and 3.3 kPa for 10 days, gas exchange responses to increased VPD showed differences among the genotypes with and without TR limitation at high VPD (Figure 2). There were slight differences between the VPD maintained in the growth chamber and the VPD obtained in the Li-Cor while taking the measurements. The obtained VPD in Li-Cor for among replicates in each growth chamber also slightly varied. The VPD maintained in each chamber and the range of VPD obtained with Li-Cor for three replicates in that chamber were listed in the Supplementary Table S1. The genotypes ‘OL220’ and ‘LKT 57,’ both with limited TR at high VPD, displayed a limitation in stomatal conductance (gs,), transpiration (E) and photosynthesis (A) after about 2.2–2.5 kPa (Figures 2A,B). This coincides with the whole TR response measured at different VPDs (0.74–3.02 kPa). The other three genotypes 06-46-153P, CS 50 and Siokra L23 with no limitation in its TR at high VPD also reduced their gs, E and A at high VPD of 3.3 kPa (Figures 2C–E). Nevertheless, the reduction in OL220 and LKT 57 appeared to be early than other three genotypes (Figure 2). The relation between stomatal conductance and internal CO2 was found significantly negative (P < 0.03) under high VPD conditions (VPD 2.65 and 3.26 kPa) (Supplementary Figure S1).
FIGURE 2. Gas exchange properties of the cotton genotypes with differences in their TR response to VPD. Stomatal conductance (μmol m-2 s-1), transpiration (mmol m-2 s-1) and photosynthesis (μmol m-2 s-1) of the genotypes with limited TR to high VPD (A) OL220 and (B) LKT 57. These two genotypes were well represented by two segmental linear regression to increase in VPD for their stomatal conductance, transpiration, and photosynthesis with BP at around 2.7 kPa. The genotypes in panels (C) 06-46-153 P, (D) CS 50, and (E) Siokra L23 showed the linear response for their gs, E, and A to the increase in VPD.
While there were slight differences (P < 0.05) among genotypes in the leaf water potential (Ψ) values measured, no conspicuous differences in Ψ among different VPD treatments were noticeable until 3.3 kPa (Figure 3A). The only significant differences were between the low VPD treatment 0.9 kPa and high VPD treatment 3.3 kPa (P < 0.01). All genotypes lowered their water potentials at high VPD (3.3 kPa). The leaf to air temperature differences appeared to be high with an increase in VPD and was less at high VPD (3.3 kPa) (Figure 3B). Low leaf temperatures than air temperatures were observed until 1.9 kPa in all genotypes. Apparent high leaf temperature than air temperature was found at high VPDs of 2.7 and 3.3 kPa in limited TR genotypes and at 3.3 kPa in genotypes without TR limitation.
FIGURE 3. Water potential (ψ) and leaf to air temperature differences (LT-T°C) of the five cotton genotypes to different VPD levels. (A) Water potential (Mpa) and (B) leaf to air temperature differences (°C) of the genotypes with limited TR and without limited TR to different VPD levels. The first two bars without pattern represent the genotypes OL220 and LKT 57 with limited TR. The other three bars with pattern represents the genotypes without limited TR. The genotypes were significantly different for water potential values (less than P < 0.05) at all VPD levels except at 3.3 kPa. The genotypes were significantly (less than P < 0.01) different for the leaf to air temperature differences at all VPD levels.
Stomatal and Epidermal Cell Frequency and Distribution
Stomatal traits were measured in the five selected genotypes in experiment 3 only in the plants exposed to long term VPDs of 0.9 and 3.3 kPa. The VPD influenced stomatal density (SD) on both adaxial and abaxial surface of the cotton leaf (Table 3). The SD was significantly different among the genotypes in plants acclimated to low VPD (0.90 ± 0.05) with an average stomatal density of 131 ± 10 stomata/mm2 on the adaxial surface and 250 ± 14 stomata/mm2 on the abaxial surface. Similarly, the SD was also varied across the genotypes in plants exposed to high VPD (3.3 ± 0.04) with an average SD of 227 ± 24 and 330 ± 28 on adaxial and abaxial leaf surfaces, respectively. While evaporative demand increased the number of stomata from low to high VPD on the adaxial surface in all genotypes except CS 50, the percentage of increase was found to be low in genotypes with limited TR (71 and 78%) compared to the genotypes with linear TR response (100 and 113) (Table 3). The percentage of increase in SD on abaxial was low in all genotypes compared to adaxial leaf surface, with no significant differences among genotypes. Even though there was a substantial increase in SD on adaxial surface, the genotypes, especially genotypes with limited TR response, reduced their stomatal area (-32.0 to -33.3%) more than the genotypes with linear TR response (-16.9 to -21.6%). All genotypes reduced their stomatal area with an increase in VPD. In contrast, all genotypes increased stomatal area with the increase in VPD on the abaxial surface with a low percentage of increase in genotypes with limited TR (Table 3).
TABLE 3. Adaxial and abaxial leaf surface stomatal density (mm-2), stomatal area (μm2), epidermal cell density (mm-2), epidermal cell area (μm2), stomatal aperture length (μm), and stomatal aperture width (μm) of the five cotton genotypes acclimated to low (0.9 kPa) and high (3.3 kPa) VPD environments.
All cultivars also increased their ratio of stomata from adaxial to abaxial surfaces from low to high VPD with a greater stomatal ratio in OL220 and LKT 57 than the other three genotypes whose TR response to VPD was linear (Figure 4A). The results also showed that the effect on the stomatal index (SI) was different on the adaxial and abaxial surfaces, when VPD significantly increased the SI by 13–51% on the adaxial surface (P < 0.001) but marginal or no increase on the abaxial surface (Figures 4B,C). Similarly, VPD also effected stomatal aperture length and width differently (Table 2). High VPD significantly reduced the aperture length in the genotypes OL220 and LKT 57 by approximately 45% on the adaxial surface. This was less than the genotypes with a linear TR response to increasing VPD (Table 3). The response was opposite on the abaxial surface with a slight increase in aperture length in high VPD plants than low VPD acclimated plants (Table 3). However, the elevated VPD had little effect on the aperture width with contrast responses on adaxial and abaxial surfaces. There was a slight increase in width on adaxial surface and reduction on abaxial surface in all genotypes with significantly less (P < 0.02) percentage of increase in genotypes with limited TR response on adaxial surface (Tables 3, 4).
FIGURE 4. Stomatal ratio and stomatal index of the cotton genotypes with differences in their TR response to VPD on adaxial and abaxial leaf surfaces. (A) Stomatal ratio (%), (B) adaxial stomatal index (%), and (C) abaxial stomatal index (%) of the cotton genotypes grown under low (0.9 kPa) and high (3.3 kPa) VPD environmental conditions. The first two bars in all three panels represent the genotypes with limited TR and other three bars are the data of the genotypes without limited TR. All genotypes were significantly different from each for all parameters under low VPD at P < 0.05 and in high VPD at P < 0.01.
TABLE 4. Percentage (%) of increase or decrease in stomatal traits (stomatal density, stomatal area, epidermal cell density, epidermal cell area, aperture length, and aperture width) from low (0.9 kPa) to high (3.3 kPa) VPD conditions of all the five cotton genotypes acclimated to low and high VPD environments.
In addition to stomatal properties, VPD also modified epidermal cell density (ED) and area. Some clear differences (P < 0.01) were observed on the adaxial surface in epidermal cell count from low VPD to high VPD and also among the genotypes with differences in their TR response to VPD. The genotypes OL220 and LKT 57 with limited TR to high VPD reduced their ED with -29 and -17%, respectively, when grown under high VPD conditions (Tables 3, 4). All other three genotypes increased their ED marginally (Tables 3, 4). There were no noticeable differences across low and high VPD treatments for their ED on abaxial surface in genotypes OL220 and Siokra L23. The other three genotypes improved their ED with increased VPD (Table 3). The response of epidermal cell area (EA) on adaxial leaf surface to VPD differences was contrast to the course of ED with increased EA in genotypes OL220 and LKT 57 (Table 3). The other three genotypes reduced their EA on both leaf surfaces (Tables 3, 4).
Correlation of Gas Exchange With Stomatal Properties
A correlation of stomatal ratio with stomatal conductance and transpiration was not observed in plants exposed to high and low VPD. Leaf stomatal index of both adaxial and abaxial surface was not correlated with stomatal conductance and hence with transpiration in plants grown in low VPD environment (Supplementary Figures S2A–D). Nevertheless, there was a significant negative correlation of leaf stomatal index of adaxial and abaxial leaf surface with stomatal conductance and transpiration (P < 0.01–0.03) (Supplementary Figures S2A–D) in high VPD treatment. However, the relationship between stomatal index and photosynthesis was scattered with no significant association.
This study displayed variation among 17 cotton genotypes tested for their transpiration response to the increasing VPD. Even though most of the studied genotypes had consistently increasing TR with increasing VPD, nearly half showed a distinct response by limiting their TR when VPD reached about 2 kPa. The BP, point where TR was limited observed in eight of the cotton genotypes studied is comparable to the VPD for the BP of genotypes in other species in which a limitation in TR to high VPD has been reported (Fletcher et al., 2007; Devi et al., 2010; Gholipoor et al., 2010; Kholová et al., 2010; Schoppach and Sadok, 2013).
The maximum stomatal conductance can be obtained under low evaporative demand conditions (Drake et al., 2013; Dow et al., 2014) which were observed in the present investigation. The genotypes OL220 and LKT 57 with limited TR at high VPD, associated with reduced stomatal conductance also constrained their photosynthetic rate (Figure 2). Often, this reduction in transpiration and photosynthesis may be counterbalanced by conserving the water with limited TR to high VPD for use later in the growing season (Sinclair et al., 2017). However, in this study, the genotypes with linear VPD- TR response 06-46-153P, CS 50 and Siokra L23 also reduced their stomatal conductance and photosynthesis at VPD of 3.3 kPa. Relatively, the reduction in genotypes with limited TR was at early VPD (2.7) than genotypes without TR limitation (3.3 kPa). The VPD in growing season on most of the days may possibly go beyond 2 kPa. At this time, the major consequence of decreasing stomatal conductance and subsequently the transpiration in the midday under high VPD is conservation of soil water. This would effectively maximize the transpiration efficiency by minimizing the percentage of the transpiration that usually arises during periods of high VPD (Sinclair et al., 2005). Reduction in photosynthesis followed abruptly by the decrease in stomatal conductance because of the fundamental link between water vapor and CO2 exchange in leaves (Tanner and Sinclair, 1983). Nevertheless, if there should be any incidence of late season drought, the genotypes with limited TR have the possibility of utilizing the water conserved by limiting the TR which would assist in improving yield than the genotypes without limited transpiration. In fact, in simulation studies of sorghum production in Australia and soybean in the United States showed a yield increase of about 75% of the seasons as a result of the limited transpiration at high VPD (Sinclair et al., 2005, 2010). However, the trait resulted in yield loss under wet conditions confirming the benefit of the trait under mild and severe drought conditions (Sinclair et al., 2010). Future simulation/crop modeling studies should be conducted in cotton to assess the utility of this trait and yield improvement in environments where mild or severe stress occurs.
Regulation of stomata to VPD is a process, by which plants adjusts their transpiration during daytime. Despite several studies, the mechanism for the stomatal response to VPD are poorly understood (McAdam and Brodribb, 2015). Some studies have suggested limited TR due to hydraulic limitation and the involvement of water channel proteins (Sinclair et al., 2017). Stomatal regulation of leaf water balance has been proposed to be controlled by active metabolic processes along with passive hydraulic process (McAdam and Brodribb, 2014). Recent studies had observed the expression of ABA synthesizing genes in guard cells and regulating stomatal responses to VPD (Waadt et al., 2014; McAdam and Brodribb, 2015). Apparently, these metabolic processes are prompted by low leaf water potential due to ABA (McAdam et al., 2016) and might result in a limitation in the TR. There were only marginal differences in leaf water potential with an increase in VPD until 3.3 kPa. The insensitivity to any decreases in leaf water potential is consistent with the observations in pea where leaves had to be exposed to more than -1 Mpa of pressure to induce a sufficient loss in turgor (McAdam and Brodribb, 2014). In this investigation, all genotypes lowered their Ψ at 3.3 kPa to the levels that were non-damaging. This is possibly due to feed forward response of the stomata to high evaporative demand (Bunce, 1997; Buckley, 2005). Reduction in stomatal conductance to high VPD caused in maintaining leaf Ψ to lower the damage to leaf.
A change in relative humidity and the concomitant change in TR often results in leaf temperature changes (Talbott et al., 2003). As VPD increased until 2.7 in limited TR genotypes and 3.3 in other genotypes, leaf temperature to air temperature (LT-T) ranged from a fraction of 1 to negative, probably because of increased transpiration cooling. However, at high VPD the difference between leaf to air temperature was positive, presumably due to the reduction in stomatal conductance and transpiration. The decrease in TR might have resulted in an increase of leaf temperature that caused the positive leaf to air differences (Pallas et al., 1967). An increase in temperature will have an impact on photorespiration causing temporary build up in internal CO2 Ci. Increased stomatal sensitivity to Ci was considered as a primary cause of reduced stomatal conductance under high VPD in wet soils (Bunce, 1997, 2003, 2007). The adjustments in the stomatal conductance are mainly to keep Ci/Ca at nearly a constant ratio, which results in a steady or reduced TR (Bunce, 1997, 2003, 2007), based on the adjustments. In the current study under high VPD conditions, where the TR was limited, a negative relation was observed between Ci and stomatal conductance (Supplementary Figure S1). As the data points in the study are small, the theory requires further examination including many data points.
The possibility of the contribution of stomatal traits in the expression of limited transpiration to high VPD was studied. Earlier studies have suggested the modification of stomatal properties by improving the density and shrinking the cell size in response to water deficit as an adaptation strategy (Martínez et al., 2007). In previous studies, in addition to other environmental factors, the effect of VPD on stomatal density and area was also investigated (Luomala et al., 2005; Carins Murphy et al., 2014). An increase in stomatal density and decrease in stomatal area from low to high VPD was observed especially on adaxial leaf surface in genotypes with limited transpiration response (Table 2). While stomatal density drives the partitioning of conductance among leaf sides, in general, stomatal size links to the gs with larger pores promoting gs (Fanourakis et al., 2015). Therefore, the decrease in stomatal size with the effect of high VPD corresponds to the decrease in gs and subsequently reduced TR in the current study.
Plants in response to environmental stress not only alter stomatal frequency (Xu and Zhou, 2008), but also change stomatal aperture size anatomically (Casson and Hetherington, 2010). All genotypes reduced their adaxial and abaxial stomatal length to high VPD but increased their width. The reduction in the stomatal aperture length to control transpiration at elevated temperature in maize leaf and different environmental variables in Arabidopsis was reported (Lomax et al., 2009; Zheng et al., 2013). Nevertheless, all genotyped reduced their stomatal aperture length to increase in VPD, with the high percentage of reduction in limited TR genotypes. Additionally, high VPD also modified stomatal ratio by increasing the ratio with no observed relation with any of the gas exchange parameters measured.
The previous study in maize had reported that the elevated VPD enhanced epidermal cell expansion in leaves but had little effect on the epidermal cell division (Tardieu et al., 2000). This suggests that the decrease of epidermal cell density in the high VPD than the low VPD treatment especially in TR limiting genotypes in the current study was mainly due to the greater expansion of the individual epidermal cells, rather than greater differentiation to stomatal cells under high VPD treatment. Our results also confirmed that the average size of the individual epidermal cells was significantly larger in the high VPD treatment than in the low VPD (Table 2) on the adaxial leaf surface. Correspondingly, some studies have identified that the plants adapt to various environmental variables by regulating their stomatal densities through modifying their epidermal cell density rather than stomatal numbers (Driscoll et al., 2006). The alteration of the stomatal index to high VPD also well related with the transpiration and stomatal conductance (Supplementary Figure S2).
The results of the present investigation indicated a range in the response of TR to VPD in cotton genotypes. Several of the genotypes exhibited a fairly limited TR at high VPD and while some demonstrated increasing TR. The limitation in TR at high VPD would contribute to water saving in the soil profile when water is available. This conserved water would be useful to the plants at the time of maturity and therefore for the yield improvement when late season drought develops. The results here can also be useful in the breeding program for the development of genotypes with optimized limited TR to utilize in water-limited environments. Field validation is required before these breeding strategies are applied to the genotypes limiting TR studied here. Future work needs the confirmation of the limited TR of the tested cotton genotypes in the current study in the field in a wide range of VPD. The genotypes also need to be tested for water conservation trait for the period of flowering through maturity.
The study here also indicates that it may be possible to adapt the transpiration response at elevated VPD by modifying their stomatal properties. The genotypes with limited exposed to the range of VPD (0.9–3.3 kPa) for longer time also retained their limited TR even after acclimatization. The genotypes with limited TR acclimated to low and high VPD reformed their stomatal traits differently. Stomatal and epidermal cells have a higher plasticity in response to a larger range of VPD, and these parameters are clearly associated with differences in transpiration response to high VPD in cotton genotypes resulting limited transpiration. The substantial changes in stomatal development might have facilitated to maintain transpiration homeostasis in cotton genotypes grown under contrasting VPD environments.
MD and VR conceived and designed the experiments. MD conducted the experiments, collected and analyzed the data. MD and VR wrote the manuscript.
Conflict of Interest Statement
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.
We thank Dr. Paxton Payton for providing the cotton seed. We also thank Jackson Fisher and Darryl Baxman for their help in maintaining the growth chambers.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01572/full#supplementary-material
Aasamaa, K., and Sõber, A. (2011). Stomatal sensitivities to changes in leaf water potential, air humidity, CO2 concentration and light intensity, and the effect of abscisic acid on the sensitivities in six temperate deciduous tree species. Environ. Exp. Bot. 71, 72–78. doi: 10.1016/j.envexpbot.2010.10.013
Aliniaeifard, S., and van Meeteren, U. (2014). Natural variation in stomatal response to closing stimuli among Arabidopsis thaliana accessions after exposure to low VPD as a tool to recognize the mechanism of disturbed stomatal functioning. J. Exp. Bot. 65, 6529–6542. doi: 10.1093/jxb/eru370
Attia, Z., Domec, J. C., Oren, R., Way, D. A., and Moshelion, M. (2015). Growth and physiological responses of isohydric and anisohydric poplars to drought. J. Exp. Bot. 66, 4373–4381. doi: 10.1093/jxb/erv195
Ben-Asher, J., Garcia, A., Flitcroft, I., and Hoogenboom, G. (2013). Effect of atmospheric water vapor on photosynthesis, transpiration and canopy conductance: a case study in corn. Plant Soil Environ. 59, 549–555. doi: 10.17221/346/2013-PSE
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
Bunce, J. A. (2003). Effects of water vapor pressure difference on leaf gas exchange in potato and sorghum at ambient and elevated carbon dioxide under field conditions. Field Crops Res. 82, 37–47. doi: 10.1016/S0378-4290(03)00004-2
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
Conaty, W. C., Mahan, J. R., Neilsen, J. E., and Constable, G. A. (2014). Vapour pressure deficit aids the interpretation of cotton canopy temperature response to water deficit. Funct. Plant Biol. 41, 535–546. doi: 10.1071/FP13223
Dow, G. J., Berry, J. A., and Bergmann, D. C. (2014). The physiological importance of developmental mechanisms that enforce proper stomatal spacing in Arabidopsis thaliana. New Phytol. 201, 1205–1217. doi: 10.1111/nph.12586
Driscoll, S. P., Prins, A., Olmos, E., Kunert, K. J., and Foyer, C. H. (2006). Specification of adaxial and abaxial stomata, epidermal structure and photosynthesis to CO2 enrichment in maize leaves. J. Exp. Bot. 57, 381–390. doi: 10.1093/jxb/erj030
Duursma, R. A., Payton, P., Bange, M. P., Broughton, K. J., Smith, R. A., Medlyn, B. E., et al. (2013). Near-optimal response of instantaneous transpiration efficiency to vapour pressure deficit, temperature and [CO2] in cotton (Gossypium hirsutum L.). Agric. For. Meteorol. 168(Suppl. C), 168–176. doi: 10.1016/j.agrformet.2012.09.005
Fanourakis, D., Giday, H., Milla, R., Pieruschka, R., Kjaer, K. H., Bolger, M., et al. (2015). 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
Fletcher, A. L., Sinclair, T. R., and Allen, L. H. (2007). Transpiration responses to vapor pressure deficit in well watered ‘slow-wilting’ and commercial soybean. Environ. Exp. Bot. 61, 145–151. doi: 10.1016/j.envexpbot.2007.05.004
Galmés, J., Flexas, J., Savé, R., and Medrano, H. (2007). Water relations and stomatal characteristics of Mediterranean plants with different growth forms and leaf habits: responses to water stress and recovery. Plant Soil 290, 139–155. doi: 10.1007/s11104-006-9148-6
Gholipoor, M., Prasad, P. V. V., Mutava, R. N., and Sinclair, T. R. (2010). Genetic variability of transpiration response to vapor pressure deficit among sorghum genotypes. Field Crops Res. 119, 85–90. doi: 10.1016/j.fcr.2010.06.018
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
Kholová, J., Hash, C. T., Kakkera, A., Koèová, M., and Vadez, V. (2010). Constitutive water-conserving mechanisms are correlated with the terminal drought tolerance of pearl millet [Pennisetum glaucum (L.) R. Br.]. J. Exp. Bot. 61, 369–377. doi: 10.1093/jxb/erp314
Lomax, B. H., Woodward, F. I., Leitch, I. J., Knight, C. A., and Lake, J. A. (2009). Genome size as a predictor of guard cell length in Arabidopsis thaliana is independent of environmental conditions. New Phytol. 181, 311–314. doi: 10.1111/j.1469-8137.2008.02700.x
Luomala, E. M., Laitinen, K., Sutinen, S., KellomÄKi, S., and Vapaavuori, E. (2005). Stomatal density, anatomy and nutrient concentrations of Scots pine needles are affected by elevated CO2 and temperature. Plant Cell Environ. 28, 733–749. doi: 10.1111/j.1365-3040.2005.01319.x
Martínez, J. P., Silva, H., Ledent, J. F., and Pinto, M. (2007). Effect of drought stress on the osmotic adjustment, cell wall elasticity and cell volume of six cultivars of common beans (Phaseolus vulgaris L.). Eur. J. Agron. 26, 30–38. doi: 10.1016/j.eja.2006.08.003
Martins, S. C., Galmés, J., Cavatte, P. C., Pereira, L. F., Ventrella, M. C., and DaMatta, F. M. (2014). Understanding the low photosynthetic rates of sun and shade coffee leaves: bridging the gap on the relative roles of hydraulic, diffusive and biochemical constraints to photosynthesis. PLoS One 9:e95571. doi: 10.1371/journal.pone.0095571
McAdam, S. A. M., Sussmilch, F. C., and Brodribb, T. J. (2016). Stomatal responses to vapour pressure deficit are regulated by high speed gene expression in angiosperms. Plant Cell Environ. 39, 485–491. doi: 10.1111/pce.12633
Ottosen, C. O., Mortensen, L. M., and Gislerod, H. R. (2002). Effect of relative air humidity on gas exchange, stomatal conductance and nutrient uptake in miniature potted roses. Eur. J. Hortic. Sci. 67, 143–147.
Pallas, J. E., Michel, B. E., and Harris, D. G. (1967). Photosynthesis, transpiration, leaf temperature, and stomatal activity of cotton plants under varying water potentials. Plant Physiol. 42, 76–88. doi: 10.1104/pp.42.1.76
Parida, A. K., Dagaonkar, V. S., Phalak, M. S., Umalkar, G. V., and Aurangabadkar, L. P. (2007). Alterations in photosynthetic pigments, protein and osmotic components in cotton genotypes subjected to short-term drought stress followed by recovery. Plant Biotechnol. Rep. 1, 37–48. doi: 10.1007/s11816-006-0004-1
Schoppach, R., and Sadok, W. (2013). Transpiration sensitivities to evaporative demand and leaf areas vary with night and day warming regimes among wheat genotypes. Funct. Plant Biol. 40, 708–718. doi: 10.1071/FP13028
Sermons, S. M., Seversike, T. M., Sinclair, T. R., Fiscus, E. L., and Rufty, T. W. (2012). Temperature influences the ability of tall fescue to control transpiration in response to atmospheric vapour pressure deficit. Funct. Plant Biol. 39, 979–986. doi: 10.1071/FP12172
Shekoofa, A., Sinclair, T. R., Messina, C. D., and Cooper, M. (2016). Variation among maize hybrids in response to high vapor pressure deficit at high temperatures. Crop Sci. 56, 392–396. doi: 10.2135/cropsci2015.02.0134
Shibuya, T., Kano, K., Endo, R., and Kitaya, Y. (2017). Photosynthetic properties and response to drought in cucumber seedlings acclimatized to different vapor-pressure-deficit levels. Hort. J. 86, 334–339. doi: 10.2503/hortj.MI-154
Sinclair, T. R., Devi, J., Shekoofa, A., Choudhary, S., Sadok, W., Vadez, V., et al. (2017). Limited-transpiration response to high vapor pressure deficit in crop species. Plant Sci. 260, 109–118. doi: 10.1016/j.plantsci.2017.04.007
Sinclair, T. R., Hammer, G. L., and van Oosterom, E. J. (2005). Potential yield and water-use efficiency benefits in sorghum from limited maximum transpiration rate. Funct. Plant Biol. 32, 945–952. doi: 10.1071/FP05047
Sinclair, T. R., Messina, C. D., Beatty, A., and Samples, M. (2010). Assessment across the United States of the benefits of altered soybean drought traits. Agron. J. 102, 475–482. doi: 10.2134/agronj2009.0195
Tanner, C., and Sinclair, T. (1983). “Efficient water use in crop production: research or re-search?” in Limitations to Efficient Water Use in Crop Production, eds H. M. Taylor, W. R. Jordan, and T. R. Sinclair (Madison, WI: ASA, CSSA and SSSA), 1–27.
Tardieu, F., Reymond, M., Hamard, P., Granier, C., and Muller, B. (2000). Spatial distributions of expansion rate, cell division rate and cell size in maize leaves: a synthesis of the effects of soil water status, evaporative demand and temperature. J. Exp. Bot. 51, 1505–1514. doi: 10.1093/jexbot/51.350.1505
Vicente-Serrano, S. M., Beguería, S., Lorenzo-Lacruz, J., Camarero, J. J., López-Moreno, J. I., Azorin-Molina, C., et al. (2012). Performance of drought indices for ecological, agricultural, and hydrological applications. Earth Interact. 16, 1–27. doi: 10.1175/2012ei000434.1
Waadt, R., Hitomi, K., Nishimura, N., Hitomi, C., Adams, S. R., Getzoff, E. D., et al. (2014). FRET-based reporters for the direct visualization of abscisic acid concentration changes and distribution in Arabidopsis. eLife 3:e01739. doi: 10.7554/eLife.01739
Keywords: cotton, gas exchange, stomata, transpiration, vapor pressure deficit, water potential
Citation: Devi MJ and Reddy VR (2018) Transpiration Response of Cotton to Vapor Pressure Deficit and Its Relationship With Stomatal Traits. Front. Plant Sci. 9:1572. doi: 10.3389/fpls.2018.01572
Received: 06 March 2018; Accepted: 09 October 2018;
Published: 30 October 2018.
Edited by:Robert Edward Sharwood, Australian National University, Australia
Reviewed by:Claudio Lovisolo, Università degli Studi di Torino, Italy
David Tissue, Western Sydney University, Australia
Copyright © 2018 Devi and Reddy. 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.